Search NASA⌕ Search

SEARCH · Search NASA

Results for “water availability”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

METSiCloPs in Cooling Tower Water, Final Presentation: Water Treatment at FSI [Slides]

Motivation: According to current projections, FSI could require anywhere from 83-288 MGY of water. Reducing excessive dependence on external water supply by making operations as water-lean as possible in the face of water availability fluctuations seasonally. The amount of Silica occurring naturally in New Mexico is several multiples higher than most of the country, with average values of 90 ppm, but reaching as high as 113 ppm. Since its saturation concentration at ambient temperatures is about 150 ppm, this corresponds to cycles of concentration as low as 1-2 for untreated water in cooling towers (i.e., the water is circulated through the towers only once or twice before they are dumped into the outfall), thereby making operations highly water intensive. There is, therefore, there is a strong need and internal motivation to remove Silica from the water that is circulated in the cooling towers.

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

SPRUCE Precipitation Chemistry and Bulk Atmospheric Deposition Beginning in 2013

This data set reports the chemistry of precipitation water that was collected at the SPRUCE site located in the S1 bog and calculated values of bulk atmospheric deposition of solutes in precipitation. Sample collection and analyses started during December of 2013 and will continue for the duration of the experiment. Precipitation samples were collected year round on an event-basis from up to three individual collectors. Samples were either composited or individually analyzed, depending on the total volume of water available relative to the amount needed to complete all laboratory analyses. Funnel/bottle collectors were used to collect rainfall and buckets were used to collect snowfall and mixed precipitation. Samples were analyzed for pH, specific conductivity, nutrient concentrations (ammonium, nitrate, soluble reactive phosphorus, total nitrogen), cation concentrations (calcium, magnesium), anion concentrations (chloride, sulfate), total organic carbon concentrations, and natural-abundance stable isotopes of water (δ18O-H2O and δD-H2O). Bulk deposition was calculated by multiplying rainfall amount by solute concentration for ammonium, nitrate, soluble reactive phosphorus (SRP), total nitrogen, calcium, magnesium, chloride, sulfate, and total organic carbon. All measurements were made at the 8.1-ha S1 bog forest site in northern Minnesota, 40 km north of Grand Rapids, in the USDA Forest Service Marcell Experimental Forest (MEF).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Optimization of leaf morphology in relation to leaf water status: A theory

The leaf economic traits such as leaf area, maximum carbon assimilation rate, and venation are all correlated and related to water availability. Furthermore, leaves are often broad and large in humid areas and narrower in arid/semiarid and hot and cold areas. We use optimization theory to explain these patterns. We have created a constrained optimization leaf model linking leaf shape to vein structure that is integrated into coupled transpiration and carbon assimilation processes. The model maximizes net leaf carbon gain (NPPleaf) over the loss of xylem water potential. Modeled relations between leaf traits are consistent with empirically observed patterns. As the results of the leaf shape-venation relation, our model further predicts that a broadleaf has overall higher NPPleaf compared to a narrowleaf. In addition, a broadleaf has a lower stomatal resistance compared to a narrowleaf under the same level of constraint. With the same leaf area, a broadleaf will have, on average, larger conduits and lower total leaf xylem resistance and thus be more efficient in water transportation but less resistant to cavitation. By linking venation structure to leaf shape and using water potential as the constraint, our model provides a physical explanation for the general pattern of the covariance of leaf traits through the safety-efficiency trade-off of leaf hydraulic design.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Climate-Water Impacts on Interconnection-Scale Electricity System Planning

A growing literature emphasizes the importance of integrating climate change impacts into electricity system planning. Rising average temperatures can increase and shift electricity demand while reducing generator and transmission efficiency. Changes to water availability and quality can reduce the output of thermally cooled generators and hydropower. Electric power grids across the US and globally are undergoing transformational changes that present new opportunities and challenges to reliability assurance. However, electric utilities and system operators have limited internal capabilities to incorporate these effects into planning practices. This work addresses gaps in utility and system planner practices by integrating climate-water-electricity expertise from universities and U.S. Department of Energy National Laboratories with electricity system planners and stakeholders in the Western Electricity Coordinating Council (WECC). Using a highly collaborative approach, global climate model data, high-resolution hydrology models, and long-term electric sector capacity expansion tools are employed to analyze a range of climate outcomes for future electricity scenarios aligned with recent WECC planning studies. Doing so allows WECC to expand its climate-agnostic planning assessments to consider how future temperature and precipitation patterns could influence generation and transmission planning. We explore how changes to climate-water conditions can affect power plant investment and operation, system economics, and environmental impacts, providing an expanded perspective on interconnection-wide decision making under climate uncertainty.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Burning questions for a warming and changing world: 15 unknowns in plant abiotic stress

Abstract We present unresolved questions in plant abiotic stress biology as posed by 15 research groups with expertise spanning eco-physiology to cell and molecular biology. Common themes of these questions include the need to better understand how plants detect water availability, temperature, salinity, and rising carbon dioxide (CO2) levels; how environmental signals interface with endogenous signaling and development (e.g. circadian clock and flowering time); and how this integrated signaling controls downstream responses (e.g. stomatal regulation, proline metabolism, and growth versus defense balance). The plasma membrane comes up frequently as a site of key signaling and transport events (e.g. mechanosensing and lipid-derived signaling, aquaporins). Adaptation to water extremes and rising CO2 affects hydraulic architecture and transpiration, as well as root and shoot growth and morphology, in ways not fully understood. Environmental adaptation involves tradeoffs that limit ecological distribution and crop resilience in the face of changing and increasingly unpredictable environments. Exploration of plant diversity within and among species can help us know which of these tradeoffs represent fundamental limits and which ones can be circumvented by bringing new trait combinations together. Better defining what constitutes beneficial stress resistance in different contexts and making connections between genes and phenotypes, and between laboratory and field observations, are overarching challenges.

Biochemistry & Molecular Biology↗
DOE OSTI2020

Crop Modeling Application to Improve Irrigation Efficiency in Year-Round Vegetable Production in the Texas Winter Garden Region

Given a rising demand for quality assurance, rather than solely yield, supplemental irrigation plays an important role to ensure the viability and profitability of vegetable crops from unpredictable changes in weather. However, under drought conditions, agricultural irrigation is often given low priority for water allocation. This reduced water availability for agriculture calls for techniques with greater irrigation efficiency, that do not compromise crop quality and yield, and that provide economic benefit for producers. This study developed vegetable growing models for eight different vegetable crops (bush bean, green bean, cabbage, peppermint, spearmint, yellow straight neck squash, zucchini, and bell pepper) based on data from several years of field research. The ALMANAC model accurately simulated yields and water use efficiency (WUE) of all eight vegetables. The developed vegetable models were used to evaluate the effects of various irrigation regimes on vegetable growth and production in several locations in the Winter Garden Region of Texas, under variable weather conditions. Based on our simulation results from 960 scenarios, optimal irrigation amounts that produce high yield as well as reasonable economic profit to producers were determined for each vegetable crop. Overall, yields for all vegetables increased as irrigation amounts increased. However, irrigation amounts did not have a sustainable impact on vegetable yield at high irrigation treatments, and the WUEs of most vegetables were not significantly different among various irrigation regimes. When vegetable yields were compared with water cost, the rate decreased as irrigation amounts increased. Thus, producers will not receive economic benefits when vegetable irrigation water demand is too high.

Kim, Sumin↗
DOE OSTI2025

Data and scripts associated with “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” (v3)

This data package is associated with the publication “Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics within Water-Tagging Enabled Hydrologic Models” submitted to Journal of Advances in Modeling Earth Systems (Butler et al. 2025). This study developed the Sequential Precipitation Input Tagging (SPIT) framework to tag input precipitation and estimate water transit times and hydrologic tracers. SPIT tags all precipitation events at regular intervals over an extended period (monthly tags over seven years) in a hydrologic model from 2016-2022. SPIT is applied at six National Ecological Observatory Network (NEON) sites across the continental United States to calculate transit time distributions (TTD) and derive from these mean transit times (MTT), fractions of young water (Fyw), and hydrologic tracer concentrations in stream water (δ18O) within a water-tagging enabled version of the Weather Research and Forecast (WT-WRF-Hydro) model with national water model (NWM) configurations. We go on to validate WT-WRF-Hydro estimates against Butler et al. (2023), who analyzed the same NEON sites using stable water isotope data to estimate water transit times. This new tracking method provides a detailed picture of water movement and helps improve predictions about water availability in the future. This data package was originally published in January 2025. It was updated May 2025 (v2; new and modified files) and October 2025 (v3; 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. This data package contains the data and scripts used to develop the SPIT framework WT-WRF-Hydro (Water Tagging Weather Research and Forecasting Hydrologic) model and is associated with the following GitHub repository: https://github.com/zbutler33/SPIT-Framework. This data package contains five parent folders: (1) “Manipulated_outputs”, (2) “Metadata”, (3) “Observed”, (4) “Outputs”, and (5) “Scripts”. Each of these parent folders contains additional subfolders and files. Please see the FLMD (“v*_Butler_2024_WT_WRF_Hydro_flmd.csv”) for a list of all the files contained in this data package and descriptions for each. See the data dictionary (“v*_Butler_2024_WT_WRF_Hydro_dd.csv”) for definitions and units of all of the tabular (files ending in “.csv” and ".tsv") column headers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Large Divergence in Tropical Hydrological Projections Caused by Model Spread in Vegetation Responses to Elevated CO 2

Increasing atmospheric CO 2 and associated global warming are expected to alter the global hydrological cycle, thereby posing widespread threats to freshwater availability. However, future hydrological projections differ greatly between models, particularly over the tropical regions. The large difference between model projections directly limits policy planning efforts, and the responsible modeling processes remain unclear. Here, we identify the primary processes accounting for model differences in tropical hydrological changes using multiple CO 2 sensitivity experiments in the Coupled Model Intercomparison Project. We show that differences in projected changes to tropical evapotranspiration, precipitation, and surface water availability mainly arise from model representations of vegetation cover and stomatal conductance responses to elevated CO 2 and associated changes in atmospheric moisture and circulation. Atmospheric responses to sea surface warming contribute additionally to the divergence in hydrological projections. Given the importance of vegetation responses to elevated CO 2 and associated atmosphere feedbacks, our results underscore the need to improve representations of the vegetation physiological response to rising CO 2 and its coupling to the atmosphere, to provide reliable tropical hydrological projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

GCAM-USA: IM3 Phase 2 Official Simulations

Overview GCAM-USA simulates 21st century energy-water-land interactions in global regions with sub-national detail in the United States (US). The IM3 project is focused on US multisector dynamics and has created eight GCAM-USA simulations by combining two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight GCAM-USA simulations are named: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 Detailed Assumptions The SSPs describe alternative socioeconomic futures and their implications for both the US and global environmental challenges. As implemented in GCAM-USA, the primary difference between SSP3 and SSP5 in the US is their population growth projections, with SSP3 foreseeing lower growth and SSP5 projecting higher growth. The four climate projections explore moderate (RCP4.5) and high (RCP8.5) greenhouse gas emission (GHG) trajectories, each with a hotter and cooler climate variation based on the range of global climate model simulations in CMIP6s. These four climate projections are used in precursor models to develop climate impacts on building heating and cooling degree-hours, water availability, and agricultural yields that are provided as inputs to GCAM-USA. GCAM-USA reproduces the RCP4.5 and RCP8.5 GHG emissions trajectories consistent with each climate projection, using dynamic carbon prices to achieve RCP4.5. When the four climate projections are combined with the two SSP projections, the resulting eight scenarios allow GCAM-USA to explore how a wide, yet plausible, range of climate and socioeconomic impacts affects the interconnected energy, water, and land systems across the contiguous United States. Technical Details IM3 uses a version of GCAM-USA v5.3 that has an enhanced electricity sector and an updated data system needed to represent regional to local scale dynamics.

IM3↗
DOE OSTI2025

Baseflow Identification via Explainable AI With Kolmogorov‐Arnold Networks

Abstract Hydrological models often involve constitutive laws that may not be optimal in every application. We propose to replace such laws with the Kolmogorov‐Arnold networks (KANs), a class of neural networks designed to identify symbolic expressions. We demonstrate KAN's potential on the problem of baseflow identification, a notoriously challenging task plagued by significant uncertainty. KAN‐derived functional dependencies of the baseflow components on the aridity index outperform their original counterparts; they demonstrate that water availability, rather than potential evapotranspiration, drives baseflow by constraining actual evapotranspiration under arid conditions. On a test set, they increase the Nash‐Sutcliffe efficiency (NSE) by 65%, decrease the root mean squared error by 29%, and increase the Kling‐Gupta efficiency by 34%. This superior performance is achieved while reducing the number of fitting parameters from three to two. Next, we use data from 378 catchments across the continental United States to refine the water‐balance equation at the mean‐annual scale. The KAN‐derived equations based on the refined water balance outperform both the current aridity index model, with up to a 105% increase in NSE, and the KAN‐derived equations based on the original water balance. While the performance of our model and tree‐based machine learning methods is similar, KANs offer the advantage of simplicity and transparency and require no specific software or computational tools. This case study focuses on the aridity index formulation, but the approach is flexible and transferable to other hydrological processes. Plain Language Summary Equations used in hydrologic model are often suboptimal, resulting in reduced prediction accuracy and efficiency. We implemented Kolmogorov‐Arnold networks (KAN), a machine learning algorithm for deriving symbolic formulations, to estimate groundwater recharge and showed that it outperforms an existing state‐of‐the‐art semi‐empirical formulation. In hydrology, Nash‐Sutcliffe efficiency (NSE), root mean squared error (RMSE), and Kling‐Gupta efficiency (KGE) are commonly used to evaluate model performance. Higher NSE and KGE values indicate better performance, while lower RMSE values are preferable. Our results show that NSE increased by 71%, RMSE decreased by 32%, and KGE improved by 25%. In addition, KAN identifies an optimal functional form and can be used to derive new analytical formulas using the prior knowledge. The KAN‐inspired equation outperformed the original formulation and reduced the fitting parameters. Furthermore, we refined the water‐balance equation at the mean‐annual scale and showed that, based on the new water‐balance equation, KAN can derive new formulations that are superior to the original aridity index formulations (up to 105% increase in NSE) and KAN‐derived equations based on the original water balance. These findings highlight the significant potential of KAN to advance the scientific understanding of a wide range of hydrologic processes. Key Points Kolmogorov‐Arnold networks (KANs) enhance interpretability of machine‐learned hydrological models KAN‐derived symbolic formulations outperform state‐of‐the‐art semi‐empirical aridity indices KAN‐identified functional form yields an analytical index with fewer fitting parameters and improved performance

baseflow↗
DOE OSTI2022

Preparing the Power Sector to Navigate Climate and Water Risks

As the climate changes, so will water availability and other ambient conditions, threatening the reliability of today's power system and tomorrow's clean energy grid. That's why researchers at the National Renewable Energy Laboratory (NREL) - along with those at Sandia National Laboratories, Oak Ridge National Laboratory, the City University of New York, the Electric Power Research Institute, and the National Energy Technology Laboratory - are studying both regional and national climate and hydrologic changes to provide a comprehensive assessment of climate and water impacts and risks to the U.S. power grid.

ENERGY PLANNING, POLICY, AND ECONOMY,HYDRO ENERGY↗
DOE OSTI2023

Root exudates induced coupled carbon and phosphorus cycling in a soil with low phosphorus availability

The amount and type of root exudates can influence P availability in the rhizosphere directly by desorption or dissolution of soil minerals, or indirectly by decomposition of soil organic matter (SOM). Here this study aimed to determine the mechanisms by which specific root exudates influence the distribution and availability of P in soils with low P availability. Water, glucose, alanine, and oxalate were delivered through a simulated root into soils for 15 days. Zymography and planar optodes were used to image potential phosphatase activity, and O2 and pH distribution, respectively. Soils were analyzed for resin extractable inorganic P (Pi), dissolved organic C (DOC), water soluble Fe, and Al, and microbial community structure. Characterization of SOM and P were conducted using ultra-high resolution mass spectrometry and 31P solution nuclear magnetic resonance (NMR), respectively. The addition of oxalate resulted in the greatest resin extractable Pi, DOC, and water-soluble Fe, and Al compared to the other exudates suggesting destabilization of mineral associated organic matter (MAOM) and release of organic P (Po). Both 31P solution NMR and ultra-high resolution mass spectrometry analysis provided evidence of mineralization of Po released from the destabilization of MAOM. The study demonstrates the important role microbial and plant-derived metal chelating ligands play in destabilizing MAOM, releasing SOM and importantly Po, that when mineralized may contribute to increasing Pi availability in soils with low P availability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Identifying Hydropower Operational Flexibilities in Presence of Streamflow and Net-load Uncertainty (Final Technical report)

In the existing operations, hydropower contributions to future system flexibility are generally modeled while maintaining traditional operating rules and constraints in supporting grid operation, such as the balancing of variable renewable energy production. Moreover, operation of large scale hydropower systems on major rivers has been investigated for decades, utilizing various systems engineering approaches, with the evolving electric grid, as the result of renewable resources integration, compounded by the changing climate (variability of river flows, intensification of hydrologic cycle resulting in more frequent extreme events) affecting water availability, the need for more advanced stochastic modeling and effective uncertainty analysis approaches have become necessary. The research results supported by this funding and presented in this report provide a new look at hydropower operational flexibility enforced by the changes identified above. Understanding how hydropower operates in response to the underlying uncertainties with respect to the system constraints is crucial in identifying its operational flexibility potentials. In this project, the flexibility of the operating hydropower facility is described by capturing uncertainties in both water and power system and formulating the operations as a multistage stochastic optimization problem. The proposed approach supports short- to seasonal-term operations and planning decision horizons.

13 HYDRO ENERGY↗
DOE OSTI2024

Osmosensor-mediated control of Ca 2+ spiking in pollen germination

Higher plants survive terrestrial water deficiency and fluctuation by arresting cellular activities (dehydration) and resuscitating processes (rehydration). However, how plants monitor water availability during rehydration is unknown. Although increases in hypo-osmolarity-induced cytosolic Ca 2+ concentration (HOSCA) have long been postulated to be the mechanism for sensing hypo-osmolarity in rehydration, the molecular basis remains unknown. Because osmolarity triggers membrane tension and the osmosensing specificity of osmosensing channels can only be determined in vivo, these channels have been classified as a subtype of mechanosensors. Here we identify bona fide cell surface hypo-osmosensors in Arabidopsis and find that pollen Ca 2+ spiking is controlled directly by water through these hypo-osmosensors—that is, Ca 2+ spiking is the second messenger for water status. We developed a functional expression screen in Escherichia coli for hypo-osmosensitive channels and identified OSCA2.1, a member of the hyperosmolarity-gated calcium-permeable channel (OSCA) family of proteins. We screened single and high-order OSCA mutants, and observed that the osca2.1/osca2.2 double-knockout mutant was impaired in pollen germination and HOSCA. OSCA2.1 and OSCA2.2 function as hypo-osmosensitive Ca 2+ -permeable channels in planta and in HEK293 cells. Decreasing osmolarity of the medium enhanced pollen Ca 2+ oscillations, which were mediated by OSCA2.1 and OSCA2.2 and required for germination. OSCA2.1 and OSCA2.2 convert extracellular water status into Ca 2+ spiking in pollen and may serve as essential hypo-osmosensors for tracking rehydration in plants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Temporal, Spatial, and Temperature Controls on Organic Carbon Mineralization and Methanogenesis in Arctic High-Centered Polygon Soils

Warming temperatures in continuous permafrost zones of the Arctic will alter both hydrological and geochemical soil conditions, which are strongly linked with heterotrophic microbial carbon (C) cycling. Heterogeneous permafrost landscapes are often dominated by polygonal features formed by expanding ice wedges: water accumulates in low centered polygons (LCPs), and water drains outward to surrounding troughs in high centered polygons (HCPs). These geospatial differences in hydrology cause gradients in biogeochemistry, soil C storage potential, and thermal properties. Presently, data quantifying carbon dioxide (CO 2 ) and methane (CH 4 ) release from HCP soils are needed to support modeling and evaluation of warming-induced CO 2 and CH 4 fluxes from tundra soils. This study quantifies the distribution of microbial CO 2 and CH 4 release in HCPs over a range of temperatures and draws comparisons to previous LCP studies. Arctic tundra soils were initially characterized for geochemical and hydraulic properties. Laboratory incubations at −2, +4, and +8°C were used to quantify temporal trends in CO 2 and CH 4 production from homogenized active layer organic and mineral soils in HCP centers and troughs, and methanogen abundance was estimated from mcrA gene measurements. Results showed that soil water availability, organic C, and redox conditions influence temporal dynamics and magnitude of gas production from HCP active layer soils during warming. At early incubation times (2–9 days), higher CO 2 emissions were observed from HCP trough soils than from HCP center soils, but increased CO 2 production occurred in center soils at later times (>20 days). HCP center soils did not support methanogenesis, but CH 4 -producing trough soils did indicate methanogen presence. Consistent with previous LCP studies, HCP organic soils showed increased CO 2 and CH 4 production with elevated water content, but HCP trough mineral soils produced more CH 4 than LCP mineral soils. HCP mineral soils also released substantial CO 2 but did not show a strong trend in CO 2 and CH 4 release with water content. Knowledge of temporal and spatial variability in microbial C mineralization rates of Arctic soils in response to warming are key to constraining uncertainties in predictive climate models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Global and regional hydrological impacts of global forest expansion

Abstract. Large-scale reforestation, afforestation, and forest restoration schemes have gained global support as climate change mitigation strategies due to their significant carbon dioxide removal (CDR) potential. However, there has been limited research into the unintended consequences of forestation from a biophysical perspective. In the Community Earth System Model version 2 (CESM2), we apply a global forestation scenario, within a Paris Agreement-compatible warming scenario, to investigate the land surface and hydroclimate response. Compared to a control scenario where land use is fixed to present-day levels, the forestation scenario is up to 2 °C cooler at low latitudes by 2100, driven by a 10 % increase in evaporative cooling in forested areas. However, afforested areas where grassland or shrubland are replaced lead to a doubling of plant water demand in some tropical regions, causing significant decreases in soil moisture (∼ 5 % globally, 5 %–10 % regionally) and water availability (∼ 10 % globally, 10 %–15 % regionally) in regions with increased forest cover. While there are some increases in low cloud and seasonal precipitation over the expanded tropical forests, with enhanced negative cloud radiative forcing, the impacts on large-scale precipitation and atmospheric circulation are limited. This contrasts with the precipitation response to simulated large-scale deforestation found in previous studies. The forestation scenario demonstrates local cooling benefits without major disruption to global hydrodynamics beyond those already projected to result from climate change, in addition to the cooling associated with CDR. However, the water demands of extensive forestation, especially afforestation, have implications for its viability, given the uncertainty in future precipitation changes.

Environmental Sciences & Ecology↗
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 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↗