Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fine root productivity”

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.

48 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Rooting for function: community‐level fine‐root traits relate to many ecosystem functions

Humans are driving biodiversity change, which also alters community functional traits. However, how changes in the functional traits of the community alter ecosystem functions—especially belowground—remains an important gap in our understanding of the consequences of biodiversity change. We test hypotheses for how the root traits of the root economics space (composed of the collaboration and conservation gradients) are associated with proxies for ecosystem functioning across grassland and forest ecosystems in both observational and experimental datasets from 810 plant communities. First, we assessed whether community-weighted means of the root economics space traits adhered to the same trade-offs as species-level root traits. Then, we examined the relationships between community-weighted mean root traits and aboveground biomass production, root standing biomass, soil fauna biomass, soil microbial biomass, decomposition of standard and plot-specific material, ammonification, nitrification, phosphatase activity, and drought resistance. We found evidence for a community collaboration gradient but not for a community conservation gradient. Yet, links between community root traits and ecosystem functions were more common than we expected, especially for aboveground biomass, microbial biomass, and decomposition. These findings suggest that changes in species composition, which alter root trait means, will in turn affect critical ecosystem functions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Transgenic Mixed‐Linkage‐Glucan Enhancement Affects Root Characteristics and Decomposition in Soils of Contrasting Vegetation History

ABSTRACT Development of transgenic bioenergy sorghum [ Sorghum bicolor (L.) Moench] with increased contents of mixed‐linkage (1,3;1,4)‐β‐glucan (MLG) is an important step towards enhancing quality of bioenergy feedstocks. Since MLG‐enhancement leads to greater biomass digestibility, our overarching hypothesis is that root residues of MLG‐enhanced plants may be more readily decomposed in the soil, potentially creating new opportunities for optimizing soil carbon (C) sequestration, nutrient cycling, and overall agricultural sustainability. The study examined morphological, chemical, and enzymatic characteristics of fine and coarse roots of four bioenergy sorghum genotypes. Then, we incubated the roots within soils with contrasting vegetation histories while measuring C mineralization, microbial biomass C (MBC), and activity of hydrolytic enzymes and calculating vector length and vector angle enzymatic stoichiometry parameters. The results indicated that MLG‐enhancing transformations increased root total nitrogen (N) contents, decreased C/N ratios, and were associated with higher MLG concentrations in fine than in coarse roots. Incubations with transgenic roots led to 16%–38% higher MBC and 19%–41% lower microbial metabolic quotient (qCO 2 ). While enzyme activity differed markedly among the studied genotypes, it did not directly respond to MLG levels in root tissues. The increase in MBC without concurrent increases in C mineralization or hydrolytic enzyme activities in transgenic genotypes suggests that MLG enhancement promoted microbial anabolic retention of root‐derived C rather than stimulating catabolic decomposition. Enzymatic vector results indicated that these parameters reflect a variety of drivers behind microbial enzyme production, including availability of specific substrates, such as MLG here, and/or deficiency in specific nutrients, such as phosphorus (P). The study confirms the positive impacts from the roots of engineered MLG‐enhanced bioenergy plants on soil microbial activity and highlights the interactive influences on the MLG‐enhancement effects from root size and inherent soil properties.

Mahmoodabadi, Majid [Department of Plant, Soil, an↗
DOE OSTI2023

Belowground Biomass and Root:Shoot Ratios of Three Willow Cultivars at Two Sites

Belowground biomass is an important but less studied component of energy crop systems that is essential in understanding the greenhouse gas benefits of these systems. In this study, a complete above- and belowground biomass inventory (foliage, stems, stools, coarse and fine roots) was performed on three cultivars of short-rotation willow biomass crops at two sites. Mixed models were used to analyze the proportion of biomass allocated to each component and the ratios between different components. The root:shoot (R:S) ratio, defined here as the stable unharvested biomass (stool and coarse roots) divided by the shoot biomass, averaged 0.63 (SE: +0.04). Though the portion of the plant where the willows distributed their belowground biomass varied, the R:S ratio was not significantly different across sites (p = 0.8970), cultivars (p = 0.2834), nor in the site × cultivar interaction (p = 0.8481). These results may be associated with the consistently good growth across sites and limited differences in site conditions. However, the R:S ratios were affected by the overall productivity of the stand (p = 0.0978), with higher producing stands having moderately lower ratios. This information on biomass allocation between components is essential for understanding and estimating the carbon balance of these systems and breeding and selection programs.

09 BIOMASS FUELS↗
DOE OSTI2019

Embryogenic cell suspensions for high-capacity genetic transformation and regeneration of switchgrass (Panicum virgatum L.)

Abstract Background Switchgrass ( Panicum virgatum L.), a North American prairie grassland species, is a potential lignocellulosic biofuel feedstock owing to its wide adaptability and biomass production. Production and genetic manipulation of switchgrass should be useful to improve its biomass composition and production for bioenergy applications. The goal of this project was to develop a high-throughput stable switchgrass transformation method using Agrobacterium tumefaciens with subsequent plant regeneration. Results Regenerable embryogenic cell suspension cultures were established from friable type II callus-derived inflorescences using two genotypes selected from the synthetic switchgrass variety ‘Performer’ tissue culture lines 32 and 605. The cell suspension cultures were composed of a heterogeneous fine mixture culture of single cells and aggregates. Agrobacterium tumefaciens strain GV3101 was optimum to transfer into cells the pANIC-10A vector with a hygromycin-selectable marker gene and a pporRFP orange fluorescent protein marker gene at an 85% transformation efficiency. Liquid cultures gave rise to embryogenic callus and then shoots, of which up to 94% formed roots. The resulting transgenic plants were phenotypically indistinguishable from the non-transgenic parent lines. Conclusion The new cell suspension-based protocol enables high-throughput Agrobacterium -mediated transformation and regeneration of switchgrass in which plants are recovered within 6–7 months from culture establishment.

Ondzighi-Assoume, Christine A.↗
DOE OSTI2023

Can upscaling ground nadir SIF to eddy covariance footprint improve the relationship between SIF and GPP in croplands?

Ground solar-induced chlorophyll fluorescence (SIF) is important for the mechanistic understanding of the dynamics of vegetation gross primary production (GPP) at fine spatiotemporal scales. However, eddy covariance (EC) observations generally cover larger footprint areas than ground SIF observations (a bare fiber with nadir), and this footprint mismatch between nadir SIF and GPP could complicate the canopy SIF-GPP relationships. Here, we upscaled nadir SIF observations to EC footprint and investigated the change in SIF-GPP relationships after the upscaling in cropland. We included 13 site-years data in our study, with seven site-years corn, four site-years soybeans, and two site-years miscanthus, all located in the US Corn Belt. All sites’ crop nadir SIF observations collected from the automated FluoSpec2 system (a hemispheric-nadir system) were upscaled to the GPP footprint-based SIF using vegetation indices (VIs) calculated from high spatiotemporal satellite reflectance data. We found that SIF-GPP relationships were not substantially changed after upscaling nadir SIF to GPP footprint at our crop sites planted with corn, soybean, and miscanthus, with R 2 change after the upscaling ranging from -0.007 to 0.051 and root mean square error (RMSE) difference from -0.658 to 0.095 umol m –2 s –1 relative to original nadir SIF-GPP relationships across all the site-years. The variation of the SIF-GPP relationship within each species across different site-years was similar between the original nadir SIF and upscaled SIF. Different VIs, EC footprint models, and satellite data led to marginal differences in the SIF-GPP relationships when upscaling nadir SIF to EC footprint. Furthermore, our study provided a methodological framework to correct this spatial mismatch between ground nadir SIF and GPP observations for croplands and potentially for other ecosystems. Our results also demonstrated that the spatial mismatch between ground nadir SIF and GPP might not significantly affect the SIF-GPP relationship in cropland that are largely homogeneous.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Root Foraging Alters Global Patterns of Ecosystem Legacy From Climate Perturbations

The response of terrestrial ecosystems to climate perturbations typically persist longer than the timescale of the forcing, a phenomenon broadly referred to as legacy. Understanding the strength of legacy is critical for predicting ecosystem sensitivity to climate extremes and the extent that persistent changes in surface-atmosphere exchange might feedback onto climate. The cause of ecosystem legacy has been associated with myriad factors such as changes in aboveground biomass, however, few studies have tested how changes in the depth distribution of fine roots in response to perturbation might alter an ecosystem's recovery time. We explore this question using an Earth System Model that includes a dynamic root module where vegetation can forage for water and nutrients by altering their root profiles. The global simulations presented here show that in response to water stress events most ecosystems deepen their root profiles. In semi-arid ecosystems, this response expedites recovery (i.e., less legacy) relative to simulations without dynamics roots because access to deeper water pools after the initial event remains favorable. In wetter ecosystems, the development of deeper root profiles slows down the recovery timescale (i.e., more legacy) because the deeper root profile reduces access to nutrients and is thus unfavorable. Here, the simulations show that while the inclusion of dynamic roots might only minimally affect global patterns of Gross Primary Productivity and Evapotranspiration, the shift in root profile alters the timescale of recovery. Studies interested in the sustained response of land surfaces fluxes to climate disturbances should consider models that include dynamic root capability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Predicting weather impacts on corn production in a data-limited region using a transfer learning approach

The stability of food supply and prices may depend more on annual changes in yields from year-to-year variability in weather than on longer-term average changes from changing climatic conditions. However, the absence of high-quality data on crop yields at fine spatial resolutions in many regions of the world makes it challenging to statistically model their response to interannual variability in weather patterns. Therefore, there is a need for empirical methods that can project annual crop yield changes even in limited data regions. Here, we propose a transfer learning algorithm that uses high spatial resolution data from one region to project yields in another region with more limited data. The goal of our work is to understand what data types can be beneficial for transferring learning from a source region to a very different target region with more limited data. We utilize Long Short-Term Memory to develop a transfer learning model that is trained on historical county-level corn yield in the United States and predicts district-level corn yield variations in India. Even using smaller amounts of data in India, simulating a data-scarce region, we achieve an average root mean square error of 0.48 bu acre−1 in predicting interannual yield variations. Using Shapley values to interpret results, we explore the contribution of the different weather parameters to interannual yield variability and find a larger influence of precipitation-related variables. Our study demonstrates the usefulness of this method for transferring models of weather impacts on crop yields trained on a data-rich country to one with more limited data. It suggests the potential of applying the transfer learning model to mitigate the need for extensive raw data globally.

Vishwakarma, Srishti [ORNL] (ORCID:000000031674419↗
DOE OSTI2026

Unlocking plant-microbial interactions in deep Mollisols in the Midwestern US: Linking depth gradients in roots, microbial activity, and soil carbon in agroecosystems

Deep-rooted plants may build soil carbon (C) stocks, but most research has focused on shallow soils, leaving gaps in our understanding of how shifts in the balance between decomposition and C inputs drive soil C accumulation with depth. Thus, our objectives were to: (1) link depth gradients in root biomass with microbial activity and soil C stocks down to 1 m, and (2) examine the potential of simple C inputs to prime soil C across depths. To this end, we dug 5 quantitative soil pits in Argiudolls under mature perennial miscanthus plots in the SoyFACE Farm (Champaign-Urbana, IL). We added 13 C labeled glucose to our soils to determine the fate of simple C inputs with depth. We found that fine root biomass, total soil C, mineral-associated organic C (MAOC), particulate organic C (POC), and microbial activity (as measured by potential enzyme activity) declined with depth. POC declined more rapidly than MAOC, resulting in an increase in the ratio of MAOC-to-POC. Root biomass, enzyme activity (either acid phosphatase or n-acetyl-glucosaminadase) activity, and microbial respiration explained 74% and 38% of the variability in soil total C and MAOC, respectively, while POC was dependent on root biomass and microbial respiration (47%). Although the incorporation of simple 13 C inputs into MAOC was similar across depths, these inputs led to greater net MAOC losses in shallow soils than in deeper soils between 50 and 100 cm. The divergent impact of simple C inputs across depths may suggest that MAOC in shallow soils is more susceptible to priming losses, while C inputs into deep soils may instead be more persistent. Collectively, our results suggest that depth gradients in soil C stocks represents a balance between inputs, decomposition, and microbial necromass production and that increases in root C inputs by deep-rooted plants may have the potential to build stable MAOC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Water foraging with dynamic roots in E3SM; The role of roots in terrestrial ecosystem memory on intermediate timescales (Final Technical Report)

Terrestrial ecosystems can show sustained responses to stress events that may last for years after the initial perturbation. This phenomenon is referred to as legacy or memory and it emerges from a set of ecosystem processes that are poorly captured by Earth System and Land Surface Models. Consequently, these models struggle to predict the impact that extreme climate events have on surface energy, carbon and hydrological exchange once a stressor such as drought has relaxed or in response to repeated exposure to stress. In this project, we set out to test how the addition of dynamic root profiles in land surface models impact the capacity for Earth System Models, such as the Department of Energy’s E3SM, to capture realistic legacy effects. Observations have shown that vegetation shifts their root profiles during stress events to forage for water and these altered root profiles may take years to relax back to the initial state. We hypothesized that the alteration of the belowground root structures may be a source of terrestrial ecosystem legacy that is missing from models. To test this idea we undertook three core activities. (1) We developed a global-scale analysis of the impact of dynamic roots on ecosystem legacy building from a recently developed dynamic root module. (2) We developed new root dynamics in the DOE’s Energy Land Model (ELM) that allowed not only root profiles to shift their profile but also simultaneously alter carbon allocation to fine root pools. We then tested the impact of these new dynamics with intensive sensitivity analysis at four long-term AmeriFlux sites. (3) We undertook detailed isotopic analysis of tree rings from these AmeriFlux sites to assess ecophysiological and ecohydrological legacy to stress events that can be used to benchmark the sensitivity experiments with ELM. From these core activities, we report the following key findings. Firstly, on a global scale, the addition of dynamic roots led to chronically water-stressed ecosystems recovering faster to climate stress while wetter ecosystem showed enhanced legacy. This is because across ecosystems, stress events were almost universally associated with water shortages that led to the development of deeper root profiles. These deeper root profiles proved beneficial for recovery from drought stress. While the root dynamics did not universally improve the modeled representation of legacy it showed complex transient dynamics that emerge from the addition of dynamic roots. Secondly, the sensitivity analysis illustrated long term shifts in rooting depths away from the prescribed default profile suggesting that initializing of root profiles could benefit from spin-up simulations that converge on locally optimized root profiles. In addition, by enabling dynamic allocation some of the sites predicted unrealistically low allocation to roots (and high allocation to leaves). While these changes did not dramatically alter modeled gross primary productivity, they illustrated that without more sophisticated root processes in models, there is little penalty to dramatically disinvest in roots. Thirdly, the isotopic analysis of tree rings showed highly distinct legacy responses across species and sites. For example, T. canadensis showed reduced transpiration the year after stress events illustrated by sustained elevated $\delta^{18}$O. In contrast, A saccharum displayed elevated $\delta^{13}$C associated with reduced stomatal conductance in response to the previous years’ stress event. These geochemical signatures of legacy were present despite tree growth returning to normal the year after the stress. The results show the importance of species-level dynamics in legacy that are absent in modeling that assumes common traits within plant functional types. In summary, the work here established new avenues to explore root dynamics in models while illustrating how additional root processes are needed before dynamic carbon allocation can be implemented. Lastly, this project provided mentorship to a postdoctoral fellow, training for an early career scientist and multiple undergraduate students recruited from a minority serving institution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Divide and conquer: using RhizoVision Explorer to aggregate data from multiple root scans using image concatenation and statistical methods

Roots are important in agricultural and natural systems for determining plant productivity and soil carbon inputs. Sometimes, the amount of roots in a sample is too much to fit into a single scanned image, so the sample is divided among several scans, and there is no standard method to aggregate the data. Here, we describe and validate two methods for standardizing measurements across multiple scans: image concatenation and statistical aggregation. We developed a Python script that identifies which images belong to the same sample and returns a single, larger concatenated image. These concatenated images and the original images were processed with RhizoVision Explorer, a free and open-source software. An R script was developed, which identifies rows of data belonging to the same sample and applies correct statistical methods to return a single data row for each sample. These two methods were compared using example images from switchgrass, poplar, and various tree and ericaceous shrub species from a northern peatland and the Arctic. Most root measurements were nearly identical between the two methods except median diameter, which cannot be accurately computed by statistical aggregation. We believe the availability of these methods will be useful to the root biology community.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

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↗