USGS Science⌕ Search

SEARCH · USGS Science

Results for “Groundwater”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64Linked to original sources

Riparian groundwater and baseflow studies in the Upper Colorado River Basin

Executive summary As part of an ongoing effort to understand baseflow in the Upper Colorado River Basin (UCRB) and implications for stream-dependent ecosystems, we conducted a brief review of literature related to groundwater and baseflow in the UCRB. We included primary literature, federal and state resources, databases and gray literature studies on groundwater, baseflow, and springs in the UCRB region. This review provides examples of the types of groundwater and baseflow studies published for the UCRB with sections on whole-basin studies, catchment and reach studies and their locations, water quality studies, studies adjacent to the UCRB, state and federal groundwater resources, and finally a discussion of potential further directions. Despite the limited nature of the review, we summarize numerous studies related to groundwater in the UCRB which will be valuable to researchers interested in groundwater and baseflow dynamics in the region.

Arizona, Colorado, New Mexico, Utah, Wyoming↗

Regional land subsidence caused by the compaction of susceptible aquifer systems accompanying groundwater extraction

Land subsidence includes both gentle downwarping and sudden sinking of segments of the land surface. Major anthropogenic causes of land subsidence are extraction of fluids including water, oil, and gas. Measurement and detec- tion of land subsidence include both ground-based and remotely sensed air- borne and space-based methods. Methods for measurement of subsidence at points include differential leveling, global positioning system surveys, and extensometers. Satellite-borne differential interferometric synthetic aperture radar and airborne LiDAR techniques can detect land-surface movement over wide areas of interest. Aquifer-system compaction and subsidence owing to groundwater extraction typically occurs in areas of unconsolidated alluvial or basin-fill aquifer systems comprising aquifers and aquitards. Approaches to analyzing and modeling deformation of aquifer systems follow from the basic relations between head, stress, compressibility, and groundwater flow. Analysis and simulation of aquifer-system compaction have been addressed primarily using either an approach based on conventional groundwater flow theory or an approach based on linear poroelasticity theory. Both approaches rely on the principle of effective stress outlined by Karl Terzaghi in 1925. In the approach based on conventional groundwater flow theory, an aquitard drainage model explains the compaction of fine grained material using the principle of effective stress and theory of hydrodynamic lag. Packages for the widely-used MODFLOW groundwater model are available to simulate aqui- fer-system compaction and land subsidence using the aquitard-drainage approach. Poroelasticity theory describes the more fully coupled processes of groundwater flow and three-dimensional deformation of aquifer systems. The general theory accounts for compressible fluid, porous matrix and solid grains. Simulation codes using the poroelastic theory include some commer- cial software products and a few research codes.

Book chapter↗

Depletion mapping and constrained optimization to support managing groundwater extraction

Groundwater models often serve as management tools to evaluate competing water uses including ecosystems, irrigated agriculture, industry, municipal supply, and others. Depletion potential mapping—showing the model-calculated potential impacts that wells have on stream baseflow—can form the basis for multiple potential management approaches in an oversubscribed basin. Specific management approaches can include scenarios proposed by stakeholders, systematic changes in well pumping based on depletion potential, and formal constrained optimization, which can be used to quantify the tradeoff between water use and stream baseflow. Variables such as the maximum amount of reduction allowed in each well and various groupings of wells using, for example, K-means clustering considering spatial proximity and depletion potential are considered. These approaches provide a potential starting point and guidance for resource managers and stakeholders to make decisions about groundwater management in a basin, spreading responsibility in different ways. We illustrate these approaches in the Little Plover River basin in central Wisconsin, United States—home to a rich agricultural tradition, with farmland and urban areas both in close proximity to a groundwater-dependent trout stream. Groundwater withdrawals have reduced baseflow supplying the Little Plover River below a legally established minimum. The techniques in this work were developed in response to engaged stakeholders with various interests and goals for the basin. They sought to develop a collaborative management plan at a watershed scale that restores the flow rate in the river in a manner that incorporates principles of shared governance and results in effective and minimally disruptive changes in groundwater extraction practices.

Groundwater↗

A Monte-Carlo chemical budget approach to assess ambient groundwater flow in bedrock open boreholes

In low-permeability rocks, ambient groundwater flow in open boreholes may go undetected using conventional borehole-flowmeter tools and alternative approaches may be needed to identify flow. Understanding ambient flow in open boreholes is important for tracking of cross contamination in groundwater. Chlorinated volatile organic compound (CVOC) concentrations from three open boreholes set in a crystalline-rock aquifer (two of three open boreholes) and a siltstone aquifer (one of three open boreholes) were examined using a new approach and associated software program called the AFCE (Aqueous-Flow-Concentration-Estimator). The program allows comparison of coupled chemical datasets through a Monte-Carlo simulation and a chemical-budget approach to assess ambient groundwater flow in open boreholes. The coupled datasets required for the comparison include aqueous CVOC concentrations from groundwater samples from (1) discrete fractures, such as those measured from temporary deployment of straddle-borehole packer assemblies; and (2) the concentration of the open borehole (wellbore) water, as measured by a vertical profile of passive samplers from within the same open borehole. Because results from the passive samplers represent a composite mixture of the results from the discrete samples under ambient groundwater-flow conditions, potentially at unknown proportions, the comparison between coupled datasets affords the ability to discern likely water contributions of CVOC from discrete fractures (or fracture zones), and which fractures may be dominating the water chemistry of the open borehole.

Groundwater Monitoring and Remediation↗

Evaluating geothermal and hydrogeologic controls on regional groundwater temperature distribution

A one-dimensional (1-D) analytic solution is developed for heat transport through an aquifer system where the vertical temperature profile in the aquifer is nearly uniform. The general anisotropic form of the viscous heat generation term is developed for use in groundwater flow simulations. The 1-D solution is extended to more complex geometries by solving the equation for piece-wise linear or uniform properties and boundary conditions. A moderately complex example, the Eastern Snake River Plain (ESRP), is analyzed to demonstrate the use of the analytic solution for identifying important physical processes. For example, it is shown that viscous heating is variably important and that heat conduction to the land surface is a primary control on the distribution of aquifer and spring temperatures. Use of published values for all aquifer and thermal properties results in a reasonable match between simulated and measured groundwater temperatures over most of the 300 km length of the ESRP, except for geothermal heat flow into the base of the aquifer within 20 km of the Yellowstone hotspot. Previous basal heat flow measurements (∼110 mW/m 2 ) made beneath the ESRP aquifer were collected at distances of >50 km from the Yellowstone Plateau, but a higher basal heat flow of 150 mW/m 2 is required to match groundwater temperatures near the Plateau. The ESRP example demonstrates how the new tool can be used during preliminary analysis of a groundwater system, allowing efficient identification of the important physical processes that must be represented during more-complex 2-D and 3-D simulations of combined groundwater and heat flow.

Idaho↗

Combined use of thermal methods and seepage meters to efficiently locate, quantify, and monitor focused groundwater discharge to a sand-bed stream

Quantifying flow of groundwater through streambeds often is difficult due to the complexity of aquifer-scale heterogeneity combined with local-scale hyporheic exchange. We used fiber-optic distributed temperature sensing (FO-DTS), seepage meters, and vertical temperature profiling to locate, quantify, and monitor areas of focused groundwater discharge in a geomorphically simple sand-bed stream. This combined approach allowed us to rapidly focus efforts at locations where prodigious amounts of groundwater discharged to the Quashnet River on Cape Cod, Massachusetts, northeastern USA. FO-DTS detected numerous anomalously cold reaches one to several m long that persisted over two summers. Seepage meters positioned upstream, within, and downstream of 7 anomalously cold reaches indicated that rapid groundwater discharge occurred precisely where the bed was cold; median upward seepage was nearly 5 times faster than seepage measured in streambed areas not identified as cold. Vertical temperature profilers deployed next to 8 seepage meters provided diurnal-signal-based seepage estimates that compared remarkably well with seepage-meter values. Regression slope and R 2 values both were near 1 for seepage ranging from 0.05 to 3.0 m d −1 . Temperature-based seepage model accuracy was improved with thermal diffusivity determined locally from diurnal signals. Similar calculations provided values for streambed sediment scour and deposition at subdaily resolution. Seepage was strongly heterogeneous even along a sand-bed river that flows over a relatively uniform sand and fine-gravel aquifer. FO-DTS was an efficient method for detecting areas of rapid groundwater discharge, even in a strongly gaining river, that can then be quantified over time with inexpensive streambed thermal methods.

Massachusetts↗

Generation of 3-D hydrostratigraphic zones from dense airborne electromagnetic data to assess groundwater model prediction error

We present a new methodology to combine spatially dense high-resolution airborne electromagnetic (AEM) data and sparse borehole information to construct multiple plausible geological structures using a stochastic approach. The method developed allows for quantification of the performance of groundwater models built from different geological realizations of structure. Multiple structural realizations are generated using geostatistical Monte Carlo simulations that treat sparse borehole lithological observations as hard data and dense geophysically derived structural probabilities as soft data. Each structural model is used to define 3-D hydrostratigraphical zones of a groundwater model, and the hydraulic parameter values of the zones are estimated by using nonlinear regression to fit hydrological data (hydraulic head and river discharge measurements). Use of the methodology is demonstrated for a synthetic domain having structures of categorical deposits consisting of sand, silt, or clay. It is shown that using dense AEM data with the methodology can significantly improve the estimated accuracy of the sediment distribution as compared to when borehole data are used alone. It is also shown that this use of AEM data can improve the predictive capability of a calibrated groundwater model that uses the geological structures as zones. However, such structural models will always contain errors because even with dense AEM data it is not possible to perfectly resolve the structures of a groundwater system. It is shown that when using such erroneous structures in a groundwater model, they can lead to biased parameter estimates and biased model predictions, therefore impairing the model's predictive capability.

Water Resources Research↗

Predicting redox-sensitive contaminant concentrations in groundwater using random forest classification

Machine learning techniques were applied to a large (n > 10,000) compliance monitoring database to predict the occurrence of several redox-active constituents in groundwater across a large watershed. Specifically, random forest classification was used to determine the probabilities of detecting elevated concentrations of nitrate, iron, and arsenic in the Fox, Wolf, Peshtigo, and surrounding watersheds in northeastern Wisconsin. Random forest classification is well suited to describe the nonlinear relationships observed among several explanatory variables and the predicted probabilities of elevated concentrations of nitrate, iron, and arsenic. Maps of the probability of elevated nitrate, iron, and arsenic can be used to assess groundwater vulnerability and the vulnerability of streams to contaminants derived from groundwater. Processes responsible for elevated concentrations are elucidated using partial dependence plots. For example, an increase in the probability of elevated iron and arsenic occurred when well depths coincided with the glacial/bedrock interface, suggesting a bedrock source for these constituents. Furthermore, groundwater in contact with Ordovician bedrock has a higher likelihood of elevated iron concentrations, which supports the hypothesis that groundwater liberates iron from a sulfide-bearing secondary cement horizon of Ordovician age. Application of machine learning techniques to existing compliance monitoring data offers an opportunity to broadly assess aquifer and stream vulnerability at regional and national scales and to better understand geochemical processes responsible for observed conditions.

Water Resources Research↗

Regional variability of nitrate fluxes in the unsaturated zone and groundwater, Wisconsin, USA

Process-based modeling of regional NO3− fluxes to groundwater is critical for understanding and managing water quality, but the complexity of NO3− reactive transport processes make implementation a challenge. This study introduces a regional vertical flux method (VFM) for efficient estimation of reactive transport of NO3− in the vadose zone and groundwater. The regional VFM was applied to 443 well samples in central-eastern Wisconsin. Chemical measurements included O2, NO3−, N2 from denitrification, and atmospheric tracers of groundwater age including carbon-14, chlorofluorocarbons, tritium, and tritiogenic helium. VFM results were consistent with observed chemistry, and calibrated parameters were in-line with estimates from previous studies. Results indicated that (1) unsaturated zone travel times were a substantial portion of the transit time to wells and streams (2) since 1945 fractions of applied N leached to groundwater have increased for manure-N, possibly due to increased injection of liquid manure, and decreased for fertilizer-N, and (3) under current practices and conditions, approximately 60% of the shallow aquifer will eventually be affected by downward migration of NO3−, with denitrification protecting the remaining 40%. Recharge variability strongly affected the unsaturated zone lag times and the eventual depth of the NO3− front. Principal components regression demonstrated that VFM parameters and predictions were significantly correlated with hydrogeochemical landscape features. The diverse and sometimes conflicting aspects of N management (e.g. limiting N volatilization versus limiting N losses to groundwater) warrant continued development of large-scale holistic strategies to manage water quality and quantity.

Wisconsin↗

The Hawai'i groundwater recharge tool

This article discusses the design and implementation of the Hawai’i Groundwater Recharge Tool, an application for providing data and analyses of the impacts of land-cover modifications and changes in precipitation on groundwater-recharge rates for the island of O’ahu. This application uses simulation data based on a set of 29 land-cover types and 2 precipitation conditions to provide users with real-time recharge calculations for interactively defined land-cover modifications. The tool provides two visualizations, representing the land cover for the island and the resultant groundwater-recharge rates, and a set of metrics indicating the changes to groundwater recharge for relevant areas to present a set of easily interpretable outcomes based on user-defined scenarios. Users have varying degrees of control over the granularity of data input and output, allowing for the quick production of a roughly defined scenario, or more precise land-cover definitions. These modifications can be exported for further analysis. Heuristics are used to provide a responsive user interface and performant integration with the database containing the full set of simulation data. This tool is designed to provide user-friendly access to the information on the impacts of land-cover and precipitation changes on groundwater-recharge rates needed to assist in making data-driven decisions.

Hawai'i↗

Characterization of surface and ground water δ18O seasonal variation and its use for estimating groundwater residence times

18 O is an ideal tracer for characterizing hydrological processes because it can be reliably measured in several watershed hydrological compartments. Here, we present multiyear isotopic data, i.e. 18 O variations (δ 18 O), for precipitation inputs, surface water and groundwater in the Shingobee River Headwaters Area (SRHA), a well-instrumented research catchment in north-central Minnesota. SRHA surface waters exhibit δ 18 O seasonal variations similar to those of groundwaters, and seasonal δ 18 O variations plotted versus time fit seasonal sine functions. These seasonal δ 18 O variations were interpreted to estimate surface water and groundwater mean residence times (MRTs) at sampling locations near topographically closed-basin lakes. MRT variations of about 1 to 16 years have been estimated over an area covering about 9 km 2 from the basin boundary to the most downgradient well. Estimated MRT error (±0·3 to ±0·7 years) is small for short MRTs and is much larger (±10 years) for a well with an MRT (16 years) near the limit of the method. Groundwater transit time estimates based on Darcy's law, tritium content, and the seasonal δ 18 O amplitude approach appear to be consistent within the limits of each method. The results from this study suggest that use of the δ 18 O seasonal variation method to determine MRTs can help assess groundwater recharge areas in small headwaters catchments.

Minnesota↗

Global synthesis of groundwater recharge in semiarid and arid regions

Global synthesis of the findings from ∼140 recharge study areas in semiarid and arid regions provides important information on recharge rates, controls, and processes, which are critical for sustainable water development. Water resource evaluation, dryland salinity assessment (Australia), and radioactive waste disposal (US) are among the primary goals of many of these recharge studies. The chloride mass balance (CMB) technique is widely used to estimate recharge. Average recharge rates estimated over large areas (40–374 000 km 2 ) range from 0·2 to 35 mm year −1 , representing 0·1–5% of long-term average annual precipitation. Extreme local variability in recharge, with rates up to ∼720 m year −1 , results from focussed recharge beneath ephemeral streams and lakes and preferential flow mostly in fractured systems. System response to climate variability and land use/land cover (LU/LC) changes is archived in unsaturated zone tracer profiles and in groundwater level fluctuations. Inter-annual climate variability related to El Niño Southern Oscillation (ENSO) results in up to three times higher recharge in regions within the SW US during periods of frequent El Niños (1977–1998) relative to periods dominated by La Niñas (1941–1957). Enhanced recharge related to ENSO is also documented in Argentina. Climate variability at decadal to century scales recorded in chloride profiles in Africa results in recharge rates of 30 mm year −1 during the Sahel drought (1970–1986) to 150 mm year −1 during non-drought periods. Variations in climate at millennial scales in the SW US changed systems from recharge during the Pleistocene glacial period (≥10 000 years ago) to discharge during the Holocene semiarid period. LU/LC changes such as deforestation in Australia increased recharge up to about 2 orders of magnitude. Changes from natural grassland and shrublands to dryland (rain-fed) agriculture altered systems from discharge (evapotranspiration, ET) to recharge in the SW US. The impact of LU change was much greater than climate variability in Niger (Africa), where replacement of savanna by crops increased recharge by about an order of magnitude even during severe droughts. Sensitivity of recharge to LU/LC changes suggests that recharge may be controlled through management of LU. In irrigated areas, recharge varies from 10 to 485 mm year −1 , representing 1–25% of irrigation plus precipitation. However, irrigation pumpage in groundwater-fed irrigated areas greatly exceeds recharge rates, resulting in groundwater mining. Increased recharge related to cultivation has mobilized salts that accumulated in the unsaturated zone over millennia, resulting in widespread groundwater and surface water contamination, particularly in Australia. The synthesis of recharge rates provided in this study contains valuable information for developing sustainable groundwater resource programmes within the context of climate variability and LU/LC change.

Hydrological Processes↗

Age and growth of Freshwater Drum and Gizzard Shad occupying two reservoir-river complexes with different groundwater contributions

Restoring groundwater flow is a management option that improves water temperature regimes and benefits fishes. Although this strategy applies more readily to river systems, the thermal character of reservoirs is heavily influenced by inflowing rivers. We examined differences in age, structure, and growth of both Freshwater Drum Aplodinotus grunniens and Gizzard Shad Dorosoma cepedianum that occupy catchments with varying groundwater contributions in the south‐central United States. Seepage run data indicated that the Kiamichi River was losing surface water to groundwater in summer 2016, whereas groundwater inflows were apparent in the Elk River basin. Summer 2016 data showed that the Elk River had cooler water temperatures than the Kiamichi River and Grand Lake O’ the Cherokees water temperatures were similar to those in the incoming Elk River. We found higher densities of older Freshwater Drum and Gizzard Shad (maximums of 32 and 8 years old, respectively) in samples from the Grand basin than among fish that were sampled from the Kiamichi River basin (21 and 6 years old, respectively). Freshwater Drum grew at similar rates in both basins even though they reached larger maximum lengths in the Grand basin (649 mm TL) than in the Kiamichi River basin (600 mm). The average asymptotic length was greater for the Kiamichi population ( L ∞ = 613 mm) than for the Grand population ( L ∞ = 557 mm). Gizzard Shad from the Grand basin were larger than those from the Kiamichi River basin, though the latter population grew faster initially (Brody growth coefficient: K = 0.787 versus K = 0.179, respectively), but they had smaller asymptotic length ( L ∞ = 206 mm versus L ∞ = 343 mm). The role that groundwater plays in temperature regulation in these basins partially explains the observed differences. Our results suggest that the metabolic theory of ecology can be applied to fisheries management at a finer spatial scale.

Oklahoma↗

Bridging groundwater models and decision support with a Bayesian network

Resource managers need to make decisions to plan for future environmental conditions, particularly sea level rise, in the face of substantial uncertainty. Many interacting processes factor in to the decisions they face. Advances in process models and the quantification of uncertainty have made models a valuable tool for this purpose. Long-simulation runtimes and, often, numerical instability make linking process models impractical in many cases. A method for emulating the important connections between model input and forecasts, while propagating uncertainty, has the potential to provide a bridge between complicated numerical process models and the efficiency and stability needed for decision making. We explore this using a Bayesian network (BN) to emulate a groundwater flow model. We expand on previous approaches to validating a BN by calculating forecasting skill using cross validation of a groundwater model of Assateague Island in Virginia and Maryland, USA. This BN emulation was shown to capture the important groundwater-flow characteristics and uncertainty of the groundwater system because of its connection to island morphology and sea level. Forecast power metrics associated with the validation of multiple alternative BN designs guided the selection of an optimal level of BN complexity. Assateague island is an ideal test case for exploring a forecasting tool based on current conditions because the unique hydrogeomorphological variability of the island includes a range of settings indicative of past, current, and future conditions. The resulting BN is a valuable tool for exploring the response of groundwater conditions to sea level rise in decision support.

Water Resources Research↗

Ecohydrology and Its Relation to Integrated Groundwater Management

In the twentieth century, groundwater characterization focused primarily on easily measured hydraulic metrics of water storage and flows. Twenty-first century concepts of groundwater availability, however, encompass other factors having societal value, such as ecological well-being. Effective ecohydrological science is a nexus of fundamental understanding derived from two scientific disciplines: (1) ecology, where scale, thresholds, feedbacks and tipping points for societal questions form the basis for the ecologic characterization, and (2) hydrology, where the characteristics, magnitude, and timing of water flows are characterized for a defined system of interest. In addition to ecohydrology itself, integrated groundwater management requires input from resource managers to understand which areas of the vast world of ecohydrology are important for decision making. Expectations of acceptable uncertainty, or even what ecohydrological outputs have utility, are often not well articulated within societal decision making frameworks, or within the science community itself. Similarly, “acceptable levels of impact” are difficult to define. Three examples are given to demonstrate the use of ecohydrological considerations for long-term sustainability of groundwater resources and their related ecosystem function. Such examples illustrate the importance of accommodating ecohydrogeological aspects into integrated groundwater management of the twenty-first century, regardless of society, climate, or setting.

Book chapter↗

A GIS-based groundwater travel time model to evaluate stream nitrate concentration reductions from land use change

Excessive nitrate-nitrogen (nitrate) loss from agricultural watersheds is an environmental concern. A common conservation practice to improve stream water quality is to retire vulnerable row croplands to grass. In this paper, a groundwater travel time model based on a geographic information system (GIS) analysis of readily available soil and topographic variables was used to evaluate the time needed to observe stream nitrate concentration reductions from conversion of row crop land to native prairie in Walnut Creek watershed, Iowa. Average linear groundwater velocity in 5-m cells was estimated by overlaying GIS layers of soil permeability, land slope (surrogates for hydraulic conductivity and gradient, respectively) and porosity. Cells were summed backwards from the stream network to watershed divide to develop a travel time distribution map. Results suggested that groundwater from half of the land planted in prairie has reached the stream network during the 10 years of ongoing water quality monitoring. The mean travel time for the watershed was estimated to be 10.1 years, consistent with results from a simple analytical model. The proportion of land in the watershed and subbasins with prairie groundwater reaching the stream (10-22%) was similar to the measured reduction of stream nitrate (11-36%). Results provide encouragement that additional nitrate reductions in Walnut Creek are probable in the future as reduced nitrate groundwater from distal locations discharges to the stream network in the coming years. The high spatial resolution of the model (5-m cells) and its simplicity may make it potentially applicable for land managers interested in communicating lag time issues to the public, particularly related to nitrate concentration reductions over time. ?? 2007 Springer-Verlag.

Environmental Geology↗

Use of 234U and 238U isotopes to evaluate contamination of near-surface groundwater with uranium-mill effluent: A case study in south-central Colorado, U.S.A.

The 234U/238U alpha activity ratio (AR) was determined in 47 samples of variably uraniferous groundwater from the vicinity of a uranium mill near Canon City, Colorado. The results illustrate that uranium isotopes can be used to determine the distribution of uranium contamination in groundwater and to indicate processes such as mixing and chemical precipitation that affect uranium concentrations. Highly to moderately contaminated groundwater samples collected from the mill site and land immediately downgradient from the mill site contain more than 100 ??g/l of dissolved uranium and typically have AR values in the narrow range of 1.0-1.06. Other samples from the shallow alluvial aquifer farther downgradient from the mill contain 10-100 ??g/1 uranium and plot along a broad trend of increasing AR (1.06-1.46) with decreasing uranium concentration. The results are consistent with mixing of liquid mill waste (AR ??? 1.0) with alluvial groundwater of small, but variable, uranium concentrations and AR of 1.31.5. In the alluvial aquifer, the spatial distribution of wells with AR values less than 1.3 is consistent with previous estimates of the probable distribution of contamination, based on water chemistry and hydrology. Wells more distant from the area of probable contamination have AR values that are consistently greater than 1.3 and are indicative of little or no contamination. The methodology of this study can be extended usefully to similar sites of uranium mining, milling, or processing provided that local geohydrologic settings promote uranium mobility and that introduced uranium contamination is isotopically distinct from that of local groundwater.

Environmental Geology↗

Groundwater recharge and agricultural contamination

Agriculture has had direct and indirect effects on the rates and compositions of groundwater recharge and aquifer biogeochemistry. Direct effects include dissolution and transport of excess quantities of fertilizers and associated materials and hydrologic alterations related to irrigation and drainage. Some indirect effects include changes in water–rock reactions in soils and aquifers caused by increased concentrations of dissolved oxidants, protons, and major ions. Agricultural activities have directly or indirectly affected the concentrations of a large number of inorganic chemicals in groundwater, for example NO 3 – , N 2 , Cl, SO 4 2– , H + , P, C, K, Mg, Ca, Sr, Ba, Ra, and As, as well as a wide variety of pesticides and other organic compounds. For reactive contaminants like NO 3 – , a combination of chemical, isotopic, and environmental-tracer analytical approaches might be required to resolve changing inputs from subsequent alterations as causes of concentration gradients in groundwater. Groundwater records derived from multi-component hydrostratigraphic data can be used to quantify recharge rates and residence times of water and dissolved contaminants, document past variations in recharging contaminant loads, and identify natural contaminant-remediation processes. These data indicate that many of the world's surficial aquifers contain transient records of changing agricultural contamination from the last half of the 20th century. The transient agricultural groundwater signal has important implications for long-term trends and spatial heterogeneity in discharge.

Hydrogeology Journal↗