Multi-process herbicide transport in structured soil columns: experiments and model analysis

Model predictions of pesticide transport in structured soils are complicated by multiple processes acting concurrently. In this study, the hydraulic, physical, and chemical nonequilibrium (HNE, PNE, and CNE, respectively) processes governing herbicide transport under variably saturated flow conditions were studied. Bromide (Br-), isoproturon (IPU, 3-(4-isoprpylphenyl)-1,1-dimethylurea) and terbuthylazine (TER, N2-tert-butyl-6-chloro-N4-ethyl-1,3,5-triazine-2,4-diamine) were applied to two soil columns. An aggregated Ap soil column and a macroporous, aggregated Ah soil column were irrigated at a rate of 1 cm h(-1) for 3 h. Two more irrigations at the same rate and duration followed in weekly intervals. Nonlinear (Freundlich) equilibrium and two-site kinetic sorption parameters were determined for IPU and TER using batch experiments. The observed water flow and Br- transport were inversely simulated using mobile-immobile (MIM), dual-permeability (DPM), and combined triple-porosity (DP-MIM) numerical models implemented in HYDRUS-1D, with improving correspondence between empirical data and model results. Using the estimated HNE and PNE parameters together with batch-test derived equilibrium sorption parameters, the preferential breakthrough of the weakly adsorbed IPU in the Ah soil could be reasonably well predicted with the DPM approach, whereas leaching of the strongly adsorbed TER was predicted less well. The transport of IPU and TER through the aggregated Ap soil could be described consistently only when HNE, PNE, and CNE were simultaneously accounted for using the DPM. Inverse parameter estimation suggested that two-site kinetic sorption in inter-aggregate flow paths was reduced as compared to within aggregates, and that large values for the first-order degradation rate were an artifact caused by irreversible sorption. Overall, our results should be helpful to enhance the understanding and modeling of multi-process pesticide transport through structured soils during variably saturated water flow.     

Experimental variability in characterization of cyfluthrin sorption to soil

To assess the risk of a pesticide to leach to groundwater or to run off to surface water after application, it is necessary to characterize the sorption of the pesticide to soil. For pyrethroids, their hydrophobicity, strong sorption to various materials, and low solubility make it difficult to accurately characterize sorption processes. The objective of this research was to evaluate the variability in cyfluthrin ((RS)-alpha -cyano-4-fluoro-3-phenoxybenzyl (1RS,3RS;1RS,3SR)-3-(2,2-dichlorovinyl)-2,2-dimethylcyclopropanecarboxylate) sorption to soil as affected by experiment conditions. To minimize cyfluthrin sorption on the walls of glass, silanized-glass, stainless steel, and PTFE centrifuge tubes, cyfluthrin solution was added to aqueous soil slurries or directly to soil, after which it was equilibrated with aqueous solution. Depending on the soil, variation in sorption coefficients, Koc, obtained using different experimental methodologies with one soil can be comparable to the variation in Koc values obtained for soils with different physical and chemical properties using one method. Koc values for cyfluthrin ranged from 56,000 to 300,000 L kg-1. Sorption methodology needs to be evaluated before sorption coefficients are used in predictive transport models.     

Rigorous uncertainty assessment in contaminant transport inverse modelling: a case study of fluoride diffusion through clay liners

Inverse methods used in assessing landfill liner design have not yet taken advantage of current developments in inverse procedures. Here, a method for inverting contaminant transport models is presented including a general error model and procedures for differentially weighted multiple response regression. General error models are employed in cases where the residuals are heteroscedastic and correlated, and lead to valid inference on model parameter and predictive uncertainty. The Shuffled Complex Evolution algorithm is used to optimise model parameters. Model parameter uncertainty is assessed by exploring the posterior probability distribution with the Metropolis algorithm, a Markov chain Monte Carlo sampling method. The inverse method is applied to simultaneously determine the sorption and diffusion parameters from laboratory diffusion cell experiments. In these experiments, fluoride migration through kaolin clays was measured by sampling the source and collector cells over time. To uniquely determine the transport model parameters, it was necessary to simultaneously fit the observed data from two independent diffusion cell experiments with different initial concentrations. The jointly fitted transport model parameters compared well with those fitted to independent batch experiments.     

Impact of varying soil structure on transport processes in different diagnostic horizons of three soil types

When soil structure varies in different soil types and the horizons of these soil types, it has a significant impact on water flow and contaminant transport in soils. This paper focuses on the effect of soil structure variations on the transport of pesticides in the soil above the water table. Transport of a pesticide (chlorotoluron) initially applied on soil columns taken from various horizons of three different soil types (Haplic Luvisol, Greyic Phaeozem and Haplic Cambisol) was studied using two scenarios of ponding infiltration. The highest infiltration rate and pesticide mobility were observed for the Bt₁ horizon of Haplic Luvisol that exhibited a well-developed prismatic structure. The lowest infiltration rate was measured for the Bw horizon of Haplic Cambisol, which had a poorly developed soil structure and a low fraction of large capillary pores and gravitational pores. Water infiltration rates were reduced during the experiments by a soil structure breakdown, swelling of clay and/or air entrapped in soil samples. The largest soil structure breakdown and infiltration decrease was observed for the Ap horizon of Haplic Luvisol due to the low aggregate stability of the initially well-aggregated soil. Single-porosity and dual-permeability (with matrix and macropore domains) flow models in HYDRUS-1D were used to estimate soil hydraulic parameters via numerical inversion using data from the first infiltration experiment. A fraction of the macropore domain in the dual-permeability model was estimated using the micro-morphological images. Final soil hydraulic parameters determined using the single-porosity and dual-permeability models were subsequently used to optimize solute transport parameters. To improve numerical inversion results, the two-site sorption model was also applied. Although structural changes observed during the experiment affected water flow and solute transport, the dual-permeability model together with the two-site sorption model proved to be able to approximate experimental data.     

Transport modes and pathways of the strongly sorbing pesticides glyphosate and pendimethalin through structured drained soils

Leaching of the strongly sorbing pesticides glyphosate and pendimethalin was evaluated in an 8-month field study focussing on preferential flow and particle-facilitated transport, both of which may enhance the leaching of such pesticides in structured soils. Glyphosate mainly sorbs to mineral sorption sites, while pendimethalin mainly sorbs to organic sorption sites. The two pesticides were applied in equal dosage to a structured, tile-drained soil, and the concentration of the pesticides was then measured in drainage water sampled flow-proportionally.The leaching pattern of glyphosate resembled that of pendimethalin, suggesting that the leaching potential of pesticides sorbed to either the inorganic or organic soil fractions is high in structured soils. Both glyphosate and pendimethalin leached from the root zone, with the average concentration in the drainage water being 3.5 and 2.7 μg L⁻¹, respectively. Particle-facilitated transport (particles >0.24 μm) accounted for only a small proportion of the observed leaching (13-16% for glyphosate and 16-31% for pendimethalin). Drain-connected macropores located above or in the vicinity of the drains facilitated very rapid transport of pesticide to the drains. That the concentration of glyphosate and pendimethalin in the drainage water remained high (>0.1 μg L⁻¹) for up to 7 d after a precipitation event indicates that macropores between the drains connected to underlying fractures were able to transport strongly sorbing pesticides in the dissolved phase. Lateral transport of dissolved pesticide via such discontinuities implies that strongly sorbing pesticides such as glyphosate and pendimethalin could potentially be present in high concentrations (>0.1 μg L⁻¹) in both water originating from the drainage system and the shallow groundwater located at the depth of the drainage system.     

Spatially distributed pesticide exposure assessment in the Central Valley, California, USA

Field runoff is an important transport mechanism by which pesticides move into the hydrologic environment of intensive agricultural regions such as California's Central Valley. This study presents a spatially explicit modeling approach to extend Pesticide Root Zone Model (PRZM), a field-scale pesticide transport model, into basin level. The approach was applied to simulate chlorpyrifos use in the Central Valley during 2003-2007. The average value of loading as percent of use (LAPU) is 0.031%. Results of this study provide strong evidence that surface runoff generation and pesticide application timing are the two influencing factors on the spatial and temporal variability of chlorpyrifos sources from agricultural fields. This is one of the first studies in coupling GIS and field-scale models and providing simulations for the dynamics of pesticides over an agriculturally dominated landscape. The demonstrated modeling approach may be useful for implementations of best management practice (BMP) and total maximum daily load (TMDL).     

Experimental Variability in Characterization of Cyfluthrin Sorption to Soil

To assess the risk of a pesticide to leach to groundwater or to run off to surface water after application, it is necessary to characterize the sorption of the pesticide to soil. For pyrethroids, their hydrophobicity, strong sorption to various materials, and low solubility make it difficult to accurately characterize sorption processes. The objective of this research was to evaluate the variability in cyfluthrin ((RS)-α -cyano-4-fluoro-3-phenoxybenzyl (1RS,3RS;1RS,3SR)-3-(2,2-dichlorovinyl)-2,2-dimethylcyclopropanecarboxylate) sorption to soil as affected by experiment conditions. To minimize cyfluthrin sorption on the walls of glass, silanized-glass, stainless steel, and PTFE centrifuge tubes, cyfluthrin solution was added to aqueous soil slurries or directly to soil, after which it was equilibrated with aqueous solution. Depending on the soil, variation in sorption coefficients, K _oc, obtained using different experimental methodologies with one soil can be comparable to the variation in K _oc values obtained for soils with different physical and chemical properties using one method. K _oc values for cyfluthrin ranged from 56,000 to 300,000 L kg^?1. Sorption methodology needs to be evaluated before sorption coefficients are used in predictive transport models.     

Monte Carlo analysis of field water flow comparing uni- and bimodal effective hydraulic parameters for structured soil

Soil structure critically affects the hydrological behaviour of soils. In this paper, we examined the impact of areal heterogeneity of hydraulic properties of a structured soil on soil ensemble behaviour for various soil water flow processes with different top boundary conditions (redistribution and drainage plus evaporation and infiltration). Using a numerical solution of the Richards' equation in a stochastic framework, the ensemble characteristics and flow dynamics were studied for drying and wetting processes observed during a time interval of ten days when a series of relatively intense rainfall events occurred. The effects of using unimodal and bimodal interpretative models of hydraulic properties on the ensemble hydrological behaviour of the soil were illustrated by comparing predictions to mean water contents measured over time in several sites at field scale. Although the differences between unimodal and bimodal fitting are not significant in terms of goodness of fit, the differences in process predictions are considerable with the bimodal soil simulating water content measurements much better than unimodal soil. We also investigated the relative contribution of the soil variability of each parameter on the variance of the water contents obtained as the main output of the stochastic simulations. The variability of the structural parameter, weighting the two pore space fractions in the bimodal interpretative model, has the largest contribution to water content variance. The contribution of each parameter depends only partly on the coefficient of variation, much more on the sensitivity of the model to the parameters and on the flow process being observed. We observed that the contribution of the retention parameters to uncertainty increases during drainage processes; the opposite occurs with the hydraulic conductivity parameters.     

Formulation of a soil-pesticide transport model based on a compartmental approach

A semianalytical soil-pesticide transport model is formulated based on a compartmental approach to determine spatial and temporal variations of pesticide residues across a soil profile. The compartmental model is implemented by drawing an analogy between a series of continuous-flow stirred tank reactors and a soil horizon that consists of multiple perfectly mixed compartments. The analogy is strengthened by exploiting a relation between the compartment series and the conventional convective-dispersive equation (CDE) for vertical transport in the soil. Consequently, the number of compartments in the model formulation is not free, but dictated as a function of transport parameters. The model formulation allows consideration of arbitrary boundary value specifications and also, for some cases, spatially varying initial concentration profiles. Sorption kinetics is represented via a two-site model that involves a linear sorption isotherm and a first-order irreversible sorption or a radial diffusive penetrating model. For these three cases, analysis of the compartmental model allows the resultant concentration profiles to be expressed in terms of the Poisson distribution. When a nonlinear kinetic sorption model is used to simulate the sorption processes, an analytical solution is not found and a numerical approach is required.     

Dynamic modeling of organophosphate pesticide load in surface water in the northern San Joaquin Valley watershed of California

The hydrology, sediment, and pesticide transport components of the Soil and Water Assessment Tool (SWAT) were evaluated on the northern San Joaquin Valley watershed of California. The Nash-Sutcliffe coefficients for monthly stream flow and sediment load ranged from 0.49 to 0.99 over the watershed during the study period of 1992-2005. The calibrated SWAT model was applied to simulate fate and transport processes of two organophosphate pesticides of diazinon and chlorpyrifos at watershed scale. The model generated satisfactory predictions of dissolved pesticide loads relative to the monitoring data. The model also showed great success in capturing spatial patterns of dissolved diazinon and chlorpyrifos loads according to the soil properties and landscape morphology over the large agricultural watershed. This study indicated that curve number was the major factor influencing the hydrology while pesticide fate and transport were mainly affected by surface runoff and pesticide application and in the study area.     

Modeling field-scale multiple tracer injection at a low-level waste disposal site in fractured rocks: effect of multiscale heterogeneity and source term uncertainty on conceptual understanding of mass transfer processes

Multiple factors may affect the scale-up of laboratory multi-tracer injection into structured porous media to the field. Under transient flow conditions and with multiscale heterogeneities in the field, previous attempts to scale-up laboratory experiments have not answered definitely the questions about the governing mechanisms and the spatial extent of the influence of small-scale mass transfer processes such as matrix diffusion. The objective of this research is to investigate the effects of multiscale heterogeneity, mechanistic and site model conceptualization, and source term density effect on elucidating and interpreting tracer movement in the field. Tracer release and monitoring information previously obtained in a field campaign of multiple, conservative tracer injection under natural hydraulic gradients at a low-level waste disposal site in eastern Tennessee, United States, is used for the research. A suite of two-pore-domain, or fracture-matrix, groundwater flow and transport models are calibrated and used to conduct model parameter and prediction uncertainty analyses. These efforts are facilitated by a novel nested Latin-hypercube sampling technique. Our results verify, at field scale, a multiple-pore-domain, multiscale mechanistic conceptual model that was used previously to interpret only laboratory observations. The results also suggest that, integrated over the entire field site, mass flux rates attributable to small-scale mass transfer are comparable to that of field-scale solute transport. The uncertainty analyses show that fracture spacing is the most important model parameter and model prediction uncertainty is relatively higher at the interface between the preferred flow path and its parent bedrock. The comparisons of site conceptual models indicate that the effect of matrix diffusion may be confined to the immediate neighborhood of the preferential flow path. Finally, because the relatively large amount of tracer needed for field studies, it is likely that source term density effect may exaggerate or obscure the effect of matrix diffusion on the movement of tracers from the preferred flow path into the bedrock.     

Sorption and leaching potential of acidic herbicides in Brazilian soils

Leaching of acidic herbicides (2,4-D, flumetsulam, and sulfentrazone) in soils was estimated by comparing the original and modified AF (Attenuation Factor) models for multi-layered soils (AFi). The original AFi model was modified to include the concept of pH-dependence for Kd (sorption coefficient) based on pesticide dissociation and changes in the accessibility of soil organic functional groups able to interact with the pesticide. The original and modified models, considering soil and herbicide properties, were applied to assess the leaching potential of selected herbicides in three Brazilian soils. The pH-dependent Kd values estimated for all three herbicides were observed to be always higher than pH-independent Kd values calculated using average Koc data, and therefore the original AFi model overestimated the overall leaching potential for the soils studied.     

Limitations on the role of incorporated organic matter in reducing pesticide leaching.

The use of organic amendments has been suggested as a method of controlling pesticide leaching through soils. The enarenados soils of the intensive horticulture of the Almeria province of southern Spain contain buried organic matter horizons above a soil layer amended with clay. This region is ideal for understanding the potential for and limitations of organic amendments in preventing pesticide pollution. This study measured the sorption and degradation potential of carbofuran in this soil system and the hydrological behaviour of the soil horizons. The sorption of carbofuran was controlled by the organic carbon content, the degradation was strongly pH-dependent and the acidic organic layer protected the sorbed carbofuran against degradation. Hydrologically, the soil system is dominated by ponding above an amended clay layer and by the presence of macropores that can transport water through this clay. A simple model is proposed on this basis and shows that although high levels of dissolved organic carbon can be released by buried organic horizons, the major control on re-release of sorbed pesticide is the potential for sorption hysteresis in this organic layer. A comparison of sorption and degradation data for carbamate insecticides used in the region with groundwater observations for these compounds shows that no amount of incorporated organic would protect against pollution from highly water-soluble compounds.     

Capability and limitations of first-order and diffusion approaches to describe long-term sorption of chlortoluron in soil

This paper compares the capability of a first-order and a spherical diffusion model to describe and predict long-term sorption and desorption processes of chlortoluron in two soils. Chlortoluron sorption was investigated at different time scales utilizing one rate experiment (120 days) and two sorption/desorption experiments. Experimental periods for sorption and desorption were set to 1 day (five desorption steps) and 30 days (three desorption steps), respectively. Upon fitting, the two models satisfactorily described the whole set of data. The spherical diffusion model performed better than the first-order model. We then tested the predictive capability of the models by predicting 30-day sorption/desorption data using kinetic parameters fitted on 1-day sorption/desorption data only. While the spherical diffusion model was able to predict the 30-day data set, the first-order model failed completely. Fitting both models to subsets of the data corresponding to different experimental time scales revealed that the rate parameter as well as the Freundlich coefficient of the first-order model are strongly time-dependent--a property that is not shared by parameters of the spherical diffusion model. The apparent stability of the spherical diffusion model with regard to time dependency of its parameters indicates that sorptive uptake may be diffusion-controlled. This also explains the models greater predictive power across different time scales compared to the first-order model. Finally, we investigate the suitability of solute class specific log-linear relationships between the first-order rate parameter and the Freundlich coefficient presented by earlier researchers in the light of the time dependency observed for the parameters of the first-order model.     

A lysimeter experiment to investigate temporal changes in the availability of pesticide residues for leaching

Leaching of three pesticides (isoproturon, chlorotoluron and triasulfuron) and a tracer (bromide) were determined in four contrasting soils ranging in texture from sandy loam to clay. The compounds were applied to undisturbed columns of soil and four columns for each soil were randomly selected and leached with 24-mm equivalent of water at prescribed time intervals (3, 9, 24, 37 and 57 d after application). A rapid decline in leached loads of isoproturon and chlorotoluron as time from application to irrigation increased was observed in all soils. In contrast, triasulfuron and bromide loads only decreased rapidly in the clay soil. Bromide losses decreased with decreasing clay contents of the soil and therefore with a decrease in structural development. Magnitudes of pesticide losses varied from soil to soil, depending on structural development and the organic carbon content. Pesticide degradation experiments on disturbed and undisturbed soil samples showed that the rapid decline of leached loads with time was faster than could be explained by degradation alone. Five physico-chemical processes are put forward to explain the different patterns of pesticide leached loads observed in the soils: (1) relative extent of preferential flow, (2) sorption capacity of the compounds to the different soils, (3) extent of degradation of the compounds in the soil, (4) variation in sorption kinetics between compounds associated with pesticide diffusion into soil aggregates, and (5) protection of the compounds by a combination of intra-aggregate diffusion and the presence of preferential flow pathways.     

Sorption of metolachlor and atrazine in fly ash amended soils: comparison of optimized isotherm models

Adsorption of metolachlor and atrazine was studied in the fly ash (Inderprastha and Badarpur)- amended Inceptisol and Alfisol soils using batch method. Results indicated that sorption of both the herbicides in soil+fly ash mixtures was highly nonlinear and sorption decreased with a higher herbicide concentration in the solution. Also, nonlinearity increased with an increase in the level of fly ash amendment from 0-5%. Three two-parameter monolayer isotherms viz. Langmuir, Temkin, Jovanovic and one imperical Freundlich models were used to fit the experimental data. Data analysis and comparison revealed that the Temkin and the Freundlich isotherms were best-suited to explain the sorption results and the observed and the calculated adsorption coefficient values showed less variability. The study suggested that sorption mechanism of metolachlor and atrazine involved the physical association at the sorbate surface and the nonlinearity in the sorption at higher pesticide or fly ash concentration was due to a decrease in the heat of adsorption and higher binding energy.