Mapping of glacial lakes using Sentinel-1 and Sentinel-2 data and a random forest classifier: Strengths and challenges

Glacial lakes pose a serious threat to downstream areas and significantly impact glacier melt. The number and area of lakes has grown in most regions during the last decades due to the ongoing atmospheric warming and retreating glaciers. It is therefore important to identify and monitor these lakes. However, mapping of glacial lakes in alpine regions is challenged by many factors. These factors include the small size of glacial lakes, cloud cover in optical satellite images, cast shadows from mountains and clouds, seasonal snow in satellite images, varying degrees of turbidity amongst glacial lakes, and frozen glacial lake surfaces. In our study, we have developed a fully automated method for mapping glacial lakes across alpine regions including the Python package called GLakeMap. The method uses multi-source data such as Sentinel-1 Synthetic Aperture Radar and Sentinel-2 Multispectral Instrument data, a digital elevation model, and a random forest classifier model. We use multi-source datasets as inputs for rule-based segmentation of images, mainly aiming at extracting glacial lake objects from satellite images using a set of rules. Segmented objects are then classified either as glacial lake or non-glacial lake objects by the random forest classifier model. The method was tested in eight sites across alpine regions mainly located in High Mountain Asia but also in the Alps and the Andes. We show that the proposed method overcomes a majority of the aforementioned challenges to detect and delineate glacial lakes. The method performs efficiently irrespective of geographic, geologic, and climatic conditions of glacial lakes.