Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Information Sharing”

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

wavess 1.2: presenting an HLA-aware within-host virus sequence simulation framework

Motivation Understanding how virus sequences are shaped by selection can inform vaccine design and transmission inference. Modeling within-host evolution to interrogate these questions requires a detailed mechanistic framework that accurately captures sequence diversification. The CD8 + cytotoxic T-lymphocyte (CTL) response plays an important role in immune-mediated selection and can leave strong signatures in virus sequences; however, existing sequence-based within-host virus modeling frameworks do not explicitly include a human leukocyte antigen (HLA)-aware CTL response. Results We extended our previously published within-host sequence evolution simulator, wavess, to include an explicit CTL response, and share a method for identifying HLA-specific CTL epitopes given a founder virus sequence. We also updated the model to permit a variable recombination rate, which allows for modeling non-adjacent genes, segmented genomes, and recombination hotspots. These extensions to wavess allow for more accurate simulation of viruses and virus genes, particularly in regions of the genome where the immune response is dominated by CTLs (rather than antibodies). It also provides the foundation for investigations of how these newly-added biological mechanisms influence within-host evolution. Availability and implementation The core of wavess is written in Python 3, with helper functions written in R. It is available at https://github.com/MolEvolEpid/wavess.

60 APPLIED LIFE SCIENCES↗

State-Selective Double Photoionization of Atomic Carbon and Neon

Double photoionization (DPI) allows for a sensitive and direct probe of electron correlation, which governs the structure of all matter. For atoms, much of the work in theory and experiment that informs our fullest understanding of this process has been conducted on helium, and efforts continue to explore many-electron targets with the same level of detail to understand the angular distributions of the ejected electrons in full dimensionality. Expanding on previous results, we consider here the double photoionization of two 2p valence electrons of atomic carbon and neon and explore the possible continuum states that are connected by dipole selection rules to the coupling of the outgoing electrons in 3P, 1D, and 1S initial states of the target atoms. Carbon and neon share these possible symmetries for the coupling of their valence electrons. Results are presented for the energy-sharing single differential cross section (SDCS) and triple differential cross section (TDCS), further elucidating the impact of the initial state symmetry in determining the angular distributions that are impacted by the correlation that drives the DPI process.

Yip, Frank L. (ORCID:000000025409495X)↗

Nuclear Criticality Safety Integral Experiment Covariance Determination

Integral benchmarks for criticality safety and nuclear data validation require expensive uncertainty quantification studies. Commonly, the uncertainty quantification ignores correlations between experiments that share components. Experiments such as the TEX (Thermal/Epithermal eXperiments) campaigns consist of many shared parts, such as fuel, which create a strong correlation in their uncertainties. While these correlations are known to exist, they are often not estimated due to the complexity of such calculations. This paper describes a software package that uses an intuitive method of determining the covariance for each of the experimental components, providing a correlation matrix for each family of parts across the multiple cases examined within a benchmark. The code uses the TEX-HEU campaign as a proof of concept, and we show that the correlations can be calculated with information commonly found in ICSBEP (International Criticality Safety Benchmark Evaluation Project) benchmarks. The estimated covariances are used in χ 2 trending studies to evaluate their impact on nuclear data validation. Without covariances, χ 2 per degree of freedom was calculated as 2.203 and with covariances it was 1.179. The difference shows that omitting covariance information may cause overly pessimistic bias quantifications. The covariance determination code can be easily integrated into current benchmark evaluations as well as reevaluating legacy benchmark uncertainties. Uncertainty correlation calculations should become the baseline for criticality safety integral experiment benchmarks and can now be easily calculated with the described software package.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Exploiting CRIC to streamline the configuration management of GlideinWMS factories for CMS support

GlideinWMS is a workload management and provisioning system that allows sharing computing resources distributed over independent sites. Based on the requests made by GlideinWMS frontends, a dynamically sized pool of resources is created by GlideinWMS pilot factories via pilot job submission to resource sites’ CEs. More than 400 CEs are currently serving more than ten virtual organizations through GlideinWMS, with CMS being the biggest user with 230 CEs. The complex configurations of the parameters defining resource requests, as submitted to those CEs, have been historically managed by manually editing a set of different XML files. New possibilities arise with CMS adopting the CRIC, an information system that collects, aggregates, stores, and exposes, among other things, computing resource data coming from various data providers. The paper will describe the challenges faced when CMS started to use CRIC to automatically generate the GlideinWMS factory configurations. The architecture of the prototype, and the ancillary tools developed to ease this transition, will be discussed. Finally, future plans and milestones will be outlined.

Dost, Jeffrey↗

Science and Solutions for a Planet Under Pressure: Synthetic Research, Analytics, and Community Engagement on Five Priority Theme Areas

Support from the Department of Energy proved instrumental in supporting the work of Global Council for Science and the Environment (GCSE) on synthetic research, analysis, and community engagement, oriented around the themes and discussions from the GCSE Drawdown 2021 Conference — Research to Action: Science and Solutions for a Planet Under Pressure. These efforts helped inform and enhance the societal relevance of the knowledge that has been generated through the conference and, at the same time, advance the scientific insights acquired through synthesis and scientific translation efforts. The work centered on three main objectives: 1) documentation of findings from the completed GCSE 2021, culminating in the publication of a synthesis report; 2) development and vetting of integrated use cases that illustrate synthesis methods and enable scientific outreach, communications, and translational discussions among universities and the broader research community; and 3) strengthening the scientific basis for energy education research by sharing findings in university-educator fora.

54 ENVIRONMENTAL SCIENCES↗

Small Hydropower Interconnections: Best Practices

Small hydropower projects have been the predominant source of capacity growth of U.S. hydropower for more than a decade, and they present the most cost-effective and environmentally permissible avenues for hydropower growth (DOE 2016; Johnson et al. 2018). However, interconnection to electricity distribution and transmission grids is a persistent barrier due to cost surprises and schedule overruns. As a culmination to research into the status and requirements of small hydropower interconnection across the United States, this paper presents the best practices for setting interconnection standards that can improve the process for small hydropower developers. As part of the analysis, the interconnection costs are compared between small hydropower, solar, and wind. The analysis of the small hydropower interconnection landscape across the United States was carried out by Pacific Northwest National Laboratory (PNNL) and Oak Ridge National Laboratory (ORNL) with support from the U.S. Department of Energy Water Power Technologies Office. The research team was guided by a Technical Advisory Group (TAG) and gleaned data from publicly available sources, such as the HydroSource database (ORNL 2020) and interconnection queues hosted by utilities, balancing authorities, independent system operators (ISOs), and regional transmission organizations (RTOs). The results of this work are shared in a series of papers detailing the state of small hydropower in the United States (“Small Hydropower Interconnections: Small Hydropower in the United States”), the variety of state interconnection processes to connect power generators with the grid (“Small Hydropower Interconnections: State Interconnection Processes”), and an analysis of the interconnection processes (“Small Hydropower Interconnections: Analysis of Interconnection Processes”). In this, the final paper in the series, best practices for interconnection processes (“Small Hydropower Interconnections: Best Practices”) are identified from the solar energy and distributed wind energy industries that are transferrable to small hydropower development. This information will help overcome barriers to future small hydropower development.

13 HYDRO ENERGY↗

Nuclear Criticality Safety Integral Experiment Covariance Determination

Integral benchmarks for criticality safety and nuclear data validation require expensive uncertainty quantification studies. Commonly, the uncertainty quantification ignores correlations between experiments that share components. Experiments such as the TEX (Thermal/Epithermal eXperiments) campaigns consist of many shared parts between experiments, such as fuel, which creates a strong correlation in their errors. While these correlations are known to exist, they are often not estimated due to the complexity of such calculations. This paper describes a software package that uses an intuitive method of determining the covariance for each of the experimental components, providing a correlation matrix for each family of parts across the multiple cases examined within a benchmark. The code uses the TEX-HEU campaign as a proof of concept, and we show that the correlations can be calculated with information commonly found in ICSBEP (International Criticality Safety Benchmark Evaluation Project) benchmarks. The estimated covariances are used in χ 2 trending studies to evaluate their impact on nuclear data validation. The covariance determination code can be easily integrated into current benchmark evaluations as well as reevaluating legacy benchmark uncertainties. Uncertainty correlation calculations should become the baseline for criticality safety integral experiment benchmarks and can now be easily calculated with the described software package.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hot Springs and Geysers: Exploring Historical and Modern Impacts of Geothermal Energy Production on Associated Natural Surface Systems and Standardizing Management Practices

Surface thermal features, most notably hot springs and geysers are increasingly being recognized for their importance to ecosystems, indigenous cultures, and in some cases agriculture, recreation, and tourism. Geothermal project development poses a potential risk to these natural features but current regulatory requirements for assessing and managing these risks during exploration, permitting and monitoring are somewhat inconsistent and unpredictable across different geothermal fields. This has resulted in uncertainty and increases in exploration risk for geothermal energy developers that have led to costly project delays, cancellations, or hesitation to commit. Varying regulatory requirements may also influence public perception, fostering confusion, distrust and ultimately opposition to geothermal projects, further contributing to project delays or cancellations. At a time when there is an increasing urgency for reliable baseload clean energy, geothermal is a net-zero, renewable solution that additionally provides access to more equitable and environmentally just clean power. Continued integration of geothermal energy into the national energy roadmap can be facilitated through consistent and predictable permitting, providing regulators the framework they need, developers a clear path forward, and transparency that the public deserves. This project, currently in its beginning phases, seeks to address this important issue by providing a technical basis from which to build a preliminary protocol for assessing and managing potential impacts from new or existing geothermal energy projects to surface thermal features and their associated ecosystems. Development of this preliminary protocol will be informed by (1) a literature review of well-documented case studies in the western U.S. and New Zealand to understand the range of conditions that exemplify geothermal-surface thermal systems; (2) development of generic illustrative conceptual-numerical models to quantify, understand, and predict the first-order controls (e.g., pressure and permeability) on surface flows; and (3) additional independent and scientifically rigorous evaluations of geothermal-surface thermal system case studies from the Basin and Range Province that incorporate publicly available data as well as data provided by industry through data-sharing agreements. Learning from the successes of the process used to develop the Induced Seismicity Management Protocol (ISMP), we ultimately aim to use these initial efforts as a springboard for establishing a surface thermal feature management working group that will work collaboratively to finalize the protocol as well as co-create recommended best practices for implementation. We envision that the working group will primarily be composed of representatives from regulatory entities, government agencies, Tribes, academia, national laboratories, and industry, and will include early and regular engagement with community organizations and environmental groups. This will help ensure broad acceptance and implementation of the protocol, which will facilitate a more consistent, predictable, and standardized regulatory process, and help to ensure that geothermal energy continues to provide a reliable source of clean energy, and a pathway to achieving greater energy equity in the U.S.

Best Practices↗

Modeling direct air carbon capture and storage in a 1.5 °C climate future using historical analogs

Limiting the rise in global temperature to 1.5 °C will rely, in part, on technologies to remove CO 2 from the atmosphere. However, many carbon dioxide removal (CDR) technologies are in the early stages of development, and there is limited data to inform predictions of their future adoption. Here, we present an approach to model adoption of early-stage technologies such as CDR and apply it to direct air carbon capture and storage (DACCS). Our approach combines empirical data on historical technology analogs and early adoption indicators to model a range of feasible growth pathways. We use these pathways as inputs to an integrated assessment model (the Global Change Analysis Model, GCAM) and evaluate their effects under an emissions policy to limit end-of-century temperature change to 1.5 °C. Adoption varies widely across analogs, which share different strategic similarities with DACCS. If DACCS growth mirrors high-growth analogs (e.g., solar photovoltaics), it can reach up to 4.9 GtCO 2 removal by midcentury, compared to as low as 0.2 GtCO 2 for low-growth analogs (e.g., natural gas pipelines). For these slower growing analogs, unabated fossil fuel generation in 2050 is reduced by 44% compared to high-growth analogs, with implications for energy investments and stranded assets. Residual emissions at the end of the century are also substantially lower (by up to 43% and 34% in transportation and industry) under lower DACCS scenarios. The large variation in growth rates observed for different analogs can also point to policy takeaways for enabling DACCS.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Future Electric Power Industry and Grids—Now What Again is Our Destination?

Having a vision that others agree to support and work toward is highly desirable but hard to achieve. We seem to lack a common and shared understanding of the vision—or worse, multiple visions (vivid mental images or documented statements) with varying areas of focus and details: • some appear to be similar but have differing underlying goals and characteristics, or • some reflect differing viewpoints as to effects on various stakeholders. These desirable and undesirable situations apply to realizing visions for enterprise and industry, including the electricity sector. Multiple industry stakeholder groups have developed goals, industry vision statements, and characterizations of the future. The viewpoints are promoted, discussed, refined by their stakeholder group, and often published to promote broad understanding and to inform or influence others. The GridWise Architecture Council asked itself how well-aligned these characterizations of the future are. If these publications collectively set the overall direction for the industry, it is useful to identify their answers to questions such as, where is the electric industry headed, guided by what objectives, and with what role(s) for the customers, electric utilities, and other stakeholders? Are these goals, visions, and future states moving toward a common vision, do they provide value for the stakeholders, and are they likely to meet the objective stated? This paper addresses these questions via an assessment and characterization of a sampling of stakeholder groups’ publicly available vision and future state reports for the electricity industry. Identified electric power grid architectural topic areas needing further work are described, along with GridWise Architecture Council analysis and observation, potential collaborative work efforts, and suggested next steps.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Phenome‐to‐genome insights for evaluating root system architecture in field studies of maize

Abstract Understanding the genetic basis of root system architecture (RSA) in crops requires innovative approaches that enable both high‐throughput and precise phenotyping in field conditions. In this study, we evaluated multiple phenotyping and analytical frameworks for quantifying RSA in mature, field‐grown maize in three field experiments. We used forward and reverse genetic approaches to evaluate >1700 maize root crowns, including a diversity panel, a biparental mapping population, and maize mutant and wild‐type alleles at two known RSA genes,DEEPER ROOTING 1(DRO1) andRootless1(Rt1). We show the utility of increasing the dimensionality of traditional two‐dimensional (2D) techniques, referred to as the “2D multi‐view” method, to improve the capture of whole root system information for mapping genetic variation influencing RSA. Comparison of univariate and multivariate genome‐wide association study (GWAS) approaches revealed that multivariate traits were effective at dissecting complex RSA phenotypes and identifying pleiotropic quantitative trait loci (QTLs). Overall, three‐dimensional (3D) root models generated from X‐ray computed tomography and digital phenotyping captured a larger proportion of RSA trait variations compared to other methods of root phenotyping, as evidenced by both genome‐wide and single‐gene analyses. Among the individual root traits, root pulling force emerged as a highly heritable estimate of RSA that identified the largest number of shared QTLs with 3D phenotypes. Our study shows that integrating complementary phenotyping technologies helps to provide a more comprehensive understanding of the genetic architecture of RSA in field‐grown maize.

Genetics & Heredity↗

Modified Andronov-Hopf Oscillator-Based Grid-Forming Converter with Emulated Virtual Cable for Enhanced Power Sharing Performance

Nonlinear oscillator-based grid-forming converters offer superior dynamic and steady-state performance, making them an attractive solution for interconnecting renewable resources. This paper proposes a novel modified Andronov-Hopf oscillator to enhance the operating spectrum and facilitate the integration of renewable energy sources. An inner loop controller based on the Lyapunov energy function is implemented to achieve robust stability and performance, while a virtual cable emulation strategy enables seamless parallel operation. Comprehensive modeling and simulation studies validate the effectiveness of the proposed system, demonstrating its capabilities in addressing diverse operating scenarios, including grid faults, renewable energy fluctuations, and parallel operation. The proposed solution exhibits fast transient response, robust stability, and flexible operation, making it a valuable contribution to the field of renewable energy integration. The results of this study can be used to inform the design and implementation of next-generation grid-forming converters, enabling a more sustainable and reliable energy future. Additionally, the proposed system's ability to operate in both grid-connected and islanded modes makes it an ideal candidate for remote and off-grid renewable energy applications. The proposed solution's scalability and modularity also make it suitable for large-scale renewable energy integration. The proposed system is verified through MATLAB/Simulink and PLECS simulations, demonstrating its effectiveness in ensuring robust and efficient operation.

Andronov-Hopf Oscillator (AHO)↗

Experiences and Lessons from Field Demonstration of Grid-forming Inverter in An AC Microgrid

The validation of GFM control strategies through simulation and hardware demonstration is important before their large-scale deployments in the real grid. Considering the importance of testing and validation, several works have explored the GFM inverter’s capability to blackstart a microgrid, synchronize and share loads, and interact with various types of generation sources and loads present in the grid. Herein, this paper complements the existing works by presenting the results and analysis of a field demonstration in an actual AC microgrid. The capability of a three-level neutral point clamped (NPC) GFM inverter equipped with a recursive feedback type of non-linear device level control to operate with PV source on its DC input and off-the-shelf PVGFL inverters and EV chargers of different kinds on the AC side is explored. The analysis and conclusions drawn would inform the readers to make better decisions during the field demonstration process.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Coal-dependent Communities in Transition: Identifying Best Practices to Ensure Equitable Outcomes

The U.S. coal industry is experiencing a sharp increase in the numbers of retired and/or decommissioned coal-fired power plants across the U.S. In the years between 2010-2019, around 102 gigawatts (GW) of coal-fired generating capacity has been announced to be decommissioned, representing more than 546 coal-fired power plant units, and an additional 17 GW is planned to be decommissioned by 2025. This change in the energy production landscape presents an impact on the social, environmental, and economic prospects of coal-dependent communities. This report examined the role of communities in the coal power plant decommissioning process and provided community-identified best practices to ensure an equitable process. The experiences of four coal-dependent communities—Wise County, VA, Muskegon, MI, Anderson County, TN, and Becker, MN—are presented as case studies to understand the impacts of the decommissioning process, and associated best practices, from the communities’ perspective. The report results highlight the need to recognize that the decommissioning decision-making process must be community-based to be equitable. Each community’s input is key to the transition away from coal power because there is no one-size-fits-all development plan. In other words, each community’s trajectory through the decommissioning process—from the retirement decision-making stage to the final site redevelopment phase—is unique because each community has distinct needs and wants from the energy transition. What is best for one community may not be suited for another. Ultimately, the framework for site development and community revitalization post-decommissioning cannot be universal because each community’s profile—from a social, cultural, and economic perspective—is different. Community impacts of power plant decommissioning are not limited to job and revenue losses. Communities are likely to be impacted culturally, socially, environmentally, and have long-term health-based impacts that should be acknowledged and addressed in post-retirement plans. Commonly identified decommissioning best practices include: Early and continued engagement throughout, with a number of mediums for communication and feedback (e.g., in-person sessions, virtual meetings, written comment opportunities); Early planning of post-decommissioning projects to replace lost jobs, revenue, and economic activity; Recognition (and mitigation, if possible) of social impacts on the community due to plant closure; Transparency throughout the process, with trusted information being provided about the decommissioning process and timeline; potential impacts on the workforce, economy, and environment; and the feasibility of alternative site uses; Identification of funding sources, technical experts, and/or strategic partnerships to support decommissioning and the affected communities upfront; and, Acknowledgment of communities as stakeholders who have a role in the conversation and right to determine their futures. Three key areas for assisting coal-dependent communities affected by the energy transition: Technical assistance: assessment of site feasibility for alternative uses or to repower with new technologies; Cross-partnership engagement and collaboration: facilitate knowledge-sharing of “lessons learned” about the decommissioning process between communities and provide guidance for decision-making processes; and, Financial assistance: access to grant and/or loan programs to assist with redevelopment survey, bolster community economic security through job creation, and cover environmental clean-up costs. Technical assistance, cross-partnership engagement, and financial aid can be mobilized to help communities throughout various stages of the decommissioning process, including the retirement decision, the site reclamation phase, and eventual revitalization of the site and surrounding community.

01 COAL, LIGNITE, AND PEAT↗

Developing Bloom Filters for Web Archives’ Holdings (Final Project Report)

The main goal of the project was to develop a framework for web archives to create Bloom filters based on their holdings of archived web resources. A Bloom filter (BF) is a data structure that, in our scenario, contains hash values of all (or a subset of) URLs, of which an archive has one or more archival copies. Two main use cases fall in scope for this project and are supported by a BF implementation: 1) Sharing of an archive’s holdings (URLs) and 2) Querying the holdings of one or more archives. Since URL strings are hashed before ingested into the BF, the index of an archive is not shared in plain text when the BF is shared with trusted parties. On the other hand, a BF implementation does allow queries for a URL to confirm if an archive indeed has one or more archival copies of that URL. This collaborative project between Los Alamos National Laboratory (LANL) and the Croatian Web Archive (HAW), from the National and University Library in Zagreb (NSK), developed by NSK and University of Zagreb University Computing Center SRCE), aimed at developing software to create BFs, evaluate the scalability of the approach, pilot a search service based on BFs, and design a framework for archives to share their filters with trusted parties. In the remainder of this document, we will report on the work completed with respect to the individual deliverables, outline aspects of future work, and conclude with our recommendations for the use of BFs for the web archiving community.

97 MATHEMATICS AND COMPUTING↗

Pooled Rideshare in the U.S.: An Exploratory Study of User Preferences

Pooled ridesharing offers on-demand, one-way, cost-effective transportation for passengers traveling in similar directions via a shared vehicle ride with others they do not know. Despite its potential benefits, the adoption of pooled rideshare remains low in the United States. This exploratory study aims to evaluate potential service improvements and features that may increase users’ willingness to adopt the service. The study analyzed transportation behaviors, rideshare preferences, and willingness to adopt pooled rideshare services among 8296 U.S. participants in 2025, building on findings from a 2021 nationwide survey of 5385 U.S. participants. The study incorporated 77 actionable items developed from the results of the 2021 survey to assess whether addressing specific user-generated topics such as safety, reliability, convenience, and privacy can improve pooled rideshare use. A side-by-side comparison of the 2021 and 2025 data revealed shifts in transportation behavior, with personal rideshare usage increasing from 22% to 28%, public transportation from 21% to 27%, and pooled rideshare from 6% to 8%, while personal vehicle (79%) use remained dominant. Participants rated features such as driver verification (94%), vehicle information (93%), peak time reliability (93%), and saving time and money (92–93%) as most important for improving rideshare services. A pre-to-post analysis of willingness to use pooled rideshare utilizing the actionable items as per respondents’ preferences showed improvement: “definitely will” increased from 15.9% to 20.1% and “probably will” rose from 35.6% to 47.7%. These results suggest that well-targeted service improvements may meaningfully enhance pooled rideshare acceptance. This study offers practical guidance for Transportation Network Companies (TNCs) and policymakers aiming to improve pooled rideshare as well as potential future research opportunities.

Transportation Network Companies (TNCs)↗

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↗

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↗