Extended sorption partitioning models for pesticide leaching risk assessments: Can we improve upon the koc concept?

Assessing the suitability of leachability-based screening levels for per- and polyfluoroalkyl substances (PFAS) risk assessment

In recent years, soil screening levels have been adopted by regulatory agencies for certain per- and polyfluoroalkyl substances (PFAS) to assess the risk of groundwater contamination through leaching. These soil screening levels, determined using an established equilibrium-based partitioning equation, have high variability among regulatory groups largely attributed to the diverse reported partitioning coefficients in the literature. This variability between reported partitioning coefficients, and subsequently soil screening levels, is due to the complex leaching behavior of PFAS not being predicted well by the standard equilibrium-based model. This has led one regulatory group to require batch leaching to assess risk rather than setting default soil screening levels based on partitioning equations. In this work, we conducted leaching experiments on five field-sampled soils impacted by aqueous film-forming foams (AFFF), following Leaching Environmental Assessment Framework (LEAF) Method 1316 and compared the results to expected leaching utilizing an equilibrium-based partitioning equation commonly employed by regulatory agencies to establish soil screening levels. Our analysis found among the six PFAS detected in the soils, which have regulatory leaching thresholds established, the partitioning values assumed by the U.S. EPA exhibited the highest accuracy in predicting leachate concentrations. These partitioning values predicted actual leaching within a ± 20 % margin of error for approximately 50 % of sample points, highlighting limitations in relying solely on equilibrium-based partitioning values as predictors of leaching behavior. This discrepancy between predicted and actual leaching has implications for site managers and regulatory entities overseeing PFAS-contaminated sites, suggesting that soil screening level determinations for PFAS might need to be revised to account for the unique transport characteristics of PFAS.

Investigating residue dissolution of insensitive high explosives in two sandy soil types: A predictive modelling approach

The demand for munitions that are less likely to detonate accidentally has led to an increased use of Insensitive High Explosives (IHE), which contain substances like 2,4-dinitroanisole (DNAN) and 5-nitro-1,2,4-triazol-3-one (NTO). These substances have different properties compared to traditional explosives, and their potential environmental impact is not well understood. When these explosives are used in live-fire training exercises, their residues end up in the soil. It is important to determine how these residues dissolve and enter the soil. This study aimed to experimentally measure the rate at which an IHE formulation dissolves when exposed to rainwater with pH levels of 5.0 and 6.5, and to simulate how these residues dissolve and move through two different soil types. The dissolution rates were determined by conducting experiments in which IHE particles (30-60 mg) were exposed to water with varying pH levels and temperatures. The results showed that the dissolution rate of NTO did not vary with pH, while the dissolution rate of DNAN and RDX decreased with decreasing pH. Specifically, the dissolution rate of DNAN decreased from 18 ± 40 μg min-1 at pH 6.5 to 6 ± 4 μg min-1 at pH 5.0, while the dissolution rate of RDX decreased from 8 ± 4 to 3 ± 1 μg min-1. These findings were used to develop a stochastic model that successfully simulated the concentration of IHE in the leachate from soil columns over time. A sensitivity analysis revealed that while dissolution rates determined the amount of mass entering the soil, they did not significantly regulate the amount of mass that migrated through the soil and leached out of the soil columns.

Mobility in the context of exposure-based assessment of chemicals for drinking water resource protection

In order to protect European Union (EU) drinking water resources from chemical contamination, criteria for identifying persistent, mobile, and toxic (PMT) chemicals and very persistent and very mobile (vPvM) chemicals under the EU REACH Regulation were proposed by the German Environment Agency (Umweltbundesamt-UBA). Additionally, new hazard classes for PMT and vPvM substances in the revised EU classification, labeling, and packaging (CLP Regulation) are intended. Therefore, a reliable approach in the identification of potential drinking water resource contaminants is needed. The scientific basis of the property-based PMT/vPvM criteria, focusing on mobility, which dictates the migration of chemical drinking water sources, was evaluated, and a critical analysis of the deviation of sorption metrics from simple behavior was carried out. Based on our evaluation, a K_oc may be used for nonionic substances on a screening level only, requiring a higher tier assessment. It is considered inappropriate for hydrophilic and ionizable chemicals, particularly for soils with low organic carbon contents. The nonextractable residue formation is complex and not well understood but remains significant in limiting the mobility of chemicals through soils and sediments. In order to inform the EU commission's work on the introduction of new hazard classes for PMT and vPvM substances into the European legislation, the derivation of a tiered approach is proposed, which utilizes the weight of evidence available, with adoption of appropriate higher tier models commensurate with the nature of the substance and the data available. Integr Environ Assess Manag 2023;19:775-791. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Pesticides groundwater modelling relies on input data characterised by a high intrinsic variability: Is the resulting risk for groundwater credible?

In the framework of Regulation (EC) 1107/2009, concerning the placing of plant protection products (PPP) on the market, FOCUS models are used to predict active substances concentration in groundwater. The predicted environmental concentration in groundwater (PEC_GW) are influenced by active substance specific parameters, namely DT₅₀, K_OM and Freundlich coefficient (1/n), whose minimal variation in certain combinations of intervals significantly affects PEC_GW output. Considering that minimal variation are intrinsic in all laboratory studies, this approach may lead to not acceptable variations in the results for regulatory purposes. In the present article, PEC_GW were calculated for all maize crop scenarios, using 808 dummy active substances with different combinations of DT₅₀, K_OM and 1/n values, in order to quantify the influence of each single parameter on the final result of PEARL and PELMO models. The results obtained were used to create a classification system for the input parameters K_OM and DT₅₀ in order to minimise the input uncertainty effects. Even if this approach is scientifically viable yet, due to its conservative nature, it cannot be considered suitable in the regulatory framework, where acceptability of an active substance is strictly related to the limit value of 0.1 μg/L. Nevertheless, this classification system could represent an important screening or preliminary assessment to plan pesticide monitoring programmes. Based on the results of this analysis, it is believed that the assessment of pesticide leaching into groundwater should be revised to take into account this variability. Considering that both PEARL and PELMO FOCUS models deal with interaction between a chemical and a complex system like soil and weather, the selection of input data cannot pretend to rely on single specific number. Considering that intrinsic uncertainty cannot be eliminated from experimental work, a revision of the criteria used to identify the proper input data and a thorough revision of the actual groundwater modelling is recommended.

Mechanisms of Tebuconazole Adsorption in Profiles of Mineral Soils

The study attempted to identify the soil components and the principal adsorption mechanisms that bind tebuconazole in mineral soils. The K_F values of the Freundlich isotherm determined in 18 soils from six soil profiles in batch experiments after 96 h of shaking ranged from 1.11 to 16.85 μg^1-1/n (mL)^1/n g^-1, and the exponent 1/n values from 0.74 to 1.04. The adsorption of tebuconazole was inversely correlated with the soil pH. Both neutral and protonated forms of this organic base were adsorbed mainly on the fraction of humins. The adsorption of the protonated form increased in the presence of hydrogen cations adsorbed in the soil sorption sites. Fourier transform infrared spectroscopy coupled with the molecular modeling studies and partial least squares regression analysis indicated that the tebuconazole molecule is bound in the organic matter through the formation of hydrogen bonds as well as hydrophobic and π-π interactions. Ion exchange was one of the adsorption mechanisms of the protonated form of this fungicide. The created mathematical model, assuming that both forms of tebuconazole are adsorbed on the organic matter and adsorption of the protonated form is affected by the potential acidity, described its adsorption in soils well.

Spatial distribution of sorption and desorption process of 14C-radiolabelled hexazinone and tebuthiuron in tropical soil

The understanding of the interaction between soil physicochemical attributes and herbicide behavior is fundamental for optimizing the efficient use of PRE-emergence herbicides in a more sustainable approach. However, it is still a poorly studied area within precision agriculture. Thus, the objective of this research was to evaluate the correlation of soil physicochemical attributes with the sorption and desorption processes of hexazinone and tebuthiuron to support application maps considering the field level variability. Soil samples from an agricultural area had their physicochemical attributes analyzed and were submitted to sorption and desorption studies of ¹⁴C-tebuthiuron and ¹⁴C-hexazinone using the batch equilibrium method. The values of sorption and desorption apparent coefficients (K_d), sorption and desorption percentage and bioavailability were correlated with soil attributes by Pearson's correlation. The K_d values of tebuthiuron and hexazinone sorption ranged from 1.2 to 2.9 mL g^-1 and 0.4-0.6 mL g^-1, respectively. For desorption of tebuthiuron and hexazinone, K_d values ranged from 3.4 to 4.4 mL g^-1 and 2.6-3.0 mL g^-1, respectively. A positive correlation among clay content, soil organic matter (OM), and tebuthiuron and hexazinone sorption K_d values were found. Both herbicides had variable retention according to geographic position in the area. The recommendation of application of PRE herbicides, such as tebuthiuron and hexazinone, observing the physicochemical attributes of the soil is an alternative to increase efficiency in weed control and decrease the risk of environmental contamination.

Differences in the sorption kinetics of various non-ionisable pesticides in a limited number of agricultural soils from the Mediterranean basin

Adsorption in soil of organic contaminants, such as pesticides, is a time-dependent process, which can be relevant for understanding and predicting the potential pollution risk of different water sources. The adsorption behavior of six different pesticides with a wide range of physicochemical properties (log K_OW 1.26-5.8) was evaluated in up to three different soils with low organic carbon (OC) content (≤1.2%). Pesticide sorbed amounts were fitted to several mathematical models to unravel the mechanisms involved in the adsorption process. The linear distribution constants revealed that pendimethalin and the pyrethroid insecticides were strongly retained in soil, whereas the other three compounds were moderately or weakly adsorbed. In the three soils, the pseudo second order model described more accurately the sorption kinetics of all the contaminants. The more hydrophobic pesticides (log K_OW ≥ 4.6) presented lower kinetic rates as compared with the other compounds under study. Both Elovich and intraparticle diffusion models reflected a strong contribution of a rapid initial adsorption on soil surface for thiacloprid, dimethenamid and fenarimol. For the hydrophobic pesticides this contribution was moderate according to the intraparticle diffusion model. Therefore, slower diffusion into the soil micropores was more relevant for the more hydrophobic compounds and for the bigger molecules, and less significant for the more polar pesticides because almost 90% of the total amount adsorbed was achieved in the rapid initial stage.

Adsorption-desorption of dimethenamid and fenarimol onto three agricultural soils as affected by treated wastewater and fresh sewage sludge-derived dissolved organic carbon

The use of treated wastewaters (TWW) in agriculture is widening in areas suffering drought, such as southern Europe, to preserve freshwater supply for human consumption. The composition of TWW, especially concerning their organic carbon (OC) content, has been demonstrated to influence the processes governing the behavior of non-ionic pesticides in soils. Three OC-poor agricultural soils (SV, RM1 and RM3) from the province of Granada (Spain) were chosen for the assessment of the adsorption and desorption of the herbicide dimethenamid (DIM) and the fungicide fenarimol (FEN). TWW and sewage sludge extracts at different dissolved OC (DOC) concentrations (30, 90 and 300 mg L^-1) were considered to evaluate their effect on pesticide adsorption-desorption. As expected by their properties, DIM and FEN were weakly and moderately adsorbed to the soils, respectively. Soil OC seemed to be the major factor controlling FEN adsorption, whereas the mineral fraction played a key role in DIM adsorption, especially in RM1 with high clay:OC ratio. Although TWW did not significantly modify the adsorption of pesticides, it enhanced DIM desorption from the three soils. Adsorption of FEN to SV and RM3 was directly related to the concentration of DOC, possibly due to co-sorption phenomena. Hysteretic desorption was found in all cases, indicating partially reversible adsorption. While FEN desorption was not altered by the solutions, the use of sludge extracts at the highest DOC concentration (300 mg L^-1) enhanced DIM desorption as occurred with TWW. Interactions with DOC in solution seemed to predominate for this less hydrophobic compound, thus increasing the risk of natural waters contamination if TWW will be used.

Meta-analysis of pesticide sorption in subsoil

Models used to assess leaching risks generally use organic carbon partition coefficient (k_OC ) values derived from batch experiments on topsoil samples to estimate pesticide sorption in subsoils of much smaller organic carbon contents. This can introduce significant errors in leaching risk calculations, because inorganic sorbents can play an important role for sorption in subsoil. The objectives of the present study were therefore to summarize the available literature data on pesticide sorption in subsoils and to test whether a simple alternative model could improve on the standard k_OC approach used in risk assessment models for pesticide leaching. This model describes the sorption constant as a power law function of the organic carbon content. A database with the results of batch sorption experiments was collated from published studies that emphasized measurements in subsoils. This database contains 1029 data entries from 36 published studies with data for 29 active substances (11 nonionic compounds, 10 weak acids, 6 weak bases, one cation, and one zwitterion). The results show that whereas the constant k_OC model proved to be an adequate model for 17 of the 63 individual datasets, the power law model gave acceptable fits (p < 0.05) for 60 of these cases. The exponent in the power law model varied over a wide range, from slightly negative to near unity. It also differed significantly (p = 0.015) for ionized and nonionized compounds, with median values of 0.25 and 0.55, respectively. It is concluded that the power law model could be used to parameterize subsoil sorption in regulatory leaching models, because it has widespread applicability and is simple enough for this purpose. Suitable ways of incorporating this approach in risk assessment procedures are discussed. Environ Toxicol Chem 2018;37:755-761. © 2017 SETAC.