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 1,693 records · Page 94Linked to original sources

USGS2004

Conceptualization and simulation of the Edwards aquifer, San Antonio region, Texas

A new numerical ground-water-flow model (Edwards aquifer model) that incorporates important components of the latest information and plausible conceptualization of the Edwards aquifer was developed. The model includes both the San Antonio and Barton Springs segments of the Edwards aquifer in the San Antonio region, Texas, and was calibrated for steady-state (1939–46) and transient (1947–2000) conditions, excluding Travis County. Transient simulations were conducted using monthly recharge and pumpage (withdrawal) data. The model incorporates conduits simulated as continuously connected (other than being separated in eastern Uvalde and southwestern Medina Counties), one-cell-wide (1,320 feet) zones with very large hydraulic-conductivity values (as much as 300,000 feet per day). The locations of the conduits were based on a number of factors, including major potentiometric-surface troughs in the aquifer, the presence of sinking streams, geochemical information, and geologic structures (for example, faults and grabens). The simulated directions of flow in the Edwards aquifer model are most strongly influenced by the presence of simulated conduits and barrier faults. The simulated flow in the Edwards aquifer is influenced by the locations of the simulated conduits, which tend to facilitate flow. The simulated subregional flow directions generally are toward the nearest conduit and subsequently along the conduits from the recharge zone into the confined zone and toward the major springs. Structures simulated in the Edwards aquifer model influencing groundwater flow that tend to restrict flow are barrier faults. The influence of simulated barrier faults on flow directions is most evident in northern Medina County. A water budget is an accounting of inflow to, outflow from, and storage change in the aquifer. For the Edwards aquifer model steady-state simulation, recharge (from seepage losses from streams and infiltration of rainfall) accounts for 93.5 percent of the sources of water to the Edwards aquifer, and inflow through the northern and northwestern model boundaries contributes 6.5 percent. The largest discharges are springflow (73.7 percent) and ground-water withdrawals by wells (25.7 percent). The principal source of water to the Edwards aquifer for the Edwards aquifer model transient simulation was recharge, constituting about 60 percent of the sources of water (excluding change in storage) to the Edwards aquifer during 1956, a drought period, and about 97 percent of the sources (excluding change in storage) during 1975, a period of above-normal rainfall and recharge. The principal discharges from the Edwards aquifer for the transient simulation were springflow and withdrawals by wells. During 1956, representing drought conditions, the change in storage (net water released from storage) was much greater than recharge, accounting for 75.9 percent of the total flow compared to 14.5 percent for recharge. Conversely, during 1975, representing above-normal rainfall and recharge conditions, recharge constituted 79.9 percent of the total flow, compared to 7.1 percent for the change in storage (net water added to storage). A series of sensitivity tests was made to ascertain how the model results were affected by variations greater than and less than the calibrated values of input data. Simulated hydraulic heads in the Edwards aquifer model were most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, and specific yield and were comparatively insensitive to spring-orifice conductance, northern boundary inflow, and specific storage. Simulated springflow in the Edwards aquifer model was most sensitive to recharge, withdrawals, hydraulic conductivity of the conduit segments, specific yield, and increases in northern boundary inflow and was comparatively insensitive to spring-orifice conductance and specific storage.

Texas↗

USGS2020

Trends in recent historical and projected climate data for the Colorado River Basin and potential effects on groundwater availability

Understanding recent historical and projected trends in precipitation and temperature in the Colorado River Basin, and estimating what the projected changes in these climate parameters may mean for groundwater resources in the region, is important for water managers and policymakers to sustainably manage water resources in the basin. Historical (1896–2019) precipitation and temperature data for the upper and lower Colorado River Basins were analyzed to better understand recent trends in climate data that may affect groundwater resources in the area. Historical data indicate multidecadal-scale cyclical patterns in precipitation in both the upper and lower basins. Although upper basin precipitation had no statistical trend over the recent historical period, the lower basin had a weak negative trend over this period. Multidecadal-scale cyclical patterns in temperature also are observed in historical climate data in both the upper and lower basins, at least until the early 1970s. Beginning at that time, both the upper and lower basins experienced strong, monotonic positive trends in temperature. Basic principles of hydrology indicate that periods of decreasing precipitation as well as increasing temperature would have a negative effect, that is, reduction in groundwater infiltration and hence, reduced recharge of aquifer systems. Projected climate data from 97 Coupled Model Intercomparison Project phase 5 (CMIP5) ensemble members across the full range of Representative Concentration Pathway (RCPs) from water years 1951 through 2099 were evaluated to understand what current global climate models are projecting about future conditions in the Colorado River Basin, and what this might mean for groundwater systems in the region. Precipitation in the upper basin is projected to increase throughout the rest of the century, rising to 6 percent above the 1951–2015 historical period by mid-century and to 9 percent above the historical period by the end of the century. Temperature in the upper basin also is projected to be above the recent historical median throughout the rest of the century, with steady warming in decadal average temperatures expected until the last quarter of this century. In contrast to projected precipitation in the upper basin, precipitation in the lower basin is projected to be the same as, or slightly less than, the historical period throughout most of the rest of this century. Like projected temperature in the upper basin, temperature in the lower basin also is projected to be above the recent historical median throughout the rest of the century. Comparing median projections for all future decades with median results from all historical decades, future precipitation is expected to be greater than that of the past in the upper basin, though no significant difference is projected for precipitation in the lower basin. Significant increases (p-value<0.05) are expected in temperature in both the upper and lower basins. To estimate the effects of projected precipitation and temperature on groundwater systems in the region, results from the 97 member CMIP5 climate projection ensemble were used as input in a Soil-Water Balance (SWB) groundwater infiltration model for the Colorado River Basin. SWB simulation results indicate that the upper Colorado River Basin is expected to experience decades of above-historical-average groundwater infiltration through the end of the century. For the lower Colorado River Basin, simulated groundwater infiltration is projected to be consistently less than the recent (1951–2015) historical period for most of the remaining century. A comparison of the distribution of all median simulated groundwater infiltration results between recent historical and future periods indicates projected groundwater infiltration in the upper basin is significantly (p-value<0.05) greater over the combined 2020–2099 future period than the recent (1951–2015) historical period. Moreover, in 41 of 71 (58 percent) possible future decades in this century, groundwater infiltration is projected to be greater than the 75th percentile of historical simulated groundwater infiltration. Projected groundwater infiltration in the lower Colorado River Basin across all future decades is significantly less than in the historical period. Of the 71 future decades in the century, projected groundwater infiltration in the lower basin is expected to be less than the 25th percentile of historical infiltration in 55 (77 percent) of the 10-year periods. Important differences in projected precipitation between the upper (increasing precipitation) and lower (decreasing precipitation) basins largely drive the different responses of simulated groundwater infiltration in the upper (increasing infiltration) and lower (decreasing infiltration) basins. It will be useful to revisit projections in groundwater infiltration in the Colorado River Basin when more up-to-date projections of precipitation become available from the next Coupled Model Intercomparison Project phases or by using climate input developments through Regional Climate Modeling efforts and stochastic weather generators.

Arizona, California, Colorado, Nevada, New Mexico,↗

USGS2022

Hydrology of the Yucaipa groundwater subbasin: Characterization and integrated numerical model, San Bernardino and Riverside Counties, California

Executive Summary Water management in the Santa Ana River watershed in San Bernardino and Riverside Counties in southern California is a complex task with various water purveyors navigating geographic, geologic, hydrologic, and political challenges to provide a reliable water supply to stakeholders. As the population has increased throughout southern California, so has the demand for water. The Yucaipa groundwater subbasin (hereafter referred to as “Yucaipa subbasin”), one of nine groundwater subbasins in what the California Department of Water Resources (DWR) refers to as the Upper Santa Ana Valley groundwater basin (California Department of Water Resources, 2016; the DWR naming convention is used within this report), is no exception; steady population growth since the 1940s and changes in water use has forced local water purveyors to regularly adapt their water infrastructure to meet demand. Groundwater has historically been the dominant source of water in the Yucaipa subbasin although recently, imported water via the California State Water Project has augmented the total water supply. Despite the influx of imported water, overall demand for groundwater continues to rise, and there is concern by local water managers that groundwater levels may adversely impact water supply and (or) decline to a point where it will be uneconomical to produce water, severely limiting the ability of local agencies to meet water-supply demand. To better understand the hydrogeology and water resources in the Yucaipa subbasin, the U.S. Geological Survey (USGS) and the San Bernardino Valley Municipal Water District initiated a cooperative study to understand the hydrogeologic system of the Yucaipa subbasin and in the encompassing Yucaipa Valley watershed (YVW). A three-dimensional hydrogeologic framework model was constructed to quantify the structure and extent of hydrogeologic units. Historical and present-day groundwater conditions were characterized to evaluate the groundwater-flow system. Lastly, the Yucaipa Integrated Hydrological Model (YIHM) was developed to simulate the integrated surface-water and groundwater systems, including natural and anthropogenic (that is, human influenced) recharge and discharge throughout the study area from 1947 to 2014. The Yucaipa subbasin is an inland groundwater basin located about 12 miles (mi) southeast of the City of San Bernardino and about 75 mi east of Los Angeles, California. The subbasin encompasses about 39 square miles (mi 2 ), including the City of Yucaipa. The geographic extent of the Yucaipa subbasin was established by the California Department of Water Resources, who defined the boundaries of the subbasin based on hydrogeologic transitions between crystalline rock and basin-fill sediments, active fault strands, surface-water drainage divides, and a portion of an adjudicated groundwater management boundary. Two groundwater subbasins of the Upper Santa Ana Valley groundwater basin are adjacent to the Yucaipa subbasin, the San Bernardino groundwater subbasin to the west and the San Timoteo groundwater subbasin to the south. The Yucaipa subbasin is encompassed by the YVW, which is in turn comprised of three sub-watersheds that represent surface-water flow across and within the Yucaipa subbasin. Although the Yucaipa subbasin is the specific area of interest for this study, the entire YVW was considered for the purposes of characterizing the hydrogeology of the Yucaipa subbasin and for development of the YIHM. The purposes of this report are to (1) describe the hydrologic and hydrogeologic settings of the Yucaipa subbasin and aquifer system, (2) describe the construction and calibration of the fully coupled groundwater and surface-water flow model for the Yucaipa subbasin and the encompassing YVW, referred to as the YIHM, and (3) present numerical results, including water budgets and hydraulic heads, and the effect of pumping and climate stresses (precipitation and temperature) on water-budget components.

California↗

USGS1996

Geohydrology and water quality of stratified-drift aquifers in the middle Connecticut River basin, west-central New Hampshire

A study was done by the U.S. Geological Survey, in cooperation with the New Hampshire Department of Environmental Services, Water Resources Division, to describe the geohydrology and water quality of stratified-drift aquifers in the Middle Connecticut River Basin, west-central New Hampshire Stratified-drift aquifers discontinuously underlie 123 mi2 (square miles) of the Middle Connecticut River Basin, which has a total drainage area of 987 mi 2. Saturated thicknesses of stratified drift in the study area are locally greater than 500 feet but generally are less than 100 feet. Aquifer transmissivity locally exceeds 4,000 ft2/d (feet squared per day) but is generally less than 1,000 ft2/d. In only 17.2 mi2 of the study area are the aquifers identified as having a transmissivity greater than 1,000 ft2/d. As of 1990, total groundwater withdrawals from stratified drift for municipal supply were about 1.5 Mgal/d (million gallons per day) in the study area. Many of the stratified-drift aquifers underlying the study area are not developed to their fullest potential. The geohydrologic investigation of the stratified-drift aquifers focused on aquifer properties, including aquifer boundaries; recharge, discharge, and direction of ground-water flow; saturated thickness and storage; and transmissivity. Surficial-geologic mapping assisted in the determination of aquifer boundaries. Data from more than 1,000 wells, test borings, and springs were used to prepare maps of water-table altitude, saturated thickness, and transmissivity of stratified drift. More than 11 miles of seismic-refraction profiling at 95 sites was used in the preparation of the water-table-altitude and saturated-thickness maps. Seismic-reflection data collected along 1.6 miles of Mascoma Lake also were used in preparation of the saturated-thickness maps. Four stratified-drift aquifers in the towns of Franconia, Haverhill, and Lisbon were analyzed to estimate the water availability on the basis of analytical ground-water model simulation based on the Theis confined-flow equation adjusted to account for boundary effects commonly associated with stratified-drift aquifers. Conservative estimates of water availability during a 180-day period of no recharge were estimated to be 1.9 Mgal/d for the Meadow Brook aquifer; 1.8 Mgal/d for the Ham Branch Brook aquifer; 1.5 Mgal/d for the Salmon Hole aquifer; and 1.4 Mgal/d for the Haverhill-French Pond aquifer. Water-availability estimates would be higher if periods of recharge were accounted for and if less conservative boundary conditions were used in the model. Results of analysis of water samples from 26 observation wells, 3 municipal water-supply wells, and 1 public-supply spring show that, with the exception of dissolved iron and manganese in some samples, water in the stratified-drift aquifers generally meets the U.S. Environmental Protection Agency's primary and secondary drinking-water standards.

New Hampshire↗

USGS2011

Modeling the potential impact of seasonal and inactive multi-aquifer wells on contaminant movement to public water-supply wells

Wells screened across multiple aquifers can provide pathways for the movement of surprisingly large volumes of groundwater to confined aquifers used for public water supply (PWS). Using a simple numerical model, we examine the impact of several pumping scenarios on leakage from an unconfined aquifer to a confined aquifer and conclude that a single inactive multi-aquifer well can contribute nearly 10% of total PWS well flow over a wide range of pumping rates. This leakage can occur even when the multi-aquifer well is more than a kilometer from the PWS well. The contribution from multi-aquifer wells may be greater under conditions where seasonal pumping (e.g., irrigation) creates large, widespread downward hydraulic gradients between aquifers. Under those conditions, water can continue to leak down a multi-aquifer well from an unconfined aquifer to a confined aquifer even when those multi-aquifer wells are actively pumped. An important implication is that, if an unconfined aquifer is contaminated, multi-aquifer wells can increase the vulnerability of a confined-aquifer PWS well.

Journal of the American Water Resources Associatio↗

USGS2010

Summary of hydrologic testing of the Floridan aquifer system at Hunter Army Airfield, Chatham County, Georgia

A 1,168-foot deep test well was completed at Hunter Army Airfield in the summer of 2009 to investigate the potential of using the Lower Floridan aquifer as a source of water supply to satisfy increased needs as a result of base expansion and increased troop levels. The U.S. Geological Survey conducted hydrologic testing at the test site including flowmeter surveys, packer-slug tests, and aquifer tests of the Upper and Lower Floridan aquifers. Flowmeter surveys were completed at different stages of well construction to determine the depth and yield of water-bearing zones and to identify confining beds that separate the main producing aquifers. During a survey when the borehole was open to both the upper and lower aquifers, five water-bearing zones in the Upper Floridan aquifer supplied 83.5 percent of the total pumpage, and five water-bearing zones in the Lower Floridan aquifer supplied the remaining 16.5 percent. An upward gradient was indicated from the ambient flowmeter survey: 7.6 gallons per minute of groundwater was detected entering the borehole between 750 and 1,069 feet below land surface, then moved upward, and exited the borehole into lower-head zones between 333 and 527 feet below land surface. During a survey of the completed Lower Floridan well, six distinct water-producing zones were identified; one 17-foot-thick zone at 768-785 feet below land surface yielded 47.9 percent of the total pumpage while the remaining five zones yielded between 2 and 15 percent each. The thickness and hydrologic properties of the confining unit separating the Upper and Lower Floridan aquifers were determined from packer tests and flowmeter surveys. This confining unit, which is composed of rocks of Middle Eocene age, is approximately 160 feet thick with horizontal hydraulic conductivities determined from four slug tests to range from 0.2 to 3 feet per day. Results of two separate slug tests within the middle confining unit were both 2 feet per day. Aquifer testing indicated the Upper Floridan aquifer had a transmissivity of 40,000 feet squared per day, and the Lower Floridan aquifer had a transmissivity of 10,000 feet squared per day. An aquifer test conducted on the combined aquifer system, when the test well was open from 333 to 1,112 feet, gave a transmissivity of 50,000 feet squared per day. Additionally, during the 72-hour test of the Lower Floridan aquifer, a drawdown response was observed in the Upper Floridan aquifer wells.

Georgia↗

USGS2015

Potentiometric surface, 2012, and water-level differences, 2005-12, of the Sparta Aquifer in north-central Louisiana

The Sparta aquifer is used in 15 parishes in north-central Louisiana, primarily for public supply and industrial purposes. Of those parishes, eight (Bienville, Claiborne, Jackson, Lincoln, Ouachita, Union, Webster, and Winn) rely on the Sparta aquifer as their principal source of groundwater. In 2010, withdrawals from the Sparta aquifer in Louisiana totaled 63.11 million gallons per day (Mgal/d), a reduction of more than 11 percent from 1995, when the highest rate of withdrawals (71.32 Mgal/d) from the Sparta aquifer were documented. The Sparta aquifer provides water for a variety of purposes which include public supply (34.61 Mgal/d), industrial (25.60 Mgal/d), rural domestic (1.50 Mgal/d), and various agricultural (1.40 Mgal/d). Of the 13 major aquifers or aquifer systems in Louisiana, the Sparta aquifer is currently (2012) the sixth most heavily pumped. The Sparta aquifer is the second most heavily pumped aquifer in Arkansas, which borders Louisiana to the north. In 2005, 170 Mgal/d were withdrawn from the Sparta aquifer in eastern and southern Arkansas; of that total, about 15.55 Mgal/d were withdrawn from the aquifer in Union County, which borders Claiborne and Union Parishes to the north. By 1997, a large cone of depression (a cone-shaped depression in the potentiometric surface caused by and centered on a pumping well or wells) in the Sparta aquifer centered over Union County had merged with the cone of depression at West Monroe. In 2004, the rate of withdrawal from the Sparta aquifer in Union County began to decline and water levels in the aquifer began to rise in nearby areas of Arkansas and Louisiana.

Louisiana↗

USGS2017

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

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

Oregon↗

USGS2019

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

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

Utah↗

USGS2020

Spectral analysis to quantify the response of groundwater levels to precipitation — Northwestern United States

Persistent atmospheric patterns that lead to wet and dry seasons and droughts over periods of months to decades and longer-term climate change over periods of decades to millennia affect groundwater resources. Changes in groundwater storage and the resulting groundwater discharge from most aquifers is relatively slow and steady compared to the variability of daily precipitation. The response of groundwater levels to precipitation can be complex because of a combination of processes that include evapotranspiration, surface runoff, and infiltration of net recharge from precipitation through the vadose zone. Typically, this response is delayed and results in a change in groundwater storage reflected in a time series of groundwater levels. Understanding the relations between variations in precipitation and changes in groundwater storage is essential to water resources planning. The objectives of this study were to (1) characterize the relation between precipitation and responses in groundwater levels at seasonal to decadal scales and (2) to develop methods that are transferable on a continental scale to any groundwater-level record. Spectral analysis was applied to daily precipitation and groundwater levels for eleven monitoring wells in the northwestern Unites States with records ranging in length from 5.9 to 23.9 years. The analysis provided a quantitative characterization for each monitoring well that met both objectives and indicated that maximum and minimum precipitation rates generally occurred in December and August, respectively. Maximum groundwater levels occurred from February to August, and minimum values occurred from January to December. The lag in the annual peak response of groundwater to peak precipitation ranged from 2.2 to 8.8 months, with a median value of 5.3 months. Groundwater responses to wet and dry seasons were evident in the relatively high amplitude frequencies of 10 and 20 cycles per decade (cpdec). A high amplitude frequency at 1 cpdec represents a drought cycle that resulted in larger groundwater level changes than typical seasonal water level fluctuations.

Oregon, Washington↗

USGS1997

Hydrogeology and ground-water flow in the Edwards-Trinity aquifer-system, west-central, Texas

Two finite-element ground-water flow models were developed for the Edwards–Trinity aquifer system, west-central Texas, to gain a better understanding of the flow system; one ground-water flow model was developed at a large scale to simulate the regional system and contiguous, hydraulically connected units, and one model was constructed at a smaller more detailed scale to simulate the most active areas of the system. The study area is divided into four geographic subareas: the Trans-Pecos (9,750 square miles), the Edwards Plateau (23,750 square miles), the Hill Country (5,500 square miles), and the Balcones fault zone (3,000 square miles). The major aquifers within the study area are the Edwards–Trinity aquifer underlying the Trans-Pecos and Edwards Plateau, the Trinity aquifer underlying the Hill Country, and the Edwards aquifer in the Balcones fault zone. Hydraulically connected aquifers include the High Plains aquifer north of the Edwards Plateau, and the Cenozoic Pecos alluvium aquifer adjacent to both the Trans-Pecos and the Edwards Plateau along the Pecos River. Minor contiguous aquifers include the Dockum, Ellenburger– San Saba, Marble Falls, Hickory, and Lipan, which is adjacent to the Colorado River in Tom Green and Concho Counties, Texas. The ground-water flow equations solved by the finite-element method are based on conservation of mass and energy. The equation for ground-water flow assumes laminar flow through a porous media. In places, the Edwards–Trinity aquifer system is a fractured karst system in which ground water flows through caverns and other features of secondary porosity development. The regional and subregional models were constructed to synthesize the known hydrologic boundaries and geologic structures into a heterogeneous continuum model of the karst ground-water flow system, rather than simulate the flow through specific fractures and caverns. A heterogeneous continuum or equivalent porous media approach uses an effective transmissivity and anisotropy for each element of the models. The models are calibrated both on water levels (representing potential energy) and estimates of recharge and discharge (for a realistic mass balance). A two-dimensional one-layer large-scale model (55,600 square miles) was developed for the Edwards–Trinity aquifer system and contiguous, hydraulically connected units, in westcentral Texas. A quasi-three-dimensional, multilayer more detailed scale ground-water flow model (12,300 square miles) was applied to the major aquifers of the Edwards–Trinity aquifer system in the Hill Country and the Balcones fault zone, and in part of the Edwards Plateau. The ground-water flow system in most of the study area within the Trans-Pecos and Edwards Plateau can be approximated with a one-layer regional model under steady-state conditions. Regionally, the Edwards–Trinity aquifer system in the Trans-Pecos and Edwards Plateau has been relatively static. Potentiometric maps from predevelopment and postdevelopment (winter 1974–75) indicate small differences in water levels. In local areas in the Trans-Pecos (in Pecos and Reeves Counties), ground-water withdrawals have exceeded recharge resulting in more than 300 feet of drawdown. Measurable differences between the 1974 and predevelopment potentiometric surfaces have been observed in small areas in the Trans-Pecos and in the northwestern part of the Edwards Plateau. The largest water-level declines in the Trans-Pecos have been observed in Pecos and Reeves Counties, and declines greater than 300 feet have been measured in Reeves County. Comparison of pre- and postdevelopment water budgets for the regional model indicates that the increase in groundwater withdrawals has captured 20 percent of the water that would have naturally discharged to streams, and 30 percent of the natural discharge to springs after ground-water development. Induced recharge from streams to the ground-water system increased by 12 percent in the postdevelopment simulation compared to the predevelopment simulation.

Texas↗

USGS2010

Nitrate Loads and Concentrations in Surface-Water Base Flow and Shallow Groundwater for Selected Basins in the United States, Water Years 1990-2006

Hydrograph separation was used to determine the base-flow component of streamflow for 148 sites sampled as part of the National Water-Quality Assessment program. Sites in the Southwest and the Northwest tend to have base-flow index values greater than 0.5. Sites in the Midwest and the eastern portion of the Southern Plains generally have values less than 0.5. Base-flow index values for sites in the Southeast and Northeast are mixed with values less than and greater than 0.5. Hypothesized flow paths based on relative scaling of soil and bedrock permeability explain some of the differences found in base-flow index. Sites in areas with impermeable soils and bedrock (areas where overland flow may be the primary hydrologic flow path) tend to have lower base-flow index values than sites in areas with either permeable bedrock or permeable soils (areas where deep groundwater flow paths or shallow groundwater flow paths may occur). The percentage of nitrate load contributed by base flow was determined using total flow and base flow nitrate load models. These regression-based models were calibrated using available nitrate samples and total streamflow or base-flow nitrate samples and the base-flow component of total streamflow. Many streams in the country have a large proportion of nitrate load contributed by base flow: 40 percent of sites have more than 50 percent of the total nitrate load contributed by base flow. Sites in the Midwest and eastern portion of the Southern Plains generally have less than 50 percent of the total nitrate load contributed by base flow. Sites in the Northern Plains and Northwest have nitrate load ratios that generally are greater than 50 percent. Nitrate load ratios for sites in the Southeast and Northeast are mixed with values less than and greater than 50 percent. Significantly lower contributions of nitrate from base flow were found at sites in areas with impermeable soils and impermeable bedrock. These areas could be most responsive to nutrient management practices designed to reduce nutrient transport to streams by runoff. Conversely, sites with potential for shallow or deep groundwater contribution (some combination of permeable soils or permeable bedrock) had significantly greater contributions of nitrate from base flow. Effective nutrient management strategies would consider groundwater nitrate contributions in these areas. Mean annual base-flow nitrate concentrations were compared to shallow-groundwater nitrate concentrations for 27 sites. Concentrations in groundwater tended to be greater than base-flow concentrations for this group of sites. Sites where groundwater concentrations were much greater than base-flow concentrations were found in areas of high infiltration and oxic groundwater conditions. The lack of correspondingly high concentrations in the base flow of the paired surface-water sites may have multiple causes. In some settings, there has not been sufficient time for enough high-nitrate shallow groundwater to migrate to the nearby stream. In these cases, the stream nitrate concentrations lag behind those in the shallow groundwater, and concentrations may increase in the future as more high-nitrate groundwater reaches the stream. Alternatively, some of these sites may have processes that rapidly remove nitrate as water moves from the aquifer into the stream channel. Partitioning streamflow and nitrate load between the quick-flow and base-flow portions of the hydrograph coupled with relative scales of soil permeability can infer the importance of surface water compared to groundwater nitrate sources. Study of the relation of nitrate concentrations to base-flow index and the comparison of groundwater nitrate concentrations to stream nitrate concentrations during times when base-flow index is high can provide evidence of potential nitrate transport mechanisms. Accounting for the surface-water and groundwater contributions of nitrate is crucial to effective management and remediat

Scientific Investigations Report↗

USGS2015

Simulation of the regional groundwater-flow system of the Menominee Indian Reservation, Wisconsin

A regional, two-dimensional, steady-state groundwater-flow model was developed to simulate the groundwater-flow system and groundwater/surface-water interactions within the Menominee Indian Reservation. The model was developed by the U.S. Geological Survey (USGS), in cooperation with the Menominee Indian Tribe of Wisconsin, to contribute to the fundamental understanding of the region’s hydrogeology. The objectives of the regional model were to improve understanding of the groundwater-flow system, including groundwater/surface-water interactions, and to develop a tool suitable for evaluating the effects of potential regional water-management programs. The computer code GFLOW was used because of the ease with which the model can simulate groundwater/surface-water interactions, provide a framework for simulating regional groundwater-flow systems, and be refined in a stepwise fashion to incorporate new data and simulate groundwater-flow patterns at multiple scales. Simulations made with the regional model reproduce groundwater levels and stream base flows representative of recent conditions (1970–2013) and illustrate groundwater-flow patterns with maps of (1) the simulated water table and groundwater-flow directions, (2) probabilistic areas contributing recharge to high-capacity pumped wells, and (3) estimation of the extent of infiltrated wastewater from treatment lagoons. The groundwater-flow model described in this report simulates the major hydrogeologic features of the modeled area, including surficial unconsolidated aquifers, groundwater/surface-water interactions, and groundwater withdrawals from existing high-capacity production wells. Areas contributing recharge to pumped high-capacity wells on the Menominee Indian Reservation were delineated by tracking simulated water particles from the water table to wells in combination with Monte Carlo techniques, and maps of the probability of capture for each well nest were produced. Groundwater-agebased areas contributing recharge to wells were simulated by using the calibrated set of parameters and porosity values adjusted to account for bias in simulated saturated thickness. Simulations were performed for current (2013) pumping rates. The simulations show a range in sensitivity of the simulated areas contributing recharge to wells given the parameters evaluated through the Monte Carlo analysis. The areas contributing recharge to supply wells for the villages of Zoar and Neopit are long and narrow, with a sharp gradation from high to low probability of capture. The areas contributing recharge to supply wells for Middle Village and the village of Keshena exhibit a sharp gradation from high to low probability over a relatively small area between the well and a local groundwater mound. The highest probability areas contributing recharge to the supply wells for the Villages of Onekewat and Redwing are in the immediate vicinity of the wells. These wells also have an extensive area with low probability for capturing water that is likely due to a locally low hydraulic gradient and the large degree of uncertainty associated with the lakebed resistance parameters that control interaction between groundwater and local lakes. Additional field investigations and associated local model refinements would facilitate further reductions in uncertainty associated with simulated areas contributing recharge to the wells. The likely extent of the Neopit wastewater plume was simulated by using the groundwater-flow model and Monte Carlo techniques to evaluate the sensitivity of predictive simulations to a range of model parameter values. Wastewater infiltrated from the currently operating lagoons flows predominantly south toward Tourtillotte Creek. Some of the infiltrated wastewater is simulated as having a low probability of flowing beneath Tourtillotte Creek to the nearby West Branch Wolf River. Results for the probable extent of the wastewater plume are considered to be qualitative because the method only considers advective flow and does not account for processes affecting contaminant transport in porous media. Therefore, results for the probable extent of the wastewater plume are sensitive to the number of particles used to represent flow from the lagoon and the resolution of a synthetic grid used for the analysis. Nonetheless, it is expected that the qualitative results may be of use for identifying potential downgradient areas of concern that can then be evaluated using the quantitative “area contributing recharge to wells” method or traditional contaminant-transport simulations.

Wisconsin↗

USGS2001

Effects of ground-water withdrawals on flow in the Sauk River Valley Aquifer and on streamflow in the Cold Spring area, Minnesota

The U. S. Geological Survey, in cooperation with the city of Cold Spring, Minnesota conducted a study during 1998-99 to: (1) determine the contributing areas of groundwater flow to high-capacity wells, (2) delineate the 10-, 20-, and 30-year time-of-travel zones to high-capacity wells, and (3) determine changes in streamflow in the Sauk River due to ground-water withdrawals. Surficial aquifers underlie a portion of the uplands in the north-central part of the study area and the Sauk River Valley. The upland surficial aquifer is unconfined, with saturated thicknesses ranging from 6 to 37 feet. The Sauk River Valley aquifer consists of hydraulically connected sand and gravel units: (1) the generally continuous surficial unconfined unit underlying the Sauk River Valley and (2) buried confined units that are overlain by till or lake clays in the upland areas. Maximum saturated thickness of the Sauk River Valley aquifer is about 50 feet. Thicknesses of the buried units are generally 10 to 20 feet. Ground-water flow in the upland surficial aquifer is (1) to the south and southeast toward the margins of the aquifer and (2) toward an unnamed creek valley in the northwestern part of the study area. Ground-water flow in the Sauk River Valley aquifer is from upland areas toward the Sauk River and the connected chain of lakes. Based on hydrograph analysis, ground-water recharge rates to the aquifers during 1999 ranged from 5.3 to 8.6 inches per year, with an average value of 7.0 inches per year. Streamflow measurements indicated net gains of 44.8 and 25.8 cubic feet per second for the Sauk River from Cold Spring to Rockville during October 1998 and August 1999, respectively. In general, the Sauk River probably is a gaining stream in all reaches, except for the reach adjacent to the Gold'n Plump Poultry Processing Plant well field. Ground-water withdrawals from the well field induce infiltration of water from the Sauk River to the underlying aquifer. The measured gains in streamflow for Brewery Creek were 2.75 and 2.25 cubic feet per second during October 1998 and August 1999, respectively. A numerical ground-water-flow model was constructed and calibrated for steady-state conditions. Based on the calibrated model, areal recharge accounts for 51.5 percent of the sources of water to the aquifers in the Cold Spring area and inflow through constant-head boundaries contributes 45.8 percent. The largest discharges from the aquifers are leakage from the Sauk River Valley aquifer to the Sauk River and Brewery Creek (53.7 percent) and outflow through constant-head boundaries (33.1 percent). The simulated contributing areas for selected watersupply wells in the Cold Spring area generally extend to and possibly beyond the model boundaries to the north and to the southeast. The contributing areas for the Gold'n Plump Poultry Processing Plant supply wells extend: (1) to the Sauk River, (2) to the north to and possibly beyond to the northern model boundary, and (3) to the southeast to and possibly beyond the southeastern model boundary. The primary effects of projected increased ground-water withdrawals of 0.23 cubic feet per second (7.5 percent increase) were to: (1) decrease outflow from the Sauk River Valley aquifer through constant-head boundaries and (2) decrease leakage from the valley unit of the Sauk River Valley aquifer to the streams. No appreciable differences were discernible between the simulated steady-state contributing areas to wells with 1998 pumpage and those with the projected pumpage.

Minnesota↗

USGS2020

A comparison of groundwater sampling technologies, including passive diffusion sampling, for radionuclide contamination

Using traditional high-flow purge methods for long-term water quality monitoring of deep groundwater wells can be expensive, affect contaminant migration, and produce excessive volumes of discharge water that can be difficult to manage. The use of low-flow pumping methods and depth discrete bailers (DDBs) can reduce the cost of sampling deep groundwater wells. In general, using different pumping methods to obtain reproducible and representative groundwater can be challenging. Passive diffusion samplers (PDSs) have successfully been used in long-term monitoring for volatile organic compounds, major ions, and trace-elements, but application has been limited for stable and radioactive isotopes. Beginning in 2018, the United States Geological Survey (USGS) conducted three sampling events to test the ability of PDSs and DDBs to obtain reproducible and representative groundwater samples. The sampled well is completed in a regional, permeable carbonate aquifer and is one in a set of wells that have historically been used for tritium tracer testing. All samples were obtained at 180 m (590 ft) below land surface (bls) in a 13.97 cm (5.5 inch) uncased well that has a total depth of 202 m (662 ft) bls. The first sampling event deployed a regenerated cellulose dialysis membrane (RCDM) PDS for 14 days with deionized water as the blank. The second sampling event deployed a RCDM for 27 days also with deionized water as the blank. The third sampling event deployed a Dual Membrane (DM) PDS for 65 days using a blank of tritium-dead carbonate water. The DM PDS was used to assess the effect of longer-term deployment on tritium concentrations and address whether or not the PDSs reached equilibrium with ambient groundwater. The day after each of the three passive samplers were retrieved a DDB was used to obtain discrete non-integrated groundwater samples. For each DDB sampling day, the bailer was lowered into the well 10 consecutive times to determine if the water chemistry changed from the first to the last bailed sample. Quality assurance samples including blanks and duplicates were obtained during all three sampling events. All blank waters had tritium concentrations less than 21±33 pCi/L. Major ion (e.g. calcium, chloride, sodium, and sulfate) results were compared between all samples obtained with RCDM and DDB. Major ion concentrations showed a coefficient of variation of less than 6% between all RCDM and DDB samples; however, the coefficient of variation between the different deployment times and the two different methods for trace-element concentrations was much larger, particularly for manganese (82%), lead (41%), and zinc (33%). Stable isotope values were compared between the RCDM and DDB samples. The DDB sample results all fell within analytical uncertainty and were considered representative of the formation groundwater. The stable isotope values from the RCDM samples indicated that a longer deployment time was necessary to gain equilibrium and to obtain representative groundwater samples. Tritium results from groundwater samples obtained from the RCDM, DM, and DDB indicate that groundwaters obtained with PDSs produced tritium concentrations 2 to 3 times higher (between 1,426 and 3,060±87 pCi/L) than groundwaters obtained with a DDB (479 to 1,219±51 pCi/L). The longer the passive diffusion samplers were deployed, the higher the tritium concentration, suggesting that equilibrium with tritium was not reached within a 27-day deployment. DDB samples showed tritium results declining from the first to the last bailed sample for all three sampling events. This research suggests that tritium results from groundwater samples obtained from PDSs are more reproducible than samples obtained from DDBs. Also, PDSs likely do not accumulate isotopes of water but rather equilibrate with the ambient groundwater. On the other hand, both PDSs and DDBs were able to provide representative groundwater samples for major ions and have the potential to produce.

Conference Paper↗

USGS2024

Predicting redox conditions in groundwater at a national scale using random forest classification

Redox conditions in groundwater may markedly affect the fate and transport of nutrients, volatile organic compounds, and trace metals, with significant implications for human health. While many local assessments of redox conditions have been made, the spatial variability of redox reaction rates makes the determination of redox conditions at regional or national scales problematic. In this study, redox conditions in groundwater were predicted for the contiguous United States using random forest classification by relating measured water quality data from over 30,000 wells to natural and anthropogenic factors. The model correctly predicted the oxic/suboxic classification for 78 and 79% of the samples in the out-of-bag and hold-out data sets, respectively. Variables describing geology, hydrology, soil properties, and hydrologic position were among the most important factors affecting the likelihood of oxic conditions in groundwater. Important model variables tended to relate to aquifer recharge, groundwater travel time, or prevalence of electron donors, which are key drivers of redox conditions in groundwater. Partial dependence plots suggested that the likelihood of oxic conditions in groundwater decreased sharply as streams were approached and gradually as the depth below the water table increased. The probability of oxic groundwater increased as base flow index values increased, likely due to the prevalence of well-drained soils and geologic materials in high base flow index areas. The likelihood of oxic conditions increased as topographic wetness index (TWI) values decreased. High topographic wetness index values occur in areas with a propensity for standing water and overland flow, conditions that limit the delivery of dissolved oxygen to groundwater by recharge; higher TWI values also tend to occur in discharge areas, which may contain groundwater with long travel times. A second model was developed to predict the probability of elevated manganese (Mn) concentrations in groundwater (i.e., ≥50 μg/L). The Mn model relied on many of the same variables as the oxic/suboxic model and may be used to identify areas where Mn-reducing conditions occur and where there is an increased risk to domestic water supplies due to high Mn concentrations. Model predictions of redox conditions in groundwater produced in this study may help identify regions of the country with elevated groundwater vulnerability and stream vulnerability to groundwater-derived contaminants.

contiguous United States↗

USGS1993

Hydrogeology and water quality of the Galena-Platteville aquifer at the Parson's Casket Hardware Superfund site, Belvidere, Illinois, 1991

The U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, investigated the hydrogeology of the Galena-Platteville aquifer and its relation to contaminant migration at the Parson's Casket Hardware Superfund site in Belvidere, Ill. This report presents the results of the second phase of the investigation, which lasted from March through October 1991. The uppermost bedrock units beneath the study site are the Galena and Platteville Groups1; these bedrock units immediately underlie a glacial drift aquifer. The Galena and Platteville Groups, which consist predominantly of dolomite, compose the Galena-Platteville aquifer, and extend from about 40 to 320 feet below land surface. The unconfined Galena-Platteville aquifer is partitioned into five hydrogeologic units. The uppermost unit, the weathered surface of the bedrock, has a horizontal hydraulic conductivity that ranges from about 1 to 200 feet per day. The four underlying units have hydraulic conductivities that range from about 0.01 to 1 foot per day. Vertical hydraulic gradients in the aquifer are typically downward. Horizontal groundwater flow generally is southward to southeastward from the site toward the Kishwaukee River. Three notable bedding-plane solution fissures and three fractures that crosscut the bedding planes are identified within the dolomite bedrock. The inclined fractures are assumed to function as conduits that connect high conductivity horizontal fissures, thus allowing more rapid vertical movement of ground water and contaminants than would be expected in the generally low conductivity dolomite matrix. A multiple-well, constant-discharge aquifer test confirms the heterogeneity and anisotropy of the dolomite aquifer. The hydraulic characteristics of the uppermost part of the bedrock aquifer are somewhat different than the characteristics of the deeper part(s) of the aquifer. This is because the principal conduits for water movement are in the deeper part(s) of the aquifer. Ground-water flow in the Galena-Platteville aquifer beneath the study site is affected by pumping of Belvidere Municipal Wells No. 4 and No. 6. Water levels in wells screened near the vertical midpoint of the aquifer are affected by the municipal-well pumping substantially more than are water levels in wells screened near the top of the aquifer.

Open-File Report↗

USGS2021

Groundwater dynamics at Kīlauea Volcano and vicinity, Hawaiʻi

Kīlauea Volcano, on the Island of Hawaiʻi, is surrounded and permeated by active groundwater systems that interact dynamically with the volcanic system. A generalized conceptual model of Hawaiian hydrogeology includes high-level dike-impounded groundwater, very permeable perched and basal aquifers, and a transition (mixing) zone between freshwater and saltwater. Most high-level groundwater is associated with the low-permeability intrusive complexes that underlie volcanic rift zones and calderas and also act to compartmentalize the groundwater system. Hydrogeologic studies of Kīlauea in recent decades, accompanied by deep research drilling, have shown that high-level groundwater is more widespread than once understood, that permeability decreases dramatically at depth, particularly in rift zones, and that freshwater can occur at depths of as much as several kilometers below the local water table. Copious groundwater recharge causes near-surface conductive heat flow to be near zero over much of Kīlauea. Approximately 95 percent of groundwater discharge occurs offshore, accompanied by approximately 99 percent of the approximately 6,000 megawatts of heat supplied by magmatic intrusion. Here, we summarize current understanding of the groundwater system of Kīlauea Volcano and describe transient changes during the decade or more preceding the 2018 eruption sequence. The changes in groundwater chemistry and thermal structure beneath Kīlauea summit hold implications for volcanic-volatile transport and the potential for explosive volcanism. Between 2008 and 2018, the magma conduit beneath the lava lake likely created an adjacent zone of very hot rock that significantly delayed liquid groundwater inflow to the draining magma conduit. Sulfate concentrations in groundwater beneath Kīlauea summit, sampled at the National Science Foundation-funded drill hole 1.5 kilometers south-southwest of the lava lake, declined substantially between 2010 and present. This decline likely reflects, at least in part, the decreased effectiveness of volatile condensation and solution into groundwater (scrubbing). The vent opening in 2008 presumably focused volatile flux into the vicinity of the vent, and progressive drying of the surroundings further restricted interaction with the groundwater system. The decrease in sulfate concentrations in the drill hole between 2010 and 2018 likely reflects decreased effectiveness of scrubbing.

Hawaii↗