Method and analysis

Drifting snow was generated in a cold wind tunnel. Measurements of mass-flux profile and observations of snow surface change were carried out. The compacted snow was broken up and sieved onto the floor of the test section, and the surface was made flat. Three kinds of snow cover, loose, semi-hard and hard, were prepared for the experiments, with mean hardness of 23, 47 and 65 kPa, respectively. The loose snow cover corresponded to that just after sieving, the semi-hard snow cover to moderately sintered snow, and the hard snow cover to highly sintered snow. The snow temperature was about –15°C. A snow seeder was installed near the upwind end, supplying snow particles into the airflow to trigger the bounding motion of snow particles. Details of the wind tunnel are shown in Reference Sato, Kosugi and SatoSato and others (2001).

The mass-flux profile was measured using a snow particle counter (SPC) located about 11 m downwind from the snow seeder and on the center line of the wind tunnel. The SPC was fixed on a traverse device and moved from z = 2 cm to z = 30 cm, stopping at each of ten heights for 5 s. The snowdrift transport rate, Q, was obtained from

where q is the horizontal mass flux. The theoretical expression for the mass-flux profile by Reference KawamuraKawamura (1948) was fitted to the measured profile, and the surface mass flux was obtained regressively (Reference Sato, Kosugi and SatoSato and others, 2001). This surface value was used to calculate the first term on the righthand side of Equation (1). The second term on the righthand side was calculated from the measured mass-flux values. Since the mass flux decreases markedly with increasing height, as shown in Sato and others (Reference Sato, Kosugi and Sato2001), Equation (1) approximates the total snowdrift transport rate.

Drifting snow is in a state of equilibrium if the numbers of impacted and ejected snow particles are the same. This means that the equilibrium state is governed by the splash process at the snow surface. On the other hand, the airflow cannot carry an infinite number of aeolian snow particles. The situation where the airflow involves the maximum number of snow particles at a certain wind speed is called saturation, which is controlled by the wind field near the snow surface. If the drifting snow does not change along the fetch, the drifting snow is in an equilibrium state and/or saturated. In some literature, this is called a stationary state or a steady state.

In the case of loose snow cover, part of the surface was eroded and the rest was covered with accumulated snow between the snow seeder and the point of SPC measurement. This suggests that the drifting snow was almost saturated. In the case of hard snow cover, erosion was prevented by the strong bonding of snow particles, nor did accumulation occur. Therefore, the drifting snow was considered to be in an equilibrium state but not saturated. In the case of semi-hard snow cover, only erosion occurred and the drifting snow developed along the fetch by incorporating newly ejected snow particles. Thus, the drifting snow was neither in an equilibrium state nor saturated.

Here, a virtual group of snow particles representing actual bounding snow particles is introduced. It bounds on the snow surface and its saltation length is denoted by L. If the snow cover is not hard, the mass of the virtual group of snow particles will increase during every impact as it erodes the snow surface. The rate of mass increase at an impact is expressed by the ejection factor, f_E, in this study. After the virtual group of snow particles impacts n times while it travels a distance, X, from the starting point, the snowdrift transport rate is given by

where N = X/Land_Q0 is the snowdrift transport rate at the starting point, whichcorresponds to the mass supplied by the snow seeder. As shown by Reference KobayashiKobayashi (1972), the mean saltation length can be estimated using a box-type drift gauge.

Results of experiments

The mean saltation length on semi-hard snow cover is plotted against the friction velocity, u_*, in Figure 1. The saltation length is longer on hard snow cover and shorter on loose snow cover compared to that on semi-hard snow cover, which is shown in Kosugi and others (in press, fig. 9), where their wind speed U can be converted to u_* by multiplying 0.04. The mean saltation length was substituted for L and their regression line was obtained as

Fig. 1. Relationship between mean saltation length on semi-hard snow cover and friction velocity. The regression line, Equation (3), is also shown.

The snowdrift transport rate measured on loose snow cover is shown in Figure 2 together with the values observed in the field by previous researchers. Although the observed values were originally expressed as functions of wind speed, friction velocity is used as a parameter in Figure 2, which was converted assuming the log-law of wind profile with a roughness length, z₀, of 10^–4m except that z₀ = 2.5 × 10^-4 m was used for Reference Budd, Dingle and RadokBudd and others (1966). Since the measurements are close to the uppermost values in the field, they can be regarded as the saturated-snowdrift transport rate, Q_sat, which is expressed as a function of friction velocity by

Fig. 2. Relationship between snowdrift transport rate and friction velocity. The three large solid circles are the measurements on loose snow cover, each of which is a mean of several measurements. Relationships obtained in the field by previous researchers are also shown.

From the experiments on semi-hard snow cover, the ejection factor was obtained using Equations (1–3). Results are shown in Figure 3. The ejection factor is slightly larger than unity for friction velocities of less than about 0.3 ms^–1, and it increases with increasing friction velocity.

Fig. 3. Relationship between ejection factor and friction velocity for semi-hard snow. Mean values are connected by a solid line.

Method

The development of a saltation layer along the fetch can be simulated by taking the ejection factor into account. Two kinds of snow cover, semi-hard and hard, were considered. Sublimation of aeolian snow particles during travel was neglected, which corresponds to low-temperature and/or high-humidity conditions. For a given wind speed at z = 10 m, U₁₀, the corresponding friction velocity was obtained assuming the log-law of wind profile and z₀ = 10^-4 m. The U₁₀ value was specified as 5.5, 8 and 12 m s^–1 for semi-hard snow cover and 12ms^–1 for hard snow cover. In the case of semi-hard snow cover, the saltation length, L, was calculated using Equation (3), and the ejection factor, f_E, was given from Figure 3. For hard snow cover, L = 2.0m (Kosugi and others, in press) and f_E = 1.02 were used. The fetch distance, X, was calculated from X = nL, where n is the frequency of impact. Since various disturbances can trigger drifting snow, it is difficult to specify the snowdrift transport rate at the starting point. In this study, the snowdrift transport rate at X = L, the first landing point, was assumed to be given by

where F is the snowfall intensity (water equivalent). Equation (5) means that the virtual group of snow particles arriving at X = L results from snowfall alone. If the drifting snow without snowfall is considered, the snowfall intensity is set to be zero for X > L, and the snowdrift transport rate after n-times impact, Q_n, can be written as

If the drifting snow is accompanied by snowfall, Q_n can be written as

where the second term represents the contribution of snowfall to drifting snow. The drifting snow will become saturated as the virtual group of snow particles repeats impact. The saturated-snowdrift transport rate, Q_sat, was substituted for Q_n if the calculated value exceeded Q_sat.

Results of simulation

Figure 4 shows the change of snowdrift transport rate along the fetch for the drifting snow without snowfall. The snowdrift transport rate increases with the fetch distance and attains a saturated value. If the snowdrift transport rate at the starting point is large, namely F is large, the distance required to attain saturation is short. Hereafter this distance is designated as the developing distance, X_sat. The developing distance is determined by saltation length, ejection factor and saturated-snowdrift transport rate, all of which depend on wind speed. In the case of semi-hard snow cover, X_sat increases as U₁₀ increases from 5.5 ms^–1 to 8 ms^–1, and then decreases as U₁₀ increases to 12ms^–1. In the former range of U₁₀, the wind-speed dependence of X_sat is substantially determined by both the increase in saturated-snowdrift transport rate and the increase in saltation length because the ejection factor is almost constant. In the latter range of U₁₀, however, the effect of increase in the ejection factor with increasing wind speed surpasses other effects. Within the range of U₁₀ considered, X_sat is 100–200m at F = 0.01mmh^–1 and 50–100m at F = 1mmh^–1. Compared with semi-hard snow cover, X_sat on hard snow cover is about one order longer due to its smaller ejection factor and longer saltation length. X_sat is about 1000 m for F = 0.01mmh^–1 and 500 m for F = 1 mmh^–1.

Fig. 4. Change of snowdrift transport rate with fetch distance for drifting snow without snowfall.

Figure 5 is the same as Figure 4 except that the drifting snow is accompanied by snowfall. Since the snowfall is incorporated into aeolian snow particles during travel, X_sat is shorter compared to the drifting snow without snowfall. The dependence of X_sat on wind speed at F = 0.01mmh^–1 is similar to that for the drifting snow without snowfall. The difference in X_sat between two types of snow cover is qualitatively the same. However, the hardness dependence of X_sat is not so distinct as that for drifting snow without snowfall. This is particularly marked when snowfall intensity is high.

Fig. 5. Change of snowdrift transport rate with fetch distance for drifting snow accompanied by snowfall.

Reference KobayashiKobayashi (1972) reported that X_sat was 30–60m under wind-speed conditions of 5–13ms^–1. Reference TakeuchiTakeuchi (1980) observed that X_sat was 150–300m for a wind speed of 5– 8ms^–1. Although precise comparison is difficult because there is insufficient knowledge of the snow and meteorological conditions at the time, their values are within a range of the simulated developing distance on semi-hard snow cover. The simulated snowdrift transport rate increases rapidly as it approaches the saturated value, in contrast to the gentle increase observed in the field. Aeolian snow particles absorb the momentum of the airflow, causing the wind speed to decrease, especially when the drifting snow approaches the saturated state. Since the present simulation does not consider this effect, the resultant behavior nearX = X_sat might be somewhat different from the actual one.