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.

SourceResearch connectionFrom yearThrough year

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

At least 865 records · Page 48Linked to original sources

USGS2022

Groundwater, surface-water, and water-chemistry data, Black Mesa area, northeastern Arizona—2018–2019

The Navajo (N) aquifer is an extensive aquifer and the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Water availability is an important issue in the Black Mesa area because of the arid climate, past industrial water use, and continued water requirements for municipal use by a growing population. Precipitation in the area typically ranges from less than 6 to more than 16 inches per year depending on location. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents results of data collected as part of the monitoring program in the Black Mesa area from calendar year 2019, and additionally uses streamflow statistics from November and December 2018. The monitoring program includes measurements of (1) groundwater withdrawals (pumping), (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater chemistry. In calendar year 2019, total groundwater withdrawals were estimated to be 3,070 acre-feet (acre-ft), industrial withdrawals were 670 acre-ft, and municipal withdrawals were estimated to be 2,400 acre-ft. Total withdrawals during 2019 were about 58 percent less than total withdrawals in 2005 because of Peabody Western Coal Company’s discontinued use of water to transport coal in a coal slurry pipeline after 2005 and cessation of mining operations in 2019. Water levels measured in 2019 from wells completed in the unconfined areas of the N aquifer within the Black Mesa area showed a decline in 10 of 16 wells when compared with water levels from the prestress period (prior to 1965). The changes in water levels across all 16 wells ranged from +8.2 feet (ft) to −40.0 ft, and the median change was −1.7 ft. Water levels also showed decline in 16 of 18 wells measured in the confined area of the aquifer when compared to the prestress period. The median change for the confined area of the aquifer was −38.8 ft, with changes across all 18 wells ranging from +12.9 ft to −185.0 ft. Spring flow was measured at four springs in 2019. Flow fluctuated during the period of record for Burro Spring and Pasture Canyon Spring, but a decreasing trend was statistically significant (p<0.05) at Moenkopi School Spring and Unnamed Spring near Dennehotso. Discharge at Burro Spring has remained relatively constant since it was first measured in the 1980s and discharge at Pasture Canyon Spring has fluctuated for the period of record. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976 to 2019), Dinnebito Wash near Sand Springs 09401110 (1993 to 2019), Polacca Wash near Second Mesa 09400568 (1994 to 2019), and Pasture Canyon Springs 09401265 (2004 to 2019). Median winter flows (November through February) of each winter were used as an estimate of the amount of groundwater discharge at the above-named sites. For the period of record, the median winter flows have generally remained constant at Polacca Wash and Pasture Canyon Springs, whereas a decreasing trend was indicated at Moenkopi Wash and Dinnebito Wash. In 2019, water samples collected from four springs and three wells in the Black Mesa area were analyzed for selected chemical constituents. Results from the four springs were compared with previous analyses from the same springs. Concentrations of dissolved solids, chloride, and sulfate increased at Moenkopi School Spring during the more than 30 years of record at that site. Concentrations of dissolved solids, chloride, and sulfate at Pasture Canyon Spring have not varied significantly (p>0.05) since the early 1980s, and there is no increasing or decreasing trend in those data. Concentrations of dissolved solids, chloride, and sulfate at Unnamed Spring near Dennehotso have varied for the period of record, but there is no statistical trend in the data. Concentrations of dissolved solids and chloride at Burro Spring have varied for the period of record, but there is no statistical trend in the data; however, concentrations of sulfate from Burro Spring now show a trend towards lower concentrations. No statistical trend tests were performed for the three wells sampled in 2019 since less historical water-quality data were available for comparison.

Arizona↗

USGS1984

Groundwater-flow parameter estimation and quality modeling of the Equus Beds aquifer in Kansas, U.S.A.

The salinity problems created in the Burrton area as a result of poor oil-field brine disposal practices of the past continue to be a major concern to the area depending on the Equus Beds aquifer for water, including the City of Wichita, Kansas. In this paper, an attempt is made to predict where and how fast the brine plume will move in this area, and what the average chloride concentrations in different parts of the aquifer are. In order to make such predictions, it was necessary to get a calibrated model of the groundwater-flow velocity field. Multiple regression analysis is used for parameter estimation of the steady-state groundwater-flow equation applied in the most critical area of the Equus Beds aquifer. Results of such an analysis produced a correlation coefficient of 0.992 between calculated and observed values of hydraulic head. A chloride transport modeling effort is then carried out despite some serious data deficiencies, the significance of which are evaluated through sensitivity analysis. Thus, starting with the quasi steady-state conditions of the early 1940's, it was possible to match the present chloride distribution satisfactorily. Chloride concentration predictions made for the year 2000 indicate that the quality of the Wichita well-field waters will not generally deteriorate from their present condition by that time.

Kansas↗

USGS2024

Groundwater, surface-water, and water-chemistry data, Black Mesa area, northeastern Arizona—2019–2021

The Navajo (N) aquifer is an extensive aquifer and the primary source of groundwater in the 5,400-square-mile Black Mesa area in northeastern Arizona. Water availability is an important issue in the Black Mesa area because of the arid climate, past industrial water use, and continued water requirements for municipal use by a growing population. Precipitation in the area typically ranges from less than 6 to more than 16 inches per year, depending on location. The U.S. Geological Survey water-monitoring program in the Black Mesa area began in 1971 and provides information about the long-term effects of groundwater withdrawals from the N aquifer for industrial and municipal uses. This report presents the results of data collected as part of the monitoring program in the Black Mesa area from calendar years 2020–2021 and, additionally, uses streamflow statistics from November and December 2019. The monitoring program includes measurements of (1) groundwater withdrawals (pumping), (2) groundwater levels, (3) spring discharge, (4) surface-water discharge, and (5) groundwater chemistry. In calendar year 2020, total groundwater withdrawals were estimated to be 2,680 acre-feet (acre-ft), and, in 2021, total withdrawals were estimated to be 2,570 acre-ft. Total withdrawals during 2021 were about 65 percent less than total withdrawals in 2005 because the Peabody Western Coal Company discontinued its use of water to transport coal in a coal slurry pipeline after 2005 and ceased mining operations in 2019. Owing to Navajo Nation and Hopi Reservation access restrictions during the Coronavirus pandemic, water levels were not collected from municipal wells in 2020 or 2021. Water levels measured in 2021 from wells completed in the unconfined areas of the N aquifer within the Black Mesa area showed a decline in 7 of 13 wells when compared with water levels from the prestress period (prior to 1965). The changes in water levels across all 13 wells ranged from +8.4 feet (ft) to −42.4 ft, and the median change was −0.4 ft. Water levels also showed decline in 11 of 12 wells measured in the confined area of the aquifer when compared to the prestress period. The median change for the confined area of the aquifer was −25.9 ft, with changes across all 12 wells ranging from +17.3 ft to −133.7 ft. Spring flow was measured at four springs between 2020 and 2021. Flow fluctuated during the period of record for Burro Spring and Pasture Canyon Spring, but a decreasing trend was statistically significant (p<0.05) at Moenkopi School Spring and Unnamed Spring near Dennehotso, Arizona. Discharge at Burro Spring has remained relatively constant since it was first measured in the 1980s, and discharge at Pasture Canyon Spring has fluctuated for the period of record. Continuous records of surface-water discharge in the Black Mesa area were collected from streamflow-gaging stations at the following sites: Moenkopi Wash at Moenkopi 09401260 (1976–2021), Dinnebito Wash near Sand Springs 09401110 (1993–2020), Polacca Wash near Second Mesa 09400568 (1994–2020), and Pasture Canyon Springs 09401265 (2004–2021). Median winter flows (November through February) of each winter were used as an estimate of the amount of groundwater discharge at the above-named sites. For the period of record, the median winter flows have generally remained constant at Polacca Wash and Pasture Canyon Springs, whereas a decreasing trend was observed at Moenkopi Wash and Dinnebito Wash. In 2020 and 2021, water samples were collected from a total of four springs in the Black Mesa area and analyzed for selected chemical constituents. Results from the four springs were compared with previous analyses from the same springs. Dissolved solids, chloride, and sulfate concentrations increased at Moenkopi School Spring during the more than 30 years of record at that site. Concentrations of dissolved solids and sulfate at Pasture Canyon Spring have not varied significantly (p>0.05) since the early 1980s, and there is no increasing or decreasing trend in those data. However, concentrations of chloride from Pasture Canyon Spring show a diminishing trend. Concentrations of dissolved solids, chloride, and sulfate at Unnamed Spring near Dennehotso have varied for the period of record, but there is no statistical trend in the data. Concentrations of dissolved solids at Burro Spring have varied for the period of record, but there is no statistical trend in the data. However, concentrations of chloride and sulfate from Burro Spring show a trend towards lower concentrations.

Arizona↗

USGS2014

The Effect of modeled recharge distribution on simulated groundwater availability and capture

Simulating groundwater flow in basin-fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin-fill aquifers by direct infiltration and transport through faults and fractures in the high-elevation areas, by flowing overland through high-elevation areas to infiltrate at basin-fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin-fill aquifers by calibrating a groundwater-flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady-state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.

Arizona↗

USGS2025

Hydrogeologic investigation, framework, and conceptual flow model of the Antlers aquifer, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. Every 20 years, the Oklahoma Water Resources Board is required to update the hydrologic investigation on which the maximum annual yield determinations were based. The maximum annual yield allocated per acre of land is used to set the equal-proportionate share pumping rate. The maximum annual yield of 5,913,600 acre-feet per year and equal-proportionate-share of 2.1 acre-feet per acre per year currently (2025) in place for the Antlers aquifer were issued by the Oklahoma Water Resources Board on February 14, 1995. Because more than 20 years have elapsed since the 1995 final order for the Antlers aquifer was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, completed an in-depth hydrologic study that included a hydrogeologic framework and conceptual groundwater-flow model for the 1980–2022 study period. The results of an analysis of land use, long-term climate patterns, streamflow and base-flow patterns, historical groundwater use, as well as groundwater-level fluctuations across the Antlers aquifer are described. In addition, groundwater quality was analyzed for total dissolved solids concentrations and major ions for the Antlers aquifer. An updated hydrogeologic framework was developed that included refining the aquifer boundary in Oklahoma, the creation of new potentiometric surface and saturated thickness of fresh groundwater maps, one multiple-well aquifer test, slug tests, and an analysis of lithologic logs across the aquifer. A conceptual groundwater flow model and water budget were developed by incorporating estimates of recharge from precipitation, saturated-zone evapotranspiration, streambed seepage, lateral groundwater flows, vertical leakage, and withdrawals from groundwater wells.

Oklahoma, Texas↗

USGS2008

Distribution and variability of redox zones controlling spatial variability of arsenic in the Mississippi River Valley alluvial aquifer, southeastern Arkansas

Twenty one of 118 irrigation water wells in the shallow (25-30??m thick) Mississippi River Valley alluvial aquifer in the Bayou Bartholomew watershed, southeastern Arkansas had arsenic (As) concentrations (< 0.5 to 77????g/L) exceeding 10????g/L. Sediment and groundwater samples were collected and analyzed from the sites of the highest, median, and lowest concentrations of As in groundwater in the alluvial aquifers located at Jefferson County, Arkansas. A traditional five-step sequential extraction was performed to differentiate the exchangeable, carbonate, amorphous Fe and Mn oxide, organic, and hot HNO3-leachable fraction of As and other compounds in sediments. The Chao reagent (0.25??M hydroxylamine hydrochloride in 0.25??M HCl) removes amorphous Fe and Mn oxides and oxyhydroxides (present as coatings on grains and amorphous minerals) by reductive dissolution and is a measure of reducible Fe and Mn in sediments. The hot HNO3 extraction removes mostly crystalline metal oxides and all other labile forms of As. Significant total As (20%) is complexed with amorphous Fe and Mn oxides in sediments. Arsenic abundance is not significant in carbonates or organic matter. Significant (40-70????g/kg) exchangeable As is only present at shallow depth (0-1??m below ground surface). Arsenic is positively correlated to Fe extracted by Chao reagent (r = 0.83) and hot HNO3 (r = 0.85). Arsenic extracted by Chao reagent decreases significantly with depth as compared to As extracted by hot HNO3. Fe (II)/Fe (the ratio of Fe concentration in the extracts of Chao reagent and hot HNO3) is positively correlated (r = 0.76) to As extracted from Chao reagent. Although Fe (II)/Fe increases with depth, the relative abundance of reducible Fe decreases noticeably with depth. The amount of reducible Fe, as well as As complexed to amorphous Fe and Mn oxides and oxyhydroxides decreases with depth. Possible explanations for the decrease in reducible Fe and its complexed As with depth include historic flushing of As and Fe from hydrous ferric oxides (HFO) by microbially-mediated reductive dissolution and aging of HFO to crystalline phases. Hydrogeochemical data suggests that the groundwater in the area falls in the mildly reducing (suboxic) to relatively highly reducing (anoxic) zone, and points to reductive dissolution of HFO as the dominant As release mechanism. Spatial variability of gypsum solubility and simultaneous SO42- reduction with co-precipitation of As and sulfide is an important limiting process controlling the concentration of As in groundwater in the area. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Contaminant Hydrology↗

USGS2011

Groundwater conditions and studies in the Augusta–Richmond County area, Georgia, 2008–2009

Groundwater studies and monitoring efforts conducted during 2008–2009, as part of the U.S. Geological Survey (USGS) Cooperative Water Program with the City of Augusta in Richmond County, Georgia, provided data for the effective management of local water resources. During 2008–2009 the USGS completed: (1) installation of three monitoring wells and the collection of lithologic and geophysical logging data to determine the extent of hydrogeologic units, (2) collection of continuous groundwater-level data from wells near Well Fields 2 and 3, (3) collection of synoptic groundwater-level measurements and construction of potentiometric-surface maps in Richmond County to establish flow gradients and groundwater-flow directions in the Dublin and Midville aquifer systems, (4) completion of a 24-hour aquifer test to determine hydraulic characteristics of the lower Dublin aquifer, and upper and lower Midville aquifers in Well Field 2, and (5) collection of groundwater samples from selected wells in Well Field 2 for laboratory analysis of volatile organic compounds and groundwater tracers to assess groundwater quality and estimate the time of groundwater recharge. Potentiometric-surface maps of the Dublin and Midville aquifer systems for 2008–2009 indicate that the general groundwater flow direction within Richmond County is eastward toward the Savannah River, with the exception of the area around Well Field 2, where pumping interrupts the eastward flow of water toward the Savannah River and causes flow lines to bend toward the center of pumping. Results from a 24-hour aquifer test conducted in 2009 within the upper and lower Midville aquifers at Well Field 2 indicated a transmissivity and storativity for the upper and lower Midville aquifers, combined, of 4,000 feet-squared per day and 2x10 -4 , respectively. The upper and lower Midville aquifers and the middle lower Midville confining unit, which is 85-feet thick in this area, yielded horizontal hydraulic conductivity and specific storage values of about 45 feet per day and 2x10 -6 ft -1 , respectively. Results from the 24-hour aquifer test also indicate a low horizontal hydraulic conductivity for the lower Dublin aquifer of less than 1 foot per day. Of the 35 volatile organic compounds (VOCs) analyzed in 23 groundwater samples during 2008–2009, only six were detected above laboratory reporting limits in samples from eight wells. No concentration in groundwater samples collected during 2008–2009 exceeded drinking water standards. Trichloroethene had the maximum VOC concentration (1.9 micrograms per liter) collected from a water sample during 2008–2009. Water-quality sampling of several wells near Well Field 2 indicate that, while in operation, the northernmost production well might have diverted groundwater, containing low levels of trichloroethene from at least two other production wells. Analysis of sulfur hexafluoride data indicate the average year of recharge ranges between 1981 and 1984 for water samples from five wells open to the upper and lower Midville aquifers, and 1991 for a water sample from one shallow well open to the lower Dublin aquifer. All of these ages suggest a short flow path and nearby source of contamination. The actual source of low levels of VOCs at Well Field 2 remains unknown. Three newly installed monitoring wells indicate that hydrogeologic units beneath Well Fields 2 and 3 are composed of sand and clay layers. Hydrogeologic units, encountered at Well Field 2, in order of increasing depth are the lower Dublin confining unit, lower Dublin aquifer, upper Midville confining unit, upper Midville aquifer, lower Midville confining unit, and lower Midville aquifer. West of Well Field 3, hydrogeologic units, in order of increasing depth are the Upper Three Runs aquifer, Gordon confining unit, Gordon aquifer, lower Dublin confining unit, lower Dublin aquifer, upper Midville confining unit, upper Midville aquifer, lower Midville confining unit, and lower Midville aquifer.

Georgia↗

USGS2023

Updates to the regional groundwater-flow model of the New Jersey Coastal Plain, 1980–2013

A 21-layer three-dimensional transient groundwater-flow model of the New Jersey Coastal Plain was developed and calibrated by the U.S. Geological Survey (USGS) in cooperation with the New Jersey Department of Environmental Protection to simulate groundwater-flow conditions during 1980–2013, incorporating average annual groundwater withdrawals and average annual groundwater recharge. This model is the third version of the New Jersey Coastal Plain regional groundwater-flow model that was initially developed as part of the USGS Regional Aquifer System Analysis (RASA) program. The model simulates groundwater flow in 11 aquifers and 10 intervening confining units of the New Jersey Coastal Plain to provide a regional overview of groundwater conditions. Averaged groundwater withdrawal data for 1980 to 2013 were used in the model. The 11 aquifers in New Jersey are, from shallowest to deepest, the Holly Beach water-bearing zone and the confined Cohansey aquifer in Cape May County; the Rio Grande water-bearing zone; the Atlantic City 800-foot sand; the Piney Point, Vincentown, and Wenonah-Mount Laurel aquifers; the Englishtown aquifer system; and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy (PRM) aquifer system. The model was developed with the MODFLOW–2005 numerical code and the UCODE parameter estimation technique and calibrated using water-level and base-flow observations. A total of 3,453 water-level observations from 392 wells in New Jersey and 48 wells in Delaware from 1983 to 2013 were used in model calibration, which includes historical water-level trends for 29 wells in New Jersey during 1980–2013 presented in time-series hydrographs. In addition, derived observations also were included by calculating the vertical gradient at 33 pairs of nested observation wells in New Jersey, for a total of 210 observations. Changes in water levels over time were calculated for 134 wells in New Jersey and four wells in Delaware where water levels had varied substantially (approximately 10 ft) over the 30-year span of synoptic water-level measurements, for a total of 767 observations. A total of 1,485 base-flow observations in 47 surface-water basins in New Jersey from 1980 to 2013 were used in model calibration. Updates to the groundwater-flow model include the conversion to a fully three-dimensional model from the previous quasi-three-dimensional model. The new model will allow for potential future uses such as particle tracking or simulation of variable-density groundwater flow that could not be accomplished with earlier versions of the model. Spatially and temporally variable recharge estimated by using a soil-water balance model resulted in a spatially and temporally finer discretization. The Rio Grande water-bearing zone was added to the model as an aquifer layer to refine estimates of simulated flow in Atlantic and Cape May Counties, New Jersey. Hydrogeologic parameters were updated to include the confining units in New Jersey and corresponding hydrogeologic units in Delaware and eastern Maryland. The simulated water levels for the New Jersey Coastal Plain aquifers were compared to water-level measurements made during 1980–2013. The average residual for 4,243 water-level observations for New Jersey (simulated water levels minus measured water levels) is 1.5 feet. The simulated water-level contours for the confined aquifers for 2013 were compared to potentiometric surfaces produced from water levels measured during 2013. Simulated water levels generally matched the 2013 potentiometric surfaces of the confined aquifers in the areas of large withdrawals. Hydrographs of wells in the confined Coastal Plain aquifers of New Jersey show that simulated water levels generally match the magnitude and seasonal variation of the observed water levels. Hydrographs of base flow for the 47 streamgaging stations in New Jersey indicate that most of the simulated and estimated data match reasonably well. Groundwater withdrawals are an important resource for water supply, agricultural, industrial, and commercial needs in the New Jersey Coastal Plain. Groundwater withdrawals from the New Jersey Coastal Plain aquifers have resulted in persistent, regionally extensive cones of depression in the Englishtown aquifer system and Wenonah-Mount Laurel aquifer in Ocean and Monmouth Counties; Wenonah-Mount Laurel and upper, middle, and lower PRM aquifers in Camden County; and Atlantic City 800-foot sand in Atlantic County. Because hydrologic stresses and water-management needs change with time, periodic updates to the groundwater-flow model are required to provide current information about hydrologic conditions in the New Jersey Coastal Plain and to maintain its usefulness as a tool to manage water resources and develop water-resource strategies. The current updates will support the continued application of this model as a tool for evaluating the regional effects of changes in groundwater withdrawals and of current and potential future water-management strategies on groundwater levels in the New Jersey Coastal Plain.

New Jersey Coastal Plain↗

USGS2016

Hydrogeological framework, numerical simulation of groundwater flow, and effects of projected water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma

This report describes a study of the hydrology, hydrogeological framework, numerical groundwater-flow models, and results of simulations of the effects of water use and drought for the Beaver-North Canadian River alluvial aquifer, northwestern Oklahoma. The purpose of the study was to provide analyses, including estimating equal-proportionate-share (EPS) groundwater-pumping rates and the effects of projected water use and droughts, pertinent to water management of the Beaver-North Canadian River alluvial aquifer for the Oklahoma Water Resources Board. The Beaver-North Canadian River alluvial aquifer consists of unconsolidated sand, gravel, silt, and clay in varying proportions that underlies the Beaver and North Canadian River Valleys for approximately 175 miles (mi) from the Oklahoma Panhandle to the western edge of Oklahoma City in central Oklahoma. The aquifer as delineated for this study varies from 4 to 12 mi wide and is as thick as 308 feet (ft) in the northwest where the aquifer includes the Ogallala Formation. There are two distinct but in most areas hydraulically connected alluvial units that compose the Beaver-North Canadian River alluvial aquifer: a Quaternary-age topographically higher terrace deposit and a topographically lower, younger alluvium along the active river channel that includes active and Quaternary-age alluvium. The Beaver River composes the headwaters of the North Canadian River, which begins at the confluence of the Beaver River and Wolf Creek. The aquifer is divided for water management into two geographic areas: Reach I upstream from Canton Dam and Reach II downstream from Canton Dam. Reach I covers an area of approximately 874 square miles (mi 2 ), and Reach II covers an area of approximately 371 mi 2 . The Beaver-North Canadian River alluvial aquifer crosses several climatic zones, from semiarid in the west to continental subhumid in the east. Mean annual precipitation varies from 23.5 inches (in.) in the western part of this aquifer to 35.7 in. in the east. Surface-water demands were met through numerous temporary and permanent surface-water diversions from the Beaver and North Canadian Rivers during the period of study. During the study period, seven diversions removed a mean annual 2,000 acre-feet (acre-ft) of water from Reach I. There were 14 diversions from Reach II with a mean annual permitted volume of approximately 81,000 acre-ft, including diversion into the Lake Hefner Canal for the Oklahoma City public water supply. During the period of this study, 17 temporary surface-water diversion permits were active in Reach I, with total permitted volumes of 2,000 acre-ft, and 41 diversions were active in Reach II, with total permitted volumes of 38,000 acre-ft. The total water use for each temporary permit was assumed to be taken over the 3-month period allotted to temporary withdrawal permits. The groundwater-use analysis full period of record, 1967–2011, was divided into two sub-intervals because of varying water use, 1970–80 and 1981–2011. Groundwater use in Reach I and Reach II was substantially greater from 1970 to 1980 compared to the rest of the period, and the sub-period 1981–2011 was used because this period includes recent population growth and modern irrigation methods. The total mean annual groundwater use in Reach I was 15,309 acre-feet per year (acre-ft/yr) during 1967–2011; 20,724 acre-ft/yr during 1970–80, and 13,739 acre-ft/yr during 1981–2011. Total mean annual groundwater use in Reach II was similar but slightly less than in Reach I, with 14,098 acre-ft/yr during 1967–2011; 19,963 acre-ft/yr during 1970–80; and 12,285 acre-ft/yr during 1981–2011. Irrigation composed 72 percent of groundwater use in Reach I and 48 percent of groundwater use in Reach II during the 1967–2011 period. Public water supply was a much smaller proportion of total groundwater use in Reach I (15 percent) than in Reach II (39 percent). The proportion of groundwater use for power was 10 percent in Reach I and 5.2 percent in Reach II. All other water-use categories in Reach I only composed 2.2 percent of groundwater use in Reach I. In Reach II, industrial, mining, and commercial categories combined accounted for 4.4 percent of groundwater use; recreation, fish, and wildlife groundwater use accounted for 2.3 percent; and nonirrigated agriculture accounted for 1.5 percent of groundwater use. Permian-age bedrock underlies the Beaver-North Canadian River alluvial aquifer. In the east, the Dog Creek Shale, the Duncan Sandstone, and the Blaine and Chickasha Formations, none of which are notable sources of groundwater in the study area, underlie the Beaver-North Canadian River alluvial aquifer. In the northwestern part of Reach I, bedrock is composed of the Rush Springs and Marlow Formations, which are productive aquifers in some areas. The Cloud Chief Formation is not a source of groundwater. One hydrogeological unit was delineated in the Beaver-North Canadian River alluvial aquifer, composed of the terrace deposits and alluvium, with limited flow between this unit and bedrock units. Groundwater in this aquifer generally flows from northwest to southeast and across the aquifer toward the Beaver and North Canadian Rivers. Groundwater recharge from precipitation was estimated for the entire Beaver-North Canadian River alluvial aquifer and then itemized for both reaches by using a soil-water-balance (SWB) model. At two locations in Reach I, a water-table fluctuation method was used to estimate local recharge. Total mean annual groundwater recharge from the soil-water-balance method was estimated to be approximately 136,400 acre-ft in Reach I and 82,400 acre-ft in Reach II; the mean annual recharge for both reaches combined was approximately 218,800 acre-ft. Two sites in Reach I located at observation wells with continuous water-level measurements and nearby streamflow-gaging stations with precipitation gages were used to estimate the percentage of precipitation that becomes groundwater recharge. The Woodward site was located at observation well OW-4 near the Woodward, Okla. (07237500), streamflow-gaging station. Total precipitation and recharge for the Woodward and Seiling sites were calculated for the water year 2013. The Woodward site had a total of 14.18 in. of precipitation and 6.3 in. of recharge was calculated, equaling 44 percent of precipitation. The mean percentage of precipitation that was estimated to become recharge in the SWB model for the period 1980–2011 at that location was 9.2 percent, although adjacent SWB-model cells were as high as 20 percent of precipitation. The Seiling site had a total of 26.84 in. of precipitation during the water year 2013, and a total of 6.9 in. of recharge was estimated, equaling 25.9 percent of precipitation. At the Seiling site, the mean percentage of precipitation that became recharge in the SWB model for the period 1980–2011 was 23.0 percent. The principal inflow to the Beaver-North Canadian River alluvial aquifer was estimated to be surface recharge from precipitation, and plant evapotranspiration was estimated to be the greatest discharge, followed by stream and lake base flow, groundwater pumping, and flow to seeps and springs along the eastern margin of the aquifer. Reach I also included inflow from the High Plains aquifer as lateral inflow of groundwater, though this flow was estimated to be a very minor component of the total water budget. Most of the Beaver and North Canadian Rivers were determined to be gaining streamflow from groundwater, but several reaches in Reach I upstream from Wolf Creek were determined to be losing streamflow through infiltration to the aquifer. Aquifer hydrogeologic characteristics were estimated from borehole lithologic logs, well-construction information, and published aquifer tests and during numerical model calibration. The maximum saturated aquifer thickness in Reach I was estimated to be 308 ft, and the mean thickness was estimated to be 36 ft. The maximum saturated thickness in Reach II was estimated to be 86 ft, and the mean thickness was estimated to be 29 ft. Mean hydraulic conductivity of Reach I was estimated to be 70 feet per day (ft/d) with a range of 7–279 ft/d. Mean hydraulic conductivity in Reach II was estimated to be 92 ft/d with a range of 4–279 ft/d. Both reach models were calibrated manually by using trial-and-error adjustment of recharge, hydraulic conductivity, specific yield, and conductance of boundary conditions. The Reach I model used 28 head observations during the steady-state period of 1980 and 487 head observations during the transient period of 1981–2011. The root-mean-square error of head residuals (observed minus simulated head) was 3.86 ft, and 83 percent of head residuals were between -5 and 5 ft. The Reach II model was calibrated to 75 steady-state head observations and 134 head observations during the transient period. The root-mean-square error of head residuals for that reach was 3.58 ft, and similar to Reach I, 85 percent of residuals were between -5 and 5 ft. Several analyses were performed by using the numeric groundwater-flow models as predictive tools, including estimating the EPS pumping rate for both reaches. The EPS is defined by the Oklahoma Water Resources Board as an annual per-acre groundwater-pumping rate that will reduce saturated thickness in half of the aquifer to 5 ft or less over a period of 20 years; additional estimates were made for periods of 40 and 50 years. Other analyses included using models to estimate the effects of groundwater pumping and a prolonged drought on groundwater in storage and streamflow and lake storage of water. The EPS pumping rate was found to be approximately 0.57 acre-feet per acre per year ([acre-ft/acre]/yr) in Reach I and 0.73 (acre-ft/acre)/yr in Reach II for a 20-year period. For a 40-year period, the annual EPS pumping rate was determined to be 0.54 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. For a 50-year period, the EPS pumping rate was determined to be 0.53 (acre-ft/acre)/yr in Reach I and 0.61 (acre-ft/acre)/yr in Reach II. Groundwater pumping at the 2011 rate for 50 years resulted in a 3.6-percent decrease in the amount of water in groundwater storage in Reach I and a decrease of 2.5 percent in the amount of groundwater in storage in Reach II. A cumulative 32-percent increase in pumping greater than the 2011 rate over a period of 50 years caused a decrease in groundwater storage of 4.0 percent in Reach I and 3.3 percent in Reach II. A hypothetical severe drought was simulated by using aquifer recharge flow rates during the drought year of 2011 for a period of 10 years. All other flows including evapotranspiration and groundwater pumping were set at estimated 2011 rates. The hypothetical drought caused a decrease in water in aquifer storage by about 7 percent in Reach I and 7 percent in Reach II. Another analysis of the effects of hypothetical drought estimated the effects of drought on streamflow and lake storage. The hypothetical drought was simulated by decreasing recharge by 75 percent for a selected 10-year period (1994–2004) during the 1980–2011 simulation. In Reach I, the amounts of water stored in Canton Lake and streamflow at the Seiling, Okla., streamflow-gaging station were analyzed. Streamflow at the Seiling station decreased by a mean of 75 percent and was still diminished by 10 percent after 2011. In Reach II, the effect of drought on the streamflow at the Yukon, Okla., streamflow-gaging station was examined. The greatest mean streamflow decrease was approximately 60 percent during the simulated drought, and after 2011, the mean decrease in streamflow was still about 5 percent. Canton Lake storage decreased by as much as 83 percent during the simulated drought and did not recover by 2011.

Oklahoma↗

USGS2002

Putting aquifers into atmospheric simulation models: An example from the Mill Creek Watershed, Northeastern Kansas

Aquifer-atmosphere interactions can be important in regions where the water table is shallow (<2 m). A shallow water table provides moisture for the soil and vegetation and thus acts as a source term for evapotranspiration to the atmosphere. A coupled aquifer-land surface-atmosphere model has been developed to study aquifer-atmosphere interactions in watersheds, on decadal timescales. A single column vertically discretized atmospheric model is linked to a distributed soil-vegetation-aquifer model. This physically based model was able to reproduce monthly and yearly trends in precipitation, stream discharge, and evapotranspiration, for a catchment in northeastern Kansas. However, the calculated soil moisture tended to drop to levels lower than were observed in drier years. The evapotranspiration varies spatially and seasonally and was highest in cells situated in topographic depressions where the water table is in the root zone. Annually, simulation results indicate that from 5-20% of groundwater supported evapotranspiration is drawn from the aquifer. The groundwater supported fraction of evapotranspiration is higher in drier years, when evapotranspiration exceeds precipitation. A long-term (40 year) simulation of extended drought conditions indicated that water table position is a function of groundwater hydrodynamics and cannot be predicted solely on the basis of topography. The response time of the aquifer to drought conditions was on the order of 200 years indicating that feedbacks between these two water reservoirs act on disparate time scales. With recent advances in the computational power of massively parallel supercomputers, it may soon become possible to incorporate physically based representations of aquifer hydrodynamics into general circulation models (GCM) land surface parameterization schemes. ?? 2002 Elsevier Science Ltd. All rights reserved.

Advances in Water Resources↗

USGS1980

Water-supply assessment of the Laramie-Fox Hills aquifer in parts of Adams, Boulder, Jefferson, and Weld Counties, Colorado

Groundwater in the Laramie-Fox Hills aquifer is a potential source of supplemental municipal water supplies for the communities of Erie, Lafayette, Louisville, and Superior in Colorado. The present water supplies for these communities are not always adequate to meet current demands. The U.S. Geological Survey made a water-supply assessment of the Laramie-Fox Hills aquifer for the U.S. Bureau of Reclamation, which is investigating and evaluating alternative sources of water for the communities. Recharge to the aquifer is mostly in the western and southwestern parts of the study area. Groundwater movement is generally from the southwest to northeast. Groundwater discharge in the study area is primarily by pumping wells. Since 1961, this pumping has caused water-level declines of about 250 to 300 feet from Broomfield to east of Erie, Colorado. Generally, water levels in other parts of the area have remained the same. The aggregate sand thickness determined from well logs ranges from 42 to 360 feet and the mean thickness is 229 feet. The volume of groundwater in storage in the study area is about 5 million acre-feet. Reported yields from 93 wells ranged from 1 to 90 gallons per minute and averaged 22 gallons per minute. Well yields tended to be larger in the areas where aggregate sand thickness is the greatest. The water generally changes from a sodium calcium bicarbonate type to a sodium calcium sulfate type as it moves through the aquifer away from the recharge areas. The maximum limit established by the U.S. Environmental Protection Agency for nitrite plus nitrate in public-water supplies was exceeded in water from three wells, the maximum limit for fluoride was exceeded in water from two wells, and the maximum limit for selenium was exceeded in water from three wells.

Colorado↗

USGS2017

The magnitude and origin of groundwater discharge to eastern U.S. and Gulf of Mexico coastal waters

Fresh groundwater discharge to coastal environments contributes to the physical and chemical conditions of coastal waters, but the role of coastal groundwater at regional to continental scales remains poorly defined due to diverse hydrologic conditions and the difficulty of tracking coastal groundwater flow paths through heterogeneous subsurface materials. We use three-dimensional groundwater flow models for the first time to calculate the magnitude and source areas of groundwater discharge from unconfined aquifers to coastal waterbodies along the entire eastern U.S. We find that 27.1 km 3 /yr (22.8–30.5 km 3 /yr) of groundwater directly enters eastern U.S. and Gulf of Mexico coastal waters. The contributing recharge areas comprised ~175,000 km 2 of U.S. land area, extending several kilometers inland. This result provides new information on the land area that can supply natural and anthropogenic constituents to coastal waters via groundwater discharge, thereby defining the subterranean domain potentially affecting coastal chemical budgets and ecosystem processes.

Gulf of Mexico↗

USGS2018

Rebuttal to “The case of the Biscayne Bay and aquifer near Miami, Florida: density-driven flow of seawater or gravitationally driven discharge of deep saline groundwater?” by Weyer (Environ Earth Sci 2018, 77:1–16)

A recent paper by Weyer (Environ Earth Sci 2018, 77:1–16) challenges the widely accepted interpretation of groundwater heads and salinities in the coastal Biscayne aquifer near Miami, Florida, USA. Weyer (2018) suggests that the body of saltwater that underlies fresh groundwater just inland of the coast is not a recirculating wedge of seawater, but results instead from upward migration of deep saline groundwater driven by regional flow. Perhaps more significantly, Weyer (2018) also asserts that established hydrologic theory is fundamentally incorrect with respect to buoyancy. Instead of acting along the direction of gravity (that is, vertically), Weyer (2018) claims, buoyancy acts instead along the direction of the pressure gradient. As a result, Weyer (2018) considers currently available density-dependent groundwater flow and transport modeling codes, and the analyses based on them, to be in error. In this rebuttal, we clarify the inaccuracies in the main points of Weyer’s (2018) paper. First, we explain that Weyer (2018) has misinterpreted observed equivalent freshwater heads in the Biscayne aquifer and that his alternative hypothesis concerning the source of the saltwater does not explain the observed salinities. Then, we review the established theory of buoyancy to identify the problem with Weyer’s (2018) alternative theory. Finally, we present theory and cite successful benchmark simulations to affirm the suitability of currently available codes for modeling density-dependent groundwater flow and transport.

Environmental Earth Sciences↗

USGS2024

Generalized potentiometric maps of the Fort Union, Hell Creek, and Fox Hills aquifers within the Standing Rock Reservation

Generalized potentiometric surfaces of the Fort Union, Hell Creek, and Fox Hills aquifers were constructed to assess the groundwater resources of the Standing Rock Reservation. Additionally, this information can provide water managers with tools and data to effectively manage water resources in the future. Previous studies that mapped the geology and hydrogeology of the area at differing scales were used to confirm in which aquifer the study wells were completed. Water-level data from wells are provided by the U.S. Geological Survey Groundwater System Inventory database, the South Dakota Department of Agriculture and Natural Resources, and the North Dakota Department of Water Resources. Hydrographs were constructed for five selected observation wells to evaluate historical water-level fluctuations and trends. Hydrographs for the Hell Creek aquifer showed a flat trend with a rise in 2020. Hydrographs for the deeper Fox Hills aquifer showed that water levels fluctuated in response to climatic conditions and demonstrated an increasing trend in water-level elevations starting in 2010. Hydrographs were not constructed for any wells completed in the Fort Union Formation because none of the wells had continuous long-term measurements. Generalized potentiometric surfaces, constructed from interpolating water-level elevations, gave insight into groundwater flow directions. Groundwater in the Fort Union aquifer likely flows radially outward from the northwest part of the study area to the northeast and south-southeast parts. Groundwater in the Hell Creek aquifer generally flows from higher elevations in the northwest towards lower areas, where surface-water tributaries have incised into the aquifer. Groundwater in the Fox Hills aquifer likely flows from higher elevations in the west, southwest, and central parts of the study area towards the valleys of the Grand River and Missouri River. Most wells used for constructing potentiometric maps had only one recorded water-level measurement from drillers at the time of well construction. These measurements are often subject to error because the well is still recovering and because of spatial limitations of data availability. Also, because single water-level measurements were recorded at different points in time, additional uncertainty is introduced by fluctuating climatic conditions effect on water levels. Potentiometric map interpretation limitations are a result of areas with sparse data. Limitations also arise from potentiometric surfaces generalizing a complex and dynamic hydrogeologic system; however, the generalized potentiometric surface maps can be used to assist water managers and can help prioritize locations for future monitoring in areas with high uncertainty from sparse existing data.

North Dakota, South Dakota↗

USGS2013

Development of a geodatabase for springs within and surrounding outcrops of the Trinity aquifer in northern Bexar County, Texas, 2010-11

The Trinity aquifer is an important source of groundwater in central Texas, including Bexar County, where population growth has resulted in an increased demand for water (Ashworth, 1983; Mace and others, 2000). Numerous springs issue from rock outcrops within and surrounding the Trinity aquifer in northern Bexar County (fig. 1). The effects of increased groundwater withdrawals from the Trinity aquifer on springflow in the area are not well documented, but because the total amount of water entering, leaving, and being stored in a groundwater system must be conserved, increased groundwater withdrawals will result in decreases in springflow (Alley and others, 1999). Documenting the location, discharge, and basic water-quality information of the springs in northern Bexar County can provide a baseline assessment for comparison to future conditions. Accordingly, the U.S. Geological Survey (USGS), in cooperation with the Trinity Glen Rose Groundwater Conservation District, the Edwards Aquifer Authority, and the San Antonio River Authority, developed a geodatabase populated with data associated with springs within and surrounding outcrops of the Trinity aquifer in northern Bexar County during 2010–11. A geodatabase provides a framework for organizing spatial and tabular data (such as the geographic location and water-quality characteristics, respectively) in a relational database environment, making it easier and more intuitive to evaluate changes over time. Data for 141 springs within and surrounding the Trinity aquifer outcrops in northern Bexar County were compiled from existing reports and databases. These data were augmented with selected data collected onsite, including the location, discharge, and water-quality characteristics of selected springs, and were entered into the geodatabase. The Trinity aquifer in central Texas is commonly divided into the upper, middle, and lower Trinity aquifers; all of the information that was compiled pertaining to the aquifer is for the upper and middle Trinity aquifers.

Texas↗

USGS1987

Effect of reduced industrial pumpage on the migration of dissolved nitrogen in an outwash aquifer at Olean, Cattaraugus County, New York

A quasi-three-dimensional digital groundwater flow model of a shallow outwash aquifer system at Olean, New York, was developed to study the effects of several pumping alternatives on groundwater flow and stream seepage. Nitrogen compounds have contaminated the aquifer in an industrial park in North Olean. Pumping from seven industrial production wells and a purge well has created a cone of depression within which all nitrogen compounds are contained, thus preventing their migration to nearby private, municipal, and industrial-supply wells. A simulated total shutdown of an industrial well field and the purge well indicates that groundwater flowing laterally southward to southwestward from the area of contamination would require about 5 yr to reach a municipal well field along the Allegheny River. Simulation of a partial shutdown with only three main production wells and the purge well or just the three production wells indicate that all nitrogen-bearing groundwater would be captured. Pumping from only the purge well would allow contaminated groundwater along the southwestern edge of the site to escape and, in 8 to 9 yr, reach the municipal well field. (USGS)

Water-Resources Investigations Report↗

USGS2026

Hydrogeology, groundwater salinity distributions, and assessment of the effect of oil-production activities on groundwater in the Midway Valley area, western Kern County, San Joaquin Valley, California

This study seeks to determine the effects of oil field produced water disposal operations and well mechanical integrity issues on groundwater quality in oil fields in the southwest San Joaquin Valley, California. Whereas previous studies used groundwater wells to study shallow aquifers outside the oil fields, this study demonstrates that future approaches may use oil well geophysical logs to map groundwater head gradients, create salinity profiles and document changes in salinity over time in oil field areas with sparse groundwater well data and at depths greater than 330 m. We also incorporate an analysis of well histories to determine potential effects of compromised wellbore seals on changes in aquifer quality that cannot be explained by water disposal practices. Water quality in the aquifers is naturally brackish across most of the area, with better quality groundwater occurring in the eastern part. Geophysical logs are used to determine salinity variations within aquifers including the depth at which TDS exceeds 10,000 mg/L. This depth ranges from 366 m in the northwest to approximately 1,500 m in the southeast. Oil well porosity logs are used to determine water table elevations. These logs indicate the water table slopes south-southeast, showing the predominant groundwater flow direction is from oil field disposal areas toward better quality groundwater east of the oil fields. Geophysical logs show formation resistivity near some disposal facilities has decreased over time, indicating the salinity of the aquifer has increased due to disposal of saline produced water in injection wells and ponds. Oil well history analysis suggests that increased salinity over time in water-saturated sand intervals >1.5 km from disposal facilities may be caused by mechanical failures and/or incomplete borehole seals in poorly constructed or abandoned wellbores prevalent throughout the study area—particularly wells drilled prior to 1930.

California↗

USGS2010

Groundwater as a nonpoint source of atrazine and deethylatrazine in a river during base flow conditions

Alluvial groundwater adjacent to the main stem river is the principal nonpoint source of atrazine and deethylatrazine in the Cedar River of Iowa after the river has been in base flow conditions for 5 days. Between two sites along a 116-km reach of the Cedar River, tributaries contributed about 25% of the increase in the atrazine and deethylatrazine load, whereas groundwater from the alluvial aquifer contributed at least 75% of the increase in load. Within the study area, tributaries aggregate almost all of the discharge from tile drains, and yet the tributaries still only contribute 25% of the increase in loads in the main stem river. At an unfamned study site adjacent to the Cedar River, the sources of atrazine and deethylatrazine in the alluvial groundwater are bank storage of river water and groundwater recharge from areas distant from the river. Atrazine and deethylatrazine associated with bank storage water will provide larger concentrations to the river during early base flow conditions. After the depletion of bank storage, stable and smaller concentrations of atrazine and deethylatrazine, originating from groundwater recharge, continue to be discharged from the alluvial aquifer to the river; thus these results indicate that alluvial aquifers are an important nonpoint source of atrazine and deethylatrazine in rivers during base flow.

Iowa↗