Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydraulic Conductivity”

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 37 records · Page 2

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↗

Agricultural soil moisture experiment, Colby, Kansas 1978: Measured and predicted hydrological properties of the soil

Predictive procedures for developing soil hydrologic properties (i.e., relationships of soil water pressure and hydraulic conductivity to soil water content) are presented. Three models of the soil water pressure-water content relationship and one model of the hydraulic conductivity-water content relationship are discussed. Input requirements for the models are indicated, and computational procedures are outlined. Computed hydrologic properties for Keith silt loam, a soil typer near Colby, Kansas, on which the 1978 Agricultural Soil Moisture Experiment was conducted, are presented. A comparison of computed results with experimental data in the dry range shows that analytical models utilizing a few basic hydrophysical parameters can produce satisfactory data for large-scale applications.

Arya, L. M.↗

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↗

Temporal persistence of postfire flood hazards under present and future climate conditions in southern Arizona, USA

Changes to soil hydraulic properties that reduce infiltration capacity following fire can increase flash flood risks. These risks are exacerbated by rainfall intensification associated with a warming climate. However, the potential effects of climate-change-driven rainfall intensification on postfire floods remain largely unexplored. Using rainfall and runoff observations from a 49.4 km 2 watershed in southern Arizona, USA, and a hydrologic model (KINEROS2), we examined the temporal evolution following a historic fire of three crucial hydrologic parameters: soil saturated hydraulic conductivity (K sp ), net capillary drive (G p ), and hydraulic roughness (n c ). We explored how the effect of fire on these parameters may influence peak flow under future climate scenarios derived from CMIP6, specifically the medium emissions scenario (SSP2-4.5) and high emissions scenario (SSP5-8.5). Results demonstrate an increase in K sp from 11 mm h −1 in the first postfire year to 60 mm h −1 in postfire year 3. G p similarly increased from 19 mm in the first postfire year to 30 mm in the third, while nc was relatively constant. The highest simulated Q p occurred in the first postfire year. Under the SSP2-4.5 scenario, the likelihood of a 100-year flood is projected to be twice as large by the middle of the century relative to its historical magnitude. Simulations further indicate that the maximum expected discharge associated with a postfire flood, as derived from historical data, could be triggered by a 10-year rainstorm under the SSP5-8.5 scenario by the late century. Simulations also demonstrate that rainfall intensification will lead to greater persistence of elevated flood hazards following fire by the late century under both the SSP2-4.5 and the SSP5-8.5 scenarios.

54 ENVIRONMENTAL SCIENCES↗

Hydraulic constraints to stomatal conductance in flooded trees

Stomatal closure is a pervasive response among trees exposed to flooded soil. We tested whether this response is caused by reduced hydraulic conductance in the soil-to-leaf hydraulic continuum (k total ), and particularly by reduced root hydraulic conductance (k root ), which has been widely hypothesized. We tracked stomatal conductance at the leaf level (g s ) and canopy scale (G s ) along with physiological conditions in two temperate tree species, Magnolia grandiflora and Quercus virginiana, that were subjected to flood and control conditions in a greenhouse experiment. Flooding reduced g s , G s , k root and k total . Path analysis showed strong support for direct effects of k total on g s and for flood duration on k total , but not k root on k total . A process-based model that accounted for the k total reduction predicted the timeseries of G s in flood and control treatment trees reasonably well (predicted versus observed G s R 2 = 0.80 and 0.51 for M. grandiflora and Q. virginiana, respectively). However, accounting only for k root reduction in flooded trees was insufficient for predicting observed G s reduction. Together, these results suggest that hydraulic constraints were not limited to roots and highlight the need to account for flooding effects on k total when projecting forest ecosystem function using process-based models.

Plant stress↗

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗

Spatial variation of corn canopy temperature as dependent upon soil texture and crop rooting characteristics

A soil plant atmosphere model for corn (Zea mays L.) together with the scaling theory for soil hydraulic heterogeneity are used to study the sensitivity of spatial variation of canopy temperature to field averaged soil texture and crop rooting characteristics. The soil plant atmosphere model explicitly solves a continuity equation for water flux resulting from root water uptake, changes in plant water storage and transpirational flux. Dynamical equations for root zone soil water potential and the plant water storage models the progressive drying of soil, and day time dehydration and night time hydration of the crop. The statistic of scaling parameter which describes the spatial variation of soil hydraulic conductivity and matric potential is assumed to be independent of soil texture class. The field averaged soil hydraulic characteristics are chosen to be representative of loamy sand and clay loam soils. Two rooting characteristics are chosen, one shallow and the other deep rooted. The simulation shows that the range of canopy temperatures in the clayey soil is less than 1K, but for the sandy soil the range is about 2.5 and 5.0 K, respectively, for the shallow and deep rooted crops.

Choudhury, B. J.↗

Water Relations, Gas Exchange, and Nutrient Response to a Long Term Constant Water Deficit

Wheat plants (Triticum aestivum) were grown for 43 days in a micro-porous tube nutrient delivery system. Roots were unable to penetrate the microporous tube, but grew on the surface and maintained capillary contact with the nutrient solution on the inside of the tube through the 5-micron pores of the porous tube. Water potential in the system was controlled at -0.4, -0.8, and -3.0 kPa by adjusting the applied pressure (hydrostatic head) to the nutrient solution flowing through the microporous tubes. A relatively small decrease in applied water potential from -0.4 to -3.0 kPa resulted in a 34% reduction of shoot growth but only a moderate reduction in the midday leaf water potential from -1.3 to -1.7 MPa. Carbon dioxide assimilation decreased and water use efficiency increased with the more negative applied water potentials, while intercellular CO2 concentration remained constant. This was associated with a decrease in stomatal conductance to water vapor from 1.90 to 0.98 mol/(sq m sec) and a decrease in total apparent hydraulic conductance from 47 to 12 (micro)mol/(sec MPa). Although the applied water potentials were in the -0.4 to -3.0 kPa range, the actual water potential perceived by the plant roots appeared to be in the range of -0.26 to -0.38 MPa as estimated by the leaf water potential of bagged plants. The amount of K, Ca, Mg, Zn, Cu, and B accumulated with each unit of transpired water increased as the applied water potential became less negative. The increase in accumulation ranged from 1.4-fold for K to 2.2-fold for B. The physiological responses observed in this study in response to small constant differences in applied water potentials were much greater than expected from either the applied water potential or the observed plant water potential. Even though the micro-porous tube may not represent natural conditions and could possibly introduce morphological and physiological artifacts, it enables a high degree of control of water potential that facilitates the investigation of many aspects of water relations not practical with other experimental systems.

Berry, Wade L.↗

Thermal Hydraulic Disaggregation of SMAP Soil Moisture Over the Continental United States

Thermal Hydraulic disaggregation of Soil Moisture (THySM) algorithm was implemented to downscale NASA’s Soil Moisture Active Passive (SMAP) Enhanced soil moisture (SM) product to 1 km over the continental United States (CONUS). This algorithm was developed by combining thermal inertia theory with a soil hydraulic-based approach that considers fine-scale SM spatial distribution driven by both heat fluxes and hydraulic conductivity in soils. Relative soil wetness values were estimated using land surface temperature and normalized difference vegetation index for the thermal inertia model and using soil properties for the hydraulic model. The relative soil wetness values at 1 km from both models were then combined by using weighting functions whereby the spatial distribution of SM was governed more by thermal fluxes during times of strong heat transport and infiltration during moisture abundant soil conditions. THySM values were evaluated using in situ SM measurements from SMAP Core Validation Sites (CVS), the USDA Soil Climate Analysis Network, and the NOAA Climate Reference Network over CONUS. THySM shows higher accuracy than the SMAP / Sentinel-1 (SPL2SMAP_S) 1 km SM product when compared to in situ measurements. The accuracy of THySM is 0.048 m3/m3 based on unbiased root mean square error (ubRMSE), outperforming SPL2SMAP_S by 0.01-0.02 m3/m3. The ubRMSE of THySM 1km SM over the SMAP grassland/rangeland-dominated CVS sites is better than 0.04 m3/m3, which meets the SMAP mission SM accuracy requirement applied at 9 and 36 km.

SMAP↗

Near-term lander experiments for growing plants on Mars: requirements for information on chemical and physical properties of Mars regolith

In order to support humans for long-duration missions to Mars, bioregenerative Advanced Life Support (ALS) systems have been proposed that would use higher plants as the primary candidates for photosynthesis. Hydroponic technologies have been suggested as the primary method of plant production in ALS systems, but the use of Mars regolith as a plant growth medium may have several advantages over hydroponic systems. The advantages for using Mars regolith include the likely bioavailability of plant-essential ions, mechanical support for plants, and easy access of the material once on the surface. We propose that plant biology experiments must be included in near-term Mars lander missions in order to begin defining the optimum approach for growing plants on Mars. Second, we discuss a range of soil chemistry and soil physics tests that must be conducted prior to, or in concert with, a plant biology experiment in order to properly interpret the results of plant growth studies in Mars regolith. The recommended chemical tests include measurements on soil pH, electrical conductivity and soluble salts, redox potential, bioavailability of essential plant nutrients, and bioavailability of phytotoxic elements. In addition, a future plant growth experiment should include procedures for determining the buffering and leaching requirements of Mars regolith prior to planting. Soil physical tests useful for plant biology studies in Mars regolith include bulk density, particle size distribution, porosity, water retention, and hydraulic conductivity.

NASA Center JSC↗

Air Sparge Pilot Study in the DNAPL Source Zone at Launch Complex 34

Multiple releases of trichloroethene (TCE) occurred at Launch Complex 34 (LC34) between the late 1950s and 1968. A 2007 conceptual site model estimated a 2-acre dense non-aqueous phase liquid (DNAPL) source area with mass in excess of 90,000 pounds, nearly 40 years after termination of launch activities. Located on a barrier island, LC34 currently has no complete exposure pathways and a historical groundwater flow radial from the DNAPL source zone (DSZ), the focus area of this study. The geology at LC34 is classified as Layers 1 through 9 with each layer representing a different lithology. These lithologies contribute to large variations in hydraulic conductivity (1x10-3 cm/sec to 1x10-8 cm/sec) with notable fine-grained units at Layer 4 (sandy clay) and Layer 7 (fine silty sand). Historically, technologies implemented at LC34 have been split vertically based on technology limitations and lithology. Remedial technologies have been evaluated to control and/or remediate the DSZ in the past; however, more aggressive technologies required significant cost and would likely leave considerable mass. Thus, an adaptive site management strategy has been implemented that adopts a treatment train approach which began in 2009 with hydraulic containment via pump and treat to control mass discharge from the DSZ while removing mass as a secondary benefit. Concurrently, hot spot areas in the larger dissolved plume are being treated by air sparging to reduce overall mass and the plume footprint. In 2019, a re-characterization of the DSZ was completed to update conditions in support of implementing more aggressive technologies as part of the treatment train approach. The results showed that the DSZ remained relatively the same size but with a slightly different morphology. Moreover, data showed that the fine-grained units (Layers 4 and 7) are storing most of the remaining mass, with TCE concentrations suggestive of DNAPL extending into Layer 7 (approximately 80 to 100 ft below land surface). Based on the current conditions, more aggressive technologies are still cost prohibitive; therefore, a different approach is warranted to determine the next implementable step in the treatment train. Air sparging is being proposed as a technology alternative based on several factors including lower treatment costs. A pilot study was conducted to test the feasibility of air sparging in the DSZ as the next step in the treatment train for that area.

Launch Complex 34↗

Climatic Versus Biotic Constraints on Carbon and Water Fluxes in Seasonally Drought-affected Ponderosa Pine Ecosystems

We investigated the relative importance of climatic versus biotic controls on gross primary production (GPP) and water vapor fluxes in seasonally drought-affected ponderosa pine forests. The study was conducted in young (YS), mature (MS), and old stands (OS) over 4 years at the AmeriFlux Metolius sites. Model simulations showed that interannual variation of GPP did not follow the same trends as precipitation, and effects of climatic variation were smallest at the OS ( 50%), and intermediate at the YS (<20%). In the young, developing stand, interannual variation in leaf area has larger effects on fluxes than climate, although leaf area is a function of climate in that climate can interact with age-related shifts in carbon allocation and affect whole-tree hydraulic conductance. Older forests, with well-established root systems, appear to be better buffered from effects of seasonal drought and interannual climatic variation. Interannual variation of net ecosystem exchange (NEE) was also lowest at the OS, where NEE is controlled more by interannual variation of ecosystem respiration, 70% of which is from soil, than by the variation of GPP, whereas variation in GPP is the primary reason for interannual changes in NEE at the YS and MS. Across spatially heterogeneous landscapes with high frequency of younger stands resulting from natural and anthropogenic disturbances, interannual climatic variation and change in leaf area are likely to result in large interannual variation in GPP and NEE.

Schwarz, P. A.↗

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 2. Evaluating Controls on Flow Persistence in an Urbanized Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in an urbanized catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, distributed temperature sensing (DTS), continuous self-potential (SP) monitoring, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Field_Application subfolder contains the ATS XML input scripts, data files, output data for the SP site. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. The flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.m can only be used with COMSOL with MATLAB) is executed using the ATS output data to simulate the potential field. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) DTS Contains collated DTS data including raw Stokes and anti-Stokes measurement (provided as .h5 file). It also includes DTS processing.ipynb, a Jupyter notebook for calibrating the DTS data using dts_calibration Python package. cooler_calibration.csv is the DTS calibration CSV used in the calibration sequence. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion. 6) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 7) SP Contains the SP data collected in field at the SP sites (provided as CSV files). 8) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). Note: Code files (.ipynb, .py, .xml) can be opened in any standard code editor, .exo file can be viewed using Paraview, .h5 files can be opened using HDFView software and h5py Python package, and .resipy file can be opened with the open-source ResIPy software.

ATS↗

JAMES BUTTLE REVIEW: Interflow, subsurface stormflow and throughflow: A synthesis of field work and modelling

Interflow, throughflow and subsurface stormflow are interchangeable terms that refer to the lateral subsurface flow above a restricting layer of lower hydraulic conductivity that occurs during and following storm events. Interflow (used here) is a more dominant process in steeper catchments with high infiltration capacity soils overlying a more impermeable soil or geologic layer. Interflow as a runoff process was first recognised in the early 1900s, yet hydrologists still struggle to predict its occurrence, persistence, importance, interaction with other streamflow generation processes, and potential to connect to valleys and streams during and following storms. We review the history of interflow research and address some of the challenges in understanding its role in runoff production. We argue that characterising the controls on interflow initiation and occurrence relies on detailed field observations of subsurface properties, which exist only in limited experimental settings. This data shortcoming contributes to our inability to predict interflow or determine its contribution to streamflow more broadly. There remain many opportunities to advance our understanding of interflow that include both modelling and experimental or observational approaches in hydrology.

hillslope hydrology↗

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗

The Role of Bedrock Circulation Depth and Porosity in Mountain Streamflow Response to Prolonged Drought

Quantitative understanding is lacking on how the depth of active groundwater circulation in bedrock affects mountain streamflow response to a multi-year drought. We use an integrated hydrological model to explore the sensitivity of a variety of streamflow metrics to bedrock circulation depth and porosity under a plausible extreme drought scenario lasting up to 5 years. Endmember depth versus hydraulic conductivity relationships and porosity values for fractured crystalline rock are simulated. With drought, a deeper circulation system with higher drainable porosity more effectively buffers minimum flow and significantly limits perennial stream loss in comparison to a shallow circulation system. Streamflow buffering is accomplished through extensive groundwater storage loss. However, deeper circulation systems experience prolonged recovery from drought in comparison to storage-limited shallow systems. Research highlights the importance of characterizing the deeper bedrock hydrogeology in mountainous watersheds to better understand and predict drought impacts on stream ecosystem health and water resource sustainability.

54 ENVIRONMENTAL SCIENCES↗

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 2. Quantitative Moisture Dynamics Estimation Model

Abstract The long‐term containment of high‐level radioactive waste in geological disposal repositories relies on Engineered Barrier Systems (EBS), with bentonite clay emerging as a candidate material due to its unique properties. Understanding moisture dynamics within bentonite buffers is crucial for EBS performance, as it directly influences the material's swelling capacity, thermal and hydraulic conductivity, mechanical properties, and long‐term evolution under complex thermal‐hydrological‐mechanical (THM) processes. This study develops an advanced Electrical Resistivity Tomography (ERT)‐based framework to quantitatively monitor moisture dynamics under THM conditions. Our framework extends the Waxman‐Smits model to incorporate the coupled effects of temperature, water content, fluid chemistry, and mechanical changes on bentonite's electrical properties. Utilizing HotBENT‐Lab data from our companion paper, which includes electrical conductivity, CT density, and thermocouple measurements, this study offers a novel methodological framework bridging different scales of the model. Our results show that the extended model can estimate water content from ERT data, capturing spatial and temporal variations in moisture distribution within bentonite columns. However, the model tends to overestimate water content compared to CT density‐derived measurements. We address this discrepancy by incorporating a simplified swelling effect model, which improves agreement between ERT and CT density‐based water content estimates. We also discuss model limitations, including simplified treatment of swelling and micropore effects, and propose a conceptual framework for transitioning from laboratory to field applications, addressing challenges such as parameter scalability, field validation methods, and integration of diverse data sources. This ERT‐based framework can potentially advance real‐world moisture monitoring of bentonite‐based EBS in nuclear waste repositories. Plain Language Summary Safely containing high‐level radioactive waste depends on barriers made from materials like bentonite clay, which is effective because it swells and seals in the waste. To ensure these barriers work well over time, it's important to understand how moisture moves through the clay. Our study developed a new method using ERT to monitor moisture levels in bentonite under conditions that mimic those in actual storage sites, including changes in temperature, water content, and mechanical stress. This study improved an existing model to better account for how these factors affect the clay, allowing us to create more accurate moisture maps. Initially, the proposed model overestimated the amount of water in the clay, but its accuracy was improved by factoring in how the clay swells when wet. This study also identified some limitations of the model and suggested ways to adapt it for use in real‐world waste storage sites. This new approach could lead to better monitoring and safety checks for nuclear waste storage systems, helping to ensure long‐term containment. Key Points This work develops an ERT‐based framework extending the Waxman‐Smits model to monitor bentonite moisture dynamics during coupled THM processes The extended model accurately estimates water content from Electrical Resistivity Tomography data, incorporating swelling effects to improve precision This work proposes a conceptual framework for transitioning from laboratory to field applications, advancing EBS monitoring in nuclear waste repositories

Chen, Hang↗