Abstract: 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.

Abstract: Croplands and grasslands are agricultural systems that contribute to land–atmosphere exchanges of carbon (C). We evaluated and compared gross primary production (GPP), ecosystem respiration (RECO), net ecosystem exchange (NEE) of CO2, and two derived outputs – C use efficiency (CUE=-NEE/GPP) and C emission intensity (IntC= -NEE/Offtake [grazed or harvested biomass]). The outputs came from 23 models (11 crop-specific, eight grassland-specific, and four models covering both systems) at three cropping sites over several rotations with spring and winter cereals, soybean and rapeseed in Canada, France and India, and two temperate permanent grasslands in France and the United Kingdom. The models were run independently over multi-year simulation periods in five stages (S), either blind with no calibration and initialization data (S1), using historical management and climate for initialization (S2), calibrated against plant data (S3), plant and soil data together (S4), or with the addition of C and N fluxes (S5). Here, we provide a framework to address methodological uncertainties and contextualize results. Most of the models overestimated or underestimated the C fluxes observed during the growing seasons (or the whole years for grasslands), with substantial differences between models. For each simulated variable, changes in the multi-model median (MMM) from S1 to S5 was used as a descriptor of the ensemble performance. Overall, the greatest improvements (MMM approaching the mean of observations) were achieved at S3 or higher calibration stages. For instance, grassland GPP MMM was equal to 1632 g C m−2 yr-1 (S5) while the observed mean was equal to 1763 m-2 yr-1 (average for two sites). Nash-Sutcliffe modelling efficiency coefficients indicated that MMM outperformed individual models in 92.3 % of cases. Our study suggests a cautious use of large-scale, multi-model ensembles to estimate C fluxes in agricultural sites if some site-specific plant and soil observations are available for model calibration. The further development of crop/grassland ensemble modelling will hinge upon the interpretation of results in light of the way models represent the processes underlying C fluxes in complex agricultural systems (grassland and crop rotations including fallow periods).

Abstract: Even brief periods of high temperatures occurring around flowering and during grain filling can severely reduce grain yield in cereals. Recently, ecophysiological and crop models have begun to represent such phenomena. Most models use air temperature (Tair) in their heat stress responses despite evidence that crop canopy temperature (Tc) better explains grain yield losses. Tc can deviate significantly from Tair based on climatic factors and the crop water status. The broad objective of this study was to evaluate whether simulation of Tc improves the ability of crop models to simulate heat stress impacts on wheat under irrigated conditions. Nine process-based models, each using one of three broad approaches (empirical, EMP; energy balance assuming neutral atmospheric stability, EBN; and energy balance correcting for the atmospheric stability conditions, EBSC) to simulate Tc, simulated grain yield under a range of temperature conditions. The models varied widely in their ability to reproduce the measured Tc with the commonly used EBN models performing much worse than either EMP or EBSC. Use of Tc to account for heat stress effects did improve simulations compared to using only Tair to a relatively minor extent, but the models that additionally use Tc on various other processes as well did not have better yield simulations. Models that simulated yield well under heat stress had varying skill in simulating Tc. For example, the EBN models had very poor simulations of Tc but performed very well in simulating grain yield. These results highlight the need to more systematically understand and model heat stress events in wheat.