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.

378 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cometary Activity Begins at Kuiper Belt Distances: Evidence from C/2017 K2

We study the development of activity in the incoming long-period comet C/2017 K2 over the heliocentric distance range 9 ≲ r {sub H} ≲ 16 au. The comet continues to be characterized by a coma of submillimeter-sized and larger particles ejected at low velocity. In a fixed co-moving volume around the nucleus we find that the scattering cross section of the coma, C, is related to the heliocentric distance by a power law, C∝r{sub H}{sup -s}, with heliocentric index s = 1.14 ± 0.05. This dependence is significantly weaker than the r {sub H} {sup -2} variation of the insolation as a result of two effects. These are, first, the heliocentric dependence of the dust velocity and, second, a lag effect due to very slow-moving particles ejected long before the observations were taken. A Monte Carlo model of the photometry shows that dust production beginning at r {sub H} ~ 35 au is needed to match the measured heliocentric index, with only a slight dependence on the particle size distribution. Mass-loss rates in dust at 10 au are of order 10{sup 3} kg s{sup -1}, while loss rates in gas may be much smaller, depending on the unknown dust to gas ratio. Consequently, the ratio of the nongravitational acceleration to the local solar gravity, α', may, depending on the nucleus size, attain values of ~10{sup -7} ≲ α' ≲ 10{sup -5}, comparable to values found in short-period comets at much smaller distances. Nongravitational acceleration in C/2017 K2 and similarly distant comets, while presently unmeasured, may limit the accuracy with which we can infer the properties of the Oort cloud from the orbits of long-period comets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Constraining hydrodynamic models of inertial confinement fusion implosions using capsule surrogate experiments

We conduct capsule surrogate experiments at the National Ignition Facility to calibrate radiation hydrodynamic simulations to infer hydrodynamic conditions that are not observable in indirect drive ignition implosions. We tune the simulations by applying laser power and cross beam energy transfer (CBET) saturation multipliers to match the observables from capsule surrogate experiments. Shock timing, velocity, and symmetry are measured in liquid D 2 filled Keyhole capsule surrogate experiments and implosion trajectory, stagnation time, and shape time history are measured in in-flight 2D backlit x-ray radiography experiments (“2DConA”) of D 2 gas filled capsule implosions. Calibrated simulations suggest that the N210808 ignition implosion (fusion target gain = 0.7) had a shell mass remaining at stagnation of less than the nominal %5 (3.8%) and resulted in less confinement. For N221204, the shell was made 5.75 μm thicker to trade implosion velocity for increased confinement and resulted in a target gain = 1.5 with a shell mass remaining of 5.7%. Furthermore, a single adjusted model can reproduce all shock timing data as changes are made to shell thickness (79–85 μm) and laser wavelength separation (1.8–4.0 Å). However, for the 2DConA implosions, a 5% variation in the peak power laser multipliers and a 30% variation in late-time CBET between experiments are needed to match the observed stagnation times, in-flight $P_2$ shape, and hot-spot $P_2$ shape. While progress is being made to improve the models in simulations using focused experiments, capsule surrogate experiments will continue to be needed to optimize future ignition designs.

Lasers↗
DOE OSTI2025

The Atacama Cosmology Telescope: DR6 constraints on extended cosmological models

We use new cosmic microwave background (CMB) primary temperature and polarization anisotropy measurements from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) to test foundational assumptions of the standard cosmological model, ΛCDM, and set constraints on extensions to it. We derive constraints from the ACT DR6 power spectra alone, as well as in combination with legacy data from the Planck mission. To break geometric degeneracies, we include ACT and Planck CMB lensing data and baryon acoustic oscillation data from DESI Year-1. To test the dependence of our results on non-ACT data, we also explore combinations replacing Planck with WMAP and DESI with BOSS, and further add supernovae measurements from Pantheon+ for models that affect the late-time expansion history. We verify the near-scale-invariance (running of the spectral index dn s /d ln k = 0.0062 ± 0.0052) and adiabaticity of the primordial perturbations. Neutrino properties are consistent with Standard Model predictions: we find no evidence for new light, relativistic species that are free-streaming (N eff = 2.86 ± 0.13, which combined with astrophysical measurements of primordial helium and deuterium abundances becomes N eff = 2.89 ± 0.11), for non-zero neutrino masses (∑m ν < 0.089 eV at 95% CL), or for neutrino self-interactions. We also find no evidence for self-interacting dark radiation (N idr < 0.134), or for early-universe variation of fundamental constants, including the fine-structure constant (α EM /α EM,0 = 1.0043 ± 0.0017) and the electron mass (m e /m e,0 = 1.0063 ± 0.0056). Our data are consistent with standard big bang nucleosynthesis (we find Y p = 0.2312 ± 0.0092), the COBE/FIRAS-inferred CMB temperature (we find T CMB = 2.698 ± 0.016 K), a dark matter component that is collisionless and with only a small fraction allowed as axion-like particles, a cosmological constant (w = -0.986 ± 0.025), and the late-time growth rate predicted by general relativity (γ = 0.663 ± 0.052). We find no statistically significant preference for a departure from the baseline ΛCDM model. In fits to models invoking early dark energy, primordial magnetic fields, or an arbitrary modified recombination history, we find H 0 = 69.9 +0.8 -1.5 , 69.1 ± 0.5, or 69.6 ± 1.0 km/s/Mpc, respectively; using BOSS instead of DESI BAO data reduces the central values of these constraints by 1–1.5 km/s/Mpc while only slightly increasing the error bars. In general, models introduced to increase the Hubble constant or to decrease the amplitude of density fluctuations inferred from the primary CMB are not favored over ΛCDM by our data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Development of a framework for sequential Bayesian design of experiments: Application to a pilot-scale solvent-based CO 2 capture process

In this paper, a methodology is developed for sequential design of experiments (SDoE) for process systems and applied to a solvent-based CO 2 capture system. In this approach, the prior knowledge of the system is used to prioritize process data collection at specific operating conditions. These data are then incorporated into a Bayesian inference methodology for updating a stochastic model by refining estimations of its underlying parameters, and the updated model is then used to generate the next set of test runs. Thus, the new knowledge obtained from the data is used to guide subsequent iterations of the experimental runs, ensuring that the overall data collection is maximally informative given that most experimental campaigns, especially at pilot or higher-scale plants, are costly, time-consuming, and resource-limited. The test run objective for this work was to minimize the maximum model prediction uncertainty for key output variables, but the methodology is generic and can be readily applied to other test run objectives. This methodology is applied to an aqueous monoethanolamine (MEA) pilot plant campaign at the National Carbon Capture Center (NCCC) in Wilsonville, Alabama, USA. The SDoE framework was utilized for two iterations, while collecting 18 sets of data representing different process conditions, and this resulted in an overall average reduction in uncertainty of approximately 50% in the prediction of CO 2 capture percentage. Moreover, 11 additional data sets were obtained with variation of absorber packing height for further model validation. This work shows the capability of the SDoE framework to maximize learning given limited resources, allowing for the reduction of model uncertainty, which is of great importance for many applications including reduction of technical risk associated with scale-up and economic analysis.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Internal Microstructure Dictates Interactions of Polymer-grafted Nanoparticles in Solution

Understanding the effects of polymer brush architecture on particle interactions in solution is requisite to enable the development of functional materials based on self-assembled polymer-grafted nanoparticles (GNPs). Static and dynamic light scattering of polystyrene-grafted silica particle solutions in toluene reveals that the pair interaction potential, inferred from the second virial coefficient, A 2 , is strongly affected by the grafting density, σ, and degree of polymerization, N, of tethered chains. In the limit of intermediate σ (~0.3 to 0.6 nm -2 ) and high N, A 2 is positive and increases with N. This confirms the good solvent conditions and can be qualitatively rationalized on the basis of a pair interaction potential derived for grafted (brush) particles. In contrast, for high σ > 0.6 nm -2 and low N, A 2 displays an unexpected reversal to negative values, thus indicating poor solvent conditions. These findings are rationalized by means of a simple analysis based on a coarse-grained brush potential, which balances the attractive core-core interactions and the excluded volume interactions imparted by the polymer grafts. The results suggest that the steric crowding of polymer ligands in dense GNP systems may fundamentally alter the interactions between brush particles in solution and highlight the crucial role of architecture (internal microstructure) on the behavior of hybrid materials. The effect of grafting density also illustrates the opportunity to tailor the physical properties of hybrid materials by altering geometry (or architecture) rather than a variation of the chemical composition.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Exposing new taxonomic variation with inflammation — a murine model-specific genome database for gut microbiome researchers

The murine CBA/J mouse model widely supports immunology and enteric pathogen research. This model has illuminated Salmonella interactions with the gut microbiome since pathogen proliferation does not require disruptive pretreatment of the native microbiota, nor does it become systemic, thereby representing an analog to gastroenteritis disease progression in humans. Despite the value to broad research communities, microbiota in CBA/J mice are not represented in current murine microbiome genome catalogs. Here we present the first microbial and viral genomic catalog of the CBA/J murine gut microbiome. Using fecal microbial communities from untreated and Salmonella-infected, highly inflamed mice, we performed genomic reconstruction to determine the impacts on gut microbiome membership and functional potential. From high depth whole community sequencing (~ 42.4 Gbps/sample), we reconstructed 2281 bacterial and 4516 viral draft genomes. Salmonella challenge significantly altered gut membership in CBA/J mice, revealing 30 genera and 98 species that were conditionally rare and unsampled in non-inflamed mice. Additionally, inflamed communities were depleted in microbial genes that modulate host anti-inflammatory pathways and enriched in genes for respiratory energy generation. Our findings suggest decreases in butyrate concentrations during Salmonella infection corresponded to reductions in the relative abundance in members of the Alistipes. Strain-level comparison of CBA/J microbial genomes to prominent murine gut microbiome databases identified newly sampled lineages in this resource, while comparisons to human gut microbiomes extended the host relevance of dominant CBA/J inflammation-resistant strains. This CBA/J microbiome database provides the first genomic sampling of relevant, uncultivated microorganisms within the gut from this widely used laboratory model. Using this resource, we curated a functional, strain-resolved view on how Salmonella remodels intact murine gut communities, advancing pathobiome understanding beyond inferences from prior amplicon-based approaches. Salmonella-induced inflammation suppressed Alistipes and other dominant members, while rarer commensals like Lactobacillus and Enterococcus endure. The rare and novel species sampled across this inflammation gradient advance the utility of this microbiome resource to benefit the broad research needs of the CBA/J scientific community, and those using murine models for understanding the impact of inflammation on the gut microbiome more generally.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Machine Learning-Based Classification of Lignocellulosic Biomass from Pyrolysis-Molecular Beam Mass Spectrometry Data

High-throughput analysis of biomass is necessary to ensure consistent and uniform feedstocks for agricultural and bioenergy applications and is needed to inform genomics and systems biology models. Pyrolysis followed by mass spectrometry such as molecular beam mass spectrometry (py-MBMS) analyses are becoming increasingly popular for the rapid analysis of biomass cell wall composition and typically require the use of different data analysis tools depending on the need and application. Here, the authors report the py-MBMS analysis of several types of lignocellulosic biomass to gain an understanding of spectral patterns and variation with associated biomass composition and use machine learning approaches to classify, differentiate, and predict biomass types on the basis of py-MBMS spectra. Py-MBMS spectra were also corrected for instrumental variance using generalized linear modeling (GLM) based on the use of select ions relative abundances as spike-in controls. Machine learning classification algorithms e.g., random forest, k-nearest neighbor, decision tree, Gaussian Naïve Bayes, gradient boosting, and multilayer perceptron classifiers were used. The k-nearest neighbors (k-NN) classifier generally performed the best for classifications using raw spectral data, and the decision tree classifier performed the worst. After normalization of spectra to account for instrumental variance, all the classifiers had comparable and generally acceptable performance for predicting the biomass types, although the k-NN and decision tree classifiers were not as accurate for prediction of specific sample types. Gaussian Naïve Bayes (GNB) and extreme gradient boosting (XGB) classifiers performed better than the k-NN and the decision tree classifiers for the prediction of biomass mixtures. The data analysis workflow reported here could be applied and extended for comparison of biomass samples of varying types, species, phenotypes, and/or genotypes or subjected to different treatments, environments, etc. to further elucidate the sources of spectral variance, patterns, and to infer compositional information based on spectral analysis, particularly for analysis of data without a priori knowledge of the feedstock composition or identity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Quantifying the Variation in the Number of Donors in Quantum Dots Created Using Atomic Precision Advanced Manufacturing

Atomic-precision advanced manufacturing enables unique silicon quantum electronics built on quantum dots fabricated from small numbers of phosphorus dopants. The number of dopant atoms comprising a dot plays a central role in determining the behavior of charge and spin confined to the dots and thus overall device performance. Here in this work, we use both theoretical and experimental techniques to explore the combined impact of lithographic variation and stochastic kinetics on the number of P incorporations in quantum dots made using these techniques and how this variation changes as a function of the size of the dot. Using a kinetic model of PH3 dissociation augmented with novel reaction barriers, we demonstrate that for a 2 × 3 silicon dimer window the probability that no donor incorporates goes to zero, allowing for certainty in the placement of at least one donor. However, this still comes with some uncertainty in the precise number of incorporated donors (either one or two), and this variability may still impact certain applications. We also examine the impact of the size of the initial lithographic window, finding that the incorporation fraction saturates to δ-layer-like coverage as the circumference-to-area ratio decreases. We predict that this incorporation fraction depends strongly on the dosage of the precursor and that the standard deviation of the number of incorporations scales as ~√n, as would be expected for a sequence of largely independent incorporation events. Finally, we characterize an array of 36 experimentally prepared multidonor 3 × 3 nm lithographic windows with scanning tunneling microscopy, measuring the fidelity of the lithography to the desired array and the final location of PH x fragments within these lithographic windows. We use our kinetic model to examine the expected variability due to the observed lithographic error, predicting a negligible impact on incorporation statistics. We find good agreement between our model and the inferred incorporation locations in these windows from scanning tunneling microscope measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Experimentally-validated correlation analysis reveals new anaerobic methane oxidation partnerships with consortium-level heterogeneity in diazotrophy

Abstract Archaeal anaerobic methanotrophs (“ANME”) and sulfate-reducing Deltaproteobacteria (“SRB”) form symbiotic multicellular consortia capable of anaerobic methane oxidation (AOM), and in so doing modulate methane flux from marine sediments. The specificity with which ANME associate with particular SRB partners in situ, however, is poorly understood. To characterize partnership specificity in ANME-SRB consortia, we applied the correlation inference technique SparCC to 310 16S rRNA amplicon libraries prepared from Costa Rica seep sediment samples, uncovering a strong positive correlation between ANME-2b and members of a clade of Deltaproteobacteria we termed SEEP-SRB1g. We confirmed this association by examining 16S rRNA diversity in individual ANME-SRB consortia sorted using flow cytometry and by imaging ANME-SRB consortia with fluorescence in situ hybridization (FISH) microscopy using newly-designed probes targeting the SEEP-SRB1g clade. Analysis of genome bins belonging to SEEP-SRB1g revealed the presence of a complete nifHDK operon required for diazotrophy, unusual in published genomes of ANME-associated SRB. Active expression of nifH in SEEP-SRB1g within ANME-2b—SEEP-SRB1g consortia was then demonstrated by microscopy using hybridization chain reaction (HCR-) FISH targeting nifH transcripts and diazotrophic activity was documented by FISH-nanoSIMS experiments. NanoSIMS analysis of ANME-2b—SEEP-SRB1g consortia incubated with a headspace containing CH4 and 15N2 revealed differences in cellular 15N-enrichment between the two partners that varied between individual consortia, with SEEP-SRB1g cells enriched in 15N relative to ANME-2b in one consortium and the opposite pattern observed in others, indicating both ANME-2b and SEEP-SRB1g are capable of nitrogen fixation, but with consortium-specific variation in whether the archaea or bacterial partner is the dominant diazotroph.

58 GEOSCIENCES↗
DOE OSTI2024

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗
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

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Recombination smooths the time-signal disrupted by latency in within-host HIV phylogenies

Within-host HIV evolution involves several features that may disrupt standard phylogenetic reconstruction. One important feature is re-activation of latently integrated provirus, which has the potential to disrupt the temporal signal, leading to variation in the branch lengths and apparent evolutionary rates in a tree. Yet, real within-host HIV phylogenies tend to show clear, ladder-like trees structured by the time of sampling. Another important feature is recombination, which violates the fundamental assumption that evolutionary history can be represented by a single bifurcating tree. Thus, recombination complicates the within-host HIV dynamic by mixing genomes and creating evolutionary loop structures that cannot be represented in bifurcating trees. In this paper, we develop a coalescent-based simulator of within-host HIV evolution that includes latency, recombination, and effective population size dynamics that allows us to study the relationship between the true, complex genealogy of within-host HIV evolution, encoded as an Ancestral Recombination Graph (ARG), and the observed phylogenetic tree. To compare our ARG results to the familiar phylogeny format, we calculate the expected bifurcating tree after decomposing the ARG into all unique site trees, their combined distance matrix, and the overall corresponding bifurcating tree. While latency and recombination separately disrupt the phylogenetic signal, remarkably, we find that recombination recovers the temporal signal of within-host HIV evolution caused by latency by mixing fragments of old, latent genomes into the contemporary population. In effect, recombination averages over extant heterogeneity, whether it stems from mixed time-signals or population bottlenecks. Further, we establish that the signals of latency and recombination can be observed in phylogenetic trees despite being an incorrect representation of the true evolutionary history. Using an Approximate Bayesian Computation method, we develop a set of statistical probes to tune our simulation model to nine longitudinally-sampled within-host HIV phylogenies. Because ARGs are exceedingly difficult to infer from real HIV data, our simulation system allows investigating effects of latency, recombination, and population size bottlenecks by matching decomposed ARGs to real data as observed in standard phylogenies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

The Simons Observatory: Combining cross-spectral foreground cleaning with multitracer B -mode delensing for improved constraints on inflation

The Simons Observatory (SO), due to start full science operations in early 2025, aims to set tight constraints on inflationary physics by inferring the tensor-to-scalar ratio r from measurements of cosmic microwave background (CMB) polarization B-modes. Its nominal design including three small-aperture telescopes (SATs) targets a precision σ⁡(r = 0) ≤ 0.003 without delensing. Achieving this goal and further reducing uncertainties requires a thorough understanding and mitigation of other large-scale B-mode sources such as Galactic foregrounds and weak gravitational lensing. We present an analysis pipeline aiming to estimate r by including delensing within a cross-spectral likelihood, and demonstrate it for the first time on SO-like simulations accounting for various levels of foreground complexity, inhomogeneous noise and partial sky coverage. As introduced in an earlier SO delensing paper, lensing B-modes are synthesized using internal CMB lensing reconstructions as well as Planck-like cosmic infrared background maps and LSST-like galaxy density maps. We then extend SO’s power-spectrum-based foreground-cleaning algorithm to include all auto- and cross-spectra between the lensing template and the SAT B-modes in the likelihood function. This allows us to constrain r and the parameters of our foreground model simultaneously. Within this framework, we demonstrate the equivalence of map-based and cross-spectral delensing and use it to motivate an optimized pixel-weighting scheme for power spectrum estimation. We start by validating our pipeline in the simplistic case of uniform foreground spectral energy distributions. In the absence of primordial B-modes, we find that the 1⁢σ statistical uncertainty on r, σ⁡(r), decreases by 37% as a result of delensing. Tensor modes at the level of r = 0.01 are successfully detected by our pipeline. Even when using more realistic foreground models including spatial variations in the dust and synchrotron spectral properties, we obtain unbiased estimates of r both with and without delensing by employing the moment-expansion method. In this case, uncertainties are increased due to the higher number of model parameters, and delensing-related improvements range between 27% and 31%. These results constitute the first realistic assessment of the delensing performance at SO’s nominal sensitivity level.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Measurement of Volatile Compounds for Real-Time Analysis of Soil Microbial Metabolic Response to Simulated Snowmelt

Snowmelt dynamics are a significant determinant of microbial metabolism in soil and regulate global biogeochemical cycles of carbon and nutrients by creating seasonal variations in soil redox and nutrient pools. With an increasing concern that climate change accelerates both snowmelt timing and rate, obtaining an accurate characterization of microbial response to snowmelt is important for understanding biogeochemical cycles intertwined with soil. However, observing microbial metabolism and its dynamics non-destructively remains a major challenge for systems such as soil. Microbial volatile compounds (mVCs) emitted from soil represent information-dense signatures and when assayed non-destructively using state-of-the-art instrumentation such as Proton Transfer Reaction-Time of Flight-Mass Spectrometry (PTR-TOF-MS) provide time resolved insights into the metabolism of active microbiomes. In this study, we used PTR-TOF-MS to investigate the metabolic trajectory of microbiomes from a subalpine forest soil, and their response to a simulated wet-up event akin to snowmelt. Using an information theory approach based on the partitioning of mutual information, we identified mVC metabolite pairs with robust interactions, including those that were non-linear and with time lags. The biological context for these mVC interactions was evaluated by projecting the connections onto the Kyoto Encyclopedia of Genes and Genomes (KEGG) network of known metabolic pathways. Simulated snowmelt resulted in a rapid increase in the production of trimethylamine (TMA) suggesting that anaerobic degradation of quaternary amine osmo/cryoprotectants, such as glycine betaine, may be important contributors to this resource pulse. Unique and synergistic connections between intermediates of methylotrophic pathways such as dimethylamine, formaldehyde and methanol were observed upon wet-up and indicate that the initial pulse of TMA was likely transformed into these intermediates by methylotrophs. Increases in ammonia oxidation signatures (transformation of hydroxylamine to nitrite) were observed in parallel, and while the relative role of nitrifiers or methylotrophs cannot be confirmed, the inferred connection to TMA oxidation suggests either a direct or indirect coupling between these processes. Overall, it appears that such mVC time-series from PTR-TOF-MS combined with causal inference represents an attractive approach to non-destructively observe soil microbial metabolism and its response to environmental perturbation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗