Search NASA⌕ Search

SEARCH · Search NASA

Results for “soil carbon model”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Genomic fingerprints of the world’s soil ecosystems

Despite the explosion of soil metagenomic data, we lack a synthesized understanding of patterns in the distribution and functions of soil microorganisms. These patterns are critical to predictions of soil microbiome responses to climate change and resulting feedbacks that regulate greenhouse gas release from soils. To address this gap, we assay 1,512 manually curated soil metagenomes using complementary annotation databases, read-based taxonomy, and machine learning to extract multidimensional genomic fingerprints of global soil microbiomes. Our objective is to uncover novel biogeographical patterns of soil microbiomes across environmental factors and ecological biomes with high molecular resolution. We reveal shifts in the potential for (i) microbial nutrient acquisition across pH gradients; (ii) stress-, transport-, and redox-based processes across changes in soil bulk density; and (iii) greenhouse gas emissions across biomes. We also use an unsupervised approach to reveal a collection of soils with distinct genomic signatures, characterized by coordinated changes in soil organic carbon, nitrogen, and cation exchange capacity and in bulk density and clay content that may ultimately reflect soil environments with high microbial activity. Genomic fingerprints for these soils highlight the importance of resource scavenging, plant-microbe interactions, fungi, and heterotrophic metabolisms. Across all analyses, we observed phylogenetic coherence in soil microbiomes—more closely related microorganisms tended to move congruently in response to soil factors. Collectively, the genomic fingerprints uncovered here present a basis for global patterns in the microbial mechanisms underlying soil biogeochemistry and help beget tractable microbial reaction networks for incorporation into process-based models of soil carbon and nutrient cycling.

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

Stoichiometrically coupled carbon and nitrogen cycling in the MIcrobial-MIneral Carbon Stabilization model version 1.0 (MIMICS-CN v1.0)

Explicit consideration of microbial physiology in soil biogeochemical models that represent coupled carbon–nitrogen dynamics presents opportunities to deepen understanding of ecosystem responses to environmental change. The MIcrobial-MIneral Carbon Stabilization (MIMICS) model explicitly represents microbial physiology and physicochemical stabilization of soil carbon (C) on regional and global scales. Here we present a new version of MIMICS with coupled C and nitrogen (N) cycling through litter, microbial, and soil organic matter (SOM) pools. The model was parameterized and validated against C and N data from the Long-Term Inter-site Decomposition Experiment Team (LIDET; six litter types, 10 years of observations, and 13 sites across North America). The model simulates C and N losses from litterbags in the LIDET study with reasonable accuracy (C: R 2 =0.63; N: R 2 =0.29), which is comparable with simulations from the DAYCENT model that implicitly represents microbial activity (C: R 2 =0.67; N: R 2 =0.30). Subsequently, we evaluated equilibrium values of stocks (total soil C and N, microbial biomass C and N, inorganic N) and microbial process rates (soil heterotrophic respiration, N mineralization) simulated by MIMICS-CN across the 13 simulated LIDET sites against published observations from other continent-wide datasets. We found that MIMICS-CN produces equilibrium values in line with measured values, showing that the model generates plausible estimates of ecosystem soil biogeochemical dynamics across continental-scale gradients. MIMICS-CN provides a platform for coupling C and N projections in a microbially explicit model, but experiments still need to identify the physiological and stoichiometric characteristics of soil microbes, especially under environmental change scenarios.

58 GEOSCIENCES↗
DOE OSTI2023

CORPSE model with litter decomposition parameters derived from the LIDET dataset

This is a version of the CORPSE model (Carbon, Organisms, Rhizosphere and Protection in the Soil Environment, Sulman et al. 2014) that uses litter decomposition parameters derived from a modified Monte Carlo simulation using the LIDET litter decomposition dataset (Long-term Intersite Decomposition Experiment Team, Harmon 2013). The code also includes the Baseline parameters, and the eight other best parameter sets identified in a modified Monte Carlo simulation. 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:The folder "Input Files" contains one folder for each LIDET site with data necessary to run the model. Note that "(site)" in the filenames below indicates where the LIDET site code appears (see Table 1 for site codes). Data streams include: CORPSE_full_spinup_litter.csv, CORPSE_full_spinup_rhizo.csv, CORPSE_full_spinup_bulk.csv, litterbag_init_100g_6spp.csv: initial C and N (kg C or N/m2) pool values for each soil layer, the litterbag_init_100_6spp.csv file is for the litterbag layer and is the same file for all sites. All initial C and N 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 soilT (site) DOY274start.csv: Average daily soil temperature (oC) interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). soilT (site) DOY274start.csv: Average daily soil volumetric water content (VWC) scalar interpolated from previously calculated monthly values used in DayCent LIDET simulations (Bonan et al., 2013). litter production.csv: Average daily litter production values for each site, data sources listed in Table S3 of related publication. litter (site) CN.csv: C:N ratio for each species from LIDET dataset (Table 2, Harmon 2013). (site).csv: Table indicating number of observations for each species decomposed at each site. Instructions: Save the model code ("CORPSE_LIDET.R") and "Input Files" folder in the same folder. Also make a folder for the model output (e.g., "results_Baseline") in the same folder. Set the working directory (setwd) in the model code to the folder with the files saved in step #1. Select the parameter set to use for the litter and litterbag compartments, comment out all other parameter sets. Run code. Output will be saved in the folder made in step 1. Output destination can be changed as necessary in code section called "Running the model." Table 1 LIDET sites and site codes used in model files. Site Code - Site AND - H.J. Andrews Experimental Forest BNZ - Bonanza Creek Experimental Forest BSF - Blodgett Research Forest CDR - Cedar Creek Natural History Area CPR - Central Plains Experimental Range HBR - Hubbard Brook Experimental Forest HFR - Harvard Forest JUN - Juneau KBS - Kellogg Biological Station KNZ - Konza Prairie Research Natural Area NWT - Niwot Ridge/Green Lakes Valley OLY - Olympic National Park OLY Conifer forest SEV - Sevilleta National Wildlife Refuge SMR - Santa Margarita Ecological Reserve UFL - University of Florida VCR - Virginia Coast Reserve Table 2 LIDET species and species codes used in model files (6 common species). Species - Species Code Sugar maple (Acer saccharum) - ACSA Drypetes (Drypetes glauca) - DRGL Red pine (Pinus resinosa) - PIRE Chestnut oak (Quercus prinus) - QUPR Western redcedar (Thuja plicata) - THPL Wheat (Triticum aestivum) - TRAE References:Bonan, G. B., Hartman, M. D., Parton, W. J., & Wieder, W. R. (2013). Evaluating litter decomposition in earth system models with long-term litterbag experiments: an example using the Community Land Model version 4 (CLM4). Global Change Biology, 19(3), 957-974. https://doi.org/https://doi.org/10.1111/gcb.12031 Harmon, M. (2013). LTER Intersite Fine Litter Decomposition Experiment (LIDET), 1990 to 2002. Long-Term Ecological Research. Forest Science Data Bank, Corvallis, OR. [Data set]. Accessed http://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023. https://doi.org/10.6073/pasta/f35f56bea52d78b6a1ecf1952b4889c5. Sulman, B. N., Phillips, R. P., Oishi, A. C., Shevliakova, E., & Pacala, S. W. (2014). Microbe-driven turnover offsets mineral-mediated storage of soil carbon under elevated CO2. Nature Climate Change, 4, 1099 - 1102. https://doi.org/10.1038/nclimate2436

Juice, Stephanie↗
DOE OSTI2021

Active microbial biomass decreases, but microbial growth potential remains similar across soil depth profiles under deeply-vs. shallow-rooted plants

Climate-smart land management practices that replace shallow-rooted annual crop systems with deeply-rooted perennial plants can contribute to soil carbon sequestration. However, deep soil carbon accrual may be influenced by active microbial biomass and their capacity to assimilate fresh carbon at depth. Incorporating active microbial biomass, dormancy, and growth in microbially-explicit models can improve our ability to predict soil's capacity to store carbon. But, so far, the microbial parameters that are needed for such modeling are poorly constrained, especially in deep soil layers. Here, we used a lab incubation experiment and growth kinetics model to estimate how microbial parameters vary along 240 cm of soil depth in profiles under shallow- (soy) and deeply-rooted (switchgrass) plants 11 years after plant cover conversion. We also assessed resource origin and availability (total organic carbon, 14 C, extractable organic carbon, specific UV absorbance of K 2 SO 4 extractable organic C, total nitrogen, total dissolved nitrogen) along the soil profiles to examine associations between soil chemical and biological parameters. Even though root biomass was greater and rooting depth was deeper under switchgrass than soy, resource availability and microbial growth parameters were generally similar between vegetation types. Instead, depth significantly influenced soil chemical and biological parameters. For example, resource availability and total and relative active microbial biomass decreased with soil depth. Decreases in the relative active microbial biomass coincided with increased lag time (response time to external carbon inputs) along the soil profiles. Even at a depth of 210–240 cm, microbial communities were activated to grow by added resources within a day. Maximum specific growth rate decreased to a depth of 90 cm and then remained consistent in deeper layers. Our findings show that >10 years of vegetation and rooting depth changes may not be long enough to alter microbial growth parameters, and suggest that at least a portion of the microbial community in deep soils can grow rapidly in response to added resources. Our study determined microbial growth parameters that can be used in microbially-explicit models to simulate carbon dynamics in deep soil layers.

14C↗
DOE OSTI2021

Leveraging observed soil heterotrophic respiration fluxes as a novel constraint on global-scale models

Microbially-explicit models may improve understanding and projections of carbon dynamics in response to future climate change, but their fidelity in simulating global-scale soil heterotrophic respiration (RH), a stringent test for soil biogeochemical models, has never been evaluated. We used statistical global RH products, as well as 7,821 daily site-scale RH measurements, to evaluate the spatio-temporal performance of one first-order decay model (CASA-CNP) and two microbially-explicit biogeochemical models (CORPSE and MIMICS) that were forced by two different climate datasets. CORPSE and MIMICS did not provide any measurable performance improvement; instead, the models were highly sensitive to the meteorological input data used to drive them. Spatial RH variability was generally well simulated except in the northern middle latitudes (~50°N) and arid regions; models captured the seasonal variability of RH well, but showed more divergence in tropic and arctic regions. Our results demonstrate that the next generation of biogeochemical models shows promise, but also needs to be improved for realistic spatio-temporal variability of RH. Finally, we emphasize the importance of net primary production, soil moisture, and soil temperature inputs, and that jointly evaluating soil models for their spatial (global scale) and temporal (site scale) performance provides crucial benchmarks for improving biogeochemical models.

Jian, Jinshi↗
DOE OSTI2022

Multi-modelling predictions show high uncertainty of required carbon input changes to reach a 4‰ target

Soils store vast amounts of carbon (C) on land, and increasing soil organic carbon (SOC) stocks in already managed soils such as croplands may be one way to remove C from the atmosphere, thereby limiting subsequent warming. The main objective of this study was to estimate the amount of additional C input needed to annually increase SOC stocks by 4‰ at 16 long-term agricultural experiments in Europe, including exogenous organic matter (EOM) additions. We used an ensemble of six SOC models and ran them under two configurations: (1) with default parametrization and (2) with parameters calibrated site-by-site to fit the evolution of SOC stocks in the control treatments (without EOM). We compared model simulations and analysed the factors generating variability across models. The calibrated ensemble was able to reproduce the SOC stock evolution in the unfertilised control treatments. We found that, on average, the experimental sites needed an additional 1.5 ± 1.2 Mg C ha -1 year -1 to increase SOC stocks by 4‰ per year over 30 years, compared to the C input in the control treatments (multi-model median ± median standard deviation across sites). That is, a 119% increase compared to the control. While mean annual temperature, initial SOC stocks and initial C input had a significant effect on the variability of the predicted C input in the default configuration (i.e., the relative standard deviation of the predicted C input from the mean), only water-related variables (i.e., mean annual precipitation and potential evapotranspiration) explained the divergence between models when calibrated. Here, our work highlights the challenge of increasing SOC stocks in agriculture and accentuates the need to increasingly lean on multi-model ensembles when predicting SOC stock trends and related processes. To increase the reliability of SOC models under future climate change, we suggest model developers to better constrain the effect of water-related variables on SOC decomposition.

4 per 1000 initiative↗
DOE OSTI2022

Optimizing process-based models to predict current and future soil organic carbon stocks at high-resolution

From hillslope to small catchment scales (< 50 km 2 ), soil carbon management and mitigation policies rely on estimates and projections of soil organic carbon (SOC) stocks. Here we apply a process-based modeling approach that parameterizes the MIcrobial-MIneral Carbon Stabilization (MIMICS) model with SOC measurements and remotely sensed environmental data from the Reynolds Creek Experimental Watershed in SW Idaho, USA. Calibrating model parameters reduced error between simulated and observed SOC stocks by 25%, relative to the initial parameter estimates and better captured local gradients in climate and productivity. The calibrated parameter ensemble was used to produce spatially continuous, high-resolution (10 m 2 ) estimates of stocks and associated uncertainties of litter, microbial biomass, particulate, and protected SOC pools across the complex landscape. Here, subsequent projections of SOC response to idealized environmental disturbances illustrate the spatial complexity of potential SOC vulnerabilities across the watershed. Parametric uncertainty generated physicochemically protected soil C stocks that varied by a mean factor of 4.4 × across individual locations in the watershed and a – 14.9 to + 20.4% range in potential SOC stock response to idealized disturbances, illustrating the need for additional measurements of soil carbon fractions and their turnover time to improve confidence in the MIMICS simulations of SOC dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Sticky roots--implications of widespread, cryptic, viral infection of plants in natural and managed ecosystems for soil carbon processing in the rhizosphere

Plants strongly influence soil properties through rhizodeposition, in which exudates diffuse from roots, additional secretions are actively released, and root cells are sloughed into the soil. This contribution by plants of carbon compounds belowground is at the core of soil health, water holding capacity, and the soil carbon storage that pulls carbon dioxide out of the atmosphere. Once in soil, organic matter can bind with minerals such as iron hydroxides, where it can be protected from microbial attack for millenia, preserving very large terrestrial soil carbon pools. However, those same compounds contributed by roots to soil may also destabilize the long-term protective associations of SOM with minerals, making that soil organic matter (SOM) more vulnerable to microbial attack and decomposition. Plant roots thus influence both the buildup and breakdown of soil carbon pools. DOE’s E3SM Land Model (ELM) includes a representation of soil carbon storage on minerals, but the potential vulnerability of SOM–mineral associations to effects of rhizodeposition is not yet represented in ELM. To begin testing for this effect of rhizodeposition on soil carbon storage and decomposition, we worked to develop a novel approach during this TES Exploratory project DE-SC0019142 – we harnessed the power of plant viral infection. We examined whether plant virus infection can serve as a tool to intensify rhizodeposition at the root surface, and therefore possibly intensify mobilization of SOM from minerals making it visible to our analytical techniques. Viral infection is widespread in terrestrial ecosystems; 25-70% of plants have virus infection, yet the influence of such infection on root traits and terrestrial soil carbon dynamics remains largely unexplored. We used two plant hosts: the annual Avena sativa (oats) and the genetically tractable, model grass Brachypodium distachyon. These grasses were infected with the broad host range virus Barley Yellow Dwarf Virus (BYDV) via aphids (Rhopalosiphum padi). BYDV infects at least 150 grass species in agricultural and natural ecosystems, and in previous experiments, oats infected with BYDV had roots that were very sticky to the touch, strongly suggesting that infection altered rhizodeposition. We developed this new experimental approach mostly in a one virus (Barley Yellow Dwarf Virus)–one plant (Avena sativa) system. (Several effects of infection in a Brachypodium-BYDV system were similar in nature to effects on Avena sativa, but were more variable.) In the BYDV-Avena system, we developed protocols for consistently infecting target plants (and avoiding infection of control plants) using aphid caging on leaves. We measured that infected plants exhibited reduced photosynthesis, plant (including root) biomass, and root:shoot ratio, as well as simplified root system architecture. We established procedures for sampling the organic compounds carried specifically in phloem (vascular tissue) of leaves and roots, using aphid stylectomy. We used FTICR-MS, Orbitrap GC-MS, and LC-MS/MS to analyze organic compounds in phloem, liquid around roots of plants grown hydroponically, and pore water around roots in soil, and found differences in the compounds in solution bathing roots when infected and uninfected plants were grown hydroponically. Finally, we synthesized isotopically-labeled mineral–organic matter (MAOM) associations in the lab and developed assays using them in solution and in soil. Assays quantified the extent and rate of mineralization of labeled MAOM that was mobilized by functionally distinct rhizodeposits and then attacked by microbes. Two mechanisms for MAOM mobilization emerged, with distinct dynamics. During “direct” mobilization, rhizodeposits such as the strong ligand oxalic acid could drive rapid dissolution of minerals, mobilizing MAOM. During “indirect” mobilization, rhizodeposits such as the simple sugar glucose did not attack minerals directly but instead intensified microbial activity, which led to mobilization via changes in e.g. pH, Eh, and microbial metabolite production (Li et al. 2021). Mechanistic understanding derived from these data and our ongoing experiments using these techniques will inform future development of ELM. Plant roots not only contribute newly fixed organic compounds to soils, but also root activities can drive mineralization of the carbon and nutrients mobilized off minerals via “indirect” or “direct” mechanisms. Using viral infection as a new tool, ongoing combined experimentation and modeling will explore the strength and larger-scale significance of the cascade of processes from rhizodeposition to MAOM mobilization for soil carbon storage and nutrient cycling in terrestrial ecosystems. And if viral infection leads quite generally to “sticky roots”, our perception of the potential importance of prevalent virus infection in terrestrial landscapes will be transformed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Advancing quantitative models of soil microbiology, ecology, and biochemistry

Rapid advances in soil biology are increasingly reflected in mathematical models. From molecular to global scales, models contribute to fundamental understanding and prediction of critical soil processes, such as carbon sequestration. Recent models explicitly incorporate microbial control over soil carbon and nutrient cycling, an approach that offers greater biological realism but raises new challenges in model parameterization and stability. Although they have proliferated at community and ecosystem scales, microbial-explicit models have yet to be incorporated into most Earth system models. Applying soil models across scales requires integration with data to support parameterization and validation using approaches such as Bayesian data assimilation. In conclusion, to catalyze future progress, we recommend breaking down barriers between modeling and empirical disciplines while broadening access to the computational infrastructure that supports soil biological modeling.

Xie, Wally↗
DOE OSTI2020

Dynamics of Fungal and Bacterial Biomass Carbon in Natural Ecosystems: Site-level Applications of the CLM-Microbe Model

Explicitly representing microbial processes has been recognized as a key improvement to Earth system models for the realistic projections of soil carbon (C) and climate dynamics. The CLM-Microbe model builds upon the CLM4.5 and explicitly represents two major soil microbial groups, fungi and bacteria. Based on the compiled time-series data of fungal (FBC) and bacterial (BBC) biomass C from nine biomes, we parameterized and validated the CLM-Microbe model, and further conducted sensitivity analysis and uncertainty analysis for simulating C cycling. The model performance was evaluated with mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2) for relative change in FBC and BBC. The CLM-Microbe model is able to reasonably capture the seasonal dynamics of FBC and BBC across biomes, particularly for tropical/subtropical forest, temperate broadleaf forest, and grassland, with MAE < 0.49 for FBC and <0.36 for BBC and RMSE <0.52 FBC and <0.39 for BBC, while R2 values are relatively smaller in some biomes (e.g., shrub) due to small sample sizes. We found good consistencies between simulated and observed FBC (R 2 =0.70, P<0.001) and BBC (R 2 =0.26, P<0.05) on average across biomes, but the model is not able to fully capture the large variation in observed FBC and BBC. Sensitivity analysis shows the most critical parameters are turnover rate, carbon-to-nitrogen ratio of fungi and bacteria, and microbial assimilation efficiency. This study confirms that the explicit representation of soil microbial mechanisms enhances model performance in simulating C variables such as heterotrophic respiration and soil organic C density. The further application of the CLM-Microbe model would deepen our understanding of microbial contributions to the global C cycle.

54 ENVIRONMENTAL 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↗
DOE OSTI2024

Model uncertainty obscures major driver of soil carbon

Understanding the formation and stabilization mechanisms of soil organic carbon (SOC) is important for managing land carbon (C) and mitigating climate change. Tao et al. reported that microbial C use efficiency (CUE) is the primary determinant of global SOC storage and that the relative impact of plant C inputs on SOC is minor. Although soil microbes undoubtedly play an important role in SOC cycling, we are concerned about the robustness of the approach taken by Tao et al. In conclusion, the potential biases in their analyses may lead to misleading, model-dependent results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A new bioenergy model that simulates the impacts of plant‐microbial interactions, soil carbon protection, and mechanistic tillage on soil carbon cycling

Abstract Advancing our predictive understanding of bioenergy systems is critical to design decision tools that can inform which feedstock to plant, where to plant it, and how to manage its production to provide both energy and ecosystem carbon (C) benefits. Here, we lay the foundation for that advancement by integrating recent developments in the science of belowground processes in shaping the C cycle into a new bioenergy model, FUN‐BioCROP (Fixation and Uptake of Nitrogen‐Bioenergy Carbon, Rhizosphere, Organisms, and Protection). We show that FUN‐BioCROP can approximate the historical trajectory of soil C dynamics as natural ecosystems were successively converted into intensive agriculture and bioenergy systems. This ability relies in part on a novel tillage representation that mechanistically models tillage as a process that increases microbial access to C. Importantly, the impacts of tillage and feedstock choice also influence FUN‐BioCROP simulations of warming responses with no‐till perennial feedstocks, miscanthus, and switchgrass, having more C that is unprotected and susceptible to warming than tilled annual feedstocks like corn–corn–soybean. However, this susceptibility to warming is balanced by a greater potential for increases in belowground C allocation to enhance soil C stocks in perennial systems. Collectively, our model results highlight the importance of belowground processes in evaluating the ecosystem C benefits of bioenergy production.

09 BIOMASS FUELS↗
DOE OSTI2019

Plant roots stimulate the decomposition of complex, but not simple, soil carbon

Roots release carbon into soil and can alleviate energy limitation of microbial organic matter decomposition. We know little about the effects of roots on microbial decomposition of different organic matter substrates, despite the importance for soil carbon stocks and turnover. Through implementing root-microbe interactions, the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model was previously shown to represent dynamics of total soil carbon in temperate forest field experiments. However, the model permits alternative hypotheses concerning microbial-substrate affinity. We investigated how root inputs affect decomposition of soil organic carbon (SOC) with variable decomposability. We simulated SOC stocks in CORPSE and compared microbial degradation of two substrates types with varying root-microbe interactions under two alternative hypotheses that varied in microbial-substrate affinity. We compared our modeled hypotheses to a forest field experiment where we quantified decomposition of isotopically-labeled starch and leaf tissues in soils with manipulated root access to microbes. Here, we tested the hypothesis that decomposition of leaves would be more sensitive to root inputs than decomposition of starch, corresponding to the alternative model hypothesis. In the field study, leaf decomposition increased with root density while starch decomposition was unchanged by root density. Microbial biomass and enzyme activity consistently increased with root inputs in CORPSE and the field study. Our field experiment supported the CORPSE simulations with high microbial-substrate affinity. Roots stimulated microbial growth and enzyme production, which increased degradation of more complex substrates such as leaf tissues. Substrates that were easily decomposed, such as starch, may already be degrading at a maximum rate in the absence of rhizosphere influence because their decomposition rate was unchanged by root inputs. We found that the degree to which roots stimulate microbial decomposition depends on the substrate being decomposed, and that root-microbe interactions influenced SOC stocks in both our model and field experiment. Environmental changes that alter root-microbe interactions could, therefore, alter soil C stocks and biogeochemical cycling, and models of these interactions should incorporate differential influence of rhizosphere inputs on different substrates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Soil-related developments of the Biome-BGCMuSo v6.2 terrestrial ecosystem model

Terrestrial biogeochemical models are essential tools to quantify climate–carbon cycle feedback and plant–soil relations from local to global scale. In this study, a theoretical basis is provided for the latest version of the Biome-BGCMuSo biogeochemical model (version 6.2). Biome-BGCMuSo is a branch of the original Biome-BGC model with a large number of developments and structural changes. Earlier model versions performed poorly in terms of soil water content (SWC) dynamics in different environments. Moreover, lack of detailed nitrogen cycle representation was a major limitation of the model. Since problems associated with these internal drivers might influence the final results and parameter estimation, additional structural improvements were necessary. In this paper the improved soil hydrology as well as the soil carbon and nitrogen cycle calculation methods are described in detail. Capabilities of the Biome-BGCMuSo v6.2 model are demonstrated via case studies focusing on soil hydrology, soil nitrogen cycle, and soil organic carbon content estimation. Soil-hydrology-related results are compared to observation data from an experimental lysimeter station. The results indicate improved performance for Biome-BGCMuSo v6.2 compared to v4.0 (explained variance increased from 0.121 to 0.8 for SWC and from 0.084 to 0.46 for soil evaporation; bias changed from -0.047 to -0.007 m 3 m -3 for SWC and from -0.68 to -0.2 mm d -1 for soil evaporation). Simulations related to nitrogen balance and soil CO 2 efflux were evaluated based on observations made in a long-term field experiment under crop rotation. The results indicated that the model is able to provide realistic nitrate content estimation for the topsoil. Soil nitrous oxide (N 2 O) efflux and soil respiration simulations were also realistic, with overall correspondence with the observations (for the N 2 O efflux simulation bias was between -0.13 and -0.1 mgNm -2 d -1 , and normalized root mean squared error (NRMSE) was 32.4 %–37.6 %; for CO 2 efflux simulations bias was 0.04–0.17 gCm -2 d -1 , while NRMSE was 34.1 %–40.1 %). Sensitivity analysis and optimization of the decomposition scheme are presented to support practical application of the model. The improved version of Biome-BGCMuSo has the ability to provide more realistic soil hydrology representation as well as nitrification and denitrification process estimation, which represents a major milestone.

58 GEOSCIENCES↗
DOE OSTI2024

Modeling microbial carbon fluxes and stocks in global soils from 1901 to 2016

Abstract. Soil microbes play a crucial role in the carbon (C) cycle; however, they have been overlooked in predicting the terrestrial C cycle. We applied a microbial-explicit Earth system model – the Community Land Model-Microbe (CLM-Microbe) – to investigate the dynamics of soil microbes during 1901 to 2016. The CLM-Microbe model was able to reproduce the variations of gross (GPP) and net (NPP) primary productivity, heterotrophic (HR) and soil (SR) respiration, microbial (MBC) biomass C in fungi (FBC) and bacteria (BBC) in the top 30 cm and 1 m, and dissolved (DOC) and soil organic C (SOC) in the top 30 cm and 1 m during 1901–2016. During the study period, simulated C variables increased by approximately 12 PgC yr−1 for HR, 25 PgC yr−1 for SR, 1.0 PgC for FBC and 0.4 PgC for BBC in 0–30 cm, and 1.2 PgC for FBC and 0.7 PgC for BBC in 0–1 m. Increases in microbial C fluxes and pools were widely found, particularly at high latitudes and in equatorial regions, but we also observed their decreases in some grids. Overall, the area-weighted averages of HR, SR, FBC, and BBC in the top 1 m were significantly correlated with those of soil moisture and soil temperature in the top 1 m. These results suggested that microbial C fluxes and pools were jointly governed by vegetation C input and soil temperature and moisture. Our simulations revealed the spatial and temporal patterns of microbial C fluxes and pools in response to environmental change, laying the foundation for an improved understanding of soil microbial roles in the global terrestrial C cycle.

He, Liyuan↗