USGS Science⌕ Search

SEARCH · USGS Science

Results for “Groundwater”

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

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

At least 307 records · Page 17Linked to original sources
USGS2015

Status and understanding of groundwater quality in the Cascade Range and Modoc Plateau study unit, 2010: California GAMA Priority Basin Project

Groundwater quality in the Cascade Range and Modoc Plateau study unit was investigated as part of the California State Water Resources Control Board’s Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project. The study was designed to provide a statistically unbiased assessment of untreated groundwater quality in the primary aquifer system. The depth of the primary aquifer system for the Cascade Range and Modoc Plateau study unit was delineated by the depths of the screened or open intervals of wells in the State of California’s database of public-supply wells. Two types of assessments were made: a status assessment that described the current quality of the groundwater resource, and an understanding assessment that made evaluations of relations between groundwater quality and potential explanatory factors representing characteristics of the primary aquifer system. The assessments characterize the quality of untreated groundwater, not the quality of treated drinking water delivered to consumers by water distributors. The status assessment was based on water-quality data collected in 2010 by the U.S. Geological Survey from 90 wells and springs (USGS-grid wells) and on water-quality data compiled from the State of California’s regulatory compliance database for samples collected from 240 public-supply wells between September 2007 and September 2010. To provide context, the water-quality data discussed in this report were compared to California and Federal drinking-water regulatory and non-regulatory benchmarks for treated drinking water. Groundwater quality is defined in terms of relative concentrations (RCs), which are calculated by dividing the concentration of a constituent in groundwater by the concentration of the benchmark for that constituent. The RCs for inorganic constituents (major ions, trace elements, nutrients, and radioactive constituents) were classified as “high” (the RC is greater than 1.0, indicating that the concentration is above the benchmark), “moderate” (the RC is from 1.0 to greater than 0.5), or “low” (the RC is less than or equal to 0.5). For organic constituents (volatile organic compounds and pesticides) and special-interest constituents (perchlorate), the boundary between moderate and low RCs was set at 0.1. All benchmarks used for organic constituents were health-based. For inorganic constituents, health-based and aesthetic-based benchmarks were used. Constituents without benchmarks were not considered in the status assessment . The primary metric used for quantifying regional-scale groundwater quality was the aquifer-scale proportion—the areal percentages of the primary aquifer system with high, moderate, and low RCs for a given constituent or class of constituents. The study unit was divided into six study areas on the basis of geologic differences (Eastside Sacramento Valley, Honey Lake Valley groundwater basin, Cascade Range and Modoc Plateau Low Use Basins, Quaternary Volcanic Areas, Shasta Valley and Mount Shasta Volcanic Area, and Tertiary Volcanic Areas), and each study area was divided into equal-area grid cells. Aquifer-scale proportions were calculated for individual constituents and constituent classes for each of the six study areas and for the study unit as a whole by using grid-based (one well per cell) and spatially weighted (many wells per cell) statistical methods. The status assessment showed that inorganic constituents were present at high and moderate RCs in greater proportions of the Cascade Range and Modoc Plateau study unit than were organic constituents. One or more inorganic constituents with health-based benchmarks were present at high RCs in 9.4 percent, and at moderate RCs in 14.7 percent of the primary aquifer system. Arsenic was present at high RCs in approximately 3 percent of the primary aquifer system; boron, molybdenum, uranium, and vanadium each were present at high RCs in approximately 2 percent of the primary aquifer system. One or more inorganic constituents with aesthetic-based benchmarks were present at high RCs in 15.1 percent of the primary aquifer system and at moderate RCs in 4.9 percent. Manganese, iron, and total dissolved solids were present at high RCs in approximately 12 percent, 5 percent, and 2 percent, respectively, of the primary aquifer system. Organic constituents were not detected at high or moderate RCs in the primary aquifer system, and one or more organic constituents were detected at low RCs in approximately 40 percent of the primary aquifer system. Two classes of organic constituents were detected in more than 10 percent of the primary aquifer system: trihalomethanes (chloroform only) and herbicides. The special interest constituent perchlorate was not detected at high RCs, but was detected at moderate RCs in approximately 2 percent of the primary aquifer system. The understanding assessment relied on statistical tests to evaluate relations between concentrations of constituents and values of potential explanatory factors representing geology, land use, well construction, hydrologic conditions, groundwater age, and geochemical conditions. The majority of the high and moderate RCs of arsenic, boron, molybdenum, uranium, and total dissolved solids were in samples from the Honey Lake Valley groundwater basin study area. Groundwater mixing with hydrothermal fluids present in the study area, evaporative concentration of groundwater in the Honey Lake playa, presence of uranium-bearing sediment derived from the adjacent Sierra Nevada, and release of arsenic and other trace elements from sediments under high pH and low dissolved oxygen conditions all appeared to contribute to these elevated concentrations. Thermal springs are in many parts of the Cascade Range and Modoc Plateau study unit and could account for locally elevated concentrations of arsenic, boron, molybdenum, and total dissolved solids in samples from the other study areas. Vanadium concentrations were greater in oxic samples than in anoxic samples, but were not correlated with pH, contrary to expectations from previous studies. Organic constituents were not detected at high or moderate RCs, and the occurrence of low organic constituents at low RCs ranged from 27 percent to 73 percent of the primary aquifers system in the six study areas. The Shasta Valley and Mount Shasta Volcanic study area had significantly greater occurrence of low RCs of herbicides compared to all of the other study areas, which could reflect the greater prevalence of modern groundwater in the Shasta Valley and Mount Shasta Volcanic study area and the presence of potential sources of herbicides, including applications to timberlands and roadside rights-of-way. The Eastside Sacramento Valley study area had the greatest occurrence of low concentrations of chloroform, and chloroform occurrence was most strongly associated with the combination of septic-tank density greater than two tanks per square kilometer and urban land use greater than 10 percent within a radius of 500 meters of the well. These conditions were most prevalent in the Eastside Sacramento Valley study area. The detection frequency of low concentrations of perchlorate was consistent with the probability of occurrence expected under natural conditions, except in the Eastside Sacramento Valley study area, where detection frequencies were much higher than expected and could not be explained by known anthropogenic sources of perchlorate.

California↗
USGS2017

Simulated groundwater flow paths, travel time, and advective transport of nitrogen in the Kirkwood-Cohansey aquifer system, Barnegat Bay–Little Egg Harbor Watershed, New Jersey

Elevated concentrations of nitrogen in groundwater that discharges to surface-water bodies can degrade surface-water quality and habitats in the New Jersey Coastal Plain. An analysis of groundwater flow in the Kirkwood-Cohansey aquifer system and deeper confined aquifers that underlie the Barnegat Bay–Little Egg Harbor (BB-LEH) watershed and estuary was conducted by using groundwater-flow simulation, in conjunction with a particle-tracking routine, to provide estimates of groundwater flow paths and travel times to streams and the BB-LEH estuary. Water-quality data from the Ambient Groundwater Quality Monitoring Network, a long-term monitoring network of wells distributed throughout New Jersey, were used to estimate the initial nitrogen concentration in recharge for five different land-use classes—agricultural cropland or pasture, agricultural orchard or vineyard, urban non-residential, urban residential, and undeveloped. Land use at the point of recharge within the watershed was determined using a geographic information system (GIS). Flow path starting locations were plotted on land-use maps for 1930, 1973, 1986, 1997, and 2002. Information on the land use at the time and location of recharge, time of travel to the discharge location, and the point of discharge were determined for each simulated flow path. Particle-tracking analysis provided the link from the point of recharge, along the particle flow path, to the point of discharge, and the particle travel time. The travel time of each simulated particle established the recharge year. Land use during the year of recharge was used to define the nitrogen concentration associated with each flow path. The recharge-weighted average nitrogen concentration for all flow paths that discharge to the Toms River upstream from streamflow-gaging station 01408500 or to the BB-LEH estuary was calculated. Groundwater input into the Barnegat Bay–Little Egg Harbor estuary from two main sources— indirect discharge from base flow to streams that eventually flow into the bay and groundwater discharge directly into the estuary and adjoining coastal wetlands— is summarized by quantity, travel time, and estimated nitrogen concentration. Simulated average groundwater discharge to streams in the watershed that flow into the BB-LEH estuary is approximately 400 million gallons per day. Particle-tracking results indicate that the travel time of 56 percent of this discharge is less than 7 years. Fourteen percent of the groundwater discharge to the streams in the BB-LEH watershed has a travel time of less than 7 years and originates in urban land. Analysis of flow-path simulations indicate that approximately 13 percent of the total groundwater flow through the study area discharges directly to the estuary and adjoining coastal wetlands (approximately 64 million gallons per day). The travel time of 19 percent of this discharge is less than 7 years. Ten percent of this discharge (1 percent of the total groundwater flow through the study area) originates in urban areas and has a travel time of less than 7 years. Groundwater that discharges to the streams that flow into the BB-LEH, in general, has shorter travel times, and a higher percentage of it originates in urban areas than does direct groundwater discharge to the Barnegat Bay–Little Egg Harbor estuary. The simulated average nitrogen concentration in groundwater that discharges to the Toms River, upstream from streamflow-gaging station 01408500 was computed and compared to summary concentrations determined from analysis of multiple surface-water samples. The nitrogen concentration in groundwater that discharges directly to the estuary and adjoining coastal wetlands is a current data gap. The particle tracking methodology used in this study provides an estimate of this concentration."

New Jersey↗
USGS2017

Status and understanding of groundwater quality in the North San Francisco Bay Shallow Aquifer study unit, 2012; California GAMA Priority Basin Project (ver. 1.1, February 2018)

Groundwater quality in the North San Francisco Bay Shallow Aquifer study unit (NSF-SA) was investigated as part of the Priority Basin Project of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program. The study unit is in Marin, Mendocino, Napa, Solano, and Sonoma Counties and included two physiographic study areas: the Valleys and Plains area and the surrounding Highlands area. The NSF-SA focused on groundwater resources used for domestic drinking water supply, which generally correspond to shallower parts of aquifer systems than that of groundwater resources used for public drinking water supply in the same area. The assessments characterized the quality of untreated groundwater, not the quality of drinking water. This study included three components: (1) a status assessment , which characterized the status of the quality of the groundwater resources used for domestic supply for 2012; (2) an understanding assessment , which evaluated the natural and human factors potentially affecting water quality in those resources; and (3) a comparison between the groundwater resources used for domestic supply and those used for public supply. The status assessment was based on data collected from 71 sites sampled by the U.S. Geological Survey for the GAMA Priority Basin Project in 2012. To provide context, concentrations of constituents measured in groundwater were compared to U.S. Environmental Protection Agency (EPA) and California State Water Resources Control Board Division of Drinking Water regulatory and non-regulatory benchmarks for drinking-water quality. The status assessment used a grid-based method to estimate the proportion of the groundwater resources that has concentrations of water-quality constituents approaching or above benchmark concentrations. This method provides statistically unbiased results at the study-area scale and permits comparisons to other GAMA Priority Basin Project study areas. In the NSF-SA study unit as a whole, inorganic constituents with human-health benchmarks were detected at high relative concentrations (RCs) in 27 percent of the shallow aquifer system, and inorganic constituents with secondary maximum contaminant levels (SMCL) were detected at high RCs in 24 percent of the system. The inorganic constituents detected at high RCs were arsenic, boron, fluoride, manganese, nitrate, iron, sulfate, and total dissolved solids (TDS). Organic constituents with human-health benchmarks were detected at high RCs in 1 percent of the shallow aquifer system. Of the 148 organic constituents analyzed, 30 constituents were detected, although only 1, chloroform, had a detection frequency greater than 10 percent. Natural and anthropogenic factors that could affect the groundwater quality were evaluated by using results from statistical testing of associations between constituent concentrations and values of potential explanatory factors. Groundwater age class (modern, mixed, or pre-modern), redox class (oxic or anoxic), aquifer lithology class (metamorphic, sedimentary, or volcanic), and dissolved oxygen concentrations were the explanatory factors that explained distribution patterns of most of the inorganic constituents best. Groundwater classified primarily as pre-modern or mixed in age was associated with higher concentrations of arsenic and fluoride than waters classified as modern. Anoxic or mixed redox conditions were associated with higher concentrations of boron, fluoride, and manganese. Similar patterns of association with explanatory variables were seen for inorganic constituents with aesthetic-based benchmarks detected at high concentrations. Nitrate and perchlorate had higher concentrations in oxic than in the anoxic redox class and were positively correlated with urban land use. The NSF-SA water-quality results were compared to those of the GAMA North San Francisco Bay Public-Supply Aquifer study unit (NSF-PA). The NSF-PA was sampled in 2004 and covers much of the same area as the NSF-SA, but focused on the deeper public-supply aquifer system. The comparison of the NSF-PA to the NSF-SA showed that there were more differences between the Valleys and Plains study areas of the two study units than between the Highlands study areas of the two study units. As expected from the shallower depth of wells, the NSF-SA Valleys and Plains study area had a lesser proportion of pre-modern age groundwater and greater proportion of modern age groundwater than the NSF-PA Valleys and Plains study area. In contrast, well depths and groundwater ages were not significantly different between the two Highlands study areas. Arsenic, manganese, and nitrate were present at high RCs, and perchlorate was detected in greater proportions of the NSF-SA Valleys and Plains study area than the NSF-PA Valleys and Plains study area.

California↗
USGS2017

Simulation of groundwater and surface-water flow in the upper Deschutes Basin, Oregon

This report describes a hydrologic model for the upper Deschutes Basin in central Oregon developed using the U.S. Geological Survey (USGS) integrated Groundwater and Surface-Water Flow model (GSFLOW). The upper Deschutes Basin, which drains much of the eastern side of the Cascade Range in Oregon, is underlain by large areas of permeable volcanic rock. That permeability, in combination with the large annual precipitation at high elevations, results in a substantial regional aquifer system and a stream system that is heavily groundwater dominated. The upper Deschutes Basin is also an area of expanding population and increasing water demand for public supply and agriculture. Surface water was largely developed for agricultural use by the mid-20th century, and is closed to additional appropriations. Consequently, water users look to groundwater to satisfy the growing demand. The well‑documented connection between groundwater and the stream system, and the institutional and legal restrictions on streamflow depletion by wells, resulted in the Oregon Water Resources Department (OWRD) instituting a process whereby additional groundwater pumping can be permitted only if the effects to streams are mitigated, for example, by reducing permitted surface-water diversions. Implementing such a program requires understanding of the spatial and temporal distribution of effects to streams from groundwater pumping. A groundwater model developed in the early 2000s by the USGS and OWRD has been used to provide insights into the distribution of streamflow depletion by wells, but lacks spatial resolution in sensitive headwaters and spring areas. The integrated model developed for this project, based largely on the earlier model, has a much finer grid spacing allowing resolution of sensitive headwater streams and important spring areas, and simulates a more complete set of surface processes as well as runoff and groundwater flow. In addition, the integrated model includes improved representation of subsurface geology and explicitly simulates the effects of hydrologically important fault zones not included in the previous model. The upper Deschutes Basin GSFLOW model was calibrated using an iterative trial and error approach using measured water-level elevations (water levels) from 800 wells, 144 of which have time series of 10 or more measurements. Streamflow was calibrated using data from 21 gage locations. At 14 locations where measured flows are heavily influenced by reservoir operations and irrigation diversions, so called “ naturalized ” flows, with the effects of reservoirs and diversion removed, developed by the Bureau of Reclamation, were used for calibration. Surface energy and moisture processes such as solar radiation, snow accumulation and melting, and evapotranspiration were calibrated using national datasets as well as data from long-term measurement sites in the basin. The calibrated Deschutes GSFLOW model requires daily precipitation, minimum and maximum air temperature data, and monthly data describing groundwater pumping and artificial recharge from leaking irrigation canals (which are a significant source of groundwater recharge). The calibrated model simulates the geographic distribution of hydraulic head over the 5,000 ft range measured in the basin, with a median absolute residual of about 53 ft. Temporal variations in head resulting from climate cycles, pumping, and canal leakage are well simulated over the model area. Simulated daily streamflow matches gaged flows or calculated naturalized flows for streams including the Crooked and Metolius Rivers, and lower parts of the mainstem Deschutes River. Seasonal patterns of runoff are less well fit in some upper basin streams. Annual water balances of streamflow are good over most of the model domain. Model fit and overall capabilities are appropriate for the objectives of the project. The integrated model results confirm findings from other studies and models indicating that most streamflow in the upper Deschutes Basin comes directly from groundwater discharge. The integrated model provides additional insights about the components of streamflow including direct groundwater discharge to streams, interflow, groundwater discharge to the land surface (Dunnian flow), and direct runoff (Hortonian flow). The new model provides improved capability for exploring the timing and distribution of streamflow capture by wells, and the hydrologic response to changes in other external stresses such as canal operation, irrigation, and drought. Because the model uses basic meteorological data as the primary input; and simulates surface energy and moisture balances, groundwater recharge and flow, and all components of streamflow; it is well suited for exploring the hydrologic response to climate change, although no such simulations are included in this report. The model was developed as a tool for future application; however, example simulations are provided in this report. In the example simulations, the model is used to explore the influence of well location and geologic structure on stream capture by pumping wells. Wells were simulated at three locations within a 12-mi area close to known groundwater discharge areas and crossed by a regional fault zone. Simulations indicate that the magnitude and timing of stream capture from pumping is largely controlled by the geographic location of the wells, but that faults can have a large influence on the propagation of pumping stresses.

Oregon↗
USGS2019

Evaluation of groundwater resources in the Spanish Valley Watershed, Grand and San Juan Counties, Utah

Groundwater resources in the Spanish Valley watershed in southern Utah were quantified for the first time since the early 1970s. The primary objectives of this study were (1) to better understand sources of recharge to, groundwater flow directions within, and discharge points for both the valley-fill and Glen Canyon Group aquifers (VFA and GCGA), and (2) to quantify groundwater budget components of the combined VFA and GCGA, including both recharge and discharge. Based on both groundwater chemistry (stable isotopes, major ions, and noble gases) and environmental tracers in vadose-zone pore water of the Navajo Sandstone outcrop along Sand Flats Road, most recharge to the GCGA occurs high in the La Sal Mountains, and not on the sandstone outcrop area. The same groundwater chemistry and environmental tracer evidence from the saturated zone indicates that Pack Creek, rather than GCGA groundwater, is the primary source of recharge to the VFA. Groundwater recharge in the study area occurs mostly from infiltration of precipitation (in the form of snowmelt) at high altitudes. Additional recharge occurs from the infiltration of runoff along losing reaches of stream channels, or as unconsumed surface-water and groundwater irrigation. Average annual recharge to the Moab-Spanish Valley watershed part of the Spanish Valley study area was estimated to be between 9,550 and 30,000 acre-feet. Based on water-levels collected in the current study, groundwater in both the GCGA and the VFA generally moves downgradient parallel to the topographic slope of the watershed towards the Colorado River. Groundwater discharge measurements, and hydraulic-flux estimates at the lower end of Spanish Valley, provide a more robust estimate of the groundwater budget than evaluating recharge. The primary base-flow discharge components in the study area include groundwater discharge to gaining reaches of streams, groundwater discharge to springs, and well withdrawals. Based on 3 years of measurements (2014–16) and hydraulic-flux calculations at the lower end of Spanish Valley, total groundwater discharge was estimated to be 14,000 to 16,000 acre-feet per year (acre-ft/yr) for the entire watershed, or 13,000 to 15,000 acre-ft/yr, excluding the watershed areas of Grandstaff (formerly Negro Bill) and Ice Box Canyons (compared to the 1971 Sumsion estimate of 22,000 acre-ft/yr). The primary difference is this study’s estimate of subsurface outflow to the Colorado River of only 300 to 1,000 acre-ft/yr, compared to 11,000 acre-ft/yr estimated by Sumsion. Because the study period (2014–16) experienced above average precipitation for 2 of the 3 years, the discharge estimates may be slightly higher than long-term average annual discharge from the groundwater system.

Utah↗
USGS2020

Numerical simulation of groundwater availability in central Moloka‘i, Hawai‘i

Since the 1990s, increased chloride concentrations of water pumped from wells (much of which is used for drinking water) and the effects of withdrawals on groundwater-dependent ecosystems have led to concerns over groundwater availability on the island of Molokaʻi, Hawaiʻi. An improved understanding of the hydrologic effects of proposed groundwater withdrawals is needed to ensure effective management of the groundwater resources of Molokaʻi, plan for possible growth, and accommodate cultural, social, and economic concerns. To address the information needs of managers and community stakeholders on Molokaʻi, the U.S. Geological Survey developed a numerical groundwater model capable of simulating salinity change and reduction in groundwater discharge in coastal areas of central and southern Molokaʻi. Estimates of groundwater recharge needed as input to the numerical groundwater model were made using a daily water budget for each decade during 1940−2012 (the period 2000−12 spanned 13 years) and the most current available data, including the distributions of monthly rainfall and potential evapotranspiration. Total island recharge during the decadal periods ranged from a low of about 189 Mgal/d during the 1970s to a high of 278 Mgal/d during the 1960s. These recharge estimates were used to develop an island-wide numerical groundwater model with simplifying assumptions (sharp interface between freshwater and saltwater; two-dimensional flow). The island-wide model provided estimates of groundwater inflows to the main area of interest simulated with a three-dimensional numerical groundwater model. Simulated withdrawal scenarios were selected in consultation with water managers and stakeholders and consisted of: (1) a baseline scenario using average recharge (1978−2007 rainfall and 2010 land cover) and average 2016−17 withdrawals; (2) a scenario using average recharge and withdrawals from existing wells at pending (as of January 2019) water-use permit rates; (3) six scenarios using average recharge and selected withdrawals from existing and proposed wells; and (4) a scenario using reduced recharge and selected withdrawals from existing and proposed wells. Results of the simulated withdrawal scenarios indicate that wells may be capable of producing groundwater with chloride concentrations below 250 mg/L at withdrawal rates exceeding average 2016−17 rates. However, the quality of water withdrawn from production wells is dependent on the rate and distribution of the withdrawals. For all nonbaseline scenarios, simulated groundwater discharge to the nearshore environment is reduced relative to the baseline scenario. Areas of discharge reduction may correspond to areas used for cultural or subsistence purposes. The three-dimensional numerical groundwater model developed for this study utilizes the latest available hydrologic and geologic information and is a useful tool for understanding the hydrologic effects of additional groundwater withdrawals in central Molokaʻi. The model has several limitations, including its nonuniqueness and inability to account for local-scale heterogeneities.

Hawaii↗
USGS2020

Geochemical assessment of groundwater in the Big Chino subbasin, Arizona, 2011–18

A geochemical characterization of groundwater in the Big Chino subbasin of Arizona was conducted by the U.S. Geological Survey, in cooperation with the City of Prescott, the Town of Prescott Valley, and the Salt River Project, to understand groundwater evolution through the study area and the source of water to springs along the gaining reach of the Verde River just downstream from its confluence with Granite Creek. Samples were collected between 2011 and 2018 in groundwater wells completed in basin-fill and carbonate aquifers and at selected springs, including two discrete springs discharging along the aforementioned stretch of the Verde River. Five newly installed monitoring wells completed in the carbonate aquifer were sampled in 2018. Water-quality results obtained from these samples include the first known geochemical data for carbonate groundwater beneath the basin-fill in the Big Chino subbasin downgradient from Walnut Creek near Paulden, Arizona, as well as other parts of the study area without previous data. Groundwater samples were collected and analyzed for major ions, arsenic, nutrients, stable isotopes of oxygen and hydrogen (δ 18 O and δ 2 H), strontium isotopes ( 87 Sr/ 86 Sr), carbon-14, isotopes of carbon (δ 13 C), and noble gases. Significant differences in groundwater geochemistry between the basin-fill and carbonate aquifers were driven primarily by higher pH, tritium, and δ 18 O and δ 2 H in the basin-fill aquifer samples and higher specific conductance and higher concentrations of calcium, sodium, bicarbonate, fluoride, and arsenic in the carbonate aquifer samples. All but one sample from the carbonate aquifer and two samples from the basin-fill aquifer exceeded the U.S. Environmental Protection Agency (EPA) drinking water standard for arsenic of 10 micrograms per liter. One basin-fill aquifer sample exceeded the EPA drinking water standard for fluoride of 4 milligrams per liter, and one carbonate aquifer sample exceeded the EPA secondary drinking water standard for fluoride of 2 milligrams per liter. A component of modern groundwater recharged following aboveground nuclear testing beginning in the mid-1950s is present in some basin-fill and spring groundwater from this study. Groundwater that can be dated using radiocarbon decay is also present in the study area, with four groundwater samples indicating possible recharge during the Pleistocene with groundwater ages ranging from approximately 34,600 to 13,300 years before present. Other groundwater sampled during this study that can dated using radiocarbon decay ranged in age from about 7,500 to 1,100 years before present, indicating possible recharge during the Holocene. The gaining reach of the Verde River downstream from the confluence with Granite Creek shows areal changes in temperature, pH, and specific conductance, indicating multiple zones of groundwater input. Surface-water samples for analyses of δ 18 O and δ 2 H have been collected at the Verde River near Paulden, Ariz. streamgage (09503700) during discharge measurements since 2009, and a trend analysis of the δ 18 O and δ 2 H data indicated no significant trend exists for the 10-year period of record. Additional groundwater samples from the carbonate aquifer beneath the basin-fill upgradient and downgradient from Walnut Creek would provide valuable information to understand groundwater evolution along the Big Chino subbasin.

Arizona↗
USGS2021

Status and understanding of groundwater quality in the northern Sierra Nevada foothills domestic-supply aquifer study units, 2015–17—California GAMA Priority Basin Project

Groundwater quality in the northern Sierra Nevada foothills region of California was investigated as part of California State Water Resources Control Board (SWRCB) Groundwater Ambient Monitoring Assessment Priority Basin Project (GAMA-PBP). The region was divided into two study units: the Yuba-Bear watersheds (YBW) study unit and the American-Cosumnes-Mokelumne watersheds (ACMW) study unit. The GAMA-PBP made a spatially unbiased assessment of aquifer systems used for domestic drinking-water supply in the study region, which are predominantly composed of fractured, hard-rock aquifers of varying lithology. These assessments characterized the quality of raw groundwater to evaluate ambient conditions in the domestic-supply aquifer and not the quality of treated drinking water. The study included three components: (1) a status assessment , which characterized the quality of groundwater resources used for domestic drinking-water supply in the YBW and ACMW study units; (2) an understanding assessment , which evaluated natural and anthropogenic explanatory factors that could potentially affect groundwater quality in the study region; and (3) a comparative assessment between the groundwater resources used for domestic and public drinking-water supply in the study region. The status assessment was based on data collected by the GAMA-PBP from 74 sites in the YBW study unit during 2015–16 and 67 sites in the ACMW study unit from 2016 to 2017. To contextualize water-quality results, concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. The status assessment used a grid-based method to estimate proportions of groundwater resources with concentrations approaching or exceeding benchmark thresholds. This method provides spatially unbiased results and allows inter-comparability with similar groundwater-quality assessments. Inorganic constituents with health-based benchmarks were present at high relative concentration (RC), meaning they exceeded the benchmark threshold, in 5.4 and 10 percent of domestic-supply aquifer systems in the YBW and ACMW study units, respectively. Inorganic constituents with aesthetic-based benchmarks were detected at high-RCs in 20 and 28 percent of the YBW and ACMW study units, respectively. The inorganic constituents present at high RC were arsenic, barium, boron, molybdenum, strontium, nitrate, adjusted gross-alpha particle activity, chloride, total dissolved solids, specific conductance, iron, manganese, and hardness. Groundwater samples were tested for presence or absence of three microbial indicators (total coliform, Escherichia coli , and Enterococci ). At least one microbial indicator was present in 26 and 28 percent of the YBW and ACMW study units, respectively. At least one organic constituent was detected in 30 and 42 percent of the YBW and ACMW study units, respectively. Organic constituents were not present at high RC, but tetrachloroethene (PCE), trichloroethene (TCE), and toluene were detected in the YBW study unit at moderate RC (between the benchmark concentration and one-tenth of the benchmark concentration). Methyl tert -butyl ether (MTBE) and chloroform were present at low RC (less than one-tenth of the benchmark concentration) in the YBW and ACMW study units with detection frequencies greater than 10 percent. Perchlorate, a constituent of special interest, was detected in 31 and 41 percent of the YBW and ACMW study units, respectively, at either low or moderate RCs. Relations among select water-quality constituents and potential explanatory factors were evaluated using statistical and graphical approaches. Nitrate, microbial indicators, and perchlorate were all correlated to elevation-dependent variables relating to climate, land use, and recharge condition. Isotopic and dissolved noble-gas tracers indicated these water-quality constituents are associated with recharge conditions associated with irrigation during the summer dry-season, which is common in areas of rural-residential or agricultural land uses. Higher concentrations of iron and manganese were primarily associated with anoxic groundwater in aquifers of metasedimentary lithology. Increased hardness was primarily associated with anoxic groundwater in aquifers of mafic-ultramafic or metavolcanics lithologies at lower elevations in the study region in the Melones fault zone. Chloroform and MTBE were associated with shallow groundwater (wells depths less than 130 m) under oxic and anoxic redox conditions, respectively. The comparative assessment evaluated differences between the aquifer systems used for domestic- and public-supply in study region based on (1) well-construction characteristics, and (2) water quality. Analysis of over 60,000 well-completion reports in the study region showed that although domestic-supply wells span the deepest depth zones in regional aquifers, median depths for public-supply wells were significantly greater than those of domestic-supply wells in both study units. Water-quality data from more than 300 public-supply wells in the study region were assessed using a spatially weighted method for calculation aquifer-scale proportions and compared with the domestic-supply assessment results. Detections of inorganic constituents at high RC and detection frequencies for organic constituents were generally similar between the domestic- and public-supply aquifer systems in both study units, with a few notable exceptions in the ACMW study unit: nitrate was greater for the public- compared to domestic-supply aquifer system and both manganese, hardness, and MTBE were greater in the domestic- compared to public-supply aquifer system. These differences are likely related to contrasting land uses, aquifer lithologies, landscape positions, and depths characterizing domestic- and public-supply wells in the ACMW study unit. Overall, fewer samples from domestic-supply wells in the northern Sierra Nevada foothills exceeded health-based benchmarks compared to aesthetic-based benchmarks for groundwater quality. Exceedences of health-based benchmarks were primarily caused by nitrate and coliform bacteria, which were associated with recharge from diverted surface water used primarily for irrigation. Exceedences of aesthetic-based benchmarks were primarily caused by iron, managanese, and hardness, which were associated with geologic factors. Regional irrigation practices and aquifer lithology can affect groundwater quality in fractured-rock aquifers in the northern Sierra Nevada foothills used for domestic drinking-water supply.

California↗
USGS2021

Groundwater hydrology and chemistry of Jamestown Island, Virginia—Potential effects of tides, storm surges, and sea-level rise on archaeological, cultural, and ecological resources

As the site of the first permanent English settlement in North America in 1607, Jamestown Island, Colonial National Historical Park (COLO), Virginia, contains a rich archaeological record that extends from the Paleoindian period (15,000 to 8,000 years ago) through the 20th century. The island is located on the lower James River near the mouth of Chesapeake Bay. Jamestown Island vegetation is dominated by upland forests surrounded by tidal, freshwater-to-oligohaline marshes. Along the Virginia coast, relative sea-level rise was more than 2.5 times the global average during the 20th century. Consequently, the National Park Service (NPS) has identified COLO as one of the 25 national parks most threatened by climate change. Surface waters across the island are hydraulically connected to the laterally continuous Surficial aquifer. The land-surface altitude of the island is low, with two-thirds of the island less than 5 feet (ft) above the North American Vertical Datum of 1988 (NAVD 88). Consequently, sea-level rise, combined with tides and storm surges, threatens the island and its resources as surface-water and groundwater levels rise, saltwater enters the Surficial aquifer, and groundwater chemistry changes. The impact of sea-level rise on the island’s surface-water resources has been well studied, but groundwater effects have been largely ignored. Quantifying the effects of tides, storm surges, and sea-level rise on groundwater levels and chemistry is essential to developing an effective strategy for managing climate-induced changes. The first step in developing a response strategy includes a parkwide general risk assessment for archaeological sites on the island, so that sites can be prioritized for management actions. The U.S. Geological Survey and the NPS began a study in 2015 to develop a long-term groundwater-monitoring program to evaluate this risk and to develop an updated management strategy. The groundwater-monitoring program consists of 45 wells and piezometers in two individual clusters and three transects across the island in different hydrologic and chemical settings. Samples for water quality were collected from the wells and piezometers from October 2015 through September 2018 at variable time intervals. Results of the monitoring identified disparate hydrologic and chemical responses to saltwater intrusion across the island. Specific conductance (an indicator of salinity) of groundwater beneath several marshes responded differently to changes in James River salinity. Groundwater response to changes in James River specific conductance appeared to be controlled by land-surface altitude and slope, differences in lateral and vertical sediment characteristics, distance from surface waters, and the degree of surface water/groundwater connectivity between channels and the aquifer. Groundwater chemistry data from monitoring wells at Black Point, a low-altitude, upland setting, are in contrast with conditions observed in Island House observation wells, a high-altitude, upland setting. Specific conductance (less than 200 microsiemens per centimeter [μS/cm]) and pH (greater than 5.0) of groundwater beneath much of the uplands that characterize the Island House observation wells are typical of groundwater in noncarbonate sedimentary aquifers recharged by precipitation. At Black Point, specific conductance ranged from 2,490 to 15,200 μS/cm, and pH ranged from 3.1 to 6.6 standard units. At the Black Point observation wells, the most saline and dense water was at the water table rather than deeper in the aquifer, causing a density inversion that persisted throughout the study. The density inversion likely resulted from differences in permeability between the shallow clay and fine-grained sands and the deeper coarse-grained sand and gravel. Groundwater with the lowest pH was at the water table. As saline groundwater flows through organic sediment beneath the marshes, bacterial biodegradation of organic matter creates anoxic conditions. Continued biodegradation concomitantly reduces iron-oxide minerals in the sediment and sulfate in saline water. When oxygen is reintroduced into groundwater, iron and sulfur can reoxidize to form sulfuric acid, locally lowering the pH of the water. This report describes the groundwater monitoring network design, rationale for site selection, monitoring approach, and results of monitoring from October 2015 through September 2018. Maps of inundation at selected water-level altitudes are included to identify the risk to archaeological, cultural, and ecological resources. The monitoring results of the hydrology and chemistry data are interpreted, and the different hydrologic and chemical settings are described. The implications of the study results for management decisions are presented, and suggestions for improving the monitoring network are included.

Virginia↗
USGS2023

Updates to the regional groundwater-flow model of the New Jersey Coastal Plain, 1980–2013

A 21-layer three-dimensional transient groundwater-flow model of the New Jersey Coastal Plain was developed and calibrated by the U.S. Geological Survey (USGS) in cooperation with the New Jersey Department of Environmental Protection to simulate groundwater-flow conditions during 1980–2013, incorporating average annual groundwater withdrawals and average annual groundwater recharge. This model is the third version of the New Jersey Coastal Plain regional groundwater-flow model that was initially developed as part of the USGS Regional Aquifer System Analysis (RASA) program. The model simulates groundwater flow in 11 aquifers and 10 intervening confining units of the New Jersey Coastal Plain to provide a regional overview of groundwater conditions. Averaged groundwater withdrawal data for 1980 to 2013 were used in the model. The 11 aquifers in New Jersey are, from shallowest to deepest, the Holly Beach water-bearing zone and the confined Cohansey aquifer in Cape May County; the Rio Grande water-bearing zone; the Atlantic City 800-foot sand; the Piney Point, Vincentown, and Wenonah-Mount Laurel aquifers; the Englishtown aquifer system; and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy (PRM) aquifer system. The model was developed with the MODFLOW–2005 numerical code and the UCODE parameter estimation technique and calibrated using water-level and base-flow observations. A total of 3,453 water-level observations from 392 wells in New Jersey and 48 wells in Delaware from 1983 to 2013 were used in model calibration, which includes historical water-level trends for 29 wells in New Jersey during 1980–2013 presented in time-series hydrographs. In addition, derived observations also were included by calculating the vertical gradient at 33 pairs of nested observation wells in New Jersey, for a total of 210 observations. Changes in water levels over time were calculated for 134 wells in New Jersey and four wells in Delaware where water levels had varied substantially (approximately 10 ft) over the 30-year span of synoptic water-level measurements, for a total of 767 observations. A total of 1,485 base-flow observations in 47 surface-water basins in New Jersey from 1980 to 2013 were used in model calibration. Updates to the groundwater-flow model include the conversion to a fully three-dimensional model from the previous quasi-three-dimensional model. The new model will allow for potential future uses such as particle tracking or simulation of variable-density groundwater flow that could not be accomplished with earlier versions of the model. Spatially and temporally variable recharge estimated by using a soil-water balance model resulted in a spatially and temporally finer discretization. The Rio Grande water-bearing zone was added to the model as an aquifer layer to refine estimates of simulated flow in Atlantic and Cape May Counties, New Jersey. Hydrogeologic parameters were updated to include the confining units in New Jersey and corresponding hydrogeologic units in Delaware and eastern Maryland. The simulated water levels for the New Jersey Coastal Plain aquifers were compared to water-level measurements made during 1980–2013. The average residual for 4,243 water-level observations for New Jersey (simulated water levels minus measured water levels) is 1.5 feet. The simulated water-level contours for the confined aquifers for 2013 were compared to potentiometric surfaces produced from water levels measured during 2013. Simulated water levels generally matched the 2013 potentiometric surfaces of the confined aquifers in the areas of large withdrawals. Hydrographs of wells in the confined Coastal Plain aquifers of New Jersey show that simulated water levels generally match the magnitude and seasonal variation of the observed water levels. Hydrographs of base flow for the 47 streamgaging stations in New Jersey indicate that most of the simulated and estimated data match reasonably well. Groundwater withdrawals are an important resource for water supply, agricultural, industrial, and commercial needs in the New Jersey Coastal Plain. Groundwater withdrawals from the New Jersey Coastal Plain aquifers have resulted in persistent, regionally extensive cones of depression in the Englishtown aquifer system and Wenonah-Mount Laurel aquifer in Ocean and Monmouth Counties; Wenonah-Mount Laurel and upper, middle, and lower PRM aquifers in Camden County; and Atlantic City 800-foot sand in Atlantic County. Because hydrologic stresses and water-management needs change with time, periodic updates to the groundwater-flow model are required to provide current information about hydrologic conditions in the New Jersey Coastal Plain and to maintain its usefulness as a tool to manage water resources and develop water-resource strategies. The current updates will support the continued application of this model as a tool for evaluating the regional effects of changes in groundwater withdrawals and of current and potential future water-management strategies on groundwater levels in the New Jersey Coastal Plain.

New Jersey Coastal Plain↗
USGS2024

An update of hydrologic conditions and distribution of selected constituents in water, eastern Snake River aquifer and perched groundwater zones, Idaho National Laboratory, Idaho, emphasis 2019–21

Since 1952, wastewater discharged to infiltration 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 (DOE), maintains groundwater-monitoring networks at the INL to determine hydrologic trends and to delineate the movement of radiochemical and chemical wastes in both the aquifer and perched groundwater zones. This report presents an analysis of water-level and water-quality data collected from the ESRP aquifer and perched groundwater wells from the USGS groundwater monitoring networks during 2019–21. From March–May 2018 to March–May 2021, water levels in wells completed in the ESRP aquifer increased in the northern part of the INL and decreased in the southwestern part. Water-level increases ranged from 0.02 to 1.04 feet in the northern part and decreases ranged from 0.03 to 2.94 feet in the southwestern part of the INL. Detectable concentrations of radiochemical constituents in water samples from wells in the ESRP aquifer at the INL generally decreased or remained constant during 2019–21. Decreases in concentrations were attributed to radioactive decay, changes in waste-disposal methods, and dilution from recharge and underflow. In 2021, tritium was detected above reporting levels in water samples collected from 46 of 105 aquifer wells and ranged from 150±50 to 4,280±150 picocuries per liter (pCi/L). Tritium concentrations from eight wells completed in deep perched groundwater near the Advanced Test Reactor Complex (ATRC) generally were greater than or equal to the reporting level during at least one sampling event during 2019–21, and concentrations ranged from 160±50 to 2,097±107 pCi/L. Concentrations of strontium-90 in water from 12 of 45 aquifer wells sampled in 2021 exceeded the reporting level, and concentrations ranged from 2.5±0.7 to 299±6 pCi/L. During 2021, concentrations of strontium-90 from five wells completed in deep perched groundwater at the ATRC equaled or exceeded the reporting levels, and concentrations ranged from 3±0.9 pCi/L to 27.8±1.3 pCi/L. Concentrations of cesium-137 were less than the reporting level in all but one aquifer well, and concentrations of plutonium-238, plutonium-239, -240 (undivided), and americium-241 were less than the reporting level in water samples from all aquifer wells sampled during this study period. Dissolved chromium concentrations in water samples from 64 ESRP aquifer wells ranged from less than (<) 0.5 to 76.4 micrograms per liter (μg/L). During 2019–21, dissolved chromium was detected in water from wells completed in deep perched groundwater above the ESRP aquifer at the ATRC, and concentrations ranged from <1 to 82.1 μg/L. In 2021, concentrations of dissolved sodium in water from most ESRP aquifer wells in the southern part of the INL were greater than the western tributary groundwater background concentration of 8.3 milligrams per liter (mg/L). During 2021, dissolved sodium concentrations in water from 15 wells completed in deep perched groundwater ranged from 11.7 to 122.5 mg/L. Variations in sodium concentrations in aquifer wells and perched groundwater zones are attributed to either migration of remnant water from the former chemical-waste ponds or disposal volume and composition variability in percolation ponds installed in 2008. In 2021, concentrations of chloride in most water samples from ESRP aquifer wells south of the Idaho Nuclear Technology and Engineering Center (INTEC) and at the Central Facilities Area (CFA) exceeded background concentrations. Chloride concentrations in water from wells south of the INTEC have generally decreased because of discontinued chloride disposal to the legacy percolation ponds since 2002 when the discharge of wastewater was discontinued. During 2019–21, dissolved chloride concentrations in deep perched groundwater above the ESRP aquifer from 18 wells at the ATRC ranged from 8.15 to 231 mg/L. In 2021, sulfate concentrations in water samples from ESRP aquifer wells in the south-central part of the INL that exceeded the background concentration of sulfate, ranged from 21 to 141 mg/L. The greater-than-background concentrations in water from these wells are attributed to sulfate disposal at the ATRC infiltration ponds or the legacy INTEC percolation ponds. In 2021, sulfate concentrations in water samples from aquifer wells near the Radioactive Waste Management Complex (RWMC) were mostly greater than background concentrations. The maximum dissolved sulfate concentration in shallow perched groundwater near the ATRC was 575 mg/L in 2021. During 2021, dissolved sulfate concentrations in water from wells completed in deep perched groundwater near the cold waste ponds at the ATRC ranged from 22.3 to 519 mg/L. In 2021, concentrations of nitrate in water from most ESRP aquifer wells at and near the INTEC exceeded the western tributary groundwater background concentration of 0.655 mg/L. Concentrations of nitrate in aquifer wells southwest of INTEC and farther away from the influence of disposal areas and the Big Lost River, in intermittent source of surface water recharge to the aquifer, show a general decrease in nitrate concentration over time. Two aquifer wells south of INTEC show increasing trends that could result from wastewater beneath the INTEC tank farm being mobilized to the aquifer. During 2019–21, water samples from several ESRP aquifer wells were collected and analyzed for volatile organic compounds (VOCs). Twelve VOCs were detected, and 1–4 VOCs were detected in water samples from 10 wells. The most frequently detected VOCs include carbon tetrachloride (tetrachloromethane), trichloromethane, tetrachloroethene, 1,1,1-trichloroethane, and trichloroethene. In 2019–21, concentrations for all VOCs were less than their respective maximum contaminant levels (MCLs) for drinking water, except carbon tetrachloride in one well, trichloroethene in two wells, and vinyl chloride in one well. During 2019–21, variability and bias were evaluated from 34 replicate and 14 blank quality-assurance samples. Results from replicate analyses were investigated to evaluate sample variability. Constituents with acceptable reproducibility were major ions, trace elements, nutrients, and VOCs. All radiochemical constituents including gross alpha- and beta- radioactivity, strontium-90, cesium-137, and tritium, had acceptable reproducibility. Bias from sample contamination was evaluated from equipment, field, and source-solution blanks. Chloride and sulfate were detected slightly above their respective method detection limits in equipment and field blanks, but at concentrations well below the co-collected sample for that well. These chloride and sulfate detections in the field and equipment blanks were inconsequential because they weren’t detected above the analysis-specific variability for those constituents as determined by replicate sample result evaluation. None of the detections of nutrients and trace inorganic constituents were high enough to indicate environmental sample or analytical procedure bias.

Idaho↗
USGS2024

Groundwater model of the Harney Basin, southeastern Oregon

Groundwater development, mainly for large-scale irrigation, has increased substantially in the Harney Basin of southeastern Oregon since 2010. Concurrently, some areas of the basin experienced groundwater-level declines of more than 100 feet, and some shallow wells have gone dry. The Oregon Water Resources Department has limited new groundwater development in the basin until an improved understanding of the groundwater-flow system is available. The groundwater resources report by Gingerich and others (2022, U.S. Geological Survey Scientific Investigations Report 2021–5103, https://doi.org/10.3133/sir20215103 ) provides that understanding. This report describes the development of a numerical groundwater-flow model that can be used as a tool to help improve that understanding. The Harney Basin Groundwater Model was developed using the finite-difference groundwater-modeling software U.S. Geological Survey modular finite-difference groundwater-flow model (MODFLOW 6) and associated Python pre- and post-processing routines. The groundwater model encompasses the entire 5,240-square-mile Harney Basin and adjacent areas and is calibrated to the hydrologic conditions from 1930 to 2018. The model has a uniform grid consisting of 78,064 nearly square cells, each covering 2,005 by 2,007 feet (about 92 acres) and has 10 layers (780,640 total cells) representing the vertical distribution of hydrogeologic units. The results from the calibrated model simulations indicate that groundwater pumpage exceeded recharge since about the mid-1980s, resulting in an estimated net cumulative depletion of groundwater storage (discharge minus recharge) of about 840,000 acre-feet and also indicated declines in groundwater evapotranspiration and spring and stream discharge. Model simulations show as much as 100 feet of groundwater-level decline in some areas and more than 40 feet of decline in widespread areas in recent decades. Model simulations are consistent with field observations of groundwater levels through time.

Oregon↗
USGS2026

Groundwater budget for the Mountain Home area, southern Idaho, 2022–23

The U.S. Geological Survey, with funding from the Idaho Department of Water Resources, developed a groundwater budget for the Mountain Home area in southern Idaho for irrigation year 2023 (November 1, 2022–October 31, 2023). This study focused on the water balance across the Cinder Cone Butte Critical Groundwater Area (CGWA), Mountain Home Groundwater Management Area (GWMA), and the rest of the study area (RoSA), compiling data from various sources, including precipitation records, groundwater level measurements, metered groundwater pumpage data, surface water diversions and evapotranspiration (ET) estimates derived from remote sensing satellite imagery, and ground-based reference data. Key inflow components included recharge from applied surface water irrigation (which incorporates incidental recharge from irrigation practices and conveyance losses), estimated tributary streamflow, and estimated mountain block recharge. The key outflow components were groundwater pumpage for irrigation, municipal, industrial, and domestic uses, and ET. Recharge from applied irrigation and mountain block recharge were the largest inflows, and groundwater pumpage for irrigation was the largest outflow. The CGWA had a positive groundwater budget residual of 2,170 acre-feet (acre-ft), which contrasts with observed long-term groundwater level declines and historical trends of storage depletion. This positive residual is likely associated with unquantified outflows, including lateral groundwater flow out of the subregion, or other complexities, such as overestimated tributary contributions relative to the actual recharge for the 2023 water budget. The GWMA exhibited a positive residual of 56,563 acre-ft, primarily owing to recharge from applied surface water irrigation and areal recharge during a wetter-than-average year, which allowed irrigation entities to deliver more water from in-basin and out-of-basin reservoirs. The RoSA showed a large positive residual of 124,933 acre-ft. The interpretation of these positive residuals must account for significant uncertainties, including estimations of areal recharge, tributary streamflow (particularly losses and diversions), ET, the volume of surface water loss to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and the unquantified groundwater discharge to the Snake River. These uncertainties, in combination with the complex hydrogeologic controls on water movement and limitations of remotely sensed data, directly affect the accuracy of water availability assessments. Future data collection efforts would help reduce these uncertainties and support water resource management decisions in the Mountain Home area. Key efforts could include installing additional streamflow gaging stations (particularly to quantify tributary losses and gains and surface water losses to the Snake River), improving groundwater pumpage metering, and validating remotely sensed ET data with ground-based measurements. Furthermore, to better quantify unrepresented or highly uncertain fluxes, focused investigations on groundwater discharge to the Snake River, lateral groundwater flows between subregions and across study area boundaries, and a more robust determination of the actual influence and volume of mountain block recharge would help refine future water availability assessments for the Mountain Home area.

Idaho↗
USGS2019

Characterizing groundwater/surface-water interaction using hydrograph-separation techniques and groundwater-level data throughout the Mississippi Delta, USA

The Mississippi Delta, located in northwest Mississippi, is an area dense with industrial-level agriculture sustained by groundwater-dependent irrigation supplied by the Mississippi River Valley Alluvial aquifer (alluvial aquifer). The Delta provides agricultural commodities across the United States and around the world. Observed declines in groundwater altitudes and streamflow contemporaneous with increases in irrigation have raised concerns about future groundwater availability and the effects of groundwater withdrawals on streamflow. To quantify the impacts of groundwater withdrawals on streamflow and increase understanding of groundwater and surface-water interaction, hydrograph-separation techniques were used to estimate baseflow and identify statistical streamflow trends. The analysis was conducted using the U.S. Geological Survey Groundwater Toolbox open-source software and daily hydrologic data provided by a spatially-distributed network of paired groundwater wells and streamgaging sites. This study found that effects of groundwater withdrawals on streamflow were observed as statistically significant reductions in baseflow in areas with substantial groundwater-altitude declines. Hydrograph-separation and trend analyses may be applicable to assess the impacts of groundwater withdrawals in altered environments and streamflow may be used as a proxy for changes in groundwater availability. Characterizing and defining hydrologic relations between groundwater and surface water will help scientists and water-resource managers refine a regional groundwater-flow model that includes the Mississippi Delta that will be used to aid water-resource managers in future decisions concerning the alluvial aquifer.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗
USGS2007

New isotopic evidence for the origin of groundwater from the Nubian Sandstone Aquifer in the Negev, Israel

The geochemistry and isotopic composition (H, O, S, O sulfate , C, Sr) of groundwater from the Nubian Sandstone (Kurnub Group) aquifer in the Negev, Israel, were investigated in an attempt to reconstruct the origin of the water and solutes, evaluate modes of water–rock interactions, and determine mean residence times of the water. The results indicate multiple recharge events into the Nubian sandstone aquifer characterized by distinctive isotope signatures and deuterium excess values. In the northeastern Negev, groundwater was identified with deuterium excess values of ∼16‰, which suggests local recharge via unconfined areas of the aquifer in the Negev anticline systems. The δ 18 O H2O and δ 2 H values (−6.5‰ and −35.4‰) of this groundwater are higher than those of groundwater in the Sinai Peninsula and southern Arava valley (−7.5‰ and −48.3‰) that likewise have lower deuterium excess values of ∼10‰. Based on the geochemical differences between groundwater in the unconfined and confined zones of the aquifer, a conceptual geochemical model for the evolution of the groundwater in the Nubian sandstone aquifer has been reconstructed. The isotopic composition of shallow groundwater from the unconfined zone indicates that during recharge oxidation of pyrite to SO 4 ( δ 34 S SO4 ∼−13‰; δ 18 O SO4 ∼+7.7‰) and dissolution of CaCO 3 ( 87 Sr/ 86 Sr ∼0.70787; δ 13 C DIC = −3.7‰) occur. In the confined zone of the aquifer, bacterial SO 4 reduction removes a significant part of dissolved SO 4 2 - "> SO42- , thereby modifying its isotopic composition ( δ 34 S SO4 ∼−2‰; δ 18 O SO4 ∼+8.5‰) and liberating dissolved inorganic C that contains little or no radiocarbon ( 14 C-free) with low δ 13 C DIC values (<−12‰). In addition to local recharge, the Sr and S isotopic data revealed contribution of external groundwater sources to the Nubian Sandstone aquifer, resulting in further modifications of the groundwater chemical and isotopic signatures. In the northeastern Negev, it is shown that SO 4 -rich groundwater from the underlying Jurassic aquifer contributes significantly to the salt budget of the Nubian Sandstone aquifer. The unique chemical and isotopic composition of the Jurassic groundwater ( δ 34 S SO4 ∼ +14‰; δ 18 O SO4 ∼ 14‰; 87 Sr/ 86 Sr ∼0.70764) is interpreted as reflecting dissolution of Late Triassic marine gypsum deposits. In the southern Arava Valley the authors postulate that SO 4 -rich groundwater with distinctively high Br/Cl (3 × 10 −3 ) low 87 Sr/ 86 Sr (0.70734), and high δ 34 S SO4 values (+15‰) is derived from mixing with underlying brines from the Paleozoic units. The radiocarbon measurements reveal low 14 C activities (0.2–5.8 pmc) in both the northeastern Negev and southern Arava Valley. Taking into account dissolution of carbonate rocks and bacterial SO 4 reduction in the unconfined area, estimated mean residence times of groundwater in the confined zone in the northeastern Negev are on the order of 21–38 ka, which suggests recharge predominantly during the last glacial period. The 14 C signal in groundwater from the southern Arava Valley is equally low but due to evidence for mixing with external water sources the residence time estimates are questionable.

Applied Geochemistry↗
USGS2014

Assessing the solubility controls on vanadium in groundwater, northeastern San Joaquin Valley, CA

The solubility controls on vanadium (V) in groundwater were studied due to concerns over possible harmful health effects of ingesting V in drinking water. Vanadium concentrations in the northeastern San Joaquin Valley ranged from <3 μg/L to 70 μg/L with a median of 21 μg/L. Concentrations of V were highest in samples collected from oxic groundwater (49% > 25 μg/L) and lowest in samples collected from anoxic groundwater (70% < 0.8 μg/L). In oxic groundwater, speciation modeling (SM) using PHREEQC predicted that V exists primarily as the oxyanion H 2 VO 4 − . Adsorption/desorption reactions with mineral surfaces and associated oxide coatings were indicated as the primary solubility control of V 5+ oxyanions in groundwater. Environmental data showed that V concentrations in oxic groundwater generally increased with increasing groundwater pH. However, data from adsorption isotherm experiments indicated that small variations in pH (7.4–8.2) were not likely as an important a factor as the inherent adsorption capacity of oxide assemblages coating the surface of mineral grains. In suboxic groundwater, accurate SM modeling was difficult since Eh measurements of source water were not measured in this study. Vanadium concentrations in suboxic groundwater decreased with increasing pH indicating that V may exist as an oxycationic species [e.g. V(OH) 3 + ]. Vanadium may complex with dissolved inorganic and organic ligands under suboxic conditions, which could alter the adsorption behavior of V in groundwater. Speciation modeling did not predict the existence of V-inorganic ligand complexes and organic ligands were not collected as part of this study. More work is needed to determine processes governing V solubility under suboxic groundwater conditions. Under anoxic groundwater conditions, SM predicts that aqueous V exists as the uncharged V(OH) 3 molecule. However, exceedingly low V concentrations show that V is sparingly soluble in anoxic conditions. Results indicated that V may be precipitating as V 3+ - or mixed V 3+ /Fe 3+ -oxides in anoxic groundwater, which is consistent with results of a previous study. The fact that V appears insoluble in anoxic (Fe reducing) redox conditions indicates that the behavior of V is different than arsenic (As) in aquifer systems where the reductive dissolution of Fe-oxides with As adsorbed to the surface is a well-documented mechanism for increasing As concentrations in groundwater. This hypothesis is supported by the relation of V to As concentrations in oxic versus anoxic redox conditions. Sequential extraction procedures (SEP) applied to aquifer material showed that the greatest amount of V was recovered by the nitric acid (HNO 3 ) extract (37–71%), followed by the oxalate-ascorbic acid extract (19–60%) and the oxalate extract (3–14%). These results indicate that V was not associated with the solid phase as an easily exchangeable fraction. Although the total amount of V recovered was greatest for the HNO 3 extract that targets V adsorbed to sorption sites of crystalline Al, Fe and Mn oxides, the greatest V saturation of sorption sites appeared to occur on the amorphous and poorly crystalline oxide solid phases targeted by the oxalate and oxalate-ascorbic acid extracts respectively. Adsorption isotherm experiments showed no correlation between V sorption and any of the fractions identified by the SEP. This lack of correlation indicates the application of an SEP alone is not adequate to estimate the sorption characteristics of V in an aquifer system.

California↗
USGS2014

Dynamics of submarine groundwater discharge and associated fluxes of dissolved nutrients, carbon, and trace gases to the coastal zone (Okatee River estuary, South Carolina)

Multiple techniques, including thermal infrared aerial remote sensing, geophysical and geological data, geochemical characterization and radium isotopes, were used to evaluate the role of groundwater as a source of dissolved nutrients, carbon, and trace gases to the Okatee River estuary, South Carolina. Thermal infrared aerial remote sensing surveys illustrated the presence of multiple submarine groundwater discharge sites in Okatee headwaters. Significant relationships were observed between groundwater geochemical constituents and 226 Ra activity in groundwater with higher 226 Ra activity correlated to higher concentrations of organics, dissolved inorganic carbon, nutrients, and trace gases to the Okatee system. A system-level radium mass balance confirmed a substantial submarine groundwater discharge contribution of these constituents to the Okatee River. Diffusive benthic flux measurements and potential denitrification rate assays tracked the fate of constituents in creek bank sediments. Diffusive benthic fluxes were substantially lower than calculated radium-based submarine groundwater discharge inputs, showing that advection of groundwater-derived nutrients dominated fluxes in the system. While a considerable potential for denitrification in tidal creek bank sediments was noted, in situ denitrification rates were nitrate-limited, making intertidal sediments an inefficient nitrogen sink in this system. Groundwater geochemical data indicated significant differences in groundwater chemical composition and radium activity ratios between the eastern and western sides of the river; these likely arose from the distinct hydrological regimes observed in each area. Groundwater from the western side of the Okatee headwaters was characterized by higher concentrations of dissolved organic and inorganic carbon, dissolved organic nitrogen, inorganic nutrients and reduced metabolites and trace gases, i.e. methane and nitrous oxide, than groundwater from the eastern side. Differences in microbial sulfate reduction, organic matter supply, and/or groundwater residence time likely contributed to this pattern. The contrasting features of the east and west sub-marsh zones highlight the need for multiple techniques for characterization of submarine groundwater discharge sources and the impact of biogeochemical processes on the delivery of nutrients and carbon to coastal areas via submarine groundwater discharge.

South Carolina↗
USGS2024

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia↗