# Causes and implications of the unforeseen 2016 extreme yield loss in the breadbasket of France

Tamara Ben-Ari, Julien Boé, Philippe Ciais, Rémi Lecerf, Marijn van Der Velde, David Makowski

► To cite this version:

Tamara Ben-Ari, Julien Boé, Philippe Ciais, Rémi Lecerf, Marijn van Der Velde, et al. Causes and implications of the unforeseen 2016 extreme yield loss in the breadbasket of France. Nature Communications, 2018, 9 (1), pp.1627. 10.1038/s41467-018-04087-x. hal-01806699

HAL Id: hal-01806699

https://hal.science/hal-01806699

Submitted on 17 Sep 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
ARTICLE

DOI: 10.1038/s41467-018-04087-x

OPEN

# Causes and implications of the unforeseen 2016 extreme yield loss in the breadbasket of France

Tamara Ben-Ari<fnref n="1" />, Julien Boé<fnref n="1" />,², Philippe Ciais³, Remi Lecerf⁴, Marijn Van der Velde⁴ & David Makowski<fnref n="1" />

In 2016, France, one of the leading wheat-producing and wheat-exporting regions in the world suffered its most extreme yield loss in over half a century. Yet, yield forecasting systems failed to anticipate this event. We show that this unprecedented event is a new type of compound extreme with a conjunction of abnormally warm temperatures in late autumn and abnormally wet conditions in the following spring. A binomial logistic regression accounting for fall and spring conditions is able to capture key yield loss events since 1959. Based on climate projections, we show that the conditions that led to the 2016 wheat yield loss are projected to become more frequent in the future. The increased likelihood of such compound extreme events poses a challenge: farming systems and yield forecasting systems, which often support them, must adapt.

<footnote n="1">¹INRA, AgroParisTech, Université Paris-Saclay, UMR 211 Agronomie, Thiverval-Grignon 78850, France. ²CECI, Université de Toulouse, CERFACS/CNRS, Toulouse 31100, France. ³Laboratoire des Sciences du Climat et de l'Environnement, 91191 Gif-sur-Yvette, France. ⁴European Commission, Joint Research Centre (JRC), Via E. Fermi 2749, Ispra, VA 21027, Italy. Correspondence and requests for materials should be addressed to T.B-A. (email: tamara.ben-ari@inra.fr)</footnote>
Crop yield forecasting systems rely on a combination of expert knowledge, data mining and analysis, and mechanistic and/or statistical models¹⁻³. The increased extreme weather conditions in recent years⁴, have challenged these systems. Heat waves and drought events take the largest toll on production⁵, and are anticipated to increase in frequency and severity in the northern mid-latitudes⁶. Extreme yield losses also can occur, however, due to events resulting from insidious constellations of climate variables forming a compound extreme. The 2016 extreme loss of wheat harvest in the breadbasket region of France is one such example, and we present here an in-depth analysis of this event and its implications for wheat yield forecasting.

France ranks fifth in the global league table of national wheat production, and despite its limited arable area, produces more wheat than any other country in the European Union (EU). This is achieved because of very high yields. Recent French yields equal about 7.4 t ha⁻¹. In comparison, China, India, Russia and the United States, the world's four largest wheat producers, harvest about 5, 2.5, and 3 t ha⁻¹, respectively⁷. Between 2000 and 2013, France was the EU's main grain exporter, exporting about 17 million tonnes of wheat mainly to North Africa, where local production covers only 10–50% of the demand⁷.

The 2016 winter-wheat harvest was disastrous. Yields in the breadbasket region dropped on average by 27.7% compared to trend expectations (Fig. 1) and by 39.5% compared to 2015. This equates to a shortfall of about 8 million tonnes compared to the 24.5 million tonnes usually harvested in this region or of about 11 million tonnes compared to the record 27.5 tonnes harvested in 2015⁸. These extremely low yields combined with lower exchange prices on international markets compared to 2015 induced a substantial income loss for farmers and about 2.3 billion dollars loss for France's trade balance⁹.

None of the public forecasting systems anticipated the magnitude of this loss. Even just before the disastrous harvest, forecasts predicted average yields of 7.23 t ha⁻¹, close to the 5-year average¹⁰, overestimating the actual value by about 2 t ha⁻¹. Towards the end of the growing season, there were concerns among regional experts about heavy rainfalls leading to flooding

and saturated agricultural soils in the Seine river basin, and about high incidences of foliar diseases. High observed wheat biomass at the very end of winter on the one hand and a strong confidence on the effectiveness of fungicides on the other, possibly explain that close-to-average yield values were anticipated by most experts until harvest started unfolding in mid-July. Although technical experts have since investigated the possible causes of this extreme yield loss, no quantitative study has characterized the precise climatic conditions that led to it—we still have little understanding of why yield forecasts failed by such a large amount in 2016. Here, we analyse long-term yield and climate time series at the scale of departments (administrative units) from 1959 to 2016 and address the following research questions: (i) how exceptional were climate conditions, individually and combined, in the breadbasket region during the 2015–2016 growing season? (ii) Can 2016 help us improve forecasting systems? And (iii) will such events become more frequent in the future?

We first search for single and compound climatic extremes that occurred during the 2015–2016 growing season. We then perform a statistical analysis to model severe and extreme wheat yield losses since 1959 based upon these climate variables. We show that the huge wheat yield loss in 2015–2016 can be predicted from a combination of climate variables related to higher temperature in autumn and wetter conditions in spring. Finally, based on future climate warming projections, we show that the specific conditions that led to the 2016 wheat yield loss are projected to become more frequent towards the end of the century.

## Results

**Extreme loss and unprecedented weather conditions.** We focus this study on the breadbasket region, comprising 27 departments, and which together accounts for more than 67% of France's total wheat production (average since 1959). In 2015, these departments for example produced more wheat than all of Ukraine and slightly less than all of Canada. All 27 departments suffered extreme yield losses in 2016, leading to the spatio-temporal

<figure><img src="image_2.png" /><figcaption>**Fig. 1** Spatio-temporal pattern of the 2016 extreme yield loss. **a** Wheat yield anomalies in 2016 relative to expected values defined in each department by the long-term yield trend (e.g., −0.1 corresponds to a loss of 10% compared to expectation, see Methods). The breadbasket region of France is delineated in bold black contours. Note that a similar figure is presented for each year in the data set in the Supplementary Movie 1. The map was generated with R based on the yield data used in the analyses. **b** Boxplot of the distribution of anomalies in the breadbasket (1959–2016). Decades are separated with thin dotted lines. Anomalies in 2016 are highlighted in red. Yearly median anomalies are presented in Supplementary Fig. 1b</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Fig. 2** Maximum temperatures and precipitation over the 1959–2016 wheat growing seasons. Boxplot of **a** daily maximum temperature and **b** daily precipitation averaged each year over the study area for the period 1959–2015. Whiskers extend to maximum and minimum values. Values in the 2015–2016 growing season are presented as red dots (for other climate variables see Supplementary Fig. 2)</figcaption></figure>

pattern shown in Fig. 1 and Supplementary Fig. 1. Before 2016, yields never dropped by more than 15% of the long-term mean, on average over the study area, except during the extreme drought of 1976 when the loss was about 16%. By contrast, in 2016, yields were between 17 and 45% below their expected values. The year 2016 thus suffered the single most extreme loss experienced over the past five decades (Fig. 1b, Supplementary Fig. 1b, Supplementary Movie 1).

To characterize climate conditions in the breadbasket region, we present monthly and seasonal distributions for seven growing-season climate variables from the autumn of 1958 to the 2016 harvest using data from the SAFRAN analysis¹¹. The 2015–2016 growing season is singled out and indicated by red dots (Fig. 2, Supplementary Fig. 2). Both maximum and minimum temperatures in the late autumn (here November and December) of 2015 were exceptionally high (Fig. 2a, Supplementary Fig. 2a), resulting in a dramatic reduction of the number of vernalizing days (see Supplementary Fig. 2b). Vernalization is a critical process that controls the development of wheat through exposure to cold temperatures. Temperatures were also high in January and February of 2016 (i.e., close to the third quartile) and, overall, the winter of 2016 was relatively wet compared to the average over 1959–2015. The data in Fig. 2 also show that the spring (here April to July inclusive) of 2016 was extremely wet, with mean precipitation of 2.66 mm day⁻¹, conditions associated with abnormally low potential evapotranspiration (in particular in May and June see Supplementary Fig. 2c). We found only three other years with whole spring conditions within 10% of spring 2016 values (i.e., 1983, 2000, and 2012). Note that seasonal values hide between-months variability. In April 2016, precipitation was close to average, but in May it reached a record high of 4.39 mm day⁻¹, which had never occurred since 1959. June 2016, the month preceding harvest, was also characterized by the persistence of high precipitation (Fig. 2b, Supplementary Fig. 2d). Daily average solar radiation was low during most of the spring in 2016, with a record low value of 160 W m⁻² in June 2016 (Supplementary Fig. 2e).

Overall, the 2015–2016 growing season was characterized by a unique combination of abnormally warm temperatures in the late autumn and abnormally high precipitation, with concurrent low radiation and potential evapotranspiration, in the spring. These variables were outside the 95th percentile of their distributions. The joint probability of all the 2015–2016 growing season climate conditions was null, which makes this event a compound extreme.

**Relating yield loss to autumn and spring conditions.** We designed an ensemble of statistical models to diagnose yield loss occurrence as a function of time series of climate variables in each department since autumn 1958. We considered probabilities of yield losses below −10 and −15% (respectively corresponding to a loss below one standard deviation and below the 10th percentile). In the rest of the manuscript we refer to these levels as severe and extreme yield losses and present results for net yield losses (i.e., negative yield anomalies) in the supplement.

The influential parsimonious subset of climate variables was selected using the Bayesian Information Criterion (BIC) independently for each yield loss level. We considered both monthly and seasonal climate variables with autumn defined as October–December, winter as January–March and spring as April–July, inclusive. We first computed BIC for each variable independently and select a subset of four best covariates relying on the extreme climate events characterizing the 2015–2016 growing season. Climate variables in the selected three best models are consistent across the levels of loss considered. These are: (1) the number of days with maximum temperatures between 0 and 10 °C in December, (2) precipitation in November, (3) precipitation in spring (or in May), and (4) minimum temperature in June. Interactions among these variables are considered (see full description of each model in Supplementary Table 2). The same variables are selected when the models are fitted to normalized anomalies (See Supplementary discussion and Supplementary Fig. 3).
<figure><img src="image_4.png" /><figcaption>**Fig. 3** Modeled probabilities. Estimated probabilities of **a** severe and **b** extreme yield loss (yields more than 10 and 15% lower than the expected value) from generalized linear regression models based on climate predictors fitted to the full yield time series (1959–2016) in the breadbasket region. The dotted blue line is the prior probability (the empirical proportion of severe or extreme yield loss events in the dataset). As an indication, vertical red lines correspond to median yield loss occurrences across all departments larger than –15% (i.e., in 1976, 2003, and 2016). Note that the figure for net yield losses is presented in the Supplementary Fig. 3</figcaption></figure>

The influential climate variables statistically associated with yield losses over 1959–2016 include those identified by the French academy of agriculture, extension services or specialized journalists after the 2016 harvest (summarized in Supplementary Table 1). The selected best models, combining the four-abovementioned variables and between-seasons interactions, trained on the entire data set, are used to assess the risk of net, severe, and extreme yield loss in the breadbasket over all harvest years (Fig. 3, Supplementary Fig. 4).

The models successfully capture a series of key yield-loss events, e.g., extreme losses in 1976, 2003, and 2016 and severe losses in 1966, 1970, and 2007. As an example, the estimated probability of extreme yield loss reaches a median (10th–90th percentiles) across departments of 0.39 (0.21–0.74) during the severe drought year of 1976, the second worst yield on record and of 0.33 (0.15–0.59) for the heatwaves of 2003 (Fig. 3). Note that probabilities computed by our statistical models need to be interpreted as a departure from prior probabilities defined as the frequency of actual losses in the data set (Supplementary Table 3). For example, the prior value for extreme losses is 0.081, computed probabilities thus correspond to an increase of roughly 5 and 4-folds in 1976 and 2003 respectively. Notably, the years 1976, 2003, and 2016 are highly contrasted in terms of autumn and spring conditions (Supplementary Fig. 5).

We find that the selected statistical models all include an interaction between the number of days between 0 and 10 °C in December and spring precipitation (p = 0.0027, Supplementary Table 2) suggesting that the strength of the relationship between high precipitation in the spring and the probability of yield loss increases with temperatures in the late autumn. For example, according to our model, if the number of days between 0 and 10 °C in December drops from 20 to 10 and the following spring is characterized by average precipitation levels, the probability of severe yield loss increases from 0.12 to 0.2. But when precipitation in the following spring is above one standard deviation this probability increases from 0.2 to 0.51 (Fig. 4a). In other words, the effect of excess precipitation in the spring is stronger if it follows warmth in late autumn. Such an effect is also observed for extreme (Fig. 4b) or net losses (Supplementary

Fig. 6). Selected models also include an interaction between temperatures in June (i.e., flowering period) and precipitation over the spring for all yield loss levels (see discussion section).

Our results illustrate the importance of taking into account the joint occurrence of multiple yet specific climate predictors during the growing season¹² if we are to successfully predict the impacts of compounded extremes like the one of 2016. Although we did not attempt to use our models for neighboring countries, we note that similar weather anomalies in autumn and spring were recorded in Belgium, which also suffered an extreme yield loss of about 24%¹³. Our results could thus probably be extrapolated to other similar agro-ecological areas, provided wheat-cropping systems use similar varieties and agricultural practices.

**On the prediction of the 2016 extreme yield loss.** Figure 5 presents odds ratios computed from estimated probabilities of severe and extreme yield loss over each unit of the breadbasket from logistic models trained in the study area on a dataset excluding 2016 (out-of-sample procedure). Odds ratios indicate the relative chance of severe or extreme yield losses. According to our statistical models, the odds of an extreme yield loss in 2016 were 35 times higher than expected (i.e., from prior values, Supplementary Table 3). This is equivalent to a risk ratio of about 11. In other words, our statistical models estimates a probability of extreme yield loss 11 times higher than a priori expected. Our models also predict between 1.8 and 4.6 more chance of losing yield severely in 2016 than in an average year. Estimated probabilities of severe and extreme yield loss in 2016—based on two separate models excluding 2016 from the training datasets—are on average of 0.46 and 0.71, respectively (Supplementary Table 3). Note that these values hide important local disparities (Supplementary Fig. 7); both the confidence intervals and the inter-unit ranges are larger for extreme yield losses reflecting the smaller number of occurrences of such events in our data set. The results obtained are robust to a change in the training and test data sets (see Methods section).
<figure><img src="image_5.png" /><figcaption>**Fig. 4** Autumn-spring interaction. Modulation by spring precipitation of the relationship between the number of days between 0 and 10 °C in December and the probability of severe (a) and extreme wheat yield loss (b). Interactions are shown here for normal precipitation (daily AMJJ precipitation equal to the median over 1959–2016; bold line), and for wet and dry spring conditions (daily AMJJ precipitation equal to +/− one standard deviation; dotted lines). Vertical grey line: median number of days with maximum temperature between 0 and 10 °C over 1958–2016; red dotted lines: modeled value for 2016</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 5** Odds of yield loss in 2016. Odds ratios of severe (a) and extreme (b) yield loss in the study area from the best models trained on yield data excluding 2016. Data are shown in Supplementary Table 3. Models include temperature and precipitation conditions in the late autumn and spring. We rely on the odds ratio classification by⁶⁸ with values from 1 to 3.2 not worth more than a bare mention, between 3.2 and 10 suggesting substantial evidence, from 10 to 100 strong and above 100 decisive. The maps were generated with R based on the yield data used in the analyses</figcaption></figure>

Thus, providing that all weather data were available in time, the use of our statistical model before the 2016 harvest would have indicated strong evidence of a severe to extreme yield loss. Unlike an actual forecast, however, the influential set of climate covariates were selected from a data set including the 2015–2016 growing season. The implication of this result lies in the key information provided by 2016 extreme wheat yield loss: outlining the necessity of considering late-autumn climatic conditions for winter wheat forecasting.

**What 2016 presages for the future.** Based on climate projections from the Coupled Model Intercomparison Project phase 5 (CMIP5), we assessed how regional climate change will impact the likelihood of climate conditions similar to the ones experienced during the 2015–2016 wheat growing season. All four

climate variables identified as influential were extracted from model outputs. Very small changes in spring and November or spring precipitation over northern France are projected under the RCP2.6 scenario, i.e., the frequency of 2016-like precipitation anomalies hardly changes over the twenty-first century (Fig. 6a, b). On the other hand, as a result of the warming trend under the RCP2.6 scenario, the number of days with Tmax between 0 and 10 °C in December is projected to decline by on average 5 days (around 22%) by the end of the twenty-first century compared to the 1950s (Fig. 6c). The likelihood of a temperature anomaly at least as severe as during the 2015–2016 growing season is extremely small in the 1950s (~2%), slightly greater in the current decade (~6%), and increases moderately by the end of the twenty-first century (~12%). The ensemble-mean June minimum temperature is projected to increase by 1.5 °C during the twenty-first century under the RCP2.6 scenario (Fig. 6d). The relative
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 6** Time evolution towards 2100 of the four climate variables identified as key in our statistical models. **a** Precipitation in spring (mm per day), **b** Precipitation in November (mm per day), **c** Number of days with Tmax between 0 and 10 °C and **d** Minimum temperature in June (°C). The anomalies relative to the 1959–1988 reference period are shown in Supplementary Fig. 8. The blue lines show the multi-model signal derived from an ensemble of 13 CMIP5 models and estimated by locally weighted scatterplot smoothing (loess) (solid line: RCP2.6 scenario. Dashed line: RCP8.5 scenario).</figcaption></figure>

SAFRAN reanalysis data are shown with dark grey lines (thick dashed line: loess results, thin line: yearly values). The light (dark) blue shaded areas show the 20-80% (5-95%) range of the model outputs distribution assuming the independence of models, for the historical and RCP26 simulations. The associated delimiting quantiles are estimated with a local polynomial quantile regression. The grey points show the yearly model outputs

cumulative frequency of the positive June 2016 minimum temperature anomaly in the CMIP5 ensemble decreases from around 90% in 1951–1980 to 70% in the 30-year period centered on 2016 and 50% at the end of the twenty-first century.

Despite the near-absence of projected change in springtime precipitation (here including July), given the strong projected twenty-first century increase in warm minimum temperature anomalies in June and the decrease in the number of days between 0 and 10 °C in December, climate conditions favorable to yield loss are projected to become more frequent under the RCP2.6 scenario. Our results also suggest that the warming trend observed in France over the last several decades—and partly attributable to anthropogenic forcings¹⁴,¹⁵—have already increased the probability of 2016-like climate conditions occurring.

The potential benefits of the aggressive mitigation policies required to follow a low warming scenario like RCP2.6 are obvious when the previous results are compared to the ones obtained with the RCP8.5 scenario (dashed blue lines in Fig. 6—see inter-model spread in Supplementary Fig. 8). In this intensive

warming scenario, a number of vernalizing days as small as in December 2015 becomes the norm by 2070 and, even if associated with a slight drying of spring months, would drastically increase the probability of a 2016-like event. Also, June 2016 minimum temperature, unusually warm relatively to the mid-twentieth century climatology, would become characteristic of an extremely cold June by the end of the twenty-first century (Supplementary Fig. 8c).

Note that we did not address the possible effects on wheat yields of an increase in atmospheric CO₂ concentration. CO₂ effects are expected to manifest through increased leaf-level carboxylation rates and stomatal closure. Both processes interact with each other, and stomatal closure has the dual consequence of saving soil water and increasing surface temperature by reducing transpiration. It is thus expected that elevated CO₂ would improve yields under dry condition¹⁶,¹⁷.

To summarize, there are two key climatic factors associated with the 2016 loss: autumn temperature for which we know with great confidence that 2016 is going to move closer to an average year and spring precipitation for which we do not detect a
noticeable time trend. Qualitatively, 2016-like years are very likely to become more frequent in the future since one of the two extreme factors will no longer be extreme anymore. However, a robust quantification of future yield loss probability from a 2016-like growing season is beyond the scope of our study and possibly hampered by the high uncertainty associated with projected changes in precipitation.

Discussion

The experts who have been analysing the possible cause of the extreme yield loss, all recognized a posteriori the abnormal pre-cipitation and radiation conditions in the spring of 2016 and, some of them, the warmth of the preceding late autumn. Cited mechanisms most often include lodging, and more widespread occurrence of pests and diseases (Supplementary Table 1).

A decrease in the number of days between 0 and 10 °C in December (i.e., vernalizing days¹⁸)—corresponding to warmth in a usually cold month—may have direct and/or indirect effects on wheat yield. The lack of sufficient vernalizing days can modify subsequent phenological development¹⁸. Vernalization for example affects the number of leaves and tillers, floral initiation time or flowering phenology¹⁹. Warmth during the vernalization period can delay the onset of the reproductive stage increasing the risk of exposure to high temperatures during anthesis²⁰. This effect alone is unlikely, however, to explain the loss observed, at least in 2016, because most soft wheat varieties cultivated in France require only about 40 vernalizing days²¹, and this requirement was still fulfilled by the end of the winter, despite warm temperatures (Supplementary Fig. 2b). This is corroborated by the independent observation that spring varieties—with no vernalization requirement—were also strongly affected in 2016 (i.e., about -20% at the national scale⁸). But, winter warmth is also known to shift the phenology of pests and diseases, causing earlier colonization of crops and earlier spreading of vector-borne viruses²²⁻²⁴ resulting in more frequent and severe infections²⁵. In addition, positive precipitation anomalies or persistent moisture facilitates the development and spread of fungal diseases in the spring²⁵. A plausible hypothesis to our findings is that the combination of a mild autumn/winter favors a build-up of parasites and a persistence of inoculum, which subsequently leads to large-scale disease prevalence in the field provided conducive spring precipitation conditions would occur²⁶. 2016 was marked by very high to abnormal precipitation levels in the spring; conditions indeed favorable for the spread of diseases. Moreover, localized extreme precipitation events—a phenomenon observed in the field in 2016—may also induce flooding which subsequently leads to anoxia and lodging²⁷,²⁸.

The 2016 wheat harvest assessment revealed lower grain numbers and very low grain weight, suggesting impaired grain filling²⁹. Minimum temperatures around June determine the length of the grain-filling period, with higher temperatures lowering kernel weights³⁰,³¹. In 2016, the minimum temperatures in June were abnormally high (Supplementary Fig. 2a). In our models, minimum June temperature indeed appears to be an important co-factor influencing yield loss. Additionally, we find a significant interaction between June temperatures and spring precipitation: the latter modulate the probability of yield loss from June heat. We find that increased spring precipitation can slightly downplay the importance of high minimum temperatures in June in the models whereas dry conditions increase loss probability from high temperatures in June (e.g., for net yield losses in Supplementary Fig. 6). This interaction more generally reflects the impact of heat and water stress on photosynthetic activity. Finally, and consistent with earlier analyses for wheat and maize yield³²,³³, the best statistical models also found a small

positive relationship between November precipitation (i.e., between sowing and emergence) and yield loss. This could suggest a negative effect of waterlogging on root growth³⁴ and/or point to enhanced survival or growth of soil-borne diseases³⁵. Note that November 2015 precipitation was close to the 50-year average (Fig. 2b).

The failure of yield forecasts in 2016 needs to be understood in the context of a difficulty in simulating winter crops compared to spring crops³⁶ yields. This is perhaps because of the wide range of growth drivers and limiting factors in wheat³³,³⁷, whereas maize yields are for example more evidently driven by water and heat stresses³⁸. An additional deficiency of deterministic wheat forecasts is the difficulty in simulating development stages in coincidence with climate events. A common example is the impact of heat stress during³⁹⁻⁴¹ or after⁴² anthesis. Complex and localized phenomena such as flooding, lodging, or the prevalence of pests and diseases, which can take a large toll on production, are ignored in both process-based and statistical models⁴³. Spatially-explicit reliable information of initial soil water conditions, rooting depth or soil drainage and soil water-holding capacity⁴⁴⁻⁴⁷ should be included in assessing the risk-benefit balance of wet years such as 2016. Models could also benefit from region-specific parameterization of agro-management practices or onset-adaptation strategies. To overcome the shortcomings of using a single crop model, the use of model ensembles is arguably the way forward⁴⁸,⁴⁹, but this strategy is not yet routine in seasonal crop forecasting and would imply implementing the above-mentioned mechanisms. Finally, there is readily available information that could be harnessed to improve forecast systems. A regional plant health bulletin, for example, is published each month in France⁵⁰, and hydrological anomalies are regularly updated⁵¹. These observational data could complement forecast estimates by improving an analyst's judgment. Turning available information from local sources into harmonized data sets at regional scale, updated on an operational basis should probably be a key priority to improve wheat yield forecasts in Europe. Early-warning procedures making real-time information available to farmers in exchange for targeted field observations could also help improve forecasts (e.g., early yield estimates collected during harvest through social media⁵²).

Other abnormal climate conditions have affected primary production in the past. These most often occurred during droughts: for example, the 2003 heatwave, which caused tremendous damage to vegetation in France⁵³, had an enormous societal impact⁵⁴. This was also true in 1976, a growing season with similar limiting factors. Those visible impacts of water and heat stresses in the spring may actually have hidden other multivariate climatic events with similar or higher impacts on wheat harvests (e.g., in 1970, 1987, 2007, or 2016).

Based on long-term wheat yield and the department-scale climate time series, we show that the compound interaction between temperature in the late autumn/early winter and precipitation in the spring is the key to understanding the severe yield drop of 2016. A series of red flags were identified that could have enabled experts to anticipate the event. Our results also show that probabilistic approaches can be very helpful for anticipating yield losses provided that they are properly trained on the right combination of climate variables. Depending on the nature of decision makers' demand, crop yield analysts may consider combining a deterministic approach with probabilistic analyses.

Wheat growing in the breadbasket region is overwhelmingly composed of highly-mechanized wheat monocultures heavily relying on the use of fertilizers and pesticides⁵⁵. Despite the steady use of fungicides⁵⁶ over the last decade, and intensified use of chemical inputs in 2016, a widespread disease occurrence severely impacted wheat yields. There are ecological arguments
suggesting that monocultures are less resilient to abnormal climate events⁵⁷ and more sensitive to disease outbreaks than more complex cropping systems or landscapes⁵⁸⁻⁶⁰. On the other hand, is the call for a loosening of government restrictions on the use of herbicides and pesticides in order to deal with yield fluctuations⁶¹. These opposing recommendations will not only shape future wheat production in France but also its transformative path toward climate change adaptation.

## Methods

**Yield and climate data.** We analysed winter-wheat yield time series in France for the period 1958–2016 at the spatial scale of departments (administrative units known in French as départements) based on official survey data⁸. In each department we applied a detrending method to crop-yield data to remove the long-term effect of technological improvements within the study period. Relative yield anomalies $\bar{a}_{i,t}$ are defined as:

$$ \bar{a}_{i,t} = \frac{(Y_{i,t} - \mu_{i,t})}{\mu_{i,t}}, \qquad (1) $$

where $Y_{i,t}$ is the yield value and $\mu_{i,t}$ the expected yield value in the $i$th unit at year $t$. Expected yield values, $\mu_{i,t}$, are estimated using local regressions (loess). We define severe and extreme yield loss from relative anomalies below −10 and −15% of expected yields. We also define net yield loss as having all negative relative anomalies.

Climate conditions are spatially uniform at the scale of these administrative units (typically 30–100 km across), thus ensuring coherence between climate data sets and the impacted wheat yield. The studied area is also relatively homogenous in terms of wheat production systems⁵⁵ topography or weather conditions. No crop mask was used because initial tests showed no difference in climate between the entire territory of administrative units and their cropland-covered fraction (not shown). The winter-wheat growing season starts with sowing in October, undergoes a vernalizing period in winter, and ends at harvest in the following July. We refer to the growing season using harvest years (e.g., 2016 encompasses October 2015 to July 2016). Our data set covers October 1958 to July 2016 for harvests occurring from 1959 to 2016. Input climate data are from the SAFRAN reanalysis¹¹,⁶² updated by the French weather service from October 1958 to July 2016. The data cover France on an 8 × 8 km grid on a daily time step. We computed monthly values in each grid cell for the following variables: average maximum (Tmax) and minimum temperatures (Tmin) (°C), average precipitation (mm d⁻¹), average solar radiation (W m⁻²), average Penman-Monteith potential evapotranspiration (mm d⁻¹), number of days with Tmax between 0 and 10 °C (i.e., vernalizing days¹⁸) average number of rainy days with precipitation per month, from October to July (the wheat growing season). Monthly data were also averaged during October–November–December (OND) and April–May–June–July (AMJJ) henceforth called autumn and spring. Climate data were then aggregated over the territory of each administrative unit.

We calculated the frequency of occurrence of climate conditions more extreme than those of the 2015–2016 growing season as follows. Let $X_{i,t}$ be the value of a climate variable (monthly or seasonal mean) during the $i$th growing season (1959–2015) in the $i$th administrative unit. $X_{i,2016}$ is the value of the same variable during the 2015–2016 growing season

$$ I_{i,t} = 0 \quad \text{if } X_{i,t} \le X_{i,2016} $$

$$ I_{i,t} = 1 \quad \text{if } X_{i,t} > X_{i,2016} $$

The frequency of a value X strictly superior to the one of the 2015–2016 growing season is given by

$$ \text{Freq}_i = \frac{\sum_{t=1}^{N} I_{i,t}}{N}, \qquad (2) $$

where N is the number of growing seasons (N = 58). We then averaged Freq_i over all administrative units (n = 27) within the study area into Freq_i and identified extreme regional variables during the 2015–2016 growing season as those with $\bar{Freq}_i<0.05$ or $\bar{Freq}_i>0.95$ (i.e., the average of occurrence frequencies of X across all the administrative units is lower than 0.05 or higher than 0.95).

The number of years with maximum (alt. minimum) temperatures in December exceeding the value of December 2015 is null in all of the 27 administrative units. November was also extremely warm (Fig. 2). Temperatures during the autumn of 2015 have a frequency of occurrence of 0–0.05. For May precipitation, May potential evapotranspiration and radiation in June, the frequency of years exceeding the value of 2016 is also below 5%, Similarly, no more than 5% of the years have a minimum temperature over the growing season higher than that of 2016.

Considering the conjunction of October–November–December temperatures and April–May–June–July (spring) precipitation averaged over the study area, the

year 2016 is a single outlier, in conditions opposed to the ones of the drought year 1976 (the second most important yield drop on the time period considered, Supplementary Fig. 5). On average over the study area, the closest years to 2016 for autumn temperatures and spring rainfall are 2007 followed by 2012, 2001, and 1995. The years 1987, 2013, and 2001 were those with the highest number of administrative units close to 2016 for of spring precipitation. The year 1995 was the closest to 2016 for the total (low) number of vernalizing days. Years 2012, 2007, and 1995 were the closest to 2016 for autumn maximum temperatures.

**Statistical analyses.** We used a suite of binomial logit regression models trained using climate and yield data to estimate the probability of net, severe and extreme yield loss from climate inputs. Both monthly and seasonal averages of climate data were used to construct different models. Among all possible models including a single input variable, we selected the most parsimonious ones, according to their Bayesian Information Criteria (BIC), to identify the most influential predictor climate variable. We then used a stepwise selection procedure to identify the best combinations of input variables, with and without interactions. The model with the lowest BIC was finally selected for each level of yield loss independently (i.e., three best models, one per level of yield loss, Supplementary Table 2). During the model selection process, each model was fitted by maximum likelihood to binary data indicating occurrence of net, severe and extreme yield losses (see SI). Computations were done with the functions glm (family = binomial), predict.glm, and step.glm of R (Version 3.1.0).

The variables found to be influential for predicting yields are consistent over all three best models (Supplementary Table 2). These variables are:

* The number of days with Tmax between 0 and 10 °C in December

* November precipitation or average number of rainy days in November

* Minimum June temperature

* AMJJ or June precipitation.

Two interactions are also selected, namely between the number of days between 0 and 10 °C and AMJJ precipitation, and between minimum June temperatures and AMJJ precipitation. The estimated parameter values of the selected models are presented in Supplementary Table 2. Probabilities of yield loss in 1959–2016 are computed in each of the 27 units of the breadbasket using each selected model. To test the robustness of our results to the definition of the study area, we computed the probability of loss over a larger number of administrative units (including those outside of the breadbasket, i.e., 35 and 45 administrative units, see Supplementary Fig. 7) with the model trained on the study area. To predict the occurrence of yield loss in 2016, the selected models were fitted to the time series excluding 2016, i.e., from 1959 to 2015, and the predictive probabilities of yield loss in 2016 were computed as above. Results for severe and extreme yield losses are presented in Fig. 3 and Supplementary Fig. 4 for net losses.

Reported probabilities must be interpreted as a departure from a prior probability with this probability set to the frequency of yield losses in the samples. The odds obtained with the prior probabilities and with the statistical model are defined by:

$$ O_{\text{prior}} = \frac{P_{\text{prior}}}{1 - P_{\text{prior}}} \quad \text{and} \quad O_{\text{stat}} = \frac{P_{\text{stat}}}{1 - P_{\text{stat}}} \qquad (3) $$

with $P_{\text{prior}}$ equal to 0.47 (alt. 0.46, 2016 excluded) for net yield loss; 0.16 (alt. 0.14, 2016 excluded) for severe yield loss, and 0.083 (alt., 0.077 2016 excluded) for extreme yield loss. The probabilities $P_{\text{stat}}$ are the ones computed from the statistical model (with or without the 2016 data). The odds ratio thus corresponds to the ratio of $O_{\text{stat}}$ to $O_{\text{prior}}$. We also refer to the risk ratio $r$, defined by

$$ r = \frac{P_{\text{stat}}}{P_{\text{prior}}} \qquad (4) $$

**Climate projections.** Climate projections from the Coupled Model Intercomparison Project phase 5 (CMIP5)⁶³ were used for the 1951–2100 time period. Before 2005, we used the so-called historical simulations, in which the climate models are forced by the historical evolution of the main natural and anthropogenic forcings. After 2005, the results of two Representative Concentration Pathways (RCP) scenarios⁶⁴ are contrasted.

The RCP2.6 scenario (RCP26) assumes aggressive mitigations policies to likely limit global warming since the pre-industrial period to 2 °C⁶⁵. This scenario is therefore close to the objectives of the Paris Agreement on Climate. The RCP8.5 scenario (RCP85) assumes a “business as usual” approach to the climate change issue, and results in a global warming close to 4 °C at the end of the 21st century⁶⁶.

Climate projections from a subset of 13 models (bcc-csm1-1-m, BNU-ESM, CanESM2, CCSM4, CNRM-CM5, CSIRO-Mk3-6-0, GFDL-CM3, HadGEM2-ES, IPSL-CM5A-MR, MIROC5, MPI-ESM-MR, MRI-CGCM3, and NorESM1-M) were analyzed. Only the models with both the RCP26 and RCP85 scenarios, and all the variables necessary for our study have been selected. We also only used one climate model by modeling center, to limit the lack of independence within the ensemble, given the strong similarity that generally exists between same center
models⁶⁷. We then assume the independence of the selected models, as most studies to date and the IPCC report⁶⁸.

Different climate variables identified as influential for yield loss are extracted from raw model outputs on their native grids. Then, the spatial averages are computed on a domain that encompasses the 27 units of interest and whose boundaries are: 45.5°N, 51.5°N, -1.5°E, 8°E. Only the grid points with a fraction of land greater than 75% are used. The same domain is used to compute the spatial averages for SAFRAN climate indices.

Anomalies relative the 1959–1988 reference period are analyzed to deal with the potential mean climatological biases in the CMIP5 projections. We then make the implicit hypothesis that the model distributions for the inter-annual anomalies are realistic. This assumption is reasonable given the purpose of our analyses and is not crucial for our conclusions. For instance, the relative cumulative frequencies of the 2016 anomalies on the 1959–2016 period are close in the models and SAFRAN (96.5%, 85%, 40%, 3% in the models versus 93%, 90%, 43%, 2% in SAFRAN for, respectively, spring precipitation, June minimum temperature, November precipitation and the number of vernalizing days in December).

**Data availability.** The yield data that support the findings of this study are available from the corresponding author upon reasonable request. SAFRAN data are accessible here: http://mistrals.sedoo.fr/HyMeX/Data-Access-Registration/? project_name=HyMeX. CMIP data are accessible here: https://esgf-node.llnl.gov/projects/esgf-llnl/

Received: 23 July 2017 Accepted: 4 April 2018

Published online: 24 April 2018

## References

1. MacDonald, R. B. & Hall, F. G. Global crop forecasting. *J. Chem. Inf. Model.* **208**, 670–679 (1980).

2. Chipanshi, A. et al. Evaluation of the integrated Canadian crop yield forecaster (ICCYF) model for in-season prediction of crop yield across the Canadian agricultural landscape. *Agric. For. Meteorol.* **206**, 137–150 (2015).

3. Bussay, A., van der Velde, M., Fumagalli, D. & Seguini, L. Improving operational maize yield forecasting in Hungary. *Agric. Syst.* **141**, 94–106 (2015).

4. IPCC. Managing the risks of extreme events and disasters to advance climate change adaptation. In *A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change* (ed. Field, C. B., et al.) 582 (Cambridge University Press, Cambridge and New York, NY, 2012).

5. Lesk, C., Rowhani, P. & Ramankutty, N. Influence of extreme weather disasters on global crop production. *Nature* **529**, 84–87 (2016).

6. Vautard, R. et al. The European climate under a 2°C global warming. *Environ. Res. Lett.* **9**, 34006 (2014).

7. FAOSTAT. Statistical database, http://faostat.fao.org (2013).

8. Ministère de l'agriculture de l'agroalimentaire et de la foret. Agreste. https://stats.agriculture.gouv.fr/disar/. Accessed 1 Jan 2017.

9. AJG Simoes, C. H. The economic complexity observatory. Workshops at the twenty-fifth AAAI conference on artificial intelligence. https://atlas.media.mit.edu/en/. Accessed 1 Nov 2017 (2011).

10. Agri4cast, J. Crop monitoring in Europe. http://ies-webarchive-ext.jrc.org/mars/mars/Bulletins-Publications.html. Accessed 1 August 2017 (2016).

11. Quintana-Segui, P. et al. Analysis of near-surface atmospheric variables: validation of the SAFRAN analysis over France. *J. Appl. Meteorol. Climatol.* **47**, 92–107 (2008).

12. Trnka, M. et al. Adverse weather conditions for European wheat production will become more frequent with climate change. *Nat. Clim. Change* **4**, 637–643 (2014).

13. Eurostat. Statistical Office of the European Communities. http://ec.europa.eu/eurostat/fr/web/gisco/geodata/reference-data. Accessed 1 May 2017 (2016).

14. Ribes, A., Corre, L., Gibelin, A. & Dubuisson, B. Issues in estimating observed change at the local scale—a case study: the recent warming over France. *Int. J. Climatol.* **36**, 3794–3806 (2016).

15. Ribes, A., Azais, J.-M. & Planton, S. A method for regional climate change detection using smooth temporal patterns. *Clim. Dyn.* **35**, 391–406 (2010).

16. McGrath, J. M. & Lobell, D. B. An independent method of deriving the carbon dioxide fertilization effect in dry conditions using historical yield data from wet and dry years. *Glob. Change Biol.* **17**, 2689–2696 (2011).

17. Fitzgerald, G. J. et al. Elevated atmospheric [CO2] can dramatically increase wheat yields in semi-arid environments and buffer against heat waves. *Glob. Change Biol.* **22**, 2269–2284 (2016).

18. Robertson, M. J., Brooking, I. R. & Ritchie, J. T. Temperature response of vernalization in wheat: modelling the effect on the final number of mainstem leaves. *Ann. Bot.* **78**, 371–381 (1996).

19. Wu, X., Liu, H., Li, X., Tian, Y. & Mahecha, M. D. Responses of winter wheat yields to warming-mediated vernalization variations across temperature Europe. *Front. Ecol. Evol.* **5**, 126 (2017).

20. Rosenzweig, C. & Tubiello, F. N. Effects of changes in minimum and maximum temperature on wheat yields in the central US. A simulation study. *Agric. For. Meteorol.* **80**, 215–230 (1996).

21. FranceAgrimer. Données économiques agricoles et alimentaires. Enquête répartition variétale. https://visionet.franceagrimer.fr/. Accessed 1 May 2017.

22. Lamichhane, J. R. et al. Robust cropping systems to tackle pests under climate change. A review. *Agron. Sustain. Dev.* **35**, 443–459 (2015).

23. Chakraborty, S. et al. Potential impact of climate change on plant diseases of economic significance to Australia. *Australas. Plant Pathol.* **27**, 15 (1998).

24. Coakley, S. M. Climate variability in the Pacific Northwest and its effect on stripe rust disease of winter wheat. *Clim. Change* **2**, 33–51 (1979).

25. Luck, J. et al. Climate change and diseases of food crops. *Plant Pathol.* **60**, 113–121 (2011).

26. West, J. S., Townsend, J. A., Stevens, M. & Fitt, B. D. L. Comparative biology of different plant pathogens to estimate effects of climate change on crop diseases in Europe. *Eur. J. Plant Pathol.* **133**, 315–331 (2012).

27. Olgun, M., Metin Kumlay, A., Cemal Adiguzel, M. & Caglar, A. The effect of waterlogging in wheat (T. aestivum L.). *Acta Agric. Scand. Sect. B* **58**, 193–198 (2008).

28. Shao, G. C. et al. Photosynthesis and growth of winter wheat in response to waterlogging at different growth stages. *Photosynthetica* **51**, 429–437 (2013).

29. Deswarte, J.-C. Récolte des blés en 2016 une année complètement atypique. http://www.franceagrimer.fr/fam/content/download/47241/452162/file/DOSSIER%20DE%20PRESSE%20ARVALIS-Institut%20du%20v%C3%A9g%C3%A9tal%20%BILAN%20DE%20CAMPAGNE%20BLES%2014%20septembre%202016.pdf. Arvalis, Institut du végétal. Conférence de presse du 14 septembre 2016. Accessed 20 March 2018 (2016).

30. Porter, J. R. & Gawith, M. Temperatures and the growth and development of wheat: a review. *Eur. J. Agron.* **10**, 23–36 (1999).

31. Wiegand, C. L. & Cuellar, J. A. Duration of grain filling and kernel weight of wheat as affected by temperature. *Crop Sci.* **21**, 95–101 (1981).

32. Ben-Ari, T. et al. Identifying indicators for extreme wheat and maize yield losses. *Agric. For. Meteorol.* **220**, 130–140 (2016).

33. Gouache, D., Bouchon, A.-S., Jouanneau, E. & Le Bris, X. Agrometeorological analysis and prediction of wheat yield at the departmental level in France. *Agric. For. Meteorol.* **209**, 1–10 (2015).

34. Malik, A. I., Colmer, T. D., Lambers, H., Setter, T. L. & Schortemeyer, M. Short-term waterlogging has long-term effects on the growth and physiology of wheat. *New Phytol.* **153**, 225–236 (2002).

35. Bailey, K. L. & Lazarovits, G. Suppressing soil-borne diseases with residue management and organic amendments. *Soil Tillage Res.* **72**, 169–180 (2003).

36. Pirttioja, N. et al. A crop model ensemble analysis of temperature and precipitation effects on wheat yield across a European transect using impact response surfaces. FACCE MACSUR Reports, [S.l.], v. 6, p. D-C4.4.3. http://ojs.macsur.eu/index.php/Reports/article/view/D-C4.4.3. Accessed 27 Mar 2018 (2015).

37. Ceglar, A., Toreti, A., Lecerf, R., Van der Velde, M. & Dentener, F. Impact of meteorological drivers on regional inter-annual crop yield variability in France. *Agric. For. Meteorol.* **216**, 58–67 (2016).

38. Daryanto, S., Wang, L. & Jacinthe, P.-A. Global synthesis of drought effects on maize and wheat production. *PloS ONE* **11**, e0156362 (2016).

39. Stratonovitch, P. & Semenov, M. A. Heat tolerance around flowering in wheat identified as a key trait for increased yield potential in Europe under climate change. *J. Exp. Bot.* **66**, 3599–3609 (2015).

40. Semenov, M. A. & Shewry, P. R. Modelling predicts that heat stress, not drought, will increase vulnerability of wheat in Europe. *Sci. Rep.* **1**, 66 (2011).

41. Talukder, A. S. M. H. M., McDonald, G. K. & Gill, G. S. Effect of short-term heat stress prior to flowering and early grain set on the grain yield of wheat. *Field Crops Res.* **160**, 54–63 (2014).

42. Asseng, S., Foster, I. A. N. & Turner, N. C. The impact of temperature variability on wheat yields. *Glob. Change Biol.* **17**, 997–1012 (2011).

43. Lobell, D. B. & Asseng, S. Comparing estimates of climate change impacts from process-based and statistical crop models. *Environ. Res. Lett.* **12**, 1–12 (2017).

44. Asseng, S. et al. Uncertainty in simulating wheat yields under climate change. *Nat. Clim. Change* **3**, 827–832 (2013).

45. Moeller, C., Asseng, S., Berger, J. & Milroy, S. P. Plant available soil water at sowing in Mediterranean environments—Is it a useful criterion to aid nitrogen fertiliser and sowing decisions? *Field Crops Res.* **114**, 127–136 (2009).

46. Wong, M. T. F. & Asseng, S. Yield and environmental benefits of ameliorating subsoil constraints under variable rainfall in a Mediterranean environment. *Plant Soil* **297**, 29–42 (2007).

47. Rosenzweig, C., Tubiello, F. N., Goldberg, R., Mills, E. & Bloomfield, J. Increased crop damage in the US from excess precipitation under climate change. *Glob. Environ. Change* **12**, 197–202 (2002).
48. Palosuo, T. et al. Simulation of winter wheat yield and its variability in different climates of Europe: a comparison of eight crop growth models. *Eur. J. Agron.* **35**, 103–114 (2011).

49. Martre, P. et al. Multimodel ensembles of wheat growth: many models are better than one. *Glob. Change Biol.* **21**, 911–925 (2015).

50. BSV. Bulletin de santé du végétal. http://agriculture.gouv.fr/bulletins-de-sante-du-vegetal (2016).

51. EauFrance. Bulletin national de situation hydrologique. http://www.eaufrance.fr/documents/documents-213/bulletin-national-de-situation/

52. Semences de France. Moisson 2016: les premières tendances. http://www.semencesdefrance.com/actualite-semences-de-france/moisson-2016-rendement-ble-tendre-premières-tendances/. Accessed 5 May 2017 (2016).

53. Ciais, P. et al. Europe-wide reduction in primary productivity caused by the heat and drought in 2003. *Nature* **437**, 529 (2005).

54. Le Tertre, A. et al. Impact of the 2003 heatwave on all-cause mortality in 9 French cities. *Epidemiology* **17**, 75–79 (2006).

55. Lechenet, M., Makowski, D., Py, G. & Munier-Jolain, N. Profiling farming management strategies with contrasting pesticide use in France. *Agric. Syst.* **149**, 40–53 (2016).

56. Hossard, L., Guichard, L., Pelosi, C. & Makowski, D. Lack of evidence for a decrease in synthetic pesticide use on the main arable crops in France. *Sci. Total Environ.* **575**, 152–161 (2017).

57. Altieri, M. A., Nicholls, C. I., Henao, A. & Lana, M. A. Agroecology and the design of climate change-resilient farming systems. *Agron. Sustain. Dev.* **35**, 869–890 (2015).

58. Costanzo, A. & Bärberi, P. Functional agrobiodiversity and agroecosystem services in sustainable wheat production. A review. *Agron. Sustain. Dev.* **34**, 327–348 (2014).

59. Rottstock, T., Joshi, J., Kummer, V. & Fischer, M. Higher plant diversity promotes higher diversity of fungal pathogens, while it decreases pathogen infection per plant. *Ecology* **95**, 1907–1917 (2014).

60. Claflin, S. B., Jones, L. E., Thaler, J. S. & Power, A. G. Crop-dominated landscapes have higher vector-borne plant virus prevalence. *J. Appl. Ecol.* **54**, 1190–1198 (2017).

61. Académie d'agriculture de France. La production de céréales à paille en France en 2016-qualité sanitaire et technologique des grains, quelques éléments d'explication à partir de l'exemple du blé tendre. Note Académique. https://www.academie-agriculture.fr/publications/academie-communique/lemensuel/ndeg15-oct-2016 (2016).

62. Vidal, J., Martin, E., Franchistéguy, L., Baillon, M. & Soubeyroux, J. A 50-year high-resolution atmospheric reanalysis over France with the Safran system. *Int. J. Climatol.* **30**, 1627–1644 (2010).

63. Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An overview of CMIP5 and the experiment design. *Bull. Am. Meteorol. Soc.* **93**, 485–498 (2012).

64. Van Vuuren, D. P. et al. The representative concentration pathways: an overview. *Clim. Change* **109**, 5 (2011).

65. Schleussner, C.-F. et al. Differential climate impacts for policy-relevant limits to global warming: the case of 1.5 °C and 2 °C. *Earth Syst. Dyn. Discuss.* **6**, 2 (2015).

66. Collins, M. et al. In: *Climate Change 2013: the physical science basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change* (eds. Stocker, T. F. et al.) (Cambridge University Press, 2013).

67. Leduc, M., Laprise, R., de Elia, R. & Šeparović, L. Is institutional democracy a good proxy for model independence? *J. Clim.* **29**, 8301–8316 (2016).

68. Kass, R. E. & Raftery, A. E. Bayes factors. *J. Am. Stat. Assoc.* **90**, 773–795 (1995).

## Acknowledgements

This work was supported by the CLAND convergence institute funded by the French National Research Agency (ANR). T.B.-A. and D.M. acknowledge the INRA-CIRAD meta-program GloFoods (PrevSaison). P.C. acknowledges the support of the European Research Council Synergy grant ERC-2013-SyG-610028 IMBALANCE-P. The authors thank Météo France (Direction de la Climatologie et des Services Climatiques) for providing the SAFRAN data set. We acknowledge the World Climate Research Program's Working Group on Coupled Modeling, which is responsible for CMIP, and we thank the climate modeling groups that developed the models listed in the method section of this paper for producing and making available their model output. For CMIP the U.S. Department of Energy's Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portals. The authors would also like to thank Phil Tajitsu Nash for editing, Emmanuelle Gourdain for helping us retrieve some of the wheat data, and Marie Launay, Antoine Gardarin, Gilles Grandeanu and Dominique Le Floch for insightful discussions.

## Author contributions

T.B.-A., J.B., P.C. and D.M. framed the study. T.B.-A. and J.B. completed the data sets. T.B.-A. and D.M. performed statistical analyses of historical data. J.B. analysed CMIP5 data. All authors interpreted the results. J.B., R.L. & M.V.d.V. contributed to the writing of specific sections of the manuscript; T.B.-A., P.C. and D.M. wrote the paper.

## Additional information

Supplementary Information accompanies this paper at https://doi.org/10.1038/s41467-018-04087-x.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018