
Recently I generated one year of snow maps from Sentinel-2 in the Canadian Rockies for a talented colleague who is working on the numerical simulation of the snow cover at high-resolution with an exciting new hydrological model. The area is not covered by Theia, hence I used Start_Maja to generate the L2A products and then LIS to generate the snow masks on the CNES supercomputer (thanks!).

This area is quite challenging for snow optical remote sensing: the terrain is steep and there are a lot of forests. After processing this area, I also found some unexpected issues in the LIS processor, which need to be fixed, like turbid rivers detected as snow, or wildfires smoke detected as snow. I tried to use Pekel’s water mask to remove the rivers and lakes pixels but there are some glaciers that are misclassified as water in this product hence I simply masked out areas below 2000 m. This eliminates most of the water surfaces but not all the forests, hence I used Hansen’s global forest product to mask out pixels with a tree cover density larger than 50%.

As a sanity check I often look if areas with high snow cover duration match the glacierized areas. Here I zoomed on Haig Glacier (50°43’N, 115°18’W) in Kananaskis Country and plotted the glaciers outlines from the Randolph Glacier Inventory 6.0:

As noted by Marshall (2014), « Haig Glacier is the largest outlet of a 3.3 km² ice field that straddles the North American continental divide », which means that melt water from Haig Glacier flows towards the Atlantic Ocean!