USGS Science⌕ Search

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.

At least 991 records · Page 55Linked to original sources

Description and evaluation of the effects of urban and agricultural development on the surficial aquifer system, Palm Beach County, Florida

The surficial aquifer system in Palm Beach County was studied during 1982-85 to determine the effects of increased urban and agricultural development on groundwater levels, flow directions, and quality. The surficial aquifer system and its geologic matrix are divisible into three zones on the bases of relative permeabilities and lithologic characteristics. The two greatest water users in the county, public supply utilities and agricultural irrigators, increased total water withdrawals by 123 and 50%, respectively, during 1970-80. By 1980, 76% of public supply withdrawals were from zones I and II of the surficial aquifer system, whereas groundwater pumpage for irrigation decreased to 9% of the total irrigation water used. Increases in groundwater withdrawals for public supply were greatest in the southeast and central coastal parts of the county and served as an indicator for potential changes of flow directions and water quality in the surficial aquifer system. Residual seawater, emplaced in the aquifer system during the Pleistocene Epoch, is still prevalent in the central and western parts of Palm Beach County where low permeabilities in the geologic matrix have retarded its dilution. Chemical analyses of canal-water and groundwater samples collected in April 1984 were used to evaluate the effects of groundwater/surface water exchange on the quality of water during canal conveyance across the area containing residual seawater.

Florida↗

Noble gas loss may indicate groundwater flow across flow barriers in southern Nevada

Average calculated noble gas temperatures increase from 10 to 22 o C in groundwater from recharge to discharge areas in carbonate-rock aquifers of southern Nevada. Loss of noble gases from groundwater in these regional flow systems at flow barriers is the likely process that produces an increase in recharge noble gas temperatures. Emplacement of low permeability rock into high permeability aquifer rock and the presence of low permeability shear zones reduce aquifer thickness from thousands to tens of meters. At these flow barriers, which are more than 1,000 m lower than the average recharge altitude, noble gases exsolve from the groundwater by inclusion in gas bubbles formed near the barriers because of greatly reduced hydrostatic pressure. However, re-equilibration of noble gases in the groundwater with atmospheric air at the low altitude spring discharge area, at the terminus of the regional flow system, cannot be ruled out. Molecular diffusion is not an important process for removing noble gases from groundwater in the carbonate-rock aquifers because concentration gradients are small.

Conference Paper↗

Water-level altitudes 2004 and water-level changes in the Chicot, Evangeline, and Jasper aquifers and compaction 1973-2003 in the Chicot and Evangeline aquifers, Houston-Galveston region, Texas

This report is one in an annual series of reports that depicts water-level altitudes and water-level changes in the Chicot, Evangeline, and Jasper aquifers, and compaction in the Chicot and Evangeline aquifers in the Houston-Galveston region. The Houston-Galveston region comprises Harris, Galveston, Fort Bend, Waller, and Montgomery Counties and adjacent parts of Brazoria, Grimes, Walker, San Jacinto, Liberty, and Chambers Counties. The report was prepared in cooperation with the Harris-Galveston Coastal Subsidence District, the City of Houston, the Fort Bend Subsidence District, and the Lone Star Groundwater Conservation District. For the Chicot and Evangeline aquifers, maps show approximate water-level altitudes in 2004, water-level changes from 2003 to 2004, approximate water-level changes from 1977 to 2004, and approximate water-level changes from 1990 to 2004 (figs. 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 ). For the Jasper aquifer, maps show approximate water-level altitudes in 2004 and water-level changes from 2003 to 2004 and 2000 to 2004 (figs. 9 , 10 , 11 ). The report also contains a map showing borehole extensometer (well equipped with compaction monitor) site locations ( fig. 12 ) and graphs showing measured compaction of subsurface material at these sites from 1973 or later to 2003 ( fig. 13 ). The U.S. Geological Survey (USGS) has published annual reports of water-level altitudes and water-level changes for the Chicot and Evangeline aquifers in the Houston-Galveston region since 1979; and annual reports of same for the Fort Bend subregion (Fort Bend County and adjacent areas) since 1990. The USGS published its first water-level-altitude map for the Jasper aquifer in the greater Houston area (primarily Montgomery County) in 2001. The 2004 water-level-altitude and water-level-change maps for the three aquifers are included in this report.

Texas↗

Water-level altitudes 2005 and water-level changes in the Chicot, Evangeline, and Jasper aquifers and compaction 1973-2004 in the Chicot and Evangeline aquifers, Houston-Galveston region, Texas

This report is one in an annual series of reports that depicts water-level altitudes and water-level changes in the Chicot, Evangeline, and Jasper aquifers, and compaction in the Chicot and Evangeline aquifers in the Houston-Galveston region. The Houston-Galveston region comprises Harris, Galveston, Fort Bend, Waller, and Montgomery Counties and adjacent parts of Brazoria, Grimes, Walker, San Jacinto, Liberty, and Chambers Counties. The report was prepared in cooperation with the Harris-Galveston Coastal Subsidence District, the City of Houston, the Fort Bend Subsidence District, and the Lone Star Groundwater Conservation District. For the Chicot and Evangeline aquifers, maps show approximate water-level altitudes in 2005, water-level changes from 2004 to 2005, and approximate water-level changes from 2000 to 2005, from 1990 to 2005, and from 1977 to 2005 (figs. 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 , and 10 ). For the Jasper aquifer, maps show approximate water-level altitudes in 2005 and water-level changes from 2004 to 2005 and 2000 to 2005 (figs. 11 , 12 , and 13 ). The report also contains a map showing borehole extensometer (well equipped with compaction monitor) site locations ( fig. 14 ) and graphs showing measured compaction of subsurface material at these sites from 1973 or later to 2004 ( fig. 15 ). The U.S. Geological Survey (USGS) has published annual reports of water-level altitudes and water-level changes for the Chicot and Evangeline aquifers in the Houston-Galveston region since 1979; and annual reports of same for the Fort Bend subregion (Fort Bend County and adjacent areas) since 1990. The USGS published its first water-level-altitude map for the Jasper aquifer in the greater Houston area (primarily Montgomery County) in 2001. The 2005 water-level-altitude and water-level-change maps for the three aquifers are included in this report.

Texas↗

Water-level altitudes 2006 and water-level changes in the Chicot, Evangeline, and Jasper aquifers and compaction 1973-2005 in the Chicot and Evangeline aquifers, Houston-Galveston region, Texas

This report is one in an annual series of reports that depicts water-level altitudes and water-level changes in the Chicot, Evangeline, and Jasper aquifers, and compaction in the Chicot and Evangeline aquifers in the Houston-Galveston region. The Houston-Galveston region comprises Harris, Galveston, Fort Bend, Waller, and Montgomery Counties and adjacent parts of Brazoria, Grimes, Walker, San Jacinto, Liberty, and Chambers Counties. The report was prepared in cooperation with the Harris-Galveston Coastal Subsidence District, the City of Houston, the Fort Bend Subsidence District, and the Lone Star Groundwater Conservation District. For the Chicot and Evangeline aquifers, maps show approximate water-level altitudes in 2006, water-level changes from 2005 to 2006, and approximate water-level changes from 2001 to 2006, from 1990 to 2006, and from 1977 to 2006 (figs. 1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 , 9 , 10 ). For the Jasper aquifer, maps show approximate water-level altitudes in 2006 and water-level changes from 2005 to 2006 and 2000 to 2006 (figs. 11 , 12 , 13 ). The report also contains a map showing borehole extensometer (well equipped with compaction monitor) site locations ( fig. 14 ) and graphs showing measured compaction of subsurface material at these sites from 1973 or later to 2005 ( fig. 15 ). The U.S. Geological Survey (USGS) has published annual reports of water-level altitudes and water-level changes for the Chicot and Evangeline aquifers in the Houston-Galveston region since 1979; and annual reports of same for the Fort Bend subregion (Fort Bend County and adjacent areas) since 1990. The USGS published its first water-level-altitude map for the Jasper aquifer in the greater Houston area (primarily Montgomery County) in 2001. The 2006 water-level-altitude and water-level-change maps for the three aquifers are included in this report.

Texas↗

Simulation of groundwater-level changes from projected groundwater withdrawals in the Truxton basin, northwestern Arizona

A three-dimensional, numerical groundwater flow model of the Hualapai Plateau and Truxton basin was developed to assist water-resource managers in understanding the potential effects of projected groundwater withdrawals on groundwater levels and storage in the basin. The Truxton Basin Hydrologic Model (TBHM) is a transient model that simulates the hydrologic system for the years 1976 through 2139, including hypothetical low-, medium-, and high-groundwater withdrawal scenarios beginning in 2020. The simulated effects of these withdrawal scenarios are presented as groundwater-level changes from the year 2020 to 2070, and from 2020 to 2140. Hydrologic properties in the TBHM are derived from calibration of a steady-state model of the predevelopment (before 1976) groundwater system. The future pumping scenarios are each simulated with three different interpretations of basin depth supported by geophysical data. For each of the resulting nine transient models, a Monte Carlo approach is used to produce a range of possible and probable groundwater-level changes at points throughout the basin given probabilistic ranges of hydrologically reasonable aquifer property values supported by the model calibration results. The ensemble of models that simulate the future pumping scenarios include pumping from the existing well field (three wells) plus additional pumping from a proposed new well. Simulated high future pumping increases progressively to 1,840 acre-feet per year in 2120 and produces a range of drawdowns between 20 and 39 feet (ft) near the pumping center, with a median drawdown of 28 ft. The low future pumping scenario, which increases progressively to 650 acre-ft per year in 2120, produces a range of drawdowns between 5 and 15 ft, with a median drawdown of 10 ft at the same location over the same period of time.

Arizona↗

A model of transmissivity and hydraulic conductivity from electrical resistivity distribution derived from airborne electromagnetic surveys of the Mississippi River Valley Alluvial Aquifer, Midwest USA

Groundwater-flow models require the spatial distribution of the hydraulic conductivity parameter. One approach to defining this spatial distribution in groundwater-flow model grids is to map the electrical resistivity distribution by airborne electromagnetic (AEM) survey and establish a petrophysical relation between mean resistivity calculated as a nonlinear function of the resistivity layering and thicknesses of the layers and aquifer transmissivity compiled from historical aquifer tests completed within the AEM survey area. The petrophysical relation is used to transform AEM resistivity to transmissivity and to hydraulic conductivity over areas where the saturated thickness of the aquifer is known. The US Geological Survey applied this approach to a gain better understanding of the aquifer properties of the Mississippi River Valley alluvial aquifer. Alluvial-aquifer transmissivity data, compiled from 160 historical aquifer tests in the Mississippi Alluvial Plain (MAP), were correlated to mean resistivity calculated from 16,816 line-kilometers (km) of inverted resistivity soundings produced from a frequency-domain AEM survey of 95,000 km 2 of the MAP. Correlated data were used to define petrophysical relations between transmissivity and mean resistivity by omitting from the correlations the aquifer-test and AEM sounding data that were separated by distances greater than 1 km and manually calibrating the relation coefficients to slug-test data. The petrophysical relation yielding the minimum residual error between simulated and slug-test data was applied to 2,364 line-km of AEM soundings in the 1,000-km 2 Shellmound (Mississippi) study area to calculate hydraulic property distributions of the alluvial aquifer for use in future groundwater-flow models.

Mississippi River Valley Alluvial Aquifer↗

Infiltration and hydraulic connections from the Niagara River to a fractured-dolomite aquifer in Niagara Falls, New York

The spatial distribution of hydrogen and oxygen stable-isotope values in groundwater can be used to distinguish different sources of recharge and to trace groundwater flow directions from recharge boundaries. This method can be particularly useful in fractured-rock settings where multiple lines of evidence are required to delineate preferential flow paths that result from heterogeneity within fracture zones. Flow paths delineated with stable isotopes can be combined with hydraulic data to form a more complete picture of the groundwater flow system. In this study values of ??D and ??18O were used to delineate paths of river-water infiltration into the Lockport Group, a fractured dolomite aquifer, and to compute the percentage of fiver water in groundwater samples from shallow bedrock wells. Flow paths were correlated with areas of high hydraulic diffusivity in the shallow bedrock that were delineated from water-level fluctuations induced by diurnal stage fluctuations in man-made hydraulic structures. Flow paths delineated with the stable-isotope and hydraulic data suggest that fiver infiltration reaches an unlined storm sewer in the bedrock through a drainage system that surrounds aqueducts carrying river water to hydroelectric power plants. This finding is significant because the storm sewer is the discharge point for contaminated groundwater from several chemical waste-disposal sites and the cost of treating the storm sewer's discharge could be reduced if the volume of infiltration from the river were decreased.The spatial distribution of hydrogen and oxygen stable-isotope values in groundwater can be used to distinguish different sources of recharge and to trace groundwater flow directions from recharge boundaries. This method can be particularly useful in fractured-rock settings where multiple lines of evidence are required to delineate preferential flow paths that result from heterogeneity within fracture zones. Flow paths delineated with stable isotopes can be combined with hydraulic data to form a more complete picture of the groundwater flow system. In this study values of ??D and ??18O were used to delineate paths of river-water infiltration into the Lockport Group, a fractured dolomite aquifer, and to compute the percentage of river water in groundwater samples from shallow bedrock wells. Flow paths were correlated with areas of high hydraulic diffusivity in the shallow bedrock that were delineated from water-level fluctuations induced by diurnal stage fluctuations in man-made hydraulic structures. Flow paths delineated with the stable-isotope and hydraulic data suggest that river infiltration reaches an unlined storm sewer in the bedrock through a drainage system that surrounds aqueducts carrying river water to hydroelectric power plants. This finding is significant because the storm sewer is the discharge point for contaminated groundwater from several chemical waste-disposal sites and the cost of treating the storm sewer's discharge could be reduced if the volume of infiltration from the river were decreased.

Journal of Hydrology↗

Changes in groundwater flow and volatile organic compound concentrations at the Fischer and Porter Superfund Site, Warminster Township, Bucks County, Pennsylvania, 1993-2009

The 38-acre Fischer and Porter Company Superfund Site is in Warminster Township, Bucks County, Pa. Historically, as part of the manufacturing process, trichloroethylene (TCE) degreasers were used for parts cleaning. In 1979, the Bucks County Health Department detected TCE and other volatile organic compounds (VOCs) in water from the Fischer and Porter on-site supply wells and nearby public-supply wells. The Fischer and Porter Site was designated as a Superfund Site and placed on the National Priorities List in September 1983. A 1984 Record of Decision for the site required the Fischer and Porter Company to pump and treat groundwater contaminated by VOCs from three on-site wells at a combined rate of 75 gallons per minute to contain groundwater contamination on the property. Additionally, the Record of Decision recognized the need for treatment of the water from two nearby privately owned supply wells operated by the Warminster Heights Home Ownership Association. In 2004, the Warminster Heights Home Ownership Association sold its water distribution system, and both wells were taken out of service. The report describes changes in groundwater levels and contaminant concentrations and migration caused by the shutdown of the Warminster Heights supply wells and presents a delineation of the off-site groundwater-contamination plume. The U.S. Geological Survey (USGS) conducted this study (2006-09) in cooperation with the U.S. Environmental Protection Agency (USEPA). The Fischer and Porter Site and surrounding area are underlain by sedimentary rocks of the Stockton Formation of Late Triassic age. The rocks are chiefly interbedded arkosic sandstone and siltstone. The Stockton aquifer system is comprised of a series of gently dipping lithologic units with different hydraulic properties. A three-dimensional lithostratigraphic model was developed for the site on the basis of rock cores and borehole geophysical logs. The model was simplified by combining individual lithologic units into generalized units representing upward fining sedimentary cycles capped by a siltstone bed. These cycles were labeled units 1 through 8 and are called stratigraphic units in this report. Groundwater in the unweathered zone mainly moves through a network of interconnecting secondary openings--bedding-plane fractures and joints. Groundwater generally is unconfined in the shallower part of the aquifer and confined or semiconfined in the deeper part of the aquifer. The migration of VOCs from the Fischer and Porter Site source area is influenced by geologic and hydrologic controls. The hydrologic controls have changed with time. Stratigraphic units 2 and 3 crop out beneath the former Fischer and Porter plant. VOCs originating at the plant source area entered these stratigraphic units and moved downdip to the northwest. When the wells at and in the vicinity of the site were initially sampled in 1979-80, three public-supply wells (BK-366, BK-367, MG-946) and three industrial-supply wells (BK-368, BK-370, and BK-371) were pumping. Groundwater contaminated with VOCs flowed downdip and then northeast along strike toward well BK-366, downdip toward well BK-368, and downdip and then west along strike toward well MG-946. The long axis of the TCE plume is oriented about N. 18° W. in the direction of dip. In 1979-80, the leading edge of the plume was about 3,500 feet wide. With the cessation of pumping of the supply wells in 2004, the size of the plume has decreased. In 2007-09, the plume was approximately 2,000 feet long and 2,000 feet wide at the leading edge. On the western side of the site, TCE and tetrachloroethylene (PCE) appear to be moving downdip though stratigraphic unit 3. The downdip extent of TCE and PCE migration extended approximately 550 feet off-site to the northwest and 750 feet off-site to the north. TCE concentrations in water samples from wells at the western site boundary increased from 1996 to 2007. On the northern side of the site, TCE and PCE appeared to be moving downward and laterally though stratigraphic units 2, 3, and 4. Groundwater-flow directions shifted to the northwest in the intermediate and deep zones after cessation of pumping of well BK-366 in 2004. The shutdown of the Warminster Heights wells had little effect on the direction of groundwater flow in the shallow zone. In 2007, TCE concentrations measured in water samples from the three remediation wells by the USGS ranged from less than 340 to 3,000 µg/L, and PCE concentrations ranged from less than 8.4 to 51 µg/L. TCE concentrations in water samples from the source-area remediation wells have decreased with time but remain highly variable. From 2001 to 2008, the TCE and PCE concentrations in water samples from wells BK-370 and BK-371 showed a linear decreasing trend. TCE and PCE concentrations in water samples from well BK-1324 showed an exponentially decreasing trend. In 2007, TCE concentrations measured in water samples from shallow wells ranged from less than 0.1 to 14,000 µg/L, and PCE concentrations ranged from less than 0.1 to 340 µg/L. The TCE and PCE plumes followed the hydraulic gradient in the shallow zone. In 2007, TCE concentrations measured in water samples from on-site intermediate-depth monitor wells ranged from less than 0.1 to 500 µg/L, and PCE concentrations ranged from 1.3 to 28 µg/L. The TCE and PCE plumes followed the hydraulic gradient in the intermediate zone and extended off-site to the north and northwest of the source area. Concentrations of TCE in water samples north and west of the source area increased from 1996 to 2007. In 2007, the TCE concentrations measured in water samples from on-site monitor wells in the deep zone ranged from 1.1 to 86 µg/L, and PCE concentrations ranged from less than 0.1 to 8.4 µg/L. The TCE and PCE plumes generally followed the hydraulic gradient in the deep zone and extended off-site to the northwest of the source area. In general, concentrations of TCE in water samples from monitor wells outside the source area increased between 1996 and 2005 and decreased between 2005 and 2007; concentrations were less in 2007 than in 1996.

Pennsylvania↗

Shallow groundwater quality and geochemistry in the Fayetteville Shale gas-production area, north-central Arkansas, 2011

The Mississippian Fayetteville Shale serves as an unconventional gas reservoir across north-central Arkansas, ranging in thickness from approximately 50 to 550 feet and varying in depth from approximately 1,500 to 6,500 feet below the ground surface. Primary permeability in the Fayetteville Shale is severely limited, and successful extraction of the gas reservoir is the result of advances in horizontal drilling techniques and hydraulic fracturing to enhance and develop secondary fracture porosity and permeability. Drilling and production of gas wells began in 2004, with a steady increase in production thereafter. As of April 2012, approximately 4,000 producing wells had been completed in the Fayetteville Shale. In Van Buren and Faulkner Counties, 127 domestic water wells were sampled and analyzed for major ions and trace metals, with a subset of the samples analyzed for methane and carbon isotopes to describe general water quality and geochemistry and to investigate the potential effects of gas-production activities on shallow groundwater in the study area. Water-quality analyses from this study were compared to historical (pregas development) shallow groundwater quality collected in the gas-production area. An additional comparison was made using analyses from this study of groundwater quality in similar geologic and topographic areas for well sites less than and greater than 2 miles from active gas-production wells. Chloride concentrations for the 127 groundwater samples collected for this study ranged from approximately 1.0 milligram per liter (mg/L) to 70 mg/L, with a median concentration of 3.7 mg/L, as compared to maximum and median concentrations for the historical data of 378 mg/L and 20 mg/L, respectively. Statistical analysis of the data sets revealed statistically larger chloride concentrations (p-value <0.001) in the historical data compared to data collected for this study. Chloride serves as an important indicator parameter based on its conservative transport characteristics and relatively elevated concentrations in production waters associated with gas extraction activities. Major ions and trace metals additionally had lower concentrations in data gathered for this study than in the historical analyses. Additionally, no statistical difference existed between chloride concentrations from water-quality data collected for this study from 94 wells located less than 2 miles from a gas-production well and 33 wells located 2 miles or more from a gas-production well; a Wilcoxon rank-sum test showed a p-value of 0.71. Major ion chemistry was investigated to understand the effects of geochemical and reduction-oxidation (redox) processes on the shallow groundwater in the study area along a continuum of increased rock-water interaction represented by increases in dissolved solids concentration. Groundwater in sandstone formations is represented by a low dissolved solids concentration (less than 30 mg/L) and slightly acidic water type. Shallow shale aquifers were represented by dissolved solids concentrations ranging upward to 686 mg/L, and water types evolving from a dominantly mixed-bicarbonate and calcium-bicarbonate to a strongly sodium-bicarbonate water type. Methane concentration and carbon isotopic composition were analyzed in 51 of the 127 samples collected for this study. Methane occurred above a detection limit of 0.0002 mg/L in 32 of the 51 samples, with concentrations ranging upward to 28.5 mg/L. Seven samples had methane concentrations greater than or equal to 0.5 mg/L. The carbon isotopic composition of these higher concentration samples, including the highest concentration of 28.5 mg/L, shows the methane was likely biogenic in origin with carbon isotope ratio values ranging from -57.6 to -74.7 per mil. Methane concentrations increased with increases in dissolved solids concentrations, indicating more strongly reducing conditions with increasing rock-water interaction in the aquifer. As such, groundwater-quality data collected for this study indicate that groundwater chemistry in the shallow aquifer system in the study area is a result of natural processes, beginning with recharge of dilute atmospheric precipitation and evolution of observed groundwater chemistry through rock-water interaction and redox processes.

Arkansas↗

Status of water-level altitudes and long-term and short-term water-level changes in the Chicot and Evangeline (undifferentiated) and Jasper aquifers, greater Houston area, Texas, 2023

Since the early 1900s, groundwater withdrawn from the primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers—has been an important source of water in the greater Houston area, Texas. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of water-level altitudes and water-level changes in these aquifers in the greater Houston area. In this report, the Chicot and Evangeline aquifers are treated as a single aquifer for the purposes of providing annual assessments of regional-scale water-level altitudes and water-level changes over time. In 2023, shaded depictions of water-level altitudes for the Chicot and Evangeline aquifers (undifferentiated) ranged from about 286 feet (ft) below the North American Vertical Datum of 1988 (NAVD 88) to about 169 ft above NAVD 88. The largest decline in water-level altitudes indicated by the 1977–2023 long-term water-level-change map was in south-central Montgomery County southeast of The Woodlands. In comparison, the 1990–2023 long-term water-level-change map depicts the largest declines in water-level altitudes in localized areas at or near certain wells in parts of northwestern Harris County and south-central Montgomery County. The largest rise in water-level altitudes for 1977–2023 is depicted in a relatively large area in southeastern Harris County, whereas the largest rise in water-level altitudes for 1990–2023 is depicted in a relatively large area in central Harris County. The 5-year short-term water-level-change map depicts the largest declines at three wells in northern Fort Bend County, one well in western Harris County, and three wells in south-central Montgomery County and the largest rise at one well in central Harris County. The 1-year short-term water-level-change map depicts the largest declines at one well in northern Fort Bend County and two wells in southwestern Harris County and the largest rises at one well in northern Brazoria County and one well in south-central Montgomery County. In 2023, shaded depictions of water-level altitudes for the Jasper aquifer ranged from about 242 ft below NAVD 88 to about 218 ft above NAVD 88. The 2000–23 long-term water-level-change map depicts water-level declines throughout the study area where water-level-measurement data from the aquifer were collected, with the largest declines in north-central Harris County and south-central Montgomery County south of The Woodlands. The 5-year short-term water-level-change map depicts the largest declines at two wells in central Montgomery County near Conroe and two wells in south-central Montgomery County southeast of The Woodlands and the largest rise at one well in western Montgomery County. The 1-year short-term water-level-change map depicts the largest declines at four wells in south-central Montgomery County southeast of The Woodlands and one well in central Montgomery County near Conroe and the largest rises at two wells in western Montgomery County.

Texas↗

Recharge of shallow aquifers through two ephemeral-stream channels in northeastern Wyoming, 1982-83

Quantifying the recharge from ephemeral streams to alluvial and bedrock aquifers will help evaluate the effects of surface mining on alluvial valley floors in Wyoming. Two stream reaches were chosen for study in the Powder River basin. One reach was located along the North Fork Dry Fork Cheyenne River near Glenrock, Wyoming, and the other reach was located along Black Thunder Creek near Hampshire, Wyoming. The reach along the North Fork Dry Fork Cheyenne River was instrumented with 3 gaging stations to measure streamflow and with 6 observation wells to measure groundwater level fluctuations in alluvial and bedrock aquifers in response to streamflow. The 3 streamflow gaging stations were located within the 2.5-mi study reach to measure the approximate gain or loss of discharge along the reach. Computed streamflow losses ranged from 0.43 acre-ft/mi on July 9 , 1982, to 1.44 acre-ft/mi on August 9, 1982. The observation wells completed only in the alluvial aquifer were dry during flow in the North Fork Dry Fork Cheyenne River, whereas water levels in half of the observation wells completed in the bedrock aquifers or the alluvial and bedrock aquifers rose in response to flow in the North Fork Dry Fork Cheyenne River. Groundwater recharge on August 9, 1982, was calculated using a convolution technique using groundwater levels at the upstream site and was estimated to be 26.5 acre-ft/mi. The reach along Black Thunder Creek was instrumented with one gaging station to measure streamflow and with 4 observation wells to measure water level response in alluvial and bedrock aquifers to streamflow. Recharge to the alluvial aquifer from flow in Black Thunder Creek ranged from 3.56 to 12.4 acre-ft/mi. The recharge was estimated using the convolution technique using water level measurements in the observation wells completed in the alluvial aquifer. Water level measurements in the observation wells indicated water level rises in the alluvial and bedrock aquifers in response to flow in Black Thunder Creek.

Wyoming↗

A combined radio- and stable-isotopic study of a California coastal aquifer system

Stable and radioactive tracers were utilized in concert to characterize geochemical processes in a complex coastal groundwater system and to provide constraints on the kinetics of rock/water interactions. Groundwater samples from wells within the Dominguez Gap region of Los Angeles County, California were analyzed for a suite of major cations (Na + , K + , Mg 2+ , Ca 2+ ) and anions (Cl − , SO 4 2− ), silica, alkalinity, select trace elements (Ba, B, Sr), dissolved oxygen, stable isotopes of hydrogen (δD), oxygen (δ 18 O), dissolved inorganic carbon (δ 13 C DIC ), and radioactive isotopes ( 3 H, 222 Rn and 223,224,226,228 Ra). In the study area, groundwater may consist of a complex mixture of native groundwater, intruded seawater, non-native injected water, and oil-field brine water. In some wells, Cl − concentrations attained seawater-like values and in conjunction with isotopically heavier δ 18 O values, these tracers provide information on the extent of seawater intrusion and/or mixing with oil-field brines. Groundwater 3 H above 1 tritium unit (TU) was observed only in a few select wells close to the Dominguez Gap area and most other well groundwater was aged pre-1952. Based on an initial 14 C value for the study site of 90 percent modern carbon (pmc), groundwater age estimates likely extend beyond 20 kyr before present and confirm deep circulation of some native groundwater through multiple aquifers. Enriched values of groundwater δ 13 C DIC in the absence of SO 4 2− imply enhanced anaerobic microbial methanogenesis. While secular equilibrium was observed for 234 U/ 238 U (activity ratios ~1) in host matrices, strong isotopic fractionation in these groundwater samples can be used to obtain information of adsorption/desorption kinetics. Calculated Ra residence times are short, and the associated desorption rate constant is about three orders of magnitude slower than that of the adsorption rate constant. Combined stable- and radio-isotopic results provide unique insights into aquifer characteristics, such as geochemical cycling, rock/water interactions, and subsurface transport and mixing.

California↗

Source, movement and age of groundwater in the upper part of the Mojave River Basin, California, USA

Water samples from wells were collected and analysed for oxygen-18, deuterium, tritium, carbon-14, and carbon-13 to determine the source, movement and age of groundwater in the upper part of the Mojave River basin. Water in the alluvial aquifer has a median deuterium composition of -66??? and contains tritium, and was recently recharged by water from the Mojave River. Water in the regional aquifer near the Mojave River, near Summit Valley, and underlying several small washes has deuterium compositions heavier than -60???. Although some water in the regional aquifer near the Mojave River contains tritium, most of this water does not contain tritium. Carbon-14 data indicate that this water was recharged less than 2400 years ago. Water in the remainder of the regional aquifer has a median deuterium composition of -84???, which is as much as 20??? lighter than the volume-weighted deuterium composition of present-day precipitation. These data show that this water was recharged under climatic conditions different from average conditions today. Carbon-14 data indicate that some water in the regional aquifer was recharged more than 20 000 years ago.Water samples from wells were collected and analyzed for oxygen-18, deuterium, tritium, carbon-14, and carbon-13 to determine the source, movement and age of groundwater in the upper part of the Mojave River basin. Water in the alluvial aquifer has a median deuterium composition of -66qq and contains tritium, and was recently recharged by water from the Mojave River. Water in the regional aquifer near the Mojave River, near Summit Valley, and underlying several small washes has deuterium compositions heavier than -60qq. Although some water in the regional aquifer near the Mojave River contains tritium, most of this water does not contain tritium. Carbon-14 data indicate that this water was recharged less than 2400 years ago. Water in the remainder of the regional aquifer has a median deuterium composition of -84qq, which is as much as 20qq lighter than the volume-weighted deuterium composition of present-day precipitation. These data show that this water was recharged under climatic conditions different from average conditions today. Carbon-14 data indicate that some water in the regional aquifer was recharged more than 20 000 years ago.

IAHS-AISH Publication↗

Groundwater modeling

The state of the science and practice in groundwater modeling brings to mind highly sophisticated computer models that are running in parallel on many multi-processor machines. These models are expected to incorporate many different processes of both saturated and unsaturated groundwater flow and transport and possibly the media to which it connects, like surface waters and the atmosphere. We are increasingly aware we cannot study groundwater flow in isolation if we are to make useful predictions of, for instance, the impacts of climate change on the groundwater regime. We have come a long way. Today we are no longer limited to equations for flow toward a well, perhaps near an infinitely long straight canal (method of images), to sandbox models in the laboratory, or to simple steady state models of flow in a single aquifer. We now have computer models that solve groundwater flow and transport in multi-aquifer settings under transient conditions and with a user-friendly graphical user interface that allows widespread use. Additionally, multi-media models are now leaving the research environment and becoming available to mainstream consultants. So in that sense the science of groundwater modeling has matured. The practice of groundwater modeling, however, has also matured. We have come to realize that model output, being a necessary simplification of an unknowably complex natural world, has inherent limitations. That is, a model of reality is not reality itself. There is uncertainty associated with all facets of our model—parameterization, aquifer geometry and discretization, boundary conditions, and future hydrologic drivers such as future pumping regimes and climates. Today a model is now more appropriately seen as a tool that provides a quantitative framework to make supportable forecasts rather than an oracle that gives us all the answers. In this chapter we set out to briefly review the state of the science and practice in modeling. In doing so, we augment existing assessments from the journal Groundwater (e.g., Hunt and Zheng 2012; Langevin and Panday 2012; Molz 2017a,b; White 2017), specifically in terms of modeling approach. An effective modeling approach is critical. If a modeler does not decompose the societal problem correctly, the model will not be fit-for-purpose, no matter how sophisticated the code’s capabilities. Moreover, capabilities of codes will be ever improving; good modeling practices have a timelessness that is more robust. How best to decompose the problem and provide models that are accepted? We lay out here some approaches for today’s applied groundwater modeling. Specifically, we suggest: (1) a step-wise modeling process; (2) including a two-dimensional areal model within this process; (3) keeping abreast of industry standards; and (4) ways to increase acceptance of the models we produce.

Book chapter↗

Long-term groundwater depletion in the United States

The volume of groundwater stored in the subsurface in the United States decreased by almost 1000&thinsp;km 3 during 1900&ndash;2008. The aquifer systems with the three largest volumes of storage depletion include the High Plains aquifer, the Mississippi Embayment section of the Gulf Coastal Plain aquifer system, and the Central Valley of California. Depletion rates accelerated during 1945&ndash;1960, averaging 13.6&thinsp;km 3 /year during the last half of the century, and after 2000 increased again to about 24&thinsp;km 3 /year. Depletion intensity is a new parameter, introduced here, to provide a more consistent basis for comparing storage depletion problems among various aquifers by factoring in time and areal extent of the aquifer. During 2001&ndash;2008, the Central Valley of California had the largest depletion intensity. Groundwater depletion in the United States can explain 1.4% of observed sea-level rise during the 108-year study period and 2.1% during 2001&ndash;2008. Groundwater depletion must be confronted on local and regional scales to help reduce demand (primarily in irrigated agriculture) and/or increase supply.

Groundwater↗

Improved method for simulating groundwater inundation using the MODFLOW 6 Lake Transport Package

Groundwater inundation due to sea level rise can affect island and coastal freshwater resources by exposing water tables to direct, continuous evaporation. Numerical simulations of groundwater inundation effects on coastal and island aquifers have been limited by an inability to simulate solute transport and variable density flow between the aquifer and lakes formed by groundwater inundation. Consequently, we contributed to the development of a new tool, the Lake Transport Package, for MODFLOW 6 that can calculate solute concentrations within lakes and allows for variable density flow between lakes and aquifers. Here we use groundwater inundation as an example application to showcase the functionality of the Lake Transport Package and the advantages of using this tool over past methods of representing groundwater inundation. We developed hypothetical island simulations based on hydrogeological characteristics of the Bahamas. Multiple sea level rise and lake evaporation rates were simulated to evaluate the effects of groundwater inundation on freshwater lens size for different climates. The results demonstrate the ability of the Lake Transport Package to calculate the solute concentration of the lake for transient simulations, including hypersaline concentrations. Higher sea level rise and greater lake evaporation rates lead to a greater loss of the freshwater lens and higher lake salinity. The formation of a lake and corresponding expansion due to groundwater inundation increases the loss of freshwater by 6–36%, depending on the lake evaporation rate. These simulations validate the performance and demonstrate usefulness of the Lake Transport Package as a tool in representing groundwater inundation.

Groundwater↗

Modeling future groundwater depletion to evaluate sustainability goals set under the Sustainable Groundwater Management Act in the critically overdrafted basins of the Central Valley, California, USA (2020–2070)

In 2014, California's Sustainable Groundwater Management Act (SGMA) mandated local agencies to devise and implement groundwater sustainability plans to address critically overdrafted conditions throughout the state's aquifers. However, the feasibility of these agencies' sustainability goals has not previously been assessed through a regional-scale, integrative lens. Here, we develop and analyze a novel, basin-wide database of 936 sustainability indicator wells located within Central Valley subbasins designated as critically overdrafted, most of which lie in the San Joaquin Valley. Our database shows 2040 groundwater elevation goals vary widely from 60 m above to 80 m below 2020 levels, with variability within and between adjacent subbasins. To evaluate the feasibility of achieving these goals, we coupled the database with a regional hydrologic model (Central Valley Hydrologic Model version 2) and simulated multiple future pumping scenarios. Results show that under increased groundwater demand, 60%–70% of indicator wells may fail to meet their 2040 goals. Even a 50% reduction from 2020 demand levels leaves nearly 40% of wells failing to meet their sustainability thresholds by 2040. Baseline models show that by 2070, up to 70% of wells could fail to meet their goals due to large-scale, spatially connected regions of groundwater depletion. This integrated framework, linking the first region-wide compilation of SGMA indicator wells with a regional groundwater model, demonstrates that many local sustainability goals may be unattainable with substantial (up to 50%) reductions in pumping. Additional management interventions, such as expanded recharge or coordinated demand reductions, may help achieve sustainability goals.

California↗