Geology topics

Groundwater & aquifers

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

Follow this topic with RSS

At least 181 records · Page 10Linked to original sources

USGS2023

Hydrogeology and simulated groundwater availability in reaches 3 and 4 of the Washita River aquifer, southern Oklahoma, 1980–2017

The 1973 Oklahoma Groundwater Law (Oklahoma Statutes §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s aquifers to determine the maximum annual yield for each groundwater basin. Because more than 20 years have elapsed since the final order was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an updated hydrologic investigation and evaluated the effects of potential groundwater withdrawals on groundwater flow and availability in reaches 3 and 4 of the Washita River aquifer in southern Oklahoma for a study period spanning 1980–2017. A hydrogeologic framework and conceptual model were developed to guide the construction and calibration of a numerical model of the Washita River aquifer. The numerical model was calibrated to water-table-altitude observations at selected wells, base-flow observations at selected U.S. Geological Survey streamgages, and the conceptual-model recharge. Three types of groundwater-availability scenarios were run using the calibrated numerical model. These scenarios were used to (1) estimate equal-proportionate-share pumping rates, (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. With Washita River main-stem inflows, the 20-, 40-, and 50-year equal-proportionate-share pumping rates under normal recharge conditions were about 3.08 acre-feet per acre per year for reach 3 and about 3.80 acre-feet per acre per year for reach 4. Projected 50-year pumping scenarios were used to simulate the effects of modified well withdrawal rates. Because well withdrawals were less than 1 percent of the calibrated numerical-model water budget, changes to the well pumping rates had little effect on Washita River base flows and groundwater storage in the Washita River aquifer. A hypothetical 10-year drought scenario was used to simulate the potential effects of a prolonged period of reduced recharge on groundwater storage. Groundwater storage at the end of the drought period was 4.6 percent less than the groundwater storage of the calibrated numerical model at the end of the drought period.

Oklahoma↗

USGS2010

Simulated groundwater flow in the Ogallala and Arikaree aquifers, Rosebud Indian Reservation area, South Dakota – Revisions with data through water year 2008 and simulations of potential future scenarios

The Ogallala and Arikaree aquifers are important water resources in the Rosebud Indian Reservation area and are used extensively for irrigation, municipal, and domestic water supplies. Drought or increased withdrawals from the Ogallala and Arikaree aquifers in the Rosebud Indian Reservation area have the potential to affect water levels in these aquifers. This report documents revisions and recalibration of a previously published three-dimensional, numerical groundwater-flow model for this area. Data for a 30-year period (water years 1979 through 2008) were used in steady-state and transient numerical simulations of groundwater flow. In the revised model, revisions include (1) extension of the transient calibration period by 10 years, (2) the use of inverse modeling for steady-state calibration, (3) model calibration to base flow for an additional four surface-water drainage basins, (4) improved estimation of transient aquifer recharge, (5) improved delineation of vegetation types, and (6) reduced cell size near large capacity water-supply wells. In addition, potential future scenarios were simulated to assess the potential effects of drought and increased groundwater withdrawals. The model comprised two layers: the upper layer represented the Ogallala aquifer and the lower layer represented the Arikaree aquifer. The model’s grid had 168 rows and 202 columns, most of which were 1,640 feet (500 meters) wide, with narrower rows and columns near large water-supply wells. Recharge to the Ogallala and Arikaree aquifers occurs from precipitation on the outcrop areas. The average recharge rates used for the steady-state simulation were 2.91 and 1.45 inches per year for the Ogallala aquifer and Arikaree aquifer, respectively, for a total rate of 255.4 cubic feet per second (ft 3 /s). Discharge from the aquifers occurs through evapotranspiration, discharge to streams as base flow and spring flow, and well withdrawals. Discharge rates for the steady-state simulation were 171.3 ft 3 /s for evapotranspiration, 74.4 ft 3 /s for net outflow to streams and springs, and 11.6 ft 3 /s for well withdrawals. Estimated horizontal hydraulic conductivity used for the numerical model ranged from 0.2 to 84.4 feet per day (ft/d) in the Ogallala aquifer and from 0.1 to 4.3 ft/d in the Arikaree aquifer. A uniform vertical hydraulic conductivity value of 4.2x10 -4 ft/d was estimated for the Ogallala aquifer. Vertical hydraulic conductivity was estimated for five zones in the Arikaree aquifer and ranged from 8.8x10 -5 to 3.7 ft/d. Average rates of recharge, maximum evapotranspiration, and well withdrawals were included in the steady-state simulation, whereas the time-varying rates were included in the transient simulation. Inverse modeling techniques were used for steady-state model calibration. These methods were designed to estimate parameter values that are, statistically, the most likely set of values to result in the smallest differences between simulated and observed hydraulic heads and base-flow discharges. For the steady-state simulation, the root mean square error for simulated hydraulic heads for all 383 wells was 27.3 feet. Simulated hydraulic heads were within ±50 feet of observed values for 93 percent of the wells. The potentiometric surfaces of the two aquifers calculated by the steady-state simulation established initial conditions for the transient simulation. For the transient simulation, the difference between the simulated and observed means for hydrographs was within ±40 feet for 98 percent of 44 observation wells. A sensitivity analysis was used to examine the response of the calibrated steady-state model to changes in model parameters including horizontal and vertical hydraulic conductivity, evapotranspiration, recharge, and riverbed conductance. The model was most sensitive to recharge and maximum evapotranspiration and least sensitive to riverbed and spring conductances. To simulate a potential future drought scenario, a synthetic recharge record was created, the mean of which was equal to 64 percent of the average estimated recharge rate for the 30-year calibration period. This synthetic recharge record was used to simulate the last 20 years of the calibration period under drought conditions. Compared with results of the calibrated model, decreases in hydraulic-head values for the drought scenario at the end of the simulation period were as much as 39 feet for the Ogallala aquifer. To simulate the effects of potential increases in pumping, well withdrawal rates were increased by 50 percent from those estimated for the 30-year calibration period for the last 20 years of the calibration period. Compared with results of the calibrated model, decreases in hydraulic-head values for the scenario of increased pumping at the end of the simulation period were as much as 13 feet for the Ogallala aquifer. This numerical model is suitable as a tool to help understand the flow system, to help confirm that previous estimates of aquifer properties were reasonable, and to estimate aquifer properties in areas without data. The model also is useful to help assess the effects of drought and increases in pumping by simulations of these scenarios, the results of which are not precise but may be considered when making water management decisions.

South Dakota↗

USGS2023

Availability of groundwater from the volcanic aquifers of the Hawaiian Islands

The islands of Hawaiʻi were built by basaltic shield volcanoes in the Pacific Ocean. These volcanoes formed aquifers that supply hundreds of millions of gallons of fresh water per day to the islands’ residents and diverse industries. Groundwater discharge from the volcanic aquifers to streams and the coast also supports traditional practices and ecosystems. The aquifers' capacity to yield fresh groundwater is limited, however, by limitations placed on the consequences of groundwater withdrawal, such as water-table depression, rise of saltwater, and reduction of groundwater to streams, springs, and the ocean. The magnitude of these consequences will change if groundwater recharge changes as a result of climate and land-use changes. This fact sheet summarizes groundwater availability and its relation to the consequences of groundwater withdrawals in the volcanic aquifers of Hawaiʻi. This topic was addressed in detail in the Hawaiʻi Volcanic Aquifer Study ( https://doi.org/10.3133/pp1876 ).

Hawaii↗

USGS2017

Groundwater model of the Great Basin carbonate and alluvial aquifer system version 3.0: Incorporating revisions in southwestern Utah and east central Nevada

The groundwater model described in this report is a new version of previously published steady-state numerical groundwater flow models of the Great Basin carbonate and alluvial aquifer system, and was developed in conjunction with U.S. Geological Survey studies in Parowan, Pine, and Wah Wah Valleys, Utah. This version of the model is GBCAAS v. 3.0 and supersedes previous versions. The objectives of the model for Parowan Valley were to simulate revised conceptual estimates of recharge and discharge, to estimate simulated aquifer storage properties and the amount of reduction in storage as a result of historical groundwater withdrawals, and to assess reduction in groundwater withdrawals necessary to mitigate groundwater-level declines in the basin. The objectives of the model for the area near Pine and Wah Wah Valleys were to recalibrate the model using new observations of groundwater levels and evapotranspiration of groundwater; to provide new estimates of simulated recharge, hydraulic conductivity, and interbasin flow; and to simulate the effects of proposed groundwater withdrawals on the regional flow system. Meeting these objectives required the addition of 15 transient calibration stress periods and 14 projection stress periods, aquifer storage properties, historical withdrawals in Parowan Valley, and observations of water-level changes in Parowan Valley. Recharge in Parowan Valley and withdrawal from wells in Parowan Valley and two nearby wells in Cedar City Valley vary for each calibration stress period representing conditions from March 1940 to November 2013. Stresses, including recharge, are the same in each stress period as in the steady-state stress period for all areas outside of Parowan Valley. The model was calibrated to transient conditions only in Parowan Valley. Simulated storage properties outside of Parowan Valley were set the same as the Parowan Valley properties and are not considered calibrated. Model observations in GBCAAS v. 3.0 are groundwater levels at wells and discharge locations; water-level changes; and discharge to springs, evapotranspiration of groundwater, rivers, and lakes. All observations in the model outside of Parowan Valley are considered to represent steady-state conditions. Composite scaled sensitivities indicate the observations of discharge to rivers and springs provide more information about model parameters in the model focus area than do water-level observations. Water levels and water-level changes, however, provide the only information about specific yield and specific storage parameters and provide more information about recharge and withdrawals in Parowan Valley than any other observation group. Comparisons of simulated water levels and measured water levels in Parowan Valley indicated that the model fits the overall trend of declining water levels and provides reasonable estimates of long-term reduction in storage and of storage changes from 2012 to 2013. The conceptual and simulated groundwater budgets for Parowan Valley from November 2012 to November 2013 are similar, with recharge of about 20,000 acre-feet and discharge of about 45,000 acre-feet. In the simulation, historical withdrawals averaging about 28,000 acre-feet per year (acre-ft/yr) cause major changes in the groundwater system in Parowan Valley. These changes include the cessation of almost all natural discharge in the valley and the long-term removal of water from storage. Simulated recharge in Pine Valley of 11,000 acre-ft/yr and in Wah Wah Valley of 3,200 acre-ft/yr is substantially less in GBCAAS v. 3.0 than that simulated by previous model versions. In addition, the valleys have less simulated inflow from and outflow to other hydrographic areas than were simulated by previous model versions. The effects of groundwater development in these valleys, however, are independent of the amount of water recharging in and flowing through the valleys. Groundwater withdrawals in Pine and Wah Wah Valleys will decrease groundwater storage (causing drawdown) until discharge in surrounding areas and mountain springs around the two valleys is reduced by the rate of withdrawal. The model was used to estimate that reducing withdrawals in Parowan Valley from 35,000 to about 22,000 acre-ft/yr would likely stabilize groundwater levels in the valley if recharge varies as it did from about 1950 to 2012. The model was also used to demonstrate that withdrawals of 15,000 acre-ft/yr from Pine Valley and 6,500 acre-ft/yr from Wah Wah Valley could ultimately cause long-term steady-state water-level declines of about 1,900 feet near the withdrawal wells and of more than 5 feet in an area of about 10,500 square miles. The timing of drawdown and capture and the ultimate amount of drawdown are dependent on the proximity to areas of simulated natural groundwater discharge, simulated transmissivity, and simulated storage properties. The model projections are a representation of possible effects.

Nevada, Utah↗

USGS2013

California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project--shallow aquifer assessment

The California State Water Resources Control Board’s (SWRCB) GAMA Program is a comprehensive assessment of statewide groundwater quality in California. From 2004 to 2012, the GAMA Program’s Priority Basin Project focused on assessing groundwater resources used for public drinking-water supplies. More than 2,000 public-supply wells were sampled by U.S. Geological Survey (USGS) for this effort. Starting in 2012, the GAMA Priority Basin Project began an assessment of water resources in shallow aquifers in California. These shallow aquifers provide water for domestic and small community-supply wells, which are often drilled to shallower depths in the groundwater system than public-supply wells. Shallow aquifers are of interest because shallow groundwater may respond more quickly and be more susceptible to contamination from human activities at the land surface, than the deeper aquifers. The SWRCB’s GAMA Program was developed in response to the Groundwater Quality Monitoring Act of 2001 (Water Code sections 10780-10782.3): a public mandate to assess and monitor the quality of groundwater resources used for drinking-water supplies, and to increase the availability of information about groundwater quality to the public. The U.S. Geological Survey is the technical lead of the Priority Basin Project. Stewardship of California’s groundwater resources is a responsibility shared between well owners, communities, and the State. Participants and collaborators in the GAMA Program include Regional Water Quality Control Boards, Department of Water Resources, Department of Public Health, local and regional groundwater management entities, county and local water agencies, community groups, and private citizens. Well-owner participation in the GAMA Program is entirely voluntary.

California↗

USGS2023

Predevelopment water levels, groundwater recharge, and selected hydrologic properties of aquifer materials, Hinkley and Water Valleys, California

Hydrologic and geophysical data were collected to support updates to an existing groundwater-flow model of Hinkley Valley, California, in the Mojave Desert about 80 miles northeast of Los Angeles, California. These data provide information on predevelopment (pre-1930) water levels, groundwater recharge, and selected hydrologic properties of aquifer materials. A predevelopment groundwater-level map, drawn using water-level measurements from 48 wells collected as early as 1918, showed groundwater movement from recharge areas along the Mojave River to evaporative discharge areas near the margin of Harper (dry) Lake in Water Valley. During predevelopment conditions, depth to water ranged from near land surface along the Mojave River to above land surface near Harper (dry) Lake, consistent with flowing wells in Water Valley at that time. Depths to water in much of Hinkley Valley downgradient from the Lockhart fault were less than 20 feet below land surface. By 2017, water-level declines as a result of agricultural pumping, were as much as 60 feet near the Hinkley compressor station. Areal recharge from infiltration of precipitation on the valley floor is negligible. Average annual recharge as infiltration of runoff from upland drainages to Hinkley and Water Valleys averages 64.7 acre-feet per year. In most years recharge does not occur; in years when it occurs, recharge to Hinkley Valley is typically about 296 acre-feet. In contrast, average recharge as infiltration of streamflow from the Mojave River from 1931 to 2015 was between 13,400 and 17,100 acre-feet per year; in some years recharge from the Mojave River exceeded 100,000 acre-ft. Estimates of predevelopment groundwater movement through Hinkley Gap and groundwater discharge to Harper (dry) Lake ranged from 570 to 1,900 and 820 to 2,460 acre-feet per year, respectively; at the time of this study in 2017, groundwater movement through Hinkley Gap was estimated to be about 83 acre-feet per year. Hydraulic-conductivity values estimated from slug-test data for 95 monitoring wells ranged from less than 0.1 to 680 feet per day (ft/d); values generally decreased with depth. Median hydraulic-conductivity values calculated from nuclear magnetic resonance (NMR) data for Mojave River alluvium and near-shore lake deposits were 73 and 11 ft/d, respectively; median hydraulic-conductivity values for locally derived alluvium and weathered bedrock were 6 and 2 ft/d, respectively. Hydraulic-conductivity values, estimated from NMR data for formerly saturated deposits overlying the 2017 water table, were as high as 300 ft/d near the Hinkley compressor station. Downgradient from the Hinkley compressor station, formerly saturated deposits had hydraulic-conductivity values of about 150 ft/d, which were higher than values in saturated material. Coarse-textured, permeable material in formerly saturated deposits above the 2017 water table may have allowed groundwater, released from the Hinkley compressor station that may have contained Cr(VI), to move rapidly downgradient. The Lockhart fault is an impediment to groundwater flow within Hinkley Valley. Groundwater-flow directions from horizontal point-velocity probe data were deflected to the west on the upgradient side of the fault compared to the nominal direction of groundwater flow estimated from water-level data. Younger groundwater was present on the upgradient and downgradient sides of the fault, and older groundwater with unadjusted carbon-14 ages as old as 5,650 years before present was in water from wells within splays of the Lockhart fault, consistent with limited groundwater movement across the fault. As a result, groundwater and Cr(VI) released from the Hinkley compressor station moved to the northwest along the downgradient side of the fault. Coupled well-bore flow and depth-dependent water-quality data show water from wells C-01 and IW-03 within the Q4 2015 (October–December 2015) regulatory Cr(VI) plume was yielded from thin layers within the aquifer that are composed of well-sorted lake-margin (beach) deposits that likely have high lateral and longitudinal connectivity. Collectively, data show highly permeable deposits above the regional water table and thin permeable deposits within saturated portions of the upper aquifer that may have conducted groundwater and Cr(VI) downgradient when releases from the Hinkley compressor station first occurred.

California↗

USGS2016

Groundwater-flow model of the northern High Plains aquifer in Colorado, Kansas, Nebraska, South Dakota, and Wyoming

The High Plains aquifer is a nationally important water resource underlying about 175,000 square miles in parts of eight states: Colorado, Kansas, Oklahoma, Nebraska, New Mexico, South Dakota, Texas, and Wyoming. Droughts across much of the Northern High Plains from 2001 to 2007 have combined with recent (2004) legislative mandates to elevate concerns regarding future availability of groundwater and the need for additional information to support science-based water-resource management. To address these needs, the U.S. Geological Survey began the High Plains Groundwater Availability Study to provide a tool for water-resource managers and other stakeholders to assess the status and availability of groundwater resources. A transient groundwater-flow model was constructed using the U.S. Geological Survey modular three-dimensional finite-difference groundwater-flow model with Newton-Rhapson solver (MODFLOW–NWT). The model uses an orthogonal grid of 565 rows and 795 columns, and each grid cell measures 3,281 feet per side, with one variably thick vertical layer, simulated as unconfined. Groundwater flow was simulated for two distinct periods: (1) the period before substantial groundwater withdrawals, or before about 1940, and (2) the period of increasing groundwater withdrawals from May 1940 through April 2009. A soil-water-balance model was used to estimate recharge from precipitation and groundwater withdrawals for irrigation. The soil-water-balance model uses spatially distributed soil and landscape properties with daily weather data and estimated historical land-cover maps to calculate spatial and temporal variations in potential recharge. Mean annual recharge estimated for 1940–49, early in the history of groundwater development, and 2000–2009, late in the history of groundwater development, was 3.3 and 3.5 inches per year, respectively. Primary model calibration was completed using statistical techniques through parameter estimation using the parameter estimation suite of software with Tikhonov regularization. Calibration targets for the groundwater model included 343,067 groundwater levels measured in wells and 10,820 estimated monthly stream base flows at streamgages. A total of 1,312 parameters were adjusted during calibration to improve the match between calibration targets and simulated equivalents. Comparison of calibration targets to simulated equivalents indicated that, at the regional scale, the model correctly reproduced groundwater levels and stream base flows for 1940–2009. This comparison indicates that the model can be used to examine the likely response of the aquifer system to potential future stresses. Mean calibrated recharge for 1940–49 and 2000–2009 was smaller than that estimated with the soil-water-balance model. This indicated that although the general spatial patterns of recharge estimated with the soil-water-balance model were approximately correct at the regional scale of the Northern High Plains aquifer, the soil-water-balance model had overestimated recharge, and adjustments were needed to decrease recharge to improve the match of the groundwater model to calibration targets. The largest components of the simulated groundwater budgets were recharge from precipitation, recharge from canal seepage, outflows to evapotranspiration, and outflows to stream base flow. Simulated outflows to irrigation wells increased from 7 percent of total outflows in 1940–49 to 38 percent of 1970–79 total outflows and 49 percent of 2000–2009 total outflows.

Colorado, Kansas, Nebraska, South Dakota, Wyoming↗

USGS2021

Natural and anthropogenic geochemical tracers to investigate residence times and groundwater–surface-water interactions in an urban alluvial aquifer

A multi-component geochemical dataset was collected from groundwater and surface-water bodies associated with the urban Fountain Creek alluvial aquifer, Colorado, USA, to facilitate analysis of recharge sources, geochemical interactions, and groundwater-residence times. Results indicate that groundwater can be separated into three distinct geochemical zones based on location within the flow system and proximity to surface water, and these zones can be used to infer sources of recharge and groundwater movement through the aquifer. Rare-earth-element concentrations and detections of wastewater-indicator compounds indicate the presence of effluent from wastewater-treatment plants in both groundwater and surface water. Effluent presence in groundwater indicates that streams in the area lose to groundwater in some seasons and are a source of focused groundwater recharge. Distributions of pharmaceuticals and wastewater-indicator compounds also inform an understanding of groundwater–surface-water interactions. Noble-gas isotopes corroborate rare-earth-element data in indicating geochemical evolution within the aquifer from recharge area to discharge area and qualitatively indicate variable groundwater-residence times and mixing with pre-modern groundwater. Quantitative groundwater-residence times calculated from 3 H/ 3 He, SF 6 , and lumped-parameter modeling generally are less than 20 years, but the presence of mixing with older groundwater of an unknown age is also indicated at selected locations. Future investigations would benefit by including groundwater-age tracers suited to quantification of mixing for both young (years to decades) and old (centuries and millennia) groundwater. This multi-faceted analysis facilitated development of a conceptual model for the investigated groundwater-flow system and illustrates the application of an encompassing suite of analytes in exploring hydrologic and geochemical interactions in complex systems.

Colorado↗

USGS2021

Occurrence and distribution of PFAS in sampled source water of public drinking-water supplies in the surficial aquifer in Delaware, 2018; PFAS and groundwater age-dating results

The U.S. Geological Survey, in cooperation with the Delaware Department of Natural Resources and Environmental Control and the Delaware Geological Survey, conducted a groundwater-quality investigation to (1) describe the occurrence and distribution of PFAS, and (2) document any changes in groundwater quality in the Columbia aquifer public water-supply wells in the Delaware Coastal Plain between 2000 and 2008 and between 2008 and 2018. Thirty public water-supply wells located throughout the Columbia aquifer of the Delaware Coastal Plain were sampled from August through November 2018. Groundwater collected from the wells was analyzed for the occurrence and distribution of 18 per- and polyfluorinated alkyl substances (PFAS) as well as groundwater age. Descriptive statistical analyses were performed to assess PFAS analytical results within the well network and the combined perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) concentrations were compared to the U.S. Environmental Protection Agency’s (EPA) health advisory level (HAL) for informational purposes only and not for evidence of compliance or noncompliance with Federal regulations. The EPA’s HAL is a health-based reference level for public drinking water as supplied to customers and is not applied to source (raw) water. Groundwater-age data were compared for sites sampled in 2000, 2008, and 2018 to document any changes. All samples were analyzed for 18 PFAS using EPA Method 537 (modified). Forty-four percent of the analyzed PFAS were detected in the study well network. Sixteen of the sampled wells have one or more PFAS detections, and as many as eight different PFAS were found in a single sample. Wells with a higher number of PFAS detected (five or more) were in New Castle and Sussex Counties. The PFAS most frequently detected were PFOA, with 47 percent detection; perfluorohexanoic acid (PFHxA), with 33 percent detection; and PFOS and perfluorohexane sulfonate (PFHxS), with 27 percent detection each. PFAS concentrations were below 1,000 parts per trillion (ppt). Two wells exceeded the EPA’s lifetime-drinking water health advisory level of 70 ppt for combined concentrations of PFOA and PFOS. The average age of groundwater entering the screens of the supply wells sampled in 2018 ranged from 8.2 to 45.8 years, with a median groundwater age of 25.7 years. Groundwater age was positively correlated with well depth and negatively correlated with dissolved oxygen. Groundwater age and PFAS concentrations were negatively correlated in the Columbia aquifer. Data from the 23 resampled wells indicate a significant positive difference in the average modeled groundwater-sample-age results. The average groundwater age from samples collected in 2018 was generally 5 years older than the average groundwater age from samples collected in 2008. The same pattern was found during cycle two (2008) of this study, where the 2008 groundwater age was on average 7 years older than the samples collected in 2000. The distribution of groundwater sample ages among the 17 trend wells and during the three study cycles (2000, 2008, and 2018) indicates that sample-age medians were statistically different from zero; well-water sample-age data show a slight increase in groundwater sample age.

Delaware↗

USGS2010

Effects of Changes in Irrigation Practices and Aquifer Development on Groundwater Discharge to the Jobos Bay National Estuarine Research Reserve near Salinas, Puerto Rico

Since 1990, about 75 acres of black mangroves have died in the Jobos Bay National Estuarine Research Reserve near Salinas, Puerto Rico. Although many factors can contribute to the mortality of mangroves, changes in irrigation practices, rainfall, and water use resulted in as much as 25 feet of drawdown in the potentiometric surface of the aquifer in the vicinity of the reserve between 1986 and 2002. To clarify the issue, the U.S. Geological Survey, in cooperation with the Puerto Rico Department of Natural and Environmental Resources, conducted a study to ascertain how aquifer development and changes in irrigation practices have affected groundwater levels and groundwater flow to the Mar Negro area of the reserve. Changes in groundwater flow to the mangrove swamp and bay from 1986 to 2004 were estimated in this study by developing and calibrating a numerical groundwater flow model. The transient simulations indicate that prior to 1994, high irrigation return flows more than offset the effect of reduced groundwater withdrawals. In this case, the simulated discharge to the coast in the modeled area was 19 million gallons per day. From 1994 through 2004, furrow irrigation was completely replaced by micro-drip irrigation, thus eliminating return flows and the simulated average coastal discharge was 7 million gallons per day, a reduction of 63 percent. The simulated average groundwater discharge to the coastal mangrove swamps in the reserve from 1986 to 1993 was 2 million gallons per day, compared to an average simulated discharge of 0.2 million gallons per day from 1994 to 2004. The average annual rainfall for each of these periods was 38 inches. The groundwater discharge to the coastal mangrove swamps in the Jobos Bay National Estuarine Research Reserve was estimated at about 0.5 million gallons per day for 2003-2004 because of higher than average annual rainfall during these 2 years. The groundwater flow model was used to test five alternatives for increasing groundwater discharge to the coastal mangrove swamps to approximately 1.4 million gallons per day: (1) artificially recharging the aquifer with injection wells or (2) by increasing irrigation return flow by going back to furrow irrigation; (3) termination of groundwater withdrawals near the mangroves; (4) reduction of groundwater withdrawals at irrigation wells by 50 percent; and (5) a combination of alternatives 2 and 4 increasing irrigation return flows and decreasing irrigation withdrawals. Each alternative assumed average climatic conditions and groundwater withdrawals at 2004 rates. Alternative 1 required 1.5 million gallons per day of injected water. Alternative 2 required flooding 958 acres with a rate of 1.84 million gallons per day if no crops are grown. Alternative 3 required the termination of 2.44 million gallons per day of withdrawals to achieve 1.34 million gallons per day of discharge to the mangroves. Alternative 4 did not achieve the objective with only 0.80 million gallons per day simulated discharge to the mangroves, while requiring a 1.26 million gallon per day reduction in groundwater withdrawals. Alternative 5 required flooding fields with additional 1.13 million gallons of day and the same reduction in groundwater withdrawals, but did achieve the objective of about 1.4 million gallons per day discharge to the mangroves. Alternative 1, incorporating injection wells near the reserve required the least amount of water to raise groundwater levels and maintain discharge of 1.4 million gallons per day through the mangroves.

Scientific Investigations Report↗

USGS2010

Estimating 14C groundwater ages in a methanogenic aquifer

This paper addresses the problem of 14 C age dating of groundwaters in a confined regional aquifer affected by methanogenesis. Increasing CH 4 concentrations along the groundwater flow system and 13 C and 14 C isotopic data for dissolved inorganic carbon, dissolved organic carbon, and CH 4 clearly show the effect of methanogenesis on groundwater chemistry. Inverse reaction path modeling using NETPATH indicates the predominant geochemical reactions controlling the chemical evolution of groundwater in the aquifer are incongruent dissolution of dolomite, ion exchange, methanogenesis, and oxidation of sedimentary organic matter. Modeling of groundwater 14 C ages using NETPATH indicates that a significant part of groundwater in the Alliston aquifer is less than 13,000 years old; however, older groundwater in the range of 15,000–23,000 years is also present in the aquifer. This paper demonstrates that 14C ages calculated using NETPATH, incorporating the effects of methanogenesis on the carbon pools, provide reasonable groundwater ages that were not possible by other isotopic methods.

Water Resources Research↗

USGS2011

Using Cl/Br ratios and other indicators to assess potential impacts on groundwater quality from septic systems: A review and examples from principal aquifers in the United States

A detailed review was made of chemical indicators used to identify impacts from septic tanks on groundwater quality. Potential impacts from septic tank leachate on groundwater quality were assessed using the mass ratio of chloride–bromide (Cl/Br), concentrations of selected chemical constituents, and ancillary information (land use, census data, well depth, soil characteristics) for wells in principal aquifers of the United States. Chemical data were evaluated from 1848 domestic wells in 19 aquifers, 121 public-supply wells in 6 aquifers, and associated monitoring wells in four aquifers and their overlying hydrogeologic units. Based on previously reported Cl/Br ratios, statistical comparisons between targeted wells (where Cl/Br ratios range from 400 to 1100 and Cl concentrations range from 20 to 100 mg/L) and non-targeted wells indicated that shallow targeted monitoring and domestic wells (<20 m depth below land surface) had a significantly ( p < 0.05) higher median percentage of houses with septic tanks (1990 census data) than non-targeted wells. Higher ( p = 0.08) median nitrate–N concentration (3.1 mg/L) in oxic (dissolved oxygen concentrations >0.5 mg/L) shallow groundwater from target domestic wells, relative to non-target wells (1.5 mg/L), corresponded to significantly higher potassium, boron, chloride, dissolved organic carbon, and sulfate concentrations, which may also indicate the influence of septic-tank effluent. Impacts on groundwater quality from septic systems were most evident for the Eastern Glacial Deposits aquifer and the Northern High Plains aquifer that were associated with the number of housing units using septic tanks, high permeability of overlying sediments, mostly oxic conditions, and shallow wells. Overall, little or no influence from septic systems were found for water samples from the deeper public-supply wells. The Cl/Br ratio is a useful first-level screening tool for assessing possible septic tank influence in water from shallow wells (<20 m) with the range of 400–1100. The use of this ratio would be enhanced with information on other chloride sources, temporal variability of chloride and bromide concentrations in shallow groundwater, knowledge of septic-system age and maintenance, and the use of multiple tracers (combination of additional chemical and microbiological indicators).

Journal of Hydrology↗

USGS2013

Analysis and simulation of regional subsidence accompanying groundwater abstraction and compaction of susceptible aquifer systems in the USA

Regional aquifer-system compaction and land subsidence accompanying groundwater abstraction in susceptible aquifer systems in the USA is a challenge for managing groundwater resources and mitigating associated hazards. Developments in the assessment of regional subsidence provide more information to constrain analyses and simulation of aquifer-system compaction. Current popular approaches to simulating vertical aquifer-system deformation (compaction), such as those embodied in the aquitard drainage model and the MODFLOW subsidence packages, have proven useful from the perspective of regional groundwater resources assessment. However, these approaches inadequately address related local-scale hazards-ground ruptures and damages to engineered structures on the land surface arising from tensional stresses and strains accompanying groundwater abstraction. This paper presents a brief overview of the general approaches taken by the U.S. Geological Survey toward understanding aquifer-system compaction and subsidence with regard to a) identifying the affected aquifer systems; b) making regional assessments; c) analyzing the governing processes; and d) simulating historical and future groundwater flow and subsidence conditions. Limitations and shortcomings of these approaches, as well as future challenges also are discussed.

Boletín de la Sociedad Geológica Mexicana↗

USGS2015

Hydrogeology and simulation of groundwater flow in fractured-rock aquifers of the Piedmont and Blue Ridge Physiographic Provinces, Bedford County, Virginia

An annual groundwater budget was computed as part of a hydrogeologic characterization and monitoring effort of fractured-rock aquifers in Bedford County, Virginia, a growing 764-square-mile (mi 2 ) rural area between the cities of Roanoke and Lynchburg, Virginia. Data collection in Bedford County began in the 1930s when continuous stream gages were installed on Goose Creek and Big Otter River, the two major tributaries of the Roanoke River within the county. Between 2006 and 2014, an additional 2 stream gages, 3 groundwater monitoring wells, and 12 partial-record stream gages were operated. Hydrograph separation methods were used to compute base-flow recharge rates from the continuous data collected from the continuous stream gages. Mean annual base-flow recharge ranged from 8.3 inches per year (in/yr) for the period 1931–2012 at Goose Creek near Huddleston (drainage area 188 mi 2 ) to 9.3 in/yr for the period 1938–2012 at Big Otter River near Evington (drainage area 315 mi 2 ). Mean annual base-flow recharge was estimated to be 6.5 in/yr for the period 2007–2012 at Goose Creek at Route 747 near Bunker Hill (drainage area 125 mi 2 ) and 8.9 in/yr for the period 2007–2012 at Big Otter River at Route 221 near Bedford (drainage area 114 mi 2 ). Base-flow recharge computed from the partial-record data ranged from 5.0 in/yr in the headwaters of Goose Creek to 10.5 in/yr in the headwaters of Big Otter River. A steady-state groundwater-flow simulation for Bedford County was developed to test the conceptual understanding of flow in the fractured-rock aquifers and to compute a groundwater budget for the four major drainages: James River, Smith Mountain and Leesville Lakes, Goose Creek, and Big Otter River. Model results indicate that groundwater levels mimic topography and that minimal differences in aquifer properties exist between the Proterozoic basement crystalline rocks and Late Proterozoic-Cambrian cover crystalline rocks. The Big Otter River receives 40.8 percent of the total daily groundwater outflow from fractured-rock aquifers in Bedford County; Goose Creek receives 25.8 percent, the James River receives 18.2 percent, and Smith Mountain and Leesville Lakes receive 15.2 percent. The remaining percentage of outflow is attributed to pumping from the aquifer (consumptive use).

Virginia↗

USGS2024

Assessment of the interconnection between Tampa Bay and the Floridan aquifer system: Historical groundwater data compilation and analysis, 1976–2022

The U.S. Geological Survey used existing data collected after the last major navigational channel modification in the mid-1980s to investigate groundwater levels and chloride concentrations in wells in the Floridan aquifer system and other aquifers beneath and near Tampa Bay. Tampa Bay is located on the west-central coast of Florida and provides access for commercial shipping. In 2021, the U.S. Army Corps of Engineers began to investigate alternatives to improve the efficiency of the deep-draft navigation channels within Tampa Bay, reduce costs, improve safety, and analyze the implications of modifying navigational channels. The Floridan aquifer system underlies Tampa Bay and is the primary source of public water supply in the region. Modifications to the channels have the potential to expose the Floridan aquifer system to the saltwater in Tampa Bay, with the potential to increase salinity in wells in the region. Other factors affecting the groundwater levels and the location of the freshwater/saltwater interface were also examined, including changes in sea level, groundwater extraction, and variations in climate. Groundwater levels and well-construction reports were used to identify whether different aquifer units are well-connected. Twelve wells had available data before and after the last major channel modifications, which took place in the 1980s, with six datasets of chloride concentration available in areas along the northern and eastern coastline of Tampa Bay, which is nearest to historical dredging activities. Of these six, Kendall’s τ and p -values indicated increasing trends in chloride concentration for three datasets (TR 11-2, TR 10-2, 51), no trend in chloride concentration for two datasets (TR 9-3, 50), and a decreasing trend in chloride concentration for one dataset (TR 9-1). The upward trends in chloride concentration observed for TR 10-2 and 51 are likely the result of changes in local groundwater withdrawals. Well TR 11-2 had a gradual increasing trend in chloride concentration, fresh groundwater throughout the period of record, and a 3- to 4-foot increase in hydraulic head during the period of record, possibly caused by the construction and control of the Tampa Bypass Canal, resulting in changes to the regional potentiometric surface.

Florida↗

USGS2010

Sources of groundwater based on Helium analyses in and near the freshwater/saline-water transition zone of the San Antonio segment of the Edwards Aquifer, South-Central Texas, 2002-03

This report evaluates dissolved noble gas data, specifically helium-3 and helium-4, collected by the U.S. Geological Survey, in cooperation with the San Antonio Water System, during 2002-03. Helium analyses are used to provide insight into the sources of groundwater in the freshwater/saline-water transition zone of the San Antonio segment of the Edwards aquifer. Sixty-nine dissolved gas samples were collected from 19 monitoring wells (categorized as fresh, transitional, or saline on the basis of dissolved solids concentration in samples from the wells or from fluid-profile logging of the boreholes) arranged in five transects, with one exception, across the freshwater/saline-water interface (the 1,000-milligrams-per-liter dissolved solids concentration threshold) of the Edwards aquifer. The concentration of helium-4 (the dominant isotope in atmospheric and terrigenic helium) in samples ranged from 63 microcubic centimeters per kilogram at standard temperature (20 degrees Celsius) and pressure (1 atmosphere) in a well in the East Uvalde transect to 160,587 microcubic centimeters per kilogram at standard temperature and pressure in a well in the Kyle transect. Helium-4 concentrations in the 10 saline wells generally increase from the western transects to the eastern transects. Increasing helium-4 concentrations from southwest to northeast in the transition zone, indicating increasing residence time of groundwater from southwest to northeast, is consistent with the longstanding conceptualization of the Edwards aquifer in which water recharges in the southwest, flows generally northeasterly (including in the transition zone, although more slowly than in the fresh-water zone), and discharges at major springs in the northeast. Excess helium-4 was greater than 1,000 percent for 60 of the 69 samples, indicating that terrigenic helium is largely present and that most of the excess helium-4 comes from sources other than the atmosphere. The helium data of this report cannot be used to identify sources of groundwater in and near the transition zone of the Edwards aquifer in terms of specific geologic (stratigraphic) units or hydrogeologic units (aquifers or confining units). However, the data indicate that the source or sources of the helium, and thus the water in which the helium is dissolved, in the transition zone are mostly terrigenic in origin rather than atmospheric. Whether most helium in and near the transition zone of the Edwards aquifer originated either in rocks outside the transition zone and at depth or in the adjacent Trinity aquifer is uncertain; but most of the helium in the transition zone had to enter the transition zone from the Trinity aquifer because the Trinity aquifer is the hydrogeologic unit immediately beneath and laterally adjacent to the transition zone of the Edwards aquifer. Thus the helium data support a hypothesis of sufficient hydraulic connection between the Trinity and Edwards aquifers to allow movement of water from the Trinity aquifer to the transition zone of the Edwards aquifer.

Texas↗

USGS2026

Status and understanding of groundwater quality in the San Joaquin Valley Kern County subbasin domestic-supply aquifer study unit, 2022—California GAMA Priority Basin Project

The quality of water accessed by domestic wells (here referred to as domestic groundwater resources) in the San Joaquin Valley Kern County subbasin (basin number 5-022.14) was assessed as part of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project (GAMA-PBP), in cooperation with the California State Water Resources Control Board. Kern County is at the southern end of the San Joaquin Valley in California, and about 30,000 residents are estimated to use privately owned domestic wells for drinking water. Domestic wells typically draw from shallower parts of the aquifer system than public-supply wells and can be more vulnerable to effects from surface activities. Kern County is host to a highly productive agricultural industry, with Bakersfield as the main urban center. The Kern River runs through Bakersfield from the southern Sierra Nevada and intersects the Kern Water Bank, one of the largest groundwater banking operations in California, at the Kern River Intertie. The section of the Kern River running through the Kern Water Bank is dry most years. Kern County also encompasses some of the most productive oil and gas basins in California, with extensive underground and surface disposal of oil-field wastewater. This study was based on data collected from 33 sites sampled by the U.S. Geological Survey for the GAMA-PBP in 2022. To provide context for the water quality assessment, measured concentrations were compared to regulatory and non-regulatory health-based and aesthetic benchmarks. A grid-based method was used to estimate the proportions of the groundwater resources used for domestic-supply wells that have water-quality constituents below (low relative concentration), approaching (moderate relative concentration), or above (high relative concentration) benchmark concentrations. At least one measured constituent with a regulatory benchmark was categorized as having a high relative concentration in 72 percent of the aquifer area used for domestic groundwater resources. Inorganic constituents were detected at high concentrations in 45 percent of the domestic groundwater resources, and the constituents detected above regulatory benchmarks were arsenic, nitrate, and uranium. At least one organic constituent was detected at high concentrations in 41 percent of the domestic groundwater resources, and the constituents exceeding regulatory benchmarks were the fumigants 1,2,3-trichloropropane (1,2,3-TCP), 1,2-dibromo-3-chloropropane (dibromochloropropane [DBCP]), 1,2-dibromoethane (EDB), and the per-and polyfluoroalkyl substance (PFAS) perfluorooctanesulfonate. The disinfection by-product chloroform, the fumigant 1,2-dichloropropane, the herbicides atrazine and hexazinone, and the herbicide degradates 2-chloro-6-ethylamino-4-amino-s-triazine, 2-chloro-4,6-diamino-s-triazine, 4-hydroxychlorothalonil, and metolachlor sulfonic acid were detected in more than 10 percent of domestic groundwater resources, but concentrations did not exceed regulatory benchmarks. Land use, groundwater age (fraction of modern water and mean age), and geochemical environment (oxic or anoxic conditions, pH, alkalinity) were associated with the distribution of high relative concentrations of inorganic and organic constituents. Young, oxygenated water is recharged along the Kern River and adjacent recharge ponds, or as irrigation water in the agricultural areas. High concentrations of nitrate and volatile organic compounds occurred in the oxic water in urban and agricultural areas. The fumigants 1,2,3-TCP, DBCP, and EDB were reported throughout the agricultural areas, whereas chloroform, tetrachloroethene, and PFAS were associated with urban land use. High uranium concentrations were associated with young, modern groundwater in agricultural areas with low pH and high bicarbonate. Total dissolved solids increased with distance from the Kern River, as the contributions of fresh, oxic water decreased. High concentrations of arsenic were present in older anoxic or alkaline groundwater away from areas of recharge. Overall, groundwater age, redox conditions, and the source of recharge as a result of different land uses contribute to large aquifer-scale portions of domestic groundwater resources that exceed health-based benchmarks for nitrate, uranium, and fumigant concentrations.

California↗

USGS2019

Groundwater quality of a public supply aquifer in proximity to oil development, Fruitvale Oil Field, Bakersfield, California

Due to concerns over the effects of oil production activities on groundwater quality in California, chemical, isotopic, dissolved gas and age-dating tracers were analyzed in samples collected from public-supply wells and produced-water sites in the Fruitvale oil field (FVOF). A combination of newly collected and historical data was used to determine whether oil formation fluids have mixed with groundwater used for public supply and what the potential pathways for the migration of oil formation fluids into groundwater may be. Stable isotopes of water (δ 2 H and δ 18 O) and age dating ( 3 H, 3 He trit , SF 6 and 14 C) tracers in groundwater samples were consistent with the Kern River being the main source of recharge to aquifers. The distribution of major ion concentrations and pH with distance from the Kern River indicate that natural processes were the primary controls on groundwater salinity. Two of 14 groundwater samples had δ 13 C-DIC values (−2.4 to +1.9 per mil) consistent with mixtures of <1 to about 9 percent oil-field water. Concentrations of TDS in groundwater samples were generally much lower (129–1,200 milligrams per liter (mg/l), median 216 mg/l) than produced water samples (586–24,930 mg/l, median 2,717 mg/l), suggesting that any mixing of oil-field water with groundwater has not significantly affected groundwater salinity. Trace concentrations of thermogenic methane were detected in three groundwater samples that did not have dissolved inorganic or isotopic indicators consistent with mixing of oil-field water, suggesting that stray gases may have migrated from the subsurface via preferential pathways such as leaky well bores into groundwater aquifers. Low concentrations of petroleum hydrocarbons were detected in samples that also contained anthropogenic VOCs and components of post- and pre-1950s recharge, indicating that petroleum hydrocarbons could have come from subsurface and/or surface sources. Overall, the results of this study indicated that groundwater currently used for public supply in the FVOF was of good quality with little, if any, effects from oil production activities. This may be due in part to the relatively rapid flushing of the aquifer system by recharge from the Kern River.

California↗