Spatially and Temporally Resolved Monitoring of Glacial Lake Changes in Alps During the Recent Two Decades

Climate warming is intensifying the melting of glaciers and the growth of glacial lakes in the Alps, which has a profound impact on the management of water resources and high-mountain hydropower in this region. However, the research on the spatial distribution and temporal evolution of the Alps glacial lakes of various types still lacks a holistic view. In this study, we developed an inventory of Alps glacial lakes of different types and then obtained the annual areas of these lakes from 2000 to 2019 using JRC Global Surface Water and Global Land Analysis and Discovery data at a resolution of 30 m. A total of 498 glacial lakes (>0.01 km2) with the net area of 33.77 ± 6.94 km2 were identified in the Alps in 2019 and are mainly distributed in the western and central Alps. These Alps glacial lakes, with the area ranging 0.01–1.59 km2, are generally dominated by small-sized ones. The comparison of lakes of different types indicated that ice-uncontacted lakes are dominant in number and area, accounting for 59.4 and 58.4%, respectively. In terms of the elevation distribution, almost half of the lakes are concentrated at the altitude of 2,250–2,750 m (a.s.l.). Meanwhile, the mean altitude of small glacial lakes is higher than that of large lakes. The distribution of ice-contacted lakes and supraglacial lakes were more concentrated, and the mean altitude was higher. During the study period, the number, area, and water volume of glacial lakes were increasing, but the expansion varied between different periods. The changing trends of the glacial lake area and volume were consistent and presents in three stages, as the glacial lake expanded rapidly in the first 5 years and in the last 7 years and remained relatively stable between 2005 and 2012. The number and area of glacier-fed lakes increased rapidly, while the non-glacier-fed lakes were relatively stable. The area change rate of supraglacial lakes was the largest (+47%). This study provides a spatially-complete and temporally-consecutive picture of glacial lake changes in the Alps and can be greatly helpful for future research on climate-glacier-lake interactions, glacial lake outburst floods, and freshwater resources in this region.


INTRODUCTION
With global warming, the temperature rising is particularly obvious in high-altitude regions, and glaciers show accelerated retreat or thinning (Zemp et al., 2006;Fischer et al., 2015;Zemp et al., 2019). Persistently glacier retreat and meltwater supply promotes the circulation of surface water and increases the abundance of glacial lake (Yao et al., 2018). As the source of precious water resources and mountain flood in high mountain regions, glacial lakes are not only the carrier of environmental information, but also the indicator of climate change (Che et al., 2014). At the same time, glacial lakes promote the melting of alpine glaciers and landscape evolution and acts vital role in surface water cycle and glacier disaster evolution effect Song et al., 2017). On the one hand, the glacial lake retains considerable glacier meltwater, a feature which plays a temporary role in storing freshwater resources and maintaining sea level (Wang et al., 2020). On the other hand, in view of the global expansion of glacial lakes and the instability of many glacial lakes (Shugar et al., 2020), their unexpected outburst will lead to catastrophic floods which seriously threaten downstream residents, infrastructure, and regional ecological environment security. Therefore, the consecutive monitoring of glacial lake changes is particularly important and urgently needed (Allen et al., 2019;Dubey and Goyal, 2020;Veh et al., 2020).
Given the rapid development of remote sensing technology, varieties of high-resolution satellites and sensors have shown outstanding performance for the identification, extraction, and monitoring of glacial lakes (Petrov et al., 2017;Begam and Sen, 2019). For example, the near five-decade Landsat archive launched by NASA has become the main data source for monitoring glacial lake dynamics. On the basis of remote sensing data, many scholars or institutions have combined research needs to perform glacial lake cataloging, especially for High Mountain Asia (HMA) (Chen et al., 2020;Wang et al., 2020). For example, remote sensing imagery has been employed to map glacial lake inventories for regions of Bhutan (Veettil et al., 2016), Nepal (Khadka et al., 2018), China , and Patagonia (Wilson et al., 2018) and provides an effective way to reveal the characteristics of glacial lake temporal and spatial changes and understand the response of glacial lakes to climate change, as well as offer a data source for subsequent glacial lake water volume estimation and outburst floods research.
So far, the researches on Alps glacial lakes are sporadic and regional (Buckel et al., 2018), few studies explore the spatial distribution and dense time-series changes of glacial lakes for the entire Alps area, and this shortage leads to the lack of detailed understanding of the changes of glacial lakes in this region. Recent glacier studies indicate a rapid retreat of Alps glaciers (-39 km 2 a −1 ) and ice thickness variations (-0.5 m a −1 to -0.9 m a −1 ) (Huss, 2012;Zekollari et al., 2018;Sommer et al., 2020). Rapidly melt of Alps glaciers may cause substantial changes in glacial lakes during the past several decades. Besides, lakes connected to glaciers provide positively feedback on the glacier system, resulting in changes in mass balance, and enhanced melting (Song et al., 2017). In order to deepen the understanding of the impact of glacial lakes on mountain environment and the response to glacier change, detailed data on glacial lake evolution and their spatially and temporally heterogeneous changing characteristics are needed. Here, we aim to develop an inventory of Alps glacial lakes with the maximum inundation distribution at 30 m resolution and further characterize the annual variations of glacial lakes from 2000 to 2019. The specific objectives of this research are as follows: 1) to map the glacial lake inventory of annual time series in the Alps region over the past 20 years to fill the gap in a complete glacial lake inventory of the entire region; 2) to reveal the heterogeneous changes of glacial lakes of different types, scales, and elevation distributions in the Alps; 3) to provide open access to the resulting inventory data set so as to offer geographic information on glacial lakes and improve awareness of the Alps glacial lakes.

Study Area
The Alps is the highest mountain range in Europe and is located in southern Europe and known as the "European water tower" (Figure 1). The Alps accounts for only 11% of Europe's total area, but provides more than 90% of the water supply of Europe, especially in arid regions and in summer (Vanham, 2012). The mountain ranges arc in shape and extend from west to east, starting from the Mediterranean coast near Nice in southeast France in the West and ending in the Vienna basin in Austria in the East (Paul et al., 2020). The mountain ranges are 1,200 km long and 130-260 km wide, with a mean altitude of approximately 3,000 m and a total area of up to 220000 km 2 (Auer et al., 2007;Zemp et al., 2015). The highest peak is Mont Blanc (4,808 m).
Located between temperate and subtropical latitudes, the Alps is the boundary between the temperate continental humid climate of central Europe and the summer dry climate of the subtropical southern Europe. The peak is cold all year round, and the annual mean temperature is 0°C at the elevation of 2000 m, and the annual precipitation is up to 1,200-2000 mm (Frei et al., 2003). The annual precipitation in the edge area varies significantly from that in the interior of the mountains. The elevation of 3,000 m is the maximum precipitation zone, and the annual precipitation in the high mountains exceeds 2,500 mm (Vanham, 2012;Ács et al., 2018). The area is home to more than 1,000 modern glaciers covering a total area of 3,600 km 2 . These glaciers were mostly shrinking and thinning rapidly, and substantial glacial meltwater provide continuous supply for the development of glacial lakes (Paul et al., 2020).

Data Source
We developed our glacial lake inventory of the Alps according to the JRC Global Surface Water (GSW) and the Global Land Analysis & Discovery (GLAD) data stored in the Google Earth Engine (GEE, https://earthengine.google.com/) cloud computing platform for free to global users. The GSW is a spatial and temporal distribution map of the global surface water generated by Pekel et al. (2016) using 4,185,439 scenes of Landsat images acquired between 1984 and 2019 and provides spatial raster data on the water extent and variations. Similar to the GSW data set, the GLAD data (https://www.glad.umd.edu/ dataset/global-surface-water-dynamics) is generated by Pickens et al. (2020) using a classification tree to categorize the land and water in all Landsat scenes from 1999 to 2020 and perform a timeseries analysis to produce maps that characterize inter-annual and intra-annual open surface water dynamics. The spatial resolution of both datasets is 30 m. In this study, the annual data layers of the GSW (JRC Yearly Water Classification History, v1.2) as the primary data source were used to map the annual areas of glacial lakes in the Alps. To ensure the accuracy of the glacial lake inventory, we use GLAD annual data as supplementary information for manually inspecting and editing the mapping results from the GSW.
In addition, the Shuttle Radar Topography Mission (SRTM V3) digital elevation data provided by NASA JPL (http://gdex.cr. usgs.gov/gdex/) with a spatial resolution of up to 1 arc-second are used to extract the elevation of the glacial lakes. RGI 6.0 glacier data (https://www.glims.org/RGI/) are utilized to determine the buffer range of glacial lake data distribution.

Glacial Lake Area Extraction
At present, the extraction methods of glacial lakes are mainly divided into two categories: automatic mapping and manual visual interpretation. Automatic extraction algorithms have the advantages of being relatively efficient and less labor-consuming, but manual visual interpretation has higher accuracy and ease of implementation (Huang et al., 2018); moreover, the latter is widely used in different regions for delineating the extent of glacial lakes (Wang et al., 2013a;Zhang et al., 2015;Wang et al., 2020). The public release of the GSW and GLAD data provide a great convenience for extracting information on glacial lakes. However, the dataset maps all surface water bodies of different types, e.g., lakes, rivers, floodplains, and paddy water bodies. Therefore, many non-lake water bodies were excluded from glacial lakes in the study area. Manual inspection can also largely reduce influences of mountain shadows and exclude other types of water bodies (such as rivers, and tectonic lakes).
The research on glacial lake has been carried out over a long period, but a consensus on the remote sensing-based specific standards of glacial lake categorization is yet to be reached. Researchers typically combine their research purposes and the geographical conditions of the study area to establish some thresholds to differentiate the glacial lakes from other natural or artificial lakes (Li et al., 2020). Prior studies emphasize the relationship between glacial lakes and glaciers and suggest that lakes within 3, 5, or 10 km from the margin of glaciers are regarded as glacial lakes (Wang et al., 2013b;Worni et al., 2013;Veh et al., 2018). Others define glacial lakes by setting altitude thresholds and other parameters (Raj and Kumar, 2016;Wang et al., 2017). The extraction of glacial lakes based on altitude threshold is not universal, and the approach of establishing distance buffer is still a common method to extract and assist glacial lake recognition. According to the visual judgement of numerous glacial lake cases on the geographical location of Alps glaciers and surrounding lakes, we determine that it is reasonable to map glacial lakes within the distance of 5 km from the glacier termini. We also conducted a procedure of quality control on the mistaken or missed glacial lakes around the buffering distance by manual inspection. In addition, the glacial lake inventory aims to provide basic data for the management of freshwater resources in alpine regions and the research of glacial lake outburst flood (GLOFs), glacial lakes with an area of more than 0.01 km 2 are considered to be a kind of lakes that pose a threat to downstream facilities when dam break occurs (Worni et al., 2013;Veh et al., 2020)). Hence, we established a 5-km buffer of distance to upstream glaciers and set a threshold of 0.01 km 2 as the minimum area of Alps glacial lakes to exclude other water bodies in the region. In other words, we defined glacial lake as natural lake formed by glaciation, with an area of greater than 0.01 km 2 within 5 km from the terminal of the glacier and then employ the high-resolution Google Earth imagery to assist in the manual identification for improved accuracy of the glacial lake inventory.
To estimate the long-term change rates of area time series, a simple and robust strategy is to calculate the annual change rate by referring to the observed values in starting year and ending year. Specifically, the area change rate of the glacial lake in each period is obtained via the difference between the starting and ending lake area divided by the number of years between them, and the corresponding standard error is obtained by combining the annual glacial lake area with the multi-year average area.

Glacial Lake Classification
Prior studies categorize glacial lakes according to their physical characteristics, including the glacial lake size, the geometrical relationship with modern glaciers, the geographical location of the lake, and whether or not the lake is being fed by glacial meltwater (Lizong et al., 2011;Yao et al., 2018). Different researchers classify glacial lakes according to their research aims or focuses. The common approaches can be summarized as two systems. One method is implemented on the basis of the relationship between the glacial lake and glacier. For example, to examine the inter-annual variation characteristics of glacial lakes in HMA, Chen et al. (2020) classified glacial lakes into four categories: 1) pro-glacial lake, 2) supraglacial lake, 3) unconnected glacial lake, and 4) ice-marginal lake. The other method is in accordance with the formation mechanism of glacial lakes (Wilson et al., 2018), which enables them to be divided into 1) glacial erosion lakes, 2) moraine-dammed lakes, and 3) iceblocked lake, among others. Evaluating the formation mechanism of glacial lakes on the basis of remote sensing images only is generally difficult. The classification system of glacial lakes used in this study is developed according to the relationship between glacial lakes and glaciers, as this mechanism reflects the important role that glaciers play in the formation of glacial lakes. Wang et al.
(2020) provides a practical reference for our classification.
In this research, we divided glacial lakes into four categories: 1) glacier-fed lakes (GFL, lakes fed by modern glacier meltwater) which were further divided into 1) ice-uncontacted lakes (IUL, lakes not in contact with glaciers but fed by glacier meltwater), 2) ice-contacted lakes (ICL, lakes that are in close contact with a glacier terminal), and 3) supraglacial lakes (SGL, lakes developed on glacier surfaces); and 2) non-glacier-fed lakes (NGFL, lakes are formed due to the historical glacial movements and are not fed by modern glacial meltwater). To obtain reliable classification results, the glacial lakes are manually distinguished by using high-resolution imagery from Google and Bing and the glacier data of RGI 6.0 (Pfeffer et al., 2014).

Uncertainty Estimation
Given the impact of remote sensing image quality (such as spatial resolution, mountain shadow, cloud cover, and water extraction algorithm, etc.,), the interpreters' prior knowledge, and the definition of glacial lake, errors inevitably occur when extracting glacial lake information (Zhang et al., 2015;Li et al., 2020). In this study, we used a semiautomatic method to map the glacial lake inventory, and every lake was manually validated and edited. However, no appropriate and quantifiable uncertainty estimation exists for the error caused by artificial correction process and the definition of glacial lake. To improve the credibility of the inventory, recheck it is necessary. Previous studies have shown that the mixed pixel is the key factor of error sources (Mixed pixel refers to the existence of different types of ground objects in a pixel, mainly at the boundary of ground types. The existence of mixed pixels is one of the main factors affecting the recognition and classification accuracy) (Hanshaw and Bookhagen, 2014). In other words, the error is mainly caused by the spatial resolution of the image (Wang et al., 2020). Using the buffer zone of a half-pixel distance outside the boundary is optimal to evaluate the uncertainty (Fujita et al., 2009;Salerno et al., 2012). Therefore, the uncertainty of a single lake area within 1 standard deviation (1σ) can be expressed as (Hanshaw and Bookhagen, 2014): where Error is the area error of each lake, G is the image spatial resolution (the resolution of the dataset is 30 m), P (m) is the perimeter of each lake, 0.6872 is the revised coefficient under 1σ (i.e., approximately 69% of pixels conform to the error distribution), R is the uncertainty of the lake area, and A is the lake area. Table 1 shows that the uncertainty of the glacial lake inventory is inversely proportional to the area size: as the size increases, the uncertainty gradually decreases. In addition, the uncertainty of the glacial lake area extraction is approximately 50% when the lake area is 0.01 km 2 . Therefore, the minimum area of 0.01 km 2 was set as the threshold in this work, and this approach improved the accuracy of the inventory to a certain extent. Overall, the mean uncertainty of the inventory in the Alps is approximately 20.5% (The uncertainty in the study region is dominated by small lakes, yet the result is acceptable).

Characteristics of the Spatial Distribution of Glacial Lakes in the Alps
A total of 498 glacial lakes (>0.01 km 2 ) were identified in 2019 by using remote sensing data products and cover a net area of 33.77 ± 6.94 km 2 . Glacial lakes are widely distributed across the Alps, and several of the largest glacial lakes (>1 km 2 ) in the study area are mainly concentrated in the western and central regions (Figure 1). On the whole, the Alps are dominated by small glacial lakes, and all glacial lakes are in the size range of 0.01-1.59 km 2 . Glacial lakes less than 0.1 km 2 account for 86.3% of all the glacial lakes of the inventory in number, but their area is only 41.8% of the total. Only 68 glacial lakes exceed 0.1 km 2 , yet they account for 58.2% of the total area. Therefore, although many small glacial lakes exist in the Alps, the total area of glacial lakes in the region is still dominated by relatively large lakes. Figure 2 presents the spatial distribution of glacial lakes of different types. GFL are dominant in number and area, and NGFL are in the minority. Specifically, the number of IUL take the leading position, accounting for approximately 59.4%. NGFL, ICL and SGL account for 20.4, 11.4 and 8.8%, respectively. In terms of area, the four types (IUL, NGFL, ICL, and SGL) of glacial lakes account for 58.4, 12.6, 20.4, 8.6, and 12.6%, respectively. Compared with other lakes, ICL have larger sizes, probably because they are connected with the glacier termini and can evolve with the retreating of glaciers. Compared with ICL, SGL do not exist all year round and their area fluctuation is evident.
The elevation bins by the step of 250 m were defined to analyze the dependence of glacial lake distribution on elevation. Figure 3 shows the distribution of the Alps' glacial lakes in different elevation zones. As the altitude increases, the quantity and area of glacial lakes are approximately normally distributed. The altitudes of all glacial lakes mapped in 2019 ranged from 1,060 m a.s.l. (above sea level) to 3,133 m a.s.l., of which 72% are distributed at the altitude belt of 2,250-2,750 m a.s.l., with a mean of 2,410 m. The lakes between 2000 and 2,500 m also have the largest proportion in area, accounting for 70.3%. Further analysis reveals that in different elevation bins, the distribution of glacial lakes at different sizes is also uneven and vary greatly ( Figure 3C). Glacial lakes predominantly range from 2000 m to 2,750 m, and most are small lakes (≤0.04 km 2 ), with the larger lakes distributed in relatively low altitude regions  mainly because the land is flat in the area with low altitude, thereby providing a good water storage environment for some large NGFL and a few IUL.
In addition, the distribution of glacial lakes varies with the types of glacial lakes given the influence of glacial movement and action (Figure 4). Among the lakes, ICL and SGL are closely related to glaciers and therefore have a higher and more concentrated elevation distribution, with a mean elevation of approximately 2,600 m a.s.l. However, IUL and NGFL have more discrete elevation distributions, with lower mean elevations (2,400 m and 2,100 m, respectively). NGFL are particularly prominent as products of pre-glacial movement.

Number and Area Changes of Glacial Lakes in the Alps
We further obtained the annual glacial lake areas (>0.01 km 2 ) in the Alps from 2000 to 2019. During that period, the number and area of glacial lakes present an ascending tendency ( Figure 5), yet the increasing rates varied between different periods. There were 485 and 498 glacial lakes in 2000 and 2019, respectively. The total area of the glacial lakes in the study region increased from 29.78 ± 6.48 km 2 in 2000 to 33.77 ± 6.94 km 2 in 2019 and expanded by approximately 3.99 ± 0.46 km 2 (by 13%). During the sub-period of 2000-2004, the quantity and total area of glacial lakes presented relatively more rapid increases, with an increase of 1.74 ± 0.21 km 2 (5.8%) and an annual growth rate of 0.44 ± 0.07 km 2 a −1 . During 2005-2012, the number and area of glacial lakes were relatively stable. From 2013 to 2019, the number of glacial lakes changed gradually and the area increased rapidly, increasing by 1.68 ± 0.20 km 2 (5.2%) with an annual growth rate of 0.28 ± 0.05 km 2 a −1 .
Overall, the majority of glacial regions experienced glacial lake growth during the study period (Shugar et al., 2020), particularly in the central Alps. However, some separate districts recorded a decrease in glacial lake area (light-blue grid cell in Figure 5C) and were mainly affected by some large NGFL.
Globally, the number of available lake inventories has gradually increased recently, the glacial lake inventory from individual basin (Kumar et al., 2020) to long-time series and global scale (Shugar et al., 2020) deepens the understanding of the evolution differences of glacial lakes in different regions, so as to better constrain and simulate the hydrological response of glaciers to climate change (Wood et al., 2021). we discover that, compared with those of other regions with glacial lake distributions such as the Andes (the total area of glacial lake increased by 14% from 2000 to 2016, excluding Lake Argentino, Lake Viedma, Lake O'Higgins/San Martín) (Wilson et al., 2018) and Nepal (increased by approximately 17% from 1997 to 2017) (Haritashya et al., 2018), the change rate of Alps' glacial lakes in the recent years seems to be relatively slow. However, compared with Alaska (Shugar et al., 2020) and Himalayas (increased by 5% from 2000 to 2015) , the change rate seems to be relatively fast (Due to the difference of glacial lake identification and statistical years, the actual results may be a little biased.).

Changes in the Water Volume of Glacial Lakes
The "real" water volume of a glacial lake can be calculated by combining the water level and shoreline data and accurate lake bathymetry. However, most glacial lakes are located in harsh Frontiers in Earth Science | www.frontiersin.org August 2021 | Volume 9 | Article 723386 6 environments and remote areas, making the field measurements of glacial lake depths difficult and infeasible. Therefore, the empirical model fitted by the relationship of glacial lake area-depth-volume is a commonly-used method for estimating the glacial lake water volume on regional scale. The water volume estimation model (Eq. 3) of glacial lakes proposed by Cook and Quincey, (2015) has strong applicability to the Alps.
Where V and A represent the water volume and area of each lake, respectively. Substantial retreat and downwasting of glaciers have occurred in the Alps since the Little Ice Age due to global warming and is expected to exacerbate. Accelerated glacier shrinkage has been reported over the past few decades. In the first 5 years of the 21st Century, the mass-change rate was close to −1 m w.e.a −1 (Huss, 2012); from 2000 to 2012, it was -0.71 m w.e.a −1 (Sommer et al., 2020). The glacier mass  Frontiers in Earth Science | www.frontiersin.org August 2021 | Volume 9 | Article 723386 7 change rate decreased to -0.5 m w.e.a −1 from 2006 to 2016 (Zemp et al., 2019). Then, the glacier mass entered a period of rapid retreat (Zekollari et al., 2018;Carrivick et al., 2020), this outcome is consistent with research findings which show several years of record-breaking volume loss and glacier retreat after 2012 (2015 and 2018) (Zekollari et al., 2020). Referring to the glacier mass changes reported recently suggest that the changing tendency of glacial lakes in the Alps is generally consistent with the temporal patterns of negative glacier mass balance. The volume change of a glacial lake was similar to its area change and generally behaves in a three-stage evolution trajectory ( Figure 6). During the study period, lake volumes increased from 0.57 to 0.66 km 3 mainly in the two sub-periods of 2000-2004 and 2013-2019, with growth rates of 5 and 8%, respectively. From 2005 to 2012, the total water volume of glacial lakes was relatively stable.

Heterogeneous Changes of Glacial Lakes by Different Types
Differences in the area changes of the four types of glacial lakes were further examined as illustrated in Figure 7. During the study period, the number and area of GFL increased rapidly because of the influence of glacier retreating and glacial meltwater supply, while that of NGFL were relatively stable (Figure 7). The number of GFL increased from 384 in 2000 to 397 in 2019, and their area expanded from 25.85 ± 5.39 km 2 in 2000 to 29.52 ± 5.81 km 2 in 2019, with an increase of 3.67 km 2 (+14.2%). By contrast, the differences in the numbers (+0) and area (+0.32 km 2 , +8%) of NGFL during 2000-2019 varied slightly given the absence of the influence of modern glaciation.
For the past 2 decades, SGL were closely related to glacier changes and were vulnerable to glacier ablation and precipitation given their small sizes. As a result, they have the highest rate of change in area (+47%), as well as large inter-annual fluctuations in area. Except for the case with SGL, the glacier terminus retreats faster and the glacial meltwater increases due to the rising temperature, thereby providing a favorable condition for the expansions of IUL and ICL (their area growth percentages were 11 and 12%, respectively). A feedback relationship also occurs between glacial lakes and glaciers. First, as the product of glacier movement, glacial lakes receive energy and materials transported by glaciers and can highly respond to the changes of glaciers. Second, the expansion of glacial lakes connected with glaciers accelerated the melting of glaciers while also indirectly promoting the appearance of new and higher elevation glacial lakes (mainly ICL) at the end of the glacier and the conversion between glacial lakes (transformation from SGL to ICL or ICL to IUL).

Comparison of Inventoried Glacial Lakes With Previous Study
We compared our glacial lake inventory with that of Buckel et al. (2018) for the same period over the approximately same  Frontiers in Earth Science | www.frontiersin.org August 2021 | Volume 9 | Article 723386 extent. Our inventory includes 84 glacial lakes (≥0.02 km 2 ) with a total area of 6.33 ± 1.28 km 2 , while the artificially interpreted inventory of Buckel et al. (2018) contains 105 glacial lakes (≥0.02 km 2 ) with a total area of 6.26 km 2 . The differences of glacial lakes between these two inventories are 21 and 0.07 km 2 in the total number and area, respectively. Some of reasons for the missed glacial lakes in our inventory are the interference of some not-so-great observations (shade or snow), identification criteria for glacial lakes, and a few glacial lakes that have dried up. The major reason for the shortage of glacial lake inventory is the difference of extraction methods. Different from our glacial lake inventory, Buckel mapped Austria glacial lakes by manual digitization using multi-source remote sensing image, such as Google Earth, Esri World Imagery, and Orthofoto Österreich.
In short, the semi-automatic method shows unprecedented advantages in glacial lake extraction efficiency. However, in terms of extraction accuracy, it is slightly insufficient compared with the manual digitization method. It cannot be ignored that the integrity and availability of data sources fundamentally affect the extraction of glacial lakes. Lake mapping using satellite images and orthophotos from multiple sources from different years is also a vital channel to map a complete glacial lake inventory.

CONCLUSION
Glacial lakes in the alpine regions are highly sensitive to glacier mass evolution and climate change. However, the research on the spatial distribution and temporal evolution of various types of glacial lakes in the Alps still lacks a holistic view. Through a semi-automatic approach, we used remote sensing products (the GSW and GLAD) to map annual glacial lakes and compiled an inventory of these lakes in the Alps. Our results show that 498 glacial lakes (>0.01 km 2 ) were identified in 2019 using remote sensing data within the study area, covering a total area of 33.77 ± 6.94 km 2 . The glacial lakes are widely distributed and concentrated mainly in the western and central regions of the Alps. On the whole, the Alps are dominated by small glacial lakes in number, with all glacial lakes in the size range of 0.01-1.59 km 2 . The differences among different types of glacial lakes are also significant. Glacier-fed lakes are dominant in both number and total area, and non-glacier-fed lakes are in the minority. Among all the lakes, ice-uncontacted lakes occupy the leading position in number and area, accounting for 59.4 and 58.4%, respectively. In terms of the elevation distribution, almost half of the lakes are concentrated between 2,250 and 2,750 m above sea level. Meanwhile, the mean elevation of the small glacial lakes is higher, and the larger lakes are distributed in relatively low altitude regions. The distributions of ice-contacted lakes and supraglacial lakes were more concentrated, and their mean elevations were higher. However, ice-uncontacted lakes and non-glacier-fed lakes have more discrete altitude distributions with lower mean elevations.
Referring to the glacier mass changes reported recently indicate that the changing tendency of glacial lakes in Alps is generally consistent with the temporal patterns of negative glacier mass balance. Glacier retreating and additional meltwater supply could play a leading role in the area expansion of most glacial lakes. During the study period, the quantity and area of glacial lakes show an increasing tendency, while the expansion varied between different periods. The area shows a three-stage trend whereby the glacial lake expansion is fast in the first 5 years and the last 7 years (with an annual growth rate of 0.44 ± 0.07 km 2 a −1 and 0.28 ± 0.05 km 2 a −1 ) and remained relatively stable between 2005 and 2012. The water volume and area of the glacial lake in the region are similar and reveal a consistent trend of change.
The number and area of glacier-fed lakes increased rapidly during the study period because of the influence of glacier retreat and glacial meltwater supply, while those of non-glacier-fed lakes were relatively stable. For the past 2 decades, the area change rate of supraglacial lakes is the largest percentage (+47%). Given the rapid melting of glaciers, several glacial lakes are undergoing irreversible changes.
In summary, this study attempted to produce a complete glacial lake inventory in the Alps and examined the dense time-series of lake area changes linked to glacier retreating through detailed classification of glacial lakes. This work is expected to provide basic data and information for cryospheric hydrology research, water resource utilization and management, and glacier hazard assessment in the Alps.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
RGI. We are also grateful to YW for numerous fruitful discussions, advice, and comments on the present article. Finally, we are grateful to the scientific editor and the reviewers for their careful review of the article used to improve the study.