Abstract: Crop growth simulation models can differ greatly in their treatment of key processes and hence in their response to environmental conditions. Here, we used an ensemble of 26 process-based wheat models applied at sites across a European transect to compare their sensitivity to changes in temperature (−2 to +9°C) and precipitation (−50 to +50%). Model results were analysed by plotting them as impact response surfaces (IRSs), classifying the IRS patterns of individual model simulations, describing these classes and analysing factors that may explain the major differences in model responses.   The model ensemble was used to simulate yields of winter and spring wheat at sites in Finland, Germany and Spain. Results were plotted as IRSs that show changes in yields relative to the baseline with respect to temperature and precipitation. IRSs of 30-year means and selected extreme years were classified using two approaches describing their pattern.   The expert diagnostic approach (EDA) combines two aspects of IRS patterns: location of the maximum yield (nine classes, Figure 1) and strength of the yield response with respect to climate (four classes), resulting in a total of 36 combined classes defined using criteria pre-specified by experts. The statistical diagnostic approach (SDA) groups IRSs by comparing their pattern and magnitude, without attempting to interpret these features. It applies a hierarchical clustering method, grouping response patterns using a distance metric that combines the spatial correlation and Euclidian distance between IRS pairs. The two approaches were used to investigate whether different patterns of yield response could be related to different properties of the crop models, specifically their genealogy, calibration and process description.   Although no single model property across a large model ensemble was found to explain the integrated yield response to temperature and precipitation perturbations, the application of the EDA and SDA approaches revealed their capability to distinguish: (i) stronger yield responses to precipitation for winter wheat than spring wheat; (ii) differing strengths of response to climate changes for years with anomalous weather conditions compared to period-average conditions; (iii) the influence of site conditions on yield patterns; (iv) similarities in IRS patterns among models with related genealogy; (v) similarities in IRS patterns for models with simpler process descriptions of root growth and water uptake compared to those with more complex descriptions; and (vi) a closer correspondence of IRS patterns in models using partitioning schemes to represent yield formation than in those using a harvest index.   Such results can inform future crop modelling studies that seek to exploit the diversity of multi-model ensembles, by distinguishing ensemble members that span a wide range of responses as well as those that display implausible behaviour or strong mutual similarities.   The full manuscript of this study is currently under revision (Fronzek et al. 2017).

Abstract: Highlights • The larger simulated attainable yield for a specific crop season, the larger the yield gap. • Average size of the yield gap is not affected by the inter-annual variability of attainable yield. • Technology levels (resource input and accessibility) determine average yield gap. • To reduce yield gaps in rainfed environments, farmers need to improve season-specific crop management. Abstract Provision of food security in the face of increasing global food demand requires narrowing of the gap between actual farmer’s yield and maximum attainable yield. So far, assessments of yield gaps have focused on average yield over 5–10 years, but yield gaps can vary substantially between crop seasons. In this study we hypothesized that climate-induced inter-annual yield variability and associated risk is a major barrier for farmers to invest, i.e. increase inputs to narrow the yield gap. We evaluated the importance of inter-annual attainable yield variability for the magnitude of the yield gap by utilizing data for wheat and maize at ten sites representing some major food production systems and a large range of climate and soil conditions across the world. Yield gaps were derived from the difference of simulated attainable yields and regional recorded farmer yields for 1981 to 2010. The size of the yield gap did not correlate with the amplitude of attainable yield variability at a site, but was rather associated with the level of available resources such as labor, fertilizer and plant protection inputs. For the sites in Africa, recorded yield reached only 20% of the attainable yield, while for European, Asian and North American sites it was 56–84%. Most sites showed that the higher the attainable yield of a specific season the larger was the yield gap. This significant relationship indicated that farmers were not able to take advantage of favorable seasonal weather conditions. To reduce yield gaps in the different environments, reliable seasonal weather forecasts would be required to allow farmers to manage each seasonal potential, i.e. overcoming season-specific yield limitations.

Abstract: Highlights • The parameters to be fixed were consistent across sites. • Model calibration must be performed separately for each specific case. • Possible to reduce model parameters from 66 to 45. • Strong model reductions must be avoided. • The error term for the training data were characterised by timing (phase shift). Abstract Proper parameterisation and quantification of model uncertainty are two essential tasks in improvement and assessment of model performance. Bayesian calibration is a method that combines both tasks by quantifying probability distributions for model parameters and outputs. However, the method is rarely applied to complex models because of its high computational demand when used with high-dimensional parameter spaces. We therefore combined Bayesian calibration with sensitivity analysis, using the screening method by Morris (1991), in order to reduce model complexity by fixing parameters to which model output was only weakly sensitive to a nominal value. Further, the robustness of the model with respect to reduction in the number of free parameters were examined according to model discrepancy and output uncertainty. The process-based grassland model BASGRA was examined in the present study on two sites in Norway and in Germany, for two grass species (Phleum pratense and Arrhenatherum elatius). According to this study, a reduction of free model parameters from 66 to 45 was possible. The sensitivity analysis showed that the parameters to be fixed were consistent across sites (which differed in climate and soil conditions), while model calibration had to be performed separately for each combination of site and species. The output uncertainty decreased slightly, but still covered the field observations of aboveground biomass. Considering the training data, the mean square error for both the 66 and the 45 parameter model was dominated by errors in timing (phase shift), whereas no general pattern was found in errors when using the validation data. Stronger model reduction should be avoided, as the error term increased and output uncertainty was underestimated.