Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bulk properties”

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

Learning interpretable surface elasticity properties from bulk properties via neural network equation learners

Surface elasticity is central to understanding the mechanics and stability of surfaces and interfaces. It is characterized by quantities such as surface tension, residual surface stress, and surface stiffness. However their analytical expressions are typically difficult to derive from atomistic data, and depend strongly on modeling choices. This work presents a neural network-based equation learner which combines customized activation functions and connection-based pruning to discover parsimonious, closed-form equations for surface elasticity from atomistic simulations. Applying the method to seven face-centered cubic (FCC) metals, our equation learner uncovers interpretable equations that describe both low-Miller index and high-Miller index surface properties, capturing long-tail property distributions accurately. The discovered expressions are decoupled into two components: a universal, geometry-driven orientation function, and material-specific baseline coefficients. We find that lower-order properties such as surface tension are fundamentally geometry dependent, while higher-order properties such as surface stress and elasticity show more complex geometry and material dependence. We also relate material dependent coefficients to bulk properties, forming a clear map from bulk material properties to surface elasticity. Overall, this approach demonstrates that interpretable neurosymbolic machine learning can bridge the gap between atomistic simulations and physical laws, enabling the discovery of generalizable structure–property relationships for materials science phenomena such as surface elasticity.

Equation learning↗
DOE OSTI2025

Dynamics of hot QCD matter 2024 — Bulk properties

The second Hot QCD Matter 2024 conference at IIT Mandi focused on various ongoing topics in high-energy heavy-ion collisions, encompassing theoretical and experimental perspectives. This proceedings volume includes 19 contributions that collectively explore diverse aspects of the bulk properties of hot QCD matter. The topics encompass the dynamics of electromagnetic fields, transport properties, hadronic matter, spin hydrodynamics, and the role of conserved charges in high-energy environments. These studies significantly enhance our understanding of the complex dynamics of hot QCD matter, the quark–gluon plasma (QGP) formed in high-energy nuclear collisions. Advances in theoretical frameworks, including hydrodynamics, spin dynamics and fluctuation studies, aim to improve theoretical calculations and refine our knowledge of the thermodynamic properties of strongly interacting matter. Experimental efforts, such as those conducted by the ALICE and STAR collaborations, play a vital role in validating these theoretical predictions and deepening our insight into the QCD phase diagram, collectivity in small systems, and the early-stage behavior of strongly interacting matter. In conclusion, combining theoretical models with experimental observations offers a comprehensive understanding of the extreme conditions encountered in relativistic heavy-ion and proton-proton collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Mesoscale Linear Elastic Modeling and Homogenization of Marine Energy Composites

The design of fiber-reinforced composite (FRC)-based components for marine energy applications necessitates a fundamental understanding of material properties and the resulting geometry to predict long-term performance. In this work, we present a modeling workflow to predict linear elastic and diffusive bulk properties at the mesoscale for an idealized geometry based on knowledge of fiber and resin properties. A parametric study was performed to identify the key model input parameters that influence bulk properties. Furthermore, we demonstrate how bulk properties can be leveraged in high-fidelity image-based simulations, where imperfections in tow geometry and voids captured during X-ray computed tomography imaging are explicitly represented within the simulation. Bulk properties of interest include moduli, Poisson’s ratios, hygroscopic swelling, diffusivity, and moisture uptake, which are key parameters for characterizing FRC performance within marine environments. Modeling predictions agreed well with experimental data, except for estimating swelling coefficients, likely due to crack accumulation as a function of moisture uptake. The mesoscale modeling workflow ultimately highlights a versatile framework for understanding the influence of material and geometric properties, which can be leveraged to rapidly assess new FRC-based components.

computational mechanics↗
DOE OSTI2024

Multiscale Modeling of Silicon Carbide Cladding for Nuclear Applications: Thermal Performance Modeling

The complex multiscale and anisotropic nature of silicon carbide (SiC) ceramic matrix composite (CMC) makes it difficult to accurately model its performance in nuclear applications. The existing models for nuclear grade composite SiC do not account for the microstructural features and how these features can affect the thermal and structural behavior of the cladding and its anisotropic properties. In addition to the microstructural features, the properties of individual constituents of the composites and fiber tow architecture determine the bulk properties. Models for determining the relationship between the individual constituents’ properties and the bulk properties of SiC composites for nuclear applications are absent, although empirical relationships exist in the literature. Here, a hierarchical multiscale modeling approach was presented to address this challenge. This modular approach addressed this difficulty by dividing the various aspects of the composite material into separate models at different length scales, with the evaluated property from the lower-length-scale model serving as an input to the higher-length-scale model. The multiscale model considered the properties of various individual constituents of the composite material (fiber, matrix, and interphase), the porosity in the matrix, the fiber volume fraction, the composite architecture, the tow thickness, etc. By considering inhomogeneous and anisotropic contributions intrinsically, our bottom-up multiscale modeling strategy is naturally physics-informed, bridging constitutive law from micromechanics to meso-mechanics and structural mechanics. The effects that these various physical attributes and thermo-physical properties have on the composite’s bulk thermal properties were easily evaluated and demonstrated through the various analyses presented herein. Since silicon carbide fiber-reinforced SiC CMCs are also promising thermal–structural materials with a broad range of high-end technology applications beyond nuclear applications, we envision that the multiscale modeling method we present here may prove helpful in future efforts to develop and construct reinforced CMCs and other advanced composite nuclear materials, such as MAX phase materials, that can service under harsh environments of ultrahigh temperatures, oxidation, corrosion, and/or irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Effect of metal encapsulation on bulk superconducting properties of niobium thin films used in qubits

Niobium metal occupies nearly 100% of the volume of a typical 2D transmon device. While the aluminum Josephson junction is of utmost importance, maintaining quantum coherence across the entire device means that pair-breaking in Nb leads, capacitive pads, and readout resonators can be a major source of decoherence. The established contributors are surface oxides and hydroxides, as well as absorbed hydrogen and oxygen. Metal encapsulation of freshly grown surfaces with non-oxidizing metals, preferably without breaking the vacuum, is a successful strategy to mitigate these issues. While the positive effects of encapsulation are undeniable, it is important to understand its impact on the macroscopic behavior of niobium films. We present a comprehensive study of the bulk superconducting properties of Nb thin films encapsulated with gold and palladium/gold, and compare them to those of bare Nb films. Magneto-optical imaging, magnetization, resistivity, and London and Campbell penetration depth measurements reveal significant differences in encapsulated samples. Both sputtered, and epitaxial Au-capped films exhibit the highest residual resistivity ratio and superconducting transition temperature, as well as the lowest upper critical field, London penetration depth, and critical current. These results are in good agreement with the microscopic theory of anisotropic normal and superconducting states of Nb. We conclude that pair-breaking in the bulk of niobium films, driven by disorder throughout the film rather than just at the surface, is a significant source of quantum decoherence in transmons. We also conclude that gold capping not only passivates the surface but also affects the properties of the entire film, significantly reducing the scattering rate due to defects likely induced by surface diffusion if the film is not protected immediately after fabrication.

Datta, Amlan [Ames Lab; Iowa State U.]↗
DOE OSTI2025

Creation of an Fe 3 P Schreibersite Density Functional Tight Binding Model for Astrobiological Simulations

The mineral schreibersite, e.g., Fe 3 P, is commonly found in iron-rich meteorites and could have served as an abiotic phosphorus source for prebiotic chemistry. However, atomistic calculations of its degradation chemistry generally require quantum simulation approaches, which can be too computationally cumbersome to study sufficient time and length scales for this process. In this regard, we have created a computationally efficient semiempirical quantum density functional tight binding (DFTB) model for iron and phosphorus-containing materials by adopting an existing semiautomated workflow that represents many-body interactions by linear combinations of Chebyshev polynomials. Here, we have utilized a relatively small training set to optimize a DFTB model that is accurate for schreibersite physical and chemical properties, including its bulk properties, surface energies, and water absorption. We then show that our model shows strong transferability to several iron phosphide solids as well as multiple allotropes of iron metal. Our resulting DFTB parametrization will allow us to interrogate schreibersite aqueous decomposition at longer time and length scales than standard quantum approaches, providing for more detailed investigations of its role in prebiotic chemistry on early Earth.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Synthesis and Properties of Bulk Mg 3 WN 4 in a Wurtzite-Derived Structure

Experimental synthesis of theoretically predicted materials with controlled elemental coordination environments can lead to the realization of useful properties, such as facile ion transport or ferroelectric switching. Among such materials are ternary nitrides in the Mg−W−N composition space, where several stable and metastable compounds have been predicted and synthesized recently in bulk and film forms. Here, we report for the first time on the bulk synthesis of Mg 3 WN 4 in a wurtzite-derived crystal structure via a solid-state metathesis reaction. In situ synchrotron powder X-ray diffraction shows how the ion exchange proceeds from Li 6 WN 4 + MgCl 2 precursors to Mg 3 WN 4 + LiCl products, with the reaction starting slowly near 380 °C and completing by 600 °C, including the presence of a competing disordered rocksalt-derived phase (Mg,W)N above 440 °C. The follow-up ex situ powder synthesis at 400 °C for 0.5 h with 10% excess MgCl 2 reveals the cation-ordered nature of the wurtzite-derived Mg 3 WN 4 structure with polar symmetry confirmed by second-harmonic generation measurements. Optical absorption spectra, chemical composition analysis, and electron microscopy imaging suggest that bulk Mg 3 WN 4 obtained via metathesis reaction is prone to defect formation. Overall, this study shows that selective ex situ synthesis of the phase-pure ternary nitrides, informed by in situ measurements, is possible by carefully controlling the thermal budget of the reaction and paves the way toward property characterization of the Mg 3 WN 4 wurtzite-derived material.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthesis and properties of bulk cobalt-carbide permanent magnets

Nanocrystalline bulk Co-C permanent magnets were produced by rapid hot consolidation of ball milled powder precursors. Bulk magnets with an average grain size of 20 nm and densities up to 94 percent of the theoretical density were successfully processed at pressures exceeding 1 GPa for the first time. Their microstructural and magnetic properties along with their magnetocrystalline anisotropy values are reported. Presence of iron contamination that results from the use steel milling media in milling was confirmed by elemental mapping of the bulk magnets. It is concluded that iron bonds in the form of hematite and does not adversely affect coercivity values. Here, this study shows that the Co 3 C compound does possess a sizable magnetocrystalline anisotropy, making it a suitable candidate as a rare earth free permanent magnet.

Co3 C↗
DOE OSTI2025

Characterizing the Oscillatory Properties of Bulk Electric Systems

This paper presents a process for characterizing the oscillatory dynamics of a large bulk power system. As a demonstration, the process is applied to the Western Interconnection of North America. Several complementary analysis approaches, both new and existing, are employed to provide a comprehensive understanding of the oscillatory properties of the system. Established modal analysis techniques based on ringdown and mode-meter algorithms are utilized. In addition, we derive and apply methods based on spectral correlation analysis to identify modal frequencies, distinguish between modes that are closely spaced in frequency, and determine locations at which the modes are observable. Critical interarea modes are identified and characterized using actual-system synchrophasor measurements taken over several years of operation in concert with industry-standard simulation models. This includes 145 hours of PMU data and two planning base cases.

42 ENGINEERING↗
DOE OSTI2025

Electrocatalytic grafting of polyvinyl chloride plastics

Polyvinyl chloride (PVC) plastics require high loadings of plasticizers and stabilizers to achieve commercially useful bulk properties. However, these non-covalent additives leach from PVC over time, resulting in the loss of their tailored functionality. This work details the electrocatalytic functionalization of PVC to covalently graft plasticizing additives directly onto the polymer backbone. Here, mechanistic insights guided the design of electrocatalysts capable of modifying C–Cl bonds of PVC under mild conditions with high selectivity while suppressing side reactions such as elimination and chain scission. Functional groups that mimic PVC plasticizers are covalently installed into the backbone of PVC to create new materials with distinct bulk properties from the original polymer. The degree of polymer grafting is easily controlled by simply changing the redox capacity that is passed during electrolysis. This strategy is employed to create chemically- and leach-resistant PVC materials by directly electrolyzing mixtures of consumer PVC products

Zackasee, Jordan L.S.↗
DOE OSTI2025

Probing Alkylpolyglycoside Hydrogen Bonding and Its Destabilization by Sulfonate Hydrotropes with X‑ray and Vibrational Spectroscopy

It has been proposed that intersurfactant H-bonding networks can produce highly stable foams, specifically those generated from solutions containing alkylpolyglycoside surfactants. In this work, we aim to characterize the presence and destabilization of these networks by introducing a hydrotrope (sodium p-toluene sulfate) at various concentrations into an alkylpolyglycoside (Glucopon 225 DK) surfactant containing solution. Solution surface properties are probed by using aerosol velocity map imaging X-ray photoelectron spectroscopy (A-VMI-XPS) and tensiometer measurements. Bulk properties below the surface are probed using Fourier transform infrared spectroscopy (FTIR) and C-edge near-edge X-ray fine structure spectroscopy (NEXAFS). The surface measurements provide the framework to describe the destabilization of the solution H-bonding network, while the bulk solution measurements provide hints about the disruption of the hydrogen bonding network upon hydrotrope addition. The collected data support the hypothesis that the destabilization of the intersurfactant H-bonding network frees surfactant molecules from the bulk, increasing surfactant population at the air/water interface. This was quantified through an increase in peak area and width in XPS measurements as well as a decrease in surfactant critical micelle concentration. Three regimes with increasing amounts of hydrotrope addition are suggested, described as (1) a hydrotrope affecting only surfactant surface properties, (2) a hydrotrope affecting surfactant surface and bulk properties, and (3) hydrotrope-dominated surface and bulk properties. Future studies will characterize foam stability across the hydrotrope concentration regimes to better define correlations between intersurfactant H-bonding networks and foam stability.

Molecules↗
DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
DOE OSTI2025

Amorphous complexion-aided sintering enables scalable processing of bulk nanocrystalline Cu-Zr with high strength and compressive plasticity

Nanocrystalline alloys can have exceptional strengths, yet due to limited microstructural stability it is difficult to fabricate bulk pieces through traditional processing routes that retain nanosized grains. In this study, centimeter-sized Cu-Zr alloy pellets were fabricated via a simple and improved powder metallurgy processing route. Different consolidation temperatures and times were employed to investigate the effect of amorphous grain boundary complexions on densification and the resulting mechanical properties. Bulk compression tests were carried out, with the samples that were hot pressed at 900 °C for 10 h exhibiting an excellent combination of average yield strength of 722 ± 45 MPa and average failure strain of 25.3 ± 2.4 %. Furthermore, we find that a powder processing route which enables amorphous complexion-assisted sintering leads to samples that (1) reach full density without requiring quenching treatments or other complex processing, (2) demonstrate appreciable plasticity, and (3) have strength that competes with commercially available high-strength Cu alloys.

Bulk nanocrystalline alloy↗
DOE OSTI2026

Insulating ground state and 2−𝑘 magnetic structure of candidate Weyl hydrogen-atom K2⁢Mn3⁢(AsO4)3

The ideal Weyl "hydrogen-atom" semimetal exhibits only a single pair of Weyl nodes and no other trivial states at the Fermi energy. Such a material would be a panacea in the study of Weyl quasi-particles, allowing direct unambiguous observation of their topological properties. The alluaudite-like K2⁢Mn3⁢(AsO4)3 compound was recently proposed as such a material. Here, we use comprehensive experimental work and first-principle calculations to assess this prediction. We find K2⁢Mn3⁢(AsO4)3 crystallizes in the 𝐶⁢2/𝑐 symmetry with a quasi-one-dimensional Mn sublattice, growing as small needle-like crystals. Bulk property measurements reveal magnetic transitions at ≈8 and ≈4 K, which neutron scattering experiments show correspond to two distinct magnetic orders, first a partially ordered ferrimagnetic 𝐤𝟏=(0,0,0) structure at 8 K and a second transition of 𝐤𝟐=(1,0,0) at 4 K to a fully ordered state. Below the second transition, both ordering vectors are necessary to describe the complex magnetic structure with modulated spin magnitudes. Both of the best-fit magnetic structures in this work are found to break the symmetry necessary for the generation of Weyl nodes, though one of the magnetic structures allowed by 𝐤𝟏 does preserve this symmetry. However, the crystals are optically transparent and ellipsometry measurements reveal a large band gap, undermining expectations of semimetallic behavior. Density functional theory calculations predict an insulating antiferromagnetic ground state, in contrast to previous reports, and suggest potential frustration on the magnetic sublattice. Given the wide tunability of the alluaudite structure, we consider ways to push the system closer to a semimetallic state.

Taddei, Keith M [Argonne National Laboratory]↗
DOE OSTI2025

Across the Scales of the Nucleus: Understanding Short Range Correlations from Medium Modification to Probe Independence

The atomic nucleus presents an intricate system due to the non-linear forces described by Quantum Chromodynamics (QCD) that govern its structure. The range of scales involved is remarkable; the most massive nuclei weigh approximately five orders of magnitude more than the quarks that compose them. The nucleus can be analyzed at various levels, from quarks to hadrons to the nucleus as a whole. Short-Range Correlations (SRCs) within the nucleus play a significant role that spans these diverse scales. At the most fundamental level, SRCs influence the interaction between nucleons. The nucleon-nucleon (NN) interaction, arising from QCD, is crucial in determining nuclear properties. SRCs serve as valuable probes for measuring this NN interaction, as the nucleons within SRCs become effectively decoupled from the rest of the nucleus. Multiple experimental techniques, including electron scattering, have been employed to investigate the NN interaction through SRCs. However, our first project demonstrates that inclusive measurements alone are inadequate to constrain this interaction fully. Moving to the scale of the nucleus, SRCs contribute to the high-momentum tail of the nuclear spectral function. While the low-momentum region is characterized by nucleons exhibiting bulk properties, nucleons begin to pair into SRCs at higher momenta. Our research aims to bridge the understanding between the mean-field portion of the nucleus and its high-momentum SRC components. Additionally, SRCs affect the quark structure of protons, as evidenced by the EMC effect, which indicates that quarks behave differently when protons are embedded within a nucleus—an effect referred to as medium modification. This thesis explores the correlation between SRCs and medium modification across various experimental setups. Finally, we seek to establish an interpretation of the nuclear ground state. Accomplishing this requires demonstrating that our SRC observables are independent of the probe’s scale and scheme. The concluding project of this thesis illustrates how we utilize triple coincidence quasi-elastic scattering across a range of (Q2) values to develop a model-dependent framework for understanding SRC distributions within the nucleus’s ground-state wavefunction.

Denniston, A. W. Denniston [University of Tel Aviv↗
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↗