Search USGS⌕ Search

Geology topics

Christopher T. Green

Publications and source records attributed to Christopher T. Green.

29 records · Page 2Linked to original sources

Factors controlling nitrate fluxes in groundwater in agricultural areas

The impact of agricultural chemicals on groundwater quality depends on the interactions of biogeochemical and hydrologic factors. To identify key processes affecting distribution of agricultural nitrate in groundwater, a parsimonious transport model was applied at 14 sites across the U.S. Simulated vertical profiles of NO 3 - , N 2 from denitrification, O 2 , Cl - , and environmental tracers of groundwater age were matched to observations by adjusting the parameters for recharge rate, unsaturated zone travel time, fractions of N and Cl - inputs leached to groundwater, O 2 reduction rate, O 2 threshold for denitrification, and denitrification rate. Model results revealed important interactions among biogeochemical and physical factors. Chloride fluxes decreased between the land surface and water table possibly because of Cl - exports in harvested crops (averaging 22% of land-surface Cl - inputs). Modeled zero-order rates of O 2 reduction and denitrification were correlated. Denitrification rates at depth commonly exceeded overlying O 2 reduction rates, likely because shallow geologic sources of reactive electron donors had been depleted. Projections indicated continued downward migration of NO 3 - fronts at sites with denitrification rates <0.25 mg-N L -1 yr -1 . The steady state depth of NO 3 - depended to a similar degree on application rate, leaching fraction, recharge, and NO 3 - and O 2 reaction rates. Steady state total mass in each aquifer depended primarily on the N application rate. In addition to managing application rates at land surface, efficient water use may reduce the depth and mass of N in groundwater because lower recharge was associated with lower N fraction leached. Management actions to reduce N leaching could be targeted over aquifers with high-recharge and low-denitrification rates.

Water Resources Research↗

The fate and transport of nitrate in shallow groundwater in northwestern Mississippi, USA

Agricultural contamination of groundwater in northwestern Mississippi, USA, has not been studied extensively, and subsurface fluxes of agricultural chemicals have been presumed minimal. To determine the factors controlling transport of nitrate-N into the Mississippi River Valley alluvial aquifer, a study was conducted from 2006 to 2008 to estimate fluxes of water and solutes for a site in the Bogue Phalia basin (1,250 km 2 ). Water-quality data were collected from a shallow water-table well, a vertical profile of temporary sampling points, and a nearby irrigation well. Nitrate was detected within 4.4 m of the water table but was absent in deeper waters with evidence of reducing conditions and denitrification. Recharge estimates from 6.2 to 10.9 cm/year were quantified using water-table fluctuations, a Cl – tracer method, and atmospheric age-tracers. A mathematical advection-reaction model predicted similar recharge to the aquifer, and also predicted that 15% of applied nitrogen is leached into the saturated zone. With current denitrification and application rates, the nitrate-N front is expected to remain in shallow groundwater, less than 6–9 m deep. Increasing application rates resulting from intensifying agricultural demands may advance the nitrate-N front to 16–23 m, within the zone of groundwater pumping.

Mississippi↗

Relations of hydrogeologic factors, groundwater reduction-oxidation conditions, and temporal and spatial distributions of nitrate, Central-Eastside San Joaquin Valley, California, USA

In a 2,700-km 2 area in the eastern San Joaquin Valley, California (USA), data from multiple sources were used to determine interrelations among hydrogeologic factors, reduction-oxidation (redox) conditions , and temporal and spatial distributions of nitrate (NO 3 ), a widely detected groundwater contaminant. Groundwater is predominantly modern, or mixtures of modern water, with detectable NO 3 and oxic redox conditions , but some zones have anoxic or mixed redox conditions . Anoxic conditions were associated with long residence times that occurred near the valley trough and in areas of historical groundwater discharge with shallow depth to water. Anoxic conditions also were associated with interactions of shallow, modern groundwater with soils. NO 3 concentrations were significantly lower in anoxic than oxic or mixed redox groundwater , primarily because residence times of anoxic waters exceed the duration of increased pumping and fertilizer use associated with modern agriculture. Effects of redox reactions on NO 3 concentrations were relatively minor. Dissolved N 2 gas data indicated that denitrification has eliminated gt;5 mg/L NO 3 -N in about 10% of 39 wells. Increasing NO 3 concentrations over time were slightly less prevalent in anoxic than oxic or mixed redox groundwater . Spatial and temporal trends of NO 3 are primarily controlled by water and NO 3 fluxes of modern land use.

California↗

Mixing effects on apparent reaction rates and isotope fractionation during denitrification in a heterogeneous aquifer

Gradients in contaminant concentrations and isotopic compositions commonly are used to derive reaction parameters for natural attenuation in aquifers. Differences between field‐scale (apparent) estimated reaction rates and isotopic fractionations and local‐scale (intrinsic) effects are poorly understood for complex natural systems. For a heterogeneous alluvial fan aquifer, numerical models and field observations were used to study the effects of physical heterogeneity on reaction parameter estimates. Field measurements included major ions, age tracers, stable isotopes, and dissolved gases. Parameters were estimated for the O 2 reduction rate, denitrification rate, O 2 threshold for denitrification, and stable N isotope fractionation during denitrification. For multiple geostatistical realizations of the aquifer, inverse modeling was used to establish reactive transport simulations that were consistent with field observations and served as a basis for numerical experiments to compare sample‐based estimates of “apparent” parameters with “true“ (intrinsic) values. For this aquifer, non‐Gaussian dispersion reduced the magnitudes of apparent reaction rates and isotope fractionations to a greater extent than Gaussian mixing alone. Apparent and true rate constants and fractionation parameters can differ by an order of magnitude or more, especially for samples subject to slow transport, long travel times, or rapid reactions. The effect of mixing on apparent N isotope fractionation potentially explains differences between previous laboratory and field estimates. Similarly, predicted effects on apparent O 2 threshold values for denitrification are consistent with previous reports of higher values in aquifers than in the laboratory. These results show that hydrogeological complexity substantially influences the interpretation and prediction of reactive transport.

Water Resources Research↗

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.

Journal of Environmental Quality↗

Sustainability of natural attenuation of nitrate in agricultural aquifers

Increased concentrations of nitrate in groundwater in agricultural areas, coinciding with increased use of chemical and organic fertilizers, have raised concern because of risks to environmental and human health. At some sites, these problems are mitigated by natural attenuation of nitrate as a result of microbially mediated reactions. Results from U.S. Geological Survey (USGS) research under the National Water-Quality Assessment (NAWQA) program show that reactions of dissolved nitrate with solid aquifer minerals and organic carbon help lower nitrate concentrations in groundwater beneath agricultural fields. However, increased fluxes of nitrate cause ongoing depletion of the finite pool of solid reactants. Consumption of the solid reactants diminishes the capacity of the aquifer to remove nitrate, calling into question the long-term sustainability of these natural attenuation processes.

Fact Sheet↗

Unintended consequences of biofuels production?The effects of large-scale crop conversion on water quality and quantity

In the search for renewable fuel alternatives, biofuels have gained strong political momentum. In the last decade, extensive mandates, policies, and subsidies have been adopted to foster the development of a biofuels industry in the United States. The Biofuels Initiative in the Mississippi Delta resulted in a 47-percent decrease in cotton acreage with a concurrent 288-percent increase in corn acreage in 2007. Because corn uses 80 percent more water for irrigation than cotton, and more nitrogen fertilizer is recommended for corn cultivation than for cotton, this widespread shift in crop type has implications for water quantity and water quality in the Delta. Increased water use for corn is accelerating water-level declines in the Mississippi River Valley alluvial aquifer at a time when conservation is being encouraged because of concerns about sustainability of the groundwater resource. Results from a mathematical model calibrated to existing conditions in the Delta indicate that increased fertilizer application on corn also likely will increase the extent of nitrate-nitrogen movement into the alluvial aquifer. Preliminary estimates based on surface-water modeling results indicate that higher application rates of nitrogen increase the nitrogen exported from the Yazoo River Basin to the Mississippi River by about 7 percent. Thus, the shift from cotton to corn may further contribute to hypoxic (low dissolved oxygen) conditions in the Gulf of Mexico.

Open-File Report↗

Spatial and temporal trends in nitrate concentrations in the eastern San Joaquin Valley regional aquifer and implications for fertilizer management

Ground-water withdrawals in the San Joaquin Valley totaled 64 million m 3 /day (19 million ac-ft) in 2000, supplying about 45% of agricultural irrigation demand and about 80% of municipal supply (Hutson et al., 2004). Most of the population and ground-water use are in the eastern San Joaquin Valley, where reliance on ground water is expected to increase as a result of rapid population growth and limited surface water supplies. Protection of ground-water quality for future use requires monitoring and understanding the mechanisms controlling the long-term quality of ground water in the regional aquifer system. Nitrate has been widely detected above background concentrations in ground water in the eastern San Joaquin Valley. Nitrate concentrations (reported as nitrogen in this paper) were above the MCL of 10 mg/L in 24% of domestic wells screened in the shallow part of the aquifer that were sampled during 1993–95 (Dubrovsky et al., 1998) and the Central Valley is one of the top three regions in the state in terms of the number of public drinking-water wells exceeding the MCL for nitrate (California State Water Resources Control Board, 2002). To assess spatial and temporal trends in nitrate concentrations in the eastern San Joaquin Valley and to evaluate the long-term effects of nitrogen fertilizer use on ground-water quality in this region, data were evaluated at multiple spatial scales. Data from regional-scale monitoring networks were used to map the regional occurrence of nitrate and to determine whether shallow ground water containing elevated nitrate is migrating to deeper parts of the aquifer system. At the local scale, mean ground-water ages from analysis of age-dating tracers were combined with concentrations of nitrate to reconstruct nitrate inputs in recharge through time and to compare with estimated nitrogen applications. Ground-water flow and transport simulations of a typical public-supply well screened from about 100 to 400 ft below the water table were used to evaluate long-term concentrations beneath agricultural areas under different nitrogen management scenarios.

California↗

Simulation of Multiscale Ground-Water Flow in Part of the Northeastern San Joaquin Valley, California

The transport and fate of agricultural chemicals in a variety of environmental settings is being evaluated as part of the U.S. Geological Survey (USGS) National Water-Quality Assessment Program. One of the locations being evaluated is a 2,700-km2 (square kilometer) regional study area in the northeastern San Joaquin Valley surrounding the city of Modesto, an area dominated by irrigated agriculture in a semi-arid climate. Ground water is a key source of water for irrigation and public supply, and exploitation of this resource has altered the natural flow system. The aquifer system is predominantly alluvial, and an unconfined to semiconfined aquifer overlies a confined aquifer in the southwestern part of the study area; these aquifers are separated by the lacustrine Corcoran Clay. A regional-scale 16-layer steady-state model of ground-water flow in the aquifer system in the regional study area was developed to provide boundary conditions for an embedded 110-layer steady-state local-scale model of part of the aquifer system overlying the Corcoran Clay along the Merced River. The purpose of the local-scale model was to develop a better understanding of the aquifer system and to provide a basis for simulation of reactive transport of agricultural chemicals. The heterogeneity of aquifer materials was explicitly incorporated into the regional and local models using information from geologic and drillers? logs of boreholes. Aquifer materials were differentiated in the regional model by the percentage of coarse-grained sediments in a cell, and in the local model by four hydrofacies (sand, silty sand, silt, and clay). The calibrated horizontal hydraulic conductivity values of the coarse-grained materials in the zone above the Corcoran Clay in the regional model and of the sand hydrofacies used in the local model were about equal (30?80 m/d [meter per day]), and the vertical hydraulic conductivity values in the same zone of the regional model (median of 0.012 m/d), which is dominated by the finer-grained materials, were about an order of magnitude less than that for the clay hydrofacies in the local model. Data used for calibrating both models included long-term hourly water-level measurements in 20 short-screened wells installed by the USGS in the Modesto and Merced River areas. Additional calibration data for the regional model included water-level measurements in 11 wells upslope and 17 wells downslope from these areas. The root mean square error was 2.3 m (meter) for all wells in the regional model and 0.8 m for only the USGS wells; the associated average errors were 0.9 m and 0.3 m, respectively. The root mean square error for the 12 USGS wells along a transect in the local model area was 0.08 m; the average error was 0.0 m. Particle tracking was used with the local model to estimate the concentration of an environmental tracer, sulfur hexafluoride, in 10 USGS transect wells near the Merced River that were sampled for this constituent. Measured and estimated concentrations in the mid-depth and deepest wells, which would be most sensitive to errors in hydraulic conductivity estimates, were consistent. The combined results of particle tracking and sulfur hexafluoride analysis suggest that most water sampled from the transect wells was recharged less that 25 years ago.

Scientific Investigations Report↗

Percolation and transport in a sandy soil under a natural hydraulic gradient

Unsaturated flow and transport under a natural hydraulic gradient in a Mediterranean climate were investigated with a field tracer experiment combined with laboratory analyses and numerical modeling. Bromide was applied to the surface of a sandy soil during the dry season. During the subsequent rainy season, repeated sediment sampling tracked the movement of bromide through the profile. Analysis of data on moisture content, matric pressure, unsaturated hydraulic conductivity, bulk density, and soil texture and structure provides insights into parameterization and use of the advective‐dispersive modeling approach. Capturing the gross features of tracer and moisture movement with model simulations required an order‐of‐magnitude increase in laboratory‐measured hydraulic conductivity. Wetting curve characteristics better represented field results, calling into question the routine estimation of hydraulic characteristics based only on drying conditions. Measured increases in profile moisture exceeded cumulative precipitation in early winter, indicating that gains from dew drip can exceed losses from evapotranspiration during periods of heavy (“Tule”) fog. A single‐continuum advective‐dispersive modeling approach could not reproduce a peak of bromide that was retained near the soil surface for over 3 years. Modeling of this feature required slow exchange of solute at a transfer rate of 0.5–1 × 10 −4 d −1 with an immobile volume approaching the residual moisture content.

Water Resources Research↗

Evaluation of unsaturated-zone solute-transport models for studies of agricultural chemicals

Seven unsaturated-zone solute-transport models were tested with two data sets to select models for use by the Agricultural Chemical Team of the U.S. Geological Survey's National Water-Quality Assessment Program. The data sets were from a bromide tracer test near Merced, California, and an atrazine study in the White River Basin, Indiana. In this study the models are designated either as complex or simple based on the water flux algorithm. The complex models, HYDRUS2D, LEACHP, RZWQM, and VS2DT, use Richards' equation to simulate water flux and are well suited to process understanding. The simple models, CALF, GLEAMS, and PRZM, use a tipping-bucket algorithm and are more amenable to extrapolation because they require fewer input parameters. The purpose of this report is not to endorse a particular model, but to describe useful features, potential capabilities, and possible limitations that emerged from working with the model input data sets. More rigorous assessment of model applicability involves proper calibration, which was beyond the scope of this study. Uncalibrated ("cold") simulations were run using all seven models to predict the transport of bromide (Merced) and the transport and fate of atrazine and three of its transformation products (White River Basin). Among the complex models, HYDRUS2D successfully predicted both the surface retention and accumulation of bromide at depth at the Merced site, whereas RZWQM and VS2DT predicted only the latter. RZWQM predictions of atrazine were closest to observed values at the White River Basin site, where preferential flow has been observed. LEACHP predicted smaller solute concentrations than observed at both the Merced and White River Basin sites. Among the simple models, CALF predicted the highest values of atrazine and deethylatrazine at the measurement depth of 1.5 meters. CALF includes the Addiscott flow option for preferential flow, and also accepts user-specified dispersivity. PRZM underpredicted solute concentrations, probably because control of dispersion is a problem with this model. GLEAMS has a maximum simulation depth of 1.5 meters, which is limiting for mass-balance purposes because it creates a potential disconnect between unsaturated-zone transport and the water table. Of the models tested, RZWQM, HYDRUS2D, VS2DT, GLEAMS and PRZM had graphical user interfaces. Extensive documentation was available for RZWQM, HYDRUS2D, and VS2DT. RZWQM can explicitly simulate water and solute flux in macropores, and both HYDRUS2D and VS2DT can simulate water and solute flux in two dimensions. The version of RZWQM tested had a maximum simulation depth of 3 meters. The complex models simulate the formation, transport, and fate of degradates of up to three to five compounds including the parent, with the exception of VS2DT, which simulates the transport and fate of a single compound.

California, Indiana↗