Search NASA⌕ Search

SEARCH · Search NASA

Results for “soil carbon:nitrogen”

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

Versatile soil gas concentration and isotope monitoring: optimization and integration of novel soil gas probes with online trace gas detection

Abstract. Gas concentrations and isotopic signatures can unveil microbial metabolisms and their responses to environmental changes in soil. Currently, few methods measure in situ soil trace gases such as the products of nitrogen and carbon cycling or volatile organic compounds (VOCs) that constrain microbial biochemical processes like nitrification, methanogenesis, respiration, and microbial communication. Versatile trace gas sampling systems that integrate soil probes with sensitive trace gas analyzers could fill this gap with in situ soil gas measurements that resolve spatial (centimeters) and temporal (minutes) patterns. We developed a system that integrates new porous and hydrophobic sintered polytetrafluoroethylene (sPTFE) diffusive soil gas probes that non-disruptively collect soil gas samples with a transfer system to direct gas from multiple probes to one or more central gas analyzer(s) such as laser and mass spectrometers. Here, we demonstrate the feasibility and versatility of this automated multiprobe system for soil gas measurements of isotopic ratios of nitrous oxide (δ18O, δ15N, and the 15N site preference of N2O), methane, carbon dioxide (δ13C), and VOCs. First, we used an inert silica matrix to challenge probe measurements under controlled gas conditions. By changing and controlling system flow parameters, including the probe flow rate, we optimized recovery of representative soil gas samples while reducing sampling artifacts on subsurface concentrations. Second, we used this system to provide a real-time window into the impact of environmental manipulation of irrigation and soil redox conditions on in situ N2O and VOC concentrations. Moreover, to reveal the dynamics in the stable isotope ratios of N2O (i.e., 14N14N16O, 14N15N16O, 15N14N16O, and 14N14N18O), we developed a new high-precision laser spectrometer with a reduced sample volume demand. Our integrated system – a tunable infrared laser direct absorption spectrometry (TILDAS) in parallel with Vocus proton transfer reaction mass spectrometry (PTR-MS), in line with sPTFE soil gas probes – successfully quantified isotopic signatures for N2O, CO2, and VOCs in real time as responses to changes in the dry–wetting cycle and redox conditions. Broadening the collection of trace gases that can be monitored in the subsurface is critical for monitoring biogeochemical cycles, ecosystem health, and management practices at scales relevant to the soil system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Stimulation of soil respiration by elevated CO 2 is enhanced under nitrogen limitation in a decade-long grassland study

Whether and how CO 2 and nitrogen (N) availability interact to influence carbon (C) cycling processes such as soil respiration remains a question of considerable uncertainty in projecting future C–climate feedbacks, which are strongly influenced by multiple global change drivers, including elevated atmospheric CO 2 concentrations (eCO 2 ) and increased N deposition. However, because decades of research on the responses of ecosystems to eCO 2 and N enrichment have been done largely independently, their interactive effects on soil respiratory CO 2 efflux remain unresolved. Here, we show that in a multifactor free-air CO 2 enrichment experiment, BioCON (Biodiversity, CO 2 , and N deposition) in Minnesota, the positive response of soil respiration to eCO 2 gradually strengthened at ambient (low) N supply but not enriched (high) N supply for the 12-y experimental period from 1998 to 2009. In contrast to earlier years, eCO 2 stimulated soil respiration twice as much at low than at high N supply from 2006 to 2009. In parallel, microbial C degradation genes were significantly boosted by eCO 2 at low but not high N supply. Further, incorporating those functional genes into a coupled C–N ecosystem model reduced model parameter uncertainty and improved the projections of the effects of different CO 2 and N levels on soil respiration. If our observed results generalize to other ecosystems, they imply widely positive effects of eCO 2 on soil respiration even in infertile systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Applying NIR and MIR spectroscopy for C and soil property prediction in northern cold-region ecosystems. Which approach works better?

Here, developing reliable predictions of soil attributes is necessary to understand northern cold-region climate-soil feedback. Calibration models using near-infrared (NIR) and mid-infrared (MIR) spectroscopy were developed to predict eight commonly measured soil properties for 119 soil samples representing a range of vegetation types, parent materials, and soil types spanning >23° of latitude from southeast Alaska to the Canadian high Arctic. In order to obtain a more accurate prediction, this study compared the performance of linear and non-linear calibration techniques, including lasso regression (Lasso), support vector machine (SVM), random forest (RF) and classic partial least squares (PLS) to predict different soil properties of these soils. Comparing the four models, we noticed that their performance was quite similar for MIR overall, while NIR achieved better results with a PLS model for our dataset. PLS coupled with MIR showed a better performance for soil parameters, such as total organic carbon (TOC), total nitrogen (TN), cation exchange capacity (CEC) and clay (R-squared of 0.9, 0.81, 0.80, and 0.84) when compared with NIR (R-squared of 0.85, 0.72, 0.81 and 0.68). However, using either MIR or NIR spectroscopy, PLS predictions for bulk density (BD) and sand content were not accurate. The variable importance analysis based on the PLS model successfully estimated the relative contribution of wavelengths influencing soil property predictions most. Overall, TOC, TN, CEC and clay mineral predictions are closely related to the occurrence of specific spectral bands in the MIR region. For example, wavelengths at 2978 and 1761 cm -1 for TOC and TN, as well as at 3064 cm -1 for CEC, were selected as the most influential predictor variables. We demonstrated that MIR spectroscopy is a powerful tool for more extensive monitoring in soils of the northern cold climate region; however, NIR could be utilized for rapid estimates when the highest accuracy is not essential.

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↗
NASA NTRS1977

What's new on the moon. II

Apollo missions and returned lunar samples have provided new information about the moon, the earth, the sun, and the universe. Analyses show that all the planets were formed by the rapid accumulation of small bodies into larger ones about 4.6 billion years ago. The existence of simple molecules formed by reactions between the soil particles and atoms of carbon, oxygen, and nitrogen that have come from the sun, suggests that the basic ingredients for life are common in the universe. The ratio of hydrogen to helium in the solar wind reaching the moon is found to be 20 to 1, whereas the earth-based measurements show this ratio to be 10 to 1. Although the moon does not have any magnetic field at the present, the analyses revealed the existence of such a field three billion years ago. The understanding of the reasons for the disappearance of this field is vital for understanding planetary magnetic fields. The determination of the chemical composition of the whole moon, the explanation of the moon's observed asymmetry, and the understanding of the nature of moon's interior will have to be achieved by future, possibly unmanned, missions.

French, B. M.↗
DOE OSTI2021

Simulating measurable ecosystem carbon and nitrogen dynamics with the mechanistically defined MEMS 2.0 model

Abstract. For decades, predominant soil biogeochemical models have used conceptual soil organic matter (SOM) pools and only simulated them to a shallow depth in soil. Efforts to overcome these limitations have prompted the development of the new generation SOM models, including MEMS 1.0, which represents measurable biophysical SOM fractions, over the entire root zone, and embodies recent understanding of the processes that govern SOM dynamics. Here we present the result of continued development of the MEMS model, version 2.0. MEMS 2.0 is a full ecosystem model with modules simulating plant growth with above- and belowground inputs, soil water and temperature by layer, decomposition of plant inputs and SOM, and mineralization and immobilization of nitrogen (N). The model simulates two commonly measured SOM pools – particulate and mineral-associated organic matter (POM and MAOM, respectively). We present results of calibration and validation of the model with several grassland sites in the US. MEMS 2.0 generally captured the soil carbon (C) stocks (R2 of 0.89 and 0.6 for calibration and validation, respectively) and their distributions between POM and MAOM throughout the entire soil profile. The simulated soil N matches measurements but with lower accuracy (R2 of 0.73 and 0.31 for calibration and validation of total N in SOM, respectively) than for soil C. Simulated soil water and temperature were compared with measurements, and the accuracy is comparable to the other commonly used models. The seasonal variation in gross primary production (GPP; R2 = 0.83), ecosystem respiration (ER; R2 = 0.89), net ecosystem exchange (NEE; R2 = 0.67), and evapotranspiration (ET; R2 = 0.71) was well captured by the model. We will further develop the model to represent forest and agricultural systems and improve it to incorporate new understanding of SOM decomposition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Long-term integrated soil-crop management improves soil microbial community structure to reduce GHG emission and increase yield

Integrated soil-crop management (ISCM) has been shown as an effective strategy to increase efficiency and yield while its soil microbial community structure and function remain unclear. We evaluated changes in soil physicochemical factors, bacterial community structure responses, and the contributions of soil properties and bacterial communities to summer maize-winter wheat yield and GHG emissions through an ISCM experiment [T1 (local smallholder farmers practice system), T2 (improved management system), T3 (high–yield production system), and T4 (optimized management system)], which could provide scientific guidance for sustainable development of soil in summer maize-winter wheat rotation system. The results showed that the optimized ISCM could improve the soil quality, which significantly changed the soil bacterial community structure to reduce GHG emissions and increase yield. The co-occurrence network density of T3 was increased significantly. The Acidobacteria (class) and OM190 (class) were enriched in T2 and T4. The Frankiales (order) and Gaiellales (order) were enriched in T3. However, the changes in different crop growth stages were different. At the wheat jointing stage and maize mature stage, T4 could enhance carbon-related functional groups, such as aromatic hydrocarbon degradation and hydrocarbon degradation, to increase the soil organic carbon content. And at the maize tasseling stage, T4 could enhance nitrogen-related functional groups. And soil bacteria structure and function indirectly affected annual yield and GHG emission. T2 and T4 exhibited a similar soil microbial community. However, the yield and nitrogen use efficiency of T2 were reduced compared to those of T4. The yield of T3 was the highest, but the GHG emission increased and soil pH and nitrogen use efficiency decreased significantly. Therefore, T4 was a suitable management system to improve soil quality and soil bacterial community structure and function to decrease GHG emissions and increase the yield of the summer maize-winter wheat rotation system.

Yu, Ningning↗
DOE OSTI2022

Dynamics of organic matter molecular composition under aerobic decomposition and their response to the nitrogen addition in grassland soils

Grassland soils store a substantial proportion of the global soil carbon (C) stock. The transformation of C in grassland soils with respect to chemical composition and persistence strongly regulate the predicted terrestrial-atmosphere C flux in global C biogeochemical cycling models. In addition, increasing atmospheric nitrogen (N) deposition alters C chemistry in grassland soils. However, there remains controversy about the importance of mineralogical versus biochemical preservation of soil C, as well as uncertainty regarding how grassland soil C chemistry responds to elevated N. This study used grassland soils with diverse soil organic matter (SOM) chemistries in an 8-month aerobic incubation experiment to evaluate whether the chemical composition of SOM converged across sites over time, and how SOM persistence responded to the N addition. This study demonstrates that over the course of incubation, the richness of labile compounds decreased in soils with less ferrihydrite content, whereas labile compounds were more persistent in ferrihydrite rich soils. In contrast, we found that the richness of more complex compounds increased over the incubation in most sites, independent of soil mineralogy. Moreover, we demonstrate the extent to which the diverse chemical composition of SOM converged among sites in response to microbial decomposition. N fertilization decreased soil respiration and inhibited the convergence of molecular composition across ecosystems by altering N demand for microbial metabolism and chemical interactions between minerals and organic molecules. This study provides original evidence that the decomposition and metabolism of labile organic molecules were largely regulated by soil mineralogy (physicochemical preservation), while the metabolism of more complex organic molecules was controlled by substrate complexity (biochemical preservation) independent to mineral-organic interactions. This study advanced our understanding of the dynamic biogeochemical cycling of C by unveiling that N addition dampened C respiration and diminished the convergence of SOM chemistry across diverse grassland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Terraforming - Making an earth of Mars

The possibility of creating a habitable environment on Mars via terraforming is discussed. The first step is to determine the amount, distribution, and chemical state of water, carbon dioxide, and nitrogen. The process of warming Mars and altering its atmosphere naturally divides into two steps: in the first step, the planet would be heated by a warm thick carbon dioxide atmosphere, while the second step would be to convert the atmospheric carbon dioxide and soil nitrates to the desired oxygen and nitrogen mixture. It is concluded that life will play a major role in any terraforming of Mars, and that terraforming will be a gradual evolutionary process duplicating the early evolution of life on earth.

Mckay, Christopher P.↗
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 OSTI2026

Soil nitrogen mineralization rates, nutrient stocks, stable isotopes, and water volumetric measurements across terrestrial-aquatic interfaces from three wetlands at the Tanglewood Biological Station, Alabama

This dataset supports a broader study investigating wetland hydrologic and biogeochemical responses to inundation events. Soil samples were collected across four sampling events along terrestrial-aquatic gradients at three wetland sites located within the Tanglewood Biological Station in Alabama from April 2024 to June 2025. The contents in this data package include soil in-situ nitrogen mineralization rates (measurements of net nitrification, net ammonification, and net mineralization), nutrient stocks (total carbon, total nitrogen, and organic matter), stable isotopes (carbon and nitrogen), and water volumetric measurements (water-filled pore space). Water level data related to each wetland location can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/2530253 (Kirker et al., 2024), related water geochemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3001967 (Forbes et al., 2025), and related surface water sediment chemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3377325 (Molina Serpas et al., 2026). In addition to this readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata and international generic sample numbers (IGSNs); (4) readme; (5) the field protocol; and (6) a subfolder with sample data. The sample data subfolder contains (1) net nitrification rate, (2) net ammonification rate, (3) areal net mineralization rate, (4) percent organic matter, (5) water-filled pore space, (6) total carbon content, (7) total nitrogen content, (8) stable carbon isotope (delta carbon-13), and (9) stable nitrogen isotope (delta nitrogen-15), and (10) methods codes. All files are .csv or .pdf.

13-C↗
DOE OSTI2025

Quantifying spatial and vertical variations in soil C:N relationships in permafrost-affected landscapes

Permafrost regions are experiencing rapid changes that affect carbon (C) and nitrogen (N) cycles, with implications for vegetation dynamics and gas exchanges with the atmosphere. Soil C:N ratio is a key indicator of organic matter quality, yet spatial estimates of N stocks and C:N ratios lag behind those for C. We used quantile regression forests to compare direct and indirect digital soil mapping approaches for predicting soil C:N ratios at 0–30, 30–60, and 60–100 cm depths across a latitudinal transect in Alaska. The indirect approach – deriving C:N from separately predicted C and N stocks – outperformed direct mapping for the surface layer (0–30 cm), while direct mapping was marginally better at greater depths. However, prediction accuracy decreased with depth for both methods. Temperature and topography were the most important predictors. Both approaches overestimated low and underestimated high C:N ratios, with direct mapping showing greater bias. Our results underscore the challenges of modeling C:N ratios in heterogeneous, data-sparse permafrost soils, but also suggest that indirect mapping holds promise if supported by more extensive datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Reviews and syntheses: Ironing out wrinkles in the soil phosphorus cycling paradigm

Soil phosphorus (P) management remains a critical challenge for agriculture worldwide, and yet we are still unable to predict soil P dynamics as confidently as that of carbon (C) or nitrogen (N). This is due to both the complexity of inorganic P (P i ) and organic P (P o ) cycling and the methodological constraints that have limited our ability to trace P dynamics in the soil–plant system. In this review, we describe the challenges for building parsimonious, accurate, and useful biogeochemical models that represent P dynamics and explore the potential of new techniques to usher P biogeochemistry research and modeling forward. We conclude that research efforts should focus on the following: (1) updating the McGill and Cole (1981) model of P o mineralization by clarifying the role and prevalence of biochemical and biological P o mineralization, which we suggest are not mutually exclusive and may co-occur along a continuum of P o substrate stoichiometry; (2) further understanding the dynamics of phytate, a six C compound that can regulate the poorly understood stoichiometry of soil P; (3) exploring the effects of C and P o saturation on P sorption and P o mineralization; and (4) resolving discrepancies between hypotheses about P cycling and the methods used to test these hypotheses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Managing Plant Microbiomes for Sustainable Biofuel Production

The development of environmentally sustainable, economical, and reliable sources of energy is one of the great challenges of the 21st century. Large-scale cultivation of cellulosic feedstock crops (henceforth, bioenergy crops) is considered one of the most promising renewable sources for liquid transportation fuels. However, the mandate to develop a viable cellulosic bioenergy industry is accompanied by an equally urgent mandate to deliver not only cheap reliable biomass but also ecosystem benefits, including efficient use of water, nitrogen, and phosphorous; restored soil health; and net negative carbon emissions. Thus, sustainable bioenergy crop production may involve new agricultural practices or feedstocks and should be reliable, cost effective, and minimal input, without displacing crops currently grown for food production on fertile land. In this editorial perspective for the Phytobiomes Journal Focus Issue on Phytobiomes of Bioenergy Crops and Agroecosystems, we consider the microbiomes associated with bioenergy crops, the effects beneficial microbes have on their hosts, and potential ecosystem impacts of these interactions. We also address outstanding questions, major advances, and emerging biotechnological strategies to design and manipulate bioenergy crop microbiomes. This approach could simultaneously increase crop yields and provide important ecosystem services for a sustainable energy future.

09 BIOMASS FUELS↗
DOE OSTI2025

Trade-offs among restored ecosystem functions are context-dependent in Mediterranean-type regions

Global biodiversity hotspots, including Mediterranean-type ecosystems worldwide, are highly threatened by global change that alters biodiversity, ecosystem functions, and services. Some restoration activities enhance ecosystem functions by reintroducing plant species based on known relationships between plant traits and ecosystem processes. Achieving multiple functions across different site conditions, however, requires understanding how abiotic factors like climate and soil, along with plant assemblages, influence ecosystem functions, including their trade-offs and synergies. We used the ModEST ecosystem simulation model, which integrates carbon, water, and nutrient processes with plant traits, to assess the relationships between restored plant assemblages and ecosystem functions in Mediterranean-type climates and soils. We investigated whether maximised carbon increment, water use efficiency, and nitrogen use efficiency, along with their trade-offs and synergies, varied across different abiotic contexts. Further, we asked whether assemblages that maximised functions varied across environments and among these functions. We found that maximised ecosystem carbon increment and nitrogen use efficiency occurred under moist, warm conditions, while water use efficiency peaked under drier conditions. Generally, the assemblage that maximised one function differed from those for other maximised functions. Synergies were rare, except between water and nitrogen use efficiencies in loam soils across most climates. Trade-offs among maximised functions were common, varying in strength with abiotic context and plant assemblages, and were more pronounced in sandy loam soils compared to clay-rich soils. Our findings suggest that due to variation in abiotic conditions within and across Mediterranean-type regions at the global scale, site-specific plant assemblages are required to maximise ecosystem functions. Thus, lessons from a single site cannot be transferred to another site, even where the same plant functional types are available for restoration. Our simulation results offer valuable insights into potential ecosystem performance under specific abiotic conditions following restoration with particular plant functional types, thereby informing local restoration efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Switchgrass Rhizosphere Metabolite Chemistry Driven by Nitrogen Availability

Plants and soil microorganisms interact closely in the rhizosphere where plants may exchange carbon (C) for functional benefits from the microbial community. For example, the bioenergy crop, switchgrass (Panicum virgatum) is thought to exchange root-exuded C for nitrogen (N) fixed by diazotrophs (free-living N-fixers). However, this interaction is not well characterized and it is not known how or if switchgrass responds to diazotrophs or their activity. To explore this question, we assessed rhizosphere metabolite chemistry of switchgrass grown in a hydroponic system under two N levels and under inoculated or uninoculated conditions. Plants were grown with the inoculum Azotobacter vinelandii DJ for 3 days before harvest. We found switchgrass root exudate chemistry to be driven by N availability. Total metabolite concentrations were generally greater under high N versus low N and unaffected by inoculation. Examination of rhizosphere chemical fingerprints indicates metabolite chemistry was also driven strongly by N availability with a greater relative abundance of carbohydrates under high N and greater relative abundance of organic acids under low N. We also found evidence of changes in rhizosphere chemical fingerprints by inoculation treatment. However, we found little evidence of N treatment and inoculation interaction effects that suggests this response is not directly mediated by N availability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

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

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

Juice, Stephanie↗