Modelling hazardous surface hoar layers across western Canada with a coupled weather and snow cover model

•We modelled surface hoar formation and evolution at sites across western Canada.•Widespread formation events were modelled, but over-predicted.•Hazardous layers were modelled and identified with a structural stability index.•Regional-scale patterns between major mountain ranges were resolved.

Destructive snow avalanches in western Canada are often caused by failure in buried surface hoar layers. Numerical snow cover models can simulate the formation and evolution of these layers, which could help avalanche forecasters assess the location and timing of avalanches. To investigate this application, we compared modelled surface hoar layers with snow and avalanche observations from the Coast, Columbia, and Rocky Mountains of western Canada. Surface hoar formation and evolution was modelled by forcing the snow cover model SNOWPACK with data from a high-resolution numerical weather prediction model. Surface hoar formation was verified with daily snow surface observations at 88 observation sites over two winters, and the evolution of buried layers was verified with avalanche observations and persistent weak layer assessments from 67 avalanche forecast regions. The frequency of surface hoar formation was over-predicted by 40%, although since modelled crystal sizes were moderately correlated with observed sizes, the more hazardous layers were often distinguished. A structural stability index in SNOWPACK often identified surface hoar layers during storm slab avalanche activity, but identifying layers during persistent slab avalanche activity was more difficult. Model limitations included uncertain meteorological inputs, errors in SNOWPACK's snow surface energy balance, and representing the necessary spatial scales. Despite these limitations, the coupled model resolved differences between major Canadian mountain ranges and could improve avalanche forecasts in data-sparse regions.