High‐resolution snowline delineation from Landsat imagery to infer snow cover controls in a Himalayan catchment

The snow cover dynamics of High Mountain Asia are usually assessed at spatial resolutions of 250 m or greater, but this scale is too coarse to clearly represent the rugged topography common to the region. Higher-resolution measurement of snow-covered area often results in biased sampling due to cloud cover and deep shadows. We therefore develop a Normalized Difference Snow Index-based workflow to delineate snow lines from Landsat Thematic Mapper/Enhanced Thematic Mapper+ imagery and apply it to the upper Langtang Valley in Nepal, processing 194 scenes spanning 1999 to 2013. For each scene, we determine the spatial distribution of snow line altitudes (SLAs) with respect to aspect and across six subcatchments. Our results show that the mean SLA exhibits distinct seasonal behavior based on aspect and subcatchment position. We find that SLA dynamics respond to spatial and seasonal trade-offs in precipitation, temperature, and solar radiation, which act as primary controls. We identify two SLA spatial gradients, which we attribute to the effect of spatially variable precipitation. Our results also reveal that aspect-related SLA differences vary seasonally and are influenced by solar radiation. In terms of seasonal dominant controls, we demonstrate that the snow line is controlled by snow precipitation in winter, melt in premonsoon, a combination of both in postmonsoon, and temperature in monsoon, explaining to a large extent the spatial and seasonal variability of the SLA in the upper Langtang Valley. We conclude that while SLA and snow-covered area are complementary metrics, the SLA has a strong potential for understanding local-scale snow cover dynamics and their controlling mechanisms.