• Energy Research

Abstract This paper introduces a spatially-explicit bioeconomic model for the study of potential cellulosic biomass supply. For biomass crops to begin to replace current crops, farmers must earn more from them than from current crops. Using weather, topographic and soil data, the terrestrial ecosystem model, EPIC, dynamically simulates multiple cropping systems that vary by crop rotation, tillage, fertilization and residue removal rate. EPIC generates predicted crop yield and environmental outcomes over multiple watersheds. These EPIC results are used to parameterize a regional profit-maximization mathematical programming model that identifies profitable cropping system choices. The bioeconomic model is calibrated to 2007–09 crop production in a 9-county region of southwest Michigan. A simulation of biomass supply in response to rising biomass prices shows that cellulosic residues from corn stover and wheat straw begin to be supplied at minimum delivered biomass:corn grain price ratios of 0.15 and 0.18, respectively. At the mean corn price of $162.6/Mg ($4.13 per bushel) at commercial moisture content during 2007–2009, these ratios correspond to stover and straw prices of $24 and $29 per dry Mg. Perennial bioenergy crops begin to be supplied at price levels 2–3 times higher. Average biomass transport costs to the biorefinery plant range from $6 to $20/Mg compared to conventional crop production practices in the area, biomass supply from annual crop residues increased greenhouse gas emissions and reduced water quality through increased nutrient loss. By contrast, perennial cellulosic biomass crop production reduced greenhouse gas emissions and improved water quality.

AbstractThe current or “conventional” paradigm for producing process energy in a biorefinery processing cellulosic biomass is on‐site energy recovery through combustion of residual solids and biogas generated by the process. Excess electricity is then exported, resulting in large greenhouse gas (GHG) credits. However, this approach will cause lifecycle GHG emissions of biofuels to increase as more renewable energy sources (wind, solar, etc.) participate in grid‐electricity generation, and the GHG credits from displacing fossil fuel decrease. To overcome this drawback, a decentralized (depot‐based) biorefinery can be integrated with a coal‐fired power plant near a large urban area. In an integrated, decentralized, depot‐based biorefinery (IDB), the residual solids are co‐fired with coal either in the adjacent power plant or in coal‐fired boilers elsewhere to displace coal. An IDB system does not rely on indirect GHG credits through grid‐electricity displacement. In an IDB system, biogas from the wastewater treatment facility is also upgraded to biomethane and used as a transportation biofuel. The GHG savings per unit of cropland in the IDB systems (2.7–2.9 MgCO2/ha) are 1.5–1.6 fold greater than those in a conventional centralized system (1.7–1.8 MgCO2/ha). Importantly, the biofuel selling price in the IDBs is lower by 28–30 cents per gasoline‐equivalent liter than in the conventional centralized system. Furthermore, the total capital investment per annual biofuel volume in the IDB is much lower (by ~80%) than that in the conventional centralized system. Therefore, utilization of biomethane and residual solids in the IDB systems leads to much lower biofuel selling prices and significantly greater GHG savings per unit of cropland participating in the biorefinery system compared to the conventional centralized biorefineries.

Climate change is one of the most compelling modern issues and has important implications for almost every aspect of natural and human systems. The Soil and Water Assessment Tool (SWAT) model has been applied worldwide to support sustainable land and water management in a changing climate. However, the inadequacies of the existing carbon algorithm in SWAT limit its application in assessing impacts of human activities on CO2 emission, one important source of greenhouse gasses (GHGs) that traps heat in the earth system and results in global warming. In this research, we incorporate a revised version of the CENTURY carbon model into SWAT to describe dynamics of soil organic matter (SOM)-residue and simulate land-atmosphere carbon exchange. We test this new SWAT-C model with daily eddy covariance (EC) observations of net ecosystem exchange (NEE) and evapotranspiration (ET) and annual crop yield at six sites across the U.S. Midwest. Results show that SWAT-C simulates well multi-year average NEE and ET across the spatially distributed sites and capture the majority of temporal variation of these two variables at a daily time scale at each site. Our analyses also reveal that performance of SWAT-C is influenced by multiple factors, such as crop management practices (irrigated vs. rainfed), completeness and accuracy of input data, crop species, and initialization of state variables. Overall, the new SWAT-C demonstrates favorable performance for simulating land-atmosphere carbon exchange across agricultural sites with different soils, climate, and management practices. SWAT-C is expected to serve as a useful tool for including carbon flux into consideration in sustainable watershed management under a changing climate. We also note that extensive assessment of SWAT-C with field observations is required for further improving the model and understanding potential uncertainties of applying it across large regions with complex landscapes.

Abstract. Statistical emulation allows combining advantageous features of statistical and process-based crop models for understanding the effects of future climate changes on crop yields. We describe here the development of emulators for nine process-based crop models and five crops using output from the Global Gridded Model Intercomparison Project (GGCMI) Phase 2. The GGCMI Phase 2 experiment is designed with the explicit goal of producing a structured training dataset for emulator development that samples across four dimensions relevant to crop yields: atmospheric carbon dioxide (CO2) concentrations, temperature, water supply, and nitrogen inputs (CTWN). Simulations are run under two different adaptation assumptions: that growing seasons shorten in warmer climates, and that cultivar choice allows growing seasons to remain fixed. The dataset allows emulating the climatological-mean yield response of all models with a simple polynomial in mean growing-season values. Climatological-mean yields are a central metric in climate change impact analysis; we show here that they can be captured without relying on interannual variations. In general, emulation errors are negligible relative to differences across crop models or even across climate model scenarios; errors become significant only in some marginal lands where crops are not currently grown. We demonstrate that the resulting GGCMI emulators can reproduce yields under realistic future climate simulations, even though the GGCMI Phase 2 dataset is constructed with uniform CTWN offsets, suggesting that the effects of changes in temperature and precipitation distributions are small relative to those of changing means. The resulting emulators therefore capture relevant crop model responses in a lightweight, computationally tractable form, providing a tool that can facilitate model comparison, diagnosis of interacting factors affecting yields, and integrated assessment of climate impacts.

Ensembles of process-based crop models are increasingly used to simulate crop growth for scenarios of temperature and/or precipitation changes corresponding to different projections of atmospheric CO2 concentrations. This approach generates large datasets with thousands of simulated crop yield data. Such datasets potentially provide new information but it is difficult to summarize them in a useful way due to their structural complexities. An associated issue is that it is not straightforward to compare crops and to interpolate the results to alternative climate scenarios not initially included in the simulation protocols. Here we demonstrate that statistical models based on random-coefficient regressions are able to emulate ensembles of process-based crop models. An important advantage of the proposed statistical models is that they can interpolate between temperature levels and between CO2 concentration levels, and can thus be used to calculate temperature and [CO2] thresholds leading to yield loss or yield gain, without re-running the original complex crop models. Our approach is illustrated with three yield datasets simulated by 19 maize models, 26 wheat models, and 13 rice models. Several statistical models are fitted to these datasets, and are then used to analyze the variability of the yield response to [CO2] and temperature. Based on our results, we show that, for wheat, a [CO2] increase is likely to outweigh the negative effect of a temperature increase of +2°C in the considered sites. Compared to wheat, required levels of [CO2] increase are much higher for maize, and intermediate for rice. For all crops, uncertainties in simulating climate change impacts increase more with temperature than with elevated [CO2].

Projected global warming and population growth will reduce future water availability for agriculture. Thus, it is essential to increase the efficiency in using water to ensure crop productivity. Quantifying crop water use (WU; i.e. actual evapotranspiration) is a critical step towards this goal. Here, sixteen wheat simulation models were used to quantify sources of model uncertainty and to estimate the relative changes and variability between models for simulated WU, water use efficiency (WUE, WU per unit of grain dry mass produced), transpiration efficiency (Teff, transpiration per kg of unit of grain yield dry mass produced), grain yield, crop transpiration and soil evaporation at increased temperatures and elevated atmospheric carbon dioxide concentrations ([CO2]). The greatest uncertainty in simulating water use, potential evapotranspiration, crop transpiration and soil evaporation was due to differences in how crop transpiration was modelled and accounted for 50% of the total variability among models. The simulation results for the sensitivity to temperature indicated that crop WU will decline with increasing temperature due to reduced growing seasons. The uncertainties in simulated crop WU, and in particularly due to uncertainties in simulating crop transpiration, were greater under conditions of increased temperatures and with high temperatures in combination with elevated atmospheric [CO2] concentrations. Hence the simulation of crop WU, and in particularly crop transpiration under higher temperature, needs to be improved and evaluated with field measurements before models can be used to simulate climate change impacts on future crop water demand.