Predicting unsaturated zone nitrogen mass balances in agricultural settings of the United States.

Unsaturated zone N fate and transport were evaluated at four sites to identify the predominant pathways of N cycling: an almond [Prunus dulcis (Mill.) D.A. Webb] orchard and cornfield (Zea mays L.) in the lower Merced River study basin, California; and corn-soybean [Glycine max (L.) Merr.] rotations in study basins at Maple Creek, Nebraska, and at Morgan Creek, Maryland. We used inverse modeling with a new version of the Root Zone Water Quality Model (RZWQM2) to estimate soil hydraulic and nitrogen transformation parameters throughout the unsaturated zone; previous versions were limited to 3-m depth and relied on manual calibration. The overall goal of the modeling was to derive unsaturated zone N mass balances for the four sites. RZWQM2 showed promise for deeper simulation profiles. Relative root mean square error (RRMSE) values for predicted and observed nitrate concentrations in lysimeters were 0.40 and 0.52 for California (6.5 m depth) and Nebraska (10 m), respectively, and index of agreement (d) values were 0.60 and 0.71 (d varies between 0 and 1, with higher values indicating better agreement). For the shallow simulation profile (1 m) in Maryland, RRMSE and d for nitrate were 0.22 and 0.86, respectively. Except for Nebraska, predictions of average nitrate concentration at the bottom of the simulation profile agreed reasonably well with measured concentrations in monitoring wells. The largest additions of N were predicted to come from inorganic fertilizer (153-195 kg N ha(-1) yr(-1) in California) and N fixation (99 and 131 kg N ha(-1) yr(-1) in Maryland and Nebraska, respectively). Predicted N losses occurred primarily through plant uptake (144-237 kg N ha(-1) yr(-1)) and deep seepage out of the profile (56-102 kg N ha(-1) yr(-1)). Large reservoirs of organic N (up to 17,500 kg N ha(-1) m(-1) at Nebraska) were predicted to reside in the unsaturated zone, which has implications for potential future transfer of nitrate to groundwater.

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.

Pesticide fate and transport throughout unsaturated zones in five agricultural settings, USA.

Pesticide transport through the unsaturated zone is a function of chemical and soil characteristics, application, and water recharge rate. The fate and transport of 82 pesticides and degradates were investigated at five different agricultural sites. Atrazine and metolachlor, as well as several of the degradates of atrazine, metolachlor, acetochlor, and alachlor, were frequently detected in soil water during the 2004 growing season, and degradates were generally more abundant than parent compounds. Metolachlor and atrazine were applied at a Nebraska site the same year as sampling, and focused recharge coupled with the short time since application resulted in their movement in the unsaturated zone 9 m below the surface. At other sites where the herbicides were applied 1 to 2 yr before sampling, only degradates were found in soil water. Transformations of herbicides were evident with depth and during the 4-mo sampling time and reflected the faster degradation of metolachlor oxanilic acid and persistence of metolachor ethanesulfonic acid. The fraction of metolachlor ethanesulfonic acid relative to metolachlor and metolachlor oxanilic acid increased from 0.3 to >0.9 at a site in Maryland where the unsaturated zone was 5 m deep and from 0.3 to 0.5 at the shallowest depth. The flux of pesticide degradates from the deepest sites to the shallow ground water was greatest (3.0-4.9 micromol m(-2) yr(-1)) where upland recharge or focused flow moved the most water through the unsaturated zone. Flux estimates based on estimated recharge rates and measured concentrations were in agreement with fluxes estimated using an unsaturated-zone computer model (LEACHM).

Variations in pesticide leaching related to land use, pesticide properties, and unsaturated zone thickness.

Pesticide leaching through variably thick soils beneath agricultural fields in Morgan Creek, Maryland was simulated for water years 1995 to 2004 using LEACHM (Leaching Estimation and Chemistry Model). Fifteen individual models were constructed to simulate five depths and three crop rotations with associated pesticide applications. Unsaturated zone thickness averaged 4.7 m but reached a maximum of 18.7 m. Average annual recharge to ground water decreased from 15.9 to 11.1 cm as the unsaturated zone increased in thickness from 1 to 10 m. These point estimates of recharge are at the lower end of previously published values, which used methods that integrate over larger areas capturing focused recharge in the numerous detention ponds in the watershed. The total amount of applied and leached masses for five parent pesticide compounds and seven metabolites were estimated for the 32-km2 Morgan Creek watershed by associating each hectare to the closest one-dimensional model analog of model depth and crop rotation scenario as determined from land-use surveys. LEACHM parameters were set such that branched, serial, first-order decay of pesticides and metabolites was realistically simulated. Leaching is predicted to be greatest for shallow soils and for persistent compounds with low sorptivity. Based on simulation results, percent parent compounds leached within the watershed can be described by a regression model of the form e(-depth) (a ln t1/2-b ln K OC) where t1/2 is the degradation half-life in aerobic soils, K OC is the organic carbon normalized sorption coefficient, and a and b are fitted coefficients (R2 = 0.86, p value = 7 x 10(-9)).

Use of 87Sr/86Sr and delta11B to identify slag-affected sediment in southern Lake Michigan.

Slag is a ubiquitous byproduct of the iron-smelting industry and influences geochemistry and water quality in adjacent geologic units, ground and surface water. Despite extensive slag deposition along the Indiana shoreline of Lake Michigan, definitive evidence that slag has affected lakebed sediments has not been established. Concerns for the protection of water and ecosystem resources in the Great Lakes motivated this study to determine if strontium and boron isotopes could be used to identify and delineate slag-affected bed sediment in Lake Michigan. Sixty-five samples of bed sediment were acquired from the southern lobe of Lake Michigan and analyzed for 87Sr/86Sr and deltat11B. Samples immediately offshore from Indiana steel mills and slag-disposal sites contained higher median 87Sr/86Sr values (0.70881) than shoreline sediments collected elsewhere in the basin (0.70847) and uniquely decreased with increasing distance from the shoreline. The highest delta11B values occurred in sediments from the Indiana shoreline (+12.9 to 16.4/1000) but were also elevated in sediments collected offshore from three Lake Michigan cities (+11.7 to 12.7/1000). Contoured isotope data indicated that 82-154 km2 of bed sediment along the Indiana shoreline had elevated 87Sr/86Sr and delta11B values relative to shoreline sediments elsewhere in southern Lake Michigan.