Have precipitation extremes and annual totals been increasing in the world’s dry regions over the last 60 years?
Daily rainfall extremes and annual totals have increased in large parts of the global land area over the last decades. These observations are consistent with theoretical considerations of a warming climate. However, until recently these global tendencies have not been shown to consistently affect land regions with limited moisture availability. A recent study, published by Donat et al. (2016), now identified rapid increases in globally aggregated dry region daily extreme rainfall and annual rainfall totals. Here, we reassess the respective analysis and find that a) statistical artifacts introduced by the choice of the reference period prior to data standardization lead to an overestimation of the reported trends by up to 40%, and also that b) the definition of ‘dry regions of the globe’ affect the reported globally aggregated trends in extreme rainfall. Using the same observational dataset, but accounting for the statistical artifacts and using alternative, well-established dryness definitions, we find no significant increases in heavy precipitation in the world’s dry regions. Adequate data pre-processing approaches and accounting for uncertainties regarding the definition of dryness are crucial to the quantification of spatially aggregated trends in the world’s dry regions. In view of the high relevance of the question to many potentially affected stakeholders, we call for a cautionary consideration of specific data processing methods, including issues related to the definition of dry areas, to guarantee robustness of communicated climate change relevant findings.
2]Jakob Zscheischler 1]Martin Heimann 3]Holger Lange 1,4,5]Miguel D. Mahecha 6]Geert Jan van Oldenborgh 7]Friederike E. L. Otto 1,4,5]Markus Reichstein
1]Max Planck Institute for Biogeochemistry, Jena, Germany 2]Institute for Atmospheric and Climate Science, ETH Zürich, Zürich, Switzerland 3]Norwegian Institute of Bioeconomy Research, Ås, Norway 4]German Centre for Integrative Biodiversity Research (iDiv), Leipzig, Germany 5]Michael Stifel Center Jena for Data-Driven and Simulation Science, Jena, Germany 6]Weather and Climate Modeling, Koninklijk Nederlands Meteorologisch Instituut, De Bilt, Netherlands 7]Environmental Change Institute, University of Oxford, South Parks Road, Oxford, United Kingdom
Sebastian Sippel (firstname.lastname@example.org)
Daily rainfall extremes are expected to increase roughly by 6-7% per C of warming following the Clausius-Clapeyron equation (Allen and Ingram, 2002), if there is enough moisture available. Quantifying and predicting changes in precipitation characteristics due to climate change is crucial for water availability assessments and adaptation to climate change (IPCC, 2012; Greve et al., 2014). On a global scale, daily precipitation extremes have been observed to intensify (Donat et al., 2013; Westra et al., 2013; O’Gorman, 2015), consistent with global model simulations (Fischer and Knutti, 2015), and coincide with a global-scale increase in observed annual rainfall totals (Donat et al., 2013). However, there is little consensus to date on how precipitation characteristics have changed in the past over dry land areas and how they will change in the future. Donat et al. (2016) investigate whether and to what extent daily precipitation extremes (Rx1d) and annual precipitation totals (PRCPTOT) have increased over the last 60 years using observational data. The authors identify rapid increases in rainfall extremes over dry regions, which strongly outpace the corresponding increases over wet areas, and find a similar pattern for annual rainfall.
The question whether precipitation extremes increase in dry regions is highly relevant in the context of climate change adaptation, as generally dry areas may be less prepared to deal with precipitation extremes (Ingram, 2016). Consequently, the recent report on increasing precipitation extremes in dry areas was highlighted in major Science journals (including Nature (Tollefson, 2016) and Nature Climate Change (Ingram, 2016)) and received major media coverage 111http://www.huffingtonpost.com/entry/global-warming-will-bring-extreme-rain-and-flooding-study-finds_us_56e081c7e4b0860f99d796ab 222https://www.theguardian.com/environment/2016/mar/08/hotter-planet-spells-harder-rains-to-come-study 333https://www.sciencedaily.com/releases/2016/03/160308105625.htm 444http://phys.org/news/2016-03-global-world-driest-areas.html 555http://www.abc.net.au/news/2016-03-08/climate-change-could-bring-more-rain-to-deserts-study/7229236 666http://www.asce.org/magazine/20160412-climate-change-to-cause-more-precipitation-in-dry-regions,-researchers-say/, which indicates the importance of this topic for the scientific community, the public and in the context of climate-related decision-making.
However, scrutinizing the findings by Donat et al. (2016) reveals two major issues of concern: Firstly, the applied statistical approach introduces two systematic biases that leads to a substantial overestimation of the increase in precipitation totals and extremes of up to 40% in dry regions. Wet regions, by contrast, are only affected to a limited degree due to a cancellation of errors in trend estimates. Secondly, the definition of dryness used in Donat et al. (2016) based on precipitation alone does not reflect the common understanding of a dry region and thus induces considerable uncertainty. If we test alternative but well-established definitions of a ‘dry region’ (based on water supply and demand, either implicitly or explicitly, cf. Köppen, 1900; Greve et al., 2014) and reapply the appropriate statistical tools, we find no significant increases in precipitation extremes or annual totals in the world’s dry regions. The increases in wet regions are slightly underestimated by Donat et al. (2016) due to a cancellation of errors. These results are of high relevance in the context of making climate change adaptation decisions in dry regions.
1 On data pre-processing based on a time-invariant reference period
As a first step in the analysis of Donat et al. (2016), the authors normalize the 60-year time series in the gridded HadEX2 dataset (Donat et al., 2013) for each grid point with the sample mean of a 30-year reference period (1951–1980), which is a widespread procedure in climate science (Zhang et al., 2005; Sippel et al., 2015). However, this procedure artificially increases the mean and spread of the spatial distribution in the out-of-base period (1981-2010) in all investigated time series, simply because variability in the sample means inflates the signal in the latter period (Sippel et al., 2015) (Fig. 1). To illustrate this point, consider two hypothetical climate regions of the same size: In region one, the mean of a rainfall quantity increases between two periods (from 100 to 200mm, say), for example due to a few large extremes, whereas it decreases by exactly the same amount in region two (i.e. from 200 to 100mm). Consequently, in both time periods the spatial average and the spread of the two regions would be statistically indistinguishable. However, normalizing by the mean of the first time period would imply that the spatial average across both regions for the second period is 1.25 (the average of 0.5 and 2), i.e. a spurious increase of 25% between both periods. In addition, the normalization induces considerable spread in the spatial distribution in the second period.
An additional statistical bias stems from the choice of the world’s 30% wettest and 30% driest regions based on the climatology of PRCPTOT and Rx1d in the reference period (1951-1980). Because 30 years are fairly short to derive a robust climatology of the tails of the precipitation distribution, the computed changes in wet and dry regions are distorted by the "regression to the mean" phenomenon (Galton, 1886; Barnett et al., 2005). In other words, selecting from the dry (wet) end of the spatial distribution in one subset of the dataset will result in a higher probability for wetter (drier) conditions in the remaining years (Figure in Appendix A, Tables in Appendix B).
Hence, the chosen normalization approach combined with the spatial point selection method gives a bias toward precipitation totals and extremes increasing at a faster rate in dry regions compared to wet regions. Yet over dry regions, these two effects lead to an overestimation of the trends in rainfall totals and extremes by 40.3% and 33.2%, respectively (Fig. 2, Tables in Appendix B, reference: 1951-2010). In contrast, in wet regions these effects roughly cancel each other out in the case of extremes (+8.7%) and lead to an underestimation of the increase in total rainfall (-13.7%).
2 On the definition of a dry region
Climatological dryness is generally not determined by water supply alone as implied by the definition in Donat et al. (2016) but also depends on atmospheric water demand, that is the ability to evaporate water from the land surface (Köppen, 1900). This means that "we cannot tell whether a climate is moist or dry by knowing precipitation alone; we must know whether precipitation is greater or less than potential evapotranspiration", as Charles Warren Thornthwaite put it in a landmark paper (Thornthwaite, 1948); a statement that is clearly mirrored in present-day literature (e.g. Greve et al., 2014). For the analysis of precipitation extremes (Rx1d), regions in northern Europe such as parts of Scandinavia or the Netherlands fall in the ‘dry’ class because of relatively small annual maximum daily rainfall (Figure in Appendix A and Fig. 1c in Donat et al. (2016)). This is in contrast to what is commonly understood by the term ‘dry’ and consequently induces considerable confusion regarding increases in rainfall extremes and totals in globally ‘dry’ regions. To clarify this issue, we also test the sensitivity of the reported increases in Rx1d to the choice of dryness definition by using a variety of different dryness definitions (Figure in Appendix A). Hence, we evaluate trends and period increments in Rx1d and PRCPTOT in
regions that fall below the global 30% quantile in HadEX2 in the respective diagnostic (Rx1d or PRCPTOT), following Donat et al. (2016),
dry regions as identified from an aridity-based definition of dryness (Greve et al., 2014), and
dry and transitional regions combined from the latter definition (Greve et al., 2014).
In addition, we test uncertainties related to the temporal coverage of the dataset by relying on time series with 90% coverage (cf. Donat et al. (2016)) and in addition we analyse only time series without missing values (100% coverage).
Our results show that, based on a two-sided trend test, no significant increases in precipitation totals and extremes can be detected in the world’s dry regions as defined based on either Greve et al. (2014) or Köppen (1900), or in dry and transitional regions combined or if time series with incomplete temporal coverage are removed (Fig. 3, Table 1-2). These results are consistent with earlier studies that report no or modest changes in Rx1d and PRCPTOT in (predominantly dry) subsidence regions based on model simulations (Kharin et al., 2007; Fischer and Knutti, 2015) and in observations for individual subtropical regions such as Australia or the Mediterranean (Westra et al., 2013; Lehmann et al., 2015). In contrast, if ‘the world’s dry regions’ are defined based on falling below a global 30% threshold in Rx1d or PRCPTOT in the HadEX2 dataset (Donat et al., 2016), we indeed confirm robust increases in both Rx1d and PRCPTOT. Thus, the reported increases in both diagnostics are highly sensitive to the definition of a ‘dry region’, and appear to stem from regions with relatively moderate extreme (Rx1d) or average (PRCPTOT) rainfall, such as regions in Northern Europe (Rx1d, Figure in Appendix A) or North-East Siberia (PRCPTOT, Figure in Appendix A) that are not commonly understood as ‘dry’. Hence, we conclude extreme rainfall (Rx1d) or annual totals (PRCPTOT) do not show significant increases in the world’s dry (i.e. arid, or arid-transitional) regions in conventional definitions.
An accurate quantification of spatially aggregated trends in a rapidly changing Earth system is of highest relevance to decision-makers in the context of climate change adaptation (Fischer and Knutti, 2015). Therefore, short reference periods that are defined on a subset of the available dataset for normalization or data pre-processing purposes (Zhang et al., 2005; Sippel et al., 2015) should be avoided, as this procedure inevitably introduces biases. Furthermore, the definition of a ‘dry region’ induces considerable uncertainty in quantifying changes in rainfall extremes or totals. We find no evidence for changes in observed annual rainfall or heavy precipitation in the world’s dry or dry-transitional regions if the notion of dryness is based on water supply and demand (i.e. aridity) and biases due to short reference periods are avoided. Thus, understanding and disentangling the discrepancy in rainfall trends between regions of moderate extreme or total rainfall (Donat et al., 2016) and data-scarce arid regions, and relating both to model simulations, remains a research priority of high relevance.
S.S. and J.Z. conceived the study. All authors contributed to writing the paper.
Acknowledgements.S.S. and M.D.M. are grateful to the European Commission for funding the BACI project (grant agreement No 640176) and to the European Space agency for funding the STSE project CAB-LAB.
- Allen and Ingram (2002) Allen, M. R. and Ingram, W. J.: Constraints on future changes in climate and the hydrologic cycle, Nature, 419, 224–232, 2002.
- Barnett et al. (2005) Barnett, A. G., van der Pols, J. C., and Dobson, A. J.: Regression to the mean: what it is and how to deal with it, International journal of epidemiology, 34, 215–220, 2005.
- Donat et al. (2013) Donat, M., Alexander, L., Yang, H., Durre, I., Vose, R., Dunn, R., Willett, K., Aguilar, E., Brunet, M., Caesar, J., et al.: Updated analyses of temperature and precipitation extreme indices since the beginning of the twentieth century: the HadEX2 dataset, Journal of Geophysical Research: Atmospheres, 118, 2098–2118, 2013.
- Donat et al. (2016) Donat, M. G., Lowry, A. L., Alexander, L. V., O’Gorman, P. A., and Maher, N.: More extreme precipitation in the world’s dry and wet regions, Nature Climate Change, 2016.
- Fischer and Knutti (2015) Fischer, E. M. and Knutti, R.: Anthropogenic contribution to global occurrence of heavy-precipitation and high-temperature extremes, Nature Climate Change, 5, 560–564, 2015.
- Galton (1886) Galton, F.: Regression towards mediocrity in hereditary stature., The Journal of the Anthropological Institute of Great Britain and Ireland, 15, 246–263, 1886.
- Greve et al. (2014) Greve, P., Orlowsky, B., Mueller, B., Sheffield, J., Reichstein, M., and Seneviratne, S. I.: Global assessment of trends in wetting and drying over land, Nature Geoscience, 7, 716–721, 2014.
- Ingram (2016) Ingram, W.: Extreme precipitation: Increases all round, Nature Climate Change, 6, 443–444, 2016.
- IPCC (2012) IPCC: Summary for Policymakers, in: Managing the risks of extreme events and disasters to advance climate change adaptation: special report of the intergovernmental panel on climate change, edited by Field, C., Barros, V., Stocker, T., Dahe, Q., Dokken, D., Ebi, K., Mastrandrea, M., Mach, K., Plattner, G., Allen, S., Tignor, M., and Midgley, P., Cambridge University Press, 2012.
- Kharin et al. (2007) Kharin, V. V., Zwiers, F. W., Zhang, X., and Hegerl, G. C.: Changes in temperature and precipitation extremes in the IPCC ensemble of global coupled model simulations, Journal of Climate, 20, 1419–1444, 2007.
- Köppen (1900) Köppen, W.: Versuch einer Klassifikation der Klimate, vorzugsweise nach ihren Beziehungen zur Pflanzenwelt, Geographische Zeitschrift, 6, 593–611, 1900.
- Kottek et al. (2006) Kottek, M., Grieser, J., Beck, C., Rudolf, B., and Rubel, F.: World map of the Köppen-Geiger climate classification updated, Meteorologische Zeitschrift, 15, 259–263, 2006.
- Lehmann et al. (2015) Lehmann, J., Coumou, D., and Frieler, K.: Increased record-breaking precipitation events under global warming, Climatic Change, 132, 501–515, 2015.
- O’Gorman (2015) O’Gorman, P. A.: Precipitation extremes under climate change, Current climate change reports, 1, 49–59, 2015.
- Sippel et al. (2015) Sippel, S., Zscheischler, J., Heimann, M., Otto, F. E., Peters, J., and Mahecha, M. D.: Quantifying changes in climate variability and extremes: Pitfalls and their overcoming, Geophysical Research Letters, 42, 9990–9998, 2015.
- Thornthwaite (1948) Thornthwaite, C. W.: An approach toward a rational classification of climate, Geographical review, 38, 55–94, 1948.
- Tollefson (2016) Tollefson, J.: Global warming already driving increases in rainfall extremes, Nature, pp. 443–444, doi:10.1038/nature.2016.19508, 2016.
- Westra et al. (2013) Westra, S., Alexander, L. V., and Zwiers, F. W.: Global increasing trends in annual maximum daily precipitation, Journal of Climate, 26, 3904–3918, 2013.
- Zhang et al. (2005) Zhang, X., Hegerl, G., Zwiers, F. W., and Kenyon, J.: Avoiding inhomogeneity in percentile-based indices of temperature extremes, Journal of Climate, 18, 1641–1651, 2005.
|\tophlineDry Region Definition||Ref. Period||Temporal Coverage (%)||Period Inc. (%)||Trend Slope (% per year)||two-sided p-value (one-sided)||Sample size|
|\middlehlineDonat et al. (2016), global 30% quantile in Rx1d||1951-1980||90%||()||132|
|Donat et al. (2016), global 30% quantile in Rx1d||1951-1980||100%||()||57|
|Donat et al. (2016), global 30% quantile in Rx1d||1951-2010||90%||()||132|
|Donat et al. (2016), global 30% quantile in Rx1d||1951-2010||100%||()||57|
|Köppen (1900), dry climates (‘B-climates’)||1951-2010||90%||()||71|
|Köppen (1900), dry climates (‘B-climates’)||1951-2010||100%||()||36|
|Greve et al. (2014), dry regions||1951-2010||90%||()||73|
|Greve et al. (2014), dry regions||1951-2010||100%||()||40|
|Greve et al. (2014), dry+transitional regions||1951-2010||90%||()||118|
|Greve et al. (2014), dry+transitional regions||1951-2010||100%||()||61|
|\tophlineDry Region Definition||Ref. Period||Temporal Coverage (%)||Period Inc. (%)||Trend Slope (% per year)||two-sided p-value (one-sided)||Sample size|
|\middlehlineDonat et al. (2016), global 30% quantile in PRCPTOT||1951-1980||90%||()||299|
|Donat et al. (2016), global 30% quantile in PRCPTOT||1951-1980||100%||()||108|
|Donat et al. (2016), global 30% quantile in PRCPTOT||1951-2010||90%||()||299|
|Donat et al. (2016), global 30% quantile in PRCPTOT||1951-2010||100%||()||108|
|Köppen (1900), dry climates (‘B-climates’)||1951-2010||90%||()||183|
|Köppen (1900), dry climates (‘B-climates’)||1951-2010||100%||()||119|
|Greve et al. (2014), dry regions||1951-2010||90%||()||193|
|Greve et al. (2014), dry regions||1951-2010||100%||()||120|
|Greve et al. (2014), dry+transitional regions||1951-2010||90%||()||296|
|Greve et al. (2014), dry+transitional regions||1951-2010||100%||()||205|