Engineering PapersSearch

SEARCH · Engineering Papers

Results for “AWARE campaign”

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.

ARM FY2025 Radar Plan

The U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility is dedicated to delivering high-quality radar data that significantly advance our understanding of cloud and precipitation processes, improving climate models. With radars operating across diverse frequencies, scanning modes, and global climate conditions, extensive staffing is essential for effective management. Due to current staffing constraints, achieving the expected level of operational excellence requires a more strategic approach. To address this challenge, ARM has developed an operational radar plan for the upcoming Fiscal Year 2025 (FY25) based on budget and staffing considerations. This report summarizes the Fiscal Year 2024 (FY24) ARM radar-related activities and outlines the strategic plan for FY25. It presents a comprehensive radar plan that includes detailed activities, priorities, and a projected timeline aligned with the ARM radar roadmap. Key tasks, detailed in Table 1, cover various radar operational stages and coordinated efforts among ARM teams.

54 ENVIRONMENTAL SCIENCES

Deep Gaussian process-based cost-aware batch Bayesian optimization for complex materials design campaigns

The accelerating pace and expanding scope of materials discovery demand optimization frameworks that efficiently navigate vast design spaces with complex response surfaces while judiciously allocating limited evaluation resources. We present a cost-aware, batch Bayesian optimization scheme powered by deep Gaussian process (DGP) surrogates and a heterotopic querying strategy. Our DGP surrogate, formed by stacking GP layers, models complex hierarchical relationships among high-dimensional compositional features and captures correlations across multiple target properties, propagating uncertainty through successive layers. We integrate evaluation cost into an upper-confidence-bound acquisition extension, which, together with heterotopic querying, proposes small batches of candidates in parallel, balancing exploration of under-characterized regions with exploitation of high-mean, low-variance predictions across correlated properties. Applied to refractory high-entropy alloys for high-temperature applications, our framework converges to optimal formulations in fewer iterations with cost-aware queries than conventional GP-based BO, highlighting the value of deep, uncertainty-aware, cost-sensitive strategies in materials campaigns.

36 MATERIALS SCIENCE

Bayesian Exploration and Surrogate Emulation of Nonlinear Beam-Response Geometry in the LBNF Beamline

Next-generation long-baseline neutrino experiments aim to achieve multi-MW proton beam power while reducing accelerator-induced systematic uncertainties. At Fermilab, the LBNF beamline is designed for 1.2 MW operation with PIP-II and is upgradeable to 2.4 MW. DUNE will probe the three-flavor neutrino paradigm and search for CP violation, requiring precise neutrino-flux normalization and improved control of accelerator-related uncertainties. Within the LBNF beamline, the System for On-Axis Neutrino Detection (SAND) will constrain flux uncertainties using precision near-detector measurements, while the Muon Monitor System (MuMS) will provide beamline diagnostics sensitive to the proton beam, target, and horn configuration. However, the pion phase space relevant for DUNE depends simultaneously on many correlated parameters, including beam centroid, beam width, horn current and alignment, target position, optics shifts, and radiation-induced changes. Consequently, MuMS observables exhibit nonlinear and coupled responses that are difficult to characterize using traditional one-parameter scans. To address this challenge, we are developing a Bayesian Exploration framework coupled to physics-informed surrogate emulators trained on Geant4 beamline simulations. Gaussian-process emulators provide both fast predictions and uncertainty estimates, enabling adaptive selection of new simulation points in beam-parameter space. As an initial demonstration, we construct surrogate emulators for MuMS response observables using a verified simulation campaign spanning proton-beam steering conditions. The emulators reproduce the simulated dependence of MuMS centroid and gradient observables while providing predictive uncertainties, and serve as the foundation for future multidimensional exploration including beam width, horn current, and additional beamline parameters. This work establishes a framework for uncertainty-aware beam monitoring, adaptive simulation campaigns, and rapid beam-response inference for future DUNE operations.

Ganguly, Sudeshna [Fermilab] (ORCID:00000003163482

Southern Hemisphere Polar Infrared Radiation Experiment (SHIRE) Field Campaign Report

An extended-range atmospheric emitted radiance interferometer (AERI) instrument to be deployed to Scott Base (78 S, 166 E) in McMurdo Sound, Antarctica for the purpose of synergizing with the National Aeronautics and Space Administration (NASA) Polar Radiant Energy in the Far-InfraRed Experiment (PREFIRE) satellite mission (https://prefire.ssec.wisc.edu/), examining the closure of atmospheric thermal emission, and providing validation for PREFIRE. The PREFIRE mission is the first satellite mission to map Earth’s spectrally resolved thermal emission, from 5 to 50 microns, over seasonal to annual timescales. These new observations represent the beginning of an intensive examination of the nature of thermal emission from Earth (emission from both the surface and atmosphere). A primary requirement of the PREFIRE mission is to provide validation of the retrieved variables. Numerous field sites in the Arctic (such as the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) North Slope of Alaska observatory and Summit Station, Greenland) have existing instrumentation to facilitate validation, leaving an observational gap in Antarctica. Through SHIRE, we propose to establish a temporary, intensive measurement site in Antarctica to ensure that both poles are suitably represented, and benefit from this state-of-the art science. This will also allow Arctic and Antarctic atmospheric radiative states to be compared and contrasted. For SHIRE specifically, we also propose further objectives to examine the spectrally resolved thermal radiative closure of the polar atmosphere and the historical changes in Antarctic radiative forcing, using previous Antarctic field campaigns, such as the ARM West Antarctic Radiation Experiment (AWARE).

54 ENVIRONMENTAL SCIENCES

Early Stages in the Lifecycle of Polar Liquid-Bearing Clouds

Stratiform liquid-bearing clouds are ubiquitous over the polar regions, where they are predominantly mixed-phase. These polar clouds induce substantial radiative forcing on the surface and continuously modify the atmospheric thermodynamic budget, with direct implications for the polar ice pack resilience. However, the physical representation of these clouds is still a major challenge for climate models. A significant part of polar liquid-bearing cloud lifecycle is often manifested in a quasi-steady self-sustaining, persistent, and turbulent state, which is driven by longwave cloud radiative cooling and can last for multiple days. This self-sustaining cloud lifecycle stage has been thoroughly investigated in numerous studies, though some of its aspects such as precipitation still lack robust quantification and evaluation. The preceding cloud lifecycle stages, which may last up to several hours, have nonetheless remained widely overlooked. These preceding stages initiate at cloud formation, often in a stable and non-turbulent atmospheric layer and serve as a key junction between cloud persistence and cloud dissipation. These two contrasting cloud lifecycle trajectories pose the question if general circulation models (GCMs) can capture the full lifecycle accurately for the real physical reasons, a necessary condition to improve our confidence in climate projections given the changing polar climate. The purpose of this project was to improve the characterization and understanding of these early stages in the lifecycle of polar stratiform liquid-bearing cloud and to aid their representation in GCMs. This research relied on measurements from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) field campaign, as well as the observations from the ARM West Antarctic Radiation Experiment (AWARE) and the permanent ARM site at Utqiagvik, North Slope of Alaska (NSA).

58 GEOSCIENCES

Advancing the Understanding of Cloud Microphysical Processes and Aerosol Indirect Effects in High-Latitude Mixed-Phase Clouds by Linking ARM Measurements with Climate Model Simulations (Final Report)

The key objectives of this project were to advance our understanding of cloud microphysical characteristics and aerosol indirect effects on mixed-phase clouds in high latitudes. To improve the representation of ice and mixed-phase clouds in Earth System Models (ESMs), we propose an integrated observation and modeling study of cloud macro- and microphysical properties, including spatial heterogeneities, mass partitioning between ice crystals and supercooled liquid water, effects of ice nucleating particles (INPs), and efficiency of secondary ice production (SIP), etc. Specifically, we took four main approaches in this project: (1) examining macro- and microphysical properties of ice and mixed-phase clouds based on in-situ and ground-based observations from multiple field campaigns funded by the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program, including the Mixed-Phase Arctic Cloud Experiment (M-PACE), Indirect and Semi-Direct Aerosol Campaign (ISDAC), Ice Nucleating Particle Sources at Oliktok Point (INPOP), ARM West Antarctic Radiation Experiment (AWARE), Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS), and Macquarie Island Cloud and Radiation Experiment (MICRE); (2) evaluating the DOE Energy Exascale Earth System Model (E3SM) simulations based on observations, particularly for ice and mixed-phase cloud microphysical properties; (3) examining the impacts of INPs on ice and mixed-phase clouds. Specifically, a series of comparisons were conducted using observations over the Arctic, Southern Ocean, and Antarctica, including comparisons between the lower and higher southern latitudes as well as comparisons between the northern and southern hemispheres. In addition, aerosol indirect effects from distinct sources of dust particles were examined; and (4) investigating the impacts of SIP. Ultimately, these results helped to improve cloud microphysics and aerosol-cloud interaction parameterizations in the E3SM model. Overall, the project provided improved understanding regarding various factors, including thermodynamic, dynamic, and aerosol conditions, on the micro- and macrophysical properties of ice and mixed-phase clouds in the high latitudes. Resulting analysis helped to provide an improved physical basis for refining the current cloud microphysics parameterizations related to ice and mixed-phase clouds in E3SM.

54 ENVIRONMENTAL SCIENCES

Uncertainty-Aware Machine Learning for Small-Angle X-ray Scattering Analysis in Autonomous Experimentation

Small-angle X-ray scattering (SAXS) is a powerful high-throughput characterization tool for probing nanoscale structure in native sample environments, providing real-time morphological information such as nanoparticle size and shape during synthesis. However, automated SAXS data analysis for extracting meaningful structural parameters is non-trivial and remains a bottleneck in closed-loop experimentation towards autonomous materials discovery, which demands fast, reliable, and uncertainty-aware data analysis. Here, we develop a machine-learning approach for automated SAXS analysis tailored to closed-loop nanoparticle synthesis. A Random Forest (RF) regression model is trained on 100,000 synthetic SAXS curves generated from polydisperse spherical nanoparticles with realistic background contributions. Using normalized one-dimensional SAXS intensity profiles as input, the RF model directly predicts nanoparticle radius, size polydispersity, and background parameters, while the ensemble standard deviation across trees provides built-in uncertainty quantification (UQ). On synthetic data, we show that combining fit-quality metrics (R 2 , MAE) with thresholds on prediction uncertainty reliably identifies accurate parameter estimates without access to ground truth. We then apply the trained model to 365 experimental SAXS profiles of citrate-reduced gold nanoparticles synthesized using an automated droplet-flow microreactor with in situ SAXS at a synchrotron beamline, classifying the results into high- and low-confidence subsets based on UQ metrics. Finally, we integrate RF-based SAXS analysis into a simulated closed-loop optimization campaign using Gaussian process Bayesian optimization to minimize nanoparticle polydispersity, benchmarking against conventional automated Levenberg–Marquardt fitting. The RF-guided campaign exhibits substantially faster convergence and lower relative opportunity cost (∼0.07 vs ∼0.3), demonstrating that uncertainty-aware machine-learning SAXS analysis significantly enhances the efficiency and robustness of autonomous nanomaterials synthesis workflows.

Bayesian optimization

Modular Autonomous Experimentation for Biological Applications (Full Report)

The Modular Autonomous Research System (MARS) was developed to address the pressing need for faster, more reliable, and more adaptable scientific discovery. Traditional experimentation is limited by manual labor, long cycle times, and fragmented data streams, which constrain the ability to explore complex chemical and materials design spaces. To overcome these limitations, we created an integrated, modular platform that combines laboratory robotics, diverse measurement instruments, and a central data infrastructure with artificial intelligence–driven decision-making. The system links liquid handling robots, robotic arms, and optical plate readers into a closed loop where experiments are executed automatically, data is analyzed in real time, and subsequent experimental conditions are adaptively chosen to maximize information gain. Over the course of the project, MARS was validated on two primary test cases—spectroscopic metal–ligand binding assays and peptide-directed mineralization—which highlighted the system’s ability to handle uncertainty and variability in experimental measurements. To further demonstrate modularity and extensibility, we also established additional testbeds in electrochemistry for catalyst discovery and electrolyte formulation for advanced batteries. The results show that MARS can reliably conduct autonomous campaigns with minimal human intervention, adapt to distinct scientific domains, and provide a scalable model for future self-driving laboratories. This work establishes new capabilities for modular, uncertainty-aware automation and directly supports the need for advanced, data-driven research platforms capable of accelerating discovery across a wide range of scientific and national security missions.

59 BASIC BIOLOGICAL SCIENCES

Development and Validation of Smart Building Technology Modules for Academic and Professional Education (Final Technical Report)

Smart building technologies can improve building energy efficiency and resilience, reduce carbon emissions, and provide load flexibility to the grid. However, in both college curricula and building professionals’ continuing education, there is a lack of systematic instruction on smart building technologies. Slipstream, partnering with Texas A&M University (TAMU), the Society of Building Science Educators (SBSE), and the National Institute of Building Sciences (NIBS), developed a semester-long smart building curriculum for college students and 16 training videos for building professionals and the general public. The education and training cover the drivers and benefits of smart building technologies, key building energy systems, the latest sensor technologies and IoT devices, and focus on topics related to smart building controls (i.e., energy management information systems, smart building control platforms, cybersecurity, grid-interactive-efficient buildings [GEBs], smart building control methods, and occupant-centric control). The smart building curriculum for college students was taught at TAMU in the Spring semester of 2024 as part of the validation process. Student feedback was collected and summarized in a validation report by TAMU. The curriculum material was also reviewed by SBSE faculty who are interested in teaching smart building technology-related courses. Suggestions on revisions and better adoption of the materials by other faculty across the architectural, engineering, and construction (AEC) domains were compiled in a distinct validation report by SBSE. The SBSE validation report was used to create structured subsets of the curriculum material for adoption at different levels in different sub-disciplines. These subsets are categorized and offered on the SBSE website (https://www.sbse.org/courses/Smart-Building-Technologies). The 16 training videos for building professionals and the general public were previewed by 17 industry experts, and feedback and suggested changes were incorporated into the final version of these videos. The videos are organized into a smart building technology training course and published on the Whole Building Design Guide website (https://www.wbdg.org/ce/doe/bto/sbtt), which is hosted by the National Institute of Building Sciences (NIBS). Project team members created marketing materials to promote the awareness of these free, publicly available education and training resources. Outreach and marketing activities included creating short promotional videos, building project webpages, making project announcements on social media, conducting an email campaign, and directly reaching out to faculties and building professionals. This report describes the project approach, provides outlines of the training materials, along with links to resources, and identifies lessons learned in creating the content. We also suggest ways to scale the instruction of smart building concepts to empower the workforce to accelerate the adoption of smart building technologies in the real world.

99 GENERAL AND MISCELLANEOUS

SpecSims: A Scalable Speculative Tree-based Simulation Cloning Framework for Finite Memory Machines

Simulation cloning is a technique in which cloned simulations whose state spaces differ partially from their parent simulation due to intervening events are spawned at runtime and concurrently advanced. It is a powerful method to carry out what-if analysis by speculatively exploring and evaluating the impact of various permutations of intervening cascade of events. Due to the exponential growth in the number of possible clones even for a small number of distinct intervening events, the practical efficacy of the approach is often severely limited by the maximum available memory of the computing host. In this paper, we introduce a novel speculative simulation cloning framework that executes a simulation cloning campaign capable of efficiently exploring an exponentially large space of clone simulations created by permutation of intervening events under a finite memory constraint. We provide a theoretical analysis of the runtime characteristics of our proposed approach and highlight its novel advantages such as memory-aware and as-long-as-needed execution. Furthermore, in support of our analytical findings and to demonstrate its practical feasibility, we implement a prototype of the cloning framework on a shared memory system and report its performance characteristics in the context of a heat diffusion simulation, and a power grid simulation subject to cascading disruptions from geomagnetic disturbances.

Simulation framework

5G-TSN Architecture Capable of Providing Real-time Situational Awareness to Fossil-Energy (FE) Generation Systems (Final Technical Report)

This final report highlights the comprehensive achievements of the project focused on developing and validating a 5G-Time Sensitive Networking (TSN) architecture tailored for real-time operational awareness in fossil energy systems. The initiative successfully advanced through a series of technical milestones, including the integration of EMI-aware network models, deployment of advanced simulation frameworks, and real-world performance characterization at key sites such as UTEP and Fabens. Through the strategic use of NetSim® software, the team created and validated network configurations for wired and wireless environments, tested under varying congestion conditions, and verified network slicing implementations for URLLC-specific applications. Major accomplishments include the migration of simulation tools to the latest NetSim® version to support accurate modeling of TSN and network slicing, extensive EMI measurement campaigns, and the development of a robust simulation model for end-to-end SCADA system integration. Simulations compared both TDD and FDD duplexing modes, revealing insights into their performance under congested conditions. The wireless network was benchmarked for throughput, jitter, and delay metrics, aligning with 3GPP Release 15/16 and IEEE 802.1-TSN standards. A peer-reviewed conference paper was accepted and published, contributing to the broader academic and industrial discourse on 5G-TSN integration in energy systems, in addition to a journal article. Despite minor delays due to software limitations, the project achieved its objectives and delivered validated architecture ready for deployment in advanced energy network environments.

01 COAL, LIGNITE, AND PEAT

Measurements and Analyses to Enable Science for the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE)

Coastal cities provide the opportunity to characterize the substantial effects of manmade particles on marine cloud properties and processes. La Jolla lies to the north of San Diego, California, but it is often about a day directly downwind of the major pollution sources located in the ports of Los Angeles and Long Beach. The large dynamic range of aerosol particle concentrations combined with the multi-hour to multi-day persistence of stratocumulus cloud layers makes the site ideal for investigating the seasonal changes in cloud and aerosol properties as well as the quantitative relationships between cloud and aerosol properties. The Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) characterized the extent, radiative properties, aerosol interactions, and precipitation characteristics of stratocumulus clouds in the Eastern Pacific across all four seasons at two coastal sites in La Jolla. This project was designed to enhance and expand the scientific uses of the ARM AMF1 measurements from EPCAPE. The goal was to ensure meaningful observations were collected that would enable science. The project included the following objectives: (1) Reviewing the ARM AMF1 measurements [ARM, 2021a; b] and circulating a summary of the measurements each week, (2) Comparing AMF1 measurements to those provided by collaborators (including filter measurements at the pier), (3) Collecting and analyzing filter samples from Scripps Pier by Fourier Transform Infrared spectroscopy (FTIR) and X-ray Fluorescence (XRF), (4) Assisting in operations of instruments provided by Guest PIs when possible, and (5) Providing an initial compilation of EPCAPE aerosol and cloud seasonal differences. The expected outcomes of these objectives were enhanced proposals and publications using EPCAPE measurements by helping to identify instrument issues, expanded data access and awareness by distributing weekly plots and related summaries, improved source-related attribution of aerosols with elemental tracers, additional observations provided by Guest PIs, and accelerated ACI studies enabled by the compiled seasonal summaries of aerosol and cloud properties. Three examples of the science enabled by this project are findings that (i) aerosol and cloud aqueous production contributes more than half of sulfate particle mass concentration, (ii) upwind sources make chemical composition very similar at nearby sites despite local differences in meteorology, and (iii) most of the large mass concentration of semi-volatile organic components is co-emitted and co-evaporated with nitrate. Together these findings illustrate how ARM extended field campaigns in coastal regions can be used to constrain ACI processes with direct observations. By using the unique ARM suite of cloud radiative products in addition to the measured aerosol properties at a coastal location, we were able to address the more specific question of which aerosol particles cause how much of the effects on clouds. Identifying this signature in coastal areas provides an opportunity to test the representation of aerosol sources by global models in a range of clean and urban-influenced conditions.

Russell, Lynn [Univ. of California, San Diego, CA

IRMA

IRMA (In)elastic Representation of Materials As S(α,β) evaluations IRMA turns one phonon model into three outputs that usually require three separate tool chains: an evaluated nuclear-data file, predicted neutron-scattering spectra, and scattering kernels for Monte Carlo transport. The three outputs draw on a single, consistent description of the material, so the evaluation, the spectroscopy that can validate it, and the transport that uses it always agree about the physics. Nuclear data. IRMA writes ENDF-6 File 7 thermal scattering evaluations on automatically constructed (α, β) grids. This part reimplements and generalizes NJOY's LEAPR: the classic kernels reproduce freshly generated NJOY2016 tapes digit for digit and published reference tapes to about 1e-4, and the generalized paths add the exact coherent one-phonon term, anisotropic Debye-Waller tensors, coherent elastic for arbitrary crystals, and a per-species partition for polyatomic materials. The tapes feed NJOY, AMPX, FUDGE, and every transport code downstream of them. Neutron spectroscopy. The irma.spectra forward model projects the same physics onto an instrument's kinematics and resolution: INS spectra for VISION and generic indirect geometries, and 2-D S(Q,E) powder maps for direct-geometry spectrometers, from a phonopy model or straight from a phonon DOS. It can be used to predict a proposed measurement before beam time; in analysis, it supplies the calculated single-scattering counterpart of a measured spectrum, from the same material description the evaluation was built from. Monte Carlo transport. The irma.ncrystal exporter writes per-temperature scattering kernels for the companion NCrystal plugin, so McStas, OpenMC, and other NCrystal-aware codes sample the same physics. The exported kernels carry the per-site anisotropic Debye-Waller tensors, keeping directional coherent-elastic physics that NCrystal's standard scalar treatment does not represent. With the same physics inside a transport code, an entire beamline becomes a virtual experiment: IRMA's end-to-end validation ran a custom McStas implementation of the ARCS spectrometer, assembled from the existing McVine and McStas models, against measured data. From a bare crystal structure. The irma mlip front end builds the phonon model itself: a structure file and a choice of potential are enough. Nine pretrained machine-learned interatomic potentials are supported, on a laptop CPU, with no first-principles calculation; an approximate phonon model for a new material costs minutes, not a DFT campaign, and the build emits prefilled inputs for all three outputs. The result is a good starting point rather than a finished evaluation: survey-quality physics with every parameter exposed for review. A converged atomistic calculation enters the same way, as a phonopy model, when higher fidelity is needed.

Ramic, Kemal [Oak Ridge National Laboratory (ORNL)

Impacts of Synoptic-Scale Dynamics on Clouds and Radiation in High Southern Latitudes

High-latitudinal mixed-phase clouds significantly affect Earth's radiative balance. Observations of cloud and radiative properties from two field campaigns in the Southern Ocean and Antarctica were compared with two global climate model simulations. A cyclone compositing method was used to quantify “dynamics-cloud-radiation” relationships relative to the extratropical cyclone centers. Observations show larger asymmetry in cloud and radiative properties between western and eastern sectors at McMurdo compared with Macquarie Island. Most observed quantities at McMurdo are higher in the western (i.e., post-frontal) than the eastern (frontal) sector, including cloud fraction, liquid water path (LWP), net surface shortwave and longwave radiation (SW and LW), except for ice water path (IWP) being higher in the eastern sector. Here, the two models were found to overestimate cloud fraction and LWP at Macquarie Island but underestimate them at McMurdo Station. IWP is consistently underestimated at both locations, both sectors, and in all seasons. Biases of cloud fraction, LWP, and IWP are negatively correlated with SW biases and positively correlated with LW biases. The persistent negative IWP biases may have become one of the leading causes of radiative biases over the high southern latitudes, after correcting the underestimation of supercooled liquid water in the older model versions. By examining multi-scale factors from cloud microphysics to synoptic dynamics, this work will help increase the fidelity of climate simulations in this remote region.

58 GEOSCIENCES

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