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 217 records · Page 12Linked to original sources

USGS2021

Groundwater and surface-water data from the C-aquifer monitoring program, Northeastern Arizona, 2012–2019

The Coconino aquifer (C aquifer) is a regionally extensive multiple-aquifer system supplying water for municipal, agricultural, and industrial use in northeastern Arizona, northwestern New Mexico, and southeastern Utah. This report focuses on the C aquifer in the arid to semi-arid area between St. Johns, Ariz., and Flagstaff, Ariz., along the Interstate-40 corridor where an increase in groundwater withdrawals coupled with ongoing drought conditions increase the potential for substantial water-level decline within the aquifer. The U.S. Geological Survey (USGS) C-aquifer Monitoring Program began in 2005 to establish baseline groundwater and surface-water conditions and to quantify physical and water-chemistry responses to pumping stresses and climate. This report presents data previously reported in Brown and Macy (2012) that extend back as far as the 1950s, along with new data collected from the USGS C-aquifer Monitoring Program since that publication, from water years 2012 to 2019. Water levels in 17 wells are measured quarterly as part of the C-aquifer Monitoring Program, and five of those are continuously monitored at 15-minute intervals. Water levels in an additional 18 wells in the study area are measured periodically by the USGS or other agencies. The largest historical change in water level in the study area was a decrease of 81.20 feet in Lake Mary 1 Well near Flagstaff between 1962 and 2018. Changes in water levels were greatest around major pumping centers and in the eastern extent of the study area. Surface-water water-quality parameters (pH, water temperature, specific conductance, and dissolved oxygen) and streamflow discharge measurements were collected and analyzed along perennial, groundwater-fed reaches of Clear Creek, Chevelon Creek, and the Little Colorado River during nine baseflow investigations of varying extent between 2005 and 2019. Both Clear Creek and Chevelon Creek gain in flow from the beginning of their perennial reaches to their outflow into the Little Colorado River. The Little Colorado River has relatively steady streamflow in the reach between where the two tributaries enter the river. Chevelon Creek showed an increase in median specific conductance during all baseflow investigations of nearly 4,000 microsiemens per centimeter (μS/cm) from near the headwaters to the confluence with the Little Colorado River; Clear Creek also showed an increase in median specific conductance of almost 5,000 μS/cm from headwaters to confluence. Water temperature, dissolved oxygen, and pH do not show substantial trends along the reaches of Clear Creek, Chevelon Creek, or the Little Colorado River.

Arizona↗

USGS2010

Steady-state and transient models of groundwater flow and advective transport, Eastern Snake River Plain aquifer, Idaho National Laboratory and vicinity, Idaho

Three-dimensional steady-state and transient models of groundwater flow and advective transport in the eastern Snake River Plain aquifer were developed by the U.S. Geological Survey in cooperation with the U.S. Department of Energy. The steady-state and transient flow models cover an area of 1,940 square miles that includes most of the 890 square miles of the Idaho National Laboratory (INL). A 50-year history of waste disposal at the INL has resulted in measurable concentrations of waste contaminants in the eastern Snake River Plain aquifer. Model results can be used in numerical simulations to evaluate the movement of contaminants in the aquifer. Saturated flow in the eastern Snake River Plain aquifer was simulated using the MODFLOW-2000 groundwater flow model. Steady-state flow was simulated to represent conditions in 1980 with average streamflow infiltration from 1966-80 for the Big Lost River, the major variable inflow to the system. The transient flow model simulates groundwater flow between 1980 and 1995, a period that included a 5-year wet cycle (1982-86) followed by an 8-year dry cycle (1987-94). Specified flows into or out of the active model grid define the conditions on all boundaries except the southwest (outflow) boundary, which is simulated with head-dependent flow. In the transient flow model, streamflow infiltration was the major stress, and was variable in time and location. The models were calibrated by adjusting aquifer hydraulic properties to match simulated and observed heads or head differences using the parameter-estimation program incorporated in MODFLOW-2000. Various summary, regression, and inferential statistics, in addition to comparisons of model properties and simulated head to measured properties and head, were used to evaluate the model calibration. Model parameters estimated for the steady-state calibration included hydraulic conductivity for seven of nine hydrogeologic zones and a global value of vertical anisotropy. Parameters estimated for the transient calibration included specific yield for five of the seven hydrogeologic zones. The zones represent five rock units and parts of four rock units with abundant interbedded sediment. All estimates of hydraulic conductivity were nearly within 2 orders of magnitude of the maximum expected value in a range that exceeds 6 orders of magnitude. The estimate of vertical anisotropy was larger than the maximum expected value. All estimates of specific yield and their confidence intervals were within the ranges of values expected for aquifers, the range of values for porosity of basalt, and other estimates of specific yield for basalt. The steady-state model reasonably simulated the observed water-table altitude, orientation, and gradients. Simulation of transient flow conditions accurately reproduced observed changes in the flow system resulting from episodic infiltration from the Big Lost River and facilitated understanding and visualization of the relative importance of historical differences in infiltration in time and space. As described in a conceptual model, the numerical model simulations demonstrate flow that is (1) dominantly horizontal through interflow zones in basalt and vertical anisotropy resulting from contrasts in hydraulic conductivity of various types of basalt and the interbedded sediments, (2) temporally variable due to streamflow infiltration from the Big Lost River, and (3) moving downward downgradient of the INL. The numerical models were reparameterized, recalibrated, and analyzed to evaluate alternative conceptualizations or implementations of the conceptual model. The analysis of the reparameterized models revealed that little improvement in the model could come from alternative descriptions of sediment content, simulated aquifer thickness, streamflow infiltration, and vertical head distribution on the downgradient boundary. Of the alternative estimates of flow to or from the aquifer, only a 20 percent decrease in

Scientific Investigations Report↗

USGS2012

Simulation of groundwater and surface-water interaction and effects of pumping in a complex glacial-sediment aquifer, east central Massachusetts

The effects of groundwater pumping on surface-water features were evaluated by use of a numerical groundwater model developed for a complex glacial-sediment aquifer in northeastern Framingham, Massachusetts, and parts of surrounding towns. The aquifer is composed of sand, gravel, silt, and clay glacial-fill sediments up to 270 feet thick over an irregular fractured bedrock surface. Surface-water bodies, including Cochituate Brook, the Sudbury River, Lake Cochituate, Dudley Pond, and adjoining wetlands, are in hydraulic connection with the aquifer and can be affected by groundwater withdrawals. Groundwater and surface-water interaction was simulated with MODFLOW-NWT under current conditions and a variety of hypothetical pumping conditions. Simulations of hypothetical pumping at reactivated water supply wells indicate that captured groundwater would decrease baseflow to the Sudbury River and induce recharge from Lake Cochituate. Under constant (steady-state) pumping, induced groundwater recharge from Lake Cochituate was equal to about 32 percent of the simulated pumping rate, and flow downstream in the Sudbury River decreased at the same rate as pumping. However, surface water responded quickly to pumping stresses. When pumping was simulated for 1 month and then stopped, streamflow depletions decreased by about 80 percent within 2 months and by about 90 percent within about 4 months. The fast surface water response to groundwater pumping offers the potential to substantially reduce streamflow depletions during periods of low flow, which are of greatest concern to the ecological integrity of the river. Results indicate that streamflow depletion during September, typically the month of lowest flow, can be reduced by 29 percent by lowering the maximum pumping rates to near zero during September. Lowering pumping rates for 3 months (July through September) reduces streamflow depletion during September by 79 percent as compared to constant pumping. These results demonstrate that a seasonal or streamflow-based groundwater pumping schedule can reduce the effects of pumping during periods of low flow.

Massachusetts↗

USGS2024

Streamflow, base flow, and precipitation trends and simulated effects of Rush Springs aquifer groundwater withdrawals on base flows upgradient from Fort Cobb Reservoir, western Oklahoma

To better understand the relation between groundwater use in the Rush Springs aquifer and inflows to the Fort Cobb Reservoir, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used a previously published numerical groundwater-flow model and historical streamflow records to evaluate four scenarios to investigate how changing groundwater withdrawals could affect base flows in streams that flow into Fort Cobb Reservoir. These scenarios consisted of observing simulated base-flow response by (1) scaling the 20-year equal-proportionate-share groundwater-withdrawal rate by various percentages over a 50-year period; (2) scaling the historical groundwater-withdrawal rates by various percentages across the entire Rush Springs aquifer; (3) scaling the historical groundwater-withdrawal rates within various subareas (zones) of the Fort Cobb Reservoir surface watershed; and (4) simulating a base-flow-depletion scenario. Cobb, Lake, and Willow Creeks are the major streams upgradient from the Fort Cobb Reservoir (listed from highest to lowest mean annual base flow). The results of scenarios 1 and 2 indicated that Willow Creek is the most susceptible to drying, but Cobb Creek was the most likely to have reduced base flow. Scenarios 3 and 4 indicated that groundwater withdrawals affect Cobb Creek base flows over a broader watershed area compared to Lake and Willow Creeks. In scenario 4, Cobb Creek base-flow depletion was higher across a larger area than Lake Creek and Willow Creek. Groundwater withdrawals in the Cobb Creek watershed tended to affect total inflows into Fort Cobb Reservoir more than other areas in the extent of the Rush Springs aquifer.

Oklahoma↗

USGS2024

Use of a numerical groundwater-flow model and projected climate scenarios to simulate the effects of future climate conditions on base flow for reach 1 of the Washita River alluvial aquifer and Foss Reservoir storage, western Oklahoma

To better understand the relation between climate variability and future groundwater resources in reach 1 of the Washita River alluvial aquifer and Foss Reservoir in western Oklahoma, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used a previously published numerical groundwater-flow model and climate-model data to investigate changes in base flow and reservoir storage by evaluating three scenarios. The three projected climate scenarios were (1) a central-tendency scenario, (2) a warmer/drier scenario, and (3) a less-warm/wetter scenario. To estimate future base flow and groundwater availability in western Oklahoma, specifically in reach 1 of the Washita River alluvial aquifer, downscaled climate-model data from 231 Coupled Model Intercomparison Project phase 5 (CMIP5) projections coupled with a previously published numerical groundwater-flow model were used to compare the effects of different climate scenarios on the aquifer. Changes in base flow and groundwater-level elevations during a 30-year baseline scenario (1985–2014) and the three 30-year projected climate scenarios (2050–79) under central-tendency, warmer/drier, and less-warm/wetter climatic conditions were assessed by using the calibrated model. In the simulations, the amount of base flow and reservoir storage declined in the central-tendency and warmer/drier scenarios compared to the amount of base flow and reservoir storage under historical climatic conditions (baseline scenario). Mean annual change in reservoir storage decreased from the baseline scenario the most in the warmer/drier scenario, followed by the central-tendency scenario, but increased in the less-warm/wetter scenario compared to the baseline scenario. At the end of the simulation period (2079), the largest magnitude differences in groundwater-level elevations in all three projected climate scenarios relative to the baseline scenario occurred upstream from Foss Reservoir. Results from incorporating downscaled climate projections into localized numerical groundwater-flow models can highlight potential future changes in and implications for groundwater resources and availability.

Oklahoma↗

USGS2004

Strontium isotope geochemistry of groundwater in the central part of the Dakota (Great Plains) aquifer, USA

The Dakota aquifer of the central and eastern Great Plains of the United States is an important source of water for municipal supplies, irrigation and industrial use. Although the regional flow system can be characterized generally as east to northeasterly from the Rocky Mountains towards the Missouri River, locally the flow systems are hydrologically complex. This study uses Sr isotopic data from groundwater and leached aquifer samples to document the complex subsystems within the Dakota aquifer in Nebraska and Kansas. The interaction of groundwater with the geologic material through which it flows has created spatial patterns in the isotopic measurements that are related to: long-term water-rock interaction, during which varying degrees of isotopic equilibrium between water and rock has been achieved; and the alteration of NaCl fluids by water-rock interaction. Specifically, Sr isotopic data distinguish brines from Kansas and western Nebraska from those in eastern Nebraska: the former are interpreted to reflect interaction with Permian rocks, whereas the latter record interaction with Pennsylvanian rocks. The Sr isotopic composition of groundwater from other parts of Nebraska and Kansas are a function of the dynamic interaction between groundwater and unlithified sediments (e.g., glacial till and loess), followed by interaction with oxidized and unoxidized sediments within the Dakota Formation. This study illustrates the power of combining Sr chemistry with more conventional geochemical data to obtain a more complete understanding of groundwater flow systems within regional aquifer systems where extensive monitoring networks do not exist. ?? 2003 Elsevier Ltd. All rights reserved.

Applied Geochemistry↗

USGS2014

Simulation of groundwater flow in the Edwards-Trinity and related aquifers in the Pecos County region, Texas

The Edwards-Trinity aquifer, a major aquifer in the Pecos County region of western Texas, is a vital groundwater resource for agricultural, industrial, and public supply uses. Resource managers would like to better understand the future availability of water in the Edwards-Trinity aquifer in the Pecos County region and the effects of the possible increase or temporal redistribution of groundwater withdrawals. To that end, the U.S. Geological Survey (USGS), in cooperation with the Middle Pecos Groundwater Conservation District, Pecos County, City of Fort Stockton, Brewster County, and Pecos County Water Control and Improvement District No. 1, completed a comprehensive, integrated analysis of available hydrogeologic data to develop a groundwater-flow model of the Edwards-Trinity and related aquifers in parts of Brewster, Jeff Davis, Pecos, and Reeves Counties. Following calibration, the model was used to evaluate the sustainability of recent (2008) and projected water-use demands on groundwater resources in the study area.

Texas↗

USGS2012

Bibliography of groundwater resources of the glacial aquifer systems in Washington, Idaho, and northwestern Montana, 1905-2011

The U.S. Geological Survey Groundwater Resources Program is undertaking a series of regional groundwater availability studies to improve our understanding of groundwater availability in major aquifers across the Nation. One of the objectives of the Glacial Principal Aquifers study (proposed) is to provide information on the occurrence of groundwater in glacial aquifers in the United States, an area that includes parts of the northern continental States and much of Alaska. Toward this effort, a literature search was conducted to identify readily available documents that describe the occurrence of groundwater in glacial aquifers in the United States. This bibliography provides citations for documents, as well as codes indicating types of information available in each, for Washington, Idaho, and northwestern Montana—an area corresponding approximately to the southern extent of the Cordilleran ice sheet.

Washington;Idaho;Montana↗

USGS2012

Simulation of groundwater flow in the shallow aquifer system of the Delmarva Peninsula, Maryland and Delaware

Estimating future loadings of nitrogen to the Chesapeake Bay requires knowledge about the groundwater flow system and the traveltime of water and chemicals between recharge at the water table and the discharge to streams and directly to the bay. The Delmarva Peninsula has a relatively large proportion of its land devoted to agriculture and a large associated nitrogen load in groundwater that has the potential to enter the bay in discharging groundwater. To better understand the shallow aquifer system with respect to this loading and the traveltime to the bay, the U.S. Geological Survey constructed a steady-state groundwater flow model for the region. The model is based on estimates of recharge calculated using recently developed regression equations for evapotranspiration and surface runoff. The hydrogeologic framework incorporated into the model includes unconfined surficial aquifer sediments, as well as subcropping confined aquifers and confining beds down to 300 feet below land surface. The model was calibrated using 48 water-level measurements and 24 tracer-based ages from wells located across the peninsula. The resulting steady-state flow solution was used to estimate ages of water in the shallow aquifer system through the peninsula and the distribution and magnitude of groundwater traveltime from recharge at the water table to discharge in surface-water bodies (referred to as return time). Return times vary but are typically less than 10 years near local streams and greater than 100 years near the stream divides. The model can be used to calculate nitrate transport parameters in various local watersheds and predict future trends in nitrate loadings to Chesapeake Bay for different future nitrogen application scenarios.

Delaware;Maryl↗

USGS2011

Estimated 2008 groundwater potentiometric surface and predevelopment to 2008 water-level change in the Santa Fe Group aquifer system in the Albuquerque area, central New Mexico

The water-supply requirements of the Albuquerque metropolitan area of central New Mexico have historically been met almost exclusively by groundwater withdrawal from the Santa Fe Group aquifer system. Previous studies have indicated that the large quantity of groundwater withdrawal relative to recharge has resulted in water-level declines in the aquifer system throughout the metropolitan area. Analysis of the magnitude and pattern of water-level change can help improve understanding of how the groundwater system responds to withdrawals and variations in the management of the water supply and can support water-management agencies' efforts to minimize future water-level declines and improve sustainability. This report, prepared by the U.S. Geological Survey in cooperation with the Albuquerque Bernalillo County Water Utility Authority, presents the estimated groundwater potentiometric surface during winter (from December to March) of the 2008 water year and the estimated changes in water levels between predevelopment and water year 2008 for the production zone of the Santa Fe Group aquifer system in the Albuquerque and surrounding metropolitan and military areas. Hydrographs from selected wells are included to provide details of historical water-level changes. In general, water-level measurements used for this report were measured in small-diameter observation wells screened over short intervals and were considered to best represent the potentiometric head in the production zone-the interval of the aquifer, about 300 feet below land surface to 1,100 feet or more below land surface, in which production wells generally are screened. Water-level measurements were collected by various local and Federal agencies. The 2008 water year potentiometric surface map was created in a geographic information system, and the change in water-level elevation from predevelopment to water year 2008 was calculated. The 2008 water-level contours indicate that the general direction of groundwater flow is from the Rio Grande towards clusters of production wells in the east, north, and west. Water-level changes from predevelopment to 2008 are variable across the area. Hydrographs from piezometers on the east side of the river generally indicate a trend of decline in the annual highest water level through most of the period of record. Hydrographs from piezometers in the valley near the river and on the west side of the river indicate spatial variability in water-level trends.

New Mexico↗

USGS2014

Estimated 2012 groundwater potentiometric surface and drawdown from predevelopment to 2012 in the Santa Fe Group aquifer system in the Albuquerque metropolitan area, central New Mexico

Historically, the water-supply requirements of the Albuquerque metropolitan area of central New Mexico were met almost exclusively by groundwater withdrawal from the Santa Fe Group aquifer system. In response to water-level declines, the Albuquerque Bernalillo County Water Utility Authority (ABCWUA) began diverting water from the San Juan-Chama Drinking Water Project in December 2008 to reduce the use of groundwater to meet municipal demand. Modifications in the demand for water and the source of the supply of water for the Albuquerque metropolitan area have resulted in a variable response in the potentiometric surface of the production zone (the interval of the aquifer, from within about 200 feet below the water table to 900 feet or more, in which supply wells generally are screened) of the Santa Fe Group aquifer system. Analysis of the magnitude and spatial distribution of water-level change can help improve the understanding of how the groundwater system responds to withdrawals and variations in the management of the water supply and can support water-management agencies’ efforts to minimize future water-level declines and improve sustainability. The U.S. Geological Survey (USGS), in cooperation with the ABCWUA, has developed an estimate of the 2012 potentiometric surface of the production zone of the Santa Fe Group aquifer system in the Albuquerque metropolitan area. This potentiometric surface is the latest in a series of reports depicting the potentiometric surface of the area. This report presents the estimated potentiometric surface during winter (from December to March) of water year 2012 and the estimated changes in potentiometric surface between predevelopment (pre-1961) and water year 2012 for the production zone of the Santa Fe Group aquifer system in the Albuquerque metropolitan area. Hydrographs from selected piezometers are included to provide details of historical water-level changes. In general, water-level measurements used for this report were collected in small-diameter observation wells screened over short intervals near the middle of the production zone and were considered to best represent the potentiometric head in the production zone. The water-level measurements were collected by various local and Federal agencies. The water year 2012 potentiometric surface map was created in a geographic information system, and the change in water-level altitude from predevelopment to water year 2012 was calculated. The 2012 potentiometric surface indicates that the general direction of groundwater flow is from the Rio Grande towards clusters of supply wells in the east, north, and west. Water-level changes from predevelopment to 2012 were variable across the Albuquerque metropolitan area. Estimated drawdown from 2008 was spatially variable across the Albuquerque metropolitan area. Hydrographs from piezometers on the east side of the river indicate an increase in the annual highest water-level measurement from 2008 to 2012. Hydrographs from piezometers in the northwest part of the study area indicate either steady decline of the water-level altitude over the period of record or recently variable trends in which water-level altitudes increased for a number of years but have declined since water year 2012.

New Mexico↗

USGS2017

An update of hydrologic conditions and distribution of selected constituents in water, eastern Snake River Plain aquifer and perched groundwater zones, Idaho National Laboratory, Idaho, emphasis 2012-15

Since 1952, wastewater discharged to in ltration ponds (also called percolation ponds) and disposal wells at the Idaho National Laboratory (INL) has affected water quality in the eastern Snake River Plain (ESRP) aquifer and perched groundwater zones underlying the INL. The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy, maintains groundwater-monitoring networks at the INL to determine hydrologic trends and to delineate the movement of radiochemical and chemical wastes in the aquifer and in perched groundwater zones. This report presents an analysis of water-level and water-quality data collected from the ESRP aquifer, multilevel monitoring system (MLMS) wells in the ESRP aquifer, and perched groundwater wells in the USGS groundwater monitoring networks during 2012-15.

Idaho↗

USGS2025

Estimated hydrogeologic, spatial, and temporal distribution of self-supplied domestic groundwater withdrawals for aquifers of the Virginia Coastal Plain

Water use from private-domestic wells accounts for nearly 40 percent of total groundwater withdrawals in the Virginia Coastal Plain Physiographic Province (henceforth called the Virginia Coastal Plain). However, because self-supplied domestic water use generally falls below the Virginia Department of Environmental Quality (VDEQ) reporting and management threshold of 300,000 gallons per month, quantifying these withdrawals is challenging. This report builds upon the foundation of previous U.S. Geological Survey investigations by providing revised techniques to improve estimates of the aquifer source, spatial distribution, and monthly magnitude of these groundwater withdrawals. The aquifer sources of private-domestic wells in the Virginia Coastal Plain were estimated by cross-referencing 8,264 well records from the VDEQ and the Virginia Department of Health to a digital model of the Virginia Coastal Plain hydrogeologic framework. This analysis highlights the regional importance of the Yorktown-Eastover, Potomac, and surficial aquifers. Collectively, these three aquifers account for 80 percent of self-supplied domestic groundwater withdrawals. The population using self-supplied domestic water was estimated using census blocks, well-use ratios, building footprints, and land-use and land-cover data to produce a high-resolution, disaggregated, raster-based dataset. This approach improves upon previous models at the census-block or road-network scale by reducing the low-density spread of the self-supplied domestic population across undeveloped areas and concentrating the population and its corresponding water use in the areas where it is most likely to occur. Results show that an estimated 475,332 people comprise the 2020 self-supplied domestic population of the Virginia Coastal Plain, an increase of 5.7 percent since 2010, and the greatest concentrations of self-supplied domestic population surround large cities. Estimates could be further refined with the addition of current and complete spatial data on public water-system service areas. The quantity of water used by the self-supplied domestic population was estimated by modifying published state per-capita water-use coefficients with the corresponding monthly variability assessed from Virginia Coastal Plain public water-system withdrawal data. This analysis estimates an average increase of 12 percent from June through August and an average decrease of 8 percent from December through March from the baseline annual average of 80 gallons per day per capita, which generally matches similar studies in the eastern United States. The application of these revised methodologies for the estimation of private-domestic wells and the self-supplied domestic population improves understanding of domestic groundwater use in the Virginia Coastal Plain across hydrogeologic, spatial, and temporal scales. These revisions help better inform water-resource managers and decision makers and support higher resolution groundwater modeling. Furthermore, these methods are transferrable to other areas where self-supplied domestic water withdrawals are important to the overall water budget.

Virginia↗

USGS2020

Comparison of groundwater storage changes from GRACE satellites with monitoring and modeling of major U.S. aquifers

GRACE satellite data are widely used to estimate groundwater storage (GWS) changes in aquifers globally; however, comparisons with GW monitoring and modeling data are limited. Here we compared GWS changes from GRACE over 15 yr (2002–2017) in 14 major U.S. aquifers with groundwater-level (GWL) monitoring data in ~23,000 wells and with regional and global hydrologic and land surface models. Results show declining GWS trends from GRACE data in the six southwestern and south-central U.S. aquifers, totaling −90 km 3 over 15 yr, related to long-term (5–15 yr) droughts, and exceeding Lake Mead volume by ~2.5×. GWS trends in most remaining aquifers were stable or slightly rising. GRACE-derived GWS changes agree with GWL monitoring data in most aquifers (correlation coefficients, R = 0.52–0.95), showing that GRACE satellites capture groundwater (GW) dynamics. Regional GW models (eight models) generally show similar or greater GWS trends than those from GRACE. Large discrepancies in the Mississippi Embayment aquifer, with modeled GWS decline approximately four times that of GRACE, may reflect uncertainties in model storage parameters, stream capture, pumpage, and/or recharge rates. Global hydrologic models (2003–2014), which include GW pumping, generally overestimate GRACE GWS depletion (total: approximately −172 to −186 km 3 ) in heavily exploited aquifers in southwestern and south-central U.S. by ~2.4× (GRACE: −74 km 3 ), underscoring needed modeling improvements relative to anthropogenic impacts. Global land surface models tend to track GRACE GWS dynamics better than global hydrologic models. Intercomparing remote sensing, monitoring, and modeling data underscores the importance of considering all data sources to constrain GWS uncertainties.

Water Resources Research↗

USGS2020

Water-quality comparison of the Gulf Coast aquifer system at various scales in Texas from National Water-Quality Assessment groundwater studies, 2013–15

One of the objectives of the U.S. Geological Survey National Water-Quality Assessment (NAWQA) Project is to assess groundwater quality in aquifers that are important sources of drinking water such as the coastal lowlands aquifer system, which is often referred to in Texas as the “Gulf Coast aquifer system.” The Gulf Coast aquifer system extends from Louisiana to Mexico and is a source of groundwater for several cities including Houston, Tex. The NAWQA groundwater studies in Texas in 2013–15 that assessed the Gulf Coast aquifer system included Principal Aquifer Surveys (PAS), Major Aquifer Studies (MAS), and Land Use Studies (LUS). These three study types are based on sampling networks of wells distributed in an area of interest. The PAS networks typically consist of public-supply wells that are relatively deep, the MAS networks typically consist of domestic-supply wells that are intermediate in depth, and the LUS networks typically consist of monitoring wells that are relatively shallow.

Texas↗

USGS2018

Hydrogeology of, simulation of groundwater flow in, and potential effects of sea-level rise on the Kirkwood-Cohansey aquifer system in the vicinity of Edwin B. Forsythe National Wildlife Refuge, New Jersey

The Edwin B. Forsythe National Wildlife Refuge encompasses more than 47,000 acres of New Jersey coastal habitats, including salt marshes, freshwater wetlands, tidal wetlands, barrier beaches, woodlands, and swamps. The refuge is along the Atlantic Flyway and provides breeding habitat for fish, migratory birds, and other wildlife species. The refuge area may be threatened by global climate change, including sea-level rise (SLR). The Kirkwood-Cohansey aquifer system underlies the Edwin B. Forsythe National Wildlife Refuge. Groundwater is an important source of freshwater flow into the refuge, but information about the interaction of surface water and groundwater in the refuge area and the potential effects of SLR on the underlying aquifer system is limited. The U.S. Geological Survey (USGS), in cooperation with the U.S. Fish and Wildlife Service (USFWS), conducted a hydrologic assessment of the refuge in New Jersey and developed a groundwater flow model to improve understanding of the geohydrology of the refuge area and to serve as a tool to evaluate changes in groundwater-level altitudes that may result from a rise in sea level. Groundwater flow simulations completed for this study include a calibrated baseline simulation that represents 2005–15 hydraulic conditions and three SLR scenarios―20, 40, and 60 centimeters (cm) (0.656, 1.312, and 1.968 feet, respectively). Results of the three SLR simulations indicate that the water table in the unconfined Kirkwood-Cohansey aquifer system in the refuge area will rise, resulting in increased discharge of fresh groundwater to freshwater wetlands and streams. As sea level rises, simulated groundwater discharge to the salt marsh, bay, and ocean is projected to decrease. Flow from the salt marsh, bay, and ocean to the overlying surface water is projected to increase as sea level rises. The simulated movement of the freshwater-seawater interface as sea level rises depends on the hydraulic-head gradient. In the center of the Forsythe model area, topographic relief is 23 feet (ft) and the hydraulic-head gradient is 0.0033. In the center of the Forsythe model area, the simulated interface moved inland about 600 ft and downward about 15 ft from the baseline simulation to scenario 3 as a result of a SLR of 60 cm. In the southern part of the Forsythe model area, the topography is flatter (relief of 8 ft) and the hydraulic-head gradient is smaller (0.001). In the southern part of the Forsythe model study area, the simulated interface in this area is projected to move inland about 200 ft from the baseline simulation to scenario 3 and does not move downward.

New Jersey↗

USGS2025

Groundwater response to managed aquifer recharge at the Southeast Houghton Artificial Recharge Project in Tucson, Arizona

Managed aquifer recharge is a widespread practice for storing water in the subsurface as groundwater. At a managed aquifer recharge facility in southern Arizona, groundwater-level and repeat microgravity data were collected to monitor aquifer response. These data were used to inform parameter identification for an unsaturated-zone flow model used to simulate the recharge process. The facility, the Southeast Houghton Artificial Recharge Project (SHARP), consists of 3 surface basins (about 27,600 square meters [6.8 acres] total surface area) where recycled water is distributed in recharge cycles lasting several months, with dry periods in between. During the study period, December 2020–December 2022, Tucson Water (the City of Tucson’s water utility) reported 6.56×10 6 cubic meters of water (5,320 acre-feet) recharged. Monitoring included groundwater-level observations at 3 monitoring wells and repeat microgravity measurements at as many as 22 locations (some stations were destroyed between surveys). Six gravity surveys were carried out using absolute- and relative-gravity meters. Large gravity increases, more than 250 microgals, were observed during the first repeat survey, 3.5 months after the start of recharge, but only in the immediate vicinity of the recharge basins. Data show that water moved downward to the water table, and storage changes in the unsaturated zone away from the facility were likely minimal. Gravity decreased at stations more than 1 kilometer from the facility, consistent with regional groundwater-level changes. Groundwater-level increases in wells adjacent to the recharge basins began 2 months after the second repeat gravity survey, and 5.5 months after recharge began. Unsaturated-zone flow modeling was carried out using software that simulates water movement and parameter estimation. Model calibration was carried out by minimizing an objective function calculated from the differences between simulated and observed groundwater levels, and between simulated and observed repeat microgravity data. Including repeat microgravity data in the objective function reduced the uncertainty in estimated parameter values for saturated hydraulic conductivity and saturated water content. Modeling indicated that the unsaturated zone between the recharge basins and the water table does not become saturated even after 685 days of simulated infiltration. This gradual wetting may account for increasing infiltration rates over time, as hydraulic conductivity increases with increasing water content. Unsaturated-zone water content decreased rapidly between recharge cycles. Model-simulated groundwater mounding extended about 1 kilometer from the center of SHARP after the 685-day period following the onset of recharge.

Arizona↗

USGS2023

Per- and polyfluoroalkyl substances in groundwater from the Great Miami buried-valley aquifer, southwestern Ohio, 2019–20

Groundwater samples were collected during 2019 and 2020 from 23 wells in the Great Miami buried-valley aquifer (GM-BVA) in southwestern Ohio by the U.S. Geological Survey, in cooperation with the Miami Conservancy District, Dayton, Ohio, to determine concentrations of selected per- and polyfluoroalkyl substances (PFAS). The GM-BVA is a glacial outwash and alluvial fill aquifer that is the sole source of water supply for much of the region. Wells had total depths that ranged from 21 to 101 feet below land surface, and groundwater levels that ranged from 1.39 to 52.15 feet below land surface before sampling in 2019. Groundwater and related quality-control samples were sequentially collected from 22 of the 23 wells and analyzed for 24 different PFAS by 2 methods that used proprietary isotope-dilution based adaptations of U.S. Environmental Protection Agency (EPA) method 537.1, termed methods 1 and 2. Method 2 had smaller reporting limits (RL) for 22 of 24 PFAS analyzed and smaller detection limits (DLs) for all 24 PFAS analyzed compared with method 1, which made method 2 more sensitive to detect PFAS. Concentrations of perfluorooctanesulfonate (PFOS) in a groundwater (GW)-method 2 sample from well CL–275 of 1.9 nanograms per liter (ng/L) and perfluorooctanoate (PFOA) in a GW-method 2 sample from well BU–1106 of 2.1 ng/L were greater than their EPA interim health advisory guidances for drinking water (as of June 2022) by about 9,500 and 52,500 percent, respectively. The EPA interim health advisory guidances for PFOS (0.02 ng/L) and PFOA (0.004 ng/L) were also 65 and 215 times less, respectively, than the smallest method 2 DLs for PFOS (1.3 ng/L) and PFOA (0.86 ng/L). Other PFAS were either not detected in GM-BVA groundwater samples or were detected in concentrations less than Ohio action levels or Federal health-risk-based guidance. The most detected PFAS in groundwater was perfluorobutanesulfonate (PFBS), which had concentrations in samples from eight wells that ranged from 1.0 to 8.0 ng/L or from 0.05 to 0.4 percent of its EPA health advisory of 2,000 ng/L for drinking water. The similarity of PFBS (7.8 ng/L), perfluoropentanesulfonate (PFPeS; 8.1 ng/L), and perfluorohexanesulfonate (PFHxS; 14 ng/L) concentrations yielded from the GW-method 1 sample from well CL–275 on July 9, 2019, to those of PFBS (8.0 ng/L), PFPeS (7.8 ng/L), and PFHxS (16 ng/L) from the paired GW-method 2 sample demonstrated the capability of both methods to reproduce PFAS concentrations that were greater than their respective DLs. Non-detection of these PFAS in follow-up GW-method 1 and sequential replicate (Rep–GW-method 1) samples from CL–275 on April 21, 2020, indicated that the 2019 results represented a transient detection in groundwater. Eleven of twenty-three wells sampled in 2019 had from 1 to 4 PFAS detected in one or more groundwater samples or in a paired replicate sample: PFBS in 8 wells and 9 samples; PFHxS in 4 wells and 5 samples; and PFPeS, PFOS, perfluorobutanoate, perfluoropentanoate, PFOA, and perfluorooctanesulfonamide in 1 well and 1 sample each. More PFAS were detected in GW-method 2 samples than GW-method 1 samples because method 2 had smaller RLs and DLs. Results indicate benefits from the analysis of paired samples, sequential replicate samples, and other quality-control samples using analytical methods with sensitive RLs and DLs to verify PFAS concentrations in groundwater. Groundwater-age estimates indicate that water produced from all sampled wells had infiltrated to the water table within the 1947–present (2022) period of PFAS use or environmental presence. Eight wells with detectable PFBS in groundwater from 2019 samples also had groundwater-recharge dates that ranged from 1991 to 2016. Those ages coincided with the possible environmental presence of PFBS as a PFAS byproduct or use as an alternative to PFOS after about 2002. Two wells that had detections of PFHxS in 2019 groundwater samples also had post-2000 groundwater-recharge dates that coincided with the period of use of PFHxS as an alternative to PFOS. Six of nine wells with more than 66-percent of urban land use that was within 0.3 miles of each well, as of 2012, also had 1 to 4 PFAS detected in one of their groundwater samples. Seven of nine wells that produced groundwater in 2019 with an oxic redox category also had one or more PFAS detected in a sample.

Ohio↗