USGS Science⌕ Search

SEARCH · USGS Science

Results for “groundwater aquifer”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 433 records · Page 24Linked to original sources
USGS2019

Potentiometric surface of the Mississippi River Valley alluvial aquifer, spring 2016

A potentiometric surface map for spring 2016 was created for the Mississippi River Valley alluvial (MRVA) aquifer using selected available groundwater-altitude data from wells and surface-water-altitude data from streamgages. Most of the wells were measured annually or one time after installation, but some wells were measured more than one time or continually; streamgages are typically operated continuously. Personnel from the Arkansas Natural Resources Commission, Arkansas Department of Health, Arkansas Geological Survey, Illinois Department of Agriculture, Illinois State Water Survey, Louisiana Department of Natural Resources, Louisiana Department of Transportation and Development, Mississippi Department of Environmental Quality, Yazoo Mississippi Delta Joint Water Management District, U.S. Department of Agriculture–Natural Resources Conservation Service, and the U.S. Geological Survey (USGS) routinely collect groundwater data from wells screened in the MRVA aquifer. The USGS and the U.S. Army Corps of Engineers routinely collect data on river stage and discharge for the rivers overlying the MRVA aquifer. The potentiometric surface map for 2016 was created using existing data as part of the USGS Water Availability and Use Science Program to support investigations that characterize the MRVA aquifer. Sufficient groundwater-altitude data were available to create a potentiometric-surface map for spring 2016 for about 81 percent of the aquifer area. The potentiometric contours ranged from 10 to 340 feet. The regional direction of groundwater flow in the MRVA aquifer was generally towards the south-southwest, except in areas of groundwater-altitude depressions, where groundwater flows into the depressions, and near rivers, where groundwater flow generally parallels the flow in the rivers. There are large depressions in the potentiometric surface of the MRVA aquifer in the lower half of the Cache region and in most of the Grand Prairie and Delta regions.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗
USGS2016

Application of hydrogeology and groundwater-age estimates to assess the travel time of groundwater at the site of a landfill to the Mahomet Aquifer, near Clinton, Illinois

The U.S. Geological Survey used interpretations of hydrogeologic conditions and tritium-based groundwater age estimates to assess the travel time of groundwater at a landfill site near Clinton, Illinois (the “Clinton site”) where a chemical waste unit (CWU) was proposed to be within the Clinton landfill unit #3 (CLU#3). Glacial deposits beneath the CWU consist predominantly of low-permeability silt- and clay-rich till interspersed with thin (typically less than 2 feet in thickness) layers of more permeable deposits, including the Upper and Lower Radnor Till Sands and the Organic Soil unit. These glacial deposits are about 170 feet thick and overlie the Mahomet Sand Member of the Banner Formation. The Mahomet aquifer is composed of the Mahomet Sand Member and is used for water supply in much of east-central Illinois. Eight tritium analyses of water from seven wells were used to evaluate the overall age of recharge to aquifers beneath the Clinton site. Groundwater samples were collected from six monitoring wells on or adjacent to the CLU#3 that were open to glacial deposits above the Mahomet aquifer (the upper and lower parts of the Radnor Till Member and the Organic Soil unit) and one proximal production well (approximately 0.5 miles from the CLU#3) that is screened in the Mahomet aquifer. The tritium-based age estimates were computed with a simplifying, piston-flow assumption: that groundwater moves in discrete packets to the sampled interval by advection, without hydrodynamic dispersion or mixing. Tritium concentrations indicate a recharge age of at least 59 years (pre-1953 recharge) for water sampled from deposits below the upper part of the Radnor Till Member at the CLU#3, with older water expected at progressively greater depth in the tills. The largest tritium concentration from a well sampled by this study (well G53S; 0.32 ± 0.10 tritium units) was in groundwater from a sand deposit in the upper part of the Radnor Till Member; the shallowest permeable unit sampled by this study. That result indicated that nearly all groundwater sampled from well G53S entered the aquifer as recharge before 1953. Tritium was detected in a trace concentration in one sample from a second monitoring well open to the upper part of the Radnor Till Member (well G07S; 0.11 ± 0.09 tritium units), and not detected in samples collected from two monitoring wells open to a sand deposit in the lower part of the Radnor Till Member, from two samples collected from two monitoring wells open to the Organic Soil unit, and in two samples collected from a production well screened in the middle of the Mahomet aquifer (a groundwater sample and a sequential replicate sample). The lack of tritium in five of the six groundwater samples collected from the shallow permeable units beneath CLU#3 site and the two samples from the one Mahomet aquifer well indicates an absence of post-1952 recharge. Groundwater-flow paths that could contribute post-1952 recharge to the lower part of the Radnor Till Member, the Organic Soil unit, or the Mahomet aquifer at the CLU#3 are not indicated by these data. Hypothetical two-part mixtures of tritium-dead, pre-1953 recharge water and decay-corrected tritium concentrations in post-1952 recharge were computed and compared with tritium analyses in groundwater sampled from monitoring wells at the CLU#3 site to evaluate whether tritium concentrations in groundwater could be represented by mixtures involving some post-1952 recharge. Results from the hypothetical two-part mixtures indicate that groundwater from monitoring well (G53S) was predominantly composed of pre-1953 recharge and that if present, younger, post-1955 recharge, contributed less than 2.5 percent to that sample. The hypothetical two-part mixing results also indicated that very small amounts of post-1952 recharge composing less than about 2.5 percent of the sample volume could not be distinguished in groundwater samples with tritium concentrations less than about 0.15 TU. The piston-flow based age of recharge determined from the tritium concentration in the groundwater sample from monitoring well G53S yielded an estimated maximum vertical velocity from the land surface to the upper part of the Radnor Till Member of 0.85 feet per year or less. This velocity, ifassumed to apply to the remaining glacial till deposits above the Mahomet aquifer, indicates that recharge flows through the 170 feet of glacial deposits between the base of the proposed chemical waste unit and the top of the Mahomet aquifer in a minimum of 200 years or longer. Analysis of hydraulic data from the site, constrained by a tritium-age based maximum groundwater velocity estimate, computed minimum estimates of effective porosity that range from about 0.021 to 0.024 for the predominantly till deposits above the Mahomet aquifer. Estimated rates of transport of recharge from land surface to the Mahomet aquifer for the CLU#3 site computed using the Darcy velocity equation with site-specific data were about 260 years or longer. The Darcy velocity-based estimates were computed using values that were based on tritium data, estimates of vertical velocity and effective porosity and available site-specific data. Solution of the Darcy velocity equation indicated that maximum vertical groundwater velocities through the deposits above the aquifer were 0.41 or 0.61 feet per year, depending on the site-specific values of vertical hydraulic conductivity (laboratory triaxial test values) and effective porosity used for the computation. The resulting calculated minimum travel times for groundwater to flow from the top of the Berry Clay Member (at the base of the proposed chemical waste unit) to the top of the Mahomet aquifer ranged from about 260 to 370 years, depending on the velocity value used in the calculation. In comparison, plausible travel times calculated using vertical hydraulic conductivity values from a previously published regional groundwater flow model were either slightly less than or longer than those calculated using site data and ranged from 230 to 580 years. Tritium data from 1996 to 2011 USGS regional sampling of groundwater from domestic wells in the confined part of the Mahomet aquifer—which are 2.5 to about 40 miles from the Clinton site—were compared with site-specific data from a production well at the Clinton site. Tritium-based groundwater-age estimates indicated predominantly pre- 1953 recharge dates for USGS and other prior regional samples of groundwater from domestic wells in the Mahomet aquifer. These results agreed with the tritium-based, pre-1953 recharge age estimated for a groundwater sample and a sequential replicate sample from a production well in the confined part of the Mahomet aquifer beneath the Clinton site. The regional tritium-based groundwater age estimates also were compared with pesticide detections in samples from distal domestic wells in the USGS regional network that are about 2.5 to 40 miles from the Clinton site to identify whether very small amounts of post-1952 recharge have in places reached confined parts of the Mahomet aquifer at locations other than the Clinton site in an approximately 2,000 square mile area of the Mahomet aquifer. Very small amounts of post-1952 recharge were defined in this analysis as less than about 2.5 percent of the total recharge contributing to a groundwater sample, based on results from the two-part mixing analysis of tritium data from the Clinton site. Pesticide-based groundwater-age estimates based on 22 detections of pesticides (13 of these detections were estimated concentrations), including atrazine, deethylatrazine (2-Chloro-4-isopropylamino-6-amino- s-triazine), cyanazine, diazinon, metolachlor, molinate, prometon, and trifluralin in groundwater samples from 10 domestic wells 2.5 to about 40 miles distant from the Clinton site indicate that very small amounts of post-1956 to post-1992 recharge can in places reach the confined part of the Mahomet aquifer in other parts of central Illinois. The relative lack of tritium in these samples indicate that the amounts of post-1956 to post-1992 recharge contributing to the 10 domestic wells were a very small part of the overall older groundwater sampled from those wells. The flow process by which very small amounts of pesticide-bearing groundwater reached the screened intervals of the 10 domestic wells could not be distinguished between well-integrity related infiltration and natural hydrogeologic features. Potential explanations include: (1) infiltration through man-made avenues in or along the well, (2) flow of very small amounts of post-1956 to post-1992 recharge through sparsely distributed natural permeable aspects of the glacial till and diluted by mixing with older groundwater, or (3) a combination of both processes. Presuming the domestic wells sampled by the USGS in 1996–2011 in the regional study of the confined part of the Mahomet aquifer are adequately sealed and produce groundwater that is representative of aquifer conditions, the regional tritium and pesticide-based groundwater-age results indicate substantial heterogeneity in the glacial stratigraphy above the Mahomet aquifer. The pesticide-based groundwater-age estimates from the domestic wells distant from the Clinton site also indicate that parts of the Mahomet aquifer with the pesticide detections can be susceptible to contaminant sources at the land surface. The regional pesticide and tritium results from the domestic wells further indicate that a potential exists for possible contaminants from land surface to be transported through the glacial drift deposits that confine the Mahomet aquifer in other parts of central Illinois at faster rates than those computed for recharge at the Clinton site, including CLU#3. This analysis indicates the potential value of sub-microgram-per-liter level concentrations of land-use derived indicators of modern recharge to indicate the presence of very small amounts of modern, post-1952 age recharge in overall older, pre-1953 age groundwater.

Illinois↗
USGS2019

Groundwater-flow model and analysis of groundwater and surface-water interactions for the Big Sioux aquifer, Sioux Falls, South Dakota

The city of Sioux Falls, in southeastern South Dakota, is the largest city in South Dakota. The U.S. Geological Survey (USGS), in cooperation with the city of Sioux Falls, completed a groundwater-flow model to use for improving the understanding of groundwater-flow processes, estimating hydrogeologic properties, and analyzing groundwater and surface-water interactions for the Big Sioux aquifer in the model area. The model area includes the Big Sioux aquifer and the underlying hydrogeologic units from Dell Rapids, South Dakota, to the confluence of the Big Sioux River and the outlet of the Sioux Falls Diversion Channel in eastern Sioux Falls, S. Dak. The Big Sioux aquifer is the primary aquifer in the model area and the focus of the groundwater-flow model. The Big Sioux River is the largest stream in the model area and is in hydraulic connection with the Big Sioux aquifer. A conceptual model for the area was constructed and includes a characterization of the hydrogeologic framework, analysis and construction of potentiometric surfaces, and summary of estimated water budget components in the model area. The primary hydrogeologic units in the model area consist of (1) the Big Sioux aquifer, (2) a glacial till confining unit, and (3) bedrock aquifers (Split Rock Creek and Sioux Quartzite aquifers). Sources of groundwater recharge included infiltration of precipitation, stream seepage, and groundwater exchanges among the hydraulically connected Big Sioux aquifer, glacial till confining unit, and bedrock aquifers. Groundwater losses included evapotranspiration, groundwater discharge to streams, and groundwater withdrawal to supply water-use needs. A numerical groundwater-flow model (numerical model) was constructed and was used to simulate all aspects of the conceptual model for predevelopment (steady-state) and time-varying (transient) monthly conditions for 1950–2017. The numerical model was constructed using the USGS modular hydrologic simulation program, MODFLOW–6, and was calibrated using the Parameter ESTimation software, PEST++. The transient numerical model was calibrated for steady-state and transient monthly conditions for 1950–2017. Calibration targets were observations of hydraulic head, changes in hydraulic head, monthly mean streamflow (as a rate), and cumulative monthly stream discharge (as a volume). Parameters adjusted during model calibration were horizontal and vertical hydraulic conductivity, specific storage, specific yield, recharge and evapotranspiration multipliers, and streambed hydraulic conductivity. Horizontal and vertical hydraulic conductivity were estimated at pilot points distributed within the model area; specific storage and specific yield were assigned to uniform values in each layer in the model area; recharge and evapotranspiration multipliers were assigned uniformly for every stress period in the numerical model; and streambed hydraulic conductivity values were assigned uniformly between stream confluences. The final calibrated parameter values of horizontal and vertical hydraulic conductivity, specific yield, specific storage, streambed hydraulic conductivity, recharge, and evapotranspiration were considered reasonable for the hydrogeologic materials and conditions in the model area for 1950–2017. Overall, simulated hydraulic head altitudes had a linear regression coefficient of determination (R 2 ) of 0.48. Hydraulic head altitude residuals for the glacial till confining unit and bedrock aquifers were typically greater in magnitude when compared to residuals in the Big Sioux aquifer, but simulated hydraulic head altitudes in the Big Sioux aquifer compared favorably with mean observed hydraulic head altitudes and had a linear regression R 2 of 0.93. Simulated streamflow hydrographs matched the general trends of observed increases and decreases in streamflow for USGS streamgages 06482000 (Big Sioux River at Sioux Falls, S. Dak.) and 06482020 (Big Sioux River at North Cliff Avenue at Sioux Falls, S. Dak.), but larger streamflows were overestimated at the first streamgage and underestimated at the second streamgage. The numerical model reasonably estimated cumulative monthly stream discharge for the first 10–15 years of available streamflow records at both USGS streamgages. After the first 10–15 years of available streamflow record, cumulative monthly stream discharge was closely estimated for USGS streamgage 06482000 and underestimated at USGS streamgage 06482020. Composite sensitivities without regularization were calculated by PEST++ for the calibrated numerical model parameters and were averaged by parameter group. The parameter group with the highest mean composite sensitivity was the recharge multiplier parameter group. Model simplifications, assumptions, and limitations were necessary for construction of the conceptual and numerical models and for calibration efficiency. Spatial simplification of hydraulic properties could cause the numerical model to misrepresent reactions to changes in localized stresses, such as additional demands for groundwater withdrawal. The numerical model was temporally discretized into monthly periods and required scaling daily rates into representative monthly rates for model input and calibration targets. Based on the comparison between the observed and simulated groundwater levels, monthly mean streamflow and cumulative monthly stream discharge, and general groundwater distribution and flow, the numerical model favorably simulated the flow in the Big Sioux aquifer. Eventual capture was calculated in the model area using a steady-state numerical groundwater-flow model. The eventual capture map shows areas of higher streamflow capture adjacent to the Big Sioux River north of the city of Sioux Falls and along the lower part of the Sioux Falls Diversion Channel, and areas of lower streamflow capture along aquifer boundaries and near the southern Sioux Quartzite barrier. The timing of capture was determined using a transient numerical groundwater-flow model to determine the likely captured water sources for 30 years of groundwater withdrawal at three hypothetical wells using three continuous withdrawal rates (112.5, 450.0, and 900.0 gallons per minute). Supply for all three hypothetical wells became capture-dominated after only a short period of continuous withdrawal. Capture stabilized after about 10–15 years for well A, and after 20–25 years for well B, and after about 10–15 years for well C. The groundwater-flow model is a suitable tool to use for improving the understanding of groundwater-flow processes, estimating hydrogeologic properties, and analyzing groundwater and surface-water interactions for the Big Sioux aquifer near Sioux Falls, S. Dak. The numerical model can be used to simulate hydrologic scenarios, advance understanding of groundwater budgets, compute system response to stress, and determine likely sources of water supplied to wells.

South Dakota↗
USGS2018

Comparing groundwater quality in public-supply and shallow aquifers in the Monterey Bay and Salinas Valley Basins, California

Groundwater provides more than 40 percent of California’s drinking water. To protect this vital resource, the State of California created the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The Priority Basin Project of the GAMA Program (GAMA-PBP) provides a comprehensive assessment of the State’s groundwater quality and increases public access to groundwater-quality information.

California↗
USGS2014

Hydrogeology and simulation of groundwater flow in the Central Oklahoma (Garber-Wellington) Aquifer, Oklahoma, 1987 to 2009, and simulation of available water in storage, 2010–2059

The Central Oklahoma (Garber-Wellington) aquifer underlies about 3,000 square miles of central Oklahoma. The study area for this investigation was the extent of the Central Oklahoma aquifer. Water from the Central Oklahoma aquifer is used for public, industrial, commercial, agricultural, and domestic supply. With the exception of Oklahoma City, all of the major communities in central Oklahoma rely either solely or partly on groundwater from this aquifer. The Oklahoma City metropolitan area, incorporating parts of Canadian, Cleveland, Grady, Lincoln, Logan, McClain, and Oklahoma Counties, has a population of approximately 1.2 million people. As areas are developed for groundwater supply, increased groundwater withdrawals may result in decreases in long-term aquifer storage. The U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, investigated the hydrogeology and simulated groundwater flow in the aquifer using a numerical groundwater-flow model. The purpose of this report is to describe an investigation of the Central Oklahoma aquifer that included analyses of the hydrogeology, hydrogeologic framework of the aquifer, and construction of a numerical groundwater-flow model. The groundwater-flow model was used to simulate groundwater levels and for water-budget analysis. A calibrated transient model was used to evaluate changes in groundwater storage associated with increased future water demands.

Oklahoma↗
USGS2022

Groundwater budgets for the Big Lost River Basin, south-central Idaho, 2000–19

The Big Lost River Basin, located in parts of Butte and Custer Counties in south-central Idaho, supports the communities surrounding the cities of Arco, Leslie, Mackay, and Moore and provides for agricultural resources that depend on a sustainable supply of surface water from the Big Lost River and its tributaries and groundwater from an unconfined aquifer. The aquifer, situated in a structurally controlled intermontane valley, is composed of unconsolidated alluvium, consolidated sedimentary and volcanic rocks, and younger interbedded volcanic rocks. This report presents two separate groundwater budgets for the aquifer, one above and one below Mackay Dam, as well as a combined groundwater budget for the aquifer within the entire Big Lost River Basin. The budgets span a 20-year period (2000–19), characterizing average conditions, a dry year (2014), and a wet year (2017). The groundwater budgets will help address questions regarding the availability of groundwater supply in the Big Lost River Basin and inform future groundwater modeling. The Idaho Geological Survey has prepared the groundwater budgets as part of a larger hydrogeologic investigation completed by the U.S. Geological Survey and the Idaho Geological Survey in cooperation with the Idaho Department of Water Resources during 2018–21. Other reports describe the hydrogeologic framework and several streamflow-measurement events to evaluate gains and losses on the Big Lost River. Collectively, these reports provide an updated characterization of groundwater resources in the Big Lost River Basin which will help address water resources challenges. A groundwater budget is a conceptual and numerical accounting of inflow (recharge) to groundwater and outflow (discharge) from groundwater. The predominant sources of recharge to the aquifer include losing river reaches (33 percent), areal recharge (as precipitation less evapotranspiration and surface runoff, comprising about 23 percent of the total inflow), tributary canyon underflow from higher altitudes (20 percent), canal seepage (13 percent), recharge through applied irrigation on fields below the root zone and other minor sources (11 percent), and Mackay Reservoir seepage (less than 1 percent). The primary sources of discharge from the aquifer are groundwater pumpage to meet irrigation demand, domestic supply, and municipal supply (76 percent) and gaining river reaches (24 percent). The positive or negative difference between the sum of all inflows and outflows is regarded as the residual, representing the change in groundwater storage, groundwater outflow from the basin or subbasins, and uncertainty and errors in the budget. In the Big Lost River Basin, groundwater outflow is at the mouth of the basin below Arco into the eastern Snake River Plain aquifer. The total mean annual estimated recharge to the Big Lost River Basin was 439,100 acre-feet per year (acre-ft/yr; 607 cubic feet per second [ft 3 /s]) for 2000–19, 373,900 acre-ft/yr (516 ft 3 /s) in 2014, and 762,100 acre-ft/yr (1,053 ft 3 /s) in 2017. The mean annual estimated groundwater discharge from the aquifer was about 112,300 acre-ft/yr (155 ft 3 /s) for 2000–19, 153,500 acre-ft/yr (212 ft 3 /s) in 2014, and 53,400 acre-ft/yr (74 ft 3 /s) in 2017. The estimated mean annual groundwater residual was 326,800 acre-ft/yr (451 ft 3 /s) for 2000–19, 220,400 acre-ft/yr (304 ft 3 /s) in 2014, and 708,700 acre-ft/yr (979 ft 3 /s) in 2017. The mean annual residual above Mackay Dam was 100,400 acre-ft/yr (2000-19), 96,700 acre-ft (2014), and 248,300 acre-ft (2017). The mean annual residual contribution below Mackay Dam, minus any groundwater-flow above Mackay Dam, was 226,400 acre-ft/yr (2000-19), 123,700 acre-ft (2014), and 460,400 acre-ft (2017). These results are highly sensitive to assumptions about certain budget inflow parameters. In particular, the magnitude of the budget residuals during especially dry and wet periods is amplified by the groundwater-budget terms tributary streamflow and tributary underflow that contribute appreciable recharge but also have high uncertainty. The results of the groundwater-budget evaluation describe an interconnected and complex hydrologic response throughout the basin to various climatic and water-use trends. The part of the basin above Mackay Dam typically has a positive groundwater residual derived from snowmelt recharge to tributary canyons and areal recharge in excess of groundwater pumpage for irrigation demand. This supply is used to meet irrigation demand above Mackay Dam and to provide for water supply below Mackay Dam. On average, groundwater inflow from above Mackay Dam to below Mackay Dam, assuming negligible reservoir storage effects, accounts for about 25 percent of the total groundwater recharge below Mackay Dam. Considerable recharge to groundwater below Mackay Dam occurs through seepage from the Big Lost River and canals and ditches. Most groundwater discharge from the aquifer is through irrigation pumping. The water supply below Mackay Dam is highly dependent on available upstream surface-water flows, the magnitude of the groundwater residual from above Mackay Dam, and annual variability in local groundwater conditions.

Idaho↗
USGS2013

Hydrogeologic framework, arsenic distribution, and groundwater geochemistry of the glacial-sediment aquifer at the Auburn Road landfill superfund site, Londonderry, New Hampshire

Leachate continues to be generated from landfills at the Auburn Road Landfill Superfund Site in Londonderry, New Hampshire. Impermeable caps on the three landfills at the site inhibit direct infiltration of precipitation; however, high water-table conditions allow groundwater to interact with landfill materials from below, creating leachate and ultimately reducing conditions in downgradient groundwater. Reducing conditions can facilitate arsenic transport by allowing it to stay in solution or by liberating arsenic adsorbed to surfaces and from geologic sources, such as glacial sediments and bedrock. The site occupies a 180-acre parcel of land containing streams, ponds, wetlands, and former gravel pits located in glacial sediment. Four areas, totaling 14 acres, including three landfills and one septage lagoon, were used for waste disposal. The site was closed in 1980 after volatile organic compounds associated with industrial waste dumping were detected. The site was added to the U.S. Environmental Protection Agency National Priority List in 1982, and the landfills were capped in 1996. Although volatile organic compound concentrations in groundwater have declined substantially, some measurable concentrations remain. Temporally variable and persistent elevated arsenic concentrations have been measured in groundwater affected by the landfill leachate. Microbial consumption of carbon found in leachate is a driver of reducing conditions that liberate arsenic at the site. In addition to sources of carbon in landfill leachate, wetland areas throughout the site also could contribute carbon to groundwater, but it is currently unknown if any of the wetland areas have downward or reversing gradients that could allow the infiltration of surface water to groundwater. Red-stained sediments and water indicate iron-rich groundwater discharge to surface water and are also associated with elevated concentrations of arsenic in sediment and groundwater. Ironrich groundwater seeps have been observed in the wetland, streams, and pond downgradient of the landfills. Piezometers were installed in some of these locations to confirm groundwater discharge, measure vertical-flow gradients, and to provide a way to sample the discharging groundwater. Understanding the movement of leachate in groundwater is complicated by the presence of preferential flow paths through aquifer materials with differing hydraulic properties; these preferential flow paths can affect rates of recharge, geochemical conditions, and contaminant fluxes. In areas adjacent to the three capped landfills, infiltration of precipitation containing oxygenated water through permeable deltaic sediments in the former gravel pit area causes increases in dissolved oxygen concentrations and decreases in arsenic concentrations. Layered deltaic sediments produce anisotropic hydraulic characteristics and zones of high hydraulic conductivity. The glacial-sediment aquifer also includes glaciolacustrine sediments that have low permeability and limit infiltration at the surface Discharge of leachate-affected groundwater may be limited in areas of organic muck on the bottom of Whispering Pines Pond because the muck may act as a semiconfining layer. Geophysical survey results were used to identify several areas with continuous beds of muck and an underlying highresistivity layer on top of a layer of low resistivity that may represent leachate-affected groundwater. The high-resistivity layer is likely groundwater associated with oxygenated recharge, which would cause arsenic to adsorb onto aquifer sediments and reduce concentrations of dissolved arsenic in groundwater. Surface and borehole geophysical data collected in 2011 were used to identify potentially high-permeability or contaminated zones in the aquifer (preferential flowpaths) as well as low-permeability zones that may promote contamination through back diffusion. Some groundwater in parts of the glacial-sediment aquifer where the leachate plumes were present had low electrical resistivity, low dissolved oxygen, and high concentrations of arsenic. Low-resistivity zones in the underlying bedrock were associated with fractures that also may contain leachate. Although surveying the fractured bedrock was not a specific objective of this study, the results suggest that such a survey would help to determine if leachate and associated concentrations of arsenic are migrating downward into the fractured-bedrock-aquifer system. An uncalibrated, one-dimensional, reactive-transport model was used to assess several conditions that affect arsenic mobility. The results indicate that reductive dissolution and desorption from glacial sediments control dissolved arsenic concentrations. Parameter sensitivity analysis was used to identify key data that are needed in order to accurately assess the time required for arsenic concentrations to fall to levels below the maximum contaminant level at the site. Quantifying this time will require accurate characterization of carbon, sediment-surface sorption sites, and groundwater fluxes at the site.

New Hampshire↗
USGS2023

Modeling groundwater-level responses to multiple stresses using transfer-function models and wavelet analysis in a coastal aquifer system

In coastal aquifers, dynamic stresses such as climate forcings, groundwater withdrawals, and ocean tidal fluctuations cause nonlinear responses to groundwater levels. Such responses to the stresses impact groundwater resources and related flooding and infrastructure risks at multiple scales. We used time-series models such as transfer-function models and wavelet analysis to quantify the relative contribution of these stresses to groundwater-level fluctuation in wells from the unconfined and confined aquifers in an Atlantic coastal aquifer. Climate forcings, such as precipitation and temperature, explained most of the groundwater-level variation for wells in the unconfined aquifer, whereas groundwater withdrawals were the dominant driver of groundwater levels for wells in the confined aquifer. The impact of groundwater withdrawals also was detected in several wells in the unconfined aquifer. Although the influence of ocean tides on groundwater levels commonly is observed in coastal aquifers, we found that daily groundwater withdrawals can obscure the semi-diurnal coherence signal of the two series. The magnitude of groundwater-level fluctuation that could be explained solely by tides was minor compared to that explained by climate or withdrawal stresses. Transfer-function modeling showed seasonal withdrawals from wells in confined aquifers had a significant, yet heterogeneous influence on groundwater levels in coastal aquifers, which highlights climate and withdrawals as key compounding stresses in coastal hydrology. This study demonstrates the value of time-series approaches to advance characterization of groundwater systems in areas with limited hydrogeologic parameter information.

Virginia↗
USGS2014

Groundwater-quality data in the North San Francisco Bay Shallow Aquifer study unit, 2012: results from the California GAMA Program

Groundwater quality in the 1,850-square-mile North San Francisco Bay Shallow Aquifer (NSF-SA) study unit was investigated by the U.S. Geological Survey (USGS) from April to August 2012, as part of the California State Water Resources Control Board (SWRCB) Groundwater Ambient Monitoring and Assessment (GAMA) Program’s Priority Basin Project (PBP). The GAMA-PBP was developed in response to the California Groundwater Quality Monitoring Act of 2001 and is being conducted in collaboration with the SWRCB and Lawrence Livermore National Laboratory (LLNL). The NSF-SA study unit was the first study unit to be sampled as part of the second phase of the GAMA-PBP, which focuses on the shallow aquifer system. The GAMA NSF-SA study was designed to provide a spatially unbiased assessment of untreated-groundwater quality in the shallow aquifer systems and to facilitate statistically consistent comparisons of untreated-groundwater quality throughout California. The shallow aquifer system in the NSF-SA study unit was defined as the part of the aquifer system that is used by many private domestic wells and is shallower than the primary aquifer system used by many public-supply wells. In the NSF-SA study unit located in Marin, Mendocino, Napa, Solano, and Sonoma Counties, groundwater samples were collected from 71 wells. Seventy of the wells were selected by using a spatially distributed, randomized grid-based method to provide statistical representation of the study unit (grid wells), and one well was selected to aid in evaluation of water-quality issues (understanding well). The groundwater samples were analyzed for organic constituents (volatile organic compounds [VOCs], pesticides, and pesticide degradates); constituents of special interest (perchlorate and 1,2,3-trichloropropane [1,2,3-TCP]); naturally occurring inorganic constituents (trace elements, nutrients, major and minor ions, silica, and total dissolved solids [TDS]); and radioactive constituents (radon-222 and gross alpha and gross beta radioactivity). Naturally occurring isotopes (stable isotopes of hydrogen, oxygen, boron, strontium, and inorganic carbon in water, tritium activities, and carbon-14 abundances) were measured to help identify the sources and ages of the sampled groundwater. In total, 207 constituents and water-quality indicators were measured. Three types of quality-control samples (blanks, replicates, and matrix spikes) were collected at up to 13 percent of the wells in the NSF-SA study unit, and the results for these samples were used to evaluate the quality of the data for the groundwater samples. Blanks rarely contained detectable concentrations of any constituent, suggesting that contamination from sample-collection procedures was not a significant source of bias in the data for the groundwater samples. Replicate samples generally were within the limits of acceptable analytical reproducibility. Matrix-spike recoveries were within the acceptable range (70 to 130 percent) for approximately 91 percent of the compounds. Most of the wells sampled for this study were private domestic wells. Private domestic wells are not regulated in California, and groundwater from these wells is rarely analyzed for water-quality constituents. Although regulatory benchmarks for drinking-water quality do not apply to private domestic wells, to provide some context for the results, concentrations of constituents measured in the untreated groundwater were compared with regulatory and non-regulatory health-based benchmarks established by the U.S. Environmental Protection Agency (USEPA) and California Department of Public Health (CDPH), to non-regulatory health-based benchmarks established by the USGS in cooperation with the USEPA, and to non-regulatory benchmarks established for aesthetic concerns by the CDPH. Comparisons between data collected for this study and benchmarks for drinking water are for illustrative purposes only and are not indicative of compliance or non-compliance with those benchmarks. Most of the organic and inorganic constituents that were detected in groundwater samples from the 70 grid wells in the NSF-SA study unit were detected at concentrations less than drinking-water benchmarks. Of the 149 organic and special-interest constituents analyzed for in groundwater samples, 31 were detected; concentrations of most detected constituents were less than regulatory and non-regulatory health-based benchmarks. One VOC, benzene, and one insecticide, dieldrin, were detected at concentrations above their respective health-based benchmarks. In total, VOCs were detected in 40 percent of the grid wells sampled, pesticides and pesticide degradates were detected in 13 percent, and perchlorate was detected in 27 percent of the 70 grid wells sampled. Groundwater samples from 70 grid wells were analyzed for trace elements, major and minor ions, nutrients, and radioactive constituents; most detected concentrations were less than health-based benchmarks. Exceptions are 12 detections of manganese greater than the USGS Health-Based Screening Level (HBSL), 7 detections of arsenic greater than the USEPA maximum contaminant level (MCL-US) of 10 micrograms per liter (μg/L), 2 detections of boron greater than the HBSL of 6,000 μg/L, 2 detections of fluoride greater than the CDPH maximum contaminant level (MCL-CA) of 2 milligrams per liter (mg/L), 2 detections of nitrate greater than the MCL-US of 10 mg/L, and two detections of radon-222 greater than the proposed MCL-US of 4,000 picocuries per liter. Results for constituents with non-regulatory benchmarks set for aesthetic concerns from the grid wells showed that iron concentrations greater than the CDPH secondary maximum contaminant level (SMCL-CA) of 300 μg/L were detected in 13 grid wells. Chloride was detected at a concentration greater than the SMCL-CA recommended benchmark of 250 mg/L in two grid wells. Sulfate concentrations greater than the SMCL-CA recommended benchmark of 250 mg/L were measured in two grid wells, and the concentration in one of these wells was also greater than the SMCL-CA upper benchmark of 500 mg/L. TDS concentrations greater than the SMCL-CA recommended benchmark of 500 mg/L were measured in 15 grid wells, and concentrations in 4 of these wells were also greater than the SMCL-CA upper benchmark of 1,000 mg/L.

California↗
USGS1998

The chemical evolution of Kurnub Group palcowater in the Sinai-Negev province - A mass balance approach

The chemical evolution of the Kurnub Group paleowater was studied starting from rainwater in recharge areas of the Sinai and along groundwater flowpaths leading to the natural outlets of this regional aquifer. This was achieved by investigating the chemical composition of groundwater, ionic ratios, degrees of saturation with common mineral species, normative analysis of dissolved salts and by modeling of rock/water interaction and mixing processes occurring along groundwater flow paths. The initial groundwater composition used is from the Nakhel well in Sinai. It evolves from desert rainwater percolating through typical Kurnub Group lithology in Sinai. This rainwater dissolves mainly gypsum, halite and dolomite together with smaller amounts of marine aerosol and K-feldspar. At the same time it precipitates calcite, SiO2, smectite and degasses CO2. Between the area of Nakhel and the northern Negev the chemistry of Kurnub Group waters is influenced by dissolution of halite and lesser amounts of gypsum of surficial origin in recharge areas, small amounts of feldspars and of dolomite cement in sandstones eroded from the Arabo-Nubian igneous massif of Sinai and organic degradation-derived CO2. Concomitantly, there is precipitation of calcite, smectite, SiO2 and probably analcime characteristic of sediments in continental closed basins. North of the Negev, the Kurnub Group fluids are diluted and altered by mixing with Judea Group aquifer groundwaters. On the E there is mixing with residual brines from the water body ancestral to the Dead Sea, prior to discharge into the Arava valley. Rock/water interaction indicated by NETPATH and PHREEQC modeling is in agreement with lithology and facies changes previously observed in the Kurnub Group sequence.

Applied Geochemistry↗
USGS2018

Assessment of water resources and the potential effects from oil and gas development in the Bureau of Land Management Tri-County planning area, Sierra, Doña Ana, and Otero Counties, New Mexico

The U.S. Geological Survey (USGS), in cooperation with the Bureau of Land Management, conducted a study to assess the water resources and potential effects on the water resources from oil and gas development in the Tri-County planning area, Sierra, Doña Ana, and Otero Counties, New Mexico. Publicly available data were used to assess these resources and effects and to identify data gaps in the Tri-County planning area. The Tri-County planning area includes approximately 9.3 million acres and is within the eastern extent of the Basin and Range Province, which consists of mountain ranges and low elevation basins. Three specific areas of interest within the Tri-County planning area are the Jornada del Muerto, Tularosa Basin, and Otero Mesa, which is adjacent to the Salt Basin. Surface-water resources are limited in the Tri-County planning area, with the Rio Grande as the main perennial river flowing from north to south through Sierra and Doña Ana Counties. The Tularosa Creek is an important surface-water resource in the Tularosa Basin. The Sacramento River, which flows southeast out of the Sacramento Mountains, is an important source of recharge to aquifers in the Salt Basin. Groundwater resources vary in aquifer type, depth to water, and water quality. For example, the Jornada del Muerto, Tularosa Basin, and Salt Basin each have shallow and deep aquifer systems, and water can range from freshwater, with less than 1,000 milligrams per liter (mg/L) of total dissolved solids, to brine, with greater than 35,000 mg/L of total dissolved solids. Water quality in the Tri-County planning area is affected by the dissolution of salt deposits and evaporation which are common in arid regions such as southern New Mexico. The potential for oil and gas development exists in several areas within the Tri-County area. As many as 81 new conventional wells and 25 coalbed natural gas wells could be developed by 2035. Conventional oil and gas well construction in the Tri-County planning area is expected to require 1.53 acre-feet (acre-ft) (500,000 gallons) of water per well, similar to requirements in the nearby Permian Basin of New Mexico, while construction of unconventional wells is expected to require 7.3 acre-ft of water per well. Produced waters in the Permian Basin have high total dissolved solids, in the brackish to brine range. Data gaps identified in this study include the limited detailed data on surface-water resources, the lack of groundwater data in areas of interest, and the lack of water chemistry data related to oil and gas development issues. Surface waters in the Tri-County planning area are sparse; some streams are perennial, and most are ephemeral. A more detailed study of the ephemeral channels and their interaction with groundwater could provide a better understanding of the importance of these surface-water resources. Groundwater data used in this study are from the USGS National Water Information System, which does not have continuous water-level depth data at many of the sites in the Tri-County planning area. On Otero Mesa, no recurrent groundwater-level data are available at any one site. The water-quality data compiled in this study provide a good overview of the general chemistry of groundwater in the Tri-County planning area. To fully understand the groundwater resources, it would be helpful to have more wells in specific areas of interest for groundwater-level and water-quality measurements.

New Mexico↗
USGS2021

Numerical modeling of groundwater flow in the crystalline-rock aquifer in the vicinity of the Savage Municipal Water-Supply Well Superfund site, Milford, New Hampshire

In 2010, tetrachloroethylene (PCE), a chlorinated volatile organic compound, was detected in groundwater from deep (more than 300 feet below land surface) fractures in monitoring wells tapping a crystalline-rock aquifer. The aquifer underlies the Milford-Souhegan glacial-drift aquifer, a high water-producing aquifer, and the Savage Municipal Water-Supply Well Superfund site in Milford, New Hampshire. Between 30 and 40 residential water-supply wells are near (0.25 mile north of) the PCE-contaminated monitoring wells. Some of the residential water-supply wells are likely installed in similar rock types and formations as those of the monitoring wells installed as part of the Superfund site. As of 2020, periodic sampling by the U.S. Environmental Protection Agency and New Hampshire Department of Environmental Services (cooperative partners for this study) since 1996 had not detected PCE in groundwater from the residential water-supply wells. Nevertheless, understanding the vulnerability of the residential water wells to capture PCE contaminated groundwater was of concern. A numerical groundwater flow model was developed by the U.S. Geological Survey to assess groundwater flow and advective transport of PCE-contaminated groundwater in the crystalline-rock aquifer of the Milford area. The model (called the area-wide model) encompasses a 26.5-square mile area to allow for more accurate computation of water fluxes near the PCE-contaminated monitoring wells and the residential water wells. Simulations with the area-wide model show that, with the 2016 configuration of residential wells, capture of PCE by the residential water wells appears unlikely for hydrologic conditions typical of 2010 based on steady-state, advective transport modeling. However, simulations also show that adding residential water wells to the north of the PCE-contaminated monitoring wells could affect the transport of PCE. Groundwater withdrawals at other adjacent wells in the overlying Milford-Souhegan glacial-drift aquifer affect advective transport in the crystalline-rock aquifer. Therefore, the potential for future changes in withdrawals in the area, as well as changes in hydrologic conditions, including groundwater recharge and streamflow amounts, should be considered in the remedial assessment process. The development of the area-wide model and linkages established by this study with previously developed Milford-Souhegan glacial-drift aquifer transport models will help facilitate the development of remedial strategies for this Superfund site.

New Hampshire↗
USGS2010

Comment on "Horizontal aquifer movement in a theis-theim confined system" by Donald C. Helm

In a recent paper, Helm [1994] presents an analysis of horizontal aquifer movement induced by groundwater withdrawal from a confined aquifer in which fluid and grains are incompressible. The analysis considers the aquifer in isolation (ignoring overlying and underlying strata) and assumes that the aquifer deforms purely in the horizontal direction (with no vertical movement). Helm's solution for grain displacement is obtained through introduction of a quantity known as bulk flux, q b , defined as q b = n v w + (1 - n )v s where n is porosity, v w is velocity of water, and v s is the velocity of the solid grains. On the basis of the bulk flux concept, Helm develops an explanation for the driving force on the bulk material. It is our view that Helm's analysis is subject to four limitations. First, Helm's assumption of zero vertical displacement is not supported by field observations and could result in over- estimation of radial displacement. Second, in ignoring the role of overlying and underlying strata, Helm's solution does not yield reliable estimates of aquifer deformation. Third, Helm's solution method works only for problems that involve one spatial coordinate (for example, x or r) but does not generally work for problems involving three-dimensional flow and de- formation. Fourth, Helm's explanation of the driving force on the bulk material is faulty for general three-dimensional problems. The purpose of our comment is to discuss these four issues.

Water Resources Research↗
USGS2017

Groundwater-level trends in the U.S. glacial aquifer system, 1964-2013

The glacial aquifer system in the United States is a major source of water supply but previous work on historical groundwater trends across the system is lacking. Trends in annual minimum, mean, and maximum groundwater levels for 205 monitoring wells were analyzed across three regions of the system (East, Central, West Central) for four time periods: 1964-2013, 1974-2013, 1984-2013, and 1994-2013. Trends were computed separately for wells in the glacial aquifer system with low potential for human influence on groundwater levels and ones with high potential influence from activities such as groundwater pumping. Generally there were more wells with significantly increasing groundwater levels (levels closer to ground surface) than wells with significantly decreasing levels. The highest numbers of significant increases for all four time periods were with annual minimum and/or mean levels. There were many more wells with significant increases from 1964 to 2013 than from more recent periods, consistent with low precipitation in the 1960s. Overall there were low numbers of wells with significantly decreasing trends regardless of time period considered; the highest number of these were generally for annual minimum groundwater levels at wells with likely human influence. There were substantial differences in the number of wells with significant groundwater-level trends over time, depending on whether the historical time series are assumed to be independent, have short-term persistence, or have long-term persistence. Mean annual groundwater levels have significant lag-one-year autocorrelation at 26.0% of wells in the East region, 65.4% of wells in the Central region, and 100% of wells in the West Central region. Annual precipitation across the glacial aquifer system, on the other hand, has significant autocorrelation at only 5.5% of stations, about the percentage expected due to chance.

Journal of Hydrology↗
USGS2021

Hydrogeology and model-simulated groundwater availability in the Salt Fork Red River aquifer, southwestern Oklahoma, 1980–2015

The 1973 Oklahoma Water Law (82 OK Stat § 82-1020.5) requires that the Oklahoma Water Resources Board (OWRB) conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin (hereinafter referred to as an “aquifer”). The maximum annual yield allocated per acre of land is known as the equal-proportionate-share (EPS) pumping rate. At present (2021), the OWRB has not yet established a maximum annual yield and EPS pumping rate for the Salt Fork Red River aquifer. To provide updated information to the OWRB that could support evaluation and determination of an appropriate maximum annual yield, the U.S. Geological Survey (USGS), in cooperation with the OWRB, conducted a hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater availability in the Salt Fork Red River aquifer. The Salt Fork Red River aquifer in Greer, Harmon, and Jackson Counties of southwestern Oklahoma is composed of about 274.5 square miles of alluvium and terrace deposits associated with the Salt Fork Red River. The mean annual recharge rate to the Salt Fork Red River aquifer for the period 1980–2015 was estimated to be about 2.94 inches per year, or 10.0 percent of the mean annual precipitation for the same period (29.4 inches per year). This 1980–2015 mean annual recharge rate is equivalent to a mean annual recharge rate of about 38,000 acre-feet per year (acre-ft/yr) for the Salt Fork Red River aquifer excluding about 19,764 acres comprising the Mulberry Creek and Horse Creek terraces. The mean annual recharge rates upgradient and downgradient from USGS streamgage 07300500 Salt Fork Red River at Mangum, Okla. (hereinafter referred to as the “Mangum gage”), apportioned by aquifer area (41.5 and 58.5 percent, respectively), were about 16,000 and 22,000 acre-ft/yr, respectively. Mean annual groundwater use for the study period (1980–2015) was 3,532.7 acre-ft/yr; about 77 percent of that groundwater use was for irrigation, and about 23 percent was for public supply. Most groundwater use for irrigation was associated with wells in the Martha terrace. A hydrogeologic framework was developed for the Salt Fork Red River aquifer and included a definition of the aquifer extent and potentiometric surface, as well as a description of the textural and hydraulic properties of aquifer materials. The hydrogeologic framework was used in the construction of the numerical groundwater-flow model of the Salt Fork Red River aquifer described in this report. A conceptual model for the Salt Fork Red River aquifer that reasonably represents the groundwater-flow system was developed to constrain the construction and calibration of the numerical model. The conceptual-model water budget estimated mean annual inflows to, and outflows from, the Salt Fork Red River aquifer for the period 1980–2015 and included a subaccounting of mean annual inflows and outflows for the portions of the aquifer that were upgradient and downgradient from the Mangum gage. The numerical groundwater-flow model of the Salt Fork Red River aquifer was constructed by using MODFLOW-2005 with the Newton formulation solver. The model of the Salt Fork Red River aquifer was spatially discretized into 1,050 rows, 1,125 columns, about 170,000 active cells measuring 200 by 200 feet (ft), and a single convertible layer. The model was temporally discretized into 432 monthly transient stress periods (each with two time steps to improve model stability). An initial steady-state stress period represented mean annual inflows to, and outflows from, the aquifer and produced a solution that was used as the initial condition for subsequent transient stress periods as well as some groundwater-availability scenarios. The model was calibrated to water-table-altitude observations at selected wells and base-flow observations at selected streamgages. The simulated saturated thickness of the Salt Fork Red River aquifer was determined by subtracting the altitude of the aquifer base from the simulated water-table altitude at the end of the numerical-model period (2015). The simulated saturated thickness was more than 75 ft in a paleochannel in the Dodson terrace near the Texas border. The mean aquifer thickness (sum of saturated and unsaturated) was 49.62 ft, and the mean saturated thickness was 28.55 ft. A simulated mean transmissivity of 1,024 feet squared per day was computed from the calibrated hydraulic conductivity and saturated thickness of each cell. The simulated available water in storage at the end of the numerical-model period (2015) was 526,117 acre-feet (acre-ft); about 42 percent of that total was available upgradient from the Mangum gage, and about 58 percent of that total was available downgradient from the Mangum gage (including the Mangum terrace). Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate the EPS pumping rate that ensures a minimum 20-, 40-, and 50-year life of the aquifer, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought on base flow and groundwater storage. The 20-, 40-, and 50-year EPS pumping rates under normal recharge conditions were about 0.51, 0.48, and 0.48 acre-foot per acre per year, respectively. Given the 155,929-acre modeled aquifer area, these rates correspond to annual yields of about 78,800, 74,900, and 74,700 acre-ft/yr, respectively. For the 20-year EPS scenario, decreasing and increasing recharge by 10 percent resulted in a 6-percent change in the EPS pumping rate in both cases; for the 40- and 50-year EPS scenarios, decreasing and increasing recharge by 10 percent resulted in a 7-percent change in the EPS pumping rate in both cases. Projected 50-year pumping scenarios were used to simulate the effects of selected well withdrawal rates on groundwater storage of the Salt Fork Red River aquifer and base flows in the Salt Fork Red River. The effects of well withdrawals were evaluated by quantifying differences in groundwater storage and base flow in four 50-year scenarios, which applied (1) no groundwater pumping, (2) mean pumping rates for the study period (1980–2015), (3) 2015 pumping rates, and (4) increasing demand pumping rates at simulated wells. The increasing demand pumping rates assumed a cumulative 20.4-percent increase in pumping over 50 years based on 2010–60 demand projections for southwestern Oklahoma. Groundwater storage after 50 years with no pumping was 535,000 acre-ft, or 8,900 acre-ft (1.7 percent) greater than the initial groundwater storage; this groundwater storage increase is equivalent to a mean water-table-altitude increase of 0.48 ft. Groundwater storage after 50 years of pumping at the mean rate for the study period (1980–2015) was 519,900 acre-ft, or 6,200 acre-ft (1.2 percent) less than the initial groundwater storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.34 ft. Groundwater storage at the end of the 50-year period with 2015 pumping rates was 513,100 acre-ft, or 13,000 acre-ft (2.5 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.71 ft. Groundwater storage at the end of the 50-year period with increasing demand pumping rates was 509,700 acre-ft, or 16,500 acre-ft (3.1 percent) less than the initial storage; this groundwater storage decrease is equivalent to a mean water-table-altitude decline of 0.89 ft. A hypothetical 10-year drought scenario was used to simulate the effects of a prolonged period of reduced recharge on groundwater storage. The period January 1983–December 1992 was chosen as the simulated drought period. Drought effects were quantified by comparing the results of the drought scenario to those of the calibrated numerical model (no drought) at the end of the simulated drought period (1992). To simulate the hypothetical drought, recharge in the calibrated numerical model was reduced by 50 percent during the simulated drought period (1983–92). Upstream inflows from the Salt Fork Red River, Turkey Creek, and Bitter Creek were reduced by 75 percent. Groundwater storage at the end of the drought period (1992) was 479,200 acre-ft, or 53,200 acre-ft (10.0 percent) less than the groundwater storage of the calibrated numerical model at the end of the drought period. This decrease in groundwater storage is equivalent to a mean water-table-altitude decline of 2.9 ft. At the end of the 10-year hypothetical drought period, simulated base flows at the Mangum gage and USGS streamgage 07301110 Salt Fork Red River near Elmer, Okla., had decreased by about 80 and 70 percent, respectively.

Oklahoma↗
USGS2011

Regional groundwater-flow model of the Redwall-Muav, Coconino, and alluvial basin aquifer systems of northern and central Arizona

A numerical flow model (MODFLOW) of the groundwater flow system in the primary aquifers in northern Arizona was developed to simulate interactions between the aquifers, perennial streams, and springs for predevelopment and transient conditions during 1910 through 2005. Simulated aquifers include the Redwall-Muav, Coconino, and basin-fill aquifers. Perennial stream reaches and springs that derive base flow from the aquifers were simulated, including the Colorado River, Little Colorado River, Salt River, Verde River, and perennial reaches of tributary streams. Simulated major springs include Blue Spring, Del Rio Springs, Havasu Springs, Verde River headwater springs, several springs that discharge adjacent to major Verde River tributaries, and many springs that discharge to the Colorado River. Estimates of aquifer hydraulic properties and groundwater budgets were developed from published reports and groundwater-flow models. Spatial extents of aquifers and confining units were developed from geologic data, geophysical models, a groundwater-flow model for the Prescott Active Management Area, drill logs, geologic logs, and geophysical logs. Spatial and temporal distributions of natural recharge were developed by using a water-balance model that estimates recharge from direct infiltration. Additional natural recharge from ephemeral channel infiltration was simulated in alluvial basins. Recharge at wastewater treatment facilities and incidental recharge at agricultural fields and golf courses were also simulated. Estimates of predevelopment rates of groundwater discharge to streams, springs, and evapotranspiration by phreatophytes were derived from previous reports and on the basis of streamflow records at gages. Annual estimates of groundwater withdrawals for agriculture, municipal, industrial, and domestic uses were developed from several sources, including reported withdrawals for nonexempt wells, estimated crop requirements for agricultural wells, and estimated per capita water use for exempt wells. Accuracy of the simulated groundwater-flow system was evaluated by using observational control from water levels in wells, estimates of base flow from streamflow records, and estimates of spring discharge. Major results from the simulations include the importance of variations in recharge rates throughout the study area and recharge along ephemeral and losing stream reaches in alluvial basins. Insights about the groundwater-flow systems in individual basins include the hydrologic influence of geologic structures in some areas and that stream-aquifer interactions along the lower part of the Little Colorado River are an effective control on water level distributions throughout the Little Colorado River Plateau basin. Better information on several aspects of the groundwater flow system are needed to reduce uncertainty of the simulated system. Many areas lack documentation of the response of the groundwater system to changes in withdrawals and recharge. Data needed to define groundwater flow between vertically adjacent water-bearing units is lacking in many areas. Distributions of recharge along losing stream reaches are poorly defined. Extents of aquifers and alluvial lithologies are poorly defined in parts of the Big Chino and Verde Valley sub-basins. Aquifer storage properties are poorly defined throughout most of the study area. Little data exist to define the hydrologic importance of geologic structures such as faults and fractures. Discharge of regional groundwater flow to the Verde River is difficult to identify in the Verde Valley sub-basin because of unknown contributions from deep percolation of excess surface water irrigation.

Scientific Investigations Report↗
USGS2023

Where the past meets the present: Connecting nitrogen from watersheds to streams through groundwater flowpaths

Groundwater discharge to streams is a nonpoint source of nitrogen (N) that confounds N mitigation efforts and represents a significant portion of the annual N loading to watersheds. However, we lack an understanding of where and how much groundwater N enters streams and watersheds. Nitrogen concentrations at the end of groundwater flowpaths are the culmination of biogeochemical and physical processes from the contributing land area where groundwater recharges, within the aquifer system, and in the near-stream riparian area where groundwater discharges to streams. Our research objectives were to quantify the spatial distribution of N concentrations at groundwater discharges throughout a mixed land-use watershed and to evaluate how relationships among contributing and riparian land cover, modeled aquifer characteristics, and groundwater discharge biogeochemistry explain the spatial variation in groundwater discharge N concentrations. We accomplished this by integrating high-resolution thermal infrared surveys to locate groundwater discharge, biogeochemical sampling of groundwater, and a particle tracking model that links groundwater discharge locations to their contributing area land cover. Groundwater N loading from groundwater discharges within the watershed varied substantially between and within streambank groundwater discharge features. Groundwater nitrate concentrations were spatially heterogeneous ranging from below 0.03–11.45 mg-N/L, varying up to 20-fold within meters. When combined with the particle tracking model results and land cover metrics, we found that groundwater discharge nitrate concentrations were best predicted by a linear mixed-effect model that explained over 60% of the variation in nitrate concentrations, including aquifer chemistry (dissolved oxygen, Cl − , SO 4 2− ), riparian area forested land cover, and modeled physical aquifer characteristics (discharge, Euclidean distance). Our work highlights the significant spatial variability in groundwater discharge nitrate concentrations within mixed land-use watersheds and the need to understand groundwater N processing across the many spatiotemporal scales within groundwater cycling.

Connecticut↗
USGS1999

The role of groundwater chemistry in the transport of bacteria to water-supply wells

Static mini-columns and in situ injection and recovery tests were used to assess the effects of modest changes in groundwater chemistry upon the pH-dependence of bacterial attachment, a primary determinant of bacterial mobility in drinking water aquifers. In uncontaminated groundwater (<1 mg l-1 dissolved organic carbon, DOC), bacterial attachment to aquifer grain surfaces declined steadily from 93 to 20% in response to an increase in pH from 5.8 to 7.8. However, bacterial attachment in modestly-contaminated groundwater (4 mg l-1 DOC) was relatively insensitive to pH change from pH 3.5 to pH 8, as was bacterial attachment in uncontaminated groundwater amended with only ~3 mg l-1 of purified humic acid. Destruction by UV-oxidation of the DOC in contaminated groundwater partially restored the pH-dependence of bacterial attachment. Results from static column tests and from a small-scale (3.6 m) natural-gradient injection and recovery study suggest that low concentrations of surfactants can also substantively alter the attraction of groundwater bacteria for grain surfaces and, therefore can alter the transport of bacteria to water-supply wells. This phenomenon was pH-sensitive and dependent upon the nature of the surfactant. At pH 7.6, 200 mg l-1 of the non-ionic surfactant, Imbentin, caused a doubling of fractional bacterial attachment in aquifer-sediment columns, but had little effect under slightly acidic conditions (e.g. at pH 5.8). In contrast, 1 mg l-1 of linear alkylbenzene sulphonate (LAS) surfactant, a common sewage-derived contaminant, decreased the fractional bacterial attachment by more than 30% at pH 5.8, but had little effect at pH 7.3.Static mini-columns and in situ injection and recovery tests were used to assess the effects of modest changes in groundwater chemistry upon the pH-dependence of bacterial attachment, a primary determinant of bacterial mobility in drinking water aquifers. In uncontaminated groundwater (<1 mg l-1 dissolved organic carbon, DOC), bacterial attachment to aquifer grain surfaces declined steadily from 93 to 20% in response to an increase in pH from 5.8 to 7.8. However, bacterial attachment in modestly-contaminated groundwater (4 mg l-1 DOC) was relatively insensitive to pH change from pH 3.5 to pH 8, as was bacterial attachment in uncontaminated groundwater amended with only approx. 3 mg l-1 of purified humic acid. Destruction of UV-oxidation of the DOC in contaminated groundwater partially restored the pH-dependence of bacterial attachment. Results from the static column tests and from a small-scale (3.6 m) natural-gradient injection and recovery study suggest that low concentrations of surfactants can also substantively alter the attraction of groundwater bacteria for grain surfaces and, therefore can alter the transport of bacteria to water-supply wells. This phenomenon was pH-sensitive and dependent upon the nature of the surfactant. At pH 7.6, 200 mg l-1 of the non-ionic surfactant, Imbentin, caused a doubling of fractional bacterial attachment in aquifer-sediment columns, but had little effect under slightly acidic conditions (e.g. at pH 5.8). In contrast, 1 mg l-1 of linear alkylbenzene sulphonate (LAS) surfactant, a common sewage-derived contaminant, decreased the fractional bacterial attachment by more than 30% at pH 5.8, but had little effect at pH 7.3.

Conference Paper↗