## Introduction

A snow load is given by the water equivalent of a snow-pack. Data from snow gauging networks are frequently insufficient for establishing the snow accumulation characteristics for different regions and altitudes. In particular in most countries the direct measurements of the water equivalent of snow are incomplete although their expediency and reliability is gradually improved by new instruments. On the other hand, daily observations of the total depth of snow cover are generally available and can be used for an indirect reconstruction of data.

## Settling Curves of Snow

When the maximum load of a winter is reached, snow has a varying density depending on the preceding build-up of the snow-pack. At Weissfluhjoch, 2 540 m a.s.l., a density of 490 kg m^{-3} was measured at this point in the winter 1975 while in 1949 it was only 290 kg m^{-3} and in lower stations it can be less than 200 kg m^{-3} These examples illustrate the errors to be expected if one converts the snow depth at the point of maximum accumulation into the water equivalent using a constant ratio.

The densification process is reflected in a gradual decrease of snow depth. A settling curve derived following Haefeli (Bader and others, 1939, p. 80-81)

where *H* is the initial snow depth, *H* is the snow depth after the elapse of time *t*, *ρ*_{s}
is the initial density of snow, *ρ*_{i}
is the density of ice, *α, β* are parameters related to the temperature, vertical load and type of snow, leads logically to the final depth *H*
_{so} (*ρ*
_{s}/*ρ*
_{i}) corresponding to ice. The settling curve can assume varying shapes which can be simulated by laboratory experiments and described by a more physical approach (de (Quervain, 1945, 1946).

When there is an immediate practical need of estimating the settling of a natural snow cover, it is not possible to take into account the variation of temperature, wind velocity, load of superimposed snow layers and partly the intermittent melting or rain, especially if the pertinent data are not available. In order to meet the demands, a simplified settling equation with a time function can be derived from measurements of snow layers in conditions of a natural snow cover, such as

where *H*
_{so} is the depth of new snow measured on the first morning after snowfall and *H*_{sn}
is the depth of snow after *n* days.

The validity of such an equation is limited to snow conditions and to the time period in which it was derived. The theoretical merit of converging to a depth corresponding to ice after a very long time is not fulfilled (such extrapolation is anyhow questionable due to the increasing effect of evaporation and is of no practical interest for a seasonal snow cover). The aim is to obtain a good estimate of snow density attained at the point of the maximum snow accumulation (in terms of the water equivalent) of a winter.

Equation (2) was originally derived in the Krkonoše mountains in Central Europe ((Martinec, 1956). In Figure 1 it is compared with an example of a settling curve obtained with regulated temperature conditions in a laboratory.

Fig. 1. Comparison of settling curves, 1. de Quervain (1945), T = -2.5°C, ρ_{a} = 115 kg m^{-3}, 2. (Martinec (1956), natural snow cover 1 410 ma.s.l.

## Simulation of the Water Equivalent of Snow

The mean density of snow is related to its water equivalent *H*_{w}
and depth *H*_{s}
by the equation

where *ρ*
_{w} is the density of water and, substituting in Equation (2),

Neglecting the evaporation, *H*_{wn}
equals *H*_{wo}
and Equation (3) becomes

In order to calculate *ρ* and *H*
_{wn} a reasonable value of *ρ*_{0}
is thus needed.

When an assessment of snow loads in Switzerland was requested for a revision of building standard specifications, Equation (2) was found also to approximate the average conditions in the Alps if *ρ*_{0}
= 100 kg m^{-3} was substituted equally as in previous applications. By a computerized procedure, the sequence of daily total depths was converted into daily water equivalents so that the maximum could be determined for each year. Figure 2 shows this simulation compared with direct measurements carried out twice a month in two Alpine stations.

Fig. 2. Snow water-equivalents computed from snow depths compared with direct measurements.

An approximation which is satisfactory in different localities seems surprising, since the procedure depends on the substitution of an initial snow density which is variable. An effect of temperature was established by (Diamond and Lowry 1953, p. 1-3) and the wind velocity as well as the type of the falling snow crystals are also involved. According to Equation (4) a range of density *ρ*_{0}
from 80 kg m^{-3} to 120 kg m^{-3} would cause an uncertainty of ±20% in computing *ρ*_{n}
after *n* days. The actual scatter is far better probably for the following reason: Equation (2) indicates the initial rate of settling for the density of 100 kg m^{-3}. If the new snow has a greater density, it can be assumed that this stage, although caused by other effects, is equivalent to a certain time of deposition after which the rate of settling according to Equation (2) is already reduced. A snow density of 120 kg m^{-3} then corresponds to a fictitious time shift of 0.85 d and 80 kg m^{-3} to ̶ 0.5 d (backwards from the starting time of Equation (2)).

If the settling and density of snow are computed for example for *t* = 60 d the following deviations result:

Thus, when variations of the initial snow density ranging from 80 kg m^{-3} to 120 kg m^{-3} are ignored by assuming always *ρ*
_{0} = 100 kg m^{-3}, the resulting deviations after 60 d are within ±0.5%. Recalling that the actual snow depths are available, *H*_{s 60.85}
is only used to determine *ρ*_{60.85}
(which can be also computed by Equation (4)). It is further understood that if the starting snow depth for a layer with *ρ*
_{0} = 120 kg m^{-3} is 1 m, *H*_{s}
for *ρ* = 100 kg m^{-3 } would have been 1.2 m. The actual snow depths are given by the settling curve starting with these values but taken up at the time *n* = 0, *ρ*_{0}
= 120 kg m^{-3}, *H*_{s0}
= 1 m. Consequently the actual snow depth after 60 d is

and the resulting water equivalent of the snow

The error for *ρ*_{0}
= 80 kg m^{-3} is by analogy

In the given example, *H*
_{w 60} is computed as

In a real case, the slower decrease of snow depths can also result from subsequent snowfalls. Therefore

where *ᵀ* is a medium age of all new snow layers (*ᵀ* ˂ 60 d). The diminution of error with time is due to the characteristic shape of settling curves.

Since it takes usually several weeks or months to reach the maximum snow accumulation of a winter, the simple Equation (2) can be applied to estimate the maximum water equivalent of snow although it is not able to take into account the variable starting conditions of each snow layer.

Figure 3 shows two different patterns of snow accumulation with the computed water equivalents and densities. In converting a sequence of daily snow depths to water equivalents, the decrease of snow depth was attributed to the increase of the snow density neglecting a possible small effect of evaporation. On the other hand, an intermittent melting of snow could have affected the results especially in lower stations. Therefore, whenever the snow depth decreased more than indicated by the settling equation, it was taken as a signal of snow melt and the corresponding amount of water was subtracted.

Fig. 3. Different patterns of snow accumulation, station Disentis, 1 173 m a.s.l.

By using the snow depth instead of the elapsed time as an index of the snow density, as for example in the formula (Schriever and Otstavnov, 1968, p. 14)

where *ρ* is in kg m^{-3}, *H*
_{s} is in cm and *T* is the average temperature in °C for the three coldest months of the winter, results contradictory to the age of snow are obtained:

## Frequency Analysis of Water Equivalent of Snow

The maximum water equivalent of each winter’s snow was computed for 35 stations with daily observations of the total depth of snow cover only. Together with direct measurements, 76 stations ranging from 276 m a.s.l. to 2 540 m a.s.l. with a total of 1 482 years of record were available for assessing the expectancy of snow loads. The data agreed well with the Fisher-Tippett extremal distribution

where *P* is the probability that an extreme value *X* is equal to or smaller than the variate *x* and *a* and *c* are parameters.

The parameters have been evaluated (Chow, 1964) as

where *γ* = 0.57721...(Euler’s constant), *μ* is the mean value and *σ* is the standard deviation.

In Figure 4, the maximum water equivalents of each winter’s snow measured at Weiss-fluhjoch are plotted against their respective return periods *T*_{r}
on the Gumbel-Powell probability paper designed for the Fisher-Tippett distribution, with (Weibull, 1939, p. 29)

where *N* is the number of items, *m* is the order number of items arranged in descending magnitude, and

where *T*_{r}
= 1/(1-*P*).

Fig. 4. Maximum yearly snow water-equivalents at Weissfluhjoch, 2 540 m a.s.l., related to return periods.

The relation between the average return period and the expected yearly maximum of the water equivalent of snow at Weissfluhjoch was evaluated by the least-squares method as

Similar equations have been determined for all other stations. Exceptionally, the highest points departed upwards from the general path. This was considered as a sign of a recurrence interval greater than the length of record and such data were excluded from the fitting.

The satisfactory performance of the Fisher-Tippett distribution in this study does not prove its general validity for the frequency analysis of the water equivalent of snow. The Log-normal and the Log-Pearson Type 3 distributions, for example, have been found suitable in other cases (U.S. Weather Bureau, 1964; McAndrew, 1973).

As shown in Figure 5, the relative increase of the expected water equivalent with increasing return period generally decreases with the altitude. The frequency analysis enables the water equivalents to be determined for any desired return period within the limits of a reasonable extrapolation so that values of a comparable probability of occurrence become available for further evaluation.

Fig. 5. Increase of the snow water-equivalent with the growing return period, average rates for the respective altitude ranges.

## Snow Accumulation Related to the Altitude

In mountainous Switzerland the altitude is the dominant factor in snow accumulation. Figure 6 shows polynomial regressions

where *H*
_{w} is the water equivalent in mm for *T*
_{r} = 10 years and *h* is the elevation in m a.s.l., of which the equation

was found satisfactory for practical purposes. The curves in Figure 7 have been obtained for return periods ranging from 2 years to 100 years. All available data from 76 stations have been used. Curves in Figure 8 are based on selected stations with the biggest snow accumulation in the respective elevation zones.

Fig. 6. Snow accumulation in Switzerland related to the altitude.

Fig. 7. Average relations between the maximum yearly snow water-equivalent and the altitude for various return periods.

Fig. 8. Extreme curves of the maximum yearly snow water-equivalent related to the altitude for various return periods.

With regard to the irregular distribution of rare events, the risk of a value characterized by a design return period *T*
_{r} to be exceeded in a given period *N* (the desired lifetime of a structure) is obtained as

where *N* is the number of years.

Probabilities for various alternatives are given in Figure 9. For *T*
_{r} = 10 years and *N* = 10 years the risk is 65%. A value characterized by *T*
_{r} = 100 years is needed to reduce this risk to 10%.

Fig. 9. Calculated-risk nomograph. (*T*
_{r} in years.)

The scatter of points in Figure 6 reflects local and regional deviations in the accumulation of snow. According to an earlier study (Zingg, 1968), the regional effect in Switzerland is less important than that of the altitude but it can become a major factor in vast areas with a changing climate influenced for example by the Great Lakes in North America (U.S. Weather Bureau, 1964) or by the distance from the sea ((Schaerer, 1971).

## Conversion of Snow-Cover Data to Snow Loads

By determining snow water-equivalents to be expected on the ground, the problem of snow loads on roofs is not yet solved. Factors involved in the different accumulation of snow on roofs have been summarized (Tobiasson and Redfield, unpublished, p. 15) as follows:

where *p*
_{r} is the roof snow load, *p*
_{g} is the ground snow load, *c*
_{r} is the regional ground-to-roof conversion factor, *c*
_{e} is the exposure of the structure, and *c*
_{t} is the thermal characteristic of the roof.

A variety of effects can be considered, for example ((Martinec, 1975)

where *p*
_{r} is the roof snow load [in Pa = Nm^{-2}], *H*
_{w} is the snow water-equivalent on the ground [in mm of water = kg m^{-2}], *c*
_{w} is the effect of wind, *c*
_{f} is a coefficient of the roof form, *B* is the snow melt from the ground [in mm of water = kg m^{-2}], *M* is the snows melt by the heat transfer from the building [mm of water = kg m^{-2}], *R* is the incident rainwater [in mm of water = kg m^{-2}], *A* is the load from sliding snow or avalanche [in mm of water = kg m^{-2}], and *g* = 9.81 m s^{-2}.

The aerodynamic effects have been found very variable ((Schriever and others, 1970). At the same time, the snow water-equivalent on the ground might have been underestimated by assuming a snow density of 192 kg m ^{-3}. It is difficult to reconcile maximum snow loads on the ground in British Columbia corresponding to 600 mm water-equivalent ((Peter and others, 1963) with ground load-altitude relations ((Schaerer, 1971) or with snow surveys ((Snow Survey Bulletin, 1972) whereby snow water-equivalents exceeding 2 000 mm or even 3 000 mm and snow densities around 500 kg m^{-3} have been measured.

The uncertainty in the ground-to-roof conversion is reflected in Building Codes valid in various countries ((Schriever and Otstavnov, 1968, p. 14). In most cases the altitude serves as an index of snow load. Countries like Italy understandably emphasize local climatic conditions. The form factor is accounted for by the slope of the roof resulting in considerable reductions of the snow load.

## Conclusion

In order to improve data for the structural design in Switzerland, the snow water-equivalent was indirectly determined in the absence of measurements and a frequency analysis was carried out to ensure a comparability of values and a choice of design return periods.

Foreign studies demonstrated great differences between the roof load and ground load depending on specific conditions. These effects, together with regional anomalies of the snow accumulation, have to be considered as further steps leading to more realistic indications of roof loads, especially in high areas and exposed localities.

## References

Bader,, H.,
*and others*. 1939. Der Schnee und seine Metamorphose, Bader,, von H.
Haefeli,, R.
Bucher,, E.
Neher,, J.
Eckel,, O.
Thams,, C.
Niggli., P.
*Beiträge zur Geologie der Schweiz. Geotechnische Serie. Hydrologie* Lief. 3. [English translation: *U.S. Snow, Ice and Permafrost Research Establishment. Translation* 14, 1954.]

Chow,, V. T.
1964. Statistical and probability analysis of hydrologie data. *In*
Chow,, V. T., ed. Handbook of applied hydrology. A compendium of water-resources technology. New York, etc.,McGraw-Hill, p. 8-1–8-97.

Diamond,, M.
*and*
Lowry,, W. P.
1953. Correlation of density of new snow with 700 mb temperature. U.S. Snow, Ice and Permafrost Research Establishment. Research Paper 1.

McAndrew,, D. W.
1973. Snow load analysis and recreational uses of snow data. *In*
Santeford,, H. S.,
*and*
Smith,, J. L., comp. Advanced concepts and techniques in the study of snow and ice resources... Washington, D.C.,National Academy of Sciences, p. 11–21.

Martinec,, J.
1956. Zimni prognosys pouzitim radioisotopü [Winter forecasts with the use of radioisotopes]. *In*
*Vltavska kaskada [The Vltava reservoir system]*. Prague-Podbaba, Vyzkumny Üstav Vodohospodafsky, p. 45–60.

Martinec,, J.
1975. Periodizität der Schneelasten in der Schweiz. *Interner Bericht des Eidg. Institutes für Schnee- und Lawinenforschung* Nr. 545.

Peter,, B. G. W.,
*and others*. 1963. Variation of snow loads on roofs, by Peter,, B. G. W.
Dalgliesh, W. A. and Schriever, W. R.. Transactions of the Engineering Institute of Canada Vol. 6, No. A-1, p. 1–12
.

Quervain,, M. R. de.
1945. Die Setzung der Schneedecke. Interner Bericht des Eidg. Institutes für Schnee- und Lawinenforschung Nr. 11.

Quervain,, M. R. de.
1946. Mechanisch-kristallographische Untersuchungen, 4. Zur Setzung der Schneedecke. Interner Bericht des Eidg. Institutes für Schnee- und Lawinenforschung Nr. 26.

Schaerer,, P. A.
1971. Variations of ground snow loads in British Columbia. Proceedings of the Western Snow Conference 38th annual meeting, p. 44–48
.

Schriever,, W. R.
*and*
Otstavnov,, V. A.
1968. Snow loads. *Canada. National Research Council. Division of Building Research. Research Paper* No. 366, p. 13–33
.

Schriever,, W. R.
*and others*. 1970. Snow loads, by Schriever, W. R. and Lutes, D. A. and Peter., B. G. W.
*In Canadian structural design manual. Supplement No. 4 to the National Building Code of Canada*. Ottawa, Canada. National Research Council. Associate Committee on the National Building Code, p. 567–78
. (NRC No. 11530))

Snow Survey Bulletin. 1972. [Issued by British Columbia. Dept. of Lands, Forests, and Water Resources. Water Resources Service. Water Investigation Branch.]

Tobiasson,, W.
*and*
Redfield,, R. Unpublished. Alaskan snow loads. [Paper presented at 24th Alaska Science Conference,9 1973. Obtainable from CRREL. Hanover, N.H., U.S.A.]

U.S. Weather Bureau. 1964. Frequency of maximum water equivalent of March snow cover in north central United States. *U.S. Weather Bureau. Technical Paper* No. 50, p. 1–24
.

Weibull,, W.
1939. A statistical theory of the strength of materials. Ingeniörsvetenskapsakademiens Handlingar (Stockholm), Nr. 151.

Zingg,, T.
1968. Maximale Schneelasten in der Schweiz. *Schweizerische Bauzeitung* 86. Jahrg., Ht. 31, p. 555–57
.

## Discussion

L. LLIBOUTRY: Did you find that precipitation continues to increase with altitude? We found in the Western Alps that it goes through a maximum which depends upon the altitude of the surrounding summits and then declines with further increase in altitude.

J. MARTINEC: The highest station included in the relation between altitude and water equivalent of snow was 2 540 m a.s.l. For the practical use of this relation, it should be limited to 2 000 m a.s.l. I agree that at the highest altitudes the increase in precipitation with altitude can cease or even be reversed.