USGS Science⌕ Search

Geology topics

Groundwater & aquifers

Find USGS groundwater and aquifer publications. Search an aquifer name or place, compare source-supplied dates and open reports for hydrogeology, recharge and water-level studies.

At least 865 records · Page 48Linked to original sources

Measuring basin-scale aquifer storage change and mapping specific yield in Albuquerque, New Mexico, USA, with repeat microgravity data

Study Region The groundwater basin underlying the City of Albuquerque, New Mexico, USA. Study Focus The study focuses on changes in groundwater storage and how those changes relate to groundwater-level changes. Groundwater storage change was measured using repeat microgravity at 35 stations from 2016 to 2021. Usually, storage is monitored by converting groundwater-level changes to storage changes using the aquifer storage coefficient, a difficult property to measure. With gravity, storage change can be measured directly. The storage coefficient, or specific yield in an unconfined aquifer, was estimated using the gravity method at individual sites, and a map created showing how this property varies over the region. New Hydrological Insights for the Region For the first time, a map of specific yield was produced based on interpolated maps of gravity-derived storage change and groundwater-level change, allowing inference of this property over a broad region even without collocated monitoring wells and gravity stations. Gravity data indicate aquifer drawdown and recovery mainly occurs in the central part of the basin where pumping is greatest. This spatial distribution is not captured by the more widely spaced monitoring well network. Because the aquifer is recovering from historically greater-magnitude pumping, net-neutral storage change (inflow=outflow) occurs when pumping from the central part of the well field is about 1.54 × 10 7 cubic meters (12,500 acre-feet) per year.

New Mexico↗

Application of multiple isotopic and geochemical tracers for investigation of recharge, salinization, and residence time of water in the Souss-Massa aquifer, southwest of Morocco

Groundwater and surface water in Souss-Massa basin in the west-southern part of Morocco is characterized by a large variation in salinity, up to levels of 37 g L-1. The high salinity coupled with groundwater level decline pose serious problems for current irrigation and domestic water supplies as well as future exploitation. A combined hydrogeologic and isotopic investigation using several chemical and isotopic tracers such as Br/Cl, ??18O, ??2H, 3H, 87Sr/86Sr, ??11B, and 14C was carried out in order to determine the sources of water recharge to the aquifer, the origin of salinity, and the residence time of water. Stable isotope, 3H and 14C data indicate that the high Atlas mountains in the northern margin of the Souss-Massa basin with high rainfall and low ??18O and ??2H values (-6 to -8??? and -36 to -50???) is currently constitute the major source of recharge to the Souss-Massa shallow aquifer, particularly along the eastern part of the basin. Localized stable isotope enrichments offset meteoric isotopic signature and are associated with high nitrate concentrations, which infer water recycling via water agricultural return flows. The 3H and 14C data suggest that the residence time of water in the western part of the basin is in the order of several thousands of years; hence old water is mined, particularly in the coastal areas. The multiple isotope analyses and chemical tracing of groundwater from the basin reveal that seawater intrusion is just one of multiple salinity sources that affect the quality of groundwater in the Souss-Massa aquifer. We differentiate between modern seawater intrusion, salinization by remnants of seawater entrapped in the middle Souss plains, recharge of nitrate-rich agricultural return flow, and dissolution of evaporate rocks (gypsum and halite minerals) along the outcrops of the high Atlas mountains. The data generated in this study provide the framework for a comprehensive management plan in which water exploitation should shift toward the eastern part of the basin where current recharge occurs with young and high quality groundwater. In contrast, we argued that the heavily exploited aquifer along the coastal areas is more vulnerable given the relatively longer residence time of the water and salinization processes in this part of the aquifer. ?? 2008 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

Metabolic adaptation and in situ attenuation of chlorinated ethenes by naturally occurring microorganisms in a fractured dolomite aquifer near Niagara Falls, New York

A combination of hydrogeological, geochemical, and microbiological methods was used to document the biotransformation of trichloroethene (TCE) to ethene, a completely dechlorinated and environmentally benign compound, by naturally occurring microorganisms within a fractured dolomite aquifer. Analyses of groundwater samples showed that three microbially produced TCE breakdown products (cis-1,2-dichloroethene, vinyl chloride, and ethene) were present in the contaminant plume. Hydrogen (H2) concentrations in groundwater indicated that iron reduction was the predominant terminal electron-accepting process in the most contaminated geologic zone of the site. Laboratory microcosms prepared with groundwater demonstrated complete sequential dechlorination of TCE to ethene. Microcosm assays also revealed that reductive dechlorination activity was present in waters from the center but not from the periphery of the contaminant plume. This dechlorination activity indicated that naturally occurring microorganisms have adapted to utilize chlorinated ethenes and suggested that dehalorespiring rather than cometabolic, microbial processes were the cause of the dechlorination. The addition of pulverized dolomite to microcosms enhanced the rate of reductive dechlorination, suggesting that hydrocarbons in the dolomite aquifer may serve as electron donors to drive microbially mediated reductive dechlorination reactions. Biodegradation of the chlorinated ethenes appears to contribute significantly to decontamination of the site.A combination of hydrogeological, geochemical, and microbiological methods was used to document the biotransformation of trichloroethene (TCE) to ethene, a completely dechlorinated and environmentally benign compound, by naturally occurring microorganisms within a fractured dolomite aquifer. Analyses of groundwater samples showed that three microbially produced TCE breakdown products (cis-1,2-dichloroethene, vinyl chloride, and ethene) were present in the contaminant plume. Hydrogen (H2) concentrations in groundwater indicated that iron reduction was the predominant terminal electron-accepting process in the most contaminated geologic zone of the site. Laboratory microcosms prepared with groundwater demonstrated complete sequential dechlorination of TCE to ethene. Microcosm assays also revealed that reductive dechlorination activity was present in waters from the center but not from the periphery of the contaminant plume. This dechlorination activity indicated that naturally occurring microorganisms have adapted to utilize chlorinated ethenes and suggested that dehalorespiring rather than cometabolic, microbial processes were the cause of the dechlorination. The addition of pulverized dolomite to microcosms enhanced the rate of reductive dechlorination, suggesting that hydrocarbons in the dolomite aquifer may serve as electron donors to drive microbially mediated reductive dechlorination reactions. Biodegradation of the chlorinated ethenes appears to contribute significantly to decontamination of the site.

Environmental Science & Technology↗

Radiogenic 4He as a conservative tracer in buried‐valley aquifers

The accumulation of 4 He in groundwater can be a powerful tool in hydrogeologic investigations. However, the use of 4 He often suffers from disagreement or uncertainty related to in situ and external sources of 4 He. In situ sources are quantified by several methods, while external sources are often treated as calibration parameters in modeling. We present data from direct laboratory measurements of 4 He release from sediments and field data of dissolved 4 He in the Mahomet Aquifer, a well‐studied buried‐valley aquifer in central Illinois. The laboratory‐derived accumulation rates (0.13–0.91 μcm 3 STP kg water −1 yr −1 ) are 1–2 orders of magnitude greater than the accumulation rates based on the U and Th concentrations of the sediments (0.004–0.009 μcm 3 STP kg water −1 yr −1 ). The direct measurement of accumulation rates are more consistent with dissolved concentrations of 4 He in the groundwater. We suggest that the direct measurement method is applicable in a variety of hydrogeologic settings. The patterns of accumulation of 4 He are consistent with the conceptual model of flow in the aquifer based on hydraulic and geochemical evidence and show areas where in situ production and external sources of 4 He are dominant. In the southwestern part of the study area, Ne concentrations are less than atmospheric solubility, indicating gases have been lost from the groundwater. Available evidence indicates that the gases are lost as groundwater passes by pockets of CH 4 in glacial deposits overlying the aquifer. However, the external flux from the underlying bedrock appears to dominate the accumulation of radiogenic 4 He in the aquifer in the southwestern part of the study area, and the loss or gain of helium as groundwater passes through the overlying sediments is minor in comparison.

Water Resources Research↗

Selecting remediation goals by assessing the natural attenuation capacity of groundwater systems

Remediation goals for the source areas of a chlorinated ethene‐contaminated groundwater plume were identified by assessing the natural attenuation capacity of the aquifer system. The redox chemistry of the site indicates that sulfate‐reducing (H 2 ∼ 2 nanomoles [nM]) per liter conditions near the contaminant source grade to Fe(III)‐reducing conditions (H 2 ∼ 0.5 nM) downgradient of the source. Sulfate‐reducing conditions facilitate the initial reduction of perchloroethene (PCE) to trichloroethene (TCE), cis ‐dichloroethene ( cis ‐DCE), and vinyl chloride (VC). Subsequently, the Fe(III)‐reducing conditions drive the oxidation of cis ‐DCE and VC to carbon dioxide and chloride. This sequence gives the aquifer a substantial capacity for biodegrading chlorinated ethenes. Natural attenuation capacity (the slope of the steady‐state contaminant concentration profile along a groundwater flowpath) is a function of biodegradation rates, aquifer dispersive characteristics, and groundwater flow velocity. The natural attenuation capacity at the Kings Bay, Georgia site was assessed by estimating groundwater flowrates (∼0.23 ± 0.12 m/d) and aquifer dispersivity (∼1 m) from hydrologic and scale considerations. Apparent biodegradation rate constants (PCE and TCE ∼ 0.01 d −1 ; cis ‐DCE and VC ∼ 0.025 d −1 ) were estimated from observed contaminant concentration changes along aquifer flowpaths. A boundary‐value problem approach was used to estimate levels to which contaminant concentrations in the source areas must be lowered (by engineered removal), or groundwater flow velocities lowered (by pumping) for the natural attenuation capacity to achieve maximum concentration limits (MCLs) prior to reaching a predetermined regulatory point of compliance.

Bioremediation Journal↗

Effects of climate and irrigation on GRACE-based estimates of water storage changes in major US aquifers

Understanding climate and human impacts on water storage is critical for sustainable water-resources management. Here we assessed climate and human drivers of total water storage (TWS) variability from Gravity Recovery and Climate Experiment (GRACE) satellites compared with drought severity and irrigation water use in 14 major aquifers in the United States. Results show that long-term variability in TWS tracked by GRACE satellites is dominated by interannual variability in most of the 14 major US aquifers. Low TWS trends in the humid eastern U.S. are linked to low drought intensity. Although irrigation pumpage in the humid Mississippi Embayment aquifer exceeded that in the semi-arid California Central Valley, a surprising lack of TWS depletion in the Mississippi Embayment aquifer is attributed to extensive streamflow capture. Marked storage depletion in the semi-arid southwestern Central Valley and south-central High Plains totaled ∼90 km 3 , about three times greater than the capacity of Lake Mead, the largest U.S. reservoir. Depletion in the Central Valley was driven by long-term droughts (⩽5 yr) amplified by switching from mostly surface water to groundwater irrigation. Low or slightly rising TWS trends in the northwestern (Columbia and Snake Basins) US are attributed to dampening drought impacts by mostly surface water irrigation. GRACE satellite data highlight synergies between climate and irrigation, resulting in little impact on TWS in the humid east, amplified TWS depletion in the semi-arid southwest and southcentral US, and dampened TWS deletion in the northwest and north central US Sustainable groundwater management benefits from conjunctive use of surface water and groundwater, inefficient surface water irrigation promoting groundwater recharge, efficient groundwater irrigation minimizing depletion, and increasing managed aquifer recharge. This study has important implications for sustainable water development in many regions globally.

Environmental Research Letters↗

Nitrate in aquifers beneath agricultural systems

Research from several regions of the world provides spatially anecdotal evidence to hypothesize which hydrologic and agricultural factors contribute to groundwater vulnerability to nitrate contamination. Analysis of nationally consistent measurements from the U.S. Geological Survey's NAWQA program confirms these hypotheses for a substantial range of agricultural systems. Shallow unconfined aquifers are most susceptible to nitrate contamination associated with agricultural systems. Alluvial and other unconsolidated aquifers are the most vulnerable and also shallow carbonate aquifers that provide a substantial but smaller contamination risk. Where any of these aquifers are overlain by permeable soils the risk of contamination is larger. Irrigated systems can compound this vulnerability by increasing leaching facilitated by additional recharge and additional nutrient applications. The system of corn, soybean, and hogs produced significantly larger concentrations of groundwater nitrate than all other agricultural systems because this system imports the largest amount of N-fertilizer per unit production area. Mean nitrate under dairy, poultry, horticulture, and cattle and grains systems were similar. If trends in the relation between increased fertilizer use and groundwater nitrate in the United States are repeated in other regions of the world, Asia may experience increasing problems because of recent increases in fertilizer use. Groundwater monitoring in Western and Eastern Europe as well as Russia over the next decade may provide data to determine if the trend in increased nitrate contamination can be reversed. If the concentrated livestock trend in the United States is global, it may be accompanied by increasing nitrogen contamination in groundwater. Concentrated livestock provide both point sources in the confinement area and intense non-point sources as fields close to facilities are used for manure disposal. Regions where irrigated cropland is expanding, such as in Asia, may experience the greatest impact of this practice on groundwater nitrate. ?? USDA 2007.

Conference Paper↗

Groundwater geochemical and selected volatile organic compound data, Operable Unit 1, Naval Undersea Warfare Center, Division Keyport, Washington, June 2009

Previous investigations indicate that natural attenuation and biodegradation of chlorinated volatile organic compounds (VOCs) are substantial in groundwater beneath the 9-acre former landfill at Operable Unit 1 (OU 1), Naval Undersea Warfare Center, Division Keyport, Washington. Phytoremediation combined with ongoing natural attenuation processes was the preferred remedy selected by the U.S. Navy, as specified in the Record of Decision for the site. The U.S. Navy planted two hybrid poplar plantations on the landfill in spring 1999 to remove and to control the migration of chlorinated VOCs in shallow groundwater. The U.S. Geological Survey (USGS) has continued to monitor groundwater geochemistry to ensure that conditions remain favorable for contaminant biodegradation as specified in the Record of Decision. This report presents groundwater geochemical and selected VOC data collected at OU 1 by the USGS during June 15-17, 2009, in support of long-term monitoring for natural attenuation. For 2009, groundwater samples were collected from 13 wells and 9 piezometers. Samples from all wells and piezometers were analyzed for redox sensitive constituents, and samples from 10 of 18 upper-aquifer wells and piezometers and 3 of 4 intermediate-aquifer wells also were analyzed for chlorinated VOCs. Concentrations of redox sensitive constituents measured in 2009 were consistent with previous years, with dissolved hydrogen (H2) concentrations ranging from less than 0.1 to 1.8 nanomolar (nM), dissolved oxygen concentrations all at 0.6 milligram per liter or less; little to no detectable nitrate; abundant dissolved manganese, iron, and methane; and commonly detected sulfide. The reductive declorination byproducts-methane, ethane, and ethene-were not detected in samples collected from the upgradient wells in the landfill or the upper aquifer beneath the northern phytoremediation plantation. Chlorinated VOC concentrations in 2009 at most piezometers were similar to or slightly less than chlorinated VOC concentrations measured in previous years. In 2009, concentrations of reductive dechlorination byproducts ethane and ethene were less than those measured in 2008 at most northern plantation wells and piezometers. For the upper aquifer beneath the southern phytoremediation plantation, chlorinated VOC concentrations in 2009 at the piezometers were extremely high and continued to vary considerably over space and between years. At piezometer P1-9, the total chlorinated VOC concentration increased from 25,000 micrograms per liter in 2008 to more than 172,000 micrograms per liter in 2009. At piezometer P1-7 in 2009, the concentrations of trichloroethene and cis-1,2-dichloroethene (cis-DCE) were the highest to date. The reductive dechlorination byproducts ethane and ethene were detected at all wells and piezometers in the southern plantation with the exception of piezometer P1-8, although the measured concentrations were not consistently high. For the intermediate aquifer, concentrations of redox sensitive constituents and VOCs in 2009 at wells MW1-25, MW1-28, and MW1-39 were consistent with concentrations measured in previous years. Concentrations of the reductive dechlorination byproducts ethane and ethene at wells MW1-25 and MW1-28 were equal to or greater than previously measured concentrations.

Washington↗

Hydraulic, geochemical, and thermal monitoring of an aquifer system in the vicinity of Mammoth Lakes, Mono County, California, 2015–17

Since 2014, the U.S. Geological Survey has been working in cooperation with the Bureau of Land Management, Mono County, Ormat Technologies, Inc., and the Mammoth Community Water District to design and implement a groundwater-monitoring program for the proposed Casa Diablo IV Geothermal Power Project in Long Valley Caldera, California, to characterize baseline groundwater-level, water-temperature, and water-chemistry conditions at dedicated monitoring wells and municipal supply wells. The publicly available data and the analyses provided here represent quality-assured and peer-reviewed information to help with the management of the thermal and non-thermal water resources beneath and in the vicinity of the town of Mammoth Lakes, California. The methods of data collection for continuous water levels and quarterly water-temperature profiles for two 600-foot-deep monitoring wells during 2016 through 2017 are discussed. Also discussed are the methods of water-sample collection and characterizations of the water chemistry in numerous wells in the multilayered aquifer system beneath Mammoth Lakes. Additionally, the methodology used to develop digital (mathematical) filters to remove or reduce the effects of barometric pressure and solid Earth tides on the continuous water-level records is discussed. Digitally filtered water levels for a 2017 flow test of a deep geothermal production well are described, and various aquifer responses observed during the flow test are discussed. These are further considered in a companion evaluation of potential physical and chemical influences on the water-level data collected during the flow test. The digitally filtered water-level data indicated that some hydraulic communication exists between the deep geothermal aquifer and shallow groundwater aquifer at the location of the flow test, northeast of Mammoth Lakes. Groundwater-chemistry data from three wells indicated that shallow groundwater naturally mixes with a small component of geothermal water along the northern periphery of the shallow aquifer system at Mammoth Lakes.

California↗

Transmissivity of the Upper Floridan aquifer in Florida and parts of Georgia, South Carolina, and Alabama

The Floridan aquifer system (FAS) covers an area of approximately 100,000 square miles in Florida and parts of Georgia, South Carolina, Alabama, and Mississippi. Groundwater wells for water supply were first drilled in the late 1800s and by the year 2000, the FAS was the primary source of drinking water for about 10 million people. One of the methods for assessing groundwater availability is the development of regional or subregional groundwater flow models of the aquifer system that can be used to develop water budgets spatially and temporally, as well as evaluate the groundwater resource change over time. Understanding the distribution of transmissivity within the FAS is critical to the development of groundwater flow models. The map presented herein differs from previously published maps of the FAS in that it is based on interpolation of 1,487 values of transmissivity. The transmissivity values in the dataset range from 8 to 9,000,000 feet squared per day (ft 2 /d) with the majority of the values ranging from 10,000 to 100,000 ft 2 /d. The wide range in transmissivity (6 orders of magnitude) is typical of carbonate rock aquifers, which are characterized by a wide range in karstification. Commonly, the range in transmissivity is greatest in areas where groundwater flow creates conduits in facies that dissolve more readily or areas of high porosity units that have interconnected vugs, with diameters greater than 0.1 foot. These are also areas where transmissivity is largest. Additionally, first magnitude springsheds and springs are shown because in these springshed areas, the estimates of transmissivity from interpolation may underestimate the actual range in transmissivity. Also shown is an area within the Gulf Trough in Georgia where high yielding wells are unlikely to be developed in the Upper Floridan aquifer. The interpolated transmissivity ranges shown on this map reflect the geologic structure and karstified areas. Transmissivity is large in the areas where the system is unconfined, such as west-central Florida and southwest Georgia just northwest of the Gulf Trough. Transmissivity is small along the Gulf Trough and Southwest Georgia Embayment (referred to as Apalachicola Embayment in some reports). Transmissivity is also small in the thin, updip part of the system near its northern boundary. Another area of large transmissivity coincides with the Southeast Georgia Embayment.

Alabama, Florida, Georgia, South Carolina↗

A water-budget model and assessment of groundwater recharge for the Island of Hawai'i

Concern surrounding increasing demand for groundwater on the Island of Hawaiʻi, caused by a growing population and an increasing reliance on groundwater as a source for municipal and private water systems, has prompted a study of groundwater recharge on the island using the most current data and accepted methods. For this study, a daily water-budget model for the entire Island of Hawaiʻi was developed and used to estimate mean recharge for various land-cover and rainfall conditions, and a submodel for the Kona area was developed and used to estimate historical groundwater recharge in the Kona area during the period 1984–2008. Estimated mean annual recharge on the Island of Hawaiʻi is 6,594 million gallons per day, which is about 49 percent of mean annual rainfall. Recharge is highest on the windward slopes of Mauna Loa, below the tradewind inversion, and lowest on the leeward slopes of Kohala and Mauna Kea. Local recharge maxima also occur on (1) windward Kohala, with the exception of the northern tip, (2) windward Mauna Kea below the tradewind inversion, (3) windward Kīlauea, (4) the middle elevations of southeastern Mauna Loa, and (5) the lower-middle elevations of leeward Mauna Loa and southwestern Hualālai, in the Kona area. Local recharge minima also occur on (1) Mauna Kea and Mauna Loa, above the tradewind inversion, (2) the northern tip of Kohala, (3) leeward Kīlauea, (4) the southern tip of Mauna Loa, and (5) the northwestern slopes of Mauna Loa and Hualālai. In 18 of the 24 aquifer systems on the island, estimated mean annual recharge for baseline conditions was higher than the recharge estimates used in the 2008 State of Hawaiʻi Water Resource Protection Plan (2008 WRPP). Baseline conditions for this study were 2008 land cover and mean annual rainfall from the period 1916–1983. Estimates of recharge for the Māhukona, Waimea, and Hāwī aquifer systems, however, were between 29 and 38 percent lower than the 2008 WRPP estimates, mainly because of much higher evapotranspiration estimates in this study compared to the 2008 WRPP. For the drought simulation (1991–95 rainfall), the estimates of recharge for these three aquifer systems were only 15 to 33 percent of the sustainable yields (maximum allowable pumping rates) set by the 2008 WRPP. This may be cause for concern, as these areas are experiencing a rapid growth in development and a related growth in water demand. Recent projections of change in rainfall owing to effects of ongoing climate change generally indicate a slight increase in islandwide rainfall, and estimates of annual recharge in the late 21st century are higher than baseline estimates for every aquifer system, except ʻAnaehoʻomalu. On average, these aquifer-system recharge estimates are higher by about 8 percent compared to baseline estimates. In the Kona area, estimated groundwater recharge during the period 1984–2008 was highest during 2004–8 and lowest during 1999–2003, with the 1999–2003 recharge being about 50 percent of the 2004–8 recharge. These extremes in recharge coincided with the periods of lowest and highest mean rainfall, respectively. No seasonal pattern in recharge is discernible. Spatially, the highest recharge occurred in a belt about 4 miles wide running parallel to the coast about 2 miles inland. The sensitivity of recharge estimates to input parameters is related to the climate and land-cover conditions of the particular area of study. For the wet, forested areas characteristic of the windward side of the island, recharge was most sensitive to the ratio of runoff to rainfall. For the dry, grassland areas characteristic of the northwestern leeward side of the island, recharge was most sensitive to root depth. For the Kona area, characterized by moderate rainfall and a wide variety of land cover, recharge was most sensitive to the pan coefficient and canopy-evaporation rates in

Hawai'i↗

Updated numerical model with uncertainty assessment of 1950-56 drought conditions on brackish-water movement within the Edwards aquifer, San Antonio, Texas

In 2010, the U.S. Geological Survey, in cooperation with the San Antonio Water System, began a study to assess the brackish-water movement within the Edwards aquifer (more specifically the potential for brackish-water encroachment into wells near the interface between the freshwater and brackish-water transition zones, referred to in this report as the transition-zone interface) and effects on spring discharge at Comal and San Marcos Springs under drought conditions using a numerical model. The quantitative targets of this study are to predict the effects of higher-than-average groundwater withdrawals from wells and drought-of-record rainfall conditions of 1950–56 on (1) dissolved-solids concentration changes at production wells near the transition-zone interface, (2) total spring discharge at Comal and San Marcos Springs, and (3) the groundwater head (head) at Bexar County index well J-17. The predictions of interest, and the parameters implemented into the model, were evaluated to quantify their uncertainty so the results of the predictions could be presented in terms of a 95-percent credible interval. The model area covers the San Antonio and Barton Springs segments of the Edwards aquifer; the history-matching effort was focused on the San Antonio segment. A previously developed diffuse-flow model of the Edwards aquifer, which forms the basis for the model in this assessment, is primarily based on a conceptualization in which flow in the aquifer is predominately through a network of numerous small fractures and openings. Primary updates to this model include an extension of the active area downdip, a conversion to an 8-layer SEAWAT variable-density flow and transport model to simulate dissolved-solids concentration effects on water density, history matching to 1999–2009 conditions, and parameter estimation in a highly parameterized context using automated methods in PEST (a model-independent Parameter ESTimation code). In addition to the best-fit parameter values derived from history matching, the uncertainty of model parameters was also estimated by using linear uncertainty analysis. Comparison of “prior” (before history matching) and “posterior” (after history matching) variances of parameters indicate that the information within the observation dataset used for history matching informs many parameters. The concentration threshold parameters were well-informed by the observation dataset as their posterior distributions were much narrower than their prior distributions. The transition-zone scaling parameters of hydraulic conductivity, effective porosity, and specific storage were all informed by the observation dataset, as evidenced by the difference between the prior and posterior variances. Saline-zone scaling parameters, alternatively, were not informed by the observation dataset for effective porosity and specific storage. Resulting posterior drier-month, wetter-month, and annual recharge multiplier parameter variances are important to understanding how well recharge is estimated and implemented within the model. The shifts of the posterior distributions left and right indicate that there were zones where less or more water was needed in the model. The widths of the distributions were not decreased substantially, indicating that many of the best-fit recharge parameters are not statistically different from the initial values specified in the history-matching effort. Recharge from rainfall is the driving force behind groundwater flow and heads in the aquifer; therefore, an increase in understanding of this process would benefit model development by potentially decreasing the uncertainty of this parameter. The history-matching effort was most helpful in informing the parameters in the model that control discharge at springs, namely, the spring orifice (drain) altitude and drain conductance parameters for each spring. The uncertainty assessment of the predictive model (a hypothetical recurrence of 1950–56 drought conditions and higher-than-average groundwater withdrawals from wells) provided insights into the potential effects of these conditions on dissolved-solids concentration changes at production wells near the transition-zone interface, discharges at Comal and San Marcos Springs, and heads at Bexar County index well J-17. Results at the 25 production wells near the transition-zone interface indicate that the uncertainty of model input parameters based on expert knowledge yielded an upper bound of the 95-percent credible interval of dissolved-solids concentrations that exceeds the secondary drinking water standards of 1,000 milligrams per liter (mg/L) of the Texas Commission on Environmental Quality (TCEQ) for many wells. However, the history-matching process provided key information to inform prediction-sensitive model parameters and therefore, contributed to a substantial decrease of the upper bound of the 95-percent credible interval to below the secondary drinking water standards. Reductions in dissolved-solids concentration changes were on the order of 400 mg/L to 1,300 mg/L. The reduction in uncertainty in regards to this prediction implies that this prediction of dissolved-solids concentration change can be made with some certainty using this current model and that those parameters that control this prediction are informed by the observation dataset. Even though predictive uncertainty was reduced for this prediction, dissolved-solids concentration changes were still greater than zero, indicating a minimal increase in concentration at these 25 production wells during the 7-year simulation period is likely. However, this minimal concentration increase indicates a small potential for movement of the brackish-water transition zone near these wells during the 7-year simulation period of drought-ofrecord (1950–56) rainfall conditions with higher-than-average groundwater withdrawals by wells. Predictive results of total spring discharge during the 7-year period, as well as head predictions at Bexar County index well J-17, were much different than the dissolved-solids concentration change results at the production wells. These upper bounds are an order of magnitude larger than the actual prediction which implies that (1) the predictions of total spring discharge at Comal and San Marcos Springs and head at Bexar County index well J-17 made with this model are not reliable, and (2) parameters that control these predictions are not informed well by the observation dataset during historymatching, even though the history-matching process yielded parameters to reproduce spring discharges and heads at these locations during the history-matching period. Furthermore, because spring discharges at these two springs and heads at Bexar County index well J-17 represent more of a cumulative effect of upstream conditions over a larger distance (and longer time), many more parameters (with their own uncertainties) are potentially controlling these predictions than the prediction of dissolved-solids concentration change at the prediction wells, and therefore contributing to a large posterior uncertainty.

Texas↗

Chemical characteristics of water in the surficial aquifer system, Dade County, Florida

Geohydrologic test drilling was conducted throughout Dade County to describe the chemical characteristics of water from geohydrologic units in the surficial aquifer system. Water quality analysis of samples collected from the test wells completed in central Dade County indicates that the Biscayne aquifer (unit A), the upper clastic unit of the Tamiami Formation (unit B), and parts of a limestone, sandstone, and sand unit of the Tamiami Formation (unit C) have been effectively flushed of residual seawater, and contain primarily calcium bicarbonate groundwater. The lower parts of unit C and the lower clastic unit of the Tamiami Formation (unit D) primarily contain calcium-sodium bicarbonate or sodium bicarbonate type water. An analysis of variance test indicates that water in units A and B is not significantly (0.05 probability level) different in composition. Mean concentrations of dissolved solids, sodium, and chloride are significantly different between units A and C, whereas mean concentrations of dissolved solids, calcium, sodium, and chloride are significantly different between units A and D. Groundwater in the surficial aquifer system in northwestern Dade County is more mineralized than groundwater that occurs elsewhere in the county (except in coastal areas affected by saltwater). An analysis of variance test indicates that mean concentrations of dissolved solids, sodium, and chloride at sites in northwestern Dade County, developed in units A and C, are significantly different from the mean concentrations of these constituents at sites throughout the rest of the county. Water in this part of the surficial aquifer system is similar to highly mineralized water found in western Broward County. However, in northwestern Dade County, the groundwater has been diluted to a greater extent by less mineralized recharge water. Overall, groundwater in the four major geohydrologic units of the surficial aquifer system of Dade County is suitable for most uses. Maximum concentrations of sodium, chloride, color, fluoride, iron, and dissolved solids in some areas of Dade County, however, exceed maximum contaminant levels established by the Florida Primary and Secondary Drinking Water Regulations standards. (Author 's abstract)

Florida↗

Sustainable aquifer management for food security

In aquifer-dependent regions, balancing aquifer protection, desalination, economic development, agricultural irrigation, and food security can be better managed through discovery and development of sources of sustainable groundwater pumping. Aquifer desalination for irrigation to protect food security can mitigate pressure on local freshwater aquifers. Despite its importance, little peer reviewed work to date has identified the economic capacity to pay for aquifer desalination for irrigation to mitigate freshwater aquifer drawdown. The novel contribution of this work is the development and application of an innovative method to assess the economic capacity to pay for aquifer desalination for irrigation for a recently discovered large saline aquifer. It develops an original framework to assess the capacity to pay for aquifer desalination, the results of which can help guide policymakers on efficient and sustainable pumping approaches across users, aquifers, and time periods. A mathematical programming model is developed to economically analyze the 200 billion cubic meter Lotikipi Aquifer, discovered in 2013 in northern Kenya using modern remote sensing methods. While initial pumping of the Lotikipi Aquifer was halted due to high groundwater salinity levels, interest remains strong in assessing the economic capacity to pay for groundwater desalination because of its potential role in protecting regional food security generated by aquifer pumping for irrigation. The model is formulated by calibrating optimized pumping patterns in two existing freshwater aquifers to replicate observed historical pumping levels. Based on that exercise, a second model is developed to identify a least cost set of pumping restrictions that return each of three regional aquifers to starting conditions over a seven-year time period. A third model extends the second by adding a constraint of a minimum required level of food grain security supported by irrigation pumping from the aquifer system. Results show that the economic capacity to pay for aquifer desalination for irrigated agriculture lies in the range of $0.08 - $0.18 USD per cubic meter under current economic conditions and desalination technologies available. While this economic capacity to pay is lower than its current cost in most places, the future could be more optimistic. Advances in desalination technology, higher crop prices, technical advance in agriculture, and development of drought-resistant crops can all contribute to a future capacity to economically justify the expense of desalination.

Agricultural Systems Journal↗

Relation of hydrologic processes to groundwater and surface-water levels and flow directions in a dune-beach complex at Indiana Dunes National Lakeshore and Beverly Shores, Indiana

The potential for high groundwater levels to cause wet basements (groundwater flooding) is of concern to residents of communities in northwestern Indiana. Changes in recharge from precipitation increases during 2006-9, water-level changes from restoration of nearby wetlands in the Great Marsh in 1998-2002, and changes in recharge due to the end of groundwater withdrawals for water supply since 2005 in a community at Beverly Shores, Ind., were suspected as factors in increased groundwater levels in an unconfined surficial aquifer beneath nearby parts of a dune-beach complex. Results of this study indicate that increased recharge from precipitation and snowmelt was the principal cause of raised water levels in the dune-beach complex from 2006 to 2009. Annual precipitation totals in 2006-9 ranged from 43.88 to 55.75 inches per year (in/yr) and were substantially greater than the median 1952-2009 precipitation of 36.35 in/yr. Recharge to groundwater from precipitation in 2006-9 ranged from 13.5 to 22 in/yr; it was higher than the typical 11 in/yr because of large precipitation events and precipitation amounts received during non-growing-season months. An estimated increase in net recharge from reduced groundwater use in Beverly Shores since 2005 ranged from 1.6 in/yr in 2006 to 1.9 in/yr in 2009. Surface-water levels in the wetland were as much as about 1.1 feet higher in 2007-9 (after the 1998-2002 wetland restoration) than during seasonally wet periods in 1979-89. Similar surface-water levels and ponded water were likely during winter and spring wet periods before and after wetland restoration. High water levels similar to those in 2009 were measured elsewhere in the dune-beach complex near a natural wetland during the spring months in 1991 and 1993 after receipt of near record precipitation. Recharge from similarly high precipitation amounts in 2008-9 was also a likely cause of high groundwater levels in other parts of the dune-beach complex, such as at Beverly Shores. Perennial mounding of the water table in the surficial aquifer indicates that the recharge that created the water-table mound originates within the dune-beach complex and not through flow from the adjacent hydrologic boundaries: the restored wetland, Lake Michigan, and Derby Ditch. Infiltrating precipitation causes most seasonal and episodic rises in groundwater levels beneath the dune-beach complex. Groundwater-level fluctuations lasting days to weeks in the dune-beach complex in 2008-9 were superimposed on a seasonal high water-table altitude that began with the recharge from snowmelt and rain in February 2009 and maintained through July 2009. Increases in water-table-mound altitude under the dune-beach complex recurred in 2008-9 in response to the largest rain events of 1 inch or more and to snowmelt. Smaller, shorter-term rises in water level after individual rain events persisted over hours to less than 1 week. Groundwater-level fluctuations varied over a relatively narrow range of about 2 to 3 feet, with no net fluctuations greater than 4 feet. Groundwater levels in or near low parts of the dune-beach complex were frequently within 0 to 6 feet of the land surface and indicate the potential for groundwater flooding. Groundwater-level gradients from the water-table mound to wells next to surface-water discharges increase after rainfall and snowmelt events and recede slowly as groundwater discharges from the aquifer. Evapotranspiration is responsible for part of the general pattern of decreasing water-table altitudes observed from May to August 2009. Rapid water-level rises in the restored wetland after precipitation do not likely have an effect on groundwater flooding elsewhere in the dune-beach complex. Surface-water-level fluctuations during this study generally varied over a narrower range, approximately from 1 to 1.5 feet, as compared with groundwater fluctuations, except after a very large, 10.77-inch rainfall. Time-delayed and smaller groundwater-level

Indiana↗

Conceptual models of groundwater flow in the Grand Canyon region, Arizona

The conceptual models of groundwater flow outlined herein synthesize what is known and hypothesized about the groundwater-flow systems that discharge to the Grand Canyon of Arizona. These models interpret the hydrogeologic characteristics and hydrologic dynamics of the physical systems into a framework for understanding key aspects of the physical systems as they relate to groundwater flow and contaminant transport. This report describes five individual groundwater-flow systems draining to the Grand Canyon: Kaibab, Uinkaret-Kanab, Marble-Shinumo, Cataract, and Blue Spring. These systems are present in the saturated parts of the lower Paleozoic carbonate section exposed on the walls of the Grand Canyon; specifically, the Mississippian Redwall Limestone down through the Cambrian Muav Limestone of Tonto Group. Together, the systems described in this report compose the regional groundwater-flow system. Local to subregional flow systems in the sedimentary units of the overlying Permian section could provide transport pathways from the land surface to the regional flow system. Despite the potential importance of the local systems, the focus of this report is on the systems present in the lower Paleozoic section because all major springs in the Grand Canyon discharge from those units. The most important hydrogeologic characteristics include system boundaries imposed by major tectonic structures, and the degree to which karstification influences the magnitude and direction of flow in each system. Important hydrologic dynamics include locations and rates of potential groundwater recharge, vertical pathways to the regional aquifer, and the locations, magnitude, geochemical signature, and hydrostratigraphic setting of groundwater discharge from springs. Unknown properties or conditions that represent the greatest uncertainties in our current understanding of the regional groundwater-flow system are identified for additional consideration. Groundwater data are sparse owing to geographic remoteness and extreme depth to water throughout much of the study area. This paucity of information was diminished with the development of a structural contour map of the top and bottom surfaces of the regional aquifer, and a Soil-Water-Balance model that produces spatial distributions of rates of potential recharge. Investigation of the five groundwater-flow systems reveals important, though mostly qualitative, characteristics controlling the rates and directions of groundwater flow. Karstification has produced dissolution-enhanced conduit flow pathways to various degrees in each of the systems. Parts of each system exhibit relative structural uplift or downdropping of the hydrostratigraphic units of the regional aquifer, with some uplifted sections dipping inward toward the Grand Canyon and others dipping outward. The Kaibab groundwater system is archetypical of an uplifted, inward-dipping karst system, whereas the Blue Spring groundwater system and most of the Cataract groundwater system are representative instances of a downdropped or basin karst system. The Uinkaret-Kanab groundwater-flow system is structurally similar to the basin karst systems but karstification has not progressed to nearly the same degree. The Marble-Shinumo groundwater system does not fall cleanly into either category and its boundaries are the most uncertain of all the groundwater systems.

Arizona↗

Hydrogeology of sand and gravel aquifers in the Owasco Inlet watershed, Cayuga and Tompkins Counties, New York

This study is a continuation of a series of hydrogeologic appraisals that have been conducted since 1980, as part of a cooperative, long-term, detailed aquifer mapping program by the U.S. Geological Survey and the New York State Department of Environmental Conservation. These appraisals provide a foundation for wellhead protection programs, water-resource management and planning decisions, and groundwater remediation in upstate New York. The Owasco Inlet watershed drains north directly to Owasco Lake, one of New York’s Finger Lakes. The watershed is similar in form to other watersheds of the Finger Lakes region of New York State and is characterized by broad, smooth uplands punctuated by the 18-mile-long Owasco Inlet valley and the secondary valleys of Decker, Dresserville, and Hemlock Creeks. All these streams occupy “through valleys”—the valleys are continuous into adjacent watersheds such that watershed divides are high points within the valley bottoms, rather than in uplands. Most glacial deposits in the watershed were deposited during the Valley Heads Readvance and subsequent retreat (with at least one minor readvance). Estimates are that the Valley Heads Readvance likely peaked in the Cayuga basin about 17,000 (calendar) years ago based on dates from western New York and the Mohawk Valley in east-central New York. It was the last major ice advance of the Pleistocene Epoch in this region. This readvance covered the entire Owasco watershed. Valley Heads ice in the Owasco watershed was part of the Cayuga ice lobe, which flowed parallel to the Owasco Inlet valley as far south as Moravia, but which became progressively more eastward south of Locke, effectively raking across the valley from the west. The Owasco Inlet watershed drains north, like other Finger Lakes watersheds, and during deglaciation, northward meltwater drainage was largely blocked by the ice, which resulted in development of a series of regional proglacial lakes in which fine lacustrine sediments predominated. Lacustrine sediments constitute much of the Owasco Inlet valley fill, but it has been pointed out that there are no obvious outlet channels exiting the Owasco Inlet watershed, and it has been proposed that there was northward drainage of meltwater into the ice within the larger, neighboring Cayuga Trough. The lacustrine deposits in the Owasco Inlet valley are commonly underlain by thinner coarse-grained stratified deposits (subaqueous fans and eskers) of variable sorting and permeability. An exception occurs at Groton, N.Y., where retreating ice paused long enough for coarse-grained sediments to fill much of the valley. Smaller, higher elevation valleys have a wider variety of glacial valley-fill deposits, ranging from fine lacustrine sediments to sand and gravel to till. Groundwater is the sole source of water supply in the area; glacial sand and gravel aquifers are the primary water source in the Owasco Inlet and Decker Creek valleys, and fractured bedrock aquifers typically supply domestic wells in the remaining valleys and upland areas. Municipal supplies tap a variety of aquifer types in the Owasco Inlet valley and nearby uplands. The hamlet of Locke taps the extensive confined aquifer beneath fine-grained lacustrine deposits in the valley. The Village of Moravia taps a semiconfined aquifer that overlies proglacial lake deposits and is partly confined by overlying recent lake deposits and alluvium. Withdrawals from this aquifer may also induce water from the Owasco Inlet into the aquifer. The Village of Groton draws from two aquifers: (1) a thin upland unconfined sand and gravel aquifer tapped by an infiltration gallery and (2) the local unconfined or semiconfined aquifer in the Owasco Inlet valley. Another aquifer with potential for municipal supply is the unconfined aquifer in the Decker Creek valley near Wilson Corners. The confined aquifer in the lower Dresserville Creek valley may have some water-resource potential, but it is largely untested. Unconfined aquifers are the most susceptible to contamination from activities at land surface directly above the aquifer because precipitation and subsequent recharge can transport contaminants directly to the water table. Adjacent upland areas can also contribute contaminants. Confined aquifers are less susceptible to contamination from overlying land surface areas because confining units largely prevent downward movement of water. Recharge occurs elsewhere at unconfined upvalley locations and along valley walls where alluvial fans, ice-contact deposits, or stream incision into the valley wall may provide pathways for downward movement of groundwater.

New York↗

Spatial variability of the response to climate change in regional groundwater systems -- examples from simulations in the Deschutes Basin, Oregon

We examine the spatial variability of the response of aquifer systems to climate change in and adjacent to the Cascade Range volcanic arc in the Deschutes Basin, Oregon using downscaled global climate model projections to drive surface hydrologic process and groundwater flow models. Projected warming over the 21st century is anticipated to shift the phase of precipitation toward more rain and less snow in mountainous areas in the Pacific Northwest, resulting in smaller winter snowpack and in a shift in the timing of runoff to earlier in the year. This will be accompanied by spatially variable changes in the timing of groundwater recharge. Analysis of historic climate and hydrologic data and modeling studies show that groundwater plays a key role in determining the response of stream systems to climate change. The spatial variability in the response of groundwater systems to climate change, particularly with regard to flow-system scale, however, has generally not been addressed in the literature. Here we simulate the hydrologic response to projected future climate to show that the response of groundwater systems can vary depending on the location and spatial scale of the flow systems and their aquifer characteristics. Mean annual recharge averaged over the basin does not change significantly between the 1980s and 2080s climate periods given the ensemble of global climate models and emission scenarios evaluated. There are, however, changes in the seasonality of groundwater recharge within the basin. Simulation results show that short-flow-path groundwater systems, such as those providing baseflow to many headwater streams, will likely have substantial changes in the timing of discharge in response changes in seasonality of recharge. Regional-scale aquifer systems with flow paths on the order of many tens of kilometers, in contrast, are much less affected by changes in seasonality of recharge. Flow systems at all spatial scales, however, are likely to reflect interannual changes in total recharge. These results provide insights into the possible impacts of climate change to other regional aquifer systems, and the streams they support, where discharge points represent a range of flow system scales.

Oregon↗