Next Article in Journal
Detection of Ships Cruising in the Azimuth Direction Using Spotlight SAR Images with a Deep Learning Method
Next Article in Special Issue
Recent Seasonal Spatiotemporal Variations in Alpine Glacier Surface Elevation in the Pamir
Previous Article in Journal
An Efficient Information-Reinforced Lidar Deep Completion Network without RGB Guided
Previous Article in Special Issue
Large-Scale Monitoring of Glacier Surges by Integrating High-Temporal- and -Spatial-Resolution Satellite Observations: A Case Study in the Karakoram
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of Long-Time Series Variation of Glacial Lakes in Southwestern Tibet: A Case Study in the Nyalam County

1
School of Earth Science, Chengdu University of Technology, Chengdu 610059, China
2
State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China
3
Department of Land Surveying and Geo-Informatics, The Hong Kong Polytechnic University, Hong Kong 999077, China
4
State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu 610065, China
5
College of Hydraulic and Hydroelectric Engineering, Sichuan University, Chengdu 610065, China
6
Northwest Engineering Corporation Limited, Xi’an 710065, China
7
Southwest Branch of China Petroleum Engineering Construction Co., Ltd., Chengdu 610095, China
8
Department Geosciences, College of Sciences and Mathematics, Auburn University, Auburn, AL 36849, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(19), 4688; https://doi.org/10.3390/rs14194688
Submission received: 31 July 2022 / Revised: 13 September 2022 / Accepted: 16 September 2022 / Published: 20 September 2022

Abstract

:
Glacial lakes are important freshwater resources in southern Tibet. However, glacial lake outburst floods have significantly jeopardized the safety of local residents. To better understand the changes in glacial lakes in response to climate change, it is necessary to conduct a long-term evaluation on the areal dynamics of glacial lakes, assisted with local observations. Here, we propose an innovative method of classification and stacking extraction to accurately delineate glacial lakes in southwestern Tibet from 1990 to 2020. Based on Landsat images and meteorological data, we used geographic detectors to detect correlation factors. Multiple regression models were used to analyze the driving factors of the changes in glacier lake area. We combined bathymetric data of the glacial lakes with the changes in climatic variables and utilized HEC-RAS to determine critical circumstances for glacial lake outbursts. The results show that the area of glacial lakes in Nyalam County increased from 27.95 km2 in 1990 to 52.85 km2 in 2020, and eight more glacial lakes were observed in the study area. The glacial lake area expanded by 89.09%, where we found significant growth from 2015 to 2020. The correlation analysis between the glacial lake area and climate change throughout the period shows that temperature and precipitation dominate the expansion of these lakes from 1990 to 2020. We also discover that the progressive increase in water volume of glacial lakes can be attributed to the constant rise in temperature and freeze–thaw of surrounding glaciers. Finally, the critical conditions for the glacial lake’s outburst were predicted by using HEC-RAS combined with the changes in the water volume and climatic factors. It is concluded that GangxiCo endures a maximum water flow of 4.3 × 108 m3, and the glacial lake is in a stable changing stage. This conclusion is consistent with the field investigation and can inform the prediction of glacial lake outbursts in southwestern Tibet in the future.

Graphical Abstract

1. Introduction

Glacial lakes, which are formed by the convergence of meltwater, are essential freshwater resources of local residents [1]. They are also hot spots of natural disasters in the Qinghai–Tibet Plateau [2]. In recent years, the warming of the Qinghai–Tibet Plateau has gradually led to significant changes in the number and size of glacial lakes [3,4,5,6]. Most of the glaciers in southeastern Tibet began to shrink and melt in the mid-19th century [7]. Especially after 2000, the melting of glaciers and the loss of materials have intensified [8,9], resulting in the increased area and frequency of changes in glacial lakes [10,11]. The glaciers in this area account for about 81.6% of China’s total reserves, and it is one of the areas with the largest ice avalanches in China [12]. These glaciers are the source of many rivers [13], such as PumQu, Boqu, Tamba Kosi, and Lapchekhun Khola. Still, they are also essential in glacial lake disasters, including many steep mountains and fault zones, which are more likely to cause ice avalanches [14]. The abnormal changes of glacial lakes inevitably lead to a significant number of natural disasters, such as the outbursts of glacial lakes [15,16]. Some scholars have found that glacial lakes continue to develop [17,18], and more in-depth research is needed referring to the impact of climate change on the changes and eruptions of glacial lakes [19].
To effectively map glacier lakes across a mountain region [20], optical remote sensing imagery shows obvious advancement in estimating ice phenology of lakes and rivers around the globe over in situ measurement methods [21]. Optical sensors such as MODIS (moderate-resolution imaging spectroradiometer) have relatively low resolution and relatively short observation periods. The spatial resolution of Landsat TM, ETM+, and OLI images can reach 30 m, which is more conducive to observing glacial lakes and long-term trend analysis [22]. Most imagery is susceptible to cloud cover due to the Indian Ocean monsoon. However, the collection of no-cloud remote sensing imagery is often susceptible to cloud-covered imagery, further affecting the accuracy of contouring [23]. In southwestern Tibet, most glacial lakes are heavily covered by fragments of rock fragments and surface moraine [24]. The surrounding area of the lakeshore is covered by debris, mainly composed of rock fragments and surface moraines [17,25]. These glacial lakes are difficult to distinguish from hillsides due to the spectral similarity between surface debris and their surroundings [26]. To improve the accuracy of extracting glacial lakes, conventional approaches focus on combining multisource data to enhance recognition features beyond spectral features [27], improving or proposing new classification algorithms (e.g., support vector machines, random forest, and k-nearest neighbor), and relying directly on manual interpretation for complex domains [28]. To improve the accuracy of extracting glacial lake information based on the water body index method and digital elevation model, we proposed a classification and stacking extraction method to extract glacial lake boundaries under different states. The results show that the proposed threshold combined with the classification overlay extraction method can extract glacial lake boundaries in three states: unfrozen, semifrozen, and fully frozen, and can reduce the influence of glacier and hillshade and improve the utilization of optical images. In addition, few studies link a glacial lake’s water level with climate to simulate the outburst of glacial lakes [29]. The characteristics of glacial lakes illustrate significant differences with respect to spatial and temporal scales [30]. How climate change affects the hydrological changes in southwest Tibet requires a comprehensive examination, remaining challenging to the existing research [31]. Currently, the methods of correlation analysis of glacial lakes are mostly Pearson correlation coefficient [32], principal component analysis, geographically weighted regression analysis, etc. Geographic detectors are mainly used to analyze the driving forces and influencing factors of various phenomena and multifactor interactions, which can quantify the strength of the interaction among influencing factors. This paper innovatively uses geographic detectors to study the driving factors of glacial lake changes, combined with multiple linear regression models and the Mann–Kendall trend test and estimation of Sen’s slopes to accurately measure the impact of climate change on glacial lake changes. Considering the potential impact of glacial lakes on the local population in the southwestern Tibetan Plateau [33], it is imperative to study the changes in glacial lakes over the past 30 years [34]. A glacial lake outburst flood (GLOF) is a typical outburst flood caused by the failure of a dam containing a glacial lake, which consists of glacial ice or a terminal moraine [35]. Several GLOF studies have found that glacial retreat was changing water flow patterns, affecting the incidence of glacial lake outburst floods and increasing the risk of flooding and water scarcity associated with future climate change [36]. However, to accurately predict the changes of glacial lakes, it is necessary to predict their distribution characteristics and the underlying laws to understand the dynamic response of climate change to hydrological processes [37,38] and to study the eruption disasters of glacial lakes [39]. In contrast to empirical models, physical models consider the interaction between fluid and particles and flow behaviors such as turbulence during the simulation for glacial lake outburst simulations. Currently, commonly used physical numerical simulation models, including FLDWAN [40], BASEMENT [41], HEC-RAS [42], FLO-2D [43], etc., have been applied to modeling and analyzing glacial lake outburst floods in different regions of the world. This paper uses the HEC-RAS V.5 model to simulate and analyze the outburst of glacial lakes by combining the observation of water volume change information and bathymetric data. At the same time, according to the influence of climatic factors on the glacial lake and the trend of the long sequence change of the glacial lake, the prediction and analysis of the hazard of glacial lake outburst are carried out.
Here, we investigate the spatial and temporal changes in the size and location of several glacial lakes in Nyalam County from 1990 to 2020 using Landsat images and analyze the spatial distribution and characteristics of different types of glacial lakes in the study area from the perspectives of elevation, watershed, and mountain system. Furthermore, the effect of elevation on the changes of the glacial lake was estimated by analyzing the digital elevation model (DEM) and the changing pattern of the glacial lakes. The impacts of climate change and the corresponding drivers are further assessed. The monitoring and evaluation of the outburst hazards of glacial lakes are finally carried out by combining water volume change and bathymetric, can help provide essential data support for scientific assessment and prevention of glacial lake disaster risks.

2. Study Area and Data

2.1. Study Area

The study area is located in the Qinghai–Tibet Plateau in southwestern China (Figure 1). This study selects the upper reaches of Poqu, Nyalam County, Shigatse Region, and Tibet Autonomous Region to understand further the various characteristics of glacial lakes in this region. The coordinate of the study area is 86°00′00″E~86°15′00″E and 28°00′00″N~28°20′00″N. Nyalam County is located in the southwestern part of the Tibet Autonomous Region, near the edge of the Himalayan Mountains and close to the Lhari Gangri Range, with a 2100 km2 surface area. In 2020, 119 glaciers and 366 glacial lakes were recorded in Nyalam County. These glaciers are mainly located in the northeast and south-central regions of the county, with 98.63 km2 [44]. Glacial lakes are primarily concentrated in the northeastern part of Nyalam County, and the junction of the Yala township glacial lake distribution area is relatively large. Most of them are moraine lakes or glacier blockage lakes, such as PacuCo, Chama QudanCo, DareCo, etc.
The geomorphology of the study area presents an alternating pattern of alpine valleys and plateau lake basins. Among them, large glacial lakes are formed in a concentrated distribution near the alpine glaciers, and small glacial trough lakes and ice lakes are scattered and distributed in the valleys. The region is mainly a subarctic semiarid and humid subtropical climate, with a distinct rainy season. There is an excellent possibility that the glacial lake will burst due to continuous rainfall.

2.2. Data Source and Preprocessing

Landsat multispectral images from 1990 to 2020 were used to extract the number and area of glacial lakes. As the glacial lake morphology is relatively stable in autumn and winter, this paper selects the images from September to November each year. This is also the season with low cloud cover and lake ice cover. In addition, less cloudy distribution images were chosen to avoid cloudiness’s influence on the extraction boundary. Other auxiliary data such as DEM, meteorological data, and zoning data were applied to study the distribution of the number and area of glacial lakes and the relationship with meteorological factors (Table 1).

3. Methods

To investigate the spatial and temporal changes of glacial lakes in Nyalam County and their response to climate change, GangxiCo was selected as the study area. This paper used classification and superposition methods to extract glacial lake information automatically. Differences in area and distribution characteristics were sorted and analyzed for 30 years, from 1990 to 2020. The response mechanism of glacial lakes to other factors was then investigated. The multivariate linear simulation was then used to analyze the importance of climate change on the glacial lakes and reveal the changes in the glacial lake. The calculation of water level–flow of the glacial lake employed Manning’s equation. The HEC-RAS was applied to determine the future trends of the glacial lake. The specific technical method proposed in this study was illustrated in the following subsections. Figure 2 shows the flow chart of our methodology.

3.1. Glacial Lake Boundary Extraction

3.1.1. Unfrozen Frozen Lake Boundary Extraction

This paper uses the calibration coefficients in the Landsat satellite header file. First, the radiometric calibration was performed using the software ENVI 5.1. Then, the surface reflectance values of each Landsat image were obtained by multispectral atmospheric correction conducted by the FLAASH model. Finally, this paper combines the normalized difference water index (NDWI) with the threshold method to extract the unfrozen glacial lake boundary. The NDWI uses the ratio of the green band to the NIR band to highlight the image’s water body information and suppress vegetation’s interference information to the maximum extent. Its calculation Equation (1) is as follows:
NDWI = ( BandGreen BandNIR ) / ( BandGreen + BandNIR )
where BandGreen is the reflectance in the green band (Landsat TM, Landsat7 ETM+ band2, Landsat8 OLI band3) and BandNIR is the reflectance value in the NIR band (Landsat TM, Landsat7 ETM+ band4, Landsat8 OLI band5). The whole image was subjected to NDWI operation, and the brightness of different features was divided. Extraction of the unfrozen lake was implemented through a threshold filter. This paper establishes a significant threshold value (0.28) for extraction. In this way, all ground objects except the water of the unfrozen glacial lake were filtered as the background, leaving only the unfrozen glacial lake. Finally, the boundary of the unfrozen part of the glacial lake water body was obtained.

3.1.2. Single-Band Threshold Boundary Extraction

After obtaining the water boundary of the glacial lake in the unfrozen part through boundary extraction, the boundary was obtained based on the spectral characteristics of the glacial lake. Because the difference between the reflectance of the glacial lake in the red band and the NIR band is more significant than a certain threshold, it was usually more remarkable than a specific value in the red band. Thus, Equations (2) and (3) were developed according to these two spectral characteristics of frozen glacial lakes. Considering that frozen glacial lakes have similar spectral characteristics to glaciers, this paper combined the regional attributes of Nyalam County. An experimental plan was set up to magnify the reflectance 10,000 times, highlight it, and then calculate it [45]. Finally, the most suitable threshold for the study area was obtained according to the results extracted from the glacial lake and the Google Earth map. Through continuous analysis of the glacial lakes in Nyalam County, it was found that when a = 1800 and b = 9100, the boundaries of the extracted frozen glacial lakes are relatively complete. Most glacial lakes were formed by the convergence of meltwater produced by the retreat of glaciers and are distributed around the glaciers. The glacier represents strong reflection in the red band, and an appropriate threshold is set in the red light band to eliminate the influence of the glacier in the research area. Using Equation (4), we find that the impact of glaciers in the extraction process of glacial lakes can be effectively reduced when c is 12,000 through continuous analysis of glacial lakes in Nyalam County.
BandRed BandNIR > a
BandRed > b
BandRed < c
where BandRed is the reflectance of the red band (Landsat TM, Landsat7 ETM+ band3, Landsat8 OLI band4) and BandNIR is the reflectance of the NIR band (Landsat TM, Landsat7 ETM+ band4, Landsat8 OLI band5). Since the surrounding environment of the glacial lake is generally relatively gentle, the topographic relief is slight, and the slope is usually tiny, nearly 0. The hillshade is attached to the backside of the mountain ridgelines, the topographic relief is significant, and the topographic slope varies greatly. This paper uses digital elevation model (DEM) data to distinguish lake and water pixels according to the above characteristics.

3.1.3. Precision Inspection

To further test the extraction method’s performance and the new index’s validity, the confusion matrix and overall classification accuracy were calculated, and the manual visual interpretation results of glacial lake boundary extraction were imported into Google Earth for rectification. The confusion matrix is constructed using ENVI 5.1 and calculates overall classification accuracy. The extraction accuracy of typical glacial lakes is evaluated by the confusion matrix’s Kappa coefficient and the misclassification error and omission error. The Kappa coefficient Equation (5) is as follows:
K = P 0 P ε 1 P ε
where P 0 is the sum of the number of correctly classified samples in each category divided by the total number of samples. Suppose the number of actual samples in each category is a 1 , a 2 , a 3 , … a f , the number of predicted samples in each category is b 1 , b 2 , b 3 , … b f , and the total number of samples is n, then there are:
P ε = a 1 × b 1 + a 2 × b 2 + a f × b f n × n

3.2. Response Analysis of Glacial Lake Changes to Climate Change

3.2.1. GeogDetector Model

The GeogDetector model quantifies the strength of various interactions between different ecological factors. The model follows an approach based on the assumption of spatial heterogeneity and is suitable for environmental research fields with high spatial heterogeneity. This paper uses the GeogDetector model to determine the spatial correlation between ecological factors and glacial lake changes. There are four detectors in the GeogDetector model, namely risk, factor, ecological, and interaction [46]. In this study, the differential factor and interaction detectors were used to quantitatively investigate the influence of environmental factors on the changes in the glacial lake. In the process of divergence and factor detection, detecting the spatial variation of Y and detecting how much of the spatial divergence of attribute Y are explained by a specific factor, X. Using the q-value metric, the expression is:
q = 1 h = 1 L N h σ h 2 N σ 2 = 1 SSW SST
SSW = h = 1 L N h σ h 2   ,   SST = N σ 2
where h = 1, …, L is the stratum of variable Y or factor X, that is, the classification or partition; Nh and N are the number of cells in stratum h and the whole area, respectively; σ h 2 and σ 2 are the variance of Y values in stratum h and the entire region, respectively. SSW represents the within sum of squares, and SST represents the total sum of squares. The value of q is in the range [0, 1], and the larger the value, the more pronounced the spatial differentiation of Y. If the independent variable X generates the stratification, then the larger value of q indicates the stronger explanatory power of the independent variable X on the attribute Y, and the weaker the opposite. Interaction detectors were then implemented to assess the combined effects of pairwise glacial lake changes.

3.2.2. Multivariate Linear Regression

This paper uses the main climate factors screened by the GeogDetector model analysis as independent variables. The multiple regression model is used to calculate the trend and annual change of the long-term series data of the glacial lake. The multiple regression model combines the relationship between the response and predictor variables. The general expression form of the multiple linear regression computational model is:
Y = a 0 + a 1 x 1 + a 2 x 2 + a n x n + C
where Y is the expected response value, a 0 is the model intercept, x is the set independent variable, n is the number of independent variables x, and C is the random error after removing the influence of n independent variables on Y.

3.2.3. Mann–Kendall Trend Test and Estimation of Sen’s Slopes

This paper applied the Mann–Kendall test to find significant trends for each climate variable, with separate Mann–Kendall tests for each season or month. The Sen slope, which estimates the median rate of linear change over time, is also calculated to quantify the magnitude of the change. Finally, we verified the significance of lake area trends in southwestern Tibet. To examine the responses of trends to geographic and lake morphology predictors, we used multiple linear regression, using temperature, precipitation, humidity, insolation, and wind speed as model predictions and Sen’s slope values as responses to identify trend patterns in glacial lake variables. In this study, the statistics Z-value, p-value, and Sen slope estimator slope value of the M-K test were calculated by R version 4.1.3.

3.3. Risk Analysis of Glacial Lake Outburst

3.3.1. Water Flow Calculation

The research estimate is based on the field survey of GangxiCo glacial lake on 24 May 2019, to study the reflection of a glacial lake with climate and the glacial lake outburst and the water level–flow relationship of a glacial lake using Manning’s equation to provide water body data for the later evaluation of glacial lake hazard. Moreover, joint glacial lake remote sensing data were used for outburst simulation. For the water level–flow estimates of natural rivers, the Manning Equations (10)–(12) was used as shown below:
Q = A 1 n R 2 3 S 1 2
R = A L
Q = S 1 2 1 n A 5 3 L 2 3
where A is the overwater cross-sectional area, n is the channel roughness, R is the hydraulic radius, S is the slope drop, Q is the flow quantity, and L is the wet perimeter.
The water level–flow equation derived from Manning’s equation is the relationship equation of the downstream river channel of GangxiCo. Combined with remote sensing images and downstream longitudinal section measurement results, the estimated river channel slope drop is S = 0.048. The channel sediments are mainly coarse sand and gravel formed by water scouring after the weathering of granite gneiss [47], referring to the “Hydraulic Calculation Manual” (Second Edition). The roughness ratio of n = 0.04 is comprehensively selected according to the cross-sectional shape of the riverbed and the beach structure [48]. Finally, combined with the “Tibet Autonomous Region Shigatse Glacial Lake Surveying and Mapping Project,” the water level–flow formula of the outlet section and the elevation obtained by the Manning formula is the relationship formula of the channel in the lower reaches of GangxiCo.
When the water surface elevation ≤5181.01 m, the water depth h ≤ 0.15 m, which can be shown as Equations (13)–(15):
A = 57.29 h 2
L = 114.60 h
Q = 197.66 h 3
When the water surface elevation is between 5181.01 m~5185.51 m, that is, 0.15 < h < 4.65 m, it can be concluded that:
A = 1.29 + 0.5 ( 17.19 2 + 20.94 h ) h
L = 17.19 + 21.04 h
Q = 5.48 × ( 1.29 + 17.19 h + 10.47 h 2 ) 5 3 ( 17.19 + 21.04 h ) 2 3

3.3.2. HEC-RAS V.5 Model

The HEC-RAS is one of the most popular open-source models for glacier hazard studies. The glacial lake burst model of HEC-RAS combines a parametric burst model and a physical model. The basic principle is to use a parametric model to predict the parameters of the rupture and introduce a simplified physical dam rupture mechanism, which is used for the development of the crack to the final burst and to solve the process of the downstream propagation of the ruptured flood. The outburst process in this study means that when the flood volume increases to a specific value, it overflows the boundary of the glacial lake, scours the downstream to form an undercut, and gradually expands until the water flow continues to erode down both banks through the breach. This paper applies the new HEC-RAS V.5 model to solve the complete 2D diffuse wave equation given as:
δ τ + p x + q y = 0
p τ + x ( p 2 h ) + y ( p q h ) = n 2 p q p 2 + q 2 h 2 g h δ x + p f + ρ x ( h τ x x ) + ρ y ( h τ x y )
where h is the water depth (m), p and q are the specific flow rates in the x and y directions (m/s), δ is the surface elevation (m), g is the acceleration of gravity (m/s2), n is the Manning resistance, ρ is the water density (kg/m3), τ x x and τ x y are the practical shear stress components, and f is the Coriolis force (s−1). The inertia term of the momentum equation when the diffusion wave is chosen.

4. Results

4.1. Glacial Lake Extraction Results

To verify the extraction accuracy of the glacial lakes, the glacial lakes in three states of unfrozen, semifrozen, and fully frozen in 2020 were extracted as typical glacial lakes for accuracy analysis. The proposed classification extraction superposition algorithm in this paper first needs to set the threshold value for NDWI results to remove unfrozen glacial lakes. Then, a threshold is set for the difference between the red and NIR bands. Next, the frozen lake is extracted, and the red band is used to reduce the influence of glaciers. Finally, the unfrozen lake boundary is superimposed with the frozen lake boundary to obtain the complete glacial lake boundary. Figure 3a–c are remote sensing images, and Figure 3a1–c1 are the glacial lake boundaries by superimposition. In the figure, a is the remote sensing image of the unfrozen glacial lake, a1 is the extraction result of the unfrozen glacial lake, b is the remote sensing image of the not-fully frozen glacial lake, b1 is the extraction result of the fully frozen glacial lake boundary, c is the remote sensing image of the semifrozen glacial lake, and c1 is the extraction result of the unfrozen part of the semifrozen glacial lake.

4.1.1. General Distribution Characteristics of Glacial Lakes

  • Number and area changes
The above method obtained information on glacial lakes in 1990, 1995, 2000, 2005, 2010, 2015, and 2020 (Figure 4a). Statistical analysis was carried out on the research results, and the characteristics of changes in the number and area of glacial lakes from 1990 to 2020 were obtained (Table 2). It was found that during the 30 years from 1990 to 2020, the number of glacial lakes in Nyalam County increased by 7. The area increased by 24.90 km2. Among them, the number growth rate was 0.12%, and the area growth rate was 89.12%. The data indicated that although the area of glacial lakes increased dramatically in the last three decades, the number of glacial lakes has not grown substantially simultaneously. Between 2000 and 2010, the number of lakes decreased, primarily due to the disappearance of small glacial lakes near steep mountains and the merging of broken glacial lakes due to the increase in glacial meltwater. Therefore, while the number decreased, the area continued to increase [49]. After 2010, smaller glaciers were more sensitive to climate change because smaller glaciers reduced their area at a higher rate than more enormous glaciers [50]. According to the current situation, glaciers will further decrease, and the area of glacial lakes will continue to increase [51].
The area time series of six glacial lakes with rapid changes in the area between 1990 and 2020 (Figure 4b), shows continued lake growth. According to the overall analysis of glacial lake changes, the area of lakes in this region continues to grow. At the same time, the number of lakes has also grown to a lesser extent. From the overall situation of the changing trend of glacial lakes, the number and area of glacial lakes are increasing. Among them, the area of glacial lakes grew the fastest from 2015 to 2020, with a growth rate of 0.26 km2·(5a)−1. From 2000 to 2005, the area of glacial lakes grew the slowest, with a growth rate of 0.08 km2·(5a)−1. From 1995 to 2000, the area of glacial lakes increased moderately (More details in Supplementary Materials).
2.
Altitude analysis
Historically, the Nyalam region has experienced four glacial periods, forming moraine landforms and lakes at different elevations. Therefore, the elevation distribution reflects the changes in lakes in different periods (Figure 5). The elevation difference of the lakes in the area is 1240 m, ranging from 4340 m at the lowest to 5580 m at the highest. The study area was divided into six elevation ranges to account for the impact of various Pleistocene glacial periods: <5000 m, 5000 m–5100 m, 5100 m–5200 m, 5200 m–5300 m, 5300 m–5400 m, and <5400 m. Table 3 lists the number, area, area increment, and the number of lakes with or without area increase for each type of lake by elevation range. The distribution of different types of lakes and their elevations are also shown in Table 3. The moraine-dammed lakes are evenly distributed in all elevation ranges, and the cirque lakes are concentrated below 5000 m. U-type valley lakes cover a range of 5000 m to 5300 m, with most distributed between 5100 and 5200 m. Glacial erosion lakes are distributed in two fields, below 5000 m and between 5100 m and 5300 m. Lateral moraine lakes are concentrated between 5000 m and 5100 m. Figure 6 and Figure 7 show the changes in the number and area of lakes at different elevations. As the elevation increases, it is observed that the area forms a single-peak distribution and reaches a maximum between 5000 m and 5400 m. The number of lakes in the area also increased with the elevation.

4.1.2. Characteristics of Glacial Lake Changes in Different Periods

Figure 6 shows that the glacial lakes in Nyalam County have expanded significantly in the past three decades. Among them, Figure 6a–c is the highlight of the local area of GongCuo Lake. The surrounding glaciers have gradually receded with time, the number of glacial lakes has increased, and the area has developed. Most of the increased glacial lakes are located at an elevation of 5000 m, all of which are glacial moraine lakes formed by the retreat of glaciers. At the same time, it can be seen from the changes that glacial retreat is the direct cause of the increase in the number and area of glacial lakes [52].
Based on historical glacial lake outburst data and visual interpretation [53,54], the distribution characteristics and types of glacial lake outbursts in Nyalam County were extracted and analyzed. Among them, Figure 7a,b is the highlight of the local area of GangxiCo and YinreCo Lake. The loss of glacier mass balance can quantify the information of glacial avalanche hazards [55]. The primary trend in glaciers in the Himalayas and the central and southern Tibetan Plateau is overall mass loss [56]. The southern mountain range’s edge was where Nyalam County’s ice avalanche disasters were mainly concentrated. Many glaciers in the region create conditions for glacial avalanche disasters [57]. Glacial lakes more prone to outbursts are glacial moraine lakes that end close to supply glaciers [35]. (More details in Supplementary Materials Table S2 [58])
High-resolution remote sensing images in 2018 and 2019 were extracted (Figure 8a,b). Based on the above remote sensing images of GangxiCo glacier lake and its surroundings in Nyalam County and data from the GF-1 satellite, compared with 2018, the area of glacial lakes in 2019 increased by 35,642 m2. Compared with the 107% expansion (0.34 km2/year) of GangxiCo Lake from 1974–2014 [59], the area of the GangxiCo glacier lake expanded faster. When the storage capacity of the glacial lake was calculated according to the relationship between the area and volume, it was found as V = 0.104A1.42. The water storage capacity of the glacial lake in 2018 was 301,724,326 m3, and in 2019 the water storage capacity of the glacial lake was 305,041,629 m3, an increase of 3,317,303 m3 compared to 2018.

4.1.3. Precision Test Results

The calculation result of the Kappa coefficients is −1~1, but the results are usually between 0 and 1. They can be divided into five groups to indicate different levels of consistency: 0.0~0.20 (very low consistency), 0.21~0.4 (average consistency), 0.41~0.6 (moderate consistency), 0.61~0.80 (high consistency), and 0.81~1 (almost completely consistent). Matlab calculations in this paper show (Table 4) that the classification extraction and superposition algorithm are used to extract the glacial lakes in the three states. The misclassification and omission errors of the extracted glacial lake boundary under the three states are less than 2%. The kappa coefficient is more than 0.81, which indicates that the consistency degree of the glacial lake boundary extracted by the classification extraction and superposition algorithm achieves the best results.

4.2. Analysis Results of Climate Response to Glacial Lake Changes

The expansion of the glacial lake in Nyalam is closely related to external climate change. This paper analyzes the glacial lakes’ evolution characteristics statistically sorted by periods. The average statistics of climate factors are calculated in each decade as a period, and the climate change rate in different periods is calculated by linear fitting based on the annual data. By sorting out and calculating the changes in the glacial lake area and the changes in the corresponding climatic factors, it is found that the GangxiCo in Nyalam County has expanded by a total of 1.88 km. Among them, the area of glacial lakes developed by 0.699 km2, 1.051 km2, and 0.13 km2, respectively, in 1990–2000, 2000–2010, and 2010–2020. The statistics of the growth area and the average value of meteorological factors in the three periods are shown in Table 5. The nearest observed meteorological station in the Nyalam area to the glacial lake, Nyalam station, was selected (85°58′E, 28°11′N).

4.2.1. GeogDetector Model Analysis Results

GeogDetector conducted an environmental climate impact factor analysis on the occurrence (Y) of glacial lake areas in Nyalam County from 1990 to 2020. The climate impact factors or proxy variables (X) include sunshine hours, average temperature, wind speed, relative humidity, and precipitation. Table 6(a) shows the risk factor interaction detection results, where rows 2~5 show the q-values after the interaction of two variables. The results show that the exchange of any two variables on the spatial distribution of the glacial lake is greater than that of the first variable alone, reflecting the results of the risk area detection for a single risk (climate impact factor).
In the first row of the table in Table 6(b), “temperature”, “precipitation”, and “humidity” are the names of each partition of this environmental factor, which are the type quantities. The numerical quantities in the second row are the change in the area within each type zone. The following rows 3~4 are for determining whether there is a statistical difference between the changes of a glacial lake in each type of partition, using the t-test with a significance level of 0.05, and “Y” indicating significant differences between “temperature”, “precipitation” and “humidity”.
Table 6(b) shows the calculation results of q-values for all risk factors. The results indicate that precipitation variables, temperature, and humidity have high q-values, meaning that rivers among these variables are the main environmental factors that determine the glacial lakes’ spatial pattern.

4.2.2. Multivariate Linear Regression Analysis Results

Changes in glacial lakes are closely related to global climate change. Glaciers recede when the weather warms up, supplying high lakes with meltwater as the lakes grow in size [60,61]. Seepage lubricates the glacier’s bottom during severe melting, which may lead to frontal collapse if the glacier extends into the lake, thus triggering lake eruption. The expansion of TP glacial lakes can be explained by the trend of increased precipitation and accelerated melting of glaciers associated with rising temperatures [62]. Many studies have provided substantial evidence for lake and glacier dynamics on the TP under climate fluctuations [63]. Some researchers suggest that increased regional precipitation between 1990 and 2010 led to a significant expansion of enclosed lakes on the TP [64,65]. At the same time, others argue that the expansion of lakes fed by glaciers is primarily caused by rising temperatures that accelerate glacier shrinkage [52,66]. Based on the above research experience, this paper selects these climate-influencing factors for correlation analysis. This paper uses data from 1990 to 2020 from meteorological stations in Nyalam County based on five climate impact factors: annual average precipitation, annual average temperature, annual average sunshine hours, annual average relative humidity, and annual average wind speed. Figure 9 shows the statistical data, which contain the annual average temperature of the Nyalam region and 10-year and 30-year average data, and a linear trend line composed of extremely monthly averages. The figure shows that the temperature does not vary much yearly, and the trend is slowly increasing. The temperature increased by 0.91 °C, as demonstrated by the 30-year average temperature of 4.04 °C, and the temperature rapidly increased from 2005 to 2010, and the average temperature was significantly higher than the overall average temperature. The average temperature shows a general warming trend, with a temperature increase of 1.5 °C and a temperature increase rate of 0.6 °C (10 years), which is higher than the national temperature increase rate ((0.26 ± 0.032 °C)/10 years). Precipitation in different years is significantly different, with an overall upward trend and an increase of 81.17 mm, and a precipitation growth rate of 27.06 mm (10 years). The average annual precipitation varies: 30-year average precipitation is 595.51 mm. The average precipitation from 1990 to 2000 is 551.76 mm, which is 43.75 mm less than the overall observations. After 2000, the value began to fluctuate around the linear trend line.
A tolerance of less than 0.1 indicates severe multicollinearity in the covariance statistics, and a variance expansion factor VIF, which was more significant than 10, indicates severe collinearity. When the VIF value is less than 10, we believe that the data conform to the multiple linear analysis, and there is no multicollinearity among the independent variables. It can be seen from the table that the VIFs of sunshine hours, average temperature, average wind speed, average relative humidity, and precipitation are all less than 10, and their tolerances are all close to 0.8. It can be stated that there is no multicollinearity problem in this model. Analysis of variance (ANOVA) is a significant test invented by R.A. Fisher for the difference in the mean of more than one sample. This paper uses ANOVA to judge whether the regression equation is essential. The corresponding p-value in the ANOVA model, which is the sig in the table, is significant and indicates the size of the difference between the control and the experimental group. In the table, p = 0.011 < 0.05 suggests that the null hypothesis is supported. That is, the linear regression equation is significant. The corresponding T-test in Table 7 is to judge the importance of each variable of the regression equation.
From the above calculation, the optimized regression equation is:
Y = 0.010 x 1 2.462 x 2 6.992 x 3 0.961 x 4 0.006 x 5 + 163.1260
The goodness of fit R2 is 0.689, and the regression coefficients of all independent variables passed the t-test, with an excellent overall fit and no multicollinearity. The criterion of standardized residuals by regression illustrates that this study’s polynomial linear simulation model meets the regression requirements better.
Figure 10 shows that the significance ranking of the variation of a glacial lake in the Nyalam County area is mean relative humidity > mean temperature > sunshine hours > precipitation > mean wind speed. Overall, the average annual climate is mainly influenced by the winter temperature. The yearly average temperature from 1990 to 2020 shows a slowly increasing trend with a climate tendency rate of 0.34 °C (10 years). The annual precipitation shows a weak (insignificant) increasing trend, and the climate trend rate is 5.64 mm/(10 years). The yearly temperature is positively correlated with the change in the glacial lake area, and the annual precipitation is negatively correlated with the difference in the glacial lake area (r = −0.23). However, the glacial lake area change is highly associated with summer precipitation (r = −0.57), indicating that the area varies with the seasons. Precipitation is also one of the important reasons for the change of glacial lakes. Since the humidity is highly correlated with the area change (r = 0.56), the variance of the average relative humidity and the sunshine change obtained by the univariate linear model is slight, and the annual difference is not apparent. Therefore, the changes in Nyalam County’s glacial lakes are mainly affected by temperature and precipitation changes. From 1990 to 2020, the temperature and precipitation in Nyalam County show an increasing trend. Temperature and precipitation are meteorological factors with high correlation. The continuous increase in temperature is the dominant factor in the melting of glaciers, and the melting of glaciers has caused an increase in the area of glacial lakes. Continued increases in temperature cause glaciers to retreat and produce meltwater, increasing the size of glacial lakes. At the same time, it is also affected by the topography and the temperature rise in summer and winter to varying degrees, and the humidity and sunshine hours also change accordingly. Global warming is essential in the growth of glacial lakes in Nyalam County (more details in Supplementary Materials Figure S2).

4.2.3. Mann–Kendall Trend Test and Estimation of Sen’s Slopes Results

As mentioned before, in addition to the correlation analysis, we also estimated the slope of Sen to capture the magnitude of change in each lake with climatic variables. This paper uses the Mann–Kendall test method for significant trend analysis. At the same time, the glacial lake area data from 1990 to 2020 were imported, and the standardized test statistic Z was defined by calculation. When the absolute values of Z were more significant than 1.65, 1.96, and 2.58, the trend passed the significance test with 90%, 95%, and 99% reliability, respectively. After testing, Z = 2.73 > 2.58, which indicates that there is an apparent trend change in the data. In addition, p = 0.007 < 0.05, meaning that the time series data have a trend. In addition, the difference in the Sen slope indicates that the glacial lake area is gradually increasing (Sen’s slope = 4.16). Finally, the Mann–Kendall test examined the importance of trends. All trends related to the glacial lake area were significant, and the Mann–Kendall test shows important trends (p ≤ 0.05) for wind speed, air temperature, and humidity. In addition, the change in Sen’s slope shows that depending on the lake, although the temperature has increased significantly (Sen’s slope = 0.34), the volume of some lakes may decline or remain the same. This difference may imply that the rates of lake area change are not the same under the influence of climatic factors, so it is worth investigating the relationship between changes in lake variables and geography, and other lake characteristics. Furthermore, there are many factors to consider, such as lake morphology, changes in frozen glacier soil, and surrounding glaciers.
To investigate the impact of temperature and glacier-frozen soil on the area change of glacial lakes, this paper sealed the soil samples of the glacial lake and moisture content, wet density, dry density, porosity ratio, saturation, soil particle specific gravity, and particle size analysis on the soil samples. First, the essential engineering–geological parameters of the samples, such as the maximum, minimum, average value, standard deviation, coefficient of variation, etc., were calculated according to the routine analysis of the laboratory test results. Then, according to the test analysis results, the particle size curves of each soil sample were drawn, as shown in Figure 11a–d. After arranging the data in units of test pits, the frequency map of each test pit sample’s particle size distribution, as shown in Figure 11e,f, was obtained. It can be seen that the particle size of glacial lake sediment particles is distributed in multipeaks between 0.005~0.075 mm and 2~10 mm, with poor sorting and extremely poor rounding, which are typical glacial sands.
Since the water level of the glacial lake between 2019 and 2020 was highly similar to the temperature trend, it directly shows an increase in precipitation and a decrease in lake evaporation, indicating that glacial melting was likely the leading cause of glacial lake expansion in the region. Because of the strong correlation between temperature and water volume, most glacial lake sediment particles belong to typical glacial sands. Therefore, this paper infers that the increase in temperature not only melts the glacier but also causes the glacier’s frozen soil to thaw.

4.3. Simulation Results of a Glacial Lake Outburst Event

Based on DEM in the simulation, the physical input parameters of HEC-GeoRAS were obtained through ArcGIS analysis. The range of calibrated Manning resistivity n in the main channel was obtained using the nominal curve values exercise. The simulation results were exported to ArcGIS to depict the water surface in the floodplain. The maximum flow values from the observatory were used for the upstream flood input to simulate a glacial lake outbursts scenario. The measured water surface elevation was 5181.01 m, the contact length between the water crossing section and the water body was L = 17.19 m, and the water depth was h = 0.15 m. Taking these parameters into the Manning formula, the flow rate of GangxiCo is about 0.67 m3/s during this measurement period. Compares the images of the glacial lake before and after the breach (Figure 12). At this time, the maximum storage volume can be calculated as 431,042 × 103 m3. The flood inundation area of HEC-RAS in the simulated river reach (402.53 m) is about 42.87 m2. In the top-left diagram, the line inside the lake is the isobath, and the line outside the lake is the Contour line.
The overall pattern of glacial lake extension over the last three decades is obvious. The GangxiCo lake’s current change rate has been quite gradual, and the estimated value of the lake evolution rate under climate change is 0.52% (0.02 km2·a−1). GangxiCo Lake is now in equilibrium and showing a slightly growing tendency [67,68]. Combined with the stable changes in the dam structure and the nearby glaciers, there is no risk of dyke breach in the GangxiCo glacial lake. However, the glacial lake in Nyalam County is susceptible to changes in climate factors. It is not only affected by highly correlated climatic factors such as precipitation, humidity, and temperature, but the surface water temperature and wind speed are important factors affecting the evaporation of the glacial lake. These factors have different degrees of influence in various areas of Nyalam County. The monitoring and simulation of surface water temperature in Namtso and GangxiCo show that the surface water temperature of the glacial lake increase with the increase in temperature and long-wave radiation. We can infer the corresponding climate factors that cause the changes in glacial lakes through the changes in glacial lake water volume. Hence, future estimates are further analyzed according to the predicted precipitation: from 2016 to 2025, the rate of climate change may continue to the current rate, and the lake area will continue to increase, while from 2026 to 2035, the rate of evolution of the lake may increase as the climate is expected to become moister and the risk of glacial lake outbursts can gradually increase [69].

5. Discussion

The particular geographical environment of the Qinghai–Tibet Plateau leads to the freezing of glacial lakes in most cases. In contrast, most existing methods aim to extract nonfrozen lake boundaries. The frozen lake’s spectral features are similar to those of glaciers [70]. Therefore, this paper uses the classification extraction superposition method to extract the glacial lake boundaries. According to the geographical conditions and glacial lake distribution in Nyalam County, the classification method of classification extraction superposition of glacial lake information is proposed under a stepwise iterative extraction algorithm for remote sensing information. This method overcomes the interference of common unfavorable factors such as glaciers and mountain shadows on glacial lake information extraction. It effectively improves glacial lake information extraction detection speed and quality. The glacial lake is classified as frozen and unfrozen, and the NDWI method is used to extract the unfrozen glacial lake boundary. The frozen glacial lake boundary is extracted using the difference between the frozen surface in the red band and NIR, which are more significant than the appropriate threshold and the red bar is more effective than a specific value. Finally, the unfrozen and frozen glacial lake boundaries are superimposed to obtain the complete glacial lake boundary. Some studies used object-oriented image-processing methods to extract glacial lakes. This paper proposes an innovative analysis method for classification and stacking extraction, which improves the extraction accuracy (Kappa coefficient reaches 90%). This method can perform spatiotemporal analysis on a cloud platform with high speed and efficiency. It can overcome the shortcomings of incomplete information, long processing time, and repetitive tasks and improve the extraction efficiency of long-sequence glacial lakes. At the same time, the classification stack extraction method proposed in this study has the potential to be extended to other debris-covered glacial environment observations. However, since this paper’s extraction method for unfrozen glacial lakes is based on the traditional NDWI method, the extraction of unfrozen glacial lakes is limited mainly by this method, especially in areas with severe snowfall in winter. When snow accumulates on the glacial lake, it is difficult to distinguish the specific boundary of the glacial lake, and the threshold needs to be adjusted according to the regional characteristics. The method is subsequently improved to obtain a better way of extracting the boundaries of glacial lakes.
Previous studies reported that glacial lakes began to expand with climate warming but did not specifically measure the changes and attribute the associated drivers [71]. In this paper, the characteristics of climate change and glacial lakes’ temporal and spatial changes are coupled for analysis. When analyzing the driving factors of glacial lake area changes, it is considered that there may be correlations between influencing factors and the interaction of influencing factors on dependent variables. The geographic detector measured the climatic factors affecting the area change of the glacial lake, and then the multiple linear regression analysis was carried out. The response mechanism of glacial lakes in the Nyalam area to climate change and the differences in response mechanisms between different regions are revealed.
Additionally, a combined analysis of the changes in glacial lake area, water volume, and climate conditions was conducted to inform the future changes. Based on the monitoring of the water level, the flow of the glacial lake in Nyalam County, and the analysis of the corresponding climate, it was found that the lake showed a significant expansion between 1999 and 2010. Because the change in the glacial lake area depends on the balance of water revenue and expenditure [72], it is closely related to climate change [73]. The warm-season temperature increase will replenish glacier and seasonal snow meltwater to the glacial lake and promote the evaporation dissipation of the lake water when the evaporation of the glacial lake is less than the recharge of the glacial lake, making the glacial lake area decrease. The increase in cold-season temperature may cause the ice and snow temperature to increase to some extent, accelerating the melting rate of ice and snow during the warm season. The increase in precipitation also means the recharge of the glacial lake directly or indirectly by snowmelt through slope runoffs. Among them, the water level of the glacial lake from 2019 to 2020 is highly similar to the temperature trend, which means that the temperature change is likely to cause the melting of the glaciers. This melting is the main reason for the region’s expansion of the glacial lake because of the strong correlation between temperature and water volume and the multipeak distribution of particle size in glacial lake sediment. This paper infers that the temperature increase is not only the leading cause of glacial melting but also the thawing of frozen soil. This provides a preliminary research conclusion with specific reference significance for studying glacial lakes in the area. In terms of glacial lakes’ spatiotemporal characteristics and their response to climate, this study mainly used mathematical analysis methods such as linear fitting, numerical comparison, and multivariate joint analysis for analysis and comparison. Methods are complex, and it is also necessary to introduce confidential data to reveal temporal and spatial evolution of the glacial lakes from a deeper dimension.
Recent studies on the hazard/risk assessment of glacial lakes in the Himalayas and some distinct regions have demonstrated a high degree of consistency and consensus in evaluating high-risk glacial lakes [74], even though assessment protocols vary [75,76]. Current research on glacial lakes mainly focuses on the distribution characteristics of contemporary or potential glacial lakes and their accompanying hazards or risks [77]. However, research on prevention and mitigation measures for GLOF from high-hazard/-risk glacial lakes is still lacking [78,79]. This paper provides a new model for evaluating GLOF risk from water volume variability and associated climate impact. In addition, it fills the gap in information about glacial lake outburst data in the Nyalam region using HEC-RAS software. The two-dimensional unsteady flow numerical model proposed in this paper provides a quantitative method for studying the characteristics of glacial lake outburst floods. The primary control equation is based on a complete two-dimensional water balance. The HEC-RAS software applied a numerical solution that captures the topography and automatically handles complex boundaries.
Furthermore, the numerical model has good performance in handling complex topography for simulating glacial lake outburst flood processes [80]. The HEC-RAS uses a parametric model to predict the break parameters. Simultaneously, it introduces a simplified physical dam break mechanism for the development of the break to the final break while solving the process of the break outflow. The model can be applied to the outburst simulation of small and medium glacial lakes by introducing a physical outburst model. This paper analyzes the impact of climate conditions on the GLOF outburst mechanism, combined with glacial lake changes and related climate changes over the past three decades, including the characteristics of the outburst formation process, duration, and flood peak flow. The threshold is set to provide experimental data for the future prevention and control of the area, which is of great significance for the areas at risk of glacial lake outbursts with the increasing medium–low flow. However, this paper mainly monitors and examines factors such as quantity, area, and water volume to analyze glacial lakes’ spatial structure. Considering the lack of analysis of the overall system formed by the glacial lake and nearby glaciers, more methods need to be introduced to analyze the structural stability of the whole glacial lake system.

6. Conclusions

The study of glacial lakes is significant due to their sensitivity and particularity to inform climate change. However, remote sensing of glacial lakes in southern Tibet faces the challenges of covered glacial environments and hillsides. In this study, we proposed an improved estimation method for the area distribution and development of glacial lakes and an innovative approach for extracting glacial lakes by classification and stacking, reducing the impact of hillsides and surrounding glaciers on extraction precision. We also analyzed glacial lake changes and outbursts using long-term imageries that combine meteorological and bathymetric data. Our studies illustrated that: (1) The area of glacial lakes increased sharply from 1990 to 2020, but the number of glacial lakes did not increase simultaneously. The increase in the area of glacial lakes has a unimodal distribution, mainly concentrated in the range of 5100–5400 m above sea level. (2) The continuous increase in temperature leads to the long-term retreat of glaciers and the production of meltwater, which increases the area of glacial lakes and has a far-reaching impact on the growth of glacial lakes in Nyalam County. The study found that glacial lake sediments are typical of glacial sands, and rising temperatures help melt glaciers, which is the primary reason for the thawing of frozen soil. (3) Through mathematical simulation, the maximum water volume that the GangxiCo can withstand is 431,042 × 103 m3. The introduction of the classification extraction stacking algorithm and the outburst analysis of the HEC-RAS model provides a new perspective for remonitoring the environmental dynamics of glacial lakes in climate change. This method may be extended to the outburst simulation of other glacial lakes. The future trend can be projected by comparing the increase and maximum water volume of glacial lakes with glacial meltwater and precipitation in the future, providing accurate and reasonable data for regional flood control and disaster reduction.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/rs14194688/s1, Figure S1: Glacial lake area change; Figure S2: Standard P-P plot of regression standardized residuals; Table S1: Glacial lake change Information; Table S2: Since the 20th century, the outburst of the Glacier Lake in Nyalam County, Tibet, China; Table S3: The Confusion Matrix for the extraction accuracy of Glacial Lake.

Author Contributions

Conceptualization, X.D.; data curation, G.Q., J.C., W.L. (Weile Li) and Z.Y.; funding acquisition, X.D. and W.L. (Weile Li); methodology, G.Q., X.D. and J.C.; resources, G.Q., X.D., W.L. (Weile Li), Z.Y. and Y.S.; software, G.Q. and W.L. (Wenxin Liu); validation, G.Q., W.L. (Wenxin Liu) and Y.S.; visualization, G.Q., X.D., Y.S. and J.R.; writing—original draft, G.Q., X.D. and B.Z.; writing—review and editing, G.Q., J.C., M.W., W.L. (Wenxin Liu), Z.Y., H.L., Y.W. and M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China, grant No. 2021YFC3000401; the Research Center for Human Geography of Tibetan Plateau and Its Eastern Slope (Chengdu University of Technology), grant No. RWDL2021-ZD003; Key Research Bases of Humanities and Social Sciences in Higher Education in Sichuan Province, Sichuan Center for Disaster Economic Research, grant No. ZHJJ2021-ZD001 and the National Key Research and Development Program of China under Grant 2017YFB0503601.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge any support given which is not covered by the author’s contributions or funding sections. Includes administrative and technical support that is not covered, as well as data materials for experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pritchard, H.D. Asia’s shrinking glaciers protect large populations from drought stress. Nature 2019, 569, 649–654. [Google Scholar] [CrossRef] [PubMed]
  2. Cook, S.J.; Quincey, D.J. Estimating the volume of Alpine glacial lakes. Earth Surf. Dyn. 2015, 3, 559–575. [Google Scholar] [CrossRef]
  3. Bookhagen, B.; Burbank, D.W. Toward a complete Himalayan hydrological budget: Spatiotemporal distribution of snowmelt and rainfall and their impact on river discharge. J. Geophys. Res. Earth Surf. 2010, 115. [Google Scholar] [CrossRef]
  4. Kraaijenbrink, P.D.; Bierkens, M.; Lutz, A.; Immerzeel, W. Impact of a global temperature rise of 1.5 degrees Celsius on Asia’s glaciers. Nature 2017, 549, 257–260. [Google Scholar] [CrossRef] [PubMed]
  5. Immerzeel, W.W.; Van Beek, L.; Konz, M.; Shrestha, A.; Bierkens, M. Hydrological response to climate change in a glacierized catchment in the Himalayas. Clim. Chang. 2012, 110, 721–736. [Google Scholar] [CrossRef]
  6. Kaser, G.; Großhauser, M.; Marzeion, B. Contribution potential of glaciers to water availability in different climate regimes. Proc. Natl. Acad. Sci. USA 2010, 107, 20223–20227. [Google Scholar] [CrossRef]
  7. Zalazar, L.; Ferri, L.; Castro, M.; Gargantini, H.; Gimenez, M.; Pitte, P.; Ruiz, L.; Masiokas, M.; Costa, G.; Villalba, R. Spatial distribution and characteristics of Andean ice masses in Argentina: Results from the first National Glacier Inventory. J. Glaciol. 2020, 66, 938–949. [Google Scholar] [CrossRef]
  8. Rounce, D.R.; Hock, R.; Shean, D.E. Glacier mass change in High Mountain Asia through 2100 using the open-source python glacier evolution model (PyGEM). Front. Earth Sci. 2020, 7, 331. [Google Scholar] [CrossRef]
  9. Bolch, T.; Kulkarni, A.; Kääb, A.; Huggel, C.; Paul, F.; Cogley, J.G.; Frey, H.; Kargel, J.S.; Fujita, K.; Scheel, M. The state and fate of Himalayan glaciers. Science 2012, 336, 310–314. [Google Scholar] [CrossRef]
  10. 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]
  11. Gardelle, J.; Arnaud, Y.; Berthier, E. Contrasted evolution of glacial lakes along the Hindu Kush Himalaya mountain range between 1990 and 2009. Glob. Planet. Chang. 2011, 75, 47–55. [Google Scholar] [CrossRef]
  12. Liu, C.; Shi, Y.; Wang, Z.; Xie, Z. Glacier resources and their distributive characteristics in China, a review on Chinese Glacier Inventory. J. Glaciol. Geocryol. 2000, 22, 106–112. [Google Scholar]
  13. Immerzeel, W.W.; Lutz, A.; Andrade, M.; Bahl, A.; Biemans, H.; Bolch, T.; Hyde, S.; Brumby, S.; Davies, B.; Elmore, A. Importance and vulnerability of the world’s water towers. Nature 2020, 577, 364–369. [Google Scholar] [CrossRef]
  14. Racoviteanu, A.E.; Armstrong, R.; Williams, M.W. Evaluation of an ice ablation model to estimate the contribution of melting glacier ice to annual discharge in the Nepal Himalaya. Water Resour. Res. 2013, 49, 5117–5133. [Google Scholar] [CrossRef]
  15. Carrivick, J.L.; Tweed, F.S. A global assessment of the societal impacts of glacier outburst floods. Glob. Planet. Chang. 2016, 144, 1–16. [Google Scholar] [CrossRef]
  16. Schwanghart, W.; Worni, R.; Huggel, C.; Stoffel, M.; Korup, O. Uncertainty in the Himalayan energy–water nexus: Estimating regional exposure to glacial lake outburst floods. Environ. Res. Lett. 2016, 11, 074005. [Google Scholar] [CrossRef]
  17. Scherler, D.; Bookhagen, B.; Strecker, M.R. Spatially variable response of Himalayan glaciers to climate change affected by debris cover. Nat. Geosci. 2011, 4, 156–159. [Google Scholar] [CrossRef]
  18. Tsutaki, S.; Fujita, K.; Nuimura, T.; Sakai, A.; Sugiyama, S.; Komori, J.; Tshering, P. Contrasting thinning patterns between lake-and land-terminating glaciers in the Bhutanese Himalaya. Cryosphere 2019, 13, 2733–2750. [Google Scholar] [CrossRef]
  19. Ragettli, S.; Immerzeel, W.W.; Pellicciotti, F. Contrasting climate change impact on river flows from high-altitude catchments in the Himalayan and Andes Mountains. Proc. Natl. Acad. Sci. USA 2016, 113, 9222–9227. [Google Scholar] [CrossRef] [PubMed]
  20. Huggel, C.; Kääb, A.; Haeberli, W.; Teysseire, P.; Paul, F. Remote Sensing Based Assessment of Hazards from Glacier Lake Outbursts: A Case Study in the Swiss Alps; NRC Research Press: Ottawa, ON, Canada, 2002; p. 39. [Google Scholar] [CrossRef]
  21. Schubert, A.; Faes, A.; Kääb, A.; Meier, E. Glacier surface velocity estimation using repeat TerraSAR-X images: Wavelet-vs. correlation-based image matching. ISPRS J. Photogramm. Remote Sens. 2013, 82, 49–62. [Google Scholar] [CrossRef]
  22. Quincey, D.; Lucas, R.; Richardson, S.; Glasser, N.; Hambrey, M.; Reynolds, J. Optical remote sensing techniques in high-mountain environments: Application to glacial hazards. Prog. Phys. Geogr. 2005, 29, 475–505. [Google Scholar] [CrossRef]
  23. Pepin, N.; Bradley, R.S.; Diaz, H.; Baraër, M.; Caceres, E.; Forsythe, N.; Fowler, H.; Greenwood, G.; Hashmi, M.; Liu, X. Elevation-dependent warming in mountain regions of the world. Nat. Clim. Chang. 2015, 5, 424–430. [Google Scholar] [CrossRef]
  24. Sakai, A.; Fujita, K. Formation conditions of supraglacial lakes on debris-covered glaciers in the Himalaya. J. Glaciol. 2010, 56, 177–181. [Google Scholar] [CrossRef]
  25. Sakai, A.; Nishimura, K.; Kadota, T.; Takeuchi, N. Onset of calving at supraglacial lakes on debris-covered glaciers of the Nepal Himalaya. J. Glaciol. 2009, 55, 909–917. [Google Scholar] [CrossRef]
  26. Paul, F.; Bolch, T.; Kääb, A.; Nagler, T.; Nuth, C.; Scharrer, K.; Shepherd, A.; Strozzi, T.; Ticconi, F.; Bhambri, R. The glaciers climate change initiative: Methods for creating glacier area, elevation change and velocity products. Remote Sens. Environ. 2015, 162, 408–426. [Google Scholar] [CrossRef]
  27. Lippl, S.; Vijay, S.; Braun, M. Automatic delineation of debris-covered glaciers using InSAR coherence derived from X-, C-and L-band radar data: A case study of Yazgyl Glacier. J. Glaciol. 2018, 64, 811–821. [Google Scholar] [CrossRef]
  28. Du, P.; Bai, X.; Tan, K.; Xue, Z.; Samat, A.; Xia, J.; Li, E.; Su, H.; Liu, W. Advances of four machine learning methods for spatial data handling: A review. J. Geovisualization Spat. Anal. 2020, 4, 13. [Google Scholar] [CrossRef]
  29. Jing, L.J.; Lan, C.Z.; Cheng, S.P. The relationship between air temperature fluctuation and Glacial Lake Outburst Floods in Tibet, China. Quat. Int. 2014, 321, 78–87. [Google Scholar] [CrossRef]
  30. Rounce, D.R.; McKinney, D.C.; Lala, J.M.; Byers, A.C.; Watson, C.S. A new remote hazard and risk assessment framework for glacial lakes in the Nepal Himalaya. Hydrol. Earth Syst. Sci. 2016, 20, 3455–3475. [Google Scholar] [CrossRef]
  31. 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]
  32. Narama, C.; Daiyrov, M.; Tadono, T.; Yamamoto, M.; Kääb, A.; Morita, R.; Ukita, J. Seasonal drainage of supraglacial lakes on debris-covered glaciers in the Tien Shan Mountains, Central Asia. Geomorphology 2017, 286, 133–142. [Google Scholar] [CrossRef]
  33. Hewitt, K. Glaciers of the Karakoram Himalaya: Glacial Environments, Processes, Hazards and Resources; Springer Science & Business Media: Berlin, Germany, 2013. [Google Scholar]
  34. Anand, J.; Gosain, A.; Khosa, R.; Srinivasan, R. Regional scale hydrologic modeling for prediction of water balance, analysis of trends in streamflow and variations in streamflow: The case study of the Ganga River basin. J. Hydrol. Reg. Stud. 2018, 16, 32–53. [Google Scholar] [CrossRef]
  35. Richardson, S.D.; Reynolds, J.M. An overview of glacial hazards in the Himalayas. Quat. Int. 2000, 65–66, 31–47. [Google Scholar] [CrossRef]
  36. Nie, Y.; Pritchard, H.D.; Liu, Q.; Hennig, T.; Wang, W.; Wang, X.; Liu, S.; Nepal, S.; Samyn, D.; Hewitt, K. Glacial change and hydrological implications in the Himalaya and Karakoram. Nat. Rev. Earth Environ. 2021, 2, 91–106. [Google Scholar] [CrossRef]
  37. Masood, M.; Yeh, P.-F.; Hanasaki, N.; Takeuchi, K. Model study of the impacts of future climate change on the hydrology of Ganges–Brahmaputra–Meghna basin. Hydrol. Earth Syst. Sci. 2015, 19, 747–770. [Google Scholar] [CrossRef]
  38. Su, F.; Zhang, L.; Ou, T.; Chen, D.; Yao, T.; Tong, K.; Qi, Y. Hydrological response to future climate changes for the major upstream river basins in the Tibetan Plateau. Glob. Planet. Chang. 2016, 136, 82–95. [Google Scholar] [CrossRef]
  39. Westoby, M.J.; Glasser, N.F.; Brasington, J.; Hambrey, M.J.; Quincey, D.; Reynolds, J. Modelling outburst floods from moraine-dammed glacial lakes. Earth-Sci. Rev. 2014, 134, 137–159. [Google Scholar] [CrossRef]
  40. Bajracharya, B.; Shrestha, A.B.; Rajbhandari, L. Glacial lake outburst floods in the Sagarmatha region. Mt. Res. Dev. 2007, 27, 336–344. [Google Scholar] [CrossRef]
  41. Worni, R.; Stoffel, M.; Huggel, C.; Volz, C.; Casteller, A.; Luckman, B. Analysis and dynamic modeling of a moraine failure and glacier lake outburst flood at Ventisquero Negro, Patagonian Andes (Argentina). J. Hydrol. 2012, 444, 134–145. [Google Scholar] [CrossRef]
  42. Wang, W.; Gao, Y.; Anacona, P.I.; Lei, Y.; Xiang, Y.; Zhang, G.; Li, S.; Lu, A. Integrated hazard assessment of Cirenmaco glacial lake in Zhangzangbo valley, Central Himalayas. Geomorphology 2018, 306, 292–305. [Google Scholar] [CrossRef]
  43. Petrakov, D.; Tutubalina, O.; Aleinikov, A.; Chernomorets, S.; Evans, S.; Kidyaeva, V.; Krylenko, I.; Norin, S.; Shakhmina, M.; Seynova, I. Monitoring of Bashkara Glacier lakes (Central Caucasus, Russia) and modelling of their potential outburst. Nat. Hazards 2012, 61, 1293–1316. [Google Scholar] [CrossRef]
  44. Liu, Z.; Su, Z. Reasources and distrubution of glaciers on the Tibetan Plateau. Resour. Sci. 2000, 22, 49–52. [Google Scholar] [CrossRef]
  45. Han, D.; Li, C. Based on Landsat8 OLI Remote Sensing Image Glacial Lake Boundary Classification and Stacking Extraction; Beijing Surveying and Mapping: Beijing, China, 2020. [Google Scholar] [CrossRef]
  46. Carling, P.; Villanueva, I.; Herget, J.; Wright, N.; Borodavko, P.; Morvan, H. Unsteady 1D and 2D hydraulic models with ice dam break for Quaternary megaflood, Altai Mountains, southern Siberia. Glob. Planet. Chang. 2010, 70, 24–34. [Google Scholar] [CrossRef]
  47. Hazra, P.; Krishna, A.P. Assessment of Proglacial Lakes in Sikkim Himalaya, India for Glacial Lake Outburst Flood (GLOF) Risk Analysis using HEC-RAS and Geospatial Techniques. J. Geol. Soc. India 2022, 98, 344–352. [Google Scholar] [CrossRef]
  48. Li, Y.; Liu, C.; Zhang, D.; Liang, K.; Li, X.; Dong, G. Reduced runoff due to anthropogenic intervention in the Loess Plateau, China. Water 2016, 8, 458. [Google Scholar] [CrossRef]
  49. Jain, S.K.; Mir, R.A. Glacier and glacial lake classification for change detection studies using satellite data: A case study from Baspa basin, western Himalaya. Geocarto Int. 2019, 34, 391–414. [Google Scholar] [CrossRef]
  50. 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]
  51. Li, D.-S.; Cui, B.-L.; Wang, Y.; Xiao, B.; Jiang, B.-F. Glacier extent changes and possible causes in the Hala Lake Basin of Qinghai-Tibet Plateau. J. Mt. Sci. 2019, 16, 1571–1583. [Google Scholar] [CrossRef]
  52. Yao, T.; Pu, J.; Lu, A.; Wang, Y.; Yu, W. Recent glacial retreat and its impact on hydrological processes on the Tibetan Plateau, China, and surrounding regions. Arct. Antarct. Alp. Res. 2007, 39, 642–650. [Google Scholar] [CrossRef] [Green Version]
  53. Yin, Z.-Y.; Zhang, X.; Liu, X.; Colella, M.; Chen, X. An assessment of the biases of satellite rainfall estimates over the Tibetan Plateau and correction methods based on topographic analysis. J. Hydrometeorol. 2008, 9, 301–326. [Google Scholar] [CrossRef]
  54. Wang, S.-J.; Zhang, T. Glacial lakes change and current status in the central Chinese Himalayas from 1990 to 2010. J. Appl. Remote Sens. 2013, 7, 073459. [Google Scholar] [CrossRef]
  55. Bazai, N.A.; Cui, P.; Carling, P.A.; Wang, H.; Hassan, J.; Liu, D.; Zhang, G.; Jin, W. Increasing glacial lake outburst flood hazard in response to surge glaciers in the Karakoram. Earth-Sci. Rev. 2021, 212, 103432. [Google Scholar] [CrossRef]
  56. Che, Y.; Zhang, M.; Li, Z.; Li, H.; Wang, S.; Sun, M.; Zha, S. Glacier mass-balance and length variation observed in China during the periods 1959–2015 and 1930–2014. Quat. Int. 2017, 454, 68–84. [Google Scholar] [CrossRef]
  57. Jain, S.K.; Sinha, R.K.; Chaudhary, A.; Shukla, S. Expansion of a glacial lake, Tsho Chubda, Chamkhar Chu Basin, Hindukush Himalaya, Bhutan. Nat. Hazards 2015, 75, 1451–1464. [Google Scholar] [CrossRef]
  58. Yao, X.; Liu, S.; Sun, M.; Zhang, X. Study on the glacial lake outburst flood events in Tibet since the 20th century. Ziran Ziyuan Xuebao 2014, 29, 1377–1390. [Google Scholar] [CrossRef]
  59. Zhang, G.; Bolch, T.; Allen, S.; Linsbauer, A.; Chen, W.; Wang, W. Glacial lake evolution and glacier–lake interactions in the Poiqu River basin, central Himalaya, 1964–2017. J. Glaciol. 2019, 65, 347–365. [Google Scholar] [CrossRef]
  60. Nepal, S. Impacts of climate change on the hydrological regime of the Koshi river basin in the Himalayan region. J. Hydro-Environ. Res. 2016, 10, 76–89. [Google Scholar] [CrossRef]
  61. Qian, H.; Xiao’ai, D.; Shiqi, C. Assessing the effects of vegetation and precipitation on soil erosion in the Three-River Headwaters Region of the Qinghai-Tibet Plateau, China. J. Arid. Land 2020, 12, 865–886. [Google Scholar] [CrossRef]
  62. Lei, Y.; Yang, K.; Wang, B.; Sheng, Y.; Bird, B.W.; Zhang, G.; Tian, L. Response of inland lake dynamics over the Tibetan Plateau to climate change. Clim. Chang. 2014, 125, 281–290. [Google Scholar] [CrossRef]
  63. Sun, J.; Zhou, T.; Liu, M.; Chen, Y.; Shang, H.; Zhu, L.; Shedayi, A.A.; Yu, H.; Cheng, G.; Liu, G. Linkages of the dynamics of glaciers and lakes with the climate elements over the Tibetan Plateau. Earth-Sci. Rev. 2018, 185, 308–324. [Google Scholar] [CrossRef]
  64. Huang, L.; Liu, J.; Shao, Q.; Liu, R. Changing inland lakes responding to climate warming in Northeastern Tibetan Plateau. Clim. Chang. 2011, 109, 479–502. [Google Scholar] [CrossRef]
  65. Yang, Z.; Dai, X.; Wang, Z.; Gao, X.; Qu, G.; Li, W.; Li, J.; Lu, H.; Wang, Y. The dynamics of Paiku Co lake area in response to climate change. J. Water Clim. Chang. 2022, 13, 2725–2746. [Google Scholar] [CrossRef]
  66. Bian, D.; Bian, B.; La, B.; Wang, C.; Chen, T. The response of water level of Selin Co to climate change during 1975–2008. Acta Geogr. Sin. 2010, 65, 313–319. [Google Scholar] [CrossRef]
  67. Li, W.; Zhao, B.; Xu, Q.; Scaringi, G.; Lu, H.; Huang, R. More frequent glacier-rock avalanches in Sedongpu gully are blocking the Yarlung Zangbo River in eastern Tibet. Landslides 2022, 19, 589–601. [Google Scholar] [CrossRef]
  68. Su, P.; Liu, J.; Li, Y.; Liu, W.; Wang, Y.; Ma, C.; Li, Q. Changes in glacial lakes in the Poiqu River basin in the central Himalayas. Hydrol. Earth Syst. Sci. 2021, 25, 5879–5903. [Google Scholar] [CrossRef]
  69. Radić, V.; Bliss, A.; Beedlow, A.C.; Hock, R.; Miles, E.; Cogley, J.G. Regional and global projections of twenty-first century glacier mass changes in response to climate scenarios from global climate models. Clim. Dyn. 2014, 42, 37–58. [Google Scholar] [CrossRef]
  70. Jiang, S.; Nie, Y.; Liu, Q.; Wang, J.; Liu, L.; Hassan, J.; Liu, X.; Xu, X. Glacier change, supraglacial debris expansion and glacial lake evolution in the Gyirong river basin, central Himalayas, between 1988 and 2015. Remote Sens. 2018, 10, 986. [Google Scholar] [CrossRef]
  71. Sorg, A.; Bolch, T.; Stoffel, M.; Solomina, O.; Beniston, M. Climate change impacts on glaciers and runoff in Tien Shan (Central Asia). Nat. Clim. Chang. 2012, 2, 725–773. [Google Scholar] [CrossRef]
  72. Hirabayashi, Y.; Kanae, S.; Emori, S.; Oki, T.; Kimoto, M. Global projections of changing risks of floods and droughts in a changing climate. Hydrol. Sci. J. 2008, 53, 754–772. [Google Scholar] [CrossRef]
  73. Harrison, S.; Kargel, J.S.; Huggel, C.; Reynolds, J.; Shugar, D.H.; Betts, R.A.; Emmer, A.; Glasser, N.; Haritashya, U.K.; Klimeš, J. Climate change and the global pattern of moraine-dammed glacial lake outburst floods. Cryosphere 2018, 12, 1195–1209. [Google Scholar] [CrossRef]
  74. Gurung, D.R.; Khanal, N.R.; Bajracharya, S.R.; Tsering, K.; Joshi, S.; Tshering, P.; Chhetri, L.K.; Lotay, Y.; Penjor, T. Lemthang Tsho glacial Lake outburst flood (GLOF) in Bhutan: Cause and impact. Geoenviron. Disasters 2017, 4, 17. [Google Scholar] [CrossRef]
  75. Wang, S.; Che, Y.; Xinggang, M. Integrated risk assessment of glacier lake outburst flood (GLOF) disaster over the Qinghai–Tibetan Plateau (QTP). Landslides 2020, 17, 2849–2863. [Google Scholar] [CrossRef]
  76. Zhang, T.; Wang, W.; Gao, T.; An, B.; Yao, T. An integrative method for identifying potentially dangerous glacial lakes in the Himalayas. Sci. Total Environ. 2022, 806, 150442. [Google Scholar] [CrossRef]
  77. Miles, E.S.; Watson, C.S.; Brun, F.; Berthier, E.; Esteves, M.; Quincey, D.J.; Miles, K.E.; Hubbard, B.; Wagnon, P. Glacial and geomorphic effects of a supraglacial lake drainage and outburst event, Everest region, Nepal Himalaya. Cryosphere 2018, 12, 3891–3905. [Google Scholar] [CrossRef]
  78. Emmer, A.; Vilímek, V.; Zapata, M. Hazard mitigation of glacial lake outburst floods in the Cordillera Blanca (Peru): The effectiveness of remedial works. J. Flood Risk Manag. 2018, 11, S489–S501. [Google Scholar] [CrossRef]
  79. Bhambri, R.; Hewitt, K.; Kawishwar, P.; Kumar, A.; Verma, A.; Tiwari, S.; Misra, A. Ice-dams, outburst floods, and movement heterogeneity of glaciers, Karakoram. Glob. Planet. Chang. 2019, 180, 100–116. [Google Scholar] [CrossRef]
  80. 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]
Figure 1. Study area and distribution of glacial lake area change.
Figure 1. Study area and distribution of glacial lake area change.
Remotesensing 14 04688 g001
Figure 2. Flow chart of the methodology.
Figure 2. Flow chart of the methodology.
Remotesensing 14 04688 g002
Figure 3. Remote sensing extraction map of the glacial lake area (16 November 2018).
Figure 3. Remote sensing extraction map of the glacial lake area (16 November 2018).
Remotesensing 14 04688 g003
Figure 4. Changes in the number and area of glacial lakes from 1990 to 2020. (a) Changes in the number and area of glacial lakes from 1990 to 2020; (b) Changes in the number of area increases in glacial lakes from 1990 to 2020.
Figure 4. Changes in the number and area of glacial lakes from 1990 to 2020. (a) Changes in the number and area of glacial lakes from 1990 to 2020; (b) Changes in the number of area increases in glacial lakes from 1990 to 2020.
Remotesensing 14 04688 g004
Figure 5. (a,b) Thirty-year comparison of glacial lake area changes in Nyalam County.
Figure 5. (a,b) Thirty-year comparison of glacial lake area changes in Nyalam County.
Remotesensing 14 04688 g005
Figure 6. Thirty-year area change of glacial lake in Nyalam County ((ac) GongCuo Lake).
Figure 6. Thirty-year area change of glacial lake in Nyalam County ((ac) GongCuo Lake).
Remotesensing 14 04688 g006
Figure 7. Thirty-year area change of glacial lake in Nyalam County ((a,b) GangxiCo and YinreCo Lake).
Figure 7. Thirty-year area change of glacial lake in Nyalam County ((a,b) GangxiCo and YinreCo Lake).
Remotesensing 14 04688 g007
Figure 8. (a). Interpretation map of GangxiCo glacial lake and surroundings in 2018. (b). Interpretation map of GangxiCo glacial lake and surroundings in 2019.
Figure 8. (a). Interpretation map of GangxiCo glacial lake and surroundings in 2018. (b). Interpretation map of GangxiCo glacial lake and surroundings in 2019.
Remotesensing 14 04688 g008aRemotesensing 14 04688 g008b
Figure 9. The thirty-year trend of climate factors. (A), year, (B), sunshine hours, (C), temperature, (D), wind speed, (E), average relative humidity, (F), precipitation.
Figure 9. The thirty-year trend of climate factors. (A), year, (B), sunshine hours, (C), temperature, (D), wind speed, (E), average relative humidity, (F), precipitation.
Remotesensing 14 04688 g009
Figure 10. Correlation of climate factors.
Figure 10. Correlation of climate factors.
Remotesensing 14 04688 g010
Figure 11. Particle size profile of glacial lake soil sample (a), SNH01-1, (b), SNH01-2, (c), SNH02-1, (d), SNH02-2, (e), SNH01, (f), SNH02.
Figure 11. Particle size profile of glacial lake soil sample (a), SNH01-1, (b), SNH01-2, (c), SNH02-1, (d), SNH02-2, (e), SNH01, (f), SNH02.
Remotesensing 14 04688 g011
Figure 12. Changes before and after the glacial lake outburst.
Figure 12. Changes before and after the glacial lake outburst.
Remotesensing 14 04688 g012
Table 1. Data of the study area.
Table 1. Data of the study area.
Data ClassificationNameResolutionData Source
Remote sensing imagery dataLandsat4/5 TM1990.9.15USGS
(https://glovis.usgs.gov/ (accessed on 1 January 2021))
Landsat4/5 TM1995.10.13
Landsat7 ETM+2000.11.4
Landsat5 TM2005.11.12
Landsat8 OLI2010.11.16
Landsat8 OLI2015.11.24
Landsat8 OLI2020.10.9
GF-12018.10.10Geospatial Data Cloud
(http://www.gscloud.cn/ (accessed on 1 January 2021))
GF-12019.10.5
SRTMDEM dataUSGS (https://glovis.usgs.gov/ (accessed on 1 January 2021))
Meteorological dataMeteorological stationForms for reportMeteorological Data Center of China Meteorological Administration (http://cdc.cma.gov.cn/home.do (accessed on 1 January 2021))
ZoneAdministrative boundary vectorVector dataResource and Environment
Science and Data Center
(https://www.resdc.cn/ (accessed on 1 January 2021))
Table 2. Changes in the number and area of glacial lakes from 1990 to 2020.
Table 2. Changes in the number and area of glacial lakes from 1990 to 2020.
Time (Year)Area (km2)Number (pcs)Change in Area from the Previous Period (km2)Change in Number from the Previous Period (pcs)
199027.94265853
199529.210104511.267446−2
200032.739735553.5296314
200533.861283521.121548−3
201037.055987493.194704−3
201541.941051564.8850647
202052.8461726310.905127
Table 3. Change of glacial lakes for different elevations.
Table 3. Change of glacial lakes for different elevations.
Altitude (m)Current (2020) Total Area (km2)Moraine-Dammed LakeU-type Valley LakeCirque LakeGlacial Erosion LakeLateral Moraine LakeTotal Area of Increase (km2)Number of Lakes Showing Areal IncreasesNumber of Lakes Showing No Increase in Area
<50001.446211 0.31437
5000–51007.9131 113.01751
5100–52003.45532 0.10675
5200–53007.5112 1 2.384103
5300–540014.3114 1 12.263142
>54001.1098 0.00626
Table 4. Evaluation of the accuracy of the glacial lake in different states.
Table 4. Evaluation of the accuracy of the glacial lake in different states.
StatusMisclassification Error (%)Omission Error (%)Kappa Coefficient
Half icing0.671.400.86
Total icing0.551.080.90
Unfrozen0.230.860.92
Table 5. Thirty-year annual change in Nyalam County.
Table 5. Thirty-year annual change in Nyalam County.
1990–20002000–20102010–20201990–2020
Glacial lake growth area (km2)0.6991.0510.131.88
Annual temperature (°C)3.694.3124.1584.042
Summer temperature (°C)10.291910.606710.419110.441
Winter temperature (°C)−2.897−1.673−1.987−2.2
Annual precipitation (mm)551.76613.08632.93595.51
Summer precipitation (mm)203.63198.86220.95206.36
Precipitation in winter (mm)99.72102.79141.1111.6
Sunshine hours (h)2515.492477.952507.762499.58
Average wind speed (m/s)4.6864.0634.0694.29
Average relative humidity (%)67.3967.7162.966.34
Table 6. (a). Results of cross-detection. (b). Ecological detection and risk factor detection.
Table 6. (a). Results of cross-detection. (b). Ecological detection and risk factor detection.
(a) Interaction_detector
Sunshine hoursAverage temperatureWind speedRelative humidityPrecipitation
Sunshine hours0.082
Average temperature0.2700.073
Wind speed0.4920.2480.193
Relative humidity0.5590.5570.6060.451
Precipitation0.5310.5580.7320.6770.342
(b) Sig. F test: 0.05
TemperaturePrecipitationHumidity
Temperature
PrecipitationY
HumidityYY
q statistic0.9977031170.9994941560.994113327
Table 7. Table of coefficients.
Table 7. Table of coefficients.
Coefficients a
ModelUnstandardized CoefficientStandardized CoefficienttSig.Covariance Statistics
BStandard ErrorTrial VersionToleranceVIF
1(Constant)163.12646.453 3.5120.002
Sunshine hours−0.0100.012−0.134−0.8250.4190.9261.080
Average temperature−2.4622.180−0.239−1.1300.2710.5461.833
Average wind speed−6.9223.013−0.513−2.2980.0320.4902.040
Average relative humidity−0.9610.256−0.660−3.7580.0010.7911.264
Precipitation−0.0060.009−0.122−0.6740.5080.7461.341
a Dependent variables: area.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Qu, G.; Dai, X.; Cheng, J.; Li, W.; Wang, M.; Liu, W.; Yang, Z.; Shan, Y.; Ren, J.; Lu, H.; et al. Characterization of Long-Time Series Variation of Glacial Lakes in Southwestern Tibet: A Case Study in the Nyalam County. Remote Sens. 2022, 14, 4688. https://doi.org/10.3390/rs14194688

AMA Style

Qu G, Dai X, Cheng J, Li W, Wang M, Liu W, Yang Z, Shan Y, Ren J, Lu H, et al. Characterization of Long-Time Series Variation of Glacial Lakes in Southwestern Tibet: A Case Study in the Nyalam County. Remote Sensing. 2022; 14(19):4688. https://doi.org/10.3390/rs14194688

Chicago/Turabian Style

Qu, Ge, Xiaoai Dai, Junying Cheng, Weile Li, Meilian Wang, Wenxin Liu, Zhichong Yang, Yunfeng Shan, Jiashun Ren, Heng Lu, and et al. 2022. "Characterization of Long-Time Series Variation of Glacial Lakes in Southwestern Tibet: A Case Study in the Nyalam County" Remote Sensing 14, no. 19: 4688. https://doi.org/10.3390/rs14194688

APA Style

Qu, G., Dai, X., Cheng, J., Li, W., Wang, M., Liu, W., Yang, Z., Shan, Y., Ren, J., Lu, H., Wang, Y., Zeng, B., & Atasoy, M. (2022). Characterization of Long-Time Series Variation of Glacial Lakes in Southwestern Tibet: A Case Study in the Nyalam County. Remote Sensing, 14(19), 4688. https://doi.org/10.3390/rs14194688

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