Journal cover Journal topic
The Cryosphere An interactive open-access journal of the European Geosciences Union
Journal topic
The Cryosphere, 13, 1325–1347, 2019
https://doi.org/10.5194/tc-13-1325-2019
The Cryosphere, 13, 1325–1347, 2019
https://doi.org/10.5194/tc-13-1325-2019

Research article 24 Apr 2019

Research article | 24 Apr 2019

# Winter tourism under climate change in the Pyrenees and the French Alps: relevance of snowmaking as a technical adaptation

Winter tourism under climate change in the Pyrenees and the French Alps: relevance of snowmaking as a technical adaptation
Pierre Spandre1,2, Hugues François1, Deborah Verfaillie2,3, Marc Pons4, Matthieu Vernay2, Matthieu Lafaysse2, Emmanuelle George1, and Samuel Morin2 Pierre Spandre et al.
• 1Univ. Grenoble Alpes, Irstea, UR LESSEM, Grenoble, France
• 2Univ. Grenoble Alpes, Université de Toulouse, Météo-France, CNRS, CNRM, Centre d'Études de la Neige, 38000 Grenoble, France
• 3Barcelona Supercomputing Center, Barcelona, Spain
• 4Snow and Mountain Research Center of Andorra, IEA, Sant Julià de Lòria, Andorra

Correspondence: Samuel Morin (samuel.morin@meteo.fr)

Abstract

Climate change is increasingly regarded as a threat for winter tourism due to the combined effect of decreasing natural snow amounts and decreasing suitable periods for snowmaking. The present work investigated the snow reliability of 175 ski resorts in France (Alps and Pyrenees), Spain and Andorra under past and future conditions using state-of-the-art snowpack modelling and climate projections using Representative Concentration Pathways RCP2.6, RCP4.5 and RCP8.5. The natural snow reliability (i.e. without snowmaking) elevation showed a significant spatial variability in the reference period (1986–2005) and was shown to be highly impacted by the ongoing climate change. The reliability elevation using snowmaking is projected to rise by 200 to 300 m in the Alps and by 400 to 600 m in the Pyrenees in the near future (2030–2050) compared to the reference period for all climate scenarios. While 99 % of ski lift infrastructures exhibit adequate snow reliability in the reference period when using snowmaking, a significant fraction (14 % to 25 %) may be considered in a critical situation in the near future. Beyond the mid-century, climate projections highly depend on the scenario with either steady conditions compared to the near future (RCP2.6) or continuous decrease in snow reliability (RCP8.5). Under RCP8.5, our projections show that there would no longer be any snow-reliable ski resorts based on natural snow conditions in the French Alps and Pyrenees (France, Spain and Andorra) at the end of the century (2080–2100). For this time period and this scenario, only 24 resorts are projected to remain reliable with snowmaking, all being located in the Alps.

1 Introduction

The ongoing evolution of natural snow conditions related to climate change is increasingly regarded as a major threat for winter tourism . This prompts the question of how climate change affects ski resorts and the relevance of snowmaking as an adaptation measure . Initial studies in the early 2000s quantified the snow reliability of ski resorts based on the “100 days” rule, later considered as the reference approach for investigations of climate-induced impacts on the winter tourism . This rule states that a ski resort is snow reliable if the snow depth exceeds 30 cm during 100 d or more, which provides objective information when comparing distinct periods (past and future) or locations . The snow reliability line is defined as the elevation above which these conditions are met, allowing the assessment of the reliability of a ski resort by comparing its elevation to the snow reliability line .

Most investigations based on the 100 days rule used single-point representations of ski slopes to assess the snow and meteorological conditions of a given ski resort, often using the median elevation of a ski resort defined as the average of summit and base elevations . and used the “highest”, “middle” and “lowest” elevations of the study area while mixed various approaches by considering either a single point or three distinct elevations for each ski resort. Alternatively, considered the lowest and highest elevations, between which 75 % of the ski slope surface area was concentrated. These remain coarse representations limiting the analysis of the situation of a ski resort to a binary conclusion of reliable or unreliable . and in Switzerland and later in the rest of the European Alps based their analysis on the natural snow conditions. reviewed the existing literature to address the snow reliability line for regions of Europe (Austria, Italy, Germany, Switzerland and France) based on distinct methods and reference periods . They concluded that 91 % of the 666 ski resorts in the European Alps were snow reliable around 2005. Significant spatial variations in the snow reliability line were shown, ranging from 1050 to 1500 m above sea level (a.s.l.) with consequences on local reliability of ski resorts: 69 % of ski resorts were snow reliable in Germany and up to 97 % were snow reliable in Switzerland and France. similarly addressed the impact of climate change on the snow reliability line and concluded that under a +1C warming compared to present only 75 % of European Alps ski resorts would remain reliable and respectively 61 % and 30 % for +2 and +4C warming compared to the present would remain reliable. These investigations were limited to the analysis of natural snow using average conditions over large regions. later showed by the analysis of 52 climate stations in Austria over the 1981–2001 period that an elevation of 1200 m a.s.l. could not be confirmed as snow reliable for all regions of Tyrol (Austria). Using natural snow conditions to assess the snow reliability of ski resorts has also been questioned, due to the strong role of snow management, in particular grooming and snowmaking .

Recent studies have increasingly taken into account snow grooming and snowmaking . developed a simple modelling approach accounting for a required snow depth of 50 cm for skiing activities and computed snowmaking requirements based on this target. This method provided consistent season durations for the 1961–1990 reference period in the southern Ontario region (Canada) which were shown to significantly decrease under projected climate conditions despite an increasing need for snowmaking. later used this modelling approach and a 60 cm snow base depth requirement in the Quebec region (Canada). applied this method in Tyrol (Austria) and concluded that snowmaking could guarantee snow reliability at elevations above 1000 m a.s.l. for the 1971–2000 reference period and would remain a suitable adaptation method until the 2050s with a significant increase in water and energy requirements (Steiger2010). Similar investigations were conducted to assess the impact of climate change on the ski season duration and the snowmaking requirements so as to compensate for the loss in regions of Austria , Germany , Switzerland , Andorra , the Pyrenees , the northeastern US , New Zealand and Australia . Major limitations remain. First, little investigation was undertaken in France, still a major area for winter tourism . Second, meteorological and snow input data considered for the analysis were aggregated over large regions where high spatial variability can be observed . Third, snow conditions were often simulated using simplified degree day modelling approaches and neglected the differences between natural snow and groomed or machine-made snow properties .

The present work aims at producing snow reliability investigations of a wide range of ski resorts in France (Alps and Pyrenees), Spain and Andorra under past and future conditions using state-of-the-art snowpack modelling. We accounted for snow grooming and snowmaking using a detailed snowpack model and used adjusted and downscaled climate projections from the EURO-CORDEX dataset to compute snow reliability elevations with distinct levels of snow reliability requirements. The mean elevation of residential population in a ski resort and the mean elevation of ski lifts were compared to the snow reliability line. We defined seven distinct categories for ski resorts based on their natural snow reliability, their degree of dependence on snowmaking to achieve reliability and whether snowmaking may be a technically efficient method to guarantee snow reliability under present and future climate conditions.

2 Method

## 2.1 Ski resort definition and features

### 2.1.1 Definition of relevant elevations of ski resorts

All data on the geographical location and technical data on ski resorts were extracted from the “BD Stations” database . Ski lift installation and operation in France are supervised by the STRMTG (“Services Techniques de Remontées Mécaniques et Transports Guidés”). The STRMTG is a public service in charge of the safety control of French ski lifts providing authorizations for ski lift operations. The STRMTG manages a database (CAIRN: CAtalogue Informatisé des Remontées Mécaniques Nationales) dedicated to ski lifts which includes technical characteristics of each ski lift such as the ski lift power. The ski lift power is an indicator of the size of a ski lift, defined as the product of the elevation difference between the bottom and the top of a ski lift (in kilometres) and its capacity, i.e. the flow of persons per hour (persons h−1), expressed in persons km h−1. Ski lift infrastructures in France have a total ski lift power of 977 000 persons km h−1, 94 % of which are included in the present study (Appendix B). These data are complemented with geographical information from the database BD TOPO (25 m of resolution) developed by the French Geographical Institute (IGN, “Institut Géographique National”). The following elevations were used to be compared with the snow reliability line.

• The mean ski lift elevation of a ski resort is defined as the average of top and bottom elevations of each ski lift weighted by its ski lift power, referred to as the mean elevation of the ski resort .

• The village elevation of a ski resort is defined as the mean elevation of tourism housing infrastructure, where tourists stay during their ski holidays. It is computed using IGN data on the location and characteristics of buildings. Buildings located within 300 m from the bottom of the ski lifts are selected, and the selection procedure continues by iterations, using a 200 m radius around each identified building and so on, until no more buildings are found. We then compute the net floor surface area of each selected building (taking account the number of floors, based on building height), which is used to compute the weighted mean elevation of the built area associated with each resort, weighted by their net floor area .

Figure 1Relationship between OpenStreetMap (OSM) data on French ski resorts and ski lift power (SLP) and elevations of the ski resort (min, mean, max), used to estimate similar indicators in ski resorts in Spain and Andorra. In panel (b), “BD Stations” refers to the database of French ski resorts used in this study (see text for details).

Data for computing the mean ski lift elevation and village elevation of ski resorts using the method described above are only available for France. Another approach was required for addressing the characteristics of ski resorts in Andorra and Spain. Based on the OpenStreetMap (OSM) project (http://www.openstreetmap.org/, last access: 12 April 2019), we estimated the main features for the Spanish and Andorran ski resorts (village elevation, ski lift mean elevation and ski lift power). However, ski lift capacity is not included in OSM, and building height data are incomplete, which hampers our ability to proceed using the method developed for the French ski resorts. To circumvent this issue, we extracted all ski slopes from OSM for all ski resorts in France, Spain and Andorra, and computed linear regressions between information extracted from OSM and independent estimates for French ski resorts, which were then used to compute the indicators for ski resorts in Spain and Andorra (Fig. 1).

• The linear model of the ski lift power versus the OSM surface area (Fig. 1a) had a correlation coefficient R2=0.87 (p value $<{\mathrm{10}}^{-\mathrm{15}}$), proving relevant to estimate the ski lift power based on the OSM surface area.

• Elevations derived from the OSM spatial representation also proved significantly correlated to data from the BD Stations (Fig. 1b):

• All elevations together had a RMSD = 149 m and mean difference = 15 m, and the linear model of slope was 0.97 (R2=0.91, p value $<{\mathrm{10}}^{-\mathrm{15}}$).

• Mean elevation had a RMSD = 154 m and mean difference = 51 m, and the linear model of slope was 0.82 (R2=0.83, p value $<\phantom{\rule{0.125em}{0ex}}{\mathrm{10}}^{-\mathrm{15}}$).

• The village elevation proved significantly correlated to the mean elevation derived from OSM spatial representations (slope 0.64, intercept 326 m, R2=0.62, p value $<\phantom{\rule{0.125em}{0ex}}{\mathrm{10}}^{-\mathrm{15}}$). The linear model was applied to estimate the village elevation from the OSM mean elevation and compared to the BD Stations data: RMSD = 179 m; mean difference $<\phantom{\rule{0.125em}{0ex}}{\mathrm{10}}^{-\mathrm{12}}$ (Fig. 1b).

### 2.1.2 Study area

A sample of 175 ski resorts in the French Alps (n= 129), the French Pyrenees (n= 28), the Spanish Pyrenees (n= 14) and Andorra (n= 4) were included in the present study (Fig. 2, Appendix B). The French ski resorts included in this study (n= 157) represent 94 % of the national ski lift infrastructures. For Andorra, our study accounts for 100 % of the ski tourism infrastructures. For Spain, there are a total of 30 ski resorts, 14 of which are in the Spanish Pyrenees and considered in this study (note that in our study, the ski resorts La Molina and Masella were considered together). In terms of skiers, the Spanish Pyrenees represent around 63 % of the total ski market in Spain.

Figure 2(a) The 175 ski resorts covered by the present study and the 44 massifs from the SAFRAN reanalysis and (b) distribution of ski resort elevations depending on their location: northern Alps, southern Alps, French Pyrenees, and Spanish and Andorran Pyrenees (“Sp.–And. Pyrenees”). See Appendix B.

## 2.2 Definition and computation of the snow reliability line

### 2.2.1 Snowpack modelling

The “Crocus Resort” version of the multilayer snowpack model SURFEX/ISBA-Crocus was used in the present study . Crocus Resort allows us to take into account the effect of grooming and snowmaking on snow properties so as to provide simulations of snow conditions on ski slopes . The impacts of grooming are simulated and machine-made snow can be added to the snowpack specifying the precipitation rate ($\mathrm{1.2}×{\mathrm{10}}^{-\mathrm{3}}$ kg m−2 s−1; ) and conditions for triggering the production (wet-bulb temperature threshold −2C, target quantity or target snow depth). The production of snow was based on the following rules, dividing the winter season into distinct periods :

• Between 1 November and 15 December, a 30 cm deep “base layer” (snow mass of 150 kg m−2, for a typical snow density of 500 kg m−3) is produced, weather conditions permitting, regardless of natural snowfalls during the period.

• Between 15 December and 28 February, snow is produced, if meteorologically possible, so as to maintain a total snow depth of 60 cm.

• After 1 March, no more snow is produced.

### 2.2.2 Climate forcing data

The meteorological system SAFRAN provides meteorological data (temperature, precipitations, etc.) for mountain areas of an approximate 1000 km2 surface referred to as a “massif”, covering the French Alps and Pyrenees, including the Spanish and Andorran Pyrenees (Fig. 2). Within each massif, the meteorological conditions are considered to be homogeneous and to depend only on the elevation (by steps of 300 m) with a time resolution of 1 h. SAFRAN forcing data are available for the 1958–2015 period . Computations of snow conditions over the reference period using SAFRAN forcing data are further referred to as “SAFRAN” and can be considered as the reference observational dataset.

Table 1EURO-CORDEX GCM–RCM combinations used in this study (rows: RCMs; columns: GCMs), with the time period available for the HIST and RCP4.5 and RCP8.5 scenarios (RCPs). Model combinations additionally using RCP2.6 are displayed in bold. Contributing institutes are indicated inside parentheses – CLMcom: Climate Limited-area Modelling Community with contributions by BTU, DWD, ETHZ, UCD,WEGC; CNRM: Météo France; IPSL-INERIS: Institut Pierre Simon Laplace, CNRS, France – Laboratoire des Sciences du Climat et de l'Environnement, IPSL, CEA/CNRS/UVSQ – Institut National de l'Environnement Industriel et des Risques, Verneuil en Halatte, France; KNMI: Kingdom of Netherlands Meteorological Institute, Ministry of Infrastructure and the Environment; MPI-CSC: Max Planck Institute for Meteorology, Climate Service Center, Hamburg, Germany; SMHI: Swedish Meteorological and Hydrological Institute, Rossby Centre, Norrköping, Sweden.

This study uses the EURO-CORDEX dataset for climate projections consisting of six regional climate models (RCMs) forced by five different global climate models (GCMs) from the CMIP5 ensemble over Europe, for the historical, RCP2.6, RCP4.5 and RCP8.5 scenarios . All EURO-CORDEX data were adjusted using the ADAMONT method using the SAFRAN data as the reference observation dataset . Historical runs generally cover the period 1950–2005 and climate projections (RCPs) cover the period 2006–2100 (Table 1). Continuous hourly resolution meteorological time series derived from RCM output by the ADAMONT statistical adjustment method are then used as input of the SURFEX/ISBA-Crocus snowpack model .

### 2.2.3 Snow indicators

The snow reliability line was computed from the simulated snow conditions for the reanalysis and all GCM–RCM pairs and scenarios. The snow reliability line was based on the 100 days rule and defined for a given season as the elevation above which a minimum quantity of 100 kg m−2 of snow (i.e. 20 cm of snow at 500 kg m−3 density) was simulated during at least 100 d between 15 December and 15 April . The use of snow mass instead of snow depth appeared more relevant for our study, considering the differences between natural snow properties and machine-made snow or groomed snow . Based on the season length computed for SAFRAN massif elevations (300 m step), a linear interpolation was used to compute the snow reliability line meeting the 100 d threshold. In cases where the season length at the minimum (respectively maximum) elevation was longer (respectively shorter) than 100 d, the snow reliability line was set to half the altitudinal step (150 m) below (respectively above) the minimum (respectively maximum) elevation for a given massif. We further computed for each massif the snow reliability line by considering distinct periods, climate scenarios, snow requirements and snow management, providing 48 distinct values of the snow reliability elevation resulting from the combination of these parameters (Tables A1, A2, A3 in Appendix). Eight periods and scenario configurations are based on the reference period (1986–2005) using the SAFRAN reanalysis and available GCM–RCM pairs (HIST), the near future (2030–2050) and the end of the century (2080–2100), using climate scenarios RCP2.6, RCP4.5 and RCP8.5 for all available GCM–RCM pairs (Table 1). Three distinct levels of snow reliability requirements were defined as the elevation where the season length reached 100 d one season out of two (50 % percentile of annual values), seven seasons out of 10 (70 % percentile of annual values) and nine seasons out of 10 (90 % percentile of annual values). Last, we considered the groomed snow conditions (no snowmaking) and including snowmaking (two configurations). We did not compute indicators based on unmanaged natural snow conditions alone, i.e. without grooming or snowmaking.

## 2.3 Definition of snow reliability categories

Seven snow reliability categories have been designed with respect to the natural snow reliability and the relevance of snowmaking as an efficient adaptation method to reduce the effect of snow variability and scarcity, in line with previous investigations . Following , we considered a strict threshold of nine winters out of 10 for snowmaking reliability (90 % percentile of annual values), considering that snowmaking facilities are an investment for the operations of ski resorts and should therefore target a high level of reliability. The following categories were defined to characterize the snow reliability of ski resorts, depending on the relationship between village elevation and mean ski lift elevation, and the reliability lines with and without snowmaking. The village elevation is critical because this corresponds to the entry point of skiers to the ski slopes from their tourism housing infrastructure. This often corresponds to the lower elevation of the major ski lift infrastructure, which is a key area for snow managers because snow reliability there is both challenging and a strong asset for the ski resort operations . Categories are ordered by decreasing levels of natural and managed snow reliability. For each ski resort, its category corresponds to the first one for which the criterion is fulfilled, from category 1 until category 7. A ski resort fulfilling the condition of category N-1 also fulfills the condition of category N. Ski resorts in category N fulfill the condition of category N but not the condition of category N-1.

• Category 1: village elevation above the 90 % groomed snow reliability line.

• Category 2: village elevation above the 70 % groomed snow reliability line and village elevation above the 90 % snowmaking reliability line.

• Category 3: mean ski lift elevation above the 70 % groomed snow reliability line and village elevation above the 90 % snowmaking reliability line.

• Category 4: mean ski lift elevation above the 50 % groomed snow reliability line and village elevation above the 90 % snowmaking reliability line.

• Category 5: village elevation above the 90 % snowmaking reliability line.

• Category 6: mean ski lift elevation above the 90 % snowmaking reliability line.

• Category 7: mean ski lift elevation below the 90 % snowmaking reliability line.

Categories 1, 2 and 3 illustrate ski resorts where natural snow conditions are generally reliable . Snowmaking is generally employed only at the lowest elevations, and it makes a difference only for a minority of seasons when natural snow conditions are too scarce. Categories 4 and 5 illustrate ski resorts where natural snow conditions may not be considered as reliable as the previous categories, but snowmaking can generally guarantee the reliability in all elevations of the resort. In these two categories, snowmaking is useful and efficient in reducing natural snow scarcity at all elevations of the resort . Categories 6 and 7 illustrate ski resorts where natural snow conditions are generally not considered reliable and snowmaking is not efficient in reducing natural snow scarcity at the lowest elevations of the resort.

Figure 3Spatial variability between massifs and evolution for the reference period, the near future (2030–2050) and the end of the century (2080–2100) of the snow reliability line based on RCP2.6, RCP4.5 and RCP8.5 for the main areas covered in the present study (northern and southern Alps, French and Spanish–Andorran Pyrenees).

3 Results

## 3.1 Snow conditions and snow reliability line

### 3.1.1 Past climate conditions

Figure 3 shows that a significant spatial variability of the snow reliability line can be observed for the reference period (1986–2005). The median elevation of the 70 % groomed snow reliability ranges between 1750 m a.s.l. in the northern Alps, 2000 m a.s.l. in the French Pyrenees, 2250 m a.s.l. in the southern Alps, and up to 2300 m a.s.l. in the Spanish and Andorran Pyrenees (HIST, Fig. 3). Although a deviation can be observed, the spatial variability is consistent between climate model computations over the reference period (HIST) and the reference dataset (SAFRAN). The 90 % snow reliability using snowmaking is significantly lower than the 70 % groomed snow reliability line (Fig. 3). Due to snowmaking the median reliability elevation increases between 700 m in the French Pyrenees, 900 m in the Spanish and Andorran Pyrenees, 1000 m in the northern Alps, and up to 1200 m in the southern Alps. This results in a snowmaking reliability line significantly lower in the southern Alps compared to the Pyrenees despite poorer natural snow conditions (Fig. 3). Although the improvement of snow conditions thanks to snowmaking is lower in the Pyrenees compared to the Alps, the annual snowmaking requirements are higher, with 400 to 550 kg m−2 machine-made snow produced at the snow reliability line in the northern and southern Alps (10 %–90 % percentiles of annual values) and 400 to 700 kg m−2 in the French, Spanish and Andorran Pyrenees (HIST). Such production is equivalent to 80 cm to 1.1 m of snow in the Alps and 80 cm to 1.4 m of snow in the Pyrenees at the snow reliability line (using a machine-made snow density value of 500 kg m−3).

### 3.1.2 Future change in the near future (2030–2050)

Natural snow conditions are projected to be significantly affected by climate change in the near future (2030–2050) with similar evolution between climate scenarios (Fig. 3). The median 70 % groomed snow reliability line is projected to range between

• 1850 and 2000 m a.s.l. in the northern Alps (100 to 250 m above the reference period);

• 2500 and 2650 m a.s.l. in the southern Alps (200 to 400 m above the reference period);

• 2250 and 2300 m a.s.l. in the French Pyrenees (300 to 350 m above the reference period);

• 2550 and 2650 m a.s.l. in the Spanish and Andorran Pyrenees (300 to 350 m above the reference period).

Due to the combined effect of decreasing natural snow conditions and decreasing suitable conditions for snowmaking, the 90 % snow reliability line using snowmaking is projected to rise by 200 to 300 m in the northern Alps, 300 m in the southern Alps and up to 400 to 600 m in the Pyrenees compared to the reference period. In the near future the median elevation of the snowmaking reliability is projected to range between 950 and 1050 m a.s.l. in the northern Alps, 1350 m a.s.l. in the southern Alps, 1700 to 1850 m a.s.l. in the French Pyrenees, and 1750 to 1900 m a.s.l. in the Spanish and Andorran Pyrenees. The production of machine-made snow at the snow reliability line is projected to remain steady or to decrease in the Pyrenees, up to 15 % compared to the reference period. In the Alps, the production of machine-made snow is projected to increase for all scenarios up to 15 %. This highlights the higher suitability of climate conditions for snowmaking in the Alps compared to the Pyrenees and increases the gap in the elevation of the snowmaking reliability between these areas (Fig. 3).

### 3.1.3 Future change at the end of the century (2080–2100)

The impact of climate change on natural snow conditions beyond the mid-century is projected to be highly dependent on the climate scenario. Conditions at the end of the century (2080–2100) are projected to remain similar to those in the near future, for RCP2.6, with the median 70 % groomed snow reliability line ranging between 200 and 300 m above the elevation for the reference period. Under RCP8.5, this elevation at the end of the century would be 850 m higher than the value for the reference period in the northern and southern Alps, 900 m in the Spanish and Andorran Pyrenees, and up to 1050 m in the French Pyrenees.

The snowmaking reliability elevation is projected to suffer from the decrease in periods suitable for snowmaking. The median elevation at the end of the century is projected to be 200 m (northern Alps) to 450 m (French Pyrenees) higher than the value for the reference period for RCP2.6 and up to 1100 m (northern and southern Alps) to 1450 m (French Pyrenees) higher for RCP8.5. The median elevation of the reliability using snowmaking for RCP8.5 is projected to range at the end of the century between 1850 m a.s.l. in the northern Alps, 2150 m a.s.l. in the southern Alps, and 2700 m a.s.l. in the French, Spanish and Andorran Pyrenees (Fig. 3).

In the Pyrenees, the production of machine-made snow is projected to decrease by 15 % to 35 % in the French Pyrenees and 10 % to 20 % in the Spanish and Andorran Pyrenees (10 %–90 % percentiles) compared to the reference period due to the lack of suitable conditions. In the Alps, snowmaking is projected to remain relatively steady at the snow reliability elevation compared to the near future with higher requirements compared to the reference period up to 10 %.

## 3.2 Snow reliability of ski resorts

### 3.2.1 Past climate conditions

Figures 3 and 4 and Table 2 show a deviation between the SAFRAN reference dataset and results derived from climate models (HIST) for the reference period. We therefore focus our analysis on the comparison of snow conditions computed by climate models for the reference and future periods. Based on climate models, ski lift infrastructures were reliable during the reference period (1986–2005), with either natural snow conditions (50 % in categories 1, 2 and 3 altogether; Table 2, HIST) or with snowmaking (49 % in categories 4 and 5 altogether). Natural snow conditions in larger ski resorts were more reliable than in the smaller ones with 44 resorts representing 50 % of the ski lift power being reliable using natural snow only and 129 ski resorts also representing 49 % of the ski lift power being only reliable with snowmaking (Table 2). Categories 6 and 7 include resorts where 90 % snowmaking reliability can not be achieved at the elevation of the village (category 6) or at the mean ski lift elevation (category 7). These categories represent the situation of a marginal fraction of ski resorts in the reference period: less than 1 % unreliable facilities (two resorts in these categories) and might therefore be considered in a critical situation in terms of snow conditions. Figures 4 and 5 also illustrate a significant geographical pattern with most naturally snow-reliable ski resorts being located in the northern Alps and central Pyrenees. This can be related to the lower elevation of the snow reliability line in the northern Alps compared to the southern Alps or the Pyrenees (Fig. 3, Appendix A1) and the higher elevation of larger ski resorts, most of them being located in the northern Alps and central Pyrenees (Figs. 2 and 5). The variability is particularly high between the northern Alps (a majority of ski resorts were naturally snow reliable: 67 % of ski lift power) and the southern Alps (89 % were reliable with snowmaking) highlighting a higher dependence of southern Alps ski resorts on snowmaking in the reference period (only 12 % of ski lift power was naturally snow reliable). The situation of the Pyrenees ski resorts lies in between (Fig. 5).

Figure 4Distribution of the total ski lift power (%) within reliability categories for distinct periods and scenarios (Table 2).

Table 2Distribution of the total ski lift power (%) within reliability categories for distinct periods and scenarios (Fig. 4) with the corresponding number of ski resorts (n).

Figure 5Fraction of ski lift power (%) for a given category (Sect. 2.3). Categories 1, 2 and 3 illustrate ski resorts where natural snow conditions are reliable. Categories 4 and 5 illustrate ski resorts where snow conditions are reliable with snowmaking. Categories 6 and 7 illustrate ski resorts where snowmaking is no longer efficient in reducing the effect of natural snow scarcity at the lowest elevations of the resort.

### 3.2.2 Future change in the near future (2030–2050)

In the near future (2030–2050) and depending on the RCP, only 14 to 24 ski resorts (21 % to 32 % of ski lift power) are projected to remain snow reliable based on natural conditions, all being located in the northern Alps except one in the central Pyrenees (Table 2). An additional 83 to 116 resorts (representing 49 % to 64 % of ski lift power) are projected to remain reliable with snowmaking. Overall, a majority of ski resorts would remain reliable, either with snowmaking or under natural snow conditions (75 % to 86 % of ski lift power). A significant fraction of 45 to 75 ski resorts (14 % to 25 % of ski lift power) would however turn either into category 6 (12 % to 18 % of ski lift power) or even category 7 (2 % to 7 % of ski lift power) where 90 % snowmaking reliability can not be achieved at the elevation of the village (category 6) or at the mean ski lift elevation (category 7). The geographical pattern identified for past climate conditions is projected to remain in the near future. Even though there would not be any naturally snow-reliable ski resort in the southern Alps, snow conditions are projected to remain reliable with snowmaking for most resorts (reduction from 100 % to 89 % of reliable with snowmaking ski lift power), displaying a consistent distribution between reliability categories compared to the reference period (Fig. 4). Conversely, the projected impact on the Pyrenees ski resorts is significant, particularly in the French Pyrenees. There would remain a single naturally snow reliable resort, but more important is the fraction of resorts turning into category 6 (45 % to 58 % of ski lift power in the French Pyrenees and 32 % to 59 % in the Spanish and Andorran Pyrenees) or even category 7 (12 % to 42 % of ski lift power in the French Pyrenees and 7 % to 20 % in the Spanish and Andorran Pyrenees).

### 3.2.3 Future change at the end of the century (2080–2100)

Beyond the near future, the evolution of snow conditions strongly depends on the climate scenario, due to both the evolution of natural snow conditions and the availability of suitable periods for snowmaking (Fig. 3). According to the scenario RCP2.6, snow reliability is projected to remain similar or even improve at the end of the century (2080–2100) compared to the near future (2030–2050). Figures 3 and 4 and Table 2 illustrate the significant impact of climate change on the snow conditions and ski resort reliability for RCP8.5 compared to the two other scenarios. Our projections indicate that there would not remain any ski resort with reliable natural snow conditions based on RCP8.5, with only 24 ski resorts (28 % of ski lift power) benefiting from snowmaking reliability (Table 2), all of them being located in the Alps. Figure 4 illustrates a strong geographical pattern within the Alps with higher snow reliability in eastern central Alps compared to external and southern massifs. End of century RCP8.5 ski resorts reliable with snowmaking are projected to be located in Vanoise (n= 7), Haute-Tarentaise (n= 5), Maurienne (n= 5) and Haute-Maurienne (n= 3) in the northern Alps, and Thabor (n= 1), Pelvoux (n= 1), Queyras (n= 1) and Champsaur (n= 1) in the southern Alps.

4 Discussion

A number of limitations remain in our approach and should be carefully considered in the interpretation of our results. Concerning the modelling of the snowpack evolution under past and future climate conditions, meteorological forcing data are aggregated at the scale of a massif (an approximately 1000 km2 surface area) and by elevation steps of 300 m, which is a significant improvement compared to previous investigations although local effects are still neglected. The snow melting rate is probably underestimated in the model leading to somewhat optimistic results . The main reason for this is the one-dimensional assumption in the snowpack model neglecting the snow–ground partitioning, particularly when the natural snow melts out and leaves the ski slope as an isolated snow patch in grass or rock fields . This situation is likely to be more frequent under future climate conditions, resulting in increasingly optimistic results compared to the reference period. Additionally, all results computed based on the observational reference dataset and climate models exhibit differences in the reference period (Figs. 3 and 4 and Table 2). Discrepancies may be due to potential biases of the multivariate distribution of the meteorological variables produced by the adjustment and downscaling method . This could result in potential nonlinear effects due to multiple dependencies, especially on temperature, relative humidity, precipitation and wind speed.

Beyond the modelling of snow conditions, the main limitations pertain to the snow reliability line approach. Single-point representations are considered on a flat field, i.e. neglecting the aspect and slope angles of a given ski area, which is of high importance in the seasonal evolution of the snowpack and might highly differ from one resort to another. These representations also neglect that all slopes are not covered by snowmaking facilities, hampering any detailed investigation of the evolution of water requirements (results are limited to values per unit surface area). Modelling chains including spatial representations of ski resorts may overcome such weaknesses of the snow reliability line approach . Additionally, even though snowmaking may appear to be an efficient method to technically reduce the impacts of natural snow scarcity, the attractiveness of a given resort may be damaged due to either the lack of snow in parts of the ski resort not equipped with facilities or even the lack of natural snow (landscapes, winter spirit).

We provided information beyond a binary assessment of reliable or unreliable by creating reliability categories, although economic implications should be specifically investigated with a more detailed approach. For example, the relative economic importance of specific periods (Christmas and Winter school holidays) is also neglected in this approach, similarly to previous uses of the snow reliability line. More importantly, our study highlights ski resorts which, under present climate conditions, exhibit challenging snow reliability indicators (category 5 ski resorts in outer northern Alps regions, the southernmost southern Alps, and the eastern and western parts of the Pyrenees), but are currently operational. This indicates that snow reliability is only one factor of the socio-economic performance of ski resorts. This corroborates that the assessment of the sustainability of winter ski tourism destinations must encompass dimensions other than snow conditions alone, consistent with earlier findings .

5 Conclusion

State-of-the-art snowpack modelling and climate projections were used in the present investigation to provide a snow reliability assessment of a large sample of 175 ski resorts in the French Alps and Pyrenees (France, Andorra and Spain) under past and future climate conditions. We report on a significant spatial variability in snow reliability, with or without snowmaking. The northern Alps showed the best natural snow conditions for the reference period (1986–2005) and under future climate conditions. Snowmaking appears to be an efficient method to improve the snow reliability with 99 % of ski lift facilities reliable with snowmaking for the reference period. This is particularly true in the southern Alps where snowmaking leads to a lower elevation of the snowmaking reliability compared to the Pyrenees, while the natural snow reliability line is higher. This situation is projected to remain in future climate conditions and snow reliability elevation is projected to significantly rise due to the decrease in natural snow conditions and the suitable conditions for snowmaking. The difference between projected deviation between climate scenarios is very low in the near future (2030–2050). Depending on the RCP, 21 % to 32 % of ski lift infrastructures would remain reliable based on natural snow conditions while another 14 % to 25 % might be considered in a critical situation, i.e. for which snowmaking reliability can not be achieved. Significant snowmaking requirements are projected to be necessary at the snow reliability line ranging between 400 and 700 kg m−2, i.e. an equivalent 80 to 140 cm of machine-made snow production. Deviations between climate scenarios only appear after the mid-century with limited changes compared to the near future (RCP2.6) or continuous decrease in the snow reliability (RCP8.5). At the end of the century and for RCP8.5, our projections indicate that there would not remain any reliable resort based on natural snow conditions and only 24 resorts (28 % of ski lift facilities) would benefit from snowmaking reliability, all being located in the Alps.

The past and future snow reliability of ski resorts in the French Alps and Pyrenees is highly variable, and the present investigation illustrates the relevance of considering local situations rather than drawing general conclusions. We believe that our results might be substantial material for discussions of the relevance of snowmaking as a technical adaptation and the decision-making regarding investments in these facilities. Management implications and economic issues might also be derived from this approach, which should be extended to mid-elevation areas in France (Jura, Vosges, Massif Central). This also bears potential for wider extension including at the European scale taking advantage of the fact that the method does not require complex data to characterize ski resorts (village and mean ski lift elevation) and could be applied to simulations of natural and managed snow at the European scale .

Assessing the impact of climate change on the ski tourism economy requires not only an estimate of future changes of natural and managed snow conditions, which we provided here, but also additional information on water requirements for snowmaking and how it affects the environmental context and business model of the ski industry. While until the mid-21st century snowmaking appears to be an efficient adaptation option to reduce the climate change hazard to ski resort operating conditions, their environmental footprint and socio-economic functioning may be altered, thereby increasing the vulnerability dimension of the socio-ecological risk if the mountain tourism business model remains unchanged. Towards the end of the century, under high-emission climate change scenarios (RCP8.5), the snow reliability will severely be questioned for most ski resorts currently operating in the Pyrenees and the French Alps, with and without snowmaking, with increased climate change risk under the current mountain tourism business model. Regardless of the time period of interest, future studies are required to analyse and assess all dimensions of climate change impacts and risk to this key mountain economic sector, in France and many other places on Earth.

Data availability
Data availability.

The climate projections used in this study are available on the Drias portal: http://drias-climat.fr/ (last access: 24 April 2019). Model results are available upon request to the corresponding author.

Appendix A: Snow reliability elevation

## A1 Reference period (1986–2005)

Table A1Snow reliability elevation for the reference period (1986–2005) for the 42 massifs distributed over the northern Alps, southern Alps, French Pyrenees, and Spanish and Andorran Pyrenees, computed by the reference dataset (SAFRAN) and the climate models (HIST) for three distinct reliability requirements (50 %, 70 % and 90 %).

## A2 Near future (2030–2050)

Table A2Snow reliability elevation for the near future (2030–2050) for the 42 massifs distributed over the northern Alps, southern Alps, French Pyrenees, and Spanish and Andorran Pyrenees, computed by climate models for the RCP2.6, RCP4.5 and RCP8.5 and for three reliability requirements (50 %, 70 % and 90 %).

## A3 End of the century (2080–2100)

Table A3Snow reliability elevation for the end of the century (2080–2100) for the 42 massifs distributed over the northern Alps, southern Alps, French Pyrenees, and Spanish and Andorran Pyrenees, computed by climate models for the RCP2.6, RCP4.5 and RCP8.5 and for three reliability requirements (50 %, 70 % and 90 %).

Appendix B: Detailed features of individual ski resorts

Table B1Main features of the 175 ski resorts included in the present work grouped by massifs and major areas (northern and southern Alps, French and Spanish and Andorran Pyrenees).

Author contributions
Author contributions.

SM and EG designed the research; PS developed the model, carried out the experiments, and produced the data and most figures with support from co-authors; HF contributed to produce the data and mapping figures; DV produced the adjusted climate projections; ML contributed to the development of the model and production of climate forcing data; MV contributed to the production of the reanalysis forcing data; MP provided input for the Spanish and Andorran resorts; all authors contributed to the analysis and interpretation of the results; PS wrote the paper, based on input and feedback from all co-authors.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Acknowledgements
Acknowledgements.

This study was funded by Région Rhône-Alpes (PhD grant of Pierre Spandre), and benefited from funding from the French Ministry for Ecology (MTES) to the ADAMONT project through GICC and ONERC, from the Interreg project POCTEFA/Clim'Py (OPCC2), from the IDEX Univ. Grenoble Alpes Cross Disciplinary Project “Trajectories”, and from the Conseil Départemental de l'Isère – Isère Tourisme. The CGET – Comité de Massif Alpes funded the creation of the BD Stations database. CNRM/CEN and Irstea are part of LabEX OSUG@2020 (ANR10 LABX56). Deborah Verfaillie's work has been funded by the European project EUCP (H2020-SC5-2016-776613). We acknowledge useful comments and suggestions from O. Cenk Demiroglu, Robert Steiger, one anonymous reviewer and the editor Jürg Schweizer.

Review statement
Review statement.

This paper was edited by Jürg Schweizer and reviewed by Cenk Demiroglu, Robert Steiger, and one anonymous referee.

References

Abegg, B., Agrawala, S., Crick, F., and de Montfalcon, A.: Climate change impacts and adaptation in winter tourism, in: Climate Change in the European Alps, edited by: Agrawala, S., OECD Paris, 25–60, https://doi.org/10.1787/9789264031692-en, 2007. a, b, c, d, e, f, g, h, i, j

Beniston, M., Farinotti, D., Stoffel, M., Andreassen, L. M., Coppola, E., Eckert, N., Fantini, A., Giacona, F., Hauck, C., Huss, M., Huwald, H., Lehning, M., López-Moreno, J.-I., Magnusson, J., Marty, C., Morán-Tejéda, E., Morin, S., Naaim, M., Provenzale, A., Rabatel, A., Six, D., Stötter, J., Strasser, U., Terzago, S., and Vincent, C.: The European mountain cryosphere: a review of its current state, trends, and future challenges, The Cryosphere, 12, 759–794, https://doi.org/10.5194/tc-12-759-2018, 2018. a

Breiling, M. and Charamza, P.: The impact of global warming on winter tourism and skiing: a regionalised model for Austrian snow conditions, Reg. Environ. Change, 1, 4–14, https://doi.org/10.1007/s101130050003, 1999. a, b

Brun, E., David, P., Sudul, M., and Brunot, G.: A numerical model to simulate snow-cover stratigraphy for operational avalanche forecasting, J. Glaciol., 38, 13–22, 1992. a

Damm, A., Greuell, W., Landgren, O., and Prettenthaler, F.: Impacts of +2C global warming on winter tourism demand in Europe, Climate Services, 7, 31–46, https://doi.org/10.1016/j.cliser.2016.07.003, 2017. a, b

Dawson, J. and Scott, D.: Managing for climate change in the alpine ski sector, Tourism Manage., 35, 244–254, https://doi.org/10.1016/j.tourman.2012.07.009, 2013. a, b, c

Durand, Y., Brun, E., Mérindol, L., Guyomarc'h, G., Lesaffre, B., and Martin, E.: A meteorological estimation of relevant parameters for snow models, Ann. Glaciol., 18, 65–71, 1993. a

Durand, Y., Giraud, G., Laternser, M., Etchevers, P., Mérindol, L., and Lesaffre, B.: Reanalysis of 47 Years of Climate in the French Alps (1958–2005): Climatology and Trends for Snow Cover, J. Appl. Meteorol. Clim., 48, 2487–2512, https://doi.org/10.1175/2009JAMC1810.1, 2009a. a

Durand, Y., Giraud, G., Laternser, M., Etchevers, P., Mérindol, L., and Lesaffre, B.: Reanalysis of 44 Yr of Climate in the French Alps (1958–2002): Methodology, Model Validation, Climatology, and Trends for Air Temperature and Precipitation., J. Appl. Meteorol. Clim., 48, 429–449, https://doi.org/10.1175/2008JAMC1808.1, 2009b. a, b

Durand, Y., Giraud, G., Goetz, D., Maris, M., and Payen, V.: Modeled Snow Cover in Pyrenees Mountains and Cross-Comparisons Between Remote-Sensed and Land-Based Observation Data, in: Proceedings of the International Snow Science Workshop, 16–21 September 2012, Anchorage, Alaska, USA, vol. 9981004, 2012. a

Elsasser, H. and Bürki, R.: Climate change as a threat to tourism in the Alps, Clim. Res., 20, 253–257, https://doi.org/10.3354/cr020253, 2002. a, b, c, d

Falk, M. and Vanat, L.: Gains from investments in snowmaking facilities, Ecol. Econ., 130, 339–349, https://doi.org/10.1016/j.ecolecon.2016.08.003, 2016. a, b

François, H., Morin, S., Lafaysse, M., and 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, https://doi.org/10.1016/j.coldregions.2014.08.005, 2014. a, b, c, d, e

Gilaberte-Burdalo, M., Lopez-Martin, F., M. R. Pino-Otin, M., and Lopez-Moreno, J.: Impacts of climate change on ski industry, Environ. Sci. Pol., 44, 51–61, https://doi.org/10.1016/j.envsci.2014.07.003, 2014. a

Gilaberte-Búrdalo, M., López-Moreno, J., Morán-Tejeda, E., Jerez, S., Alonso-González, E., López-Martín, F., and Pino-Otín, M.: Assessment of ski condition reliability in the Spanish and Andorran Pyrenees for the second half of the 20th century, Appl. Geogr., 79, 127–142, https://doi.org/10.1016/j.apgeog.2016.12.013, 2017. a, b, c, d

Hanzer, F., Marke, T., and Strasser, U.: Distributed, explicit modeling of technical snow production for a ski area in the Schladming region (Austrian Alps), Cold Reg. Sci. Technol., 108, 113–124, https://doi.org/10.1016/j.coldregions.2014.08.003, 2014. a, b, c

Hendrikx, J. and Hreinsson, E.: The potential impact of climate change on seasonal snow in New Zealand: industry vulnerability and future snowmaking potential, Theor. Appl. Climatol., 110, 619–630, https://doi.org/10.1007/s00704-012-0713-z, 2012. a, b

Hennessy, K., Whetton, P., Walsh, K., Smith, I., Bathols, J., Hutchinson, M., and Sharples, J.: Climate change effects on snow conditions in mainland Australia and adaptation at ski resorts through snowmaking, Clim. Res., 35, 255, https://doi.org/10.3354/cr00706, 2007. a, b

Hoegh-Guldberg, O., Jacob, D., Taylor, M., Bindi, M., Brown, S., Camilloni, I., Diedhiou, A., Djalante, R., Ebi, K., Engelbrecht, F., Guiot, J., Hijioka, Y., Mehrotra, S., Payne, A., Seneviratne, S., Thomas, A., Warren, R., and Zhou, G.: Impacts of 1.5 C Global Warming on Natural and Human Systems, in: Global Warming of 1.5 C. An IPCC Special Report on the impacts of global warming of 1.5 C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty, edited by: Masson-Delmotte, V., Zhai, P., Portner, H.-O., Roberts, D., Skea, J., Shukla, P. R., Pirani, A., Moufouma-Okia, W., Pean, C., Pidcock, R., Connors, S., Matthews, J. B. R., Chen, Y., Zhou, X., Gomis, M. I., Lonnoy, E., Maycock, T., Tignor, M., and Waterfield, T., in press, 2019. a

Jacob, D., Petersen, J., Eggert, B., Alias, A., Christensen, O. B., Bouwer, L. M., Braun, A., Colette, A., Deque, M., Georgievski, G., Georgopoulou, E., Gobiets, A., Menut, L., Nikulin, G., Haensler, A., Hempelmann, N., Jones, C., Keuler, K., Kovats, S., Kroner, N., Kotlarski, S., Kriegsmann, A., Martin, E., Meijgaard, E. V., Moseley, C., Pfeifer, S., Preuschmann, S., Radermacher, C., Radtke, K., Rechid, D., Rounsevell, M., Samuelsson, P., Somot, S., Soussana, J.-F., Teichmann, C., Valentini, R., Vautard, R., Weber, B., and Yiou, P.: EURO-CORDEX: new high-resolution climate change projections for European impact research, Reg. Environ. Change, 14, 563–578, https://doi.org/10.1007/s10113-013-0499-2, 2014. a

Koenig, U. and Abegg, B.: Impacts of climate change on winter tourism in the Swiss Alps, J. Sustain. Tour., 5, 46–58, https://doi.org/10.1080/09669589708667275, 1997. a, b, c, d

Kotlarski, S., Keuler, K., Christensen, O. B., Colette, A., Déqué, M., Gobiet, A., Goergen, K., Jacob, D., Lüthi, D., van Meijgaard, E., Nikulin, G., Schär, C., Teichmann, C., Vautard, R., Warrach-Sagi, K., and Wulfmeyer, V.: Regional climate modeling on European scales: a joint standard evaluation of the EURO-CORDEX RCM ensemble, Geosci. Model Dev., 7, 1297–1333, https://doi.org/10.5194/gmd-7-1297-2014, 2014. a

Laternser, M. and Schneebeli, M.: Long-term snow climate trends of the Swiss Alps (1931–99), Int. J. Climatol., 23, 733–750, https://doi.org/10.1002/joc.912, 2003. a

Luthe, T., Wyss, R., and Schuckert, M.: Network governance and regional resilience to climate change: empirical evidence from mountain tourism communities in the Swiss Gotthard region, Reg. Environ. Change, 12, 839–854, https://doi.org/10.1007/s10113-012-0294-5, 2012. a

Maris, M., Giraud, G., Durand, Y., Navarre, J.-P., and Mérindol, L.: Results of 50 years of climate reanalysis in the French Pyrenees (1958–2008) using the SAFRAN and CROCUS models, in: Proceedings of the International Snow Science Workshop, 27 September–2 October 2009, Davos, Switzerland, 219–223, 2009. a

Marke, T., Strasser, U., Hanzer, F., Stötter, J., Wilcke, R., and Gobiet, A.: Scenarios of future snow conditions in Styria (Austrian Alps), J. Hydrometeorol., 16, 261–277, https://doi.org/10.1175/JHM-D-14-0035.1, 2014. a, b, c

Matulla, C., Auer, I., Böhm, R., Ungersböck, M., Schöner, W., Wagner, S., and Zorita, E.: Outstanding past decadal-scale climate events in the Greater Alpine Region analysed by 250 years data and model runs, GKSS-Forschungszentrum, Geesthacht, Germany, 2005. a

Morin, S., Abegg, B., Demiroglu, O. C., Pons, M., Weber, F., Amacher, A., François, H., George, E., Soubeyroux, J. M., Samacoits, R., Lafaysse, M., Franklin, S., Clifford, D., Cauchy, A., and Dubois, G.: The “Mountain Tourism” component of the Copernicus Climate Change Services – Sectoral Information Service “European Tourism”: towards pan-European analysis and projections of natural and managed snow conditions, in: Proceedings of the International Snow Science Workshop, 7–12 October 2018, Innsbruck, Austria, 542–547, 2018. a

Moss, R., Edmonds, J., Hibbard, K., Manning, M., Rose, S., Van Vuuren, D., Carter, T., Emori, S., Kainuma, M., Kram, T., Meehl, G. A., Mitchell, J. F. B., Nakicenovic, N., Riahi, K., Smith, S. J., Stouffer, R. J., Thomson, A. M., Weyant, J. P., and Wilbanks, T. J.: The next generation of scenarios for climate change research and assessment, Nature, 463, 747–756, https://doi.org/10.1038/nature08823, 2010. a

Mott, R., Daniels, M., and Lehning, M.: Atmospheric Flow Development and Associated Changes in Turbulent Sensible Heat Flux over a Patchy Mountain Snow Cover, J. Hydrometeorol., 16, 1315–1340, https://doi.org/10.1175/JHM-D-14-0036.1, 2015. a

Pons, M., López-Moreno, J., Rosas-Casals, M., and Jover, E.: The vulnerability of Pyrenean ski resorts to climate-induced changes in the snowpack, Climatic Change, 131, 591–605, https://doi.org/10.1007/s10584-015-1400-8, 2015. a, b, c, d, e, f, g, h, i

Pons-Pons, M., Johnson, P. A., Rosas Casals, M., Sureda Carbonell, B., and Jover Comas, E.: Modeling climate change effects on winter ski tourism in Andorra, Clim. Res., 54, 197–207, https://doi.org/10.3354/cr01117, 2012. a, b, c

Rixen, C., Teich, M., Lardelli, C., Gallati, D., Pohl, M., Pütz, M., and Bebi, P.: Winter tourism and climate change in the Alps: an assessment of resource consumption, snow reliability, and future snowmaking potential, Mt. Res. Dev., 31, 229–236, https://doi.org/10.1659/MRD-JOURNAL-D-10-00112.1, 2011. a, b

Schmidt, P., Steiger, R., and Matzarakis, A.: Artificial snowmaking possibilities and climate change based on regional climate modeling in the Southern Black Forest, Meteorol. Z., 21, 167–172, https://doi.org/10.1127/0941-2948/2012/0281, 2012. a, b

Scott, D., McBoyle, G., and Mills, B.: Climate change and the skiing industry in Southern Ontario (Canada): exploring the importance of snowmaking as a technical adaptation, Clim. Res., 23, 171–181, https://doi.org/10.3354/cr023171, 2003. a, b, c, d, e

Scott, D., McBoyle, G., Minogue, A., and Mills, B.: Climate change and the sustainability of ski-based tourism in eastern North America: A reassessment, J. Sustain. Tour., 14, 376–398, https://doi.org/10.2167/jost550.0, 2006. a, b

Spandre, P., François, H., Morin, S., and George-Marcelpoil, E.: Snowmaking in the French Alps. Climatic context, existing facilities and outlook, Rev. Geogr. Alp., 103-2, https://doi.org/10.4000/rga.2913, 2015. a

Spandre, P., François, H., George-Marcelpoil, E., and Morin, S.: Panel based assessment of snow management operations in French ski resorts, Journal of Outdoor Recreation and Tourism, 16, 24–36, https://doi.org/10.1016/j.jort.2016.09.002, 2016a. a, b, c

Spandre, P., Morin, S., Lafaysse, M., George-Marcelpoil, E., François, H., and Lejeune, Y.: Integration of snow management in a detailed snowpack model, Cold Reg. Sci. Technol., 125, 48–64, https://doi.org/10.1016/j.coldregions.2016.01.002, 2016b. a, b, c, d, e

Spandre, P., François, Verfaillie, D., Lafaysse, M., Déqué, M., Eckert, N., George, E., and Morin, S.: Climate change impacts on the snow reliability of French Alps ski resorts, in: Proceedings of the International Snow Science Workshop, 7–12 October 2018, Innsbruck, Austria, 512–516, 2018. a

Steiger, R.: The impact of climate change on ski season length and snowmaking requirements in Tyrol, Austria, Clim. Res., 43, 251–262, https://doi.org/10.3354/cr00941, 2010. a, b, c, d, e, f, g

Steiger, R. and Mayer, M.: Snowmaking and climate change: Future options for snow production in Tyrolean ski resorts, Mt. Res. Dev., 28, 292–298, https://doi.org/10.1659/mrd.0978, 2008. a, b, c

Steiger, R., Scott, D., Abegg, B., Pons, M., and Aall, C.: A critical review of climate change risk for ski tourism, Curr. Issues Tour., 1–37, https://doi.org/10.1080/13683500.2017.1410110, 2017. a, b, c, d, e, f

Taylor, K., Stouffer, R., and Meehl, G.: An overview of CMIP5 and the experiment design, B. Am. Meteorol. Soc., 93, 485–498, https://doi.org/10.1175/BAMS-D-11-00094.1, 2012.  a

Verfaillie, D., Déqué, M., Morin, S., and Lafaysse, M.: The method ADAMONT v1.0 for statistical adjustment of climate projections applicable to energy balance land surface models, Geosci. Model Dev., 10, 4257–4283, https://doi.org/10.5194/gmd-10-4257-2017, 2017. a, b, c, d

Verfaillie, D., Lafaysse, M., Déqué, M., Eckert, N., Lejeune, Y., and Morin, S.: Multi-component ensembles of future meteorological and natural snow conditions for 1500 m altitude in the Chartreuse mountain range, Northern French Alps, The Cryosphere, 12, 1249–1271, https://doi.org/10.5194/tc-12-1249-2018, 2018. a, b, c

Vionnet, V., Brun, E., Morin, S., Boone, A., Faroux, S., Le Moigne, P., Martin, E., and Willemet, J.-M.: The detailed snowpack scheme Crocus and its implementation in SURFEX v7.2, Geosci. Model Dev., 5, 773–791, https://doi.org/10.5194/gmd-5-773-2012, 2012. a

Wielke, L. M., Haimberger, L., and Hantel, M.: Snow cover duration in Switzerland compared to Austria, Meteorol. Z., 13, 13–17, https://doi.org/10.1127/0941-2948/2004/0013-0013, 2004. a