Search NASA⌕ Search

SEARCH · Search NASA

Results for “CORPSE 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 19 records
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗
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

Hicks Pries, Sulman, et al 2018. In situ incubation of 13C-labeled litter in three soil pits at Blodgett Forest, CA from 2013 to 2016

Even though over half of the world's soil organic carbon (SOC) is stored in subsoils (>20 cm deep), and the old ages of subsoil OC indicate its cycling differs from surface SOC, there are few studies examining in situ decomposition processes in deep soils. The purpose of this dataset is to elucidate these processes. In this study we added 13C-labeled fine roots to 15, 55, and 95 cm depths of a well-characterized coniferous forest Alfisol and monitored the amount of root-derived C remaining over 6, 12, and 30 months. We recovered the root-derived C in microbial phospholipid fatty acids (PLFAs) after 6 months and in coarse (>2 mm) particulate, fine (<2 mm) particulate, and dense, mineral-associated pools after 6, 12, and 30 months. Overall, root decomposition in the first 6 months was similar among all depths but significantly diverged at 30 months with faster decomposition at 15 cm than at 95 cm. There were more fungal and Gram negative-associated PLFAs at 15 cm than at 95 cm, and 13C analysis revealed those microbial groups preferred the added root carbon to native SOC. Mineral-associations were not the cause of slower decomposition at depth because similar amounts of applied root C was recovered in the dense fraction at all depths. The largest difference among depths was in the amount of root C recovered in the coarse particulate fraction, which was greater at 95 cm (50%) than at 15 cm (15%). Slower decomposition of the particulate pool at depth likely contributed to the increase in C:N ratios and depletion of δ13C values below 60 cm depth in our soil profiles. Simulations of these soils using the CORPSE model, which incorporates microbial priming effects and mineral stabilization of SOC, reproduced patterns of particulate and mineral-associated SOC over both time and depth and suggested that a lack of priming by root exudates at depth could account for the slower decomposition rate of particulate root material. Decomposition of deep particulate SOC may increase if root exudation or dissolved OC transport to depth increases.This dataset includes (1) characterization of the soil pits prior to litter addition, (2) characterization of how the litter changed after 6, 12, and 30 months, and (3) the CORPSE model structure, parameterization, and output.

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 OSTI2025

Testing mechanisms of how mycorrhizal associations affect forest soil carbon and nitrogen cycling (Final Technical Report)

Trees are in a symbiotic partnership with mycorrhizal fungi in which they provide the fungi with carbon from photosynthesis and the fungi provide the trees with nutrients and water. In temperate forests, the vast majority of trees form symbioses with one of two types of mycorrhizal fungi—arbuscular mycorrhizal (AM) fungi or ectomycorrhizal (EcM) fungi. These fungi differ in their morphology, hyphal length, and nutrient acquisition strategies. Many studies have found systematic differences in soil organic matter and nitrogen availability between forest stands dominated by AM-associating trees versus EcM-associating trees. For instance, there is a larger proportion of organic matter that is mineral-associated, more available nitrogen, and lower soil carbon to nitrogen ratios in AM forest stands relative to EcM forests stands. However, the mechanisms driving these patterns are not known, which complicates our ability to model soil organic matter dynamics in forested ecosystems. The main objective of this research was to understand the degree to which the observed differences in soil C and N dynamics between AM and EcM dominated forests are driven by tree traits like litter decomposability and root exudation versus mycorrhizal fungal nutrient acquisition strategies. We investigated these mechanisms using observations and targeted experiments and incorporated this knowledge into a process-based soil organic matter model. The observational studies compared the importance of leaf litter decomposability versus fungal identity on soil organic matter processes. We found that often fungal identity and traits were more important drivers of soil organic matter patterns than leaf litter decomposability. We ran two novel experiments: 1) a growth chamber experiment across four EcM and four AM tree species using a 13 C-labeled atmosphere to trace seedling-derived C into hyphae, the rhizosphere, and soil; and 2) an in situ decomposition experiment of six different 13 C and 15 N labeled litters that ranged in decomposability incubated across a gradient of EcM dominance at three sites that capture important variation in climate, soils, and forest species composition. In the first experiment, we found no significant differences of seedling mycorrhizal association on soil carbon sequestration over a growing season, but we did find that mycorrhizal association affected rhizodeposition with EcM-associating seedlings depositing more carbon in response to increased nitrogen availability. The decomposition experiment is still ongoing, but thus far, we have found slower litter decomposition in only one of three EcM-dominated forests which suggests that differences between AM- and EcM-dominated forests depend on the environmental context and identity of the EcM fungi. Lastly, we explicitly incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model creating Myco-CORPSE. By including the different nutrient acquisition strategies of AM and EcM fungi, we explored the conditions under which EcM fungi can slow decomposition rates and lead to greater soil organic carbon accumulation compared to AM fungi. We found that the effect of EcM fungi was highly context dependent and that EcM fungi decreased decomposition in colder forests with recalcitrant litter inputs and when they produced oxidases and necromass-degrading enzymes. Our research highlights the importance of fungal nutrient acquisition in driving soil organic matter patterns and the need to move beyond the AM-EcM dichotomy to consider the identity and traits of the specific fungi participating in the symbiosis. Overall, this research has resulted in six, peer-reviewed published papers in journals such as Global Change Biology, Ecology (2), Soil Biology and Biochemistry, and Ecosystems (2). There are at least two more papers in progress on this research including one that was recently submitted to Global Change Biology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment

It has been proposed that competition between ectomycorrhizal (ECM) fungi and free-living saprotrophs for resources like nitrogen (N) slows decomposition and increases the soil carbon storage in ECM ecosystems compared to arbuscular (AM) ecosystems. However, empirical evidence for the generality of such ECM effects is equivocal, and confounding mechanisms have been proposed that affect the magnitude and direction of ECM effects on soil carbon. Here, for this work, we conduct a theoretical modeling experiment, where we explicitly incorporate mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model. We use the model to explore the conditions under which ECM N acquisition processes can induce stronger saprotrophic N limitation and result in slower decomposition rates and greater soil organic carbon accumulation compared to AM processes. We found that the ECM fungi more strongly inhibited decomposition when litter inputs were N-depleted and relatively recalcitrant and when ECM fungi possessed a strong capacity to mine N from both recalcitrant soil organic matter and microbial necromass. Climate and seasonality also played a role as the ECM competition effect was strongest at low mean annual temperatures and when litterfall peaked seasonally. Priming effects driven by high root exudation rates in ECM-dominated systems could overwhelm the competition effect and reduce soil carbon under some circumstances. The ECM effect on decomposition in our simulations was highly context dependent. Based on our model results, we expect to see a strong ECM competition effect in temperate deciduous and boreal forests with relatively recalcitrant litter inputs, and with ECM fungi that produce oxidases and necromass-degrading enzymes. However, even a relatively strong ECM competition effect on decomposition only increased soil organic carbon accumulation by ~10%.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

TES: Modelling Microbes to Predict Post-Fire Carbon Cycling in the Boreal Forest across Burn Severities (Exploratory Proposal)

Wildfires in boreal ecosystems represent a globally significant ecological process that is expected to be sensitive to changing climate as well as forest management strategies, but is currently insufficiently understood and represented in models. We sought to determine whether linking belowground microbial community composition, size, and activity to aboveground properties of burn severity and plant community composition (upland jack pine and upland spruce) would allow us to better model post-fire soil CO 2 fluxes using the Carbon, Organisms, Rhizosphere and Protection in the Soil Environment (CORPSE) model. We used an integrated modelling experimental approach, with mechanistic laboratory experiments designed to inform models.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2021

Soil organic matter and plant carbon allocated to nitrogen acquisition simulated by the FUN-BioCROP model

This data package contains the model input, results, and validation data from Juice et al (citation below). The FUN-BioCROP model (Fixation and Uptake of Nitrogen- Bioenergy Carbon, Rhizosphere, Organisms, and Protection) advances the field of bioenergy modeling by integrating new empirical paradigms of the role of belowground processes in shaping coupled carbon (C) and nitrogen (N) cycles. It was developed by modifying the FUN-CORPSE model (Fixation and Uptake of Nitrogen- Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment, Sulman et al. 2017 Ecology Letters) for use in bioenergy systems by including mechanistic tillage, organic matter addition, nitrogen fertilization, harvest, and feedstock-specific parameters, and to be driven by DayCent plant productivity and biomass data.

biofuel sustainability↗
DOE OSTI2025

Modeling Mycorrhizal Carbon Costs in Temperate Forests: The Impacts of Functional Diversity and Global Change Factors

Mycorrhizal fungi form symbiotic relationships with most plant species, facilitating nutrient acquisition while consuming a significant fraction of the plant's photosynthetic carbon (C), which we define as the mycorrhizal C cost. Drivers of the mycorrhizal C cost, which is crucial for predicting environmental impacts on plant productivity, remain under-explored and difficult to quantify. Ecosystem models that incorporate mycorrhizae can offer insights into mycorrhizal C cost dynamics, but their predictions have rarely been validated against empirical data. Here, in this study, we used the Myco-CORPSE model, which explicitly simulates mycorrhizal processes alongside soil carbon and nitrogen cycling, to investigate the drivers of mycorrhizal C cost in temperate forests. Applying this model to over 1,800 forest inventory plots across the eastern United States, we found that the simulations matched published data, showing higher C allocation to ectomycorrhizal (ECM) fungi (16.0% of net primary production (NPP)) compared to arbuscular mycorrhizal (AM) fungi (5.8% of NPP). Further analysis showed that mixed forests, co-dominated by both AM and ECM trees, allocated less C to mycorrhizal fungi compared to forests dominated by either AM or ECM fungi alone, due to complementary nutrient acquisition strategies. Elevated Nitrogen (N) deposition and higher temperatures reduce mycorrhizal C costs, favoring AM strategies. Conversely, elevated CO 2 (eCO 2 ) increased plant N demand and mycorrhizal C costs, favoring ECM strategies that access organic N sources. These findings underscore the critical role of mycorrhizal functional diversity in plant nutrient acquisition and C dynamics, providing new insights into how mycorrhizal symbioses respond to global change.

Shao, Siya [Dartmouth College, Hanover, NH (United↗
DOE OSTI2024

Reparameterizing Litter Decomposition Using a Simplified Monte Carlo Method Improves Litter Decay Simulated by a Microbial Model and Alters Bioenergy Soil Carbon Estimates

Litter decomposition determines soil organic matter (SOM) formation and plant-available nutrient cycles. Therefore, accurate model representation of litter decomposition is critical to improving soil carbon (C) projections of bioenergy feedstocks. Soil C models that simulate microbial physiology (i.e., microbial models) are new to bioenergy agriculture, and their parameterization is often based on small datasets or manual calibration to reach benchmarks. Here, we reparameterized litter decomposition in a microbial soil C model (CORPSE - Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment) using the continental-scale Long-term Inter-site Decomposition Experiment Team (LIDET) dataset which documents decomposition across a range of litter qualities over a decade. Here, we conducted a simplified Monte Carlo simulation that constrained parameter values to reduce computational costs. The LIDET-derived parameters improved modeled C and nitrogen (N) remaining, decomposition rates, and litter mean residence times as compared to Baseline parameters. We applied the LIDET litter decomposition parameters to a microbial bioenergy model (Fixation and Uptake of Nitrogen – Bioenergy Carbon, Rhizosphere, Organisms, and Protection) to examine soil C estimates generated by Baseline and LIDET parameters. LIDET parameters increased estimated soil C in bioenergy feedstocks, with even greater increases under elevated plant inputs (i.e., by increasing residue, N fertilization). This was due to the integrated effects of plant litter quantity, quality, and agricultural practices (tillage, fertilization). Collectively, we developed a simple framework for using large-scale datasets to inform the parameterization of microbial models that impacts projections of soil C for bioenergy feedstocks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Root litter decomposition slows with soil depth

Even though over half of the world's soil organic carbon (SOC) is stored in subsoils (>20 cm deep), and the old ages of subsoil OC indicate its cycling differs from surface SOC, there are few studies examining in situ decomposition processes in deep soils. Here, we added 13C-labeled fine roots to 15, 55, and 95 cm depths of a well-characterized coniferous forest Alfisol and monitored the amount of root-derived C remaining over 6, 12, and 30 months. We recovered the root-derived C in microbial phospholipid fatty acids (PLFAs) after 6 months and in coarse (>2 mm) particulate, fine (<2 mm) particulate, and dense, mineral-associated pools after 6, 12, and 30 months. Overall, root decomposition in the first 6 months was similar among all depths but significantly diverged at 30 months with faster decomposition at 15 cm than at 95 cm. There were more fungal and Gram negative-associated PLFAs at 15 cm than at 95 cm, and13C analysis revealed those microbial groups preferred the added root carbon to native SOC. Mineral-associations were not the cause of slower decomposition at depth because similar amounts of applied root C was recovered in the dense fraction at all depths. The largest difference among depths was in the amount of root C recovered in the coarse particulate fraction, which was greater at 95 cm (50%) than at 15 cm (15%). Slower decomposition of the particulate pool at depth likely contributed to the increase in C:N ratios and depletion of δ13C values below 60 cm depth in our soil profiles. Simulations of these soils using the CORPSE model, which incorporates microbial priming effects and mineral stabilization of SOC, reproduced patterns of particulate and mineral-associated SOC over both time and depth and suggested that a lack of priming by root exudates at depth could account for the slower decomposition rate of particulate root material. Decomposition of deep particulate SOC may increase if root exudation or dissolved OC transport to depth increases.

13C↗
DOE OSTI2023

Soil Respiration Responses to Throughfall Exclusion Are Decoupled From Changes in Soil Moisture for Four Tropical Forests, Suggesting Processes for Ecosystem Models

Abstract Climatic drying is predicted for many tropical forests yet models remain poorly parameterized for these ecosystems, hampering predictions of forest‐climate interactions. We applied an integrated model–experiment approach, parameterizing an ecosystem model with tropical forest observational data and comparing model predictions to a field drying manipulation. We hypothesized that drying suppresses soil CO 2 fluxes (i.e., respiration) in relatively dry tropical forests but increases CO 2 fluxes in wetter tropical forests by alleviating anaerobiosis. We measured soil CO 2 fluxes during wet‐dry cycles from 2015 to 2022 in four Panamanian forests that vary in rainfall and soil fertility. Measured soil CO 2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in lower soil moisture optima for respiration than previously modeled. We then parameterized the model using field data and the new moisture‐respiration response functions. The updated model predicted increased soil CO 2 fluxes with drying in wetter and fertile forests and suppressed fluxes in drier, infertile forests. In contrast to model predictions, a chronic throughfall exclusion experiment initially suppressed soil respiration across forests, with sustained suppression for four years in the wettest forest only (−28% ± 4% during the dry season). In the fertile forest, drying eventually elevated CO 2 fluxes over this period (+75% ± 28% during the late wet season). The unexpected negative drying effect in the wettest, infertile forest could have resulted from reduced vertical flushing of nutrients into soils. Including hydro‐nutrient interactions in ecosystem models could improve predictions of tropical forest‐climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Carbon Organisms Rhizosphere and Protection in Soil Environment model script and input data for soil moisture-respiration responses in tropical forests

Objectives: Climatic drying is predicted for many tropical forests, yet models remain poorly parameterized for tropical forests, hampering predictions of forest-climate feedbacks. We applied an integrated model–experiment approach, parameterizing an ecosystem model Carbon Organisms Rhizosphere and Protection in the Soil Environment (CORPSE) with tropical forest observational data, and comparing model predictions with a field drying manipulation. We hypothesized that drying would suppress soil CO2 fluxes (i.e., respiration) in already-drier tropical forests, but increases CO2 fluxes in wetter tropical forests by alleviating anaerobiosis. We measured soil CO2 fluxes, soil moisture, soil temperature, and forest floor biomass during wet-dry cycles (2015 – 2022) in four Panamanian forests that vary in rainfall and soil fertility. We used the field data to parameterize and run tests in the model.Results: Measured CO2 fluxes declined in the dry season and peaked in the early wet season ahead of peak soil moisture, resulting in a lower soil moisture optimum for respiration than previously modeled. We used this data to parameterize the model, which then predicted increased soil CO2 fluxes in wetter and fertile forests with drying, and decreased fluxes in drier, infertile forests. In contrast to model predictions, a chronic throughfall exclusion experiment in the forests initially suppressed soil CO2 fluxes across forests, with sustained suppression after four years in the wettest forest only (-28 ± 4% during the dry season), but elevated soil CO2 fluxes in a fertile forest after four years (+75 ± 28% during the late wet season), as predicted by the model. The unexpected negative drying effect in the wettest, most infertile forest could have resulted from reduced vertical flushing of nutrients into soils. Including hydro-nutrient interactions in ecosystem models could improve predictions of tropical forest-climate feedbacks (results presented in Cusack et al. 2023). Datasets included: Code files:CORPSE_array.py: Defines the equations of the CORPSE modelCORPSE_solvers: Functions for running the CORPSE model using either iterative or ordinary differential equation (ODE) solversrun_Panama_sims.py: Read in datasets and run the model simulations for this studyInput data:PanamaGradientEcosystemChem_BT_CPools_20152016CO2_DC_20190615.xlsx: Plot characteristics used in running model simulationsLiCor compiled surface flux only to 2020_03 DC_20200825.xlsx: Surface gas exchange fluxes used in model-data comparisonsPARCHED litterfall data for Ben Sulman LD 20200902.xlsx: Litterfall data used to drive model simulationsInitialization data:state_500y_20190823.csv: Initial state of model pools based on previous spinup runsOutput data:Outputs/prev_moisture_response.csv: Simulations of multiple sites using original model moisture response function.Outputs/updated_moisture_response.csv: Simulations of multiple sites using updated model moisture response function.Outputs/dry15_prev_moisture_response.csv: Simulations with soil moisture reduced by 15%, using original moisture response function.Outputs/dry15_updated_moisture_response.csv: Simulations with soil moisture reduced by 15%, using updated moisture response function.Outputs/dry30_prev_moisture_response.csv: Simulations with soil moisture reduced by 30%, using original moisture response function.Outputs/dry30_updated_moisture_response.csv: Simulations with soil moisture reduced by 30%, using updated moisture response function.Outputs/latestart_prev_moisture_response.csv: Simulations with extended dry season, using original moisture response function.Outputs/latestart_updated_moisture_response.csv: Simulations with extended dry season, using updated moisture response function.Outputs/[site name]_oneyear.csv: One-year simulation for each site in expanded site list using original moisture response function.Outputs/[site name]_oneyear_dried.csv: One-year simulation for each site in expanded site list using original moisture response function, with soil moisture reduced by 25%.Outputs/[site name]_oneyear_updated_moisture_response.csv: One-year simulation for each site in expanded site list using updated moisture response function.Outputs/[site name]_oneyear_updated_moisture_response_dried.csv: One-year simulation for each site in expanded site list using updated moisture response function, with soil moisture reduced by 25%.Field plot location data:There is also a .kml file that includes coordinates for all 32 plots included in the study of four forests (n = 4 throughfall reduction and n = 4 control plots per site).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Myco-CORPSE simulations assessing mycorrhizal carbon allocation across U.S. forests and global change scenarios

Plants allocate a substantial portion of their fixed carbon belowground to mycorrhizal fungi in exchange for nutrients and other benefits. However, most current ecosystem models omit mycorrhizal processes, limiting our ability to predict plant–soil carbon dynamics under environmental change. To address this gap, we used a mycorrhiza-explicit soil biogeochemical model, Myco-CORPSE (Mycorrhizal Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment), to simulate tree carbon allocation to arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) fungi in temperate forests.The dataset includes outputs from two sets of model simulations:1. Perturbation experiments: Simulations across gradients of ECM dominance (0–100%), nitrogen deposition, soil temperature, and net primary productivity (NPP) to test how these factors affect mycorrhizal C allocation and nutrient cycling.2. FIA-based simulations: Model applications to over 1,800 U.S. forest sites using site-specific data from the U.S. Forest Inventory and Analysis (FIA) program, including vegetation composition, mycorrhizal type, climate, litter traits, soil properties, and N deposition.Model outputs include simulated mycorrhizal carbon allocation and related biogeochemical variables, such as soil and microbial carbon and nitrogen stocks. Data are provided in CSV format and organized by experiment type (in separate ZIP files). Python scripts for running simulations, plotting, and spatial mapping are also included and organized similarly. No proprietary software is required. These outputs support a peer-reviewed study and were used to generate figures and tables in the associated publication.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Benchmarking and Improving Microbial-explicit Soil Biogeochemistry Models

Soils represent the largest terrestrial carbon pool on Earth. Yet emerging theories regarding stabilization of soil organic matter remain poorly represented in global-scale models; thus, underestimating the true uncertainty associated with potential terrestrial carbon cycle – climate feedbacks. With this project we created a testbed to compare and evaluate soil biogeochemical models (W. Wieder at the Univ. of CO, Boulder), data synthesis (R. Phillips, Indiana Univ.) and parameter estimation (S Allison, Univ. CA Irvine) to quantify uncertainty associated with potential soil C changes. This project resulted in 14 Presentations, 14 Publications and helped support the training of three graduate students. A significant accomplishment with this work was the creation of a soil biogeochemical testbed, which forces three soil carbon models (CASA, MIMICS, and CORPSE) with common environmental drivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Leveraging the signature of heterotrophic respiration on atmospheric CO 2 for model benchmarking

Spatial and temporal variations in atmospheric carbon dioxide (CO 2 ) reflect large-scale net carbon exchange between the atmosphere and terrestrial ecosystems. Soil heterotrophic respiration (HR) is one of the component fluxes that drive this net exchange, but, given observational limitations, it is difficult to quantify this flux or to evaluate global-scale model simulations thereof. Here, we show that atmospheric CO 2 can provide a useful constraint on large-scale patterns of soil heterotrophic respiration. We analyze three soil model configurations (CASA-CNP, MIMICS, and CORPSE) that simulate HR fluxes within a biogeochemical test bed that provides each model with identical net primary productivity (NPP) and climate forcings. We subsequently quantify the effects of variation in simulated terrestrial carbon fluxes (NPP and HR from the three soil test-bed models) on atmospheric CO 2 distributions using a three-dimensional atmospheric tracer transport model. Our results show that atmospheric CO 2 observations can be used to identify deficiencies in model simulations of the seasonal cycle and interannual variability in HR relative to NPP. In particular, the two models that explicitly simulated microbial processes (MIMICS and CORPSE) were more variable than observations at interannual timescales and showed a stronger-than-observed temperature sensitivity. Our results prompt future research directions to use atmospheric CO 2 , in combination with additional constraints on terrestrial productivity or soil carbon stocks, for evaluating HR fluxes.

54 ENVIRONMENTAL SCIENCES↗