Regional Geomorphological Conditions Related to Recent Changes of Glacial Lakes in the Issyk-Kul Basin , Northern Tien Shan

To assess the current state of glacial lakes, we examine the seasonal lake-area changes of 339 glacial lakes in the Teskey and Kungoy Ranges of the Issyk-Kul Basin, Kyrgyzstan, during 2013–2016 based on optical satellite images (Landsat7/ETM+ and 8/OLI). The glacial lakes are classified into six types based on their seasonal variations in area: stable, increasing, decreasing, appearing, vanishing, and short-lived. We then track the number of each type in a given year and examine how each number changes from one year to the next. We find that many appearing, vanishing, and short-lived types occurred in both mountain ranges, having a large variability in number that is not directly related to the local short-term summer temperature anomaly, nor to precipitation or glacier recession. However, those in the Teskey Range vary significantly more than those in the Kungoy Range. To determine if the changing number and distribution of the various lake types may be due to changes in ground ice, we apply differential interferometric synthetic aperture radar (DInSAR) analysis using ALOS-2/PALSAR-2 for the debris landforms behind which glacial lakes appear. In the Teskey Range, ground ice occurs in 413 out of a total of 930 debris landforms, whereas in the Kungoy Range, ground ice occurs in 71 out of 180. In zones with predominant glacier-retreat during 1971–2010 (from Corona KH-4B and ALOS/PRISM), the Teskey Range had 180 new lake depressions as potential lake-basins, whereas the Kungoy Range had just 22. Existing depressions also expanded when melting ice produced subsidence. Such subsidence, together with debris landforms containing ground ice and ice tunnels, appear to cause the observed large number variability. In particular, the deposition of ice and debris by tunnel collapse or the freezing of storage water in a debris landform may close-off an ice tunnel, causing a lake to appear. Subsequent re-opening via melting of such blockage would produce either a vanishing or a short-lived type. In this way, the large variability in the number of each lake type and the distribution of types over this four-year period arises from regional geomorphological conditions and not directly from the local short-term summer temperature anomaly and precipitation or glacier recession.

The Issyk-Kul Basin in northwestern Kyrgyzstan is a region with frequent disasters related to glacial lakes.In the Teskey Range, floods from the Angisay glacial lake occurred on 14 June 1974, 17 June 1975, 25 June 1980, and 14 July 1980 [6,7].The two floods in 1980 had particularly large water volumes, but resulted in no downstream casualties.Other floods from glacial lakes in this range occurred on 7 August 1952 in Kichi-Kyzylsuu and 12 August 1964 in Tamga River [6,19].Between 2006 and 2014, this mountain range experienced four large drainages from glacial lakes [4,10].Here, we refer to them simply as the "short-lived type" of glacial lakes that drastically grow and drain within a few months [4,10].For instance, on 24 July 2008, a large drainage occurred at the western Zyndan glacial lake in the northern Teskey Range [4].This lake appeared in just two and a half months and then suddenly drained.The flood caused extensive damage, killing three people and many livestock, as well as destroying a bridge, a road, two houses, and crops.Three other large drainages also damaged a bridge, irrigation channel, roads, monuments (tomb), and agriculture fields [10].Also, in this case the lake water drained through an ice tunnel within a debris landform containing ground ice.In the Kungoy Range, a flood from the Choktal glacial lake occurred on 27 September 1978 [6].
The recent area changes of glacial lakes are influenced by climatic conditions and glacier changes.For instance, Kapitsa and others [22] showed that the evolution of glacial lakes at the glacier front in the Djungarskiy Alatau Range is controlled by climate change, glacier melt, melt of buried ice, as well as local topography and moraine stability.Another study linked the changes of glacial lakes over the Tien Shan Mountains to climate warming and glacier shrinkages [23].In particular, the glacier shrinkage is thought to be important for the areal expansion of glacial lakes in the northern Tien Shan [19,24].
Compared to the "moraine outburst type" drainages of proglacial lakes, which occur by incision of a debris dam [22][23][24][25], few reports exist of the geomorphological conditions and mechanisms of "tunnel type" drainages of short-lived lakes through an underground channel.Is the area change of lakes with tunnel-type drainage directly related to glacier recessions or climate conditions?To help address such issues, we use here satellite-data analysis to examine the current situation of glacial lakes in the Teskey and Kungoy Ranges.Then we use this analysis plus field observations to examine their geomorphological conditions.Finally, we compare the geomorphological conditions of the two ranges.

Study Area
Our study area is the Issyk-Kul Basin, a large intermountain depression between the Teskey and Kungoy Ranges.This basin includes the northern part of the Teskey and the southern part of the Kungoy Ranges (Figure 1).The Teskey Range traverses 354 km in the E-W direction, with the highest point on Karakol (5216 m a.s.l.).The Kungoy Range traverses 285 km in the E-W direction, with the highest point on Chok-Tal (4770 m a.s.l.).The Teskey and Kungoy mountain ranges consist mainly of crystalline basement rocks of Archean to Middle Paleozoic age, covered by volcano-sedimentary and sedimentary strata of Devonian-Carboniferous age [26].Along the ridge line of the eastern part of the Teskey Range, the rocks originate from granites and granitoids of the Paleozoic, whereas the central part has Late Ordovician granitoids and the western part has Paleozoic granitoids.Along the ridge line of the eastern part of the Kungoy Range, the rocks originate from Paleozoic granitoids and the western part is Late Ordovician granitoids [27].Thus, the lithology around the glacial areas are essentially the same in the two mountain ranges.
The climate is continental, but due to the topography, it has much spatial variability [28].part has Late Ordovician granitoids and the western part has Paleozoic granitoids.Along the ridge line of the eastern part of the Kungoy Range, the rocks originate from Paleozoic granitoids and the western part is Late Ordovician granitoids [27].Thus, the lithology around the glacial areas are essentially the same in the two mountain ranges.
The climate is continental, but due to the topography, it has much spatial variability [28].For example, the average annual precipitation values during 1981-1990 were 594 mm at the Chong-Kyzylsuu station (2855 m) (data for Chong-Kyzylsuu station averaged only for 1981-1987), 297 mm at the Tien Shan station (3614 m), 307 mm at the Cholpon-Ata station (1645 m), and 625 mm at the Kyrchyn station (2100 m).See Figure 1 for locations.Over the same period, the annual average temperature was 0.2 °C at Chong-Kyzylsuu, −7.6 °C at Tien Shan, 8.1 °C at Cholpon-Ata, and 2.7 °C at Kyrchyn.
The basin region contains 834 glaciers, each with areal coverage ranging from 0.1 to 11 km 2 .Based on Landsat8/OLI images (2015-2016), the total glacier coverage area was 121 km 2 in the northern part of the Teskey Range, but 77.1 km 2 in the southern part of the Kungoy Range.Glacier areas decreased by 12% in the Kungoy Range, and 8-12% in the Teskey Range from 1970-2000 [28][29][30][31].

Recent Variations of Glacial Lakes Using Satellite Data
To investigate the current state of glacial lakes in the Issyk-Kul Basin, we examined the seasonal changes of the glacial-lake area using Landsat7/Enhanced Thematic Mapper Plus (ETM+, SLC-off) and Landsat8/Operational Land Imager (OLI) images (in total 128 images, 68 for Teskey, 60 for Kungoy) taken in June, July, August, September, and October during 2013-2016.Here we examine glacial lakes of area exceeding 0.0005 km 2 that exist on debris landforms at glacier fronts.The glaciallake area was delineated manually as a polygon on pan-sharpened images.To map more recent glacier coverages and their debris landforms, we use Landsat8/OLI images from 2015-2016.
We classify the glacial lakes into six types according to their seasonal change in a given year.These types are "stable", "increasing", "decreasing", "appearing", "vanishing", and "short-lived".For a given lake, we determine its type according to its change in June, July, August, September, and October of each year."Stable" indicates that the lake area did not change seasonally from June to October of each year."Increasing" means that the area increased within a season, whereas

Recent Variations of Glacial Lakes Using Satellite Data
To investigate the current state of glacial lakes in the Issyk-Kul Basin, we examined the seasonal changes of the glacial-lake area using Landsat7/Enhanced Thematic Mapper Plus (ETM+, SLC-off) and Landsat8/Operational Land Imager (OLI) images (in total 128 images, 68 for Teskey, 60 for Kungoy) taken in June, July, August, September, and October during 2013-2016.Here we examine glacial lakes of area exceeding 0.0005 km 2 that exist on debris landforms at glacier fronts.The glacial-lake area was delineated manually as a polygon on pan-sharpened images.To map more recent glacier coverages and their debris landforms, we use Landsat8/OLI images from 2015-2016.
We classify the glacial lakes into six types according to their seasonal change in a given year.These types are "stable", "increasing", "decreasing", "appearing", "vanishing", and "short-lived".For a given lake, we determine its type according to its change in June, July, August, September, and October of each year."Stable" indicates that the lake area did not change seasonally from June to October of each year."Increasing" means that the area increased within a season, whereas "decreasing" means a decrease in area."Appearing" and "vanishing" types are lakes that appeared as a new lake, or disappeared due to shrinkage within a season, respectively.A "short-lived" type is a temporary lake that appeared or grew in the area, and then disappeared or shrank within the same year.A "short-lived lake" is a transient lake that experiences a sudden expansion and then, a few months later, a drainage.The short-lived type includes the "recurring" type that repeatedly appeared and disappeared at the same location during 2013-2016.Although the number of short-lived lakes is small in the study area, drainage of this type can be particularly dangerous.Examples include the Kashkasuu (2006), western Zyndan (2008), Jeruy (2013), and Karateke (2014) lakes [4,10].For the "recurring" type, we investigate the area variations of three glacial lakes (Kichi-Kyzylsuu, Koltor, and Tyup lakes in Figure 1) in the Teskey Range during 2013-2016 using Landsat8/OLI and Sentinel-2 images.As a result, we show the number of each lake type and the distribution of lake types.When we compare the number of lakes in a given type over different years, we refer to the changes as "number variability".The "distribution" means the pattern across all lake types within one range and one year, hereafter "type variability".
We compare this year-to-year number variability and yearly type variability with the summer temperature (May-October) anomaly and precipitation during 2013-2016 at the Tien Shan and Cholpon-Ata meteorological stations (Figure 1).We use 11-day running means of summer temperature anomalies.

Geomorphological Conditions of Glacial Lakes
To better understand the role of the regional geomorphological conditions that cause such large number and type variabilities, we closely examine the geomorphological changes in debris landforms at Jeruy Glacier, an area that was recently the source of a large drainage event.Our method involves a digital surface model (DSM) made with input from air photos (1979) and images from an unmanned aerial vehicle (UAV) in 2016.For georeferences, the DSMs use ground-control points (GCPs) from the Advanced Land Observing Satellite/Panchromatic Remote-sensing Instrument for Stereo Mapping (ALOS/PRISM) DSM data in 2007 (Universal Transverse Mercator (UTM) 43 with World Geodetic System (WGS) 84).The DSM differences of stable areas outside of the debris landforms are −3.7 ± 5.8 m (UAV DSM minus air photo DSM).The UAV image data comes from a Phantom 4 (DJI, Shenzhen, China).These data were converted to an orthoimage and DSM using the software Pix4D mapper with GCPs.The PRISM DSMs (2.5 m resolution) were processed by Japan Aerospace Exploration Agency (JAXA) Earth Observation Research Center (EORC) as a high-level product.The standard deviations of the PRISM DSM height errors (PRISM DSM without GCPs, minus reference DSM) are between 4.9 and 8.7 m [32,33].
Outlines of glacial lakes, and elevation of water levels were also measured during 2008-2017 using a Global Navigation Satellite System (GNSS) survey based on Trimble GeoExplorer 6000 (cm-edition).During this same period, we observed buried ice and ice tunnels inside of some debris landforms in the field.In our region, a short-lived glacial lake appears typically in a depression on a debris landform with buried ice [4].The existence of a depression (>0.01 km 2 ) as a potential lake basin is crucial for the formation of a short-lived lake [10].We identified such depressions (>0.001 km 2 ) on debris landforms in the northern part of the Teskey Range and southern part of the Kungoy Range using a filling model based on ALOS/PRISM DSMs.The same method was also applied to other glacial lakes [10,34].We used 14 scenes from 2007-2010 of ALOS/PRISM DSMs to determine depressions (Supplementary Table S1).We then examined the surface gradient and size of the depression.The surface gradient was calculated using the elevation (based on ALOS/PRISM DSMs) of glacier termini as derived from orthoimages, produced from Corona KH-4B (1971) and ALOS/PRISM (2007-2010) images.The Corona KH-4B data were converted to orthoimages using GCPs from ALOS/PRISM DSMs.For the GCPs, we used valley bottoms and big rocks.
According to whether or not the depression was in contact with the parent glacier, we classify it as either a depression (1) with glacier contact or (2) without glacier contact.

Radar Interferometry of Debris Landforms
A short-lived lake forms typically in a depression in a debris landform with buried ice [4,10].Mergili and others [35] also found that it is important to analyze the ground-ice distribution within the debris area.To understand the distribution of these debris landforms and their potential ground ice content, we first extracted areas of debris landforms at glacier fronts within ArcGIS using Landsat8/OLI and Google Earth data.Second, we ran differential interferometric SAR (DInSAR) analysis using PALSAR-2 (phased array type L-band synthetic aperture radar-2) data from the ALOS-2 satellite (Supplementary Table S2).The PALSAR-2 data came from an L-band (wavelength 23.8 cm) active 1.2-GHz radar.This dataset is well-suited for long-interval DInSAR analysis such as that spanning over two years, and is not significantly influenced by the small X-band-scale (wavelength 3.1 cm) vegetation cover in this region.Therefore, the PALSAR-2 enables DInSAR processing with more stable coherence than the X-band.The observation repetition of ALOS-2 is 14 days [36,37].The complex SAR data that we used has an original resolution of 9.1 m × 5.3 m (range × azimuth direction, nominal value).
We identify zones with ground ice as zones that show displacement fringes over a summer.For the Kyrgyz Range, the DInSAR data was verified using GNSS measurements [38].Previous, and similar, applications of DInSAR include the study of motion of rock glaciers in the Swiss Alps, Sierra Nevada (USA), and Iceland [39][40][41].Typically, displacement fringes on debris landforms indicate surface motion due to vertical subsidence and lateral motion.In contrast, bedrock, and grasslands or pasture areas in the valleys appear stable without any displacement.Here, the DInSAR analysis involves the GAMMA SAR software [42].The phase of the interferogram depends on orbital parameters, such as the length of the baseline and the imaging angle, as well as surface topography and surface displacement [43].We used image pairs with less than a 1000-m baseline.To remove the topographic phase, we used the SRTM DEM.Phase noise within the differential interferogram associated with temporal and spatial decorrelation was reduced by applying an adaptive filtering algorithm [44] implemented in the GAMMA software.We used 27 ALOS-2/PALSAR-2 images taken from 2014 to 2016 (Teskey Range: 20 images; Kungoy Range: 7 images).As a result, we defined debris landforms containing ground ice according to the existence of motion fringes in the DInSAR data.

Recent Changes of Glacial Lakes in the Issyk-Kul Basin
In total, we investigated 339 glacial lakes and found that the number of glacial lakes of each type differed.In the Teskey Range, 31 appearing and 40 increasing lakes occurred over just 2014 through 2015 alone (Figure 2A), whereas the Kungoy Range had only five appearing and six increasing lakes over all four years (Figure 2B).For the vanishing lakes, the Teskey Range had 42 in 2013, 16 in 2014, 19 in 2015, and 30 in 2016, whereas in the Kungoy Range, these lakes numbered just 12 in 2013, 7 in 2014, 10 in 2015, and 12 in 2016.For the short-lived lakes, 56 occurred in the Teskey Range and five occurred in the Kungoy Range over all four years.Many lakes changed their type every year through 2013-2016.Although recurring short-lived types appeared at the same location, the appearing, vanishing, and short-lived types appeared both at the same and different locations.
The type variability also differed between the two ranges.For example, in the Teskey Range, the year 2013 had a high proportion of decreasing (46 of 71) and vanishing lakes (42 of 107), but years 2014 and 2015 instead had a high proportion of appearing (31 of 51) and increasing (40 of 53) lakes, with a slightly higher percentage of these types in 2015 (Figure 2A).Although 2016 had only three increasing lakes, it had 22 short-lived cases, which is a higher percentage of this type compared to other years.
Compared to the Teskey Range, the Kungoy Range had a higher proportion of vanishing lakes in all four years (Figure 2B).Also, unlike the Teskey Range, the Kungoy Range had years without any increasing and years without any decreasing types.In contrast to the other types, these increasing and decreasing types are closely related to glacier variations.For example, a lake shrinks when the glacier separates from the lake due to glacier recession.When a glacial lake continues to have contact with the glacier under recession, the lake expands gradually.Four appearing lakes occur in 2014 and only one in 2016, whereas vanishing lakes total 41 through all four years.Of the short-lived lakes, four occur in 2013 and one in 2015.Over all four years, the total number of these other three types exceeds the total number of increasing and decreasing types.with the glacier under recession, the lake expands gradually.Four appearing lakes occur in 2014 and only one in 2016, whereas vanishing lakes total 41 through all four years.Of the short-lived lakes, four occur in 2013 and one in 2015.Over all four years, the total number of these other three types exceeds the total number of increasing and decreasing types.For these lake types, how do the number and type variabilities depend on summer temperature (May-October)?For the years 2013-2016, both the Tien Shan and Cholpon-Ata meteorological stations had similar trends of temperature anomalies, defined as temperature difference from the average over all 24 months (Figure 2C).Although the summer temperature anomaly in 2013 was comparable or lower than that in 2015, the distribution of types is more uniform in 2015 in the Teskey Range.In addition, the Teskey Range has many decreasing lakes when the summer temperature is average in 2013, whereas the Kungoy Range has none decreasing.Compared to 2014 and 2015, the type variability in the Teskey Range is about the same even though the summer temperature in 2015 is higher than that in 2014.In contrast, the type variability in the Kungoy Range is different in these years.
Also, despite both ranges having higher-than-average summer precipitation in both 2013 and 2015 (compared to 2014 and 2016, not shown here), the number of increasing lake type is relatively high only in the Teskey Range in 2015, whereas the number of decreasing lake type is relatively high in 2013.The numbers of appearing lakes show no clear trend with year.The same holds for the other types.Moreover, the distribution of lake types is not consistent from year to year.Thus, the results show that the large number and type variabilities over this short period are not directly related to the local short-term summer temperature anomaly, the precipitation, and the glacier recession.
The short-lived lake type, though much less common than the other types, is nevertheless significant.This type occurred 56 times in the Teskey Range over all four years versus just five in the For these lake types, how do the number and type variabilities depend on summer temperature (May-October)?For the years 2013-2016, both the Tien Shan and Cholpon-Ata meteorological stations had similar trends of temperature anomalies, defined as temperature difference from the average over all 24 months (Figure 2C).Although the summer temperature anomaly in 2013 was comparable or lower than that in 2015, the distribution of types is more uniform in 2015 in the Teskey Range.In addition, the Teskey Range has many decreasing lakes when the summer temperature is average in 2013, whereas the Kungoy Range has none decreasing.Compared to 2014 and 2015, the type variability in the Teskey Range is about the same even though the summer temperature in 2015 is higher than that in 2014.In contrast, the type variability in the Kungoy Range is different in these years.
Also, despite both ranges having higher-than-average summer precipitation in both 2013 and 2015 (compared to 2014 and 2016, not shown here), the number of increasing lake type is relatively high only in the Teskey Range in 2015, whereas the number of decreasing lake type is relatively high in 2013.The numbers of appearing lakes show no clear trend with year.The same holds for the other types.Moreover, the distribution of lake types is not consistent from year to year.Thus, the results show that the large number and type variabilities over this short period are not directly related to the local short-term summer temperature anomaly, the precipitation, and the glacier recession.
The short-lived lake type, though much less common than the other types, is nevertheless significant.This type occurred 56 times in the Teskey Range over all four years versus just five in the Kungoy Range.In the Teskey Range, they occurred every year, with 22 such lakes observed in 2016.Using Landsat8/OLI and Sentinel-2 images, we also identified three lakes of recurring type (for location see Figure 1) among these short-lived lakes from 2013 and 2016 (Figure 3A).This type changes size and shape every year.See Figure 3B for examples of the annual areal variations of the Koltor lake.The area changes of three lakes (Figure 3) also do not directly correlate with temperature variation.Lakes of the recurring type have different variations (timing and area) each year.Kungoy Range.In the Teskey Range, they occurred every year, with 22 such lakes observed in 2016.Using Landsat8/OLI and Sentinel-2 images, we also identified three lakes of recurring type (for location see Figure 1) among these short-lived lakes from 2013 and 2016 (Figure 3A).This type changes size and shape every year.See Figure 3B for examples of the annual areal variations of the Koltor lake.The area changes of three lakes (Figure 3) also do not directly correlate with temperature variation.Lakes of the recurring type have different variations (timing and area) each year.

Recent Surface Changes of Debris Landform at Jeruy Glacier Front
To consider another possible cause of the large variability in the number of each type and the distribution of types, we investigated the geomorphological conditions in the study area.First, we examined the surface changes of debris landforms at Jeruy Glacier front, as one of the short-lived type lakes, from 1979 (airphoto) through 2016 (UAV image).The debris landforms have no clear moraine ridges in either year.Between air photo and UAV images, several landforms have changed (Figure 4): many small ridges appeared on the debris landform, a depression on the left bank of the glacier front became unclear, the glacier terminus retreated by 150 m, a large depression formed where the glacier retreated.(A short-lived glacial lake appeared there and caused a large drainage on 15 August 2013 [10].)Also, an outlet channel expanded, with a longer section of exposed ice and small kettles (sinkholes) formed around the channel.We also observed much ice and an ice tunnel at some exposed parts in 2014.
Considering the profiles of the debris landform at the Jeruy Glacier front from 1979 to 2016, it can be seen that overall, the surface elevation steadily declined, with the depression expanding and lowering by up to 20 m at the glacier front (Figure 5A).The channel cleared by down-wasting and thawing as shown on the longitudinal profile along a-a' (Figure 5A) and on the transverse profile of the lake dam along b-b' (Figure 5B).

Recent Surface Changes of Debris Landform at Jeruy Glacier Front
To consider another possible cause of the large variability in the number of each type and the distribution of types, we investigated the geomorphological conditions in the study area.First, we examined the surface changes of debris landforms at Jeruy Glacier front, as one of the short-lived type lakes, from 1979 (airphoto) through 2016 (UAV image).The debris landforms have no clear moraine ridges in either year.Between air photo and UAV images, several landforms have changed (Figure 4): many small ridges appeared on the debris landform, a depression on the left bank of the glacier front became unclear, the glacier terminus retreated by 150 m, a large depression formed where the glacier retreated.(A short-lived glacial lake appeared there and caused a large drainage on 15 August 2013 [10].)Also, an outlet channel expanded, with a longer section of exposed ice and small kettles (sinkholes) formed around the channel.We also observed much ice and an ice tunnel at some exposed parts in 2014.
Considering the profiles of the debris landform at the Jeruy Glacier front from 1979 to 2016, it can be seen that overall, the surface elevation steadily declined, with the depression expanding and lowering by up to 20 m at the glacier front (Figure 5A).The channel cleared by down-wasting and thawing as shown on the longitudinal profile along a-a' (Figure 5A) and on the transverse profile of the lake dam along b-b' (Figure 5B).

Analysis of Depressions on Debris Landforms
Using ALOS/PRISM DSMs (2007-2010), we identified a total of 372 depressions (>0.001 km 2 ), with 328 in the Teskey Range and 44 in the Kungoy Range.By comparing them to glacier positions in 1971, we determined that 202 of the 372 had appeared due to glacier recession between 1971 and 2007-2010.Of these, 180 (of 328 in total) were in the Teskey Range and the remaining 22 (of 44 in total) in the Kungoy Range.Although the distributions of depression area by area class are roughly the same in both ranges (particularly for those without glacier contact), the Teskey Range hosts many more depressions than the Kungoy Range (Figure 6A).Another difference between the two ranges is that in the Teskey Range, nearly half (45%) of the depressions are of the glacier-contact type, whereas in the Kungoy Range, nearly all (80%) are of the type without glacier contact (Figure 6A).This difference is a result of the glacier shrinkage being larger in the Kungoy Range.Among all the

Analysis of Depressions on Debris Landforms
Using ALOS/PRISM DSMs (2007-2010), we identified a total of 372 depressions (>0.001 km 2 ), with 328 in the Teskey Range and 44 in the Kungoy Range.By comparing them to glacier positions in 1971, we determined that 202 of the 372 had appeared due to glacier recession between 1971 and 2007-2010.Of these, 180 (of 328 in total) were in the Teskey Range and the remaining 22 (of 44 in total) in the Kungoy Range.Although the distributions of depression area by area class are roughly the same in both ranges (particularly for those without glacier contact), the Teskey Range hosts many more depressions than the Kungoy Range (Figure 6A).Another difference between the two ranges is that in the Teskey Range, nearly half (45%) of the depressions are of the glacier-contact type, whereas in the Kungoy Range, nearly all (80%) are of the type without glacier contact (Figure 6A).This difference is a result of the glacier shrinkage being larger in the Kungoy Range.Among all the

Analysis of Depressions on Debris Landforms
Using ALOS/PRISM DSMs (2007-2010), we identified a total of 372 depressions (>0.001 km 2 ), with 328 in the Teskey Range and 44 in the Kungoy Range.By comparing them to glacier positions in 1971, we determined that 202 of the 372 had appeared due to glacier recession between 1971 and 2007-2010.Of these, 180 (of 328 in total) were in the Teskey Range and the remaining 22 (of 44 in total) in the Kungoy Range.Although the distributions of depression area by area class are roughly the same in both ranges (particularly for those without glacier contact), the Teskey Range hosts many more depressions than the Kungoy Range (Figure 6A).Another difference between the two ranges is that in the Teskey Range, nearly half (45%) of the depressions are of the glacier-contact type, whereas in the Kungoy Range, nearly all (80%) are of the type without glacier contact (Figure 6A).This difference is a result of the glacier shrinkage being larger in the Kungoy Range.Among all the depressions, 45 had existing lakes, but most of them were less than 0.01 km 2 , with the Teskey Range having only 13 larger lakes and the Kungoy Range having no larger lakes.We consider now the surface gradient at all glacier termini at which depressions appeared between 1971 and 2007-2010.In general, the smaller the gradient, the greater the number of depressions, with no depressions found on slopes of 15° or more (Figure 6B).Although all three surface-gradient bins include depression areas, the smaller surface gradients have more depressions.In the Teskey Range, many depressions occur on slopes less than 5°, including larger size depressions (>0.05 km 2 ).In the Kungoy Range, none of these larger depressions occur on slopes within 15°.Thus, the formation of depressions seems related to (i) the surface changes from melting ice on the debris landforms, (ii) the glacier shrinkage, and (iii) the glacier termini slope.

Distribution of Debris Landforms with Ground Ice
Ground ice occurrences were estimated using DInSAR images of the Jeruy Glacier front based on one-year interferograms between 4 October 2015 and 30 October 2016 (Figure 7A).The results show the displacement fringes on the debris landform.A displacement fringe indicates surface motion, which in turn, indicates ground ice.The stable regions (same color, no fringes) showing no motion are outside the debris landform area.These regions are riverbed, talus landform, and bedrock areas.The steep slope area at right-bank of the debris landform is labeled as the foreshortening part of the DInSAR images.The motion area seems mainly influenced by the melting of ground ice and subsequent subsidence as suggested by the vertical surface decline shown in Figure 5. Lateral movement due to gravity-driven deformation might also be occurring in the motion area.We consider now the surface gradient at all glacier termini at which depressions appeared between 1971 and 2007-2010.In general, the smaller the gradient, the greater the number of depressions, with no depressions found on slopes of 15 • or more (Figure 6B).Although all three surface-gradient bins include depression areas, the smaller surface gradients have more depressions.In the Teskey Range, many depressions occur on slopes less than 5 • , including larger size depressions (>0.05 km 2 ).In the Kungoy Range, none of these larger depressions occur on slopes within 15 • .Thus, the formation of depressions seems related to (i) the surface changes from melting ice on the debris landforms; (ii) the glacier shrinkage; and (iii) the glacier termini slope.

Distribution of Debris Landforms with Ground Ice
Ground ice occurrences were estimated using DInSAR images of the Jeruy Glacier front based on one-year interferograms between 4 October 2015 and 30 October 2016 (Figure 7A).The results show the displacement fringes on the debris landform.A displacement fringe indicates surface motion, which in turn, indicates ground ice.The stable regions (same color, no fringes) showing no motion are outside the debris landform area.These regions are riverbed, talus landform, and bedrock areas.The steep slope area at right-bank of the debris landform is labeled as the foreshortening part of the DInSAR images.The motion area seems mainly influenced by the melting of ground ice and subsequent subsidence as suggested by the vertical surface decline shown in Figure 5. Lateral movement due to gravity-driven deformation might also be occurring in the motion area.Using one-year pairs and short-period pairs between 2014 and 2016 (ALOS-2/PALSAR-2), we analyzed a total of 1110 debris landforms (Teskey Range: 930; Kungoy Range: 180) at glacier termini that show displacement fringes (Figure 7B).Among them, the area class of 0.1-1.0km 2 of debris landforms accounted for 80% of the total in the Teskey and 70% in the Kungoy Range.Overall, nearly half of the debris landforms at glacier fronts have this indication of ground ice; specifically, 413 of 930 in the Teskey Range and 71 of 180 in the Kungoy Range.

Regional Geomorphological Conditions of the Glacial Lakes in the Issyk-Kul Basin
Previous studies [22][23][24] stated that expansion of glacial lakes in the Djungarskiy Alatau Range and Tien Shan Mountains is related to climatic variations and glacier recession.Although an overall increase in lake area is influenced by glacial recession and climate changes over a decades-long timescale, our data shows that every year has large number and type variabilities.These seasonal and annual variations in the Teskey and Kungoy Ranges do not directly coincide with variations in summer temperature anomaly and precipitation, and the number of each type and the distribution of types also differ between the mountain ranges (Figure 2).Instead, the main reason for the occurrence of the large year-to-year variability in number and type appears to be local geomorphological conditions.These conditions include the existence of debris landforms containing ground ice and the formation and development of depressions.
An important influence on glacial lakes in the study region is their drainage through ice tunnels and ice dams.For example, Narama and others [4,10] reported that ice-tunnel drainage channels of four large drainages (2006-2014) close during winter and then open again in summer.They showed that this closure occurs by either freezing or by debris/ice deposition via channel collapse.Also, Mergili and others [35] pointed out that the behavior of glacial lakes in the region depends on whether the dam contains ice.In addition, other environmental factors affect these glacial lakes.For example, as shown for the Jeruy Glacier front (Figure 4), substantial down-wasting of debris landforms with buried ice causes surface changes.The subsidence and down-wasting of the debris landform changes how a depression forms and changes the ice-tunnel conditions inside of the landform.These surface changes might lead to the blockage and closure of ice tunnels due to deposition of ice and debris by tunnel collapse.Using one-year pairs and short-period pairs between 2014 and 2016 (ALOS-2/PALSAR-2), we analyzed a total of 1110 debris landforms (Teskey Range: 930; Kungoy Range: 180) at glacier termini that show displacement fringes (Figure 7B).Among them, the area class of 0.1-1.0km 2 of debris landforms accounted for 80% of the total in the Teskey and 70% in the Kungoy Range.Overall, nearly half of the debris landforms at glacier fronts have this indication of ground ice; specifically, 413 of 930 in the Teskey Range and 71 of 180 in the Kungoy Range.

Regional Geomorphological Conditions of the Glacial Lakes in the Issyk-Kul Basin
Previous studies [22][23][24] stated that expansion of glacial lakes in the Djungarskiy Alatau Range and Tien Shan Mountains is related to climatic variations and glacier recession.Although an overall increase in lake area is influenced by glacial recession and climate changes over a decades-long time-scale, our data shows that every year has large number and type variabilities.These seasonal and annual variations in the Teskey and Kungoy Ranges do not directly coincide with variations in summer temperature anomaly and precipitation, and the number of each type and the distribution of types also differ between the mountain ranges (Figure 2).Instead, the main reason for the occurrence of the large year-to-year variability in number and type appears to be local geomorphological conditions.These conditions include the existence of debris landforms containing ground ice and the formation and development of depressions.
An important influence on glacial lakes in the study region is their drainage through ice tunnels and ice dams.For example, Narama and others [4,10] reported that ice-tunnel drainage channels of four large drainages (2006-2014) close during winter and then open again in summer.They showed that this closure occurs by either freezing or by debris/ice deposition via channel collapse.Also, Mergili and others [35] pointed out that the behavior of glacial lakes in the region depends on whether the dam contains ice.In addition, other environmental factors affect these glacial lakes.For example, as shown for the Jeruy Glacier front (Figure 4), substantial down-wasting of debris landforms with buried ice causes surface changes.The subsidence and down-wasting of the debris landform changes how a depression forms and changes the ice-tunnel conditions inside of the landform.These surface changes might lead to the blockage and closure of ice tunnels due to deposition of ice and debris by tunnel collapse.
Degradation of a debris-covered glacier also causes large surface subsidence [44].For example, ground ice melting by thermokarst processes leads to the development of various landforms [45].Also, thermal erosion by lake water might enlarge a lake basin [46].In the Issyk-Kul Basin, these regional cryospheric conditions cause the observed large year-to-year variability in number and type (Figure 2).In short, the large number and type variabilities are due to regional geomorphological conditions involving the existence of debris landforms containing ground ice, the formation of depressions resulting from glacier recession, vertical subsidence on debris landforms, and surface changes of debris landforms due to ice melting.In addition, surface changes or the freezing of storage water might lead to the closure of an ice tunnel inside a debris landform due to the deposition of ice and debris by tunnel collapse.As a result of geomorphological conditions, tunnel closure and opening cause the observed appearing, vanishing, and short-lived lakes.
However, over decades-long periods, climate change and glacier recession affect the lake extension of stable or increasing lakes.Also, although the appearance of glacial lakes may be linked to climate and glacier recession, some of these lakes change to another type.In addition, the short-lived lake type occurs 56 times in the Teskey Range over all four years versus just five times in the Kungoy Range.However, only two large drainages occurred from the Jeruy (2013) and Karateke (2014) glacial lakes during 2013-2016.These large drainages might be influenced by the tunnel conditions.

Differences of Regional Geomorphological Conditions between the Teskey and Kungoy Ranges
Over the past decade, some large drainages have occurred from short-lived glacial lakes in the Teskey Range, but not the Kungoy Range.In the Kungoy Range, the most recent flood occurred from Choktal glacial lake on 27 September 1978, which drained through a dam incision [6].As a possible explanation for this difference between the two ranges, we will now describe differences between their regional geomorphological conditions.
Along the ridge line of both ranges, the rocks have similar origins, originating mainly from granites of the Paleozoic and Ordovician.In the eastern part of the Teskey Range, the rocks of the debris landforms partly originate from Paleozoic and Proterozoic gneisses and marbles.However, lake variability in the numbers of each type and their distributions are large in the same lithological area.This means that the rock types do not significantly influence either kind of lake variability.
Concerning the local short-term summer temperature anomaly and precipitation, both ranges have a similar summer temperature anomaly.Large year-to-year variability in their number and type occurred every year in both ranges.Thus, we argue that such lake variability over this short period (1-4 years) is instead strongly related to local geomorphological conditions, and is not directly influenced by the local short-term summer temperature anomaly and precipitation or glacier recession.
Consider the elevation ranges.Figure 8A shows the Teskey Range with glacial lakes at an elevation range of 3200-4100 m a.s.l. and debris landforms containing ice at 3200-4000 m a.s.l., giving an 800-m-elevation overlap between these geomorphic forms.For the Kungoy Range, the lakes occur at 3400-3900 m a.s.l. and the debris landforms with ice at 3100-3700 m a.s.l., giving a much smaller 300-m-elevation overlap (Figure 8B).The differences in altitudinal overlap reflect the differences of environmental condition such as climate and landform.In addition, the differences in the overlap also reflect the differences in the number of glacial lakes and depressions, which in turn are influenced by the debris landforms containing ice.The number of debris landforms with ice also differ, being 413 in the Teskey Range and 71 in the Kungoy Range.Thus, in the Teskey Range, the glacial lakes (depressions) form at the same elevations as the debris landforms with ice, whereas in the Kungoy Range, the glacial lakes have higher elevations than the debris landforms with ice.In addition, the number of debris landforms (black bar in Figure 8) at each overlap-range elevation is larger in the Teskey Range than in the Kungoy Range (49 vs. 4 at 3700 m and 80 vs. 7 at 3600 m).This shows that the Teskey Range has many places in which appearing, vanishing and shortlived lakes can occur.In fact, the number of appearing and short-lived lakes is quite large in the Teskey Range.This shows that the Teskey Range has favorable conditions for ice-tunnel closure.
The glacier shrinkage in both ranges influences the conditions of the depressions.Due to the smaller glacier shrinkage in the Teskey Range, a larger fraction of the depressions have contact with the parent glacier than that in the Kungoy Range (Figure 6A).Specifically, in the Teskey Range, 37% of all depressions are the glacier-contact type, but of just the middle-sized depressions (0.01-0.1 km 2 ), 60% are of the glacier-contact type.However, in the Kungoy Range, due to the greater glacial shrinkage, most depressions are separated from glacier termini.Contact with the glacier may be important for the formation of short-lived glacial lakes.In particular, we found that recent new depressions with glacier contact have a high probability of forming a short-lived glacial lake because such depressions have an inflow channel from the glacier and thus water supply to form lakes.In contrast, the behavior of depressions without glacier contact depends on the connection to meltwater channels from the glacier [10].But despite the relatively large amount of glacier shrinkage in the Kungoy Range, this range has a relatively small number of the increasing type.
Looking at other cases, Emmer and others [47] reported that proglacial lakes (i.e., glacier-contact type) turn to a glacier-detached phase (i.e., without glacier-contact type) in the Andes Mountains.The different scale of glacier shrinkage influences the formation or extension of depressions in our study region.For example, in the Teskey Range, a small glacier recession leads to an extension of the depression with glacier contact, whereas a large recession causes many separations between depressions and a larger separation from the glacier.As a result, and consistent with our findings, the Teskey Range has geomorphological conditions that are more favorable for depressions.
Thus, compared to the Kungoy Range, the Teskey Range has (1) more debris landforms with buried ice that developed along glacier fronts, (2) more large depressions with glacier contact that appeared within the past 40 years where a glacier retreated slowly, and (3) an elevation range of glacial lakes that coincides more with that of the debris landforms with buried ice.Although Narama and others [10] pointed out the geomorphological characteristics in which short-lived lakes appear, the present study stresses the distribution of geomorphological condition in the Issyk-Kul Basin, and the importance of differences in geomorphological conditions on regional scales between the Teskey and Kungoy Ranges.
As a result, the regional glacier and geomorphological changes found here impact lake variability in both number and distribution.In particular, the geomorphological changes over a short period (1-4 years) may change tunnel conditions and depression formation.These regional special In addition, the number of debris landforms (black bar in Figure 8) at each overlap-range elevation is larger in the Teskey Range than in the Kungoy Range (49 vs. 4 at 3700 m and 80 vs. 7 at 3600 m).This shows that the Teskey Range has many places in which appearing, vanishing and short-lived lakes can occur.In fact, the number of appearing and short-lived lakes is quite large in the Teskey Range.This shows that the Teskey Range has favorable conditions for ice-tunnel closure.
The glacier shrinkage in both ranges influences the conditions of the depressions.Due to the smaller glacier shrinkage in the Teskey Range, a larger fraction of the depressions have contact with the parent glacier than that in the Kungoy Range (Figure 6A).Specifically, in the Teskey Range, 37% of all depressions are the glacier-contact type, but of just the middle-sized depressions (0.01-0.1 km 2 ), 60% are of the glacier-contact type.However, in the Kungoy Range, due to the greater glacial shrinkage, most depressions are separated from glacier termini.Contact with the glacier may be important for the formation of short-lived glacial lakes.In particular, we found that recent new depressions with glacier contact have a high probability of forming a short-lived glacial lake because such depressions have an inflow channel from the glacier and thus water supply to form lakes.In contrast, the behavior of depressions without glacier contact depends on the connection to meltwater channels from the glacier [10].But despite the relatively large amount of glacier shrinkage in the Kungoy Range, this range has a relatively small number of the increasing type.
Looking at other cases, Emmer and others [47] reported that proglacial lakes (i.e., glacier-contact type) turn to a glacier-detached phase (i.e., without glacier-contact type) in the Andes Mountains.The different scale of glacier shrinkage influences the formation or extension of depressions in our study region.For example, in the Teskey Range, a small glacier recession leads to an extension of the depression with glacier contact, whereas a large recession causes many separations between depressions and a larger separation from the glacier.As a result, and consistent with our findings, the Teskey Range has geomorphological conditions that are more favorable for depressions.
Thus, compared to the Kungoy Range, the Teskey Range has (1) more debris landforms with buried ice that developed along glacier fronts; (2) more large depressions with glacier contact that appeared within the past 40 years where a glacier retreated slowly; and (3) an elevation range of glacial lakes that coincides more with that of the debris landforms with buried ice.Although Narama and others [10] pointed out the geomorphological characteristics in which short-lived lakes appear, the present study stresses the distribution of geomorphological condition in the Issyk-Kul Basin, and the importance of differences in geomorphological conditions on regional scales between the Teskey and Kungoy Ranges.
As a result, the regional glacier and geomorphological changes found here impact lake variability in both number and distribution.In particular, the geomorphological changes over a short period (1-4 years) may change tunnel conditions and depression formation.These regional special conditions may induce short-lived lake appearances due to glacier/debris/ice instability.Thus, short-lived lakes that cause large drainages should be closely monitored.

Conclusions
For the Teskey and Kungoy Ranges of the Issyk-Kul Basin, we described the seasonal area changes of 339 glacial lakes.These lakes were classified into six types with respect to area changes: stable, increasing, decreasing, appearing, vanishing, and short-lived.Lakes in the Teskey Range varied by number and type significantly more than those in the Kungoy Range.In the Teskey Range, 31 appearing and 40 increasing lakes occurred in just 2014 and 2015, whereas only five appearing and six increasing lakes in the Kungoy Range were identified over four years.For the vanishing lakes, the Teskey Range had 42, 16, 19, and 30 in each of the four years, respectively, whereas the Kungoy Range had just 12, 7, 10, and 12.For the short-lived lakes, the Teskey Range had 56 over the four years, whereas the Kungoy Range had just five.The large year-to-year variabilities in both number and distribution were not directly related to local short-term summer temperature anomaly, nor to precipitation or glacier recession during 2013-2016.
We investigated the regional geomorphological conditions that caused these large number and type variabilities using field observation and satellite data analysis.We argued that these variations arise from regional geomorphological conditions involving the existence of debris landforms containing ground ice, the formation of depressions resulting from glacier recession, vertical subsidence on debris landforms, and surface changes of debris landforms due to melting of the ice content.In addition, surface changes or the freezing of storage water might lead to the closure of an ice tunnel inside a debris landform due to the deposition of ice and debris by tunnel collapse.Three types of appearing, vanishing, and short-lived lakes (including recurring ones) are related to these regional geomorphological conditions, and not directly related to the local short-term summer temperature anomaly, the precipitation, and glacier recession.
Compared to the Kungoy Range, the Teskey Range has recent large drainage events.This latter range has glacial lakes/depressions occurring at the same elevations as debris landforms containing ground ice, a condition leading to a large ratio of depressions with glacier contact as well as a large number of depressions and debris landforms containing ground ice.These forms and conditions are more conducive to appearing, vanishing, and short-lived lakes, which we suggest explains the greater number observed in the region.
For example, the average annual precipitation values during 1981-1990 were 594 mm at the Chong-Kyzylsuu station (2855 m) (data for Chong-Kyzylsuu station averaged only for 1981-1987), 297 mm at the Tien Shan station (3614 m), 307 mm at the Cholpon-Ata station (1645 m), and 625 mm at the Kyrchyn station (2100 m).See Figure 1 for locations.Over the same period, the annual average temperature was 0.2 • C at Chong-Kyzylsuu, −7.6 • C at Tien Shan, 8.1 • C at Cholpon-Ata, and 2.7 • C at Kyrchyn.The basin region contains 834 glaciers, each with areal coverage ranging from 0.1 to 11 km 2 .Based on Landsat8/OLI images (2015-2016), the total glacier coverage area was 121 km 2 in the northern part of the Teskey Range, but 77.1 km 2 in the southern part of the Kungoy Range.Glacier areas decreased by 12% in the Kungoy Range, and 8-12% in the Teskey Range from 1970-2000 [28-31].

Figure 1 .
Figure 1.The Teskey and Kungoy Ranges in the Issyk-Kul Basin, Kyrgyzstan, and locations of existing glacial lakes (2016).

Figure 1 .
Figure 1.The Teskey and Kungoy Ranges in the Issyk-Kul Basin, Kyrgyzstan, and locations of existing glacial lakes (2016).

Figure 2 .
Figure 2. Number of lakes in the six types and the summer temperature anomaly in the Issyk-Kul Basin during 2013-2016.(A) Teskey Range.(B) Kungoy Range.The value next to the grey bar shows the number of "stable" lakes.(C) Eleven-day running means of summer temperature (May-October) anomaly during 2013-2016 at the Tien Shan (Teskey Range) and Cholpon-Ata (Kungoy Range) meteorological stations (for locations, see Figure 1).For these six months, the average temperatures of the Tien Shan meteorological station were 2.4 °C in 2013, 1.5 °C in 2014, 2.5 °C in 2015, and 2.1 °C in 2016, whereas those for the Cholpon-Ata meteorological station were 15.8 °C in 2013, 15.5 °C in 2015, and 15.3 °C in 2016.

Figure 2 .
Figure 2. Number of lakes in the six types and the summer temperature anomaly in the Issyk-Kul Basin during 2013-2016.(A) Teskey Range.(B) Kungoy Range.The value next to the grey bar shows the number of "stable" lakes.(C) Eleven-day running means of summer temperature (May-October) anomaly during 2013-2016 at the Tien Shan (Teskey Range) and Cholpon-Ata (Kungoy Range) meteorological stations (for locations, see Figure 1).For these six months, the average temperatures of the Tien Shan meteorological station were 2.4 • C in 2013, 1.5 • C in 2014, 2.5 • C in 2015, and 2.1 • C in 2016, whereas those for the Cholpon-Ata meteorological station were 15.8 • C in 2013, 15.5 • C in 2015, and 15.3 • C in 2016.

Figure 3 .
Figure 3. Recurring-type of short-lived lakes in the Teskey Range.(A) Inter-annual areal variations of the Kichi-Kyzylsuu, Koltor, and Tyup lakes (for locations, see Figure 1) between 2013 and 2016.(B) Images of the Koltor lake in 2010 (photo from fieldwork by authors), 2014, 2015, and 2017 (photos from helicopter by authors) with the boundaries highlighted.

Figure 3 .
Figure 3. Recurring-type of short-lived lakes in the Teskey Range.(A) Inter-annual areal variations of the Kichi-Kyzylsuu, Koltor, and Tyup lakes (for locations, see Figure 1) between 2013 and 2016.(B) Images of the Koltor lake in 2010 (photo from fieldwork by authors), 2014, 2015, and 2017 (photos from helicopter by authors) with the boundaries highlighted.

Figure 4 .
Figure 4. Surface changes of the debris landform at Jeruy Glacier between 1979 and 2016.Data based on orthoimages from air photos (1979) and digital images taken by a UAV (2016).

Figure 5 .
Figure 5. Elevation profiles from 1979 and 2016 of debris landforms at the Jeruy Glacier.(A) Longitudinal profile along a-a' marked in the 2016 panel of Figure 4. (B) Transverse profile along bb' from the same panel.Based on DSMs generated from air photos (1979) and UAV images (2016).

Figure 4 .
Figure 4. Surface changes of the debris landform at Jeruy Glacier between 1979 and 2016.Data based on orthoimages from air photos (1979) and digital images taken by a UAV (2016).

Figure 4 .
Figure 4. Surface changes of the debris landform at Jeruy Glacier between 1979 and 2016.Data based on orthoimages from air photos (1979) and digital images taken by a UAV (2016).

Figure 5 .
Figure 5. Elevation profiles from 1979 and 2016 of debris landforms at the Jeruy Glacier.(A) Longitudinal profile along a-a' marked in the 2016 panel of Figure 4. (B) Transverse profile along bb' from the same panel.Based on DSMs generated from air photos (1979) and UAV images (2016).

Figure 5 .
Figure 5. Elevation profiles from 1979 and 2016 of debris landforms at the Jeruy Glacier.(A) Longitudinal profile along a-a' marked in the 2016 panel of Figure 4. (B) Transverse profile along b-b' from the same panel.Based on DSMs generated from air photos (1979) and UAV images (2016).

Geosciences 2018, 8 ,
x FOR PEER REVIEW 9 of 16 depressions, 45 had existing lakes, but most of them were less than 0.01 km 2 , with the Teskey Range having only 13 larger lakes and the Kungoy Range having no larger lakes.

Figure 6 .
Figure 6.Characteristics of depressions in the Teskey and Kungoy Ranges.(A) The number of depressions with and without glacier contact of the various area classes.(B) The number of depressions of various surface-gradient ranges for each area class.

Figure 6 .
Figure 6.Characteristics of depressions in the Teskey and Kungoy Ranges.(A) The number of depressions with and without glacier contact of the various area classes.(B) The number of depressions of various surface-gradient ranges for each area class.

Figure 7 .
Figure 7. Surface motion around the Jeruy Glacier front.(A) Surface displacement map by DInSAR analysis using ALOS-2/PALSAR-2 data.(B) The resulting number of debris landforms with and without ground ice (or only little ground ice without significant effect on surface deformation) over the study region.

Figure 7 .
Figure 7. Surface motion around the Jeruy Glacier front.(A) Surface displacement map by DInSAR analysis using ALOS-2/PALSAR-2 data.(B) The resulting number of debris landforms with and without ground ice (or only little ground ice without significant effect on surface deformation) over the study region.

Figure 8 .
Figure 8. Elevation ranges for the zones of glacial lake depressions and debris landforms with buried ice in the Teskey (A) and Kungoy (B) Ranges.

Figure 8 .
Figure 8. Elevation ranges for the zones of glacial lake depressions and debris landforms with buried ice in the Teskey (A) and Kungoy (B) Ranges.