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 199 records · Page 11Linked to original sources

USGS2026

Water-withdrawal and groundwater-level fluctuations, physical and infiltration characteristics of Ohio River sediment, and isotope-based estimates of water sources relative to a riverbank filtration system in an alluvial aquifer, Louisville, Kentucky, 2000–20

The U.S. Geological Survey and the Louisville Water Company, from 2000 to 2020, cooperatively investigated physical and hydrologic conditions in the Ohio River alluvial aquifer and adjacent riverbed sediments in the northeast part of Louisville, Kentucky, to increase knowledge of processes affecting water sources to collector wells during operation of a riverbank filtration system (RBF). The RBF, in 2020, included an initial collector well (CW1) completed in 1999 and a tunnel and collector well system (TCWS) composed of four collector wells (CW2–CW5) completed in late 2010 that were connected by a tunnel that routed water to a pumping station and a treatment plant. Daily mean water-withdrawal rates typically ranged from about 10.4 to 20.7 million gallons per day (Mgal/d) from 2005 to 2007, then were cyclic from 2008 to 2010 and ranged from as low as about 9 to 13 Mgal/d in winter-spring months up to about 19 to 20 Mgal/d in summer-autumn months. From 2011 through 2017, peak daily mean water-withdrawal rates from CW2–CW5 fluctuated between winter-spring lows of about 20–30 Mgal/d to summer-autumn peak rates of about 41.6 Mgal/d in 2013 to about 48.2 Mgal/d in 2012. From 2017 to 2020, TCWS peak daily mean water-withdrawal rates decreased from about 45.7 Mgal/d in 2017 to about 37.5 Mgal/d in 2020 and peak daily mean rates from CW1 ranged from about 20.7 Mgal/d in 2020 to about 21.6 Mgal/d in 2018. Periods in 2019–20 with tunnel pumping station daily mean water-withdrawal rates smaller than about 30 Mgal/d had generally stable water levels in CW2 and indicated the sustainable yield of the TCWS collector wells under the prevailing hydrologic conditions. Drawdown during July–August of 2008 and 2009 was associated with CW1, the sole operating collector well. Groundwater-flow directions outside the CW1 cone of depression during July–August of 2008 and 2009 were from areas near the alluvial valley margin toward the Ohio River and in 2009 toward Harrods Creek. Collector wells CW2 through CW5 derived water during 2011–19 surveys from surface-water infiltration, as indicated by drawdown that extended beneath the Ohio River and Harrods Creek. Increases in drawdown toward the alluvial valley margins in July–August of 2013 and 2015, compared with July 2011, indicated less available recharge from those flow directions and decreased yield of groundwater from the aquifer after progressively longer TCWS water withdrawals. Drawdown at CW2–CW5 and at observation wells between them lessened in October 2019 and indicated that onshore groundwater levels can recover when collector well withdrawals decrease. Overlapping cones of depression in 2011–15 potentiometric surfaces centered around CW2, CW3, CW4, and CW5, and drawdown along divides between those wells indicated well interference between collector wells. Well-interference effects increased during 2013 and 2015 survey periods, despite decreased water-withdrawal rates compared with July 2011, and expanded to include CW1 when its use resumed. Groundwater-level fluctuations in observation wells near the Ohio River were consistent with gradual depletion of groundwater yield during tunnel pumping station operation at river elevations of about 420 feet above the National Geodetic Vertical Datum of 1929 (NGVD 29) and enhanced recharge of the aquifer when river elevations were above flood stage of about 430 feet above NGVD 29. Groundwater-level declines in observation wells near the Ohio River were most apparent in summer and some autumn months between 2012 and 2017 during periods when tunnel pumping station mean daily water-withdrawal rates were about 35 Mgal/d or more. Periods of relatively stable groundwater levels during 2019–20 in observation wells between CW2–CW5, when tunnel pumping station water-withdrawal rates were about 30–33 Mgal/d or less, also indicated that the prevailing sustainable yield was about 30 Mgal/d or less. Wells associated with CW3 and CW4 had the largest increases in slopes of best-fit lines relating water-level difference to water-withdrawal rate between 2012 and 2020, an indication of decreased efficiencies of those wells to supply water or of changes in groundwater yield around the well. Increased slopes of best-fit lines relating water-level difference to water-withdrawal rate from 2019 to 2020, when groundwater levels were otherwise relatively stable, also indicated that the specific capacity of one or more TCWS collector wells may have declined during that period. Riverbed surface elevation, roughness, and hardness changes between hydroacoustic surveys indicated active accretion and erosion of bed sediment offshore from the riverbank filtration system in 2002–4. Riverbed hydraulic conductivity computed from 2013 to 2015 in tests at 11 sites ranged from 0.2 to 230 feet per day. Ten temporary piezometer sites with adjacent operating collector wells had drawdowns that ranged from about 4 to 11.7 feet and hammer blows greater than or equal to 49 counts, which indicated a relatively great density of alluvial sediments and the potential for diminished riverbed infiltration rates. Contributions of Ohio River infiltration to collector wells estimated from deuterium and oxygen-18 isotopic ratios in representative Ohio River and groundwater samples were considerably larger in warmer months than cooler months and ranged from 80–90 percent in CW1–CW5 in May 2020 to 100 percent in May 2019 at CW2, CW4, and CW5. Groundwater contributions to water withdrawals by collector wells increased in cooler autumn months and ranged from 40 percent at CW3 and CW5 in November 2020 to 70 percent at CW2, CW4, and CW5 in November–December 2019. November 2019 and November 2020 isotopic compositions of water from the collector wells had enriched deuterium and oxygen-18 isotopic ratios consistent with mixtures of evaporation-modified water from the Ohio River and meteoric-origin groundwater.

Indiana, Kentucky↗

USGS2018

Comparison of groundwater age models for assessing nitrate loading, transport pathways, and management options in a complex aquifer system

In an aquifer system with complex hydrogeology, mixing of groundwater with different ages could occur associated with various flow pathways. In this study, we applied different groundwater age estimation techniques (lumped parameter model, and numerical model) to characterize groundwater age distributions and the major pathways of nitrate contamination in the Gosan agricultural field, Jeju Island. According to the lumped parameter model, groundwater age in the study area could be explained by the binary mixing of the young groundwater (4-33 years) and the old water component (>60 years). The complex hydrogeologic regimes and local heterogeneity observed in the study area (multi-layered aquifer, well leakage hydraulics) were particularly well reflected in the numerical model. The numerical model predicted that the regional aquifer of Gosan responded to the fertilizer applications more rapidly (mean age: 9.7-22.3 years) than as estimated by other models. Our study results demonstrated that application and comparison of multiple age estimation methods can be useful to understand better the flow regimes and the mixing characteristics of groundwater with different ages (pathways), hence, to reduce the risk of improper groundwater management plan arising from the aquifer heterogeneity.

Jeju Island↗

USGS2009

Groundwater availability of the Central Valley Aquifer, California

California's Central Valley covers about 20,000 square miles and is one of the most productive agricultural regions in the world. More than 250 different crops are grown in the Central Valley with an estimated value of $17 billion per year. This irrigated agriculture relies heavily on surface-water diversions and groundwater pumpage. Approximately one-sixth of the Nation's irrigated land is in the Central Valley, and about one-fifth of the Nation's groundwater demand is supplied from its aquifers. The Central Valley also is rapidly becoming an important area for California's expanding urban population. Since 1980, the population of the Central Valley has nearly doubled from 2 million to 3.8 million people. The Census Bureau projects that the Central Valley's population will increase to 6 million people by 2020. This surge in population has increased the competition for water resources within the Central Valley and statewide, which likely will be exacerbated by anticipated reductions in deliveries of Colorado River water to southern California. In response to this competition for water, a number of water-related issues have gained prominence: conservation of agricultural land, conjunctive use, artificial recharge, hydrologic implications of land-use change, and effects of climate variability. To provide information to stakeholders addressing these issues, the USGS Groundwater Resources Program made a detailed assessment of groundwater availability of the Central Valley aquifer system, that includes: (1) the present status of groundwater resources; (2) how these resources have changed over time; and (3) tools to assess system responses to stresses from future human uses and climate variability and change. This effort builds on previous investigations, such as the USGS Central Valley Regional Aquifer System and Analysis (CV-RASA) project and several other groundwater studies in the Valley completed by Federal, State and local agencies at differing scales. The principal product of this new assessment is a tool referred to as the Central Valley Hydrologic Model (CVHM) that accounts for integrated, variable water supply and demand, and simulates surface-water and groundwater-flow across the entire Central Valley system. The development of the CVHM comprised four major elements: (1) a comprehensive Geographic Information System (GIS) to compile, analyze and visualize data; (2) a texture model to characterize the aquifer system;(3) estimates of water-budget components by numerically modeling the hydrologic system with the Farm Process (FMP); and (4) simulations to assess and quantify hydrologic conditions.

Professional Paper↗

USGS2021

Machine-learning predictions of redox conditions in groundwater in the Mississippi River Valley alluvial and Claiborne aquifers, south-central United States

Machine-learning models developed by the U.S. Geological Survey were used to predict iron concentrations and the probability of dissolved oxygen (DO) concentrations exceeding a threshold of 1 milligram per liter (mg/L) in groundwater in aquifers of the Mississippi embayment physiographic region. DO and iron concentrations are driven by and reflect the oxidation-reduction (redox) conditions in groundwater. Predictions from boosted regression trees, a type of machine-learning model, of iron concentration and DO threshold probability were used to categorize redox zones in the Mississippi River Valley alluvial aquifer (MRVA), middle Claiborne aquifer (MCAQ), and lower Claiborne aquifer (LCAQ). Model predictions indicated that DO concentrations greater than 1 mg/L are uncommon across the MRVA. DO events (where the predicted probability was greater than 0.5) tended to occur on the margins of the MRVA and in upland areas where MCAQ and LCAQ units crop out at the surface or are at shallow depth. Predicted iron concentrations were higher in the MRVA than in the MCAQ and LCAQ. Uncertainty in predicted iron concentrations tended to be high in areas where measured concentrations were also high, resulting in small areas (encompassing less than 1.5 percent of the areal extent of the MRVA) of predicted iron concentrations that exceeded 100,000 micrograms per liter. Despite the large magnitude of overpredicted iron concentrations, the general proportion and spatial distribution of predicted iron concentrations reflected observed concentrations in groundwater wells. Where the probability of exceeding a DO concentration of 1 mg/L was 0.8 or more and the iron concentration was less than 1,000 micrograms per liter, aquifers were categorized as oxic. Oxic conditions were mostly in the uplands where MCAQ and LCAQ units crop out at the margins of the modeled area. The MRVA was mostly anoxic, which was controlled by DO threshold probabilities less than 0.1. The predictions and redox zones support conceptual models of redox conditions in the Mississippi embayment. The MRVA is predominantly anoxic with high iron concentrations. In the Claiborne aquifers (including the MCAQ and LCAQ), groundwater flows along regional flow paths toward the axis of the Mississippi embayment (the approximate location of the Mississippi River), the residence time in the aquifer increases, DO is consumed, and iron concentrations generally increase. Elevated concentrations of trace elements, such as manganese and arsenic, are often associated with reducing conditions in anoxic and mixed anoxic zones, but other factors such as sediment mineralogy affect the occurrence and distribution of these constituents.

Alabama, Arkansas, Kentucky, Louisiana, Mississipp↗

USGS2018

Groundwater-level analysis of selected wells in the Hoosic River Valley near Hoosick Falls, New York, for aquifer framework and properties

The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, analyzed groundwater levels, drilling record logs, and field water-quality data from selected wells, and the surficial geology in the Hoosic River valley south of the village of Hoosick Falls, New York, to provide information about the framework and properties of a confined aquifer. The aquifer, which consists of ice-contact sand and gravel overlain by lacustrine clay and silt, was evaluated by the New York State Department of Environmental Conservation as part of their investigation of alternate water supplies for the village whose wellfield has been affected by perfluorooctanoic acid. Wells inventoried in the study area were classified as confined, water table, or transitional between the two aquifer conditions. Groundwater levels in three confined-aquifer wells and a transitional-aquifer well responded to pumping of a test production well finished in the confined aquifer. Groundwater levels in a water-table well showed no detectable water-level change in response to test-well pumping. Analysis of drawdown and recovery data from the three confined-aquifer wells and a transitional-aquifer well through the application of the Theis type-curve method provided estimates of aquifer properties. Representation of a constant-head boundary in the analysis where an unnamed pond and fluvial-terrace deposits abut the valley wall resulted in satisfactory matches of the Theis type curves with the observed water-level responses. Aquifer transmissivity estimates ranged from 1,160 to 1,370 feet squared per day. Aquifer storativity estimates ranged from 5.2×10 –5 to 1.1×10 –3 and were consistent with the inferred degree of confinement and distance from the represented recharge boundary.

New York↗

USGS2009

Groundwater-Quality Survey of the South Coast Aquifer of Puerto Rico, April 2 through May 30, 2007

The increased potential for variability of groundwater quality in the South Coast aquifer of Puerto Rico due to saline water encroachment from the Caribbean Sea and from deep parts of the aquifer has become a major concern of water planners and managers. In an effort to determine the extent and sources of this encroachment, the U.S. Geological Survey (USGS) and the Puerto Rico Department of Natural and Environmental Resources conducted a synoptic groundwater-quality survey from April 2 through May 30, 2007, for the South Coast aquifer between Ponce and Arroyo (fig. 1). Groundwater resources in this aquifer extend 150 square miles in south-central Puerto Rico and provide an estimated 44.2 million gallons per day (Mgal/d) or about 61 percent of the total water needs. This amount includes: 15.3 Mgal/d for irrigation, 27.4 Mgal/d for public supply, and 1.5 Mgal/d for industrial and other uses (W.L. Molina-Rivera, U.S. Geological Survey, written commun., 2007). Since 1980 when most of the south coastal plain was intensively cultivated for sugarcane, total groundwater withdrawals have declined about 32 Mgal/d with the greatest decline occurring in irrigation (37.2 Mgal/d) and the greatest increase occurring in public supply (5.5 Mgal/d). Although withdrawals have declined substantially, a major concern is that aquifer recharge provided by irrigation return flow from surface-water irrigation canals has essentially dropped to zero because of the large-scale implementation of groundwater drip irrigation systems.

Scientific Investigations Map↗

USGS2017

Challenges for creating a site-specific groundwater-use record for the Ozark Plateaus aquifer system (central USA) from 1900 to 2010

Hydrologic budgets to determine groundwater availability are important tools for water-resource managers. One challenging component for developing hydrologic budgets is quantifying water use through time because historical and site-specific water-use data can be sparse or poorly documented. This research developed a groundwater-use record for the Ozark Plateaus aquifer system (central USA) from 1900 to 2010 that related county-level aggregated water-use data to site-specific well locations and aquifer units. A simple population-based linear model, constrained to 0 million liters per day in 1900, provided the best means to extrapolate groundwater-withdrawal rates pre-1950s when there was a paucity of water-use data. To disaggregate county-level data to individual wells across a regional aquifer system, a programmatic hierarchical process was developed, based on the level of confidence that a well pumped groundwater for a specific use during a specific year. Statistical models tested on a subset of the best-available site-specific water-use data provided a mechanism to bracket historic groundwater use, such that groundwater-withdrawal rates ranged, on average, plus or minus 38% from modeled values. Groundwater withdrawn for public supply and domestic use accounted for between 48 and 74% of total groundwater use since 1901, highlighting that groundwater provides an important drinking-water resource. The compilation, analysis, and spatial and temporal extrapolation of water-use data remain a challenging task for water scientists, but is of paramount importance to better quantify groundwater use and availability.

Arkansas, Kansas, Missouri, Oklahoma↗

USGS2014

Groundwater movement, recharge, and perchlorate occurrence in a faulted alluvial aquifer in California (USA)

Perchlorate from military, industrial, and legacy agricultural sources is present within an alluvial aquifer in the Rialto-Colton groundwater subbasin, 80 km east of Los Angeles, California (USA). The area is extensively faulted, with water-level differences exceeding 60 m across parts of the Rialto-Colton Fault separating the Rialto-Colton and Chino groundwater subbasins. Coupled well-bore flow and depth-dependent water-quality data show decreases in well yield and changes in water chemistry and isotopic composition, reflecting changing aquifer properties and groundwater recharge sources with depth. Perchlorate movement through some wells under unpumped conditions from shallower to deeper layers underlying mapped plumes was as high as 13 kg/year. Water-level maps suggest potential groundwater movement across the Rialto-Colton Fault through an overlying perched aquifer. Upward flow through a well in the Chino subbasin near the Rialto-Colton Fault suggests potential groundwater movement across the fault through permeable layers within partly consolidated deposits at depth. Although potentially important locally, movement of groundwater from the Rialto-Colton subbasin has not resulted in widespread occurrence of perchlorate within the Chino subbasin. Nitrate and perchlorate concentrations at the water table, associated with legacy agricultural fertilizer use, may be underestimated by data from long-screened wells that mix water from different depths within the aquifer.

California↗

USGS2018

Radium mobility and the age of groundwater in public-drinking-water supplies from the Cambrian-Ordovician aquifer system, north-central USA

High radium (Ra) concentrations in potable portions of the Cambrian-Ordovician (C-O) aquifer system were investigated using water-quality data and environmental tracers ( 3 H, 3 He trit , SF 6 , 14 C and 4 He rad ) of groundwater age from 80 public-supply wells (PSWs). Groundwater ages were estimated by calibration of tracers to lumped parameter models and ranged from modern (<50 yr) in upgradient, regionally unconfined areas to ancient (>1 Myr) in the most downgradient, confined portions of the potable system. More than 80 and 40 percent of mean groundwater ages were older than 1000 and 50,000 yr, respectively. Anoxic, Fe-reducing conditions and increased mineralization develop with time in the aquifer system and mobilize Ra into solution resulting in the frequent occurrence of combined Ra (Ra c = 226 Ra + 228 Ra) at concentrations exceeding the USEPA MCL of 185 mBq/L (5 pCi/L). The distribution of the three Ra isotopes comprising total Ra (Ra t = 224 Ra + 226 Ra + 228 Ra) differed across the aquifer system. The concentrations of 224 Ra and 228 Ra were strongly correlated and comprised a larger proportion of the Ra t concentration in samples from the regionally unconfined area, where arkosic sandstones provide an enhanced source for progeny from the 232 Th decay series. 226 Ra comprised a larger proportion of the Ra t concentration in samples from downgradient confined regions. Concentrations of Ra t were significantly greater in samples from the regionally confined area of the aquifer system because of the increase in 226 Ra concentrations there as compared to the regionally unconfined area. 226 Ra distribution coefficients decreased substantially with anoxic conditions and increasing ionic strength of groundwater (mineralization), indicating that Ra is mobilized to solution from solid phases of the aquifer as adsorption capacity is diminished. The amount of 226 Ra released from solid phases by alpha-recoil mechanisms and retained in solution increases relative to the amount of Ra sequestered by adsorption processes or co-precipitation with barite as adsorption capacity and the concentration of Ba decreases. Although 226 Ra occurred at concentrations greater than 224 Ra or 228 Ra, the ingestion exposure risk was greater for 228 Ra owing to its greater toxicity. In addition, 224 Ra added substantial alpha-particle radioactivity to potable samples from the C-O aquifer system. Thus, monitoring for Ra isotopes and gross-alpha-activity (GAA) is important in upgradient, regionally unconfined areas as downgradient, and GAA measurements made within 72 h of sample collection would best capture alpha-particle radiation from the short-lived 224 Ra.

Applied Geochemistry↗

USGS2018

Simulation of potential groundwater recharge for the glacial aquifer system east of the Rocky Mountains, 1980–2011, using the Soil-Water-Balance Model

An understanding of the spatial and temporal extent of groundwater recharge is critical for many types of hydrologic assessments involving water quality, contaminant transport, ecosystem health, and sustainable use of groundwater. Annual potential groundwater recharge was simulated at a 1-kilometer resolution with the Soil-Water-Balance (SWB) model for the glacial aquifer system east of the Rocky Mountains, from central Montana east to Maine, for calendar years 1980–2011. The SWB model used high resolution meteorological, land cover, and soil hydrology datasets that are nationally consistent and publicly available. The SWB model computed daily potential groundwater recharge as precipitation in excess of interception, runoff, evapotranspiration, and soil-water storage capacity. Daily potential recharge values within each year of the simulation were summed to produce annual potential recharge rates. Potential recharge as described in this report is water that infiltrates vertically below the plant rooting zone and is assumed to reach the water table. The calibrated SWB model in this report is called the glacial SWB model. Model calibration assumed that the area contributing to groundwater discharge equaled the surface watershed. The model was calibrated to stream base flows from 39 watersheds throughout the model domain that had hydrologic conditions appropriate for hydrograph separation. Base flows were calculated from daily streamflow records with the HYSEP local minimum hydrograph separation method The glacial SWB model reproduced the mean annual base-flow calibration targets well; the Nash-Sutcliffe efficiency coefficient was 0.94, and the root mean squared error was 1.28 inches per year. The glacial SWB model provides insight into the spatial and temporal variability in potential annual recharge across the glacial aquifer system. About 20 percent of the active model area had an average potential recharge rate of less than 1 inch per year. Total precipitation, total recharge, and recharge as a percentage of precipitation increased from west to east. A substantial amount of the recharge water (39 percent) entering the glacial aquifer system travels through developed (urbanized) and agricultural landscapes, which are known to cause water-quality impairments. Regional climatic events, such as the 1988 to 1989 drought, are apparent in the potential recharge time series. Potential recharge generally increased across the glacial aquifer system between 2001 and 2011. A comparison of the potential recharge from the glacial SWB model to previous broad-scale recharge estimates reveals several important considerations for future SWB modeling applications. Shifts in the overall distribution of potential recharge between separate models can be explained by methods used to generate base-flow calibration target datasets. Spatial patterns in potential recharge simulated by SWB models are strongly dependent on the data and assumptions used to assign model cells to hydrologic soil groups. A review of several SWB models used to estimate groundwater recharge (and not surface runoff) revealed that model results are most sensitive to input climatic data, followed by surface runoff (curve number) and root-zone depth parameters.

Scientific Investigations Report↗

USGS2024

Groundwater flow model for the Des Moines River alluvial aquifer near Des Moines, Iowa

Des Moines Water Works (DMWW) is a regional municipal water utility that provides residential and commercial water resources to about 600,000 customers in Des Moines, Iowa, and surrounding municipalities in central Iowa. DMWW has identified a need for increased water supply and is exploring the potential for expanding groundwater production capabilities in the Des Moines River alluvial aquifer, where it operates two radial collector wells (RCWs). The U.S. Geological Survey, in cooperation with DMWW, completed a study of the Des Moines River alluvial aquifer and interactions of the RCWs with the aquifer; no previously published model has included the existing well locations, which is the focus of this model. A conceptual and numerical groundwater flow model have been developed to characterize the Des Moines River alluvial aquifer under existing conditions, to simulate water levels observed in the RCWs, and to provide publicly accessible hydrologic data and research that advance understanding of the regional hydrologic system and can potentially be used in the future to evaluate groundwater production scenarios. Model performance was assessed by comparing observed and simulated groundwater levels that included water level elevations, water level changes, water level inequality observations, surface water streamflow, and change in surface water volume from upstream to downstream. Water table elevation in the aquifer layers is on average slightly overestimated with average absolute value error less than 1.5 meters at both RCWs and less than 2.5 meters for all observation wells in the alluvial aquifer layers. The model also accurately simulated water tables greater than the RCW design minimum (a water level threshold at which RCW pumping is reduced) in all timesteps for which water level observation data existed. Water table elevation error was higher in other model layers that were not the focus of the study, and the model did not accurately match streamflow targets.

Iowa↗

USGS2026

Simulation of groundwater flow in the Silurian aquifer, eastern Iowa, 2020–45

The Silurian aquifer is an important water source for municipalities, industry, and rural households and communities in eastern Iowa, including Johnson County. Increasing demand for groundwater from the aquifer indicated that analytical tools may be beneficial to quantify groundwater resources and inform water-management decisions for the aquifer. The U.S. Geological Survey, in cooperation with the Johnson County Board of Supervisors, developed conceptual and numerical groundwater models to simulate water levels in the Silurian aquifer, determine groundwater budgets, and forecast changes in groundwater levels through 2045. A MODFLOW numerical model was constructed to match water levels for 2020 through 2022. Model performance was assessed by comparing simulated and observed water-level hydrographs at several pumping and monitoring wells. Hydrograph comparison indicates that the model simulates water levels accurately at monitoring wells that are not near areas of pumping and generally overpredicts water levels near pumping wells, while accurately simulating the drawdown trend over time at pumping wells. This transient model framework was then used to predict water levels through 2045 in a scenario of moderate drought and increased groundwater withdrawals. Results from the predictive model indicate as much as 13 meters of additional drawdown in 2045 from 2020 water levels in the area of pumping wells.

Iowa↗

USGS2016

Arsenic and radionuclide occurrence and relation to geochemistry in groundwater of the Gulf Coast Aquifer System in Houston, Texas, 2007–11

The U.S. Geological Survey (USGS), in cooperation with the City of Houston, began a study in 2007 to determine concentrations, spatial extent, and associated geochemical conditions that might be conducive for mobility and transport of selected naturally occurring trace elements and radionuclides in the Gulf Coast aquifer system in Houston, Texas. Water samples were collected from 91 municipal supply wells completed in the Evangeline and Chicot aquifers of the Gulf Coast aquifer system in northeastern, northwestern, and southwestern Houston; hereinafter referred to as northeast, northwest and southwest Houston areas. Wells were sampled in three phases: (1) 28 municipal supply wells were sampled during 2007–8, (2) 60 municipal supply wells during 2010, and (3) 3 municipal supply wells during December 2011. During each phase of sampling, samples were analyzed for major ions, selected trace elements, and radionuclides. At a subset of wells, concentrations of arsenic species and other radionuclides (carbon-14, radium-226, radium-228, radon-222, and tritium) also were analyzed. Selected physicochemical properties were measured in the field at the time each sample was collected, and oxidation-reduction potential and unfiltered sulfides also were measured at selected wells. The source-water (the raw, ambient water withdrawn from municipal supply wells prior to water treatment) samples were collected for assessment of aquifer conditions in order to provide community water-system operators information that could be important when they make decisions about which treatment processes to apply before distributing finished drinking water. Geochemical conditions of groundwater of the Gulf Coast aquifer system are suitable in some instances for release of arsenic and radionuclides from aquifer materials. Recent changes to the U.S. Environmental Protection Agency (EPA) primary drinking-water regulations for arsenic and a selected number of natural radionuclides have highlighted the necessity for municipal supply system managers to be aware of the occurrence and distribution of these constituents in their source water. Concentrations of arsenic ranged from 0.58 to 23.5 micrograms per liter (μg/L), with relatively low median and 75th percentile concentrations (2.7 and 3.6 μg/L, respectively). The gross alpha-particle activity completed within 72 hours after sample collection ranged from R-1.1 (nondetect where the result was below the sample specific critical level) to 39.7 picocuries per liter (pCi/L), with a median of 10.3 pCi/L. After 30 days, the gross alpha-particle activities in the 91 samples ranged from R-0.94 to 25.5 pCi/L, with a median of 5.60 pCi/L. Concentrations of uranium ranged from less than 0.02 to 42.7 μg/L, with a median value of 1.69 μg/L and a 75th-percentile value of 6.48 μg/L. The maximum concentrations of radium-226 and combined radium (sum of radium-226 plus radium-228) were 4.34 pCi/L and 3.23 pCi/L, respectively. Aquifer major-ion geochemistry was characterized and shown to contain three chemical types of water as grouped by a simplified predominant cation and anion classification system: (1) calcium- bicarbonate type, (2) sodium-bicarbonate type, and (3) sodium-chloride type. Aquifer geochemistry also was characterized into four reduction-oxidation (redox) categories: (1) oxic, (2) suboxic, (3) mixed, and (4) anoxic. Within the anoxic category, groundwater was further characterized into four presumed predominant reduction processes: (1) iron or sulfate or both [Fe(III)/SO 4 ] reducing, (2) iron [Fe(III)] reducing, (3) iron and sulfate [Fe(III)-SO 4 ] reducing, or (4) methanogenic, as defined by composition of redox species. The oxic category was associated with calcium-bicarbonate-type water, and the methanogenic-anoxic process was associated exclusively with the sodium-bicarbonate-type water. The species of arsenic and the dominant radionuclide present were associated with specific redox categories. Arsenate was associated primarily with oxic water and did not exceed 3.5 µg/L, whereas arsenite was associated with iron-reducing, anoxic water samples and, at the highest concentrations, occurred in sulfate-reducing, anoxic; methanogenic-anoxic; or both water samples. Uranium was associated exclusively with the oxic water, whereas the highest concentrations of combined radium were associated with the iron-reducing, anoxic water. The gross alpha-particle activity was greatest in the oxic waters where the source of the radioactivity was the uranium. Associated geochemical conditions conducive for mobility of arsenic and radionuclides and their spatial and vertical extent in the Gulf Coast aquifer system in Houston are important aspects to the areal management of the municipal groundwater supplies in Houston. Ongoing research is seeking to define chemical or geological factors that are the optimal indicators for elevated concentrations of these naturally occurring constituents.

Texas↗

USGS2014

Climate change impacts on the temperature and magnitude of groundwater discharge from shallow, unconfined aquifers

Cold groundwater discharge to streams and rivers can provide critical thermal refuge for threatened salmonids and other aquatic species during warm summer periods. Climate change may influence groundwater temperature and flow rates, which may in turn impact riverine ecosystems. This study evaluates the potential impact of climate change on the timing, magnitude, and temperature of groundwater discharge from small, unconfined aquifers that undergo seasonal freezing and thawing. Seven downscaled climate scenarios for 2046–2065 were utilized to drive surficial water and energy balance models (HELP3 and ForHyM2) to obtain future projections for daily ground surface temperature and groundwater recharge. These future surface conditions were then applied as boundary conditions to drive subsurface simulations of variably saturated groundwater flow and energy transport. The subsurface simulations were performed with the U.S. Geological Survey finite element model SUTRA that was recently modified to include the dynamic freeze-thaw process. The SUTRA simulations indicate a potential rise in the magnitude (up to 34%) and temperature (up to 3.6°C) of groundwater discharge to the adjacent river during the summer months due to projected increases in air temperature and precipitation. The thermal response of groundwater to climate change is shown to be strongly dependent on the aquifer dimensions. Thus, the simulations demonstrate that the thermal sensitivity of aquifers and baseflow-dominated streams to decadal climate change may be more complex than previously thought. Furthermore, the results indicate that the probability of exceeding critical temperature thresholds within groundwater-sourced thermal refugia may significantly increase under the most extreme climate scenarios.

Water Resources Research↗

USGS2015

Groundwater recharge assessment at local and episodic scale in a soil mantled perched karst aquifer in southern Italy

Groundwater recharge assessment of karst aquifers, at various spatial and temporal scales, is a major scientific topic of current importance, since these aquifers play an essential role for both socio-economic development and fluvial ecosystems. In this study, groundwater recharge was estimated at local and episodic scales in a representative perched karst aquifer in a region of southern Italy with a Mediterranean climate. The research utilized measurements of precipitation, air temperature, soil water content, and water-table depth, obtained in 2008 at the Acqua della Madonna test area (Terminio Mount karst aquifer, Campania region). At this location the aquifer is overlain by ash-fall pyroclastic soils. The Episodic Master Recession (EMR) method, an improved version of the Water Table Fluctuation (WTF) method, was applied to estimate the amount of recharge generated episodically by individual rainfall events. The method also quantifies the amount of precipitation generating each recharge episode, thus permitting calculation of the Recharge to the Precipitation Ratio (RPR) on a storm-by-storm basis. Depending on the seasonally varying air temperature, evapotranspiration, and precipitation patterns, calculated values of RPR varied between 35% and 97% among the individual episodes. A multiple linear correlation of the RPR with both the average intensity of recharging rainfall events and the antecedent soil water content was calculated. Given the relatively easy measurability of precipitation and soil water content, such an empirical model would have great hydrogeological and practical utility. It would facilitate short-term forecasting of recharge in karst aquifers of the Mediterranean region and other aquifers with similar hydrogeological characteristics. By establishing relationships between the RPR and climate-dependent variables such as average storm intensity, it would facilitate prediction of climate-change effects on groundwater recharge. The EMR methodology could further be applied to other aquifers for evaluating the relationship of recharge to various hydrometeorological and hydrogeological processes.

Journal of Hydrology↗

USGS2017

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

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

Journal of Hydrology↗

USGS1995

Transport behavior of groundwater protozoa and protozoan-sized microspheres in sandy aquifer sediments

Transport behaviors of unidentified flagellated protozoa (flagellates) and flagellate-sized carboxylated microspheres in sandy, organically contaminated aquifer sediments were investigated in a small-scale (1 to 4-m travel distance) natural-gradient tracer test on Cape Cod and in flow-through columns packed with sieved (0.5-to 1.0-mm grain size) aquifer sediments. The minute (average in situ cell size, 2 to 3 μ m) flagellates, which are relatively abundant in the Cape Cod aquifer, were isolated from core samples, grown in a grass extract medium, labeled with hydroethidine (a vital eukaryotic stain), and coinjected into aquifer sediments along with bromide, a conservative tracer. The 2- μ m flagellates appeared to be near the optimal size for transport, judging from flowthrough column experiments involving a polydispersed (0.7 to 6.2 μ m in diameter) suspension of carboxylated microspheres. However, immobilization within the aquifer sediments accounted for a log unit reduction over the first meter of travel compared with a log unit reduction over the first 10 m of travel for indigenous, free-living groundwater bacteria in earlier tests. High rates of flagellate immobilization in the presence of aquifer sediments also was observed in the laboratory. However, immobilization rates for the laboratory-grown flagellates (initially 4 to 5 μ m) injected into the aquifer were not constant and decreased noticeably with increasing time and distance of travel. The decrease in propensity for grain surfaces was accompanied by a decrease in cell size, as the flagellates presumably readapted to aquifer conditions. Retardation and apparent dispersion were generally at least twofold greater than those observed earlier for indigenous groundwater bacteria but were much closer to those observed for highly surface active carboxylated latex microspheres. Field and laboratory results suggest that 2- μ m carboxylated microspheres may be useful as analogs in investigating several abiotic aspects of flagellate transport behavior in groundwater.

Massachusetts↗

USGS2010

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

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

Open-File Report↗