Search NASA⌕ Search

SEARCH · Search NASA

Results for “wheat simulation models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling root system growth around obstacles

State-of-the-Art models of Root System Architecture (RSA) do not allow simulating root growth around rigid obstacles. Yet, the presence of obstacles can be highly disruptive to the root system. We grew wheat seedlings in sealed petri dishes without obstacle and in custom 3D-printed rhizoboxes containing obstacles. Time-lapse photography was used to reconstruct the wheat root morphology network. We used the reconstructed wheat root network without obstacle to calibrate an RSA model implemented in the R-SWMS software. The root network with obstacles allowed calibrating the parameters of a new function that models the influence of rigid obstacles on wheat root growth. Experimental results show that the presence of a rigid obstacle does not affect the growth rate of the wheat root axes, but that it does influence the root trajectory after the main axis has passed the obstacle. The growth recovery time, i.e. the time for the main root axis to recover its geotropism-driven growth, is proportional to the time during which the main axis grows along the obstacle. Qualitative and quantitative comparisons between experimental and numerical results show that the proposed model successfully simulates wheat RSA growth around obstacles. Our results suggest that wheat roots follow patterns that could inspire the design of adaptive engineering flow networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Data and Code for: Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits

This repository contains the simulation outputs and processing scripts associated with the study of winter wheat traits across the United States, utilizing the Ecosys agroecosystem model. The dataset includes model results for both rainfed and irrigated winter wheat systems, supporting the findings presented in the manuscript titled "Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits." Data includes the original Ecosys simulation outputs (archived in .db format within the compressed .zip files) and extracted analysis data (stored in .pkl files for efficient processing). Python code for data processing and figure generation is provided in a Jupyter notebook. External Observational Datasets should refer to the following official repositories for the input and validation data used in this study. The eddy covariance data from the AmeriFlux network (https://ameriflux.lbl.gov/). Climate-forcing data of NLDAS-2 from NASA LDAS (https://ldas.gsfc.nasa.gov/nldas/nldas-2-forcing-data). Soil data from the Gridded Soil Survey Geographic Database (gSSURGO), available at (https://www.nrcs.usda.gov/resources/data-and-reports/gridded-soil-survey-geographic-gssurgo-database). Crop yields, planting and harvest dates from the USDA public databases (https://quickstats.nass.usda.gov/; https://webapp.rma.usda.gov/apps/actuarialinformationbrowser/CropCriteria.aspx). Satellite-derived SLOPE GPP data from ORNL DAAC (https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1786). Land use and crop progress information from the USDA Crop Data Layer and Crop Progress and Condition Gridded Layers (https://www.nass.usda.gov/Research_and_Science/). The Ecosys model code is available online at https://github.com/jinyun1tang/ECOSYS.

Wheat↗
DOE OSTI2023

Adaptation Strategies Strongly Reduce the Future Impacts of Climate Change on Simulated Crop Yields

Abstract Simulations of crop yield due to climate change vary widely between models, locations, species, management strategies, and Representative Concentration Pathways (RCPs). To understand how climate and adaptation affects yield change, we developed a meta‐model based on 8703 site‐level process‐model simulations of yield with different future adaptation strategies and climate scenarios for maize, rice, wheat and soybean. We tested 10 statistical models, including some machine learning models, to predict the percentage change in projected future yield relative to the baseline period (2000–2010) as a function of explanatory variables related to adaptation strategy and climate change. We used the best model to produce global maps of yield change for the RCP4.5 scenario and identify the most influential variables affecting yield change using Shapley additive explanations. For most locations, adaptation was the most influential factor determining the projected yield change for maize, rice and wheat. Without adaptation under RCP4.5, all crops are expected to experience average global yield losses of 6%–21%. Adaptation alleviates this average projected loss by 1–13 percentage points. Maize was most responsive to adaptive practices with a projected mean yield loss of −21% [range across locations: −63%, +3.7%] without adaptation and −7.5% [range: −46%, +13%] with adaptation. For maize and rice, irrigation method and cultivar choice were the adaptation types predicted to most prevent large yield losses, respectively. When adaptation practices are applied, some areas are predicted to experience yield gains, especially at northern high latitudes. These results reveal the critical importance of implementing adequate adaptation strategies to mitigate the impact of climate change on crop yields.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Comparison of structurally diverse simulation models for prediction of epidemic outcomes caused by a long-distance dispersed pathogen

Long-distance dispersal (LDD) pathogens pose substantial challenges for epidemic control due to their ability to generate new infection foci at great distances. While various modeling approaches have been developed to understand and manage such outbreaks, little work has compared how models of different structures behave under shared conditions. Here, in this study, we compare four structurally distinct epidemiological models — EPIMUL, GEMF, PoPS, and Warwick — each adapted to simulate the spread of wheat stripe rust (WSR), a wind-dispersed LDD pathogen, under identical epidemiological parameters and dispersal kernel. Using data from a controlled field experiment, we evaluate the ability of each model to replicate disease prevalence under nine intervention scenarios that vary in timing and culling area. While the models differ substantially in design — ranging from spatial grid-based to network-based and raster-based frameworks — the shared dispersal kernel allowed for close alignment in their predictions. All models accurately captured general epidemic trends, particularly the strong effect of early intervention on disease suppression. We qualitatively compared their behavioral responses across scenarios and also evaluated an ensemble prediction by averaging across model outputs. Our findings highlight how integrating shared epidemiological components into distinct modeling frameworks can improve consistency and accuracy, while reinforcing the importance of early culling in managing LDD pathogen outbreaks.

Dispersal kernel↗
DOE OSTI2024

Representing cropping systems with the MEMS 2 ecosystem model

Abstract Croplands have been the focus of substantial investigation due to their considerable potential for sequestering carbon. Understanding the potential for soil organic carbon (SOC) sequestration and necessary management strategies will be enabled with accurate process‐based models. Accurately representing crop growth and agricultural practices will be critical for realistic SOC modeling. The MEMS 2 model incorporates a current understanding of SOC formation and stabilization, measurable SOC pools, and deep SOC dynamics and is seen as a highly promising tool to inform management intervention for SOC sequestration. Thus far, MEMS 2 has been developed to represent grasslands. In this study, we further developed MEMS 2 to model annual grain crops and common agricultural practices, such as irrigation, fertilization, harvesting, and tillage. Using four Ameriflux sites, we demonstrated an accurate simulation of crop growth and development. Model performance was strong for simulating aboveground biomass (index of agreement [d] range of 0.89–0.98) and green leaf area index (dfrom 0.90 to 0.96) across corn, soybean, and winter wheat. Good agreement with observations was also achieved for net ecosystem CO 2 exchange (dfrom 0.90 to 0.96), evapotranspiration (dfrom 0.91 to 0.94), and soil temperature (dof 0.96), while discrepancy with the available soil water content data remain (dfrom 0.14 to 0.81 at four depths to 100 cm). While we will continue model testing and improvement, MEMS 2 (version 2.14) has now demonstrated its ability to effectively simulate the growth of common grain crops and practices.

Agriculture↗
DOE OSTI2020

Fine-Scale Variability of Observed and Simulated Surface Albedo Over the Southern Great Plains

Surface albedo can be highly variable in both space and time. The Department of Energy’s Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field study provides a unique opportunity to characterize the variability over the Southern Great Plains of the United States using integrated tower, aircraft, and satellite observations of the surface albedo and Normalized Difference Vegetation Index (NDVI). The primary advantage of the aircraft and satellite observations is the ability to examine the spatial distribution of surface albedo and NDVI over a large area, while the tower measurements have the ability to examine both diurnal and day-to-day variability at a single location. In general, consistency was found between both the broadband (BB) albedo and NDVI measured from the surface, air and space. There was a small increase from 0.186 to 0.194 in the aircraft BB surface albedo between May and September (about 4% change), while the MODIS black-sky BB surface albedo increased from 0.151 to 0.166 over the same period (about 10% change) while the standard deviations in the aircraft and MODIS BB albedo were similar. The largest seasonal differences in both the aircraft BB albedo and NDVI were found for areas with winter wheat or forest, while areas with pasture or grasses showed a smaller seasonal diversity. The Weather Research and Forecasting (WRF) model was used to simulate the BB surface albedo. In comparison with the aircraft and satellite observations, the WRF-simulated BB surface albedo had no seasonal change and a much narrower distribution.

Berg, Larry K.↗
DOE OSTI2024

Tethys: A Spatiotemporal Downscaling Model for Global Water Demand

Humans use water for many important tasks, such as drinking, growing food, and cooling power plants. Since future water demands depend on complex global interactions between economic sectors (e.g., demand for wheat in one country causing demand for water to grow that wheat in another country), it is often modeled at coarse spatial and temporal scales as part of models that account for complex, multi-sector system dynamics. However, models that project future water availability typically simulate physical processes at much finer scales. Tethys enables integration between these kinds of models by downscaling region-scale water demand projections using sector-specific proxies and formulas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

A regional nuclear conflict would compromise global food security

A limited nuclear war between India and Pakistan could ignite fires large enough to emit more than 5 Tg of soot into the stratosphere. Climate model simulations have shown severe resulting climate perturbations with declines in global mean temperature by 1.8 °C and precipitation by 8%, for at least 5 y. Here we evaluate impacts for the global food system. Six harmonized state-of-the-art crop models show that global caloric production from maize, wheat, rice, and soybean falls by 13 (±1)%, 11 (±8)%, 3 (±5)%, and 17 (±2)% over 5 y. Total single-year losses of 12 (±4)% quadruple the largest observed historical anomaly and exceed impacts caused by historic droughts and volcanic eruptions. Colder temperatures drive losses more than changes in precipitation and solar radiation, leading to strongest impacts in temperate regions poleward of 30°N, including the United States, Europe, and China for 10 to 15 y. Integrated food trade network analyses show that domestic reserves and global trade can largely buffer the production anomaly in the first year. Persistent multiyear losses, however, would constrain domestic food availability and propagate to the Global South, especially to food-insecure countries. By year 5, maize and wheat availability would decrease by 13% globally and by more than 20% in 71 countries with a cumulative population of 1.3 billion people. In view of increasing instability in South Asia, this study shows that a regional conflict using <1% of the worldwide nuclear arsenal could have adverse consequences for global food security unmatched in modern history.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling the impact of measured and projected climate and management systems on agricultural fields: Surface runoff, soil moisture, and soil erosion

Abstract As global climate change poses a challenge to crop production, it is imperative to prioritize effective adaptation of agricultural systems based on a scientific understanding of likely impacts. In this study, we applied an integrated watershed modeling framework to examine the impacts of projected climate on runoff, soil moisture, and soil erosion under different management systems in Central Oklahoma. The proposed model uses measured climate data and three downscaled ensembles from the Coupled Model Intercomparison Project Phase 6 (CMIP6) at the water resources and erosion watershed to understand the impact of climate change and various climate conditions under three management systems: (1) continuous winter wheat (Triticum aestivum) under conventional tillage (WW‐CT; baseline system), (2) continuous winter wheat under no‐till (WW‐NT), and (3) cool and warm season forage cover crop mixes under no‐till (CC‐NT). The study indicates that the occurrence of agricultural drought is projected to increase while erosion rates will remain unchanged under the WW‐CT. In contrast, climate simulations imposed on the WW‐NT and CC‐NT systems significantly reduce runoff and sediment while preserving soil moisture levels. Especially, implementing the CC‐NT system can bolster food security and foster sustainable farming practices in Central Oklahoma in the face of a changing climate.

Environmental Sciences & Ecology↗
DOE OSTI2023

Adaptation to climate change through strategic integration of long fallow into cropping system in a dryland Mediterranean-type environment

The crop-growing region of Western Australia characterized by a Mediterranean-type climate is projected to become warmer and drier. Appropriate selection of crop sequences will be of importance to cope with these climatic changes for this largest grain-producing region of Australia. Through linking a widely used crop model (APSIM), 26 General Circulation Models (GCMs) with one Shared Socioeconomic Pathway (SSP585) and economic analysis, we explored how the climate change would affect dryland wheat cropping and whether/how long fallow (the practice of leaving a field out of production for an entire growing season) could be integrated into wheat cropping system in Western Australia. The potential adaptation of long fallow into wheat system was assessed with four fixed rotations (fallow-wheat, fallow-wheat-wheat, fallow-wheat-wheat-wheat, and fallow-wheat-wheat-wheat-wheat) and four flexible sowing rule-based rotations (the land was fallowed if sowing rule was not met), compared with continuous wheat. The simulation results at four representing locations show that climate change would have negative impacts on both yield and economic return of continuous wheat cropping in Western Australia. Wheat after fallow out-yielded and out-profited wheat after wheat under future climate. But integrating fallow into wheat cropping systems with the above fixed rotations would lead to yield and economic loss. By contrast, cropping systems in which fallowing took place when sowing condition could not be met at a certain time would achieve comparable yield and economic return to continuous wheat, with wheat yield being only 5% less than continuous wheat and the gross margin being $12ha -1 more than continuous wheat averaged across locations. We highlight strategic integration of long fallow into cropping system in a dryland Mediterranean-type environment would have a great potential to cope with future climate change. These findings can be extended into other Mediterranean-type cropping regions in Australia and beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

DEEPER: An Intergrated Platform for Deeper Roots

Crops with deeper roots would have multiple benefits, including better drought tolerance, reduced requirement for nitrogen fertilizer, and better sequestration of atmospheric CO 2 . DEEPER is an integrated platform of phenomic, genomic, and in silico technologies to generate maize lines with deeper roots. DEEPER is: LEADER (Leaf Elemental Accumulation from Deep Roots) is a breakthrough technology to nondestructively measure rooting depth by using the plant itself as a sensor. LEADER uses handheld X-ray Fluorescence spectrometry to quantify foliar accumulation of elements that are differentially distributed in the soil profile. LEADER is nondestructive and is orders of magnitude cheaper, faster, and more precise than any competing assay of rooting depth in the field. LEADER is able to distinguish deep-rooted from shallow-rooted maize lines in the field without the need for costly and noisy soil coring. RootRobot/DIRT3D, to automatically phenotype root architecture in any field, combining RootRobot, a mechatronics platform to excavate, clean, section, and image mature root crowns, with DIRT3D, software to quantify architectural traits in 3D. Anatomics, a high-throughput platform to phenotype root anatomy, combining LAT 2.0, a technology for 3D imaging of root anatomy and composition, with RootScan3D, software to automatically extract 3D anatomical and cell wall composition metrics from LAT 2.0 output. Using this platform we discovered two novel root traits, parenchyma cell wall thickness and multiseriate cortical sclerenchyma, that improve rooting depth and drought tolerance in maize and wheat. OpenSimRoot/Deep, software to simulate root interaction with hard subsoils. Using this platform we discovered novel concepts regarding how to increase crop rooting depth by modulating how individual root axes respond to hard soil. DeepGenes, a toolkit of genes, parent lines, and genomic selection strategies to enable breeding hybrids with deeper roots. We discovered 3 novel root genes that increase rooting depth in maize and wheat. DEEPER discovered novel root phenotypes for deeper rooting, and delivered validated ideotypes for deeper-rooted maize; novel technologies to rapidly assess root depth, root architecture and anatomy in field-grown plants; novel software tools for root modeling and 3D image analysis of root architecture and anatomy; and validated genes and genomic selection models to deploy traits for deeper rooting in maize breeding. Each DEEPER technology is transformative in its own right, and exceeds existing technologies. They are mutually synergistic, deployable for field-grown plants, and are ready for application. The phenotyping and modeling technologies are readily applicable to many crops, and genetic leads in maize may have utility in other grasses. Taken as a whole they represent a transformative platform to develop deeper-rooted crops, with greater drought tolerance, reduced fertilizer requirement, and greater carbon sequestration.

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↗
DOE OSTI2021

Too Good to Be True? Evaluation of Colonoscopy Sensitivity Assumptions Used in Policy Models

Background: Models can help guide colorectal cancer screening policy. Although models are carefully calibrated and validated, there is less scrutiny of assumptions about test performance. Methods: We examined the validity of the CRC-SPIN model and colonoscopy sensitivity assumptions. Standard sensitivity assumptions, consistent with published decision analyses, assume sensitivity equal to 0.75 for diminutive adenomas (<6 mm), 0.85 for small adenomas (6–10 mm), 0.95 for large adenomas (≥10 mm), and 0.95 for preclinical cancer. We also selected adenoma sensitivity that resulted in more accurate predictions. Targets were drawn from the Wheat Bran Fiber study. In the study, we examined how well the model predicted outcomes measured over a three-year follow-up period, including the number of adenomas detected, the size of the largest adenoma detected, and incident colorectal cancer. Results: Using standard sensitivity assumptions, the model predicted adenoma prevalence that was too low (42.5% versus 48.9% observed, with 95% confidence interval 45.3%–50.7%) and detection of too few large adenomas (5.1% versus 14.% observed, with 95% confidence interval 11.8%–17.4%). Predictions were close to targets when we set sensitivities to 0.20 for diminutive adenomas, 0.60 for small adenomas, 0.80 for 10- to 20-mm adenomas, and 0.98 for adenomas 20 mm and larger. Conclusions: Colonoscopy may be less accurate than currently assumed, especially for diminutive adenomas. Alternatively, the CRC-SPIN model may not accurately simulate onset and progression of adenomas in higher-risk populations. Impact: Misspecification of either colonoscopy sensitivity or disease progression in high-risk populations may affect the predicted effectiveness of colorectal cancer screening. When possible, decision analyses used to inform policy should address these uncertainties.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Isolation of Methane Enriched Bacterial Communities and Application as Wheat Biofertilizer under Drought Conditions: An Environmental Contribution

The search for methanotrophs as plant-growth-promoting rhizobacteria (PGPR) presents an important contribution to mitigating the impact of global warming by restoring the natural soil potential for consuming methane while benefiting plants during droughts. Our in silico simulations suggest that water, produced as a byproduct of methane oxidation, can satisfy the cell growth requirement. In addition to water, methanotrophs can produce metabolites that stimulate plant growth. Considering this, we proposed that applying methanotrophs as PGPR can alleviate the effect of droughts on crops, while stimulating atmospheric methane consumption. In this work, we isolated a series of methanotrophic communities from the rhizospheres of different crops, including Italian sweet pepper and zucchini, using an atmosphere enriched with pure methane gas, to determine their potential for alleviating drought stress in wheat plants. Subsequently, 23 strains of nonmethanotrophic bacteria present in the methanotrophic communities were isolated and characterized. We then analyzed the contribution of the methane-consuming consortia to the improvement of plant growth under drought conditions, showing that some communities contributed to increases in the wheat plants’ lengths and weights, with statistically significant differences according to ANOVA models. Furthermore, we found that the presence of methane gas can further stimulate the plant–microbe interactions, resulting in larger plants and higher drought tolerance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Higher than expected N 2 O emissions from soybean crops in the Pampas Region of Argentina: Estimates from DayCent simulations and field measurements

In developing countries, agriculture generally represents a large fraction of GHG emissions reported in National Inventories, and emissions are typically estimated using Tier 1 IPCC guidelines. However, field data and locally adapted simulation models can improve the accuracy of IPCC estimations. In this report we aimed to quantify anthropogenic N2O emissions from croplands of Argentina through field measurements, model simulations and IPCC guidelines. Here we measured N 2 O emissions and their controlling factors in 62 plots of the Pampas Region with corn, soybean and wheat/soybean crops and in unmanaged grasslands. We accounted for gross emissions from crops and background emissions from unmanaged grasslands to calculate net anthropogenic emissions from crops as the difference between them. We calibrated and evaluated the DayCent model and then simulated different weather and management scenarios. Finally, we applied IPCC guidelines to estimate anthropogenic N 2 O emissions at the same plots. The DayCent model accurately simulated annual N 2 O emission for all crops as compared to measured data (RMSE = 1.4 g N ha -1 day -1 ). Measured and simulated emissions in soybean crops were higher than in corn and wheat/soybean crops. Gross N 2 O emissions ranged from 1.4 to 5.1 kg N ha -1 yr -1 for current environmental (soil and weather) and management (crops and fertilizer doses) conditions. Background emissions ranged between 1.1 and 1.3 kg N ha -1 yr -1 , and therefore net anthropogenic emissions ranged from 0.3 to 4.0 kg N ha-1 yr -1 . IPCC Tier 1 emission factors underestimated N2O releases from soybean, that were on average 4.87 times greater when estimated with DayCent and observations (0.53 vs 2.47 and 2.69 kg N ha -1 yr -1 , respectively). On the contrary, IPCC estimates for corn and wheat/soybean crops were similar to modeled and measured values. Our results suggest that N 2 O emissions from the vast 15 million ha of soybean croplands in the Pampas Region may be substantially underestimated.

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

Agriculture’s Potential Regional Economic Contributions to the United States Economy When Supplying Feedstock to the Bio-Economy

The economic impact of obtaining biomass could become significant to U.S. rural economies via the establishment of a bioeconomy. In 2023, the Bioenergy Technologies Office (BETO) and Oak Ridge National Laboratory provided a road map to obtain over a billion tons of biomass for conversion to bioenergy and other products. Using information from this roadmap, this study estimates the potential positive and negative economic impacts that occur because of land use change, along with increased technological advances. This is achieved by using the input–output model, IMPLAN, and impacting 179 Bureau of Economic Analysis regions in the conterminous United States. Biomass included in the analysis comprises dedicated energy crops, crop residues, and forest residues. The analysis found that managing pastures more intensively could result in releasing land to produce dedicated energy crops on 30.8 million hectares, resulting in the production of 361 million metric tons of biomass. This, coupled with crop residues from barley, corn, oats, sorghum, and wheat (162 million metric tons), plus forest residues (41 million metric tons), provide 564 million dry metric tons of biomass. Assuming the price for biomass in 2023 dollars was USD 77 per dry metric-ton, this additional production results in an economic benefit for the nation of USD 619 billion, an increase from the Business As Is scenario (Baseline) of almost USD 100 billion per year, assuming a mature biomass industry. An additional 700,000 jobs are required to grow, harvest/collect, and transport the biomass material from the land.

ForSEAM↗