Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational inference”

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 OSTI2023

Large Ensemble Diagnostic Evaluation of Hydrologic Parameter Uncertainty in the Community Land Model Version 5 (CLM5)

Abstract Land surface models such as the Community Land Model version 5 (CLM5) seek to enhance understanding of terrestrial hydrology and aid in the evaluation of anthropogenic and climate change impacts. However, the effects of parametric uncertainty on CLM5 hydrologic predictions across regions, timescales, and flow regimes have yet to be explored in detail. The common use of the default hydrologic model parameters in CLM5 risks generating streamflow predictions that may lead to incorrect inferences for important dynamics and/or extremes. In this study, we benchmark CLM5 streamflow predictions relative to the commonly employed default hydrologic parameters for 464 headwater basins over the conterminous United States (CONUS). We evaluate baseline CLM5 default parameter performance relative to a large (1,307) Latin Hypercube Sampling‐based diagnostic comparison of streamflow prediction skill using over 20 error measures. We provide a global sensitivity analysis that clarifies the significant spatial variations in parametric controls for CLM5 streamflow predictions across regions, temporal scales, and error metrics of interest. The baseline CLM5 shows relatively moderate to poor streamflow prediction skill in several CONUS regions, especially the arid Southwest and Central U.S. Hydrologic parameter uncertainty strongly affects CLM5 streamflow predictions, but its impacts vary in complex ways across U.S. regions, timescales, and flow regimes. Overall, CLM5's surface runoff and soil water parameters have the largest effects on simulated high flows, while canopy water and evaporation parameters have the most significant effects on the water balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Revisiting the Supernova Engines in the 3C 397 and W49B Supernova Remnants

The nature of the supernova remnants (SNRs) 3C 397 and W49B has long been a subject of debate, with prior studies offering conflicting interpretations between thermonuclear and core-collapse scenarios. To help settle this debate, we present a systematic, spatially resolved, spectroscopic analysis of both remnants using XMM-Newton. By applying multicomponent thermal models, we derive key physical properties including elemental abundances, ejecta temperatures, ambient densities, and explosion energetics. We compare the inferred metal abundance ratios to a wide range of core-collapse and thermonuclear nucleosynthesis models, including new models whose explosion energies differ from the canonical value of 10 51 erg. We find that the observed Fe/Si and Ca/Si ratios in both SNRs are best matched by certain thermonuclear models. However, no model fully reproduces the complete set of observed abundance patterns. In 3C 397, high Fe enrichment and spatial abundance variations suggest interaction with a dense progenitor environment, and W49B’s composition is overall consistent with a thermonuclear origin; however, both require a low-energy (∼10 50 erg) supernova explosion. We additionally map the Fe Kα line centroid energies and find a spread, with W49B falling within the core-collapse region—highlighting both environmental complexity and the limitations of this diagnostic for supernova classification. Our results highlight the need for caution in relying on any single diagnostic or nucleosynthesis model for supernova typing, underscore the need for improved nucleosynthesis models, and motivate future high-resolution, high-throughput observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A DECADE of dwarfs: first detection of weak lensing around spectroscopically confirmed low-mass galaxies

We present the first detection of weak gravitational lensing around spectroscopically confirmed dwarf galaxies, using the large overlap between DESI DR1 spectroscopic data and DECADE/DES weak lensing catalogs. A clean dwarf galaxy sample with well-defined redshift and stellar mass cuts enables excess surface mass density measurements in two stellar mass bins ($\log \rm{M}_*=[8.2, 9.2]~M_\odot$ and $\log \rm{M}_*=[9.2, 10.2]~M_\odot$), with signal-to-noise ratios of $5.6$ and $12.4$ respectively. This signal-to-noise drops to $4.5$ and $9.2$ respectively for measurements without applying individual inverse probability (IIP) weights, which mitigates fiber incompleteness from DESI's targeting. The measurements are robust against variations in stellar mass estimates, photometric shredding, and lensing calibration systematics. Using a simulation-based modeling framework with stellar mass function priors, we constrain the stellar mass-halo mass relation and find a satellite fraction of $\simeq 0.3$, which is higher than previous photometric studies but $1.5σ$ lower than $Λ$CDM predictions. We find that IIP weights have a significant impact on lensing measurements and can change the inferred $f_{\rm{sat}}$ by a factor of two, highlighting the need for accurate fiber incompleteness corrections for dwarf galaxy samples. Our results open a new observational window into the galaxy-halo connection at low masses, showing that future massively multiplexed spectroscopic observations and weak lensing data will enable stringent tests of galaxy formation models and $Λ$CDM predictions.

To, Chun-Hao [Chicago U., Astron. Astrophys. Ctr.;↗
DOE OSTI2023

Can strong substorm-associated MeV electron injections be an important cause of large radiation belt enhancements?

It has become well-established that strong outer radiation belt enhancements are due to wave-driven electron energization by whistler-mode chorus waves. However, in this study, we examine strong MeV electron injections on 10 July 2019 and find substantial evidence that such injections may be a crucial contributor to outer radiation belt enhancement events. For such an examination, it is essential to precisely separate temporal flux changes from spatial variations observed as Van Allen Probes move along their orbits. Employing a new “hourly snapshot” analysis approach, we discover unprecedented details of electron flux evolutions that suggest that for this event, the outer belt enhancement was not continuous but instead intermittent, mostly composed of 4 large discrete injection-driven flux increases. The injections appear as sharp flux increases when observed near apogee. Otherwise, by comparing hourly snapshots for different times, we infer injections and infer temporally stable fluxes between injections, despite strong and continuous chorus emission. The fast and intermittent electron flux growth successively extending earthwards implies cumulative outer belt enhancement via a series of repetitive inward transport associated with injection-induced electric fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Field-level reconstruction from foreground-contaminated 21-cm maps

Current and upcoming 21-cm experiments will soon be able to map 21-cm spatial fluctuations in three dimensions for a wide range of redshifts. However, bright foreground contamination and the nature of radio interferometry create significant challenges, making it difficult to access rich cosmological information from the Fourier modes that lie within the “foreground wedge”. Here, in this work, we introduce two approaches aiming to reconstruct the full 21-cm density field, including the missing modes in the wedge: (a) a field-level inference under an effective field theory (EFT) framework; (b) a diffusion-based deep generative model trained on simulations. Under the EFT framework, we implement a fully differentiable forward model that maps the initial conditions of matter fluctuations to the observed, foreground-filtered 21-cm maps. This enables a gradient-based sampler to simultaneously sample the initial conditions and bias parameters, allowing a physically motivated mode reconstruction. Alternatively, we apply a variational diffusion model to perform 21-cm density reconstruction at the map level. Our model is trained on semi-numerical simulations over a wide range of astrophysical parameters. Our results from both approaches should provide improved cosmological constraints from the field level and also enable cross-correlation between experiments that have little or no overlapping modes.

cosmological perturbation theory↗
DOE OSTI2021

Unnamed Pt(Cu 0.67 Sn 0.33 ) from the Bolshoy Khailyk River, Western Sayans, Russia, and a Review of Related Compounds and Solid Solutions

We describe a potentially new species of a platinum cupride–stannide mineral (PCSM) of composition Pt(Cu 0.67 Sn 0.33 ). It occurs in a placer deposit in the River Bolshoy Khailyk, southern Krasnoyarskiy kray, Russia. A synthetic equivalent of PCSM was obtained and characterized. The PCSM occurs as anhedral or subhedral grains up to 15 μm × 30 μm in association with various platinum-group minerals, Rh–Co-rich pentlandite and magnetite, all hosted by a placer grain of Cu–Au–Pt alloy. Synchrotron micro-Laue diffraction studies indicate that the PCSM mineral is tetragonal and belongs to the inferred space-group P4/mmm (#123). Its unit-cell parameters are a = 2.838 (3) Å, c = 3.650 (4) Å, and V = 29.40 (10) Å 3 , and Z = 1. The c:a ratio calculated from the unit-cell parameters is 1.286. These characteristics are in good agreement with those obtained for specimens of synthetic Pt(Cu 0.67 Sn 0.33 ). A review on related minerals and unnamed phases is provided to outline compositional variations and extents of solid solutions in the relevant systems PtNi–PtFe–PtCu, PdCu–PdHg–PdAu, PdHg–PtHg, and AuCu–PtCu. The PCSM-bearing mineralization appears to be related genetically with an ophiolitic source-rock of the Aktovrakskiy complex of the western Sayans. The unnamed phase likely crystallized from microvolumes of a highly fractionated melt rich in Cu and Sn.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Using ultrasonic attenuation in cortical bone to infer distributions on pore size

Here, in this work we infer the underlying distribution on pore radius in human cortical bone samples using ultrasonic attenuation data. We first discuss how to formulate polydisperse attenuation models using a probabilistic approach and the Waterman Truell model for scattering attenuation. We then compare the Independent Scattering Approximation and the higher-order Waterman Truell models’ forward predictions for total attenuation in polydisperse samples. Following this, we formulate an inverse problem under the Prohorov Metric Framework coupled with variational regularization to stabilize this inverse problem. We then use experimental attenuation data taken from human cadaver samples and solve inverse problems resulting in nonparametric estimates of the probability density function on pore radius. We compare these estimates to the “true” microstructure of the bone samples determined via microCT imaging. We find that our methodology allows us to reliably estimate the underlying microstructure of the bone from attenuation data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Integrating Intermediate Traits in Phylogenetic Genotype-to-Phenotype Studies

A major goal of research in evolution and genetics is linking genotype to phenotype. This work could be direct, such as determining the genetic basis of a phenotype by leveraging genetic variation or divergence in a developmental, physiological, or behavioral trait. The work could also involve studying the evolutionary phenomena (e.g., reproductive isolation, adaptation, sexual dimorphism, behavior) that reveal an indirect link between genotype and a trait of interest. When the phenotype diverges across evolutionarily distinct lineages, this genotype-to-phenotype problem can be addressed using phylogenetic genotype-to-phenotype (PhyloG2P) mapping, which uses genetic signatures and convergent phenotypes on a phylogeny to infer the genetic bases of traits. The PhyloG2P approach has proven powerful in revealing key genetic changes associated with diverse traits, including the mammalian transition to marine environments and transitions between major mechanisms of photosynthesis. However, there are several intermediate traits layered in between genotype and the phenotype of interest, including but not limited to transcriptional profiles, chromatin states, protein abundances, structures, modifications, metabolites, and physiological parameters. Each intermediate trait is interesting and informative in its own right, but synthesis across data types has great promise for providing a deep, integrated, and predictive understanding of how genotypes drive phenotypic differences and convergence. We argue that an expanded PhyloG2P framework (the PhyloG2P matrix) that explicitly considers intermediate traits, and imputes those that are prohibitive to obtain, will allow a better mechanistic understanding of any trait of interest. Furthermore, this approach provides a proxy for functional validation and mechanistic understanding in organisms where laboratory manipulation is impractical.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Combining biomarker and virus phylogenetic models improves HIV-1 epidemiological source identification

To identify and stop active HIV transmission chains new epidemiological techniques are needed. Here, we describe the development of a multi-biomarker augmentation to phylogenetic inference of the underlying transmission history in a local population. HIV biomarkers are measurable biological quantities that have some relationship to the amount of time someone has been infected with HIV. To train our model, we used five biomarkers based on real data from serological assays, HIV sequence data, and target cell counts in longitudinally followed, untreated patients with known infection times. The biomarkers were modeled with a mixed effects framework to allow for patient specific variation and general trends, and fit to patient data using Markov Chain Monte Carlo (MCMC) methods. Subsequently, the density of the unobserved infection time conditional on observed biomarkers were obtained by integrating out the random effects from the model fit. This probabilistic information about infection times was incorporated into the likelihood function for the transmission history and phylogenetic tree reconstruction, informed by the HIV sequence data. To critically test our methodology, we developed a coalescent-based simulation framework that generates phylogenies and biomarkers given a specific or general transmission history. Testing on many epidemiological scenarios showed that biomarker augmented phylogenetics can reach 90% accuracy under idealized situations. Under realistic within-host HIV-1 evolution, involving substantial within-host diversification and frequent transmission of multiple lineages, the average accuracy was at about 50% in transmission clusters involving 5–50 hosts. Realistic biomarker data added on average 16 percentage points over using the phylogeny alone. Using more biomarkers improved the performance. Shorter temporal spacing between transmission events and increased transmission heterogeneity reduced reconstruction accuracy, but larger clusters were not harder to get right. More sequence data per infected host also improved accuracy. We show that the method is robust to incomplete sampling and that adding biomarkers improves reconstructions of real HIV-1 transmission histories. The technology presented here could allow for better prevention programs by providing data for locally informed and tailored strategies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Genome-scale phylogeny and comparative genomics of the fungal order Sordariales

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Deep Learning Prediction of Interspecies Interactions from Self-organized Spatiotemporal Patterns of Co-evolving Organisms

Microorganisms colonizing natural habits such as soils co-evolve to form specific spatial patterns through interspecies interactions. These self-organized patterns are a key ecological phenotype, which provides critical information on their interaction mechanisms. However, conventional network inference techniques that analyze species population data in bulk have yet to be extended to account for such spatial heterogeneity. Here we proposed supervised deep learning as a new network inference tool for predicting interspecies interactions from spatiotemporal patterns of microbial evolution. Due to lack of biological imaging data that can be used for training deep learning networks, we used in silico data generated from high-fidelity agent-based models to determine model structure and parameters. Even though networks were trained under simple configurations where interaction coefficients are assumed to be spatially invariant, we demonstrated that the resulting model can be utilized to successfully predict spatial variation of interactions in more complex domains (i.e., configured with a context-dependent mixture of interaction coefficients) as well as in simple domains without further training. In the further test against real biological data obtained through imaging experiments of a binary consortium (Pseudomonas fluorescens and a mutant of Escherichia coli), our model also predicted the dramatic shifts in interactions of the two organisms across different environmental contexts. Through various successful demonstrations in this work, the combined use of the agent-based model and machine learning algorithm provides a means to use new type of data - microscopic images - for extracting microbial interactions, therefore presenting itself as a useful tool for the analysis of more complex microbial community interactions.

Lee, Joon-Yong↗
DOE OSTI2020

Unifying Advective and Diffusive Descriptions of Bedform Pumping in the Benthic Biolayer of Streams

Abstract Many water quality and ecosystem functions performed by streams occur in the benthic biolayer, the biologically active upper (~5 cm) layer of the streambed. Solute transport through the benthic biolayer is facilitated by bedform pumping, a physical process in which dynamic and static pressure variations over the surface of stationary bedforms (e.g., ripples and dunes) drive flow across the sediment‐water interface. In this paper we derive two predictive modeling frameworks, one advective and the other diffusive, for solute transport through the benthic biolayer by bedform pumping. Both frameworks closely reproduce patterns and rates of bedform pumping previously measured in the laboratory, provided that the diffusion model's dispersion coefficient declines exponentially with depth. They are also functionally equivalent, such that parameter sets inferred from the 2D advective model can be applied to the 1D diffusive model, and vice versa. The functional equivalence and complementary strengths of these two models expand the range of questions that can be answered, for example, by adopting the 2D advective model to study the effects of geomorphic processes (such as bedform adjustments to land use change) on flow‐dependent processes and the 1D diffusive model to study problems where multiple transport mechanisms combine (such as bedform pumping and turbulent diffusion). By unifying 2D advective and 1D diffusive descriptions of bedform pumping, our analytical results provide a straightforward and computationally efficient approach for predicting, and better understanding, solute transport in the benthic biolayer of streams and coastal sediments.

Grant, Stanley B.↗
DOE OSTI2022

Comparison of Energetic Electron Fluxes Measured by GPS and THEMIS Spacecraft in the Inner Magnetosphere

Key elements of space weather models are energetic electron fluxes in the inner magnetosphere and the outer radiation belt. Flux depletion is driven by various loss processes: scattering into atmosphere, magnetopause shadowing. Flux enhancement is driven by various acceleration processes: local wave-particle interactions, radial transport, plasma sheet injections. Many of these processes operate on ~ hour timescales. Such mesoscale flux variations are not well traced by equatorial spacecraft with much longer orbits. Energetic electron detectors onboard the Global Positioning System (GPS) constellation provide a unique opportunity for probing such ~ hour-scale flux variations. Measurements from up to 23 identically instrumented GPS satellites cover a wide energy and L-shell range with a subhour time resolution. However, their orbits are inclined and thus all measurements at L-shell >4.3 are off-equatorial. Here, in this report, we present a comparison of equatorial THEMIS and nonequatorial GPS measurements of omnidirectional ≤600 keV electron fluxes. Such a comparison allows us to derive coefficients for using off-equatorial GPS fluxes to infer the equatorial values. These coefficients depend on particle energy and L-shell. We demonstrate a new data set derived from GPS measurements and discuss how it can be used to investigate mesoscale dynamics of energetic electron fluxes in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Study of the connected four-point correlation function of galaxies from the DESI Data Release 1 luminous red galaxy sample

We present a measurement of the non-Gaussian four-point correlation function (4PCF) from the DESI DR1 luminous red galaxy (LRG) sample. For the gravitationally induced parity-even 4PCF, we detect a signal with a significance of 14.7⁢𝜎 using our fiducial setup. We assess the robustness of this detection through a series of validation tests, including auto and cross-correlation analyses, sky partitioning across multiple patch combinations, and variations in radial scale cuts. Due to the low completeness of the sample, we find that differences in fiber assignment implementation schemes can significantly impact estimation of the covariance and introduce biases in the data vector. After correcting for these effects, all tests yield consistent results. This is one of the first measurements of the connected 4PCF on the DESI LRG sample; the good agreement between the simulation and the data implies that the amplitude of the density fluctuation inferred from the connected 4PCF is consistent with the Planck Λ⁢ CDM cosmology. The methodology and diagnostic framework established in this work provide a foundation for interpreting parity-odd 4PCF.

Cosmology↗
DOE OSTI2020

3D fault architecture controls the dynamism of earthquake swarms

The vibrant evolutionary patterns made by earthquake swarms are incompatible with standard, effectively two-dimensional (2D) models for general fault architecture. Here we leverage advances in earthquake monitoring with a deep-learning algorithm to image a fault zone hosting a 4-year-long swarm in southern California. We infer that fluids are naturally injected into the fault zone from below and diffuse through strike-parallel channels while triggering earthquakes. A permeability barrier initially limits up-dip swarm migration but ultimately is circumvented. This enables fluid migration within a shallower section of the fault with fundamentally different mechanical properties. Our observations provide high-resolution constraints on the processes by which swarms initiate, grow, and arrest. These findings illustrate how swarm evolution is strongly controlled by 3D variations in fault architecture.

58 GEOSCIENCES↗
DOE OSTI2023

Ecological insights into soil health according to the genomic traits and environment-wide associations of bacteria in agricultural soils

Abstract Soil microbiomes are sensitive to current and previous soil conditions, and bacterial ‘bioindicators’ of biological, physical, and chemical soil properties have considerable potential for soil health assessment. However, the lack of ecological or physiological information for most soil microorganisms limits our ability to interpret the associations of bioindicators and, thus, their utility for guiding management. We identified bioindicators of tillage intensity and twelve soil properties used to rate soil health using a 16S rRNA gene-based survey of farmland across North America. We then inferred the genomic traits of bioindicators and evaluated their environment-wide associations (EWAS) with respect to agricultural management practice, disturbance, and plant associations with 89 studies from agroecosystems. Most bioindicators were either positively correlated with biological properties (e.g., organic matter) or negatively correlated with physical and chemical properties. Higher soil health ratings corresponded with smaller genome size and higher coding density, while lower ratings corresponded with larger genomes and higher rrn copy number. Community-weighted genome size explained most variation in health ratings. EWAS linked prominent bioindicators with the impacts of environmental disturbances. Our findings provide ecological insights into bioindicators of soil properties relevant to soil health management, illustrating the tight coupling of microbiome and soil function.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Representation Fusion Framework for Decoupling Diagnostic Information in Multimodal Learning

Modern medicine increasingly relies on multimodal data, ranging from clinical notes to imaging and genomics, to guide diagnosis and treatment. However, integrating these heterogeneous data sources in a principled and interpretable manner remains a major challenge. We present MODES (Multi-mOdal Disentangled Embedding Space), a representation fusion framework that explicitly separates shared and modality-specific factors of variation, offering a structured latent space for multimodal information that improves both prediction and interpretability. By leveraging pre-trained unimodal foundation models, MODES mitigates the dependency on extensive paired datasets, crucial in data-scarce clinical settings. We introduce a masking strategy that optimizes representation dimensionality by eliminating low-information dimensions, to achieve compact, information-rich representations. Our framework demonstrates superior performance in predicting diagnoses and phenotypes compared to unimodal and conventional fusion models. MODES also enables robust diagnostic inference in missing data scenarios, offering an opportunity toward interpretable and efficient multimodal diagnostics in personalized healthcare.

60 APPLIED LIFE SCIENCES↗