Abstract: 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 rerunning 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 degrees 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].

Abstract: Efforts to limit global warming to below 2 degrees 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 degrees 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 degrees 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 degrees C scenario and -2.4% to 10.5% under the 2.0 degrees 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 degrees 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.

Abstract: Dynamic global vegetation models (DGVMs) are important tools for modelling impacts of global change on ecosystem services. However, most models do not take full account of human land management and land use and land cover changes (LULCCs). We integrated croplands and pasture and their management and natural vegetation recovery and succession following cropland abandonment into the LPJ-GUESS DGVM. The revised model was applied to Africa as a case study to investigate the implications of accounting for land use on net ecosystem carbon balance (NECB) and the skill of the model in describing agricultural production and reproducing trends and patterns in vegetation structure and function. The seasonality of modelled monthly fraction of absorbed photosynthetically active radiation (FPAR) was shown to agree well with satellite-inferred normalised difference vegetation index (NDVI). In regions with a large proportion of cropland, the managed land addition improved the FPAR vs. NDVI fit significantly. Modelled 1991-1995 average yields for the seven most important African crops, representing potential optimal yields limited only by climate forcings, were generally higher than reported FAO yields by a factor of 2-6, similar to previous yield gap estimates. Modelled inter-annual yield variations during 1971-2005 generally agreed well with FAO statistics, especially in regions with pronounced climate seasonality. Modelled land-atmosphere carbon fluxes for Africa associated with land use change (0.07 PgC yr(-1) release to the atmosphere for the 1980s) agreed well with previous estimates. Cropland management options (residue removal, grass as cover crop) were shown to be important to the land-atmosphere carbon flux for the 20th century.