Introduction

In British Columbia, Canada, the primary method of harvesting timber is clear-cut logging which consists of removing all the timber in large tracts of land (typically areas of 10–100 ha). This practice has been banned from mountain slopes in the Alpine countries of Europe for > 100 years and is not allowed in mountainous Scandinavia or Japan (Weir, in press). Since British Columbia is largely mountainous and since there is ample snow supply, clear-cut logging has resulted in significant environmental damage due to the effects of snow avalanche activity. Our research has shown that approximately 10 000 clear-cuts have been significantly affected by snow avalanche activity in British Columbia. The information in this paper represents the first compilation of data about avalanche terrain and snow supply characteristics for avalanches which have initiated as a result of clear-cut logging anywhere, not just in British Columbia (Weir, in press).

Our research has shown that there are two classes of environmentally destructive problems for avalanches interacting with clear-cuts: type I, where avalanches initiate in clear-cuts, and type II, where avalanches descend into clear-cuts. The information in this paper is strictly focused on type I events even though type II events can also cause significant environmental problems. Type I events are known to cause severe environmental damage and risk, including destruction of forest cover, prevention of regeneration of young seedlings, threats to facilities (dwellings, power lines and roads), stripping soil cover down to bedrock, discharge of debris into streams and initiation of debris torrents and landslides.

In this paper, I describe terrain parameters and forest cover for 76 avalanche paths which have initiated in clear-cuts (type I). The study areas comprise two parcels of land (each about 20 000 km2) in the Coast Mountains and the Columbia Mountains (eastern British Columbia), as shown in Figure 1. The database consists of 40 terrain and forest-cover parameters for each site collected by extensive air-photo searches, combined with information from maps and site visits. In addition, theoretical estimates of snow supply (mean and standard deviation of Gumbel distributions for maximum annual water equivalent) are derived for each site based on altitude and location from 70 snow-course sites maintained by British Columbia: Ministry of Environment, Lands and Parks (2000).

Fig. 1.

Map of British Columbia showing general study areas in Coast Mountains (four sub-regions) and Columbia Mountains (three sub -regions), with sub-regions defined for estimating mean maximum annual water equivalent of snow with elevation.

The results of the paper provide information about starting-zone variables collected from avalanches initiating in clear-cuts that is necessary to help characterize destructive avalanches in clear-cuts. Primary variables include average start-zone angle, start-zone surface roughness, vegetation height and coverage, aspect, slope curvature and snow-supply estimates. The objective of this study is to provide the statistical data framework as input to risk-based decisions to help minimize risk of damage for future logging plans.

Rationale for Parameter Selection and Snow-Climate Descriptions

In order to minimize risk of environmental damage in the future, the chance of avalanching (related to avalanche frequency) will have to be considered along with potential avalanche magnitude. A study of terrain parameters for forested avalanche paths at Rogers’ Pass, British Columbia, by Reference Smith and McClungSmith and McClung (1997) showed that all relevant parameters in their study which correlated significantly with avalanche frequency were directly or indirectly related to snow supply. Further, field observations combined with physical reasoning indicate that certain terrain characteristics in starting zones are important in preventing avalanche formation including: ground surface roughness, vegetation density and vegetation height. Accordingly, parameters in this study were chosen to represent those expected to be important for governing snow supply, avalanche size, avalanche frequency and parameters known to be important to inhibit avalanche formation (Reference McClung and SchaererMcClung and Schaerer, 1993).

The sites in this study are from two different snow-climate areas (Fig. 1) in British Columbia, each about 20 000 km2 in area: the Coast Mountains (maritime snow climate) and the Columbia Mountains (transitional snow climate). The maritime region of the Coast Mountains is subject to high snow supply, with rain possible any time during the winter. The transitional snow climate of the Columbia Mountains usually has less snowfall than coastal areas, with higher chance of persistent weak snow layers Reference JamiesonJamieson, 1995) which can affect the character of avalanching. See Reference McClung and SchaererMcClung and Schaerer (1993) for a more detailed comparison of maritime and transitional snow climates. The variables related to terrain and forest cover for the individual avalanche paths are described in the next section, and the estimates of snow supply (mean and standard deviation of maximum annual mm water equivalent (directly related to snow climate)) are described in a later section.

Data Collection Methods and Accuracy Limitations: Terrain and Forest-Cover Parameters

Terrain and forest-cover parameters for all avalanche sites were assessed by a combination of three methods: (1) air photos (scales ranging from 1:15 000 to 1:40 000, with most 1:20 000); (2) terrain contour maps (scale 1:20 000); and (3) site visits which included visual inspection, distance measurements with a laser range finder (accuracy ±1 m), altimeter readings (±50 m) and slope angle measurements (accuracy 1° with clinometer measurements).

Many of the sites were found from inspection of provincial air-photo collections maintained in Victoria, British Columbia, but most were found from field trips throughout the study areas. Each site was visited in summer to ensure that any disturbances noted were due to snow avalanches instead of landslides, debris flows or other damage types. Whenever possible, data were measured in the field, but when this was not possible, air photos and contour maps were used in combination to estimate parameters. Sometimes in the field, vegetation obscured the view or sites were not accessible due to impassable logging roads or terrain, so that a combination of maps and air photos was necessary to complete the database at such sites.

The combination of field measurements and estimates from contour maps and air photos implies that individual data elements have variable data accuracy across the datasets. However, it is felt this limitation will not have much effect on the use of the results, which are aimed at future application of risk-based methods based on probabilistic data trends. The advantages of the mixed accuracy data-collection methods include:

1. (1) more data are available within reasonable time and expense constraints than if precision information had been sought;

2. (2) the data precision is essentially the same as that available to potential users of the information;

3. (3) variable data accuracy from one site to another should not greatly affect the statistical trends which are most valuable for the risk-based approach intended for the future use of the information.

Individual Terrain and Forest-Cover Parameters: List, Description and Method of Measurement or Estimation

Determination of Snow-Amount Estimates: Mean and Standard Deviation of Annual Maximum Water Equivalent

No data are available for snow amounts at the locations studied. However, snow amounts have been found to be very important for estimating avalanche frequency (Reference SchaererSchaerer, 1977; Reference Smith and McClungSmith and McClung, 1997). Therefore, snow amounts were estimated theoretically from snow-course data collected at 70 sites throughout the study areas by the Ministry of Environment, Lands and Parks, Province of British Columbia (data records from 1935–99 contained on the World Wide Web (http://www.elp.gov.bc.ca/wat/snow-bulletin/archive/historic. html). Two scales of snow-amount patterns were found to be important from the data analysis: (1) snow amounts and characteristics (the coefficient of variation: ratio of standard deviation to the mean of maximum annual snow water equivalent) were significantly different for the two major snow-climate areas (Coast and Columbia Mountains); and (2) snow-amount variations with elevation were found to vary significantly within regions of each major climate area, similar to the findings of Reference Claus, Russell and SchaererClaus and others (1984). Therefore, the goal of the present study was to estimate snow amounts as a function of elevation within sub-regions for both the Coast and Columbia Mountains.

The procedure for determining the snow supply at a place involves locating the avalanche path within a sub-region (Fig. 1) and then applying a sub-region linear regression model of mean maximum annual snow with elevation to the top elevation of the start zone (Fig. 2). The standard deviation is then found by applying the coefficient of variation calculated for the general snow climate (Coast or Columbia Mountains) as described below (Fig. 3).

Fig. 2.

Mean maximum water equivalent with elevation for three sub-regions of the Columbia Mountains. Linear least-squares lines are shown with R as the correlation coefficient.

Fig 3.

Mean maximum (μ) and standard deviation (a) in mm w.c.for Gumbel distributions from 70 snow-course sites in British Columbia. The points are determined by fitting a Gumbel distribution to annual maximum water equivalent of snow for each site. Least-squares lines are fitted through the data for Coast Mountain (□) and Columbia Mountain sites (•). R is correlation coefficient.

Snow Amounts for Coast and Columbia Mountains Compared

For each snow-course site, the parameters of a Gumbel distribution were determined by linear least-squares regression between ranked data for annual maximum snow amount (mmw.e.) and the reduced variate {−ln[−ln(P_i)]} corresponding to the datum where Pi is non-exceedance probability. The non-exceedance probability, P_h is related to the ranked values by Pi = (i – 0.5)/N, where i is rank of the datum and N is number of years of record. For example, the highest value for a station with 15 years of record is related to P15 = 0.97, indicating that 97% of values are not expected to exceed the value at the location, and the reduced variate corresponding to the highest point is 3.49. All stations had at least 10 years of records.

The mean, μ, and standard deviation, a, of the Gumbel distributions for each site were then determined by fitting a least-squares line. The slope of the line gave the scale parameter, b, and the intercept gave the location parameter, u for the Gumbel distribution at the site. The Gumbel mean is then given by μ = u + γb (where 7 is Euler’s constant, 0.57721… and the Gumbel standard deviation is Per cent variance explained (R2(%), where R is Pearson correlation coefficient) by the least-squares fits was 84–99%, usually exceeding 95%.

Figure 3 shows plots of mean, μ, and standard deviation, a, for the 70 individual Gumbel distributions. Least-squares fits through the points show that μ and a are related by two different linear regression lines in Figure 3 corresponding to two different coefficients of variation for these differing snow climates (Coast or Columbia Mountains). In the analysis below, the value of μ with elevation is determined within each sub-region, and the corresponding value of a can be determined from the general relations in Figure 3.

Sub-Region Analysis: Mean Maximum Water Annual Equivalent within the Coast and Columbia Mountains

Based on Gumbel means for mean maximum annual water equivalent individual snow-course stations, the Coast and Columbia Mountains were divided into four and three subregions, respectively (see Fig. 1). For each sub-region, a linear least-squares fit of Gumbel mean maximum annual water equivalent vs elevation was determined. Figure 2 shows least-squares lines and mean values from the Columbia Mountains. The sub-regions were determined by inspection such that a nearly linear relation was determined for each sub-region. Statistics from these least-squares lines showed that the per cent variance explained (R2) was 68–95%, indicating good fits in most cases. From the plots, a top avalanche starting-zone elevation within a sub-region is used to find an approximate mean maximum water equivalent for the site and, from Figure 3, the standard deviation can be estimated for the mean in the main region. For example, for a height of 1500 m in sub-region 3 of the Columbia Mountains, mean maximum annual water equivalent is 260 mmw.e. (Fig. 2) and standard deviation is 80 mm we. (Fig. 3).

The theoretical procedure above is similar but not equivalent to that outlined by Reference Claus, Russell and SchaererClaus and others (1984) for determining height dependence of water equivalence. Claus and others did not explore climate variation of the coefficient of variation as in this study (Fig. 3). Similar to the analysis above, Claus and others determined that sub-regions within a general snow climate are important, and used parabolic rather than linear relations to estimate snow amount as a function of altitude. However, in the higher altitudes where the avalanche starting zones in this study were found, the parabolic relations found by Claus and others can be approximated as nearly linear. The study presented here is somewhat simpler, and many more stations are available along with 20 years more data.

At a location, it may be assumed that maximum annual snow amounts will obey a Gumbel distribution based on the analysis of more than 70 sites. Given the mean value μ estimated at a site (altitude and sub-region) and the estimated value of a (or coefficient of variation) for the general climate region, the parameters (u, b) of the Gumbel distribution can be determined. An approximate value of water equivalent for any return period, x_T, may be estimated from: x_T < {−ln[−ln(P_i)]} where T is return period. For analysis in this paper, only the maximum mean annual water equivalent (μ) is calculated for each site, rather than estimates for other return periods.

Method of Analysis for Terrain Variables

The variables described above represent characteristics of terrain, forest cover and snow supply for avalanche starting zones caused by clear-cut logging. The terrain variables are of two general types: categorical (e.g. wind index, down-slope and cross-slope shape, ground surface roughness) and numerical (e.g. start-zone slope angle, mean maximum annual snow water equivalent). In order to display the characteristics, categorical variables are shown as cumulative samples (cumulative distributions for discrete sample points) indicating the fraction of the data in each category or bar graphs. Numerical variables are plotted as cumulative distribution functions or probability density functions. The intention here is to highlight the characteristics of the individual variables as they relate to avalanche initiation in clear-cuts, with attention given to variables which should be related to magnitude and frequency of avalanches. All avalanches in the database released by natural processes as opposed to human triggering.

Results of Analysis

For the analysis below, Ms denotes sample mean, a_s denotes sample standard deviation and N is number of data points.

Fig. 4.

Probability plot for a Gaussian distribution (-) compared with data for average start-zone slope angle (•) for 76 avalanche sites.

Fig. 5.

Probability plot for a Gaussian distribution (––) compared with mean maximum water equivalent estimates • (Coast Mountain avalanche sites, N = 58).

Fig. 6.

Probability plot for a Gaussian distribution (––) compared with mean maximum water equivalent estimates • (Columbia Mountain avalanche sites, N = 18).

Fig. 7.

Bar graph for starting-zone wind index for 76 avalanche sites. The five categories (Reference SchaererSchaerer, 1977) are defined in the text. As the wind index increases, the starting zone is progressively more exposed to wind loading of snow.

Fig. 8.

Cumulative distribution plot for mean vegetation height in start zones. The plot is similar to a cumulative distribution probability function for continuous values.

Fig. 9.

Bar graph of downslope sha_pe for 76 start zones. Values: –2, extremely convex; 0, no curvature; +2, extremely concave. Categories are defined in the text.

Fig. 10.

Bar graph of cross-slope shape for 76 start zones. Values as in Figure 9 caption. Categories are defined in the text.

Fig. 11.

Cumulative distribution plot for avalanche size. Based on the five-part Canadian size classification estimated from destructive effects (Reference McClung and SchaererMcClung and Schaerer, 1993, p. 253). The plot is similar to a cumulative distribution probability function for continuous values.

Descriptive Statistics of Clear-Cut Area, Damaged Area and Starting-Zone Area

The statistics in this section yield information about the areas involved in this study. All areas are estimated from maps and air photos in hectares and shown in Table 1. The data in Table 1 show that, on average, the damaged area is about twice the estimated start-zone area, and that average clear-cut size far exceeds either of these areas.

Table 1.

Descriptive statistics for clear-cut area, start-zone area and damaged area

Conclusions and Discussion

The data and analyses in this paper provide information on the characteristics of destructive avalanches resulting from clear-cut logging. While the information comes from two different climate areas, the only clear distinction between characteristics for these areas is in terms of snow supply and top elevation of starting zones, both of which are consistent with what is expected: higher elevations in the Columbia Mountains and greater snow supply in the mari-time-climate area of the Coast Mountains.

Important results include: (1) no relation between avalanche size, average start-zone angle and mean maximum snow amount; (2) a prevalence of events on concave slopes, with concave features being more pronounced in the cross-slope than in the downslope direction; (3) potential importance of ground roughness, vegetation height and vegetation coverage in inhibiting events; (4) a prevalence of events for windward and leeward slopes with respect to prevailing storm direction (generally southwest in British Columbia): about 50% higher than adjacent 90° aspect sectors.

Important limitations of the present study (aside from data accuracy) include the difficulty in incorporating local terrain features which may be crucial for avalanche initiation, and the lack of on-site information about avalanche release conditions. It is expected (Reference Smith and McClungSmith and McClung, 1997) that snow supply will influence expected avalanche frequency so that general avalanche frequency due to this effect might be expected to be higher in the Coast Mountains, on average. However, the events in this study are primarily large ones which presumably initiate from deep weak layers such as surface hoar and facets or deeply buried crusts. It is possible that surface hoar and facets may be more prevalent and persist longer in the transitional snow climate of the Columbia Mountains so that a counterbalance may be in place to offset lower snow supply there to increase the frequency of large destructive avalanches. However, even though this seems physically reasonable, it cannot be proven by the data contained in this study.