Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “WATER”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Pump House in East River Watershed, Colorado 2019-2024

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Pump House at Mount Crested Butte in the East River Watershed. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format ER-X-Y, where ER refers to East River, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, ER-PHS, ER-LMC, ER-LMF, and ER-SMN are associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and ER-RBTn (upslope n=1) are sampling transects during the 2019 Rootball Campaign. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2024

Produced water sharing: Improved economics and reduced community impact – A Pennsylvania case study

Here, hydraulic fracturing for oil and gas extraction from unconventional reservoirs is water intensive. Between water sourced for fracturing purposes and water present in rock formations, operators often produce a greater volume of water than oil or gas. Historically, this surplus of produced water has mainly been disposed of via deep injection wells. Rising disposal costs, seasonally limited water availability, and concerns over induced seismicity have incentivized produced water recycling practices, where an operator uses produced water for hydraulic fracturing operations. The logistical challenges associated with produced water recycling have also encouraged operators to adopt ad-hoc water exchange practices, in which competing operators will exchange produced water for mutual cost savings. In this paper, we investigate the potential benefit from systematic produced water exchange among operators in Northeastern Pennsylvania. We leverage PARETO, a free and open-source modeling framework for produced water management optimization, to quantify the benefits of water exchange practices. In an example drawn from FracFocus data, we find that the adoption of systematic water sharing could improve produced water recycling rates from 49.2% to 99%, decreasing operating and trucking costs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

GCAM–GLORY v1.0: representing global reservoir water storage in a multi-sector human–Earth system model

Abstract. Reservoirs play a significant role in modifying the spatiotemporal availability of surface water to meet multi-sector human demands, despite representing a relatively small fraction of the global water budget. Yet the integrated modeling frameworks that explore the interactions among climate, land, energy, water, and socioeconomic systems at a global scale often contain limited representations of water storage dynamics that incorporate feedbacks from other systems. In this study, we implement a representation of water storage in the Global Change Analysis Model (GCAM) to enable the exploration of the future role (e.g., expansion) of reservoir water storage globally in meeting demands for, and evolving in response to interactions with, the climate, land, and energy systems. GCAM represents 235 global water basins, operates at 5-year time steps, and uses supply curves to capture economic competition among renewable water (now including reservoirs), non-renewable groundwater, and desalination. Our approach consists of developing the GLObal Reservoir Yield (GLORY) model, which uses a linear programming (LP)-based optimization algorithm and dynamically linking GLORY with GCAM. The new coupled GCAM–GLORY approach improves the representation of reservoir water storage in GCAM in several ways. First, the GLORY model identifies the cost of supplying increasing levels of water supply from reservoir storage by considering regional physical and economic factors, such as evolving monthly reservoir inflows and demands, and the leveled cost of constructing additional reservoir storage capacity. Second, by passing those costs to GCAM, GLORY enables the exploration of future regional reservoir expansion pathways and their response to climate and socioeconomic drivers. To guide the model toward reasonable reservoir expansion pathways, GLORY applies a diverse array of feasibility constraints related to protected land, population, water sources, and cropland. Finally, the GLORY–GCAM feedback loop allows evolving water demands from GCAM to inform GLORY, resulting in an updated supply curve at each time step, thus enabling GCAM to establish a more meaningful economic value of water. This study improves our understanding of the sensitivity of reservoir water supply to multiple physical and economic dimensions, such as sub-annual variations in climate conditions and human water demands, especially for basins experiencing socioeconomic droughts.

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 OSTI2024

Using stable isotopes to inform water resource management in forested and agricultural ecosystems

Present and future climatic trends are expected to markedly alter water fluxes and stores in the hydrologic cycle. In addition, water demand continues to grow due to increased human use and a growing population. Sustainably managing water resources requires a thorough understanding of water storage and flow in natural, agricultural, and urban ecosystems. Measurements of stable isotopes of water (hydrogen and oxygen) in the water cycle (atmosphere, soils, plants, surface water, and groundwater) can provide information on the transport pathways, sourcing, dynamics, ages, and storage pools of water that is difficult to obtain with other techniques. However, the potential of these techniques for practical questions has not been fully exploited yet. Here, we outline the benefits and limitations of potential applications of stable isotope methods useful to water managers, farmers, and other stakeholders. We also describe several case studies demonstrating how stable isotopes of water can support water management decision-making. Finally, we propose a workflow that guides users through a sequence of decisions required to apply stable isotope methods to examples of water management issues. We call for ongoing dialogue and a stronger connection between water management stakeholders and water stable isotope practitioners to identify the most pressing issues and develop best-practice guidelines to apply these techniques.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Hidden features in the OH-stretching spectra of amino acid decorated air–water interfaces

Chemical reactivity at the air–water interface is governed by the interfacial solvation of reactive species. For instance, during aqueous amino acid-based CO 2 absorption, water reorganizes around the reactive sites and couples dynamically with reaction pathways, facilitating the reaction. In this context, surface-sensitive vibrational sum-frequency generation (vSFG) spectroscopy can probe the OH stretch vibrations of interfacial water and determine the solvation structures around reactants and products, thereby furthering our understanding of the role of interfacial solvation. However, vSFG spectra of the air–water interface in the presence of charged species can be remarkably complex; key species-bound local water structures with distinct orientations may be hidden beneath prominent vSFG peaks arising from water–water hydrogen bonds and remain difficult to resolve. Here, we measure and compute vSFG spectra of the water OH stretch at air–water interfaces decorated with amino acids in their zwitterionic and anionic forms, as well as equimolar mixtures of these forms with bicarbonate. The latter represents post-CO 2 -absorption conditions. We find that computing depth- and frequency-dependent spectral densities—decomposed into contributions from water molecules hydrogen-bonded exclusively to other water molecules, exclusively to amines, exclusively to carboxylates, or shared between these polar/charged groups—is indispensable for accurate interpretation of the vSFG spectra. Key findings include orientational flip-flop in water sub-layers, strong carboxylate-water H-bonding, and water orientational ordering extending into the bulk aqueous phase induced by anionic amino acids. Here, this study provides a computational spectroscopic platform for improved understanding of interfacial solvation relevant to interfacial reactivity.

Air-water interface↗
DOE OSTI2026

The Kinetic Consequences of Water on Catalytic Methane Pyrolysis

Hydrogen production from biomass and natural gas has emerged as a prominent research area in response to the growing demand for energy from alternative sources that minimize CO 2 emissions. In this study, we investigate the impact of water, which is present in and generated from biomass-derived streams, on carbon nanotube (CNT) growth and hydrogen production during methane decomposition using Ni–Mo/MgO as a catalyst. We reveal here that the role of water on CNT growth is highly complex; its effect depends on the stage of growth at which the water is incorporated. When water is introduced at the beginning of methane decomposition ( t = 0 h), methane conversion rates are negatively impacted. We hypothesize that water inhibits the significant phase changes the Ni–Mo/MgO catalyst undergoes during catalyst carburization. In contrast, the incorporation of a small percentage of water after a stabilization period ( t = 3 h) results in methane conversion rate enhancements that scale with the introduced water partial pressure as water selectively reacts with amorphous carbon deposits that lead to catalyst deactivation, thus prolonging the lifetime of some of the most active sites. Moreover, water incorporation after stabilization significantly reduces the apparent activation energy. Density Functional Theory (DFT) calculations reveal that water preferentially interacts with carbon fragments on the catalyst surface to remove carbon deposits with a barrier lower than that required for methane activation, further supporting its role in cleaning active sites on the catalyst surface. Characterization of the resulting carbon nanotubes reveals the formation of more graphitic materials produced in the presence of water, highlighting the impact of water on nanotube properties. These results provide clarity toward the many ways in which water, or cofeeding of biomass-derived materials, may impact catalytic methane pyrolysis rates.

carbon nanotubes↗
DOE OSTI2024

A Framework Using Applied Process Analysis Methods to Assess Water Security in the Vu Gia–Thu Bon River Basin, Vietnam

The Vu Gia–Thu Bon (VG–TB) river basin is facing numerous challenges to water security, particularly in light of the increasing impacts of climate change. These challenges, including salinity intrusion, shifts in rainfall patterns, and reduced water supply in downstream areas, are of great concern. This study comprehensively assessed the current state of water security in the basin using robust statistical analysis methods such as the Process Analysis Method (PAM), SMART principle, and Analytic Hierarchy Process (AHP). This resulted in the development of a comprehensive assessment framework for water security in the VG–TB river basin. This framework identified five key dimensions, with basin development activities (0.32), the ability to meet water needs (0.24), and natural disaster resilience (0.19) being the most crucial and water resource potential being the least crucial (0.11) according to the AHP methodology. The latter also highlighted 15 indicators, four of which are particularly influential, including waste resources (0.54), flood (0.53), water storage capacity (0.45), and basin governance (0.42). Furthermore, 28 variables with high weight factors were identified. This framework aligns with the UN-Water water security definition and addresses the global water sustainability criteria outlined in Sustainable Development Goal 6 (SDG6). It enables the computation of a comprehensive Water Security Index (WSI) for specific regions, providing a strong foundation for decision-making and policy formulation. It aims to enhance water security in the context of climate change and support sustainable basin development, thereby guiding future research and policy decisions in water resource management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Beyond Leakage: Non-Revenue Water Loss and Economic Sustainability

Water loss in urban supply systems poses significant challenges for water utility companies worldwide, affecting both sustainable access to clean water and the financial viability of utility operations. This study analyzes the evolution of water losses in high-level supply systems from 2017 to 2021 in Portugal, focusing on its implications for the profitability of water utility companies across NUTs II regions. Drawing on data from various sources, including the National Information System for Water Resources, PORDATA, ERSAR, and ORBIS, this analysis identifies trends, patterns, and potential factors influencing water loss dynamics. Key components of the analysis include calculating average annual losses, examining unbilled water percentages, and conducting regression analysis to quantify the impact of water loss on profit margins. The findings contribute to the literature on water loss management and financial performance in the water utility sector, offering insights for policymakers, water utility managers, and stakeholders to enhance financial sustainability and reduce water losses.

Santos, Eleonora (ORCID:0000000346930804)↗
DOE OSTI2026

Understanding Operando Water Management in Hydroxide‐Exchange‐Membrane Fuel Cells

The water balance in hydroxide-exchange-membrane fuel cells (HEMFCs) is a key challenge for improved performance and durability, intimately linked with the various interfaces and coupled phenomena. For every 4 electrons produced, 4 water molecules are generated in the anode and 2 consumed in the cathode, while electroosmosis transports water across the HEM from the cathode to the anode. Consequently, a concentration gradient drives water back, from anode to cathode. Ineffective water management could lead to cathode dry-out, limiting reaction rate and causing ionomer degradation, or to anode flooding. To address these concerns, it is critical to measure the water transport operando . Herein, a home-built water-flux station is used to measure total water flux during cell operation with different inlet relative humidities and back pressures. Increasing the HEM thickness fourfold decreases the water flux at high current density, and utilizing microporous layers on both the anode and cathode decreases the water flux from the anode to the cathode. However, the most significant variable in changing the water flux was found by increasing the anode back pressure. Furthermore, humidity cycling significantly changed electrochemical performance without affecting the overall water fluxes. These findings can be translated to other devices utilizing an HEM.

AEMFC↗
DOE OSTI2024

Evaluating the Chemical Reactivity of DFT-Simulated Liquid Water with Hydrated Electrons via the Dual Descriptor

Modeling the various properties of liquid water, particularly its reactivity, has been a longstanding challenge for simulation methods. Recently, ab initio simulations based on density functional theory (DFT) have come to the fore as tenable methods for calculating the properties and reactivity of water, with varying degrees of success for different exchange-correlation functionals. In particular, hybrid-GGA and meta-GGA functionals have been shown to reproduce many of the structural, dynamical, and energetic properties of water to a high degree of accuracy relative to their computational cost. Here, we show that the dual descriptor (DD) measure of nucleophilicity and electrophilicity, which is sometimes used to elucidate organic chemistry reaction mechanisms, can also be used to characterize the reactivity of DFTsimulated liquid water. The DD is especially apt for understanding the reactivity of excess electrons with water as its calculation explicitly involves adding and removing an excess electron from a reference system. We use the DD to explore the reactivity of water simulated using three different DFT functionals: the LDA functional (LDA), a hybrid-GGA functional (PBE0), and a hybrid meta- GGA functional (SCAN0). Using the DD, we show that the SCAN0 functional with the standard 25% Hartree–Fock exchange produces simulated liquid water with many regions that are far more reactive than either PBE0 or LDA. To understand the implications of these highly reactive regions, we then add a strong nucleophile in the form of an excess electron and find that although PBE0 and LDA predict stable hydrated electrons, the excess electron reacts nearly instantaneously with SCAN0 water via proton abstraction to form a hydrogen atom and hydroxide ion. We show that the DD provides the ability to not only predict whether or not liquid water will react with a hydrated electron but also which particular waters will be involved solely from analyzing pure water configurations generated with each functional. We rationalize this result in terms of the known trap-seeking behavior of injected hydrated electrons, which are able to find the most electronegative region in bulk water. These results highlight the utility of the dual descriptor as a fast and interpretable method for investigating condensed-phase reactivity with excess electrons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Urban morphology and urban water demand: a case study in the land constrained Los Angeles region using urban growth modeling

The interactions between population growth, urban morphology, and water demand have important implications for water resources and supply in urban regions. Water use for irrigation comprises a significant fraction of urban water demand, and is potentially influenced by long-term changes in urban morphology. To investigate this, we used spatially explicit projections of urban land development intensity (fraction impervious area) generated from a 30 m resolution urban growth model for the Los Angeles (LA) region. Recent historical data on water use and high resolution landcover were used to establish relationships between green area, urban development intensity, and outdoor water demand. These relationships were then used to project outdoor and total water demand in 2100 using the urban growth model outputs. We considered two different population scenarios informed by the shared socioeconomic pathway (SSP) projections for the region (SSP3 and SSP5), and three scenarios of urban development intensification. Our analysis is resolved for over 80 water providers in the region, from the urban core to suburban fringe, and highlights diverse demand responses influenced by initial urban form and water demand attributes. Assumptions about outdoor water use factors based on recent water supply data were found to be nearly as influential on future outdoor demand as the urban growth scenario settings. Compared to previous studies, our work is unique in coherently linking high resolution SSP population scenarios, urban land cover evolution, and urban water demand projections, demonstrating the approach for the LA region—the largest population center in the western United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Side-Chain and Ring-Size Effects on Permeability in Artificial Water Channels

Artificial water channels (AWCs) have emerged as a promising framework for stable water permeation, with water transport rates comparable to aquaporins (3.4–40.3 × 10 8 H 2 O/channel/s). In this study, we probe the influence of ring-size and side-chain length on the water permeability observed within a class of AWCs termed ligand-appended pillar[n]arenes (LAPs) that have an adjustable ring-size (m) and side-chain length (n). Through all-atom molecular dynamics simulations, we calculate the permeability of these channels using the collective diffusion model and find their permeabilities. We characterize the mechanistic influence of pillar[n]arene ring-size and side-chain length on the channel water permeability by analyzing the characteristics of the internal permeating water-wire and the surrounding channel structure. We observe that water permeability decreases as a function of increasing ring-size due to increases in hydrophilic contacts between the permeating water-wire and the oxygen groups on the channel wall. Further, we observe an increase in water permeability as a function of side-chain length due to increased partitioning of the channel terminal groups into the hydrophilic blocks of the surrounding bilayer. For the LAP6 channel, with increase in side-chain length, the distance between terminal groups increases and leads to an increase in pore size, thereby enhancing water permeability. In the case of LAP5, as side-chain length increases, the channel displays a compensatory effect between tilt and bend angle due to the flexible side-chains. Such flexibility leads to higher terminal group partitioning in the hydrophilic blocks of the bilayer and extends the permeating water-wire. Furthermore, this increase in water-wire length and hydrophilic block access overcomes the nonmonotonic pore size trend in pillar[5]arene channels.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular Structure and Thermodynamics of CO 2 and Water Adsorption on Mica

The adsorption of CO 2 and water on clay surfaces plays a key role in applications, such as gas storage in saline aquifers and depleted hydrocarbon reservoirs, but is not yet fully understood. Here, we study the adsorption of CO 2 and water vapor using Grand Canonical Monte Carlo and molecular dynamics simulations. At a bulk pressure of 100 bar, pure CO 2 adsorbs strongly on mica and forms extensive layers next to it. CO 2 adsorption is lowered substantially if introducing water vapor above mica and is largely eliminated when the relative humidity (RH) approaches about 60%. When pure water vapor is introduced above a mica surface, a subnanometer thick liquid water film develops on it to form apparent liquid–solid and liquid–vapor interfaces simultaneously. Using the identification of truly interfacial molecules (ITIM) analysis, we delineate how individual water layers develop in this film as RH increases. We highlight that the water film is spatially heterogeneous and the true liquid–vapor interface emerges only at an RH of 60–80%. Introducing 100 bar of CO 2 into the water vapor above the mica surface modulates water adsorption nonlinearly: at RH = 0.01%, the water adsorption is reduced by ∼30%; as RH increases, the reduction is weakened, and eventually, enhancement of water adsorption by about 7% occurs at RH = 90%. These variations are attributed to the interplay of film thinning by high-pressure CO 2 , competition of mica surface sites by CO 2 molecules, and energetic and entropic stabilization of interfacial water by CO 2 molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The possibility of hydrogen-water demixing in Uranus, Neptune, K2-18 b and TOI-270 d

Context. The internal structures of Uranus and Neptune remain unknown. In addition, sub-Neptunes are now thought to be the most common type of exoplanets. Improving our understanding of the physical processes that govern the interiors of such planets is therefore essential. Phase separation between planetary constituents may occur, in particular, hydrogen-water immiscibility in cold, water-rich intermediate-mass planets. Aims. We assess whether hydrogen-water demixing could occur in Uranus, Neptune, K2-18 b and TOI-270 d, and investigate its effect on the planetary evolution and inferred internal structure. Methods. We couple planetary evolution models with recent ab initio calculations of the hydrogen-water phase diagram, allowing for temperature shifts to account for uncertainties in miscibility gaps. Results. We find that demixing may occur and could lead to a complete depletion of water in the outermost regions of Uranus and Neptune. Temperature offsets of up to 1100 K lead to a depleted region comprising as much as 16% of the planet’s mass, and an increase in planetary radius by nearly 20%. For K2-18 b, our models suggest that hydrogen-water demixing is ongoing and may explain the absence of water features in its JWST spectrum. A temperature offset of 500 K is required to get a complete depletion of water in the atmosphere of K2-18b. TOI-270d may also have experienced hydrogen-water demixing. When applying a similar temperature offset on the phase diagram as for K2-18 b, we find a partial depletion of water in the atmosphere of TOI-270 d, consistent with JWST’s detection of water. Conclusions. Hydrogen-water immiscibility may play a key role in shaping the structure and evolution of both Solar System giant planets like Uranus and Neptune, and cold/temperate exoplanets such as K2-18 b and TOI-270 d. Accounting for such internal processes is crucial to accurately interpret atmospheric observations from current (e.g., JWST) and upcoming (e.g., ARIEL) missions.

planets and satellites: composition↗
DOE OSTI2024

Gaining Real-Time Water Leak Detection

Devens Reserve Forces Training Area is a United States Army Reserve (USAR) Installation that struggles with severe water leaks, often causing significant damage to the facility and requiring major renovation. Traditional water use is highly dependent on occupancy, so it can be difficult to benchmark a facility’s water use. It can be exceptionally difficult when occupancy is transient and/or varies. Pacific Northwest National Laboratory (PNNL) collaborated with Devens to implement real-time monitoring of their water consumption by utilizing the smart meter data from their existing 23 water meters. PNNL created a simple algorithm to calculate hourly water consumption and trigger an alert to be instantly emailed to Devens’ personnel when there appears to be a water leak in any building with a smart water meter. Here, this approach is expected to save hundreds of thousands of dollars in unnecessary water consumption costs and damages from leaks and was implemented with little-to-no costs or service disruptions. Next steps for this project include slow leak detection through nighttime monitoring and to extrapolate this water leak approach to the remainder 360 water meters on USAR’s Enterprise Building Control System so USAR sites across the country can be instantly notified of potential water leaks.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗