Pesticides and groundwater quality protection

A simple approach for ranking the leaching of pesticides from surface soil is presented and tentatively calibrated with field data from an agricultural area. The approach is based on the calculation of a leaching index indicating the proportion of active ingredient, with respect to the quantity applied, leaching from a soil model in a given time interval (one year). In the selected area, 85 wells tapping an unconfined aquifer were sampled for groundwater pesticide residue analysis, in order to explore the index region between leachers and nonleachers.     

Controls on atrazine leaching through a soil-unsaturated fractured limestone sequence at Brévilles, France

The objective of this study was to identify the main controls on atrazine leaching through luvisols and calcisols overlying fissured limestone using the dual-permeability model MACRO. The model parameterisation was based on a combination of direct measurements (e.g. hydraulic properties, adsorption and degradation), literature data and calibration against bromide leaching experiments in field plots. A Monte Carlo sensitivity analysis was carried out for a typical application pattern, considering two different depths of unsaturated limestone (15 and 30 m). MACRO calibrations to the field experiments demonstrated the occurrence of strong macropore flow in the luvisol, while transport in the calcisol could be described by the advection-dispersion equation. MACRO simulations of tritium and atrazine leaching qualitatively matched tritium concentration profiles measured in the limestone and atrazine concentrations measured in piezometers and in aquifer discharge via a spring. The sensitivity analysis suggested that the thickness of the limestone, as well as the transport properties and processes occurring in the unsaturated rock (e.g. matrix vs. fissure flow) will have little significant long-term effect on atrazine leaching, mainly because degradation is very slow in the limestone. No mineralization of atrazine was detected in one-year incubations and a mean half-life of 10 years was assumed in the simulations. Instead, processes occurring in the soil exerted the main control on predicted atrazine leaching, especially variations in the degradation rate and the strength of sorption and macropore flow. However, fissure flow in unsaturated rock is expected to exert a much more significant control on groundwater contamination for compounds that degrade more readily in the deep vadose zone.     

Exposure assessment of pesticides in a shallow groundwater of the Tagus vulnerable zone (Portugal): a multivariate statistical approach (JCA)

Purpose

To assess groundwater exposure to pesticides, in agricultural areas of ??Ribatejo?? region (Portugal), and the influence of some key factors in that exposure, field, laboratory and modelling studies were carried out.

Methods

The study was performed in maize, potato, sugar beet, tomato and vegetables agricultural areas, located in a shallow aquifer, with pesticides use and, in most cases, with irrigation practices. Pesticides used in the studied agricultural areas and having leaching potential were selected, being considered also other pesticides included in priority lists, defined in Europe. Evaluation of groundwater exposure to pesticides was carried out by successively: (1) groundwater sampling in seven campaigns over the period 2004?C2006; (2) pesticide analysis [including isolation and concentration from the groundwater samples and further determination by gas chromatography?Cmass spectrometry (GC?CMS) of 14 herbicides, four insecticides and two metabolites]; and (3) analysis and discussion of the results by applying joint correspondence analysis (JCA).

Results

From the 20 pesticides and metabolites selected for the study, 11 were found in groundwater. Pesticides and metabolites most frequently detected were atrazine, alachlor, metolachlor, desethylatrazine, ethofumesate, ??-endosulfan, metribuzine, lindane and ??-endosulfan. The results showed that groundwater exposure to pesticides is influenced by local factors??either environmental or agricultural, as precipitation, soil, geology, crops and irrigation practices. Spring and autumn were more associated with the detection of pesticides being more likely to observe mixtures of these compounds in a groundwater sample in these transition seasons.

Conclusions

This work evidences the importance of models, which evaluate pesticides environmental behaviour, namely their water contamination potential (as Mackay multicompartimental fugacity model) and, specially, groundwater contamination potential (as GUS and Bacci and Gaggi leaching indices), in pesticide selection. Moreover, it reveals the importance to adapt proper statistical methods according to level of left-censored data. Using JCA was still possible to establish relations between pesticides and their temporal trend in a case study where there were more than 80% of data censored. This study will contribute to the Tagus river basin management plan with information on the patterns of pesticide occurrence in the alluvial aquifer system.     

Modelling of physical and reactive processes during biodegradation of a hydrocarbon plume under transient groundwater flow conditions

Numerical experiments of non-reactive and reactive transport were carried out to quantify the influence of a seasonally varying, transient flow field on transport and natural attenuation at a hydrocarbon-contaminated field site. Different numerical schemes for solving advective transport were compared to assess their capability to model low transversal dispersivities in transient flow fields. For the field site, it is shown that vertical plume spreading is largely inhibited, particularly if sorption is taken into account. For the reactive simulations, a biodegradation reaction module for the geochemical transport model PHT3D was developed. Results of the reactive transport simulations show that under the site-specific conditions the temporal variations in groundwater flow do, to a modest extent, affect average biodegradation rates and average total (dissolved) contaminant mass in the aquifer. The model simulations demonstrate that the seasonal variability in groundwater flow only results in significantly enhanced biodegradation rates when a differential sorption of electron donor (toluene) and electron acceptor (sulfate) is assumed.     

Influence of input uncertainty on prediction of within-field pesticide leaching risks

Previous research has suggested that pesticide losses at the field scale can be dominated by a small proportion of the field area. The objective of this study was to investigate whether site-specific applications (i.e. avoiding high-risk areas) at the field scale can contribute to a reduction of pesticide leaching despite uncertainty in the underlying model-based leaching risk map. Using a meta-model of the dual-permeability model MACRO, the annual average pesticide leaching concentrations were estimated for 162 sample sites on a 47 ha field. The procedure was repeated for different scenarios describing different patterns of spatial variation of degradation half-lives and the partition coefficient to soil organic carbon. To account for interpolation uncertainty, maps of predicted pesticide leaching risk were produced by the method of sequential Gaussian simulation. The results of the case study show that larger reductions of predicted leaching were achieved by site-specific application than by that of a comparable uniform dose reduction. Hence, site-specific-applications may be a feasible method to reduce pesticide leaching at the field-scale providing that the model approach gives reasonable estimates of the spatial pattern of pesticide leaching.     

Predicting natural attenuation of xylene in groundwater using a numerical model.

The aquifer beneath an abandoned refinery in the Lower Rhine area, Germany, was contaminated with a number of different mineral oil products. Groundwater sampling in the area around the former xylene plant revealed that a xylene plume had developed in the underlying groundwater, and moreover, that there is strong evidence for in situ microbial xylene degradation with oxygen, nitrate, sulfate and ferric iron as electron acceptors. In order to prevent further xylene spreading, three pumping wells extracting contaminated water were installed downgradient of the spill zone. The numerical reactive transport code Transport Biochemisty Chemistry (TBC) was applied to this situation to quantify the relation of microbial degradation to xylene removal by the pumping wells. It could be shown that the unamended in situ degradation was an appreciable xylene removal process that contributed to about one-third to the total xylene removal (degradation plus extraction). A further objective of the model application was to predict xylene spreading under regional flow conditions, i.e. without operation of the three pumping wells, to consider the possible effects of natural xylene attenuation. To accomplish this, the model calibrated for the situation with operating wells was transferred to the hydraulic situation of regional flow while retaining the parameters of the biochemical model. It turned out that the xylene plume that is expected to develop downgradient of the source area will be limited to an extension of not more than 1000 m. An interesting feature of the simulations results was that xylene degradation under iron-reducing conditions, which was of minor importance for the situation with operating pumping wells, becomes the dominant degradation mechanism under regional flow conditions. Moreover, iron reduction will be the key process in controlling plume evolution. The model application illustrates that multi-species reactive transport models are needed to adequately transfer reactive processes from one hydraulic situation to another, while single species models are not suited for this predictive task.     

Area-differentiated modelling of excess nutrient and contaminant leaching risks at the river-basin scale using a hydrological approach in India

Application of excess nutrients, such as fertilizers, is a significant and sometimes even major component of groundwater pollution. Diffuse inputs of nutrients and contaminants to the groundwater are related to runoff generated from precipitation on a catchment. This implies that the analysis of diffuse fluxes from the land surface to the groundwater requires an analysis of water fluxes for a catchment. This requires the simulation and modelling of total runoff, groundwater recharge, and plant-available water as a function of the regional interaction of the climate, soil, hydrogeology, topography and land-use conditions in the river basin. A model has been developed for large river basins in India, and has been applied to the Upper Yamuna basin, to quantify the exchange probability of plant-available soil water, which can be taken as a measure to determine the nutrient and contaminant leaching risk of a site. It was found that, with the available large-scale databases and methods, regional patterns of the total runoff could be simulated successfully. In this way, about 75001km² of the total 121000 km² of the Upper Yamuna basin was classified as an area sensitive to nutrient and contaminant leaching.     

Vertical small scale variations of sorption and mineralization of three herbicides in subsurface limestone and sandy aquifer

Vertical variation in sorption and mineralization potential of mecoprop (MCPP), isoproturon and acetochlor were investigated at low concentrations (μg-range) at the cm-scale in unsaturated sub-surface limestone samples and saturated sandy aquifer samples from an agricultural catchment in Brévilles, France. From two intact core drills, four heterogenic limestone sections were collected from 4.50 to 26.40m below surface (mbs) and divided into 12 sub-samples of 8-25cm length, and one sandy aquifer section from 19.20 to 19.53m depth divided into 7 sub-samples of 4-5cm length. In the sandy aquifer section acetochlor and isoproturon sorption increased substantially with depth; in average 78% (acetochlor) and 61% (isoproturon) per 5cm. Also the number of acetochlor and isoproturon degraders (most-probable-number) was higher in the bottom half of the aquifer section (93->16000/g) than in the upper half (4-71/g). One 50cm long limestone section with a distinct shift in color showed a clear shift in mineralization, number of degraders and sorption: In the two brown, uppermost samples, up to 31% mecoprop and up to 9% isoproturon was mineralized during 231 days, the numbers of mecoprop and isoproturon degraders were 1300 to >16000/g, and the sorption of both isoproturon and acetochlor was more than three times higher, compared to the two deeper, grayish samples just below where mineralization (≤4%) and numbers of degraders (1-520/g) were low for all three herbicides. In both unsaturated limestone and sandy aquifer, variations and even distinct shifts in both mineralization, number of specific degraders and sorption were seen within just 4-15cm of vertical distance. A simple conceptual model of herbicides leaching to groundwater through a 10m unsaturated limestone was established, and calculations showed that a 30cm active layer with the measured sorption and mineralization values hardly impacted the fate of the investigated herbicides, whereas a total thickness of layers of 1m would substantially increase natural attenuation.     

Stochastic analysis to assess the spatial distribution of groundwater nitrate concentrations in the Po catchment (Italy)

A large database including temporal trends of physical, ecological and socio-economic data was developed within the EUROCAT project. The aim was to estimate the nutrient fluxes for different socio-economic scenarios at catchment and coastal zone level of the Po catchment (Northern Italy) with reference to the Water Quality Objectives reported in the Water Framework Directive (WFD 2000/60/CE) and also in Italian legislation. Emission data derived from different sources at national, regional and local levels are referred to point and non-point sources. While non-point (diffuse) sources are simply integrated into the nutrient flux model, point sources are irregularly distributed. Intensive farming activity in the Po valley is one of the main Pressure factors Driving groundwater pollution in the catchment, therefore understanding the spatial variability of groundwater nitrate concentrations is a critical issue to be considered in developing a Water Quality Management Plan. In order to use the scattered point source data as input in our biogeochemical and transport models, it was necessary to predict their values and associated uncertainty at unsampled locations. This study reports the spatial distribution and uncertainty of groundwater nitrate concentration at a test site of the Po watershed using a probabilistic approach. Our approach was based on geostatistical sequential Gaussian simulation used to yield a series of stochastic images characterized by equally probable spatial distributions of the nitrate concentration across the area. Post-processing of many simulations allowed the mapping of contaminated and uncontaminated areas and provided a model for the uncertainty in the spatial distribution of nitrate concentrations.     

Dynamic factor analysis of groundwater quality trends in an agricultural area adjacent to Everglades National Park

The extensive eastern boundary of Everglades National Park (ENP) in south Florida (USA) is subject to one of the most expensive and ambitious environmental restoration projects in history. Understanding and predicting the water quality interactions between the shallow aquifer and surface water is a key component in meeting current environmental regulations and fine-tuning ENP wetland restoration while still maintaining flood protection for the adjacent developed areas. Dynamic factor analysis (DFA), a recent technique for the study of multivariate non-stationary time-series, was applied to study fluctuations in groundwater quality in the area. More than two years of hydrological and water quality time series (rainfall; water table depth; and soil, ground and surface water concentrations of N-NO3-, N-NH4+, P-PO4(3-), Total P, F-and Cl-) from a small agricultural watershed adjacent to the ENP were selected for the study. The unexplained variability required for determining the concentration of each chemical in the 16 wells was greatly reduced by including in the analysis some of the observed time series as explanatory variables (rainfall, water table depth, and soil and canal water chemical concentration). DFA results showed that groundwater concentration of three of the agrochemical species studied (N-NO3-, P-PO4(3-)and Total P) were affected by the same explanatory variables (water table depth, enriched topsoil, and occurrence of a leaching rainfall event, in order of decreasing relative importance). This indicates that leaching by rainfall is the main mechanism explaining concentration peaks in groundwater. In the case of N-NH4+, in addition to leaching, groundwater concentration is governed by lateral exchange with canals. F-and Cl- are mainly affected by periods of dilution by rainfall recharge, and by exchange with the canals. The unstructured nature of the common trends found suggests that these are related to the complex spatially and temporally varying land use patterns in the watershed. The results indicate that peak concentrations of agrochemicals in groundwater could be reduced by improving fertilization practices (by splitting and modifying timing of applications) and by operating the regional canal system to maintain the water table low, especially during the rainy periods.     

Modelling the closure-related geochemical evolution of groundwater at a former uranium mine

A newly developed reactive transport model was used to evaluate the potential effects of mine closure on the geochemical evolution in the aquifer downgradient from a mine site. The simulations were conducted for the K?nigstein uranium mine located in Saxony, Germany. During decades of operation, uranium at the former mine site had been extracted by in situ acid leaching of the ore underground, while the mine was maintained in a dewatered condition. One option for decommissioning is to allow the groundwater level to rise to its natural level, flooding the mine workings. As a result, pore water containing high concentrations of dissolved metals, radionuclides, and sulfate may be released. Additional contamination may arise due to the dissolution of minerals contained in the aquifer downgradient of the mine. On the other hand, dissolved metals may be attenuated by reactions within the aquifer. The geochemical processes and interactions involved are highly non-linear and their impact on the quality of the groundwater and surface water downstream of the mine is not always intuitive. The multicomponent reactive transport model MIN3P, which can describe mineral dissolution-precipitation reactions, aqueous complexation, and oxidation-reduction reactions, is shown to be a powerful tool for investigating these processes. The predictive capabilities of the model are, however, limited by the availability of key geochemical parameters such as the presence and quantities of primary and secondary mineral phases. Under these conditions, the model can provide valuable insight by means of sensitivity analyses.     

Influence of mass transfer characteristics for DNAPL source depletion and contaminant flux in a highly characterized glaciofluvial aquifer

The transfer of contaminant mass between the nonaqueous- and aqueous-phases is a process of central importance for the remediation of sites contaminated by dense nonaqueous-phase liquids (DNAPLs). This paper describes a comparison of the results obtained with various alternative DNAPL-aqueous-phase mass transfer models contained in the literature for predicting DNAPL source-zone depletion times in groundwater systems. These dissolution models were largely developed through laboratory column experiments. To gain insight into the implications of various representations of the local-scale kinetic as well as equilibrium DNAPL dissolution processes, aquifer heterogeneity and the complex architecture of a DNAPL source-zone, the aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are analyzed in a conditional stochastic framework. The hydrogeologic setting is a heterogeneous fluvial aquifer in Southwest Germany, referred to as the aquifer analog dataset, that was intensively characterized in three dimensions for hydrogeological parameters that include permeability, effective porosity, grain size, mineralogy and sorption coefficients. By embedding the various dissolution models into the compositional, multiphase flow model, CompFlow, the relative times predicted for complete depletion of a released DNAPL source due to natural dissolution are explored. Issues related to achieving environmental benefits through, for example, partial DNAPL-zone source removal via enhanced remedial technologies are also discussed. In this context, performance metrics in the form of peak aqueous-phase contaminant concentrations and mass fluxes arriving at a down-gradient compliance boundary are compared to each other. This is done for each of the alternative mass transfer models. A significant reduction in the fractional flux at a downstream location from the DNAPL source can be achieved by partial source-zone mass reduction; however, peak concentration levels at the same location remain much higher than the United States Environment Protection Agency (US-EPA) drinking water limits. Although groundwater quality was found to improve more rapidly for the equilibrium dissolution model, it is also shown that dissolution models that promote rapid DNAPL disappearance produce greater prediction uncertainty in the aqueous-phase flux reduction.     

Calibration and validation of an integrated nitrate transport model within a well capture zone

Groundwater contamination by nitrate was investigated in an agricultural area in southern Quebec, Canada, where a municipal well is the local source of drinking water. A network of 38 piezometers was installed within the capture zone of the municipal well to monitor water table levels and nitrate concentrations in the aquifer. Nitrate concentrations were also measured in the municipal well. A Water flow and Nitrate transport Global Model (WNGM) was developed to simulate the impact of agricultural activities on nitrate concentrations in both the aquifer and municipal well. The WNGM first uses the Agriflux model to simulate vertical water and nitrate fluxes below the root zone for each of the seventy agricultural fields located within the capture zone of the municipal well. The WNGM then uses the HydroGeoSphere model to simulate three-dimensional variably-saturated groundwater flow and nitrate transport in the aquifer using water and nitrate fluxes computed with the Agriflux model as the top boundary conditions. The WNGM model was calibrated by reproducing water levels measured from 2005 to 2007 in the network of piezometers and nitrate concentrations measured in the municipal well from 1997 to 2007. The nitrate concentrations measured in the network of piezometers, however, showed greater variability than in the municipal well and could not be reproduced by the calibrated model. After calibration, the model was validated by successfully reproducing the decrease of nitrate concentrations observed in the municipal well in 2006 and 2007. Although it cannot predict nitrate concentrations in individual piezometers, the calibrated and validated WNGM can be used to assess the impact of changes in agricultural practices on global nitrate concentrations in the aquifer and in the municipal well.     

Groundwater flow and contaminant transport modeling applications in urban area: scopes and limitations

Background

A three-dimensional groundwater flow model was used to evaluate the groundwater potential and assess the effects of groundwater withdrawal on the regional water level and flow direction in the central Beijing area. A program of groundwater modeling aimed at estimating current contaminant fluxes to the central area and site streams via groundwater was developed.

Results and discussion

The conceptual model developed for the site attempted to incorporate a complex stratigraphic profile in which groundwater flow and contaminant transport is strongly controlled by a shallow aquifer. Here, a conceptual model for groundwater flow and contaminant transport in central Beijing is presented.

Conclusion

Model simulations indicated that a sharp drop in the hydraulic head occurs at the center of the model area, which generates a cone of depression and a continuous decline of head with respect to time as a result of heavy groundwater abstraction.     

Incorporating uncertainty in a chemical leaching assessment

Pesticide leaching models are being used to assist in the regulation and management of pesticides by indicating their potential for leaching to groundwater. Uncertainty in model input data is not, regrettably, included in most pesticide leaching assessments. In the work described here, we use logarithmic transformations of the attenuation factor (AF), a simple process-based index model, to represent uncertainty in a pesticide leaching assessment. Characterization of a wide range of pesticides as `leachers' or `non-leachers' for a specific Hawaii hydrogeological setting is facilitated by comparing the log-transformed AF, designated AFR, for each chemical with two reference chemicals for which leaching behavior in Hawaii is known. Defining a mean and uncertainty interval for the AFR index of each chemical being ranked provides a practical method of incorporating data uncertainty into a regulatory protocol.     

垃圾填埋场对地下水污染的模拟研究

在对长春市某垃圾填埋场进行野外调查的基础上 ,于室内进行了垃圾淋滤模拟实验。分析了垃圾渗滤液污染组分的自然衰减规律 ,建立了垃圾填埋场地下水污染的数值模型 ,采用FEFLOW软件对其进行模拟和预报 ,并取得了较好的效果 ,最后提出防止和防治垃圾填埋场污染地下水的若干措施。     

Direct coupling of a genome-scale microbial in silico model and a groundwater reactive transport model

The activity of microorganisms often plays an important role in dynamic natural attenuation or engineered bioremediation of subsurface contaminants, such as chlorinated solvents, metals, and radionuclides. To evaluate and/or design bioremediated systems, quantitative reactive transport models are needed. State-of-the-art reactive transport models often ignore the microbial effects or simulate the microbial effects with static growth yield and constant reaction rate parameters over simulated conditions, while in reality microorganisms can dynamically modify their functionality (such as utilization of alternative respiratory pathways) in response to spatial and temporal variations in environmental conditions. Constraint-based genome-scale microbial in silico models, using genomic data and multiple-pathway reaction networks, have been shown to be able to simulate transient metabolism of some well studied microorganisms and identify growth rate, substrate uptake rates, and byproduct rates under different growth conditions. These rates can be identified and used to replace specific microbially-mediated reaction rates in a reactive transport model using local geochemical conditions as constraints. We previously demonstrated the potential utility of integrating a constraint-based microbial metabolism model with a reactive transport simulator as applied to bioremediation of uranium in groundwater. However, that work relied on an indirect coupling approach that was effective for initial demonstration but may not be extensible to more complex problems that are of significant interest (e.g., communities of microbial species and multiple constraining variables). Here, we extend that work by presenting and demonstrating a method of directly integrating a reactive transport model (FORTRAN code) with constraint-based in silico models solved with IBM ILOG CPLEX linear optimizer base system (C library). The models were integrated with BABEL, a language interoperability tool. The modeling system is designed in such a way that constraint-based models targeting different microorganisms or competing organism communities can be easily plugged into the system. Constraint-based modeling is very costly given the size of a genome-scale reaction network. To save computation time, a binary tree is traversed to examine the concentration and solution pool generated during the simulation in order to decide whether the constraint-based model should be called. We also show preliminary results from the integrated model including a comparison of the direct and indirect coupling approaches and evaluated the ability of the approach to simulate field experiment.     

Contaminant mass discharge estimation in groundwater based on multi-level point measurements: a numerical evaluation of expected errors

Different methods for the field-scale estimation of contaminant mass discharge in groundwater at control planes based on multi-level well data are numerically analysed for the expected estimation error. We consider "direct" methods based on time-integrated measuring of mass flux, as well as "indirect" methods, where estimates are derived from concentration measurements. The appropriateness of the methods is evaluated by means of modelled data provided by simulation of mass transport in a three-dimensional model domain. Uncertain heterogeneous aquifer conditions are addressed by means of Monte-Carlo simulations with aquifer conductivity as a random space function. We investigate extensively the role of the interplay between the spatial resolution of the sampling grid and aquifer heterogeneity with respect to the accuracy of the mass discharge estimation. It is shown that estimation errors can be reduced only if spatial sampling intervals are in due proportion to spatial correlation length scales. The ranking of the methods with regard to estimation error is shown to be heavily dependent on both the given sampling resolution and prevailing aquifer heterogeneity. Regarding the "indirect" estimation methods, we demonstrate the great importance of a consistent averaging of the parameters used for the discharge estimation.