Mapping ice cliffs on debris-covered glaciers using multispectral satellite images

Ice cliffs play a key role in the mass balance of debris-covered glaciers, but assessing their importance is limited by a lack of datasets on their distribution and evolution at scales larger than an individual glacier. These datasets are often derived using operator-biased and time-consuming manual delineation approaches, despite the recent emergence of semi-automatic mapping methods. These methods have used elevation or multispectral data, but the varying slope and mixed spectral signal of these dynamic features makes the transferability of these approaches particularly challenging. We develop three semi-automated and objective new approaches, based on the Spectral Curvature and Linear Spectral Unmixing of multispectral images, to map these features at a glacier to regional scale. The transferability of each method is assessed by applying it to three sites in the Himalaya, where debris-covered glaciers are widespread, with varying lithologic, glaciological and climatic settings, and encompassing different periods of the melt season. We develop the new methods keeping in mind the wide range of remote sensing platforms currently in use, and focus in particular on two products: we apply the three approaches at each site to near-contemporaneous atmospherically-corrected Pléiades (2 m resolution) and Sentinel-2 (10 m resolution) images and assess the effects of spatial and spectral resolution on the results. We find that the Spectral Curvature method works best for the high spatial resolution, four band Pléaides images, while a modification of the Linear Spectral Unmixing using the scaling factor of the unmixing is best for the coarser spatial resolution, but additional spectral information of Sentinel-2 products. In both cases ice cliffs are mapped with a Dice coefficient higher than 0.48. Comparison of the Pléiades results with other existing methods shows that the Spectral Curvature approach performs better and is more robust than any other existing automated or semi-automated approaches. Both methods outline a high number of small, sometimes shallow-sloping and thinly debris-covered ice patches that differ from our traditional understanding of cliffs but may have non-negligible impact on the mass balance of debris-covered glaciers. Overall these results pave the way for large scale efforts of ice cliff mapping that can enable inclusion of these features in debris-covered glacier melt models, as well as allow the generation of multiple datasets to study processes of cliff formation, evolution and decline.