USGS Science⌕ Search

SEARCH · USGS Science

Results for “groundwater aquifer”

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

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

At least 253 records · Page 14Linked to original sources

Hydrologic assessment of three drainage basins in the Pinelands of southern New Jersey, 2004-06

The New Jersey Pinelands is an ecologically diverse area in the southern New Jersey Coastal Plain, most of which overlies the Kirkwood-Cohansey aquifer system. The demand for groundwater from this aquifer system is increasing as local development increases. Because any increase in groundwater withdrawals has the potential to affect streamflows and wetland water levels, and ultimately threaten the ecological health and diversity of the Pinelands ecosystem, the U.S. Geological Survey, in cooperation with the New Jersey Pinelands Commission, began a multi-phase hydrologic investigation in 2004 to characterize the hydrologic system supporting the aquatic and wetland communities of the New Jersey Pinelands area (Pinelands). The current investigation of the hydrology of three representative drainage basins in the Pinelands (Albertson Brook, McDonalds Branch, and Morses Mill Stream basins) included a compilation of existing data; collection of water-level and streamflow data; mapping of the water-table altitude and depth to the water table; and analyses of water-level and streamflow variability, subsurface gradients and flow patterns, and water budgets. During 2004-06, a hydrologic database of existing and new data from wells and stream sites was compiled. Methods of data collection and analysis were defined, and data networks consisting of 471 wells and 106 surface-water sites were established. Hydrographs from 26 water-level-monitoring wells and four streamflow-gaging stations were analyzed to show the response of water levels and streamflow to precipitation and recharge with respect to the locations of these wells and streams within each basin. Water-level hydrographs show varying hydraulic gradients and flow potentials, and indicate that responses to recharge events vary with well depth and proximity to recharge and discharge areas. Results of the investigation provide a detailed characterization of hydrologic conditions, processes, and relations among the components of the hydrologic cycle in the Pinelands. In the Pinelands, recharge replenishes the aquifer system and contributes to groundwater flow, most of which moves to wetlands and surface water where natural discharge occurs. Some groundwater flow is intercepted by supply wells. Recharge rates generally are highest during the non-growing season and are inversely related to evapotranspiration. Analysis of subsurface hydraulic gradients, water-table fluctuations, and streamflow variability indicates a strong linkage between groundwater and wetlands, lakes and streams. Gradient analysis indicates that most wetlands are in groundwater discharge areas, but some wetlands are in groundwater recharge areas. The depth to the water table ranges from zero at surface-water features up to about 10 meters in topographically high areas. Depth to water fluctuates seasonally, and the magnitude of these fluctuations generally increases with distance from surface water. Variations in the permeability of the soils and sediments of the aquifer system strongly affect patterns of water movement through the subsurface and the interaction of groundwater with wetlands, lakes and streams. Mean annual streamflow during 2004-06 ranged from 83 to 106 percent of the long-term mean annual discharge, indicating that the data-collection period can be considered representative of average conditions. Measurements of groundwater levels, stream stage, and stream discharge and locations of start-of-flow are illustrated in basin-wide maps of water-table altitude, depth to the water table, and stream base flow during the period. Water-level data collected along 15 hydrologic transects that span the range of environments from uplands through wetlands to surface water were used to determine hydraulic gradients, potential flow directions, and areas of recharge and discharge. These data provide information about the localized interactions of groundwater with wetlands and surface water. Wetlands were categorized with r

New Jersey↗

Factors associated with sources, transport, and fate of volatile organic compounds and their mixtures in aquifers of the United States

Factors associated with sources, transport, and fate of volatile organic compounds (VOCs) in groundwater from aquifers throughout the United States were evaluated using statistical methods. Samples were collected from 1631 wells throughout the conterminous United States between 1996 and 2002 as part of the National Water-Quality Assessment (NAWQA) Program of the U.S. Geological Survey. Water samples from wells completed in aquifers used to supply drinking water were analyzed for more than 50 VOCs. Wells were primarily rural domestic water supplies (1184), followed by public water supplies (216); the remaining wells (231) supplied a variety of uses. The median well depth was 50 meters. Age-date information shows that about 60% of the samples had a fraction of water recharged after 1953. Chloroform, toluene, 1,2,4-trimethylbenzene, and perchloroethene were some of the frequently detected VOCs. Concentrations generally were less than 1 ??g/L. Source factors include, in order of importance, general land-use activity, septic/sewer density, and sites where large concentrations of VOCs are potentially released, such as leaking underground storage tanks. About 10% of all samples had VOC mixtures that were associated with concentrated sources; 20% were associated with dispersed sources. Important transport factors included well/screen depth, precipitation/groundwater recharge, air temperature, and various soil characteristics. Dissolved oxygen was strongly associated with VOCs and represents the fate of many VOCs in groundwater. Well type (domestic or public water supply) was also an important explanatory factor. Results of multiple analyses show the importance of (1) accounting for both dispersed and concentrated sources of VOCs, (2) measuring dissolved oxygen when sampling wells to help explain the fate of VOCs, and (3) limiting the type of wells sampled in monitoring networks to avoid unnecessary variance in the data, or controlling for this variance during data analysis.

Environmental Science & Technology↗

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

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

Virginia↗

The Maryland Coastal Plain Aquifer Information System: A GIS-based tool for assessing groundwater resources

Groundwater is the source of drinking water for ∼1.4 million people in the Coastal Plain Province of Maryland (USA). In addition, groundwater is essential for commercial, industrial, and agricultural uses. Approximately 0.757 × 109 L d –1 (200 million gallons/d) were withdrawn in 2010. As a result of decades of withdrawals from the coastal plain confined aquifers, groundwater levels have declined by as much as 70 m (230 ft) from estimated prepumping levels. Other issues posing challenges to long-term groundwater sustainability include degraded water quality from both man-made and natural sources, reduced stream base flow, land subsidence, and changing recharge patterns (drought) caused by climate change. In Maryland, groundwater supply is managed primarily by the Maryland Department of the Environment, which seeks to balance reasonable use of the resource with long-term sustainability. The chief goal of groundwater management in Maryland is to ensure safe and adequate supplies for all current and future users through the implementation of appropriate usage, planning, and conservation policies. To assist in that effort, the geographic information system (GIS)–based Maryland Coastal Plain Aquifer Information System was developed as a tool to help water managers access and visualize groundwater data for use in the evaluation of groundwater allocation and use permits. The system, contained within an ESRI ArcMap desktop environment, includes both interpreted and basic data for 16 aquifers and 14 confining units. Data map layers include aquifer and ­confining unit layer surfaces, aquifer extents, borehole information, hydraulic properties, time-series groundwater-level data, well records, and geophysical and lithologic logs. The aquifer and confining unit layer surfaces were generated specifically for the GIS system. The system also contains select groundwater-quality data and map layers that quantify groundwater and surface-water withdrawals. The aquifer information system can serve as a pre- and postprocessing environment for groundwater-flow models for use in water-supply planning, development, and management. The system also can be expanded to include features that evaluate constraints to groundwater development, such as insufficient available drawdown, degraded groundwater quality, insufficient aquifer yields, and well-field interference. Ultimately, the aquifer information system is intended to function as an interactive Web-based utility that provides a broad array of information related to groundwater resources in Maryland’s coastal plain to a wide-ranging audience, including well drillers, consultants, academia, and the general public.

Special Paper of the Geological Society of America↗

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

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

Scientific Investigations Report↗

Methanogenic pathways of coal-bed gas in the Powder River Basin, United States: The geologic factor

Coal-bed gas of the Tertiary Fort Union and Wasatch Formations in the Powder River Basin in Wyoming and Montana, U.S. was interpreted as microbial in origin by previous studies based on limited data on the gas and water composition and isotopes associated with the coal beds. To fully evaluate the microbial origin of the gas and mechanisms of methane generation, additional data for 165 gas and water samples from 7 different coal-bed methane-bearing coal-bed reservoirs were collected basinwide and correlated to the coal geology and stratigraphy. The C1/(C2 + C3) ratio and vitrinite reflectance of coal and organic shale permitted differentiation between microbial gas and transitional thermogenic gas in the central part of the basin. Analyses of methane ??13C and ??D, carbon dioxide ??13C, and water ??D values indicate gas was generated primarily from microbial CO2 reduction, but with significant gas generated by microbial methyl-type fermentation (aceticlastic) in some areas of the basin. Microbial CO2 reduction occurs basinwide, but is generally dominant in Paleocene Fort Union Formation coals in the central part of the basin, whereas microbial methyl-type fermentation is common along the northwest and east margins. Isotopically light methane ??13C is distributed along the basin margins where ??D is also depleted, indicating that both CO2-reduction and methyl-type fermentation pathways played major roles in gas generation, but gas from the latter pathway overprinted gas from the former pathway. More specifically, along the northwest basin margin gas generation by methyl-type fermentation may have been stimulated by late-stage infiltration of groundwater recharge from clinker areas, which flowed through highly fractured and faulted coal aquifers. Also, groundwater recharge controlled a change in gas composition in the shallow Eocene Wasatch Formation with the increase of nitrogen and decrease of methane composition of the coal-bed gas. Other geologic factors, such as burial, thermal and maturation history, lateral and vertical continuity, and coalification of the coal beds, also played a significant role in controlling methanogenic pathways and provided new perspectives on gas evolution and emplacement. The early-stage gas produced by CO2 reduction has mixed with transitional thermogenic gas in the deeper, central parts of the Powder River Basin to form 'old' gas, whereas along the basin margins the overprint of gas from methyl-type fermentation represents 'new' gas. Thus, a clear understanding of these geologic factors is necessary to relate the microbiological, biogeochemical, and hydrological processes involved in the generation of coal-bed gas.

International Journal of Coal Geology↗

Drinking water quality in the glacial aquifer system, northern USA

Groundwater supplies 50% of drinking water worldwide, but compromised water quality from anthropogenic and geogenic contaminants can limit usage of groundwater as a drinking water source. Groundwater quality in the glacial aquifer system, USA (GLAC), is presented in the context of a hydrogeologic framework that divides the study area into 17 hydrogeologic terranes. Results are reported at aquifer-system scale and regional (terrane) scale. This paper presents a quantitative assessment of groundwater quality in the GLAC using data from numerous sources for samples collected 2005–2013, compared to health-based and aesthetic (non-health) benchmarks, and evaluated with areal and population metrics. Concentrations above a benchmark are considered high. Trace elements are widespread across the study area, with an estimated 5.7 million people relying on groundwater with high concentrations of one or more trace elements; manganese and arsenic are most often at high concentration. Nitrate is found at high concentration in 4.0% of the study area, serving about 740 thousand people. Organic compounds including pesticides and volatile organic compounds are high in 2.0% of the assessed study area, with about 870 thousand people relying on groundwater with high concentrations of an organic compound. High arsenic and manganese concentrations occur primarily in the terranes with thick, stratigraphically complex, fine-grained glacial sediment, coincident with groundwater under reducing conditions (indicated by iron concentrations >100 μg/L); high nitrate is uncommon in those same terranes. When nitrate is high in thick, fine-grained, complex terranes, though, it is much more commonly associated with groundwater under more oxidizing conditions. Common geogenic trace elements occur at high concentration due to characteristic geologic and geochemical conditions. Conversely, anthropogenic nitrate and organic compounds are introduced at or near the land surface. High concentrations of nitrate or organic compounds are generally limited to areas in proximity where people live and use the chemicals.

Science of Total Environment↗

Irrigated agriculture and future climate change effects on groundwater recharge, northern High Plains aquifer, USA

Understanding the controls of agriculture and climate change on recharge rates is critically important to develop appropriate sustainable management plans for groundwater resources and coupled irrigated agricultural systems. In this study, several physical (total potential ( ψ T ) time series) and chemical tracer and dating ( 3 H, Cl − , Br − , CFCs, SF 6 , and 3 H/ 3 He) methods were used to quantify diffuse recharge rates beneath two rangeland sites and irrigation recharge rates beneath two irrigated corn sites along an east-west (wet-dry) transect of the northern High Plains aquifer, Platte River Basin, central Nebraska. The field-based recharge estimates and historical climate were used to calibrate site-specific Hydrus-1D models, and irrigation requirements were estimated using the Crops Simulation Model (CROPSIM). Future model simulations were driven by an ensemble of 16 global climate models and two global warming scenarios to project a 2050 climate relative to the historical baseline 1990 climate, and simulate changes in precipitation, irrigation, evapotranspiration, and diffuse and irrigation recharge rates. Although results indicate statistical differences between the historical variables at the eastern and western sites and rangeland and irrigated sites, the low warming scenario (+1.0 °C) simulations indicate no statistical differences between 2050 and 1990. However, the high warming scenarios (+2.4 °C) indicate a 25% and 15% increase in median annual evapotranspiration and irrigation demand, and decreases in future diffuse recharge by 53% and 98% and irrigation recharge by 47% and 29% at the eastern and western sites, respectively. These results indicate an important threshold between the low and high warming scenarios that if exceeded could trigger a significant bidirectional shift in 2050 hydroclimatology and recharge gradients. The bidirectional shift is that future northern High Plains temperatures will resemble present central High Plains temperatures and future recharge rates in the east will resemble present recharge rates in the western part of the northern High Plains aquifer. The reductions in recharge rates could accelerate declining water levels if irrigation demand and other management strategies are not implemented. Findings here have important implications for future management of irrigation practices and to slow groundwater depletion in this important agricultural region.

Northern High Plains Aquifer↗

Construction and calibration of a groundwater-flow model to assess groundwater availability in the uppermost principal aquifer systems of the Williston Basin, United States and Canada

The U.S. Geological Survey developed a groundwater-flow model for the uppermost principal aquifer systems in the Williston Basin in parts of Montana, North Dakota, and South Dakota in the United States and parts of Manitoba and Saskatchewan in Canada as part of a detailed assessment of the groundwater availability in the area. The assessment was done because of the potential for increased demands and stresses on groundwater associated with large-scale energy development in the area. As part of this assessment, a three-dimensional groundwater-flow model was developed as a tool that can be used to simulate how the groundwater-flow system responds to changes in hydrologic stresses at a regional scale. The three-dimensional groundwater-flow model was developed using the U.S. Geological Survey’s numerical finite-difference groundwater model with the Newton-Rhapson solver, MODFLOW–NWT, to represent the glacial, lower Tertiary, and Upper Cretaceous aquifer systems for steady-state (mean) hydrological conditions for 1981‒2005 and for transient (temporally varying) conditions using a combination of a steady-state period for pre-1960 and transient periods for 1961‒2005. The numerical model framework was constructed based on existing and interpreted hydrogeologic and geospatial data and consisted of eight layers. Two layers were used to represent the glacial aquifer system in the model; layer 1 represented the upper one-half and layer 2 represented the lower one-half of the glacial aquifer system. Three layers were used to represent the lower Tertiary aquifer system in the model; layer 3 represented the upper Fort Union aquifer, layer 4 represented the middle Fort Union hydrogeologic unit, and layer 5 represented the lower Fort Union aquifer. Three layers were used to represent the Upper Cretaceous aquifer system in the model; layer 6 represented the upper Hell Creek hydrogeologic unit, layer 7 represented the lower Hell Creek aquifer, and layer 8 represented the Fox Hills aquifer. The numerical model was constructed using a uniform grid with square cells that are about 1 mile (1,600 meters) on each side with a total of about 657,000 active cells. Model calibration was completed by linking Parameter ESTimation (PEST) software with MODFLOW–NWT. The PEST software uses statistical parameter estimation techniques to identify an optimum set of input parameters by adjusting individual model input parameters and assessing the differences, or residuals, between observed (measured or estimated) data and simulated values. Steady-state model calibration consisted of attempting to match mean simulated values to measured or estimated values of (1) hydraulic head, (2) hydraulic head differences between model layers, (3) stream infiltration, and (4) discharge to streams. Calibration of the transient model consisted of attempting to match simulated and measured temporally distributed values of hydraulic head changes, stream base flow, and groundwater discharge to artesian flowing wells. Hydraulic properties estimated through model calibration included hydraulic conductivity, vertical hydraulic conductivity, aquifer storage, and riverbed hydraulic conductivity in addition to groundwater recharge and well skin. The ability of the numerical model to accurately simulate groundwater flow in the Williston Basin was assessed primarily by its ability to match calibration targets for hydraulic head, stream base flow, and flowing well discharge. The steady-state model also was used to assess the simulated potentiometric surfaces in the upper Fort Union aquifer, the lower Fort Union aquifer, and the Fox Hills aquifer. Additionally, a previously estimated regional groundwater-flow budget was compared with the simulated steady-state groundwater-flow budget for the Williston Basin. The simulated potentiometric surfaces typically compared well with the estimated potentiometric surfaces based on measured hydraulic head data and indicated localized groundwater-flow gradients that were topographically controlled in outcrop areas and more generalized regional gradients where the aquifers were confined. The differences between the measured and simulated (residuals) hydraulic head values for 11,109 wells were assessed, which indicated that the steady-state model generally underestimated hydraulic head in the model area. This underestimation is indicated by a positive mean residual of 11.2 feet for all model layers. Layer 7, which represents the lower Hell Creek aquifer, is the only layer for which the steady-state model overestimated hydraulic head. Simulated groundwater-level changes for the transient model matched within plus or minus 2.5 feet of the measured values for more than 60 percent of all measurements and to within plus or minus 17.5 feet for 95 percent of all measurements; however, the transient model underestimated groundwater-level changes for all model layers. A comparison between simulated and estimated base flows for the steady-state and transient models indicated that both models overestimated base flow in streams and underestimated annual fluctuations in base flow. The estimated and simulated groundwater budgets indicate the model area received a substantial amount of recharge from precipitation and stream infiltration. The steady-state model indicated that reservoir seepage was a larger component of recharge in the Williston Basin than was previously estimated. Irrigation recharge and groundwater inflow from outside the Williston Basin accounted for a relatively small part of total groundwater recharge when compared with recharge from precipitation, stream infiltration, and reservoir seepage. Most of the estimated and simulated groundwater discharge in the Williston Basin was to streams and reservoirs. Simulated groundwater withdrawal, discharge to reservoirs, and groundwater outflow in the Williston Basin accounted for a smaller part of total groundwater discharge. The transient model was used to simulate discharge to 571 flowing artesian wells within the model area. Of the 571 established flowing artesian wells simulated by the model, 271 wells did not flow at any time during the simulation because hydraulic head was always below the land-surface altitude. As hydraulic head declined throughout the simulation, 68 of these wells responded by ceasing to flow by the end of 2005. Total mean simulated discharge for the 571 flowing artesian wells was 55.1 cubic feet per second (ft 3 /s), and the mean simulated flowing well discharge for individual wells was 0.118 ft 3 /s. Simulated discharge to individual flowing artesian wells increased from 0.039 to 0.177 ft 3 /s between 1961 and 1975 and decreased to 0.102 ft 3 /s by 2005. The mean residual for 34 flowing wells with measured discharge was 0.014 ft 3 /s, which indicates the transient model overestimated discharge to flowing artesian wells in the model area. Model limitations arise from aspects of the conceptual model and from simplifications inherent in the construction and calibration of a regional-scale numerical groundwater-flow model. Simplifying assumptions in defining hydraulic parameters in space and hydrologic stresses and time-varying observational data in time can limit the capabilities of this tool to simulate how the groundwater-flow system responds to changes in hydrologic stresses, particularly at the local scale; nevertheless, the steady-state model adequately simulated flow in the uppermost principal aquifer systems in the Williston Basin based on the comparison between the simulated and estimated groundwater-flow budget, the comparison between simulated and estimated potentiometric surfaces, and the results of the calibration process.

Montana, North Dakota, South Dakota, Wyoming↗

300,000 yr history of water-table fluctuations at Wind Cave, South Dakota, USA—Scale, timing, and groundwater mixing in the Madison Aquifer

Deposits of calcite coating the lower passages of Wind Cave in the southern Black Hills of South Dakota were precipitated under phreatic conditions. Data from samples associated with a new cave survey and hydrologic studies indicate that past water tables within Wind Cave reached a maximum height of 45 m above modern levels but were mostly confined to 25 m or less. Uranium-series ages for basal layers deposited on weathered wall rock indicate subaerial conditions in this part of the cave persisted between 1000 and 300 ka. Ages and elevations of wall coatings and cave rafts establish a 300,000 yr paleohydrograph indicating that water-table highstands occurred during interglacial or interstadial-to-early glacial periods and lowstands occurred during full-glacial and stadial episodes. Isotopes of Sr, U, C, and O from dated calcite samples were obtained to evaluate potential shifts in paleo-groundwater composition. For comparison, Sr and U isotopic compositions were determined for modern groundwater from 18 sites previously classified into five hydrogeologic domains. Isotope data for different domains tend to cluster in separate fields, although several fields overlap. Compositions of Calcite Lake (informal name) water reflect modern recharge to shallow aquifers. In contrast, speleothem data indicate that paleo-groundwater highstands were not supported by increased infiltration associated with local recharge, or by upwelling from deeper Proterozoic sources. Instead, cave water was similar to deeper, warmer groundwater from the Madison Aquifer discharging at modern artesian springs flanking the southern Black Hills. Highstands were likely influenced by large-scale hydraulic processes associated with recharge to the Madison Aquifer under the Laurentide ice sheet on the northeast side of the Williston Basin, causing increased hydrostatic pressures in confined aquifers on the south side of the basin.

South Dakota↗

Simulating groundwater flow in the Mississippi Alluvial Plain with a focus on the Mississippi Delta

The Mississippi Alluvial Plain has become one of the most important agricultural regions in the United States but relies heavily on groundwater for irrigation. On average, more than 12 billion gallons are withdrawn daily from the Mississippi River Valley alluvial aquifer. Declining groundwater levels, especially in the Delta region of northwest Mississippi and the Cache and Grand Prairie regions of eastern Arkansas, have led to concerns about future sustainability. The U.S. Geological Survey Mississippi Alluvial Plain Project is focused on quantifying the groundwater system in the alluvial plain and the response of groundwater resources to future development. A key objective of the project is to provide updated groundwater flow models supported by extensive data collection and analyses. MODFLOW 6, PEST++, and several open-source python packages were used to develop a simplified, faster running version of the Mississippi Embayment Regional Aquifer Study model that can provide boundary conditions for local inset models, including the Mississippi Delta model described in this report. An automated workflow was used for model construction, history matching, and development of baseline future climate scenarios. The models incorporate information from a Soil-Water-Balance code simulation of the terrestrial water balance, metering-based estimates of water use from thousands of wells, measured and estimated streamflow and stages, and the largest airborne electromagnetic survey flown to date in the United States. Baseline scenarios for the Mississippi Delta under potential future climates were constructed using recharge, surface runoff and irrigation pumping forcings from a future version of the Soil-Water-Balance model, driven by downscaled temperature and precipitation output from 10 general circulation model simulations, including high and moderate carbon emissions pathways. Results indicate a complex water balance that varies in time and space in terms of the terrestrial recharge, stream leakage, and regional groundwater flow components, which are affected by seasonal forcings, human activity, and alluvial geomorphology. The general circulation model outputs indicate a continued rise in average temperatures but no clear precipitation trend. Increased crop water demand is anticipated from the higher temperatures, resulting in increased irrigation withdrawals to sustain current levels of irrigated agriculture. Simulated drawdowns in groundwater levels at the mid-21st century vary greatly. Under moderate or wet climate scenarios, and in parts of the aquifer that are well connected to surface water, little to no additional drawdown is anticipated. Under dry or warm scenarios, drawdowns of as much as 10 meters or more are possible in parts of the aquifer that are relatively disconnected from surface water. Under dry or warm scenarios, the portion of the Delta with greater than 60 feet of saturated thickness could be reduced from near 100 percent currently (2018) to 80–90 percent by mid-century. Future simulations with the model could include alternative management scenarios to identify options for improving groundwater sustainability. The automated model construction workflows are designed to facilitate regular updating, making this a “living” framework that the Mississippi Department of Environmental Quality and other stakeholders can use for adaptive management going forward.

Arkansas, Louisiana, Mississippi↗

Water-quality comparison of the Gulf Coast aquifer system at various scales in Texas from National Water-Quality Assessment groundwater studies, 2013–15

One of the objectives of the U.S. Geological Survey National Water-Quality Assessment (NAWQA) Project is to assess groundwater quality in aquifers that are important sources of drinking water such as the coastal lowlands aquifer system, which is often referred to in Texas as the “Gulf Coast aquifer system.” The Gulf Coast aquifer system extends from Louisiana to Mexico and is a source of groundwater for several cities including Houston, Tex. The NAWQA groundwater studies in Texas in 2013–15 that assessed the Gulf Coast aquifer system included Principal Aquifer Surveys (PAS), Major Aquifer Studies (MAS), and Land Use Studies (LUS). These three study types are based on sampling networks of wells distributed in an area of interest. The PAS networks typically consist of public-supply wells that are relatively deep, the MAS networks typically consist of domestic-supply wells that are intermediate in depth, and the LUS networks typically consist of monitoring wells that are relatively shallow.

Texas↗

Groundwater hydrology and estimation of horizontal groundwater flux from the Rio Grande at selected locations in Albuquerque, New Mexico, 2003-9

The Albuquerque, New Mexico, area has two principal sources of water: groundwater from the Santa Fe Group aquifer system and surface water from the San Juan-Chama Diversion Project. From 1960 to 2002, groundwater withdrawals from the Santa Fe Group aquifer system have caused water levels to decline more than 120 feet in some places within the Albuquerque area, resulting in a great deal of interest in quantifying the river-aquifer interaction associated with the Rio Grande. In 2003, the U.S. Geological Survey in cooperation with the Bureau of Reclamation, the Middle Rio Grande Endangered Species Collaborative Program, and the U.S. Army Corps of Engineers began a detailed characterization of the hydrogeology of the Rio Grande riparian corridor in the Albuquerque, New Mexico, area to provide hydrologic data and enhance the understanding of rates of water leakage from the Rio Grande to the alluvial aquifer, groundwater flow through the aquifer, and discharge of water from the aquifer to the riverside drains. A simple conceptual model of flow indicates that the groundwater table gently slopes from the Rio Grande towards riverside drains and the outer boundaries of the inner valley. Water infiltrating from the Rio Grande initially moves vertically below the river, but, as flow spreads farther into the Rio Grande inner valley alluvial aquifer, flow becomes primarily horizontal. The slope of the water-table surface may be strongly controlled by the riverside drains and influenced by other more distal hydrologic boundary conditions, such as groundwater withdrawals by wells. Results from 35 slug tests performed in the Rio Grande inner valley alluvial aquifer during January and February 2009 indicate that hydraulic-conductivity values ranged from 5 feet per day to 160 feet per day with a median hydraulic-conductivity for all transects of 40 feet per day. Median annual horizontal hydraulic gradients in the Rio Grande inner valley alluvial aquifer ranged from 0.011 to 0.002. Groundwater fluxes through the alluvial aquifer calculated by using median slug-test results (qm slug ) and Darcy's law ranged from about 0.1 feet per day to about 0.7 feet per day. Groundwater fluxes calculated by using the Suzuki-Stallman method (qm heat ) ranged from 0.52 feet per day to 0.23 feet per day. Results from the Darcy's law and Suzuki-Stallman flux calculations were compared to discharge measured in riverside drains on both sides of the river north of the Montaño Bridge on February 26, 2009. Flow in the Corrales Riverside Drain increased by 1.4 cubic feet per second from mile 2 to mile 4, about 12 cubic feet per day per linear foot of drain. Flow in the Albuquerque Riverside Drain increased by 15 cubic feet per second between drain miles 0 and 3, about 82 cubic feet per day per linear foot of drain. The flux of water from the river to the aquifer was calculated to be 2.2 cubic feet per day per linear foot of river by using the median qm slug of 0.09 feet per day at Montaño transects west of the river. The total flux was calculated to be 6.0 cubic feet per day per linear foot of river by using the mean(qm heat of 0.24 feet per day for the Montaño transects west of the river. Assuming the Corrales Riverside Drain intercepted all of this flow, the qm slug or qm heat fluxes account for 18 to 50 percent, respectively, of the increase of flow in the drain. The flux of water from the river to the aquifer was calculated to be 15 cubic feet per day per linear foot of river by using the median qm slug of 0.30 feet per day at the Montaño transects east of the river. The flux of water from the river to the aquifer was calculated to be 17 cubic feet per day per linear foot of river by using the mean flux calculated from the Suzuki-Stallman method for the Montaño East transects of 0.34 feet per day. Assuming the Albuquerque Riverside Drain intercepted all this flow, the qm slug or (qm heat fluxes would only account for 18 to 21 percent, respectively, of the increase in flow in the drain. The comparison of these results with those of previous investigations suggests that calculated flux through the Rio Grande inner valley alluvial aquifer is strongly scale dependent and that the thickness of aquifer through which river water flows may be greater than indicated by the vertical temperature profiles.

Albuquerque, New Mexico↗

Geohydrology and water quality of the stratified-drift aquifers in West Branch Cayuga Inlet and Fish Kill Valleys, Newfield, Tompkins County, New York

From 2011 to 2016, the U.S. Geological Survey, in cooperation with the Town of Newfield and the Tompkins County Planning Department, performed a study of the stratified-drift aquifers in the West Branch Cayuga Inlet and Fish Kill Valleys in Newfield, Tompkins County, New York. Both confined and unconfined aquifers were identified, mostly in the valleys. The confined aquifer consists of a discontinuous sand and gravel layer that overlies bedrock and is commonly confined by overlying fine-grained sediments. The unconfined aquifer consists of surficial ice contact sand and gravel, alluvial silt, sand and gravel, and areas where several large tributary streams deposited alluvial fans in the valley, all of which were deposited during and after the last glacial recession. The unconfined aquifers are primarily recharged by direct infiltration of precipitation at the land surface, by surface runoff and shallow subsurface flow from adjacent hillsides, and by seepage loss from streams crossing the aquifer, especially on alluvial fans. The confined aquifers are primarily recharged by groundwater stored in the overlying sand and gravel aquifer that slowly seeps downward through the underlying confining layer. Other sources of recharge are precipitation that falls directly on the surficial confining unit and adjacent valley walls, which then slowly seeps downward and enters the confined aquifer, and groundwater flow from bordering till and bedrock and from bedrock below the valley. There may also be some recharge where confining units are absent or where parts of the confining units contain sediments with moderate permeability. The groundwater naturally discharges to the Fish Kill and West Branch Cayuga Inlet streams and to wetlands overlying the aquifer boundaries, with additional losses due to evapotranspiration. Groundwater is pumped from the aquifers by domestic, municipal, and agricultural wells. Approximately 57.9 million gallons per year was withdrawn from the stratified-drift (sand and gravel) aquifers. Groundwater samples were collected from 11 wells, and surface water samples were collected at 2 sites, one each from Fish Kill and West Branch Cayuga Inlet. None of the common ions (for example, sodium, chloride, and magnesium) exceeded existing drinking water standards at either surface water site. The concentration of nitrate plus nitrite detected was 0.4 milligram per liter as nitrogen in the West Branch Cayuga Inlet site. Total phosphorus was detected at 0.01 milligram per liter as phosphate for both sites. Of the 11 wells sampled, 8 were finished in confined sand and gravel aquifers, 1 was finished in unconfined sand and gravel, and 2 were finished in shale bedrock. Groundwater quality in the study area generally met Federal and State drinking-water standards. However, of the 11 samples taken, 2 exceeded the U.S. Environmental Protection Agency drinking water advisory taste threshold of 20 milligrams per liter for sodium, 8 exceeded the secondary maximum contaminant level of 300 micrograms per liter for iron, and 9 exceeded the secondary maximum contaminant level of 50 micrograms per liter for manganese.

New York↗

Groundwater quality in the Mississippi embayment–Texas coastal uplands aquifer system, south-central United States

Groundwater provides nearly 50 percent of the Nation’s drinking water. To help protect this vital resource, the U.S. Geological Survey (USGS) National Water-Quality Assessment (NAWQA) Project assesses groundwater quality in aquifers that are important sources of drinking water (Burow and Belitz, 2014). The Mississippi embayment–Texas coastal uplands aquifer system constitutes one of the important aquifer systems being evaluated.

Mississippi Embayment–Texas Coastal Uplands Aquife↗

The use of simulation and multiple environmental tracers to quantify groundwater flow in a shallow aquifer

Measurements of the concentrations of chlorofluorocarbons (CFCs), tritium, and other environmental tracers can be used to calculate recharge ages of shallow groundwater and estimate rates of groundwater movement. Numerical simulation also provides quantitative estimates of flow rates, flow paths, and mixing properties of the groundwater system. The environmental tracer techniques and the hydraulic analyses each contribute to the understanding and quantification of the flow of shallow groundwater. However, when combined, the two methods provide feedback that improves the quantification of the flow system and provides insight into the processes that are the most uncertain. A case study near Locust Grove, Maryland, is used to investigate the utility of combining groundwater age dating, based on CFCs and tritium, and hydraulic analyses using numerical simulation techniques. The results of the feedback between an advective transport model and the estimates of groundwater ages determined by the CFCs improve a quantitative description of the system by refining the system conceptualization and estimating system parameters. The plausible system developed with this feedback between the advective flow model and the CFC ages is further tested using a solute transport simulation to reproduce the observed tritium distribution in the groundwater. The solute transport simulation corroborates the plausible system developed and also indicates that, for the system under investigation with the data obtained from 0.9-m-long (3-foot-long) well screens, the hydrodynamic dispersion is negligible. Together the two methods enable a coherent explanation of the flow paths and rates of movement while indicating weaknesses in the understanding of the system that will require future data collection and conceptual refinement of the groundwater system.

Water Resources Research↗

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2024

Potentiometric-surface and depth-to-water maps for spring 2024 were created for the Mississippi River Valley alluvial aquifer using groundwater-altitude data from 1,151 wells completed in the Mississippi River Valley alluvial aquifer and from the altitude of the top of the water surface in area rivers from 160 U.S. Geological Survey and U.S. Army Corps of Engineers streamgages. The potentiometric-surface and depth-to-water maps for 2024 were created to support investigations to characterize the Mississippi River Valley alluvial aquifer as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the Mississippi River Valley alluvial aquifer for spring 2024 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 330 feet (ft) above the North American Vertical Datum of 1988. The regional direction of the groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient direction was into the depression, and near rivers, where the groundwater gradient direction was either from the aquifer to the river or from the river into the aquifer. Groundwater depressions in the potentiometric-surface map are in the lower one-half of the Cache region and in most of the Grand Prairie and Delta regions. Depth to water by well in the Mississippi River Valley alluvial aquifer in spring 2024 ranged from 0.27 ft above land surface to 145.65 ft below land surface.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗

Extraction and development of inset models in support of groundwater age calculations for glacial aquifers

The U.S. Geological Survey developed a regional model of Lake Michigan Basin (LMB). This report describes the construction of five MODFLOW inset models extracted from the LMB regional model and their application using the particle-tracking code MODPATH to simulate the groundwater age distribution of discharge to wells pumping from glacial deposits. The five study areas of the inset model correspond to 8-digit hydrologic unit code (HUC8) basins. Two of the basins are tributary to Lake Michigan from the east, two are tributary to the lake from the west, and one is just west of the western boundary of the Lake Michigan topographic basin. The inset models inherited many of the inputs to the parent LMB model, including the hydrostratigraphy and layering scheme, the hydraulic conductivity assigned to bedrock layers, recharge distribution, and water use in the form of pumping rates from glacial and bedrock wells. The construction of the inset models entailed modifying some inputs, most notably the grid spacing (reduced from cells 5,000 feet on a side in the parent LMB model to 500 feet on a side in the inset models). The refined grid spacing allowed for more precise location of pumped wells and more detailed simulation of groundwater/surface-water interactions. The glacial hydraulic conductivity values, the top bedrock surface elevation, and the surface-water network input to the inset models also were modified. The inset models are solved using the MODFLOW–NWT code, which allows for more robust handling of conditions in unconfined aquifers than previous versions of MODFLOW. Comparison of the MODFLOW inset models reveals that they incorporate a range of hydrogeologic conditions relative to the glacial part of the flow system, demonstrated by visualization and analysis of model inputs and outputs and reflected in the range of ages generated by MODPATH for existing and hypothetical glacial wells. Certain inputs and outputs are judged to be candidate predictors that, if treated statistically, may be capable of explaining much of the variance in the simulated age metrics. One example of a predictor that model results indicate strongly affects simulated age is the depth of the well open interval below the simulated water table. The strength of this example variable as an overall predictor of groundwater age and its relation to other predictors can be statistically tested through the metamodeling process. In this way the inset models are designed to serve as a training area for metamodels that estimate groundwater age in glacial wells, which in turn will contribute to ongoing studies, under the direction of the U.S. Geological Survey National Water Quality Assessment, of contaminant susceptibility of shallow groundwater across the glacial aquifer system.

Illinois, Indiana, Michigan, Wisconsin↗