Search NASA⌕ Search

SEARCH · Search NASA

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 91 records · Page 5

Spatial Mapping of Riverbed Grain-Size Distribution Using Machine Learning

Recent alluvial sediments in riverbeds play a significant role in controlling hydrologic exchange flows (HEFs) in river systems. The alluvial layer is usually associated with strong heterogeneity in physical properties (e.g., permeability and hydraulic conductivity), which affects local HEFs and therefore biogeochemical processes. The spatial distribution of these physical properties needs to be determined to inform the numerical models used to reveal the realistic hydro-biogeochemical behaviors. Such information can be obtained based on the intrinsic link between sediment grain-size distribution and hydraulic properties where sediment texture information is available. However, grain-size measurements are usually spatially sparse and do not have adequate coverage and resolution, particularly for a relatively large domain such as the Hanford Reach of the Columbia River. In this paper, we adopted machine learning (ML) approaches for categorizing and mapping the spatial distributions of riverbed substrate grain size and filling in missing areas of substrate data using the ML models along the reach. Such ML models for substrate size mapping were trained at 13,372 locations using measured substrate sizes along with observed and simulated attributes, including bathymetric attributes (e.g., elevation, slope, and aspect ratio) from LIDAR and bathymetric surveys, and hydrodynamic properties (e.g., water depth, velocity, shear stress, and their statistical moments). An ensemble bagging-based ML technique, Random Forest, was adopted to identify the most influential factors as predictors to develop the predictive models with over-fitting issues addressed. The models were evaluated with respect to each individual substrate size class and the lumped group, and then used to generate the final substrate size maps covering all the grid cells in the numerical modeling domain.

54 ENVIRONMENTAL SCIENCES↗

Single-Well Push–Pull Tracer Test Analyses to Determine Aquifer Reactive Transport Parameters at a Former Uranium Mill Site (Grand Junction, Colorado)

At a former uranium mill site where tailings have been removed, prior work has determined several potential ongoing secondary uranium sources. These include locations with uranium sorbed to organic carbon, uranium in the unsaturated zone, and uranium associated with the presence of gypsum. To better understand uranium mobility controls at the site, four single-well push–pull tests (with a drift phase) were completed with the goal of deriving aquifer flow and contaminant transport parameters for inclusion in a future sitewide reactive transport model. This goes beyond the traditional use of a constant sorption distribution coefficient (K d ) and allows for the evaluation of alternative remedial injection fluids, which can produce variable K d values. Dispersion was first removed from the resulting data to determine possible reactions before conducting reactive transport simulations. These initial analyses indicated the potential need to include cation exchange, uranium sorption, and gypsum dissolution. A reactive transport model using multiple layers to account for partially penetrating wells was completed using the PHT-USG reactive transport modeling code and calibrated using PEST. The model results quantify the hydraulic conductivity and dispersion parameters using the injected tracer concentrations. Uranium sorption, cation exchange, and gypsum dissolution parameters were quantified by comparing the simulated versus observed geochemistry. All simulations required some cation exchange and calcite equilibrium, and one simulation required gypsum dissolution to improve the model fit for calcium and sulfate. Uranium sorption parameters were not strongly influenced by the other parameter values but were highly influenced by uranium concentrations during the drift phase, with possible kinetic rate limitations. Thus, a future recommendation for such push–pull tests is to collect more geochemical data during the drift phase. The final uranium sorption parameters were within the range of values determined from prior column testing. The flow and transport parameters derived from these single-well push–pull tests will provide initial parameters for any future sitewide reactive transport model.

54 ENVIRONMENTAL SCIENCES↗

The uncharacterized gene EVE contributes to vessel element dimensions in Populus

The radiation of angiosperms led to the emergence of the vast majority of today’s plant species and all our major food crops. Their extraordinary diversification occurred in conjunction with the evolution of a more efficient vascular system for the transport of water, composed of vessel elements. The physical dimensions of these water-conducting specialized cells have played a critical role in angiosperm evolution; they determine resistance to water flow, influence photosynthesis rate, and contribute to plant stature. However, the genetic factors that determine their dimensions are unclear. Here we show that a previously uncharacterized gene, ENLARGED VESSEL ELEMENT ( EVE ) , contributes to the dimensions of vessel elements in Populus , impacting hydraulic conductivity. Our data suggest that EVE is localized in the plasma membrane and is involved in potassium uptake of differentiating xylem cells during vessel development. In plants, EVE first emerged in streptophyte algae, but expanded dramatically among vessel-containing angiosperms. The phylogeny, structure and composition of EVE indicates that it may have been involved in an ancient horizontal gene-transfer event.

59 BASIC BIOLOGICAL SCIENCES↗

Hydraulically‐vulnerable trees survive on deep‐water access during droughts in a tropical forest

Summary Deep‐water access is arguably the most effective, but under‐studied, mechanism that plants employ to survive during drought. Vulnerability to embolism and hydraulic safety margins can predict mortality risk at given levels of dehydration, but deep‐water access may delay plant dehydration. Here, we tested the role of deep‐water access in enabling survival within a diverse tropical forest community in Panama using a novel data‐model approach. We inversely estimated the effective rooting depth (ERD, as the average depth of water extraction), for 29 canopy species by linking diameter growth dynamics (1990–2015) to vapor pressure deficit, water potentials in the whole‐soil column, and leaf hydraulic vulnerability curves. We validated ERD estimates against existing isotopic data of potential water‐access depths. Across species, deeper ERD was associated with higher maximum stem hydraulic conductivity, greater vulnerability to xylem embolism, narrower safety margins, and lower mortality rates during extreme droughts over 35 years (1981–2015) among evergreen species. Species exposure to water stress declined with deeper ERD indicating that trees compensate for water stress‐related mortality risk through deep‐water access. The role of deep‐water access in mitigating mortality of hydraulically‐vulnerable trees has important implications for our predictive understanding of forest dynamics under current and future climates.

54 ENVIRONMENTAL SCIENCES↗

Estimating Vadose Zone Flow Properties at the 100 K-East Soil Flushing Site Using ERT Monitoring Data: 2023 Interim Report - 100 KE Soil Property Estimation

In situ soil flushing is being using at the Hanford 100 K-East (100 KE) area to transport mobile chromium contamination in the vadose zone to the water table, where it can be collected and treat through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. In practice, it is infeasible to comprehensively monitor which regions of the vadose zone are being infiltrated through direct sampling of pore water. Consequently, there can be significant uncertainty about which regions of the subsurface have been treated, especially if hydrogeologic conditions are favorable for the development of unstable flows and preferred flow pathways through the vadose zone (Jarvis, Koestel, and Larsbo 2016). Current approaches for quantitative monitoring of soil flushing performance rely on contaminant concentration measurements collected from extractions wells. There is no quantitative information on the volume of flush water delivered to targeted regions of the vadose zone at the Hanford Site, leading to significant uncertainty regarding source term removal and long-term impacts to groundwater. If the subsurface hydrogeologic properties at the 100 KE Area were adequately known, qualitative metrics of soil flushing performance could be simulated, thereby negating expenses required to obtain quantitative performance information through borehole drilling/sampling. However, estimating in situ hydrogeologic properties has long proven elusive, due primarily to a lack of sufficient information to constrain heterogeneous property estimates to a useful degree of certainty. Estimating vadose zone hydrogeologic properties is particularly challenging due the dependence of hydraulic conductivity on saturation. This report describes progress toward a first-of-its-kind demonstration using surface time-lapse 3D electrical resistivity tomography (ERT) monitoring data to estimate the hydrogeologic properties that control flush water transport at the 100 KE soil flushing site. The ultimate objective is (1) to verify sufficient information exists in the ERT monitoring data to adequately resolve vadose zone hydraulic properties, and (2) generate a “digital twin” (i.e., a numerical simulator) that can be used to simulate the amount of flush water that has been delivered to each targeted region of the vadose zone, and thereby assess the efficacy of flush water delivery. Resulting performance estimates can be used in leu of comprehensive borehole drilling and direct sampling (or wellbore logging) that would otherwise be required to obtain the same information.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Measuring a Low Horizontal Hydraulic Gradient in a High Transmissivity Aquifer

Abstract High transmissivity aquifers typically have low hydraulic gradients (i.e., a flat water table). Measuring low gradients using water levels can be problematic because measurement error may be greater than the true difference in water levels (i.e., a low signal‐to‐noise ratio). In this study, the feasibility of measuring a hydraulic gradient in the range of 10 −6 to 10 −5 m/m was demonstrated. The study was performed at a site where the depth to water from land surface ranged from 40.1 to 94.2 m and the aquifer transmissivity was estimated at 41,300 m 2 /d (hydraulic conductivity of 18,800 m/d). The goals of the study were to reduce measurement error as much as practicable and assess the importance of factors affecting water level measurement accuracy. Well verticality was the largest source of error (0.000 to 0.168 m; median of 0.014 m), and geodetic survey of casing elevations was the next most important source of error (0.002 to 0.013 m; median of 0.005 m). Variability due to barometric pressure fluctuations was not an important factor at the site. Hydraulic heads were measured to an accuracy of ±0.0065 m, and the average hydraulic gradient was estimated to be 8.0 × 10 −6 (±0.9 × 10 −6 ) m/m. The improvement in accuracy allowed for two reversals in the groundwater flow direction to be identified, after which the gradient averaged 2.5 × 10 −5 (±0.4 × 10 −5 ) m/m. This study showed it is possible to sufficiently control sources of error to measure hydraulic gradients in the 10 −6 to 10 −5 m/m range.

McDonald, John P.↗

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗

Divergent species-specific impacts of whole ecosystem warming and elevated CO 2 on vegetation water relations in an ombrotrophic peatland

Boreal peatland forests have relatively low species diversity and thus impacts of climate change on one or more dominant species could shift ecosystem function. Despite abundant soil water availability, shallowly rooted vascular plants within peatlands may not be able to meet foliar demand for water under drought or heat events that increase vapor pressure deficits while reducing near surface water availability, although concurrent increases in atmospheric CO 2 could buffer resultant hydraulic stress. We assessed plant water relations of co-occurring shrub (primarily Rhododendron groenlandicum and Chamaedaphne calyculata) and tree (Picea mariana and Larix laricina) species prior to, and in response to whole ecosystem warming (0 to +9°C) and elevated CO 2 using 12.8-m diameter open-top enclosures installed within an ombrotrophic bog. Water relations (water potential [Ψ], turgor loss point, foliar and root hydraulic conductivity) were assessed prior to treatment initiation, then Ψ and peak sap flow (trees only) assessed after 1 or 2 years of treatments. Under the higher temperature treatments, L. laricina Ψ exceeded its turgor loss point, increased its peak sap flow, and was not able to recover Ψ overnight. In contrast, P. mariana operated below its turgor loss point and maintained constant Ψ and sap flow across warming treatments. Similarly, C. calyculata Ψ stress increased with temperature while R. groenlandicum Ψ remained at pretreatment levels. The more anisohydric behavior of L. laricina and C. calyculata may provide greater net C uptake with warming, while the more conservative P. mariana and R. groenlandicum maintained greater hydraulic safety. These latter species also responded to elevated CO 2 by reduced Ψ stress, which may also help limit hydraulic failure during periods of extreme drought or heat in the future. Along with Sphagnum moss, the species-specific responses of peatland vascular communities to drier or hotter conditions will shape boreal peatland composition and function in the future.

59 BASIC BIOLOGICAL SCIENCES↗

Hydraulic Properties Data Package for the E-Area Soils, Cementitious Materials, and Waste Zones–Update

Hydraulic property estimates for the soils, cementitious materials, and waste zones associated with the E-Area low-level radioactive waste disposal units have been provided to support the Performance Assessments (PA) for the E-Area Low-Level Waste Facility (LLWF). Nominal or “best estimate” hydraulic property values for use in the deterministic modeling are provided along with representations of the hydraulic property value uncertainty for use in sensitivity and uncertainty modeling. The hydraulic properties provided for each of the E-Area materials include: porosity ( η ), dry bulk density ( ρ b ), particle density (ρ p ), saturated hydraulic conductivity ( K sat ), characteristic curves (suction head, saturation, and relative permeability), and effective diffusion coefficient ( D e ). A representation of the uncertainty associated with each property, except for the characteristic curves, is provided for each material, except for the E-Area waste zones. These nominal parameter values for each of the E-Area soils, cementitious materials, and waste zones, and, where indicated, parameter uncertainty representations, are based upon the following in order of priority: Site-specific field data, Site-specific laboratory data, Similarity to material with site-specific field or laboratory data, and Literature data. Additionally, a methodology to represent long-term concrete degradation is provided. Finally, because much of the nominal hydraulic property values and uncertainty representations for the E-Area soils, cementitious materials, and waste zones are based on similarity to other materials or literature data, a methodology to prioritize additional work to better define these values and representations is outlined. Prioritization should be based on the importance of the material and/or property to the results of deterministic, sensitivity, and uncertainty modeling. This prioritization should be established through a process of sensitivity modeling. This report supersedes the previous hydraulic properties data package. Revision 1 contains clarifications, updates and corrections to text, figures and tables in a number of places in the report along with some minor editing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

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↗

Understanding the relative importance of vertical and horizontal flow in ice-wedge polygons

Abstract. Ice-wedge polygons are common Arctic landforms. The future of these landforms in a warming climate depends on the bidirectional feedback between the rate of ice-wedge degradation and changes in hydrological characteristics. This work aims to better understand the relative roles of vertical and horizontal water fluxes in the subsurface of polygonal landscapes, providing new insights and data to test and calibrate hydrological models. Field-scale investigations were conducted at an intensively instrumented location on the Barrow Environmental Observatory (BEO) near Utqiaġvik, AK, USA. Using a conservative tracer, we examined controls of microtopography and the frost table on subsurface flow and transport within a low-centered and a high-centered polygon. Bromide tracer was applied at both polygons in July 2015 and transport was monitored through two thaw seasons. Sampler arrays placed in polygon centers, rims, and troughs were used to monitor tracer concentrations. In both polygons, the tracer first infiltrated vertically until encountering the frost table and was then transported horizontally. Horizontal flow occurred in more locations and at higher velocities in the low-centered polygon than in the high-centered polygon. Preferential flow, influenced by frost table topography, was significant between polygon centers and troughs. Estimates of horizontal hydraulic conductivity were within the range of previous estimates of vertical conductivity, highlighting the importance of horizontal flow in these systems. This work forms a basis for understanding complexity of flow in polygonal landscapes.

54 ENVIRONMENTAL SCIENCES↗

Impacts of Topography‐Driven Water Redistribution on Terrestrial Water Storage Change in California Through Ecosystem Responses

Abstract Lateral subsurface flow plays an essential role in sustaining the terrestrial ecosystem, but it is not explicitly represented in most Earth System Models. In this study, we implemented an explicit lateral saturated flow model into the E3SM land model (ELM). The model explicitly describes lateral flow in the saturated zone by representing, for each model grid, an idealized hillslope consisting of five hydrologically connected soil columns. We conducted three model experiments driven by 0.125° atmospheric forcing data during 1980–2015 over California using models of the default ELM, a modified version of ELM to enhance infiltration, and the model with the lateral saturated flow model. The simulated runoff, evapotranspiration, and terrestrial water storage anomaly (TWSA) from the three simulations were evaluated against available observations, and the model explicitly representing lateral flow performs best. The new model produces greater gridcell‐averaged evapotranspiration especially over the mountainous regions with moderate relief and seasonally dry climates. Most importantly, it improves the modeled seasonal variations, interannual variabilities, and the recent decadal decline of TWSA. Many of these improvements can be attributed to the enhanced ecosystem resilience to droughts as demonstrated by transpiration increases caused by lateral flow. Model sensitivity experiments suggest that subsurface runoff is most sensitive to the ratio between horizontal and vertical saturated hydraulic conductivity, followed by hillslope planforms (convergent, divergent, and uniform), number of columns, and lower boundary conditions. Future work should effectively characterize hillslopes in global models and explore the long‐term influences of lateral water movement on modeled biogeochemical cycle.

54 ENVIRONMENTAL SCIENCES↗

Anatomical and hydraulic responses to desiccation in emergent conifer seedlings

Premise The young seedling life stage is critical for reforestation after disturbance and for species migration under climate change, yet little is known regarding their basic hydraulic function or vulnerability to drought. Here, we sought to characterize responses to desiccation including hydraulic vulnerability, xylem anatomical traits, and impacts on other stem tissues that contribute to hydraulic functioning. Methods Larix occidentalis , Pseudotsuga menziesii , and Pinus ponderosa (all ≤6 weeks old) were imaged using x‐ray computed microtomography during desiccation to assess seedling biomechanical responses with concurrently measured hydraulic conductivity ( k s ) and water potential ( Ψ ) to assess vulnerability to xylem embolism formation and other tissue damage. Results In non‐stressed samples for all species, pith and cortical cells appeared circular and well hydrated, but they started to empty and deform with decreasing Ψ which resulted in cell tearing and eventual collapse. Despite the severity of this structural damage, the vascular cambium remained well hydrated even under the most severe drought. There were significant differences among species in vulnerability to xylem embolism formation, with 78% xylem embolism in L. occidentalis by Ψ of −2.1 MPa, but only 47.7% and 62.1% in P. ponderosa and P. menziesii at −4.27 and −6.73 MPa, respectively. Conclusions Larix occidentalis seedlings appeared to be more susceptible to secondary xylem embolism compared to the other two species, but all three maintained hydration of the vascular cambium under severe stress, which could facilitate hydraulic recovery by regrowth of xylem when stress is relieved.

Miller, Megan L.↗

Numerical study of proppant transport and settling processes in fractures

Reservoir stimulation by creating hydraulically conductive fractures is the key step for enabling enhanced geothermal systems (EGS). The effectiveness of stimulation is significantly influenced by the deposition of proppant inside induced fractures. The transportation and settling of proppant in a propagating fracture is controlled by a multitude of operational and physical parameters, including the fracturing fluid rheology, injection rate, proppant concentration, fracture length/aperture evolution, proppant size/density/shape, etc. A numerical tool that robustly and efficiently accounts for all important attributes can facilitate the design and optimization of reservoir stimulation. This study presents the novel computational tool ELK (ELectrical fracKing) developed for the numerical simulation of proppant-fluid mixture circulation in a fractured geothermal reservoir. We enriched the MOOSE-based PorousFlow module with a suite of equations to consider the fluid-proppant mixture with particle-particle/fluid interactions, which include gravitational settling, particle convection, particle hampering, and strong density and viscosity contrasts. The computational tool is validated by comparing the predicted proppant bed evolution against two different laboratory scale experiments of proppant transport in a fixed aperture channel. Further parameter studies were performed, and the modeling results show that the proppant deposition is determined by the mixing characteristics and settling of the particles from the slurry. Concentration-dependent density and viscosity lead to an inhomogeneous distribution of the proppant, particle collision, and enhanced settling at the bottom of the fractures. Preliminary coupling with dynamic fracture propagation shows promising results and will be further developed to simulate hydraulic stimulation at high fidelity.

15 GEOTHERMAL ENERGY↗

Zero-Valent Iron Permeable Reactive Barrier to Remediate Volatile Organic Compounds in Groundwater - 20209

The US DOE-SRS, the US EPA, and the South Carolina Department of Health and Environmental Control determined it was appropriate to perform a non-time critical removal action at the P-Area Groundwater Operable Unit at the SRS to reduce the mass and downgradient transport of trichloroethylene in the P-Area groundwater plume. Contaminated groundwater discharges to a nearby stream, Steel Creek, within the SRS boundaries, resulting in trichloroethylene concentrations above the maximum contaminant level. Impact to surface water is limited in areal extent and supported by recently collected characterization data. The P-Area Groundwater Operable Unit encompasses the groundwater beneath an industrial area within SRS, P Area, where the P-Reactor once operated. The boundaries of the P-Area Groundwater Operable Unit extend northwest to Steel Creek, northeast toward PAR Pond, and southeast to Meyers Branch. Groundwater in the Upper Three Runs Aquifer of the P-Area Groundwater Operable Unit has been impacted by reactor and facility operations between 1954 and 1991, including tritium and volatile organic compounds. The P-Area surface units contributing to groundwater contamination were remediated as part of the P-Area Operable Unit in 2011. The P-Reactor closure is one of the first of its kind in the DOE Complex and is one of only a few full-sized production reactors in the US to undergo completion of final closure activities. The nature and extent of groundwater contamination was determined using a variety of investigative approaches such as groundwater monitoring wells, direct-push technology, and surface water samples. Groundwater contamination associated with trichloroethylene is primarily exhibited in a narrow plume that extends from the source area at P-Reactor and west to Steel Creek. Maximum contaminant level exceedances in groundwater occur over an area of ∼6.9 hectares for trichloroethylene with concentrations as high as 7.7 milligrams per liter. To the west of the P-Area facility area, the trichloroethylene groundwater plume is controlled by a buried geologic feature, assumed to be an old stream bed, that further narrows the groundwater plume in what has been designated as the 'neck area.' This narrowing of the groundwater plume provides an ideal location for a treatment barrier. The non-time critical removal action alternative chosen is to install a zero-valent iron permeable reactive barrier within the neck area of the trichloroethylene groundwater plume, perpendicular to groundwater flow direction. This technology will provide a treatment barrier that will reduce trichloroethylene groundwater concentrations by 90% and has an anticipated useful life of at least 25 years. A pre-design investigation was performed in the neck area to confirm site lithology, hydrogeology, geochemistry, and extent of trichloroethylene contamination prior to a final design. A treatability study, conducted as part of the pre-design investigation, indicated that the subsurface and groundwater in the PArea Groundwater Operable Unit is compatible with the zero-valent iron and will not lead to excessive buildup from mineralization/precipitation or biofouling. Probabilistic modeling was conducted using field and laboratory data to determine the expected performance of the zero-valent iron permeable reactive barrier. The model simulations indicated that a 3.81-centimeter thick barrier would provide greater than 90% reduction of trichloroethylene groundwater concentration. The final design of the zero-valent iron permeable reactive barrier is a barrier that will extend 80.5 linear meters in a 'zigzag' orientation to best transect the trichloroethylene plume and account for varying groundwater flow. The barrier will be installed from 13.7 meters below ground surface to 41.1 meters below ground surface for 65.8 linear meters and from 13.7 meters below ground surface to 36.6 meters below ground surface for 14.6 linear meters, the base of which is 'keyed' into a low permeability zone. The barrier is designed to a thickness of 10.2 centimeters, which was determined to reduce the trichloroethylene groundwater concentrations by greater than 90% with a safety factor of 2.67. A total of approximately 689 metric tons of zero-valent iron will be injected through 22 injection wells spaced 3.66 meters apart, using guar to suspend the zero-valent iron. Zero-valent iron permeable reactive barrier construction will be monitored through 23 installed resistivity receivers offset 7.32 meters from the zerovalent iron permeable reactive barrier. The zero-valent iron will be energized with a low-voltage 100 Hertz signal during injection and will be monitored using the resistivity receivers to ensure complete coalescence of the zero-valent iron permeable reactive barrier. The zero-valent iron was sized to have a hydraulic conductivity greater than the natural subsurface, thus promoting groundwater flow through the barrier. As contaminated groundwater contacts the zero-valent iron, volatile organic compounds, including trichloroethylene, are immediately degraded to harmless compounds such as ethylene. The performance of the zero-valent iron permeable reactive barrier will be monitored using three upgradient monitoring well clusters, six downgradient monitoring well clusters, and four in-wall monitoring wells. The in-wall monitoring wells will indicate immediate reduction of trichloroethylene mass in the groundwater and allow analyses of the zero-valent iron permeable reactive barrier health. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Studies on the Encapsulation of Brine Generated from a Process Using Selective Electrodialysis Reversal

This topical report summarizes the laboratory testing, geochemical modeling, and mineralogical testing conducted to examine the efficacy of wastewater encapsulation for the brine generated from the Saltworks Flex electrodialysis reversal technology. Within the laboratory testing, the short- and long-term physical properties of seven separate mix designs were examined along with the leaching of each mix design. Additionally, geochemical thermodynamic modeling was used along with mineralogical analyses to help gain a fundamental understanding of the chemistry involved in the encapsulation process. Encapsulation is a process that results in the immobilization of constituents of interest by creating an engineered material that contains the waste within its matrix, both physically and chemically, in such a manner that the receiving environment is protected from adverse effects. Within the context of this study, encapsulation mix designs involved mixing the Saltworks brine with fly ash and various additives, principally lime (calcium oxide) and cement, to form a paste that eventually sets and hardens into a monolith upon placement in a landfill. Moisture conditioning fly ash with wastewater sources for on-site landfill disposal is a process that is in use at some sites and has been presented as a viable option to dispose of liquid wastes. However, because chemical reactions that can lead to sequestration of constituents of interest are not optimized, leaching can occur. Therefore, in this study, mixes were created to mimic moisture conditioning, as well as engineered pastes of varying mix designs. The study found that engineered encapsulation mixes provide significant resistance to leaching compared to moisture conditioning alone. The best performing encapsulation mix used lime as the additive, which proved superior to the use of cement. After 90 days of curing, this encapsulated mix was 5.5 times stronger than the moisture conditioned mix. Additionally, this mix exhibited an approximate 2.5 order of magnitude drop in hydraulic conductivity compared to the moisture conditioned mix. When comparing leaching, nearly all the chloride was found to evolve from the moisture conditioned mix, whereas only 6.7% of the chlorine contained in the best performing encapsulated mix leached out. Mineralogical analysis via X-ray diffraction (XRD) showed a significant concentration of the chloride sequestering mineral, Friedel’s salt, in the best performing encapsulation mix, whereas none was observed in the moisture conditioned mix, pointing to it being the primary chemical mechanism for chloride retention, along with the physical impermeability of the material. With respect to the modeling study, the model was successful in predicting the presence of most mineral phases observed in the XRD analysis.

01 COAL, LIGNITE, AND PEAT↗

Artificial intelligence models, photos, and data associated with the manuscript “Quantifying Streambed Grain Size, Uncertainty, and Hydrobiogeochemical Parameters Using Machine Learning Model YOLO” (v2)

This data package is associated with the manuscript “Quantifying Streambed Grain Size, Uncertainty, and Hydrobiogeochemical Parameters Using Machine Learning Model YOLO” published in Water Resources Research (Chen et al., 2024). This data package includes the training, validation, testing, and prediction data used by the artificial intelligence (AI) model for automated grain size and hydro-biogeochemistry quantification using streambed photos. The grain size data are extracted for each photo using You Look Only Once (YOLO), a pre-trained object detection model. This data package was originally published in October 2023. It was updated August 2025 (v2; new and modified files). File and folder names were not revised to indicate changes. See the change history section in the readme for more details. Please see flmd.csv for a list of all files contained in this data package and descriptions for each. Please see dd.csv for a data dictionary that defines the column headers of .csv files in the data package. This dataset is comprised of one data folder containing (1) file-level metadata; (2) data dictionary; (3) readme; and (4) six subfolders. Subfolders 1 to 4 include the training, validation, testing, and prediction data. Subfolder 5_Summary includes the summary results of different combinations of training, validation, testing, and prediction data. Subfolder 6_SupplementalData includes additional data downloaded from public sources (Kaufman et al., 2023a; Kaufman et al., 2023b; Garefalakis et al., 2023; Mair et al., 2024; https://github.com/river-corridors-sfa/Geospatial_variables). In total, the data package includes 110 folders and 44,283 files. These files include 9,047 .jpg photos, 1 .png photo, 3 .tif photos; 26,639 photo labels and individual grain sizes and probability from AI (.txt); 8,447 grain size distribution data (.dat); and 126 CSV files for results summary, and 14 required metadata files (.xlsx). The summary CSV files contain 68 columns and approximately 2,200 rows that represent photo names, site locations, recording time, GPS coordinates, grains sizes (D10, D50, D60, and D84), number of grains, and additional hydro-biogeochemical data such as water depth, flow velocity, Manning’s coefficient, friction factor, hydraulic conductivity, permeability, streambed interstitial velocity magnitude, mass transfer rate, and nitrate uptake velocity. The photos were obtained from 75 sites in the Yakima River Basin and the Columbia River shorelines, and other associated data from samples and sensors obtained when the photos were taken are publicly available (Fulton et al. 2022; Grieger et al. 2023). All files are .csv, .txt, .dat, .jpg, or .pdf. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected some of these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗