A fully automatic framework for sub-pixel mapping of thermokarst lakes using Sentinel-2 images

Mapping and monitoring thermokarst lakes are crucial to understanding the impact of climate change on permafrost regions and quantifying permafrost-related carbon emissions. Several automatic methods based on remote sensing images have been developed for thermokarst lake mapping. However, mixed pixels containing both land and water characteristics in the lakeshore zones pose a significant challenge to the accuracy of these methods. Furthermore, few approaches were able to fully automate the identification of thermokarst lakes without the manual training sample selection or parameter tuning. In this study, we present a fully automatic framework for thermokarst lake mapping using moderate-resolution Sentinel-2 images. The proposed method combines multidimensional hierarchical clustering and sub-pixel mapping (SPM) based on the radial basis function (RBF) interpolation and Markov random field (MRF) (referred to as RBF-then-MRF SPM), so as to achieve thermokarst lake mapping at a spatial resolution of 3.3 m. We apply the proposed method to two representative thermokarst lake distribution regions in the Northern Hemisphere and achieve a mean Kappa coefficient of 0.89 and 0.99, and a mean $Quality$ of 89.86% and 96.60% on the central Tibetan Plateau and the northern Seward Peninsula, respectively. The results demonstrate that the proposed method significantly improves the accuracy of mixed pixel extraction, and the automatic thermokarst lake mapping is applicable to diverse permafrost regions.