1. Introduction
The central arctic experiences the most occurrences of blowing snow in Canada, while the Canadian prairies have significant number of events as well (Phillips 1990). Blowing snow can be a hazard to public safety and transportation because of poor visibility associated with some events. An accurate forecast of the occurrence and resultant visibility of this phenomenon may help to mitigate this hazard. Unfortunately, predicting the occurrence of blowing snow and the associated reduction in visibility is difficult due to complexities related to snowpack state and its interaction with the lower atmosphere. Blowing snow occurs when the wind shear stress exceeds the snow particle resistance to dislocation (Li and Pomeroy 1997a). Thus, the occurrence of blowing snow will vary spatially and temporally with most of the variation attributable to differences in wind speed, air temperature, and characteristics of the snow.
Although there have been several studies devoted to blowing snow from a climatological, hydrological, or pure research perspective (e.g., Lawson 2002; Li and Pomeroy 1997a, Li and Pomeroy 1997b; Dery and Yau 2001), there remains a need for a practical forecast method of the phenomenon. Pomeroy and Male (1988) developed a scheme for estimating the reduction of visibility due to blowing snow; however, it is dependent upon parameters that are in most cases not known to the forecaster (e.g., drift density and mean particle radius). Li and Pomeroy (1997b) included agricultural stubble height and local fetch distances as predictors but these elements vary on a scale smaller than usually considered by a weather forecaster. Li and Pomeroy (1997b) defined the frequency of occurrence of blowing snow for any given wind, temperature, and snowpack age as the ratio of the number of events with blowing snow divided by the number of events without. However, the utility of a forecast can be harmed by “crying wolf” or forecasting too many false alarms. This present paper attempts to put more emphasis on the false alarm aspect by defining a credibility factor. This study is aimed at creating a practical technique for the prediction of significant blowing snow and the associated visibility reductions through available observational data.
From an operational forecasting perspective, the definition of “significant” blowing snow varies with the application. For instance, a public weather forecast is mainly concerned with visibility less than 1 km, whereas the aviation industry operates using established visibility thresholds including 0.8, 1.6, 3.2, 4.8, and 8 km (½, 1, 2, 3, and 5 mi).
Knowledge of snow state and meteorological conditions are required to forecast blowing snow. While detailed snowpack properties are not typically available to the weather forecaster in real time, a record of snow depth, time since last snowfall (proxy for snow age), and interim weather conditions should be known.
The paper first discusses the historical meteorological station data that were used in this study and the methods employed to develop the blowing snow forecast scheme (section 2). Section 3 tests the ability of the scheme to forecast blowing snow and visibility at various Canadian arctic and the Prairie provinces stations, with concluding remarks to follow in section 4.
2. Data and methods
a. Historical meteorological observations
The Environment Canada weather stations (Fig. 1) were selected based on the completeness of their observational record, their proximity to large population centers, and aviation terminal weather forecast locations. Hourly weather observations (air temperature, wind direction, mean wind speed, prevailing visibility, and precipitation) dating from 1960 to 1999 were used from each station. Prevailing visibility is defined as the maximum visibility value common to sectors comprising one-half or more of the horizon circle (Atmospheric Environment Service 1977). Wind gust data were also included in a portion of the database. Intermediate observations, or “specials,” were not available. Some stations in the arctic reported weather every 3 h in the early 1960s. Consequently, the portions of observations with incomplete records at these particular stations were omitted from this study. Observations from October to April were selected as most blowing snow events occur over these months (>95%).
Blowing snow is defined as snow particles raised by wind to sufficient heights above the ground to reduce horizontal visibility to 9.7 km (6 mi) or less (Atmospheric Environment Service 1977). Horizontal visibility is measured at 1.8 m above the ground. Visibilities recorded during periods of darkness were modified in this study. Rasmussen et al. (1998) found that visibility is reduced to a lesser degree at night compared to daylight hours for the same rate of falling snow. By analogy, this process was applied to visibilities in blowing snow. Sunrise and sunset times were calculated for each station based on date, latitude, and longitude. A daytime equivalent visibility was calculated for snow and/or blowing snow by multiplying nighttime visibilities by a factor of 2/3 (Rasmussen et al. 1998).
For illustration, a reported hourly mean wind of 20 km h−1 gusting to 40 km h−1, where it is assumed that wind cycles from 0 to 40 km h−1, has the same energy as a steady wind of 25 km h−1. While this assumption oversimplifies the chaotic nature of gusty winds, subsequent testing indicated that this adjustment showed an improvement over the unmodified mean wind speed.
Based upon previous snow work done by Jordan et al. (1999) and Hanesiak et al. (1999), the daily maxima net radiation on the snow surface are better reflected in maximum daily temperatures. The latter would therefore have more influence on the aging of the snowpack than the temperature at the time of the blowing snow event, which in turn should influence the blowing snow processes. The highest temperature since the last snowfall was used as an independent variable in our analysis.
Thus the reported wind speed may be adjusted for anemometer height, gustiness, and/or temperature. Incorporating station barometric pressure will affect air density and should therefore have some theoretical influence on blowing snow. Testing did not show any measurable difference, and would represent a complication to the forecasting process. Thus, the effect of barometric pressure was not incorporated into any of the results.
b. Analysis and forecast scheme
To test whether or not regional differences were important, all of the observations from arctic1 stations were combined and compared to a collection of reports from the Prairie provinces. Blowing snow events were divided into observations with and without concurrent precipitating snow. For observations with blowing snow with visibility 1 km or less in the Prairie provinces, 77.9% occurred with precipitating snow, whereas in the arctic this figure was just 31.9%. The blowing snow events that occurred without concurrent falling snow (ground blowing snow) were grouped by snow age. Since the likelihood of blowing snow decreased with time after the last occurrence of falling snow, the time span of each age group was set larger at greater ages (e.g., the age groups were set to 1–2, 3–5, 6–11, 12–23, . . . 192–383 h since last snowfall). It was found that this grouping scheme maintains a reasonable sample size in each snow age group.
There were a few cases of blowing snow with extremely old snowpacks (>1 month), or with snowpacks that had undergone liquid precipitation or significant melting. We suspect these cases to be blowing snow being advected from a distant location, but this could not be determined using point observational data. To reduce the effect of advected blowing snow on the results, all positive events and nonevents that occurred after a melting period with temperatures greater than 4°C (<5% of the total hits) were excluded in the calculations. Similarly, all observations that occurred after liquid precipitation were discarded. Individual reports that included fog, smoke, blowing dust, or ice pellets were removed from the dataset.
Site-specific topography may preclude blowing snow with certain wind directions at particular observation sites. A climatology of blowing snow events was derived for each station and for each wind direction (in tens of degrees) with and without concurrent precipitating snow. A similar climatology can be found in Hanesiak et al. (2003). If there were no cases of blowing snow from a certain wind direction at a particular location, null events from that direction and location were also removed from the analysis. Since blowing snow is much less common in forested areas than in the open prairies, stations such as La Ronge, Saskatchewan, contributed very little to the final statistics. A forecaster may mimic this process by checking the climatology or the characteristics of the terrain and vegetation of a particular site before predicting the occurrence of blowing snow (Hanesiak et al. 2003).
With concurrent precipitating snow, a positive event is defined as any hourly report of visibility less than or equal to a given specified threshold (public or aviation limits) with observed light snow and blowing snow. A nonevent is defined as any hourly observation with reported light snow, but with visibility greater than the specified threshold. The latter group may contain reported blowing snow provided the visibility remains greater than the threshold. Light snow is defined as falling snow giving visibility greater than or equal to 1.01 km (5/8 mi). Observations with moderate and heavy snow were excluded, as this precipitating snow would give low visibility whether or not blowing snow occurred. These definitions also apply to cases without concurrent precipitating snow except that the hourly observations will not include any falling snow.
The mean and variance of this probability function was calculated at each temperature using an automated iterative method that found the best fit to the observed CRED (Fig. 2). Critical wind speeds at the 5%, 50%, and 95% CRED were computed at each temperature with Eq. (5) using the best mean and variance. The individual critical wind speeds were replotted over the range of tested temperatures and the best-fit quadratic lines at each of the 5%, 50%, and 95% CRED were derived using standard statistical software (Fig. 3).
The thresholds represented by the 5%, 50%, and 95% CREDs wind speeds represent a historic correlation but could be used as a nomogram for predicting upcoming blowing snow events. The estimate of CRED can be obtained by applying the expected temperature and wind speed. This value of CRED could then be used as an indicator of the probability of blowing snow in the future. It also shows how sensitive the likelihood of blowing snow can be to the speed of the wind. Small differences in wind speeds within the 5%–95% limits can produce large variations in the CRED. Wind speeds outside of these limits are not so critical.
In summary, the CRED of significant blowing snow occurrence was calculated over a range of temperatures and wind speeds extracted from the quality-controlled observations. Critical wind speeds were interpolated from a probabilistic function that was best fit to the observed CRED points at each temperature. Other best-fit equations then described these critical wind speeds over a range of temperatures. The result was an entirely objective estimate of the critical winds associated with significant blowing snow. This process was repeated for every permutation of
snowpack age groups;
specified visibilities of 0.8, 1.0, 1.6, 3.2, 4.8, and 8 km;
using just the mean wind, and also by incorporating the mean wind and gust as available;
all stations aggregated together, and as two groups divided as stations from the arctic and those from the Prairie provinces; and
using TCurrent versus using TMax for ground blowing snow.
The number of combinations is large and precludes including all the graphs in this paper. Figure 3 is one example of a critical wind nomogram.
3. Results and discussion
a. Blowing snow frequency
Four arctic stations (Baker Lake, Coral Harbour, Resolute Bay, and Hall Beach) account for over half the blowing snow events of the entire test group. Some of this is explained by the length of the arctic cold season. Each site in the full test group will have local effects caused by terrain and land usage that will affect the occurrence of blowing snow. To test the overall significance of these local effects, the correlation between the incidence of blowing snow and the characteristic of the wind was calculated.
The Pearson moment correlation coefficient, r (Naiman et al. 1977), was used to determine the dependency of blowing snow reports with visibilities less than 1 km on the character of the wind. The wind characteristics tested were the mean speed, and the number of occurrences of speeds above set thresholds. An r of 0 would indicate that none of the variation in the one variable could be attributed to a second while an r of 1 would indicate that there is a complete positive linear correlation between the two sets of variables.
Table 1 shows that using the frequency of any wind from a given direction (threshold = 0; Fig. 4b) shows a higher r than using the mean wind from a given direction (Fig. 4a). In other words, there is a better correlation of the number of significant blowing snow events recorded at each station and direction to the number of wind occurrences from a particular direction, as opposed to the simple speed average of those occurrences. The number of strong wind events (threshold =7.5 m s−1; Fig. 4c) showed a higher correlation than the frequency of winds of any speed. The number of hours with speeds above the calculated level (threshold = critical speed at 5% CRED; Fig. 4d) has a correlation of over 90% to the number of significant blowing snow occurrences for the entire sample.
The result that the best correlation occurs with winds above the CRED threshold indicates that the local topographical effects of a particular station determines significant blowing snow occurrence largely through its influence on the wind speed. Therefore, it would be valid to apply the forecast nomograms derived from this set of sites to other locations provided one could correct the forecast wind speed for the local effects of that area. Some stations in the arctic do have a high number of blowing snow events but this is mainly due to having more wind, cold temperatures, and loose snow. This may also suggest that a blowing snow climatology could be approximated for any station, provided it had a good record of temperature and wind speed (Hanesiak et al. 2003).
b. Duration of blowing snow
Blowing snow will cease when the wind diminishes and/or the snow becomes too compacted or all the loose snow is scoured away. In the observed events in which snow and blowing snow has occurred together, blowing snow usually terminates when precipitating snow ends (Table 2). The remaining percentage of blowing snow events that persist after the precipitating snow has ended decreases rapidly with time (Fig. 5). However, some blowing snow events in the arctic have lasted weeks after the falling snow has ended. Note that the end of an event is defined as the first report not reporting blowing snow, or reporting an improvement in visibility above the set threshold. Blowing snow may or may not have redeveloped later.
c. Air temperature dependency of blowing snow
As discussed by Li and Pomeroy (1997a), there is a nonlinear relationship between the frequency of blowing snow and temperature. The 5% CRED critical wind for blowing snow with concurrent precipitating snow is relatively constant over the range of temperatures (Fig. 6a). Ground blowing snow typically requires a stronger wind, especially at mild and very cold temperatures. One should expect that the critical wind would increase sharply near 0°C, but this was not always captured by the resulting curves as the number of occurrences of blowing snow drops off dramatically at temperatures near melting. The minima in critical wind speed for freshly fallen snow are at a milder temperature than for well-aged snow. Older snowpacks do not necessarily require a stronger wind to generate blowing snow compared to fresh snow. The graphs even suggest less wind is required for older snowpacks at some temperatures. Keep in mind that there were a limited number of cases of blowing snow with well-aged snowpacks (>192 h since last snowfall) and even fewer of these with temperatures near 0°C. The arctic does experience long periods of cold and darkness that can contribute to a slow aging process for the snowpack. There were almost no blowing snow events with well-aged snowpacks in the record from the Prairie provinces. Figure 6b shows that blowing snow with lower specified visibilities consistently requires stronger winds through the whole range of observed temperatures.
Figure 7a shows a simple count of reports of visibility less than 1 km in blowing snow at each temperature for various snowpack ages. Figure 7b has the same distribution, but presented as a proportion of occurrences of temperatures at each snow age. In other words, Fig. 7b represents the number of significant blowing snow reports at a given temperature and snow age divided by the total number of reports of that same temperature and snow age. All hourly reports that occurred after melting conditions or after liquid precipitation were discarded. The result is the number of hits in a given set of circumstances normalized by the number of times those circumstances arise. The small number of cases of blowing snow with very old snowpacks leads to some scatter in the plots and is not shown. Blowing snow with concurrent falling snow does not show a strong dependency on temperature. These normalized frequency distributions of ground blowing snow peaks at an air temperature near minus 34°C regardless of the age of the snowpack even though the minima in the threshold wind are at milder temperatures for young snowpacks. These charts also indicate that at a given temperature, the likelihood of ground blowing snow is largely independent of the age of the snowpack.
d. Visibility distribution
One technique of forecasting visibility might begin with a method of predicting any snow movement using a critical wind speed (Li and Pomeroy 1997b) and then applying a visibility distribution in blowing snow by wind speed in excess of this threshold (Fig. 8). Instances of blowing snow are rare with wind less than 5 m s−1 and very few of these are significant. Even at the 5% CRED critical speed, a minority of blowing snow cases has low visibilities. The distribution of visibilities varies with temperature (Figs. 8a–d). Colder temperatures favor a higher proportion of low visibilities at marginal wind speeds. The same trend was observed for older snowpacks. Not only is blowing snow more likely at higher wind speeds, a greater proportion of these events have significant visibility reduction. The specified visibility approach (Fig. 6a) where the critical wind speed is determined for set visibility restrictions simplifies the forecast process by requiring one nomogram instead of two. A suggested forecast process is provided in section 4.
e. Evaluating wind thresholds for specified visibilities
The 5% CRED critical wind was objectively evaluated as a potential forecast tool by testing it as a discriminator between cases where blowing snow with visibilities less than 1 km would or would not be expected. Observations with actual wind speeds above the 5% CRED critical wind would constitute a positive “forecast” of significant blowing snow that could then be verified against the reported weather. Various 5% CRED critical winds were derived using the data manipulation criteria listed in section 2b. The value of these criteria could then be assessed by their critical winds ability to distinguish cases of significant blowing snow.
Tables 3a and 3b shows that the 5% CRED critical wind is a poorer discriminator of blowing snow events in cases with concurrent falling snow (CSI of 34.49%) than for those without (55.65%). Light falling snow on its own can have a visibility restriction close to 1 km. It would be very difficult to isolate the visibility restriction component due solely to the blowing snow. Blowing snow with falling snow does not have a strong temperature dependency, and so the 5% CRED critical wind is only marginally better than using a single wind speed threshold of 7.6 m s−1 for the whole temperature range. However, using a temperature independent speed of 9.2 m s−1 for ground blowing snow has a CSI of 44.14%, which is much poorer than using the 5% CRED critical wind (55.65%).
Generally, using thresholds calculated for arctic and Prairie province groups separately shows a small improvement over using one threshold equation derived from the two groups combined (Table 3b). There were not enough cases of ground blowing snow in the Prairie provinces in the dataset with both mean wind and gusts to justify separate statistics. Using separate thresholds based on snow age is beneficial as is using the maximum temperature since the last snowfall (Table 3b, c). One would expect these two effects to be more important with older snowpacks. However, events with aged snow are fewer in number than those with freshly fallen snow. Therefore improving the discrimination of old snow events would not be expected to make great improvements in the overall CSI. Incorporating gust information makes a more significant improvement over just using the mean wind (Table 3c) as wind gusts will affect any snow age group. The CSI with gusts is 66.10% compared to 63.00% without.
4. Summary and conclusions
This paper has examined the occurrence of significant blowing snow events on the Canadian arctic and the Prairie provinces derived from historical hourly weather observations at 32 locations. Various descriptive statistics display some blowing snow characteristics. A critical success index of 66% can be achieved in discriminating blowing snow occurrences with certain specified visibilities using wind speed thresholds derived from this dataset. The highest CSI resulted when incorporating the most input data. This includes using information of wind gust, snow age, and highest air temperature since the last snowfall. Knowledge of the frequency and circumstances of past blowing snow occurrences at each station was also beneficial.
Developing the wind speed thresholds based on specified visibility restrictions proved to be the key in the performance of this technique of distinguishing conditions that are likely to generate blowing snow. Not only does the likelihood of blowing snow increase with wind speed, an increasing proportion of these events produce low visibilities. The coupling of these two trends leads to a sharper division between conditions that give high-impact blowing snow from all other cases. While blowing snow remains a chaotic process, forecast success is more likely if one focuses the objective to be predicting events giving very restricted visibilities. It is also these events that have the highest impact on transportation and public safety.
The wind speed thresholds found in this study could be employed as a method of forecasting future events. Although this is a purely empirical technique, it can be useful to the forecaster as an objective first guess of significant blowing snow events. A suggested process is as follows.
Forecast a wind speed and direction, and air temperature. Blowing snow prediction requires a precise and accurate wind forecast. The forecaster must account for the effects of local topography on wind speed and direction. A standard wind rose may be useful in this task.
Examine a blowing snow rose for the site in question, if applicable. These diagrams are useful for point forecasts but may not apply to nearby areas if the topography at the observation site is not representative of the surrounding area.
Check and note the snowpack history, if applicable. Do not forecast ground blowing snow if the snowpack has undergone melting, or has been exposed to liquid precipitation. Be aware that blowing snow may advect in from a distant location.
Use the appropriate nomogram to find the CRED. The nomogram selection would be based on the specified visibility, snow age, day versus night, current versus maximum temperature, and mean wind versus a combination of mean wind and gusts. For example, if one were expecting a wind of 60 km h−1 at −10°C, Fig. 3 would suggest a CRED around 80%. The forecaster would know that while blowing snow will likely occur, it is not necessarily guaranteed to do so. The forecast of wind and temperature will include some amount of uncertainty. Therefore, one should find CRED for a range of these forecast parameters. The spread of the extracted CRED will be an indicator of the confidence the forecaster should have in this process. In the example above, an error range of just 1°C and 10 km h−1 could produce a CRED as low as 30% or as high as 95%. On the other hand, if the temperature was expected to be −30 instead of −10°C, the CRED exceeds 95% over a similar range of uncertainty in the forecast inputs.
Knowledge of the local terrain, vegetation, fetch, climatology, surface temperature, depth of snow on the ground, etc., may be used to subjectively adjust the forecast. The aerial occurrence of blowing snow can vary greatly due to these factors. The effects of these elements are described in Li and Pomeroy (1997b).
Acknowledgments
The authors wish to thank Frank Svistovski (Commercial Weather Services Division) for his advice on obtaining the meteorological observations; and Bevan Lawson, Ron Hopkinson and Rick Raddatz [Atmospheric and Hydrologic Sciences Division (AHSD)] for their contributions in converting wind observations to a standard reference height and helpful comments on the paper. We acknowledge Rachel Jordan (Cold Regions Research and Engineering Lab) and John Pomeroy (University of Saskatchewan) for discussions on snow and modeling.
REFERENCES
Atmospheric Environment Service, 1977: Manual of Surface Weather Observations (MANOBS). 7th ed. Atmospheric Environment, 322 pp.
Davenport, A. G., 1960: Wind loads on structures. Division of Building Research, National Research Council, Canada, NRCC-5576, 60 pp.
Dery, S. J., and Yau M. K. , 2001: Simulation of an arctic ground blizzard using a coupled blowing snow–atmosphere model. J. Hydrometeor., 2 , 579–598.
Elliot, D. L., 1986: Wind Energy Resource Atlas of the United States. Solar Technical information Program, Solar Energy Research Institute Rep. DOE/CH 10093-4, 210 pp.
Hanesiak, J. M., Barber D. G. , and Flato G. M. , 1999: The role of diurnal processes in the seasonal evolution of sea ice and its snow cover. J. Geophys. Res., 104 , C6,. 13593–13604.
Hanesiak, J. M., Fisico T. , and Carriere E. , 2003: Climatology of adverse weather events in the Canadian arctic. Centre for Earth Observation Science Tech. Rep. CEOS-Tech-2003-2, 3540 pp. [Available from Faculty of Environment, University of Manitoba, Winnipeg, MB R32 2N2, Canada.].
Jordan, R., Andreas E. A. , and Makshtas A. P. , 1999: Heat budget of snow-covered sea ice at North Pole 4. J. Geophys. Res., 104 , 7785–7806.
Lawson, B., 2002: Trends in blizzards at selected locations in the Canadian prairies. J. Natural Hazards, 29 , 2,. 123–138.
Li, L., and Pomeroy J. W. , 1997a: Estimates of threshold wind speeds for snow transport using meteorological data. J. Appl. Meteor., 36 , 205–213.
Li, L., and Pomeroy J. W. , 1997b: Probability of occurrence of blowing snow. J. Geophys. Res., 102 , 21955–21964.
Naiman, R., Rosenfield R. , and Zirkel G. , 1977: Understanding Statistics. 2d ed. McGraw-Hill, 307 pp.
Phillips, D. W., 1990: The Climates of Canada. Canadian Government Publishing Centre, 176 pp.
Pomeroy, J. W., and Male D. H. , 1988: Optical properties of blowing snow. J. Glaciol., 34 , 116,. 3–10.
Rasmussen, R. M., Vivekanandan J. , and Cole J. , 1998: The estimation of snowfall rate using visibility. J. Appl. Meteor., 38 , 1542–1563.
Pearson moment correlation coefficient of blowing snow occurrence visibility ≤ 1 km and the character of the wind for compass directions in tens of degrees
Proportion of concurrent snow and blowing snow events that end when the precipitating snow ends
Table 3a. Blowing snow with visibility ≤ 1 km and concurrent falling snow
Table 3b. Blowing snow with visibility ≤ 1 km but without concurrent falling snow
Table 3c. Blowing snow with visibility ≤ 1 km but without concurrent falling snow
Arctic stations are all north of 60° north latitude. Churchill, Manitoba, was also treated as an arctic station.