Title: Past aquifer responses to climate recorded by fossil groundwater
Billions of people rely upon groundwater for drinking water and agriculture, yet predicting how climate change may affect aquifer storage remains challenging. To gain insight beyond the short historical record, we reconstruct changes in groundwater levels in western North America during the last glacial termination (LGT, ~20 to 11 thousand years ago) using noble gas isotopes. Our reconstructions indicate remarkable stability of water table depth in a Pacific Northwest aquifer throughout the LGT despite increasing precipitation, closely matching independent Earth system model (ESM) simulations. In the American Southwest, ESM simulations and noble gas isotopes both suggest a pronounced LGT decline in water table depth in in response to decreasing precipitation, indicating distinct regional groundwater responses to climate. Despite the hydrologic simplicity of ESMs, their agreement with proxy reconstructions of past water table depth suggests that these models hold value in understanding groundwater dynamics and projecting large-scale aquifer responses to climate forcing.  more » « less
Award ID(s):
2102984 2238641
PAR ID:
10615897
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Science Advances
Date Published:
Journal Name:
Science Advances
Volume:
11
Issue:
24
ISSN:
2375-2548
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The abundance and isotopic composition of noble gases dissolved in water have many applications in the geosciences. In recent years, new analytical techniques have opened the door to the use of high-precision measurements of noble gas isotopes as tracers for groundwater hydrology, oceanography, mantle geochemistry, and paleoclimatology. These analytical advances have brought about new measurements of solubility equilibrium isotope effects (SEIEs) in water (i.e., the relative solubilities of noble gas isotopes) and their sensitivities to the temperature and salinity. Here, we carry out a suite of classical molecular dynamics (MD) simulations and employ the theoretical method of quantum correction to estimate SEIEs for comparison with experimental observations. We find that classical MD simulations can accurately predict SEIEs for the isotopes of Ar, Kr, and Xe to order 0.01‰, on the scale of analytical uncertainty. However, MD simulations consistently overpredict the SEIEs of Ne and He by up to 40% of observed values. We carry out sensitivity tests at different temperatures, salinities, and pressures and employ different sets of interatomic potential parameters and water models. For all noble gas isotopes, the TIP4P water model is found to reproduce observed SEIEs more accurately than the SPC/E and TIP4P/ice models. Classical MD simulations also accurately capture the sign and approximate magnitude of temperature and salinity sensitivities of SEIEs for heavy noble gases. We find that experimental and modeled SEIEs generally follow an inverse-square mass dependence, which implies that the mean-square force experienced by a noble gas atom within a solvation shell is similar for all noble gases. This inverse-square mass proportionality is nearly exact for Ar, Kr, and Xe isotopes, but He and Ne exhibit a slightly weaker mass dependence. We hypothesize that the apparent dichotomy between He–Ne and Ar–Kr–Xe SEIEs may result from atomic size differences, whereby the smaller noble gases are more likely to spontaneously fit within cavities of water without breaking water–water H-bonds, thereby experiencing softer collisions during translation within a solvation shell. We further speculate that the overprediction of simulated He and Ne SEIEs may result from the neglection of higher-order quantum corrections or the overly stiff representation of van der Waals repulsion by the widely used Lennard-Jones 6–12 potential model. We suggest that new measurements of SEIEs of heavy and light noble gases may represent a novel set of constraints with which to refine hydrophobic solvation theories and optimize the set of interatomic potential models used in MD simulations of water and noble gases. 
    more » « less
  2. null (Ed.)
    Abstract. Groundwater table dynamics extensively modify the volume of the hyporheic zoneand the rate of hyporheic exchange processes. Understanding the effects ofdaily groundwater table fluctuations on the tightly coupled flow and heattransport within hyporheic zones is crucial for water resourcesmanagement. With this aim in mind, a physically based model is used to explorehyporheic responses to varying groundwater table fluctuationscenarios. The effects of different timing and amplitude of groundwater tabledaily drawdowns under gaining and losing conditions are explored in hyporheiczones influenced by natural flood events and diel river temperaturefluctuations. We find that both diel river temperature fluctuations and dailygroundwater table drawdowns play important roles in determining thespatiotemporal variability of hyporheic exchange rates, temperature ofexfiltrating hyporheic fluxes, mean residence times, and hyporheicdenitrification potentials. Groundwater table dynamics present substantiallydistinct impacts on hyporheic exchange under gaining or losing conditions. Thetiming of groundwater table drawdown has a direct influence on hyporheicexchange rates and hyporheic buffering capacity on thermaldisturbances. Consequently, the selection of aquifer pumping regimes hassignificant impacts on the dispersal of pollutants in the aquifer and thermalheterogeneity in the sediment. 
    more » « less
  3. Coastal aquifers experience water table fluctuations that push and pull water and air through organic‐rich soils. This exchange affects the supply of oxygen, dissolved organic carbon (DOC), and nitrogen (N) to shallow aquifers and influences groundwater quality. To investigate the fate of N species, we used a meter‐long column containing a sequence of natural organic topsoil and aquifer sediments. A fluctuating head was imposed at the column bottom with local, nitrate‐rich groundwater (16.5 mg/L NO3‐N). We monitored in‐situ redox potential and collected pore water samples for analysis of inorganic N species and DOC over 16 days. Reactive processes were more complex than anticipated. The organic‐rich topsoil remained anaerobic, while mineral sediments beneath alternated between aerobic, when the water table dropped and sucked air across preferential flow paths, and anaerobic conditions, when the water table was high. A fluid flow and reactive transport model shows that when the water table rises into organic‐rich soils, it limits the flow of oxygen, while the soils release DOC, which stimulates the removal of nitrate from groundwater by denitrification. At the end of the experiment, we introduced seawater to the column to mimic a storm surge. Seawater mobilized N and DOC from shallow soil horizons, which could reach the aquifer if the surge is long enough. These processes are relevant for groundwater quality in developed coastal areas with anthropogenic N sources, as climate change and rising seas will drive changes in water table and flood dynamics. 
    more » « less
  4. Abstract Non-perennial rivers and streams make up over half the global river network and are becoming more widespread. Transitions from perennial to non-perennial flow are a threshold-type change that can lead to alternative stable states in aquatic ecosystems, but it is unknown whether streamflow itself is stable in either wet (flowing) or dry (no-flow) conditions. Here, we investigated drivers and feedbacks associated with regime shifts between wet and dry conditions in an intermittent reach of the Arkansas River (USA) over the past 23 years. Multiple lines of evidence suggested that these regimes represent alternative stable states, including (a) significant jumps in discharge time series that were not accompanied by jumps in flow drivers such as precipitation and groundwater pumping; (b) a multi-modal state distribution with 92% of months experiencing no-flow conditions for <10% or >90% of days, despite unimodal distributions of precipitation and pumping; and (c) a hysteretic relationship between climate and flow state. Groundwater levels appear to be the primary control over the hydrological regime, as groundwater levels in the alluvial aquifer were higher than the stream stage during wet regimes and lower than the streambed during dry regimes. Groundwater level variation, in turn, was driven by processes occurring at both the regional scale (surface water inflows from upstream, groundwater pumping) and the reach scale (stream–aquifer exchange, diffuse recharge through the soil column). Historical regime shifts were associated with diverse pressures including network disconnection caused by upstream water use, increased flow stability potentially associated with reservoir operations, and anomalous wet and dry climate conditions. In sum, stabilizing feedbacks among upstream inflows, stream–aquifer interactions, climate, vegetation, and pumping appear to create alternative wet and dry stable states at this site. These stabilizing feedbacks suggest that widespread observed shifts from perennial to non-perennial flow will be difficult to reverse. 
    more » « less
  5. Abstract Accurate groundwater representation in land surface models (LSMs) is vital for water and energy cycle studies, water resource assessments, and climate projections. Yet, many LSMs do not consider key processes including lateral groundwater flow and aquifer pumping, especially at the global scale. This study simulates these processes using an enhanced version of the Community Land Model (CLM5) and evaluates their roles at three spatial resolutions (0.5°, 0.25°, 0.1°). Results show that lateral flow strongly modulates water table depth and capillary rise at all resolutions. The magnitude of mean lateral flow increases from 25 mm/year at 0.5° to 36 mm/year at 0.25°, and 52 mm/year at 0.1° resolution, with pumping inducing lateral flow even at 0.5° (∼50 km), a typical grid size in global LSMs. Further, lateral flow alters runoff in regions with high recharge and shallow water table (e.g., eastern North America and Amazon basin), and soil moisture and ET in regions with comparatively low recharge and deeper water table (e.g., western North America, central Asia, and Australia) through enhanced capillary rise. Runoff alteration by lateral flow increases substantially with resolution, from a maximum of 15 mm/month at 0.5° to 20 mm/month and 25 mm/month at 0.25° and 0.1°, respectively; the impact of resolution on soil moisture and ET is less pronounced. While the model does not fully capture deeper water tables—warranting further enhancements—it provides valuable insights on how lateral groundwater flow impacts land surface processes, highlighting the importance of lateral groundwater flow and pumping in global LSMs. 
    more » « less