Glacial Lake Area Changes in High Mountain Asia during 1990–2020 Using Satellite Remote Sensing

Changes in a large-scale glacial lake area directly reflect the regional glacier status and climate changes. However, long time series of glacial lake dataset and comprehensive investigation of the spatiotemporal changes in the glacial lake area in the whole High Mountain Asia (HMA) region remained elusive. Satellite remote sensing provides an indispensable way for dynamic monitoring of glacial lakes over large regions. But glacial lakes are quite small and discretely distributed, and the extraction of glacial lakes is usually influenced by clouds, snow/ice cover, and terrain shadows; thus, there is a lack of an automatic method to continuously monitor the dynamic changes of glacial lakes in a large scale. In this paper, we developed a per-pixel composited method named the “multitemporal mean NDWI composite” to automatically extract the glacial lake area in HMA from 1990 to 2020 using time-series Landsat data. There were 19,294 glacial lakes covering a total area of 1471.85 ± 366.42 km(2) in 1990, and 22,646 glacial lakes with an area of 1729.08 ± 461.31 km(2) in 2020. It is noted that the glacial lake area in the whole HMA region expanded by 0.58 ± 0.21%/a over the past three decades, with high spatiotemporal heterogeneity. The glacial lake area increased at a consistent speed over time. The fastest expansion was in East Kun Lun at an average rate of 2.01 ± 0.54%/a, while in the Pamir and Hengduan Shan, they show slow increases with rates of 0.33 ± 0.08%/a and 0.39 ± 0.01%/a, respectively, during 1990–2020. The greatest increase in lake area occurred at 5000-5200 m a.s.l., which increased by about 45 km(2) (~25%). We conclude that the temperature rise and glacier thinning are the leading factors of glacial lake expansion in HMA, and precipitation is the main source of lake water increase in West Kun Lun. Using the proposed method, a large amount of Landsat images from successive years of melting seasons can be fully utilized to obtain a pixel-level composited cloud-free and solid snow/ice-free glacial lake map. The uncertainties from supraglacial ponds and glacial meltwater were also estimated to improve the reliability and comparability of glacial lake area changes among different regions. This study provides important technical and data support for regional climate changes, glacier hydrology, and disaster analysis.