Effects of soil variability and weather conditions on pesticide leaching--a farm-level evaluation

In line with European regulations, Dutch law imposes an environmental threshold of 0.1 microg L(-1) on pesticide concentrations in ground water. During registration, the risk of exceeding this threshold is assessed through simulations for one or a few standard scenarios that do not reflect spatial variability under field conditions. The introduction of precision agriculture, where soil variability is actively managed, can increase control over pesticide leaching. This study presents a step-wise evaluation of the effects of soil variability and weather conditions on pesticide leaching. The evaluation was conducted on a 100-ha arable farm and aimed at identifying opportunities for precision management. As a first step, a relative risk assessment identified pesticides presenting a relatively high risk to the environment. Second, the effect of weather conditions was analyzed through 20 years of simulations for three distinct soil profiles. Results were summarized in cumulative probability plots to provide a probabilistic characterization of historical weather data. The year matching 90% probability (1981) served as a reference to simulate pesticide leaching from 612 soil profiles. After interpolation, areas where concentrations exceeded the environmental threshold were identified. Out of a total of 19 pesticides, isoproturon [N-dimethyl-N'-(4-(1-methylethyl)phenyl)urea], metribuzin [4-amino-6-tert-butyl-3-(methylthio)-as-triazin-5(4H)-one], and bentazon [2,1,3-benzothiadiazin-4(3H)-one, 3-isopropyl-, 2,2-dioxide] showed the highest risk for leaching. Leaching was strongly affected by soil variability at the within-field, field, and farm levels. Opportunities for precision management were apparent, but depended on the scale level at which environmental thresholds were implemented. When legislation is formulated in this issue, the presented step-wise evaluation can serve as a basis for identification and precision management of high-risk pesticides.     

Degradation and sorption of metribuzin and primary metabolites in a sandy soil

Leaching to the ground water of metabolites from the herbicide metribuzin [4-amino-6-(1,1-dimethylethyl)-3-(methylthio)-1,2,4-triazin-5-one] has been measured in a Danish field experiment in concentrations exceeding the European Union threshold limit for pesticides at 0.1 microg/L. In the present work, degradation and sorption of metribuzin and the metabolites desamino-metribuzin (DA), diketo-metribuzin (DK), and desamino-diketo-metribuzin (DADK) were studied in a Danish sandy loam topsoil and subsoil from the field in question, using accelerated solvent extraction and liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis. Fast dissipation of metribuzin and the metabolites was observed in the topsoil, with 50% disappearance within 30 to 40 d. A two-compartment model described degradation of metribuzin and DA, whereas that of DADK could be described using first-order kinetics. Part of the dissipation was probably due to incorporation into soil organic matter. Degradation in subsoil occurred very slowly, with extrapolated half-lives of more than one year. Sorption in the topsoil followed the order DA > metribuzin > DK > DADK. Subsoil sorption was considerably lower, and was hardly measurable for metribuzin and DK. Abiotic degradation was considerably higher in the topsoil than the subsoil, especially concerning the de-amination step, indicating that organic matter may be related to the degradation process. The present results confirm observations of metribuzin and transformation product leaching made in the field experiment and demonstrate the need for knowledge on primary metabolites when assessing the risk for pesticide leaching.     

Effects of Forest Composition and Spatial Patterns on Storm Flows of a Small Watershed1

Abstract: The PRMS_Storm model was built as a storm event, distributed hydrological model for studying the hydrological effects of forest composition and spatial distribution on storm‐flow volume and peakflow rates in the Xiangshuixi Watershed in the Three Gorges Reservoir Area, in the Yangtze River Basin in southwestern China. We developed three simulation scenarios based on forest composition and their spatial arrangements across the watershed, including all mixed conifer‐evergreen broadleaf forests (Scenario 1), all mixed evergreen broadleaf forests (Scenario 2), and mixed conifer + evergreen broadleaf + shrub forests (Scenario 3). We examined 11 storm events observed during 2002‐2005. Compared with the existing forest covers, modeling results suggested that the amount of overland flow was reduced by 21, 23, and 22%, and the interflow increased by 16, 88, and 30%, for Scenarios 1, 2, and 3, respectively. During the same time, peakflow rates were reduced by 20.8, 9.6, and 18.9%, respectively. The reduction of peakflow rates was most significant when rainfall intensity exceeded 0.8 mm/min and events with a short duration and effect was minor when rainfall intensity was below 0.5 mm/min. In general, we found that Scenarios 1 and 3 were preferred for reducing storm‐flow volume and peakflow rates due to their higher interception rates, large soil water holding capacity, and higher soil infiltration capacity. The modeled results suggested soil properties are important in affecting the flow processes and thus forest composition and forest spatial distributions will affect storm‐flow volume and peakflow rates at the watershed scale. To maximize flood reduction functions of a watershed, high priority should be given to those forest types (Scenarios 1 and 3) in reforestation practices in the study region. This study suggests both forest composition and spatial pattern are important reforestation designs for flood reduction in the Three Gorges Reservoir Area.     

SPATIALLY DISTRIBUTED MODELING OF STREAM FLOW DURING STORM EVENTS1

ABSTRACT: The purpose of this paper is to present a new approach for the spatially distributed modeling of water flow during storm events. Distributed modeling of flow during storm events is an important basis for any environmental modeling, including turbidity or sediment transport. During the initial phase of a rainstorm, surface runoff is the main contributor of flow. To provide the spatial components for distributed hydrological modeling a Geographic Information System (GIS) was used to map and visualize contributing areas around a stream channel. Stream segments were defined using the hydrologic response unit (HRU) concept. Lateral flows were derived from GIS output for each segment of the stream and at each time interval of the rain storm and were routed using the kinematic routing equation. This approach is new in hydrological modeling and can be used to enhance many existing simulations. The model is also unique in the fine time scale (i.e., intervals are on the order of minutes). Model results showed good correlation with measured discharge values; however, further studies of contributing area behavior, its relationship with soil types and slope categories, and the influence of watershed size are needed to improve model performance. This model will be used in the future as the basis to model turbidity in streams.     

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.     

Predicting phosphorus losses with the PLEASE model on a local scale in Denmark and the Netherlands

To reduce losses from agricultural soils to surface water, mitigation options have to be implemented as a local scale. For a cost-effective implementation of these measures, an instrument to identify critical areas for P leaching is indispensable. In many countries, P-index methods are used to identify areas as risk for P losses to surface water. In flat areas, where losses by leaching are dominant, these methods have their limitations because leaching is often not described in detail, PLEASE, is a simple mechanistic model designed to stimulate P Losses by leaching at the field scale using a limited amount of local field data. In this study, PLEASE, was applied to 17 lowland sites in Denmark and 14 lowland sites in the Netherlands. Results show that the simple model simulated measured fluxes and concentrations in water from pipe drains, suction cups, and groundwater quite well. The modeling efficiency ranged from 0.92 for modeling total-P fluxes to 0.36 fr modeling concentrations in groundwater. Poor results were obtained for heavy clay soils and eutrophic peat soils, where fluxes and concentration were strongly underestimated by the model. The poot performance for the heavy clay soil can be explained by the transport of P through macropores to the drain pipes and the underestimation of overland flow on this heavy-textured soil. In the eutrophic peat soils, fluxes were underestimated due to the release of P from deep soil layers.     

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.     

Retention and runoff losses of atrazine and metribuzin in soil

Minimizing herbicide runoff and mobility in the soil and thus potential contamination of water resources is a national concern. Metribuzin [4-amino-6-(1,1-dimethylethyl)-3-(methylthio)-1,2,4-triazin-5(4H)-one] and atrazine [2-chloro-4-ethylamino-6-isopropylamino-1,3,5-triazine] dynamics in surface soils and in runoff waters were studied on six 0.2-ha sugarcane (Saccharum spp.) plots of a Commerce silt loam (fine-silty, mixed, superactive, nonacid, thermic Fluvaquentic Endoaquept) during three growing seasons under different best management practices. Metribuzin was applied in the spring as a postemergence herbicide and atrazine was applied following winter harvest. Both herbicides were applied on top of the sugarcane rows as 0.6- or 0.9-m band width application, or broadcast application, where the entire area was treated. Maximum effluent concentrations were measured from the broadcast treatment and ranged from 600 to 1100 microg L(-1) for atrazine and 250 to 450 microg L(-1) for metribuzin. Atrazine runoff losses were highest for the broadcast treatment (2.8-11% of that applied) and lowest for the 0.6-m band treatment (1.9-7.6%), with a similar trend for metribuzin losses. Measured extractable herbicides from the surface soil exhibited a sharp decrease with time and were well described with a simple first-order decay model. For atrazine, estimates for the decay rate (lambda) were higher than for metribuzin. Results based on laboratory adsorption-desorption (kinetic-batch) measurements were consistent with field observations. The distribution coefficients (Kd) for atrazine exhibited stronger retention over time in comparison with metribuzin on the Commerce soil. Moreover, discrepancies between adsorption isotherm and desorption indicated slower release and that hysteresis was more pronounced for atrazine compared with metribuzin.     

Modeling macropore flow effects on pesticide leaching: inverse parameter estimation using microlysimeters

Macropore flow is a key factor determining pesticide fate, but models accounting for this process need parameters that cannot be easily measured. This study was conducted to investigate the use of inverse techniques to estimate parameters controlling macropore flow and pesticide fate in the dual-permeability model MACRO. Undisturbed columns were sampled at three landscape positions (hilltop, slope, hollow) with contrasting texture and organic carbon content. Transient leaching experiments were performed for an anionic tracer and the herbicide MCPA (4-chloro-2methylphenoxy acetic acid) during a 4-mo period, first under natural rainfall, and then under controlled irrigation in the laboratory. The tracer breakthrough for the liner-textured soil from the hilltop showed strong evidence of macropore flow, resulting in a rapid leaching of MCPA, while leaching was minimal from the organic-rich hollow soil, since macropore flow was weaker and adsorption stronger. The MACRO model was linked to the inverse modeling program SUFI (Sequential Uncertainty Fitting) to enable calibration of nine key model parameters. Based on calculated model efficiencies, MACRO-SUFI gave generally good predictions of water movement and tracer and pesticide transport, although some errors were attributed to difficulties in simulating the effects of soil moisture on degradation and the timing of water outflows. Even after calibration, significant uncertainties remained for some key parameters controlling macropore flow. Nevertheless, the parameter estimates were significantly different between landscape positions and could also be related to basic soil properties. The posterior uncertainty ranges could probably be reduced with a more exhaustive sampling of the parameter space and improved experimental designs.     

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.     

Fungicide leaching from golf greens: effects of root zone composition and surfactant use

Soil water repellency in golf putting greens may induce preferential "finger flow," leading to enhanced leaching of surface applied fungicides. We examined the effects of root zone composition, treatment with a non-ionic surfactant, and the use of the fungicide iprodion or a combination of azoxystrobin and propiconazole on soil water repellency, soil water content distributions, fungicide leaching, and turf quality during 1 yr. Soil water repellency was measured using the water drop penetration time (WDPT) test and tension infiltrometers. Our study was made on a 3-yr-old experimental green seeded with creeping bentgrass (Agrostis stolonifera L.) 'Penn A-4' at Landvik in southeast Norway. The facility consists of 16 lysimeters with two different root zone materials: (i) straight sand (1% gravel, 96% sand, 3% silt and clay, 4 g kg(-1) organic matter) (SS) and (ii) straight sand mixed with garden compost to an organic matter content of 21 g kg(-1) (Green Mix [GM]). Surfactant treatment resulted in 96% lower average WDPTs at 1 cm depth, three times higher water infiltration rates at the soil surface, and reduced spatial variation in soil water contents. Fungicide leaching was close to zero for the GM lysimeters probably due to stronger sorption. Concentrations in the drainage water from SS lysimeters often exceeded surface water guideline values for all three fungicides, but surfactant treatment dramatically reduced fungicide leaching from these lysimeters. In autumn and winter, surfactant-treated plots were more infected with fungal diseases probably because of higher water content in the turfgrass thatch layer.     

Solute movement through an allophanic soil

Allophanic soils are widespread around the world, but little research has been done on their transport properties. This study reveals the effect of two soil water potential heads and two water-flow regimes of continuous and intermittent flow on solute transport through undisturbed soil columns of Horotiu silt loam (Typic Hapludand), an allophanic soil. Two different methods--breakthrough curves (BTCs) and time domain reflectometry (TDR)--were employed to determine the extent of preferential solute transport in the topsoil. The TDR data were also used to look at the depth dependence of the transport properties. The convection-dispersion equation (CDE) with the appropriate boundary conditions adequately described the movement of both Br and Cl under the various flow conditions. Although no preferential flow was found under the imposed unsaturated flow conditions, the flow of water and transport of solute became more uniform with depth. The results show that both Br and Cl are retarded in this allophanic soil. Retardation values range from 1.5 to 1.9, and, as the TDR data showed, increase from the depth of 5.0 to 10.0 cm. Intermittent leaching results showed that there was no effect on solute concentrations in the leachate following no-flow periods. This suggests that water and solute transport in this soil were either relatively uniform or that transverse mixing during flow was already fast enough to eliminate concentration gradients between regions of different "mobility."     

Effect of manure application timing, crop, and soil type on phosphorus leaching

Phosphorus (P) leaching losses from manure applications may be of concern when artificial drainage systems allow for hydrologic short-cuts to surface waters. This study quantified P leaching losses from liquid manure applications on two soil textural extremes, a clay loam and loamy sand soil, as affected by cropping system and timing of application. For each soil type, manure was applied at an annual rate of 93 800 L ha(-1) on replicated drained plots under maize (Zea mays L.) in early fall, late fall, early spring, and as a split application in early and late spring. Manure was applied on orchardgrass (Dactylis glomerata L.) in split applications in early fall and late spring, and early and late spring. Drain water was sampled at least weekly when lines were flowing, and outflow rate and total P content were determined. High P leaching losses were measured in the clay loam as soon as drain lines initiated flow after manure application. Flow-weighted mean P leaching losses on clay loam plots averaged 39 times higher (0.504 mg L(-1)) than those on loamy sand plots (0.013 mg L(-1)), and were above the USEPA level of concern of 0.1 mg L(-1). Phosphorus losses varied among application seasons on the clay loam soil, with highest losses generally measured for early fall applications. Phosphorus leaching patterns in clay loam showed short-term spikes and high losses were associated with high drain outflow rates, suggesting preferential flow as the main transport mechanism. Phosphorus leaching from manure applications on loamy sand soils does not pose environmental concerns as long as soil P levels remain below the saturation level.     

MODELING PERCOLATION LOSSES FROM A PONDED FIELD UNDER VARIABLE WATER-TABLE CONDITIONS1

ABSTRACT: A numerical simulation model was developed to predict the vertical and lateral percolation losses from a ponded agricultural field. The two-dimensional steady-state unsaturated/ saturated flow equation was solved using the finite-difference technique. A constant ponding depth was maintained at the soil surface with different water table conditions in an application of the model for rice fields bordered by bunds. Field experiments were conducted for two different water table depths to collect data on the spatial distribution of volumetric soil-moisture content for model verification. The measured soil-moisture content values were found to be in close agreement with those predicted by the model. The sensitivity analysis of the model with selected hydrologic conditions shows that the model is most sensitive to the values of saturated hydraulic conductivity, but relatively less sensitive to water table depth, ponding depth, and evaporation rate from the soil surface. It implies that, in a ponded rice field condition, the lateral and vertical percolation losses are mostly governed by the hydraulic conductivity of the soil. The vertical percolation losses were almost equal to the saturated hydraulic conductivity values and, in most cases, these losses increased with deeper water table depths. The lateral percolation losses also increased with deeper water table depths; however, these losses were relatively small in comparison to the vertical percolation losses. The vertical and lateral percolation losses increased with the increase in ponding depths. The lateral percolation losses through the bund decreased when the evaporation losses increased from the soil surface. The results of this study indicate that the percolation losses from a ponded field may be predicted accurately for a wide range of soil and hydrological conditions when the values of hydraulic conductivity, evaporation rate, depth of ponding, and water table depth are accurately known.     

Runoff and drainage losses of atrazine, metribuzin, and metolachlor in three water management systems

Rainfall can transport herbicides from agricultural land to surface waters, where they become an environmental concern. Tile drainage can benefit crop production by removing excess soil water but tile drainage may also aggravate herbicide and nutrient movement into surface waters. Water management of tile drains after planting may reduce tile drainage and thereby reduce herbicide losses to surface water. To test this hypothesis we calculated the loss of three herbicides from a field with three water management systems: free drainage (D), controlled drainage (CD), and controlled drainage with subsurface irrigation (CDS). The effect of water management systems on the dissipation of atrazine (6-chloro-N2-ethyl-N4-isopropyl-1,3,5-triazine-2,4-diamine), metribuzin [4-amino-6-(1,1-dimethylethyl)-3-(methylthio)-1,2,4-triazine-5(4H)-one), and metolachlor [2-chloro-N-(2-ethyl-6-methylphenyl)-N-(2-methoxy-1-methylethyl)acetamide] in soil was also monitored. Less herbicide was lost by surface runoff from the D and CD treatments than from CDS. The CDS treatment increased surface runoff, which transported more herbicide than that from D or CD treatments. In one year, the time for metribuzin residue to dissipate to half its initial value was shorter for CDS (33 d) than for D (43 d) and CD (46 d). The half-life of atrazine and metolachlor were not affected by water management. Controlled drainage with subsurface irrigation may increase herbicide loss through increased surface runoff when excessive rain is received soon after herbicide application. However, increasing soil water content in CDS may decrease herbicide persistence, resulting in less residual herbicide available for aqueous transport.     

Assessing ground water vulnerability with the type transfer function model in the San Joaquin Valley, California

The recently developed type transfer function (TTF) simulation approach was applied to generate a regional-scale nonpoint-source ground water vulnerability assessment for the San Joaquin Valley, California. The computationally comparatively inexpensive TTF approach produces quantitative estimates of contaminant concentrations for large regional scales through characteristic functions based on different soil textures and their leaching properties. The TTF simulations employed an extensive soil and recharge database to estimate atrazine (1-chloro-3-ethylamino-5-isopropylamino-2,4,6-triazine) concentrations at a compliance depth of 3 m resulting from a surface application. Two different sets of TTFs with two different levels of upscaling were used for spatially uniform and distributed recharge estimates. Results show that estimated atrazine concentrations can be related to soil survey information. Areas with high potential vulnerability to atrazine leaching were found for soils with low organic carbon content and sandy loam and loam textures. Travel times for atrazine peak concentrations to the compliance depth ranged from 350 to 730 d. The extent of areas with estimated atrazine concentrations above the maximum contaminant level was less extensive when uniform annual recharge values were used. Simulated TTF concentrations were highest for eastern Fresno County, a vulnerability pattern that is also supported by field observations. The TTF modeling approach is shown to be a useful tool for quantitative pesticide leaching estimates at regional scales significantly larger than those of previous studies.     

Soil water repellency induced by long-term irrigation with treated sewage effluent

This study describes soil water repellency developed under prolonged irrigation with treated sewage effluent in a semiarid environment. Soil surface layer (0-5 cm) and soil profile (0-50 cm) transects were sampled at a high resolution at the close of the irrigation season and rainy winter season. Samples from 0- to 5-cm transects were subdivided into 1-cm slices to obtain fine scale resolution of repellency and organic matter distribution. Extreme to severe soil water repellency in the 0- to 5-cm soil surface layer persisted throughout the 2-yr study period in the effluent-irrigated Shamouti orange [Citrus sinensis (L.) Osbeck cv. Shamouti] orchard plot. Nearby Shamouti orange plots irrigated with tap water were either nonrepellent or only somewhat repellent. Repellency was very variable spatially and with depth, appearing in vertically oriented "repellency tongues." Temporal and spatial variability in repellency in the uppermost 5-cm soil surface layer was not related to seasonality, soil moisture content, or soil organic matter content. Nonuniform distribution of soil moisture and fingered flow were observed in the soil profile after both seasons, demonstrating that the repellent layer had a persistent effect on water flow in the soil profile. A lack of correlation between bulk density and volumetric water content in the soil profile demonstrates that the observed nonuniform spatial distribution of moisture results from preferential flow and not heterogeneity in soil properties. Soil water repellency can adversely affect agricultural production, cause contamination of underlying ground water resources, and result in excessive runoff and soil erosion.     

LANDSCAPE ATTRIBUTES AS CONTROLS ON GROITHD WATER NITRATE REMOVAL CAPACITY OF RIPARIAN ZONES1

ABSTRACT: Inherent site factors can generate substantial variation in the ground water nitrate removal capacity of riparian zones. This paper examines research in the glaciated Northeast to relate variability in ground water nitrate removal to site attributes depicted in readily available spatial databases, such as SSUIRGO. Linking site‐specific studies of riparian ground water nitrate removal to spatial data can help target high‐value riparian locations for restoration or protection and improve the modeling of watershed nitrogen flux. Site attributes, such as hydric soil status (soil wetness) and geomorphology, affect the interaction of nitrate‐enriched ground water with portions of the soil ecosystem possessing elevated biogeochemical transformation rates (i.e., biologically active zones). At our riparian sites, high ground water nitrate‐N removal rates were restricted to hydric soils. Geomorphology provided insights into ground water flowpaths. Riparian sites located on outwash and organic/alluvial deposits have high potential for nitrate‐enriched ground water to interact with biologically active zones. In till deposits, ground water nitrate removal capacity may be limited by the high occurrence of surface seeps that markedly reduce the time available for biological transformations to occur within the riparian zone. To fully realize the value of riparian zones for nitrate retention, landscape controls of riparian nitrate removal in different climatic and physiographic regions must be determined and translated into available spatial databases.     

Spatial Calibration and Temporal Validation of Flow for Regional Scale Hydrologic Modeling1

Abstract: Physically based regional scale hydrologic modeling is gaining importance for planning and management of water resources. Calibration and validation of such regional scale model is necessary before applying it for scenario assessment. However, in most regional scale hydrologic modeling, flow validation is performed at the river basin outlet without accounting for spatial variations in hydrological parameters within the subunits. In this study, we calibrated the model to capture the spatial variations in runoff at subwatershed level to assure local water balance, and validated the streamflow at key gaging stations along the river to assure temporal variability. Ohio and Arkansas‐White‐Red River Basins of the United States were modeled using Soil and Water Assessment Tool (SWAT) for the period from 1961 to 1990. R² values of average annual runoff at subwatersheds were 0.78 and 0.99 for the Ohio and Arkansas Basins. Observed and simulated annual and monthly streamflow from 1961 to 1990 is used for temporal validation at the gages. R² values estimated were greater than 0.6. In summary, spatially distributed calibration at subwatersheds and temporal validation at the stream gages accounted for the spatial and temporal hydrological patterns reasonably well in the two river basins. This study highlights the importance of spatially distributed calibration and validation in large river basins.     

REGIONAL NITRATE LEACHING VARIABILITY: WHAT MAKES A DIFFERENCE IN NORTHEASTERN COLORADO

ABSTRACT: The high spatial variability of nitrate concentrations in ground water of many regions is thought to be closely related to spatially-variable leaching rates from agricultural activities. To clarify the relative roles of the different nitrate leaching controlling variables under irrigated agriculture in northeastern Colorado, we conducted an extensive series of leaching simulations with the NLEAP model using best estimates of local agricultural practices. The results of these simulations were then used with GIS to estimate the spatial variability of leachate quality for a 14,000 ha area overlying the alluvial aquifer of the South Platte River. Simulations showed that in the study area, differences in soil type might lead to 5–10 kg/ha of N variation in annual leaching rates while variability due to crop rotations was as much as 65 kg-N/ha for common rotations. Land application of manure from confined animal feeding operations may account for more than 100 kg-N/ha additional leaching. For a selected index rotation, the simulated nitrogen leaching rates across the area varied from 10 to 299 kg/ha and simulated water volumes leached ranged from 13 to 76 cm/yr depending on soil type, irrigation type, and use of manure. Resulting leachate concentrations of 3.5–140 mg/l NO₃ as N were simulated. Land application of manure was found to be the most important factor determining the mass flux of nitrate leached and the combination of sprinkler irrigation and manure application yields the highest leachate concentrations.