Future area expansion outweighs increasing drought risk for soybean in Europe.

The European Union is highly dependent on soybean imports from overseas to meet its protein demands. Individual Member States have been quick to declare self-sufficiency targets for plant-based proteins, but detailed strategies are still lacking. Rising global temperatures have painted an image of a bright future for soybean production in Europe, but emerging climatic risks such as drought have so far not been included in any of those outlooks. Here, we present simulations of future soybean production and the most prominent risk factors across Europe using an ensemble of climate and soybean growth models. Projections suggest a substantial increase in potential soybean production area and productivity in Central Europe, while southern European production would become increasingly dependent on supplementary irrigation. Average productivity would rise by 8.3% (RCP 4.5) to 8.7% (RCP 8.5) as a result of improved growing conditions (plant physiology benefiting from rising temperature and CO2 levels) and farmers adapting to them by using cultivars with longer phenological cycles. Suitable production area would rise by 31.4% (RCP 4.5) to 37.7% (RCP 8.5) by the mid-century, contributing considerably more than productivity increase to the production potential for closing the protein gap in Europe. While wet conditions at harvest and incidental cold spells are the current key challenges for extending soybean production, the models and climate data analysis anticipate that drought and heat will become the dominant limitations in the future. Breeding for heat-tolerant and water-efficient genotypes is needed to further improve soybean adaptation to changing climatic conditions.

Decline in climate resilience of European wheat.

Food security relies on the resilience of staple food crops to climatic variability and extremes, but the climate resilience of European wheat is unknown. A diversity of responses to disturbance is considered a key determinant of resilience. The capacity of a sole crop genotype to perform well under climatic variability is limited; therefore, a set of cultivars with diverse responses to weather conditions critical to crop yield is required. Here, we show a decline in the response diversity of wheat in farmers' fields in most European countries after 2002-2009 based on 101,000 cultivar yield observations. Similar responses to weather were identified in cultivar trials among central European countries and southern European countries. A response diversity hotspot appeared in the trials in Slovakia, while response diversity "deserts" were identified in Czechia and Germany and for durum wheat in southern Europe. Positive responses to abundant precipitation were lacking. This assessment suggests that current breeding programs and cultivar selection practices do not sufficiently prepare for climatic uncertainty and variability. Consequently, the demand for climate resilience of staple food crops such as wheat must be better articulated. Assessments and communication of response diversity enable collective learning across supply chains. Increased awareness could foster governance of resilience through research and breeding programs, incentives, and regulation.

Decoupling of impact factors reveals the response of German winter wheat yields to climatic changes.

Yield development of agricultural crops over time is not merely the result of genetic and agronomic factors, but also the outcome of a complex interaction between climatic and site-specific soil conditions. However, the influence of past climatic changes on yield trends remains unclear, particularly under consideration of different soil conditions. In this study, we determine the effects of single agrometeorological factors on the evolution of German winter wheat yields between 1958 and 2015 from 298 published nitrogen (N)-fertilization experiments. For this purpose, we separate climatic from genetic and agronomic yield effects using linear mixed effect models and estimate the climatic influence based on a coefficient of determination for these models. We found earlier occurrence of wheat growth stages, and shortened development phases except for the phase of stem elongation. Agrometeorological factors are defined as climate covariates related to the growth of winter wheat. Our results indicate a general and strong effect of agroclimatic changes on yield development, in particular due to increasing mean temperatures and heat stress events during the grain-filling period. Except for heat stress days with more than 31°C, yields at sites with higher yield potential were less prone to adverse weather effects than at sites with lower yield potential. Our data furthermore reveal that a potential yield levelling, as found for many West-European countries, predominantly occurred at sites with relatively low yield potential and about one decade earlier (mid-1980s) compared to averaged yield data for the whole of Germany. Interestingly, effects related to high precipitation events were less relevant than temperature-related effects and became relevant particularly during the vegetative growth phase. Overall, this study emphasizes the sensitivity of yield productivity to past climatic conditions, under consideration of regional differences, and underlines the necessity of finding adaptation strategies for food production under ongoing and expected climate change.

Modelling climate change impacts on maize yields under low nitrogen input conditions in sub-Saharan Africa.

Smallholder farmers in sub-Saharan Africa (SSA) currently grow rainfed maize with limited inputs including fertilizer. Climate change may exacerbate current production constraints. Crop models can help quantify the potential impact of climate change on maize yields, but a comprehensive multimodel assessment of simulation accuracy and uncertainty in these low-input systems is currently lacking. We evaluated the impact of varying [CO2 ], temperature and rainfall conditions on maize yield, for different nitrogen (N) inputs (0, 80, 160 kg N/ha) for five environments in SSA, including cool subhumid Ethiopia, cool semi-arid Rwanda, hot subhumid Ghana and hot semi-arid Mali and Benin using an ensemble of 25 maize models. Models were calibrated with measured grain yield, plant biomass, plant N, leaf area index, harvest index and in-season soil water content from 2-year experiments in each country to assess their ability to simulate observed yield. Simulated responses to climate change factors were explored and compared between models. Calibrated models reproduced measured grain yield variations well with average relative root mean square error of 26%, although uncertainty in model prediction was substantial (CV = 28%). Model ensembles gave greater accuracy than any model taken at random. Nitrogen fertilization controlled the response to variations in [CO2 ], temperature and rainfall. Without N fertilizer input, maize (a) benefited less from an increase in atmospheric [CO2 ]; (b) was less affected by higher temperature or decreasing rainfall; and (c) was more affected by increased rainfall because N leaching was more critical. The model intercomparison revealed that simulation of daily soil N supply and N leaching plays a crucial role in simulating climate change impacts for low-input systems. Climate change and N input interactions have strong implications for the design of robust adaptation approaches across SSA, because the impact of climate change in low input systems will be modified if farmers intensify maize production with balanced nutrient management.

Climate change impact and adaptation for wheat protein.

Wheat grain protein concentration is an important determinant of wheat quality for human nutrition that is often overlooked in efforts to improve crop production. We tested and applied a 32-multi-model ensemble to simulate global wheat yield and quality in a changing climate. Potential benefits of elevated atmospheric CO2 concentration by 2050 on global wheat grain and protein yield are likely to be negated by impacts from rising temperature and changes in rainfall, but with considerable disparities between regions. Grain and protein yields are expected to be lower and more variable in most low-rainfall regions, with nitrogen availability limiting growth stimulus from elevated CO2 . Introducing genotypes adapted to warmer temperatures (and also considering changes in CO2 and rainfall) could boost global wheat yield by 7% and protein yield by 2%, but grain protein concentration would be reduced by -1.1 percentage points, representing a relative change of -8.6%. Climate change adaptations that benefit grain yield are not always positive for grain quality, putting additional pressure on global wheat production.

Contribution of crop model structure, parameters and climate projections to uncertainty in climate change impact assessments.

Climate change impact assessments are plagued with uncertainties from many sources, such as climate projections or the inadequacies in structure and parameters of the impact model. Previous studies tried to account for the uncertainty from one or two of these. Here, we developed a triple-ensemble probabilistic assessment using seven crop models, multiple sets of model parameters and eight contrasting climate projections together to comprehensively account for uncertainties from these three important sources. We demonstrated the approach in assessing climate change impact on barley growth and yield at Jokioinen, Finland in the Boreal climatic zone and Lleida, Spain in the Mediterranean climatic zone, for the 2050s. We further quantified and compared the contribution of crop model structure, crop model parameters and climate projections to the total variance of ensemble output using Analysis of Variance (ANOVA). Based on the triple-ensemble probabilistic assessment, the median of simulated yield change was -4% and +16%, and the probability of decreasing yield was 63% and 31% in the 2050s, at Jokioinen and Lleida, respectively, relative to 1981-2010. The contribution of crop model structure to the total variance of ensemble output was larger than that from downscaled climate projections and model parameters. The relative contribution of crop model parameters and downscaled climate projections to the total variance of ensemble output varied greatly among the seven crop models and between the two sites. The contribution of downscaled climate projections was on average larger than that of crop model parameters. This information on the uncertainty from different sources can be quite useful for model users to decide where to put the most effort when preparing or choosing models or parameters for impact analyses. We concluded that the triple-ensemble probabilistic approach that accounts for the uncertainties from multiple important sources provide more comprehensive information for quantifying uncertainties in climate change impact assessments as compared to the conventional approaches that are deterministic or only account for the uncertainties from one or two of the uncertainty sources.

Multimodel ensembles of wheat growth: many models are better than one.

Crop models of crop growth are increasingly used to quantify the impact of global changes due to climate or crop management. Therefore, accuracy of simulation results is a major concern. Studies with ensembles of crop models can give valuable information about model accuracy and uncertainty, but such studies are difficult to organize and have only recently begun. We report on the largest ensemble study to date, of 27 wheat models tested in four contrasting locations for their accuracy in simulating multiple crop growth and yield variables. The relative error averaged over models was 24-38% for the different end-of-season variables including grain yield (GY) and grain protein concentration (GPC). There was little relation between error of a model for GY or GPC and error for in-season variables. Thus, most models did not arrive at accurate simulations of GY and GPC by accurately simulating preceding growth dynamics. Ensemble simulations, taking either the mean (e-mean) or median (e-median) of simulated values, gave better estimates than any individual model when all variables were considered. Compared to individual models, e-median ranked first in simulating measured GY and third in GPC. The error of e-mean and e-median declined with an increasing number of ensemble members, with little decrease beyond 10 models. We conclude that multimodel ensembles can be used to create new estimators with improved accuracy and consistency in simulating growth dynamics. We argue that these results are applicable to other crop species, and hypothesize that they apply more generally to ecological system models.

Redistribution of nitrogen to feed the people on a safer planet.

Lack of nitrogen limits food production in poor countries while excessive nitrogen use in industrial countries has led to transgression of the planetary boundary. However, the potential of spatial redistribution of nitrogen input for food security when returning to the safe boundary has not been quantified in a robust manner. Using an emulator of a global gridded crop model ensemble, we found that redistribution of current nitrogen input to major cereals among countries can double production in the most food-insecure countries, while increasing global production of these crops by 12% with no notable regional loss or reducing the nitrogen input to the current production by one-third. Redistribution of the input within the boundary increased production by 6-8% compared to the current relative distribution, increasing production in the food-insecure countries by two-thirds. Our findings provide georeferenced guidelines for redistributing nitrogen use to enhance food security while safeguarding the planet.

Understanding spatial and temporal variability of N leaching reduction by winter cover crops under climate change.

Winter cover crops are sown in between main spring crops (e.g. cash and forage crops) to provide a range of benefits, including the reduction of nitrogen (N) leaching losses to groundwater. However, the extent by which winter cover crops will remain effective under future climate change is unclear. We assess variability and uncertainty of climate change effects on the reduction of N leaching by winter oat cover crops. Field data were collected to quantify ranges of cover crop above-ground biomass (7 to 10 t DM/ha) and N uptake (70 to 180 kg N/ha) under contrasting initial soil conditions. The data were also used to evaluate the APSIM-NextGen model (R2 from 62 to 96% and RMSEr from 7 to 50%), which was then applied to simulate cover crop and fallow conditions across four key agricultural locations in New Zealand, under baseline and future climate scenarios. Cover crops reduced N leaching risks for all location/scenario combinations but with large variability in space and time (e.g. 21 to 47% of fallow) depending on the climate change scenario. For instance, end-of-century estimates for northern (warmer) locations mostly showed non-significant effects of climate change on cover crop effectiveness and N leaching. In contrast for southern (colder) locations, there was a systematic increase in N leaching risks with climate change intensity despite a concomitant, but less than proportional, increase in cover crop effectiveness (up to ~5% of baseline) due to higher winter yields and N uptake. This implies that climate change may not only modify the geography of N leaching hotspots, but also the extent by which cover crops can locally reduce pollution risks, in some cases requiring complementary adaptive measures. The patchy- and threshold-nature of leaching events indicates that fine spatio-temporal resolutions are better suited to evaluate cover crop effectiveness under climate change.

Water footprint of winter wheat under climate change: Trends and uncertainties associated to the ensemble of crop models.

The impact of climate change could undermine the future grain production as a consequence of increased temperature and drought condition or improve the crop performance owing to the increased CO2 in the atmosphere. Wheat water demand and yield are strictly related to climate conditions of the area where the plants are cropped. In this study, we assessed the future trends of grain yield and water consumption in two European regions, Germany (Continental region) and Italy (Mediterranean region) in the light of the multiple sources of uncertainty related to climate and yield forecasts. Four crop models were set up under combinations of two European climate regions, five Global Circulation Models and two Representative CO2 Concentration Pathways, 486 ppm and 540 ppm in 2050. Yield and water use were assessed under rainfed and irrigated regimes, and the water footprint of green water and total water was estimated. Our results indicated that projected yields were comparable (Mediterranean area) or even improved (+9%; Continental area) in rainfed conditions in comparison to the current trend; and water supply enhanced crop performance (+22% in Germany and +19% in Italy, as mean). Crop water consumption (both green and blue) remained stable in future projections but the water footprint was 5% lower on average in Italy and 23% in Germany when compared to the baseline. Despite the uncertainty in future predictions related to the factors analysed, our result indicated that current wheat production and its water footprint could become more favourable under climate change.

Mitigation efforts will not fully alleviate the increase in water scarcity occurrence probability in wheat-producing areas.

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.

Diverging importance of drought stress for maize and winter wheat in Europe.

Understanding the drivers of yield levels under climate change is required to support adaptation planning and respond to changing production risks. This study uses an ensemble of crop models applied on a spatial grid to quantify the contributions of various climatic drivers to past yield variability in grain maize and winter wheat of European cropping systems (1984-2009) and drivers of climate change impacts to 2050. Results reveal that for the current genotypes and mix of irrigated and rainfed production, climate change would lead to yield losses for grain maize and gains for winter wheat. Across Europe, on average heat stress does not increase for either crop in rainfed systems, while drought stress intensifies for maize only. In low-yielding years, drought stress persists as the main driver of losses for both crops, with elevated CO2 offering no yield benefit in these years.

Impact of Spatial Soil and Climate Input Data Aggregation on Regional Yield Simulations.

We show the error in water-limited yields simulated by crop models which is associated with spatially aggregated soil and climate input data. Crop simulations at large scales (regional, national, continental) frequently use input data of low resolution. Therefore, climate and soil data are often generated via averaging and sampling by area majority. This may bias simulated yields at large scales, varying largely across models. Thus, we evaluated the error associated with spatially aggregated soil and climate data for 14 crop models. Yields of winter wheat and silage maize were simulated under water-limited production conditions. We calculated this error from crop yields simulated at spatial resolutions from 1 to 100 km for the state of North Rhine-Westphalia, Germany. Most models showed yields biased by <15% when aggregating only soil data. The relative mean absolute error (rMAE) of most models using aggregated soil data was in the range or larger than the inter-annual or inter-model variability in yields. This error increased further when both climate and soil data were aggregated. Distinct error patterns indicate that the rMAE may be estimated from few soil variables. Illustrating the range of these aggregation effects across models, this study is a first step towards an ex-ante assessment of aggregation errors in large-scale simulations.