# Huge decrease of frost frequency in the Mont-Blanc Massif under climate change

Benjamin Pohl, Daniel Joly, Julien Pergaud, Jean-François Buoncristiani, Paul Soare, Alexandre Berger

► To cite this version:

Benjamin Pohl, Daniel Joly, Julien Pergaud, Jean-François Buoncristiani, Paul Soare, et al.. Huge decrease of frost frequency in the Mont-Blanc Massif under climate change. Scientific Reports, 2019, 9 (1), pp.4919. 10.1038/s41598-019-41398-5. hal-02074047

HAL Id: hal-02074047

https://hal.science/hal-02074047

Submitted on 20 Mar 2019

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
SCIENTIFIC REPORTS

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

OPEN

# Huge decrease of frost frequency in the Mont-Blanc Massif under climate change

Benjamin Pohl<fnref n="1" />, Daniel Joly², Julien Pergaud<fnref n="1" />, Jean-François Buoncristiani<fnref n="1" />, Paul Soare<fnref n="1" /> & Alexandre Berger<fnref n="1" />

Received: 22 June 2018

Accepted: 6 March 2019

Published online: 20 March 2019

Mountains are a sensitive indicator of climate change and these areas are an early glimpse of what could happen in lowland environments. Peaking at 4808 m asl, the Mont-Blanc summit, at the boundary between France and Italy, is the highest of the Alps, in Western Europe. Its Massif is world-famous for outdoor and extreme sport activities, especially since the 1924 Olympic games held in Chamonix. Here, we use a novel statistical downscaling approach to regionalize current and future climate change over the Mont-Blanc Massif at an unequalled spatial resolution of 200 m. The algorithm is applied to daily minimum and maximum temperature derived from global climate models used in the fifth assessment report of the International Panel on Climate Change (IPCC). This new high-resolution database allows for a precise quantification of frost occurrence and its evolution until 2100. In the winter season and by the end of the 21st century, under a pessimistic scenario (RCP8.5), frost frequency in the morning could decrease by 30-35 percentage points in the valley of Chamonix, and in the afternoon, similar changes could occur for elevations comprised between 2000 and 3000 m. In summertime, changes are even larger, reaching a huge drop of 45-50 points in the afternoon between 3500 and 4500 m. These changes are much reduced under an optimistic scenario. They could have huge impacts on the environment (glacier shrinking, permafrost degradation, floods, changes in the distribution of species and ecosystems) and societies (summer tourism for climbing and hiking, and winter tourism for skiing).

The vulnerability of temperate mountain ranges to global warming is partly related to their elevation<fnref n="1" />,², steep slopes and orientation towards the sun<fnref n="1" />. The Mont-Blanc massif is one of the highest in the Alps (Fig. 1): 38 km² of its surface shows elevations of at least 3500 m, 8 km² peak above 4000 m but only 1 km² above 4500 m. This raises the question of the availability of so-called climatic shelters, that is, areas where the ecosystems can migrate to remain in environmental conditions favourable to their development³. In the plains of the mid-latitudes, a rise of 0.6 °C typically corresponds to a poleward shift of one latitude degree (about 111 km). In the temperate mountain ranges, this corresponds to a vertical shift of 100 m only⁴. A warming of several Celsius degrees, as estimated by the fifth assessment report of the IPCC, could thus lead to a general migration of several hundreds of meters vertically, and cause a dramatic transformation of high-mountain environments². In high elevations, temperature remains negative almost constantly, favouring snow accumulation and permafrost in non glaciated areas⁵,⁶. This is likely to change in the future, hence the need for very high-resolution climate projections.

Regionalizing climate change⁷,⁸ in mountain regions is a challenging task, because reliefs are drastically smoothed in the coarse grids of current climate models⁹,¹⁰. At spatial resolutions of a few tens of kilometres, the Alps barely reach 2000 or 2500 m above sea level, which is not suited to estimate the interaction between the topography and climate variability and change in these regions¹¹,¹². Statistical downscaling algorithms¹³,¹⁴ (see Supplementary Information) are appropriate tools to estimate the spatial variability of variables like air temperature, using spatial predictands derived from high-resolution digital elevation models, and large-scale temperature changes (predictors) taken from climate models. Here, we use a novel algorithm, allowing for more precise and refined temperature estimations⁴, to regionalize climate change as simulated by a sample of 13 climate models (Supplementary Table 1) at a spatial resolution of 200 m. The two most contrasted greenhouse gas emission scenarios are selected, namely RCP2.6 (optimistic) and RCP8.5 (pessimistic), to consider the whole spread of possible climate evolutions, together with historical simulations conducted with observed concentrations of

<footnote n="1">¹Biogéosciences, UMR6282 CNRS/université Bourgogne Franche-Comté, Dijon, France. ²ThéMA, UMR6049 CNRS/université Bourgogne Franche-Comté, Besançon, France. Correspondence and requests for materials should be addressed to B.P. (email: benjamin.pohl@ubfc.fr)</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Presentation of the Mont-Blanc Massif. (a) The Mont-Blanc summit seen from Emosson Dam (French-Swiss border). (b) 3D view of the Mont-Blanc Massif.</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Average frost occurrence (%, see colorbar) for Tn and Tx and for summer and winter seasons under current climate conditions (period 1979–2014) according to spatially interpolated observations. The black star corresponds to the Chamonix Mont-Blanc station in France (1035 m), the red star to the Mont-Blanc summit (4810 m) and the green star to the Courmayeur station in Italy (1224 m).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Spatially-averaged temperature evolutions over the Mont-Blanc Massif for Tn and Tx and for summer and winter seasons. The spatial domain used for computation is that shown in Fig. 2. Black curves: spatially interpolated observations, period 1979–2014. Blue colours: historical (HIST) simulations, period 1850–2005. The solid curve shows the ensemble mean, the colour shading extends to ±1 standard deviation to show model uncertainties. Green colours: the same for RCP2.6 simulations, period 2006–2100. Red colours: the same for RCP8.5 simulations, period 2006–2100.</figcaption></figure>

greenhouse gases since the 1850s. This downscaling procedure is necessary to consider the combined effects of elevation, slopes, orientation and valleys, on daily minimum (Tn) and maximum (Tx) air temperature. In the resulting elevation model (Fig. 1b), the Mont-Blanc summit peaks at 4763 m instead of 4808 m in the real world, illustrating the high degree of realism reached by such hectometric resolution.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Evolution of frost occurrence (percentage points) for the summer season. Maps represent the difference between RCP and HIST simulations (averaged over period 1970–2000) for mid-century (period 2036–2065) and late century (2070–2100). Maps are multi-model averages for Tn and Tx and for RCP2.6 and RCP8.5.</figcaption></figure>

## Frost Occurrence: a Relevant Metric for Impacts on Mountain Environment and Societies

We focus here on frost frequency evolutions, since this metric is of major importance for the rain-snow line, glacier melting or ice accumulation, the hydrology of the region and rock wall erosion though cryoclasty. Variability and long-term changes in frost frequency are thus accompanied by dramatic effects on the environment³, landscape¹⁵,¹⁶, and human activities¹⁷,¹⁸. Under current climate conditions (Fig. 2), in the morning (Tn) during the winter season, the 0°C threshold is almost never exceeded, except in the bottom of the lower valleys like Chamonix or Courmayeur (<1200–1500 m). In the afternoon (Tx), the air mass warms and the occurrence of frost strongly decreases under 1500 and 2000 m asl. In summer frost is circumscribed to higher elevations,
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Boxplot representation of the evolution of frost occurrence (percentage points) for the summer period, according to elevation. Analyses are the same as the maps of Fig. 4. The boxes have lines at the lower, median and upper quartile values. The whiskers are lines extending from each end of the box to 1.5 interquartile range. Plus signs correspond to statistical outliers.</figcaption></figure>

and frost predominates above 2500–3000 m asl for Tn and 3500–4000 m asl for Tx. These estimates (obtained as multiple regression based on observed temperature measurements⁴: the “interpolation” step⁴), are associated with very weak errors and uncertainties, and can thus be considered as reliable.

The next step consists in applying the same spatial models to downscale climate change simulations (the “downscaling” *per se*⁴). In such temperature estimates, the spatial variability is thus inherited from the statistical laws deduced, under current climate conditions, from the aforementioned interpolations, while the large-scale temporal variability of air temperature is derived from the global models involved in the Coupled Model Inter-comparison Project phase 5¹⁹ (CMIP5). Using the same 200-m horizontal projected grid over the Mont-Blanc Massif and its surroundings, summer (JJA) and winter (DJF) mean temperature evolutions from 1850 to 2100 are compared to interpolated observations (Fig. 3). Our data clearly show the heat wave for summer of 2003²⁰, while the trimester of DJF 1989–1990 holds the record for winter temperature over the period. Over their common period, the difference between observations and historical simulations allows for an estimate of the typical errors of climate models over this region (Supplementary Figure 2). Models over-estimate Tn both in summer and winter, and produce much more realistic results for Tx estimates. The amplitude of the future warming is strongly dependent on the emission scenario. Optimistic RCP2.6 produces stabilized temperatures at about +1.5°C compared to pre-industrial levels: it is the only scenario that meets the terms of the Paris Agreement on climate change locally, that is, a “temperature rise this century well below 2 Celsius degrees above pre-industrial levels”. Pessimistic RCP8.5 leads to a warming of roughly +2°C in wintertime (+2.8°C in summertime) by the mid-century, these values reaching respectively +3.8°C and +5.6°C for the late century. Although all downscaled models produce qualitatively similar results, uncertainties remain large, as denoted by the inter-model spread in
Fig. 3. For the late 21st century for instance, and from one model to another, summer temperature rise is comprised between +5 and nearly +9°C for Tx (+4°C to +7.5°C for Tn). Depending on the models, the record heat wave of 2003 could become representative of a typical summer season between 2060 and 2075. The abnormally warm winter of 1989–1990 could represent the typical winters of the years 2070–2080.

## Evolution of frost occurrence by the mid- and late-century

Results are next refined at the daily timescale, through an analysis of daily frost occurrence, to represent the evolution in the summertime frost (Fig. 4). This is done under contrasted emission scenarios and for the mid- and late century, compared to current climate conditions (as described by historical simulations). Errors inherited from the CMIP5 are shown in Supplementary Figure 2, corresponding uncertainties in Supplementary Figure 3, and similar analyses applied to the winter season in Supplementary Figure 4. Projected changes are of larger magnitude for Tn, because most of the domains already experiences positive temperature in summer under current climate conditions (Fig. 2). They are much sensitive to the emission scenario, especially for the late century. Morning frost occurrence (Tn) decreases by about 25–30 pp over most of the Mont-Blanc massif in the mid-century under RCP8.5, and slightly less than that for the late century under RCP2.6. Frost occurrence experiences a huge drop of 45–50 pp for the late century for RCP8.5. The geography of these changes is strongly driven by the topography (Fig. 1). Crucially, even the Mont-Blanc summit shows obvious evolutions: up to –20 pp by the end of the century for frost occurrence in the morning (Tn), and even –35 pp in the afternoon (Tx), according to the RCP8.5 scenario. Such decrease may lead to major changes of its local climate. Under current climate conditions, it has only experienced brief thaw events²¹, like during the heat wave of summer 2003. This could favour negative mass balance for most glaciers²², due to drastically enhanced melting and decreased accumulation in high elevation environments.

The present altitudinal range of frost occurrence will show major changes during the 21st century (Fig. 5). Main changes occur between 2000 and 4000 m asl for RCP2.6 in the morning (Tn), and between 3000 and 4500 m in the afternoon (Tx), according to RCP2.6 downscaled simulations. They are much higher in elevation and magnitude for RCP8.5 (up to –45–50 pp between 3500 and 4000 m asl for Tn and –30 to –50 pp at the same elevation for Tx by the late century).

These evolutions in the climate of the Massif are likely to induce major changes in the dynamics of the snow-pack²³, cryosphere²⁴, hydrosphere²⁵ and biosphere¹⁵. They could lead to major, but complex changes in hydrological regimes, especially in the upper alpine valleys, through changes in the relative contribution of their glacial and nival components²⁶. Increased frequency in frost-thawing alternations could also enhance gelifraction, which, together with degradations in the permafrost¹⁸,²⁷, could increase rockfall occurrence²⁸⁻³⁰ in the future. These changes might also force human activities to adapt to such shifts of the isotherms and rain-snow line. Major impact are expected on winter sports, leading ski resorts³¹,³² to higher elevations in parallel with the migration of biomes³³ and environments. Given the narrowness and vulnerability of high mountain environments, this could potentially lead land-use conflicts. Increased anthropogenic pressure on these vulnerable areas has the potential to further perturb their fragile ecosystems, highlighting the need for sustainable development and planning strategies.

## Methods

A frost event is defined as a day for which Tn < 0°C for morning frost, or Tx < 0°C for afternoon (permanent) frost. Frost frequency (in %) is obtained as the ratio between the number of frost events and the total length of the period. Frost evolutions under climate change consist in differences between future and current frost frequencies and are thus expressed as percentage points (pp). The statistical downscaling⁴, as used here, combines temporal variability derived from CMIP5 simulations, and spatial variability based on statistical relationships obtained on observed time series (an original database of 60 weather stations available daily over a 35-year period). Temperature biases are inherited from both the CMIP5 model and the statistical downscaling: see Supplementary Information and methodological paper⁴ for details. Uncertainties arise from disagreements between CMIP5 models, hence the use of 13 different simulations for robust convergence assessment.

## References

1. Beniston, M. Climatic Change in Mountain Regions: A Review of Possible Impacts. *Clim. Chang.* **59**, 5–31 (2003).

2. Gobiet, A. et al. 21st century climate change in the European Alps-A review. *Sci. Total Environ.* **493**, 1138–1151 (2014).

3. Rogora, M. et al. Assessment of climate change effects on mountain ecosystems through a cross-site analysis in the Alps and Apennines. *Sci. Total Environ.* **624**, 1429–1442 (2018).

4. Joly, D. et al. Geomatic downscaling of temperatures in the Mont Blanc massif. *Int. J. Climatol.* **38**, 1846–1863 (2018).

5. Haeberli, W. Mountain permafrost - research frontiers and a special long-term challenge. *Cold Reg. Sci. Technol.* **96**, 71–76 (2013).

6. Marty, C., Schlögl, S., Bavay, M. & Lehning, M. How much can we save? Impact of different emission scenarios on future snow cover in the Alps. *Cryosph.* **11**, 517–529 (2017).

7. Giorgi, F. et al. Enhanced summer convective rainfall at Alpine high elevations in response to climate warming. *Nat. Geosci.* **9**, 584 (2016).

8. Torma, C., Giorgi, F. & Coppola, E. Added value of regional climate modeling over areas characterized by complex terrain—Precipitation over the Alps. *J. Geophys. Res. Atmos.* **120**, 3957–3972 (2015).

9. Suklitsch, M., Gobiet, A., Leuprecht, A. & Frei, C. High resolution sensitivity studies with the regional climate model CCLM in the Alpine region. *Meteorol. Zeitschrift* **17**, 467–476 (2008).

10. Montesarchio, M., Zollo, A. L., Bucchignani, E., Mercogliano, P. & Castellari, S. Performance evaluation of high-resolution regional climate simulations in the Alpine space and analysis of extreme events. *J. Geophys. Res. Atmos.* **119**, 3222–3237 (2014).

11. Salzmann, N. et al. Ground surface temperature scenarios in complex high-mountain topography based on regional climate model results. *J. Geophys. Res.* **112**, F02S12 (2007).

12. Suklitsch, M. et al. Error characteristics of high resolution regional climate models over the Alpine area. *Clim. Dyn.* **37**, 377–390 (2011).
13. Le Roux, R., Katurji, M., Zawar-Reza, P., Quénol, H. & Sturman, A. Comparison of statistical and dynamical downscaling results from the WRF model. *Environ. Model. Softw.* **100**, 67–73 (2018).

14. Roberts, D. R., Wood, W. H. & Marshall, S. J. Assessments of downscaled climate data with a high-resolution weather station network reveal consistent but predictable bias. *Int. J. Climatol*, https://doi.org/10.1002/joc.6005 (2019).

15. Theurillat, J. & Guisan, A. Potential Impact of Climate Change on Vegetation in the European Alps: A Review. *Clim. Change* **50**, 77–109 (2001).

16. Pauli, H., Gottfried, M. & Grabherr, G. Effects of climate change on the alpine and nival vegetation of the Alps. *J. Mt. Ecol.* **7**, 9–12 (2003).

17. Gilaberte-Búrdalo, M., López-Martín, F., Pino-Otín, M. R. & López-Moreno, J. I. Impacts of climate change on ski industry. *Environ. Sci. Policy* **44**, 51–61 (2014).

18. Duvillard, P.-A., Ravanel, L., Deline, P. & Dubois, L. Paraglacial Rock Slope Adjustment Beneath a High Mountain Infrastructure — The Pilatte Hut Case Study (Écrins Mountain Range, France). *Front. Environ. Sci.* **6**, 94 (2018).

19. 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).

20. Casty, C., Wanner, H., Luterbacher, J., Esper, J. & Böhm, R. Temperature and precipitation variability in the European Alps since 1500. *Int. J. Climatol.* **25**, 1855–1880 (2005).

21. Schiermeier, Q. Alpine thaw breaks ice over permafrost's role. *Nature* **424**, 712 (2003).

22. Six, D. & Vincent, C. Sensitivity of mass balance and equilibrium-line altitude to climate change in the French Alps. *J. Glaciol.* **60**, 867–878 (2014).

23. Piazza, M. et al. Projected 21st century snowfall changes over the French Alps and related uncertainties. *Clim. Change* **122**, 583–594 (2014).

24. Stoffel, M., Tiranti, D. & Huggel, C. Climate change impacts on mass movements - Case studies from the European Alps. *Sci. Total Environ.* **493**, 1255–1266 (2014).

25. Beniston, M. & Stoffel, M. Rain-on-snow events, floods and climate change in the Alps: Events may increase with warming up to 4°C and decrease thereafter. *Sci. Total Environ.* **571**, 228–236 (2016).

26. Bavay, M., Grünewald, T. & Lehning, M. Advances in Water Resources Response of snow cover and runoff to climate change in high Alpine catchments of Eastern Switzerland. *Adv. Water Resour.* **55**, 4–16 (2013).

27. Magnin, F. et al. Modelling rock wall permafrost degradation in the Mont Blanc massif from the LIA to the end of the 21st century. *Cryosphere* **11**, (2017).

28. Ravanel, L., Magnin, F. & Deline, P. Science of the Total Environment Impacts of the 2003 and 2015 summer heatwaves on permafrost-affected rock-walls in the Mont Blanc massif. *Sci. Total Environ.* **609**, 132–143 (2017).

29. Allen, S. & Huggel, C. Extremely warm temperatures as a potential cause of recent high mountain rockfall. *Glob. Planet. Change* **107**, 59–69 (2013).

30. Luigi, S. & Guzzetti, F. Earth-Science Reviews Landslides in a changing climate. *Earth Sci. Rev.* **162**, 227–252 (2016).

31. Gonseth, C. Impact of snow variability on the Swiss winter tourism sector: implications in an era of climate change. *Clim. Change* **119**, 307–320 (2013).

32. François, H., Morin, S., Lafaysse, M. & George-Marcelpoil, E. Crossing numerical simulations of snow conditions with a spatially-resolved socio-economic database of ski resorts: A proof of concept in the French Alps. *Cold Reg. Sci. Technol.* **108**, 98–112 (2014).

33. Chamberlain, D. E., Negro, M., Caprio, E. & Rolando, A. Assessing the sensitivity of alpine birds to potential future changes in habitat and climate to inform management strategies. *Biol. Conserv.* **167**, 127–135 (2013).

## Acknowledgements

We acknowledge the World Climate Research Programme's Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Supplementary Table 1 of this paper) for producing and making available their model output. Calculations were performed using HPC resources from DNUM CCUB (Centre de Calcul de l'Université de Bourgogne). This work is part of the ANR project VIP-Mont Blanc (grant ANR-14-CE03-0006-04).

## Author Contributions

B.P. designed the research; D.J. performed the statistical downscaling; P.S., J.P. and B.P. analyzed the data; J.P. and J.F.B. plotted the figures; all authors contributed to write the manuscript.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-41398-5.

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

**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 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) 2019