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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

How do soil microbes shape ecosystem biogeochemistry in the context of global change?

Earth system models disagree on the future of global soil carbon C, with projections ranging from losses of 70 Pg to gains of 250 Pg by the year 2100 (Todd-Brown et al., 2014). This disagreement is due, in part, to the high degree of uncertainty in the mechanisms that control soil C, the largest dynamic pool of organic C on Earth. It remains unclear how climate change, rising atmospheric CO2 levels, and other environmental changes will combine to affect global soil C. As key engineers of biogeochemical cycling, microbes play a crucial role in soil C responses to these environmental changes. Yet, the path toward representing microbes in soil models remains uncertain. Furthermore, we address four burning questions relevant for modelling microbial ecological, evolutionary, and biogeochemical processes in soils: (1) What is the role of soil microbes in biogeochemical cycling? (2) Does microbial community composition matter for soil carbon cycling? (3) Which ecological and evolutionary processes contribute to functional changes? (4) How should we model microbial eco-evolutionary mechanisms in soil C models?

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Hysteretic temperature sensitivity of wetland CH4 fluxes explained by substrate availability and microbial activity: Model Archive

This Modeling Archive is in support of an NGEE Arctic publication "Hysteretic temperature sensitivity of wetland CH4 fluxes explained by substrate availability and microbial activity" in the Journal Biogeosciences (https://doi.org/10.5194/bg-17-5849-2020), which includes the model data used in the publication. CH4 emissions from terrestrial systems are posited to increase, which can offset mitigation efforts and accelerate climate change. Yet, the accuracy of modeled CH4 emissions is sensitive to the prescribed CH4 production (or emission) temperature dependencies that are currently uncertain. Here, we use a comprehensive biogeochemistry model (ecosys) to investigate factors modulating CH4 production and emission rates across a permafrost thaw gradient encompassing a partly thawed bog and a fully thawed fen. We find that seasonally varying substrate availability drives lower and higher modeled methanogen biomass and activity, and thereby CH4 production, during the earlier and later periods of the thawed season, respectively. Package follows the Model data archiving guidelines with data in a *.zip file with three subfolders containing *.csv files and raw model output files; a data dictionary table (data_dictionary.csv) and two model output description tables (ecosys_plantspecies_output_notes.csv and ecosys_soil_ouput_notes.csv) to explain the format and meaning of individual output variables; and a user guide as a *.pdf. A detailed model description can be found in the supplement of (Grant, 2013). The ecosys source code is available at https://doi:10.5281/zenodo.3906642. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Targeted rRNA depletion enables efficient mRNA sequencing in diverse bacterial species and complex co-cultures

Bacterial mRNA sequencing is inefficient due to the abundance of ribosomal RNA that is challenging to deplete. While commercial kits target rRNA from common bacterial species, they are frequently inefficient when applied to divergent species, including those from environmental isolates. Similarly, other methods typically employ large probe sets that tile the entire length of rRNAs; however, such approaches are infeasible when applied to many species. Therefore, we present EMBR-seq+, which requires fewer than 10 oligonucleotides per rRNA by combining rRNA blocking primers with RNase H-mediated depletion to achieve rRNA removal efficiencies of up to 99% in diverse bacterial species. Furthermore, in more complex microbial co-cultures between Fibrobacter succinogenes strain UWB7 and anaerobic fungi, EMBR-seq+ depleted both bacterial and fungal rRNA, with a fourfold improvement in bacterial rRNA depletion compared with a commercial kit, thereby demonstrating that the method can be applied to non-model microbial mixtures. Notably, for microbes with unknown rRNA sequences, EMBR-seq+ enables rapid iterations in probe design without requiring to start experiments from total RNA. Finally, efficient depletion of rRNA enabled systematic quantification of the reprogramming of the bacterial transcriptome when cultured in the presence of the anaerobic fungi Anaeromyces robustus or Caecomyces churrovis. We observed that F. succinogenes strain UWB7 downregulated several lignocellulose-degrading carbohydrate-active enzymes in the presence of anaerobic gut fungi, suggesting close interactions between two cellulolytic species that specialize in different aspects of biomass breakdown. Thus, EMBR-seq+ enables efficient, cost-effective, and rapid quantification of the transcriptome to gain insights into non-model microbial systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population’s mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn .

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

How microscopic epistasis and clonal interference shape the fitness trajectory in a spin glass model of microbial long-term evolution

The adaptive dynamics of evolving microbial populations takes place on a complex fitness landscape generated by epistatic interactions. The population generically consists of multiple competing strains, a phenomenon known as clonal interference. Microscopic epistasis and clonal interference are central aspects of evolution in microbes, but their combined effects on the functional form of the population’s mean fitness are poorly understood. Here, we develop a computational method that resolves the full microscopic complexity of a simulated evolving population subject to a standard serial dilution protocol. Through extensive numerical experimentation, we find that stronger microscopic epistasis gives rise to fitness trajectories with slower growth independent of the number of competing strains, which we quantify with power-law fits and understand mechanistically via a random walk model that neglects dynamical correlations between genes. We show that increasing the level of clonal interference leads to fitness trajectories with faster growth (in functional form) without microscopic epistasis, but leaves the rate of growth invariant when epistasis is sufficiently strong, indicating that the role of clonal interference depends intimately on the underlying fitness landscape. The simulation package for this work may be found at https://github.com/nmboffi/spin_glass_evodyn .

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1999

Extremes of Survival Achieved by the Radiophile Deinococcus Radiodurans: A Model for Microbial Life on Mars

Deinococcus radiodurans is an extremophile in more than one defined way. First it is extreme in its resistance to freeze drying. Second it is probably uniquely extreme on Earth in its resistance to ionizing radiation. The polyextremophilic capacity of D. radiodurans will be considered. The selection pressures on Mars will then be considered in relation to D. radiodurans in order to support a hypothesis that if microbial life exists on Mars, then it likely includes polyextremophiles.

Daly, M.↗
NASA NTRSDate not supplied

The Pheno-Evo Model: Evolution of Microbial Phenotypic Diversity in 2D Space

Biologists appreciate microbes for their simplicity and predictability: we assume that a population of genetically identical cells in a uniform environment will all behave the same way. While this assumption is often useful, it is also often wrong. Not only might microbes in a clonal population act differently from one another, the differences may be categorical (growing v. non-growing; toxin-sensitive v. tolerant), and this diversity might be an evolved trait conferring increased fitness on the population. How does such phenotypic heterogeneity evolve? And how does a population find the optimal distribution of phenotypes for a given environment? Efforts at modeling microbial phenotypic heterogeneity often focus on populations with two discrete phenotypic types; phenotypes in continuous distributions remain poorly explored. To address this gap and to explore the role of spatial relationships, we use agent-based modeling to simulate a phenotypically diverse population of microbial cells evolving in the presence of periodic toxic stress. Cells on patches in a 2D grid may degrade toxin, suffer damage from toxin, switch phenotype, and reproduce. An individual's phenotype its toxin degradation rate-- for which there is a tradeoff with reproduction-- and the genotype encodes the distribution of phenotypic values in the population. The rate of toxin diffusion determines how individuals influence their neighbors' environments. We examine the effects of toxin concentration, diffusion rate, and environmental predictability on the survival success of populations with different phenotype distributions, and populations' evolutionary trajectories when phenotype distribution is allowed to evolve. We conduct all simulations on the platform NetLogo, which provides a friendly interface allowing users of any experience level to tweak parameters and run their own simulations. We have also created custom tools in R for analyzing and visualizing the results of multiple runs.

microbial evolution↗
DOE OSTI2021

Cross-Feeding of a Toxic Metabolite in a Synthetic Lignocellulose-Degrading Microbial Community

The recalcitrance of complex organic polymers such as lignocellulose is one of the major obstacles to sustainable energy production from plant biomass, and the generation of toxic intermediates can negatively impact the efficiency of microbial lignocellulose degradation. Here, we describe the development of a model microbial consortium for studying lignocellulose degradation, with the specific goal of mitigating the production of the toxin formaldehyde during the breakdown of methoxylated aromatic compounds. Included are Pseudomonas putida, a lignin degrader; Cellulomonas fimi, a cellulose degrader; and sometimes Yarrowia lipolytica, an oleaginous yeast. Unique to our system is the inclusion of Methylorubrum extorquens, a methylotroph capable of using formaldehyde for growth. We developed a defined minimal “Model Lignocellulose” growth medium for reproducible coculture experiments. We demonstrated that the formaldehyde produced by P. putida growing on vanillic acid can exceed the minimum inhibitory concentration for C. fimi, and, furthermore, that the presence of M. extorquens lowers those concentrations. We also uncovered unexpected ecological dynamics, including resource competition, and interspecies differences in growth requirements and toxin sensitivities. Finally, we introduced the possibility for a mutualistic interaction between C. fimi and M. extorquens through metabolite exchange. This study lays the foundation to enable future work incorporating metabolomic analysis and modeling, genetic engineering, and laboratory evolution, on a model system that is appropriate both for fundamental eco-evolutionary studies and for the optimization of efficiency and yield in microbially-mediated biomass transformation.

Methylorubrum extorquens↗
DOE OSTI2021

Representing methane emissions from wet tropical forest soils using microbial functional groups constrained by soil diffusivity

Tropical ecosystems contribute significantly to global emissions of methane (CH 4 ), and landscape topography influences the rate of CH 4 emissions from wet tropical forest soils. However, extreme events such as drought can alter normal topographic patterns of emissions. Here we explain the dynamics of CH 4 emissions during normal and drought conditions across a catena in the Luquillo Experimental Forest, Puerto Rico. Valley soils served as the major source of CH 4 emissions in a normal precipitation year (2016), but drought recovery in 2015 resulted in dramatic pulses in CH 4 emissions from all topographic positions. Geochemical parameters including (i) dissolved organic carbon (C), acetate, and soil pH and (ii) hydrological parameters like soil moisture and oxygen (O 2 ) concentrations varied across the catena. During the drought, soil moisture decreased in the slope and ridge, and O 2 concentrations increased in the valley. We simulated the dynamics of CH 4 emissions with the Microbial Model for Methane Dynamics-Dual Arrhenius and Michaelis–Menten (M3D-DAMM), which couples a microbial functional group CH 4 model with a diffusivity module for solute and gas transport within soil microsites. Contrasting patterns of soil moisture, O 2 , acetate, and associated changes in soil pH with topography regulated simulated CH 4 emissions, but emissions were also altered by rate-limited diffusion in soil microsites. Changes in simulated available substrate for CH 4 production (acetate, CO 2 , and H 2 ) and oxidation (O 2 and CH 4 ) increased the predicted biomass of methanotrophs during the drought event and methanogens during drought recovery, which in turn affected net emissions of CH 4 . A variance-based sensitivity analysis suggested that parameters related to aceticlastic methanogenesis and methanotrophy were most critical to simulate net CH 4 emissions. This study enhanced the predictive capability for CH 4 emissions associated with complex topography and drought in wet tropical forest soils.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Rejuvenation of Spent Media via Supported Emulsion Liquid Membranes

The overall goal of this project is to maximize the reuseability of spent fermentation media. Supported emulsion liquid membrane separation, a highly efficient extraction technique, is used to remove inhibitory byproducts during fermentation; thus, improving the yield while reducing the need for fresh water. The key objectives of this study are: Develop an emulsion liquid membrane system targeting low molecular weight organic acids which has minimal toxicity on a variety of microbial systems; Conduct mass transfer studies to allow proper modeling and design of a supported emulsion liquid membrane system; Investigate the effect of gravity on emulsion coalescence within the membrane unit; Access the effect of water re-use on fermentation yields in a model microbial system; Develop a perfusion-type fermentor utilizing a supported emulsion liquid membrane system to control inhibitory fermentation byproducts; Work for the coming year will focus on the determination of toxicity of various solvents, selection of the emulsifying agents, as well as characterizing the mass transfer of hollow-fiber contactors.

Wiencek, John M.↗
DOE OSTI2025

genomeocean: a pretrained microbial genome foundational model (genomeoceanLLM) v1.0

We present Genomeocean, a foundational genome language model that represents the microbial genome sequences from complex environmental samples. By training on a large, diverse metagenomic dataset, Genomeocean learns species-specific sequence composition and can generate long, realistic open reading frames (ORFs). Our model employs a Byte-pair-encoding (BPE) tokenization strategy, allowing it to efficiently process large genomic datasets and generate long sequences up to 50kb. We demonstrate that fine-tuning Genomeocean can generate novel gene clusters encoding biosynthetic pathways, showcasing its ability to model both fundamental and complex biological processes. Our work establishes Genomeocean as a powerful tool for understanding microbial genome biology and paves the way for its application in a range of fields, from synthetic biology to microbiome research.

Wang, Zhong [Lawrence Berkeley National Laboratory↗
DOE OSTI2021

A Bayesian nonparametric analysis for zero-inflated multivariate count data with application to microbiome study

High-throughput sequencing technology has enabled researchers to profile microbial communities from a variety of environments, but analysis of multivariate taxon count data remains challenging. Here, we develop a Bayesian nonparametric (BNP) regression model with zero inflation to analyse multivariate count data from microbiome studies. A BNP approach flexibly models microbial associations with covariates, such as environmental factors and clinical characteristics. The model produces estimates for probability distributions which relate microbial diversity and differential abundance to covariates, and facilitates community comparisons beyond those provided by simple statistical tests. We compare the model to simpler models and popular alternatives in simulation studies, showing, in addition to these additional community-level insights, it yields superior parameter estimates and model fit in various settings. The model's utility is demonstrated by applying it to a chronic wound microbiome data set and a Human Microbiome Project data set, where it is used to compare microbial communities present in different environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2002

Rejuvenation of Spent Media via Supported Emulsion Liquid Membranes

The overall goal of this project was to maximize the reuseability of spent fermentation media. Supported emulsion liquid membrane separation, a highly efficient extraction technique, was used to remove inhibitory byproducts during fermentation; thus, improve the yield while reducing the need for fresh water. The key objectives of this study were: (1) Develop an emulsion liquid membrane system targeting low molecular weight organic acids which has minimal toxicity on a variety of microbial systems. (2) Conduct mass transfer studies to allow proper modeling and design of a supported emulsion liquid membrane system. (3) Investigate the effect of gravity on emulsion coalescence within the membrane unit. (4) Access the effect of water re-use on fermentation yields in a model microbial system. and (5) Develop a perfusion-type fermentor utilizing a supported emulsion liquid membrane system to control inhibitory fermentation byproducts (not completed due to lack of funds)

Wiencek, John M.↗
DOE OSTI2024

Predictions of rhizosphere microbiome dynamics with a genome-informed and trait-based energy budget model

Abstract Soil microbiomes are highly diverse, and to improve their representation in biogeochemical models, microbial genome data can be leveraged to infer key functional traits. By integrating genome-inferred traits into a theory-based hierarchical framework, emergent behaviour arising from interactions of individual traits can be predicted. Here we combine theory-driven predictions of substrate uptake kinetics with a genome-informed trait-based dynamic energy budget model to predict emergent life-history traits and trade-offs in soil bacteria. When applied to a plant microbiome system, the model accurately predicted distinct substrate-acquisition strategies that aligned with observations, uncovering resource-dependent trade-offs between microbial growth rate and efficiency. For instance, inherently slower-growing microorganisms, favoured by organic acid exudation at later plant growth stages, exhibited enhanced carbon use efficiency (yield) without sacrificing growth rate (power). This insight has implications for retaining plant root-derived carbon in soils and highlights the power of data-driven, trait-based approaches for improving microbial representation in biogeochemical models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

FUN-BioCROP model with litter decomposition parameters derived from the LIDET dataset

This repository contains the code and data necessary to run the FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms, and Protection) model with litter decomposition parameters derived from a modified Monte Carlo simulation that used the Long-term Intersite Decomposition Experiment Team dataset. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of Files: FUNBioCROP_LIDET Study.Rmd R code with FUN-BioCROP model that can be run with 10 different sets of parameters for litter decomposition (Baseline, LIDET, or eight other best parameter sets identified in the modified Monte Carlo simulation. CORPSE Functions_Bioenergy_V2.R Code with CORPSE model functions, called by FUNBioCROP_LIDET Study.Rmd Model Input Data: bulk.csv, bulk_till.csv, rhizo.csv, rhizo_till.csv, litter.csv Initial C and N (kg C or N/m2) pool values for each soil compartment, final values from spin up. All five files have the same columns: (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 Model Input Data: FluxTower_AvgSoilT.csv: Average daily soil temperature (oC) at 10 cm depth at University of Illinois Urbana-Champaign (UIUC) Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: FluxTower_AvgSoilVWC.csv: Average daily soil volumetric water content (VWC) at 10 cm depth at UIUC Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: input_CCS_LIDET Study.csv: This file has daily data to run FUN-BioCROP (Column - Description - Units): yr - calendar year - year doy - day of year (1 to 365) (no leap year) - day anpp - aboveground NPP (DayCent) - kg C/m2/day bnpp - belowground NPP (DayCent) - kg C/m2/day aglivc - live aboveground biomass carbon (DayCent) - kg C/m2 bglivcj - live juvenile fine root biomass carbon (DayCent) - kg C/m2 bglivcm - live mature fine root biomass carbon (DayCent) - kg C/m2 aglivn - live aboveground biomass nitrogen (DayCent) - kg N/m2 bglivnj - live juvenile fine root biomass nitrogen (DayCent) - kg N/m2 bglivnm - live mature fine root biomass nitrogen (DayCent) - kg N/m2 nyr - simulation year - year cult - indicates a cultivation event (0 or 1) crop - indicates a new crop (0 or 1) fert - indicates a fertilizer event (0 or 1) frst - indicates the first day of the growing season (0 or 1) harv - indicates a harvest event (0 or 1) last - indicates the end of the growing season (0 or 1) croptype - crop type (0=none; 1=alfalfa; 2=corn; 3=grass clover pasture; 4=soybean; 5=wheat) cropsrl - crop specific root length - mm/g root cultrhizmix - fraction of rhizosphere mixed with bulk soil during cultivation (0.0-1.0) - fraction cultlitmix - fraction of litter mixed with bulk soil during cultivation (0.0-1.0) - fraction harvremov - fraction of above ground biomass removed during harvest (0.0-1.0) - fraction fertamt - fertilization amount - g N/m2 lifehist - plant life history (0 = annual, 1 = perennial) froot_turnover_c - amount of C in fine root turnover - kg C/m2 froot_turnover_n - amount of N in fine root turnover - kg N/m2 agrd_turnover_c - amount of C in aboveground biomass turnover - kg C/m2 agrd_turnover_n - amount of N in aboveground biomass turnover - kg N/m2 leaf_litter_fastfrac - Fast decomposing fraction of leaf litter (0.0-1.0) - fraction root_litter_fastfrac - Fast decomposing fraction of root litter (0.0-1.0) - fraction root_diameter - root diameter - mm root_length - root length - mm root/m2 rhizo_frac - fraction of total soil volume that is rhizosphere (0.0 - 1.0) - fraction date - date in format YYYY-MM-DD Instructions: Save the model code ("FUN-BioCROP_LIDET Study.Rmd") and accompanything files (data streams and CORPSE function code) in the same folder. In model code "Chunk 3: Load CORPSE Data Streams" set the working directory (setwd) to the folder with the files saved in step #1. In "Chunk 5: Define LIDET parameter sets" select the litter decomposition parameter set to be used in the run, and comment out all other sets. If changing any parameter values, edit them in "Chunk 6: Load parameters." Run all chunks up to and including "Chunk 10: Prepare Data for Export." In "Chunk 11: Export Output Data" edit data frames for export and filenames, as necessary. "Chunk 12: Graph Total Soil C" makes a figure of C remaining over the model run period. Description of each model chunk (in file FUN-BioCROP_LIDET Study.Rmd): Chunk 1: Remove all functions, clear memory. Removes all functions from R environment, clears the memory. Chunk 2: Load Packages. Loads packages necessary to run the code. Chunk 3: Load CORPSE Data Streams. Sets the working directory and loads the data files necessary to run CORPSE. Chunk 4: Load CORPSE Functions. Loads the R script with CORPSE functions from the working directory, "CORPSE Functions_Bioenergy_V2.R". Chunk 5: Define LIDET parameter sets. Has ten different parameter sets for litter decomposition tested in this study: Baseline parameters, LIDET parameters, and the other 8 best performing parameter sets identified in the modified Monte Carlo. To run the model, all but one parameter set must be commented out. Chunk 6: Load Parameters. Loads all fixed parameters to run the model. Data frame with definitions of parameters is in the CORPSE function script "CORPSE Functions_Bioenergy_V2.R" Chunk 7: Prepare Data Streams. Takes data streams loaded in Chunk 3 and puts them in the format necessary to run the model. The model is coded to run at least two sites at a time, so if only one site is being run it must be run in duplicate. Individual data tables of daily values are created in this chunk from the input data file. Chunk 8: Set Initial Conditions. Creates data tables of soil C and N pools for each soil compartment (rhizo_till, rhizo, bulk_till, bulk, litter) and loads initial values into the data tables. Creates lists for each soil compartment to hold model output. Chunk 9: Load FUN Data and Set Up Matrices. Uses DayCent data to calculate FUN input data: root and leaf N demand, total N demand, plant CN, leaf N available for retranslocation, and litter production. Creates matrices for FUN model outputs. Chunk 10: Run Model. Runs the model. Chunk 11: Prepare Data for Export. Combines data from each day saved as lists into data frames for each soil compartment. Adds values from all soil compartments together to calculate total soil values, creates separate data frames for each soil C and N pool (e.g., protected slow C) for the total soil value. Adds different C and N pools together to calculate total soil C and N for all layers. Creates data frame of ratio of protected to unprotected SOC. Organizes FUN data for export. Chunk 12: Export Results. Exports CSV files of model results to the working directory. Chunk 13: Graph Total Soil C. Makes figure of C remaining over time. Related Links: Original FUN-BioCROP model: https://github.com/BrzostekEcologyLab/FUN-BioCROP LIDET dataset: https://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023

Juice, Stephanie↗