Abstract: Global warming is expected to increase the frequency and intensity of severe water scarcity (SWS) events, which negatively affect rain-fed crops such as wheat, a key source of calories and protein for humans. Here, we develop a method to simultaneously quantify SWS over the world’s entire wheat-growing area and calculate the probabilities of multiple/sequential SWS events for baseline and future climates. Our projections show that, without climate change mitigation (representative concentration pathway 8.5), up to 60% of the current wheat-growing area will face simultaneous SWS events by the end of this century, compared to 15% today. Climate change stabilization in line with the Paris Agreement would substantially reduce the negative effects, but they would still double between 2041 and 2070 compared to current conditions. Future assessments of production shocks in food security should explicitly include the risk of severe, prolonged, and near- simultaneous droughts across key world wheat-producing areas.

Abstract: Modelling N transformations within cropping systems is crucial for N management optimization in order to increase N use efficiency and reduce N losses. Such modelling remains challenging because of the complexity of N cycling in soil–plant systems. In the current study, the uncertainties of six widely used process-based models (PBMs), including APSIM, CROPSYST, DAISY, FASSET, HERMES and STICS, were tested in simulating different N managements (catch crops (CC) and different N fertilizer rates) in 12-year rotations in Western Europe. Winter wheat, sugar beet and pea were the main crops, and radish was the main CC in the tested systems. Our results showed that PBMs simulated yield, aboveground biomass, N export and N uptake well with low RMSE values, except for sugar beet, which was generally less well parameterized. Moreover, PBMs provided more accurate crop simulations (i.e. N export and N uptake) compared to simulations of soil (N mineralization and soil mineral N (SMN)) and environmental variables (N leaching). The use of multi-model ensemble mean or median of four PBMs significantly reduced the mean absolute percentage error (MAPE) between simulations and observations to less than 15% for yield, aboveground biomass, N export and N uptake. Multi-model ensemble also significantly reduced the MAPE for net N mineralization and annual N leaching to around 15%, while it was larger than 20% for SMN. Generally, PBMs well simulated the CC effects on N fluxes, i.e. increasing N mineralization and reducing N leaching in both short-term and long-term, and all PBMs correctly predicted the effects of the reduced N rate on all measured variables in the study. The uncertainties of multi-model ensemble for N mineralization, SMN and N leaching were larger, mainly because these variables are influenced by plant-soil interactions and subject to cumulative long-term effects in crop rotations, which makes them more difficult to simulate. Large differences existed between individual PBMs due to the differences in formalisms for describing N processes in soil–plant systems, the skills of modelers and the model calibration level. In addition, the model performance also depended on the simulated variables, for instance, HERMES and FASSET performed better for yield and crop biomass, APSIM, DAISY and STICS performed better for N export and N uptake, STICS provided best simulation for SMN and N leaching among the six individual PBMs in the study, but all PBMs met difficulties to well predict either average or variance of soil N mineralization. Our results showed that better calibration for soil N variables is needed to improve model predictions of N cycling in order to optimize N management in crop rotations.

Abstract: Crop models can be sensitive to climate input data aggregation and this response may differ among models. This should be considered when applying field-scale models for assessment of climate change impacts on larger spatial scales or when coupling models across scales. In order to evaluate these effects systematically, an ensemble of ten crop models was run with climate input data on different spatial aggregations ranging from 1, 10, 25, 50 and 100 km horizontal resolution for the state of North Rhine-Westphalia, Germany. Models were minimally calibrated to typical sowing and harvest dates, and crop yields observed in the region, subsequently simulating potential, water-limited and nitrogen-limited production of winter wheat and silage maize for 1982-2011. Outputs were analysed for 19 variables (yield, evapotranspiration, soil organic carbon, etc.). In this study the sensitivity of the individual models and the model ensemble in response to input data aggregation is assessed for crop yield. Results show that the mean yield of the region calculated from climate time series of 1 km horizontal resolution changes only little when using climate input data of higher aggregation levels for most models. However, yield frequency distributions change with aggregation, resembling observed data better with increasing resolution. With few exceptions, these results apply to the two crops and three production situations (potential, water-, nitrogen-limited) and across models including the model ensemble, regardless of differences among models in simulated yield levels and spatial yield patterns. Results of this study improve the confidence of using crop models at varying scales.