# Complementary mechanisms stabilize national food production

Lucie Mahaut, Cyrille Violle, Delphine Renard

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

## To cite this version:

Lucie Mahaut, Cyrille Violle, Delphine Renard. Complementary mechanisms stabilize national food production. Scientific Reports, 2021, 11 (1), pp.4922. 10.1038/s41598-021-84272-z. hal-04906792

HAL Id: hal-04906792

https://hal.science/hal-04906792v1

Submitted on 22 Jan 2025

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_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Complementary mechanisms stabilize national food production

Lucie Mahaut*, Cyrille Violle & Delphine Renard

Ensuring the temporal stability of national food production is crucial for avoiding sharp drops in domestic food availability. The average stability of individual crop yields and asynchrony among crop yield fluctuations are two candidate mechanisms to stabilize national food production. However, the quantification of their respective influence on the stability of national food production is lacking, as is the identification of the factors regulating both mechanisms. Using yield data for 138 crops and 115 countries over a 50-year period, we first show that the stability of total national yield mostly relies on the fluctuations of the yield of crops covering the largest share of cropland. The average yield stability of these crops exert a stabilizing effect on national food production that is twice as important as the one of the asynchronous yield fluctuations among them. Climate variability reduces the stability of national food production by synchronizing yield fluctuations among crops and destabilizing the yield of individual crops. However, our results suggest that increasing crop diversity can counteract the synchronizing effects of climate variability by enhancing asynchronous dynamics among crops. Irrigation can promote the average stability of individual crop yields but cannot compensate for the destabilizing effect of climate variability. Considering both the response of each crop to climatic variations and the dynamics emerging from crop baskets will help agricultural policies to ensure stable food supply at the national level.

Maintaining the year-to-year stability of national food supply is a critical target for agricultural policies to support food security worldwide. The increasing frequency of extreme climate events such as drought led to strong fluctuations in food production at the national (e.g.¹) and global²⁻⁴ scales. Such fluctuations in crop production can adversely affect the stability of food supply and impede food security⁵. Increasing domestic grain reserves⁶ and importing food through international trade⁷⁻⁹ can compensate for sharp drops in domestic food availability. However, grain reserves are declining worldwide⁸ and international trade could be source of instability and threaten food sovereignty and security in low-income nations¹⁰. Moreover, with increasing synchronization of major crop failures at the global scale⁴, relying on food imports to ensure food supply seems a limited solution in the long term. In this context, increasing the stability of national food production is a promising, alternative solution to safeguard society against food shortages. Achieving this goal requires the investigation of the different levers of action and of the mechanisms leading to greater stability of national food production.

Decades of research in ecology show that the year-to-year stability of ecosystem productivity results from two distinct mechanisms: the stability of each individual species and the asynchronous variation among species that allows for losses and gains in productivity to compensate each other (i.e. asynchronous dynamics¹¹). In cropping systems, the use of agricultural inputs such as irrigation water and fertilizers is one way to buffer the effect of environmental fluctuations and stabilize individual crop yields (e.g.³,¹²). However, ecological experiments in grasslands show that nutrient enrichment reduces the stability of ecosystem productivity by synchronizing the fluctuations among plant species (e.g.¹³⁻¹⁵). Investigating whether and how agricultural inputs affect asynchronous dynamics among crops is therefore a critical issue to understand the overall effects of agricultural inputs on the stability of national food production.

Another promising strategy to stabilize total national yield is to increase the diversity of crop species or crop groups grown within a country¹⁶. However, the mechanism underlying the crop diversity-national yield stability relationship remains to be elucidated. Theories and experiments in biodiversity-ecosystem functioning research (BEF) show that increasing plant diversity enhances asynchronous dynamics (reviewed in¹⁷,¹⁸). While the role of species diversity in asynchrony has been extensively studied at the local scale, where species interact (e.g.¹⁹), recent theoretical developments predict that these dynamics can also play out across larger spatial scales²⁰,²¹. Importantly, attaining greater yield stability through higher crop diversity will also require identifying crops with asynchronous dynamics. Renard and Tilman¹⁶ provide evidences that the stabilizing effect of crop diversity on national food production depends on the number of crops representing an equal share of harvested area
('effective' diversity) and not simply on the total number of crops grown in a country. This suggests that fluctua-
tions in the yield of crops covering a large share of national cropland area have a greater impact on the stability
of total national yield than fluctuations in the yield of crops occupying only small areas²².

In this study, we quantify the respective role of the average stability of individual crop yields and asynchrony
among crop yield fluctuations on national yield stability. Following Renard and Tilman¹⁶, we compute total
national yield stability as the ratio between the mean national yield (i.e. the summed annual kilocalories har-
vested across all edible crops in a nation, divided by the total area [ha] harvested) and the year-to-year standard
deviation of this ratio. Although calories are not the most limiting nutrients in terms of nutrition at the global
scale²³, using caloric yields allow us to compare the production of the different crops grown within each nation.
We further test whether fluctuations of the yield of the most abundant crops have a greater impact on the stability
of total national yield. To do so, we calculate the average stability of yields of individual crops and asynchrony
among crop yield fluctuations using two measures of crop yield, one non-weighted, and the other for which
the yield is weighted by the proportion of total cropland it occupies (see Methods). Similarly, we quantify crop
diversity by both crop richness and the exponential of the Shannon diversity index (H') that is a measure of the
effective diversity²⁴. H' weights each crop in a nation by the proportion of total cropland it occupies so that the
crops that are most produced in a country count more than the minor crops in the calculation of diversity (see
Methods). Finally, we assess by which of these mechanisms climate variability, crop diversity and agricultural
inputs affect the stability of national food production.

Results and discussion

We used structural equation models (SEMs) to disentangle the role of two stabilizing mechanisms (i.e. the aver-
age stability of individual crop yields and asynchrony) and three factors (i.e. climate variability, crop diversity
and agricultural inputs) on national yield stability. We built two separate models, one based on crop richness,
individual crop yield stability and asynchrony and the other based on these same indices weighted by the crop's
share of total area harvested.

Our results reveal that the stability of national food production mainly depends on yield fluctuations of crops that cover a large proportion of the national harvested area. The model including the area-weighted measure of national crop yield explains 87% of the variance of national yield stability versus 48% for the model that does not account for differences in area planted for different crops (Fig. 1). This echoes recent findings showing that the fluctuations of the most largely produced crops have a larger weight on the stability of the total national food production<sup>22</sup>. Although crops with a lower share of cropland weight less in the process, their cultivation deliver other critical services for food security, for example by providing higher economic incomes or nutritional intakes<sup>23,25</sup>.

Both the average stability of individual crop yield and yield asynchrony have significant, stabilizing effects on total national yield when yields are weighted by the proportion of area harvested (Figs. 1, 2). Importantly, this finding indicates that both fluctuations in the yield of a single crop and the dynamics that occur between crop yield fluctuations must be considered to ensure stable food production. However, the influence of the average stability of individual crop yields was twice as important as the one of asynchrony (Fig. 1). Conversely, when crop's abundance is not accounted for, asynchrony has no significant effect on the stability of national yield (Fig. 1), suggesting that yield gains in crops representing low share of area harvested cannot compensate for yield losses in the most abundant crops at the national scale.

Climate variability has been identified as a main determinant of national food production instability¹⁶. Our results show that this destabilizing effect arises not only from a negative effect of precipitation and temperature variability on the average stability of individual crop yields, as previously observed¹, but also from a synchronization of individual crop yield fluctuations (Fig. 1). This brings novel evidence that climate is a main driver of the synchronization of the global production of major commodities⁴.

Previous studies have revealed the potential of crop diversity to increase the stability of crop production at
the national scale¹⁶,²². Accordingly, we find a strong, positive effect of crop diversity on asynchrony (Fig. 3) that
can largely counteract the synchronizing effect of climate variability (Fig. 1). This result indicates that asynchrony
among the most abundant crops is the main process by which crop diversity stabilizes national food production¹⁶,
confirming that the mechanisms invoked in BEF studies, classically conducted at local scales, can be extended
to larger spatial scales²⁰,²¹. However, whereas local-scale BEF studies often involve direct biotic interactions in
driving asynchronous dynamics among species (e.g., ref.¹⁹), our dataset at national scale does not allow consid-
erations of the effects of crop-crop interactions. Further research will be needed to go beyond our correlative
approach to better understand the links between crop diversity, asynchrony and yield stability at such a large
scale. In particularly, the role of spatial heterogeneity in environmental conditions in shaping the distribution
of crops and in generating spatial asynchronous dynamics among different crops merits deeper exploration²⁰,²¹.

Investigating the role of agricultural inputs shows that greater use of irrigation stabilizes total national food production by increasing the average stability of crops that occupy the largest proportion of the national harvested area (Figs. 1, 4). This probably reflects the fact that the most important crops for food consumption and trade (e.g. rice, cereals, maize) are also the most irrigated²⁶. However, the stabilizing effect of irrigation on the yields of these crops is by far lower than the destabilizing one of climate variability (Fig. 1). The selection of varieties that are less reliant on irrigation and more resistant to climate variability will consequently be important to ensure a higher stability of crop yields²⁷⁻²⁹. Finally, we find that nitrogen fertilisation weakly affects the stability of the total, national yield although it is negatively associated with the average stability of individual crop yields when crop's abundance is not accounted for (Fig. 1).

Ensuring the stability of food production at national level is becoming an increasingly important challenge
for agricultural policies. Overall, our study reveals that both the average stability of individual crop yields and
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Structural equation model representing the drivers of the stability of national food production. For a seek of clarity, the contribution of agricultural inputs (red path) and crop diversity (blue path) to the stability of national food production are represented separately from the ones of climate variability (green path) [left and right panels, respectively]. Top: Non-weighted model. Crop diversity is quantified by crop richness. Down: Abundance weighted-model. Crop diversity is quantified by the exponential of Shannon diversity index. Asynchrony and average stability are weighted by the proportion of total cropland each crop occupied. The thickness of the arrows indicates the relative contribution of each variable. Plain arrows represent positive relationships while dotted arrows represent negative relationships between two boxes. Only relationships supported by the data (*P* > 0.05) are shown. Standardized regressions weights (along arrows) and squared coefficient of regressions (r²) for the fitting model are shown. Test indicates close model-data fit (Fisher's C = 0.16, *P* = 0.92 and Fisher's C = 0.46, *P* = 0.79 for non-weighted and weighed models, respectively). National yield stability, nitrogen use intensity and the percentage of land equipped for irrigation are log-transformed.</figcaption></figure>

asynchrony in yield fluctuations among crops are important mechanisms to stabilize national food production in the face of climate variability. Promoting crop diversity at the national level might be a solution to promote multiple benefits, including greater stability of food production¹⁶,²², a diversified diet³⁰ and the reduction of the use of agricultural inputs³¹. However, not all baskets of crops can provide all of these services, and our results suggest that it is important to design crop diversification strategies in a way that promotes asynchronous yield fluctuations between crops, for example by selecting species with different responses to environmental fluctuations. Finally, while the national level is the one at which agricultural policies are made, working at such a large scale does not allow to fully capturing determinants of food supply at finer scales. Transposing our approach to smaller scales will provide a better understanding of the determinants of stable food production worldwide.

## Materials and methods

**National yield stability.** We used the FAOSTAT database (http://www.fao.org/faostat, visited in September 2019) to obtain data on annual crop production (in tons) and area harvested (in hectares) from 1961 to 2010 for 138 crops in 91 populous nations. Following Renard and Tilman¹⁶, we accounted for differences among nations in data quality and excluded five nations, namely North Korea, Guinea, Kenya, Mozambique and Zambia, for which at least 20% of the data on area harvested or production were extrapolated by the FAO (see details in¹⁶). We calculated for each nation and each year the total annual caloric yield (millions of kcal ha⁻¹). To do so, we first calculated the kcal production of each crop by multiplying the production of each crop by its commodity-specific kilocalorie conversion factor from the USDA Nutrient Database³². In doing so, we were able to compare the production of different crops. Then, we summed these kcal harvests across all crops and divided this value by the sum of harvested area for all crops. We calculated national yield stability (S) as the ratio of mean total annual caloric yield (μT) over its time-detrended standard deviation (σT) for fifty consecutive years (1961–2010). We accounted for a temporal trend of increasing total annual crop yield by implementing a loess regression between annual crop yield and years. σT corresponds to the standard deviation of the residuals of this regression. Finally, we compared this stability index (largely used in the biodiversity-ecological functioning research, e.g.¹⁴,¹⁶⁻¹⁸) with the resilience index used by Zampieri et al.²². Both indices were strongly correlated (r = 0.992), strengthening our findings.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Effects of the average stability of individual crop yield and asynchrony on total, national yield stability. Black lines show the relationships estimated by the structural equation model when the crop yields are non-weighted (a, b) and weighted (c, d) by the proportion of cropland occupied by each crop within a country. Dashed lines represents non significant relationship.</figcaption></figure>

**Individual crop yield stability and yield asynchrony.** For each country, we quantified the average stability of yields of individual crops as the mean of the inverse of the coefficient of variation of yield of each crop:

$$ \left( \sum_{i=1}^{N} \frac{\mu_i}{\sigma_i} \right) / N \qquad (1) $$

where $\mu_i$ is the temporal mean of crop's annual kcal yield and $\sigma_i$ its time-detrended standard deviation. Time-detrended crop yield was computed through a loess regression between individual, annual crop yield and years.

We computed the asynchrony between crop yield fluctuations following the index developed by Loreau and De Mazancourt¹¹:

$$ \Phi = 1 - \frac{\sigma_T^2}{\left(\sum_{i=1}^{N} \sigma_i\right)^2} \qquad (2) $$

where $\Phi$ is the asynchrony of crop species based on annual caloric yield (millions of kcal ha⁻¹) with $\sigma_T^2$ the temporal variance of the time-detrended national yield and $\sigma_i$ the time-detrended standard deviation of each crop's annual kcal yield. The value of asynchrony varies between zero (perfect synchrony) and one (perfect asynchronous temporal fluctuations).
<figure><img src="image_6.png" /><figcaption>Figure 3. Effects of crop diversity on asynchronous dynamics among crop yield fluctuations. Black lines show the relationships estimated by the structural equation model when the crop yields are non-weighted (a) and weighted (b) by the proportion of cropland occupied by each crop within a country.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Figure 4. Effects of the percentage of land equipped for irrigation on the average stability of individual crop yields. Black lines show the relationships estimated by the structural equation model when the crop yields are non-weighted (a) and weighted (b) by the proportion of cropland occupied by each crop within a country.</figcaption></figure>

To test whether yield fluctuations of the most abundant crops have a greater impact on the stability of national food production, we weighted the annual yield of each crop by the proportion of total harvested area occupied by that crop. Average stability of yields of individual crops and yield asynchrony were computed on both the non-weighted and abundance-weighted yields.

**Crop diversity.** For each country and year, we used both the total number of crop commodities (i.e. crop richness) and the Shannon information index (H') to quantify crop diversity. H' weights each crop in a nation by the proportion of total cropland it occupies ($p_i$):

$$H' = - \sum_{i=1}^{N} (p_i \ln[p_i]) \qquad (3)$$

with N being the total number of crops grown in a country each year.

The exponential form of the Shannon diversity index gives the effective crop diversity that is the number of crops representing an equal share of harvested area²⁴. In other words, the exponential of the Shannon diversity index weighs all species by their frequency, without favouring either common or rare species²⁴. We averaged
the annual effective diversity of crop across the fifty years studied to test the effect of crop diversity on national yield stability.

**Agricultural inputs.** We extracted the annual national application of nitrogen and the annual cropland area equipped for irrigation from the FAOSTAT database. Because Ireland, New Zealand and Netherlands use much of their fertilizers on pastures rather than croplands, we excluded these nations from our analysis. Similarly, we excluded Egypt because it has 100% of cropland equipped for irrigation. We calculated the annual rates of nitrogen application and irrigation per hectare by dividing their use by the total annual cropland area.

**Climate variability.** We used global gridded climatic data from the Climate Research Unit of the University of East Anglia³³ to compute the year-to-year variability of growing season precipitation and temperature for each country, both strongly affecting the stability of national food production¹⁶. From these data, we derived annual precipitation and temperature for each grid cell in a country by taking the sum of monthly precipitation and the mean of monthly temperature values weighted by the proportion of cropland in each grid cell³⁴. We then computed the year-to-year coefficient of variation of cropland-based temperature and precipitation for each country.

**Statistical analysis.** We used structural equation models (SEMs) to evaluate how irrigation, intensity of use of nitrogen fertilizers and crop diversity affected national yield stability through changes in the average stability of yields of individual crops and asynchrony of yields. SEMs represent a powerful way to disentangle complex mechanisms controlling crop diversity-stability relationships, as previously done in natural ecosystems (e.g.¹⁴,¹⁵,³⁵,³⁶). We set up two different structural equation models, one based on non-weighted indices of stability of individual crops and asynchrony, the other based on the same indices weighted by the proportion of total harvested area accounted for by each crop. We firstly considered the effects of agricultural inputs and crop diversity on the stability of national food production via the path of average yield stability. The second path quantified the indirect effects of agricultural inputs and crop diversity on national stability via their impacts on crop yield asynchrony. We also accounted for the direct effects of agricultural inputs and crop diversity on national yield stability. Finally, we controlled for the effects of climate variability on total, national yield stability, individual crop yield stability and yield asynchrony. SEMs were run with the lavaan R library³⁷. We used the standardized estimates to compare the relative importance of the different paths. The model fit was evaluated using the Fisher C-score and its associated *p* values. Because the structural equation model assumes linear relationships between predictors and the dependent variable, we also plotted the relationships between total national yield stability and both asynchrony and average stability of individual crop yield to control for linearity (Fig. 2). Similarly, we investigated the relationships between crop diversity and asynchrony (Fig. 3), as well as between irrigation rate and the average stability of individual crop yield (Fig. 4).

**Data availability**

The sources of all data used in this study are referenced in the "Materials and methods" section and all raw data are freely accessible at the URLs provided. The datasets used for the analyses are available from the corresponding author upon request.

Received: 16 October 2020; Accepted: 10 February 2021

Published online: 01 March 2021

**References**

1. Brown, J. K. M. et al. Yield instability of winter oilseed rape modulated by early winter temperature. *Sci. Rep.* **9**, 6953 (2019).

2. Lesk, C. et al. Influence of extreme weather disasters on global crop production. *Nature* **529**, 84 (2016).

3. Ray, D. K. et al. Climate variation explains a third of global crop yield variability. *Nat. Commun.* **6**, 5989 (2015).

4. Mehrabi, Z. & Ramankutty, N. Synchronized failure of global crop production. *Nat. Ecol. Evol.* **3**, 780–786 (2019).

5. Schmidhuber, J. & Tubiello, F. N. Global food security under climate change. *Proc. Natl. Acad. Sci. U. S. A.* **104**, 19703–19708 (2007).

6. Laio, F. et al. The past and future of food stocks. *Environ. Res. Lett.* **11**, 035010 (2016).

7. D'Odorico, P. et al. Feeding humanity through global food trade. *Earths Future* **2**, 458–469 (2014).

8. Marchand, P. et al. Reserves and trade jointly determine exposure to food supply shocks. *Environ. Res. Lett.* **11**, 095009 (2016).

9. Puma, M. J. et al. Assessing the evolving fragility of the global food system. *Environ. Res. Lett.* **10**, 024007 (2015).

10. d'Amour, C. B. & Anderson, W. International trade and the stability of food supplies in the Global South. *Environ. Res. Lett.* https://doi.org/10.1088/1748-9326/ab832f (2020).

11. Loreau, M. & de Mazancourt, C. Species synchrony and its drivers: neutral and nonneutral community dynamics in fluctuating environments. *Am. Nat.* **172**, E48-66 (2008).

12. Herrera, J. M. et al. Lessons from 20 years of studies of wheat genotypes in multiple environments and under contrasting production systems. *Front. Plant Sci.* **10**, 1745 (2020).

13. Liu, J. et al. Nitrogen addition reduced ecosystem stability regardless of its impacts on plant diversity. *J. Ecol.* **107**, 2427–2435 (2019).

14. Zhang, Y. et al. Nitrogen enrichment weakens ecosystem stability through decreased species asynchrony and population stability in a temperate grassland. *Glob. Change Biol.* **22**, 1445–1455 (2016).

15. Song, M.-H. et al. Nutrient-induced shifts of dominant species reduce ecosystem stability via increases in species synchrony and population variability. *Sci. Total Environ.* **692**, 441–449 (2019).

16. Renard, D. & Tilman, D. National food production stabilized by crop diversity. *Nature* https://doi.org/10.1038/s41586-019-1316-y (2019).

17. Loreau, M. & de Mazancourt, C. Biodiversity and ecosystem stability: a synthesis of underlying mechanisms. *Ecol. Lett.* **16**, 106–115 (2013).

18. Thibaut, L. M. & Connolly, S. R. Understanding diversity-stability relationships: towards a unified model of portfolio effects. *Ecol. Lett.* **16**, 140–150 (2013).
19. Tilman, D. et al. Biodiversity and ecosystem stability in a decade-long grassland experiment. *Nature* **441**, 629–632 (2006).

20. Gonzalez, A. et al. Scaling-up biodiversity-ecosystem functioning research. *Ecol. Lett.* **23**, 757–776 (2020).

21. Wang, S. & Loreau, M. Biodiversity and ecosystem stability across scales in metacommunities. *Ecol. Lett.* **19**, 510–518 (2016).

22. Zampieri, M. et al. Estimating resilience of crop production systems: From theory to practice. *Sci. Total Environ.* **735**, 139378 (2020).

23. Kc, K. B. et al. When too much isn't enough: Does current food production meet global nutritional needs? *PLoS One* **13**, e0205683 (2018).

24. Jost, L. Entropy and diversity. *Oikos* **113**, 363–375 (2006).

25. Cheng, A. et al. Diversifying crops for food and nutrition security—a case of teff. *Biol. Rev.* **92**, 188–198 (2017).

26. FAO Aquastat. http://www.fao.org/aquastat/en/geospatial-information/global-maps-irrigated-areas.

27. Harfouche, A. L. et al. Accelerating climate resilient plant breeding by applying next-generation artificial intelligence. *Trends Biotechnol.* **37**, 1217–1235 (2019).

28. Kahiluoto, H. et al. Decline in climate resilience of European wheat. *Proc. Natl. Acad. Sci. U. S. A.* **116**, 123–128 (2019).

29. Zhang, H. et al. Developing naturally stress-resistant crops for a sustainable agriculture. *Nat. Plants* **4**, 989–996 (2018).

30. Herrero, M. et al. Farming and the geography of nutrient production for human use: a transdisciplinary analysis. *Lancet Planet. Health* **1**, e33–e42 (2017).

31. Isbell, F. et al. Benefits of increasing plant diversity in sustainable agroecosystems. *J. Ecol.* **105**, 871–879 (2017).

32. FoodData Central. [Online]. Available: https://fdc.nal.usda.gov/. Accessed: 27-May-2020.

33. Harris, I. et al. Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset. *Sci. Data* **7**, 109 (2020).

34. Ramankutty, N., Evan, A. T., Monfreda, C. & Foley, J. A. Farming the planet: 1. Geographic distribution of global agricultural lands in the year 2000. *Glob. Biogeochem. Cycle* **22**, GB1003 (2008).

35. Bluethgen, N. et al. Land use imperils plant and animal community stability through changes in asynchrony rather than diversity. *Nat. Commun.* **7**, 10697 (2016).

36. Shi, Z. et al. Dual mechanisms regulate ecosystem stability under decade-long warming and hay harvest. *Nat. Commun.* **7**, 11973 (2016).

37. Rosseel, Y. et al. *lavaan: Latent Variable Analysis* (2020).

## Acknowledgements

This work was funded by the French National Research Agency under the Programme "Investissements d'Avenir" under the reference ANR 17 MPG A 0004. CV was supported by the European Research Council (ERC) Starting Grant Project 'ecophysiological and biophysical constraints on domestication in crop plants' (Grant ERC-StG-2014-639706-CONSTRAINTS). We thank Doyle McKey for helpful comments during the preparation of the manuscript and English proofreading.

## Author contributions

L.M., C.V. and D.R. designed the research question; L.M. analysed the data, prepared Fig. 1 and wrote the manuscript. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to L.M.

Reprints and permissions information is available at www.nature.com/reprints.

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

<figure><img src="image_8.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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021