Modeling depth-variant and domain-specific sorption and biodegradation in dual-permeability media

The impact of managed aquifer recharge on the fate and transport of pesticides in agricultural soils

Groundwater aquifers worldwide experience unsustainable depletion, compounded by population growth, economic development, and climate forcing. Managed aquifer recharge provides one tool to alleviate flood risk and replenish groundwater. However, concerns grow that intentional flooding of farmland for groundwater recharge, a practice known as Ag-MAR, may increase the leaching of pesticides and other chemicals into groundwater. This study employs a physically based unsaturated flow model to determine the fate and transport of residues of four pesticide in three vadose zone profiles characterized by differing fractions of sand (41 %, 61 %, and 84 %) in California's Central Valley. Here, we show that the complex heterogeneity of alternating coarse and fine-grain hydrogeologic units controls the transit times of pesticides and their adsorption and degradation rates. Unsaturated zones that contain a higher fraction of sand are more prone to support preferential flow, higher recharge rates (+8 %), and faster (42 %) water flow and pesticide transport, more flooding-induced pesticide leaching (about 22 %), as well as more salt leaching correlating with increased risks of groundwater contamination. Interestingly, considering preferential flow predicted higher degradation and retention rates despite shorter travel times, attributed to the trapping of pesticides in immobile zones where they degrade more effectively. The findings underscore the importance of considering soil texture and structure in Ag-MAR practices to minimize environmental risks while enhancing groundwater recharge. The study also highlights that selecting less mobile pesticides can reduce leaching risks in sandy areas.

Transport of bromide and pesticides through an undisturbed soil column: a modeling study with global optimization analysis

The fate of pesticides in tropical soils is still not understood as well as it is for soils in temperate regions. In this study, water flow and transport of bromide tracer and five pesticides (atrazine, imazaquin, sulfometuron methyl, S-metolachlor, and imidacloprid) through an undisturbed soil column of tropical Oxisol were analyzed using a one-dimensional numerical model. The numerical model is based on Richards' equation for solving water flow, and the advection-dispersion equation for solving solute transport. Data from a laboratory column leaching experiment were used in the uncertainty analysis using a global optimization methodology to evaluate the model's sensitivity to transport parameters. All pesticides were found to be relatively mobile (sorption distribution coefficients lower than 2 cm(3) g(-1)). Experimental data indicated significant non-conservative behavior of bromide tracer. All pesticides, with the exception of imidacloprid, were found less persistent (degradation half-lives smaller than 45 days). Three of the five pesticides (atrazine, sulfometuron methyl, and S-metolachlor) were better described by the linear kinetic sorption model, while the breakthrough curves of imazaquin and imidacloprid were more appropriately approximated using nonlinear instantaneous sorption. Sensitivity analysis suggested that the model is most sensitive to sorption distribution coefficient. The prediction limits contained most of the measured points of the experimental breakthrough curves, indicating adequate model concept and model structure for the description of transport processes in the soil column under study. Uncertainty analysis using a physically-based Monte Carlo modeling of pesticide fate and transport provides useful information for the evaluation of chemical leaching in Hawaii soils.

Fine scale spatial variability of microbial pesticide degradation in soil: scales, controlling factors, and implications

Pesticide biodegradation is a soil microbial function of critical importance for modern agriculture and its environmental impact. While it was once assumed that this activity was homogeneously distributed at the field scale, mounting evidence indicates that this is rarely the case. Here, we critically examine the literature on spatial variability of pesticide biodegradation in agricultural soil. We discuss the motivations, methods, and main findings of the primary literature. We found significant diversity in the approaches used to describe and quantify spatial heterogeneity, which complicates inter-studies comparisons. However, it is clear that the presence and activity of pesticide degraders is often highly spatially variable with coefficients of variation often exceeding 50% and frequently displays non-random spatial patterns. A few controlling factors have tentatively been identified across pesticide classes: they include some soil characteristics (pH) and some agricultural management practices (pesticide application, tillage), while other potential controlling factors have more conflicting effects depending on the site or the pesticide. Evidence demonstrating the importance of spatial heterogeneity on the fate of pesticides in soil has been difficult to obtain but modeling and experimental systems that do not include soil's full complexity reveal that this heterogeneity must be considered to improve prediction of pesticide biodegradation rates or of leaching risks. Overall, studying the spatial heterogeneity of pesticide biodegradation is a relatively new field at the interface of agronomy, microbial ecology, and geosciences and a wealth of novel data is being collected from these different disciplinary perspectives. We make suggestions on possible avenues to take full advantage of these investigations for a better understanding and prediction of the fate of pesticides in soil.

Contaminant transport in soil with depth-dependent reaction coefficients and time-dependent boundary conditions

Predicting the fate and movement of contaminant in soils and groundwater is essential to assess and reduce the risk of soil contamination and groundwater pollution. Reaction processes of contaminant often decreased monotonously with depth. Time-dependent input sources usually occurred at the inlet of natural or human-made system such as radioactive waste disposal site. This study presented a one-dimensional convection-dispersion equation (CDE) for contaminant transport in soils with depth-dependent reaction coefficients and time-dependent inlet boundary conditions, and derived its analytical solution. The adsorption coefficient and degradation rate were represented as sigmoidal functions of soil depth. Solute breakthrough curves (BTCs) and concentration profiles obtained from CDE with depth-dependent and constant reaction coefficients were compared, and a constant effective reaction coefficient, which was calculated by arithmetically averaging the depth-dependent reaction coefficient, was proposed to reflect the lumped depth-dependent reaction effect. With the effective adsorption coefficient and degradation rate, CDE could produce similar BTCs and concentration profiles as those from CDE with depth-dependent reactions in soils with moderate chemical heterogeneity. In contrast, the predicted concentrations of CDE with fitted reaction coefficients at a certain depth departed significantly from those of CDE with depth-dependent reactions. Parametric analysis was performed to illustrate the effects of sinusoidally and exponentially decaying input functions on solute BTCs. The BTCs and concentration profiles obtained from the solutions for finite and semi-infinite domain were compared to investigate the effects of effluent boundary condition. The finite solution produced higher concentrations at the increasing limb of the BTCs and possessed a higher peak concentration than the semi-infinite solution which had a slightly long tail. Furthermore, the finite solution gave a higher concentration in the immediate vicinity of the exit boundary than the semi-infinite solution. The applicability of the proposed model was tested with a field herbicide and tracer leaching experiment in an agricultural area of northeastern Greece. The simulation results indicated that the proposed CDE with depth-dependent reaction coefficients was able to capture the evolution of metolachlor concentration at the upper soil depths. However, the simulation results at deep depths were not satisfactory as the proposed model did not account for preferential flow observed in the field.

Field leaching of pesticides at five test sites in Hawaii: modeling flow and transport

BACKGROUND

Physically based tier-II models may serve as possible alternatives to expensive field and laboratory leaching experiments required for pesticide approval and registration. The objective of this study was to predict pesticide fate and transport at five different sites in Hawaii using data from an earlier field leaching experiment and a one-dimensional tier-II model. As the predicted concentration profiles of pesticides did not provide close agreement with data, inverse modeling was used to obtain adequate reactive transport parameters. The estimated transport parameters of pesticides were also utilized in a tier-I model, which is currently used by the state authorities to evaluate the relative leaching potential.

RESULTS

Water flow in soil profiles was simulated by the tier-II model with acceptable accuracy at all experimental sites. The observed concentration profiles and center of mass depths predicted by the tier-II simulations based on optimized transport parameters provided better agreements than did the non-optimized parameters. With optimized parameters, the tier-I model also delivered results consistent with observed pesticide center of mass depths.

CONCLUSION

Tier-II numerical modeling helped to identify relevant transport processes in field leaching of pesticides. The process-based modeling of water flow and pesticide transport, coupled with the inverse procedure, can contribute significantly to the evaluation of chemical leaching in Hawaii soils.

Leaching standards for mineral recycling materials--a harmonized regulatory concept for the upcoming German Recycling Decree

In this contribution we give a first general overview of results of recent studies in Germany which focused on contaminant leaching from various materials and reactive solute transport in the unsaturated soil zone to identify the key factors for groundwater risk assessment. Based on these results we developed new and improved existing methods for groundwater risk assessment which are used to derive a new regulatory concept for the upcoming "Decree for the Requirements of the Use of Alternative Mineral Building Materials in Technical Constructions and for the Amendment of the Federal Soil Protection and Contaminated Sites Ordinance" of the German Federal Ministry of Environment. The new concept aims at a holistic and scientifically sound assessment of the use of mineral recycling materials (e.g., mineral waste, excavated soils, slag and ashes, recycling products, etc.) in technical constructions (e.g., road dams) and permanent applications (e.g., backfilling and landscaping) which is based on a mechanistic understanding of leaching and transport processes. Fundamental for risk assessment are leaching standards for the mineral recycling materials. For each application of mineral recycling materials specific maximum concentrations of a substance in the seepage water at the bottom of an application were calculated. Technical boundary conditions and policy conventions derived from the "German precautionary groundwater and soil protection policy" were accounted to prevent adverse environmental effects on the media soil and groundwater. This includes the concentration decline of highly soluble substances (e.g., chloride and sulphate), retardation or attenuation of solutes, accumulation of contaminants in sub-soils and the hydraulic properties of recycling materials used for specific applications. To decide whether the use of a mineral recycling material is possible in a specific application, the leaching qualities were evaluated based on column percolation tests with various samples and compared with application-specific maximum concentrations. In the upcoming federal decree this simplified concept is realized using detailed tables which classify the leaching quality of mineral recycling materials and demonstrate potential application. A quality assurance system will be mandatory which defines specific testing programs (material properties and limit concentrations to be tested, number and schedule of testing) for the different mineral recycling materials using standardized methods (column percolation test).

Simulation of pharmaceutical and personal care product transport to tile drains after biosolids application

Pharmaceuticals and personal care products (PPCPs) carried in biosolids may reach surface waters or ground water when these materials are applied as fertilizer to agricultural land. During preferential flow conditions created by land application of liquid municipal biosolids (LMB), the residence time of solutes in the macropores may be too short for sorption equilibration. The physically based dual-permeability model MACRO is used in environmental risk assessments for pesticides and may have potential as an environmental risk assessment tool for PPCPs. The objective of this study was to evaluate MACRO and an updated version of MACRO that included non-equilibrium sorption in macropores using data from experiments conducted in eastern Ontario, Canada on the transport of three PPCPs (atenolol, carbamazepine, and triclosan), the nicotine metabolite cotinine, and the strongly sorbing dye rhodamine WT applied in LMB. Results showed that the MACRO model could not reproduce the measured rhodamine WT concentrations (Nash-Sutcliffe coefficient [NS] for the best simulation = -0.057) in drain discharge. The updated version resulted in better fits to measured data for PPCP (average NS = 0.97) and rhodamine WT (NS = 0.84) concentrations. However, it was not possible to simulate all compounds using the same set of hydraulic parameters, which indicates that the model does not fully account for all relevant processes. The results presented herein show that non-equilibrium sorption in macropores has a large impact on simulated solute transport for reactive compounds contained in LMB. This process should be considered in solute transport models that are used for environmental risk assessments for such compounds.

A review of model applications for structured soils: b) Pesticide transport

The past decade has seen considerable progress in the development of models simulating pesticide transport in structured soils subject to preferential flow (PF). Most PF pesticide transport models are based on the two-region concept and usually assume one (vertical) dimensional flow and transport. Stochastic parameter sets are sometimes used to account for the effects of spatial variability at the field scale. In the past decade, PF pesticide models were also coupled with Geographical Information Systems (GIS) and groundwater flow models for application at the catchment and larger regional scales. A review of PF pesticide model applications reveals that the principal difficulty of their application is still the appropriate parameterization of PF and pesticide processes. Experimental solution strategies involve improving measurement techniques and experimental designs. Model strategies aim at enhancing process descriptions, studying parameter sensitivity, uncertainty, inverse parameter identification, model calibration, and effects of spatial variability, as well as generating model emulators and databases. Model comparison studies demonstrated that, after calibration, PF pesticide models clearly outperform chromatographic models for structured soils. Considering nonlinear and kinetic sorption reactions further enhanced the pesticide transport description. However, inverse techniques combined with typically available experimental data are often limited in their ability to simultaneously identify parameters for describing PF, sorption, degradation and other processes. On the other hand, the predictive capacity of uncalibrated PF pesticide models currently allows at best an approximate (order-of-magnitude) estimation of concentrations. Moreover, models should target the entire soil-plant-atmosphere system, including often neglected above-ground processes such as pesticide volatilization, interception, sorption to plant residues, root uptake, and losses by runoff. The conclusions compile progress, problems, and future research choices for modelling pesticide displacement in structured soils.

Dual-continuum analysis of a cadmium tracer field experiment

Cadmium penetration into a sandy-loam soil during the field ponded infiltration at Kralovska luka site in Southern Slovakia was observed in a controlled experiment. Adsorption of cadmium was examined using the radioactive tracer techniques in laboratory batch tests. Quite deep (115m)Cd penetration during the field experiment (65 cm beneath the soil surface) gives evidence of the particle-facilitated transport of cadmium through preferential pathways. More than 40% of the applied cadmium moved deeper than 10 cm. The field experiment was analyzed using both conventional single-continuum approach (based on Richards' equation for water flow and advection-dispersion equation for Cd transport) and dual-continuum approach. The latter approach allowed us to simulate more adequately the observed movement of cadmium.

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.