Search NASA⌕ Search

SEARCH · Search NASA

Results for “density estimation”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Analyzing LF/VLF Lightning Waveforms to Estimate D-region Electron Density Profiles

Lightning waveforms in the low frequency (LF; 30-300 kHz) and the very low frequency (VLF; 3-30 kHz) bands can be exploited to produce data-driven ionospheric D-region electron density profile (EDP) estimates with significantly higher spatial and temporal coverage than previously available. The lightning waveforms used in this paper are signals detected in the LF/VLF of negative cloud-to-ground lightning by the Earth Networks Total Lightning Detection Network. Each waveform contains a ground wave and a time-delayed ionospheric reflection. The time delay between the ground wave and ionospheric reflection has previously been used to estimate a single specular reflection altitude, where LF/VLF emissions are reflected by the ionosphere. Here, we expand upon previous methods to include filtering and spectral analysis, and account for oblique propagation to produce higher-order estimates for reflection altitudes and corresponding electron densities. Once estimated, reflection altitudes and corresponding electron densities can be used to derive parameters β and h’, which define an EDP for the D-region. In this study, the lightning waveform (LW) analysis is demonstrated using a single representative 24-hour dataset over the Southeast United States, and then extended to a total of 10 separate datasets with varying locations and ionospheric conditions. The LW-derived D-region EDPs are in agreement with predictions made by the Faraday International Reference Ionosphere model, and the LW EDPs β and h’ values are consistent with previous LF/VLF-derived estimates.

D-region ionosphere↗
DOE OSTI2013

Atmospheric pressure He-air plasma jet: Breakdown process and propagation phenomenon

In this paper He-discharge (plasma jet/bullet) in atmospheric pressure air and its progression phenomenon has been studied experimentally using ICCD camera, optical emission spectroscopy (OES) and calibrated dielectric probe measurements. The repetitive nanosecond pulse has applied to a plasma pencil to generate discharge in the helium gas channel. The discharge propagation speed was measured from the ICCD images. The axial electric field distribution in the plasma jet is inferred from the optical emission spectroscopic data and from the probe measurement. The correlation between the jet velocities, jet length with the pulse duration is established. It shows that the plasma jet is not isolated from the input voltage along its propagation path. The discharge propagation speed, the electron density and the local and average electric field distribution along the plasma jet axis predicted from the experimental results are in good agreement with the data predicted by numerical simulation of the streamer propagation presented in different literatures. The ionization phenomenon of the discharge predicts the key ionization parameters, such as speed, peak electric field in the front, and electron density. The maximum local electric field measured by OES is 95 kV/cm at 1.3 cm of the jet axis, and average EF measured by probe is 24 kV/cm at the same place of the jet. The average and local electron density estimated are in the order of 10{sup 11} cm{sup -3} and it reaches to the maximum of 10{sup 12} cm{sup -3}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Stochastic multiscale modeling for quantifying statistical and model errors with application to composite materials

This paper provides a coherent and efficient computational framework for stochastic multiscale analysis of material systems in the presence of parametric uncertainties and modeling errors. Uncertainty in those model parameters that are not deduced as upscaled quantities is attributed to an uncertainty “germ”. While such parameters can appear at any scale, they are predominant at the finest analysis scale. Additional uncertainties stemming from statistical estimation, attributed to lack of data and model error, are associated with each submodel contributing to the multiscale system. Here, a robust and efficient framework based on a generalized extended polynomial chaos expansion (gEPCE) is proposed to simultaneously propagate all these uncertainties in order to provide a probabilistic representation of specific quantities of interest (QoI). We characterize the full probability distribution of the QoI and the uncertainty in the failure probability pertaining to its tails. By combining gEPCE with kernel density estimation (KDE) and directional derivatives, we construct sensitivity measures that connect these statistical metrics of QoI to the various sources of uncertainty to assess their individual and combined impacts. An illustrative problem featuring three-point bending of a composite beam is investigated to demonstrate the presented approach.

36 MATERIALS SCIENCE↗
NASA NTRS1986

Voyager 2 radio science observations of the Uranian system Atmosphere, rings, and satellites

The results of preliminary analyses of radio occultation data obtained by Voyager 2 as it passed Uranus are described. The occultations took place between 2-7 deg S latitude and yielded atmospheric temperature profiles between pressure levels of 10-900 mbar, an altitude range of 100 km. The mole fractions of hydrogen and helium in the tropopause were estimated, in conjunction with IR data, to about 0.85 and 0.10-0.20, respectively. Radio signal intensity data indicated the presence of a cloud deck of CH4 ice at a pressure level of 1300 mbar and a temperature of 81 K, implying a CH4 mole fraction of 0.02 at very low altitudes. The ionosphere extended upward, in two levels, to more than 10,000 km altitude. The ring system was different than the one around Saturn and possessed cylindrical substructures. The radio data also permitted mass density estimates for the five major moons, i.e., about 1.40 gr/cu cm, a value which rules out cometary origins.

Tyler, G. L.↗
NASA NTRS2010

Mars-GRAM: Increasing the Precision of Sensitivity Studies at Large Optical Depths

The Mars Global Reference Atmospheric Model (Mars-GRAM) is an engineering-level atmospheric model widely used for diverse mission applications. Mars-GRAM's perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL). It has been discovered during the Mars Science Laboratory (MSL) site selection process that Mars-GRAM, when used for sensitivity studies for MapYear=0 and large optical depth values such as tau=3, is less than realistic. A comparison study between Mars atmospheric density estimates from Mars-GRAM and measurements by Mars Global Surveyor (MGS) has been undertaken for locations of varying latitudes, Ls, and LTST on Mars. The preliminary results from this study have validated the Thermal Emission Spectrometer (TES) limb data. From the surface to 80 km altitude, Mars-GRAM is based on the NASA Ames Mars General Circulation Model (MGCM). MGCM results that were used for Mars-GRAM with MapYear=0 were from a MGCM run with a fixed value of tau=3 for the entire year at all locations. This has resulted in an imprecise atmospheric density at all altitudes. To solve this pressure-density problem, density factor values were determined for tau=.3, 1 and 3 that will adjust the input values of MGCM MapYear 0 pressure and density to achieve a better match of Mars-GRAM MapYear 0 with TES observations for MapYears 1 and 2 at comparable dust loading. The addition of these density factors to Mars-GRAM will improve the results of the sensitivity studies done for large optical depths.

Justh, Hilary L.↗
DOE OSTI2010

SEEN AND UNSEEN TIDAL CAUSTICS IN THE ANDROMEDA GALAXY

Indirect detection of high-energy particles from dark matter interactions is a promising avenue for learning more about dark matter, but is hampered by the frequent coincidence of high-energy astrophysical sources of such particles with putative high-density regions of dark matter. We calculate the boost factor and gamma-ray flux from dark matter associated with two shell-like caustics of luminous tidal debris recently discovered around the Andromeda galaxy, under the assumption that dark matter is its own supersymmetric antiparticle. These shell features could be a good candidate for indirect detection of dark matter via gamma rays because they are located far from the primary confusion sources at the galaxy's center, and because the shapes of the shells indicate that most of the mass has piled up near the apocenter. Using a numerical estimator specifically calibrated to estimate densities in N-body representations with sharp features and a previously determined N-body model of the shells, we find that the largest boost factors do occur in the shells but are only a few percent. We also find that the gamma-ray flux is an order of magnitude too low to be detected with Fermi for likely dark matter parameters, and about two orders of magnitude less than the signal that would have come from the dwarf galaxy that produces the shells in the N-body model. We further show that the radial density profiles and relative radial spacing of the shells, in either dark or luminous matter, is relatively insensitive to the details of the potential of the host galaxy but depends in a predictable way on the velocity dispersion of the progenitor galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

UNVEILING THE NATURE OF THE UNIDENTIFIED GAMMA-RAY SOURCES. IV. THE SWIFT CATALOG OF POTENTIAL X-RAY COUNTERPARTS

A significant fraction (∼30%) of the high-energy γ-ray sources listed in the second Fermi Large Area Telescope (LAT) catalog are still of unknown origin, having not yet been associated with counterparts at lower energies. To investigate the nature of these enigmatic sources, we present an extensive search of X-ray sources lying in the positional uncertainty region of a selected sample of these unidentified gamma-ray sources (UGSs) that makes use of all available observations performed by the Swift X-ray Telescope before 2013 March 31, available for 205 UGSs. To detect the fainter sources, we merged all the observations covering the Fermi LAT positional uncertainty region at a 95% level of confidence of each UGS. This yields a catalog of 357 X-ray sources, finding candidate X-ray counterparts for ∼70% of the selected sample. In particular, 25% of the UGSs feature a single X-ray source within their positional uncertainty region, while 45% have multiple X-ray sources. For each X-ray source, we also looked in the corresponding Swift UVOT merged images for optical and ultraviolet counterparts, also performing source photometry. We found ultraviolet-optical correspondences for ∼70% of the X-ray sources. We searched several major radio, infrared, optical, and ultraviolet surveys for possible counterparts within the positional error of the sources in the X-ray catalog to obtain additional information on their nature. Applying the kernel density estimation technique to infrared colors of Wide-Field Infrared Survey Explorer counterparts of our X-ray sources we select six γ-ray blazar candidates. In addition, comparing our results with previous analyses, we select 11 additional γ-ray blazar candidates.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

Operational Data Reduction Procedure for Determining Density and Vertical Structure of the Martian Upper Atmosphere from Mars Global Surveyor Accelerometer Measurements

The success of aerobraking by the Mars Global Surveyor (MGS) spacecraft was partly due to the analysis of MGS accelerometer data. Accelerometer data was used to determine the effect of the atmosphere on each orbit, to characterize the nature of the atmosphere, and to predict the atmosphere for future orbits. To interpret the accelerometer data, a data reduction procedure was developed to produce density estimations utilizing inputs from the spacecraft, the Navigation Team, and pre-mission aerothermodynamic studies. This data reduction procedure was based on the calculation of aerodynamic forces from the accelerometer data by considering acceleration due to gravity gradient, solar pressure, angular motion of the MGS, instrument bias, thruster activity, and a vibration component due to the motion of the damaged solar array. Methods were developed to calculate all of the acceleration components including a 4 degree of freedom dynamics model used to gain a greater understanding of the damaged solar array. The total error inherent to the data reduction procedure was calculated as a function of altitude and density considering contributions from ephemeris errors, errors in force coefficient, and instrument errors due to bias and digitization. Comparing the results from this procedure to the data of other MGS Teams has demonstrated that this procedure can quickly and accurately describe the density and vertical structure of the Martian upper atmosphere.

Cancro, George J.↗
DOE OSTI2024

Machine learning materials properties with accurate predictions, uncertainty estimates, domain guidance, and persistent online accessibility

One compelling vision of the future of materials discovery and design involves the use of machine learning (ML) models to predict materials properties and then rapidly find materials tailored for specific applications. However, realizing this vision requires both providing detailed uncertainty quantification (model prediction errors and domain of applicability) and making models readily usable. At present, it is common practice in the community to assess ML model performance only in terms of prediction accuracy (e.g. mean absolute error), while neglecting detailed uncertainty quantification and robust model accessibility and usability. Here, we demonstrate a practical method for realizing both uncertainty and accessibility features with a large set of models. We develop random forest ML models for 33 materials properties spanning an array of data sources (computational and experimental) and property types (electrical, mechanical, thermodynamic, etc). All models have calibrated ensemble error bars to quantify prediction uncertainty and domain of applicability guidance enabled by kernel-density-estimate-based feature distance measures. All data and models are publicly hosted on the Garden-AI infrastructure, which provides an easy-to-use, persistent interface for model dissemination that permits models to be invoked with only a few lines of Python code. We demonstrate the power of this approach by using our models to conduct a fully ML-based materials discovery exercise to search for new stable, highly active perovskite oxide catalyst materials.

domain of applicability↗
DOE OSTI2011

Assessment of the Species Composition, Densities, and Distribution of Native Freshwater Mussels along the Benton County Shoreline of the Hanford Reach, Columbia River, 2004

The Hanford Reach of the Columbia River is the last unimpounded section of the river and contains substrate characteristics (cobble, gravel, sand/silt) suitable for many of the native freshwater mussels known to exist in the Pacific Northwest. Information concerning the native mussel species composition, densities, and distributions in the mainstem of the Columbia River is limited. Under funding from the U.S. Department of Energy Richland Operations Office (DOE-RL), Pacific Northwest National Laboratory conducted an assessment of the near-shore habitat on the Hanford Reach. Surveys conducted in 2004 as part of the Ecological Monitoring and Compliance project documented several species of native mussels inhabiting the near-shore habitat of the Hanford Reach. Findings reported here may be useful to resource biologists, ecologists, and DOE-RL to determine possible negative impacts to native mussels from ongoing near-shore remediation activities associated with Hanford Site cleanup. The objective of this study was to provide an initial assessment of the species composition, densities, and distribution of the freshwater mussels (Margaritiferidae and Unionidae families) that exist in the Hanford Reach. Researchers observed and measured 201 live native mussel specimens. Mussel density estimated from these surveys is summarized in this report with respect to near-shore habitat characteristics including substrate size, substrate embeddedness, relative abundance of aquatic vegetation, and large-scale geomorphic/hydrologic characteristics of the Hanford Reach.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1998

Tidal Disruption of Strengthless Rubble Piles: A Dimensional Analysis

A relatively simple prescription for estimating the number of debris clumps (n) that form after a catastrophic tidal disruption event is presented. Following the breakup event, it is assumed that the individual debris particles follow keplerian orbits about the planet until the debris' gravitational contraction timescale (t(sub c)) becomes shorter than its orbital spreading timescale (t(sub s)). When the two timescales become comparable, self-gravity breaks up the debris train into n = L/D clumps, which is the debris length/diameter ratio at that instant. The clumps subsequently orbit the planet independent of each other. The predicted number of clumps n is in good agreement with more sophisticated N-body treatments of tidal breakup for parabolic encounters, and the dependence of n upon the progenitor's density as well as its orbit is also mapped out for hyperbolic encounters. These findings may be used to further constrain both the orbits and densities of the tidally disrupted bodies that struck Callisto and Ganymede. A cursory analysis shows that the Gomul and Gipul crater chains on Callisto, which have the greatest number of craters among the known chains, were formed by projectiles having comet-like densities estimated at rho(sub o) < 1 gm/cc.

Hahn, Joseph M.↗
DOE OSTI2015

Simulations of defense strategies for Bennu: Material characterization and impulse delivery

Assessments of asteroid deflection strategies depend on material characterization to reduce the uncertainty in predictions of the deflection velocity resulting from impulsive loading. In addition to strength, equation of state, the initial state of the material including its competency (i.e. fractured or monolithic) and the amount of micro- or macroscopic porosity are important considerations to predict the thermomechanical response. There is recent interest in observing near-Earth asteroid (101955) Bennu due to its classification of being potentially hazardous with close approaches occurring every 6 years. Bennu is relatively large with a nominal diameter of 492 m, density estimates ranging from 0.9-1.26 g/cm³ and is composed mainly of carbonaceous chondrite. There is a lack of data for highly porous carbonaceous chondrite at very large pressures and temperatures. In the absence of the specific material composition and state (e.g. layering, porosity as a function of depth) on Bennu we introduce a continuum constitutive model based on the response of granular materials and provide impact and standoff explosion simulations to investigate the response of highly porous materials to these types of impulsive loading scenarios. Simulations with impact speeds of 5 km/s show that the shock wave emanating from the impact site is highly dispersive and that a 10% porous material has a larger compacted volume compared with a 40% porous material with the same bulk density due to differences in compaction response.

36 MATERIALS SCIENCE↗
DOE OSTI2015

ASTEROID LIGHT CURVES FROM THE PALOMAR TRANSIENT FACTORY SURVEY: ROTATION PERIODS AND PHASE FUNCTIONS FROM SPARSE PHOTOMETRY

We fit 54,296 sparsely sampled asteroid light curves in the Palomar Transient Factory survey to a combined rotation plus phase-function model. Each light curve consists of 20 or more observations acquired in a single opposition. Using 805 asteroids in our sample that have reference periods in the literature, we find that the reliability of our fitted periods is a complicated function of the period, amplitude, apparent magnitude, and other light-curve attributes. Using the 805-asteroid ground-truth sample, we train an automated classifier to estimate (along with manual inspection) the validity of the remaining ∼53,000 fitted periods. By this method we find that 9033 of our light curves (of ∼8300 unique asteroids) have “reliable” periods. Subsequent consideration of asteroids with multiple light-curve fits indicates a 4% contamination in these “reliable” periods. For 3902 light curves with sufficient phase-angle coverage and either a reliable fit period or low amplitude, we examine the distribution of several phase-function parameters, none of which are bimodal though all correlate with the bond albedo and with visible-band colors. Comparing the theoretical maximal spin rate of a fluid body with our amplitude versus spin-rate distribution suggests that, if held together only by self-gravity, most asteroids are in general less dense than ∼2 g cm{sup −3}, while C types have a lower limit of between 1 and 2 g cm{sup −3}. These results are in agreement with previous density estimates. For 5–20 km diameters, S types rotate faster and have lower amplitudes than C types. If both populations share the same angular momentum, this may indicate the two types’ differing ability to deform under rotational stress. Lastly, we compare our absolute magnitudes (and apparent-magnitude residuals) to those of the Minor Planet Center’s nominal (G = 0.15, rotation-neglecting) model; our phase-function plus Fourier-series fitting reduces asteroid photometric rms scatter by a factor of ∼3.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

The stellar radiation field and the ionization of H and He in the local interstellar medium

Detailed photoionization calculations for the local interstellar medium (LISM) are presented based on constraints imposed by H I column density estimates derived from IUE and Copernicus data toward nearby B stars and hot white dwarfs. It is found that the nearby hot white dwarfs dominate the stellar contributions to the EUV radiation field. Considering stellar contribution to the EUV background alone, the resulting lower limits to the hydrogen and helium fractional ionization in the local diffuse cloud, at the sun, are 0.091 and 0.089, respectively. This result is insensitive to the absorbing geometry of the local cloud. The comparable ionizations of H and He near the sun are due to the lack of absorption from photospheric helium in the hot DA white dwarfs, to the difference in optical depths at wavelengths shortward of the H I and He I ionization edges, and to the fact that the hydrogen ionization controls the electron number density. The derived H and He ionization limits have important implications for the interpretation of the He I and H I backscattering results, for which uncertainties in charge-exchange interactions occurring in or near the heliopause reflect additional significant uncertainties in determining the ionization in the LISM.

Bruhweiler, Frederick C.↗
DOE OSTI2013

THE TWO-DIMENSIONAL PROJECTED SPATIAL DISTRIBUTION OF GLOBULAR CLUSTERS. I. METHOD AND APPLICATION TO NGC 4261

We present a new method for the determination of the two-dimensional (2D) projected spatial distribution of globular clusters (GCs) in external galaxies. This method is based on the K-Nearest Neighbor density estimator of Dressler, complemented by Monte-Carlo simulations to establish the statistical significance of the results. We apply this method to NGC 4261, a ''test galaxy'' where significant 2D anisotropy in the GC distribution has been reported. We confirm that the 2D distribution of GC is not azimuthally isotropic. Moreover, we demonstrate that the 2D distribution departures from the average GC radial distribution results in highly significant spiral-like or broken shell features. Overall, the same perturbations are found in ''red'' and ''blue'' GCs, but with some differences. In particular, we observe a central feature, roughly aligned with the minor axis of NGC 4261, composed of red and most luminous GCs. Blue and fainter GCs are more frequent at large radial distances and follow the spiral-like features of the overall density structure. These results suggest a complex merging history for NGC 4261.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Properties of weak lensing clusters detected on Hyper Suprime-Cam's 2.3 deg{sup 2} field

We present properties of moderately massive clusters of galaxies detected by the newly developed Hyper Suprime-Cam on the Subaru telescope using weak gravitational lensing. Eight peaks exceeding a signal-to-noise ratio (S/N) of 4.5 are identified on the convergence S/N map of a 2.3 deg{sup 2} field observed during the early commissioning phase of the camera. Multi-color photometric data are used to generate optically selected clusters using the Cluster finding algorithm based on the Multiband Identification of Red-sequence galaxies algorithm. The optical cluster positions were correlated with the peak positions from the convergence map. All eight significant peaks have optical counterparts. The velocity dispersion of clusters is evaluated by adopting the Singular Isothemal Sphere fit to the tangential shear profiles, yielding virial mass estimates, M{sub 500{sub c}}, of the clusters which range from 2.7 × 10{sup 13} to 4.4 × 10 {sup 14} M{sub ⊙}. The number of peaks is considerably larger than the average number expected from ΛCDM cosmology but this is not extremely unlikely if one takes the large sample variance in the small field into account. We could, however, safely argue that the peak count strongly favors the recent Planck result suggesting a high σ{sub 8} value of 0.83. The ratio of stellar mass to the dark matter halo mass shows a clear decline as the halo mass increases. If the gas mass fraction, f{sub g}, in halos is universal, as has been suggested in the literature, the observed baryon mass in stars and gas shows a possible deficit compared with the total baryon density estimated from the baryon oscillation peaks in anisotropy of the cosmic microwave background.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2009

Strategies to Improve the Accuracy of Mars-GRAM Sensitivity Studies at Large Optical Depths

The Mars Global Reference Atmospheric Model (Mars-GRAM) is an engineering-level atmospheric model widely used for diverse mission applications. Mars-GRAM s perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL). It has been discovered during the Mars Science Laboratory (MSL) site selection process that Mars-GRAM when used for sensitivity studies for MapYear=0 and large optical depth values such as tau=3 is less than realistic. A comparison study between Mars atmospheric density estimates from Mars- GRAM and measurements by Mars Global Surveyor (MGS) has been undertaken for locations of varying latitudes, Ls, and LTST on Mars. The preliminary results from this study have validated the Thermal Emission Spectrometer (TES) limb data. From the surface to 80 km altitude, Mars- GRAM is based on the NASA Ames Mars General Circulation Model (MGCM). MGCM results that were used for Mars-GRAM with MapYear=0 were from a MGCM run with a fixed value of tau=3 for the entire year at all locations. Unrealistic energy absorption by uniform atmospheric dust leads to an unrealistic thermal energy balance on the polar caps. The outcome is an inaccurate cycle of condensation/sublimation of the polar caps and, as a consequence, an inaccurate cycle of total atmospheric mass and global-average surface pressure. Under an assumption of unchanged temperature profile and hydrostatic equilibrium, a given percentage change in surface pressure would produce a corresponding percentage change in density at all altitudes. Consequently, the final result of a change in surface pressure is an imprecise atmospheric density at all altitudes. To solve this pressure-density problem, a density factor value was determined for tau=.3, 1 and 3 that will adjust the input values of MGCM MapYear 0 pressure and density to achieve a better match of Mars-GRAM MapYear=0 with MapYears 1 and 2 MGCM output at comparable dust loading. Currently, these density factors are fixed values for all latitudes and Ls. Results will be presented of the work underway to derive better multipliers by including possible variation with latitude and/or Ls. This is achieved by comparison of Mars-GRAM MapYear=0 output with TES limb data. The addition of these density factors to Mars-GRAM will improve the results of the sensitivity studies done for large optical depths. Answers may also be provided to the issues raised in a recent study by Desai(2008). Desai has shown that the actual landing sites of Mars Pathfinder, the Mars Exploration Rovers and the Phoenix Mars Lander have been further downrange than predicted by models prior to landing. Desai s reconstruction of their entries into the Martian atmosphere showed that the models consistently predicted higher densities than those found upon EDL. The solution of this problem would be important to the Mars Program since future exploration of Mars by landers and rovers will require more accurate landing capabilities, especially for the proposed Mars Sample Return mission.

Justh, Hilary L.↗
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↗