Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Functional Data Analysis”

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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cell-Type-Specific Proteomics Analysis of a Small Number of Plant Cells by Integrating Laser Capture Microdissection with a Nanodroplet Sample Processing Platform

Plant organs and tissues contain multiple cell types, which are well organized in 3-dimensional structure to efficiently perform physiological functions such as homeostasis, response to environmental perturbation, pathogen infection. It is critically important to perform molecular measurements at the cell-type-specific level to discover mechanisms and unique features of cell populations that govern differentiation and respond to external perturbations. Although mass spectrometry-based proteomics has been demonstrated as an enabling discovery tool to study plant physiology, conventional approaches require millions of cells to generate robust biological conclusions. Such requirements mask the cell-to-cell heterogeneities and limit the comprehensive profiling of plant proteins at spatially resolved and cell-type-specific resolutions. This protocol describes a recently-developed proteomics workflow for studying a small number of plant cells by integrating laser capture microdissection, microfluidic nanodroplet-based sample preparation, with ultrasensitive liquid chromatography-mass spectrometry. Using poplar as a model tree species, we provide detailed protocols, including plant tissue harvest, tissue preparation, cryosectioning, laser microdissection, protein digestion, mass spectrometry measurement, and data analysis. We show the workflow enables the precise identification and quantification of thousands of proteins from hundreds of isolated plant root and leaf cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Physics-guided dual implicit neural representations for source separation

Significant challenges exist in efficient data analysis of most advanced experimental and observational techniques because the collected signals often include unwanted contributions, such as background and signal distortions, that can obscure the physically relevant information of interest. To address this, we have developed a self-supervised machine-learning approach for source separation using a dual implicit neural representation framework that jointly trains two neural networks: one for approximating distortions of the physical signal of interest and the other for learning the effective background contribution. Our method learns directly from the raw data by minimizing a reconstruction-based loss function without requiring labeled data or pre-defined dictionaries. We demonstrate the effectiveness of our framework by considering a challenging case study involving large-scale simulated, as well as experimental, momentum-energy-dependent inelastic neutron scattering data in a four-dimensional parameter space, characterized by heterogeneous background contributions and unknown distortions to the target signal. The method is found to successfully separate physically meaningful signals from a complex or structured background even when the signal characteristics vary across all four dimensions of the parameter space. An analytical approach that informs the choice of the regularization parameter is presented. Our method offers a versatile framework for addressing source separation problems across diverse domains, ranging from superimposed signals in astronomical measurements to structural features in biomedical image reconstructions.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

The role of moisture in MgCl 2 salt: A multiscale approach to TES performance

This study presents a comprehensive multiscale analysis to evaluate the influence of moisture on the thermal performance of thermal energy storage systems using magnesium chloride (MgCl 2 ) as the phase-changing material. The system uses graphite foam with 90% relative density to enhance thermal conductivity. The analysis includes thermal conductivity calculations and specific heat capacity for different hydrate phases of MgCl 2 : Anhydrous, Mono, Di, Tetra, and Hexa. These properties were evaluated for the first time using the phonon density of states from Density Functional Theory simulations. Results showed that thermal conductivity decreased, while specific heat capacity increased by a factor of two as the phase changed from anhydrous to hexahydrate. Meso-scale models were created to account for the anisotropy of graphite foam and property variations of the MgCl 2 hydrate phases. Asymptotic Expansion homogenization simulations determined the anisotropic thermal conductivity for all phases. This unique methodology improved simulation accuracy, which matched experimental data for anhydrous MgCl 2 . A parametric study examined various operating conditions and their effect on TES performance. It revealed that higher charging temperatures did not enhance exergy efficiency, but increased discharge mass flow rates improved it due to better heat transfer. Thermal performance evaluated by the exergy efficiency remained consistent across hydrate systems under the tested conditions. Furthermore, the study suggests that this uniformity is linked to missing key data, particularly latent heat and melting point for different hydrates. Overall, it highlights the importance of multiscale effects and accurate material properties in designing and optimizing TES systems.

Molten salt degradation↗
DOE OSTI2020

Notes on lattice observables for parton distributions: nongauge theories

We review recent theoretical developments concerning the definition and the renormalization of equal-time correlators that can be computed on the lattice and related to Parton Distribution Functions (PDFs) through a factorization formula. We show how these objects can be studied and analyzed within the framework of a nongauge theory, gaining insight through a one-loop computation. We use scalar field theory as a playground to revise, analyze and present the main features of these ideas, to explore their potential, and to understand their limitations for extracting PDFs. We then propose a framework that would allow to include the available lattice QCD data in a global analysis to extract PDFs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Prevalence and diversity of TAL effector-like proteins in fungal endosymbiotic Mycetohabitans spp.

EndofungalMycetohabitans(formerlyBurkholderia) spp. rely on a type III secretion system to deliver mostly unidentified effector proteins when colonizing their host fungus,Rhizopus microsporus. The one known secreted effector family fromMycetohabitansconsists of homologues of transcription activator-like (TAL) effectors, which are used by plant pathogenicXanthomonasandRalstoniaspp. to activate host genes that promote disease. These ‘BurkholderiaTAL-like (Btl)’ proteins bind corresponding specific DNA sequences in a predictable manner, but their genomic target(s) and impact on transcription in the fungus are unknown. Recent phenotyping of Btl mutants of twoMycetohabitansstrains revealed that the single Btl in oneMycetohabitans endofungorumstrain enhances fungal membrane stress tolerance, while others in aMycetohabitans rhizoxinicastrain promote bacterial colonization of the fungus. The phenotypic diversity underscores the need to assess the sequence diversity and, given that sequence diversity translates to DNA targeting specificity, the functional diversity of Btl proteins. Using a dual approach to maximize capture of Btl protein sequences for our analysis, we sequenced and assembled nineMycetohabitansspp. genomes using long-read PacBio technology and also mined available short-read Illumina fungal–bacterial metagenomes. We show thatbtlgenes are present across diverseMycetohabitansstrains from Mucoromycota fungal hosts yet vary in sequences and predicted DNA binding specificity. Phylogenetic analysis revealed distinct clades of Btl proteins and suggested thatMycetohabitansmight contain more species than previously recognized. Within our data set, Btl proteins were more conserved acrossM. rhizoxinicastrains than acrossM. endofungorum, but there was also evidence of greater overall strain diversity within the latter clade. Overall, the results suggest that Btl proteins contribute to bacterial–fungal symbioses in myriad ways.

Genetics & Heredity↗
DOE OSTI2026

Mitigating Algorithmic Bias in Cancer Site Classification Models

Purpose Integrating artificial intelligence in cancer diagnostics has improved tumor classification beyond rule-based systems. Despite these advancements, these models may still encode demographic biases. We conducted a large-scale, applied bias-probing study of a deep learning–based cancer site classifier to quantify race information encoded in document embeddings. We then evaluated how performance changes when race-correlated embedding dimensions are removed in a post-training sensitivity analysis. Methods The cancer site classifier was trained using 3.5 million electronic cancer pathology reports from six of the National Cancer Institute's SEER registries. We trained a hierarchical self-attention network to generate 400-dimensional document embeddings. These embeddings were used to train two downstream, gradient-boosted decision tree classifiers: one to classify the cancer sites and another to predict racial categories. We identified overlapping features by intersecting the top 50 feature-importance rankings from the site and race models and computed their cumulative feature importance in each model. As a post hoc sensitivity analysis, we progressively pruned these overlapping dimensions, retrained the site model, and compared overall macro-F1 and accuracy, race-stratified macro-F1, and group fairness metrics on the basis of demographic parity and equalized odds before and after pruning. Results The analysis revealed minimal feature overlap between the cancer site and race prediction models, and the cumulative importance scores indicated a negligible influence of racial information on clinical predictions. Post-training pruning of overlapping features did not compromise the models' diagnostic accuracy, with a 0.07% loss in accuracy. Conclusion Our findings demonstrate that HiSAN-generated embeddings from SEER data can be used effectively in cancer site classification without significant demographic bias influencing the outcomes. Post-training pruning therefore functions as a practical audit and sensitivity check.

Shivanna, Abhishek [ORNL] (ORCID:0009000665228593)↗
DOE OSTI2020

Cylinder Test

In 1941 G.I. Taylor proposed the cylinder test as an experiment to determine the properties of explosive products. The test consists of a hollow metal cylinder containing an unreacted cylindrical charge of explosives. The charge is detonated at one end of the cylinder and, as the detonation wave propagates, the product gas expands behind it. As the product gas is produced, the pressure acting on the metal increases and begins to push the walls outward. Consider the time at which the detonation wave has reached the halfway point between the two ends of the cylinder. In front of the wave nothing has changed and the metal tube has the same initial radius. Behind the wave, the radius of the tube has expanded. The farther behind the wave a section of tube is, the larger its radius has become. The longer the wall is exposed to the pressurized gas, the more it expands, so that the radius is largest where the detonation was initiated. The metal most commonly used for this experiment is soft copper. Its material properties are well known and soft copper is also easy to shape with precision machining. During the cylinder test the radial expansion of the cylinder is measured over time and this data is used to determine the amount of work applied to the cylinder wall by the product gas as a function of time. Since volume increases in this problem, the time dependence can be seen as a volume dependence. The volume dependent work yields a path in thermodynamic space. The cylinder test experiment can be modeled using simulation. The simulation gives you the answer, but it cannot tell you why. The goal of the proposed work is to use analytic analysis methods to give insight into the simulation solution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accelerated Materials Design for Molten Salt Technologies Using Innovative High-Throughput Methods

The focus of the project is on building an innovative accelerated materials design platform for molten salt technologies using novel high-throughput methods coupled to data analytics. The main objectives is to predict a FeCrMnNi alloy compositional space with better corrosion resistance than stainless steel 316 and identify new molten salt corrosion mechanisms. The project demonstrates the feasibility to use high-throughput methods coupled to data analytics to accelerate alloy design for extreme environments applications. Using a trained and tested machine learning (ML) model, 2000 FeCrMnNi alloy corrosion rate in molten chloride salts were predicted and a compositional field with corrosion rate lower than 316 stainless steel was identified. The ML model interpretability unveiled multiple features of importance in the model prediction. Some features were expected to be of relative significance, such as work function, surface energy and alloy electronegativity, and the ML model interpretability analysis confirmed those. On the other hand, the most important feature is the diffusion coefficient of Ni in the bulk alloy which indicates that a surface diffusion mechanism plays an important role in the overall corrosion mechanism in molten salts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Application of the Metabolic Modeling Pipeline in KBase to Categorize Reactions, Predict Essential Genes, and Predict Pathways in an Isolate Genome

The DOE Systems Biology Knowledgebase (KBase) platform offers a range of powerful tools for the reconstruction, refinement, and analysis of genome-scale metabolic models built from microbial isolate genomes. In this chapter, we describe and demonstrate these tools in action with an analysis of isoprene production in the Bacillus subtilis DSM genome. Two different methods are applied to build initial metabolic models for the DSM genome, then the models are gapfilled in three different growth conditions. Next, flux balance analysis (FBA) and flux variability analysis (FVA) techniques are applied to both study the growth of these models in minimal media and classify reactions within each model based on essentiality and functionality. The models are applied with the FBA method to predict essential genes, which are then compared to an updated list of essential genes obtained for B. subtilis 168, a very similar strain to the DSM isolate. The models are also applied to simulate Biolog growth conditions, and these results are compared with Biolog data collected for B. subtilis 168. Finally, the DSM metabolic models are applied to explore the pathways and genes responsible for producing isoprene in this strain. These studies demonstrate the accuracy and utility of models generated from the KBase pipelines, as well as exploring the tools available for analyzing these models.

Allen, Benjamin↗
DOE OSTI2021

The picasso map-making code: application to a simulation of the QUIJOTE northern sky survey

ABSTRACT Map-making is an important step for the data analysis of cosmic microwave background (CMB) experiments. It consists of converting the data, which are typically a long, complex, and noisy collection of measurements, into a map, which is an image of the observed sky. We present in this paper a new map-making code named picasso (Polarization and Intensity CArtographer for Scanned Sky Observations), which was implemented to construct intensity and polarization maps from the Multi Frequency Instrument (MFI) of the QUIJOTE (Q-U-I Joint TEnerife) CMB polarization experiment. picasso is based on the destriping algorithm, and is suited to address specific issues of ground-based microwave observations, with a technique that allows the fit of a template function in the time domain, during the map-making step. This paper describes the picasso code, validating it with simulations and assessing its performance. For this purpose, we produced realistic simulations of the QUIJOTE-MFI survey of the northern sky (approximately ∼20 000 deg2), and analysed the reconstructed maps with picasso, using real and harmonic space statistics. We show that, for this sky area, picasso is able to reconstruct, with high fidelity, the injected signal, recovering all the scales with ℓ > 10 in TT, EE, and BB. The signal error is better than 0.001 per cent at 20 < ℓ < 200. Finally, we validated some of the methods that will be applied to the real wide-survey data, like the detection of the CMB anisotropies via cross-correlation analyses. Despite that the implementation of picasso is specific for QUIJOTE-MFI data, it could be adapted to other experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations

Abstract Molecular dynamics (MD) simulations are performed to derive an equation of state (EOS) for helium (He) bubbles in tungsten (W) and to study the growth of He bubbles under a W(100) surface until they burst. We study the growth as a function of the initial nucleation depth of the bubbles. During growth, successive loop-punching events are observed, accompanied by shifts in the depth of the bubble towards the surface. Subsequently, the MD data are used to derive models that describe the conditions that cause the loop punching and bursting events. Simulations have been performed at 500, 933, 1500, 2000, and 2500 K to fit the parameters in the models. To compute the pressure in the bubble at the loop punching and bursting events from the models, we derive an EOS for He bubbles in tungsten with an accompanying volume model to compute the bubble volume for a given number of vacancies ( $$N_\text {V}$$ N V ), He atoms ( $$N_\text {He}$$ N He ), and temperature ( T ). To derive the bubble EOS, we firstly derive the EOS for a free He gas. The derived free-gas EOS can accurately predict all MD data included in the analysis (which span up to 54 GPa at 2500 K). Subsequently, the bubble EOS is derived based on the free-gas EOS by correcting the gas density to account for the interaction between He and W atoms. The EOS for the bubbles is fitted to data from MD simulations of He bubbles in bulk W that span a wide range of gas density and sizes up to about 3 nm in diameter. The pressure of subsurface bubbles at the loop punching events as calculated using the bubble-EOS and the volume model agrees well with the pressure obtained directly from the MD simulations. In the loop punching model, for bubbles consisting of $$N_\text {V}$$ N V vacancies and $$N_\text {He}$$ N He helium atoms, the $$N_\text {He}/N_\text {V}$$ N He / N V ratio that causes the event, the resulting increase in $$N_\text {V}$$ N V , and the associated shift of the bubble depth are formulated as a function of $$N_\text {V}$$ N V and T . In the bursting model, a bubble must simultaneously reach a certain depth and $$N_\text {He}/N_\text {V}$$ N He / N V ratio in order to burst. The burst depth and $$N_\text {He}/N_\text {V}$$ N He / N V are also modeled as a function of $$N_\text {V}$$ N V and T . The majority of the loop punching events occur at bubble pressures between 20 and 60 GPa, depending on the bubble size and temperature. The larger the bubble and the higher the temperature, the lower the bubble pressure. Furthermore, our results indicate that at a higher temperature, a bubble can burst from a deeper region.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗
DOE OSTI2025

A structured framework for predicting sustainable aviation fuel properties using liquid-phase FTIR and machine learning

Sustainable aviation fuels have the potential to improve efficiency, reduce emissions, and enhance energy security. To help identify viable sustainable aviation fuels and accelerate research, machine learning models have been developed to predict relevant physicochemical properties. However, many models have limited applicability, leverage data from complex analytical techniques with confined spectral ranges, or use feature decomposition methods that offer limited interpretability. Using liquid-phase Fourier Transform Infrared (FTIR) spectra, this study presents a structured method for creating accurate and interpretable property prediction models for neat molecules, aviation fuels, and blends. Liquid FTIR spectra can be collected quickly and consistently, offering high reliability, sensitivity, and component specificity using less than 2 ml of sample. The method first decomposes FTIR spectra into fundamental building blocks using non-negative matrix factorization (NMF) to enable scientific analysis of FTIR spectra attributes and fuel properties. The NMF features are then used to create five ensemble models for predicting final boiling point, flash point, freezing point, density at 15°C, and kinematic viscosity at -20°C. All models were trained using experimental property data from neat molecules, aviation fuels, and blends. The models accurately predict key properties across a broad range of neat molecules and representative fuels and blends, while enabling interpretation of relationships between compositional elements, such as functional groups or chemical classes, and their resulting properties. This demonstrates strong potential to support sustainable aviation fuel research and development. The models and data are available on an interactive web tool.

Fourier transform infrared spectroscopy↗
DOE OSTI2021

A probabilistic creep model incorporating test condition, initial damage, and material property uncertainty

Uncertainty is prevalent in the creep resistance of alloys, where at elevated temperature and low pressure, rupture can range across logarithmic decades. In this study, a probabilistic continuum-damage-mechanics (CDM)-based model is derived to capture the uncertainty of creep resistance. To meet this objective, creep data for alloy 304 Stainless Steel is gathered. A constitutive model, “Sinh”, is calibrated deterministically to determine the statistical variability of the material properties. Three sources of uncertainty are injected into the model: test condition (stress and temperature), initial damage, and material properties. Probabilistic simulations are carried out by (a) calibrating probability distribution functions (pdfs) for each source of uncertainty (b) randomly sampling the pdfs using Monte Carlo methods and (c) executing simulations to replicate the uncertain creep behavior. A sensitivity analysis is performed to evaluate the relative effect of each source of uncertainty. In full probabilistic simulations, the cumulative uncertainty of creep behavior is evaluated. The probabilistic model accurately predicts the creep deformation and rupture of the available experiments. The probabilistic model is validated for interpolation but lacks extrapolation ability. Several future works are proposed to further improve the model.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Volatility Basis Set Distributions and Viscosity of Organic Aerosol Mixtures: Insights from Chemical Characterization Using Temperature-Programmed Desorption–Direct Analysis in Real-Time High-Resolution Mass Spectrometry

Quantitative assessment of gas-particle partitioning of individual components within complex atmospheric organic aerosols (OA) mixtures is critical for predicting and comprehending the formation and evolution of OA particles in the atmosphere. This investigation leverages previously documented data obtained through a temperature programmed desorption - direct analysis in real time – high resolution mass spectrometry (TPD-DART-HRMS) platform. This methodology facilitates the bottom-up construction of volatility basis set (VBS) distributions for constituents found in three biogenic secondary organic aerosol (SOA) mixtures produced through the ozonolysis of a-pinene, limonene, and ocimene. The apparent enthalpies (ΔH*, kJ mol -1 ) and saturated vapor mass concentrations (C T *, µg∙m -3 ) of individual SOA components, determined as a function of temperature (T, K), facilitated an assessment of changes in VBS distributions and gas-particle partitioning with respect to T and atmospheric total organic mass loadings (tOM, µg∙m -3 ). Further, the VBS distributions reveal distinct differences in volatilities among monomers, dimers, and trimers, enabling their categorization into separate volatility bins. At the ambient temperature of T = 298 K, only monomers efficiently partition between gas and particle phases across a broad range of atmospherically relevant total organic mass loadings (tOM) values of 1–100 µg∙m -3 . Partitioning of dimers and trimers becomes notable only at T > 360 K and T > 420 K, respectively. The viscosity of SOA mixtures is assessed using a bottom-up calculation approach, incorporating the input of elemental formulas, ΔH*, C T *, and particle-phase mass fractions of the SOA components. Through this approach, we are able to accurately estimate the variations in SOA viscosity that result from the evaporation of its components. These variations are, in turn, influenced by atmospherically relevant changes in tOM and T. Comparison of the calculated SOA viscosity and diffusivity values with literature reported experimental results shows close agreement, thereby validating the employed calculation approach. These findings underscore the significant potential for TPD-DART-HRMS measurements in enabling the untargeted analysis of organic molecules within OA mixtures. This approach facilitates quantitative assessment of their gas-particle partitioning and allows for the estimation of their viscosity and condensed-phase diffusion, thereby contributing valuable insights to atmospheric models.

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

Excited-State Dynamics and Nonlinear Optical Properties of Hyperpolarizable Chromophores Based on Conjugated Bis(terpyridyl)Ru(II) and Palladium and Platinum Porphyrinic Components: Impact of Heavy Metals upon Supermolecular Electro-Optic Properties

A new series of strongly coupled oscillators based upon (porphinato)Pd, (porphinato)Pt, and bis(terpyridyl)-ruthenium(II) building blocks is described. These RuPPd, RuPPt, RuPPdRu, and RuPPtRu chromophores feature bis- (terpyridyl)Ru(II) moieties connected to the (porphinato)metal unit via an ethyne linker that bridges the 4!-terpyridyl and porphyrin macrocycle meso-carbon positions. Pump–probe transient optical data demonstrate sub-picosecond excited singlet-to-triplet-state relaxation. The relaxed lowest-energy triplet (T 1 ) excited states of these chromophores feature absorption manifolds that span the 800–1200 nm spectral region, microsecond triplet-state lifetimes, and large absorptive extinction coefficients [ε(T 1 → T n ) > 4 × 10 4 M –1 cm –1 ]. Dynamic hyperpolarizability (β λ ) values were determined from hyper-Rayleigh light scattering (HRS) measurements carried out at several incident irradiation wavelengths over the 800–1500 nm spectral region. Relative to benchmark RuPZn and RuPZnRu chromophores which showed large βHRS values over the 1200–1600 nm range, RuPPd, RuPPt, RuPPdRu, and RuPPtRu displayed large βHRS values over the 850–1200 nm region. Generalized Thomas–Kuhn sum (TKS) rules and experimental hyperpolarizability values were utilized to determine excited state-to-excited state transition dipole terms from experimental electronic absorption data and thus assessed frequency-dependent β λ values, including two- and three-level contributions for both β zzz and β xzx tensor components to the RuPPd, RuPPt, RuPPdRu, and RuPPtRu hyperpolarizability spectra. Furthermore, these analyses qualitatively rationalize how the βzzz and βxzx tensor elements influence the observed irradiation wavelength-dependent hyperpolarizability magnitudes. The TKS analysis suggests that supermolecules related to RuPPd, RuPPt, RuPPdRu, and RuPPtRu will likely feature intricate dependences of experimentally determined βHRS values as a function of irradiation wavelength that derive from substantial singlet–triplet mixing, and complex interactions among multiple different β tensor components that modulate the long wavelength regime of the nonlinear optical response.

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

Harnessing the Power of Machine Learning and Omics to Identify Environmental Regulation on Microbial Functional Composition for Soil C, N, and P Cycling

Microbial enzyme-mediated soil organic matter (SOM) decomposition regulates many key ecosystem functions, such as elemental cycling, soil carbon sequestration, and soil fertility. However, representing microbial processes in Earth system models (ESMs) remains challenging due to a limited understanding of the spatial patterns of diverse microbial functions responsible for soil carbon (C), nitrogen (N), and phosphorus (P) cycling as well as the underlying mechanisms regulating their relative abundances across various environments. We collected published metagenomics data across the continental US (CONUS) to identify hundreds of microbial genes involved in soil C, N, and P cycling and grouped them into eight enzyme functional classes (EFCs). Each EFC represented a group of gene-encoded potential enzymes that decompose similar soil compounds. By integrating the abundances of omics-informed EFCs with the corresponding environmental information, we trained a machine learning (ML) model to identify key edaphic, climate, and vegetation factors regulating the abundances of each EFC. Quantitative analysis of effects of these factors revealed that the spatial distribution of eight EFCs for soil C, N, and P cycling across CONUS reflected potential resource optimization strategies of microbial communities under nutrient limitation, preferential organic-mineral associations, and climatological stresses. This insight, together with the interpreted ML tool and the CONUS-level benchmark for EFCs abundances, paves the way for parameterizing environmental-regulated microbial functional dynamics in biogeochemical models.

machine learning↗
DOE OSTI2020

Imaging systematics and clustering of DESI main targets

ABSTRACT We evaluate the impact of imaging systematics on the clustering of luminous red galaxies (LRG), emission-line galaxies (ELG), and quasars (QSO) targeted for the upcoming Dark Energy Spectroscopic Instrument (DESI) survey. Using Data Release 7 of the DECam Legacy Survey, we study the effects of astrophysical foregrounds, stellar contamination, differences between north galactic cap and south galactic cap measurements, and variations in imaging depth, stellar density, galactic extinction, seeing, airmass, sky brightness, and exposure time before presenting survey masks and weights to mitigate these effects. With our sanitized samples in hand, we conduct a preliminary analysis of the clustering amplitude and evolution of the DESI main targets. From measurements of the angular correlation functions, we determine power law fits $r_0 = 7.78 \pm 0.26\, h^{-1}$Mpc, γ = 1.98 ± 0.02 for LRGs and $r_0 = 5.45 \pm 0.1\, h^{-1}$Mpc, γ = 1.54 ± 0.01 for ELGs. Additionally, from the angular power spectra, we measure the linear biases and model the scale-dependent biases in the weakly non-linear regime. Both sets of clustering measurements show good agreement with survey requirements for LRGs and ELGs, attesting that these samples will enable DESI to achieve precise cosmological constraints. We also present clustering as a function of magnitude, use cross-correlations with external spectroscopy to infer dN/dz and measure clustering as a function of luminosity, and probe higher order clustering statistics through counts-in-cells moments.

79 ASTRONOMY AND ASTROPHYSICS↗