Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “microbial models”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Fungal and bacterial growth variation due to drought and nitrogen addition experimental treatments. Loma Ridge Experimental Project. 2010-2012

Terrestrial ecosystem models assume that microbial communities respond instantaneously, or are immediately resilient, to environmental change. Here we tested this assumption by quantifying the resilience of a leaf litter community to changes in precipitation or nitrogen availability. By manipulating composition within a global change experiment, we decoupled the legacies of abiotic parameters versus that of the microbial community itself. After one rainy season, more variation in fungal composition could be explained by the original microbial inoculum than the litterbag environment (18% versus 5.5% of total variation). This compositional legacy persisted for 3 years, when 6% of the variability in fungal composition was still explained by the microbial origin. In contrast, bacterial composition was generally more resilient than fungal composition. Microbial functioning (measured as decomposition rate) was not immediately resilient to the global change manipulations; decomposition depended on both the contemporary environment and rainfall the year prior. Finally, using metagenomic sequencing, we showed that changes in precipitation, but not nitrogen availability, altered the potential for bacterial carbohydrate degradation, suggesting why the functional consequences of the two experiments may have differed. Predictions of how terrestrial ecosystem processes respond to environmental change may thus be improved by considering the legacies of microbial communities. This data package includes ten csv files (five data files and their corresponding data dictionaries) and one file-level metadata excel file. Data files contains information about which plots were exposed to treatments related to drought and nitrogen, information about litter bags reciprocal transplants manipulation for water input and nitrogen, detail information about water addition, precipitation records, and litter variables collected. Data dictionary files include detail explanation for each column in the data files. The file-level metadata file describes each file mentioned above. All the analyses were done using the R software.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Small Molecule Interactions with Membranes (CRADA No. TC02321 Final Report)

This was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Procter & Gamble (P&G) to develop and use a multiscale simulation method to explore and evaluate the nature of industrially relevant compounds in the presence of models that represent microbial membranes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Small Molecule Interactions with Membranes. Final CRADA Report

Final Report CRADA No. TC02321 Date Technical Work Ended: 1/29/20. This was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS) as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Procter & Gamble (P&G) to develop and use a multiscale simulation method to explore and evaluate the nature of industrially relevant compounds in the presence of models that represent microbial membranes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Probing interspecies metabolic interactions within a synthetic binary microbiome using genome-scale modeling

Metabolic interactions within a microbial community play a key role in determining the structure, function, and composition of the community. However, due to the complexity and intractability of natural microbiomes, limited knowledge is available on interspecies interactions within a community. In this work, using a binary synthetic microbiome, a methanotroph-photoautotroph (M-P) coculture, as the model system, we examined different genome-scale metabolic modeling (GEM) approaches to gain a better understanding of the metabolic interactions within the coculture, how they contribute to the enhanced growth observed in the coculture, and how they evolve over time. Using batch growth data of the model M-P coculture, we compared three GEM approaches for microbial communities. Two of the methods are existing approaches: SteadyCom, a steady state GEM, and dynamic flux balance analysis (DFBA) Lab, a dynamic GEM. We also proposed an improved dynamic GEM approach, DynamiCom, for the M-P coculture. SteadyCom can predict the metabolic interactions within the coculture but not their dynamic evolutions; DFBA Lab can predict the dynamics of the coculture but cannot identify interspecies interactions. DynamiCom was able to identify the cross-fed metabolite within the coculture, as well as predict the evolution of the interspecies interactions over time. A new dynamic GEM approach, DynamiCom, was developed for a model M-P coculture. Constrained by the predictions from a validated kinetic model, DynamiCom consistently predicted the top metabolites being exchanged in the M-P coculture, as well as the establishment of the mutualistic N-exchange between the methanotroph and cyanobacteria. The interspecies interactions and their dynamic evolution predicted by DynamiCom are supported by ample evidence in the literature on methanotroph, cyanobacteria, and other cyanobacteria-heterotroph cocultures.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Building a genome-based understanding of bacterial pH preferences

The environmental preferences of many microbes remain undetermined. This is the case for bacterial pH preferences, which can be difficult to predict a priori despite the importance of pH as a factor structuring bacterial communities in many systems. We compiled data on bacterial distributions from five datasets spanning pH gradients in soil and freshwater systems (1470 samples), quantified the pH preferences of bacterial taxa across these datasets, and compiled genomic data from representative bacterial taxa. While taxonomic and phylogenetic information were generally poor predictors of bacterial pH preferences, we identified genes consistently associated with pH preference across environments. We then developed and validated a machine learning model to estimate bacterial pH preferences from genomic information alone, a model that could aid in the selection of microbial inoculants, improve species distribution models, or help design effective cultivation strategies. More generally, we demonstrate the value of combining biogeographic and genomic data to infer and predict the environmental preferences of diverse bacterial taxa.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Divergent carbon use efficiency-growth rate tradeoff in popular biological growth models

Carbon use efficiency (CUE) is an important trait emerging from processes regulating biological growth. CUE can be computed either based on the growth of structural biomass or total biomass divided by substrate uptake rate. Nonequilibrium thermodynamics and observations suggest that, for an exponentially growing population of cells, structural biomass CUE should first increase, then peak, and finally decrease with specific growth rate; meanwhile, total biomass CUE increases asymptotically with specific growth rate. We compared predictions from six popular models that are often used for plant and microbial growth in existing ecosystem models. We found that, for an exponentially growing population of biological cells, (1) the source-driven Pirt and Compromise models predict that structural biomass CUE increase asymptotically with growth rate; (2) the apparent sink-driven modified Droop model predicts that structural biomass CUE decreases with growth rate; and (3) the sink-driven variable internal storage model and two dynamic energy budget models predict that structural biomass CUE first increases, then peaks, and finally decreases with growth rate. Moreover, the modified Droop model predicts that total biomass CUE is constant with growth rate, while all other five models predict that total biomass CUE increases with growth rate asymptotically. For non-exponential biological growth, we show that there is no static relationship between total biomass CUE or structural biomass CUE with respect to either growth rate or temperature. Therefore, we contend that biological growth models should explicitly represent interactions between substrate acquisition, substate transformation, and maintenance respiration to better capture observed CUE dynamics, and the sink-driven model should be preferred for general ecosystem biogeochemistry modeling.

Tang, Jinyun [Lawrence Berkeley National Laborator↗
DOE OSTI2025

Toward a Unified Kinetic Model of Nitrogenase Catalysis

The microbial enzyme nitrogenase catalyzes the MgATP-dependent reduction of N 2 to 2NH 3 , a transformation central to the global nitrogen cycle. While the canonical Thorneley−Lowe (TL) kinetic model has long served as a mechanistic framework, it does not incorporate several recent insights. Here, we present an updated kinetic model for Monitrogenase that incorporates these new findings. A significant insight is that electron transfer (ET) from the reduced Fe protein to the FeMo-cofactor is gated by MgATP-dependent conformational transitions and can be described as a probabilistic event that is dependent on the ligand bound to the active-site metallocofactor. The updated kinetic model quantitatively reproduces steady-state product formation rates across a broad range of experimental conditions, yielding revised estimates for key rate constants. It is demonstrated that under N 2 turnover, the probability of productive ET to the active site decreases by ∼60%, resulting in a significant fraction of Fe protein cycles that are unproductive for electron delivery. This mechanistic feature explains the observed rate limitation in N 2 reduction and implies a revised minimum energetic cost of approximately 25 MgATP per N 2 reduced. Integrating these new features into the revised kinetic model provides a more complete and usable foundation for understanding nitrogenase catalysis.

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

A tale of two extremes: Temperature sensitivity of carbon loss from cool and hot soils

Soils represent the largest terrestrial carbon (C) pool, and the flux of carbon dioxide (CO 2 ) from soils to the atmosphere is ~ 6-10 times more than anthropogenic emissions. Understanding responses of soil CO 2 emissions to warming is crucial for evaluating feedback to ongoing environmental changes. The relationship between microbial respiration and temperature is typically modeled using a Q 10 function. Generally, observations of the apparent Q 10 of soil respiration are higher for cold vs. warm ecosystems, reflecting expected biophysical controls of Arrhenius kinetics. However, results from two field warming experiments in the tropics contradict this expectation, both observing extraordinarily high soil respiration responses to in situ warming. Our overall objective for the proposed work is to reduce uncertainty in temperature sensitivity of soil C loss by systematically synthesizing underlying mechanisms related to soil C turnover and stabilization. We are evaluating the temperature sensitivity of soil respiration in ecosystems across temperature extremes (e.g., arctic/boreal and tropical systems) by integrating data collected from field warming experiments with machine learning and biogeochemical models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of a plant carbon–nitrogen interface coupling framework in a coupled biophysical-ecosystem–biogeochemical model (SSiB5/TRIFFID/DayCent-SOM v1.0)

Plant and microbial nitrogen (N) dynamics and N availability regulate the photosynthetic capacity and capture, allocation, and turnover of carbon (C) in terrestrial ecosystems. Studies have shown that a wide divergence in representations of N dynamics in land surface models leads to large uncertainties in the biogeochemical cycle of terrestrial ecosystems and then in climate simulations as well as the projections of future trajectories. In this study, a plant C–N interface coupling framework is developed and implemented in a coupled biophysical-ecosystem–biogeochemical model (SSiB5/TRIFFID/DayCent-SOM v1.0). The main concept and structure of this plant C–N framework and its coupling strategy are presented in this study. This framework takes more plant N-related processes into account. The dynamic ratio (CNR) for each plant functional type (PFT) is introduced to consider plant resistance and adaptation to N availability to better evaluate the plant response to N limitation. Furthermore, when available N is less than plant N demand, plant growth is restricted by a lower maximum carboxylation capacity of RuBisCO (V c,max ), reducing gross primary productivity (GPP). In addition, a module for plant respiration rates is introduced by adjusting the respiration with different rates for different plant components at the same N concentration. Since insufficient N can potentially give rise to lags in plant phenology, the phenological scheme is also adjusted in response to N availability. All these considerations ensure a more comprehensive incorporation of N regulations to plant growth and C cycling. This new approach has been tested systematically to assess the effects of this coupling framework and N limitation on the terrestrial carbon cycle. Long-term measurements from flux tower sites with different PFTs and global satellite-derived products are employed as references to assess these effects. The results show a general improvement with the new plant C–N coupling framework, with more consistent emergent properties, such as GPP and leaf area index (LAI), compared to the observations. The main improvements occur in tropical Africa and boreal regions, accompanied by a decrease in the bias in global GPP and LAI by 16.3 % and 27.1 %, respectively.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Non-growing season plant nutrient uptake controls Arctic tundra vegetation composition under future climate

Abstract Plant growth and distribution in high-latitude tundra ecosystems is strongly limited by nutrient availability and is critical for quantifying centennial-scale carbon-climate interactions. However, land model representations of plant–nutrient interactions are uncertain, leading to poor comparisons with high-latitude observations. Although it has been recognized for decades in the observational community that plants continue to acquire nutrients well past when aboveground activity has ceased, most large-scale land models ignore this process. Here we address the role tundra plant nutrient acquisition during the non-growing season (NGS) has on centennial-scale vegetation growth and dynamics, with a focus on shrub expansion. We apply a well-tested mechanistic model of coupled plant, microbial, hydrological, and thermal dynamics that explicitly represents nutrient acquisition based on plant and microbial traits, thereby allowing a prognostic assessment of NGS nutrient uptake. We first show that the model accurately represents observed seasonality of NGS plant nutrient uptake in a northern Alaskan tundra site. Applying the model across the North America tundra indicates that NGS nutrient uptake is consistent with observations and ranges between ∼5% and 50% of annual uptake, with large spatial variability and dependence on plant functional type. We show that NGS plant nutrient acquisition strongly enhances modeled 21 st century tundra shrub growth and expansion rates. Our results suggest that without NGS nutrient uptake, total shrub aboveground dominance would be ∼50% lower, limited primarily by their inability to grow tall enough to maximize their inherent capacity for light competition. Evergreen shrubs would be more strongly affected because of their relatively lower capacity for nutrient remobilization and acquisition compared to deciduous shrubs. Our results highlight the importance of NGS plant and soil processes on high-latitude biogeochemistry and vegetation dynamics and motivate new observations and model structures to represent these dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Stable carbon isotope values of syndepositional carbonate spherules and micrite record spatial and temporal changes in photosynthesis intensity

Marine and lacustrine carbonate minerals preserve carbon cycle information, and their stable carbon isotope values (δ 13 C) are frequently used to infer and reconstruct paleoenvironmental changes. However, multiple processes can influence the δ 13 C values of bulk carbonates, confounding the interpretation of these values in terms of conditions at the time of mineral precipitation. Co-existing carbonate forms may represent different environmental conditions, yet few studies have analyzed δ 13 C values of syndepositional carbonate grains of varying morphologies to investigate their origins. Here, we combine stable isotope analyses, metagenomics, and geochemical modeling to interpret δ 13 C values of syndepositional carbonate spherules (>500 μm) and fine-grained micrite (<63 μm) from a ~1600-year-long sediment record of a hypersaline lake located on the coral atoll of Kiritimati, Republic of Kiribati (1.9°N, 157.4°W). Petrographic, mineralogic, and stable isotope results suggest that both carbonate fractions precipitate in situ with minor diagenetic alterations. The δ 13 C values of spherules are high compared to the syndepositional micrite and cannot be explained by mineral differences or external perturbations, suggesting a role for local biological processes. We use geochemical modeling to test the hypothesis that the spherules form in the surface microbial mat during peak diurnal photosynthesis when the δ 13 C value of dissolved inorganic carbon is elevated. In contrast, we hypothesize that the micrite may precipitate more continuously in the water as well as in sub-surface, heterotrophic layers of the microbial mat. Both metagenome and geochemical model results support a critical role for photosynthesis in influencing carbonate δ 13 C values. The down-core spherule–micrite offset in δ 13 C values also aligns with total organic carbon values, suggesting that the difference in the δ13C values of spherules and micrite may be a more robust, inorganic indicator of variability in productivity and local biological processes through time than the δ 13 C values of individual carbonate forms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

TUTORIAL: A new custom metabolic model for iron-oxidizing bacteria

In this tutorial narrative, we introduce a novel template developed to enable the creation of stoichiometric genome-scale metabolic models for iron-oxidizing bacteria. We demonstrate the development of this template by applying it to Sideroxydans lithotrophicus ES-1, and validate our model using transcriptomic data (Published in Zhou et al., 2022 AEM). Below, we further show that our template facilitates the modeling of mixotrophic iron-oxidizing bacteria and metagenome-assembled genomes (MAGs), by applying our template to the MAG of the mixotrophic iron oxidizer Leptothrix ochracea (Published in Tothero et al, 2024). This work represents the first instance of a generalized and adaptable template for modeling diverse iron-oxidizing microbial systems, expanding the accessibility and applicability of metabolic modeling in this field.

genome-scale model↗
DOE OSTI2026

Data for Zheng et al. (2025), "AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry"

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystems and soil function. Salinity affects soil carbon cycling by directly impacting microbial activity and indirectly altering soil physicochemical properties, but current models inadequately represent these complexities. This dataset contains the observational and modeling data from Zheng et al. (2025), which described a process-based modeling framework that couples soil solution chemistry with microbial carbon cycling reactions to study the impacts of soil salinization. This conceptual model is implemented numerically into the open-source geochemical program PHREEQC 3.0 (Parkhurst and Appelo, 2013). This dataset consists of: - Figure2_AquaMEND_salinity_buffer: Contains model simulation outputs to assess the impact of three different cation exchange and surface complexation processes on salinity buffering (Fig. 2 from Zheng et al. 2025). - Figure3_Salinity_function: Contains salinity function fitting for literature data (Fig. 3 from Zheng et al. 2025). - Figure4_AquaMEND_microbial_mechanisms: Contains model simulation outputs for testing various microbial process-based hypotheses related to soil salinization, including microbial mortality, carbon use efficiency (CUE), extracellular enzyme activity, and other microbial mechanisms (Fig. 4 from Zheng et al. 2025). - Figure5_AquaMEND_Redox: Contains on model simulation outputs to evaluate shifts among key redox processes, such as aerobic respiration, sulfate reduction, and methanogenesis (Fig.5 from Zheng et al. 2025). - Figure6_AquaMEND_sorption: Contains on model simulation outputs for investigating the effects of salinity on dissolved organic matter (DOM) sorption and desorption processes (Fig. 6 from Zheng et al. 2025). - Figure7_AquaMEND_process_couple: Contains on model simulation outputs for exploring coupled biotic-abiotic processes and their interactions (Fig. 7 from Zheng et al. 2025). - data: Includes datasets used to develop salinity response functions and evaluate salinity buffering capacity. Datasets for MEND model calibration. - database: Contains the `.dat` file required by PHREEQC for model execution. - README.md: A Markdown plain text file describing the computational tools and directories. Files are a mixture of plain text CSV (comma-separated value) and plain text *.dat files written by the model; no special software is required to read them.

EARTH SCIENCE > AGRICULTURE > SOILS > SOIL SALINIT↗
DOE OSTI2021

Multisubstrate DNA stable isotope probing reveals guild structure of bacteria that mediate soil carbon cycling

Significance Soil organic matter contains a large fraction of global C whose fate is largely governed by soil microbial communities. Global C-cycle models can be improved using information about microbial processes in soil, but the inner workings of soil communities remain largely hidden from view. We used stable isotopes to characterize the guild structure of soil bacteria based on how they assimilated C from different sources. Guild structure corresponded with growth dynamics and rrn copy number. We also show that this guild structure explains microbial biogeographical patterns at continental and global scales. These results suggest that microbial contributions to soil C cycling are largely governed by an anabolic process that governs energy allocation to growth, competition, and survival.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2023

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗