Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Groundwater Flux”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Data and scripts associated with a manuscript on residence time distribution simulation in two 10-kilometer long river sections

This data package is associated with the publication “On the Transferability of Residence Time Distributions in Two 10-km Long River Sections with Similar Hydromorphic Units” submitted to the Journal of Hydrology (Bao et al. 2024).Quantifying hydrologic exchange fluxes (HEFs) at the stream-groundwater interface, along with their residence time distributions (RTDs) in the subsurface, is crucial for managing water quality and ecosystem health in dynamic river corridors. However, directly simulating high-spatial resolution HEFs and RTDs can be a time-consuming process, particularly for watershed-scale modeling. Efficient surrogate models that link RTDs to hydromorphic units (HUs) may serve as alternatives for simulating RTDs in large-scale models. One common concern with these surrogate models, however, is the transferability of the relationship between the RTDs and HUs from one river corridor to another. To address this, we evaluated the HEFs and the resulting RTD-HU relationships for two 10-kilometer-long river corridors along the Columbia River, using a one-way coupled three-dimensional transient surface-subsurface water transport modeling framework that we previously developed. Applying this framework to the two river corridors with similar HUs allows for quantitative comparisons of HEFs and RTDs using both statistical tests and machine learning classification models. This data package includes the model inputs files and the simulation results data. This data package contains 10 folders. The modeling simulation results data are in the folders 100H_pt_data and 300area_pt_data, for the study domain Hanford 100H and 300 area respectively. The remaining eight folders contain the scripts and data to generate the manuscript figures. The file-level metadata file (Bao_2024_Residence_Time_Distribution _flmd.csv) includes a list of all files contained in this data package and descriptions for each. The data dictionary file (Bao_2024_Residence_Time_Distribution _dd.csv) includes column header definitions and units of all tabular files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Development of a vadose zone advanced monitoring system: Tools to assess groundwater vulnerability

Performing repeat pore-fluid sampling over long time-scales can provide valuable information on unsaturated zone contaminants and their potential flux to ground water. This information can be used to manage groundwater remedies and identify contaminants that need to be sequestered in the vadose zone to minimize flux to ground water. Pore-water samples are commonly used to obtain contaminant concentrations within the vadose zone, but existing methods are limited as they only provide a single sample at one location and time. The vadose zone advanced monitoring system (VZAMS) has been designed to integrate multiple technologies into a single down-borehole system that allows for sampling of pore fluids (liquid and gas) to provide information about contamination and hydraulic conditions at multiple depths (~0.3-m intervals) within a cased borehole. Testing has been completed at the laboratory scale to verify the sampling elements of VZAMS, including geochemical testing for representative contaminants known to exist at the Hanford Site, located in southeastern Washington State. Physical tests focused on the ability of the sampler to draw fluid under unsaturated conditions. Initial geochemical testing showed that the stainless steel material used with the porous cuff may affect the sampled concentrations of redox-sensitive contaminants under very dry conditions. Additional laboratory testing demonstrated that the VZAMS components are able to collect representative samples for substances of interest under expected field conditions. In this paper, the design and functionality of a novel instrument are demonstrated in support of subsequent testing in the field.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Structure and composition of natural ferrihydrite nano-colloids in anoxic groundwater

Fe-rich mobile colloids play vital yet poorly understood roles in the biogeochemical cycling of Fe in groundwater by influencing organic matter (OM) preservation and fluxes of Fe, OM, and other essential (micro-)nutrients. Yet, few studies have provided molecular detail on the structures and compositions of Fe-rich mobile colloids and factors controlling their persistence in natural groundwater. Here, we provide comprehensive new information on the sizes, molecular structures, and compositions of Fe-rich mobile colloids that accounted for up to 72% of aqueous Fe in anoxic groundwater from a redox-active floodplain. The mobile colloids are multi-phase assemblages consisting of Si-coated ferrihydrite nanoparticles and Fe(II)-OM complexes. Ferrihydrite nanoparticles persisted under both oxic and anoxic/sulfidic conditions, which we attribute to passivation by Si and OM. These findings suggest that mobile Fe-rich colloids generated in floodplains can persist during transport through redox-variable soils and could be discharged to surface waters. These results shed new light on their potential to transport Fe, OM, and nutrients across terrestrial-aquatic interfaces.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Global Biogeochemical Cycle of Rhenium

Here, this paper is the first comprehensive synthesis of what is currently known about the different natural and anthropogenic fluxes of rhenium (Re) on Earth's surface. We highlight the significant role of anthropogenic mobilization of Re, which is an important consideration in utilizing Re in the context of a biogeochemical tracer or proxy. The largest natural flux of Re derives from chemical weathering and riverine transport to the ocean (dissolved = 62 × 10 6 g yr -1 and particulate = 5 × 10 6 g yr -1 ). This review reports a new global average [Re] of 16 ± 2 pmol L -1 , or 10 ± 1 pmol L -1 for the inferred pre-anthropogenic concentration without human impact, for rivers draining to the ocean. Human activity via mining (including secondary mobilization), coal combustion, and petroleum combustion mobilize approximately 560 × 10 6 g yr -1 Re, which is more than any natural flux of Re. There are several poorly constrained fluxes of Re that merit further research, including: submarine groundwater discharge, precipitation (terrestrial and oceanic), magma degassing, and hydrothermal activity. The mechanisms and the main host phases responsible for releasing (sources) or sequestrating (sinks) these fluxes remain poorly understood. This study also highlights the use of dissolved [Re] concentrations as a tracer of oxidation of petrogenic organic carbon, and stable Re isotopes as proxies for changes in global redox conditions.

58 GEOSCIENCES↗
DOE OSTI2020

Evaluation of Gaseous Substrates for Microbial Immobilization of Contaminant Mixtures in Unsaturated Subsurface Sediments

Extensive vadose zone metals and organic contamination remains at many former industrial and defense manufacturing sites, and effective remedial solutions are needed to slow or prevent its migration to groundwater. In this study, the application of gaseous substrates to stimulate microbial respiratory reduction of comingled radioisotopes and nitrate under unsaturated conditions was examined for possible application at the Hanford Site, a former nuclear production facility in southeastern WA, USA. First, screening studies were performed to qualitatively measure the sediment respiratory response to 14 gaseous or volatile organic substrates at two moisture contents, 4% and 8%. Volatile substrates produced the strongest respiratory response, among them were butyrate, pentane, butyl acetate. Ethane and butane were the most effective gaseous substrates but only at 8% water content. Hanford sediment from two waste sites with distinctive chemistries were wetted to 7% moisture content, packed into columns, and treated with ethane or butane. After 4 weeks, columns were then leached to quantify retardation in the mobility of aqueous contaminant concentrations compared to no gas control columns. Treatment with both gases resulted in >80% removal of Cr from the aqueous phase. However, NO3 concentration and a waste sites exposure history to NO3 had a major effect on U and Tc reduction. Incomplete nitrate reduction outcompeted U and Tc in waste site sediments having limited prior exposure to NO3. Conversely, waste site sediments co-contaminated with NO3 were able to achieve highly reduced conditions resulting in complete denitrification of NO3, and delayed leaching of U and Tc. This implied effective reduction of both contaminants to less mobile species. This study demonstrates that unsaturated vadose sediments at Hanford waste sites have the capacity for a sustained respiratory response to gaseous substrate injection, which could potentially be deployed as part of an overall strategy to reduce the flux of long-lived radionuclides to groundwater at Hanford and other legacy waste sites.

Bagwell, Christopher E.↗
DOE OSTI2023

QA/QC-ed Groundwater Level Time Series in PLM-1 and PLM-6 Monitoring Wells, East River, Colorado (2016-2022)

This data set contains QA/QC-ed (Quality Assurance and Quality Control) water level data for the PLM1 and PLM6 wells. PLM1 and PLM6 are location identifiers used by the Watershed Function SFA project for two groundwater monitoring wells along an elevation gradient located along the lower montane life zone of a hillslope near the Pumphouse location at the East River Watershed, Colorado, USA. These wells are used to monitor subsurface water and carbon inventories and fluxes, and to determine the seasonally dependent flow of groundwater under the PLM hillslope. The downslope flow of groundwater in combination with data on groundwater chemistry (see related references) can be used to estimate rates of solute export from the hillslope to the floodplain and river. QA/QC analysis of measured groundwater levels in monitoring wells PLM-1 and PLM-6 included identification and flagging of duplicated values of timestamps, gap filling of missing timestamps and water levels, removal of abnormal/bad and outliers of measured water levels. The QA/QC analysis also tested the application of different QA/QC methods and the development of regular (5-minute, 1-hour, and 1-day) time series datasets, which can serve as a benchmark for testing other QA/QC techniques, and will be applicable for ecohydrological modeling. The package includes a Readme file, one R code file used to perform QA/QC, a series of 8 data csv files (six QA/QC-ed regular time series datasets of varying intervals (5-min, 1-hr, 1-day) and two files with QA/QC flagging of original data), and three files for the reporting format adoption of this dataset (InstallationMethods, file level metadata (flmd), and data dictionary (dd) files).QA/QC-ed data herein were derived from the original/raw data publication available at Williams et al., 2020 (DOI: 10.15485/1818367). For more information about running R code file (10.15485_1866836_QAQC_PLM1_PLM6.R) to reproduce QA/QC output files, see README (QAQC_PLM_readme.docx). This dataset replaces the previously published raw data time series, and is the final groundwater data product for the PLM wells in the East River. Complete metadata information on the PLM1 and PLM6 wells are available in a related dataset on ESS-DIVE: Varadharajan C, et al (2022). https://doi.org/10.15485/1660962. These data products are part of the Watershed Function Scientific Focus Area collection effort to further scientific understanding of biogeochemical dynamics from genome to watershed scales. 2022/09/09 Update: Converted data files using ESS-DIVE’s Hydrological Monitoring Reporting Format. With the adoption of this reporting format, the addition of three new files (v1_20220909_flmd.csv, V1_20220909_dd.csv, and InstallationMethods.csv) were added. The file-level metadata file (v1_20220909_flmd.csv) contains information specific to the files contained within the dataset. The data dictionary file (v1_20220909_dd.csv) contains definitions of column headers and other terms across the dataset. The installation methods file (InstallationMethods.csv) contains a description of methods associated with installation and deployment at PLM1 and PLM6 wells. Additionally, eight data files were re-formatted to follow the reporting format guidance (er_plm1_waterlevel_2016-2020.csv, er_plm1_waterlevel_1-hour_2016-2020.csv, er_plm1_waterlevel_daily_2016-2020.csv, QA_PLM1_Flagging.csv, er_plm6_waterlevel_2016-2020.csv, er_plm6_waterlevel_1-hour_2016-2020.csv, er_plm6_waterlevel_daily_2016-2020.csv, QA_PLM6_Flagging.csv). The major changes to the data files include the addition of header_rows above the data containing metadata about the particular well, units, and sensor description. 2023/01/18 Update: Dataset updated to include additional QA/QC-ed water level data up until 2022-10-12 for ER-PLM1 and 2022-10-13 for ER-PLM6. Reporting format specific files (v2_20230118_flmd.csv, v2_20230118_dd.csv, v2_20230118_InstallationMethods.csv) were updated to reflect the additional data. R code file (QAQC_PLM1_PLM6.R) was added to replace the previously uploaded HTML files to enable execution of the associated code. R code file (QAQC_PLM1_PLM6.R) and ReadMe file (QAQC_PLM_readme.docx) were revised to clarify where original data was retrieved from and to remove local file paths.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Field Test Bed for Vadose-Zone Monitoring Approaches - 20404

In the Central Plateau at the U.S. Department of Energy Hanford Site, a large inventory of contaminants resides in unsaturated sediments within the approximately 100-meter-thick vadose zone, posing a potential continuing risk to groundwater. Vadose zone remedies used to address these contaminants will require performance monitoring to provide feedback during implementation and for long-term verification that remedial action objectives have been met. Passive approaches may also need long-term monitoring to demonstrate that the flux of contaminants from the vadose zone to the groundwater are below thresholds established to meet groundwater protection goals. Collection of physical (e.g., groundwater or sediment) samples is a common method for identifying contaminant concentration distributions, but this approach is limited by the number of locations and the frequency with which data can be collected. In situ vadose zone measurements have evolved over the past few years to include key measurements of water content, soil water pressure, temperature, and chemical concentration. However, the current generation of sensors is designed for relatively short-term use in near-surface soils or sediments. Geophysical methods have been evolving but are also limited in that they have not been designed for the specific long-term vadose zone monitoring needs at the Hanford Site. Overall, monitoring under unsaturated conditions can be difficult due to the need to install and maintain instrumentation over a large area and depth and the need to identify preferential flow pathways due to geologic and chemical heterogeneities over long periods. Thus, a vadose zone monitoring test bed was initiated to address these challenges and identify cost-effective approaches for implementation and postclosure monitoring of the deep vadose zone. The monitoring test bed is expected to provide valuable field-scale information for the design of vadose zone monitoring systems. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Utah FORGE: Geochemical Data for Cold Groundwaters and Produced Geothermal Fluids

Geochemical data for cold groundwaters and produced geothermal fluids around the Utah FORGE site. The data is compiled into four tables in the attached Excel File. Table 1 is a compilation of compositions (anions, cations, weak acids, oxygen, hydrogen, and carbon isotopes) for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 2 is a compilation of noble gas (He, Ne, Ar) and He and Ne isotopic compositions for cold groundwaters and produced geothermal waters in the Milford valley, Utah. Table 3 provides values for calculated advective and diffusive fluxes of helium. Table 4 provides values of calculated subsurface stored heat between the Opal Mound fault and the Utah FORGE site, which are related to volumes of recently solidified magmatic heat sources.

15 GEOTHERMAL ENERGY↗
NASA NTRS2023

Irrigation in the Earth System

Irrigation accounts for ~70% of global freshwater withdrawals and ~90% of consumptive water use, driving myriad Earth system impacts. In this Review, we summarize how irrigation currently impacts key components of the Earth system. Estimates suggest that more than 3.6 million km 2 of currently irrigated land, with hot spots in the intensively cultivated US High Plains, California Central Valley, Indo-Gangetic Basin and northern China. Process-based models estimate that ~2,700 ± 540 km 3 irrigation water is withdrawn globally each year, broadly consistent with country-reported values despite these estimates embedding substantial uncertainties. Expansive irrigation has modified surface energy balance and biogeochemical cycling. A shift from sensible to latent heat fluxes, and resulting land–atmosphere feedbacks, generally reduce regional growing season surface temperatures by ~1–3 °C. Irrigation can ameliorate temperature extremes in some regions, but conversely exacerbates moist heat stress. Modelled precipitation responses are more varied, with some intensive cropping regions exhibiting suppressed local precipitation but enhanced precipitation downstream owing to atmospheric circulation interactions. Additionally, irrigation could enhance cropland carbon uptake; however, it can also contribute to elevated methane fluxes in rice systems and mobilize nitrogen loading to groundwater. Cross-disciplinary, integrative research efforts can help advance understanding of these irrigation–Earth system interactions, and identify and reduce uncertainties, biases and limitations.

Climate sciences↗
DOE OSTI2025

On the emergent scale of bedrock groundwater contribution to headwater mountain streams

We investigated the contribution of bedrock groundwater to streamflow as a function of catchment scale in a headwater stream. Synoptic surveys were conducted during hydrologically important periods of the year using multiple environmental tracers in stream water, soil water, and bedrock groundwater, along a first-order montane stream, in west-central Montana. Sampled analytes included 222 Rn, used to constrain total subsurface flux, and major and minor elements, used in end-member mixing analysis (EMMA) to identify the contributions of soil and bedrock groundwater to the stream. Partitioning between soil-derived and bedrock-derived groundwater was then analyzed as a function of the incremental and accumulated sub-catchment sizes. Radon results indicated that subsurface water contributions accounted for the majority of streamflow at all surveyed times. EMMA results revealed that the bedrock groundwater contribution to streamflow varied between 26% during peak snowmelt and 44% during late summer. Streamflow generation was dominated by soil groundwater contribution along the entire reach, but the bedrock groundwater contribution increased consistently with accumulated sub-catchment size. However, groundwater contributions were not well-correlated with incremental sub-catchment size. The scale at which increased bedrock groundwater discharge can be correlated with sub-catchment size appears to be >1 km 2 for our study. Our results are consistent with a conceptual model where streamflow is predominantly generated by a 3D subsurface nested flow system. Local subsurface heterogeneities control the stream source at local scales but begin to average out at scales >2 km 2 . Our study indicates that, while soil groundwater is the dominant source, bedrock groundwater remains an important and predictable contributor to streamflow throughout the year, even in a snow-dominated, mountainous headwater catchment.

environmental tracers↗
DOE OSTI2022

Interspecific variation in the timing and magnitude of hydraulic redistribution in a forest with distinct water sources

Trees regulate water availability among their rooting strata through a nocturnal, passive transference of water known as hydraulic redistribution (HR). This study investigates differences in HR and groundwater use among common canopy species in longleaf pine (Pinus palustris Mill., Pinaceae) woodlands and explores environmental factors influencing HR. HR was estimated by sap flux of lateral roots and main stems of three mature canopy species (P. palustris, Quercus laevis Walter., Fagaceae and Quercus margarettae Ashe., Fagaceae). We used δ18O and δD of xylem water, soil water, and groundwater to determine water source. Finally, we related HR to environmental factors (Temperature, VWC, VPD) to better understand controls of HR dynamics. Pinus palustris had higher water use than either Quercus species, and also redistributed significantly more water as a nocturnal subsidy. HR fluxes were inversely related with mean nightly temperature and independent of shallow soil moisture. Stable isotope mixing models, based on δ 18 O and δD, indicated that all species have access to groundwater, but utilized shallow soil water in differing amounts when available. In systems with strong water potential gradients among soil strata, any species with access to a groundwater source is likely capable of HR; however, the magnitude of HR varies significantly by species, even among closely related taxa.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling geogenic and atmospheric nitrogen through the East River Watershed, Colorado Rocky Mountains

There is a growing understanding of the role that bedrock weathering can play as a source of nitrogen (N) to soils, groundwater and river systems. The significance is particularly apparent in mountainous environments where weathering fluxes can be large. However, our understanding of the relative contributions of rock-derived, or geogenic, N to the total N supply of mountainous watersheds remains poorly understood. In this study, we develop the High-Altitude Nitrogen Suite of Models (HAN-SoMo), a watershed-scale ensemble of process-based models to quantify the relative sources, transformations, and sinks of geogenic and atmospheric N through a mountain watershed. Our study is based in the East River Watershed (ERW) in the Upper Colorado River Basin. The East River is a near-pristine headwater watershed underlain primarily by an N-rich Mancos Shale bedrock, enabling the timing and magnitude of geogenic and atmospheric contributions to watershed scale dissolved N-exports to be quantified. Several calibration scenarios were developed to explore equifinality using >1600 N concentration measurements from streams, groundwater, and vadose zone samples collected over the course of four years across the watershed. When accounting for recycling of N through plant litter turnover, rock weathering accounts for approximately 12% of the annual dissolved N sources to the watershed in the most probable calibration scenario (0–31% in other scenarios), and 21% (0–44% in other scenarios) when considering only “new” N sources (i.e. geogenic and atmospheric). On an annual scale, instream dissolved N elimination, plant turnover (including cattle grazing) and atmospheric deposition are the most important controls on N cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Nitrate, ammonium, and DON mass time series output for East River stream, vadose zone and groundwater subwatersheds from HAN-SoMo model

There is a growing understanding of the role that bedrock weathering can play as a source of nitrogen (N) to soils, groundwater and river systems. The significance is particularly apparent in mountainous environments where weathering fluxes can be large. However, our understanding of the relative contributions of rock-derived, or geogenic, N to the total N supply of mountainous watersheds remains poorly understood. This dataset contains model output from the High-Altitude Nitrogen Suite of Models (HAN-SoMo), a watershed-scale ensemble of process-based models to quantify the relative sources, transformations, and sinks of geogenic and atmospheric N through a mountain watershed. Our study is based in the East River Watershed (ERW) in the Upper Colorado River Basin. The East River is a near-pristine headwater watershed underlain primarily by an N-rich Mancos Shale bedrock, enabling the timing and magnitude of geogenic and atmospheric contributions to watershed scale dissolved N-exports to be quantified. Six calibration scenarios were developed to explore equifinality using >1600 N concentration measurements from streams, groundwater, and vadose zone samples collected over the course of four years across the watershed, from 2014-2018. The files in this dataset are named according to which calibration scenario they represent (1, 2, 3, NM1, NM2, and NC), and include daily time series model output of nitrate, ammonium and dissolved organic nitrogen concentrations for all six calibration scenarios described in Maavara et al., 2021, "Modeling geogenic and atmospheric nitrogen through the East River Watershed, Colorado Rocky Mountains," PLOS One (in press), for the stream, vadose zone, and groundwater in the five discretized subwatersheds: Rustlers Gulch, Copper Creek, East above Quigley (EAQ), the Middle East River (ME), and Lower Triangle (LT) from Gothic to the Pumphouse. The READ_ME.txt file includes details on units and naming conventions in each of these files. Model files are in .mat format, to be opened in MATLAB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Controls on the 36 Cl/Cl input ratio of paleo-groundwater in arid environments: New evidence from 81 Kr/Kr data

In this study, measurements of the long-lived 81 Kr and 36 Cl radioisotopes in groundwater from the Negev Desert (Israel) were used to assess the 36 Cl/Cl input ratios and Cl – contents for paleorecharge into the Nubian Sandstone Aquifer (NSA). The reconstructed Cl – content of the recharge flux was on the order of 300–400 mg/L. An initial 36 Cl/Cl ratio of 50 × 10 –15 was assessed for the groundwater replenishment in the Negev Desert since the late Pleistocene, in agreement with the 36 Cl/Cl ratios in recent local rainwater. This is despite possible changes in the climatic conditions and the 36 Cl production rates in the atmosphere over this timeframe. This similarity in values is explained by the major role played by the erosion and weathering of near-surface materials in the desert environment that dominate the hydrochemistry of rains, floods, and the consequent groundwater recharge. Spatial variation in the reconstructed initial 36 Cl/Cl ratio is accounted for by the differences in the mineral aerosol sources for specific recharge areas of the NSA. Accordingly, regional variations in the 36 Cl/Cl input in groundwater reservoirs surrounding the Mediterranean Sea indicate various processes that govern the 36 Cl/Cl system. Finally, the results of this study highlight the great advantage of integrating 81 Kr age information in evaluating the initial 36 Cl/Cl and Cl – input, which is essential for the calibration of 36 Cl radioisotope as an available long-term dating tool for a given basin.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Dominant Controls on Preferential Flow and Their Implications for Future Soil Water Fluxes

Abstract Soil water flow, particularly preferential flow (PF), is a critical control on hydrological and biogeochemical processes, including groundwater recharge, contaminant transport, and carbon cycling. However, it remains challenging to predict PF occurrence across large environmental gradients. Here, we developed a deep learning (DL) model to estimate event‐scale soil water flow velocity and the probability of PF occurrence using high‐frequency soil moisture and precipitation data from 33 sites across the National Ecological Observatory Network. The model demonstrated high skill in predicting the binary occurrence of PF (91% F1‐score; 85% accuracy) but the performance was limited in predicting soil water velocity ( R 2 = 0.31). We found that precipitation characteristics (duration, volume, and intensity) were the most important predictors for soil water velocity. Among the non‐precipitation event variables, sand content showed relatively high predictive skill, though differences among non‐event climate variables were generally modest. Lower sand content was associated with increased predicted soil water velocity, a finding that highlights the role of soil structure in producing more non‐uniform flow, which contrasts with traditional uniform flow models. Projecting a reduced DL model under both moderate and high‐emissions future climate scenarios (2060–2099 Representative Concentration Pathways 4.5 and 8.5), we found ∼7.3% increase under RCP4.5 and ∼15% under RCP8.5 of soil water velocities compared to the historical simulation, while modeled likelihood of PF changed little. These findings suggest climate change is not making PF more frequent, but it is making existing PF pathways more efficient with important consequences for associated nutrient and contaminant transport under climate change. Plain Language Summary Water movement in soil is critical for water quality. While often modeled as a uniform flow process, in reality water moves rapidly through cracks and burrows in what is called “preferential flow” (PF), which limits natural filtration and can transport pollutants. We developed a deep learning model, trained on data from 33 U.S. sites, to predict when and how fast this PF occurs based on precipitation, soil, and climate data. The model showed that precipitation characteristics (duration, intensity, volume) were the most important predictors of PF. Lower soil sand content/higher clay content was associated with faster water flow, likely due to clay soils forming aggregates and cracks that water moves through rather than infiltrating uniformly. Further analyses based on climate projections suggest that the speed at which PF occurs will become more rapid under future climate scenarios compared to historical simulation. This highlights the need to represent PF in soil water models when assessing future water quality. Key Points The effect of precipitation peak intensity on soil water velocities declined with increasing precipitation intensity Antecedent soil moisture failed to predict preferential flow (PF), contrasting the high predictive power of sand content Climate predictions suggest that soil water velocities through PF paths will increase ∼15% by 2099

Li, Bonan↗
DOE OSTI2025

Quantifying Groundwater Response and Uncertainty in Beaver‐Influenced Mountainous Floodplains Using Machine Learning‐Based Model Calibration

Abstract Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model‐data integration workflow to address uncertainty in groundwater response to beaver‐induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning‐based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down‐valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down‐valley underflow, from 2 during dry pond periods to 20 during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down‐valley underflow, minimizes the effect of beaver‐induced inundation on water quality. We emphasize the need for field‐scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response. Plain Language Summary Beavers change the flow of water in river corridors by creating ponds, expanding wetlands, and flooding floodplains. This increases surface water area, promotes plant growth, and enhances biodiversity. However, the impact of this flooding on groundwater flow is not well understood, especially in mountainous areas with gravel layers where water moves easily beneath soil. In this study, we used numerical modeling to investigate how beaver ponds influence groundwater in a mountainous floodplain of the Upper Colorado River Basin. We adapted a machine learning method to validate our numerical models using multiple field data sets. Our findings show that beaver ponds significantly increase vertical water flow from the soil to the gravel during wet periods, compared to when the ponds are fully drained. The study also highlights the importance of floodplain structure in controlling both water flow in gravel layers along the river direction and vertical flow from the soil to the gravel with the presence of beavers. To reduce uncertainty in groundwater response, we emphasize the need for more field‐scale measurements of floodplain structure, hydraulic properties, and evapotranspiration changes. Key Points Floodplain structures and hydraulic conductivities are important for groundwater response with beaver ponds in mountainous floodplains Large down‐valley underflow in permeability‐stratified floodplains reduces beaver‐induced impacts on groundwater storage and water quality Machine learning‐based model calibration methods are effective for estimating posterior distributions of groundwater model parameters

Wang, Lijing↗
DOE OSTI2024

Old-Aged groundwater contributes to mountain hillslope hydrologic dynamics

Understanding connectivity between the soil and deeper bedrock groundwater is needed to accurately predict a watershed’s response to perturbation, such as drought. Yet, the bedrock groundwater dynamics in mountainous environments are typically under-constrained and excluded from watershed hydrologic models. Here, we investigate the role of groundwater characterized with decadal and longer water ages on the hydrologic and mass-transport processes within a steep snow-dominated mountain hillslope in the Central Rocky Mountains (USA). We quantify subsurface and surface water mass-balance, groundwater flowpaths, and age distributions using the ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models, which are compared to hydrometric and environmental tracer observations. An ensemble of models with varied soil and hydrogeologic parameters reproduces observed groundwater levels and century-scale mean ages inferred from environmental tracers. The numerical models suggest soil water near the toe of the hillslope contains considerable (>60 % of the mass-flux) contributions from bedrock flowpaths characterized with water ages >10 years. Flowpath connectivity between the deeper bedrock and soil systems is present throughout the year, highlighting the potentially critical role of groundwater with old ages on processes such as evapotranspiration and streamflow generation. The coupled numerical model and groundwater age observations show the bedrock groundwater system influences the hillslope hydrodynamics and should be considered in mountain watershed conceptual and numerical models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2018

Data Assimilation of Terrestrial Water Storage to Adjust Precipitation Fluxes

The Gravity Recovery and Climate Experiment (GRACE) mission has provided unprecedented observations of terrestrial water storage (TWS) dynamics at basin to continental scales. TWS is defined as the sum of groundwater, soil moisture, snow, surface water, ice and biomass water. Data assimilation of GRACE TWS observations has been shown to improve simulation of groundwater, streamflow, and snow water equivalent, and has also proven useful for drought monitoring and identifying human impacts on the water cycle. From a modeling perspective, the TWS components are defined as "prognostic hydrological states". Existing GRACE data assimilation schemes update these prognostic states directly. In this work, we propose an alternate approach in which precipitation fluxes are adjusted in order to achieve the desired change in the hydrological prognostic states. Limitations of such an approach include the assumption that all errors in TWS originate from errors in precipitation. Nonetheless, benefits comprise (1) the water balance is maintained, as opposed to having to add increments to the water budget components, (2) the model automatically determines how to distribute the updates among the TWS prognostic states, and (3) it is not necessary to know the exact time of the observation TWS, because the TWS change timing is determined by the precipitation forcing.

Girotto, Manuela↗