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,045 records · Page 58Linked to original sources

Hillslope hydrologic connectivity controls riparian groundwater turnover: Implications of catchment structure for riparian buffering and stream water sources

Hydrologic connectivity between catchment upland and near stream areas is essential for the transmission of water, solutes, and nutrients to streams. However, our current understanding of the role of riparian zones in mediating landscape hydrologic connectivity and the catchment scale export of water and solutes is limited. We tested the relationship between the duration of hillslope‐riparian‐stream (HRS) hydrologic connectivity and the rate and degree of riparian shallow groundwater turnover along four HRS well transects within a set of nested mountain catchments (Tenderfoot Creek Experimental Forest, MT). Transect HRS water table connectivity ranged from 9 to 123 days during the annual snowmelt hydrograph. Hillslope water was always characterized by low specific conductance (∼27 μ S cm −1 ). In transects with transient hillslope water tables, riparian groundwater specific conductance was elevated during base flow conditions (∼127 μ S cm −1 ) but shifted toward hillslope signatures once a HRS groundwater connection was established. The degree of riparian groundwater turnover was proportional to the duration of HRS connectivity and inversely related to the riparian: hillslope area ratios (buffer ratio; r 2 = 0.95). We applied this relationship to the stream network in seven subcatchments within the Tenderfoot Creek Experimental Forest and compared their turnover distributions to source water contributions measured at each catchment outlet. The amount of riparian groundwater exiting each of the seven catchments was linearly related ( r 2 = 0.92) to their median riparian turnover time. Our observations suggest that the size and spatial arrangement of hillslope and riparian zones along a stream network and the timing and duration of groundwater connectivity between them is a first‐order control on the magnitude and timing of water and solutes observed at the catchment outlet.

Water Resources Research↗

Metamodeling for groundwater age forecasting in the Lake Michigan Basin

Groundwater age is an important indicator of groundwater susceptibility to anthropogenic contamination and a key input to statistical models for forecasting water quality. Numerical models can provide estimates of groundwater age, enabling interpretation of measured age tracers. However, to extend to national‐scale groundwater systems where numerical models are not routinely available, a more efficient metamodeling approach can provide a less precise but widely applicable estimate of groundwater age, trained to make forecasts based on predictor variables that can be measured independent of numerical models. We trained gradient‐boosted regression tree statistical metamodels to MODFLOW/MODPATH‐derived groundwater age estimates in five inset models in the Lake Michigan Basin, USA. Using high‐throughput computing, we explored an exhaustive range of tuning parameters and tested metamodels through cross validation, a 20% holdout, and a round robin approach among the five inset models withholding each inset model from training and testing on the held‐out inset model. Forecast skill—measured by Nash Sutcliffe efficiency—was high for age‐related responses in the 20% hold‐out case (ranging from 0.73 to 0.84). The round robin analysis provided the opportunity to explore extending to unmodeled areas and a greater range of skill indicated the need to evaluate when it is appropriate to apply a metamodel from one region to another. We further explored the ramifications of metamodel simplification achieved through removing predictor variables based on their estimated importance. We found that similar metamodel performance was achievable with a fraction of the candidate set of predictor variables with well construction variables being most important.

Water Resources Research↗

Probabilistic categorical groundwater salinity mapping from airborne electromagnetic data adjacent to California’s Lost Hills and Belridge oil fields

Growing water stress has led to emerging interest in protecting fresh and brackish groundwater as a potential supplement to water supplies and raised questions about factors that could affect the future quality of fresh and brackish aquifers. Limited well infrastructure, particularly in regions where elevated salinity has led to limited historical groundwater development, hinders traditional mapping of salinity distributions through groundwater sampling. This paper presents a quantitative salinity mapping approach of the upper 300 m using high‐resolution, regionally comprehensive resistivity models derived from Bayesian inversion of an airborne electromagnetic survey adjacent to the Lost Hills and Belridge oil fields in the southwestern San Joaquin Valley of California. Using local water quality observations as an interpretational foundation, a probabilistic approach yields maps of fresh, saline, and brackish groundwater while quantifying joint uncertainty inherited from the geophysical data and interpretational relations. Saline and fresh regions are mapped with relatively high confidence in many locations, while areas of lower confidence, particularly at depth, can be mapped as their most probable salinity category while reflecting the relative uncertainty in the interpretation. These maps identify a stratified salinity structure, where saline water commonly occurs in the surficial aquifer overlying fresher groundwater in the Tulare aquifer, separated by regional confining clay layers. Downgradient of unlined surface water diversions, recharge of imported surface water results in relatively fresh groundwater throughout the depth of investigation.

California↗

Calibrated simulation of the long-term average surficial groundwater system and derived spatial distributions of its characteristics for the contiguous United States

While the physical processes governing groundwater flow are well understood, and the computational resources now exist for solving the governing equations in three dimensions over continental-scale domains, there remains substantial uncertainty about the subsurface distribution of the properties that control groundwater flow and transport for much of the contiguous United States (CONUS). The transmissivity of the shallow subsurface is a key parameter for the simulation of water table position, shallow groundwater flow, and base-flow discharge, but is not well-characterized at large regional to continental scales. We used a process-based inversion of CONUS-extent groundwater information to generate national data sets of (a) the transmissivity of the shallow groundwater system, (b) the depth to the water table, (c) groundwater discharge as base-flow, and (d) long-term average water content in the unsaturated zone. CONUS-extent coverage was developed in the form of 75 subdomain models, with the spatial distribution of long-term average transmissivity for each subdomain model calibrated against water-levels derived from U.S. Geological Survey (USGS) observation wells, NHDPlusV2 first-order perennial streams, and National Wetlands Inventory (NWI) freshwater wetlands. Estimated transmissivities were lower in the western CONUS than the eastern CONUS, and across the CONUS both transmissivity and depth to water correlate with recharge, elevation, and topographic slope. These generated data sets provide spatially distributed, long-term average estimates of subsurface properties and hydrological states that we anticipate will complement other environmental modeling efforts as explanatory variables, boundary conditions, or transport pathways.

contiguous United States↗

Improved prediction of management-relevant groundwater discharge characteristics throughout river networks

Groundwater discharge zones connect aquifers to surface water, generating baseflow and serving as ecosystem control points across aquatic ecosystems. The influence of groundwater discharge on surface flow connectivity, fate and transport of contaminants and nutrients, and thermal habitat depends strongly on hydrologic characteristics such as the spatial distribution, age, and depth of source groundwater flow paths. Groundwater models have the potential to predict spatial discharge characteristics within river networks, but models are often not evaluated against these critical characteristics and model equifinality with respect to discharge processes is a known challenge. We quantify discharge characteristics across a suite of groundwater models with commonly used frameworks and calibration data. We developed a base model (MODFLOW‐NWT) for a 1,570‐km 2 watershed in the northeastern United States and varied the calibration data, control of river‐aquifer exchange directionality, and resolution. Most models ( n = 11 of 12) fit similarly to calibration metrics, but patterns in discharge location, flow path depth, and subsurface travel time varied substantially. We found (1) a 15% difference in the percent of discharge going to first‐order streams, (2) threefold variations in flow path depth, and (3) sevenfold variations in the subsurface travel times among the models. We recalibrated three models using a synthetic discharge location data set. Calibration with discharge location data reduced differences in simulated discharge characteristics, suggesting an approach to improved equifinality based on widespread field‐based mapping of discharge zones. Our work quantifying variation across common modeling approaches is an important step toward characterizing and improving predictions of groundwater discharge characteristics.

Connecticut, Massachusetts↗

Estimating and forecasting time-varying groundwater recharge in fractured rock: A state-space formulation with preferential and diffuse flow to the water table

Rapid infiltration following precipitation may result in groundwater contamination from surface contaminants or pathogens. In fractured rock, contaminants can migrate rapidly to points of groundwater withdrawals. In contrast to the temporal availability of groundwater quality chemical indicators, meteorological and groundwater level observations are available in real-time to estimate time-varying recharge, which can act as a surrogate to identify periods of rapid infiltration that may indicate contamination susceptibility. Estimating recharge using methods, such as base-flow recession, unsaturated infiltration models, or Water-Table Fluctuations (WTF), cannot capitalize on currently available technologies and telecommunication infrastructure to conduct real-time recharge estimation at scales relevant to characterizing rapid infiltration. We present a linear, physics-based State-Space (SS) model of one-dimensional infiltration to estimate recharge, which includes preferential and diffuse-flow to the water table. The model can take advantage of real-time data for water-table altitude, precipitation, and evapotranspiration. Model parameters are calibrated over an observation period, and the Kalman Filter (KF) is subsequently applied to continuously update the observed (water-table altitude) and unobserved (groundwater recharge) system states and predict future states as new data become available. The SS/KF algorithm is demonstrated at the Masser Groundwater Recharge Site in Pennsylvania, USA and comparisons are made with recharge estimates from WTF methods. Model results indicate that the frequency of observations (daily versus sub-daily) dictates the allocation between preferential and diffuse flow. Additionally, because infiltration processes encompass many nonlinearities, model parameters estimated from observation periods need to be updated at least seasonally to account for changing recharge conditions.

Water Resources Research↗

WellSTIC: A cost-effective sensor for performing point dilution tests to measure groundwater velocity in shallow aquifers

Many individual measurement points are required to characterize groundwater velocity within an aquifer. Groundwater velocity is most commonly measured using a network of >5 cm diameter monitoring wells, which, if not already present at a site, are expensive and labor-intensive to install. Drive-point piezometers—simple, cost-effective wells that can be installed by hand—are a common tool for sampling groundwater in shallow, alluvial aquifers, but most groundwater velocity measurement techniques require equipment that is too large for these narrow (usually <2 cm inside diameter) piezometers. In this technical note, we introduce a low-cost sensor and well packer system (<$90 USD) for performing point dilution tests in narrow piezometers. Field data show that the magnitude of groundwater velocity measured with this technique agrees with velocities computed from natural gradient tracer tests. Additionally, with proper calibration, these sensors can be used to continuously monitor in-well specific conductance, either during inter-well tracer tests with saline tracers or for water quality monitoring. This system is a viable tool for rapid assessment of the magnitude of groundwater velocity in shallow aquifers.

Water Resources Research↗

Seasonal and species-level water-use strategies and groundwater dependence in dryland riparian woodlands during extreme drought

Drought-induced groundwater decline and warming associated with climate change are primary threats to dryland riparian woodlands. We used the extreme 2012–2019 drought in southern California as a natural experiment to assess how differences in water-use strategies and groundwater dependence may influence the drought susceptibility of dryland riparian tree species with overlapping distributions. We analyzed tree-ring stable carbon and oxygen isotopes collected from two cottonwood species ( Populus trichocarpa and P . fremontii ) along the semi-arid Santa Clara River. We also modeled tree source water δ 18 O composition to compare with observed source water δ 18 O within the floodplain to infer patterns of groundwater reliance. Our results suggest that both species functioned as facultative phreatophytes that used shallow soil moisture when available but ultimately relied on groundwater to maintain physiological function during drought. We also observed apparent species differences in water-use strategies and groundwater dependence related to their regional distributions. P . fremontii was constrained to more arid river segments and ostensibly used a greater proportion of groundwater to satisfy higher evaporative demand. P . fremontii maintained ∆ 13 C at pre-drought levels up until the peak of the drought, when trees experienced a precipitous decline in ∆ 13 C. This response pattern suggests that trees prioritized maintaining photosynthetic processes over hydraulic safety, until a critical point. In contrast, P . trichocarpa showed a more gradual and sustained reduction in ∆ 13 C, indicating that drought conditions induced stomatal closure and higher water use efficiency. This strategy may confer drought avoidance for P . trichocarpa while increasing its susceptibility to anticipated climate warming.

California↗

Quantifying groundwater response and uncertainty in beaver-influenced mountainous floodplains using machine learning-based model calibration

Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model-data integration workflow to address uncertainty in groundwater response to beaver-induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning-based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down-valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down-valley underflow, from 2% during dry pond periods to 20% during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down-valley underflow, minimizes the effect of beaver-induced inundation on water quality. We emphasize the need for field-scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response.

Colorado↗

Groundwater sustainability strategies

Groundwater extraction has facilitated significant social development and economic growth, enhanced food security and alleviated drought in many farming regions. But groundwater development has also depressed water tables, degraded ecosystems and led to the deterioration of groundwater quality, as well as to conflict among water users. The effects are not evenly spread. In some areas of India, for example, groundwater depletion has preferentially affected the poor. Importantly, groundwater in some aquifers is renewed slowly, over decades to millennia, and coupled climate–aquifer models predict that the flux and/or timing of recharge to many aquifers will change under future climate scenarios. Here we argue that communities need to set multigenerational goals if groundwater is to be managed sustainably.

Nature Geoscience↗

Groundwater discharge impacts marine isotope budgets of Li, Mg, Ca, Sr, and Ba

Groundwater-derived solute fluxes to the ocean have long been assumed static and subordinate to riverine fluxes, if not neglected entirely, in marine isotope budgets. Here we present concentration and isotope data for Li, Mg, Ca, Sr, and Ba in coastal groundwaters to constrain the importance of groundwater discharge in mediating the magnitude and isotopic composition of terrestrially derived solute fluxes to the ocean. Data were extrapolated globally using three independent volumetric estimates of groundwater discharge to coastal waters, from which we estimate that groundwater-derived solute fluxes represent, at a minimum, 5% of riverine fluxes for Li, Mg, Ca, Sr, and Ba. The isotopic compositions of the groundwater-derived Mg, Ca, and Sr fluxes are distinct from global riverine averages, while Li and Ba fluxes are isotopically indistinguishable from rivers. These differences reflect a strong dependence on coastal lithology that should be considered a priority for parameterization in Earth-system models.

Nature Communications↗

Arsenic and other geogenic contaminants in global groundwater

Geogenic groundwater contaminants (GGCs) affect drinking-water availability and safety, with up to 60% of groundwater sources in some regions contaminated by more than recommended concentrations. As a result, an estimated 300–500 million people are at risk of severe health impacts and premature mortality. In this Review, we discuss the sources, occurrences and cycling of arsenic, fluoride, selenium and uranium, which are GGCs with widespread distribution and/or high toxicity. The global distribution of GGCs is controlled by basin geology and tectonics, with GGC enrichment in both orogenic systems and cratonic basement rocks. This regional distribution is broadly influenced by climate, geomorphology and hydrogeochemical evolution along groundwater flow paths. GGC distribution is locally heterogeneous and affected by in situ lithology, groundwater flow and water–rock interactions. Local biogeochemical cycling also determines GGC concentrations, as arsenic, selenium and uranium mobilizations are strongly redox-dependent. Increasing groundwater extraction and land-use changes are likely to modify GGC distribution and extent, potentially exacerbating human exposure to GGCs, but the net impact of these activities is unknown. Integration of science, policy, community involvement programmes and technological interventions is needed to manage GGC-enriched groundwater and ensure equitable access to clean water.

Nature Reviews Earth & Environment↗

Over a third of groundwater in USA public-supply aquifers is Anthropocene-age and susceptible to surface contamination

The distribution of groundwater age is useful for evaluating the susceptibility and sustainability of groundwater resources. Here, we compute the aquifer-scale cumulative distribution function to characterize the age distribution for 21 Principal Aquifers that account for ~80% of public-supply pumping in the United States. The aquifer-scale cumulative distribution function for each Principal Aquifer was derived from an ensemble of modeled age distributions (~60 samples per aquifer) based on multiple tracers: tritium, tritiogenic helium-3, sulfur hexafluoride, chlorofluorocarbons, carbon-14, and radiogenic helium-4. Nationally, the groundwater is 38% Anthropocene (since 1953), 34% Holocene (75 – 11,800 years ago), and 28% Pleistocene (>11,800 years ago). The Anthropocene fraction ranges from <5 to 100%, indicating a wide range in susceptibility to land-surface contamination. The Pleistocene fraction of groundwater exceeds 50% in 7 eastern aquifers that are predominately confined. The Holocene fraction of groundwater exceeds 50% in 5 western aquifers that are predominately unconfined. The sustainability of pumping from these Principal Aquifers depends on rates of recharge and release of groundwater stored in fine-grained layers.

Nature Communications Earth & Environment↗

Effects of stormwater infiltration on quality of groundwater beneath retention and detention basins

Infiltration of storm water through detention and retention basins may increase the risk of groundwater contamination, especially in areas where the soil is sandy and the water table shallow, and contaminants may not have a chance to degrade or sorb onto soil particles before reaching the saturated zone. Groundwater from 16 monitoring wells installed in basins in southern New Jersey was compared to the quality of shallow groundwater from 30 wells in areas of new-urban land use. Basin groundwater contained much lower levels of dissolved oxygen, which affected concentrations of major ions. Patterns of volatile organic compound and pesticide occurrence in basin groundwater reflected the land use in the drainage areas served by the basins, and differed from patterns in background samples, exhibiting a greater occurrence of petroleum hydrocarbons and certain pesticides. Dilution effects and volatilization likely decrease the concentration and detection frequency of certain compounds commonly found in background groundwater. High recharge rates in storm water basins may cause loading factors to be substantial even when constituent concentrations in infiltrating storm water are relatively low.

New Jersey↗

Evaluation of the sustainability of deep groundwater as an arsenic-safe resource in the Bengal Basin

Tens of millions of people in the Bengal Basin region of Bangladesh and India drink groundwater containing unsafe concentrations of arsenic. This high-arsenic groundwater is produced from shallow (<100 m) depths by domestic and irrigation wells in the Bengal Basin aquifer system. The government of Bangladesh has begun to install wells to depths of >150 m where groundwater arsenic concentrations are nearly uniformly low, and many more wells are needed, however, the sustainability of deep, arsenic-safe groundwater has not been previously assessed. Deeper pumping could induce downward migration of dissolved arsenic, permanently destroying the deep resource. Here, it is shown, through quantitative, large-scale hydrogeologic analysis and simulation of the entire basin, that the deeper part of the aquifer system may provide a sustainable source of arsenic-safe water if its utilization is limited to domestic supply. Simulations provide two explanations for this result: deep domestic pumping only slightly perturbs the deep groundwater flow system, and substantial shallow pumping for irrigation forms a hydraulic barrier that protects deeper resources from shallow arsenic sources. Additional analysis indicates that this simple management approach could provide arsenic-safe drinking water to >90% of the arsenic-impacted region over a 1,000-year timescale. This insight may assist water-resources managers in alleviating one of the world's largest groundwater contamination problems.

Proceedings of the National Academy of Sciences of↗

Evaluating impacts of subdivision density on shallow groundwater in Southeastern Wisconsin, USA

Using simple numerical groundwater flow models, we tested the impacts of suburban developments on groundwater levels and discharge to streams. We used lot sizes of 1, 3 and 5 acres (4000, 12,000 and 20,000 m 2 ) with one domestic well per lot that pumped water from shallow aquifers. Our modelling showed that pumping had little impact on water levels and groundwater discharge to streams if the developed area is of a moderate size. However, domestic wells had the potential to impact local groundwater levels and baseflows in large developments. In township-wide development scenarios of 1-acre (4000 m 2 ) lots, simulated drawdowns beneath developed areas ranged from 1 to 18 ft (0.3 to 5.5 m), and baseflow reductions ranged from 20 to 40%. Impacts generally were inversely proportional to lot size, recharge rate and hydraulic conductivity of the aquifer materials. Developments using individual domestic wells have the potential to impact local groundwater levels and surface water features. The impacts can range from negligible to severe, depending on local hydrogeologic conditions and on whether wastewater is recharged onsite or is removed from the basin. An assessment of groundwater impacts should be a part of the planning process for all suburban developments.

Wisconsin↗

Understanding heat and groundwater flow through continental flood basalt provinces: insights gained from alternative models of permeability/depth relationships for the Columbia Plateau, USA

Heat-flow mapping of the western USA has identified an apparent low-heat-flow anomaly coincident with the Columbia Plateau Regional Aquifer System, a thick sequence of basalt aquifers within the Columbia River Basalt Group (CRBG). A heat and mass transport model (SUTRA) was used to evaluate the potential impact of groundwater flow on heat flow along two different regional groundwater flow paths. Limited in situ permeability (k) data from the CRBG are compatible with a steep permeability decrease (approximately 3.5 orders of magnitude) at 600&ndash;900 m depth and approximately 40&deg;C. Numerical simulations incorporating this permeability decrease demonstrate that regional groundwater flow can explain lower-than-expected heat flow in these highly anisotropic (k x /k z ~ 10 4 ) continental flood basalts. Simulation results indicate that the abrupt reduction in permeability at approximately 600 m depth results in an equivalently abrupt transition from a shallow region where heat flow is affected by groundwater flow to a deeper region of conduction-dominated heat flow. Most existing heat-flow measurements within the CRBG are from shallower than 600 m depth or near regional groundwater discharge zones, so that heat-flow maps generated using these data are likely influenced by groundwater flow. Substantial k decreases at similar temperatures have also been observed in the volcanic rocks of the adjacent Cascade Range volcanic arc and at Kilauea Volcano, Hawaii, where they result from low-temperature hydrothermal alteration.

Idaho;Oregon;Washington↗

Evaluating micrometeorological estimates of groundwater discharge from Great Basin desert playas

Groundwater availability studies in the arid southwestern United States traditionally have assumed that groundwater discharge by evapotranspiration (ETg) from desert playas is a significant component of the groundwater budget. However, desert playa ETg rates are poorly constrained by Bowen Ratio energy budget (BREB) and eddy-covariance (EC) micrometeorological measurement approaches. Best attempts by previous studies to constrain ETg from desert playas have resulted in ETg rates that are within the measurement error of micrometeorological approaches. This study uses numerical models to further constrain desert playa ETg rates that are within the measurement error of BREB and EC approaches, and to evaluate the effect of hydraulic properties and salinity-based groundwater-density contrasts on desert playa ETg rates. Numerical models simulated ETg rates from desert playas in Death Valley, California and Dixie Valley, Nevada. Results indicate that actual ETg rates from desert playas are significantly below the uncertainty thresholds of BREB- and EC-based micrometeorological measurements. Discharge from desert playas likely contributes less than 2 percent of total groundwater discharge from Dixie and Death Valleys, which suggests discharge from desert playas also is negligible in other basins. Simulation results also show that ETg from desert playas primarily is limited by differences in hydraulic properties between alluvial fan and playa sediments and, to a lesser extent, by salinity-based groundwater density contrasts.

Ground Water↗