Stochastic modeling of diffuse pesticide losses from a small agricultural catchment

The objective of this study was to identify the main sources of variation in pesticide losses at field and catchment scales using the dual permeability model MACRO. Stochastic simulations of the leaching of the herbicide MCPA (4-chloro-2-methylphenoxyacetic acid) were compared with seven years of measured concentrations in a stream draining a small agricultural catchment and one year of measured concentrations at the outlet of a field located within the catchment. MACRO was parameterized from measured probability distributions accounting for spatial variability of soil properties and local pedotransfer functions derived from information gathered in field- and catchment-scale soil surveys. At the field scale, a single deterministic simulation using the means of the input distributions was also performed. The deterministic run failed to reproduce the summer outflows when most leaching occurred, and greatly underestimated pesticide leaching. In contrast, the stochastic simulations successfully predicted the hydrologic response of the field and catchment and there was a good resemblance between the simulations and measured MCPA concentrations at the field outlet. At the catchment scale, the stochastic approach underestimated the concentrations of MCPA in the stream, probably mostly due to point sources, but perhaps also because the distributions used for the input variables did not accurately reflect conditions in the catchment. Sensitivity analyses showed that the most important factors affecting MACRO modeled diffuse MCPA losses from this catchment were soil properties controlling macropore flow, precipitation following application, and organic carbon content.     

Simulating solute transport in a structured field soil: uncertainty in parameter identification and predictions

Dual-permeability models have been developed to account for the significant effects of macropore flow on contaminant transport, but their use is hampered by difficulties in estimating the additional parameters required. Therefore, our objective was to evaluate data requirements for parameter identification for predictive modeling with the dual-permeability model MACRO. Two different approaches were compared: sequential uncertainty fitting (SUFI) and generalized likelihood uncertainty estimation (GLUE). We investigated six parameters controlling macropore flow and pesticide sorption and degradation, applying MACRO to a comprehensive field data set of bromide andbentazone [3-isopropyl-1H-2,1,3-benzothiadiazin-4(3H)-one-2,2dioxide] transport in a structured soil. The GLUE analyses of parameter conditioning for different combinations of observations showed that both resident and flux concentrations were needed to obtain highly conditioned and unbiased parameters and that observations of tracer transport generally improved the conditioning of macropore flow parameters. The GLUE "behavioral" parameter sets covered wider parameter ranges than the SUFI posterior uncertainty domains. Nevertheless, estimation uncertainty ranges defined by the 5th and 95th percentiles were similar and many simulations randomly sampled from the SUFI posterior uncertainty domains had negative model efficiencies (minimum of -3.2). This is because parameter correlations are neglected in SUFI and the posterior uncertainty domains were not always determined correctly. For the same reasons, uncertainty ranges for predictions of bentazone losses through drainflow for good agricultural practice in southern Sweden were 27% larger for SUFI compared with GLUE. Although SUFI proved to be an efficient parameter estimation tool, GLUE seems better suited as a method of uncertainty estimation for predictions.     

Sensitivity and first-step uncertainty analyses for the preferential flow model MACRO

Sensitivity analyses for the preferential flow model MACRO were carried out using one-at-a-time and Monte Carlo sampling approaches. Four different scenarios were generated by simulating leaching to depth of two hypothetical pesticides in a sandy loam and a more structured clay loam soil. Sensitivity of the model was assessed using the predictions for accumulated water percolated at a 1-m depth and accumulated pesticide losses in percolation. Results for simulated percolation were similar for the two soils. Predictions of water volumes percolated were found to be only marginally affected by changes in input parameters and the most influential parameter was the water content defining the boundary between micropores and macropores in this dual-porosity model. In contrast, predictions of pesticide losses were found to be dependent on the scenarios considered and to be significantly affected by variations in input parameters. In most scenarios, predictions for pesticide losses by MACRO were most influenced by parameters related to sorption and degradation. Under specific circumstances, pesticide losses can be largely affected by changes in hydrological properties of the soil. Since parameters were varied within ranges that approximated their uncertainty, a first-step assessment of uncertainty for the predictions of pesticide losses was possible. Large uncertainties in the predictions were reported, although these are likely to have been overestimated by considering a large number of input parameters in the exercise. It appears desirable that a probabilistic framework accounting for uncertainty is integrated into the estimation of pesticide exposure for regulatory purposes.     

Mapping ground water vulnerability to pesticide leaching with a process-based metamodel of EuroPEARL

To support EU policy, indicators of pesticide leaching at the European level are required. For this reason, a metamodel of the spatially distributed European pesticide leaching model EuroPEARL was developed. EuroPEARL considers transient flow and solute transport and assumes Freundlich adsorption, first-order degradation and passive plant uptake of pesticides. Physical parameters are depth dependent while (bio)-chemical parameters are depth, temperature, and moisture dependent. The metamodel is based on an analytical expression that describes the mass fraction of pesticide leached. The metamodel ignores vertical parameter variations and assumes steady flow. The calibration dataset was generated with EuroPEARL and consisted of approximately 60,000 simulations done for 56 pesticides with different half-lives and partitioning coefficients. The target variable was the 80th percentile of the annual average leaching concentration at 1-m depth from a time series of 20 yr. The metamodel explains over 90% of the variation of the original model with only four independent spatial attributes. These parameters are available in European soil and climate databases, so that the calibrated metamodel could be applied to generate maps of the predicted leaching concentration in the European Union. Maps generated with the metamodel showed a good similarity with the maps obtained with EuroPEARL, which was confirmed by means of quantitative performance indicators.     

Developing risk models of Cryptosporidium transport in soils from vegetated, tilted soilbox experiments

Transport of Cryptosporidium parvum through macroporous soils is poorly understood yet critical for assessing the risk of groundwater contamination. We developed a conceptual model of the physics of flow and transport in packed, tilted, and vegetated soilboxes during and immediately after a simulated rainfall event and applied it to 54 experiments implemented with different soils, slopes, and rainfall rates. Using a parsimonious inverse modeling procedure, we show that a significant amount of subsurface outflow from the soilboxes is due to macropore flow. The effective hydraulic properties of the macropore space were obtained by calibration of a simple two-domain flow and transport model that accounts for coupled flow in the matrix and in the macropores of the soils. Using linear mixed-effects analysis, macropore hydraulic properties and oocyst attenuation were shown to be associated with soil bulk density and rainfall rate. Macropore flow was shown to be responsible for bromide and C. parvum transport through the soil into the underlying pore space observed during the 4-h experiments. We confirmed this finding by conducting a pair of saturated soil column studies under homogeneously repacked conditions with no macropores in which no C. parvum transport was observed in the effluent. The linear mixed-effects and logistic regression models developed from the soilbox experiments provide a basis for estimating macropore hydraulic properties and the risk of C. parvum transport through shallow soils from bulk density, precipitation, and total shallow subsurface flow rate. The risk assessment is consistent with the reported occurrence of oocysts in springs or groundwater from fractured or karstic rocks protected only by shallow overlying soils.     

Modeling Parent and Metabolite Fate and Transport in Subsurface Drained Fields With Directly Connected Macropores1

Abstract: Few studies exist that evaluate or apply pesticide transport models based on measured parent and metabolite concentrations in fields with subsurface drainage. Furthermore, recent research suggests pesticide transport through exceedingly efficient direct connections, which occur when macropores are hydrologically connected to subsurface drains, but this connectivity has been simulated at only one field site in Allen County, Indiana. This research evaluates the Root Zone Water Quality Model (RZWQM) in simulating the transport of a parent compound and its metabolite at two subsurface drained field sites. Previous research used one of the field sites to test the original modification of the RZWQM to simulate directly connected macropores for bromide and the parent compound, but not for the metabolite. This research will evaluate RZWQM for parent/metabolite transformation and transport at this first field site, along with evaluating the model at an additional field site to evaluate whether the parameters for direct connectivity are transferable and whether model performance is consistent for the two field sites with unique soil, hydrologic, and environmental conditions. Isoxaflutole, the active ingredient in BALANCE^® herbicide, was applied to both fields. Isoxaflutole rapidly degrades into a metabolite (RPA 202248). This research used calibrated RZWQM models for each field based on observed subsurface drain flow and/or edge of field conservative tracer concentrations in subsurface flow. The calibrated models for both field sites required a portion (approximately 2% but this fraction may require calibration) of the available water and chemical in macropore flow to be routed directly into the subsurface drains to simulate peak concentrations in edge of field subsurface drain flow shortly after chemical applications. Confirming the results from the first field site, the existing modification for directly connected macropores continually failed to predict pesticide concentrations on the recession limbs of drainage hydrographs, suggesting that the current strategy only partially accounts for direct connectivity. Thirty‐year distributions of annual mass (drainage) loss of parent and metabolite in terms of percent of isoxaflutole applied suggested annual simulated percent losses of parent and metabolite (3.04 and 1.31%) no greater in drainage than losses in runoff on nondrained fields as reported in the literature.     

Numerical modeling of diazinon transport through inter-row vegetative filter strips

A numerical simulation model of pesticide runoff through vegetative filer strips (PRVFS) was developed as a tool for investigating the effects of pesticide transport mechanisms on VFS design in dormant-sprayed orchard. The PRVFS model was developed applying existing theories such as kinematic wave theory and mixing zone theory for pesticide transport in the bare soil area. For VFS area, the model performs flow routing by simple mass accounting in sequential segments and the pesticide mass balance by considering pesticide washoff and adsorption processes on the leaf, vegetative litter, root zone and soil. Model sensitivity analysis indicated that pesticide transfer from surface soil to overland flow and pesticide washoff from the VFS were important mechanisms affecting diazinon transport. The VFS cover ratio and rainfall intensity can be important design parameters for controlling diazinon runoff using inter-row VFS in orchard. The PRVFS model was validated using micro-ecosystem simulation of diazinon transport for 0, 50 and 100% VFS cover conditions. The PRVFS model is shown to be a beneficial tool for evaluating and analyzing possible best management practices for controlling offsite runoff of dormant-sprayed diazinon in orchards during the rainy season.     

Propagation of uncertainties in soil and pesticide properties to pesticide leaching

In the new Dutch decision tree for the evaluation of pesticide leaching to groundwater, spatially distributed soil data are used by the GeoPEARL model to calculate the 90th percentile of the spatial cumulative distribution function of the leaching concentration in the area of potential usage (SP90). Until now it was not known to what extent uncertainties in soil and pesticide properties propagate to spatially aggregated parameters like the SP90. A study was performed to quantify the uncertainties in soil and pesticide properties and to analyze their contribution to the uncertainty in SP90. First, uncertainties in the soil and pesticide properties were quantified. Next, a regular grid sample of points covering the whole of the agricultural area in the Netherlands was randomly selected. At the grid nodes, realizations from the probability distributions of the uncertain inputs were generated and used as input to a Monte Carlo uncertainty propagation analysis. The analysis showed that the uncertainty concerning the SP90 is 10 times smaller than the uncertainty about the leaching concentration at individual point locations. The parameters that contribute most to the uncertainty about the SP90 are, however, the same as the parameters that contribute most to uncertainty about the leaching concentration at individual point locations (e.g., the transformation half-life in soil and the coefficient of sorption on organic matter). Taking uncertainties in soil and pesticide properties into account further leads to a systematic increase of the predicted SP90. The important implication for pesticide regulation is that the leaching concentration is systematically underestimated when these uncertainties are ignored.     

Modeling organic carbon dynamics and cadmium fate in long-term sludge-amended soil

A model is described that may help to resolve uncertainty and controversy over the long-term consequences of sludge applications to arable land, especially with regard to the effects of sludge adsorption characteristics on trace metal solubility and bioavailability (e.g., the sludge "time bomb" or sludge "protection" hypotheses). Mass balances of organic and inorganic material derived from sludge and crop residues are simulated. Each pool has a potentially different adsorption affinity for trace metals, and this leads to changes in the adsorption capacity of sludge-amended soil that influence leaching and crop uptake. Model simulations were compared with measured changes in organic carbon and ethylenediaminetetraacetic acid (EDTA)-extractable cadmium contents in a clay loam soil following 41 years of sludge applications. The model adequately reproduced the data, although discrepancies in the vertical distribution of Cd were attributed to the effects of macropore transport and root-uptake driven recirculation. A Monte Carlo sensitivity analysis demonstrated that the most important parameters affecting leaching and crop uptake were the Cd loading and parameters controlling adsorption, especially the partition coefficient for sludge-derived inorganic material and the exponent regulating the effect of pH on sorption. Scenario simulations show that no general conclusions can be drawn with respect to the validity of the sludge "time bomb" and sludge "protection" hypotheses. Either may occur, or neither, depending on three key system parameters: the ratio of sludge adsorption capacity to the initial adsorption capacity of the soil, the proportion of the sludge adsorption capacity contributed by the inorganic fraction, and the sludge Cd loading.     

Simulating pesticide leaching and runoff in rice paddies with the RICEWQ-VADOFT model

There is a current need to simulate leaching and runoff of pesticide from rice (Oryza sativa L.) paddies for assessing environmental impacts on a valuable agricultural system. The objective of this study was to develop a model for determining predicted environmental concentration (PEC) in soil, runoff, and ground water through the linkage of two models, rice water quality model (RICEWQ) and vadose zone transport model (VADOFT), to simulate pesticide fate and transport within a rice paddy and underlying soil profile. Model performance was evaluated with a field data set obtained from a 2-yr field experiment in 1997 and 1998 in northern Italy. The predictions of amount of pesticide running off from the paddy field and accumulating in the paddy sediment were in agreement with measured values. Leaching into the vadose zone accounted for approximately 19% of the applied dose, but only a small amount of chemical (<0.1%) was predicted to reach ground water at a 5-m depth due to sorption and transformation in the soil. The permeability of the soil and the water management practices in the paddy field were shown to have a strong influence on pesticide fate. These factors need to be well characterized in the field if model predictions are to be successful. The combined model developed in this work is an effective tool for exposure assessments for soil, surface water, and ground water, in the particular conditions of rice cultivation.     

Simulated fate and transport of metolachlor in the unsaturated zone, Maryland, USA

An unsaturated-zone transport model was used to examine the transport and fate of metolachlor applied to an agricultural site in Maryland, USA. The study site was instrumented to collect data on soil-water content, soil-water potential, ground water levels, major ions, pesticides, and nutrients from the unsaturated zone during 2002-2004. The data set was enhanced with site-specific information describing weather, soils, and agricultural practices. The Root Zone Water Quality Model was used to simulate physical, chemical, and biological processes occurring in the unsaturated zone. Model calibration to bromide tracer concentrations indicated flow occurred through the soil matrix. Simulated recharge rates were within the measured range of values. The pesticide transport model was calibrated to the intensive data collection period (2002-2004), and the calibrated model was then used to simulate the period 1984 through 2004 to examine the impact of sustained agricultural management practices on the concentrations of metolachlor and its degradates at the study site. Simulation results indicated that metolachlor degrades rapidly in the root zone but that the degradates are transported to depth in measurable quantities. Simulations indicated that degradate transport is strongly related to the duration of sustained use of metolachlor and the extent of biodegradation.     

Interrelationship of macropores and subsurface drainage for conservative tracer and pesticide transport

Macropore flow results in the rapid movement of pesticides to subsurface drains, which may be caused in part by a small portion of macropores directly connected to drains. However, current models fail to account for this direct connection. This research investigated the interrelationship between macropore flow and subsurface drainage on conservative solute and pesticide transport using the Root Zone Water Quality Model (RZWQM). Potassium bromide tracer and isoxaflutole, the active ingredient in BALANCE herbicide [(5-cyclopropyl-4-isoxazolyl) [2(methylsulfonyl)-4-(trifluoromethyl)phenyl] methanone], with average half-life of 1.7 d were applied to a 30.4-ha Indiana corn (Zea mays L.) field. Water flow and chemical concentrations emanating from the drains were measured from two samplers. Model predictions of drain flow after minimal calibration reasonably matched observations (slope = 1.03, intercept = 0.01, and R(2) = 0.75). Without direct hydraulic connection of macropores to drains, RZWQM under predicted bromide and isoxaflutole concentration during the first measured peak after application (e.g., observed isoxaflutole concentration was between 1.2 and 1.4 mug L(-1), RZWQM concentration was 0.1 mug L(-1)). This research modified RZWQM to include an express fraction relating the percentage of macropores in direct hydraulic connection to drains. The modified model captured the first measured peak in bromide and isoxaflutole concentrations using an express fraction of 2% (e.g., simulated isoxaflutole concentration increased to 1.7 mug L(-1)). The RZWQM modified to include a macropore express fraction more accurately simulates chemical movement through macropores to subsurface drains. An express fraction is required to match peak concentrations in subsurface drains shortly after chemical applications.     

Transport of di(2-ethylhexyl)phthalate (DEHP) applied with sewage sludge to undisturbed and repacked soil columns

Municipal sewage sludge is often used on arable soils as a source of nitrogen and phosphorus, but it also contains organic contaminants that may be leached to the ground water. Di(2-ethylhexyl)phthalate (DEHP) is a priority pollutant that is present in sewage sludge in ubiquitous amounts. Column experiments were performed on undisturbed soil cores (20-cm depth x 20-cm diameter) with three different soil types: a sand, a loamy sand, and a sandy loam soil. Dewatered sewage sludge was spiked with 14C-labeled DEHP (60 mg kg(-1)) and bromide (5 g kg(-1)). Sludge was applied to the soil columns either as five aggregates, or homogeneously mixed with the surface layer. Also, two leaching experiments were performed with repacked soil columns (loamy sand and sandy loam soil). The DEHP concentrations in the effluent did not exceed 1.0 microg L(-1), and after 200 mm of outflow less than 0.5% of the applied amount was recovered in the leachate in all soils but the sandy loam soil with homogeneous sludge application (up to 3.4% of the applied amount recovered). In the absence of macropore flow, DEHP in the leachate was primarily sorbed to mobilized dissolved organic macromolecules (DOM, 30.3 to 81.3%), while 2.4 to 23.6% was sorbed to mobilized mineral particles. When macropore flow occurred, this changed to 16.5 to 37.4% (DOM) and 36.9 to 40.6% (mineral particles), respectively. The critical combination for leaching of considerable amounts of DEHP was homogeneous sludge application and a continuous macropore structure.     

Turfgrass thatch effects on pesticide leaching: a laboratory and modeling study

Process-based models are frequently used to assess the water quality impacts of turfgrass management emanating from proposed or existing golf courses. Thatch complicates the prediction of pesticide transport because surface-applied pesticides must pass through an organic-rich layer before entering the soil. This study was conducted to (i) compare the use of a linear equilibrium model (LEM) and two-site nonequilibrium (2SNE) model to predict pesticide transport through soil and thatch + soil columns, and (ii) evaluate thatch effects on pesticide transport through soil columns with a volume-averaging approach. Pesticide breakthrough curves were obtained for soil and thatch + soil columns from a 1 cm h(-1) flux applied one day after applying triclopyr (3,5,6-trichloro-2-pyridinyloxyacetic acid) and carbaryl (1-napthyl-methyl carbamate). Pesticide and bromide transport parameters indicated that nonequilibrium processes were affecting pesticide transport. Columns containing zoysiagrass (Zoysia japonica Steud.) thatch had lower triclopyr and carbaryl leaching losses than did soil-only columns, although total reductions attributable to thatch did not exceed 15% of the applied pesticide. When laboratory-based retardation factors were used, the 2SNE model explained 88 to 93% of the variability for triclopyr and 70 to 94% of the variability for carbaryl. Laboratory-based retardation factors performed well in a 2SNE model to predict the peak concentration and tailing behavior of triclopyr and carbaryl with a volume-averaging approach. These results suggest that separate representation of the thatch layer in process-based models is not a prerequisite to obtain reasonable estimates of pesticide transport under steady state flow conditions.     

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

Accurate input data for leaching models are expensive and difficult to obtain which may lead to the use of "general" non-site-specific input data. This study investigated the effect of using different quality data on model outputs. Three models of varying complexity, GLEAMS, LEACHM, and HYDRUS-2D, were used to simulate pesticide leaching at a field trial near Hamilton, New Zealand, on an allophanic silt loam using input data of varying quality. Each model was run for four different pesticides (hexazinone, procymidone, picloram and triclopyr); three different sets of pesticide sorption and degradation parameters (i.e., site optimized, laboratory derived, and sourced from the USDA Pesticide Properties Database); and three different sets of soil physical data of varying quality (i.e., site specific, regional database, and particle size distribution data). We found that the selection of site-optimized pesticide sorption (Koc) and degradation parameters (half-life), compared to the use of more general database derived values, had significantly more impact than the quality of the soil input data used, but interestingly also more impact than the choice of the models. Models run with pesticide sorption and degradation parameters derived from observed solute concentrations data provided simulation outputs with goodness-of-fit values closest to optimum, followed by laboratory-derived parameters, with the USDA parameters providing the least accurate simulations. In general, when using pesticide sorption and degradation parameters optimized from site solute concentrations, the more complex models (LEACHM and HYDRUS-2D) were more accurate. However, when using USDA database derived parameters, all models performed about equally.     

Phosphorus leaching under a restored tallgrass prairie and corn agroecosystems

Most studies of phosphorus (P) movement in soil have based their conclusions on patterns of extractable soil P as a function of depth, which has led to the assumption that no substantial leaching loss occurs because of high P-fixation capacity in mineral soils. Few studies have involved high-quality leachate samples collected below the root zone; rather, most have involved tile drainage systems. Equilibrium-tension lysimeters installed at a depth of 1.4 m were used to evaluate and compare P leaching from a restored tallgrass prairie and corn (Zea mays L.) agroecosystems on Plano silt loam soil (fine-silty, mixed, superactive, mesic Typic Argiudoll) in southcentral Wisconsin during a 5-yr period. The corn agroecosystem treatments included nitrogen (N)-fertilized (f) or N-unfertilized (nf) and no-tillage (NT) or chisel-plowed (CP). Mean volume-weighted molybdate-reactive phosphorus (MRP) and total dissolved phosphorus (TDP) concentrations were similar within replicate samples, but always higher in NTf corn than in the prairie or CPf corn systems, though drainage from the CPf corn was always higher than from the NTf corn system. Water-extractable soil P concentrations at any given depth were not positively correlated with leachate concentrations, suggesting that macropore flow causes infiltrating runoff to preferentially bypass the bulk of the soil matrix. Leachate-P concentrations from the natural and managed agroecosystems exceeded 0.01 mg P L(-1) and leaching losses were significantly higher from N-fertilized corn, regardless of tillage, than from the prairie or N-unfertilized corn systems, from which leachate-P concentrations and loads were similar. Increased root growth from N fertilization could cause more macropore formation, preferential flow, and P mineralization from decaying roots compared with N-unfertilized systems, which could contribute to a N-fertilization effect on P leaching.     

Partitioning and persistence of trichlorfon and chlorpyrifos in a creeping bentgrass putting green

Golf course putting greens typically receive high pesticide applications to meet high quality demands. Research on pesticide fate in turf ecosystems is important to better understand the potential impact of pesticide use on the environment and human health. This research was conducted to evaluate the environmental fate of two commonly used insecticides--trichlorfon (dimethyl 2,2,2-trichloro-1-hydroxyethylphosphonate) and chlorpyrifos (O,O-diethyl O-3,5,6-trichloro-2-pyridylphosphorothioate)--in a creeping bentgrass (Agrostis palustris Huds.) putting green under customary field management practices at the University of California-Riverside Turf Research Facility during 1996 and 1997. The two insecticides were chosen because of their difference in water solubility, persistence, adsorption, and vapor pressure. Volatilization, clipping removal, and soil residues of the insecticides were quantified and leaching was monitored using lysimeters installed in putting green plots. Results showed trichlorfon volatilization, clipping removal, and leaching loss was insignificant (in the range of 0.0001-0.06% of applied mass) both in 1996 and 1997. No significant difference in clipping removal of trichlorfon and chlorpyrifos was observed in both years (0.06 and 0.05% of applied mass for trichlorfon and 0.15 and 0.19% of applied mass for chlorpyrifos, respectively, in 1996 and 1997), but significantly lower cumulative leaching and lower soil concentration was observed in 1997 than in 1996. Volatilization loss of chlorpyrifos was not significantly different between 1996 (2.05%) and 1997 (2.71%). Volatilization loss of trichlorfon in 1996 (0.01%) was significantly higher than in 1997 (0.008%). This study demonstrated the fraction of applied insecticides leaving the turf putting greens was minimal.     

Field-scale variation in microbial activity and soil properties in relation to mineralization and sorption of pesticides in a sandy soil

Pesticides applied to agricultural soils are subject to environmental concerns because leaching to groundwater reservoirs and aquatic habitats may occur. Knowledge of field variation of pesticide-related parameters is required to evaluate the vulnerability of pesticide leaching. The mineralization and sorption of the pesticides glyphosate and metribuzin and the pesticide degradation product triazinamin in a field were measured and compared with the field-scale variation of geochemical and microbiological parameters. We focused on the soil parameters clay and organic carbon (C) content and on soil respiratory and enzymatic processes and microbial biomass. These parameters were measured in soil samples taken at two depths (Ap and Bs horizon) in 51 sampling points from a 4-ha agricultural fine sandy soil field. The results indicated that the spatial variation of the soil parameters, and in particular the content of organic C, had a major influence on the variability of the microbial parameters and on sorption and pesticide mineralization in the soil. For glyphosate, with a co-metabolic pathway for degradation, the mineralization was increased in soils with high microbial activity. The spatial variability, expressed as the CV, was about five times higher in the Bs horizon than in the Ap horizon, and the local-scale variation within 100 m(2) areas were two to three times lower than the field-scale variation within the entire field of about 4 ha.     

Denitrification at a long-term forested land treatment system in the Piedmont of Georgia

Spray irrigation of forested land can provide an effective system for nutrient removal and treatment of municipal wastewater. Evolution of N2 + N2O from denitrifying activity is an important renovation pathway for N applied to forested land treatment systems. Federal and state guidance documents for design of forested land treatment systems indicate the expected range for denitrification to be up to 25% of applied N, and most forest land treatment systems are designed using values from 15 to 20% of applied N. However, few measurements of denitrification following long-term wastewater applications at forested land treatment sites exist. In this study, soil N2 + N2O-N evolution was directly measured at four different landscape positions (hilltop, midslope, toe-slope, and riparian zone) in a forested land treatment facility in the Georgia Piedmont that has been operating for more than 13 yr. Denitrification rates within effluent-irrigated areas were significantly greater than rates in adjacent nonirrigated buffer zones. Rates of N2 + N2O-N evolved from soil in irrigated forests ranged from 5 to 10 kg ha(-1) yr(-1) N on the three upland landscape positions and averaged 38 kg ha(-1) yr(-1) N within the riparian zone. The relationship between measured riparian zone denitrification rates and soil physical and chemical properties was poor. The best relationship was with soil temperature, with an r2 of 0.18. Overall, on a landscape position weighted basis, only 2.4% of the wastewater-applied N was lost through denitrification.