Geology topics

Groundwater & aquifers

Find USGS groundwater and aquifer publications. Search an aquifer name or place, compare source-supplied dates and open reports for hydrogeology, recharge and water-level studies.

Follow this topic with RSS

At least 811 records · Page 45Linked to original sources

USGS2017

Using groundwater age distributions to understand changes in methyl tert-butyl ether (MtBE) concentrations in ambient groundwater, northeastern United States

Temporal changes in methyl tert-butyl ether (MtBE) concentrations in groundwater were evaluated in the northeastern United States, an area of the nation with widespread low-level detections of MtBE based on a national survey of wells selected to represent ambient conditions. MtBE use in the U.S. peaked in 1999 and was largely discontinued by 2007. Six well networks, each representing specific areas and well types (monitoring or supply wells), were each sampled at 10 year intervals between 1996 and 2012. Concentrations were decreasing or unchanged in most wells as of 2012, with the exception of a small number of wells where concentrations continue to increase. Statistically significant increasing concentrations were found in one network sampled for the second time shortly after the peak of MtBE use, and decreasing concentrations were found in two networks sampled for the second time about 10 years after the peak of MtBE use. Simulated concentrations from convolutions of estimates for concentrations of MtBE in recharge water with age distributions from environmental tracer data correctly predicted the direction of MtBE concentration changes in about 65 percent of individual wells. The best matches between simulated and observed concentrations were found when simulating recharge concentrations that followed the pattern of national MtBE use. Some observations were matched better when recharge was modeled as a plume moving past the well from a spill at one point in time. Modeling and sample results showed that wells with young median ages and narrow age distributions responded more quickly to changes in the contaminant source than wells with older median ages and broad age distributions. Well depth and aquifer type affect these responses. Regardless of the timing of decontamination, all of these aquifers show high susceptibility for contamination by a highly soluble, persistent constituent.

Science of the Total Environment↗

USGS2016

Assessing groundwater depletion and dynamics using GRACE and InSAR: Potential and limitations

In the last decade, remote sensing of the temporal variation of ground level and gravity has improved our understanding of groundwater dynamics and storage. Mass changes are measured by GRACE (Gravity Recovery and Climate Experiment) satellites, whereas ground deformation is measured by processing synthetic aperture radar satellites data using the InSAR (Interferometry of Synthetic Aperture Radar) techniques. Both methods are complementary and offer different sensitivities to aquifer system processes. GRACE is sensitive to mass changes over large spatial scales (more than 100,000 km 2 ). As such, it fails in providing groundwater storage change estimates at local or regional scales relevant to most aquifer systems, and at which most groundwater management schemes are applied. However, InSAR measures ground displacement due to aquifer response to fluid-pressure changes. InSAR applications to groundwater depletion assessments are limited to aquifer systems susceptible to measurable deformation. Furthermore, the inversion of InSAR-derived displacement maps into volume of depleted groundwater storage (both reversible and largely irreversible) is confounded by vertical and horizontal variability of sediment compressibility. During the last decade, both techniques have shown increasing interest in the scientific community to complement available in situ observations where they are insufficient. In this review, we present the theoretical and conceptual bases of each method, and present idealized scenarios to highlight the potential benefits and challenges of combining these techniques to remotely assess groundwater storage changes and other aspects of the dynamics of aquifer systems.

Groundwater↗

USGS2016

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 (&sim;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↗

USGS2019

Hydrogeology of an alpine talus aquifer: Cordillera Blanca, Peru

The dramatic loss of glacial mass in low latitudes is causing shifts in downstream water availability and use during the driest months of the year. The world’s largest concentration of tropical glaciers lies in the Cordillera Blanca range of Peru, where glacial runoff is declining and regional stresses are emerging over water resources. Throughout the Cordillera Blanca, groundwater inputs from alpine meadow–talus complexes, locally known as pampas, supply proglacial streams with up to 80% of their flow during the region’s dry season. Structural knowledge of the pampa aquifers is needed to estimate their drainable groundwater storage capacity and residence time, to elucidate the role and importance of alpine groundwater storage in the regional water budget of the Cordillera Blanca. To understand the structure of these proglacial aquifers, multiple near-surface geophysical methods were implemented in a proglacial valley near dense networks of spring-fed tributaries. Geophysical results and borehole logs suggest groundwater is stored in a confined aquifer composed of buried talus deposits overlain by lacustrine clay, while deeper portions of the unit, 10–15 m in depth, are relatively clay-free and more hydraulically conductive. Based on these findings and assumptions of aquifer porosity, the pampas of the Callejon de Huaylas may store from 0.006 to 0.02 km 3 of groundwater. Furthermore, these findings suggest that the talus aquifers of the Cordillera Blanca were formed in proglacial lakes, followed by infilling with fine lacustrine sediments that confine lower units and allow for groundwater discharge to springs via macropores and preferential flow.

Hydrogeology Journal↗

USGS2014

Nitrogen speciation and trends, and prediction of denitrification extent, in shallow US groundwater

Uncertainties surrounding nitrogen cycling complicate assessments of the environmental effects of nitrogen use and our understanding of the global carbon–nitrogen cycle. In this paper, we synthesize data from 877 ambient-monitoring wells across the US to frame broad patterns of nitrogen speciation and trends. At these sites, groundwater frequently contains substantial co-occurring NO 3 − and XSN 2 (N 2 from denitrification), reflecting active/ongoing denitrification and/or a mixture of undenitrified and denitrified groundwater. NO 3 − and NH 4 + essentially do not co-occur, indicating that the dominant source of NH 4 + at these sites likely is not dissimilatory reduction of NO 3 − to NH 4 + . Positive correlations of NH 4 + with apparent age, CH 4 , dissolved organic carbon, and indicators of reduced conditions are consistent with NH 4 + mobilization from degradation of aquifer organic matter and contraindicate an anthropogenic source of NH 4 + for most sites. Glacial aquifers and eastern sand and gravel aquifers generally have lower proportions of NO 3 − and greater proportions of XSN 2 than do fractured rock and karst aquifers and western sand and gravel aquifers. NO 3 − dominates in the youngest groundwater, but XSN 2 increases as residence time increases. Temporal patterns of nitrogen speciation and concentration reflect (1) changing NO 3 − loads over time, (2) groundwater residence-time controls on NH 4 + mobilization from solid phases, and (3) groundwater residence-time controls on denitrification. A simple classification tree using readily available variables (a national coverage of soil water depth, generalized geology) or variables reasonably estimated in many aquifers (residence time) identifies categorical denitrification extent (<10%, 10–50%, and >50%) with 79% accuracy in an independent testing set, demonstrating a predictive application based on the interconnected effects of redox, geology, and residence time.

Journal of Hydrology↗

USGS2002

Denitrification in the recharge area and discharge area of a transient agricultural nitrate plume in a glacial outwash sand aquifer, Minnesota

Recharge rates of nitrate (NO 3 − ) to groundwater beneath agricultural land commonly are greater than discharge rates of NO 3 − in nearby streams, but local controls of NO 3 − distribution in the subsurface generally are poorly known. Groundwater dating (CFC, 3 H) was combined with chemical (ions and gases) and stable isotope (N, S, and C) analyses to resolve the effects of land use changes, flow patterns, and water‐aquifer reactions on the distributions of O 2 , NO 3 − , SO 4 = , and other constituents in a two‐dimensional vertical section leading from upland cultivated fields to a riparian wetland and stream in a glacial outwash sand aquifer near Princeton, Minnesota. Within this section a “plume” of oxic NO 3 − ‐rich groundwater was present at shallow depths beneath the fields and part of the wetland but terminated before reaching the stream or the wetland surface. Groundwater dating and hydraulic measurements indicate travel times in the local flow system of 0 to >40 years, with stratified recharge beneath the fields, downward diversion of the shallow NO 3 − ‐bearing plume by semiconfining organic‐rich valley‐filling sediments under the wetland and upward discharge across the valley and stream bottom. The concentrations and δ 15 N values of NO 3 − and N 2 indicate that the NO 3 − plume section was bounded in three directions by a curvilinear zone of active denitrification that limited its progress; however, when recalculated to remove the effects of denitrification, the data also indicate changes in both the concentrations and δ 15 N values of NO 3 − that was recharged in the past. Isotope data and mass balance calculations indicate that FeS 2 and other ferrous Fe phases were the major electron donors for denitrification in at least two settings: (1) within the glacial‐fluvial aquifer sediments beneath the recharge and discharge areas and (2) along the bottom of the valley‐filling sediments in the discharge area. Combined results indicate that the shape and progress of the oxic NO 3 − plume termination were controlled by a combination of (1) historical and spatial variations in land use practices, (2) contrast in groundwater flow patterns between the agricultural recharge area and riparian wetland discharge area, and (3) distribution and abundance of electron donors in both the sand aquifer and valley‐filling sediments. The data are consistent with slow migration of redox zones through the aquifer in response to recharging oxic groundwater during Holocene time, then an order‐of‐magnitude increase in the flux of electron acceptors as a result of agricultural NO 3 − contamination in the late twentieth century, to which the redox zone configuration still may be adjusting. The importance of denitrification for NO 3 − movement through formerly glaciated terrains should depend on the source areas and depositional environments of the glacial sediments, as well as geomorphology and recent stream‐valley sediment history.

Water Resources Research↗

USGS2013

Assessment of groundwater quality data for the Turtle Mountain Indian Reservation, Rolette County, North Dakota

The Turtle Mountain Indian Reservation relies on groundwater supplies to meet the demands of community and economic needs. The U.S. Geological Survey, in cooperation with the Turtle Mountain Band of Chippewa Indians, examined historical groundwater-level and groundwater-quality data for the Fox Hills, Hell Creek, Rolla, and Shell Valley aquifers. The two main sources of water-quality data for groundwater were the U.S. Geological Survey National Water Information System database and the North Dakota State Water Commission database. Data included major ions, trace elements, nutrients, field properties, and physical properties. The Fox Hills and Hell Creek aquifers had few groundwater water-quality data. The lack of data limits any detailed assessments that can be made about these aquifers. Data for the Rolla aquifer exist from 1978 through 1980 only. The concentrations of some water-quality constituents exceeded the U.S. Environmental Protection Agency secondary maximum contaminant levels. No samples were analyzed for pesticides and hydrocarbons. Numerous water-quality samples have been obtained from the Shell Valley aquifer. About one-half of the water samples from the Shell Valley aquifer had concentrations of iron, manganese, sulfate, and dissolved solids that exceeded the U.S. Environmental Protection Agency secondary maximum contaminant levels. Overall, the data did not indicate obvious patterns in concentrations.

North Dakota↗

USGS2013

Transient calibration of a groundwater-flow model of Chimacum Creek Basin and vicinity, Jefferson County, Washington: a supplement to Scientific Investigations Report 2013-5160

A steady-state groundwater-flow model described in Scientific Investigations Report 2013-5160, ”Numerical Simulation of the Groundwater-Flow System in Chimacum Creek Basin and Vicinity, Jefferson County, Washington” was developed to evaluate potential future impacts of growth and of water-management strategies on water resources in the Chimacum Creek Basin. This supplement to that report describes the unsuccessful attempt to perform a calibration to transient conditions on the model. The modeled area is about 64 square miles on the Olympic Peninsula in northeastern Jefferson County, Washington. The geologic setting for the model area is that of unconsolidated deposits of glacial and interglacial origin typical of the Puget Sound Lowlands. The hydrogeologic units representing aquifers are Upper Aquifer (UA, roughly corresponding to recessional outwash) and Lower Aquifer (LA, roughly corresponding to advance outwash). Recharge from precipitation is the dominant source of water to the aquifer system; discharge is primarily to marine waters below sea level and to Chimacum Creek and its tributaries. The model is comprised of a grid of 245 columns and 313 rows; cells are a uniform 200 feet per side. There are six model layers, each representing one hydrogeologic unit: (1) Upper Confining unit (UC); (2) Upper Aquifer unit (UA); (3) Middle Confining unit (MC); (4) Lower Aquifer unit (LA); (5) Lower Confining unit (LC); and (6) Bedrock unit (OE). The transient simulation period (October 1994–September 2009) was divided into 180 monthly stress periods to represent temporal variations in recharge, discharge, and storage. An attempt to calibrate the model to transient conditions was unsuccessful due to instabilities stemming from oscillations in groundwater discharge to and recharge from streamflow in Chimacum Creek. The model as calibrated to transient conditions has mean residuals and standard errors of 0.06 ft ±0.45 feet for groundwater levels and 0.48 ± 0.06 cubic feet per second for flows. Although the expected seasonal trends were observed in model results, the typical observed annual variation of groundwater levels of about 2 feet was not. Streamflow at the most downstream observation point was about three times larger than simulated streamflow. Because the transient version of the model proved inherently unstable, it was not used to simulate forecast conditions for alternate hydrologic or anthropogenic changes. Adaptation of alternate stream simulation packages, such as RIV, or newer versions of MODFLOW, such as MODFLOW-NWT, could possibly assist with achieving calibration to transient conditions.

Washington↗

USGS2016

Potential effects of alterations to the hydrologic system on the distribution of salinity in the Biscayne aquifer in Broward County, Florida

To address concerns about the effects of water-resource management practices and rising sea level on saltwater intrusion, the U.S. Geological Survey in cooperation with the Broward County Environmental Planning and Community Resilience Division, initiated a study to examine causes of saltwater intrusion and predict the effects of future alterations to the hydrologic system on salinity distribution in eastern Broward County, Florida. A three-dimensional, variable-density solute-transport model was calibrated to conditions from 1970 to 2012, the period for which data are most complete and reliable, and was used to simulate historical conditions from 1950 to 2012. These types of models are typically difficult to calibrate by matching to observed groundwater salinities because of spatial variability in aquifer properties that are unknown, and natural and anthropogenic processes that are complex and unknown; therefore, the primary goal was to reproduce major trends and locally generalized distributions of salinity in the Biscayne aquifer. The methods used in this study are relatively new, and results will provide transferable techniques for protecting groundwater resources and maximizing groundwater availability in coastal areas. The model was used to (1) evaluate the sensitivity of the salinity distribution in groundwater to sea-level rise and groundwater pumping, and (2) simulate the potential effects of increases in pumping, variable rates of sea-level rise, movement of a salinity control structure, and use of drainage recharge wells on the future distribution of salinity in the aquifer. Results from the simulation of historical conditions indicate that the model generally represents the observed greater westward extent of elevated salinity in the central part of the intruded area relative to the northern and southernmost parts of the intruded area. Results of sensitivity testing indicate that the extent of elevated salinity is most sensitive to pumping in areas where the source of saltwater is largely offshore, from the Atlantic Ocean, and is most sensitive to sea-level rise in areas where the source of salinity is downward leakage of brackish water from canals. Simulations of future scenarios indicate that increases in pumping near the existing interface may cause the interface to advance and decreases in pumping may cause it to retreat. Climatic effects, such as periods of prolonged drought or high precipitation, may augment or counteract long-term effects of changes in pumping on aquifer salinity at well fields. With increasing rates of sea-level rise, the freshwater-saltwater interface advances progressively inland, and flow-averaged salinities at well fields near the existing interface increase commensurately. Hypothetical southeastward (downstream) re-positioning of the existing G–54 salinity-control structure may prevent the interface from moving northwestward along and near the North New River canal, but beneficial effects are localized. Implementation of freshwater recharge wells in the city of Hallandale Beach may also have only a localized freshening effect in the aquifer and little appreciable effect on the freshwater-saltwater interface or on concentrations of salinity at well fields. Model accuracy and use are limited by uncertainty in the physical properties and boundary conditions of the system, uncertainty in historical and future conditions, and generalizations made in the mathematical relationships used to describe the physical processes of groundwater flow and transport. Because of these limitations, model results should be considered in relative rather than absolute terms. Nonetheless, model results do provide useful information on the relative scale of response of the system to changes in pumping distribution, sea-level rise, and mitigation activities.

Florida↗

USGS2012

Hydrochemical evolution of sodium-sulfate and sodium-chloride groundwater beneath the Northern Chihuahuan Desert, Trans-Pecos, Texas, USA

Groundwater beneath the northern Chihuahuan Desert, Trans-Pecos, Texas, USA, occurs in both carbonate and siliciclastic aquifers beneath a thick unsaturated zone in shallow Rio Grande alluvium. Groundwater hydrochemical evolution was investigated by analyzing soils, soils leachates, bolson-fill sediments, water from the unsaturated zone, and groundwater from three major aquifers. Ionic relations, mineral saturation states, and geochemical modeling show that groundwater compositions are controlled by reactions in the unsaturated zone, mineralogy of unsaturated sediments and aquifers, position in the groundwater flow system, and extensive irrigation. Recharge to aquifers unaffected by irrigation is initially a Ca-HCO 3 type as a result of dissolving carbonate surficial salts. With continued flow and mineral-water interaction, saturation with calcite and dolomite is maintained, gypsum is dissolved, and aqueous Ca and Mg are exchanged for adsorbed Na to produce a Na-SO 4 water. Groundwater in Rio Grande alluvium is a Na-Cl type, reflecting river-water composition and the effects of irrigation, evapotranspiration, and surficial salt recycling. These results document two hydrochemical evolution paths for groundwater in arid lands. If recharge is dilute precipitation, significant changes in water chemistry can occur in unsaturated media, ion exchange can be as important as dissolution-precipitation reactions in determining groundwater composition, and mineral-water reactions ultimately control groundwater composition. If recharge is return flow of irrigation water that already contains appreciable solutes, mineral-water reactions are less important than irrigation-water composition in determining groundwater chemistry.

Hydrogeology Journal↗

USGS2000

Nutrients in groundwaters of the conterminous United States, 1992-1995

Results of a national water quality assessment indicate that nitrate is detected in 71% of groundwater samples, more than 13 times as often as ammonia, nitrite, organic nitrogen, and orthophosphate, based on a common detection threshold of 0.2 mg/L. Shallow groundwater (typically 5 m deep or less) beneath agricultural land has the highest median nitrate concentration (3.4 mg/L), followed by shallow groundwater beneath urban land (1.6 mg/L) and deeper groundwater in major aquifers (0.48 mg/L). Nitrate exceeds the maximum contaminant level, 10 mg/L as nitrogen, in more than 15% of groundwater samples from 4 of 33 major aquifers commonly used as a source of drinking water. Nitrate concentration in groundwater is variable and depends on interactions among several factors, including nitrogen loading, soil type, aquifer permeability, recharge rate, and climate. For a given nitrogen loading, factors that generally increase nitrate concentration in groundwater include well-drained soils, fractured bedrock, and irrigation. Factors that mitigate nitrate contamination of groundwater include poorly drained soils, greater depth to groundwater, artificial drainage systems, intervening layers of unfractured bedrock, a low rate of groundwater recharge, and anaerobic conditions in aquifers.

Environmental Science & Technology↗

USGS2014

Base of principal aquifer for parts of the North Platte, South Platte, and Twin Platte Natural Resources Districts, western Nebraska

Water resources in the North and South Platte River valleys of Nebraska, including the valley of Lodgepole Creek, are critical to the social and economic health of the area, and for the recovery of threatened and endangered species in the Platte River Basin. Groundwater and surface water are heavily used resources, and uses are regulated in the study area. Irrigation is the dominant water use and, in most instances, is supplied by both groundwater and surface-water sources. The U.S. Geological Survey and its partners have collaborated to use airborne geophysical surveys for areas of the North and South Platte River valleys including the valley of Lodgepole Creek in western Nebraska. The objective of the surveys was to map the aquifers and underlying bedrock topography of selected areas to help improve the understanding of groundwater–surface-water relations to guide water-management decisions. This project was a cooperative study involving the North Platte Natural Resources District, the South Platte Natural Resources District, the Twin Platte Natural Resources District, the Conservation and Survey Division of the University of Nebraska-Lincoln, and the Nebraska Environmental Trust. This report presents the interpreted base-of-aquifer surface for part of the area consisting of the North Platte Natural Resources District, the South Platte Natural Resources District, and the Twin Platte Natural Resources District. The interpretations presented herein build on work done by previous researchers from 2008 to 2009 by incorporating additional airborne electromagnetic survey data collected in 2010 and additional test holes from separate, related studies. To make the airborne electromagnetic data useful, numerical inversion was used to convert the measured data into a depth-dependent subsurface resistivity model. An interpretation of the elevation and configuration of the base of aquifer was completed in a geographic information system that provided x, y, and z coordinates. The process of interpretation involved manually picking locations (base-of-aquifer elevations) on the displayed airborne electromagnetic-derived resistivity profile by the project geophysicist, hydrologist, and geologist. These locations, or picks, of the base-of-aquifer elevation (typically the top of the Brule Formation of the White River Group) were then stored in a georeferenced database. The pick was made by comparing the inverted airborne electromagnetic-derived resistivity profile to the lithologic descriptions and borehole geophysical logs from nearby test holes. The database of interpretive picks of the base-of-aquifer elevation was used to create primary input for interpolating the new base-of-aquifer contours. The automatically generated contours were manually adjusted based on the interpreted location of paleovalleys eroded into the base-of-aquifer surface and associated bedrock highs, many of which were unmapped before this study. When contours are overlain by the water-table surface, the saturated thickness of the aquifer can be computed, which allows an estimate of total water in storage. The contours of the base-of-aquifer surface presented in this report may be used as the lower boundary layer in existing and future groundwater-flow models. The integration of geophysical data into the contouring process facilitated a more continuous and spatially comprehensive view of the hydrogeologic framework.

Nebraska↗

USGS2012

Numerical model simulations of nitrate concentrations in groundwater using various nitrogen input scenarios, mid-Snake region, south-central Idaho

As part of the U.S. Geological Survey’s National Water Quality Assessment (NAWQA) program nitrate transport in groundwater was modeled in the mid-Snake River region in south-central Idaho to project future concentrations of nitrate. Model simulation results indicated that nitrate concentrations would continue to increase over time, eventually exceeding the U.S. Environmental Protection Agency maximum contaminant level for drinking water of 10 milligrams per liter in some areas. A subregional groundwater model simulated the change of nitrate concentrations in groundwater over time in response to three nitrogen input scenarios: (1) nitrogen input fixed at 2008 levels; (2) nitrogen input increased from 2008 to 2028 using the same rate of increase as the average rate of increase during the previous 10 years (1998 through 2008); after 2028, nitrogen input is fixed at 2028 levels; and (3) nitrogen input related to agriculture completely halted, with only nitrogen input from precipitation remaining. Scenarios 1 and 2 project that nitrate concentrations in groundwater continue to increase from 10 to 50 years beyond the year nitrogen input is fixed, depending on the location in the model area. Projected nitrate concentrations in groundwater increase by as much as 2–4 milligrams per liter in many areas, with nitrate concentrations in some areas reaching 10 milligrams per liter. Scenario 3, although unrealistic, estimates how long (20–50 years) it would take nitrate in groundwater to return to background concentrations—the “flushing time” of the system. The amount of nitrate concentration increase cannot be explained solely by differences in nitrogen input; in fact, some areas with the highest amount of nitrogen input have the lowest increase in nitrate concentration. The geometry of the aquifer and the pattern of regional groundwater flow through the aquifer greatly influence nitrate concentrations. The aquifer thins toward discharge areas along the Snake River which forces upward convergence of good-quality regional groundwater that mixes with the nitrate-laden groundwater in the uppermost parts of the aquifer, which results in lowered nitrate concentrations. A new method of inputting nitrogen to the subregional groundwater model was used that prorates nitrogen input by the probability of detecting nitrate concentrations greater than 2 mg/L. The probability map is based on correlations with physical factors, and prorates an existing nitrogen input dataset providing an estimate of nitrogen flux to the water table that accounts for new factors such as soil properties. The effectiveness of this updated nitrogen input method was evaluated using the software UCODE_2005.

Idaho↗

USGS2014

Estimates of groundwater recharge rates and sources in the East Mountain area, Eastern Bernalillo County, New Mexico, 2005-12

The U.S. Geological Survey, in cooperation with the Bernalillo County Public Works Division, has conducted a monitoring program in the East Mountain area of eastern Bernalillo County, New Mexico, since 2000 to better define the hydrogeologic characteristics of the East Mountain area and to provide scientific information that will assist in the sustainable management of water resources. This report presents estimates of groundwater recharge to the aquifers that supply water to a network of springs that discharged within the East Mountain area of eastern Bernalillo County during 2005–12. Chloride concentration, the mass ratio of chloride to bromide, and the stable isotope ratios of hydrogen and oxygen were used to estimate annual groundwater recharge rates and to identify the sources and timing of recharge to the aquifers in the East Mountain area. Groundwater recharge rates were estimated by using a chloride mass-balance (CMB) method applied to data from selected springs located in the study area. Eleven springs and four downgradient monitoring wells were sampled for this study. On the basis of chloride concentrations and the mass ratio of chloride to bromide, eight of the eleven sampled springs are considered representative of dilute groundwater recharged by meteoric water in the Sandia Mountains. Eight of the eleven springs sampled as part of this investigation are considered representative of dilute groundwater not influenced by nonmeteoric chloride sources on the basis of analysis of chloride concentrations and the mass ratio of chloride to bromide. Chloride concentrations at three of the sampled springs were likely affected by nonmeteoric chloride sources. Results of CMB calculations for the eight springs with Cl:Br ratios and chloride concentrations within the range of dilute groundwater (not influenced by nonmeteoric chloride sources) indicated that between about 5.5 and 23 percent of annual precipitation recharges the groundwater system. The variation in estimated recharge rates indicated that the mechanisms for recharge and groundwater movement in the East Mountain area are complex and that factors such as climate variability, the extent and interconnection of structural features such as faults and fractures, and potential solution enhancement of the aquifers all play important roles in the rates and timing of recharge. Stable isotope data from springs and snowpacks sampled in the East Mountain area were compared with local, regional, and global meteoric water lines and were analyzed along with values representing the stable isotope composition of winter precipitation and summer monsoonal rains. Results of the stable isotope analysis from springs in this study suggested that winter precipitation is the primary source of groundwater recharge to the aquifers supplying the springs, but there is a component of more isotopically enriched precipitation being recharged as well, likely from summer monsoonal rains. Specific conductance, groundwater-level hydrographs, snowpack chemistry, and snow-water equivalent data were used to inform the analyses and corroborate the findings of the CMB and stable isotope results.

New Mexico↗

USGS2023

Hydrologic effects of possible changes in water-supply withdrawals from, and effluent recharge to, the Kirkwood-Cohansey aquifer system, Winslow Township, Camden County, New Jersey

Winslow Township and the Camden County Municipal Utility Authority (CCMUA) developed a plan to shut down the Winslow sewage-treatment facility and associated effluent infiltration facility and transfer the effluent to the CCMUA sewage-treatment facility on the Delaware River in Camden, New Jersey. Winslow Township reduced groundwater withdrawals from the Kirkwood-Cohansey aquifer system to offset groundwater recharge lost with the cessation of effluent infiltration. The U.S. Geological Survey, in cooperation with Winslow Township and the CCMUA, collected data to evaluate conditions prior to cessation of effluent infiltration and installed two continuous-record streamflow-gaging stations. Streamflow measurements also were made at two low-flow partial-record sites, and groundwater levels were measured in 17 wells at high and low water-level periods (May and September 2010). A groundwater-flow model provides estimated changes in base flow of the Great Egg Harbor River under several groundwater-withdrawal and effluent infiltration scenarios. Water levels were measured in an observation well 480 feet (ft) from the infiltration lagoons during 1971–2010. A downward trend in water levels in the well prior to 1985 is attributed in part to increased impervious surfaces and groundwater withdrawals associated with development in the area that began in the early 1970s. From late 1985 to 2010, there was an upward trend in water levels in the well that is attributed to the construction of nearby effluent infiltration lagoons in 1985 and the increasing rate of effluent infiltration during the period. Recent and historical measurements made at the four surface-water sites were correlated with same-day discharges measured at three nearby index stations to estimate continuous low-flow record at the sites. Effects on base flow caused by reductions in groundwater withdrawals or the cessation of effluent infiltration in Winslow Township could not be ascertained from the available data with the statistical and analysis methods used. Groundwater discharge to streams (base flow) was simulated with a groundwater-flow model of the Great Egg Harbor and Mullica River Basins. Simulated monthly base flows using 2008–10 withdrawal rates and effluent recharge (Scenario 1) are generally about 1.5 million gallons per day (Mgal/d) greater than simulated base flows using 2003–07 withdrawal rates (Baseline Scenario) because of the 1.57 Mgal/d reduction in average withdrawals by Winslow Township from the Kirkwood-Cohansey aquifer system from 2003–07 to 2008–10. Simulated monthly base flows using 2008–10 withdrawals but without effluent infiltration (Scenario 2) are very similar to, but typically slightly lower than, Baseline Scenario base flows. Three hypothetical future distributions of groundwater withdrawals from existing Winslow Township wells are simulated, each without effluent infiltration and using the same groundwater withdrawal rate as Scenario 2, but with different hypothetical distributions of withdrawals among existing Winslow Township wells. The Scenario 3 and 4 base flows are greater than the Baseline Scenario base flows in all months, and the Scenario 5 base flows are less than the Baseline Scenario base flows in all months. The simulation results indicate that a reduction in average withdrawals from the Kirkwood-Cohansey aquifer system by 1.57 Mgal/d offsets the reduction of effluent infiltration by about the same rate, resulting in nearly unchanged base flows in the Great Egg Harbor River near Blue Anchor (01410820).

New Jersey↗

USGS1989

Yield and quality of ground water from stratified-drift aquifers, Taunton River basin, Massachusetts : executive summary

Water shortages are a chronic problem in parts of the Taunton River basin and are caused by a combination of factors. Water use in this part of the Boston metropolitan area is likely to increase during the next decade. The Massachusetts Division of Water Resources projects that about 50% of the cities and towns within and on the perimeter of the basin may have water supply deficits by 1990 if water management projects are not pursued throughout the 1980s. Estimates of the long-term yield of the 26 regional aquifers indicate that the yields of the two most productive aquifers equal or exceed 11.9 and 11.3 cu ft/sec, 90% of the time, respectively, if minimum stream discharge is maintained at 99.5% flow duration. Eighteen of the 26 aquifers were pumped for public water supply during 1983. Further analysis of the yield characteristics of these 18 aquifers indicates that the 1983 pumping rate of each of these 18 aquifers can be sustained at least 70% of the time. Selected physical properties and concentrations of major chemical constituents in groundwater from the stratified-drift aquifers at 80 sampling sites were used to characterize general water quality in aquifers throughout the basin. The pH of the groundwater ranged from 5.4 to 7.0. Natural elevated concentrations of Fe and Mn in water in the stratified-drift aquifers are present locally in the basin. Natural concentrations of these two metals commonly exceed the limits of 0.3 mg/L for Fe and 0.05 mg/L for Mn recommended for drinking water. Fifty-one analyses of selected trace metals in groundwater samples from stratified-drift aquifers throughout the basin were used to characterize trace metal concentrations in the groundwater. Of the 10 constituents sampled that have US EPA limits recommended for drinking water, only the Pb concentration in water at one site (60 micrograms/L) exceeded the recommended limit of 50 micrograms/L. Analyses of selected organic compounds in water in the stratified-drift aquifers at 74 locations revealed that 13 of the sample contained one or more of the following compounds: chloroform; carbon tetrachloride; dichloroethane; dichloroethylene; tetrachloroethylene; and, toluene. (Lantz-PTT)

Water-Resources Investigations Report↗

USGS2010

Death Valley regional groundwater flow system, Nevada and California: Hydrogeologic framework and transient groundwater flow model

A numerical three-dimensional (3D) transient groundwater flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the groundwater flow system and previous less extensive groundwater flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect groundwater flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley regional groundwater flow system (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the groundwater flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural groundwater discharge occurring through evapotranspiration (ET) and spring flow; the history of groundwater pumping from 1913 through 1998; groundwater recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional groundwater flow system. Groundwater flow in the Death Valley region is composed of several interconnected, complex groundwater flow systems. Groundwater flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional groundwater flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and ET are the dominant natural groundwater discharge processes. Groundwater also is withdrawn for agricultural, commercial, and domestic uses. Groundwater flow in the DVRFS was simulated using MODFLOW-2000, the U.S. Geological Survey 3D finitedifference modular groundwater flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 meters (m) on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins simulated by constant-head boundaries. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model generally fit observed heads reasonably well (residuals with absolute values less than 10 meters) with two exceptions: in most areas of nearly flat hydraulic gradient the fit is considered moderate (residuals with absolute values of 10 to 20 meters), and in areas of steep hydraulic gradient along the Eleana Range and western part of Yucca Flat, southern part of the Owlshead Mountains, southern part of the Bullfrog Hills, and the north-northwestern part of the model domain (residuals with absolute values greater than 20 meters). Groundwater discharge residuals are fairly random, with as many areas where simulated flows are less than observed flows as areas where simulated flows are greater. The highest unweighted groundwater discharge residuals occur at Death Valley, Sarcobatus Flat (northeastern area), Tecopa, and early observations at Manse Spring in Pahrump Valley. High weighted-discharge residuals were computed in Indian Springs Valley and parts of Death Valley. Most of these inaccuracies in head and discharge can be attributed to insufficient representation of the hydrogeology in the HFM and(or) discharge estimates, misrepresentation of water levels, and(or) model error associated with grid-cell size. The model represents the large and complex groundwater flow system of the Death Valley region at a greater degree of refinement and accuracy than has been possible previously. The representation of detail provided by the 3D digital hydrogeologic framework model and the numerical groundwater flow model enabled greater spatial accuracy in every model parameter. The lithostratigraphy and structural effects of the hydrogeologic framework; recharge estimates from simulated net infiltration; discharge estimates from ET, spring flow, and pumping; and boundary inflow and outflow estimates all were reevaluated, some additional data were collected, and accuracy was improved. Uncertainty in the results of the flow model simulations can be reduced by improving on the quality, interpretation, and representation of the water-level and discharge observations used to calibrate the model and improving on the representation of the HGU geometries, the spatial variability of HGU material properties, the flow model physical framework, and the hydrologic conditions.

California, Nevada↗

USGS2022

Hydrogeology and groundwater quality in the San Agustin Basin, New Mexico, 1975–2019

This report describes the findings of a U.S. Geological Survey study, completed in cooperation with the Bureau of Land Management, focused on better understanding the present-day (1975–2019) hydrogeology and groundwater quality of the San Agustin Basin in west-central New Mexico to support sustainable groundwater resource management. The basin hosts a relatively undeveloped basin-fill and alluvium aquifer system and is topographically divided into east and west subbasins by the McClure Hills. Groundwater chemistry and groundwater elevation data were compiled, collected, and interpreted in the context of groundwater flow and quality. The analyses presented in this report consider groundwater chemistry data collected within the last decade (2010–19) and groundwater elevation data collected from 1975 through 2019 to provide insight into present-day conditions. Groundwater elevations show that groundwater typically moves from the highlands to the lowlands, with a prominent east to west regional trend. Groundwater elevations were lowest in the southwestern portion of the west subbasin, where estimated flow directions suggest underflow through the local highlands into the northern East Fork Gila River watershed, which is further supported by historical groundwater elevation data from the northern East Fork Gila River watershed. Gradual groundwater elevation gradients (about 2 feet per mile) near the east and west subbasin divide suggest that groundwater slowly flows from the east subbasin to the west subbasin. Quantitative analyses of groundwater chemistry data show that groundwater in both subbasins has similar chemical characteristics. A systematic east to west groundwater evolution in water chemistry was not observed despite evidenced subbasin connectivity. The absence of this pattern suggests that groundwater mixing is regionally prevalent, sediment reactivity is low and variable, and (or) recharge conditions are comparable in both subbasins. Groundwater chemistry was generally independent of aquifer type, suggesting that the aquifers are hydrologically well connected. Corrected carbon-14 groundwater age estimates in the basin ranged from 232 to 13,916 years before present with a median of 5,409 years. A wide range of groundwater ages is therefore present in the basin, with waters commonly being thousands of years old, thereby supporting generally slow regional groundwater movement. A component of relatively young groundwater, for which estimated ages could not be accurately computed, is also present in the basin, and it may commonly mix with older waters. The spatial distribution of categorical and quantitative groundwater ages indicates that most recharge likely occurs in the highlands through mountain-block recharge and as focused recharge within arroyos, although evidence of modern (1953 and after) groundwater was minimal at sampled sites. Median annual gradients (groundwater elevation change over time) indicate that most groundwater elevations in the lowlands changed little (−0.2 to 0.2 foot per year) from 1975 through 2019. Groundwater elevations in the highlands varied more annually, which is likely due to recharge from precipitation events. These more variable groundwater elevations in the highlands compared with the lowlands, along with groundwater ages, provide further evidence that most groundwater recharge takes place in the highlands, with minimal recharge in the lowlands. Median groundwater elevation change for all sites was −0.05 foot per year. Temporal consistency of lowland groundwater elevations suggests that regional groundwater dynamics have been more or less stable through time under current climate and development conditions, although median annual gradients indicate that groundwater elevations may have slightly declined on average between 1975 and 2019.

New Mexico↗