• Energy Research

Asseng, S. et al. Crop models are essential tools for assessing the threat of climate change to local and global food production1. Present models used to predict wheat grain yield are highly uncertain when simulating how crops respond to temperature2. Here we systematically tested 30 different wheat crop models of the Agricultural Model Intercomparison and Improvement Project against field experiments in which growing season mean temperatures ranged from 15 °C to 32 °C, including experiments with artificial heating. Many models simulated yields well, but were less accurate at higher temperatures. The model ensemble median was consistently more accurate in simulating the crop temperature response than any single model, regardless of the input information used. Extrapolating the model ensemble temperature response indicates that warming is already slowing yield gains at a majority of wheat-growing locations. Global wheat production is estimated to fall by 6% for each °C of further temperature increase and become more variable over space and time. We thank the Agricultural Model Intercomparison and Improvement Project and its leaders C. Rosenzweig from NASA Goddard Institute for Space Studies and Columbia University (USA), J. Jones from University of Florida (USA), J. Hatfield from United States Department of Agriculture (USA) and J. Antle from Oregon State University (USA) for support. We also thank M. Lopez from CIMMYT (Turkey), M. Usman Bashir from University of Agriculture, Faisalabad (Pakistan), S. Soufizadeh from Shahid Beheshti University (Iran), and J. Lorgeou and J-C. Deswarte from ARVALIS—Institut du Végétal (France) for assistance with selecting key locations and quantifying regional crop cultivars, anthesis and maturity dates and R. Raymundo for assistance with GIS. S.A. and D.C. received financial support from the International Food Policy Research Institute (IFPRI). C.S. was funded through USDA National Institute for Food and Agriculture award 32011-68002-30191. C.M. received financial support from the KULUNDA project (01LL0905L) and the FACCE MACSUR project (031A103B) funded through the German Federal Ministry of Education and Research (BMBF). F.E. received support from the FACCE MACSUR project (031A103B) funded through the German Federal Ministry of Education and Research (2812ERA115) and E.E.R. was funded through the German Science Foundation (project EW 119/5-1). M.J. and J.E.O. were funded through the FACCE MACSUR project by the Danish Strategic Research Council. K.C.K. and C.N. were funded by the FACCE MACSUR project through the German Federal Ministry of Food and Agriculture (BMEL). F.T., T.P. and R.P.R. received financial support from FACCE MACSUR project funded through the Finnish Ministry of Agriculture and Forestry (MMM); F.T. was also funded through National Natural Science Foundation of China (No. 41071030). C.B. was funded through the Helmholtz project ‘REKLIM—Regional Climate Change: Causes and Effects’ Topic 9: ‘Climate Change and Air Quality’. M.P.R. and P.D.A. received funding from the CGIAR Research Program on Climate Change, Agriculture, and Food Security (CCAFS). G.O’L. was funded through the Australian Grains Research and Development Corporation and the Department of Environment and Primary Industries Victoria, Australia. R.C.I. was funded by Texas AgriLife Research, Texas A&M University. E.W. and Z.Z. were funded by CSIRO and the Chinese Academy of Sciences (CAS) through the research project ‘Advancing crop yield while reducing the use of water and nitrogen’ and by the CSIRO-MoE PhD Research Program. Peer reviewed

Abstract Conventional methods of modelling impacts of future climate change on crop yields often rely on a limited selection of projections for representing uncertainties in future climate. However, large ensembles of climate projections offer an opportunity to estimate yield responses probabilistically. This study demonstrates an approach to probabilistic yield estimation using impact response surfaces (IRSs). These are constructed from a set of sensitivity simulations that explore yield responses to a wide range of changes in temperature and precipitation. Options for adaptation and different levels of future atmospheric carbon dioxide concentration [CO2] defined by representative concentration pathways (RCP4.5 and RCP8.5) were also considered. Model-based IRSs were combined with probabilistic climate projections to estimate impact likelihoods for yields of spring barley (Hordeum vulgare L.) in Finland during the 21st century. Probabilistic projections of climate for the same RCPs were overlaid on IRSs for corresponding [CO2] levels throughout the century and likelihoods of yield shortfall calculated with respect to a threshold mean yield for the baseline (1981–2010). Results suggest that cultivars combining short pre- and long post-anthesis phases together with earlier sowing dates produce the highest yields and smallest likelihoods of yield shortfall under future scenarios. Higher [CO2] levels generally compensate for yield losses due to warming under the RCPs. Yet, this does not happen fully under the more moderate warming of RCP4.5 with a weaker rise in [CO2], where there is a chance of yield shortfall throughout the century. Under the stronger warming but more rapid [CO2] increase of RCP8.5, the likelihood of yield shortfall drops to zero from mid-century onwards. Whilst the incremental IRS-based approach simplifies the temporal and cross-variable complexities of projected climate, it was found to offer a close approximation of evolving future likelihoods of yield impacts in comparison to a more conventional scenario-based approach. The IRS approach is scenario-neutral and existing plots can be used in combination with any new scenario that falls within the sensitivity range without the need to perform new runs with the impact model. A single crop model is used for demonstration, but an ensemble IRS approach could additionally capture impact model uncertainties.

The data set includes a current representative management treatment from detailed, quality-tested sentinel field experiments with wheat from four contrasting environments including Australia, The Netherlands, India and Argentina. Measurements include local daily climate data (solar radiation, maximum and minimum temperature, precipitation, surface wind, dew point temperature, relative humidity, and vapor pressure), soil characteristics, frequent growth, nitrogen in crop and soil, crop and soil water and yield components. Simulations include results from 27 wheat models and a sensitivity analysis with 26 models and 30 years (1981-2010) for each location, for elevated atmospheric CO2 and temperature changes, a heat stress sensitivity analysis at anthesis, and a sensitivity analysis with soil and crop management variations and a Global Climate Model end-century scenario. Data access via DOI 10.17026/DANS-ZB6-6FVQ.

Adaptation of crops to climate change has to be addressed locally due to the variability of soil, climate and the specific socio-economic settings influencing farm management decisions. Adaptation of rainfed cropping systems in the Mediterranean is especially challenging due to the projected decline in precipitation in the coming decades, which will increase the risk of droughts. Methods that can help explore uncertainties in climate projections and crop modelling, such as impact response surfaces (IRSs) and ensemble modelling, can then be valuable for identifying effective adaptations. Here, an ensemble of 17 crop models was used to simulate a total of 54 adaptation options for rainfed winter wheat (Triticum aestivum) at Lleida (NE Spain). To support the ensemble building, an ex post quality check of model simulations based on several criteria was performed. Those criteria were based on the “According to Our Current Knowledge” (AOCK) concept, which has been formalized here. Adaptations were based on changes in cultivars and management regarding phenology, vernalization, sowing date and irrigation. The effects of adaptation options under changed precipitation (P), temperature (T), [CO2] and soil type were analysed by constructing response surfaces, which we termed, in accordance with their specific purpose, adaptation response surfaces (ARSs). These were created to assess the effect of adaptations through a range of plausible P, T and [CO2] perturbations. The results indicated that impacts of altered climate were predominantly negative. No single adaptation was capable of overcoming the detrimental effect of the complex interactions imposed by the P, T and [CO2] perturbations except for supplementary irrigation (sI), which reduced the potential impacts under most of the perturbations. Yet, a combination of adaptations for dealing with climate change demonstrated that effective adaptation is possible at Lleida. Combinations based on a cultivar without vernalization requirements showed good and wide adaptation potential. Few combined adaptation options performed well under rainfed conditions. However, a single sI was sufficient to develop a high adaptation potential, including options mainly based on spring wheat, current cycle duration and early sowing date. Depending on local environment (e.g. soil type), many of these adaptations can maintain current yield levels under moderate changes in T and P, and some also under strong changes. We conclude that ARSs can offer a useful tool for supporting planning of field level adaptation under conditions of high uncertainty.