Search NASA⌕ Search

SEARCH · Search NASA

Results for “uncertainty assessment”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Opinion: The importance of historical and paleoclimate aerosol radiative effects

Abstract. Estimating past aerosol radiative effects and their uncertainties is an important topic in climate science. Aerosol radiative effects propagate into large uncertainties in estimates of how present and future climate evolves with changing greenhouse gas emissions. A deeper understanding of how aerosols affected the atmospheric energy budget under past climates is hindered in part by a lack of relevant paleo-observations and in part because less attention has been paid to the problem. Because of the lack of information we do not seek here to determine the change in the radiative forcing due to aerosol changes but rather to estimate the uncertainties in those changes. Here we argue that current uncertainties from emission uncertainties (90 % confidence interval range spanning 2.8 W m−2) are just as large as model spread uncertainties (2.8 W m−2) in calculating preindustrial to present-day aerosol radiative effects. There are no estimates of radiative forcing for important aerosols such as wildfire and dust aerosols in most paleoclimate time periods. However, qualitative analysis of paleoclimate proxies suggests that changes in aerosols between different past climates are similar in magnitude to changes in aerosols between the preindustrial and present day; plus, there is the added uncertainty from the variability in aerosols and fires in the preindustrial. From the limited literature we crudely estimate a paleoclimate aerosol uncertainty for the Last Glacial Maximum relative to preindustrial of 4.8 W m−2, and we estimate the uncertainty in the aerosol feedback in the natural Earth system over the paleoclimate (Last Glacial Maximum to preindustrial) to be about 3.2 W m−2 K−1. In order to more accurately assess the uncertainty in historical aerosol radiative effects, we propose a new model intercomparison project, which would include multiple plausible emission scenarios tested across a range of state-of-the-art climate models over the historical period. These emission scenarios would then be compared to the available independent aerosol observations to constrain which are most probable. In addition, future efforts should work to characterize and constrain paleo-aerosol forcings and uncertainties. Careful propagation of aerosol uncertainties in the literature is required to ensure an accurate quantification of uncertainties in projections of future climate changes.

Environmental Sciences & Ecology↗
DOE OSTI2021

National Opacity Program: Analysis of Opacity-relevant X-ray Emission Spectra (Milestone Report ID 7121)

We have produced high energy density iron plasmas at temperatures above 1 keV and electron densities exceeding 1023 cm -3 (~ 1 g/cm 3 ) using the Orion laser at the Atomic Weapons Establishment. These plasmas were created by irradiating 50 µm diameter layered targets with frequency doubled (λ = 527 nm), 1 ps laser pulses focused to a 100 µm diameter producing an irradiance of ~ 2 x 10 18 W/cm 2 . The buried layer targets consist of 160 nm iron sulfide (FeS), 60 nm potassium chloride (KCl), and 15 nm carbon. The combined layers are tamped on both sides with 3 µm of parylene-N. The x-ray emission from the plasma was measured using two time-resolved, and four time-integrated Bragg crystal spectrometers, as well as one time-integrated imaging system. One time-resolved x-ray spectrometer measured emission from L-shell transitions in highly charged iron, the other from K-shell transitions in helium-like S 14+ and hydrogen-like S 15+ . The time-integrated spectrometers are intensity-calibrated and measured emission from K-shell transitions in sulfur, potassium, chlorine, and both K-shell and L-shell transitions in iron. The density and temperature of the plasma were determined by modeling the x-ray spectra using different spectral and hydrodynamic modeling packages. A brief overview of the uncertainties associated with the measurements and models are presented. We also give an overview of our 1-D HYDRA-DCA radiation-hydrodynamics model and improvements for future work. Our results aid in assessing experimental uncertainties associated with plasma uniformity and with x-ray emission employed as diagnostics in opacity experiments at temperatures and densities not achievable elsewhere and represent a significant step in creating and diagnosing plasmas near LTE. These results are summarized as part of the completion requirements for milestone 7121

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

3D Deep Learning Joint Inversion of Active Seismic Full Waveform and Passive Seismic Traveltime Data for Reservoir Imaging and Uncertainty Quantification

Here, we present deep learning (DL) networks for three-dimensional (3D) joint inversion of active seismic full waveform and passive seismic traveltime data to image reservoirs and their properties and quantify imaging uncertainties. Active seismic full-waveform data can provide high-resolution monitoring images but are collected only intermittently because of their high acquisition cost. In contrast, passive seismic data can be gathered at relatively low cost between regular active surveys, although their imaging quality can be compromised by factors such as low signal-to-noise ratios and limited ray coverage of the target. Although these datasets are routinely acquired together at CO 2 storage sites, their combined inversion within a 3D DL framework has not been previously demonstrated. To our knowledge, this is the first study to address this gap, combining the strength of both data types. For efficient data storage and DL training with large 3D seismic datasets, we use a 3D data matrix in which a random number of passive seismic traveltime data are stored as parabolic envelopes using one-hot encoding and a 3D full-waveform data matrix in which multiple shot gathers are summed. Two network architectures are evaluated: a single-encoder U-Net for single-data type inversion and a dual-encoder U-Net for joint inversion of active and passive seismic data. We also evaluate the single-encoder U-Net for joint inversion by concatenating full-waveform data and traveltime data. We propose a systematic approach for selecting an optimal dropout rate that balances regularization during training and Monte Carlo dropout-based uncertainty quantification during prediction by examining the correlation coefficient between standard deviation and prediction error, along with the training misfit, across a range of dropout rates. 3D DL inversion experiments include five different network configurations, with evaluations under ideal, noisy and dropout-enabled conditions. Both model and data uncertainties are assessed, as well as their combined effects. Across all conditions, the networks consistently predict accurate CO 2 saturation models with low prediction errors, such as a structural similarity index of 0.993 and CO 2 difference of 1.1%. Uncertainty estimates show strong spatial correlation with prediction errors, confirming the effectiveness of the proposed dropout selection approach. The results demonstrate that our DL approach, utilizing compact data representations and appropriate uncertainty quantification, yields accurate subsurface images under various inversion conditions and provides valuable insights into the reliability of predictions.

Um, Evan Schankee [Lawrence Berkeley National Labo↗
DOE OSTI2022

Uncertainty in land carbon budget simulated by terrestrial biosphere models: the role of atmospheric forcing

Global estimates of the land carbon sink are often based on simulations by terrestrial biosphere models (TBMs). The use of a large number of models that differ in their underlying hypotheses, structure and parameters is one way to assess the uncertainty in the historical land carbon sink. Here we show that the atmospheric forcing datasets used to drive these TBMs represent a significant source of uncertainty that is currently not systematically accounted for in land carbon cycle evaluations. We present results from three TBMs each forced with three different historical atmospheric forcing reconstructions over the period 1850–2015. We perform an analysis of variance to quantify the relative uncertainty in carbon fluxes arising from the models themselves, atmospheric forcing, and model-forcing interactions. We find that atmospheric forcing in this set of simulations plays a dominant role on uncertainties in global gross primary productivity (GPP) (75% of variability) and autotrophic respiration (90%), and a significant but reduced role on net primary productivity and heterotrophic respiration (30%). Atmospheric forcing is the dominant driver (52%) of variability for the net ecosystem exchange flux, defined as the difference between GPP and respiration (both autotrophic and heterotrophic respiration). In contrast, for wildfire-driven carbon emissions model uncertainties dominate and, as a result, model uncertainties dominate for net ecosystem productivity. At regional scales, the contribution of atmospheric forcing to uncertainty shows a very heterogeneous pattern and is smaller on average than at the global scale. We find that this difference in the relative importance of forcing uncertainty between global and regional scales is related to large differences in regional model flux estimates, which partially offset each other when integrated globally, while the flux differences driven by forcing are mainly consistent across the world and therefore add up to a larger fractional contribution to global uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Uncertainty of SW Cloud Radiative Effect in Atmospheric Models Due to the Parameterization of Liquid Cloud Optical Properties

Clouds are largely responsible for the spread of climate models predictions. Here we focus on the uncertainties in cloud shortwave radiative effect due to the parameterization of liquid cloud single scattering properties (SSPs) from liquid water content (LWC) and droplet number concentration (N), named parameterization of cloud optical properties. Uncertainties arise from not accounting for the droplet size distribution (DSD)—which affects the estimation of the effective radius (r eff ) and modulates the r eff -dependency of the SSPs—and from averaging SSPs over wide spectral bands. To assess these uncertainties a series of r eff -dependent SSPs parameterizations corresponding to various DSDs and spectral averaging methods are derived and implemented in a radiative code. Combined with the DSD-dependent estimation of r eff they are used to compute the bulk radiative properties (reflectance, transmittance, absorptance) of various clouds (defined in terms of LWC and N), including a homogeneous cloud, more realistic case studies, and outputs of a climate model. The results show that the cloud radiative forcing can vary up to 20% depending on the assumed DSD. Likewise, differences up to 20% are obtained for heating rates. The estimation of r eff is the main source of uncertainty, while the SSPs parameterization contributes to around 20% of the total uncertainty. Spectral averaging is less an issue, except for atmospheric absorption. Overall, global shortwave cloud radiative effect can vary by 6 W m –2 depending on the assumed DSD shape, which is about 13% of the best observational estimate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Final Design for Thermal/Epithermal eXperiments (TEX) with Chloride Absorbers to Provide Validation Benchmarks for Y-12 Electrorefining Facility

Uranium electrorefining operations at Y-12 require validation for chlorine absorption and reflection. Plutonium chloride solution operations at Los Alamos National Laboratory require chlorine absorption validation. The growing need for chlorine validation is evident to the wider criticality safety community, with multiple attendees at the recent TEX 2.0 meeting at Lawrence Livermore National Laboratory (May 2023) requesting validation for chlorine (Idaho National Laboratory (INL), Institute de radioprotection et de surete nucleaire (IRSN), Savannah River Nuclear Site (SRNS), LANL, and Y-12). Los Alamos National Laboratory has recently performed, and benchmarked, an experiment titled Chlorine Worth Study (CWS) for the internal operations at the lab, but due to the difficulties in precisely characterizing the material compositions of the chlorine absorbers it is advantageous to perform a complimentary study with a different chlorine-based absorber material. Furthermore, having a uranium-based vs plutonium-based experiment provides a separate and important validation basis for criticality safety and nuclear data evaluation. The original final design report for the study of chlorine absorption using the TEX-HEU experimental base was presented in 2022, but used the same chlorine-bearing materials that were found to be difficult to characterize in the LANL benchmark. A complete redesign of the experiment has been performed looking at alternative absorber materials in various forms to produce an experiment that is fully characterizable. This report presents five novel chlorine experiments using the TEX-HEU test bed which provides direct comparison to the Y-12 and INL/Terrapower application needs, utilizing sodium chloride (NaCl) absorber plates. The absorber plates will consist of granulated NaCl (≥99.5% pure), which will be fully encapsulated in aluminum tins, providing a simple but effective chlorine-based absorber material that can be completely characterized. Of the five configurations presented in this report, it is expected that two or three configurations will be down selected for the actual experiment, with the other configurations being alternates. Three of the configurations are in the standard configuration, where the absorber is placed directly on the HEU fuel plates, and two in the sandwich configuration, where the absorber is surrounded by polyethylene moderators to force additional neutron thermalization prior to reaching the absorber. There are two thicknesses of NaCl absorber plates: 3/16” and 1/4” active thicknesses (i.e. not including the encapsulation). Both variations of the absorbers have an active absorber radius of 6” and a total radius of 7.5” to match the diameter of the HEU plates, with the outer 1.5” being aluminum encapsulation. The high-density polyethylene (HDPE) moderators are of the thicknesses: 27/16”, 7/4”, 1/8”, 11/16”, and 3/4”. The final configurations have six (one sandwich and one standard configuration), eight (one sandwich and one standard configuration), and 18 total fuel layers (standard configuration). The standard and sandwich configurations were designed such that the differences in moderator thicknesses are 1” HDPE, which were already procured for the original CED-2. The proposed configurations were precisely tuned to closely match the sensitivity profiles and neutron spectra of the Y-12 upset cases and were also compared to the INL/Terrapower upset cases. The assessment of experimental uncertainties of the non-absorber components was predicted to be 0.00114 Δk eff . The assessment of uncertainties resulting from the absorbers was predicted to be 0.00029 Δk eff . This results in a total uncertainty of 0.00118 Δk eff . Many of the largest uncertainties, namely the moderator densities, may be reduced with precision dimensional inspection of the components. The 1” HDPE moderators as well as the HDPE reflectors from the original CED-2 were incorporated in the final designs presented here. Additional HDPE moderator plates must be fabricated to complete the configurations. NaCl absorber plates will by fabricated at LLNL. The total additional cost is expected to be $\$$54,250 for the remaining components. Precise inspection, including dimensional, mass, density, and impurity, is recommended for all components. LLNL estimates that these costs are around $\$$12,000. It is expected, based on previous TEX-HEU experiments, that three weeks of experimental facility time is needed to complete the experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

An Assessment of the Influence of Uncertainty in Temporally Evolving Streamflow Forecasts on Riverine Inundation Modeling

Continental-scale river forecasting platforms forecast streamflow at reaches that can be used as boundary conditions to drive a local-scale flood inundation model. Uncertainty accumulated during this process stems not only from any part of the forecasting chain but can also be caused by the daily variations in weather forcing that keeps evolving as the event advances. This work aims to examine the influence of the evolving forecast streamflow on predicting the maximum inundation for extreme floods. A diagnostic case study was made on the basis of a hindcast of Hurricane Matthew striking the eastern U.S. in 2016. The U.S. National Water Model was one-way coupled to a hydrodynamic inundation model through a developed automated workflow. Although the river forcing has significantly mismatched hydrographs versus observations, the simulated peak water surface elevations and maximum extents were validated to be comparable with the observations, which indicates that the inundation model may not be sensitive to the inherited uncertainty from the weather forcing. Moreover, the uncertainty of the forecast streamflow time series caused only one order of magnitude fewer variations in inundation prediction; this dampening effect may become clearer for extreme events with large areas inundated. In addition, the forecast total volume of stream discharge appears to be an important metric for assessing the performance of river forcing for inundation mapping, as a linear correlation between the total volume and the accuracy of the predicted peak water surface elevation and maximum extent was found, with the coefficients of determination all above 0.8. Extra best-practice experience of running similar operational tasks demonstrated the tradeoff between the cost and accuracy gain.

flood, national water model, RIFT↗
DOE OSTI2022

Derivation and assessment of regional electricity generation emission factors in the USA

Electricity production is one of the largest sources of environmental emissions - especially greenhouse gases (GHGs) - in the USA. Emission factors (EFs) vary from region to region, which requires the use of spatially relevant EF data for electricity production while performing life cycle assessments (LCAs). Uncertainty information, which is sought by LCA practitioners, is rarely supplied with available life cycle inventories (LCIs). To address these challenges, we present a method for collecting data from different sources for electricity generation and environmental emissions; discuss the challenges involved in agglomerating such data; provide relevant suggestions and solutions to merge the information; and calculate EFs for electricity generation processes from various fuel sources for different spatial regions and spatial resolutions. The EFs from the US 2016 Electricity Life Cycle Inventory (eLCI) are analyzed and explored in this study. We also explore the method of uncertainty information derivation for the EFs. We explore the EFs from different technologies across Emissions & Generation Resource Integrated Database (eGRID) regions in the USA. We find that for certain eGRID regions, the same electricity production technology may have worse emissions. This may be a result of the age of the plants in the region, the quality of fuel used, or other underlying factors. Region-wise life cycle impact assessment (LCIA) ISO 14040 impacts for total generation mix activities provide an overview of the total sustainability profile of electricity production in a particular region, rather than only global warming potential (GWP). We also find that, for different LCIA impacts, several eGRID regions are consistently worse than the US average LCIA impact for every unit of electricity generated. Here this work describes the development of an electricity production LCI at different spatial resolutions by combining and harmonizing information from several databases. The inventory consists of emissions, fuel inputs, and electricity and steam outputs from different electricity production technologies located across various regions of the USA. This LCI for electricity production in the USA will prove to be an enormous resource for all LCA researchers - considering the detailed sources of the information and the breadth of emissions covered by it.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Mobility assessment of the BCC and carbide phases in the C-Nb, C-U and Nb-U systems

Uranium carbides with refractory metal additions are considered for Gen IV nuclear reactors and nuclear thermal propulsion as fuels for their high-temperature and corrosion resistant properties. Understanding kinetic effects that dictate microstructural evolution during fabrication and operating conditions is essential to advance technological development of these fuels. This work presents the development of an atomic mobility database for C-Nb-U systems based off available experimental data supported with ab-initio methods. The mobility assessments and uncertainty quantification (using Markov chain Monte Carlo) were conducted in the Kawin software. Carbon diffusion is considered dominant, as metal diffusion is much slower, with niobium diffusion being even slower and rate limiting than uranium metal. We provide a comprehensive and self-consistent thermo-kinetic database that is validated by diffusion couple simulations through Kawin. In conclusion, this enables prediction of microstructural and phase evolution critical for the development and lifetime assessment of next generation nuclear fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Detection of atmospheric rivers with inline uncertainty quantification: TECA-BARD v1.0.1

It has become increasingly common for researchers to utilize methods that identify weather features in climate models. There is an increasing recognition that the uncertainty associated with choice of detection method may affect our scientific understanding. For example, results from the Atmospheric River Tracking Method Intercomparison Project (ARTMIP) indicate that there are a broad range of plausible atmospheric river (AR) detectors and that scientific results can depend on the algorithm used. There are similar examples from the literature on extratropical cyclones and tropical cyclones. It is therefore imperative to develop detection techniques that explicitly quantify the uncertainty associated with the detection of events. We seek to answer the following question: given a “plausible” AR detector, how does uncertainty in the detector quantitatively impact scientific results? We develop a large dataset of global AR counts, manually identified by a set of eight researchers with expertise in atmospheric science, which we use to constrain parameters in a novel AR detection method. We use a Bayesian framework to sample from the set of AR detector parameters that yield AR counts similar to the expert database of AR counts; this yields a set of “plausible” AR detectors from which we can assess quantitative uncertainty. This probabilistic AR detector has been implemented in the Toolkit for Extreme Climate Analysis (TECA), which allows for efficient processing of petabyte-scale datasets. We apply the TECA Bayesian AR Detector, TECA-BARD v1.0.1, to the MERRA-2 reanalysis and show that the sign of the correlation between global AR count and El Niño–Southern Oscillation depends on the set of parameters used.

54 ENVIRONMENTAL SCIENCES↗
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↗
DOE OSTI2020

γ -ray spectroscopy of astrophysically important states in Ca 39

Background: Nova explosions synthesize elements up to A ≃ 40 , and discrepancies exist between calculated and observed abundances of Ar and Ca created in the explosion. The K 38 ( p , γ ) Ca 39 reaction rate has been illustrated to be influential on these isotopic abundances at the endpoint of nova nucleosynthesis. The energies of the three most important resonances, corresponding to J π = 5 / 2 + excited states in the Ca 39 nucleus above the proton separation threshold, are uncertain and one has been measured with conflicting values [ E r = 679 ( 2 ) versus E r = 701 ( 2 ) keV] in previous experiments. Purpose: Reducing the uncertainties on the resonance energies would allow for a better understanding of the reaction rate. To improve these uncertainties, we searched for γ rays from the depopulation of the corresponding excited states in Ca 39 . Methods: We report a new measurement of these resonance energies via the observation of previously unobserved γ -ray transitions. These transitions were observed by studying the Ca 40 ( 3 He , α γ ) Ca 39 reaction with Gammasphere ORRUBA Dual Detectors for Experimental Structure Studies (GODDESS). The updated resonance energies were then used to calculate the K 38 ( p , γ ) Ca 39 reaction rate and assess its uncertainties. Results: In total, 23 new transitions were found, including three γ -ray transitions corresponding to the three J π = 5 / 2 + states of astrophysical interest at energies of 6156.2(16), 6268.8(22), and 6470.8(19) keV. These correspond to resonance energies in the K 38 ( p , γ ) Ca 39 reaction of 386(2), 498(2), and 701(2) keV. Conclusions: Updated K 38 ( p , γ ) Ca 39 reaction rate calculations reflect a reduced upper limit at nova temperatures. However, the lower-than-previously-measured energy of the 498-keV resonance and uncertainty in its resonance strength increases the upper limit of the rate close to previous estimates at 0.4 GK.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Stochastic parametric skeletal dosimetry model for humans: Pediatric and adult computational skeleton phantoms for internal bone marrow dosimetry

Currently, computational phantoms that simulate skeletal tissues are used in active red bone marrow (AM) internal dosimetry. Up-to-date reference computational phantoms recommended by the ICRP are based on the analysis of CT-images of cadavers. Such phantoms have significant disadvantages. One disadvantage is that the assessment of uncertainty due to the population variability of skeleton dimensions and microstructure results from the limited availability of autopsy material. Another disadvantage is the simplified modelling of cortical layer and bone microarchitecture. A method of stochastic parametric skeletal dosimetry modelling of the bone structures – SPSD modelling – has been developed as an alternative to the ICRP reference phantoms. In the framework of this approach, skeletal phantom parameters are evaluated based on extensively reviewed results of published measurements of real bones. The SPSD approach allows for the assessment of both population-average values and their variability. SPSD-phantoms of the skeleton are modelled in voxel representation. They consist of smaller phantoms of the bone sites – segments – described by simple geometric shapes with uniform microarchitecture parameters. Such segmentation makes it possible to account for non-homogeneous skeletal microarchitecture and to model the bone structure with the required voxel resolution to elaborate suitable skeletal phantoms. The current study presents the parameters of the SPSD skeletal phantoms for the following age-groups: newborn, 1-year-old, 5-year-old, 10-year-old, 15-year-old (male and female), and adult (male and female). This skeletal phantom can be used for dosimetry as an alternative to available reference phantoms for bone-seeking radionuclides. The above-mentioned age- and sex-specific skeletal phantoms are comprised of 289 unique segments. The characteristics of the SPSD phantoms do not contradict published data and are in good agreement with the measurement results of real bones.

Science & Technology - Other Topics↗
DOE OSTI2023

Bootstrapping outperforms community‐weighted approaches for estimating the shapes of phenotypic distributions

Abstract Estimating phenotypic distributions of populations and communities is central to many questions in ecology and evolution. These distributions can be characterized by their moments (mean, variance, skewness and kurtosis) or diversity metrics (e.g. functional richness). Typically, such moments and metrics are calculated using community‐weighted approaches (e.g. abundance‐weighted mean). We propose an alternative bootstrapping approach that allows flexibility in trait sampling and explicit incorporation of intraspecific variation, and show that this approach significantly improves estimation while allowing us to quantify uncertainty. We assess the performance of different approaches for estimating the moments of trait distributions across various sampling scenarios, taxa and datasets by comparing estimates derived from simulated samples with the true values calculated from full datasets. Simulations differ in sampling intensity (individuals per species), sampling biases (abundance, size), trait data source (local vs. global) and estimation method (two types of community‐weighting, two types of bootstrapping). We introduce the traitstrap R package, which contains a modular and extensible set of bootstrapping and weighted‐averaging functions that use community composition and trait data to estimate the moments of community trait distributions with their uncertainty. Importantly, the first function in the workflow, trait_fill , allows the user to specify hierarchical structures (e.g. plot within site, experiment vs. control, species within genus) to assign trait values to each taxon in each community sample. Across all taxa, simulations and metrics, bootstrapping approaches were more accurate and less biased than community‐weighted approaches. With bootstrapping, a sample size of 9 or more measurements per species per trait generally included the true mean within the 95% CI. It reduced average percent errors by 26%–74% relative to community‐weighting. Random sampling across all species outperformed both size‐ and abundance‐biased sampling. Our results suggest randomly sampling ~9 individuals per sampling unit and species, covering all species in the community and analysing the data using nonparametric bootstrapping generally enable reliable inference on trait distributions, including the central moments, of communities. By providing better estimates of community trait distributions, bootstrapping approaches can improve our ability to link traits to both the processes that generate them and their effects on ecosystems.

Maitner, Brian S.↗
DOE OSTI2025

D–MOPH–25: diverse MOF–molecule pairs for Henry’s constants prediction

Computational methods like grand-canonical Monte Carlo simulations and machine learning (ML) have accelerated metal–organic frameworks (MOF) exploration but are typically limited to a narrow range of adsorbates due to data availability and force field constraints. In this study, we introduce a dataset of diverse MOF–molecule pairs for Henry’s constant prediction, D–MOPH–25, which systematically explores a diverse chemical space by combining 113 molecular adsorbates with over 5000 MOF structures through an active learning process. D–MOPH–25 constitutes the most diverse adsorbate dataset used in any ML study of molecular adsorption in MOFs to date. Our workflow builds a benchmark for predicting Henry’s constants at 300 K, leveraging conformal prediction for uncertainty quantification. Assessment through Shannon entropy and uniform manifold approximation and projection confirms the comprehensiveness of D–MOPH–25 while highlighting the importance of robust classification to filter out unphysical data points in regression tasks. Although future enhancements in model architecture and sampling criteria could improve predictive performance, our dataset already spans the target space using only 2.31% of total possibilities. This comprehensive dataset facilitates assessment of model generalizability across adsorbate species and can establish a foundation for high-throughput MOF screening and ML-driven separation processes.

active learning↗
DOE OSTI2020

Integrated hydrogeophysical modelling and data assimilation for geoelectrical leak detection

Time-lapse electrical resistivity tomography (ERT) measurements provide indirect observations of hydrological processes in the Earth's shallow subsurface at high spatial and temporal resolution. ERT has been used in the past decades to detect leaks and monitor the evolution of associated contaminant plumes. Specifically, inverted resistivity images allow visualization of the dynamic changes in the structure of the plume. However, existing methods do not allow the direct estimation of leak parameters (e.g. leak rate, location, etc.) and their uncertainties. We propose an ensemble-based data assimilation framework that evaluates proposed hydrological models against observed time-lapse ERT measurements without directly inverting for the resistivities. Each proposed hydrological model is run through the parallel coupled hydro-geophysical simulation code PFLOTRAN-E4D to obtain simulated ERT measurements. The ensemble of model proposals is then updated using an iterative ensemble smoother. In this paper, we demonstrate the proposed framework on synthetic and field ERT data from controlled tracer injection experiments. Our results show that the approach allows joint identification of contaminant source location, initial release time, and solute loading from the cross-borehole time-lapse ERT data, alongside with an assessment of uncertainties in these estimates. We demonstrate a reduction in site-wide uncertainty by comparing the prior and posterior plume mass discharges at a selected image plane. This framework is particularly attractive to sites that have previously undergone extensive geological investigation (e.g., nuclear sites). It is well suited to complement ERT imaging and we discuss practical issues in its application to field problems.

58 GEOSCIENCES↗
DOE OSTI2023

MaNGA DynPop – I. Quality-assessed stellar dynamical modelling from integral-field spectroscopy of 10K nearby galaxies: a catalogue of masses, mass-to-light ratios, density profiles, and dark matter

ABSTRACT This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here, we present a catalogue of dynamically determined quantities for over 10 000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically aligned JAMcyl or spherically aligned JAMsph. We also make four assumptions for the models’ dark versus luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results.

Astronomy & Astrophysics↗
DOE OSTI2021

A comprehensive and synthetic dataset for global, regional, and national greenhouse gas emissions by sector 1970–2018 with an extension to 2019

To track progress towards keeping global warming well below 2 °C or even 1.5 °C, as agreed in the Paris Agreement, comprehensive up-to-date and reliable information on anthropogenic emissions and removals of greenhouse gas (GHG) emissions is required. Here we compile a new synthetic dataset on anthropogenic GHG emissions for 1970–2018 with a fast-track extension to 2019. Our dataset is global in coverage and includes CO 2 emissions, CH 4 emissions, N 2 O emissions, as well as those from fluorinated gases (F-gases: HFCs, PFCs, SF 6 , NF 3 ) and provides country and sector details. We build this dataset from the version 6 release of the Emissions Database for Global Atmospheric Research (EDGAR v6) and three bookkeeping models for CO 2 emissions from land use, land-use change, and forestry (LULUCF). We assess the uncertainties of global greenhouse gases at the 90 % confidence interval (5th–95th percentile range) by combining statistical analysis and comparisons of global emissions inventories and top-down atmospheric measurements with an expert judgement informed by the relevant scientific literature. We identify important data gaps for F-gas emissions. The agreement between our bottom-up inventory estimates and top-down atmospheric-based emissions estimates is relatively close for some F-gas species (~ 10 % or less), but estimates can differ by an order of magnitude or more for others. Our aggregated F-gas estimate is about 10 % lower than top-down estimates in recent years. However, emissions from excluded F-gas species such as chlorofluorocarbons (CFCs) or hydrochlorofluorocarbons (HCFCs) are cumulatively larger than the sum of the reported species. Using global warming potential values with a 100-year time horizon from the Sixth Assessment Report by the Intergovernmental Panel on Climate Change (IPCC), global GHG emissions in 2018 amounted to 58 ± 6.1 GtCO 2 eq. consisting of CO 2 from fossil fuel combustion and industry (FFI) 38 ± 3.0 GtCO 2 , CO 2 -LULUCF 5.7 ± 4.0 GtCO 2 , CH 4 10 ± 3.1 GtCO 2 eq., N2O 2.6 ± 1.6 GtCO 2 eq., and F-gases 1.3 ± 0.40 GtCO 2 eq. Initial estimates suggest further growth of 1.3 GtCO 2 eq. in GHG emissions to reach 59 ± 6.6 GtCO 2 eq. by 2019. Our analysis of global trends in anthropogenic GHG emissions over the past 5 decades (1970–2018) highlights a pattern of varied but sustained emissions growth. There is high confidence that global anthropogenic GHG emissions have increased every decade, and emissions growth has been persistent across the different (groups of) gases. There is also high confidence that global anthropogenic GHG emissions levels were higher in 2009–2018 than in any previous decade and that GHG emissions levels grew throughout the most recent decade. While the average annual GHG emissions growth rate slowed between 2009 and 2018 (1.2 % yr –1 ) compared to 2000–2009 (2.4 % yr –1 ), the absolute increase in average annual GHG emissions by decade was never larger than between 2000–2009 and 2009–2018. Our analysis further reveals that there are no global sectors that show sustained reductions in GHG emissions. There are a number of countries that have reduced GHG emissions over the past decade, but these reductions are comparatively modest and outgrown by much larger emissions growth in some developing countries such as China, India, and Indonesia. There is a need to further develop independent, robust, and timely emissions estimates across all gases. As such, tracking progress in climate policy requires substantial investments in independent GHG emissions accounting and monitoring as well as in national and international statistical infrastructures. The data associated with this article (Minx et al., 2021) can be found at https://doi.org/10.5281/zenodo.5566761.

54 ENVIRONMENTAL SCIENCES↗