Predicting Hydrocarbon Primary Biodegradation in Soil and Sediment Systems Using System Parameterization and Machine Learning.

Technical complexity associated with biodegradation testing, particularly for substances of unknown or variable composition, complex reaction products, or biological materials (UVCB), necessitates the advancement of non-testing methods such as quantitative structure-property relationships (QSPRs). Models for describing the biodegradation of petroleum hydrocarbons (HCs) have been previously developed. A critical limitation of available models is their inability to capture the variability in biodegradation rates associated with variable test systems and environmental conditions. Recently, the Hydrocarbon Biodegradation System Integrated Model (HC-BioSIM) was developed to characterize the biodegradation of HCs in aquatic systems with the inclusion of key test system variables. The present study further expands the HC-BioSIM methodology to soil and sediment systems using a database of 2195 half-life (i.e., degradation time [DT]50) entries for HCs in soil and sediment. Relevance and reliability criteria were defined based on similarity to standard testing guidelines for biodegradation testing and applied to all entries in the database. The HC-BioSIM soil and sediment models significantly outperformed the existing biodegradation HC half-life (BioHCWin) and virtual evaluation of chemical properties and toxicities (VEGA) quantitative Mario Negri Institute for Pharmacological Research (IRFMN) models in soil and sediment. Average errors in predicted DT50s were reduced by up to 6.3- and 8.7-fold for soil and sediment, respectively. No significant bias as a function of HC class, carbon number, or test system parameters was observed. Model diagnostics demonstrated low variability in performance and high consistency of parameter usage/importance and rule structure, supporting the generalizability and stability of the models for application to external data sets. The HC-BioSIM provides improved accuracy of Persistence categorization, with correct classification rates of 83.9%, and 90.6% for soil and sediment, respectively, demonstrating a significant improvement over the existing BioHCWin (70.7% and 58.6%) and VEGA (59.5% and 18.5%) models. Environ Toxicol Chem 2024;43:1352-1363. © 2024 Concawe. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Comprehensive Two-Dimensional Gas Chromatography with Peak Tracking for Screening of Constituent Biodegradation in Petroleum UVCB Substances.

Petroleum substances, as archetypical UVCBs (substances of unknown or variable composition, complex reaction products, or biological substances), pose a challenge for chemical risk assessment as they contain hundreds to thousands of individual constituents. It is particularly challenging to determine the biodegradability of petroleum substances since each constituent behaves differently. Testing the whole substance provides an average biodegradation, but it would be effectively impossible to obtain all constituents and test them individually. To overcome this challenge, comprehensive two-dimensional gas chromatography (GC × GC) in combination with advanced data-handling algorithms was applied to track and calculate degradation half-times (DT50s) of individual constituents in two dispersed middle distillate gas oils in seawater. By tracking >1000 peaks (representing ∼53-54% of the total mass across the entire chromatographic area), known biodegradation patterns of oil constituents were confirmed and extended to include many hundreds not currently investigated by traditional one-dimensional GC methods. Approximately 95% of the total tracked peak mass biodegraded after 64 days. By tracking the microbial community evolution, a correlation between the presence of functional microbial communities and the observed progression of DT50s between chemical classes was demonstrated. This approach could be used to screen the persistence of GC × GC-amenable constituents of petroleum substance UVCBs.

Improved closed test setup for biodegradation testing of slightly volatile substances in water-sediment systems (OECD 308).

Standardized biodegradation testing methods, like the OECD 308 Aerobic and Anaerobic Transformation in Aquatic Sediment Systems, generate data on biodegradation required during environmental risk and hazard assessment of chemicals under different European and international regulations. However, difficulties arise when applying the OECD 308 guideline for testing hydrophobic volatile chemicals. Especially the use of a co-solvent (like acetone) as a measure to facilitate the application of the test chemical in combination with a closed setup to reduce losses due to volatilization tend to deplete/restrict the amount of oxygen in the test system. The result is a low oxygen or even anoxic water column in the water-sediment system. Thus, the degradation half-lives of the chemical generated from such tests are not directly comparable to the regulatory half-life values for Persistence assessment of the test chemical. The aim of this work was to further develop the closed setup to improve and maintain aerobic conditions in the water phase of the water-sediment systems for testing slightly volatile hydrophobic test chemicals. This improvement was attained by optimizing the test system geometry and agitation technique to maintain aerobic conditions in the water phase in a closed test setup, investigating appropriate co-solvent application strategy, and trialing the resulting test setup. This study shows that when using a closed test setup for OECD 308 tests, agitation of the water phase overlaying the sediment and the test item application using low co-solvent volume is critical for maintaining an aerobic water layer.

Using weight of evidence to assess degradation potential of UVCB hydrocarbon solvents.

Hydrocarbon solvents are a diverse group of petrochemical substances that are identified as unknown or variable composition, complex reaction products, or biological materials (UVCBs) and may contain tens of thousands of individual chemical constituents. As such, it is generally not possible to analytically resolve every chemical constituent in a hydrocarbon solvent. This, along with the low water solubility and/or high vapor pressure of constituents, precludes the use of many standardized tests designed to determine biodegradation in the environment (e.g., Organization for Economic Co-operation and Development [OECD] 309). A weight of evidence approach may be needed to reduce uncertainty to an acceptable level such that a determination on the biodegradation of the substance can be drawn. Based on the OECD 2019 weight of evidence guidance, we present a framework using various lines of evidence that can be used to evaluate the biodegradation of a UVCB solvent in a weight of evidence approach. The lines of evidence include whole substance testing, data on representative constituents, quantitative structure activity relationship (QSAR) models, and biological plausibility. Using these lines of evidence, "Hydrocarbon, C11-C14, normal alkane, isoalkane, cyclic, <2% aromatics" (EC# 926-141-6) was evaluated in a case study. Data from three whole substance tests, 43 constituents (representing 152 data points), three QSAR models and evidence of microbial degradation pathways were evaluated. Based on the available data, it is concluded that the solvent for the case study is not expected to persist in the environment. This framework sets out a real-world example of how the weight of evidence can be used to evaluate hydrocarbon solvents. While focused on persistence, similar approaches can be used to evaluate other endpoints such as bioaccumulation and toxicity. Integr Environ Assess Manag 2023;19:1120-1130. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Persistence and Mobility (defined as organic­carbon partitioning) do not correlate to the detection of substances found in surface and groundwater: Criticism of the regulatory concept of Persistent and mobile substances.

The Chemical Strategy for Sustainability (CSS) includes actions to ensure the protection of drinking water resources from chemical pollution. To proactively identify potential pollutants, the German Environment Agency (UBA) proposed the Persistent and Mobile (PM) concept according to which Persistence (criteria of REACH Annex XIII) and Mobility (log Koc < 4) would be proxies for a substance's degradation potential and transport velocity, two processes believed to drive the potential for contamination of surface and groundwater as drinking water sources. Two studies identified hundreds of PM substances while three subsequent studies have selected some of these substances for monitoring in surface, ground- and/or drinking water to support the concept. In the present work, the Persistence of the aforementioned substances was reassessed based on all experimental data publicly available. Depending on the exact study examined, it was found that 15 % to 40 % of the substances were erroneously concluded as P. The reinterpretation of the data indicates that a PM substance does not have a higher likelihood to be detected in surface or groundwater than a non-PM substance. In addition, the PM properties do not have any influence on the level of contamination. Twenty-six to 75 % of the substances selected because they were identified as PM were not found in surface or ground water despite being selected for their high emission pattern. Regulations based primarily on the PM concept, like the CLP and possibly REACH and UN-GHS, are unlikely to appropriately identify substances of concern for drinking water sources. It is more likely that chemical presence in surface and groundwater is driven by emission patterns or local factors. The development of specific exposure models would better contribute to the protection of drinking water resources and consumers.

Predicting Primary Biodegradation of Petroleum Hydrocarbons in Aquatic Systems: Integrating System and Molecular Structure Parameters using a Novel Machine-Learning Framework.

Quantitative structure-property relationship (QSPR) models for predicting primary biodegradation of petroleum hydrocarbons have been previously developed. These models use experimental data generated under widely varied conditions, the effects of which are not captured adequately within model formalisms. As a result, they exhibit variable predictive performance and are unable to incorporate the role of study design and test conditions on the assessment of environmental persistence. To address these limitations, a novel machine-learning System-Integrated Model (HC-BioSIM) is presented, which integrates chemical structure and test system variability, leading to improved prediction of primary disappearance time (DT50) values for petroleum hydrocarbons in fresh and marine water. An expanded, highly curated database of 728 experimental DT50 values (181 unique hydrocarbon structures compiled from 13 primary sources) was used to develop and validate a supervised model tree machine-learning model. Using relatively few parameters (6 system and 25 structural parameters), the model demonstrated significant improvement in predictive performance (root mean square error = 0.26, R2 = 0.67) over existing QSPR models. The model also demonstrated improved accuracy of persistence (P) categorization (i.e., "Not P/P/vP"), with an accuracy of 96.8%, and false-positive and -negative categorization rates of 0.4% and 2.7%, respectively. This significant improvement in DT50 prediction, and subsequent persistence categorization, validates the need for models that integrate experimental design and environmental system parameters into biodegradation and persistence assessment. Environ Toxicol Chem 2022;41:1359-1369. © 2022 ExxonMobil Biomedical Sciences, Inc. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Selection of Representative Constituents for Unknown, Variable, Complex, or Biological Origin Substance Assessment Based on Hierarchical Clustering.

Many of the newly produced and registered substances are complex mixtures or substances of unknown or variable composition, complex reaction products, and biological materials (UVCBs). The latter often consist of a large number of constituents, some of them difficult-to-identify constituents, which complicates their (eco)toxicological assessment. In the present study, through a series of examples, different scenarios for selection of representatives via hierarchical clustering of UVCB constituents are exemplified. Hierarchical clustering allows grouping of the individual chemicals into small sets, where the constituents are similar to each other with respect to more than one criterion. To this end, various similarity criteria and approaches for selection of representatives are developed and analyzed. Two types of selection are addressed: (1) selection of the most "conservative" constituents, which could be also used to support prioritization of UVCBs for evaluation, and (2) obtaining of a small set of chemical representatives that covers the structural and metabolic diversity of the whole target UVCBs or a mixture that can then be evaluated for their environmental and (eco)toxicological properties. The first step is to generate all plausible UVCB or mixture constituents. It was found that the appropriate approach for selecting representative constituents depends on the target endpoint and physicochemical parameters affecting the endpoint of interest. Environ Toxicol Chem 2021;40:3205-3218. © 2021 SETAC.

Acritical review and weight of evidence approach for assessing the bioaccumulation of phenanthrene in aquatic environments.

Bioaccumulation (B) assessment is challenging because there are various B-metrics from laboratory and field studies, multiple criteria and thresholds for classifying bioaccumulative (B), very bioaccumulative (vB), and not bioaccumulative (nB) chemicals, as well as inherent variability and uncertainty in the data. These challenges can be met using a weight of evidence (WoE) approach. The Bioaccumulation Assessment Tool (BAT) provides a transparent WoE assessment framework that follows Organisation for Economic Co-operation and Development (OECD) principles for performing a WoE analysis. The BAT guides an evaluator through the process of data collection, generation, evaluation, and integration of various lines of evidence (LoE) (i.e., B-metrics) to inform decision-making. Phenanthrene (PHE) is a naturally occurring chemical for which extensive B and toxicokinetics data are available. A B assessment for PHE using the BAT is described that includes a critical evaluation of 74 measured in vivo LoE for fish and invertebrate species from laboratory and field studies. The number of LoE are reasonably well balanced across taxa (i.e., fish and invertebrates) and the different B-metrics. Additionally, in silico and in vitro biotransformation rate estimates and corresponding model-predicted B-metrics are included as corroborating evidence. Application of the BAT provides a consistent, coherent, and scientifically defensible WoE evaluation to conclude that PHE is not bioaccumulative (nB) because the overwhelming majority of the bioconcentration, bioaccumulation, and biomagnification metrics for both fish and invertebrates are below regulatory thresholds. An analysis of the relevant data using fugacity ratios is also provided, showing that PHE does not biomagnify in aquatic food webs. The critical review identifies recommendations to increase the consistency of B assessments, such as improved standardization of B testing guidelines, data reporting requirements for invertebrate studies, and consideration of temperature and salinity effects on certain B-metrics. Integr Environ Assess Manag 2021;17:911-925. © 2021 Concawe. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Assessing toxicity of hydrophobic aliphatic and monoaromatic hydrocarbons at the solubility limit using novel dosing methods.

Reliable delineation of aquatic toxicity cut-offs for poorly soluble hydrocarbons is lacking. In this study, vapor and passive dosing methods were applied in limit tests with algae and daphnids to evaluate the presence or absence of chronic effects at exposures corresponding to the water solubility for representative hydrocarbons from five structural classes: branched alkanes, mono, di, and polynaphthenic (cyclic) alkanes and monoaromatic naphthenic hydrocarbons (MANHs). Algal growth rate and daphnid immobilization, growth and reproduction served as the chronic endpoints investigated. Results indicated that the dosing methods applied were effective for maintaining mean measured exposure concentrations within a factor of two or higher of the measured water solubility of the substances investigated. Chronic effects were not observed for hydrocarbons with an aqueous solubility below approximately 5 µg/L. This solubility cut-off corresponds to structures consisting of 13-14 carbons for branched and cyclic alkanes and 16-18 carbons for MANHs. These data support reliable hazard and risk evaluation of hydrocarbon classes that comprise petroleum substances and the methods described have broad applicability for establishing empirical solubility cut-offs for other classes of hydrophobic substances. Future work is needed to understand the role of biotransformation on the observed presence or absence of toxicity in chronic tests.

Is the Arrhenius-correction of biodegradation rates, as recommended through REACH guidance, fit for environmentally relevant conditions? An example from petroleum biodegradation in environmental systems.

Biodegradation is a major determinant of chemical persistence in the environment and an important consideration for PBT and environmental risk assessments. It is influenced by several environmental factors including temperature and microbial community structure. According to REACH guidance, a temperature correction based on the Arrhenius equation is recommended for chemical persistence data not performed at the recommended EU mean surface water temperature. Such corrections, however, can lead to overly conservative P/vP assessments. In this paper, the relevance of this temperature correction is assessed for petroleum hydrocarbons, using measured surface water (marine and freshwater) degradation half-time (DT50) and degradation half-life (HL) data compiled from relevant literature. Stringent screening criteria were used to specifically select data from biodegradation tests containing indigenous microbes and conducted at temperatures close to their ambient sampling temperature. As a result, ten independent studies were identified, with 993 data points covering 326 hydrocarbon constituents. These data were derived from tests conducted with natural seawater, or freshwater, at temperatures ranging from 5 to 21 °C. Regressions were performed on the full hydrocarbon dataset and on several individual hydrocarbons. The results were compared to the trend as predicted by the Arrhenius equation and using the activation energy (Ea) as recommend in the REACH Guidance. The comparison shows that the correction recommended in REACH Guidance over predicts the effect of temperature on hydrocarbon biodegradation. These results contrast with temperature manipulated inocula where the test temperature is different from the ambient sampling temperature. In these manipulated systems, the effect of temperature follows the Arrhenius equation more closely. In addition, a more striking effect of temperature on the lag phase was observed with longer lag phases more apparent at lower temperatures. This indicates that the effect of temperature may indeed be even lower when considering hydrocarbon biodegradation without the initial lag phase.

Tolerance and Excretion of the Mycotoxins Aflatoxin B1, Zearalenone, Deoxynivalenol, and Ochratoxin A by Alphitobius diaperinus and Hermetia illucens from Contaminated Substrates.

This study aimed to investigate the potential accumulation of mycotoxins in the lesser mealworm (Alphitobius diaperinus, LMW) and black soldier fly (Hermetia illucens, BSF) larvae. Feed was spiked with aflatoxin B1, deoxynivalenol (DON), ochratoxin A or zearalenone, and as a mixture of mycotoxins, to concentrations of 1, 10, and 25 times the maximum limits set by the European Commission for complete feed. This maximum limit is 0.02 mg/kg for aflatoxin B1, 5 mg/kg for DON, 0.5 mg/kg for zearalenone and 0.1 mg/kg for ochratoxin A. The mycotoxins and some of their metabolites were analysed in the larvae and residual material using a validated and accredited LC-MS/MS-based method. Metabolites considered were aflatoxicol, aflatoxin P1, aflatoxin Q1, and aflatoxin M1, 3-acetyl-DON, 15-acetyl-DON and DON-3-glycoside, and α- and ß-zearalenol. No differences were observed between larvae reared on mycotoxins individually or as a mixture with regards to both larvae development and mycotoxin accumulation/excretion. None of the mycotoxins accumulated in the larvae and were only detected in BSF larvae several orders of magnitude lower than the concentration in feed. Mass balance calculations showed that BSF and LMW larvae metabolized the four mycotoxins to different extents. Metabolites accounted for minimal amounts of the mass balance, except for zearalenone metabolites in the BSF treatments, which accounted for an average maximum of 86% of the overall mass balance. Both insect species showed to excrete or metabolize the four mycotoxins present in their feed. Hence, safe limits for these mycotoxins in substrates to be used for these two insect species possibly could be higher than for production animals. However, additional analytical and toxicological research to fully understand the safe limits of mycotoxins in insect feed, and thus the safety of the insects, is required.

Long-Term Persistence of Pesticides and TPs in Archived Agricultural Soil Samples and Comparison with Pesticide Application.

For polar and more degradable pesticides, not many data on long-term persistence in soil under field conditions and real application practices exist. To assess the persistence of pesticides in soil, a multiple-compound screening method (log Kow 1.7-5.5) was developed based on pressurized liquid extraction, QuEChERS and LC-HRMS. The method was applied to study 80 polar pesticides and >90 transformation products (TPs) in archived topsoil samples from the Swiss Soil Monitoring Network (NABO) from 1995 to 2008 with known pesticide application patterns. The results reveal large variations between crop type and field sites. For the majority of the sites 10-15 pesticides were identified with a detection rate of 45% at concentrations between 1 and 330 µg/kgdw in soil. Furthermore, TPs were detected in 47% of the cases where the "parent-compound" was applied. Overall, residues of about 80% of all applied pesticides could be detected with half of these found as TPs with a persistence of more than a decade.

Effects of Milk Yield, Feed Composition, and Feed Contamination with Aflatoxin B1 on the Aflatoxin M1 Concentration in Dairy Cows' Milk Investigated Using Monte Carlo Simulation Modelling.

This study investigated the presence of aflatoxin M1 (AfM1) in dairy cows' milk, given predefined scenarios for milk production, compound feed (CF) contamination with aflatoxin B1 (AfB1), and inclusion rates of ingredients, using Monte Carlo simulation modelling. The model simulated a typical dairy farm in the Netherlands. Six different scenarios were considered, based on two lactation and three CF composition scenarios. AfB1 contamination of the CF was based on results from the Dutch national monitoring programme for AfB1 in feed materials from 2000 until 2010. Monitoring data from feed materials used in CF production for dairy cattle in the Netherlands were used. Additionally, AfB1 contamination data from an incident in maize in 2013 were used. In each scenario, five different transfer equations of AfB1 from feed to AfM1 in the milk were used, and 1000 iterations were run for each scenario. The results showed that under these six scenarios, the weekly farm concentration of AfM1 in milk was above the EC threshold in less than 1% of the iterations, with all five transfer equations considered. However, this increased substantially in weeks when concentrations from the contaminated maize batch were included, and up to 28.5% of the iterations exceeded the EC threshold. It was also observed that an increase in the milk production had a minimal effect on the exceedance of the AfM1 threshold due to an apparent dilution effect. Feeding regimes, including the composition of CF and feeding roughages of dairy cows, should be carefully considered based on the potential AfM1 contamination of the farm's milk.

Local organochlorine pesticide concentrations in soil put into a global perspective.

In this work, agricultural and background soil concentrations of p,p'-DDT, p,p'-DDE, HCB, α-, ß- and γ-HCH from 1993 to 2012 were collected from 73 peer-reviewed publications, and analysed statistically. For the period 2003-2012 and for all chemicals, the mean concentration in agricultural soil is significantly higher than the concentration in background soil. In addition to the statistical analysis, concentrations of p,p'-DDT and α-HCH in soils were calculated with a global environmental fate and transport model. A decrease in the mean soil concentration from the first decade to the second was observed with the model, but this decrease is not visible in the measured concentrations, which could result from ongoing use of p,p'-DDT and α-HCH Furthermore, modelled background soil concentrations are generally lower than measurements. This implies that background soil may have received p,p'-DDT and α-HCH through additional routes not described by the model such as spray drift.

Historical emissions of octachlorodibenzodioxin in a watershed in Queensland, Australia: estimation from field data and an environmental fate model.

An octachlorodibenzodioxin (OCDD)-dominated contamination is present along the coast of Queensland, Australia. Several findings indicate that this contamination originates from pesticide use, although due to limited information on OCDD levels in the pesticides used, estimating past and current emissions of OCDD solely from pesticide use data is unfeasible. We used all the qualitative and quantitative information available on OCDD in pesticides together with a previously validated chemical fate model for a catchment in the Queensland Wet Tropics to back-calculate the emissions of OCDD from measured soil concentrations. We estimate that under different emission scenarios an average of 2,500 kg of OCDD was emitted within the modelled 1,685 km2 (Tully river) catchment between 1950 and 2010. Because this catchment represents only approximately 0.85% of the whole coast of Queensland under a similar contamination, the total amount of OCDD released in this region is considerably larger. For all emission scenarios, we could show that the OCDD currently present in agricultural soil is a result of historical emissions, and current-day emissions are less important in comparison to past emissions. Overall 18% was lost by degradation and 62% was buried below the agricultural surface soil, as a result of facilitated transport.

Sorption kinetics and equilibrium of the herbicide diuron to carbon nanotubes or soot in absence and presence of algae.

Carbon nanotubes (CNT) are strong sorbents for organic micropollutants, but changing environmental conditions may alter the distribution and bioavailability of the sorbed substances. Therefore, we investigated the effect of green algae (Chlorella vulgaris) on sorption of a model pollutant (diuron, synonyms: 3-(3,4-Dichlorophenyl)-1,1-dimethylurea, DCMU) to CNT (multi-walled purified, industrial grade, pristine, and oxidized; reference material: Diesel soot). In absence of algae, diuron sorption to CNT was fast, strong, and nonlinear (Freundlich coefficients: 10(5.79)-10(6.24) µg/kgCNT·(µg/L)(-n) and 0.62-0.70 for KF and n, respectively). Adding algae to equilibrated diuron-CNT mixtures led to 15-20% (median) diuron re-dissolution. The relatively high amorphous carbon content slowed down ad-/desorption to/from the high energy sorption sites for both industrial grade CNT and soot. The results suggest that diuron binds readily, but - particularly in presence of algae - partially reversibly to CNT, which is of relevance for environmental exposure and risk assessment.

Diuron sorbed to carbon nanotubes exhibits enhanced toxicity to Chlorella vulgaris.

Carbon nanotubes (CNT) are more and more likely to be present in the environment, where they will associate with organic micropollutants due to strong sorption. The toxic effects of these CNT-micropollutant mixtures on aquatic organisms are poorly characterized. Here, we systematically quantified the effects of the herbicide diuron on the photosynthetic activity of the green alga Chlorella vulgaris in presence of different multiwalled CNT (industrial, purified, pristine, and oxidized) or soot. The presence of carbonaceous nanoparticles reduced the adverse effect of diuron maximally by <78% (industrial CNT) and <34% (soot) at 10.0 mg CNT/L, 5.0 mg soot/L, and diuron concentrations in the range 0.73-2990 µg/L. However, taking into account the measured dissolved instead of the nominal diuron concentration, the toxic effect of diuron was equal to or stronger in the presence of CNT by a factor of up to 5. Sorbed diuron consequently remained partially bioavailable. The most pronounced increase in toxicity occurred after a 24 h exposure of algae and CNT. All results point to locally elevated exposure concentration (LEEC) in the proximity of algal cells associated with CNT as the cause for the increase in diuron toxicity.

Describing the environmental fate of diuron in a tropical river catchment.

The use of the herbicide diuron on sugarcane fields along the river catchments of the Great Barrier Reef (GBR) in Australia is an issue of concern due to high levels of diuron reported in the GBR lagoon, and has recently led to a restriction on the use of diuron during the 2011/12 wet season. An important question in this context is how much diuron is mobilised from the agricultural area by strong rainfall and floods in the wet season and transferred to the GBR lagoon. We have set up a multimedia chemical fate model for a tropical catchment to describe the fate of diuron within the Tully River catchment, Queensland, Australia. The model includes highly variable rainfall based on meteorological data from the Tully River catchment and a flood water compartment on top of the agricultural soil that is present during times for which floods were reported. The model is driven by diuron application data estimated for the Tully River catchment and is solved for time-dependent diuron concentrations in agricultural soil and seawater. Model results show that on average 25% of the diuron applied every year is transferred to the GBR lagoon with rainwater and flood water runoff. Diuron concentrations estimated for the seawater range from 0.1 ng/L to 12 ng/L and are in good agreement with concentrations measured in the GBR lagoon. The uncertainty of the diuron concentrations estimated for seawater is approximately a factor of two and mainly derives from uncertainty in the diuron degradation half-life in soil, properties of the soil compartment such as organic matter content, and the speed of the seawater current removing diuron dissolved in seawater from the seawater compartment of the model.