National Geo-Database for Biofuel Simulations and Regional Analysis of Biorefinery Siting Based on Cellulosic Feedstock Grown on Marginal Lands

The goal of this project undertaken by GLBRC (Great Lakes Bioenergy Research Center) Area 4 (Sustainability) modelers is to develop a national capability to model feedstock supply, ethanol production, and biogeochemical impacts of cellulosic biofuels. The results of this project contribute to sustainability goals of the GLBRC; i.e. to contribute to developing a sustainable bioenergy economy: one that is profitable to farmers and refiners, acceptable to society, and environmentally sound. A sustainable bioenergy economy will also contribute, in a fundamental way, to meeting national objectives on energy security and climate mitigation. The specific objectives of this study are to: (1) develop a spatially explicit national geodatabase for conducting biofuel simulation studies and (4) locate possible sites for the establishment of cellulosic ethanol biorefineries. To address the first objective, we developed SENGBEM (Spatially Explicit National Geodatabase for Biofuel and Environmental Modeling), a 60-m resolution geodatabase of the conterminous USA containing data on: (1) climate, (2) soils, (3) topography, (4) hydrography, (5) land cover/ land use (LCLU), and (6) ancillary data (e.g., road networks, federal and state lands, national and state parks, etc.). A unique feature of SENGBEM is its 2008-2010 crop rotation data, a crucially important component for simulating productivity and biogeochemical cycles as well as land-use changes associated with biofuel cropping. ARRA support for this project and to the PNNL Joint Global Change Research Institute enabled us to create an advanced computing infrastructure to execute millions of simulations, conduct post-processing calculations, store input and output data, and visualize results. These computing resources included two components installed at the Research Data Center of the University of Maryland. The first resource was 'deltac': an 8-core Linux server, dedicated to county-level and state-level simulations and PostgreSQL database hosting. The second resource was the DOE-JGCRI 'Evergreen' cluster, capable of executing millions of simulations in relatively short periods. ARRA funding also supported a PhD student from UMD who worked on creating the geodatabases and executing some of the simulations in this study. Using a physically based classification of marginal lands, we simulated production of cellulosic feedstocks from perennial mixtures grown on these lands in the US Midwest. Marginal lands in the western states of the US Midwest appear to have significant potential to supply feedstocks to a cellulosic biofuel industry. Similar results were obtained with simulations of N-fertilized perennial mixtures. A detailed spatial analysis allowed for the identification of possible locations for the establishment of 34 cellulosic ethanol biorefineries with an annual production capacity of 5.6 billion gallons. In summary, we have reported on the development of a spatially explicit national geodatabase to conduct biofuel simulation studies and provided simulation results on the potential of perennial cropping systems to serve as feedstocks for the production of cellulosic ethanol. To accomplish this, we have employed sophisticated spatial analysis methods in combination with the process-based biogeochemical model EPIC. The results of this study will be submitted to the USDOE Bioenergy Knowledge Discovery Framework as a way to contribute to the development of a sustainable bioenergy industry. This work provided the opportunity to test the hypothesis that marginal lands can serve as sources of cellulosic feedstocks and thus contribute to avoid potential conflicts between bioenergy and food production systems. This work, we believe, opens the door for further analysis on the characteristics of cellulosic feedstocks as major contributors to the development of a sustainable bioenergy economy

Testing the EPIC Richards submodel for simulating soil water dynamics under different bottom boundary conditions

AbstractMost biogeochemical models simulate water dynamics using the tipping bucket approach, which has been often found to be too simplistic to represent vadose zone dynamics adequately under shallow groundwater conditions. Recently, a solution to the Richards equation using the Mualem–van Genuchten model (Rich‐vGM) has been added into the EPIC (Environmental Policy Integrated Climate) model to address this shortfall. Its performance was tested using lysimeters operating under free drainage (FD) and at a shallow water table (60‐ [WT60] and 120‐cm depth [WT120]). Model accuracy was also compared with the upgraded tipping bucket‐based method implemented into EPIC (the variable saturation hydraulic conductivity method [VSHC]). Soil water content (SWC) data were split into calibration and validation subsets. Model evaluation also included annual evapotranspiration (ET), percolation (PRK), and upward water movements to assess underlying soil water balance factors. The submodels provided accurate and similar results upon comparison with SWC measures under FD (Nash–Sutcliffe coefficient [NSE] = 0.26 and 0.61 using VSHC and Rich‐vGM, respectively). The Rich‐vGM model accurately reproduced observed SWC and ET (e.g., NSE = 0.70 and percentage bias [PBIAS] = −3.7% for WT120, respectively) although it slightly overestimated PRK (PBIAS = 47.8%, on average). Instead, VSHC proved unable to correctly simulate shallow groundwater conditions (e.g., NSE = −1.85 for WT60 SWC). Under shallow groundwater conditions, the Rich‐vGM method is recommended, despite the additional data required and the need to define the bottom boundary conditions according to water table fluctuations. In conclusion, the Richards solver introduced and tested in EPIC improved the model's ability to represent complex biophysical and biogeochemical processes in terrestrial ecosystems associated with the hydrological balance

National Geo-Database for Biofuel Simulations and Regional Analysis

The goal of this project undertaken by GLBRC (Great Lakes Bioenergy Research Center) Area 4 (Sustainability) modelers is to develop a national capability to model feedstock supply, ethanol production, and biogeochemical impacts of cellulosic biofuels. The results of this project contribute to sustainability goals of the GLBRC; i.e. to contribute to developing a sustainable bioenergy economy: one that is profitable to farmers and refiners, acceptable to society, and environmentally sound. A sustainable bioenergy economy will also contribute, in a fundamental way, to meeting national objectives on energy security and climate mitigation. The specific objectives of this study are to: (1) develop a spatially explicit national geodatabase for conducting biofuel simulation studies; (2) model biomass productivity and associated environmental impacts of annual cellulosic feedstocks; (3) simulate production of perennial biomass feedstocks grown on marginal lands; and (4) locate possible sites for the establishment of cellulosic ethanol biorefineries. To address the first objective, we developed SENGBEM (Spatially Explicit National Geodatabase for Biofuel and Environmental Modeling), a 60-m resolution geodatabase of the conterminous USA containing data on: (1) climate, (2) soils, (3) topography, (4) hydrography, (5) land cover/ land use (LCLU), and (6) ancillary data (e.g., road networks, federal and state lands, national and state parks, etc.). A unique feature of SENGBEM is its 2008-2010 crop rotation data, a crucially important component for simulating productivity and biogeochemical cycles as well as land-use changes associated with biofuel cropping. We used the EPIC (Environmental Policy Integrated Climate) model to simulate biomass productivity and environmental impacts of annual and perennial cellulosic feedstocks across much of the USA on both croplands and marginal lands. We used data from LTER and eddy-covariance experiments within the study region to test the performance of EPIC and, when necessary, improve its parameterization. We investigated three scenarios. In the first, we simulated a historical (current) baseline scenario composed mainly of corn-, soybean-, and wheat-based rotations as grown existing croplands east of the Rocky Mountains in 30 states. In the second scenario, we simulated a modified baseline in which we harvested corn and wheat residues to supply feedstocks to potential cellulosic ethanol biorefineries distributed within the study area. In the third scenario, we simulated the productivity of perennial cropping systems such as switchgrass or perennial mixtures grown on either marginal or Conservation Reserve Program (CRP) lands. In all cases we evaluated the environmental impacts (e.g., soil carbon changes, soil erosion, nitrate leaching, etc.) associated with the practices. In summary, we have reported on the development of a spatially explicit national geodatabase to conduct biofuel simulation studies and provided initial simulation results on the potential of annual and perennial cropping systems to serve as feedstocks for the production of cellulosic ethanol. To accomplish this, we have employed sophisticated spatial analysis methods in combination with the process-based biogeochemical model EPIC. This work provided the opportunity to test the hypothesis that marginal lands can serve as sources of cellulosic feedstocks and thus contribute to avoid potential conflicts between bioenergy and food production systems. This work, we believe, opens the door for further analysis on the characteristics of cellulosic feedstocks as major contributors to the development of a sustainable bioenergy economy

Regional scale cropland carbon budgets: Evaluating a geospatial agricultural modeling system using inventory data

Accurate quantification and clear understanding of regional scale cropland carbon (C) cycling is critical for designing effective policies and management practices that can contribute toward stabilizing atmospheric CO2 concentrations. However, extrapolating site-scale observations to regional scales represents a major challenge confronting the agricultural modeling community. This study introduces a novel geospatial agricultural modeling system (GAMS) exploring the integration of the mechanistic Environmental Policy Integrated Climate model, spatially-resolved data, surveyed management data, and supercomputing functions for cropland C budgets estimates. This modeling system creates spatiallyexplicit modeling units at a spatial resolution consistent with remotely-sensed crop identification and assigns cropping systems to each of them by geo-referencing surveyed crop management information at the county or state level. A parallel computing algorithm was also developed to facilitate the computationally intensive model runs and output post-processing and visualization. We evaluated GAMS against National Agricultural Statistics Service (NASS) reported crop yields and inventory estimated county-scale cropland C budgets averaged over 2000e2008. We observed good overall agreement, with spatial correlation of 0.89, 0.90, 0.41, and 0.87, for crop yields, Net Primary Production (NPP), Soil Organic C (SOC) change, and Net Ecosystem Exchange (NEE), respectively. However, we also detected notable differences in the magnitude of NPP and NEE, as well as in the spatial pattern of SOC change. By performing crop-specific annual comparisons, we discuss possible explanations for the discrepancies between GAMS and the inventory method, such as data requirements, representation of agroecosystem processes, completeness and accuracy of crop management data, and accuracy of crop area representation. Based on these analyses, we further discuss strategies to improve GAMS by updating input data and by designing more efficient parallel computing capability to quantitatively assess errors associated with the simulation of C budget components. The modularized design of the GAMS makes it flexible to be updated and adapted for different agricultural models so long as they require similar input data, and to be linked with socio-economic models to understand the effectiveness and implications of diverse C management practices and policies

Simulating Potential Switchgrass Production in the United States

Using results from field trials of switchgrass (Panicum virgatum L.) in the United States, the EPIC (Environmental Policy Integrated Climate) process-level agroecosystem model was calibrated, validated, and applied to simulate potential productivity of switchgrass for use as a biofuel feedstock. The model was calibrated with a regional study of 10-yr switchgrass field trials and subsequently tested against a separate compiled dataset of field trials from across the eastern half of the country. An application of the model in a national database using 8-digit watersheds as the primary modeling unit produces 30-yr average switchgrass yield estimates that can be aggregated to 18 major watersheds. The model projects average annual switchgrass productivity of greater than 7 Mg ha-1 in the Upper Mississippi, Lower Mississippi, and Ohio watersheds. The major factors limiting simulated production vary by region; low precipitation is the primary limiting factor across the western half of the country, while moderately acidic soils limit yields on lands east of the Mississippi River. Average projected switchgrass production on all crop land in the continental US is 5.6 Mg ha-1. At this level of productivity, 28.6 million hectares of crop land would be required to produce the 16 billion gallons of cellulosic ethanol called for by 2022 in the 2007 Energy Independence and Security Act. The model described here can be applied as a tool to inform the land-use and environmental consequences of switchgrass production

The Global Gridded Crop Model Intercomparison phase 1 simulation dataset

The Global Gridded Crop Model Intercomparison (GGCMI) phase 1 dataset of the Agricultural Model Intercomparison and Improvement Project (AgMIP) provides an unprecedentedly large dataset of crop model simulations covering the global ice-free land surface. The dataset consists of annual data fields at a spatial resolution of 0.5 arc-degree longitude and latitude. Fourteen crop modeling groups provided output for up to 11 historical input datasets spanning 1901 to 2012, and for up to three different management harmonization levels. Each group submitted data for up to 15 different crops and for up to 14 output variables. All simulations were conducted for purely rainfed and near-perfectly irrigated conditions on all land areas irrespective of whether the crop or irrigation system is currently used there. With the publication of the GGCMI phase 1 dataset we aim to promote further analyses and understanding of crop model performance, potential relationships between productivity and environmental impacts, and insights on how to further improve global gridded crop model frameworks. We describe dataset characteristics and individual model setup narratives. © 2019, The Author(s)

Climate and Energy-Water-Land System Interactions Technical Report to the U.S. Department of Energy in Support of the National Climate Assessment

This report provides a framework to characterize and understand the important elements of climate and energy-water-land (EWL) system interactions. It identifies many of the important issues, discusses our understanding of those issues, and presents a long-term research program research needs to address the priority scientific challenges and gaps in our understanding. Much of the discussion is organized around two discrete case studies with the broad themes of (1) extreme events and (2) regional intercomparisons. These case studies help demonstrate unique ways in which energy-water-land interactions can occur and be influenced by climate

Global wheat production with 1.5 and 2.0°C above pre‐industrial warming

Efforts to limit global warming to below 2°C in relation to the pre‐industrial level are under way, in accordance with the 2015 Paris Agreement. However, most impact research on agriculture to date has focused on impacts of warming >2°C on mean crop yields, and many previous studies did not focus sufficiently on extreme events and yield interannual variability. Here, with the latest climate scenarios from the Half a degree Additional warming, Prognosis and Projected Impacts (HAPPI) project, we evaluated the impacts of the 2015 Paris Agreement range of global warming (1.5 and 2.0°C warming above the pre‐industrial period) on global wheat production and local yield variability. A multi‐crop and multi‐climate model ensemble over a global network of sites developed by the Agricultural Model Intercomparison and Improvement Project (AgMIP) for Wheat was used to represent major rainfed and irrigated wheat cropping systems. Results show that projected global wheat production will change by −2.3% to 7.0% under the 1.5°C scenario and −2.4% to 10.5% under the 2.0°C scenario, compared to a baseline of 1980–2010, when considering changes in local temperature, rainfall, and global atmospheric CO2 concentration, but no changes in management or wheat cultivars. The projected impact on wheat production varies spatially; a larger increase is projected for temperate high rainfall regions than for moderate hot low rainfall and irrigated regions. Grain yields in warmer regions are more likely to be reduced than in cooler regions. Despite mostly positive impacts on global average grain yields, the frequency of extremely low yields (bottom 5 percentile of baseline distribution) and yield inter‐annual variability will increase under both warming scenarios for some of the hot growing locations, including locations from the second largest global wheat producer—India, which supplies more than 14% of global wheat. The projected global impact of warming <2°C on wheat production is therefore not evenly distributed and will affect regional food security across the globe as well as food prices and trade

Global gridded crop model evaluation: benchmarking, skills, deficiencies and implications

Crop models are increasingly used to simulate crop yields at the global scale, but so far there is no general framework on how to assess model performance. Here we evaluate the simulation results of 14 global gridded crop modeling groups that have contributed historic crop yield simulations for maize, wheat, rice and soybean to the Global Gridded Crop Model Intercomparison (GGCMI) of the Agricultural Model Intercomparison and Improvement Project (AgMIP). Simulation results are compared to reference data at global, national and grid cell scales and we evaluate model performance with respect to time series correlation, spatial correlation and mean bias. We find that global gridded crop models (GGCMs) show mixed skill in reproducing time series correlations or spatial patterns at the different spatial scales. Generally, maize, wheat and soybean simulations of many GGCMs are capable of reproducing larger parts of observed temporal variability (time series correlation coefficients (r) of up to 0.888 for maize, 0.673 for wheat and 0.643 for soybean at the global scale) but rice yield variability cannot be well reproduced by most models. Yield variability can be well reproduced for most major producing countries by many GGCMs and for all countries by at least some. A comparison with gridded yield data and a statistical analysis of the effects of weather variability on yield variability shows that the ensemble of GGCMs can explain more of the yield variability than an ensemble of regression models for maize and soybean, but not for wheat and rice. We identify future research needs in global gridded crop modeling and for all individual crop modeling groups. In the absence of a purely observation-based benchmark for model evaluation, we propose that the best performing crop model per crop and region establishes the benchmark for all others, and modelers are encouraged to investigate how crop model performance can be increased. We make our evaluation system accessible to all crop modelers so that other modeling groups can also test their model performance against the reference data and the GGCMI benchmark

Similar estimates of temperature impacts on global wheat yield by three independent methods

The potential impact of global temperature change on global crop yield has recently been assessed with different methods. Here we show that grid-based and point-based simulations and statistical regressions (from historic records), without deliberate adaptation or CO2 fertilization effects, produce similar estimates of temperature impact on wheat yields at global and national scales. With a 1 °C global temperature increase, global wheat yield is projected to decline between 4.1% and 6.4%. Projected relative temperature impacts from different methods were similar for major wheat-producing countries China, India, USA and France, but less so for Russia. Point-based and grid-based simulations, and to some extent the statistical regressions, were consistent in projecting that warmer regions are likely to suffer more yield loss with increasing temperature than cooler regions. By forming a multi-method ensemble, it was possible to quantify ‘method uncertainty’ in addition to model uncertainty. This significantly improves confidence in estimates of climate impacts on global food security.<br/