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 415 records · Page 23Linked to original sources

Estimation of aquifer scale proportion using equal area grids: assessment of regional scale groundwater quality

The proportion of an aquifer with constituent concentrations above a specified threshold (high concentrations) is taken as a nondimensional measure of regional scale water quality. If computed on the basis of area, it can be referred to as the aquifer scale proportion. A spatially unbiased estimate of aquifer scale proportion and a confidence interval for that estimate are obtained through the use of equal area grids and the binomial distribution. Traditionally, the confidence interval for a binomial proportion is computed using either the standard interval or the exact interval. Research from the statistics literature has shown that the standard interval should not be used and that the exact interval is overly conservative. On the basis of coverage probability and interval width, the Jeffreys interval is preferred. If more than one sample per cell is available, cell declustering is used to estimate the aquifer scale proportion, and Kish's design effect may be useful for estimating an effective number of samples. The binomial distribution is also used to quantify the adequacy of a grid with a given number of cells for identifying a small target, defined as a constituent that is present at high concentrations in a small proportion of the aquifer. Case studies illustrate a consistency between approaches that use one well per grid cell and many wells per cell. The methods presented in this paper provide a quantitative basis for designing a sampling program and for utilizing existing data.

Water Resources Research↗

Comparison of aquifer characterization approaches through steady state groundwater model validation: A controlled laboratory sandbox study

Groundwater modeling has become a vital component to water supply and contaminant transport investigations. An important component of groundwater modeling under steady state conditions is selecting a representative hydraulic conductivity (K) estimate or set of estimates which defines the K field of the studied region. Currently, there are a number of characterization approaches to obtain K at various scales and in varying degrees of detail, but there is a paucity of information in terms of which characterization approach best predicts flow through aquifers or drawdowns caused by some drawdown inducing events. The main objective of this paper is to assess K estimates obtained by various approaches by predicting drawdowns from independent cross-hole pumping tests and total flow rates through a synthetic heterogeneous aquifer from flow-through tests. Specifically, we (1) characterize a synthetic heterogeneous aquifer built in the sandbox through various techniques (permeameter analyses of core samples, single-hole, cross-hole, and flow-through testing), (2) obtain mean K fields through traditional analysis of test data by treating the medium to be homogeneous, (3) obtain heterogeneous K fields through kriging and steady state hydraulic tomography, and (4) conduct forward simulations of 16 independent pumping tests and six flowthrough tests using these homogeneous and heterogeneous K fields and comparing them to actual data. Results show that the mean K and heterogeneous K fields estimated through kriging of small-scale K data (core and single-hole tests) yield biased predictions of drawdowns and flow rates in this synthetic heterogeneous aquifer. In contrast, the heterogeneous K distribution or ?K tomogram? estimated via steady state hydraulic tomography yields excellent predictions of drawdowns of pumping tests not used in the construction of the tomogram and very good estimates of total flow rates from the flowthrough tests. These results suggest that steady state groundwater model validation is possible in this laboratory sandbox aquifer if the heterogeneous K distribution and forcing functions (boundary conditions and source/sink terms) are characterized sufficiently. ?? 2010 by the American Geophysical Union.

Water Resources Research↗

Simulated changes in potentiometric levels resulting from groundwater development for phosphate mines, west-central Florida

A digital model of two-dimensional groundwater flow was used to predict changes in the potentiometric surface of the Floridan aquifer resulting from groundwater development for proposed and existing phosphate mines during 1976-2000. The modeled area covers 15,379 km 2 in west-central Florida. In 1975, groundwater withdrawn from the Floridan aquifer for irrigation, phosphate mines, other industries and municipal supplies averaged about 28,500 1/s. Withdrawals for phosphate mines are expected to shift from Polk County to adjacent counties to the south and west, and to decline from about 7,620 1/s in 1975 to about 7,060 1/s in 2000. The model was calibrated under steady-state and transient conditions. Input parameters included aquifer transmissivity and storage coefficient; thickness, vertical hydraulic conductivity, and storage coefficient of the upper confining bed; altitudes of the water table and potentiometric surface; and groundwater withdrawals. Simulation of November 1976 to October 2000, using projected combined pumping rates for existing and proposed phosphate mines, resulted in a rise in the potentiometric surface of about 6 m in Polk County, and a decline of about 4 m in parts of Manatee and Hardee counties.

Florida↗

Status and understanding of groundwater quality in the northern Sierra Nevada foothills domestic-supply aquifer study units, 2015–17—California GAMA Priority Basin Project

Groundwater quality in the northern Sierra Nevada foothills region of California was investigated as part of California State Water Resources Control Board (SWRCB) Groundwater Ambient Monitoring Assessment Priority Basin Project (GAMA-PBP). The region was divided into two study units: the Yuba-Bear watersheds (YBW) study unit and the American-Cosumnes-Mokelumne watersheds (ACMW) study unit. The GAMA-PBP made a spatially unbiased assessment of aquifer systems used for domestic drinking-water supply in the study region, which are predominantly composed of fractured, hard-rock aquifers of varying lithology. These assessments characterized the quality of raw groundwater to evaluate ambient conditions in the domestic-supply aquifer and not the quality of treated drinking water. The study included three components: (1) a status assessment , which characterized the quality of groundwater resources used for domestic drinking-water supply in the YBW and ACMW study units; (2) an understanding assessment , which evaluated natural and anthropogenic explanatory factors that could potentially affect groundwater quality in the study region; and (3) a comparative assessment between the groundwater resources used for domestic and public drinking-water supply in the study region. The status assessment was based on data collected by the GAMA-PBP from 74 sites in the YBW study unit during 2015–16 and 67 sites in the ACMW study unit from 2016 to 2017. To contextualize water-quality results, concentrations of water-quality constituents in ambient groundwater were compared to regulatory and non-regulatory benchmarks typically used by the State of California and Federal agencies as health-based or aesthetic standards for public drinking water. The status assessment used a grid-based method to estimate proportions of groundwater resources with concentrations approaching or exceeding benchmark thresholds. This method provides spatially unbiased results and allows inter-comparability with similar groundwater-quality assessments. Inorganic constituents with health-based benchmarks were present at high relative concentration (RC), meaning they exceeded the benchmark threshold, in 5.4 and 10 percent of domestic-supply aquifer systems in the YBW and ACMW study units, respectively. Inorganic constituents with aesthetic-based benchmarks were detected at high-RCs in 20 and 28 percent of the YBW and ACMW study units, respectively. The inorganic constituents present at high RC were arsenic, barium, boron, molybdenum, strontium, nitrate, adjusted gross-alpha particle activity, chloride, total dissolved solids, specific conductance, iron, manganese, and hardness. Groundwater samples were tested for presence or absence of three microbial indicators (total coliform, Escherichia coli , and Enterococci ). At least one microbial indicator was present in 26 and 28 percent of the YBW and ACMW study units, respectively. At least one organic constituent was detected in 30 and 42 percent of the YBW and ACMW study units, respectively. Organic constituents were not present at high RC, but tetrachloroethene (PCE), trichloroethene (TCE), and toluene were detected in the YBW study unit at moderate RC (between the benchmark concentration and one-tenth of the benchmark concentration). Methyl tert -butyl ether (MTBE) and chloroform were present at low RC (less than one-tenth of the benchmark concentration) in the YBW and ACMW study units with detection frequencies greater than 10 percent. Perchlorate, a constituent of special interest, was detected in 31 and 41 percent of the YBW and ACMW study units, respectively, at either low or moderate RCs. Relations among select water-quality constituents and potential explanatory factors were evaluated using statistical and graphical approaches. Nitrate, microbial indicators, and perchlorate were all correlated to elevation-dependent variables relating to climate, land use, and recharge condition. Isotopic and dissolved noble-gas tracers indicated these water-quality constituents are associated with recharge conditions associated with irrigation during the summer dry-season, which is common in areas of rural-residential or agricultural land uses. Higher concentrations of iron and manganese were primarily associated with anoxic groundwater in aquifers of metasedimentary lithology. Increased hardness was primarily associated with anoxic groundwater in aquifers of mafic-ultramafic or metavolcanics lithologies at lower elevations in the study region in the Melones fault zone. Chloroform and MTBE were associated with shallow groundwater (wells depths less than 130 m) under oxic and anoxic redox conditions, respectively. The comparative assessment evaluated differences between the aquifer systems used for domestic- and public-supply in study region based on (1) well-construction characteristics, and (2) water quality. Analysis of over 60,000 well-completion reports in the study region showed that although domestic-supply wells span the deepest depth zones in regional aquifers, median depths for public-supply wells were significantly greater than those of domestic-supply wells in both study units. Water-quality data from more than 300 public-supply wells in the study region were assessed using a spatially weighted method for calculation aquifer-scale proportions and compared with the domestic-supply assessment results. Detections of inorganic constituents at high RC and detection frequencies for organic constituents were generally similar between the domestic- and public-supply aquifer systems in both study units, with a few notable exceptions in the ACMW study unit: nitrate was greater for the public- compared to domestic-supply aquifer system and both manganese, hardness, and MTBE were greater in the domestic- compared to public-supply aquifer system. These differences are likely related to contrasting land uses, aquifer lithologies, landscape positions, and depths characterizing domestic- and public-supply wells in the ACMW study unit. Overall, fewer samples from domestic-supply wells in the northern Sierra Nevada foothills exceeded health-based benchmarks compared to aesthetic-based benchmarks for groundwater quality. Exceedences of health-based benchmarks were primarily caused by nitrate and coliform bacteria, which were associated with recharge from diverted surface water used primarily for irrigation. Exceedences of aesthetic-based benchmarks were primarily caused by iron, managanese, and hardness, which were associated with geologic factors. Regional irrigation practices and aquifer lithology can affect groundwater quality in fractured-rock aquifers in the northern Sierra Nevada foothills used for domestic drinking-water supply.

California↗

Simulation of groundwater flow, effects of artificial recharge, and storage volume changes in the Equus Beds aquifer near the city of Wichita, Kansas well field, 1935–2008

The Equus Beds aquifer is a primary water-supply source for Wichita, Kansas and the surrounding area because of shallow depth to water, large saturated thickness, and generally good water quality. Substantial water-level declines in the Equus Beds aquifer have resulted from pumping groundwater for agricultural and municipal needs, as well as periodic drought conditions. In March 2006, the city of Wichita began construction of the Equus Beds Aquifer Storage and Recovery project to store and later recover groundwater, and to form a hydraulic barrier to the known chloride-brine plume near Burrton, Kansas. In October 2009, the U.S. Geological Survey, in cooperation with the city of Wichita, began a study to determine groundwater flow in the area of the Wichita well field, and chloride transport from the Arkansas River and Burrton oilfield to the Wichita well field. Groundwater flow was simulated for the Equus Beds aquifer using the three-dimensional finite-difference groundwater-flow model MODFLOW-2000. The model simulates steady-state and transient conditions. The groundwater-flow model was calibrated by adjusting model input data and model geometry until model results matched field observations within an acceptable level of accuracy. The root mean square (RMS) error for water-level observations for the steady-state calibration simulation is 9.82 feet. The ratio of the RMS error to the total head loss in the model area is 0.049 and the mean error for water-level observations is 3.86 feet. The difference between flow into the model and flow out of the model across all model boundaries is -0.08 percent of total flow for the steady-state calibration. The RMS error for water-level observations for the transient calibration simulation is 2.48 feet, the ratio of the RMS error to the total head loss in the model area is 0.0124, and the mean error for water-level observations is 0.03 feet. The RMS error calculated for observed and simulated base flow gains or losses for the Arkansas River for the transient simulation is 7,916,564 cubic feet per day (91.6 cubic feet per second) and the RMS error divided by (/) the total range in streamflow (7,916,564/37,461,669 cubic feet per day) is 22 percent. The RMS error calculated for observed and simulated streamflow gains or losses for the Little Arkansas River for the transient simulation is 5,610,089 cubic feet per day(64.9 cubic feet per second) and the RMS error divided by the total range in streamflow (5,612,918/41,791,091 cubic feet per day) is 13 percent. The mean error between observed and simulated base flow gains or losses was 29,999 cubic feet per day (0.34 cubic feet per second) for the Arkansas River and -1,369,250 cubic feet per day (-15.8 cubic feet per second) for the Little Arkansas River. Cumulative streamflow gain and loss observations are similar to the cumulative simulated equivalents. Average percent mass balance difference for individual stress periods ranged from -0.46 to 0.51 percent. The cumulative mass balance for the transient calibration was 0.01 percent. Composite scaled sensitivities indicate the simulations are most sensitive to parameters with a large areal distribution. For the steady-state calibration, these parameters include recharge, hydraulic conductivity, and vertical conductance. For the transient simulation, these parameters include evapotranspiration, recharge, and hydraulic conductivity. The ability of the calibrated model to account for the additional groundwater recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project was assessed by using the U.S. Geological Survey subregional water budget program ZONEBUDGET and comparing those results to metered recharge for 2007 and 2008 and previous estimates of artificial recharge. The change in storage between simulations is the volume of water that estimates the recharge credit for the aquifer storage and recovery system. The estimated increase in storage of 1,607 acre-ft in the basin storage area compared to metered recharge of 1,796 acre-ft indicates some loss of metered recharge. Increased storage outside of the basin storage area of 183 acre-ft accounts for all but 6 acre-ft or 0.33 percent of the total. Previously estimated recharge credits for 2007 and 2008 are 1,018 and 600 acre-ft, respectively, and a total estimated recharge credit of 1,618 acre-ft. Storage changes calculated for this study are 4.42 percent less for 2007 and 5.67 percent more for 2008 than previous estimates. Total storage change for 2007 and 2008 is 0.68 percent less than previous estimates. The small difference between the increase in storage from artificial recharge estimated with the groundwater-flow model and metered recharge indicates the groundwater model correctly accounts for the additional water recharged to the Equus Beds aquifer as part of the Aquifer Storage and Recovery project. Small percent differences between inflows and outflows for all stress periods and all index cells in the basin storage area, improved calibration compared to the previous model, and a reasonable match between simulated and measured long-term base flow indicates the groundwater model accurately simulates groundwater flow in the study area. The change in groundwater level through recent years compared to the August 1940 groundwater level map has been documented and used to assess the change of storage volume of the Equus Beds aquifer in and near the Wichita well field for three different areas. Two methods were used to estimate changes in storage from simulation results using simulated change in groundwater levels in layer 1 between stress periods, and using ZONEBUDGET to calculate the change in storage in the same way the effects of artificial recharge were estimated within the basin storage area. The three methods indicate similar trends although the magnitude of storage changes differ. Information about the change in storage in response to hydrologic stresses is important for managing groundwater resources in the study area. The comparison between the three methods indicates similar storage change trends are estimated and each could be used to determine relative increases or decreases in storage. Use of groundwater level changes that do not include storage changes that occur in confined or semi-confined parts of the aquifer will slightly underestimate storage changes; however, use of specific yield and groundwater level changes to estimate storage change in confined or semi-confined parts of the aquifer will overestimate storage changes. Using only changes in shallow groundwater levels would provide more accurate storage change estimates for the measured groundwater levels method. The value used for specific yield is also an important consideration when estimating storage. For the Equus Beds aquifer the reported specific yield ranges between 0.08 and 0.35 and the storage coefficient (for confined conditions) ranges between 0.0004 and 0.16. Considering the importance of the value of specific yield and storage coefficient to estimates of storage change over time, and the wide range and substantial overlap for the reported values for specific yield and storage coefficient in the study area, further information on the distribution of specific yield and storage coefficient within the Equus Beds aquifer in the study area would greatly enhance the accuracy of estimated storage changes using both simulated groundwater level, simulated groundwater budget, or measured groundwater level methods.

Kansas↗

Use of a digital model to evaluate hydrogeologic controls on groundwater flow in a fractured rock aquifer at Niagara Falls, New York, U.S.A.

The Hyde Park landfill is a 15-acre (6.1 ha) chemical waste disposal site located north of Niagara Falls, New York. Underlying the site in descending order are: (1) low-permeability glacial till and lacustrine deposits; (2) a moderately permeable fractured rock aquifer - the Lockport Dolomite; and (3) a low-permeability unit - the Rochester Shale. The site is bounded on three sides by groundwater drains; the Niagara River gorge, the Niagara Power Project canal, and the Niagara Power Project buried conduits. The mechanism by which groundwater moves through fractured rocks underlying a hazardous waste site was investigated using a digital simulation approach. Three hypotheses were tested related to flow in the fractured rocks underlying Hyde Park landfill. For this purpose we used a Galerkin finite-element approximation to solve a saturated-unsaturated flow equation. A primary focus was to investigate anisotropy in the Lockport Dolomite, that is the effectiveness of horizontal (bedding) joints vs. vertical joints as water-transmitting openings. Three hydrogeologic scenarios were set up — each with prescribed limits on the hydrologic parameters. Scenario 1 specified strongly anisotropic conditions in the Lockport Dolomite (horizontal hydraulic conductivity along bedding joints exceeds vertical conductivity by 2–3 orders of magnitude), uniform areal recharge (5 in. yr. −1 or 12.7 cm yr. −1 ) except at the landfill where there is no recharge, and no flow through the base of the Rochester Shale. Scenario 2 also specified strongly anisotropic conditions in the Lockport; however, areal recharge was 6 in. yr. −1 (15.2 cm yr. −1 ) except at the landfill where the recharge was 2 in. yr. −1 (5.1 cm yr. −1 ), and outflow from the Rochester occurred. Scenario 3 specified isotropic conditions (that is, permeability along horizontal and vertical joints is the same in the Lockport Dolomite), recharge rates were the same as in scenario 2 and outflow through Rochester occurred. Scenario 2 provided the closest agreement between the simulated and measured heads while scenario 3 provided the poorest agreement. Among the three scenarios tested, scenario 2 (with strongly anisotropic conditions in the Lockport Dolomite with added recharge through the landfill cap and limited flow through the Rocherster Shale) is considered the most realistic hydrogeologic model. Based on simulation with the hydrogeologic parameters of scenario 2, groundwater flow near the Hyde Park site can be summarized as follows: 1. (1) Specific discharge (Darcy velocity) ranges from ≈0.01 to 0.1 ft. day −1 (0.003 to 0.03 m day −1 ) in the upper unit of the Lockport Dolomite to slightly more than 0.0001 ft. day −1 (0.00003 m day −1 ) in the Rochester Shale. Real velocities are highest in the upper unit of the Lockport, ranging from ≈1 to 5 ft. day −1 (0.3 to 1.5 m day −1 ) if the average effective porosity is assumed to be 0.02. 2. (2) A groundwater divide exists east of the landfill, indicating that all groundwater originating near or flowing beneath the landfill will flow toward and discharge in the gorge. 3. (3) Highest flow velocities (and presumably greatest potential for transporting chemical contaminants) occur in the upper unit of the Lockport and part of the lower unit of the Lockport Dolomite between the landfill and the gorge. The average time required for groundwater to move from the landfill to the discharge points at the gorge along selected flow paths in the Lockport Dolomite is estimated to be 5-6 yr.

New York↗

Geochemistry of groundwater in the eastern Snake River Plain aquifer, Idaho National Laboratory and vicinity, eastern Idaho

Nuclear research activities at the U.S. Department of Energy (DOE) Idaho National Laboratory (INL) in eastern Idaho produced radiochemical and chemical wastes that were discharged to the subsurface, resulting in detectable concentrations of some waste constituents in the eastern Snake River Plain (ESRP) aquifer. These waste constituents may pose risks to the water quality of the aquifer. In order to understand these risks to water quality the U.S. Geological Survey, in cooperation with the DOE, conducted a study of groundwater geochemistry to improve the understanding of hydrologic and chemical processes in the ESRP aquifer at and near the INL and to understand how these processes affect waste constituents in the aquifer. Geochemistry data were used to identify sources of recharge, mixing of water, and directions of groundwater flow in the ESRP aquifer at the INL. The geochemistry data were analyzed from 167 sample sites at and near the INL. The sites included 150 groundwater, 13 surface-water, and 4 geothermal-water sites. The data were collected between 1952 and 2012, although most data collected at the INL were collected from 1989 to 1996. Water samples were analyzed for all or most of the following: field parameters, dissolved gases, major ions, dissolved metals, isotope ratios, and environmental tracers. Sources of recharge identified at the INL were regional groundwater, groundwater from the Little Lost River (LLR) and Birch Creek (BC) valleys, groundwater from the Lost River Range, geothermal water, and surface water from the Big Lost River (BLR), LLR, and BC. Recharge from the BLR that may have occurred during the last glacial epoch, or paleorecharge, may be present at several wells in the southwestern part of the INL. Mixing of water at the INL primarily included mixing of surface water with groundwater from the tributary valleys and mixing of geothermal water with regional groundwater. Additionally, a zone of mixing between tributary valley water and regional groundwater, trending southwesterly, extended from near the northeastern boundary of the INL to the southern boundary of the INL. Groundwater flow directions for regional groundwater were southwesterly, and flow directions for tributary groundwater were southeasterly upon entering the ESRP, but eventually began to flow southwesterly in a direction parallel with regional groundwater. Several discrepancies were identified from comparison of sources of recharge determined from geochemistry data and backward particle tracking with a groundwater-flow model. Some discrepancies observed in the particle tracking results included representation of recharge from BC near the north INL boundary, groundwater from the BC valley not extending far enough south, regional groundwater that extends too far west in the southern part of the INL, and no representation of recharge from geothermal water in model layer 1 or recharge from the BLR in the southwestern part of the INL.

Idaho↗

Temporal and spatial variability of water quality in the San Antonio segment of the Edwards aquifer recharge zone, Texas, with an emphasis on periods of groundwater recharge, September 2017–July 2019

Ongoing urbanization on the Edwards aquifer recharge zone in the greater San Antonio area raises concern about the potential adverse effects on the public water supply from development. To address this concern, the U.S. Geological Survey, in cooperation with the City of San Antonio, studied patterns of temporal and spatial changes in water quality at selected surface-water and groundwater sites in the Edwards aquifer recharge zone, with an emphasis on changes during periods of groundwater recharge. Water-quality characteristics were continuously monitored and discrete water samples were collected at two sets of paired surface-water (stream) and groundwater (well) sites during a 2-year period (2017–19) that included relatively dry conditions and a large recharge event in September 2018 when as much as 16 inches of rain fell in parts of the study area. Continuous monitoring of water-level altitude, specific conductance, and concentrations of nitrate in two wells completed in the Edwards aquifer provided high-resolution data showing detailed changes in water quality across a broad range of hydrologic conditions. Water levels in the wells responded rapidly (within hours to days) to recharge from both small and large rainfall and runoff events; changes in groundwater quality as a consequence of the influx of surface-derived recharge were indicated by changes in values of the monitored characteristics. A broad range in measured values of the stable isotopes of water expressed as delta deuterium and delta oxygen-18 in the water samples collected from two streams (Salado and West Elm Creeks), in comparison to the tight clustering of the values of these isotopes in groundwater samples, indicates that source waters (surface waters) of widely varying chemical characteristics become homogenized within the aquifer system. Concentrations of major ions, trace ions, and nutrient concentrations in stormwater runoff indicate a combination of land-derived and rainfall-derived constituents. The distribution of concentrations of nitrogen species (nitrite, nitrate, and nitrogen in ammonia) among sampling sites transitions from a more variable distribution in stormwater runoff to a more uniform distribution in groundwater in which the dominant form is nitrate. Differences in nitrate isotopic composition and concentration in groundwater across the study area are likely controlled by the relative contributions of natural and anthropogenic nitrogen (with the anthropogenic nitrogen component including a wastewater source) and by the process of nitrification. Among all measured constituents, pesticides detected in discrete stormwater-runoff samples provided the clearest indication that urbanization was adversely affecting water quality; specifically, the more urbanized surface-water site had a greater number of detections and greater variety of detected pesticides. Though temporal variability in the numbers and types of pesticides was evident, the overall proportion of pesticides was dominated by triazine herbicides including atrazine, atrazine degradates, and simazine. The observed hydrologic responses to rainfall and corresponding changes in water quality in wells are thought to result from the direct hydrologic connectivity of surface water and unconfined groundwater; however, patterns of groundwater-quality change indicate mixing from multiple sources such as ambient groundwater, recent surface-derived recharge, and possibly inflow from other aquifers. Therefore, understanding the connection between urbanization and groundwater quality cannot be inferred from the input of stormwater runoff alone as changes related to local and regional hydrologic conditions also need to be considered. It should be noted that a single study comparing the results from two site pairs is not able to support definitive conclusions about the full effect of urbanization on surface water/groundwater quality; however, this study does provide useful insights about the spatial and temporal variability of both stormwater runoff and unconfined groundwater that are consistent with expectations based on the current conceptual model that depicts the Edwards aquifer surface-water/groundwater system as a single water resource.

Texas↗

Drain tiles and groundwater resources: Understanding the relations

Executive Summary Drainage for agricultural production over the past 150 years has been an integral component of human-driven change to Minnesota’s rural landscapes. Benefits of drainage Historically, poorly drained soils across much of the State would often remain saturated or flooded after spring snowmelt, preventing timely farm operations such as tilling and planting crops (Arneman, 1963). Installation of agricultural drainage, both surface ditches and subsurface drainage, accelerated transport of water off farm fields and imparted producers higher crop yields (Beauchamp, 1987; Stoner and others, 1993). Agricultural drainage offered many other benefits such as preventing crop drown out, aerating the soil profile for improved plant growth, limiting surface runoff and soil erosion, and allowing farmers better access to croplands (Fausey and others, 1987). Without agricultural drainage on much of Minnesota’s croplands, it would have been difficult to realize high enough crop yields to remain economically viable. Environmental concerns While drainage of Minnesota croplands provided the benefits mentioned above, several environmental concerns result. These include wetland loss, degradation of downstream water quality, and reduced [potential for] recharge. Early agricultural drainage efforts (pre-20th century) led to the disappearance of much of Minnesota’s natural wetlands. Increased focus on preventing or mitigating wetland loss over the last 50 years has helped curtail further losses, even as agricultural drainage proceeds. Prior to establishment of Minnesota statehood, wetlands accounted for more than 10 million acres in Minnesota, including prairie wetlands, peatlands, and forest wetlands that comprised approximately 19 percent of the total land area (Palmer, 1915; King, 1980). In 2018, only half of Minnesota’s pre-settlement wetlands remain, mostly in parts of the State that have not experienced widespread drainage, such as northern Minnesota. Water-quality monitoring has shown that agricultural drainage, in particular the practice of subsurface drainage, provides a direct flow path for nutrient (nitrogen and soluble phosphorus) losses to surface water resources. The negative consequences of agricultural drainage on surface water quality are well documented (for example, Dinnes and others, 2002; Kladivko and others, 2004; Richards and others, 2008; Rozemeijer and others, 2010; Schottler and others, 2013). Agricultural basins with a high percentage of agricultural drainage have been implicated as part of the cause of the Gulf of Mexico hypoxia zone due to excessive nitrogen export (Goolsby and Battaglin, 2001; Randall and Mulla, 2001). The connection of hydrological effects of agricultural subsurface drainage on groundwater recharge and aquifers, on the other hand, has not been well-established. Agricultural subsurface drainage intercepts infiltrating water below croplands and directly discharges the water to nearby surface waters. However, the size of the water balance shift from drained water that would have evapotranspired or run off the land to drained water that would recharge underlying aquifers has been poorly characterized (Schuh, 2008). Drain Tiles and Groundwater Given the poor accounting of subsurface drainage effects on groundwater resources, the Minnesota Ground Water Association (MGWA) deemed it imperative that we document these effects so that groundwater resources in agricultural regions with substantial drainage can be effectively managed. This white paper documents the relations of drain tiles and groundwater resources and discusses the historical significance of agricultural drainage practices, the recognized positive benefits and potential negative consequences of agricultural drainage practices, and the gaps in understanding of the connections between agricultural drainage and groundwater resources. The major messages emerged from the findings of this white paper are: Complex history. Minnesota has a long history of agricultural drainage, spanning over 150 years. Agricultural drainage, and the eventual widespread usage of subsurface drainage, can be separated into at least four distinct periods of time: (a) early drainage to get water off the land, pre-20th century; (b) the boom and bust era (1900-1945); (c) postwar resurgence of subsurface drainage and early conservation efforts (1945-1960); and, (d) emergence of the environmental movement (1960 to present). The State’s regulatory framework that both allowed for drainage and controlled its usage during these periods is complicated and has been governed by a patchwork of both State and Federal statutes. Drainage Provinces. This white paper advances the concept of tile drainage provinces to aid in the discussion of regional differences in subsurface drainage and its overall effect on groundwater resources. Built upon the concept of groundwater provinces (Minnesota Department of Natural Resources, 2001), three distinct tile drainage provinces were conceptualized: (1) the Southeastern Province; (2) the SouthCentral Province; and, (3) the Western Province. The distinct geology and the soils that developed in these regions have implications for each region’s subsurface drainage density and the potential implications for groundwater. Knowledge gaps. Several critical knowledge gaps are identified in this paper, creating opportunities for further research to improve our understanding for better managed water resources: Extent of drainage is unknown. Direct estimates of the extent of subsurface drainage do not exist in Minnesota. However, several indirect methods have been utilized to estimate subsurface drainage, from the field-scale to countylevel through the use of geographic information system (GIS) analysis and aerial photography. Based on a 2012 U.S. Geological Survey estimate of subsurface drainage extent (Nakagaki and Wieczorek, 2016), about 21% of the land area in Minnesota has some density of subsurface drainage. Effect of drainage on underlying aquifers is unknown. A basic understanding of unconfined and confined aquifers and their recharge is necessary to connect any hydrological effects from agricultural drainage to groundwater. The basic goal of subsurface drainage to efficiently drain saturated soils clearly alters the water balance in croplands. However, its overall effect on groundwater resources has been poorly characterized, and is in large part determined by the geology below drained areas and the arrangement of underlying aquifers. Water balance shifts. An improved understanding of historical water balance shifts from pre- to post-drainage periods is necessary to understand long-term implications on net groundwater recharge. Also, more direct field-scale studies and indirect modeling studies are needed to characterize water budgets for fields with subsurface drainage.

Minnesota↗

Construction of a groundwater-flow model for the Big Sioux Aquifer using airborne electromagnetic methods, Sioux Falls, South Dakota

The city of Sioux Falls is the fastest growing community in South Dakota. In response to this continued growth and planning for future development, Sioux Falls requires a sustainable supply of municipal water. Planning and managing sustainable groundwater supplies requires a thorough understanding of local groundwater resources. The Big Sioux aquifer consists of glacial outwash sands and gravels and is hydraulically connected to the Big Sioux River, which provided about 90 percent of the city’s source-water production in 2015. Managing sustainable groundwater supplies also requires an understanding of groundwater availability. An effective mechanism to inform water management decisions is the development and utilization of a groundwater-flow model. A groundwater-flow model provides a quantitative framework for synthesizing field information and conceptualizing hydrogeologic processes. These groundwater-flow models can support decision making processes by mapping and characterizing the aquifer. Accordingly, the city of Sioux Falls partnered with the U.S. Geological Survey to construct a groundwater-flow model. Model inputs will include data from advanced geophysical techniques, specifically airborne electromagnetic methods.

South Dakota↗

Hydrogeology and groundwater quality in the Snake River alluvial aquifer at Jackson Hole Airport, Wyoming, 2011–20

The Snake River alluvial aquifer underlying the Jackson Hole Airport (JHA) in northwest Wyoming is an important source of water used for domestic, commercial, and irrigation purposes by the airport and nearby residents. The U.S. Geological Survey, in response to previously identified water-quality concerns in the area, monitored and evaluated changes in hydrogeologic characteristics and groundwater-quality conditions of the alluvial aquifer during 2011–20. During that period, the Jackson Hole Airport made several changes that potentially improved water quality at and downgradient from the airport. Well, water level, and hydrogeologic data were collected from the alluvial aquifer to identify hydrogeologic characteristic and groundwater quality changes. Additionally, results of statistical tests were applied to water-quality results to evaluate trends in selected physical properties and constituent concentrations with time. The trends of those data show that water quality did improve overall during the study period compared to previously collected data. Presumably, these trends are in response to the changes in the aircraft deicing/anti-icing fluid (ADAF) formulation used by the JHA, the many JHA infrastructure improvements made during 2011–20, the degradation of existing ADAFs in subsurface soils and groundwater, or some combination of these possibilities.

Wyoming↗

Hydrogeology, water budget, and simulated groundwater availability in the Salt Fork Arkansas River and Chikaskia River alluvial aquifers, northern Oklahoma, 1980–2020

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s aquifers to determine the maximum annual yield for each groundwater basin. The U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, conducted an updated hydrologic investigation of the Salt Fork Arkansas River and Chikaskia River alluvial aquifers in northern Oklahoma for the study period spanning 1980–2020 and evaluated the simulated effects of potential groundwater withdrawals on groundwater flow and availability in the Salt Fork Arkansas River alluvial aquifer. A hydrogeologic framework and conceptual model were developed to guide the development of a numerical model. Three groundwater-availability scenarios were evaluated by using the calibrated numerical model, which was focused on the Salt Fork Arkansas River alluvial aquifer. These scenarios were used to (1) estimate equal-proportionate-share groundwater withdrawal rates, (2) quantify the potential effects of projected well withdrawals on groundwater storage over a 50-year period, and (3) simulate the potential effects of a hypothetical 10-year drought. The 20-, 40-, and 50-year equal-proportionate-share groundwater withdrawal rates for the Salt Fork Arkansas River alluvial aquifer under normal recharge conditions were about 0.63, 0.58, and 0.57 acre-foot per acre per year, respectively. Projected 50-year groundwater withdrawal scenarios were used to simulate the effects of modified well withdrawal rates. Because well withdrawals were less than 2 percent of the calibrated numerical-model water budget, changes to the well groundwater withdrawal rates had little effect on simulated Salt Fork Arkansas River base flows and groundwater storage in the Salt Fork Arkansas River alluvial aquifer. A hypothetical 10-year drought scenario was used to simulate the potential effects of a prolonged period of reduced recharge on groundwater storage. Groundwater storage at the end of the hypothetical drought period was 14.5 percent less than the groundwater storage of the calibrated numerical model without the simulated drought.

Oklahoma↗

Evaluation of the groundwater-flow model for the Ohio River alluvial aquifer near Carrollton, Kentucky, updated to conditions in September 2010

The Ohio River alluvial aquifer near Carrollton, Ky., is an important water resource for the cities of Carrollton and Ghent, as well as for several industries in the area. The groundwater of the aquifer is the primary source of drinking water in the region and a highly valued natural resource that attracts various water-dependent industries because of its quantity and quality. This report evaluates the performance of a numerical model of the groundwater-flow system in the Ohio River alluvial aquifer near Carrollton, Ky., published by the U.S. Geological Survey in 1999. The original model simulated conditions in November 1995 and was updated to simulate groundwater conditions estimated for September 2010. The files from the calibrated steady-state model of November 1995 conditions were imported into MODFLOW-2005 to update the model to conditions in September 2010. The model input files modified as part of this update were the well and recharge files. The design of the updated model and other input files are the same as the original model. The ability of the updated model to match hydrologic conditions for September 2010 was evaluated by comparing water levels measured in wells to those computed by the model. Water-level measurements were available for 48 wells in September 2010. Overall, the updated model underestimated the water levels at 36 of the 48 measured wells. The average difference between measured water levels and model-computed water levels was 3.4 feet and the maximum difference was 10.9 feet. The root-mean-square error of the simulation was 4.45 for all 48 measured water levels. The updated steady-state model could be improved by introducing more accurate and site-specific estimates of selected field parameters, refined model geometry, and additional numerical methods. Collection of field data to better estimate hydraulic parameters, together with continued review of available data and information from area well operators, could provide the model with revised estimates of conductance values for the riverbed and valley wall, hydraulic conductivities for the model layer, and target water levels for future simulations. Additional model layers, a redesigned model grid, and revised boundary conditions could provide a better framework for more accurate simulations. Additional numerical methods would identify possible parameter estimates and determine parameter sensitivities.

Kentucky↗

Can arsenic occurrence rate in bedrock aquifers be predicted?

A high percentage (31%) of groundwater samples from bedrock aquifers in the greater Augusta area, Maine was found to contain greater than 10 μg L –1 of arsenic. Elevated arsenic concentrations are associated with bedrock geology, and more frequently observed in samples with high pH, low dissolved oxygen, and low nitrate. These associations were quantitatively compared by statistical analysis. Stepwise logistic regression models using bedrock geology and/or water chemistry parameters are developed and tested with external data sets to explore the feasibility of predicting groundwater arsenic occurrence rates (the percentages of arsenic concentrations higher than 10 μg L –1 ) in bedrock aquifers. Despite the under-prediction of high arsenic occurrence rates, models including groundwater geochemistry parameters predict arsenic occurrence rates better than those with bedrock geology only. Such simple models with very few parameters can be applied to obtain a preliminary arsenic risk assessment in bedrock aquifers at local to intermediate scales at other localities with similar geology.

Environmental Science & Technology↗

Groundwater quality of the Gulf Coast aquifer system, Houston, Texas, 2010

During March–December 2010, the U.S. Geological Survey, in cooperation with the city of Houston, collected source-water samples from 60 municipal supply wells in the Houston area. These data were collected as part of an ongoing study to determine concentrations, spatial extent, and associated geochemical conditions that might be conducive for mobility and transport of selected naturally occurring contaminants (selected trace elements and radionuclides) in the Gulf Coast aquifer system in the Houston area. In the summers of 2007 and 2008, a reconnaissance-level survey of these constituents in untreated water from 28 municipal supply wells was completed in the Houston area. Included in this report are the complete analytical results for 47 of the 60 samples collected in 2010—those results which were received from the laboratories and reviewed by the authors as of December 31, 2010. All of the wells sampled were screened in the Gulf Coast aquifer system; 22 were screened entirely in the Evangeline aquifer, and the remaining 25 wells contained screened intervals that intersected both Evangeline and Chicot aquifers. The data documented in this report were collected as part of an ongoing study to characterize source-water-quality conditions in untreated groundwater prior to drinking-water treatment. An evaluation of contaminant occurrence in source water provides background information regarding the presence of a contaminant in the environment. Because source-water samples were collected prior to any treatment or blending that potentially could alter contaminant concentrations, the water-quality results documented by this report represent the quality of the source water, not the quality of finished drinking water provided to the public. Samples were analyzed for major ions (calcium, magnesium, potassium, sodium, bromide, chloride, fluoride, silica, and sulfate), residue on evaporation (dissolved solids), trace elements (arsenic, barium, boron, chromium, iron, lithium, manganese, molybdenum, selenium, strontium, and vanadium), and selected radionuclides (gross alpha- and beta-particle activity [at 72 hours and 30 days], carbon-14, radium-226, radon-222, and uranium). Field measurements were made of selected physicochemical (relating to both physical and chemical) properties (oxidation-reduction potential, turbidity, dissolved-oxygen concentration, pH, specific conductance, water temperature, and alkalinity) and unfiltered sulfides. Similar to the results from the reconnaissance survey, physicochemical properties, major ions, and trace elements varied considerably. The ranges of selected physicochemical properties were as follows: oxidation-reduction potential ranged from -173 to 466 millivolts, dissolved oxygen ranged from less than 0.1 to 4.4 milligrams per liter, pH ranged from 7.2 to 7.8, specific conductance ranged from 439 to 724 microsiemens per centimeter at 25 degrees Celsius, and alkalinity ranged from 159 to 276 milligrams per liter as calcium carbonate. The largest ranges in concentration for filtered major ion constituents were obtained for cations sodium and calcium and for anions chloride and sulfate. Arsenic concentrations measured in samples from the 47 wells ranged from 1.6 to 23.5 micrograms per liter. The maximum concentration of arsenic (23.5 micrograms per liter) was measured in the source-water sample from well LJ-65-12-328. Quantifiable concentrations of barium, boron, lithium, molybdenum, and strontium were measured in all 47 filtered, source-water samples. Quantifiable concentrations of manganese were measured in 46 source-water samples, and an estimated concentration of manganese was measured in 1 sample. Chromium, iron, selenium, and vanadium were detected in 24 or more of the 47 source-water samples. Gross alpha-particle activities and beta-particle activities for all 47 samples were analyzed at 72 hours after sample collection and again at 30 days after sample collection, allowing for the measurement of the activity of short-lived isotopes. Gross alpha-particle activities reported in this report were not adjusted for activity contributions by radon or uranium and, therefore, are conservatively high estimates if compared to the U.S. Environmental Protection Agency National Primary Drinking Water Regulation for adjusted gross alpha-particle activity. The gross alpha-particle activities at 30 days in the samples ranged from R0.60 to 25.5 picocuries per liter and at 72 hours ranged from 2.58 to 39.7 picocuries per liter, and the "R" preceding the value of 0.60 picocuries per liter refers to a nondetected result less than the sample-specific critical level. Gross beta-particle activities measured at 30 days ranged from 1.17 to 14.4 picocuries per liter and at 72 hours ranged from 1.97 to 4.4 picocuries per liter. Filtered uranium was detected in quantifiable amounts in all of the 47 wells sampled. The uranium concentrations ranged from 0.03 to 42.7 micrograms per liter. One sample was analyzed for carbon-14, and the amount of modern atmospheric carbon was reported as 0.2 percent. Six source-water samples collected from municipal supply wells were analyzed for radium-226, and all of the concentrations were considered detectable concentrations (greater than their associated sample-specific critical level). Three source-water samples collected were analyzed for radon-222, and all of the concentrations were substantially greater than the associated sample-specific critical level.

Texas↗

Relations of specific conductance to streamflow and selected water-quality characteristics of the Arkansas River basin, Colorado

Areal, seasonal, and long-term variations in the specific conductance of surface and groundwater in the Arkansas River basin of Colorado were evaluated and relations of specific conductance to stream-flow and to concentrations of dissolved solids and major ions were determined as part of an effort to develop a comprehensive hydrologic model of the basin. Mean specific conductance of surface and groundwater was smallest in the upper basin and increased downstream. Smallest mean specific conductance occurred during summer runoff, and largest mean specific conductance occurred during spring and fall low flows. Trends in specific conductance occurred at 18 of 31 surface-water stations and in flow-adjusted specific conductance at 14 of 24 surface-water stations. Logarithmic relations of specific conductance to streamflow were determined for 69 stations. Significant seasonal differences in the relations illustrate the effect of basin characteristics on the relation of specific conductance to streamflow. Relations of specific conductance to dissolved-solids concentration were determined for 28 surface-water stations and for groundwater in alluvial aquifers along the Arkansas River. Relations of specific conductance to concentrations of major ions were determined for 26 surface-water stations and for groundwater in alluvial aquifers along the Arkansas River. (USGS)

Water-Resources Investigations Report↗

Aquifer tests in the flood-plain alluvium and Santa Fe group at the Rio Grande near Canutillo, El Paso County, Texas

An aquifer system consisting of the Rio Grande flood-plain alluvium and Santa Fe Group underlying the southern Mesilla Valley in Dona Ana County, New Mexico and El Paso County, Texas has become an important source of water for both municipal and agricultural uses. Determination of aquifer properties is essential in order to evaluate groundwater potential for increasing water demand and potential streamflow depletion of the Rio Grande due to groundwater development. The aquifer system at the Canutillo well field hydrologic section was divided into a shallow, intermediate, and deep zone based on geohydrologic characteristics. Aquifer properties of specific zones at the test site were determined from a series of multiple-well aquifer tests conducted from December 3, 1985 through January 20, 1986. The Rio Grande is hydraulically connected to the shallow flood-plain alluvium. Water generally occurs within the shallow zone under unconfined conditions, within the intermediate zone under semiconfined conditions, and within the deep zone under confined conditions. (USGS)

Water-Resources Investigations Report↗

Groundwater-level elevations in the Denver Basin bedrock aquifers and Upper Black Squirrel Creek alluvial aquifer, El Paso County, Colorado, 2021–24

El Paso County is the second-most populous county in Colorado and is projected to grow another 15 percent by 2030. Within El Paso County is the Upper Black Squirrel Creek Designated Groundwater Basin (Black Squirrel Basin), an area where surface water is scarce and water users rely primarily on groundwater from five different aquifers (the Upper Black Squirrel Creek alluvial aquifer and four bedrock aquifers within the Denver Basin aquifer system: the lower Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers) to meet their needs. Currently (2024), land within the Upper Black Squirrel Creek Basin is primarily used for rural grazing and agriculture; however, municipal development is ongoing. In 2021, the U.S. Geological Survey, in cooperation with the Upper Black Squirrel Creek Ground Water Management District, began a study to establish a baseline dataset and assess the groundwater resources of the aquifers within the Black Squirrel Basin. A network of 39 wells was established in 2021; discrete groundwater-level measurements were made bimonthly. Nine of the 39 wells were equipped with pressure transducers to record hourly groundwater-level data. Seven wells had statistically significant seasonal trends, and trends at 3 wells were negative. For the discrete data, 16 wells had a significant trend for the study period, and 4 wells had negative trends. For the time-series data, 8 wells had significant trends, and 3 wells had negative trends. Potentiometric surface maps were created for this study using discrete, static groundwater levels measured in April 2023. These maps showed the estimated groundwater flow direction from the north-northwest to the south-southeast in the alluvial aquifer and from the northwest to the east-southeast for the lower Dawson and Denver aquifer wells. This study indicates the potential benefit of monitoring wells in the areas near municipal pumping. Additional monitoring could lead to a better understanding of connectivity between aquifers and be an important tool for assessing long-term sustainability of groundwater use.

Colorado↗