Determining the drivers for snow gliding

Snow gliding is a key factor for snow-glide avalanche formation and soil erosion. This study considers atmospheric and snow variables, vegetation characteristics, and soil properties and determines their relevance for snow gliding at a test site (Wildkogel, Upper Pinzgau, Austria) during winter 2014/2015. The time-dependent data were collected at a high temporal resolution. In addition to conventional sensors, a “snow melt analyzer” was used. The analysis shows that the soil temperature 10 cm below the surface, the phytomass of mosses, the liquid water content in the snowpack, and the static friction coefficient of the glide shoes had significant influence on snow gliding during the whole winter. In the first period (October to January) the soil moisture at the surface and 1.5 cm below the surface and the length of the slope uphill of the glide shoes affected the snow gliding, too. In the second period (February to May) the soil temperature at the surface, the soil moisture 10 cm below the surface, and the slope angle had additional influence on snow gliding. The role of the vegetation in the snow-glide process is determined by the influence on the static friction coefficient caused by its composition and characteristics and by mossrich and short-stemmed canopies being seemingly more interconnected with the snowpack. In addition to the soil and snow properties, the topography and the vegetation characteristics, further investigations may be focused on the freezing and melting processes in the uppermost soil layers and at the soil surface.


Introduction
Deposited snow on the ground is in motion caused by gravity, external forces, or metamorphism.The movement inside the snowpack is called creeping, and the sliding of the entire snowpack on an inclined ground surface is referred as snow gliding (In der Gand and Zupancic, 1966).Snow gliding is favored by a smooth ground surface and a lowermost layer of wet snow (In der Gand and Zupancic, 1966).Once the glide motion turns into an avalanche movement, the process is called a glide avalanche (UNESCO, 1981).
The presence of liquid water at the bottom of the snowpack is a basic requirement for snow gliding (In der Gand, 1954;Lackinger, 1988;McClung et al., 1994;Mitterer and Schweizer, 2013).Several sources exist to provide liquid water to this location (Ceaglio et al., 2012;Ceaglio et al., 2017;Mitterer and Schweizer, 2012).Rain on the snow surface, as well as melting snow near to the surface (Koh and Jordan, 1995), can percolate the isothermal snowpack.Geothermal heat flux can provide energy to melt snow at the bottom of the snowpack (McClung and Clarke, 1987).The suction head can lift water (Mitterer and Schweizer, 2012;Ceaglio et al., 2017) which is produced by melting ice stored in the soil or it can be advected through channels in the soil (ground water outflow).
In addition to the presence of liquid water at the bottom of the snowpack, further variables influence the intensity of snow gliding.Therefore, air temperature can be used to classify the glide snow avalanches into warm-temperature events and coldtemperature events (Clarke and McClung, 1999).The viscosity of snow depends on the snow temperature (Loth et al., 1993;Morris, 1994) and snow water content (Mitterer and Schweizer, 2012;McClung and Clarke, 1987).The slope angle, the micro relief, and the hydrological properties of the slope influence the glide velocity (Ceaglio et al., 2017;McClung and Schaerer, 1999;Margreth, 2007).Friction originated by the vegetation depends on its composition and height (Höller et al., 2009).Both the vegetation and the micro relief depend on the land use, which is an input for snow glide modeling (Leitinger et al., 2008;Maggioni et al., 2016).Bartelt et al. (2012) developed a two dimensional visco-elastic continuum model to specify the start of gliding snow avalanches.Ancey and Bain (2015) summarized the knowledge concerning the formation of snow glide avalanches and its impact on obstacles in the path.They concluded that meteorological conditions and topographic features causing snow gliding are well known, but the mechanisms are poorly understood.However, the role of the vegetation and the soil conditions have not been considered in very much detail in previous studies which is also indicated by (Höller, 2014) who stated that the conditions at the snowsoil interface have to be investigated most notably.Even though Feistl et al. (2014) documented the vegetation cover, vegetation height and the terrain properties of 101 glide-snow avalanche release areas on the Dorfberg and indicated four characteristic types of vegetation.Leitinger et al. (2008) established a measure for vegetation roughness (i.e.surface roughness) and showed that this factor has a significant influence on snow-glide distances.However, detailed consideration of the soil-vegetation system in the snow-glide process is missing.This was also noted by Höller (2014) who stated that the conditions at the snow-soil interface have to be investigated most notably.
Although the impact of global change on land cover was mainly due to socio-economic drivers (Tasser et al., 2017), future impact of changing climate will accelerate changes in vegetation composition and vegetation roughness.Hence, studies on causal links and quantitative impacts are especially crucial for snow-gliding and related processes.Besides measures to simplify the complex interactions of vegetation roughness at the snow-ground interface (i.e.surface roughness, Leitinger et al., 2008), the influence of vegetation composition and liquid water to the interlocking of the soil-vegetation-snow continuum is widely unknown.
This study specifically addresses the role of the soil-vegetation system on snow gliding, with an elaborate experimental setup.
The focus was on the presence of liquid water in the snowpack, on the vegetation, the soil surface, and in the upmost soil layers, as well as vegetation composition and its consequence on snow gliding.Therefore, these key questions are addressed:  Which variables in the soil-vegetation system, the snowpack, and the lowest atmospheric boundary layer have considerable influence on snow gliding? Is it appropriate to distinguish between processes at the beginning of the winter (development of the snowpack) and the late winter (decline of the snowpack)?
 How does vegetation composition influence the snow gliding process? 2 Experimental test site and methods

Test site
The study site is located on the orographic left, south-facing slope of the upper Pinzgau Valley.From the geological point of view, it is a very homogenous area made up mainly of paragneiss and mica schist.This siliceous bedrock is responsible for the presence of cambisols on the pastures.The abandoned and unused areas are mostly based on cambic podzols.The climate at the Wildkogel can be characterized as a subalpine European climate.Long-term average annual rainfall (at 1973 m a.s.l., Schmittenhöhe) amounts to 1501 mm, with the highest monthly precipitations falling in June and August (175-200 mm per month).Long-term average annual temperature is 1.9°C, with the highest monthly average in August at around 10°C.These low temperatures, high precipitation, and the long period of snow cover impose limits on the vegetation period.The investigated slope faces SSE, with slope angles from 20° to 37°.
The area is characterized by pastures and abandoned areas in the immediate vicinity (Baumgärtner, 2016).This situation allowed a comparative approach to be used (Fig. 1).

For
The pasture is stocked with cattle between the end of June and the beginning of September.This area is dominated by grasses and has been classified as Sieversio montanae-Nardetum strictae subassociation typicum (Lüth et al., 2011) with the matgrass (Nardus stricta) as dominant species.Management of the abandoned area ceased about 10 years ago.The predominant species of the area are dwarf shrubs (e.g.Vaccinium myrtillus, V. vitis-idaea, Calluna vulgaris) and the evergreen sedge (Carex sempervirens).

Snow gliding
Snow gliding was measured with glide shoes (In der Gand and Zupancic, 1966).The glide shoes were connected to a drum with a wire.Its displacements generated rotations.A rotary switch generated pulses which were counted by HOBO H6 logger units.The date and time of each pulse was stored.One pulse represents a glide distance of 2.6 mm.A detailed description is given by Leitinger et al. (2008).Forty devices (Fig. 1) were installed at randomly selected places with different land use, topographic conditions, and vegetation characteristics in October 2014 (Baumgärtner, 2016).
The initial force required to displace each shoe was measured with a tension spring balance (Pesola Medio 1000 g).The static friction coefficients for all glide shoes were calculated as the ratio of the initial forces and the normal forces.They represent the influence of different vegetation types and different land uses on snow gliding (Leitinger et al., 2008).

Meteorology and related snow and soil properties
An automatic weather station recorded air temperature, air humidity (Rotronic MP103), snow depth (Sommer UHZ8), snow temperatures (Sommer AD592c; 0, 5, 50, 100 cm), and global radiation (Schenk 8101).It was located at the test site.The data were stored at intervals of 10 minutes by a data logger.
At the meteorological station a snow melt analyzer (SMA, Sommer) was available.It measures the dielectric coefficients with a time-domain reflectometer, using two frequencies along a flat band cable.The different dielectric properties of water and ice are used to determine the volume fractions of the LWC and the ice content (Stähli et al., 2004).The flat band cable was mounted 5 cm above the soil surface.It was aligned parallel to the surface and orientated along the fall line.The acquisitions were recorded by a data logger in 10 minute intervals.Data entries were removed in case that the snow depth was less than 5 cm.
Soil temperatures (Pt-100) and soil moistures (Decagon, ECHO®) were measured at four levels (0, 1.5, 5, 10 cm) in the pastures and the abandoned area.The data were stored at intervals of 5 minutes by a data logger (HOBO® Microstation).

Topographic features and vegetation characteristics
In order to consider the micro-relief close to the snow glide shoes, topographic features were noted at each glide shoe.The slope angle was measured directly by each glide shoe, as well as one meter uphill and one meter downhill.Along the fall line the distances where the micro relief changed were measured (uphill and downhill).The amplitudes (A) and the wavelength (ʎ) of the micro relief were determined.For that purpose, an elastic aluminum pole (length 2 m) was used, which was matched to the ground surface and resulted in a deformation of the slope.With these data, the stagnation depths were calculated according to Salm (1977) for each glide shoe position.
The parameter static friction coefficient was determined to estimate the roughness of the vegetation.For calculation, the weight of the glide shoe and the force needed to move the glide shoe on the vegetation surface was measured (Leitinger et al., 2008).
A vegetation inventory of each snow gliding measurement plot was made by a simplified phytosociological survey, according to Braun-Blanquet (1964).This involves analyzing the degree to which the important plant species are present at the position of the snow-glide shoes.
To determine phytomass pools at the sites, production analyses were carried out at the beginning of the vegetation period (end of May).Within a harvest frame (size 900 cm²), all above-ground stands were harvested destructively.The experiment consisted of 18 and 22 replicate plots for the pasture and the abandoned/agricultural unused area, respectively.
Knowledge of the absolute amounts of the different functional groups are important in order to assess qualitative vegetation composition and the resulting effects on snow gliding (Newesely et al., 2000;Leitinger et al., 2008).Therefore, the harvested phytomass was divided into the following plant functional groups: grass, herbs, dwarf shrubs, lichens, and mosses.The phytomasses were then oven-dried at 80°C until they reached a constant weight, determined as the dry weight.

Data interpretation and statistical methods
In order to identify the magnitude of the influence of the variables, the snow glide rate is defined as the dependent variable.
All other variables are interpreted as independent variables.Since snow gliding in the data set is a binary piece of information for each time step, multiple logistic regression was used to determine the relevant variables (Wilks, 1995).The magnitude of the regression parameters can be used to describe their influence on the dependent variable.
The number of independent variables should be reduced to avoid overfitting.This procedure is often called screening regression and was established by backward elimination (Wilks, 1995).The procedure starts with all potential predictors.At each step the least important predictor is removed until the termination criteria are reached (tolerance of the predictor >0.2 and variance inflation factor <10).
In about 0.5 % of the data entries snow gliding was recorded.The samples with snow gliding were subsequently weighted.
This satisfies that equal amount of 0 and 1 for snow gliding which are used for the multiple logistic regressions (period I: n = 1164096; period II: n = 1340425).A bootstrap is performed by randomly selecting a value, with replacement (i.e. a given value can be represented more than once in the sample).Each sample selected in this manner is used to calculate the regression coefficient B value.This is repeated 100 times, and the generated sample of B values is then used to estimate the standard error and the lower and upper 95% confidence interval.The bootstrapping approach is preferable to that presented by Gude et al.
The logistic regressions fit the parameters B for all variables.The magnitude of exp(B) is used to describe the intensity of its influence on snow gliding.If exp(B)>1 the effect is positive, which means that the probability of snow gliding rises with increasing values for the variable.Values below 1 have a negative effect, and the probability of snow gliding decreases if the values for the variable rises.exp(B)=1 indicates that the corresponding variable has no influence on snow gliding.
Due to the fact that liquid water at the snow-soil interface is a requirement for intense snow gliding (In der Gand and Zupancic, 1966) the measured soil moisture at 0 cm (soil surface) is analyzed in more detail.By using a multiple linear regression model, the regression coefficient was determined to identify the sign and the magnitude of the independent variables.To avoid overfitting, variables which correlate among themselves were excluded.
In order to consider the differences between the properties of a rising and a degrading snowpack, the data set was divided into two sub-periods: period I from October to January, and period II from February to May.For both periods accuracy tables are used to demonstrate how well the applied method is able to distinguishes between the two classes (gliding, no gliding).As score index the hit rate was used which is the fraction of correctly calculated data records and the sum of all data entries (Wilks, 1995).
The Whitney-Mann U-test is a nonparametric rank test (Schönwiese, 2000).It was used to determine the significance levels (p values) for selected variables.
The statistical analyses were accomplished with the software IBM SPSS Statistics (Version 21, IBM SPSS Statistics Software).

Results
The time series in Fig. 2 give an overview of the investigated period.The snow cover season started in October 2014 and ended in late May 2015.It was interrupted twice: in November and in May.In period I the soil temperatures decreased until they reached values between 0°C and 1°C.During period II the soil temperatures were nearly constant until the snow melted.
At the beginning of the winter (period I) snow gliding was recorded by all glide shoes.The LWC reached more than 4 % (volumetric percent).The soil moisture characteristics were different for pastures and abandoned areas.AtThough, at the surface, the soil moisture was close to zero until March in the pastures and in the abandoned areaboth sites (Fig. 2).
At the beginning of period II, the measured LWC values were about 2.5 %.It raised during snow melting, indicated by a rapid decrease in snow height.

Topography and vegetation
An overview of the observations and measurements at the pastures and abandoned areas is given in Tab. 1.
The frequency distributions of vegetation characteristics are L-shaped for all vegetation types (Fig. 3).This indicates that no vegetation type is dominant at the test site.The prevailing slope angle ranges from 25° to 35°.The stagnation depth was below 0.5 m, except in one case, indicating a smooth location of that glide shoe.The friction force was low, and in the majority of the cases very low.The frequency distribution of the canopy heights was between 0.01 m and 0.08 m -higher values were less frequent.The distribution of the slope lengths above and below the glide shoes were equally shaped.The distribution of the slope angles below the glide shoes had a maximum at 30°.

Snow gliding
For period I the soil temperature at 10 cm was determined as the variable with the most influence on snow gliding, followed by the LWC (Tab.2).Moderate influence was detected for soil moisture at 0 cm and soil moisture at 1.5 cm.The soil temperature at 10 cm was the most important variable for period II.A strong negative influence is indicated for the phytomass of mosses and the static friction coefficient in both periods.
The boxplots for the complete data set distinguish between snow gliding and no snow gliding for LWC andperiod I for soil moisture at 0 cm and 1.5 cm (Fig. 4a), soil temperature at 10 cm (Fig. 4b), and the 4c), the LWC (Fig. 4e), the static friction coefficient (Fig. 4g), the phytomass of mosses (Fig. 4c) for period I.4i), and the slope length uphill (Fig. 4k).The positive influence of the soil moisture at 1.5 cm and the soil temperature at 10 cm is obvious, as iswell as the negative effect of the phytomass of mosses and the static friction coefficient.The influence of LWC on snow gliding exists, but it is low.
For period I the hit rate is 85.4 %, and for period II it is 66.0 % (Tab.3).

Soil water content at 0 cm
The presence of liquid water at the bottom of the snowpack is a requirement for snow gliding (In der Gand, 1954;Mitterer and Schweizer, 2013).In order to determine the relevant variables and quantify their influence on snow gliding, a multiple linear regression was calculated for both the pastures and the abandoned area.The soil moisture at 0 cm was used as the dependent variable.The signs of the regression coefficients indicate a positive or a negative relationship (Tab.4).The magnitude represents the intensity of its influence on the soil moisture at 0 cm.For both areas, the soil moisture at 10 cm is identified as the most important variable.Negative correlations were found for soil temperature at 10 cm and snow temperature at 5 cm.
Atmospheric variables had a very low influence on the soil moisture at 0 cm.

Discussion and conclusions
Ceaglio et al. ( 2017) investigated the role of the soil in the context of snow gliding and the formation of glide cracks and avalanches.They concluded that the thermal and hydraulic processes in the soil have to be considered.Our study confirms that the soil moisture at the soil surface, and a few centimeters below the surface, are variables which influence the snow glide rates.Additionally, we found that temperatures in the soil have a significant influence on snow gliding.For example, in the second half of December 2014 a freezing process was observed where the soil moisture abruptly decreased after a cold period with no snow on the ground.Furthermore, the phytomass of mosses affects the snow glide rates at the test site.Clarke and McClung (1999) introduced the terms cold-temperature events and warm-temperature events, which indicate a correlation of glide snow avalanches with air temperatures.Since glide snow avalanches did not occur at the study site, such classification is not useful here.However, to consider different processes during the development of the snowpack and the decline of the snowpack, two sub-periods were defined (period I: October-January; period II: February-May).The soil moisture and the soil temperature had a significant influence on snow gliding in both periods.This indicates a lower viscosity of the moist snowpack and a water transport from the snowpack towards the soil surface.However, the LWC is not the predominant variable that explains the soil moisture at 0 cm (Tab.4).Dreier et al. (2016) investigated the influence of meteorological parameters on snow glide avalanches and divided the winter season into two periods.They found that warm temperature events were mostly associated with a melting snow surface, and cold temperature events are linked with hydraulic process in the basal snow layers and the uppermost soil layers.It confirms the conclusions regarding glide distances presented here.
Some topographical factors also affect snow gliding.In particular, the static friction coefficient has a negative effect on snow gliding.It seems that the friction is reduced by the vegetation, which was depressed by the weight of the snowpack.This depends on the composition and the characteristics of the vegetation (Leitinger et al., 2008).At the test site it can be concluded that dwarf shrubs are more resistant against depression than pastures.
The results also show that the vegetation has a significant effect on snow gliding.Just the phytomass of mosses had a negative influence on snow gliding in both periods.The analyses of the vegetation composition have shown that a higher percentage of mosses exists at low canopy heights (p=-0.52**).Moss-rich and short-stemmed canopies seem to be more interconnected with the snowpack, and thus contribute to a reduction in snow gliding.On the other hand, long-stemmed, grass-rich canopies can be easily felled, and they form an ideal gliding horizon.These findings are in accordance with those of Newesely et al. (2000) showing that the gliding distances are increasing from cut meadows to pastures to uncut or abandoned grasslands.Furthermore, a canopy height is positively correlated with the proportion of dwarf shrub phytomass (p = -0.73***).The predominant dwarf shrub species in the study area are Vaccinium sp. and Rhododendron ferrugineum, and are highly lignified and rigid dwarf shrubs.Such dwarf shrubs, as well as small trees, keep the snowpack back and thus reduce snow gliding (see also Newesely et al., 2000;Leitinger et al., 2008).On the other side, the canopy height is negatively correlated with the phytomass of grasses (p = -0.61***)which promotes snow-gliding (Newesely et al., 2000).
Implications for agricultural land management (Tasser and Tappeiner, 2002) are given as the type of land-use (mowing, grazing), as well as the intensity of land-use (frequency of annual mowing, fertilisation, irrigation, number of grazing animals), lead to characteristic vegetation communities.Mowing and a low level of fertilisation greatly favour the growth of herbs and high growing grasses, while Nardus stricta spreads rapidly on meadows with low land-use intensities (usually mown once a year, not fertilised).After land abandonment, Carex sp.immediately spreads, forming the climax vegetation at the higher altitudes.Below the natural timberline, however, the proliferation of dwarf shrubs and subsequent a natural reforestation are For taking place.Land-cover changes, especially the transitional forms between meadows of high land-use intensity and young forests may have crucial impact for the snow-gliding process (Newesely et al., 2000;Leitinger et al., 2008).If an adequate land-use intensity cannot be maintained, steep areas have to be reforested to shorten a critical time period of high snow-gliding activity.
These investigations on snow gliding confirmed findings from previous studies, and extended them by considering variables describing the vegetation.It seems that the use of soil moisture sensors makes sense for further investigation, which may be focused on the hydraulic processes close to the soil surface.However, upcoming measurement problems of the uppermost partially frozen soil layers must be considered.

Acknowledgements
Snow and weather data including the data of the SMA are part of the operational avalanche warning service of Salzburg (Austria).The authors thank Norbert Altenhofer for providing these data sets, and Michael Butschek for the technical support.

Georg Leitinger and Erich
In particular, they still do not say explicitly how they selected the parameters for the box-plots and the Whitney-Mann U-test (Fig. 4).In Table 2 there are 7 significant parameters for Period I and 8 parameters for Period II, but then the box-plots are made only for some of those parameters.Why?Maybe they did it and the results were not significant… no problem, they should simply state this.
We modified Figure 4. Now boxplots for all variables with significant influence on snow gliding are shown.The rule (the difference of exp(B) from 1) is the same as for table 2. This information is added in the figure caption.
At pag. 5, lines 30-33 should go to the method, where actually lines 24-27 at pag. 4 tell more or less the same.
The sentence is removed.
At page. 5, lines 1, I would change the sentence into: "The soil moisture characteristics were different for pastures and abandoned areas.Though, at the surface, the soil moisture was close to zero until March in both sites (Fig. 2)." The sentence is modified according to the suggestion of the reviewer.
And here it would be interesting to comment about freezing and melting processes in the uppermost soil layers.In fact, the authors speak about these processes in the abstract and in the conclusion, but without really using their data.But, actually, this fact appears in their data!At the end of December 2014, the soil moisture abruptly decreased after a cold period (Ta down to -15 °C) with no snow at ground… I would imagine then that this decrease was related to the freezing of the liquid water in the soil.I imagine that at this point of the manuscript the authors might give this example to support their sentences (in the abstract and at the end of the conclusion) about the importance of freezing and melting process in the uppermost soil layers.It is simply a suggestion… Thanks for this suggestion.We considered it and added a sentence in the discussions and conclusions.
We thank the editor for the helpful comments and considered two scientific articles which are related to the topic snow gliding and vegetation.

Figure 3 .
Figure 3. Histograms of topographic properties and vegetation characteristics at the glide shoes.

Figure 4 .
Figure 4. Boxplots of the most relevant variables in period I and period II (selected according to Tab. 2).Differences2, indicated by a difference >0.06 of exp(B) from 1).Significant differences between the groups are given by different letters and were determined by Whitney-Mann U-test.

Table 1 : Key characteristics of pastures and abandoned/agricultural unused areas. For each land-use type the glide distance and all topographic and vegetation factors are given (mean ± s.e.).
Tasser are part of the Interdisciplinary Research Center 'Ecology of the Alpine Region' within the major research focus 'Alpine Space -Man and Environment' at the University of Innsbruck. For

Table 2 .
Significant parameters without multi-collinearity and exp(B) of two logistic linear regressions for both periods with snow gliding as dependent variable.If exp(B)<1 then the correlation is negative, if exp(B)>1 then it is positive (bold = most relevant variables, indicated by a difference >0.0506 from 1).Bootstrap results are based on 100 bootstrap samples.B = regression coefficient B, s.e.= standard error, exp(B) = odds ratio 5