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.

SourceResearch connectionFrom yearThrough year

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

At least 1,459 records · Page 81Linked to original sources

USGS2005

Estimation of groundwater consumption by phreatophytes using diurnal water table fluctuations: A saturated‐unsaturated flow assessment

Groundwater consumption by phreatophytes is a difficult‐to‐measure but important component of the water budget in many arid and semiarid environments. Over the past 70 years the consumptive use of groundwater by phreatophytes has been estimated using a method that analyzes diurnal trends in hydrographs from wells that are screened across the water table (White, 1932). The reliability of estimates obtained with this approach has never been rigorously evaluated using saturated‐unsaturated flow simulation. We present such an evaluation for common flow geometries and a range of hydraulic properties. Results indicate that the major source of error in the White method is the uncertainty in the estimate of specific yield. Evapotranspirative consumption of groundwater will often be significantly overpredicted with the White method if the effects of drainage time and the depth to the water table on specific yield are ignored. We utilize the concept of readily available specific yield as the basis for estimation of the specific yield value appropriate for use with the White method. Guidelines are defined for estimating readily available specific yield based on sediment texture. Use of these guidelines with the White method should enable the evapotranspirative consumption of groundwater to be more accurately quantified.

Water Resources Research↗

USGS2006

Combined time-series resistivity and geochemical tracer techniques to examine submarine groundwater discharge at Dor Beach, Israel

A high-resolution, stationary geophysical and geochemical survey was conducted at Dor Beach, Israel, to examine the shallow coastal hydrogeology and its control on the exchange of submarine groundwater with the shallow Mediterranean Sea. Time-series resistivity profiles using a new 56 electrode (112-m long) marine cable produced detailed profiles of the fresh water/salt water interface and the subtle response of this interface to tidal excursions and other forcing factors. Such information, when ground-truthed with representative pore water salinities and formation resistivity values, can provide unique information of the extent and rates of submarine groundwater discharge (SGD). Time-series 222Rn measurements of the adjacent coastal water column complemented these geophysical techniques and were modeled to yield integrated advective flow rates across the sediment/water interface, which at Dor Beach ranged from about 0 to 30 cm day-1 (mean = 7.1 cm d-1), depending on the tidal range. Such results suggest that the underlying hydrogeologic framework at Dor is favorable for substantial SGD. Extrapolating these SGD estimates across a 100-m wide coastal zone implies that the Rn-derived SGD rate would equal ???7.1 m3 d-1 per m of shoreline, and that the source of this discharging groundwater is a complex mixture of fresh groundwater derived from the upland Kurkar deposits, as well as locally recycled seawater. Copyright 2006 by the American Geophysical Union.

Geophysical Research Letters↗

USGS2007

Groundwater noble gas, age, and temperature signatures in an Alpine watershed: Valuable tools in conceptual model development

Bedrock groundwater in alpine watersheds is poorly understood, mainly because of a scarcity of wells in alpine settings. Groundwater noble gas, age, and temperature data were collected from springs and wells with depths of 3–342 m in Handcart Gulch, an alpine watershed in Colorado. Temperature profiles indicate active groundwater circulation to a maximum depth (aquifer thickness) of about 200 m, or about 150 m below the water table. Dissolved noble gas data show unusually high excess air concentrations (>0.02 cm 3 STP/g, ΔNe > 170%) in the bedrock, consistent with unusually large seasonal water table fluctuations (up to 50 m) observed in the upper part of the watershed. Apparent 3 H/ 3 He ages are positively correlated with sample depth and excess air concentrations. Integrated samples were collected from artesian bedrock wells near the trunk stream and are assumed to approximate flow‐weighted samples reflecting bedrock aquifer mean residence times. Exponential mean ages for these integrated samples are remarkably consistent along the stream, four of five being from 8 to 11 years. The tracer data in combination with other hydrologic and geologic data support a relatively simple conceptual model of groundwater flow in the watershed in which (1) permeability is primarily a function of depth; (2) water table fluctuations increase with distance from the stream; and (3) recharge, aquifer thickness, and porosity are relatively uniform throughout the watershed in spite of the geological complexity of the Proterozoic crystalline rocks that underlie it.

Water Resources Research↗

USGS2007

Identifying spatial variability of groundwater discharge in a wetland stream using a distributed temperature sensor

Discrete zones of groundwater discharge in a stream within a peat‐dominated wetland were identified on the basis of variations in streambed temperature using a distributed temperature sensor (DTS). The DTS gives measurements of the spatial (±1 m) and temporal (15 min) variation of streambed temperature over a much larger reach of stream (>800 m) than previous methods. Isolated temperature anomalies observed along the stream correspond to focused groundwater discharge zones likely caused by soil pipes within the peat. The DTS also recorded variations in the number of temperature anomalies, where higher numbers correlated well with a gaining reach identified by stream gauging. Focused zones of groundwater discharge showed essentially no change in position over successive measurement periods. Results suggest DTS measurements will complement other techniques (e.g., seepage meters and stream gauging) and help further improve our understanding of groundwater–surface water dynamics in wetland streams.

Water Resources Research↗

USGS2008

Effects of turbulence on hydraulic heads and parameter sensitivities in preferential groundwater flow layers

A conduit flow process (CFP) for the Modular Finite Difference Ground‐Water Flow model, MODFLOW‐2005, has been created by the U.S. Geological Survey. An application of the CFP on a carbonate aquifer in southern Florida is described; this application examines (1) the potential for turbulent groundwater flow and (2) the effects of turbulent flow on hydraulic heads and parameter sensitivities. Turbulent flow components were spatially extensive in preferential groundwater flow layers, with horizontal hydraulic conductivities of about 5,000,000 m d −1 , mean void diameters equal to about 3.5 cm, groundwater temperature equal to about 25°C, and critical Reynolds numbers less than or equal to 400. Turbulence either increased or decreased simulated heads from their laminar elevations. Specifically, head differences from laminar elevations ranged from about −18 to +27 cm and were explained by the magnitude of net flow to the finite difference model cell. Turbulence also affected the sensitivities of model parameters. Specifically, the composite‐scaled sensitivities of horizontal hydraulic conductivities decreased by as much as 70% when turbulence was essentially removed. These hydraulic head and sensitivity differences due to turbulent groundwater flow highlight potential errors in models based on the equivalent porous media assumption, which assumes laminar flow in uniformly distributed void spaces.

Water Resources Research↗

USGS2011

Transient groundwater chemistry near a river: Effects on U(VI) transport in laboratory column experiments

In the 300 Area of a U(VI)-contaminated aquifer at Hanford, Washington, USA, inorganic carbon and major cations, which have large impacts on U(VI) transport, change on an hourly and seasonal basis near the Columbia River. Batch and column experiments were conducted to investigate the factors controlling U(VI) adsorption/desorption by changing chemical conditions over time. Low alkalinity and low Ca concentrations (Columbia River water) enhanced adsorption and reduced aqueous concentrations. Conversely, high alkalinity and high Ca concentrations (Hanford groundwater) reduced adsorption and increased aqueous concentrations of U(VI). An equilibrium surface complexation model calibrated using laboratory batch experiments accounted for the decrease in U(VI) adsorption observed with increasing (bi)carbonate concentrations and other aqueous chemical conditions. In the column experiment, alternating pulses of river and groundwater caused swings in aqueous U(VI) concentration. A multispecies multirate surface complexation reactive transport model simulated most of the major U(VI) changes in two column experiments. The modeling results also indicated that U(VI) transport in the studied sediment could be simulated by using a single kinetic rate without loss of accuracy in the simulations. Moreover, the capability of the model to predict U(VI) transport in Hanford groundwater under transient chemical conditions depends significantly on the knowledge of real-time change of local groundwater chemistry.

Washington↗

USGS2021

Measuring and interpreting multilayer aquifer-system compactions for a sustainable groundwater-system development

Ever decreasing water resources and climate change have driven the increasing use of groundwater causing land subsidence in many countries. Geodetic sensors such as InSAR, GPS and leveling can detect surface deformation but cannot measure subsurface deformation. A single‐well, single‐depth extensometer can be used to measure subsurface deformation, but it cannot delineate the depths of major compaction and provide insight about the deformation mechanism throughout a complex aquifer system, unless man extensometers at different depths are used. We present a multilayer compaction well (MLCW), installed in a borehole, that uses magnetic rings to detect stratum compaction at 25 depths as deep as 300 m below land surface. Our laboratory and field assessments indicate 1 mm precision and accuracy for one single‐depth magnetic reading. We tested the performance of MLCW by measuring aquifer‐system compaction over the proximal, middle, and distal fans of the Choushui River Alluvial Fan (CRAF) that has long experienced severe land subsidence. The MLCW measurements were used to create time‐depth diagrams of compaction, showing different compaction rates at different layers of aquifers and aquitards to identify the depths of major compactions. The elastic (reversible) and inelastic (irreversible) compactions from MLCW were used in stress‐strain analyses to estimate skeletal specific storages and the safe groundwater levels, below which groundwater extractions have caused irreversible compactions. The hydrogeological parameters derived from MLCW measurements can help governmental agencies to determine effective land‐use and water‐use policies, and ascertain the best strategy for utilizing artificial recharge to prevent land subsidence and achieve sustainable groundwater management.

Water Resources Research↗

USGS2026

Decadal shifts in groundwater age detected by environmental tracers across California, USA

Groundwater age offers important insight into recharge, storage, and contamination risk. Although models predict age changes can be driven by pumping and climate variability, direct observational evidence remains limited. Here, we analyzed paired environmental tracer suites (tritium, carbon-14, and tritiogenic helium-3) collected a decade apart from 268 wells across California to assess the prevalence of groundwater age transience. Travel-time distribution models and statistical tests indicated age transience at 29% of sites, occurring most often in agricultural regions, such as the San Joaquin Valley and Southern Coast Ranges, where large carbon-14 changes coincided with substantial nitrate and chloride shifts. Sites with tritiogenic helium-3 data showed more frequent age transience, underscoring the value of multi-tracer data sets. These results provide the first regional evidence of widespread groundwater age change and a method for detecting changing water balances with implications for groundwater sustainability and water quality.

California↗

USGS2010

Gravity-driven groundwater flow and slope failure potential: 2. Effects of slope morphology, material properties, and hydraulic heterogeneity

Hillslope morphology, material properties, and hydraulic heterogeneities influence the role of groundwater flow in provoking slope instability. We evaluate these influences quantitatively by employing the elastic effective stress model and Coulomb failure potential concept described in our companion paper (Iverson and Reid, this issue). Sensitivity analyses show that of four dimensionless quantities that control model results (i.e., Poisson's ratio, porosity, topographic profile, and hydraulic conductivity contrast), slope profiles and hydraulic conductivity contrasts have the most pronounced and diverse effects on groundwater seepage forces, effective stresses, and slope failure potentials. Gravity-driven groundwater flow strongly influences the shape of equilibrium hillslopes, which we define as those with uniform near-surface failure potentials. For homogeneous slopes with no groundwater flow, equilibrium hillslope profiles are straight; but with gravity-driven flow, equilibrium profiles are concave or convex-concave, and the largest failure potentials exist near the bases of convex slopes. In heterogeneous slopes, relatively slight hydraulic conductivity contrasts of less than 1 order of magnitude markedly affect the seepage force field and slope failure potential. Maximum effects occur if conductivity contrasts are of four orders of magnitude or more, and large hydraulic gradients commonly result in particularly large failure potentials just upslope from where low-conductivity layers intersect the ground surface.

Water Resources Research↗

USGS2010

The role of groundwater transport in aquatic mercury cycling

Mercury, which is transported globally by atmospheric pathways to remote aquatic environments, is a ubiquitous contaminant at very low (nanograms Hg per liter) aqueous concentrations. Until recently, however, analytical and sampling techniques were not available for freshwater systems to quantify the actual levels of mercury concentrations without introducing significant contamination artifacts. Four different sampling strategies were used to evaluate ground water flow as a mercury source and transport mechanism within aquatic systems. The sampling strategies employ ultraclean techniques to determine mercury concentrations in groundwater and pore water near Pallette Lake, Wisconsin. Ambient groundwater concentrations are about 2–4 ng Hg L −1 , whereas pore waters near the sediment/water interface average about 12 ng Hg L −1 , emphasizing the importance of biogeochemical processes near the interface. Overall, the groundwater system removes about twice as much mercury (1.5 g yr −1 ) as it contributes (0.7 g yr −1 ) to Pallette Lake. About three fourths of the groundwater mercury load is recycled, thought to be derived from the water column.

Wisconsin↗

USGS1996

Direct simulation of groundwater age

A new method is proposed to simulate groundwater age directly, by use of an advection-dispersion transport equation with a distributed zero-order source of unit (1) strength, corresponding to the rate of aging. The dependent variable in the governing equation is the mean age, a mass-weighted average age. The governing equation is derived from residence-time-distribution concepts for the case of steady flow. For the more general case of transient flow, a transient governing equation for age is derived from mass-conservation principles applied to conceptual “age mass.” The age mass is the product of the water mass and its age, and age mass is assumed to be conserved during mixing. Boundary conditions include zero age mass flux across all noflow and inflow boundaries and no age mass dispersive flux across outflow boundaries. For transient-flow conditions, the initial distribution of age must be known. The solution of the governing transport equation yields the spatial distribution of the mean groundwater age and includes diffusion, dispersion, mixing, and exchange processes that typically are considered only through tracer-specific solute transport simulation. Traditional methods have relied on advective transport to predict point values of groundwater travel time and age. The proposed method retains the simplicity and tracer-independence of advection-only models, but incorporates the effects of dispersion and mixing on volume-averaged age. Example simulations of age in two idealized regional aquifer systems, one homogeneous and the other layered, demonstrate the agreement between the proposed method and traditional particle-tracking approaches and illustrate use of the proposed method to determine the effects of diffusion, dispersion, and mixing on groundwater age.

Water Resources Research↗

USGS1998

Detecting influential observations in nonlinear regression modeling of groundwater flow

Nonlinear regression is used to estimate optimal parameter values in models of groundwater flow to ensure that differences between predicted and observed heads and flows do not result from nonoptimal parameter values. Parameter estimates can be affected, however, by observations that disproportionately influence the regression, such as outliers that exert undue leverage on the objective function. Certain statistics developed for linear regression can be used to detect influential observations in nonlinear regression if the models are approximately linear. This paper discusses the application of Cook's D , which measures the effect of omitting a single observation on a set of estimated parameter values, and the statistical parameter DFBETAS, which quantifies the influence of an observation on each parameter. The influence statistics were used to (1) identify the influential observations in the calibration of a three-dimensional, groundwater flow model of a fractured-rock aquifer through nonlinear regression, and (2) quantify the effect of omitting influential observations on the set of estimated parameter values. Comparison of the spatial distribution of Cook's D with plots of model sensitivity shows that influential observations correspond to areas where the model heads are most sensitive to certain parameters, and where predicted groundwater flow rates are largest. Five of the six discharge observations were identified as influential, indicating that reliable measurements of groundwater flow rates are valuable data in model calibration. DFBETAS are computed and examined for an alternative model of the aquifer system to identify a parameterization error in the model design that resulted in overestimation of the effect of anisotropy on horizontal hydraulic conductivity.

Water Resources Research↗

USGS1998

The discharge of nitrate-contaminated groundwater from developed shoreline to marsh-fringed estuary

As residential development, on-site wastewater disposal, and groundwater contamination increase in the coastal zone, assessment of nutrient removal by soil and sedimentary processes becomes increasingly important. Nitrogen removal efficiency depends largely on the specific flow paths taken by groundwater as it discharges into nitrogen-limited estuarine waters. Shoreline salinity surveys, hydraulic studies, and thermal infrared imagery indicated that groundwater discharge into the Nauset Marsh estuary (Eastham, Massachusetts) occurred in high-velocity seeps immediately seaward of the upland-fringing salt marsh. Discharge was highly variable spatially and occurred through permeable, sandy sediments during low tide. Seepage chamber monitoring showed that dissolved inorganic nitrogen (principally nitrate) traversed nearly conservatively from the aquifer through shallow estuarine sediments to coastal waters at flux rates of 1–3 mmol m −2 h −1 . A significant relationship between pore water NO 3 -N concentrations and NO 3 -N flux rates may provide a rapid method of estimating nitrogen loading from groundwater to the water column.

Massachusets↗

USGS2012

Rapid intrusion of magma into wet rock: Groundwater flow due to pore pressure increases

Analytical and numerical solutions are developed to simulate the pressurization, expansion, and flow of groundwater contained within saturated, intact host rocks subject to sudden heating from the planar surface of an igneous intrusion. For most rocks, water diffuses more rapidly than heat, assuring that groundwater is not heated along a constant-volume pressure path and that thermal expansion and pressurization adjacent to the intrusion drives a flow that extends well beyond the heated region. The forcing parameter for pressurization and flow is α Δ T , where α is a thermal expansion coefficient reflecting the overall expansion of water heated through the temperature difference Δ T between the initial ambient and intrusive values. Pore pressure increases due to heating are greatest when the intrusion is emplaced rapidly and where the intrusive contact is impervious to groundwater contained in stiff, impermeable rocks with high thermal diffusivities and porosities. The maximum velocity of water flowing in pores decays with the inverse square root of time and is insensitive to hydraulic properties of the host rocks. Pressures are lessened and flow directions are reversed with the onset of hydrothermal convection. This occurs at times ranging from hours to weeks after onset of intrusion. As magma rises into near-surface rocks, steam can be generated. Solutions indicate that pressure increases and velocities are sensitive to the overall amount of expansion rather than the behavior of the water-steam transition. Both the overall thermal expansion coefficient α and the temperature difference Δ T are greater in shallow (<1 km) environments than in deep (∼5 km) ones. Thus, for rocks with similar transport properties, pressure increases due to heating are greatest in shallow environments. Although solutions can be applied to rocks with a wide variety of properties, pressure increases are calculated for compliant quartz-rich sedimentary rocks with a porosities between 1 and 20% and permeabilities between 1 darcy and 1 μdarcy, subject to temperature increases of 500 and 1000 K at depths ranging from 0.1 to 5 km in a region of hydrostatic pressures and normal geothermal gradient. Such rocks, with porosities greater than 5%, permeabilities less than a 0.1 mdarcy, and drained hydrostatic compressibilities of 10 −4 /MPa, undergo pressure increases greater than 10 MPa (100 bars)for conditions typical of water table depths of 2.5 km and heating to 500 K above ambient. Similar rocks, but with permeabilities less than 1 mdarcy, undergo pressure increases of 10 MPa for conditions typical of 1 km water table depth. Rocks commonly considered to be good aquifers undergo pressure increases of less than 1 MPa, primarily because of their high permeability. Although these estimates neglect the effects of fracturing and brecciation that may accompany such pressure increases, calculations indicate that pressure increases due to heating of cool groundwater can lead to failure of host rocks by a phreatic mechanism.

Journal of Geophysical Research Solid Earth↗

USGS2010

Estimating groundwater exchange with lakes: 1. The stable isotope mass balance method

Groundwater inflow and outflow contributions to the hydrologic budget of lakes can be determined using a stable isotope ( 18 O/ 16 O) mass balance method. The stable isotope method provides a way of integrating the spatial and temporal complexities of the flow field around a lake, thereby offering an appealing alternative to the traditional time and labor intensive methods using seepage meters and an extensive piezometer network. In this paper the method is applied to a lake in northern Wisconsin, demonstrating that it can be successfully applied to lakes in the upper midwest where thousands of similar lakes exist. Inflow and outflow rates calculated for the Wisconsin lake using the isotope mass balance method are 29 and 54 cm/yr, respectively, which compare well to estimates, derived independently using a three-dimensional groundwater flow and solute transport model, of 20 and 50 cm/yr. Such a favorable comparison lends confidence to the use of the stable isotope method to estimate groundwater exchange with lakes. In addition, utilization of stable isotopes in studies of groundwater-lake systems lends insight into mixing processes occurring in the unsaturated zone and in the aquifer surrounding the lake and verifies assumed flow paths based on head measurements in piezometers.

Wisconsin↗

USGS2010

Estimating groundwater exchange with lakes: 2. Calibration of a three-dimensional, solute transport model to a stable isotope plume

A three-dimensional groundwater flow and solute transport model was calibrated to a plume of water described by measurements of δ 18 O and used to calculate groundwater inflow and outflow rates at a lake in northern Wisconsin. The flow model was calibrated to observed hydraulic gradients and estimated recharge rates. Calibration of the solute transport submodel to the configuration of a stable isotope ( 18 O) plume in the contiguous aquifer on the downgradient side of the lake provides additional data to constrain the model. A good match between observed and simulated temporal variations in plume configuration indicates that the model closely simulated the dynamics of the real system. The model provides information on natural variations of rates of groundwater inflow, lake water outflow, and recharge to the water table. Inflow and outflow estimates compare favorably with estimates derived by the isotope mass balance method (Krabbenhoft et al., this issue). Model simulations agree with field observations that show groundwater inflow rates are more sensitive to seasonal variations in recharge than outflow.

Wisconsin↗

USGS2016

Understanding the hydrologic impacts of wastewater treatment plant discharge to shallow groundwater: Before and after plant shutdown

Effluent-impacted surface water has the potential to transport not only water, but wastewater-derived contaminants to shallow groundwater systems. To better understand the effects of effluent discharge on in-stream and near-stream hydrologic conditions in wastewater-impacted systems, water-level changes were monitored in hyporheic-zone and shallow-groundwater piezometers in a reach of Fourmile Creek adjacent to and downstream of the Ankeny (Iowa, USA) wastewater treatment plant (WWTP). Water-level changes were monitored from approximately 1.5 months before to 0.5 months after WWTP closure. Diurnal patterns in WWTP discharge were closely mirrored in stream and shallow-groundwater levels immediately upstream and up to 3 km downstream of the outfall, indicating that such discharge was the primary control on water levels before shutdown. The hydrologic response to WWTP shutdown was immediately observed throughout the study reach, verifying the far-reaching hydraulic connectivity and associated contaminant transport risk. The movement of WWTP effluent into alluvial aquifers has implications for potential WWTP-derived contamination of shallow groundwater far removed from the WWTP outfall.

Iowa↗

USGS2024

The dominance and growth of shallow groundwater resources in continuous permafrost environments

Water is a limited resource in Arctic watersheds with continuous permafrost because freezing conditions in winter and the impermeability of permafrost limit storage and connectivity between surface water and deep groundwater. However, groundwater can still be an important source of surface water in such settings, feeding springs and large aufeis fields that are abundant in cold regions and generating runoff when precipitation is rare. Whether groundwater is sourced from suprapermafrost taliks or deeper regional aquifers will impact water availability as the Arctic continues to warm and thaw. Previous research is ambiguous about the role of deep groundwater, leading to uncertainty regarding Arctic water availability and changing water resources. We analyzed chemistry and residence times of spring, stream, and river waters in the continuous permafrost zone of Alaska, spanning the mountains to the coastal plain. Water chemistry and age tracers show that surface waters are predominately sourced from recent precipitation and have short (<50 y) subsurface residence times. Remote sensing indicates trends in the areal extent of aufeis over the last 37 y, and correlations between aufeis extent and previous year summer temperature. Together, these data indicate that surface waters in continuous permafrost regions may be impacted by short flow paths and shallow suprapermafrost aquifers that are highly sensitive to climatic and hydrologic change over annual timescales. Despite the lack of connection to regional aquifers, continued warming and permafrost thaw may promote deepening of the shallow subsurface aquifers and creation of shallow taliks, providing some resilience to Arctic freshwater ecosystems.

Proceedings of the National Academy of Sciences↗