Next Article in Journal
Multiple Images Improve Lake CDOM Estimation: Building Better Landsat 8 Empirical Algorithms across Southern Canada
Next Article in Special Issue
Evaluating the Territorial Impact of Built-Up Area Expansion in the Surroundings of Bucharest (Romania) through a Multilevel Approach Based on Landsat Satellite Imagery
Previous Article in Journal
A Deep Fusion uNet for Mapping Forests at Tree Species Levels with Multi-Temporal High Spatial Resolution Satellite Imagery
Previous Article in Special Issue
The Relative Contributions of Climate and Grazing on the Dynamics of Grassland NPP and PUE on the Qinghai-Tibet Plateau
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 13
Issue 18
10.3390/rs13183614
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rapid glacier Shrinkage and Glacial Lake Expansion of a China-Nepal Transboundary Catchment in the Central Himalayas, between 1964 and 2020

by
Yan Zhong
1,2,3,
Qiao Liu
1,3,*,
Liladhar Sapkota
1,2,
Yunyi Luo
1,2,
Han Wang
1,2,
Haijun Liao
1,2 and
Yanhong Wu
1,3
1
Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu 610041, China
2
College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
3
Branch of Sustainable Mountain Development, Kathmandu Center for Research and Education, Chinese Academy of Sciences-Tribhuvan University, Chengdu 610041, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(18), 3614; https://doi.org/10.3390/rs13183614
Submission received: 4 August 2021 / Revised: 31 August 2021 / Accepted: 8 September 2021 / Published: 10 September 2021
(This article belongs to the Special Issue Temporal Resolution, a Key Factor in Environmental Risk Assessment)
Browse Figures
Review Reports Versions Notes

Abstract

:
Climate warming and concomitant glacier recession in the High Mountain Asia (HMA) have led to widespread development and expansion of glacial lakes, which reserved the freshwater resource, but also may increase risks of glacial lake outburst floods (GLOFs) or debris floods. Using 46 moderate- and high-resolution satellite images, including declassified Keyhole and Landsat missions between 1964 and 2020, we provide a comprehensive area mapping of glaciers and glacial lakes in the Tama Koshi (Rongxer) basin, a highly glacierized China-Nepal transnational catchment in the central Himalayas with high potential risks of glacier-related hazards. Results show that the 329.2 ± 1.9 km2 total area of 271 glaciers in the region has decreased by 26.2 ± 3.2 km2 in the past 56 years. During 2000–2016, remarkable ice mass loss caused the mean glacier surface elevation to decrease with a rate of −0.63 m a−1, and the mean glacier surface velocity slowed by ~25% between 1999 and 2015. The total area of glacial lakes increased by 9.2 ± 0.4 km2 (~180%) from 5.1 ± 0.1 km2 in 1964 to 14.4 ± 0.3 km2 in 2020, while ice-contacted proglacial lakes have a much higher expansion rate (~204%). Large-scale glacial lakes are developed preferentially and experienced rapid expansion on the east side of the basin, suggesting that in addition to climate warming, the glacial geomorphological characters (aspect and slope) are also key controlling factors of the lake growing process. We hypothesize that lake expansion will continue in some cases until critical local topography (i.e., steepening icefall) is reached, but the lake number may not necessarily increase. Further monitoring should be focused on eight rapidly expanding proglacial lakes due to their high potential risks of failure and relatively high lake volumes.

1. Introduction

Glaciers play essential roles in the global cryosphere and Earth’s water cycle [1] by buffering freshwater supply [2], modifying mountain landforms [3], and regulating fluvial sediment flux [4]. Glacier and glacier-related processes are highly vulnerable to climate change [5,6]. The Tibet Plateau and its surroundings, namely the High Mountain Asia (HMA), host the largest volume of glaciers outside of the polar regions [7,8] and is referred to as the Asian water tower [9,10]. Among all, approximately 20,431 (covering 19,679 km2) glaciers and 4950 (covering 455.3 km2) glacial lakes are developed along with the Himalayas [11,12,13,14]. Nearly 800 million people in the Indus, Ganges and Yarlung Tsangpo river (i.e., Brahmaputra River) basins are benefited from the Himalayan glacial meltwaters [15,16].
With a warming global climate, Himalayan glaciers have shrunk severely [16,17,18], accelerating the development and evolution of glacial lakes in this region [19]. The existence of glacial lakes could cause a lag effect of glacial runoff that affects downstream river discharges, hydrological processes, channel erosions, and landscape evolutions [20]. They could also induce glacier-related hazards, e.g., glacial lake outburst floods [19], which is one of the main climate hazards caused by warming in high altitude regions [21]. Himalayan glacial lakes are mainly located between 4000 and 5700 m a.s.l., and expanded by approximately 14.1% from 1990 to 2015 [11]. The widespread expansion of the Himalayan glacial lakes has been the focus of attention by both scientific researchers and social communities [19,22,23,24,25,26,27,28,29,30].
Owing to frequent geo-hazards, regional development delays, poor road conditions, and low access to glaciers and glacial lakes, field investigations are difficult to conduct in the HMA. However, with the recent rapid development of remote sensing technology and cloud-based geospatial data distribution, efficient and reliable mapping methods have been widely used in cryosphere science [31,32]. To better assess and mitigate the risks of GLOFs hazard and their impact, it is crucial to investigate the evolution of the glaciers and glacial lakes and predict their future development trends under climate change [33]. A higher temporal resolution (e.g., annual) of glacial lake mapping and the expansion process of different types of glacial lakes are still limited in previous studies [34]. In this study, 46 mid- and high-resolution remote sensing images based on Landsat and Keyhole series satellite were used to provides a complete mapping (1964–2020) of glacial lakes and glaciers in the Tama Koshi (Rongxer) basin, a less reported China-Nepal transnational catchment in the central Himalayas. We aim to (1) analyze the dynamics of different types of glaciers and glacial lakes and (2) discuss the glacier-lake interactions, possible forcing mechanisms of glacial lakes expansion, their future evolution, and GLOF risk implications.

2. Study Area

The Tama Koshi basin (a sub-basin of the Koshi River Basin, 28°N, 86°E, Figure 1) is a long, narrow, and deep valley extending northeast-southwest, located on the south side of the central Himalayas. The total basin area is 4117 km2 (about 1/3 is in China), and the elevation ranges from less than 500 m a.s.l. at Beni Ghat, where the Tama Koshi joins the Sunkoshi, to above 8000 m a.s.l. adjacent to Mt. Cho Oyu (8188 m a.s.l.). Due to the enormous elevation difference (more than 6800 m), the ecological environment in the basin is unique, complex, and diversified. Data from the European Centre for Medium-Range Weather Forecasts (ECMWF, see Section 3.4) show that the mean annual temperature of the entire basin was 5.42 °C, and the mean annual precipitation is 1380 mm, and both showed increasing trends in the past 40 years, by about 1 °C and 300 mm, respectively (Figure 1B).
Glaciers and glacial lakes are widely distributed in the upper basin. According to the Randolph Glacier Inventory (RGI) 6.0 [13], a total of 279 glaciers are developed above 4600 m a.s.l., covering a total area of 334.8 km2 (8% of the entire basin). According to the latest inventory by Wang et al. [12], glacial lakes in the basin have a total number of 229 and an area of 14.6 km2 in 2018, which can be classified as five major types: cirque lake (11%), end moraine-dammed lake (54%), lateral moraine-dammed lake (5%), moraine thaw lake (1%) and supraglacial lake (16%). As reported by ICIMOD in 2020 [21], 42 lakes in the Koshi basin have been identified as potentially dangerous glacial lakes (PDGL), of which 8 are located in the Tama Koshi basin. A study by Rounce et al. [35] also noted 11 very high-risk lakes in Nepal, one of them is located in this basin. Glacier and glacial lake changes in the Tama Koshi basin have been reported previously, but most with a low spatial-temporal resolution [34,36,37]. One of the most reported glacial lakes in upper Tama Koshi is the Tsho Rolpa (Rolwaling Valley), which has been evaluated as the most dangerous glacial lake in Nepal [38]. The Tsho Rolpa has been monitored, measured, and modeled since the 1990s with comprehensive details about its expansion mechanisms [39] and failure risks [40]. Various risk prevention measurements, e.g., siphon system, automatic early warning system, and artificial channel excavation, have been carried out at the Tsho Rolpa [41]. Several studies have predicted the possible response of glaciers and glacial lakes under the climate change in the Tama Koshi basin, provided the management experience and scheme of dangerous glacial lakes and discussed the potential risks of GLOFs [37,40,42,43].

3. Data and Methods

3.1. Data Set Used

The remote sensing images used in this study were obtained from the United States Geological Survey (USGS, https://earthexplorer.usgs.gov/; accessed on 26 October 2020), Geospatial Data Cloud (http://www.gscloud.cn/; accessed on 30 October 2020) and PlanetLabs (https://www.planet.com/; accessed on 21 October 2020) between 1964 and 2020. A total of 6 KH-4A and KH-9 Keyhole (KH) images (2.7–9 m), 36 Landsat (TM/ETM+/OLI) surface reflectance images (30 m), and 4 PlanetScope (PL) Ortho Tile imagery (3.125 m) were used to manually delineated glaciers in 1964, 1980, 1990, 2000, 2010, 2018, and 2020 and glacial lakes in 1964, 1980, and automatically delineated glacial lakes in 1987–2020, respectively. Images were mainly acquired in post-monsoon and near the end of ablation seasons (August to December) with less cloud cover and seasonal snow cover. The orthorectified PL products were used as base images for the registration of KH images.
Other data sets we have used include: (1) RGI 6.0 [44], (2) Glacial Lake Inventory of High-mountain Asia in 1990 and 2018 [12], (3) ERA5 monthly aggregates from 1979 to 2019 (latest climate reanalysis produced by ECMWF/Copernicus Climate Change Service (ECCCS)), (4) the 30 m Shuttle Radar Topography Mission Digital Elevation Model (SRTM DEM), (5) 2007 ALOS PALSAR DEM (12.5 m), (6) 1975–2000, 2000–2016 glacier surface elevation change for the HMA derived by differencing KH-9 DEM (generated from optical stereo) and ASTER DEMs [18], (7) 1999–2003, 2013–2015 glacier surface flow velocity for the HMA derived by optical satellite images via feature-tracking algorithm [45], (8) LandScan: global population distribution data for 2000 and 2018 [46,47], and (9) 1:72,000 scale topographic map (base map from Esri). The glacier and glacial lake data sets were used to extract relevant attributes (for example, aspect and ID) and provide a reference for shaping delineations, and the DEM and glacier thickness change data were used to extract terrain parameters and glacier surface elevation changes. The attribute information and purpose of each data set are detailed in Supplementary Materials, Table S1.

3.2. Glacier and Glacial Lake Outline Delineation

The KH series images were characterized by high resolution, single band (panchromatic), and lack of spatial coordinates and projection information [48,49], which needs to be preprocessed before mapping. Four high-resolution orthorectified PL images (as the base image) from September to October 2020 were used for the preliminary spatial registration of KH images (as the warp image) using the ENVI Image-to-Image Registration Tool. A total of 103 (for KH-4A) and 72 (for KH-9) ground control points were selected in stable terrains such as ridgelines, river intersections, and moraine ridgelines. According to the ground control point pairs in the two images, a functional relationship is established, and the geometric correction of the distorted image is completed using coordinate transformation.
To monitor the detailed annual change of the rapidly expanding glacial lakes (Figure 1A, A–H), we automatically delineated the glacial lakes (Figure S1) in 1987–2020 (September–November) based on the Google Earth Engine (GEE) remote sensing cloud processing platform [50] Python Application Programming Interface (API) [51]. The Normalized Difference Water Index (NDVI; Equation (1)) and Normalized Difference Snow Index (NDSI; Equation (2)) were combined and applied to multi-temporal Landsat (TM/ETM+/OLI) images to reduce the errors caused by sediments and seasonally frozen lake water. Then, it is complemented with manual editing to extract glacial lake boundaries. The combination of this object-oriented image processing method and expert knowledge has shown suitable performance in glacial lake delineation in the Himalayas [52].
N D W I = g r e e n − N I R g r e e n + N I R
N D S I = g r e e n − S W I R 1 g r e e n + S W I R 1
where green, NIR, and SWIR1 are the green, near-infrared, and shortwave-infrared 1 bands, respectively. By using the Python API interactive interface, the classification thresholds can be adjusted more accurately than the Javascript API interface based on the on-screen inspection of the original images. All Landsat and PlanetScope images were atmospheric and topographic corrected products. Additionally, an error of ±0.5 pixels was estimated to calculate the area uncertainties (i.e., multiply linear error and perimeter, e.g., the linear error of Landsat images are 15 m, and of KH-4A images are ~1.35 m) [33,53]. For supraglacial lakes, only the base-level lakes with perennial water storage were retained, while the perched lakes with large seasonal variations were removed [54].

3.3. Glacier Surface Dynamics Analysis

Eight lake-terminating glaciers (Figure 1A, A–H) and six land-terminating glaciers (Figure 1A, I–N) were selected to compare their dynamics and relationships with glacio-geomorphology or lake-ice interactions. Based on the 2000 SRTM, the two periods of glacier surface elevation in 1975 and 2016 were calculated by combining the 1975–2000 and 2000–2016 glacier surface elevation change. Surface gradient was calculated based on the 2007 ALOS PALSAR DEM using the QGIS slope analysis module, by which we extracted values between the terminal moraine to the mid-point position of the glacier tongue to better explore the slope of the downstream tongue. Subsequently, equidistant points at a 30 m spacing were generated along the entire glacier center flow line to identify and visualize the elevation, flow velocity, and slope. For the elevation and velocity analysis of the ice/lake area, we extended the line to a further 300–800 m in front of the moraine dam.

3.4. Meteorological Data Extractions

Since there is no meteorological observation in the Tama Koshi basin, we used the ERA5 monthly aggregates data to analyze the impact of meteorological factors. The ERA5 data set, with a spatial resolution of 0.25 radian, is the fifth generation ECMWF atmospheric reanalysis of the global climate [55]. Its grided data set is now incorporated in the GEE data catalog, and we extracted the temperature and precipitation data inside the basin domain between 1979 and 2019 pixel by pixel.

4. Results

4.1. Glacier Distribution and Changes

4.1.1. Distribution of Glaciers

Our glacier inventories show (Figure 2A) a total of 271 glaciers with a total area of 329.2 ± 1.9 km2 in 2020 (Table 1). The largest one is Glacier M (42.6 ± 0.1 km2), which is a debris-covered glacier facing south and located in the upper Lapche, one of the western tributaries of the Tama Koshi (Figure 1, Randolph Glacier Inventory ID: RGI60-15.09675). The mean elevation of all glaciers ranges from 4760 to 6270 m a.s.l. (Figure 2B) and most (96%) glacier areas are distributed between 5200 and 5800 m a.s.l. The lowest glacier tongue, which is close to the Rongxer country (ID: RGI60-15.09760), extends to an altitude of 4475 m a.s.l. Due to the northeast (high altitude) to the southwest (low altitude) orientation of the basin, southern oriented glaciers have the maximum area coverage (84.9 ± 0.4 km2, ~26%, Figure 2C). Northwest-oriented glaciers have the second largest area of 68.2 ± 0.3 km2 (~21%), whereas the glaciers with the east orientation are the smallest (10.3 ± 0.1 km2, 3.1%).
Glacier distributions are slightly different between the basin’s east and west sides, divided by the Rongxer/Tama River corridor (Figure 1 and Table 1). On the west side of the corridor, the three largest debris-covered glaciers (L, M, and N in Figure 1) occupied ~63% of the total glacier area (133.7 km2) compared to the other 131 relatively smaller glaciers (50.0 km2). In contrast, mean glacier areas in the east (1.4 km2 for the 138 glaciers) are larger than the west (1.0 km2), with seven glaciers are >5 km2 and five are >10 km2. Glaciers on the east side distribution with a wide altitude range (mainly between 5200 and 5800 m a.s.l.), whereas the west side glacier areas are concentrated between 5600 and 5700 m a.s.l. All five lake-terminating glaciers are located on the east side, but the supraglacial lakes are more frequently found on the lower part of the three largest debris-covered glaciers.

4.1.2. Glacier Changes

From 1964 to 2020, the total area of glaciers in the Tama Koshi basin decreased by 26.2 ± 3.2 km2 (0.13% a−1) in 56 years (Table 1 and Table 2), with a mean glacier terminus retreat rate of 3.2 m a−1. Glacier area loss rates in the periods of 1980–1990 and 2018–2020 were 0.37% a−1 and 0.31% a−1, respectively, which were much higher than the rest study period (0.03–0.14% a−1, Table 1). Eight small glaciers covering an area of 0.03–0.29 km2 had disappeared between 1980 and 1990. The mean retreat rate of lake-terminating glaciers (37.7 m a−1) is about 14 times higher than those land-terminating ones (2.6 m a−1). Glacier G, which terminates in lake G, experienced the longest terminus retreat (68 m a−1) between 1964 and 2020. On the east side of the basin, glacier area loss rate (0.15% a−1) and retreat rate (4.3 m a−1) are generally higher than those on the west (0.10% a−1 and 2.08 m a−1, respectively, Table 2).
Meanwhile, glaciers also showed a trend of thinning and slowdown (Table 2). The mean change rates of glacier surface elevation during 1975–2000 and 2000–2016 were −0.32 m a−1 and −0.63 m a−1, respectively, indicating an accelerated thinning of the ice. Glacier elevation change rates on the western side (from −0.26 m a−1 in the period 1975–2000 to -0.69 m a−1 in the period 2000–2016) shows more remarkable than the east (−0.33 m a−1 in the period 1975–2000 to 0.62 m a−1 in the period 2000–2016) (Table 2). The most downwasted glacier on the west side of the basin is the Glacier L, one of the largest land-terminating glaciers (debris-covered, with a length of 12.7 km) in the upper Tama Koshi. Its mean elevation change rate increased from −0.45 m a−1 in 1975–2000 to −0.97 m a−1 in 2000–2016 (Figure 3), higher than that of Glacier E, I, and G on the east side (−0.78 m a−1, −0.60 m a−1 and −0.89 m a−1 in 2000–2016, respectively). Nevertheless, the mean elevation change rates of lake-terminating glaciers (−0.40 m a−1 in the period 1975–2000 and −0.85 m a−1 in the period 2000–2016) in the basin were higher than that of land-terminating glaciers (−0.30 m a−1 in the period 1975–2000 and −0.62 m a−1 in the period 2000–2016) during the two study periods.
The mean glacier surface velocity of the Tama Koshi basin decreased by ~25% from 5.3 m a−1 in 1999–2003 to 4.0 m a−1 in 2013–2015. Overall, mean glacier surface velocities on the west side (5.4 m a−1 in the period 1999–2003 and 4.1 m a−1 in the period 2013–2015) were slightly higher than on the east side (5.2 m a−1 in the period 1999–2003 and 3.8 m a−1 in the period 2013–2015), while both have decelerated during the two study periods. Lake-terminating glaciers (with a velocity of 6.9 m a−1 in the period 1999–2003 and 6.1 m a−1 in the period 2013–2015) moved faster than land-terminating glaciers (4.8 m a−1 in the period 1999–2003 and 3.4 m a−1 in the period 2013–2015). Even though the surface velocity declined by 11%, it was much slower than that of land-terminating glaciers (−29%).

4.2. Glacial Lake Distribution and Changes

4.2.1. Distribution of Glacial Lakes

The statistical results of glacial lake inventories between 1964 and 2020 are presented in Table 3. A total of 196 glacial lakes (more than 50% are moraine-dammed proglacial lakes) are developed in the basin, with a total area of 14.4 ± 0.3 km2 and a mean area of 0.07 km2 in 2020. These values are slightly lower than the glacial lake inventory in 2018 (see Section 2) published by Wang [12] but much higher than Wu’s study in 2005 (lake number: 28; lake area: 7.87 km2) [34]. There are 107 glacial lakes (~55% of the total number of glacial lakes) with an area of >0.02 km2 (Figure 2D), which occupy 13.3 km2 (~93%) of the total lake area. Most glacial lakes (14.2 ± 0.3 km2, ~99%) are distributed between 4500 and 5500 m a.s.l. (Figure 2B) and most are located on the east side of the corridor. The total area of the glacial lake on the east side (10.0 ± 0.1 km2, ~69%) is much larger than that of the west (4.4 ± 0.1 km2). The altitudes of the glacial lakes on the west side (5100–5500 m a.s.l.) are generally higher than the east (4500–5300 m a.s.l.). Nevertheless, there are a total of 119 glacial lakes on the western side, accounting for ~61% of the total number of our mapped glacial lakes. However, all of them are less than 0.3 km2 in area.

4.2.2. Glacial Lake Changes

From 1964 to 2020, the total number and area of glacial lakes in the Tama Koshi basin increased considerably (by ~151% in number and ~180% or 3.2% a−1 in area, Figure 2A). Lake expansion between 1980 and 1990 reached the highest rate of 3.5% a−1 (Table 3). Statistically, a total of 12 glacial lakes disappeared during 2000–2020, 36 glacial lakes were newly formed, 10 glacial lakes shrank (7 in the west), and 26 glacial lakes expanded. The glacial lakes that develop between 5000 and 5100 m a.s.l. have experienced the highest expansion rate, with a rate of 4.0% a−1.
In 1964, there were only 27 glacial lakes on the eastern side of the basin, with a total area of 3.1 ± 0.1 km2 (Table 3). This number in 2020 has increased to 77 by a rate of 3.3% a−1 (185%), which is higher than the 2.4% a−1 (133%) on the western side (from 51 in 1964 to 119 in 2020, Table 3). Their overall expansion rate (3.9% a−1) was much higher than the west side (2.2% a−1). According to Table 4, moraine-dammed ice-contacted glacial lakes show a much higher expansion rate (~204%) than non-moraine-dammed ice-contacted lakes (~175%). According to our multi-decade inventories (Table 4), the total number of moraine-dammed ice-contacted lakes was peaked around 1980 at 17, an increase from 10 in 1964. After 1980, the number of moraine-dammed ice-contacted glacial lakes gradually decreased to five in 2020; the other 12 were disconnected from the retreating mother glacier. The annual expansion rates during periods after 1980 (0.4% a−1–2.2% a−1) were also lower than the period before1980 (5.7% a−1, 1964–1980). In contrast, for the non-ice-contacted lakes, their numbers showed a continuous increase from 68 in 1964 to 191 in 2020, while their mean annual area expansion rate (3.1% a−1) is close to those moraine-dammed ice-contacted lakes (3.7% a−1).
Between 1964 and 2020, the lake size distribution has also varied with time (Figure 2D). Glacial lakes with an area of <0.01 km2 and between 0.05 and 0.2 km2 showed the most unstable, while lake sizes between 0.01 and 0.05 km2 gradually increased. For those glacial lakes with larger sizes, e.g., lakes with a size of 0.2–0.5 km2 were relatively stable, and lake size of >1 km2 showing more expansive during 1964–1990 following by relatively stable after 1990.
In this study, we found that two types of glacial lakes have experienced rapid expansion: moraine-dammed ice-contact lake and supraglacial lake (Figure 4 and Figure 5, and Table 4). Based on the results of potentially dangerous glacial lakes evaluation from ICIMOD [21], the actual changes of the glacial lakes, and the past GLOF events (see Section 5.4), we identified eight (seven (B-H) of them are same with ICIMOD’s report) rapidly expanding moraine-dammed ice-contact lakes (Lake A–H, Figure 1 and Figure 6) as PDGLs. Lake G is a moraine-dammed lake formed around in 2009 by the expanding and coalescing with several supraglacial ponds (most of which are base-level lakes) in the 1980s. The total lake area over the entire glacier basin has increased by 9.0% a−1 from 0.08 km2 in 1980 to 1.92 km2 in 2020. Lakes A–C, located on the west side of the basin, expanded rapidly before 1990 but then gradually stabilized in their size. Lake A has shrunk significantly between 2018 and 2020, with a rate of ~−36% (or 0.08 km2 a−1). It was confirmed that the lake burst in 2018 according to the comparison of satellite images (more details are given in Section 5.4 and Figure S2). The expansion rate of Lake D (Niangzongmajue) was maintained at 1.3% a−1 before 2014 and then remained relatively stable, with no evident expansion in the following six years. Lake E (Yalongcuo) has experienced steady growth (1.9% a−1) over the past 56 years. The expansion rate of Lake F increased rapidly after 2004, from 3.8% a−1 between 1990 and 2004 to 5.8% a−1 between 2004 and 2020, with an overall area growth of 52.9%. Lake H (Tsho Rolpa) expanded clearly before 1990 at a rate of ~3.5% a−1, and then its area remained virtually unchanged after the manual intervention (drainage to lower down the lake level) in the 1990s. However, after 2012, despite continued intervention, it still showed a slight expansion trend with a rate of 0.8% a−1. The above results are summarized in Table 5.

5. Discussion

5.1. Glacier and Glacial Lake Variations and Climate Changes

As most HMA lakes are marginally influenced by direct human activities, glacial lake variations are closely related to climate change. Air temperature increase is considered to be the main driving force for glacier retreat and wasting [56]. Previous studies have shown significant warming in the Himalayas. From 1979 to 2019, the temperature rising rate in the Tama Koshi basin was measured at 0.01 °C a−1, and for winter and summer is 0.002 °C a−1 and 0.02 °C a−1, respectively, according to the analysis of ERA5. The precipitation changes during the past decades in the entire Himalayas were not uniform [57]. For the Tama Koshi basin, the basin-scale annual precipitation between 1979 and 2019 showed an increasing trend with a rate of 11.0 mm a−1, and during the winter and summer seasons are 1.8 and 3.4 mm a−1, respectively. This trend was similar to the meteorological observations at a nearby Tingri Station [11]. Although both temperature and humidity will affect evapotranspiration of the lake, the increase in temperature and precipitation has elevated freezing levels in the post-monsoon and early winter periods [58] and may boost the water supply of glacial lakes due to increased glacier ablation and runoff, which cause a negative effect on evapotranspiration [59], and eventually cause the lake level to rise (for example Lake Tsho Rolpa [60]) or induce a GLOF (please see Lake A in Section 5.4).

5.2. Glacier-Lake Interactions

Compared to the land-terminating glaciers, the mean change rates of area, length, elevation, and velocity of lake-terminating glaciers are remarkably higher (Table 2). In the Tama Koshi basin, length changes of most land-terminating glaciers (mostly are debris-covered) were not obvious during the study period (e.g., Glacier I, J, K, on the eastern side and Glacier L, M, N on the western side). The decrease in ice thickness, on the other hand, was their main result of mass loss. This is contrasted with those rapid retreats of lake-terminating glaciers (Figure 3). For example, for those rapidly expanding moraine-dammed lakes (Figure 5), they are all currently connecting with glaciers or have been in contact with ice. Such significant differences in glacial retreating rates were widely observed in the HMA regions [25,61,62] since the lake-calving is suggested as one of the important mass loss processes of the lake-terminating glaciers [63,64]. When those proglacial lakes finished their rapid expansion stage and gradually separated from the mother glacier, their expansion rates become gradually decreased (e.g., Glacier A, B, C). However, the retreat rates of the ice tongue are still faster than those initial land-terminating glaciers.
In addition, glaciers terminating lakes also showed notable changes in their thicknesses and flow velocities. The surface elevations of the glaciers that are connecting with Lakes A, B, E, F, G, and H have shown notably higher lowering rates than those land-terminating debris-covered glaciers (Table 2 and Table 5, e.g., I, J, K). The largest elevation change occurred on the lake-terminating glacier connected with Lake F, reaching −1.2 m a−1 between 1975 and 2016 (Table 5). The remarkable high lowering rate (−0.64 m a−1, Table 5) occurred at the largest debris-covered land-terminating glacier (M) is likely caused by the high-density distribution of supraglacial lakes (most of them are perched lakes) and ice cliffs (Figure S3) on its lower tongue, as these features acting as ablation hot spots on debris-covered glaciers [65,66,67] are unstable and seasonal. Although both lake- and land-terminating glaciers have shown a slowdown of the velocities since 1999, the decelerating rates of lake-terminating glaciers are much lower than the land-terminating glaciers (Table 2), indicating the modulation effect of proglacial calving processes on the maintenance of high ice flow rate of the lake-terminating glaciers [61].

5.3. Evolution and Future Development of Glacial Lakes Controlled by Topography

In general, the topography where the glacial lake develops (the pink rectangle highlighted area in Figure 7) generally is relatively flat. The mean slope of these lake-terminating glaciers is around 5°, while the mother glacier of Lake E has the lowest slope (~3°), and the mother glacier of Lake F has the highest surface gradient but still no more than 8° (Table 5). When the valley between receding glacier fronts and terminal moraines provides a sufficient and achievable area for glacial lake development, e.g., Lake D–H, the lake could continually expand. However, when the glacier gradually retreats to the region with higher elevation and steeper slope (Table 5), the expansion rate of the glacial lake gradually slows down, e.g., lake A–C has been fully loaded (or have reached their maximum upper limit). For Lake F and G, although they have high-gradient ice tongues, they are surviving thanks to their moraine dams supported by adjacent glacier-moraine complexes (Figure 4). From the longitudinal profile analysis (Figure 7), the height of the moraine dam in front of Lake F is the highest (50 m), which offers a strong bearing capacity for lake growth. On the other hand, despite the evident mass loss (with the mean glacier surface elevation change rate of −0.6 m a−1), the land-terminating glaciers (I–N) are still difficult to store water due to the lack of stable dams (Figure 5) or gentle slopes (with the mean slope of 8.4°; Table 5). Supraglacial lakes on these debris-covered glaciers generally are perched lakes thus could only be filled seasonally [54], or englacial or subglacial drainage system in the glacier maybe exist that could drain the lakes periodically [68]. The number and area of supraglacial ponds on Glacier M and N are increasing slightly and may merge into a large supraglacial lake or moraine-dammed lake in the future once those perched lakes are located in the lowest part of the glacier change to base-level lakes as the glacier downwasting. Additionally, some cirque lakes may have formed after small cirque glaciers seriously downwasted or disappeared. We therefore argue that topography is not only an essential factor affecting the expansion of glacial lakes, but also a major condition for the development of glacial lakes. Moreover, the altitude and orientation of the proglacial lake may also affect the thermal condition and water balance of the glacial lake by influencing the absorption of solar radiation and ablation processes of adjacent ice.

5.4. Implications for the Risk of GLOFs

Over 50 GLOF events have been documented in the Himalayas due to glaciers recession [69], 3 of which are located in the Tama Koshi basin (Figure 8, Table 6). The first one is entitled Chubung GLOF (86°27′38″E, 27°52′37″N), located in the Rolwaling Valley in Nepal, in front of Rolpa Tsho (true northwest). This GLOF was triggered by an ice avalanche that fell into the lake on 22 July 1991 and caused the lake to completely drain [38]. The second GLOF event (86°26′49″E, 27°55′46″N), which is named Upper Langbu Tsho, is located in the upper Dagazhuoma River, Rongxer country, Tibet, China, near Lake G (true right). The outburst of this lake from September to November 1992 was possibly triggered by an ice avalanche [69]). The third GLOF event was identified in this study by a continuous satellite tracking of Lake A (86°19′16″E, 27°14′43″N, Figure S2), which is located in the upper Labujikongzangbu River. Remote sensing monitoring showed that Lake A was covered by thick clouds for 13 days (28 July–9 August 2018) before its failure. Meteorological results also showed that the summer precipitation in 2018 was remarkably higher than in other years (Figure 1B, with a total of 1234 mm between June and August). We speculate that the high level of precipitation triggered this GLOF.
With the increase in human activities (the population grew from 34 per km2 in 2000 to 59 per km2 in 2018, Figure 8) and socioeconomic development in both Nepal and Chinese parts of the Tama Koshi basin, it is demanded to construct hydropower facilities to reduce the dependence on biofuels and fossil fuels imports. At present, approximately 10 hydropower projects (the largest hydropower project has an installed capacity of 456 MW, which is also the largest hydropower project of Nepal) in the lower reaches of Tama Koshi permits have been issued, of which 6 are under construction (Figure 8). However, eight rapidly expanding glacial lakes are developed in the basin, many of which are close to the roads, village, and other infrastructures, the straight-line distance between them and the nearest hydropower station is about 15 km. Future climate warming will amplify the consequence and increase the risk of GLOFs hazard [14,40,42,70]. Meanwhile, with the improved living standard and increased public awareness of disaster prevention and mitigation (increased awareness of GLOF hazards) in the Tama Koshi basin, numerous residents have migrated from rural areas (the upper valley) to urban areas, resulting in a significant increase in population density in downstream cities (Figure 8). This avoids some small-scale GLOFs but increases the effect when large-scale GLOFs occur. For example, a recent study [60] has shown that the Tsho Rolpa lake elevation has increased in recent years. If the lake area gradually increases, as we expect (see Section 5.3), the water pressure on the dam body will be greatly increased. Once a large-scale outburst occurs, it will bring disastrous effects to the highly populated areas in Lamobagar.
Remote sensing is a vital tool in quickly monitor the glacial dynamics and hazard risk assessment [71]. While this study has obtained detailed information on glacial lakes area changes, basic data on glacial lakes (e.g., lake basin topography, type of dam construction material, the internal structure of moraine dams, etc.) are still lacking and field investigations are needed. To date, the lake basin parameters and dam characteristics were measured only at Tsho Rolpa in Nepal [39]. Detail investigations and research are relatively rarely conducted in the Chinese Tibet part, where some glacial lakes with rapid expansion and their GLOFs could lead to transnational impacts on downstream Nepal. For those lakes with high risks (e.g., Niangzongmajue and Yanongcuo in China), hazard assessment is crucial, and, if necessary, artificial drainage measurements could be adopted to reduce their GLOF risks. For example, Tsho Rolpa in Nepal has been automatically drained by the construction of an open channel and siphon since 1995, and its outburst risk has been basically controlled, and even keep stable during the Gorkha earthquake in 2015 [72].

6. Conclusions and Remarks

Using 42 moderate- and high-resolution remote sensing images between 1964 and 2020, we completed a comprehensive investigation of glaciers and glacial lakes of the entire Tama Koshi basin, central Himalayas. We mapped a total of 271 glaciers and 196 glacial lakes with coverage of 329.2 ± 1.9 km2 and 14.4 ± 0.3 km2, respectively, in 2020 across the basin. During the past 56 years, the total area of glaciers in the basin has shrunk by 26.2 ± 3.2 km2 (0.13% a−1) and thinned by ~20 m on average, with their mean velocity has gradually slowed down from 5.3 m a−1 in 1999–2003 to 4.0 m a−1 in 2013–2015. The total area of glacial lakes increased by 9.2 ± 0.4 km2 (~180%), with ice-contacted proglacial lakes showing a much higher expansion rate (~204%) than the others. Moreover, large-scale glacial lakes are developed preferentially and experienced rapid expansion on the east side of the basin. For land-terminating debris-covered glaciers, the decrease in ice thickness is the main cause of its mass loss. However, for lake-terminating glaciers, rapid calving retreat of the terminus was the dominant process. Preliminary analysis shows that the increase in temperature and precipitation caused the glacier shrinkage and glacial lake expansion in the Tama Koshi basin. Meanwhile, topography also acted as a vital factor that not only affects the expansion of glacial lakes but also controls the development of glacial lakes. We hypothesize that lake expansion will continue in some cases until the critical local topography of steepening icefall is reached, but the number of lakes may not necessarily increase. Many residents (~58 per km2 in 2018) and hydropower stations (~10 in 2020) are sited in the lower reaches of the basin. The impact of the hazard will be amplified if a transnational GLOF occurs. To better assess and reduce the GLOF risks, further detailed monitoring should be focused on the eight rapidly expanding proglacial lakes in the upper reaches in Chinese Tibet Rongxer, e.g., obtaining the basic information (such as lake bathymetry, lake area/level fluctuation, dam characteristics and downstream high-resolution terrain) of these lakes via field and in situ investigations.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/rs13183614/s1, Figure S1: An example of extracting the rapid expansion glacial lakes on Landsat images using Google Earth Engine Python API, Figure S2: Identified Lake A GLOF event by PlanetScope imagery, Figure S3: Supraglacial lakes and ice cliffs on Glacier M base on PlanetScope image, Table S1: Data used for this study.

Author Contributions

Conceptualization, Q.L.; Data curation, Y.Z.; Formal analysis, Y.Z.; Funding acquisition, Q.L. and Y.W.; Investigation, Y.Z.; Methodology, Y.Z.; Project administration, Q.L. and Y.W.; Visualization, Y.Z.; Writing—original draft, Y.Z.; Writing—review and editing, Y.Z., Q.L., L.S., Y.L., H.W., H.L. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the “Chinese Academy of Sciences Overseas Institutions Platform Project, grant number 131C11KYSB20200033” and the “IAEA Interregional Technical Cooperation Project, grant number INT/5/156”; Our research was also supported by the “National Natural Science Foundation of China Project, grant number 41871069” and the “Sichuan Science and Technology Programs, grant number 2021JDJQ0009 and 2020JDJQ0002”.

Data Availability Statement

The data sets of the glacier and glacial lake in the Tama Koshi basin and the Python code for this study can be found and downloaded from https://doi.org/10.5281/zenodo.5496959.

Acknowledgments

The authors gratefully acknowledge the PlanetLab for the provision of PlanetScope imagery and U.S. Geological Survey for Landsat and Keyhole images, and the authors would also like to thank developers in the GEE community and Qiusheng Wu from the University of Tennessee for sharing very helpful sample codes.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegert, M.J. Role of Glaciers and Ice Sheets in Climate and the Global Water Cycle. Encycl. Hydrol. Sci. 2006, 4, 164. [Google Scholar] [CrossRef]
  2. Kääb, A.; Berthier, E.; Nuth, C.; Gardelle, J.; Arnaud, Y. Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas. Nature 2012, 488, 495–498. [Google Scholar] [CrossRef] [PubMed]
  3. Carton, A.; Baroni, C. The Adamello-Presanella and Brenta Massifs, Central Alps: Contrasting High-Mountain Landscapes and Landforms. In Landscapes and Landforms of Italy; Soldati, M., Marchetti, M., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 101–112. [Google Scholar]
  4. Jansson, P.; Rosqvist, G.; Schneider, T. Glacier fluctuations, suspended sediment flux and glacio-lacustrine sediments. Geogr. Ann. Ser. A Phys. Geogr. 2005, 87, 37–50. [Google Scholar] [CrossRef]
  5. Huss, M.; Bookhagen, B.; Huggel, C.; Jacobsen, D.; Bradley, R.S.; Clague, J.J.; Vuille, M.; Buytaert, W.; Cayan, D.R.; Greenwood, G.; et al. Toward mountains without permanent snow and ice. Earth’s Future 2017, 5, 418–435. [Google Scholar] [CrossRef]
  6. Nie, Y.; Zhang, Y.; Ding, M.; Liu, L.; Wang, Z. Lake change and its implication in the vicinity of Mt. Qomolangma (Everest), central high Himalayas, 1970–2009. Environ. Earth Sci. 2012, 68, 251–265. [Google Scholar] [CrossRef]
  7. Gardner, A.S.; Moholdt, G.; Cogley, J.G.; Wouters, B.; Arendt, A.A.; Wahr, J.; Berthier, E.; Hock, R.; Pfeffer, W.T.; Kaser, G.; et al. A Reconciled Estimate of Glacier Contributions to Sea Level Rise: 2003 to 2009. Science 2013, 340, 852–857. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Wu, G.; Yao, T.; Wang, W.; Zhao, H.; Yang, W.; Zhang, G.; Li, S.; Yu, W.; Lei, Y.; Hu, W. Glacial hazards on Tibetan Plateau and surrounding alpines. Bull. Chin. Acad. Sci. 2019, 34, 1285–1292. [Google Scholar] [CrossRef]
  9. Immerzeel, W.W.; Beek, L.P.H.v.; Bierkens, M.F.P. Climate change will affect the Asian water towers. Science 2010, 328, 1382–1385. [Google Scholar] [CrossRef] [PubMed]
  10. Immerzeel, W.W.; Lutz, A.F.; Andrade, M.; Bahl, A.; Biemans, H.; Bolch, T.; Hyde, S.; Brumby, S.; Davies, B.J.; Elmore, A.C.; et al. Importance and vulnerability of the world’s water towers. Nature 2020, 577, 364–369. [Google Scholar] [CrossRef]
  11. Nie, Y.; Sheng, Y.; Liu, Q.; Liu, L.; Liu, S.; Zhang, Y.; Song, C. A regional-scale assessment of Himalayan glacial lake changes using satellite observations from 1990 to 2015. Remote Sens. Environ. 2017, 189, 1–13. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, X.; Guo, X.; Yang, C.; Liu, Q.; Wei, J.; Zhang, Y.; Liu, S.; Zhang, Y.; Jiang, Z.; Tang, Z. Glacial lake inventory of high-mountain Asia in 1990 and 2018 derived from Landsat images. Earth Syst. Sci. Data 2020, 12, 2169–2182. [Google Scholar] [CrossRef]
  13. Pfeffer, W.T.; Arendt, A.A.; Bliss, A.; Bolch, T.; Cogley, J.G.; Gardner, A.S.; Hagen, J.-O.; Hock, R.; Kaser, G.; Kienholz, C.; et al. The Randolph Glacier Inventory: A globally complete inventory of glaciers. J. Glaciol. 2017, 60, 537–552. [Google Scholar] [CrossRef] [Green Version]
  14. Nie, Y.; Pritchard, H.D.; Liu, Q.; Hennig, T.; Wang, W.; Wang, X.; Liu, S.; Nepal, S.; Samyn, D.; Hewitt, K.; et al. Glacial change and hydrological implications in the Himalaya and Karakoram. Nat. Rev. Earth Environ. 2021, 2, 91–106. [Google Scholar] [CrossRef]
  15. Yao, T.; Yu, W.; Wu, G.; Xu, B.; Yang, W.; Zhao, H.; Wang, W.; Li, S.; Wang, N.; Li, Z.; et al. Glacier anomalies and relevant disaster risks on the Tibetan Plateau and surroundings. Chin. Sci. Bull. 2019, 64, 2770–2782. [Google Scholar]
  16. Bolch, T.; Kulkarni, A.; Kääb, A.; Huggel, C.; Paul, F.; Cogley, J.G.; Frey, H.; Kargel, J.S.; Fujita, K.; Scheel, M.; et al. The State and Fate of Himalayan Glaciers. Science 2012, 336, 310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Brun, F.; Berthier, E.; Wagnon, P.; Kääb, A.; Treichler, D. A spatially resolved estimate of High Mountain Asia glacier mass balances from 2000 to 2016. Nat. Geosci. 2017, 10, 668–673. [Google Scholar] [CrossRef]
  18. Maurer, J.M.; Schaefer, J.M.; Rupper, S.; Corley, A. Acceleration of ice loss across the Himalayas over the past 40 years. Sci. Adv. 2019, 5, eaav7266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Westoby, M.J.; Glasser, N.F.; Brasington, J.; Hambrey, M.J.; Quincey, D.J.; Reynolds, J.M. Modelling outburst floods from moraine-dammed glacial lakes. Earth-Sci. Rev. 2014, 134, 137–159. [Google Scholar] [CrossRef] [Green Version]
  20. Harrison, S.; Kargel, J.S.; Huggel, C.; Reynolds, J.; Shugar, D.H.; Betts, R.A.; Emmer, A.; Glasser, N.; Haritashya, U.K.; Klimeš, J.; et al. Climate change and the global pattern of moraine-dammed glacial lake outburst floods. Cryosphere 2018, 12, 1195–1209. [Google Scholar] [CrossRef] [Green Version]
  21. Bajracharya, S.; Maharjan, S.B.; Shrestha, F.; Sherpa, T.C. Inventroy of Glacial Lakes and Identification of Potentially Dangerous Glacial Lakes in the Koshi, Gandaki, and Karnali River Basins of Nepal, the Tibet Autonomous Region of China, and India; Research Report; ICIMOD, UNDP: Kathmandu, Nepal, 2020. [Google Scholar]
  22. Cook, K.L.; Andermann, C.; Gimbert, F.; Adhikari, B.R.; Hovius, N. Glacial lake outburst floods as drivers of fluvial erosion in the Himalaya. Science 2018, 362, 53. [Google Scholar] [CrossRef] [Green Version]
  23. Wang, W.; Yao, T.; Yang, X. Variations of glacial lakes and glaciers in the Boshula mountain range, southeast Tibet, from the 1970s to 2009. Ann. Glaciol. 2017, 52, 9–17. [Google Scholar] [CrossRef] [Green Version]
  24. Yao, X.; Liu, S.; Sun, M.; Zhang, X. Study on the glacial lake outburst flood events in Tibet since the 20th Century. J. Nat. Resour. 2014, 29, 1377–1390. [Google Scholar]
  25. Liu, Q.; Guo, W.; Nie, Y.; Liu, S.; Xu, J. Recent glacier and glacial lake changes and their interactions in the Bugyai Kangri, southeast Tibet. Ann. Glaciol. 2016, 57, 61–69. [Google Scholar] [CrossRef] [Green Version]
  26. Allen, S.K.; Zhang, G.; Wang, W.; Yao, T.; Bolch, T. Potentially dangerous glacial lakes across the Tibetan Plateau revealed using a large-scale automated assessment approach. Sci. Bull. 2019, 64, 435–445. [Google Scholar] [CrossRef]
  27. Liu, J.; Zhang, J.; Gao, B.; Li, Y.; Mengyu, L.; Wujin, D.; Zhou, L. An overview of glacial lake outburst flood in Tibet, China. J. Glaciol. Geocryol. 2019, 41, 1335–1347. [Google Scholar]
  28. Fan, J.; An, C.; Zhang, X.; Li, X.; Tan, J. Hazard assessment of glacial lake outburst floods in Southeast Tibet based on RS and GIS technologies. Int. J. Remote Sens. 2019, 40, 4955–4979. [Google Scholar] [CrossRef]
  29. Cui, P.; Dang, C.; Cheng, Z.; Scott, K.M. Debris Flows Resulting from Glacial-Lake Outburst Floods in Tibet, China. Phys. Geogr. 2010, 31, 508–527. [Google Scholar] [CrossRef]
  30. Cheng, Z.L.; Liu, J.J.; Liu, J.K. Debris flow induced by glacial lake break in southeast Tibet. WIT Trans. Eng. Sci. 2010, 67, 101–111. [Google Scholar] [CrossRef] [Green Version]
  31. Avian, M.; Bauer, C.; Schlögl, M.; Widhalm, B.; Gutjahr, K.-H.; Paster, M.; Hauer, C.; Frießenbichler, M.; Neureiter, A.; Weyss, G.; et al. The Status of Earth Observation Techniques in Monitoring High Mountain Environments at the Example of Pasterze Glacier, Austria: Data, Methods, Accuracies, Processes, and Scales. Remote Sens. 2020, 12, 1251. [Google Scholar] [CrossRef] [Green Version]
  32. Chen, F.; Zhang, M.; Tian, B.; Li, Z. Extraction of Glacial Lake Outlines in Tibet Plateau Using Landsat 8 Imagery and Google Earth Engine. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 4002–4009. [Google Scholar] [CrossRef]
  33. Haritashya, U.K.; Kargel, J.S.; Shugar, D.H.; Leonard, G.J.; Strattman, K.; Watson, C.S.; Shean, D.; Harrison, S.; Mandli, K.T.; Regmi, D. Evolution and Controls of Large Glacial Lakes in the Nepal Himalaya. Remote Sens. 2018, 10, 798. [Google Scholar] [CrossRef] [Green Version]
  34. Wu, S.; Yao, Z.; Huang, H.; Liu, Z.; Liu, G. Responses of glaciers and glacial lakes to climate variation between 1975 and 2005 in the Rongxer basin of Tibet, China and Nepal. Reg. Environ. Chang. 2012, 12, 887–898. [Google Scholar] [CrossRef]
  35. Rounce, D.; Watson, C.; McKinney, D. Identification of Hazard and Risk for Glacial Lakes in the Nepal Himalaya Using Satellite Imagery from 2000–2015. Remote Sens. 2017, 9, 654. [Google Scholar] [CrossRef] [Green Version]
  36. Shilpakar, R.; Shakya, N.M.; Hiratsuka, A. Impact of Climate Change on Snowmelt Runoff: A Case Study of Tamakoshi Basin in Nepal. In Proceedings of the International Symposium on Society for Social Management Systems, SSMS, Kochi, Japan, 5–7 March 2009. SMS09-124. [Google Scholar]
  37. Khadka, D.; Babel, M.S.; Shrestha, S.; Tripathi, N.K. Climate change impact on glacier and snow melt and runoff in Tamakoshi basin in the Hindu Kush Himalayan (HKH) region. J. Hydrol. 2014, 511, 49–60. [Google Scholar] [CrossRef]
  38. Reynolds, J.M. Glacial hazard assessment at Tsho Rolpa, Rolwaling, Central Nepal. Q. J. Eng. Geol. Hydrogeol. 1999, 32, 209. [Google Scholar] [CrossRef]
  39. Sakai, A.; Chikita, K.; Yamada, T. Expansion of a moraine-dammed glacial lake, Tsho Rolpa, in Rolwaling Himal, Nepal Himalaya. Limnol. Oceanogr. 2000, 45, 1401–1408. [Google Scholar] [CrossRef]
  40. Shrestha, B.B.; Nakagawa, H.; Kawaike, K.; Zhang, H. Glacial and Sediment Hazards in the Rolwaling Valley, Nepal. Int. J. Eros. Control Eng. 2012, 5, 123–133. [Google Scholar] [CrossRef] [Green Version]
  41. Dahal, K.R.; Hagelman, R. People’s risk perception of glacial lake outburst flooding: A case of Tsho Rolpa Lake, Nepal. Environ. Hazards 2011, 10, 154–170. [Google Scholar] [CrossRef]
  42. Shrestha, S.; Bajracharya, A.R.; Babel, M.S. Assessment of risks due to climate change for the Upper Tamakoshi Hydropower Project in Nepal. Clim. Risk Manag. 2016, 14, 27–41. [Google Scholar] [CrossRef] [Green Version]
  43. Khadka, M.; Kayastha, R.B.; Kayastha, R. Future projection of cryospheric and hydrologic regimes in Koshi River basin, Central Himalaya, using coupled glacier dynamics and glacio-hydrological models. J. Glaciol. 2020, 66, 831–845. [Google Scholar] [CrossRef]
  44. Randolph Glacier Inventory Consortium. Randolph Glacier Inventory (RGI)—A Dataset of Global Glacier Outlines: Version 6.0; Technical Report; Global Land Ice Measurements from Space; Digital Media: Boulder, CO, USA, 2017. [Google Scholar] [CrossRef]
  45. Dehecq, A.; Gourmelen, N.; Gardner, A.S.; Brun, F.; Goldberg, D.; Nienow, P.W.; Berthier, E.; Vincent, C.; Wagnon, P.; Trouvé, E. Twenty-first century glacier slowdown driven by mass loss in High Mountain Asia. Nat. Geosci. 2019, 12, 22–27. [Google Scholar] [CrossRef]
  46. Rose, A.N.; McKee, J.J.; Urban, M.L.; Bright, E.A.; Sims, K.M. LandScan 2018; Oak Ridge National Laboratory: Oak Ridge, TN, USA, 2019. [Google Scholar]
  47. Bright, E.A.; Coleman, P.R. LandScan 2000; Oak Ridge National Laboratory: Oak Ridge, TN, USA, 2001. [Google Scholar]
  48. Surazakov, A.; Aizen, V. Positional accuracy evaluation of declassified hexagon KH-9 mapping camera imagery. Photogramm. Eng. Remote Sens. 2010, 76, 603–608. [Google Scholar] [CrossRef] [Green Version]
  49. Goerlich, F.; Bolch, T.; Mukherjee, K.; Pieczonka, T. Glacier mass loss during the 1960’s and 1970’s in the Ak-Shirak range (Kyrgyzstan) from multiple stereoscopic Corona and Hexagon imagery. Remote Sens. 2017, 9, 275. [Google Scholar] [CrossRef] [Green Version]
  50. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  51. Wu, Q. geemap: A Python package for interactive mapping with Google Earth Engine. J. Open Source Softw. 2020, 5, 2305. [Google Scholar] [CrossRef]
  52. Nie, Y.; Liu, Q.; Liu, S. Glacial lake expansion in the central Himalayas by Landsat images, 1990–2010. PLoS ONE 2013, 8, e83973. [Google Scholar] [CrossRef] [Green Version]
  53. Salerno, F.; Thakuri, S.; D’Agata, C.; Smiraglia, C.; Manfredi, E.C.; Viviano, G.; Tartari, G. Glacial lake distribution in the Mount Everest region: Uncertainty of measurement and conditions of formation. Glob. Planet. Chang. 2012, 92, 30–39. [Google Scholar] [CrossRef]
  54. Mertes, J.R.; Thompson, S.S.; Booth, A.D.; Gulley, J.D.; Benn, D.I. A conceptual model of supra-glacial lake formation on debris-covered glaciers based on GPR facies analysis. Earth Surf. Process. Landf. 2017, 42, 903–914. [Google Scholar] [CrossRef] [Green Version]
  55. Hersbach, H.; Bell, B.; Berrisford, P.; Biavati, G.; Horányi, A.; Muñoz Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Rozum, I.; et al. ERA5 Monthly Averaged Data on Single Levels from 1979 to Present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS). 2019. Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/10.24381/cds.f17050d7?tab=overview (accessed on 26 October 2020).
  56. Yao, T.; Thompson, L.; Yang, W.; Yu, W.; Gao, Y.; Guo, X.; Yang, X.; Duan, K.; Zhao, H.; Xu, B.; et al. Different glacier status with atmospheric circulations in Tibetan Plateau and surroundings. Nat. Clim. Chang. 2012, 2, 663–667. [Google Scholar] [CrossRef]
  57. Maussion, F.; Scherer, D.; Mölg, T.; Collier, E.; Curio, J.; Finkelnburg, R. Precipitation Seasonality and Variability over the Tibetan Plateau as Resolved by the High Asia Reanalysis. J. Clim. 2014, 27, 1910–1927. [Google Scholar] [CrossRef] [Green Version]
  58. Pelto, M.; Panday, P.; Matthews, T.; Maurer, J.; Perry, L.B. Observations of Winter Ablation on Glaciers in the Mount Everest Region in 2020–2021. Remote Sens. 2021, 13, 2692. [Google Scholar] [CrossRef]
  59. Song, C.; Huang, B.; Richards, K.; Ke, L.; Hien Phan, V. Accelerated lake expansion on the Tibetan Plateau in the 2000s: Induced by glacial melting or other processes? Water Resour. Res. 2014, 50, 3170–3186. [Google Scholar] [CrossRef] [Green Version]
  60. Thakuri, S.; Chauhan, R.; Baskota, P. Glacial Hazards and Avalanches in High Mountains of Nepal Himalaya. J. Tour. Himalaya Adventures 2020, 2, 87–104. [Google Scholar]
  61. Liu, Q.; Mayer, C.; Wang, X.; Nie, Y.; Wu, K.; Wei, J.; Liu, S. Interannual flow dynamics driven by frontal retreat of a lake-terminating glacier in the Chinese Central Himalaya. Earth Planet. Sci. Lett. 2020, 546, 116450. [Google Scholar] [CrossRef]
  62. Basnett, S.; Kulkarni, A.V.; Bolch, T. The influence of debris cover and glacial lakes on the recession of glaciers in Sikkim Himalaya, India. J. Glaciol. 2013, 59, 1035–1046. [Google Scholar] [CrossRef] [Green Version]
  63. King, O.; Dehecq, A.; Quincey, D.; Carrivick, J. Contrasting geometric and dynamic evolution of lake and land-terminating glaciers in the central Himalaya. Glob. Planet. Chang. 2018, 167, 46–60. [Google Scholar] [CrossRef]
  64. Watson, C.S.; Kargel, J.S.; Shugar, D.H.; Haritashya, U.K.; Schiassi, E.; Furfaro, R. Mass Loss from Calving in Himalayan Proglacial Lakes. Front. Earth Sci. 2020, 7, 342. [Google Scholar] [CrossRef] [Green Version]
  65. Miles, E.S.; Willis, I.; Buri, P.; Steiner, J.F.; Arnold, N.S.; Pellicciotti, F. Surface Pond Energy Absorption Across Four Himalayan Glaciers Accounts for 1/8 of Total Catchment Ice Loss. Geophys. Res. Lett. 2018, 45, 10464–10473. [Google Scholar] [CrossRef] [Green Version]
  66. Steiner, J.F.; Buri, P.; Miles, E.S.; Ragettli, S.; Pellicciotti, F. Supraglacial ice cliffs and ponds on debris-covered glaciers: Spatio-temporal distribution and characteristics. J. Glaciol. 2019, 65, 617–632. [Google Scholar] [CrossRef] [Green Version]
  67. Buri, P.; Miles, E.S.; Steiner, J.F.; Ragettli, S.; Pellicciotti, F. Supraglacial ice cliffs can substantially increase the mass loss of debris-covered glaciers. Geophys. Res. Lett. 2021, 48, e2020GL092150. [Google Scholar] [CrossRef]
  68. Bolch, T.; Buchroithner, M.F.; Peters, J.; Baessler, M.; Bajracharya, S. Identification of glacier motion and potentially dangerous glacial lakes in the Mt. Everest region/Nepal using spaceborne imagery. Nat. Hazards Earth Syst. Sci. 2008, 8, 1329–1340. [Google Scholar] [CrossRef] [Green Version]
  69. Nie, Y.; Liu, Q.; Wang, J.; Zhang, Y.; Sheng, Y.; Liu, S. An inventory of historical glacial lake outburst floods in the Himalayas based on remote sensing observations and geomorphological analysis. Geomorphology 2018, 308, 91–106. [Google Scholar] [CrossRef]
  70. Terink, W.; Immerzeel, W.W.; Lutz, A.F.; Droogers, P.; Khanal, S.; Nepal, S.; Shrestha, A.B. Hydrological and Climate Change Assessment for Hydropower Development in the Tamakoshi River Basin, Nepal. 2017. Available online: https://www.futurewater.nl/wp-content/uploads/2017/05/final_report_WT_v2.pdf (accessed on 30 October 2020).
  71. Quincey, D.J.; Lucas, R.M.; Richardson, S.D.; Glasser, N.F.; Hambrey, M.J.; Reynolds, J.M. Optical remote sensing techniques in high-mountain environments: Application to glacial hazards. Prog. Phys. Geogr. Earth Environ. 2005, 29, 475–505. [Google Scholar] [CrossRef]
  72. Kargel, J.S.; Leonard, G.J.; Shugar, D.H.; Haritashya, U.K.; Bevington, A.; Fielding, E.J.; Fujita, K.; Geertsema, M.; Miles, E.S.; Steiner, J.; et al. Geomorphic and geologic controls of geohazards induced by Nepal’s 2015 Gorkha earthquake. Science 2016, 351, 140–150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Remotesensing 13 03614 g001 550
Figure 1. (A) Distribution of glaciers and glacial lakes in the Tama Koshi basin with the background image of Landsat 8 (24 October 2018), letters indicate detailly investigated glaciers and glacial lakes. The inset (B) shows the mean annual (red), summer (blue), and winter (black) air temperature (left) and precipitation (right) extracted from ERA5. The inset (C) indicates the location of the study area.
Figure 1. (A) Distribution of glaciers and glacial lakes in the Tama Koshi basin with the background image of Landsat 8 (24 October 2018), letters indicate detailly investigated glaciers and glacial lakes. The inset (B) shows the mean annual (red), summer (blue), and winter (black) air temperature (left) and precipitation (right) extracted from ERA5. The inset (C) indicates the location of the study area.
Remotesensing 13 03614 g001
Remotesensing 13 03614 g002 550
Figure 2. Overview of glaciers and glacial lakes. (A) Number and area changes of glaciers and glacial lakes from 1964 to 2020. (B) Elevation distribution of glaciers (in 2020) and glacial lakes (between 1990 and 2020) on the west and east side of the basin. (C) Glacier distribution with aspect in the total (blue), east (red), and west (green) of the basin. (D) Frequency of glacial lakes in various size classes.
Figure 2. Overview of glaciers and glacial lakes. (A) Number and area changes of glaciers and glacial lakes from 1964 to 2020. (B) Elevation distribution of glaciers (in 2020) and glacial lakes (between 1990 and 2020) on the west and east side of the basin. (C) Glacier distribution with aspect in the total (blue), east (red), and west (green) of the basin. (D) Frequency of glacial lakes in various size classes.
Remotesensing 13 03614 g002
Remotesensing 13 03614 g003 550
Figure 3. Comparisons of elevation change rates and velocities of three typical land- and lake-terminating glaciers (Glacier E, I, and L) on the east and west sides of the basin.
Figure 3. Comparisons of elevation change rates and velocities of three typical land- and lake-terminating glaciers (Glacier E, I, and L) on the east and west sides of the basin.
Remotesensing 13 03614 g003
Remotesensing 13 03614 g004 550
Figure 4. Changes of eight (Lake A–H) rapidly expanding proglacial lakes in the past 56 years overlayed on the Esri topographic map.
Figure 4. Changes of eight (Lake A–H) rapidly expanding proglacial lakes in the past 56 years overlayed on the Esri topographic map.
Remotesensing 13 03614 g004
Remotesensing 13 03614 g005 550
Figure 5. The expansion process of eight (Lake A–H) rapidly expanding proglacial lakes with the background of KH and Landsat images.
Figure 5. The expansion process of eight (Lake A–H) rapidly expanding proglacial lakes with the background of KH and Landsat images.
Remotesensing 13 03614 g005
Remotesensing 13 03614 g006 550
Figure 6. Area changes of eight (Lake A–H) rapidly expanding proglacial lakes.
Figure 6. Area changes of eight (Lake A–H) rapidly expanding proglacial lakes.
Remotesensing 13 03614 g006
Remotesensing 13 03614 g007 550
Figure 7. Changes in elevation and velocity of the lake-terminating glaciers (A–H) and land-terminating glaciers (I–N). “E” is elevation, “V” is velocity, the pink rectangle highlighted is the area where the glacial lake develops.
Figure 7. Changes in elevation and velocity of the lake-terminating glaciers (A–H) and land-terminating glaciers (I–N). “E” is elevation, “V” is velocity, the pink rectangle highlighted is the area where the glacial lake develops.
Remotesensing 13 03614 g007
Remotesensing 13 03614 g008 550
Figure 8. GLOF events, hydropower projects, and population tendency (left 2000 and right 2018) in the basin.
Figure 8. GLOF events, hydropower projects, and population tendency (left 2000 and right 2018) in the basin.
Remotesensing 13 03614 g008
Table
Table 1. Glacier inventories between 1964 and 2020 of the Tama Koshi basin.
Table 1. Glacier inventories between 1964 and 2020 of the Tama Koshi basin.
Study Period
(Sensor)		Glacier NumberTotal Area
(km2)		Area Change
(% a−1)		West Side of BasinEast Side of Basin
NumberArea
(km2)		NumberArea
(km2)
1964 (KH-4A)282355.3 ± 1.8-141142.0 ± 0.8141213.4 ± 1.0
1980 (KH-9)283350.4 ± 5.9−0.09143141.2 ± 2.6140209.3 ± 3.4
1990 (Landsat TM)279337.5 ± 19.0−0.37139136.2 ± 8.1140201.3 ± 10.9
2000 (Landsat TM)279335.9 ± 19.0−0.05139135.5 ± 8.1140200.5 ± 10.9
2010 (Landsat ETM+)279335.0 ± 18.9−0.03139135.2 ± 8.1140199.8 ± 10.8
2018 (Landsat OLI)271331.2 ± 18.5−0.14133133.9 ± 7.9138197.4 ± 10.6
2020 (PlanetScope)271329.2 ± 1.9−0.31133133.8 ± 0.8138195.4 ± 1.1
Table
Table 2. Glacier changes between 1964 and 2020 for different glacier types and locations.
Table 2. Glacier changes between 1964 and 2020 for different glacier types and locations.
Study PeriodTotal Land-TerminatingLake-TerminatingWest Side of BasinEast Side of Basin
Area change rate (% a−1)1964–20200.130.130.140.100.15
Retreat rate (m a−1)1964–20203.22.637.72.14.3
Elevation change rate (m a−1)1975–2000−0.32−0.30−0.40−0.26−0.33
2000–2016−0.63−0.62−0.85−0.69−0.62
Mean Velocity (m a−1)1999–20035.34.86.95.45.2
2013–20154.03.46.14.13.8
Note: From 1975 to 2000, glacier-wide mean elevation changes were extracted using the glacier boundary in 1980 and from 2000 to 2016 using the boundary in 2000. From 1999 to 2003, mean velocities were extracted using glacier boundaries in 2000 and from 2013 to 2015 using glacier boundaries in 2010.
Table
Table 3. Glacial lake inventories between 1964 and 2020 and changes in meteorological conditions in the Tama Koshi basin.
Table 3. Glacial lake inventories between 1964 and 2020 and changes in meteorological conditions in the Tama Koshi basin.
Study Period
(Sensor)		Lake NumberTotal Area
(km2)		Area Change
(% a−1)		West Side of BasinEast Side of BasinTemperature Change
(°C a−1)		Precipitation Change
(mm a−1)
NumberArea
(km2)		NumberArea
(km2)
1964 (KH-4A)785.1 ± 0.1-512.0 ± 0.1273.1 ± 0.1--
1980 (KH-9)1157.6 ± 0.53.0763.0 ± 0.6394.5 ± 0.2--
1990 (Landsat TM)17110.3 ± 2.23.61083.8 ± 1.1636.5 ± 1.1--
2000 (Landsat TM)17210.7 ± 2.20.471094.0 ± 1.1636.9 ± 1.10.0330.48
2010 (Landsat ETM+)17511.7 ± 2.30.851123.8 ± 1.1637.8 ± 1.10.0386.4
2018 (Landsat OLI)19814.2 ± 2.72.81184.4 ± 1.3809.8 ± 1.40.000714.6
2020 (PlanetScope)19614.4 ± 0.30.491194.4 ± 0.17710.0 ± 0.1--
Total
Lake growth (dA%)151-180133123185217--
Lake growth rate (% a−1)2.7-3.22.42.23.33.9--
Note: Temperature and precipitation changes were obtained by calculating the difference between ten-year averages.
Table
Table 4. The expansion differences between ice-contact lake and non-ice-contact lake in the Tama Koshi basin.
Table 4. The expansion differences between ice-contact lake and non-ice-contact lake in the Tama Koshi basin.
Study Period
(Sensor)		Moraine-Dammed Ice-Contact LakeNon-Moraine-Dammed Ice-Contact Lake
NumberArea
(km2)		Expansion Rate (% a−1)NumberArea
(km2)		Expansion Rate (% a−1)
1964 (KH-4A)102.1 ± 0.1-683.0 ± 0.1-
1980 (KH-9)173.8 ± 0.25.7983.7 ± 0.31.5
1990 (Landsat TM)114.0 ± 0.40.401606.3 ± 1.86.9
2000 (Landsat TM)124.9 ± 0.62.21605.9 ± 1.7−0.65
2010 (Landsat ETM+)105.3 ± 0.50.831656.4 ± 1.90.85
2018 (Landsat OLI)55.8 ± 0.51.31938.4 ± 2.44.0
2020 (PlanetScope)56.1 ± 0.12.21918.3 ± 0.3−0.71
Total
Lake growth (dA%)−50-204181-175
Lake growth rate (% a−1)−0.9-3.73.2-3.1
Table
Table 5. Rates of changes and other related information for several lake-terminating glaciers (Glacier A-H) and land-terminating glaciers (Glacier I-N).
Table 5. Rates of changes and other related information for several lake-terminating glaciers (Glacier A-H) and land-terminating glaciers (Glacier I-N).
IDName
(Lake or Glacier ID)		(1964–2020)
Rate of Change in Lake Area
(% a−1)		(1999–2015)
Rate of Change in Glacier Velocity
(m a−1) 		(1975–2016)
Change in Glacier Elevation
(m a−1)		(2007)
Average Slope
(°)
ARGI60-15.09690−1.92.1−0.555.5
BRGI60-15.097140.334.8−0.555.8
CRGI60-15.097210.849.5−0.495.2
DNiangzongmajue1.125.4−0.405.5
EYalongcuo1.910.4−0.803.1
FDangpu Lake4.22.8−1.27.7
GRGI60-15.097719.07.4−1.15.5
HTsho Rolpa0.609.7−0.713.0
IShalong Glacier0.045.1−0.647.9
JDangpu Glacier-3.7−0.4710.5
KRGI60-15.03428-2.5−0.717.9
LBamolelingjia Glacier-5.7−0.868.2
MRGI60-15.09675-8.9−0.647.8
NRGI60-15.09729-14.7−0.268.1
Table
Table 6. Historical GLOF events of the basin.
Table 6. Historical GLOF events of the basin.
Name Date of OutburstLocation(Possible) TriggerSource
Chubung 12 July 1991Rolwaling Valley, NepalIce avalanche(Reynolds 1999)
Upper Langbu Tsho22 September 1992~17 November 1992Dagazhuoma River, China Ice avalanche(Nie et al., 2018)
Lake A30 July 2018~9 August 2018Labujikongzangbu River, China Extreme precipitationThis study
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhong, Y.; Liu, Q.; Sapkota, L.; Luo, Y.; Wang, H.; Liao, H.; Wu, Y. Rapid glacier Shrinkage and Glacial Lake Expansion of a China-Nepal Transboundary Catchment in the Central Himalayas, between 1964 and 2020. Remote Sens. 2021, 13, 3614. https://doi.org/10.3390/rs13183614

AMA Style

Zhong Y, Liu Q, Sapkota L, Luo Y, Wang H, Liao H, Wu Y. Rapid glacier Shrinkage and Glacial Lake Expansion of a China-Nepal Transboundary Catchment in the Central Himalayas, between 1964 and 2020. Remote Sensing. 2021; 13(18):3614. https://doi.org/10.3390/rs13183614

Chicago/Turabian Style

Zhong, Yan, Qiao Liu, Liladhar Sapkota, Yunyi Luo, Han Wang, Haijun Liao, and Yanhong Wu. 2021. "Rapid glacier Shrinkage and Glacial Lake Expansion of a China-Nepal Transboundary Catchment in the Central Himalayas, between 1964 and 2020" Remote Sensing 13, no. 18: 3614. https://doi.org/10.3390/rs13183614

APA Style

Zhong, Y., Liu, Q., Sapkota, L., Luo, Y., Wang, H., Liao, H., & Wu, Y. (2021). Rapid glacier Shrinkage and Glacial Lake Expansion of a China-Nepal Transboundary Catchment in the Central Himalayas, between 1964 and 2020. Remote Sensing, 13(18), 3614. https://doi.org/10.3390/rs13183614

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop