SEARCH · Geology Reports

Results for “Groundwater”

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 19 recordsLinked to original sources

USGS2022

Where groundwater seeps: Evaluating modeled groundwater discharge patterns with thermal infrared surveys at the river-network scale

Predicting baseflow dynamics, protecting aquatic habitat, and managing legacy contaminants requires explicit characterization and prediction of groundwater discharge patterns throughout river networks. Using handheld thermal infrared (TIR) cameras, we surveyed 47 km of stream length across the Farmington River watershed (1,570 km 2 ; CT and MA, USA), mapping locations of bank and waterline groundwater discharges based on their thermal signature. Using the observed groundwater discharge locations and predicted groundwater discharge rates from 6 variations of a numerical groundwater-flow model (MODFLOW-NWT), we compared 1) predicted groundwater-discharge rates in areas with and without observed groundwater discharge, 2) spatial patterns of observed and predicted groundwater discharge locations, and 3) density of observed groundwater discharge locations with predicted discharge rates. Five of six models reasonably predicted the spatial patterns of discharge locations along the 5th order mainstem, but fewer models predicted groundwater discharge patterns in smaller streams. Our results highlight 1) the feasibility of using TIR observations to evaluate groundwater models, 2) model parameters that influence discharge prediction accuracy (riverbed sediment and bedrock hydraulic conductivity and river-aquifer connections), and 3) current strengths and future opportunities for improved modeling of groundwater-discharge patterns.

Connecticut, Massachusetts↗

USGS2025

Groundwater hydrology, groundwater and surface-water interactions, water quality, and groundwater-flow simulations for the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado, 2017–19

In 2017, the U.S. Geological Survey, in cooperation with the Upper Arkansas Water Conservancy District, began a study to provide a comprehensive analysis of the Wet Mountain Valley alluvial aquifer, Custer and Fremont Counties, Colorado. The study included collection of data pertaining to groundwater hydrology, groundwater and surface-water interactions, and water quality in the alluvial aquifer. In addition to providing foundational information on the hydrology of the alluvial aquifer, a numerical groundwater-flow model was developed to estimate the potential effects of additional storage of groundwater in the alluvial aquifer. Groundwater-level elevation data from 30 wells were used to estimate groundwater-flow directions in the alluvial aquifer, which were generally from the southwest to northeast, away from the Sangre de Cristo Mountains and towards perennial streams in the center of the valley. Although some seasonal variation was apparent in groundwater-level elevation records, no statistically significant seasonal trends were indicated. Statistically significant long-term trends were indicated in groundwater-level elevation records for 8 of the 30 wells, and of these wells with statistically significant trends, all but 1 indicated a negative trend of groundwater-level elevations. Spatial evaluation of wells with statistically significant negative groundwater-level elevation trends showed many are in areas of denser well drilling for domestic or other uses, indicating increasing groundwater use could potentially be causing groundwater-level elevation declines. There were instances of wells with no statistically significant groundwater-level elevation trends also located in areas of greater density of well completions. Additional investigations may be necessary to more fully characterize the processes responsible for negative groundwater-level elevation trends. Streamflow gain or loss calculations were completed for low flow in 2017–19 and for high flow in 2018 in nine reaches of streams within the study area. Stream reaches of the upper Texas Creek, upper Grape Creek, upper-middle Grape Creek, and Taylor Creek displayed consistent streamflow loss in each period from 2017 to 2019. These stream reaches represent long-term sources of recharge to the alluvial aquifer. Streamflow gain or loss varies through time in other stream reaches (lower Texas Creek, lower-middle Grape Creek, lower Grape Creek below Westcliffe, and lower Grape Creek above DeWeese Reservoir). The temporally variable behavior indicates these stream reaches may be sources of groundwater recharge or areas of groundwater discharge, likely depending on temporal dynamics between the elevation of the water table and the stream. Water-quality samples were collected from 10 groundwater wells and 10 stream sites during September through November 2019. All groundwater and stream samples were analyzed for major and trace elements and stable isotopes of water. A subset of groundwater samples was also analyzed for the environmental tracers sulfur hexafluoride, tritium, and noble gases. Comparison of water-quality results to U.S. Environmental Protection Agency drinking water-quality standards indicated no constituents exceeded primary standards for human health. Spatial evaluation of water quality indicated the concentrations of various constituents are likely controlled by groundwater and surface-water interactions and by spatial variability in bedrock geology underlying the alluvial aquifer. Specifically, streams shown to gain from groundwater had water chemistry constituent compositions similar to groundwater, whereas streams exiting the Sangre de Cristo Mountains tended to have compositions consistent with snowmelt. Groundwater geochemistry appeared to be partially controlled by oxidation-reduction processes and by proximity to igneous rocks in the Wet Mountains. Environmental tracers used to estimate groundwater age indicated all sampled groundwater contained tracers representing modern recharge (approximately less than 65 years old) but mixing of premodern recharge (approximately more than 65 years old) also occurs. Spatial evaluation of environmental tracers indicated large faults may be conduits for upwelling of older groundwater. No trends were observed in groundwater age with well depth, indicating all sampled wells are located within the zone of active groundwater flow. The presence of modern groundwater in wells with statistically significant negative groundwater-level elevation trends indicates groundwater storage depletions may be partially offset by capture of modern recharge. Repeated sampling of groundwater age would be necessary, however, to determine if any trends in groundwater age exist, which may indicate changing groundwater recharge, storage, or discharge. Additional investigations could also consider quantifying groundwater age in deeper wells to more fully define the depth of active groundwater flow. A numerical groundwater-flow model was developed to estimate components of the water budget, simulate groundwater and surface-water interactions, and evaluate the potential effects of aquifer storage and recovery. Simulated groundwater-level elevations from the calibrated groundwater-flow model are similar to the observed pattern of groundwater-level elevations with higher elevations in the western part of the study area along the Sangre de Cristo Mountains. Simulated water-budget components indicate most of the recharge to the alluvial aquifer is derived from streamflow losses, which is consistent with observations of losing streams along the mountain front. The largest groundwater discharge component of the alluvial aquifer was to streams in the center of the valley, where observations of stream gain or loss indicated the predominance of gaining conditions. Comparison of groundwater and surface-water interactions between the calibrated groundwater-flow model for 2000-19 (the base-case model) and a simulation including additional recharge, representing potential aquifer storage and recovery operations, indicated the additional recharge distributed throughout the area had minimal effects on streamflow in the nearby Grape Creek. An analysis of subregional groundwater budgets showed approximately 54 percent of the additional recharge flowed back to nearby Grape Creek, and the other 46 percent was distributed laterally into adjacent cells in the alluvial aquifer. The comparison of simulations and subregional water budget show the additional recharge did not substantially alter groundwater-level elevations or basin wide groundwater storage. Although the analysis of additional recharge provided in the numerical groundwater-flow model considers only one of many possible recharge scenarios, the model provides a useful tool that could be modified for various scenarios to understand potential effects of managed aquifer recharge.

Colorado↗

USGS2024

Groundwater hydrology, groundwater and surface-water interactions, aquifer testing, and groundwater-flow simulations for the Fountain Creek alluvial aquifer, near Colorado Springs, Colorado, 2018–20

From 2018 through 2020, the U.S. Geological Survey, in cooperation with the Air Force Civil Engineering Center, conducted an integrated study of the Fountain Creek alluvial aquifer located near Colorado Springs, Colorado. The objective of the study was to characterize hydrologic conditions for the alluvial aquifer pertinent to the potential for transport of solutes. Specific goals of this report were to characterize the groundwater hydrology of the area, to quantify groundwater and surface-water interactions, to estimate hydraulic properties of the aquifer using aquifer testing, and to complete numerical simulations of groundwater flow. Synoptic groundwater-level elevation measurements completed throughout this study, and as part of other U.S. Geological Survey programs between 1994 and 2020, indicate groundwater-level elevations fluctuate on annual and interannual timeframes. Groundwater-level fluctuations likely were caused by temporally variable groundwater recharge and discharge components in the area, with many wells showing maximum groundwater-level elevations during the winter months (November through March). From an interannual perspective, groundwater-level fluctuations appear to have reached maximum values during 2000 to 2003, decreased during 2003 to 2006, and remained relatively constant since that time, with the exception of several wells which have displayed rising groundwater-level elevations since 2018. Spatial evaluation of groundwater-level elevations indicates groundwater flow is generally from northeast to southwest within the vicinity of several alluvial paleochannels occurring along the northeastern margin of the aquifer. Within the center of the aquifer along Fountain Creek, groundwater flow is generally from north to south, approximately paralleling surface-water flow. To quantitatively understand the potential effect of groundwater recharge and groundwater pumping on fluctuations in groundwater-level elevation, a statistical transfer-function-noise model was applied. Results of the statistical model indicate throughout most of the aquifer, fluctuations were primarily the result of recharge seasonality. In the main stem of the aquifer where groundwater pumping wells were more concentrated, however, groundwater-level elevation fluctuations were also attributable to groundwater pumping through time. Three-dimensional evaluation of the aquifer geometry near Fountain Creek was combined with synoptic streamflow measurement and accounting of stream gains and losses to evaluate groundwater and surface-water interactions in the study area. Streamflow gain or loss calculations indicate Fountain Creek both gains from and loses flow to the alluvial aquifer, and gaining or losing reaches of the stream may be partially controlled by the depth to bedrock near the stream. Reaches with streamflow gains tend to coincide with areas where the estimated depth to bedrock is decreasing, meaning the alluvial aquifer is likely thinning in these areas and groundwater-flow paths may be converging and discharging groundwater to the stream. Losing reaches tended to coincide with locally greater depth to bedrock where the alluvial aquifer is likely thicker and has greater storage potential for surface water lost from Fountain Creek. Results of aquifer testing indicate hydraulic conductivity, estimated from slug tests and single-well pumping tests, ranged from 0.32 to 1,410 feet per day (ft/d) and 4.13 to 664 ft/d, respectively. These results are similar to the range of values from previous aquifer tests in the study area. Hydraulic conductivities from aquifer testing for this study were generally greater than the estimates of previous slug tests and had a mean value less than the estimates from previous pumping tests. Spatial evaluation of aquifer testing results indicates hydraulic conductivity tends to be greater in the main stem of the alluvial aquifer and lower in paleochannels upgradient from the main stem of the aquifer. The spatial variation in hydraulic conductivity may be attributed to the geomorphologic processes that formed the alluvial aquifer. Compacted sediment in the paleochannels has not been potentially transported sufficient distance to cause grain-size sorting, resulting in a poorly sorted deposit and lower hydraulic conductivities. In the central portion of the alluvial aquifer, near Fountain Creek, the sediments have been transported farther from their source areas and are likely better sorted, removing finer grained sediments that would cause lower hydraulic conductivity. A numerical groundwater-flow model was calibrated for the Fountain Creek alluvial aquifer for 2000–19 to simulate water-budget components, groundwater-flow directions, and groundwater-flow paths. The model simulated precipitation recharge, groundwater and surface-water interactions, evapotranspiration, high-volume groundwater pumping by pumping wells, and external inflows and outflows occurring along the boundaries of the alluvial aquifer. Model calibration was completed using manual and automated approaches, the latter of which assisted in quantifying model results sensitivity to input parameters. The calibrated model corresponds well with groundwater-level elevation observations, with a mean residual (observed minus simulated groundwater-level elevation) equal to −0.60 feet. Simulated groundwater base flow to streams was typically within 10 percent of base flow estimated by independent methods. Groundwater and surface-water interactions represented the largest water-budget components of the aquifer, with the second largest groundwater discharge component coming from pumping wells. Groundwater and surface-water interactions represent both the largest gain and loss terms in the water budget, because these interactions differ spatially, meaning in some areas of the model domain groundwater is being recharged by streams, whereas in other areas, groundwater is discharged to streams. Estimates of advective groundwater-flow paths indicate pumping wells may capture groundwater recharged from losing streams and groundwater that flows into the main stem of the alluvial aquifer from paleochannels.

Colorado↗

USGS2023

Unstructured-grid approach to develop high-fidelity groundwater model to understand groundwater flow and storage responses to excessive groundwater withdrawals in the Southern Hills aquifer system in southeastern Louisiana (USA)

Study region The Southern Hills aquifer system in the Louisiana Capital Area Groundwater Conservation District (CAGCD), USA. Study focus The Southern Hills aquifer system provides abundant groundwater for public and industrial supplies in the CAGCD. Groundwater depletion, saltwater intrusion, and land subsidence are potential concerns due to prolonged excessive groundwater withdrawals. This study develops a high-fidelity groundwater flow model utilizing a complex unstructured grid to investigate groundwater flow and storage responses to excessive groundwater withdrawals for the Southern Hills aquifer system in the CAGCD. The groundwater model incorporates the Mississippi River alluvial aquifer down to the Miocene sands extending to depths around 1 km. New hydrological insights Groundwater modeling results indicate large cones of depression in the Evangeline and Jasper formations in the Baton Rouge area due to prolonged groundwater withdrawals. Low-permeability faults are inferred by significant groundwater level difference across the faults. While local groundwater storage depletion in deeper aquifers is evident, overall estimated groundwater storage changes of the Southern Hills aquifer system in the CAGCD are close to zero in the past two decades, indicating insignificant groundwater storage changes. This is attributed to dominant interactions between the major rivers and the shallower alluvial aquifer. In addition, the simulated groundwater storage changes exhibit patterns similar to those derived by the Gravity Recovery and Climate Experiment (GRACE) model that has been used in evaluation of groundwater depletion in many regional studies.

Louisiana↗

USGS2017

Simulation and assessment of groundwater flow and groundwater and surface-water exchanges in lakes of the northeast Twin Cities Metropolitan Area, Minnesota, 2003 through 2013: Chapter B of Water levels and groundwater and surface-water exchanges in lakes of the northeast Twin Cities Metropolitan Area, Minnesota, 2002 through 2015

Water levels during 2003 through 2013 were less than mean water levels for the period 1925–2013 for several lakes in the northeast Twin Cities Metropolitan Area in Minnesota. Previous periods of low lake-water levels generally were correlated with periods with less than mean precipitation. Increases in groundwater withdrawals and land-use changes have brought into question whether or not recent (2003–13) lake-water-level declines are solely caused by decreases in precipitation. A thorough understanding of groundwater and surface-water exchanges was needed to assess the effect of water-management decisions on lake-water levels. To address this need, the U.S. Geological Survey, in cooperation with the Metropolitan Council and the Minnesota Department of Health, developed and calibrated a three-dimensional, steady-state groundwater-flow model representing 2003–13 mean hydrologic conditions to assess groundwater and lake-water exchanges, and the effects of groundwater withdrawals and precipitation on water levels of 96 lakes in the northeast Twin Cities Metropolitan Area. Lake-water budgets for the calibrated groundwater-flow model indicated that groundwater is flowing into lakes in the northeast Twin Cities Metropolitan Area and lakes are providing water to underlying aquifers. Lake-water outflow to the simulated groundwater system was a major outflow component for Big Marine Lake, Lake Elmo, Snail Lake, and White Bear Lake, accounting for 45 to 64 percent of the total outflows from the lakes. Evaporation and transpiration from the lake surface ranged from 19 to 52 percent of the total outflow from the four lakes. Groundwater withdrawals and precipitation were varied from the 2003‒13 mean values used in the calibrated model (30-percent changes in groundwater withdrawals and 5-percent changes in precipitation) for hypothetical scenarios to assess the effects of groundwater withdrawals and precipitation on water budgets and levels in Big Marine Lake, Snail Lake, and White Bear Lake. Simulated lake-water levels and budgets for Snail Lake and White Bear Lake were affected by 30-percent changes in groundwater withdrawals and 5-percent changes in precipitation in the area, whereas the water level in Big Marine Lake was mainly affected by 5-percent precipitation changes. The effects of groundwater withdrawals on the lake-water levels depend on the number of wells and amount of withdrawals from wells near the lakes. Although lake-water levels are sensitive to precipitation changes, increases in groundwater withdrawals during dry periods exacerbate lake-water level declines. The calibrated, groundwater-flow model is a tool that water-resources managers can use to address future water management issues in the northeast Twin Cities Metropolitan Area.

Minnesota↗

USGS2014

Hydrologic conditions in urban Miami-Dade County, Florida, and the effect of groundwater pumpage and increased sea level on canal leakage and regional groundwater flow

The extensive and highly managed surface-water system in southeastern Florida constructed during the 20th Century has allowed for the westward expansion of urban and agricultural activities in Miami-Dade County. In urban areas of the county, the surface-water system is used to (1) control urban flooding, (2) supply recharge to production well fields, and (3) control seawater intrusion. Previous studies in Miami-Dade County have determined that on a local scale, leakage from canals adjacent to well fields can supply a large percentage (46 to 78 percent) of the total groundwater pumpage from production well fields. Canals in the urban areas also receive seepage from the Biscayne aquifer that is derived from a combination of local rainfall and groundwater flow from Water Conservation Area 3 and Everglades National Park, which are west of urban areas of Miami-Dade County. To evaluate the effects of groundwater pumpage on canal leakage and regional groundwater flow, the U.S. Geological Survey (USGS) developed and calibrated a coupled surface-water/groundwater model of the urban areas of Miami-Dade County, Florida. The model was calibrated by using observation data collected from January 1997 through December 2004. The model calibration was verified using observation data collected from January 2005 through December 2010. A 1-year warmup period (January 1996 through December 1996) was added prior to the start of the calibration period to reduce the effects of inaccurate initial conditions on model calibration. The model is designed to simulate surface-water stage and discharge in the managed canal system and dynamic canal leakage to the Biscayne aquifer as well as seepage to the canal from the aquifer. The model was developed using USGS MODFLOW–NWT with the Surface-Water Routing (SWR1) Process to simulate surface-water stage, surface-water discharge, and surface-water/groundwater interaction and the Seawater Intrusion (SWI2) Package to simulate seawater intrusion, respectively. Automated parameter estimation software (PEST) and highly parameterized inversion techniques were used to calibrate the model to observed surface-water stage, surface-water discharge, net surface-water subbasin discharge, and groundwater level data from 1997 through 2004 by modifying hydraulic conductivity, specific storage coefficients, specific yield, evapotranspiration parameters, canal roughness coefficients (Manning’s n values), and canal leakance coefficients. Tikhonov regularization was used to produce parameter distributions that provide an acceptable fit between model outputs and observation data, while simultaneously minimizing deviations from preferred values based on field measurements and expert knowledge. Analytical and simulated water budgets for the period from 1996 through 2010 indicate that most of the water discharging through the salinity control structures is derived from within the urban parts of the study area and that, on average, the canals are draining the Biscayne aquifer. Simulated groundwater discharge from the urban areas to the coast is approximately 7 percent of the total surface-water inflow to Biscayne Bay and is consistent with previous estimates of fresh groundwater discharge to Biscayne Bay. Simulated groundwater budgets indicate that groundwater pumpage in some surface-water basins ranges from 13 to 27 percent of the sum of local sources of groundwater inflow. The largest percentage of groundwater pumpage to local sources of groundwater inflow occurs in the basins that have the highest pumping rates (C–2 and C–100 Basins). The ratio of groundwater pumpage to simulated local sources of groundwater inflow is less than values calculated in previous local-scale studies. The position of the freshwater-seawater interface at the base of the Biscayne aquifer did not change notably during the simulation period (1996–2010), consistent with the similar positions of the interface in 1984, 1995, and 2011 under similar hydrologic and groundwater pumping conditions. Landward movement of the freshwater-seawater interface above the base of the aquifer is more prone to occur during relatively dry years. The model was used to evaluate the effect of increased groundwater pumpage and (or) increased sea level on canal leakage, regional groundwater flow, and the position of the freshwater-seawater interface. Permitted groundwater pumping rates, which generally exceed historical groundwater pumping rates, were used for Miami-Dade County Water and Sewer Department groundwater pumping wells in the base-case future scenario. Base-case future and increased pumping scenario results suggest seawater intrusion may occur at the Miami-Springs well field if the Miami Springs, Hialeah, and Preston well fields are operated using current permitted groundwater pumping rates. Scenario simulations also show that, in general, the canal system limits the adverse effects of proposed groundwater pumpage increases on water-level changes and saltwater intrusion. Proposed increases (up to a 7 percent increase) in groundwater pumpage do not have a notable effect on movement of the freshwater-seawater interface. Increased groundwater pumpage increased lateral groundwater inflow into basins subject to additional groundwater pumpage; however, most (55 percent) of the additional groundwater extracted from pumping wells was supplied by changes in canal seepage and leakage in urban areas of the model. Increased sea level caused increased water-table elevations in urban areas and decreased hydraulic gradients across the system; the largest increases in water-table elevations occurred seaward of the salinity control structures. The extent of flood-prone areas and the percentage of time water-table elevations in flood-prone areas were less than 0.5 foot below land surface increased with increased sea level. Increased sea level also resulted in landward migration of the freshwater-seawater interface; the largest changes in the position of the interface occurred seaward of the salinity control structures except in parts of the model area that were inundated by increased sea level. Decreased water-table gradients reduced groundwater inflow, groundwater outflow, canal exchanges, surface-water inflow, and surface-water outflow through salinity control structures. Results for the scenario that evaluated the combination of increased groundwater pumpage and increased sea level did not differ substantially from the scenario that evaluated increased sea level alone. Groundwater inflow, groundwater outflow, and canal exchanges were reduced in urban areas of the study area as a result of decreased water-table gradients across the system, although reductions were less than those in the increased sea-level scenario. The decline in groundwater levels caused by increased groundwater pumpage was less under the increased sea-level scenario than under the increased groundwater-pumpage scenario. The largest reductions in surface-water outflow from the salinity control structures occurred with increased sea level and increased groundwater pumpage. The model was designed specifically to evaluate the effect of groundwater pumpage on canal leakage at the surface-water-basin scale and thus may not be appropriate for (1) predictions that are dependent on data not included in the calibration process (for example, subdaily simulation of high-intensity events and travel times) and (or) (2) hydrologic conditions that are substantially different from those during the calibration and verification periods. The reliability of the model is limited by the conceptual model of the surface-water and groundwater system, the spatial distribution of physical properties, the scale and discretization of the system, and specified boundary conditions. Some of the model limitations are manifested in model errors. Despite these limitations, however, the model represents the complexities of the interconnected surface-water and groundwater systems that affect how the systems respond to groundwater pumpage, sea-level rise, and other hydrologic stresses. The model also quantifies the relative effects of groundwater pumpage and sea-level rise on the surface-water and groundwater systems.

Florida↗

USGS2015

Evaluation of groundwater levels in the South Platte River alluvial aquifer, Colorado, 1953-2012, and design of initial well networks for monitoring groundwater levels

The South Platte River and underlying alluvial aquifer form an important hydrologic resource in northeastern Colorado that provides water to population centers along the Front Range and to agricultural communities across the rural plains. Water is regulated based on seniority of water rights and delivered using a network of administration structures that includes ditches, reservoirs, wells, impacted river sections, and engineered recharge areas. A recent addendum to Colorado water law enacted during 2002-2003 curtailed pumping from thousands of wells that lacked authorized augmentation plans. The restrictions in pumping were hypothesized to increase water storage in the aquifer, causing groundwater to rise near the land surface at some locations. The U.S. Geological Survey (USGS), in cooperation with the Colorado Water Conservation Board and the Colorado Water Institute, completed an assessment of 60 years (yr) of historical groundwater-level records collected from 1953 to 2012 from 1,669 wells. Relations of "high" groundwater levels, defined as depth to water from 0 to 10 feet (ft) below land surface, were compared to precipitation, river discharge, and 36 geographic and administrative attributes to identify natural and human controls in areas with shallow groundwater. Averaged per decade and over the entire aquifer, depths to groundwater varied between 24 and 32 ft over the 60-yr record. The shallowest average depth to water was identified during 1983-1992, which also recorded the highest levels of decadal precipitation. Average depth to water was greatest (32 ft) during 1953-1962 and intermediate (30 ft) in the recent decade (2003-2012) following curtailment of pumping. Between the decades 1993-2002 and 2003-2012, groundwater levels declined about 2 ft across the aquifer. In comparison, in areas where groundwater levels were within 20 ft of the land surface, observed groundwater levels rose about 0.6 ft, on average, during the same period, which demonstrated preferential rise in areas with shallow groundwater. Approximately 29 percent of water-level observations were identified as high groundwater in the South Platte River alluvial aquifer over the 60-yr record. High groundwater levels were found in 17 to 33 percent of wells examined by decade, with the largest percentages occurring over three decades from 1963 to 1992. The recent decade (2003-2012) exhibited an intermediate percentage (25 percent) of wells with high groundwater levels but also had the highest percentage (30 percent) of high groundwater observations, although results by observations were similar (26-29 percent) over three decades prior, from 1963 to 1992. Major sections of the aquifer from north of Sterling to Julesburg and areas near Greeley, La Salle, and Gilcrest were identified with the highest frequencies of high groundwater levels. Changes in groundwater levels were evaluated using Kendal line and least trimmed squares regression methods using a significance level of 0.01 and statistical power of 0.8. During 2003-2012, following curtailment of pumping, 88 percent of wells and 81 percent of subwatershed areas with significant trends in groundwater levels exhibited rising water levels. Over the complete 60-yr record, however, 66 percent of wells and 57 percent of subwatersheds with significant groundwater-level trends still showed declining water levels; rates of groundwater-level change were typically less than 0.125 ft/yr in areas near the South Platte River, with greater declines along the southern tributaries. In agreement, 58 percent of subwatersheds evaluated between 1963-1972 and 2003-2012 showed net declines in average decadal groundwater levels. More areas had groundwater decline in upgradient sections to the west and rise in downgradient sections to the east, implying a redistribution of water has occurred in some areas of the aquifer. Precipitation was identified as having the strongest statistically significant correlations to river discharge over annual and decadal periods (Pearson correlation coefficients of 0.5 and 0.8, respectively, and statistical significance defined by p-values less than 0.05). Correlation coefficients between river discharge and frequency of high groundwater levels were statistically significant at 0.4 annually and 0.6 over decadal periods, indicating that periods of high river flow were often coincident with high groundwater conditions. Over seasonal periods in five of the six decades examined, peak high groundwater levels occurred after spring runoff from July to September when administrative structures were most active. Between 1993-2002 and 2003-2012, groundwater levels rose while river discharge decreased, in part from greater reliance on surface water and curtailed pumping from wells without augmentation plans. Geographic attributes of elevation and proximity to streams and rivers showed moderate correlations to high groundwater levels in wells used for observing groundwater levels (correlation coefficients of 0.3 to 0.4). Local depressions and regional lows within the aquifer were identified as areas of potential shallow groundwater. Wells close to the river regularly indicated high groundwater levels, while those within depleted tributaries tended to have low frequencies of high groundwater levels. Some attributes of administrative structures were spatially correlated to high groundwater levels at moderate to high magnitudes (correlation coefficients of 0.3 to 0.7). The number of affected river reaches or recharge areas that surround a well where groundwater levels were observed and its distance from the nearest well field showed the strongest controls on high groundwater levels. Influences of administrative structures on groundwater levels were in some cases local over a mile or less but could extend to several miles, often manifesting as diffuse effects from multiple surrounding structures. A network of candidate monitoring wells was proposed to initiate a regional monitoring program. Consistent monitoring and analysis of groundwater levels will be needed for informed decisions to optimize beneficial use of water and to limit high groundwater levels in susceptible areas. Finalization of the network will require future field reconnaissance to assess local site conditions and discussions with State authorities.

Colorado↗

USGS2016

Application of hydrogeology and groundwater-age estimates to assess the travel time of groundwater at the site of a landfill to the Mahomet Aquifer, near Clinton, Illinois

The U.S. Geological Survey used interpretations of hydrogeologic conditions and tritium-based groundwater age estimates to assess the travel time of groundwater at a landfill site near Clinton, Illinois (the “Clinton site”) where a chemical waste unit (CWU) was proposed to be within the Clinton landfill unit #3 (CLU#3). Glacial deposits beneath the CWU consist predominantly of low-permeability silt- and clay-rich till interspersed with thin (typically less than 2 feet in thickness) layers of more permeable deposits, including the Upper and Lower Radnor Till Sands and the Organic Soil unit. These glacial deposits are about 170 feet thick and overlie the Mahomet Sand Member of the Banner Formation. The Mahomet aquifer is composed of the Mahomet Sand Member and is used for water supply in much of east-central Illinois. Eight tritium analyses of water from seven wells were used to evaluate the overall age of recharge to aquifers beneath the Clinton site. Groundwater samples were collected from six monitoring wells on or adjacent to the CLU#3 that were open to glacial deposits above the Mahomet aquifer (the upper and lower parts of the Radnor Till Member and the Organic Soil unit) and one proximal production well (approximately 0.5 miles from the CLU#3) that is screened in the Mahomet aquifer. The tritium-based age estimates were computed with a simplifying, piston-flow assumption: that groundwater moves in discrete packets to the sampled interval by advection, without hydrodynamic dispersion or mixing. Tritium concentrations indicate a recharge age of at least 59 years (pre-1953 recharge) for water sampled from deposits below the upper part of the Radnor Till Member at the CLU#3, with older water expected at progressively greater depth in the tills. The largest tritium concentration from a well sampled by this study (well G53S; 0.32 ± 0.10 tritium units) was in groundwater from a sand deposit in the upper part of the Radnor Till Member; the shallowest permeable unit sampled by this study. That result indicated that nearly all groundwater sampled from well G53S entered the aquifer as recharge before 1953. Tritium was detected in a trace concentration in one sample from a second monitoring well open to the upper part of the Radnor Till Member (well G07S; 0.11 ± 0.09 tritium units), and not detected in samples collected from two monitoring wells open to a sand deposit in the lower part of the Radnor Till Member, from two samples collected from two monitoring wells open to the Organic Soil unit, and in two samples collected from a production well screened in the middle of the Mahomet aquifer (a groundwater sample and a sequential replicate sample). The lack of tritium in five of the six groundwater samples collected from the shallow permeable units beneath CLU#3 site and the two samples from the one Mahomet aquifer well indicates an absence of post-1952 recharge. Groundwater-flow paths that could contribute post-1952 recharge to the lower part of the Radnor Till Member, the Organic Soil unit, or the Mahomet aquifer at the CLU#3 are not indicated by these data. Hypothetical two-part mixtures of tritium-dead, pre-1953 recharge water and decay-corrected tritium concentrations in post-1952 recharge were computed and compared with tritium analyses in groundwater sampled from monitoring wells at the CLU#3 site to evaluate whether tritium concentrations in groundwater could be represented by mixtures involving some post-1952 recharge. Results from the hypothetical two-part mixtures indicate that groundwater from monitoring well (G53S) was predominantly composed of pre-1953 recharge and that if present, younger, post-1955 recharge, contributed less than 2.5 percent to that sample. The hypothetical two-part mixing results also indicated that very small amounts of post-1952 recharge composing less than about 2.5 percent of the sample volume could not be distinguished in groundwater samples with tritium concentrations less than about 0.15 TU. The piston-flow based age of recharge determined from the tritium concentration in the groundwater sample from monitoring well G53S yielded an estimated maximum vertical velocity from the land surface to the upper part of the Radnor Till Member of 0.85 feet per year or less. This velocity, ifassumed to apply to the remaining glacial till deposits above the Mahomet aquifer, indicates that recharge flows through the 170 feet of glacial deposits between the base of the proposed chemical waste unit and the top of the Mahomet aquifer in a minimum of 200 years or longer. Analysis of hydraulic data from the site, constrained by a tritium-age based maximum groundwater velocity estimate, computed minimum estimates of effective porosity that range from about 0.021 to 0.024 for the predominantly till deposits above the Mahomet aquifer. Estimated rates of transport of recharge from land surface to the Mahomet aquifer for the CLU#3 site computed using the Darcy velocity equation with site-specific data were about 260 years or longer. The Darcy velocity-based estimates were computed using values that were based on tritium data, estimates of vertical velocity and effective porosity and available site-specific data. Solution of the Darcy velocity equation indicated that maximum vertical groundwater velocities through the deposits above the aquifer were 0.41 or 0.61 feet per year, depending on the site-specific values of vertical hydraulic conductivity (laboratory triaxial test values) and effective porosity used for the computation. The resulting calculated minimum travel times for groundwater to flow from the top of the Berry Clay Member (at the base of the proposed chemical waste unit) to the top of the Mahomet aquifer ranged from about 260 to 370 years, depending on the velocity value used in the calculation. In comparison, plausible travel times calculated using vertical hydraulic conductivity values from a previously published regional groundwater flow model were either slightly less than or longer than those calculated using site data and ranged from 230 to 580 years. Tritium data from 1996 to 2011 USGS regional sampling of groundwater from domestic wells in the confined part of the Mahomet aquifer—which are 2.5 to about 40 miles from the Clinton site—were compared with site-specific data from a production well at the Clinton site. Tritium-based groundwater-age estimates indicated predominantly pre- 1953 recharge dates for USGS and other prior regional samples of groundwater from domestic wells in the Mahomet aquifer. These results agreed with the tritium-based, pre-1953 recharge age estimated for a groundwater sample and a sequential replicate sample from a production well in the confined part of the Mahomet aquifer beneath the Clinton site. The regional tritium-based groundwater age estimates also were compared with pesticide detections in samples from distal domestic wells in the USGS regional network that are about 2.5 to 40 miles from the Clinton site to identify whether very small amounts of post-1952 recharge have in places reached confined parts of the Mahomet aquifer at locations other than the Clinton site in an approximately 2,000 square mile area of the Mahomet aquifer. Very small amounts of post-1952 recharge were defined in this analysis as less than about 2.5 percent of the total recharge contributing to a groundwater sample, based on results from the two-part mixing analysis of tritium data from the Clinton site. Pesticide-based groundwater-age estimates based on 22 detections of pesticides (13 of these detections were estimated concentrations), including atrazine, deethylatrazine (2-Chloro-4-isopropylamino-6-amino- s-triazine), cyanazine, diazinon, metolachlor, molinate, prometon, and trifluralin in groundwater samples from 10 domestic wells 2.5 to about 40 miles distant from the Clinton site indicate that very small amounts of post-1956 to post-1992 recharge can in places reach the confined part of the Mahomet aquifer in other parts of central Illinois. The relative lack of tritium in these samples indicate that the amounts of post-1956 to post-1992 recharge contributing to the 10 domestic wells were a very small part of the overall older groundwater sampled from those wells. The flow process by which very small amounts of pesticide-bearing groundwater reached the screened intervals of the 10 domestic wells could not be distinguished between well-integrity related infiltration and natural hydrogeologic features. Potential explanations include: (1) infiltration through man-made avenues in or along the well, (2) flow of very small amounts of post-1956 to post-1992 recharge through sparsely distributed natural permeable aspects of the glacial till and diluted by mixing with older groundwater, or (3) a combination of both processes. Presuming the domestic wells sampled by the USGS in 1996–2011 in the regional study of the confined part of the Mahomet aquifer are adequately sealed and produce groundwater that is representative of aquifer conditions, the regional tritium and pesticide-based groundwater-age results indicate substantial heterogeneity in the glacial stratigraphy above the Mahomet aquifer. The pesticide-based groundwater-age estimates from the domestic wells distant from the Clinton site also indicate that parts of the Mahomet aquifer with the pesticide detections can be susceptible to contaminant sources at the land surface. The regional pesticide and tritium results from the domestic wells further indicate that a potential exists for possible contaminants from land surface to be transported through the glacial drift deposits that confine the Mahomet aquifer in other parts of central Illinois at faster rates than those computed for recharge at the Clinton site, including CLU#3. This analysis indicates the potential value of sub-microgram-per-liter level concentrations of land-use derived indicators of modern recharge to indicate the presence of very small amounts of modern, post-1952 age recharge in overall older, pre-1953 age groundwater.

Illinois↗

USGS2016

Groundwater quality, age, and susceptibility and vulnerability to nitrate contamination with linkages to land use and groundwater flow, Upper Black Squirrel Creek Basin, Colorado, 2013

The Upper Black Squirrel Creek Basin is located about 25 kilometers east of Colorado Springs, Colorado. The primary aquifer is a productive section of unconsolidated deposits that overlies bedrock units of the Denver Basin and is a critical resource for local water needs, including irrigation, domestic, and commercial use. The primary aquifer also serves an important regional role by the export of water to nearby communities in the Colorado Springs area. Changes in land use and development over the last decade, which includes substantial growth of subdivisions in the Upper Black Squirrel Creek Basin, have led to uncertainty regarding the potential effects to water quality throughout the basin. In response, the U.S. Geological Survey, in cooperation with Cherokee Metropolitan District, El Paso County, Meridian Service Metropolitan District, Mountain View Electric Association, Upper Black Squirrel Creek Groundwater Management District, Woodmen Hills Metropolitan District, Colorado State Land Board, and Colorado Water Conservation Board, and the stakeholders represented in the Groundwater Quality Study Committee of El Paso County conducted an assessment of groundwater quality and groundwater age with an emphasis on characterizing nitrate in the groundwater. Groundwater-quality samples were collected from 50 randomly selected wells between May and June 2013. The samples were analyzed for major ions, nutrients, dissolved gases, tritium ( 3 H), chlorofluorocarbons (CFC-11, CFC-12, and CFC-113), and fuel products (such as benzene, toluene, ethylbenzene, and xylenes). None of the groundwater samples exceeded the U.S. Environmental Protection Agency (EPA) National Primary Drinking Water Regulations for primary maximum contaminant levels (MCL) for major ions. Secondary maximum contaminant levels, which are not health concerns and affect mainly taste, color, or odor of the water, were observed in rare instances for pH (2 samples), chloride (1 sample), iron (3 samples), and manganese (8 samples). The secondary maximum contaminant level for total dissolved solids was also exceeded for two samples. Nitrate (nitrite plus nitrate as nitrogen in groundwater) was elevated above the estimated background concentration of natural recharge waters of 1 milligram per liter (mg/L) in 44 of the 50 wells sampled and showed a median concentration of 5.4 mg/L. Nitrate concentrations were above the MCL of 10 mg/L in 5 of the 50 wells sampled and above half of the EPA MCL (5 mg/L) in 27 of the 50 wells sampled, which included samples above the MCL. Dissolved-oxygen concentrations exceeded 0.5 mg/L in 95 percent of reported values (40 of 42 samples) and exceeded 2.0 mg/L in 90 percent of reported values (38 of 42 samples). The oxidized conditions observed in most areas indicate that nitrate from fertilizers and animal or human waste was geochemically stable and could persist in the groundwater for decades or perhaps longer. A historical analysis of median nitrate concentrations over nearly three decades showed an increase in nitrate of approximately 1 mg/L from 4.3 to 5.4 mg/L, although the increase was not determined to be significantly different using nonparametric statistical methods. Major-ion data indicate that groundwater representative of the primary aquifer was classified as calcium-sodium bicarbonate type water. Other water samples from wells located mainly along the periphery of the primary aquifer had cation-anion compositions consistent with distinct water sources, including groundwater contributions from the underlying bedrock aquifers. The areas with differentiable water sources were located mainly where alluvial deposits were thin and geologic contacts to the underlying bedrock aquifers were relatively shallow. Nitrate concentrations in the groundwater were evaluated for relations to land use. An agricultural region was defined using a sequence of land satellite imagery. Groundwater flow directions interpreted from median water-table elevations measured from 2000 to 2013 were used in conjunction with cropland locations to define the agricultural region boundaries by encompassing potential pathways of nitrate transport in the groundwater from nitrogen-based fertilizers. A statistically significant higher median nitrate concentration was observed for areas inside the agricultural region (6.7 mg/L) compared to areas outside the agricultural region (2.3 mg/L), although median concentrations in both areas were below the MCL (10 mg/L). Median nitrate concentration was also significantly greater in land parcels with septic use (4.9 mg/L) compared to nonseptic parcels (1.7 mg/L). In general, agriculture or septic use was identified as the primary source of nitrate, depending on location, while commercial, county, grazing, and residential land uses were generally secondary sources of nitrate. Apparent groundwater ages were estimated from chlorofluorocarbons (CFC-11, CFC-12, and CFC-113) and tritium ( 3 H) data using models that assumed piston flow and binary mixing (dilution of a young component with old, tracer-free water). The mean and median groundwater ages were about 30 years and the standard deviation was 6 years, indicating that most groundwater in the primary aquifer was “young” water that had recharged to the aquifer over the last few decades (post-1950s). The median fraction of young water was about 71 percent, and the standard deviation was 29 percent. The remaining water predated the 1950s, which may have originated from deeper geologic formations or may represent slow moving groundwater within the primary aquifer. Some of the oldest groundwater ages (older than 30 years) were observed in the upper reaches of the aquifer to the northwest where the primary aquifer is thin and intersects bedrock, supporting the hypothesis of geochemically distinct groundwater entering the primary aquifer from below. Groundwater that had reached the central part of the aquifer from upgradient areas of the basin was variable in age because of differences in flow paths and travel velocities. The groundwater age analysis showed that current (2013) land-use practices could affect water quality over decades to come, and that responses to remedial actions could be slow, especially for constituents, such as nitrate, that are stable under oxidized conditions. Fuel products (including acetone, benzene, diisopropyl ether, ethylbenzene, methyl acetate, methyl tertiary butyl ether (MTBE), methyl tert-pentyl ether, m- + p-xylene, o-xylene, tert-amyl alcohol, tert-butyl alcohol, tert-butyl ethyl ether, and toluene) were analyzed in groundwater from 49 of the 50 wells. Water from seven sites had detections for fuel compounds; all concentrations were below MCL. The results provided assurance of water quality and a valuable baseline to evaluate future trends of fuel constituents as the region is further developed. Probability maps were developed from logistic regression models to examine the likelihood that nitrate concentrations in groundwater exceeded specified levels. Susceptibility analysis examined relations between mid-level (5.0 mg/L) nitrate concentrations and climatic, hydrologic, and geologic variables; the significant variables were identified as depth to groundwater, soil organic matter, and soil water storage to 25-centimeter (cm) depth. The vulnerability assessments included natural factors driving susceptibility but also human factors related to land use and septic use. Vulnerability to low-level (2.5 mg/L) nitrate was related to depth to groundwater, septic zoning, and soil organic matter. The results highlighted that septic zoning affected low-level nitrate concentrations. Vulnerability to mid-level (5.0 mg/L) nitrate was examined using all 50 samples and also with two data outliers removed, which showed relatively high nitrate concentrations but also anomalous water chemistry or were located beyond the primary study area. Vulnerability to mid-level (5.0 mg/L) nitrate using all 50 samples was related to depth to groundwater, land use, septic use within a 500-meter (m) radius, soil water storage to a 25-cm depth, soil organic matter, and whether a location was within the agricultural region. The mid-level (5.0 mg/L) vulnerability model using 48 samples (two outliers removed) produced the best overall fit and was related to the same variables as when using all samples except septic use. The results for mid-level vulnerability provided additional support that septic use was associated with low levels of nitrate in the groundwater. Soil properties and land use were identified as the main drivers of moderate nitrate concentrations. Probabilities of exceeding low-level nitrate concentrations were high in most areas with the lowest probabilities usually to the northwest along thin geologic deposits in the upper part of the basin. The results of this investigation offer the foundational information needed for developing best management practices to mitigate nitrate contamination, basic concepts on water quality to aid public education, and information to guide regulatory measures if policy makers determine this is warranted. Science-based decision making will require continued monitoring and analysis of water quality in the future.

Colorado↗

USGS2018

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

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

Montana, North Dakota, South Dakota, Wyoming↗

USGS2018

Hydrogeologic controls and geochemical indicators of groundwater movement in the Niles Cone and southern East Bay Plain groundwater subbasins, Alameda County, California

Beginning in the 1970s, Alameda County Water District began infiltrating imported water through ponds in repurposed gravel quarries at the Quarry Lakes Regional Park, in the Niles Cone groundwater subbasin, to recharge groundwater and to minimize intrusion of saline, San Francisco Bay water into freshwater aquifers. Hydraulic connection between distinct aquifers underlying Quarry Lakes allows water to recharge the upper aquifer system to depths of 400 feet below land surface, and the Deep aquifer to depths of more than 650 feet. Previous studies of the Niles Cone and southern East Bay Plain groundwater subbasins suggested that these two subbasins may be hydraulically connected. Characterization of storage capacities and hydraulic properties of the complex aquifers and the structural and stratigraphic controls on groundwater movement aids in optimal storage and recovery of recharged water and provides information on the ability of aquifers shared by different water management agencies to fulfill competing storage and extraction demands. The movement of recharge water through the Niles Cone groundwater subbasin from Quarry Lakes and the possible hydraulic connection between the Niles Cone and the southern East Bay Plain groundwater subbasins were investigated using interferometric synthetic aperture radar (InSAR), water-chemistry, and isotopic data, including tritium/helium-3, helium-4, and carbon-14 age-dating techniques. InSAR data collected during refilling of the Quarry Lakes recharge ponds show corresponding ground-surface displacement. Maximum uplift was about 0.8 inches, reasonable for elastic expansion of sedimentary materials experiencing an increase in hydraulic head that resulted from pond refilling. Sodium concentrations increase while calcium and magnesium concentrations in groundwater decrease along groundwater flowpaths from the Niles Cone groundwater subbasin through the Deep aquifer to the northwest toward the southern East Bay Plain groundwater subbasin. Residual effects of pre-1970s intrusion of saline water from San Francisco Bay, including high chloride concentrations in groundwater, are evident in parts of the Niles Cone subbasin. Noble gas recharge temperatures indicate two primary recharge sources (Quarry Lakes and Alameda Creek) in the Niles Cone groundwater subbasin. Although recharge at Quarry Lakes affects hydraulic heads as far as the transition zone between the Niles Cone and East Bay Plain groundwater subbasins (about 5 miles), the effect of recharged water on water quality is only apparent in wells near (less than 2 miles) recharge sources. Groundwater chemistry from upper aquifer system wells near Quarry Lakes showed an evaporated signal (less negative oxygen and hydrogen isotopic values) relative to surrounding groundwater and a tritium concentration (2 tritium units) consistent with recently recharged water from a surface-water impoundment. Uncorrected carbon-14 activities measured in water sampled from wells in the Niles Cone groundwater subbasin range from 16 to 100 percent modern carbon (pmC). The geochemical reaction modeling software NETPATH was used to interpret carbon-14 ages along a flowpath from Quarry Lakes toward the East Bay Plain groundwater subbasin. Model results indicate that changes in groundwater chemistry are controlled by cation exchange on clay minerals and weathering of primary silicate minerals. Old groundwater (lower carbon-14 activities) is characterized by high dissolved silica and pH. Interpreted carbon-14 ages ranged from 830 to more than 7,000 years before present and are less than helium-4 ages that range from 2,000 to greater than 11,000 years before present. The average horizontal groundwater velocity along the studied flowpath, as calculated using interpreted carbon-14 ages, through the Deep aquifer of the Niles Cone groundwater subbasin is between 3 and 12 feet per year. The groundwater velocity decreases near the boundary of the transition zone to the southern East Bay Plain groundwater subbasin to about 0.5 feet per year. These changes may result from water recharged from different sources converging in flowpaths north of the transition zone, or a boundary to flow between the Niles Cone and southern East Bay Plain groundwater subbasins, likely owing to changes in lithology caused by depositional patterns.

California↗

USGS2022

Hydrogeology and simulation of groundwater flow in the Lucerne Valley groundwater basin, California

The Lucerne Valley is in the southwestern part of the Mojave Desert and is about 75 miles northeast of Los Angeles, California. The Lucerne Valley groundwater basin encompasses about 230 square miles and is separated from the Upper Mojave Valley groundwater basin by splays of the Helendale Fault. Since its settlement, groundwater has been the primary source of water for agricultural, industrial, municipal, and domestic uses. Groundwater withdrawal from pumping has exceeded the amount of water recharged to the basin, causing groundwater declines of more than 100 feet between 1917 and 2016 in the center of the basin. The continued withdrawal has resulted in an increase in pumping costs, reduced well efficiency, and land subsidence near Lucerne Lake. Although the volume of pumping has declined in recent years, there is concern that new agricultural growth and limits on imported water will continue to strain the sustainability of the groundwater system. To address these concerns, the U.S. Geological Survey entered into a cooperative agreement with the Mojave Water Agency to develop a better understanding of the Lucerne Valley hydrogeologic system and provide tools to help evaluate and manage the effects of future development in the Lucerne Valley. The objectives of this study were to (1) improve the understanding of the aquifer system, (2) improve the understanding of subsidence in the basin, and (3) incorporate the understanding into a groundwater-flow model that can be used to help manage the groundwater resources in the Lucerne Valley. The model developed for this study covers the period of 1942–2016 and can help evaluate various proposed water-management scenarios during different climatic and hydrologic conditions. The aquifer system consists of a shallow aquifer, a confining unit, and middle and lower aquifers. These layered water-bearing units were identified based on geologic units of the mostly unconsolidated sediments and hydrologic properties. These alluvial deposits consist of clay, silt, sand, and gravel; some places also contain clay and silty clay lacustrine deposits. Several faults act, at least in part, as barriers to groundwater flow on the eastern, southern, and western edges of the basin. Present-day natural recharge is primarily from the infiltration of runoff from the San Bernardino Mountains to the south; however, stable and radioactive isotopes show that groundwater from the middle of the Lucerne Valley was older than about 10,000 years and probably was recharged as infiltration from streams draining the mountains in the Mojave Desert to the north, which probably does not occur under present-day climatic conditions. The annual average natural recharge for 1942–2016, estimated by a Basin Characterization Model, was about 635 acre-feet per year; the average amount of treated wastewater effluent transferred to the Lucerne Valley for artificial recharge annually ranged from about 1,500 to 4,000 acre-feet per year during 1980–2016. Pumpage estimates for 1942–2016 ranged from about 3,000 acre-feet in 1942 to about 18,300 acre-feet in 1984. The total cumulative amount of groundwater removed from the basin by pumping between 1942 and 2016 was estimated to be about 700,000 acre-feet, which was about 10 times greater than the cumulative amount of recharge to the entire Lucerne Valley groundwater basin. Before groundwater development, the direction of groundwater flow was from the southern part of the basin northward to discharge areas near Lucerne Lake, where it discharged through springs along the Helendale Fault and by evapotranspiration. Since the early 1900s, groundwater-level declines have mostly eliminated the areas where natural discharge occurred and exceeded 100 feet in the middle of the basin between the early 1950s and mid-1990s, and as much as 25 feet near the margins from about the mid-1950s to 2000s. A decrease in the rate of pumping after the mid-1990s lessened the hydraulic stress on the middle and lower aquifers and enabled hydraulic heads in the middle of the basin to recover slightly as groundwater near the margins of the basin moved toward the pumping depression. Although trends in groundwater levels in the center of the basin have reversed since the mid-1990s, levels at the basin margins continue to decline as the movement of groundwater from the margins fills the pumping depression and gradually flattens the groundwater table throughout the basin. The long-term extraction of groundwater and associated dewatering of the fine-grained sediments present within the aquifer system has resulted in aquifer compaction and consequently land subsidence, primarily near Lucerne Lake. Analysis of interferometric synthetic aperture radar data shows that almost 11 inches of land subsidence has occurred south of Lucerne Lake between April 1992 and November 2009; less subsidence occurred elsewhere in the basin during this period. This differential land subsidence has caused fissures and cracks in the ground surface, which have buckled the pavement and undercut roads in several locations. The Lucerne Valley Hydrologic Model was developed using the finite-difference groundwater modeling software One Water Hydrologic Model to represent the hydrologic conditions and stresses during 1942–2016. The model has a uniform grid of approximately 92 acres per cell (2,000 feet by 2,000 feet) and has four layers representing the water-bearing units. The results from the calibrated model simulations indicated that groundwater pumpage exceeded recharge, resulting in an estimated net cumulative depletion of groundwater storage (discharge minus recharge) of about 465,000 acre-feet from 1942 to 2016. The model simulated as much as 7.5 feet (90 inches; 2,286 millimeters) of aquifer compaction, which indicates the extensive fine-grained deposits and measured subsidence near Lucerne Lake.

California↗

USGS2010

Death Valley regional groundwater flow system, Nevada and California: Hydrogeologic framework and transient groundwater flow model

A numerical three-dimensional (3D) transient groundwater flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the groundwater flow system and previous less extensive groundwater flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect groundwater flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley regional groundwater flow system (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the groundwater flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural groundwater discharge occurring through evapotranspiration (ET) and spring flow; the history of groundwater pumping from 1913 through 1998; groundwater recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional groundwater flow system. Groundwater flow in the Death Valley region is composed of several interconnected, complex groundwater flow systems. Groundwater flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional groundwater flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and ET are the dominant natural groundwater discharge processes. Groundwater also is withdrawn for agricultural, commercial, and domestic uses. Groundwater flow in the DVRFS was simulated using MODFLOW-2000, the U.S. Geological Survey 3D finitedifference modular groundwater flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 meters (m) on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins simulated by constant-head boundaries. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model generally fit observed heads reasonably well (residuals with absolute values less than 10 meters) with two exceptions: in most areas of nearly flat hydraulic gradient the fit is considered moderate (residuals with absolute values of 10 to 20 meters), and in areas of steep hydraulic gradient along the Eleana Range and western part of Yucca Flat, southern part of the Owlshead Mountains, southern part of the Bullfrog Hills, and the north-northwestern part of the model domain (residuals with absolute values greater than 20 meters). Groundwater discharge residuals are fairly random, with as many areas where simulated flows are less than observed flows as areas where simulated flows are greater. The highest unweighted groundwater discharge residuals occur at Death Valley, Sarcobatus Flat (northeastern area), Tecopa, and early observations at Manse Spring in Pahrump Valley. High weighted-discharge residuals were computed in Indian Springs Valley and parts of Death Valley. Most of these inaccuracies in head and discharge can be attributed to insufficient representation of the hydrogeology in the HFM and(or) discharge estimates, misrepresentation of water levels, and(or) model error associated with grid-cell size. The model represents the large and complex groundwater flow system of the Death Valley region at a greater degree of refinement and accuracy than has been possible previously. The representation of detail provided by the 3D digital hydrogeologic framework model and the numerical groundwater flow model enabled greater spatial accuracy in every model parameter. The lithostratigraphy and structural effects of the hydrogeologic framework; recharge estimates from simulated net infiltration; discharge estimates from ET, spring flow, and pumping; and boundary inflow and outflow estimates all were reevaluated, some additional data were collected, and accuracy was improved. Uncertainty in the results of the flow model simulations can be reduced by improving on the quality, interpretation, and representation of the water-level and discharge observations used to calibrate the model and improving on the representation of the HGU geometries, the spatial variability of HGU material properties, the flow model physical framework, and the hydrologic conditions.

California, Nevada↗

USGS2013

Investigations of groundwater system and simulation of regional groundwater flow for North Penn Area 7 Superfund site, Montgomery County, Pennsylvania

Groundwater in the vicinity of several industrial facilities in Upper Gwynedd Township and vicinity, Montgomery County, in southeast Pennsylvania has been shown to be contaminated with volatile organic compounds (VOCs), the most common of which is the solvent trichloroethylene (TCE). The 2-square-mile area was placed on the National Priorities List as the North Penn Area 7 Superfund site by the U.S. Environmental Protection Agency (USEPA) in 1989. The U.S. Geological Survey (USGS) conducted geophysical logging, aquifer testing, and water-level monitoring, and measured streamflows in and near North Penn Area 7 from fall 2000 through fall 2006 in a technical assistance study for the USEPA to develop an understanding of the hydrogeologic framework in the area as part of the USEPA Remedial Investigation. In addition, the USGS developed a groundwater-flow computer model based on the hydrogeologic framework to simulate regional groundwater flow and to estimate directions of groundwater flow and pathways of groundwater contaminants. The study area is underlain by Triassic- and Jurassic-age sandstones and shales of the Lockatong Formation and Brunswick Group in the Mesozoic Newark Basin. Regionally, these rocks strike northeast and dip to the northwest. The sequence of rocks form a fractured-sedimentary-rock aquifer that acts as a set of confined to partially confined layers of differing permeabilities. Depth to competent bedrock typically is less than 20 ft below land surface. The aquifer layers are recharged locally by precipitation and discharge locally to streams. The general configuration of the potentiometric surface in the aquifer is similar to topography, except in areas affected by pumping. The headwaters of Wissahickon Creek are nearby, and the stream flows southwest, parallel to strike, to bisect North Penn Area 7. Groundwater is pumped in the vicinity of North Penn Area 7 for industrial use, public supply, and residential supply. Results of field investigations by USGS at the site and results from other studies support, and are consistent with, a conceptual model of a layered leaky aquifer where the dip of the beds has a strong control on hydraulic connections in the groundwater system. Connections within and (or) parallel to bedding tend to be greater than across bedding. Transmissivities of aquifer intervals isolated by packers ranged over three orders of magnitude [from about 2.8 to 2,290 square feet per day (ft 2 /d) or 0.26 to 213 square meters per day (m 2 /d)], did not appear to differ much by mapped geologic unit, but showed some relation to depth being relatively smaller in the shallowest and deepest intervals (0 to 50 ft and more than 250 ft below land surface, respectively) compared to the intermediate depth intervals (50 to 250 ft below land surface) tested. Transmissivities estimated from multiple-observation well aquifer tests ranged from about 700 to 2,300 ft 2 /d (65 to 214 m 2 /d). Results of chemical analyses of water from isolated intervals or monitoring wells open to short sections of the aquifer show vertical differences in concentrations; chloride and silica concentrations generally were greater in shallow intervals than in deeper intervals. Chloride concentrations greater than 100 milligrams per liter (mg/L), combined with distinctive chloride/bromide ratios, indicate a different source of chloride in the western part of North Penn Area 7 than elsewhere in the site. Groundwater flow at a regional scale under steady-state conditions was simulated by use of a numerical model (MODFLOW-2000) for North Penn Area 7 with different layers representing saprolite/highly weathered rock near the surface and unweathered competent bedrock. The sedimentary formations that underlie the study area were modeled using dipping model layers for intermediate and deep zones of unweathered, fractured rock. Horizontal cell model size was 100 meters (m) by 100 meters (328 ft by 328 ft), and model layer thickness ranged from 6 m (19.7 ft) representing shallow weathered rock and saprolite up to 200 m (656 ft) representing deeper dipping bedrock. The model did not include detailed structure to account for local-scale differences in hydraulic properties, with the result that local-scale groundwater flow may not be well simulated. Additional detailed multi-well aquifer tests would be needed to establish the extent of interconnection between intervals at the local scale to address remediation of contamination at each source area. This regional groundwater-flow model was calibrated against measured groundwater levels (1996, 2000, and 2005) and base flow estimated from selected streamflow measurements by use of nonlinear-regression parameter-estimation algorithms to determine hydraulic conductivity and anisotropy of hydraulic conductivity, streambed hydraulic conductivity, and recharge during calibration periods. Results of the simulation using the calibrated regional model indicate that the aquifer appears to be anisotropic where hydraulic conductivity is greatest parallel to the orientation of bedding of the formations underlying the area and least in the cross-bed direction. The maximum hydraulic conductivity is aligned with the average regional strike of the formations, which is “subhorizontal” in the model because the altitudes of the beds and model cells vary in the strike, as well as dip, direction. Estimated subhorizontal hydraulic conductivities (in strike direction parallel to dipping beds) range from 0.001 to 1.67 meters per day (0.0032 to 5.5 feet per day). The ratio of minimum (dip direction) to maximum (strike direction) subhorizontal hydraulic conductivity ranges from 1/3.1 to 1/8.6, and the ratio of vertical to horizontal hydraulic conductivity ranges from 1/1 to 1/478. However, limited available field data precluded rigorous calibration of vertical anisotropy in the model. Estimated recharge rates corresponding to calibration periods in 1996, 2000, and 2005 are 150, 109, and 124 millimeters per year (5.9, 4.3, and 4.9 inches per year), respectively. The calibrated groundwater-flow model was used to simulate groundwater flow under steady-state conditions during periods of relatively high withdrawals (pumpage) (1990) and relatively low withdrawals (2000 and 2005). Groundwater-flow paths originating from recharge areas near known areas of soil contamination (sources) were simulated. Pumped industrial and production wells captured more groundwater from several of these sources during 1990 than after 1990 when pumping declined or ceased and greater amounts of contaminated groundwater moved away from North Penn Area 7 Superfund site to surrounding areas. Uncertainty in simulated groundwater-flow paths from contaminant sources and contributing areas, resulting from uncertainty in estimated hydraulic properties of the model, was illustrated through Monte Carlo simulations. The effect of uncertainty in the vertical anisotropy was not included in the Monte Carlo simulations. Contributing areas indicating the general configuration of groundwater flow towards production well MG-202 (L-22) in the study area also were simulated for the different time periods; as simulated, the flow paths do not pass through any identified contaminant source in North Penn Area 7. However, contributing areas to wells, such as MG-202, located near many pumped wells are particularly complex and, in some cases, include areas that contribute flow to streams that subsequently recharge the aquifer through stream loss. In these cases, water-quality constituents, including contaminants that are present in surface water may be drawn into the aquifer to nearby pumped wells. Results of a simulated shutdown of well MG-202 under steady-state 2005 conditions showed that the area contributing recharge for nearby production well MG-76 (L-17), when MG-202 is not pumping, shifts downstream and is similar to the area contributing recharge for MG-202 when both wells are pumping. Concentrations of constituents in groundwater samples collected in fall 2005 or spring 2006 were compared to simulated groundwater-flow paths for the year 2005 to provide a qualitative assessment of model results. The observed spatial distribution of selected constituents, including TCE, CFC-11, and CFC-113 in groundwater in 2005 and the chloride/bromide mass ratios in 2006, generally were consistent with the model results of the simulated 2005 groundwater-flow paths at North Penn Area 7, indicating the presence of several separate sources of contaminants within North Penn Area 7.

Pennsylvania↗

USGS2014

Simulation of the effects of rainfall and groundwater use on historical lake water levels, groundwater levels, and spring flows in central Florida

The urbanization of central Florida has progressed substantially in recent decades, and the total population in Lake, Orange, Osceola, Polk, and Seminole Counties more than quadrupled from 1960 to 2010. The Floridan aquifer system is the primary source of water for potable, industrial, and agricultural purposes in central Florida. Despite increases in groundwater withdrawals to meet the demand of population growth, recharge derived by infiltration of rainfall in the well-drained karst terrain of central Florida is the largest component of the long-term water balance of the Floridan aquifer system. To complement existing physics-based groundwater flow models, artificial neural networks and other data-mining techniques were used to simulate historical lake water level, groundwater level, and spring flow at sites throughout the area. Historical data were examined using descriptive statistics, cluster analysis, and other exploratory analysis techniques to assess their suitability for more intensive data-mining analysis. Linear trend analyses of meteorological data collected by the National Oceanic and Atmospheric Administration at 21 sites indicate 67 percent of sites exhibited upward trends in air temperature over at least a 45-year period of record, whereas 76 percent exhibited downward trends in rainfall over at least a 95-year period of record. Likewise, linear trend analyses of hydrologic response data, which have varied periods of record ranging in length from 10 to 79 years, indicate that water levels in lakes (307 sites) were about evenly split between upward and downward trends, whereas water levels in 69 percent of wells (out of 455 sites) and flows in 68 percent of springs (out of 19 sites) exhibited downward trends. Total groundwater use in the study area increased from about 250 million gallons per day (Mgal/d) in 1958 to about 590 Mgal/d in 1980 and remained relatively stable from 1981 to 2008, with a minimum of 559 Mgal/d in 1994 and a maximum of 773 Mgal/d in 2000. The change in groundwater-use trend in the early 1980s and the following period of relatively slight trend is attributable to the concomitant effects of increasing public-supply withdrawals and decreasing use of water by the phosphate industry and agriculture. On the basis of available historical data and exploratory analyses, empirical lake water-level, groundwater-level, and spring-flow models were developed for 22 lakes, 23 wells, and 6 springs. Input time series consisting of various frequencies and frequency-band components of daily rainfall (1942 to 2008) and monthly total groundwater use (1957 to 2008) resulted in hybrid signal-decomposition artificial neural network models. The final models explained much of the variability in observed hydrologic data, with 43 of the 51 sites having coefficients of determination exceeding 0.6, and the models matched the magnitude of the observed data reasonably well, such that models for 32 of the 51 sites had root-mean-square errors less than 10 percent of the measured range of the data. The Central Florida Artificial Neural Network Decision Support System was developed to integrate historical databases and the 102 site-specific artificial neural network models, model controls, and model output into a spreadsheet application with a graphical user interface that allows the user to simulate scenarios of interest. Overall, the data-mining analyses indicate that the Floridan aquifer system in central Florida is a highly conductive, dynamic, open system that is strongly influenced by external forcing. The most important external forcing appears to be rainfall, which explains much of the multiyear cyclic variability and long-term downward trends observed in lake water levels, groundwater levels, and spring flows. For most sites, groundwater use explains less of the observed variability in water levels and flows than rainfall. Relative groundwater-use impacts are greater during droughts, however, and long-term trends in water levels and flows were identified that are consistent with historical groundwater-use patterns. The sensitivity of the hydrologic system to rainfall is expected, owing to the well-drained karst terrain and relatively thin confinement of the Floridan aquifer system in much of central Florida. These characteristics facilitate the relatively rapid transmission of infiltrating water from rainfall to the water table and contribute to downward leakage of water to the Floridan aquifer system. The areally distributed nature of rainfall, as opposed to the site-specific nature of groundwater use, and the generally high transmissivity and low storativity properties of the semiconfined Floridan aquifer system contribute to the prevalence of water-level and flow patterns that mimic rainfall patterns. In general, the data-mining analyses demonstrate that the hydrologic system in central Florida is affected by groundwater use differently during wet periods, when little or no system storage is available (high water levels), compared to dry periods, when there is excess system storage (low water levels). Thus, by driving the overall behavior of the system, rainfall indirectly influences the degree to which groundwater use will effect persistent trends in water levels and flows, with groundwater-use impacts more prevalent during periods of low water levels and spring flows caused by low rainfall and less prevalent during periods of high water levels and spring flows caused by high rainfall. Differences in the magnitudes of rainfall and groundwater use during wet and dry periods also are important determinants of hydrologic response. An important implication of the data-mining analyses is that rainfall variability at subannual to multidecadal timescales must be considered in combination with groundwater use to provide robust system-response predictions that enhance sustainable resource management in an open karst aquifer system. The data-driven approach was limited, however, by the confounding effects of correlation between rainfall and groundwater use, the quality and completeness of the historical databases, and the spatial variations in groundwater use. The data-mining analyses indicate that available historical data when used alone do not contain sufficient information to definitively quantify the related individual effects of rainfall and groundwater use on hydrologic response. The knowledge gained from data-driven modeling and the results from physics-based modeling, when compared and used in combination, can yield a more comprehensive assessment and a more robust understanding of the hydrologic system than either of the approaches used separately.

Florida↗

USGS2018

Status of groundwater-level altitudes and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers, Houston-Galveston region, Texas, 2018

Since the early 1900s, most of the groundwater withdrawals in the Houston-Galveston region, Texas, have been from the three primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers. Withdrawals from these aquifers are used for municipal supply, industrial, and irrigation purposes. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of groundwater-level altitudes and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers in the Houston-Galveston region. This report contains regional-scale maps depicting approximate 2018 groundwater-level altitudes (represented by measurements made during December 2017 through March 2018) and long-term groundwater-level changes for the Chicot, Evangeline, and Jasper aquifers. In 2018, groundwater-level-altitude contours for the Chicot aquifer ranged from 200 feet (ft) below the North American Vertical Datum of 1988 (hereinafter referred to as “datum”) to 200 ft above datum. The 1977–2018 groundwater-level-change contours for the Chicot aquifer depict a large area of decline in groundwater-level altitudes (120 ft) in northwestern Harris County. The largest rise in groundwater-level altitudes in the Chicot aquifer from 1977 to 2018 (180 ft) was in southeastern Harris County. Groundwater-level-altitude contours for the Evangeline aquifer ranged from 250 ft below datum to 200 ft above datum in 2018. The 1977–2018 groundwater-level-change contours for the Evangeline aquifer depict broad areas where groundwater-level altitudes either declined or rose. The largest decline in groundwater-level altitudes (320 ft) was in southern Montgomery County. The largest rise in groundwater-level altitudes in the Evangeline aquifer from 1977 to 2018 (240 ft) was in southeastern Harris County. In 2018, groundwater-level-altitude contours for the Jasper aquifer ranged from 200 ft below datum to 200 ft above datum. The 2000–18 groundwater-level-change contours for the Jasper aquifer depict groundwater-level declines throughout most of the study area where groundwater-level-altitude data from the Jasper aquifer were collected, with the largest decline (200 ft) in southern Montgomery County.

Texas↗

USGS2019

Status of groundwater-level altitudes and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers, Houston-Galveston region, Texas, 2019

Since the early 1900s, most of the groundwater withdrawals in the Houston-Galveston region, Texas, have been from the three primary aquifers that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers. Withdrawals from these aquifers are used for municipal supply, commercial and industrial use, and irrigation. This report, prepared by the U.S. Geological Survey in cooperation with the Harris-Galveston Subsidence District, City of Houston, Fort Bend Subsidence District, Lone Star Groundwater Conservation District, and Brazoria County Groundwater Conservation District, is one in an annual series of reports depicting the status of groundwater-level altitudes and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers in the Houston-Galveston region. This report contains regional-scale maps depicting approximate 2019 groundwater-level altitudes (represented by measurements made during December 2018 through March 2019) and long-term groundwater-level changes in the Chicot, Evangeline, and Jasper aquifers. In 2019, groundwater-level-altitude contours for the Chicot aquifer ranged from 200 feet (ft) below the North American Vertical Datum of 1988 (hereinafter referred to as “datum”) to 200 ft above datum. The 1977–2019 groundwater-level-change contours for the Chicot aquifer depict a large area of decline in groundwater-level altitudes (100 ft) in northwestern Harris County. The largest rise in groundwater-level altitudes in the Chicot aquifer from 1977 to 2019 (200 ft) was in southeastern Harris County. In 2019, groundwater-level-altitude contours for the Evangeline aquifer ranged from 300 ft below datum to 200 ft above datum. The 1977–2019 groundwater-level-change contours for the Evangeline aquifer depict broad areas where groundwater-level altitudes either declined or rose. The largest decline in groundwater-level altitudes (280 ft) was in southern Montgomery and northern Harris Counties. The largest rise in groundwater-level altitudes in the Evangeline aquifer from 1977 to 2019 (240 ft) was in southeastern Harris County. In 2019, groundwater-level-altitude contours for the Jasper aquifer ranged from 200 ft below datum to 250 ft above datum. The 2000–19 groundwater-level-change contours for the Jasper aquifer depict groundwater-level declines throughout most of the study area where groundwater-level-altitude data from the Jasper aquifer were collected, with the largest decline (200 ft) in southern Montgomery County.

Texas↗

USGS2019

Groundwater-flow model and analysis of groundwater and surface-water interactions for the Big Sioux aquifer, Sioux Falls, South Dakota

The city of Sioux Falls, in southeastern South Dakota, is the largest city in South Dakota. The U.S. Geological Survey (USGS), in cooperation with the city of Sioux Falls, completed a groundwater-flow model to use for improving the understanding of groundwater-flow processes, estimating hydrogeologic properties, and analyzing groundwater and surface-water interactions for the Big Sioux aquifer in the model area. The model area includes the Big Sioux aquifer and the underlying hydrogeologic units from Dell Rapids, South Dakota, to the confluence of the Big Sioux River and the outlet of the Sioux Falls Diversion Channel in eastern Sioux Falls, S. Dak. The Big Sioux aquifer is the primary aquifer in the model area and the focus of the groundwater-flow model. The Big Sioux River is the largest stream in the model area and is in hydraulic connection with the Big Sioux aquifer. A conceptual model for the area was constructed and includes a characterization of the hydrogeologic framework, analysis and construction of potentiometric surfaces, and summary of estimated water budget components in the model area. The primary hydrogeologic units in the model area consist of (1) the Big Sioux aquifer, (2) a glacial till confining unit, and (3) bedrock aquifers (Split Rock Creek and Sioux Quartzite aquifers). Sources of groundwater recharge included infiltration of precipitation, stream seepage, and groundwater exchanges among the hydraulically connected Big Sioux aquifer, glacial till confining unit, and bedrock aquifers. Groundwater losses included evapotranspiration, groundwater discharge to streams, and groundwater withdrawal to supply water-use needs. A numerical groundwater-flow model (numerical model) was constructed and was used to simulate all aspects of the conceptual model for predevelopment (steady-state) and time-varying (transient) monthly conditions for 1950–2017. The numerical model was constructed using the USGS modular hydrologic simulation program, MODFLOW–6, and was calibrated using the Parameter ESTimation software, PEST++. The transient numerical model was calibrated for steady-state and transient monthly conditions for 1950–2017. Calibration targets were observations of hydraulic head, changes in hydraulic head, monthly mean streamflow (as a rate), and cumulative monthly stream discharge (as a volume). Parameters adjusted during model calibration were horizontal and vertical hydraulic conductivity, specific storage, specific yield, recharge and evapotranspiration multipliers, and streambed hydraulic conductivity. Horizontal and vertical hydraulic conductivity were estimated at pilot points distributed within the model area; specific storage and specific yield were assigned to uniform values in each layer in the model area; recharge and evapotranspiration multipliers were assigned uniformly for every stress period in the numerical model; and streambed hydraulic conductivity values were assigned uniformly between stream confluences. The final calibrated parameter values of horizontal and vertical hydraulic conductivity, specific yield, specific storage, streambed hydraulic conductivity, recharge, and evapotranspiration were considered reasonable for the hydrogeologic materials and conditions in the model area for 1950–2017. Overall, simulated hydraulic head altitudes had a linear regression coefficient of determination (R 2 ) of 0.48. Hydraulic head altitude residuals for the glacial till confining unit and bedrock aquifers were typically greater in magnitude when compared to residuals in the Big Sioux aquifer, but simulated hydraulic head altitudes in the Big Sioux aquifer compared favorably with mean observed hydraulic head altitudes and had a linear regression R 2 of 0.93. Simulated streamflow hydrographs matched the general trends of observed increases and decreases in streamflow for USGS streamgages 06482000 (Big Sioux River at Sioux Falls, S. Dak.) and 06482020 (Big Sioux River at North Cliff Avenue at Sioux Falls, S. Dak.), but larger streamflows were overestimated at the first streamgage and underestimated at the second streamgage. The numerical model reasonably estimated cumulative monthly stream discharge for the first 10–15 years of available streamflow records at both USGS streamgages. After the first 10–15 years of available streamflow record, cumulative monthly stream discharge was closely estimated for USGS streamgage 06482000 and underestimated at USGS streamgage 06482020. Composite sensitivities without regularization were calculated by PEST++ for the calibrated numerical model parameters and were averaged by parameter group. The parameter group with the highest mean composite sensitivity was the recharge multiplier parameter group. Model simplifications, assumptions, and limitations were necessary for construction of the conceptual and numerical models and for calibration efficiency. Spatial simplification of hydraulic properties could cause the numerical model to misrepresent reactions to changes in localized stresses, such as additional demands for groundwater withdrawal. The numerical model was temporally discretized into monthly periods and required scaling daily rates into representative monthly rates for model input and calibration targets. Based on the comparison between the observed and simulated groundwater levels, monthly mean streamflow and cumulative monthly stream discharge, and general groundwater distribution and flow, the numerical model favorably simulated the flow in the Big Sioux aquifer. Eventual capture was calculated in the model area using a steady-state numerical groundwater-flow model. The eventual capture map shows areas of higher streamflow capture adjacent to the Big Sioux River north of the city of Sioux Falls and along the lower part of the Sioux Falls Diversion Channel, and areas of lower streamflow capture along aquifer boundaries and near the southern Sioux Quartzite barrier. The timing of capture was determined using a transient numerical groundwater-flow model to determine the likely captured water sources for 30 years of groundwater withdrawal at three hypothetical wells using three continuous withdrawal rates (112.5, 450.0, and 900.0 gallons per minute). Supply for all three hypothetical wells became capture-dominated after only a short period of continuous withdrawal. Capture stabilized after about 10–15 years for well A, and after 20–25 years for well B, and after about 10–15 years for well C. The groundwater-flow model is a suitable tool to use for improving the understanding of groundwater-flow processes, estimating hydrogeologic properties, and analyzing groundwater and surface-water interactions for the Big Sioux aquifer near Sioux Falls, S. Dak. The numerical model can be used to simulate hydrologic scenarios, advance understanding of groundwater budgets, compute system response to stress, and determine likely sources of water supplied to wells.

South Dakota↗