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

Progress made in groundwater flow and transport modeling

A special session, “Automated Parameter Identification and Sensitivity Analysis in Groundwater Flow and Transport,” was held at the Spring AGU meeting last May to explore recent scientific and technological developments in this expanding research area. The problems and approaches discussed are crucial to any field of Earth science in which models are used to represent physical processes and simulated values are compared with measured quantities. Although methods for automated parameter identification and sensitivity analysis were first developed at about the same time as numerical groundwater models, they are rarely used. Presentations at this special session clearly demonstrated the potential contributions of these methods for more insightful groundwater model development and calibration.

Eos, Transactions, American Geophysical Union↗

Groundwater inflow measurements in wetland systems

Our current understanding of wetlands is insufficient to assess the effects of past and future wetland loss. While knowledge of wetland hydrology is crucial, groundwater flows are often neglected or uncertain. In this paper, groundwater inflows were estimated in wetlands in southwestern Wisconsin using traditional Darcy's law calculations and three independent methods that included (1) stable isotope mass balances, (2) temperature profile modeling, and (3) numerical water balance modeling techniques. Inflows calculated using Darcy's law were lower than inflows estimated using the other approaches and ranged from 0.02 to 0.3 cm/d. Estimates obtained using the other methods generally were higher (0.1 to 1.1 cm/d) and showed similar spatial trends. An areal map of groundwater flux generated by the water balance model demonstrated that areas of both recharge and discharge exist in what is considered a regional discharge area. While each method has strengths and weaknesses, the use of more than one method can reduce uncertainty in the estimates.

Water Resources Research↗

Inferring shallow groundwater flow in saprolite and fractured rock using environmental tracers

The Ridge and Valley Province of eastern Tennessee is characterized by (1) substantial topographic relief, (2) folded and highly fractured rocks of various lithologies that have low primary permeability and porosity, and (3) a shallow residuum of medium permeability and high total porosity. Conceptual models of shallow groundwater flow and solute transport in this system have been developed but are difficult to evaluate using physical characterization or short‐term tracer methods due to extreme spatial variability in hydraulic properties. In this paper we describe how chlorofluorocarbon 12, 3 H, and 3 He were used to infer groundwater flow and solute transport in saprolite and fractured rock near Oak Ridge, Tennessee. In the shallow residuum, fracture spacings are <0.05 m, suggesting that concentrations of these tracers in fractures and in the matrix have time to diffusionally equilibrate. The relatively smooth nature of tracer concentrations with depth in the residuum is consistent with this model and quantitatively suggests recharge fluxes of 0.2 to 0.4 m yr −1 . In contrast, groundwater flow within the unweathered rock appears to be controlled by fractures with spacings of the order of 2 to 5 m, and diffusional equilibration of fractures and matrix has not occurred. For this reason, vertical fluid fluxes in the unweathered rock cannot be estimated from the tracer data.

Water Resources Research↗

Multiphase groundwater flow near cooling plutons

We investigate groundwater flow near cooling plutons with a computer program that can model multiphase flow, temperatures up to 1200°C, thermal pressurization, and temperature‐dependent rock properties. A series of experiments examines the effects of host‐rock permeability, size and depth of pluton emplacement, single versus multiple intrusions, the influence of a caprock, and the impact of topographically driven groundwater flow. We also reproduce and evaluate some of the pioneering numerical experiments on flow around plutons. Host‐rock permeability is the principal factor influencing fluid circulation and heat transfer in hydrothermal systems. The hottest and most steam‐rich systems develop where permeability is of the order of 10 −15 m 2 . Temperatures and life spans of systems decrease with increasing permeability. Conduction‐dominated systems, in which permeabilities are ≤10 −16 m 2 , persist longer but exhibit relatively modest increases in near‐surface temperatures relative to ambient conditions. Pluton size, emplacement depth, and initial thermal conditions have less influence on hydrothermal circulation patterns but affect the extent of boiling and duration of hydrothermal systems. Topographically driven groundwater flow can significantly alter hydrothermal circulation; however, a low‐permeability caprock effectively decouples the topographically and density‐driven systems and stabilizes the mixing interface between them thereby defining a likely ore‐forming environment.

Journal of Geophysical Research B: Solid Earth↗

Interaction between stream temperature, streamflow, and groundwater exchanges in alpine streams

Four alpine streams were monitored to continuously collect stream temperature and streamflow for periods ranging from a week to a year. In a small stream in the Colorado Rockies, diurnal variations in both stream temperature and streamflow were significantly greater in losing reaches than in gaining reaches, with minimum streamflow losses occurring early in the day and maximum losses occurring early in the evening. Using measured stream temperature changes, diurnal streambed infiltration rates were predicted to increase as much as 35% during the day (based on a heat and water transport groundwater model), while the measured increase in streamflow loss was 40%. For two large streams in the Sierra Nevada Mountains, annual stream temperature variations ranged from 0° to 25°C. In summer months, diurnal stream temperature variations were 30–40% of annual stream temperature variations, owing to reduced streamflows and increased atmospheric heating. Previous reports document that one Sierra stream site generally gains groundwater during low flows, while the second Sierra stream site may lose water during low flows. For August the diurnal streamflow variation was 11% at the gaining stream site and 30% at the losing stream site. On the basis of measured diurnal stream temperature variations, streambed infiltration rates were predicted to vary diurnally as much as 20% at the losing stream site. Analysis of results suggests that evapotranspiration losses determined diurnal streamflow variations in the gaining reaches, while in the losing reaches, evapotranspiration losses were compounded by diurnal variations in streambed infiltration. Diurnal variations in stream temperature were reduced in the gaining reaches as a result of discharging groundwater of relatively constant temperature. For the Sierra sites, comparison of results with those from a small tributary demonstrated that stream temperature patterns were useful in delineating discharges of bank storage following dam releases. Direct coupling may have occurred between streamflow and stream temperature for losing stream reaches, such that reduced streamflows facilitated increased afternoon stream temperatures and increased afternoon stream temperatures induced increased streambed losses, leading to even greater increases in both stream temperature and streamflow losses.

Water Resources Research↗

A model for managing sources of groundwater pollution

The waste disposal capacity of a groundwater system can be maximized while maintaining water quality at specified locations by using a groundwater pollutant source management model that is based upon linear programing and numerical simulation. The decision variables of the management model are solute waste disposal rates at various facilities distributed over space. A concentration response matrix is used in the management model to describe transient solute transport and is developed using the U.S. Geological Survey solute transport simulation model. The management model was applied to a complex hypothetical groundwater system. Large-scale management models were formulated as dual linear programing problems to reduce numerical difficulties and computation time. Linear programing problems were solved using a numerically stable, available code. Optimal solutions to problems with successively longer management time horizons indicated that disposal schedules at some sites are relatively independent of the number of disposal periods. Optimal waste disposal schedules exhibited pulsing rather than constant disposal rates. Sensitivity analysis using parametric linear programing showed that a sharp reduction in total waste disposal potential occurs if disposal rates at any site are increased beyond their optimal values.

Water Resources Research↗

Adaptation of the Carter-Tracy water influx calculation to groundwater flow simulation

The Carter-Tracy calculation for water influx is adapted to groundwater flow simulation with additional clarifying explanation not present in the original papers. The Van Everdingen and Hurst aquifer-influence functions for radial flow from an outer aquifer region are employed. This technique, based on convolution of unit-step response functions, offers a simple but approximate method for embedding an inner region of groundwater flow simulation within a much larger aquifer region where flow can be treated in an approximate fashion. The use of aquifer-influence functions in groundwater flow modeling reduces the size of the computational grid with a corresponding reduction in computer storage and execution time. The Carter-Tracy approximation to the convolution integral enables the aquifer influence function calculation to be made with an additional storage requirement of only two times the number of boundary nodes more than that required for the inner region simulation. It is a good approximation for constant flow rates but is poor for time-varying flow rates where the variation is large relative to the mean. A variety of outer aquifer region geometries, exterior boundary conditions, and flow rate versus potentiometric head relations can be used. The radial, transient-flow case presented is representative. An analytical approximation to the functions of Van Everdingen and Hurst for the dimensionless potentiometric head versus dimensionless time is given.

Water Resources Research↗

Nonlinear-regression groundwater flow modeling of a deep regional aquifer system

A nonlinear regression groundwater flow model, based on a Galerkin finite-element discretization, was used to analyze steady state two-dimensional groundwater flow in the areally extensive Madison aquifer in a 75,000 mi 2 area of the Northern Great Plains. Regression parameters estimated include intrinsic permeabilities of the main aquifer and separate lineament zones, discharges from eight major springs surrounding the Black Hills, and specified heads on the model boundaries. Aquifer thickness and temperature variations were included as specified functions. The regression model was applied using sequential F testing so that the fewest number and simplest zonation of intrinsic permeabilities, combined with the simplest overall model, were evaluated initially; additional complexities (such as subdivisions of zones and variations in temperature and thickness) were added in stages to evaluate the subsequent degree of improvement in the model results. It was found that only the eight major springs, a single main aquifer intrinsic permeability, two separate lineament intrinsic permeabilities of much smaller values, and temperature variations are warranted by the observed data (hydraulic heads and prior information on some parameters) for inclusion in a model that attempts to explain significant controls on groundwater flow. Addition of thickness variations did not significantly improve model results; however, thickness variations were included in the final model because they are fairly well defined. Effects on the observed head distribution from other features, such as vertical leakage and regional variations in intrinsic permeability, apparently were overshadowed by measurement errors in the observed heads. Estimates of the parameters correspond well to estimates obtained from other independent sources.

Great Plains, Madison aquifer↗

A compositional multiphase model for groundwater contamination by petroleum products: 2. Numerical solution

In this paper we develop a numerical solution to equations developed in part 1 (M. Y. Corapcioglu and A. L. Baehr, this issue) to predict the fate of an immiscible organic contaminant such as gasoline in the unsaturated zone subsequent to plume establishment. This solution, obtained by using a finite difference scheme and a method of forward projection to evaluate nonlinear coefficients, provides estimates of the flux of solubilized hydrocarbon constituents to groundwater from the portion of a spill which remains trapped in a soil after routine remedial efforts to recover the product have ceased. The procedure was used to solve the one-dimensional (vertical) form of the system of nonlinear partial differential equations defining the transport for each constituent of the product. Additionally, a homogeneous, isothermal soil with constant water content was assumed. An equilibrium assumption partitions the constituents between air, water, adsorbed, and immiscible phases. Free oxygen transport in the soil was also simulated to provide an upper bound estimate of aerobic biodgradation rates. Results are presented for a hypothetical gasoline consisting of eight groups of hydrocarbon constituents. Rates at which hydrocarbon mass is removed from the soil, entering either the atmosphere or groundwater, or is biodegraded are presented. A significant sensitivity to model parameters, particularly the parameters characterizing diffusive vapor transport, was discovered. We conclude that hydrocarbon solute composition in groundwater beneath a gasoline contaminated soil would be heavily weighted toward aromatic constituents like benzene, toluene, and xylene.

Water Resources Research↗

Tectonic contraction across Los Angeles after removal of groundwater pumping effects

After the 1987 Whittier Narrows and 1994 Northridge earthquakes revealed that blind thrust faults represent a significant threat to metropolitan Los Angeles , a network of 250 continuously recording global positioning system (GPS) stations was deployed to monitor displacements associated with deep slip on both blind and surface faults. Here we augment this GPS data with interferometric synthetic aperture radar imagery to take into account the deformation associated with groundwater pumping and strike-slip faulting. After removing these non-tectonic signals, we are left with 4.4 mm yr -1 of uniaxial contraction across the Los Angeles basin, oriented N 36° E (perpendicular to the major strike-slip faults in the area). This indicates that the contraction is primarily accommodated on thrust faults rather than on the northeast-trending strike–slip faults. We have found that widespread groundwater and oil pumping obscures and in some cases mimics the tectonic signals expected from the blind thrust faults. In the 40-km-long Santa Ana basin, groundwater withdrawal and re-injection produces 12 mm yr -1 of long-term subsidence, accompanied by an unprecedented seasonal oscillation of 55 mm in the vertical direction and 7 mm horizontally.

Nature↗

Quantification of sucralose in groundwater well drinking water by silylation derivatization and gas chromatography-mass spectrometry

Sucralose is an increasingly popular artificial sweetener and has been found in the environment in groundwater, surface water, and wastewater treatment plant effluent. Its chemical properties make it strongly recalcitrant in the environment and it has been used as a conservative tracer of human wastewater in recent years. Most current methods of sucralose analysis use high performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) instrumentation. In this study, we describe an analytical method using silylation derivatization and gas chromatography mass spectrometry (GC-MS) for the quantification of sucralose in groundwater samples. This method employs a deuterium-labeled internal standard to account for reaction and sample processing imprecision. The deuterated internal standard as well as experiments using negative ion chemical ionization gas chromatography-mass spectrometry strongly indicate that sucralose is derivatized at all five hydroxyl positions with trimethyl silyl groups. A previously developed GC-MS method with derivatization of sucralose in environmental samples did not employ an internal standard for quantification. As such, this method represents a more robust methodological approach for sucralose quantification in environmental samples. The method detection limit based on a set of 15 method blanks was calculated to be 21.8 ng L −1 , which is competitive with most methods in the published literature and sufficient to detect sucralose in water with ∼0.1% wastewater contribution. The method was applied to 37 groundwater samples from drinking water wells in California's Central Valley, in which five samples (13.5%) were found to contain sucralose at concentrations greater than the 21.8 ng L −1 detection limit.

Analytical Methods↗

GIS-based method for estimating surficial groundwater levels in coastal Virginia using limited information

In many coastal areas, high water tables are present, complicating installation of some stormwater best management practices (BMPs) that rely on infiltration. Regional estimates of the seasonal high water table (SHWT) often rely on sources such as soil surveys taken over a decade ago; these data are static and do not account for groundwater withdrawals or other anthropogenic impacts. To improve estimates of the SHWT, we developed a GIS-based methodology relying on surface water elevations. Data sources included a 1.5-m (5.0 ft) resolution Lidar-derived digital elevation model (DEM), aerial imagery, and publicly available shapefiles of water boundaries. Twenty-six groundwater monitoring wells were screened to eliminate well locations influenced by pumping, yielding 22 wells. In coastal Virginia, tidal water bodies and ditches form terminal boundaries for discharge from the water-table aquifers and permit water table elevations to be fixed at the landward boundaries of surface water bodies. Water table elevations interpolated from well data and boundary elevations were used to create a triangulated irregular network representing the water table elevations for November 2012, which was the date of the DEM. An adjustment factor, calculated from the highest recorded April water table depth from long-term groundwater monitoring data, was added to estimate the SHWT elevation. SHWT elevations were subtracted from the DEM to yield SHWT depth, which was compared with long-term monitoring well data, yielding an R2 value of 0.91. Residual errors were random, although the method underpredicted the highest expected SHWT and overpredicted the median SHWT. The SHWT depth map was validated by using water table depths from 57 soil borings at 10 different sites, and consistently matched observations better than available soil survey estimates. The SHWT depth map could be useful for BMP siting and feasibility studies in similar hydrogeological settings.

Virginia↗

Management of surface water and groundwater withdrawals to maintain environmental stream flows in Michigan

In 2008, the State of Michigan enacted legislation requiring that new or increased high-capacity withdrawals (greater than 100,000 gallons per day) from either surface water or groundwater be reviewed to prevent Adverse Resource Impacts (ARI). Science- based guidance was sought in defining how groundwater or surface-water withdrawals affect streamflow and in quantifying the relation between reduced streamflow and changes in stream ecology. The implementation of the legislation led to a risk-based system based on a gradient of risk, ecological response curves, and estimation of groundwater-surface water interaction. All Michigan streams are included in the legislation, and, accordingly, all Michigan streams were classified into management types defined by size of watershed, stream-water temperature, and predicted fish assemblages. Different streamflow removal percentages define risk-based thresholds allowed for each type. These removal percentages were informed by ecological response curves of characteristic fish populations and finalized through a legislative workgroup process. The assessment process includes an on-line screening tool that may be used to evaluate new or increased withdrawals against the risk-based zones and allows withdrawals that are not likely to cause an ARI to proceed to water-use registration. The system is designed to consider cumulative impacts of high-capacity withdrawals and to promote user involvement in water resource management by the establishment of water-user committees as cumulative withdrawals indicate greater potential for ARI in the watershed.

Michigan↗

The influence of groundwater chemistry on arsenic concentrations and speciation in a quartz sand and gravel aquifer

We examined the chemical reactions influencing dissolved concentrations, speciation, and transport of naturally occurring arsenic (As) in a shallow, sand and gravel aquifer with distinct geochemical zones resulting from land disposal of dilute sewage effluent. The principal geochemical zones were: (1) the uncontaminated zone above the sewage plume [350 μ M dissolved oxygen (DO), p H 5.9]; (2) the suboxic zone (5 μ M DO, p H 6.2, elevated concentrations of sewage-derived phosphate and nitrate); and (3) the anoxic zone [dissolved iron(II) 100–300 μ M, p H 6.5–6.9, elevated concentrations of sewage-derived phosphate]. Sediments are comprised of greater than 90% quartz but the surfaces of quartz and other mineral grains are coated with nanometer-size iron (Fe) and aluminum (Al) oxides and/or silicates, which control the adsorption properties of the sediments. Uncontaminated groundwater with added phosphate (620 μ M) was pumped into the uncontaminated zone while samples were collected 0.3 m above the injection point. Concentrations of As(V) increased from below detection (0.005 μ M) to a maximum of 0.07 μ M during breakthrough of phosphate at the sampling port; As(III) concentrations remained below detection. These results are consistent with the hypothesis that naturally occurring As(V) adsorbed to constituents of the coatings on grain surfaces was desorbed by phosphate in the injected groundwater. Also consistent with this hypothesis, vertical profiles of groundwater chemistry measured prior to the tracer test showed that dissolved As(V) concentrations increased along with dissolved phosphate from below detection in the uncontaminated zone to approximately 0.07 and 70 μ M, respectively, in the suboxic zone. Concentrations of As(III) were below detection in both zones. The anoxic zone had approximately 0.07 μ M As(V) but also had As(III) concentrations of 0.07–0.14 μ M, suggesting that release of As bound to sediment grains occurred by desorption by phosphate, reductive dissolution of Fe oxides, and reduction of As(V) to As(III), which adsorbs only weakly to the Fe-oxide-depleted material in the coatings. Results of reductive extractions of the sediments suggest that As associated with the coatings was relatively uniformly distributed at approximately 1 nmol/g of sediment (equivalent to 0.075 ppm As) and comprised 20%-50% of the total As in the sediments, determined from oxidative extractions. Quartz sand aquifers provide high-quality drinking water but can become contaminated when naturally occurring arsenic bound to Fe and Al oxides or silicates on sediment surfaces is released by desorption and dissolution of Fe oxides in response to changing chemical conditions.

Geochemical Transactions↗

Winter feeding, growth and condition of brown trout Salmo trutta in a groundwater-dominated stream

Winter can be a stressful period for stream-dwelling salmonid populations, often resulting in reduced growth and survival. Stream water temperatures have been identified as a primary mechanism driving reductions in fitness during winter. However, groundwater inputs can moderate water temperature and may reduce winter severity. Additionally, seasonal reductions in prey availability may contribute to decreased growth and survival, although few studies have examined food webs supporting salmonids under winter conditions. This study employed diet, stable isotope, and mark-recapture techniques to examine winter (November through March) feeding, growth, and condition of brown trout Salmo trutta in a groundwater-dominated stream (Badger Creek, Minnesota, USA). Growth was greater for fish &le; 150 mm (mean = 4.1 mg g &minus;1 day &minus;1 ) than for those 151&ndash;276 mm (mean = 1.0 mg g &minus;1 day &minus;1 ) during the winter season. Overall condition from early winter to late winter did not vary for fish &le;150 mm (mean relative weight (Wr) = 89.5) and increased for those 151&ndash;276 mm (mean Wr = 85.8 early and 89.4 late). Although composition varied both temporally and by individual, brown trout diets were dominated by aquatic invertebrates, primarily Amphipods, Dipterans, and Trichopterans. Stable isotope analysis supported the observations of the dominant prey taxa in stomach contents and indicated the winter food web was supported by a combination of allochthonous inputs and aquatic macrophytes. Brown trout in Badger Creek likely benefited from the thermal regime and increased prey abundance present in this groundwater-dominated stream during winter.

Minnesota↗

Examining the contradiction in 'sustainable urban growth': an example of groundwater sustainability

The environmental planning literature proposes a set of 'best management practices' for urban development that assumes improvement in environmental quality as a result of specific urban patterns. These best management practices, however, often do not recognise finite biophysical limits and social impacts that urban patterns alone cannot overcome. To shed light on this debate, we explore the effects of different degrees of urban clustering on groundwater levels using a coupled land-use change and groundwater-flow model. Our simulations show that specific urban forms only slow down the impact on groundwater. As population increases, the pattern in which it is accommodated ceases to matter, and widespread depletion ensues. These results are predictable, yet current planning practice tends to take growth for granted and is reluctant to envision either no-growth scenarios or the prospect of depletion. We propose to use simulations such as those presented here to aid in policy discussions that allow decision makers to question the assumption of sustainable growth and suggest alternative forms of development.

Journal of Environmental Planning and Management↗

Linkages between temperature, macroinvertebrates, and young-of-year Coho Salmon growth in surface-water and groundwater streams

Ecological resources for fishes in stream food webs shift over space and time, providing a complex template of available resources that can be used for growth. We tracked water temperature in conjunction with young-of-year Coho Salmon size, growth, and diet in 2 streams with contrasting thermal regimes: a groundwater stream with colder temperatures and lower thermal variability all year and a surface-water stream with greater thermal variability and warmer summer temperatures more conducive to young-of-year salmon growth. We hypothesized that fry emergence would occur when rearing conditions are optimal for growth and that, all else being equal, summer fish growth will be greater in the surface-water stream. Previous work on Coho Salmon phenology in these streams showed that peak fry emergence occurred at the same time in early summer in both streams. We measured salmon fry emergence in relation to thermal variability and macroinvertebrate prey availability with subsequent tracking of somatic growth, diet, and body size during the 1 st year of life in both streams. Macroinvertebrate prey availability was highest overall in the colder and thermally-stable groundwater stream than the surface-water stream. Prey availability was particularly high in the thalweg drift during peak fry emergence in the groundwater stream. There was no difference in Coho Salmon diet composition between streams, which included invertebrates from benthic, drift, and riparian habitats. We found no differences in young-of-year Coho Salmon body size, growth, or consumption between streams. Overall, our results suggest that large differences in thermal regimes do not necessarily translate to large differences in young-of-year Coho Salmon size, growth, or diet. Many variables can influence fish growth, and there is not always a direct connection between spatial and temporal dimensions of environmental variability and their cascading effects on young-of-year Coho Salmon growth during the 1 st summer of life.

Freshwater Science↗

An assessment of radon in groundwater in New York State

Abstract: A set of 317 samples collected from wells throughout New York State (excluding Long Island) from 2003 through 2008 was used to assess the distribution of radon gas in drinking water. Previous studies have documented high concentrations of radon in groundwater from granitic and metamorphic bedrock, but there have been only limited characterizations of radon in water from sedimentary rock and unconsolidated sand-and-gravel deposits in New York. Approximately 8% of the samples from bedrock wells exceed 89 Bq L -1 (eight times the proposed regulatory limit), but only 2% of samples from sand-and-gravel wells exceed 44 Bq L -1 . Specific metamorphic and sedimentary rock formations in New York are associated with the high radon concentrations, indicating that specific areas of New York could be targeted with efforts to reduce the risk of exposure to radon in groundwater. Additionally, radon in groundwater from the sand-and-gravel aquifers was found to be directly correlated to radon in indoor air when assessed by county.

New York↗