Search NASA⌕ Search

SEARCH · Search NASA

Results for “Groundwater storage”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records

Declining groundwater storage expected to amplify mountain streamflow reductions in a warmer world

Abstract Groundwater interactions with mountain streams are often simplified in model projections, potentially leading to inaccurate estimates of streamflow response to climate change. Here, using a high-resolution, integrated hydrological model extending 400 m into the subsurface, we find groundwater an important and stable source of historical streamflow in a mountainous watershed of the Colorado River. In a warmer climate, increased forest water use is predicted to reduce groundwater recharge resulting in groundwater storage loss. Losses are expected to be most severe during dry years and cannot recover to historical levels even during simulated wet periods. Groundwater depletion substantially reduces annual streamflow with intermittent conditions predicted when precipitation is low. Expanding results across the region suggests groundwater declines will be highest in the Colorado Headwater and Gunnison basins. Our research highlights the tight coupling of vegetation and groundwater dynamics and that excluding explicit groundwater response to warming may underestimate future reductions in mountain streamflow.

Carroll, Rosemary W. H. (ORCID:0000000293028074)↗

Groundwater and Surface Water Flow (GSFLOW) model files to explore bedrock circulation depth and porosity in Copper Creek, Colorado

This data package contains integrated hydrological model input and output files for Copper Creek, Colorado (24 km2), a tributary of the East River located in the headwaters of the Upper Colorado River Basin. The model code is the U.S. Geological Survey (USGS) Groundwater and Surface Water Flow (GSFLOW) model. The model contains a 100-m grid resolution and a daily timestep. The land surface model is dynamically linked to a three-dimensional groundwater flow model that allows for streamflow gaining and losing conditions. The groundwater model contains 12 model layers and extends 400 m below land surface. The original Copper Creek model was modified to contain geologic layers representing saprolite, shallow bedrock, and deep bedrock. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. For the shallow case, median flow depths occur in the shallow saprolite at depths <8 m, while the deep case promotes a median groundwater flow depth of 100 m. With this modeling framework we compare streamflow response to a plausible worst-case drought lasting up to five years. Streamflow metrics of analysis include average streamflow, fraction of stream network that is dry, no-flow duration, average groundwater flow to streams and time to recovery following the drought. Results and implications are presented in a paper submitted to Geophysical Research Letters titled, "The role of bedrock circulation depth and porosity in mountain streamflow response to prolonged drought" by Rosemary WH. Carroll, Andrew H. Manning and Kenneth H Williams. A Readme.txt file provides instructions on how to download all model files and execute each model scenario. In addition to the GSFLOW output/prms/copper_drought.csv file containing daily basin water stores and fluxes (refer to GSFLOW manual) and the output/prms/copper_drought_statvar.dat file with output defined in the gsflow3.control file (refer to GSFLOW Manual), output files also include spatially distributed daily values of total evapotranspiration, canopy evaporation, precipitation, snowfall, infiltration, snow water equivalent, potential evapotranspiration, recharge, sublimation, soil moisture, contributing interflow, water table elevations, changes in groundwater storage, groundwater evapotranspiration, interbasin groundwater flow (limited to the alluvium below the stream outlet), and surface-groundwater exchanges within the river system.

54 ENVIRONMENTAL SCIENCES↗

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Matrix Diffusion Controls Mountain Hillslope Groundwater Ages and Inferred Storage Dynamics

Groundwater age distributions provide fundamental insights on coupled water and biogeochemical processes in mountain watersheds. Field-based studies have found mixtures of young and old-aged groundwater in mountain catchments underlain by bedrock; yet, the processes that dictate these groundwater age distributions are poorly understood. In this work, we use the coupled ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models to simulate groundwater age distributions on a lower montane hillslope in the East River Watershed, Colorado (USA). We develop a convolution-based approach to propagate fracture-matrix diffusion processes to the EcoSLIM advection-dominated age distributions. We compare observed 3 H and 4 He concentrations from two groundwater wells against model predictions that have varying advective transport times and matrix diffusion magnitudes. Based on a Monte Carlo analysis that considers uncertain matrix and fracture parameters, we find that matrix diffusion is needed to jointly predict 3 H and 4 He observations at both wells. The advection-dominated age distributions lack adequate mixing of young and old-aged water to capture the observed co-occurrence of 3 H and 4 He. The model scenario that best matches the 3 H, 4 He, and water level observations when considering both advective flowpath and matrix diffusion mixing processes has a dynamic bedrock groundwater reservoir that is susceptible to considerable storage losses during low-snow periods. This dynamic groundwater system amplifies the need to assimilate deeper bedrock groundwater into watershed hydro-biogeochemical predictions. This work further highlights the importance of considering matrix diffusion when interpreting environmental tracers in bedrock groundwater systems.

54 ENVIRONMENTAL SCIENCES↗

Evapotranspiration depletes groundwater under warming over the contiguous United States

A warmer climate increases evaporative demand. However, response to warming depends on water availability. Existing earth system models represent soil moisture but simplify groundwater connections, a primary control on soil moisture. Here we apply an integrated surface-groundwater hydrologic model to evaluate the sensitivity of shallow groundwater to warming across the majority of the US. We show that as warming shifts the balance between water supply and demand, shallow groundwater storage can buffer plant water stress; but only where shallow groundwater connections are present, and not indefinitely. As warming persists, storage can be depleted and connections lost. Similarly, in the arid western US warming does not result in significant groundwater changes because this area is already largely water limited. The direct response of shallow groundwater storage to warming demonstrates the strong and early effect that low to moderate warming may have on groundwater storage and evapotranspiration.

54 ENVIRONMENTAL SCIENCES↗

Projecting end-of-century climate extremes and their impacts on the hydrology of a representative California watershed

Abstract. In California, it is essential to understand the evolution of water resources in response to a changing climate to sustain its economy and agriculture and to build resilient communities. Although extreme conditions have characterized the historical hydroclimate of California, climate change will likely intensify hydroclimatic extremes by the end of the century (EoC). However, few studies have investigated the impacts of EoC extremes on watershed hydrology. We use cutting-edge global climate and integrated hydrologic models to simulate EoC extremes and their effects on the water-energy balance. We assess the impacts of projected driest, median, and wettest water years under Representative Concentration Pathway (RCP) 8.5 on the hydrodynamics of the Cosumnes River basin. Substantial changes to annual average temperature (>+2.5 ∘C) and precipitation (>+38 %) will characterize the EoC extreme water years compared to their historical counterparts. A shift in the dominant form of precipitation, mostly in the form of rain, is projected to fall earlier. These changes reduce snowpack by more than 90 %, increase peak surface water and groundwater storages up to 75 % and 23 %, respectively, and drive the timing of peak storage to occur earlier in the year. Because EoC temperatures and soil moisture are high, both potential and actual evapotranspiration (ET) increase. The latter, along with the lack of snowmelt in the warm EoC, causes surface water and groundwater storages to significantly decrease in summer, with groundwater showing the highest rates of decrease. These changes result in more ephemeral EoC streams with more focused flow and increased storage in the mainstem of the river network during the summer.

54 ENVIRONMENTAL SCIENCES↗

An ensemble of 48 physically perturbed model estimates of the 1/8° terrestrial water budget over the conterminous United States, 1980–2015

Terrestrial water budget (TWB) data over large domains are of high interest for various hydrological applications. Spatiotemporally continuous and physically consistent estimations of TWB rely on land surface models (LSMs). As an augmentation of the operational North American Land Data Assimilation System Phase 2 (NLDAS-2) four-LSM ensemble, this paper describes a dataset simulated from an ensemble of 48 physics configurations of the Noah LSM with multi-physics options (Noah-MP). The 48 Noah-MP physics configurations are selected to give a representative cross-section of commonly used LSMs for parameterizing runoff, atmospheric surface layer turbulence, soil moisture limitation on photosynthesis, and stomatal conductance. The dataset spans from 1980 to 2015 over the conterminous United States (CONUS) at a monthly temporal resolution and a 1/8° spatial resolution. The dataset variables include total evapotranspiration and its constituents (canopy evaporation, soil evaporation, and transpiration), runoff (the surface and subsurface components), as well as terrestrial water storage (snow water equivalent, four-layer soil water content from the surface down to 2 m, and the groundwater storage anomaly). The dataset is available at https://doi.org/10.5281/zenodo.7109816. Evaluations carried out in this study and previous investigations show that the ensemble per forms well in reproducing the observed terrestrial water storage, snow water equivalent, soil moisture, and runoff. Noah-MP complements the NLDAS models well, and adding Noah-MP consistently improves the NLDAS es timations of the above variables in most areas of CONUS. Besides, the perturbed-physics ensemble facilitates the identification of model deficiencies. The parameterizations of shallow snow, spatially varying groundwater dynamics, and near-surface atmospheric turbulence should be improved in future model versions.

54 ENVIRONMENTAL SCIENCES↗

The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach

Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Wang, Lijing [University of Connecticut, Storrs, C↗

Integrated GW Farm ABM

This Data Repository includes data used for the integrated groundwater- farm ABM model, raw model output from scenario ensemble, and processed outputs that isolate the groundwater storage depletion outcomes for the 35,000 farm cells. Model Inputs: Farm ABM Inputs: This folder contains the input data used by the integrated groundwater - farm ABM modelling script (Python file) used for the high performance computing (HPC) experiments. The sub-folder "data inputs" contains all of the farm attribute data, while the three files in the folder have the hydrogeological data lookup table (NLDAS Cost Curve Attributes.csv), a lookup table (Theis well function table.csv) for the groundwater cost curve function, and the farm indexes and corresponding NLDAS ids for all of the cells run in this experiment (nldas farms subset final.csv). NLDAS Cost curve hydrogeological data: Hydrogeological data aggregated to 1/8 degree resolution and aligned with the NLDAS grid. Parameters include: water depth below ground surface [meters], subsurface porosity [unitless], aquifer depth from ground surface to aquifer bottom [meters], annual average recharge (USGS: mm, Doll: meters), and three different hydraulic conductivity (K) values (meters/day). The three K values represent the mean value from Gleeson et al. (2018), one standard deviation above the mean from Gleeson et al. (2018), and the de Graaf et al. 2020 modifications to certain lithologies. Additional information about these datasets and their processing are documented in the supplement to Yoon et al. 2025 (in review). Output: Raw outputs: This folder contains a .zip file that has model outputs for the entire scenario ensemble. There is one csv for each farm id, using the format "farm farmid cases.csv". The relationship between the farm id and NLDAS id is defined by the "nldas farms subset final.csv" located in the Farm ABM Inputs folder. Each csv has 625 rows, corresponding to 625 combinations of different scenario parameter values. Each row (scenario) represents the outcome of a 100 year simulation. Columns define scenario settings and summary statistics for each scenario. The first four columns define the scenario settings: "hydro ratio," "econ ratio," "K scenario," and "gamma scenario." The hydro and econ ratios are values passed to the modeling script that influence multipliers for other model parameters, as documented in the supplement to Yoon et al. 2025 (in review). The gamma multiplier is a coefficient multiplier applied to the baseline gamma values (values below 1 represent lower unobserved costs compared to baseline, values above 1 represent higher costs). The K scenario names represent K values of: "low": 0.5 m/d, "int 1": 2.5 m/d, "int 2": 10 m/d, "high": 50 m/d, and "gleeson": mean Gleeson K value. "Perc vol depleted" is the fraction of groundwater depleted at the end of the 100 simulation. Processed Output: Derived depletion outcomes from raw outputs: All of the individual csv files from the Raw outputs were aggregated into a single file that has the scenario settings and fraction depletion "Perc vol depleted" for every farm cell, for every scenario. The other two files define relationships between the farm id, NLDAS id, and local and major aquifer units, used for aquifer-level depletion analysis.

Agent based modeling↗

Using Sentinel-1 and GRACE satellite data to monitor the hydrological variations within the Tulare Basin, California

Abstract Subsidence induced by groundwater depletion is a grave problem in many regions around the world, leading to a permanent loss of groundwater storage within an aquifer and even producing structural damage at the Earth’s surface. California’s Tulare Basin is no exception, experiencing about a meter of subsidence between 2015 and 2020. However, understanding the relationship between changes in groundwater volumes and ground deformation has proven difficult. We employ surface displacement measurements from Interferometric Synthetic Aperture Radar (InSAR) and gravimetric estimates of terrestrial water storage from the Gravity Recovery and Climate Experiment (GRACE) satellite pair to characterize the hydrological dynamics within the Tulare basin. The removal of the long-term aquifer compaction from the InSAR time series reveals coherent short-term variations that correlate with hydrological features. For example, in the winter of 2018–2019 uplift is observed at the confluence of several rivers and streams that drain into the southeastern edge of the basin. These observations, combined with estimates of mass changes obtained from the orbiting GRACE satellites, form the basis for imaging the monthly spatial variations in water volumes. This approach facilitates the quick and effective synthesis of InSAR and gravimetric datasets and will aid efforts to improve our understanding and management of groundwater resources around the world.

54 ENVIRONMENTAL SCIENCES↗

Warming and snow loss increase reliance on old groundwater in a Colorado River headwater

Atmospheric warming is reducing snowpack, with uncertain effects on mountainous streamflow, a crucial water resource. Despite limited historical observations of groundwater–streamflow interactions above 2,500 m, new measurements in the Upper Colorado River headwaters indicate declining groundwater storage that is dated decades to millennia old. Here we use integrated hydrologic modelling spanning water years 2015–2021 to determine whether the loss of old-age groundwater buffers streamflow during low-snow years and whether that loss is exacerbated with warming. Results show that old-groundwater contributions to streams remain relatively steady through time, unlike the more variable contributions from young groundwater. Numerical experiments of increased surface air temperatures (+2.5 °C and +4 °C) increase rain–snow fractions and evapotranspiration and decrease runoff ratio by 2–3% per degree Celsius increase. As streamflow declines with warming, the age of groundwater supporting it gets older, in part owing to intermediate-aged (1–3 year) groundwater declining twice as fast. Simulations show that water table depths at higher elevations (>3,700 m) decline disproportionately and fail to recover even during wet years. These findings suggest altered groundwater–streamflow interactions with warming and snow loss, with implications for water resources.

Siirila-Woodburn, Erica R. [Lawrence Berkeley Nati↗

Monitoring Terrestrial Water Storage, Drought and Seasonal Changes in Central Oklahoma With Ambient Seismic Noise

Abstract Significant imbalances in terrestrial water storage (TWS) and severe drought have been observed around the world as a consequence of climate changes. Improving our ability to monitor TWS and drought is critical for water‐resource management and water‐deficit estimation. We use continuous seismic ambient noise to monitor temporal evolution of near‐surface seismic velocity, dv / v , in central Oklahoma from 2013 to 2022. The derived dv / v is found to be negatively correlated with gravitational measurements and groundwater depths, showing the impact of groundwater storage on seismic velocities. The hydrological effects involving droughts and recharge of groundwater occur on a multi‐year time scale and dominate the overall derived velocity changes. The thermoelastic response to atmospheric temperature variations occurs primarily on a yearly timescale and dominates the superposed seasonal velocity changes in this study. The occurrences of droughts appear simultaneously with local peaks of dv / v , demonstrating the sensitivity of near‐surface seismic velocities to droughts.

58 GEOSCIENCES↗

The Role of Atmospheric Rivers on Groundwater: Lessons Learned From an Extreme Wet Year

Abstract In the coastal regions of the western United States, atmospheric rivers (ARs) are associated with the largest precipitation generating storms and contribute up to half of annual precipitation, but the impact of ARs on the integrated hydrologic cycle, specifically on groundwater storage and hydrodynamics, is largely unknown. To better explore the hydrologic behavior of AR versus non‐AR event precipitation, we present a novel combination of two water tracking methods (one in the atmosphere and one in the subsurface) to explicitly track the full lifecycle of water parcels generated by ARs. Simulations of northern California's Cosumnes River watershed during the record wet 2017 water year are performed via the coupling of a high‐resolution regional climate model and a land surface‐groundwater model accounting for lateral groundwater flow. Despite ARs contributing more precipitation than non‐AR storms, we find less AR water is preferentially stored in aquifers by year end. Fractionally, ARs result in 300% less snow derived groundwater‐recharged compared to non‐AR precipitation. Rain‐on‐snow (RoS) plays an important role in AR‐driven discharge, where over 50% of total discharge from ARs snow is from RoS events. Finally, despite record‐breaking annual precipitation, simulated groundwater depletion occurs by year end due to estimates of groundwater pumping activities. The results from these simulations serve as a partial analogue of future hydrologic conditions where ARs are expected to intensify and provide a greater fraction of annual precipitation due to climate change.

54 ENVIRONMENTAL SCIENCES↗

Evaluating How Climate Adaptation Measures Affect the Interconnected Water‐Energy Resource Systems of the Western United States

Abstract The Western US faces increasing water stress from the impacts of climate change, making it difficult to meet water demands for the region's cities, agriculture, and hydropower generators. Existing literature suggests that climate adaptation measures such as water conservation, cropland retirement, wastewater recycling, and managed aquifer recharge can alleviate some of these challenges. Few analyses, however, compare the relative efficacy and system‐wide effects of these adaptations under different climate projections across the entire Western United States. Here we use a Western US‐wide water systems model to evaluate, by sector and sub‐region, how the widespread implementation of these adaptive measures impacts water demands, water deliveries, and electricity use related to the water system for three different climate projections. We find that wastewater recycling has greater potential to lower unmet indoor water demands than urban indoor water conservation measures. However, when implemented at scale, indoor water conservation reduces electricity use by an average of 683 Terawatt hours while wastewater recycling increases energy use by an average of 721 Terawatt hours, cumulatively from 2020 to 2070. Cropland retirement and aquifer recharge adaptations increase the ability to meet agricultural water demand, increase groundwater storage, and reduce summertime electricity use. While most of these findings are consistent across different climate projections, the benefits of aquifer recharge are sensitive to spatial variation of precipitation. Given the limitations and tradeoffs of each individually, the results suggest that a portfolio of adaptation measures will be needed for a climate‐resilient water and energy future in the Western US. Plain Language Summary The Western US faces increased water stress from the impacts of climate change, making it difficult to meet demands for cities, agriculture, and hydropower facilities. Adaptation measures like water conservation, retiring agricultural lands, recycling wastewater, and storing water underground can address these challenges. However, there is little modeling to understand the impact of these adaptations if they were implemented across the entire Western US while also considering climate change. We use a Western US‐wide water systems model to evaluate how implementing adaptation measures impacts groundwater levels, the ability to meet water demands, and electricity use related to water in three different possible climate futures. We find that recycling wastewater to drinking water standards does a better job of meeting urban water demands than water conservation measures. However, indoor water conservation reduces energy use for water while recycling wastewater increases energy use for water. Retiring agricultural lands and storing water underground both increase the ability to meet agricultural water demands, increase underground water storage, and reduce summertime energy use. No adaptation measures provides benefits across every metric we track. Therefore, multiple adaptation measures will likely be needed to achieve a climate‐resilient future for energy and water in the Western US. Key Points Water‐related climate change adaptations have different benefits and tradeoffs for interconnected water‐energy systems in the Western U.S. Wastewater recycling most increases urban water coverage, but has a significant tradeoff in the form of greater electricity use Cropland retirement and aquifer recharge benefit water and electricity systems most in the summer when those systems are most stressed

Singhal, A↗

Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics Within Water‐Tagging Enabled Hydrologic Models

Determining the age distribution of water exiting a catchment is important for understanding groundwater storage and mixing. New water-tagging capabilities within models track precipitation events as they move through simulated storages, yet forward modeling of individual events may not systematically capture the full transit time distribution (TTD). Here, we present a “sequential precipitation input tagging” (SPIT) framework to tag all input precipitation at regular intervals during extended model simulations. Monthly tags over 7 years were applied at six National Ecological Observatory Network sites to calculate TTDs and derive mean virtual tracer age, $\overline{T_{V}}$, fractions of young water, F yw , and hydrologic tracer concentrations (water isotopes δ 18 O and δ 2 H) within a tagging enabled version of the Weather Research and Forecast hydrologic model (WRF-Hydro). Throughout seven simulation years, the fraction of simulated discharge derived from tagged events, F tag , increased each year, with the final year's F tag ranging from 66% to 100% and highlights the need to apply SPIT over many years to understand TTDs. When the F tag was >75%, simulated $\overline{T_{V}}$ ranged 179–923 days and F yw 0.6%–23.9%, with daily values exhibiting a power-law relationship with precipitation, discharge, and groundwater. Through implementation of SPIT, we find this hydrologic model configuration performs poorly in estimation of $\overline{T_{V}}$ and F yw (root mean squared error of 469 days and 14.4% respectively), suggesting it misrepresents subsurface mixing. Thus, the SPIT framework provides a reproducible approach to calculate watershed transit times within tagging enabled models and thereby assess and improve representation of hydrologic processes.

fraction of young water↗

Evaluating the effects of heatwave events on hydrological processes in the contiguous United States (2003–2022)

Extreme heat and drought conditions are affecting water availability in many regions worldwide, leading to negative impacts on human societies, agriculture, and ecosystems. However, current research lacks comprehensive spatiotemporal analysis examining the interplay between multiple hydrological factors and heatwave events, especially in the context of climate change. This research broadly pertains to understanding the dynamics of hydrological factors and their potential responses to heatwave during warm seasons across the contiguous United States for the period from 2003 to 2022. Utilizing data from the Global Land Data Assimilation System (GLDAS), we analyzed surface runoff, evapotranspiration (ET), precipitation, Groundwater Storage (GWS), Root Zone Soil Moisture (RZSM), and Total Water Storage (TWS) to discern annual patterns and the impacts of heatwave. Further, the spatial patterns of heatwave highlighted a higher occurrence in the western, central, and northeastern U.S., with longer average durations in the western and south-central regions. These events are predominantly dry, characterized by low Relative Humidity (RH), except in the southeastern U.S., where heatwave coincide with high RH levels. Post-heatwave analysis indicated a reduction in GWS, TWS, RZSM, and ET, alongside an increase in surface runoff, RH, and precipitation. An in-depth examination of rainfall and temperature dynamics during heatwave revealed weak correlations between rainfall and temperature, as well as between rainfall and heatwave duration, highlighting the complex nature of these interactions. The study also found an enhanced probability of rainfall following heatwave, particularly in the eastern regions, drawing attention to the potential for increased flood risks post-heatwave. Our findings contribute to the growing body of knowledge on the impacts of heatwave on hydrological factors, providing valuable insights for climate change adaptation and water resource management strategies.

54 ENVIRONMENTAL SCIENCES↗