Impact of data quality and model complexity on prediction of pesticide leaching

Remediation of pesticides in commercial farm soils by solarization and ozonation techniques

Soil contamination by pesticides is a growing environmental problem. Even though nowadays numerous soil remediation technologies are available, most of them have not been tested at field scale. This study attempts to demonstrate the efficiency of solarization-ozonation techniques for the removal of twelve pesticides at full scale. Initial solarization and ozonation trials were conducted in plots located in a greenhouse using freshly and aged contaminated soils under controlled pilot conditions. The combination of solarization and ozonation treatment was efficient for all the studied pesticides both in freshly and in aged contaminated soils, being the lower degradation values found for the second type. This low removal suggests that the increase of pesticides' adsorption on soil resulting from ageing decreases their availability. Once the essays were carried out at pilot scale, the solarization-ozonation applicability was evaluated in a commercial farm soil. This trial was carried out in a greenhouse whose soil had previously been contaminated with some of the pesticides studied. A significant degradation (53.8%) was observed after 40 days of treatment. Pesticides' main metabolites were identified during the different remediation experiments. In addition, the cost of the combined solarization and ozonation technology was evaluated. Finally, our results suggest that this combination of techniques could be considered a promising technology to degrade pesticides in soil.

A protocol to build soil descriptions for APSIM simulations

•

Introducing the models and user interface for characterising a soil in APSIM simulations.

•

Listing and describing the parameters needed for building soil descriptions in APSIM.

•

Providing recommendations for good practice when setting up soil parameters in APSIM.

Soil processes have a major impact on agroecosystems, controlling water and nutrient cycling, regulating plant growth and losses to the wider environment. Process-based agroecosystem simulation models generally encompass detailed descriptions of the soil, including a wide number of parameters that can be daunting to users with a limited soil science background. In this work we review and present an abridged description of the models used to simulate soil processes in the APSIM (Agricultural Production Systems sIMulator) framework. Such a resource is needed because this information is currently spread over multiple publications and some elements have become outdated. We list and briefly describe the parameters, and establish a protocol with guidelines, for building a soil description for APSIM. This protocol will promote consistency, enhancing the quality of the science done employing APSIM, and provide an easier pathway for new users. This compilation should also be of relevance to users of other models that require detailed soil information.•

This paper presents a brief description of the models for simulating soil processes in the APSIM model.

•

The method stablishes guidelines to define the parameters for building a soil description for APSIM.

Real-time forecasting of pesticide concentrations in soil

Forecasting pesticide residues in soils in real time is essential for agronomic purposes, to manage phytotoxic effects, and in catchments to manage surface and ground water quality. This has not been possible in the past due to both modelling and measurement constraints. Here, the analytical transient probability distribution (pdf) of pesticide concentrations is derived. The pdf results from the random ways in which rain events occur after pesticide application. First-order degradation kinetics and linear equilibrium sorption are assumed. The analytical pdfs allow understanding of the relative contributions that climate (mean storm depth and mean rainfall event frequency) and chemical (sorption and degradation) properties have on the variability of soil concentrations into the future. We demonstrated the two uncertain reaction parameters can be constrained using Bayesian methods. An approach to a Bayesian informed forecast is then presented. With the use of new rapid tests capable of providing quantitative measurements of soil concentrations in the field, real-time forecasting of future pesticide concentrations now looks possible for the first time. Such an approach offers new means to manage crops, soils and water quality, and may be extended to other classes of pesticides for ecological risk assessment purposes.

Environmental and Agricultural Relevance of Humic Fractions Extracted by Alkali from Soils and Natural Waters

To study the structure and function of soil organic matter, soil scientists have performed alkali extractions for soil humic acid (HA) and fulvic acid (FA) fractions for more than 200 years. Over the last few decades aquatic scientists have used similar fractions of dissolved organic matter, extracted by resin adsorption followed by alkali desorption. Critics have claimed that alkali-extractable fractions are laboratory artifacts, hence unsuitable for studying natural organic matter structure and function in field conditions. In response, this review first addresses specific conceptual concerns about humic fractions. Then we discuss several case studies in which HA and FA were extracted from soils, waters, and organic materials to address meaningful problems across diverse research settings. Specifically, one case study demonstrated the importance of humic substances for understanding transport and bioavailability of persistent organic pollutants. An understanding of metal binding sites in FA and HA proved essential to accurately model metal ion behavior in soil and water. In landscape-based studies, pesticides were preferentially bound to HA, reducing their mobility. Compost maturity and acceptability of other organic waste for land application were well evaluated by properties of HA extracted from these materials. A young humic fraction helped understand N cycling in paddy rice ( L.) soils, leading to improved rice management. The HA and FA fractions accurately represent natural organic matter across multiple environments, source materials, and research objectives. Studying them can help resolve important scientific and practical issues.

A Global Meta-Analysis to Predict Atrazine Sorption from Soil Properties

Atrazine (2-chloro-4-ethylamino-6-isopropylamino-1,3,5-triazine) is one of the most widely used herbicides worldwide, and groundwater contamination is of concern, especially in heavily used regions and in edaphic conditions prone to leaching. Soil sorption plays an essential role in atrazine environmental fate, yet consistent atrazine risk prediction remains limited. A quantitative meta-analysis was conducted to characterize the effect of soil properties on atrazine sorption, using 378 previous observations in 48 publications from 1985 to 2015 globally, which included data on soil properties and sorption parameters. A supplemental regional study was conducted to test the derived meta-analysis models. The meta-analysis indicated that percentage organic C (OC) was the most important parameter for estimating atrazine sorption, followed by percentage silt, soil pH, and percentage clay. Meta-analysis and supplemental study models were developed for Freundlich sorption coefficients () and sorption distribution coefficients () as a function of OC. The global meta-analysis models generated positive linear trends for OC with and ( = 0.197 and 0.205, respectively). Organic C was highly correlated with and in supplemental experimental study models ( = 0.93 and 0.92, respectively), indicating accurate prediction of sorption within the evaluated region. Continental models were investigated, which improved the goodness of fit. Models developed via meta-analysis may be used to predict atrazine sorption over wide ranges of data, whereas more accurate and refined prediction can be achieved by specific regional models through experimental studies. However, such models could be improved if standardized agroclimatic conditions, soil classification, and other key variables were more widely reported.

Integrated assessment of climate change impact on surface runoff contamination by pesticides

Pesticide transport by surface runoff depends on climate, agricultural practices, topography, soil characteristics, crop type, and pest phenology. To accurately assess the impact of climate change, these factors must be accounted for in a single framework by integrating their interaction and uncertainty. This article presents the development and application of a framework to assess the impact of climate change on pesticide transport by surface runoff in southern Québec (Canada) for the 1981-2040 period. The crop enemies investigated were: weeds for corn (Zea mays); and for apple orchard (Malus pumila), 3 insect pests (codling moth [Cydia pomonella], plum curculio [Conotrachelus nenuphar], and apple maggot [Rhagoletis pomonella]), 2 diseases (apple scab [Venturia inaequalis], and fire blight [Erwinia amylovora]). A total of 23 climate simulations, 19 sites, and 11 active ingredients were considered. The relationship between climate and phenology was accounted for by bioclimatic models of the Computer Centre for Agricultural Pest Forecasting (CIPRA) software. Exported loads of pesticides were evaluated at the edge-of-field scale using the Pesticide Root Zone Model (PRZM), simulating both hydrology and chemical transport. A stochastic model was developed to account for PRZM parameter uncertainty. Results of this study indicate that for the 2011-2040 period, application dates would be advanced from 3 to 7 days on average with respect to the 1981-2010 period. However, the impact of climate change on maximum daily rainfall during the application window is not statistically significant, mainly due to the high variability of extreme rainfall events. Hence, for the studied sites and crop enemies considered, climate change impact on pesticide transported in surface runoff is not statistically significant throughout the 2011-2040 period. Integr Environ Assess Managem 2016;12:559-571. © Her Majesty the Queen in Right of Canada 2015; Published 2015 SETAC.

Do lab-derived distribution coefficient values of pesticides match distribution coefficient values determined from column and field-scale experiments? A critical analysis of relevant literature

In this study, we analyzed sorption parameters for pesticides that were derived from batch and column or batch and field experiments. The batch experiments analyzed in this study were run with the same pesticide and soil as in the column and field experiments. We analyzed the relationship between the pore water velocity of the column and field experiments, solute residence times, and sorption parameters, such as the organic carbon normalized distribution coefficient ( ) and the mass exchange coefficient in kinetic models, as well as the predictability of sorption parameters from basic soil properties. The batch/column analysis included 38 studies with a total of 139 observations. The batch/field analysis included five studies, resulting in a dataset of 24 observations. For the batch/column data, power law relationships between pore water velocity, residence time, and sorption constants were derived. The unexplained variability in these equations was reduced, taking into account the saturation status and the packing status (disturbed-undisturbed) of the soil sample. A new regression equation was derived that allows estimating the values derived from column experiments using organic matter and bulk density with an value of 0.56. Regression analysis of the batch/column data showed that the relationship between batch- and column-derived values depends on the saturation status and packing of the soil column. Analysis of the batch/field data showed that as the batch-derived value becomes larger, field-derived values tend to be lower than the corresponding batch-derived values, and vice versa. The present dataset also showed that the variability in the ratio of batch- to column-derived value increases with increasing pore water velocity, with a maximum value approaching 3.5.

Toxicity of neurons treated with herbicides and neuroprotection by mitochondria-targeted antioxidant SS31

The purpose of this study was to determine the neurotoxicity of two commonly used herbicides: picloram and triclopyr and the neuroprotective effects of the mitochondria-targeted antioxidant, SS31. Using mouse neuroblastoma (N2a) cells and primary neurons from C57BL/6 mice, we investigated the toxicity of these herbicides, and protective effects of SS1 peptide against picloram and triclopyr toxicity. We measured total RNA content, cell viability and mRNA expression of peroxiredoxins, neuroprotective genes, mitochondrial-encoded electron transport chain (ETC) genes in N2a cells treated with herbicides and SS31. Using primary neurons from C57BL/6 mice, neuronal survival was studied in neurons treated with herbicides, in neurons pretreated with SS31 plus treated with herbicides, neurons treated with SS31 alone, and untreated neurons. Significantly decreased total RNA content, and cell viability in N2a cells treated with picloram and triclopyr were found compared to untreated N2a cells. Decreased mRNA expression of neuroprotective genes, and ETC genes in cells treated with herbicides was found compared to untreated cells. Decreased mRNA expression of peroxiredoxins 1-6 in N2a cells treated with picloram was found, suggesting that picloram affects the antioxidant enzymes in N2a cells. Immunofluorescence analysis of primary neurons revealed that decreased neuronal branching and degenerating neurons in neurons treated with picloram and triclopyr. However, neurons pretreated with SS31 prevented degenerative process caused by herbicides. Based on these results, we propose that herbicides--picloram and triclopyr appear to damage neurons, and the SS31 peptide appears to protect neurons from herbicide toxicity.

Agronomic and environmental implications of enhanced s-triazine degradation

Novel catabolic pathways enabling rapid detoxification of s-triazine herbicides have been elucidated and detected at a growing number of locations. The genes responsible for s-triazine mineralization, i.e. atzABCDEF and trzNDF, occur in at least four bacterial phyla and are implicated in the development of enhanced degradation in agricultural soils from all continents except Antarctica. Enhanced degradation occurs in at least nine crops and six crop rotation systems that rely on s-triazine herbicides for weed control, and, with the exception of acidic soil conditions and s-triazine application frequency, adaptation of the microbial population is independent of soil physiochemical properties and cultural management practices. From an agronomic perspective, residual weed control could be reduced tenfold in s-triazine-adapted relative to non-adapted soils. From an environmental standpoint, the off-site loss of total s-triazine residues could be overestimated 13-fold in adapted soils if altered persistence estimates and metabolic pathways are not reflected in fate and transport models. Empirical models requiring soil pH and s-triazine use history as input parameters predict atrazine persistence more accurately than historical estimates, thereby allowing practitioners to adjust weed control strategies and model input values when warranted.

The important characteristics of soil organic matter affecting 2,4-dichlorophenoxyacetic acid sorption along a catenary sequence

Variations in the characteristics of soil organic matter (SOM) at the field-scale are largely unknown, particularly in relation to observed variations in herbicide sorption. For the herbicide 2,4-D [2,4-dichlorophenoxyacetic acid], we found that its organic carbon-normalized sorption coefficient, Koc, varied by four-fold, from 76 to 315 L kg(-1), in the Ap-horizon along a slope transect in an undulating agricultural field in Manitoba, Canada. In order to explain the relatively large in-field variation in 2,4-D Koc values, techniques ranging from conventional chemical fractionation methods to solid state Cross Polarization and Magic-Angle Spinning (13)C-Nuclear Magnetic Resonance applied on whole soils, were used to derive SOM chemical, physical and structural parameters for correlation analyses with the measured 2,4-D Koc values on whole soils. Out of the 15 parameters considered, the 2,4-D Koc was significantly positively correlated with 1) the carbon (C) content of sodium hydroxide-extracted humic acids (r = 0.83, P < 0.01), a chemical parameter indicative of free form C in soil; 2) the molar absorptivity of humic acids at wavelength 280 nm (r = 0.81, P < 0.01), a physical parameter indicative of greater SOM aromaticity; and 3) the relatively intensity of aryl C (r = 0.92, P < 0.01) and O-aryl C (r = 0.93, P < 0.01) in whole soil, both structural parameters indicative of aromatic C. Consequently, the results suggest that in-field variations in 2,4-D Koc values are induced by variations in SOM aromaticity. Koc values are among the most sensitive parameters in herbicide fate models used in regulatory and environmental assessments. Currently, these herbicide fate models do not consider associations between SOM characteristics and Koc and hence revising model equations to include these associations may improve estimates of herbicide persistence, bioavailability and transport at the field-scale.

Land use and riparian effects on prairie wetland sediment properties and herbicide sorption coefficients

Sorption of commonly used herbicides by wetland sediment can provide important information for herbicide fate modeling. The influence of sediment properties on herbicide sorption as a result of different land uses in the wetland catchment is unclear. We examined the effects of land use on the physiochemical properties of wetland sediments and the associations between these sediment properties and herbicide sorption characteristics. Bottom sediments were sampled in 0- to 5- and 5- to 10-cm sections from 17 wetlands under five different land use classes: (i) ephemeral wetlands with no riparian vegetation in a cultivated catchment (ECNR), (ii) ephemeral wetlands with riparian vegetation in a cultivated catchment (ECR), (iii) ephemeral wetlands in a grassland catchment established 4 yr ago (E4G), (iv) ephemeral wetlands in a brome grass catchment established 20 yr ago (E20G), and (v) semi-permanent (SP) wetlands in a multiple-land-use catchment. Sediments were analyzed for total organic carbon (TOC), total inorganic carbon (TIC), pH, electrical conductivity, exchangeable cations (EXCAT), total cation exchangeable capacity (CEC), and percent clay (%clay). Sediment herbicide sorption partition coefficient (Kd) was measured for trifluralin, atrazine, 2,4-D, and glyphosate. The sorption of the herbicides in the sediment increased in the order of 2,4-D < atrazine < glyphosate < trifluralin. The sorption of 2,4-D, atrazine, and trifluralin was positively correlated to TOC, EXCAT, and CEC but negatively correlated to %clay. Glyphosate sorption was negatively correlated to pH, TIC, EXCAT, and %clay. Overall, wetland sediments that were recently cultivated (ECNR and E4G) had lower TOC, TIC, EC, EXCAT, CEC, and Kd values (2,4-D, trifluralin, and atrazine) than sediments that had not been recently cultivated (ECR, E20G, and SP). The ECR wetland sediments had the largest Kd for all four herbicides, suggesting that land use and riparian vegetation have a significant impact on herbicide sorption.

Characterization of nicosulfuron availability in aged soils

Sorption-desorption interactions of pesticides with soil determine their availability for transport, plant uptake, and microbial degradation. These interactions are affected by the physical-chemical properties of the pesticide and soil, and for some pesticides, their residence time in the soil. This research evaluated changes in sorption/availability of nicosulfuron (2-[[[[(4,6-dimethoxy-2-pyrimidinyl]amino]carbonyl]amino]sulfonyl]-N,N-dimethyl-3-pyridinecarboxamide) herbicide with aging in different soils, using a radiolabeled ((14)C) tracer. Aging significantly increased sorption. For instance, after the 41-day incubation, calculated K d,app increased by a factor of 2 to 3 in Mollisols from the Midwestern United States and by a factor of 5 to 9 in Oxisols from Brazil and Hawaii, as compared to freshly treated soils. In view of this outcome, potential transport of nicosulfuron would be overpredicted if freshly treated soil K d values were used to predict transport. The fact that the nicosulfuron solution concentration decreased faster than the soil concentration with time suggested that the increase in sorption was because the rate of degradation in solution and on labile sites was faster than the rate of desorption of the neutral species from the soil particles. It may have also been due to nicosulfuron anion diffusion to less accessible sites with time, leaving the more strongly bound neutral molecules for the sorption characterization. Regardless of the mechanism, these results are further evidence that increases in sorption during pesticide aging should be taken into account during the characterization of the sorption process for mathematical models of pesticide degradation and transport.

Availability of triazine herbicides in aged soils amended with olive oil mill waste

Amendments are frequently added to agricultural soils to increase organic matter content. In this study, we examined the influence of alperujo, an olive oil mill waste, on the availability of two triazine herbicides, terbuthylazine and atrazine, in two different sandy soils, one from Sevilla, Spain, and the other from Minnesota. The effect of aging on herbicide sorption and bioavailability was also studied. Soils were amended with alperujo at a rate of 3-5% (w:w) in laboratory studies. Apparent sorption coefficients for the triazine herbicides were calculated as the ratio of the concentrations of each herbicide sequentially extracted with water, followed by aqueous methanol, at each sampling time. These data showed greater sorption of terbuthylazine and atrazine in amended soils as compared to nonamended soils, and an increase in the amount of herbicide sorbed with increasing aging time in nonamended soils. The triazine-mineralizing bacterium Pseudomonas sp. strain ADP was used to characterize triazine bioavailability. Less mineralization of the herbicides by Pseudomonas sp. strain ADP was observed in soils amended with alperujo, as compared to the unamended soils, and, despite the increase in sorption with aging in unamended soils, herbicide mineralization also increased in this case. This has been attributed to Pseudomonas sp. strain ADP first using alperujo as a more readily available source of N as compared to the parent triazines. In summary, addition of alperujo to the soils studied was shown to increase triazine herbicides sorption and hence to reduce its availability and potential to leach.

Effective elution of RDX and TNT from particles of Comp B in surface soil

During live fire training exercises, large amounts of explosives are consumed. Low order detonations of high explosive payloads result in the patchy dispersal of particles of high explosive formulations over large areas of firing range soils. Dissolution of explosives from explosive formulation particles into soil pore water is a controlling factor for transport, fate, and effects of explosive compounds. We developed an empirical method to evaluate soils based on functionally defined effective dissolution rates. An automated Accelerated Solvent Extractor was used to determine the effective elution rates under controlled conditions of RDX and TNT from soil columns containing particles of Comp B. Contrived soils containing selected soil geosorbants and reactive surfaces were used to quantitatively determine the importance of these materials. Natural soils from training ranges of various soil types were also evaluated. The effects of geosorbants on effective elution rates were compound- and sorbent-specific. TNT elution was less than that of RDX and was greatly slowed by humic acid. Iron and iron-bearing clays reduced the effective elution rates of both RDX and TNT. This empirical method is a useful tool for directly generating data on the potential for explosives to leach from firing range soils, to identify general bulk soil characteristics that can be used to predict the potential, and to identify means to engineer soil treatments to mitigate potential transport.

Leaching of seven pesticides currently used in cotton crop in Mato Grosso State-Brazil

This study aimed to evaluate the leaching of pesticides and the applicability of the Attenuation Factor (AF) Model to predict their leaching. The leaching of carbofuran, carbendazim, diuron, metolachlor, alpha and beta endosulfan and chlorpyrifos was studied in an Oxisol using a field experiment lysimeter located in Dom Aquino-Mato Grosso. The samples of percolated water were collected by rain event and analyzed. Chemical and physical soil attributes were determined before pesticide application to the plots. The results showed that carbofuran was the pesticide that presented a higher leaching rate in the studied soil, so was the one representing the highest contamination potential. From the total carbofuran applied in the soil surface, around 6% leached below 50 cm. The other pesticides showed lower mobility in the studied soil. The calculated values to AF were 7.06E-12 (carbendazim), 5.08E-03 (carbofuran), 3.12E-17 (diuron), 6.66E-345 (alpha-endosulfan), 1.47E-162 (beta-endosulfan), 1.50E-06 (metolachlor), 3.51E-155 (chlorpyrifos). AF Model was useful to classify the pesticides' potential for contamination; however, that model underestimated pesticide leaching.

Sorption and predicted mobility of herbicides in Baltic soils

This study was undertaken to determine sorption coefficients of eight herbicides (alachlor, amitrole, atrazine, simazine, dicamba, imazamox, imazethapyr, and pendimethalin) to seven agricultural soils from sites throughout Lithuania. The measured sorption coefficients were used to predict the susceptibility of these herbicides to leach to groundwater. Soil-water partitioning coefficients were measured in batch equilibrium studies using radiolabeled herbicides. In most soils, sorption followed the general trend pendimethalin > alachlor > atrazine approximately amitrole approximately simazine > imazethapyr > imazamox > dicamba, consistent with the trends in hydrophobicity (log K(ow)) except in the case of amitrole. For several herbicides, sorption coefficients and calculated retardation factors were lowest (predicted to be most susceptible to leaching) in a soil of intermediate organic carbon content and sand content. Calculated herbicide retardation factors were high for soils with high organic carbon contents. Estimated leaching times under saturated conditions, assuming no herbicide degradation and no preferential water flow, were more strongly affected by soil textural effects on predicted water flow than by herbicide sorption effects. All herbicides were predicted to be slowest to leach in soils with high clay and low sand contents, and fastest to leach in soils with high sand content and low organic matter content. Herbicide management is important to the continued increase in agricultural production and profitability in the Baltic region, and these results will be useful in identifying critical areas requiring improved management practices to reduce water contamination by pesticides.