Cryosphere hazards in Ladakh: the 2014 Gya glacial lake outburst flood and its implications for risk assessment

This article attempts to reconstruct the causes and consequences of the 2014 glacial lake outburst flood (GLOF) event in Gya, Ladakh. We analyse the evolution of the Gya glacial lake using a high temporal and high spatial resolution remote sensing approach. In order to frame the case study in a larger picture, we produce a comprehensive inventory of glacial lakes for the entire Trans-Himalayan region of Ladakh. Changes in the extent and number of glacial lakes have been detected for the years 1969, 1993, 2000/02 and 2018 in order to assess the potential risk of future GLOFs in the region. The remote sensing approach was supported by field surveys between 2014 and 2019. The case study of the Gya GLOF illustrates the problem of potentially hazardous lakes being overlooked in inventories. The broader analysis of the Ladakh region and in-depth analysis of one GLOF lead us to propose an integrated approach for detecting undocumented GLOFs. This article demonstrates the necessity for using multiple methods to ensure robustness of risk assessment. The improved understanding can lead to a more accurate evaluation of exposure to cryosphere hazards and identification of alternative mechanisms and spatial patterns of GLOFs in the Himalaya.


Introduction
Cryosphere dynamics and cryosphere-related hazards are a vital component of and threat for land use development in the semi-arid Trans-Himalayan region of Ladakh. The observed and projected glacier decrease in the extended Himalayan region (Bolch et al. 2019) and specifically in Ladakh (Schmidt and Nüsser 2017) will affect livelihood security, both in the mountains and in the lowlands in the long run (Hock et al. 2019;Huss and Hock 2018;Immerzeel et al. 2010;Nüsser et al. 2019). In the short term, glacial lake outburst floods (GLOFs) are major cryosphere hazards and a potential risk for local communities 1 3 as lakes are formed and grow due to deglaciation (Carey et al. 2017;Carrivick and Tweed 2016;Emmer and Vilímek 2013;Emmer 2018;Rasul and Molden 2019).
On 6 August 2014, a glacial lake outburst flood (GLOF) hit the village of Gya (33° 39′ N, 77° 44′ E) in Ladakh. GLOF events are not a new phenomenon in the upper Indus basin: a number of such events have been reported from the Karakoram and Western Himalaya dating back to the nineteenth century (Hewitt 1982;Hewitt and Liu 2010;Iturrizaga 2019). Large and recurrent GLOF events happened in the Shyok and Nubra valley, where the river was blocked by a glacier several times during the nineteenth century, followed by destructive outburst floods (Bhambri et al. 2019(Bhambri et al. , 2020Cunningham 1854;Sheikh 2015;Sinclair 1929).
However, only after the GLOF of Dig Tsho, a glacial lake in the Khumbu region of Nepal in 1985 (Vuichard and Zimmermann 1987) was there a rise of scientific interest in glacial lake dynamics and potential risk of outburst floods in the Himalaya. Since then, various studies have analysed GLOF events (Allen et al. 2016b;Ashraf et al. 2012;Nie et al. 2018) and their potential risk (Allen et al. 2016a;Sharma et al. 2018;Wang et al. 2012) in different parts of the Hindu Kush-Karakoram-Himalaya (HKH). Contrary to the findings of Richardson and Reynolds (2000) that GLOF events have become more frequent, Komori et al. (2012) and Veh et al. (2019) could not detect an increase in GLOF events in the central and eastern Himalaya over the last few decades. However, a global analysis of GLOFs suggests that they demonstrate a lagged response to climate change and they are likely to become more frequent over the next decades, lasting into the twenty-second century . Although it is difficult to give precise figures due to differing data sets and methodology, all studies show an increase in the number and area of glacial lakes in the HKH (Bolch et al. 2019;Ives et al. 2010;Nie et al. 2017;Zhang et al. 2015). The largest absolute glacial lake growth was observed in the central and eastern Himalaya (Bajracharya and Mool 2009;Gardelle et al. 2011;Haritashya et al. 2018).
Only a few destructive GLOFs have been reported from the Trans-Himalayan region of Ladakh during recent decades (Ikeda et al. 2016;Morup 2014;Narama et al. 2012; Tabassum and Kanth 2013). One GLOF was described from Nyoma in the Changthang area in eastern Ladakh in 1971, which caused 13 to 16 fatalities. Another GLOF happened in Domkhar in 2003, which destroyed farmland and infrastructure (Ikeda et al. 2016). The most recent reported GLOF is the 2014 Gya flood (Dolma 2014). This paper attempts to reconstruct the causes and consequences of the GLOF event at Gya village and to analyse the evolution of the glacier lake using a high temporal and high spatial remote sensing approach with data from 1969 to 2019. For the entire region of Ladakh, a multi-temporal and multi-scale remote sensing approach (Corona, Landsat, Sentinel) was conducted to investigate changes of glacial lakes and GLOF events. Based on this new inventory of glacier lakes in Ladakh, probable former GLOF events and the potential risk of future GLOFs are analysed. Identification of unreported or undocumented GLOFs can lead to a more accurate assessment of exposure and risk of cryosphere hazards. The remote sensing approach was supported by field surveys between 2014 and 2019.

Study area
The Trans-Himalayan region of Ladakh is situated between the Greater Himalayan Range to the south and the Karakoram Range to the north at an altitude of over 3000 m a.s.l., with surrounding mountain ranges exceeding altitudes of 5500 m a.s.l. Rain shadow effects 1 3 result in steep hygric gradients of decreasing precipitation from west to east and south to north (Hussain Dahri et al. 2016). Despite the general cold-arid conditions, occasional torrential rainstorms may cause floods, such as the extreme cloudburst in August 2010 (Thayyen et al. 2013;Ziegler et al. 2016). The specific combination of topography and climate contributes to the small size and high altitude of glaciers, almost all of which terminate above 5200 m a.s.l. and are smaller than 0.75 km 2 . Glaciated area reduction ranges from 0.2 to 0.9% year −1 for the period between 1969 and 2016, with high variability across different watersheds (Schmidt and Nüsser 2017). Meltwater from these small glaciers and permafrost thaw (Ali et al. 2018) determines the potential for habitation and irrigated crop cultivation , only found along the Indus and Shyok rivers, on alluvial fans, or tributary valleys between 2200 m a.s.l. and 4370 m a.s.l. (Fig. 1). The region is sparsely populated and has a total of about 274,300 inhabitants (Leh and Kargil districts, since 2019 administered as a Union Territory), with the registered population expanding by 1.5% annually between 2001 and 2011 (Census of India 2011).
Gya, the case study village, is located at an altitude of 4010 m a.s.l. on the eastern slope of the Kang Yatze massif, close to the Taglang La (5328 m a.s.l.) on the main road connecting Leh and Manali. At the end of the nineteenth century, the village together with neighbouring hamlets consisted of 40 houses and a monastery and served as a halting place with grazing grounds for animals along the historical trading route (Gazetteer of Kashmir and Ladák 1890). The old settlement is located along a small ridge between a palaeo-fluvial terrace, which is used for agriculture, and the current channel of Gya rivulet (Fig. 2). A second small kernel is located above the irrigated area on the southern side of the rivulet, which divides the cultivated area into a northern and southern part. Both built-up areas have been constructed on flood-protected places. The village has grown to more than 700 inhabitants (Government of Jammu & Kashmir 2011-2012) dispersing into the cultivated area of 89 ha (Fig. 3). Only a few houses were built in the vicinity of the stream. The upper catchment of the 16-km-long Gya valley is characterized by glaciers and perennial snow fields above 5400 m a.s.l., which cover an area of about 2.9 km 2 and feed the irrigated fields of the village. The Gya glacier affected by the proglacial lake was identified as receding very fast compared to neighbouring glaciers of the Kang Yatze massif .

Lake inventory of the Trans-Himalaya
Periglacial and glacial lakes have been mapped on Corona images from 1969, Landsat imagery from 1993, 2000/2002, and Sentinel 2 data from 2018 (Table 1) in order to analyse the evolution of lakes and to assess future risk. Temporally adjoining images were added to fill data gaps caused by cloud coverage.
The delineation of lakes in radiometric calibrated and co-registered Sentinel and Landsat images was carried out using a standardized semi-automatic approach based on a Normalized Difference Water Index (NDWI) according to the formula (NIR -B)/ (NIR + B). The usage of near-infrared and blue band representing maximum and minimum reflectance for water gives good discrimination between ice and snow (Allen et al. 2016b;Bolch et al. 2012;Huggel et al. 2002;Worni et al. 2013). Misclassified pixels were removed using slope angle and shadow masks. Water bodies below 4000 m a.s.l.
were deleted because glacial lakes do not exist at these altitudes under given climatic and environmental conditions. The remaining classified lake pixels were converted to vector data for automatic calculation of lake size. All lakes smaller than 9000 m 2 were removed, because their change detection was regarded as too uncertain given the spatial resolution of Landsat images used for the years 1993 and 2002. Uncertainty in glacial lake area was estimated to be half a pixel (Kamp et al. 2011). In the final pre-processing step, all lakes were controlled by visual interpretation, remaining misclassified water bodies were removed, snow-covered and frozen lakes were added, and some lakes, which were partly ice-covered, rich in sediment load or located in shadow positions, were corrected. On orthorectified panchromatic Corona images, the lakes were digitized on screen in a geographical information system (ArcGIS 10.6) based on detailed visual image interpretation. Additionally, lakes smaller than 9000 m 2 , which can be detected on at least one of the other three time periods, were digitized in order to investigate since when these lakes have existed or if they have ever completely discharged. Finally, each lake was categorized as proglacial lake (PG), supraglacial lake, lake in a recent moraine (RM), lakes in Pleistocene moraine, cirque lake and landslide dammed lake. Change detection analyses were conducted for three distinct time spans : 1969-1993, 1993-2000/02 and 2000/02-2018. To estimate glacial lake volume, several empirical formulae exist (Cook and Quincey 2015;Muñoz et al. 2020). For this study, the water volume for all time spans was calculated according to the widely accepted empirical formula by Huggel et al. (2002): The time series analyses of Gya indicated that volume changes up to 100,000 m 3 are within the seasonal range; a threshold of more than 200,000 m 3 was used to detect potential GLOFs.

Gya lake and the 2014 Gya GLOF event
For the in-depth study of the GLOF event in Gya on 6 August 2014, high spatial resolution and multi-spectral images were used to analyse lake evolution, lake level changes before and after the event, changes of the moraine for further risk assessment, and size and damages of the flood.

Changes in the proglacial lake
In order to analyse the lake evolution, multi-temporal and multi-scale remote sensing imagery from different satellites preferably taken in summer or autumn were used to generate a long time series. Starting with Corona images from 1969, available cloud-and snow-free Landsat imagery from 1989 to 2014 and Sentinel 2, RapidEye and PlanetScope data from 2014 to 2019 were used (Table 2). Due to the growing availability of free and open satellite imagery, the temporal and spatial resolution has improved, thereby reducing uncertainty.
The outline of the proglacial lake was digitized manually on-screen on all selected images. The semi-automatic approach, which was used for the lake detection of the Trans-Himalaya, was not possible, due to the frequent ice coverage of this lake. For error estimation using the root-mean-square error, we used half a pixel size as a buffer (Kamp et al. 2011).

Investigation of the 2014 Gya GLOF event
The 2014 Gya GLOF was investigated using high spatial resolution images taken before and after the event. The Pléiades image, taken on the morning of 6 August 2014, allowed the detection of the maximum lake level a few hours before the GLOF hit Gya village.
To estimate lake level changes, the altitude of the lake before and after the event has been calculated. Therefore, on stable ground-excluding the glaciated area-vertices of the digitized outlines of the lake from 6 August 2014 and 13 October 2014 have been used to derive the elevation from TanDEM-X data, processed and delivered by German Aerospace Center (DLR) from interferometry data taken between 2011 and 2013. The calculated difference between the averaged elevations of both time steps represents lake level fluctuations. This height difference-multiplied by the size of the lake area-was used to estimate the water discharge by the GLOF. Additionally, the released water volume was derived from the difference between the estimated lake volume before and after the GLOF event according to the above-mentioned empirical formula by Huggel et al. (2002). The same approach was used to estimate the released water volume in 2019. We estimated the peak discharge (Q max ) with a water volume V (km 3 ) by the empirical formula by Walder and Costa (1996) for tunnel events: The flooded area was reconstructed conducting a change detection analyses of two radiometric calibrated RapidEye images (2013/11/26 and 2014/10/13) representing the situation before and after the GLOF event. For this purpose, an unsupervised classification of the combined Normalized Difference Vegetation Index (NDVI) derived from both images was conducted within a distance of 300 m to the stream. Ten classes were separated and later merged to differentiate flooded and non-flooded areas. The classified flooded area was corrected by visual interpretation. The flood map was used as a mask to quantify the damages in the village, for which the fields, houses and infrastructure were mapped on an orthorectified SPOT-6 image taken on 18 November 2013. During the field survey in September and October 2014, the proglacial lake was prospected in order to reconstruct the GLOF event, to map the extent of the flood and to assess the potential for further GLOF events. Furthermore, the flooded area and damages in the village were mapped. Semi-structured interviews were conducted with villagers in order to get information on the 2014 GLOF and former floods. Further field visits in 2015 and 2017 only in the village environs and a final survey in August 2019 were used to inspect the situation at the proglacial lake and recent reconstruction in the village.

Results
The results are divided into two sections, the first on the glacial lake inventory for the entire Ladakh region and the second on the specific case of the Gya GLOF of 2014.

Glacial lakes in the Trans-Himalaya of Ladakh
A total of 192 glacial lakes cover an area of 5.93 ± 0.70 km 2 with an estimated water volume of about 61.11 ± 8.5 million m 3 , including 127 proglacial (PG) and 9 supraglacial lakes and 56 lakes located on recent moraines (RM) in 2018 (Fig. 4). Apart from these, 30 lakes are located in front of recent moraines (combined area: 1.53 ± 0.14 km 2 ) and 14 lakes in cirques (combined area: 0.47 ± 0.04 km 2 ). Corresponding to the location of glaciers, the altitude of glacial lakes increases from west to east from below 4900 m a.s.l. to above 5500 m a.s.l. The total number of lakes is higher in the western part compared to the east and higher along the Ladakh range compared to the Zanskar range. Furthermore, in the western Ladakh range the highest amount of lakes in pleistocene moraines can be detected (in total 129 lakes with an area of 6.83 ± 0.55 km 2 ), which can be dangerous in the context of cascading flood events. In 2018, most glacial lakes were smaller than 50,000 m 2 (equivalent to 490,000 m 3 ) with only 10 larger than 100,000 m 2 (equal to approximately 1.3 million m 3 , Fig. 4). The area of the two largest glacial lakes is about 245,000 m 2 (equal to approximately 4.7 million m 3 ); one of them is located in eastern Ladakh, 24 km to the south of Pangong Lake, and the second one is located in the western part of the study area, close to the confluence of Shyok and Indus river.
The number and size of PG and RM lakes have increased over time. On the 1969 Corona images, at least 73 glacial lakes larger than 9000 m 2 are already detectable. An additional 42 glacial lakes smaller than 9000 m 2 have been taken into consideration. As the coverage of Corona data does not include the entire study area, detailed information on 44 glacial lakes detected in later images is not available for 1969. However, in 1993, ten of the undetectable lakes did not yet exist and only two were larger than 100,000 m 2 .
The change detection analyses indicated that 22 glacial lakes disappeared (decrease by more than 90%) between 1969 and 2018: four between 1969 and 1993, nine between 1993 and 2002, and nine from 2002 to 2018. Nine lakes can be detected which released more than 200,000 m 3 of water within each period; only two of them completely disappeared between 2002. From 1969to 1993and 1993to 2002, most of these lakes were located in the sparsely populated Changthang region in eastern Ladakh (Fig. 5).
The number of lakes with a water volume increase of more than 200,000 m 3 has doubled from 14 lakes in the period between 1993 and 2002 to 29 between 2002 and 2018. (Only one lake, which is not covered by Corona images, was larger than 200,000 m 3 in 1993.) The time series also shows that some of the former released lakes were refilled in the aftermath, such as in the Changthang area and Zanskar range. Many expanding lakes can be found in the western part of the Ladakh range, where the majority of habitations are located.

Evolution of the proglacial lake
The proglacial Gya lake has developed in front of a cirque glacier at an altitude of 5400 m a.s.l. (Fig. 6). Glacier retreat and corresponding lake expansion of about 430 m can be detected for the period between 1969 and 2019. The year of initial lake formation remains unclear. Surrounded by steep slopes, this lake was largely unknown due to its location behind a 400-m-wide and 100-m-high moraine. Only a few aged shepherds, who had to look after their yaks on the high pasture in their younger days, had knowledge of its existence (own interviews). The terminal moraine with a small thermokarst hollow, which was water-filled (100 m 2 ) in August 2014, is cut by a channel, connected with the proglacial lake by a 6-m-deep depression. The freeboard amounts to about 10 m; thus, the lake discharge occurs subsurface through the moraine and water only reappears at the surface at an altitude of 5300 m a.s.l.
The time series analyses (Fig. 7) show that the proglacial lake is regularly ice-covered even in summer and has expanded continuously from 0.03 ± 0.002 km 2 (with an estimated water volume of 260,900 ± 17,900 m 3 ) in July 1969 to 0.12 ± 0.003 km 2 (estimated water volume: 1,648,050 ± 67,800 m 3 ) in August 2019. Since 2005, a rapid increase in the lake area can be detected, which has almost doubled from 0.06 ± 0.007 km 2 to 0.12 ± 0.003 km 2 over the past 15 years. In 2014, a drastic short-term lake level rise occurred with an area increase of about 0.02 km 2 to 0.11 ± 0.002 km 2 (estimated water volume of 1,459,250 ± 30,200 m 3 ) from July 2014 to 6 August 2014, the day of the flood. After the GLOF event, the lake area decreased to almost its former size of about 0.09 ± 0.009 km 2 . The estimated released water volume, derived from the difference of the lake volume using satellite data from shortly before and after the event, using the empirical lake volume formula (Huggel et al. 2002) amounts to 360,650 ± 123,600 m 3 . The peak discharge (Q max ) amounts to 23 ± 6 m 3 s −1 . However, considering the lake area and lowering of lake level by about 5 m (derived from TanDEM-X (Fig. 8) and verified by field observation (Fig. 6), the estimated discharge ranges from 450,000 to 600,000 m 3 with an estimated peak discharge between 27 and 32 m 3 s −1 . A lake level fluctuation of a similar magnitude can also be detected in September 1998, when the lake size was reduced to 81% of its previous extent, releasing an estimated water volume of 237,600 ± 53,500 m 3 . Due to the insufficient temporal resolution of satellite images from that period, a GLOF event cannot be stated with confidence; however, according to some elderly villagers, a flood without rain in the catchment, inundating some fields, occurred three decades ago (own interviews, September 2014). Recent rapid lake level rise can be detected in August 2019, when the lake reached the same level as in 2014 before it decreased to its former level between August and September (Fig. 8). The estimated released water volume was 155,900 ± 9200 m 3 (according to the empirical formula by Huggel et al. 2002) and 204,000 to 219,000 m 3 according to lake level differences; however, a flood was not observed in the village (personal communication Stangzin Dorjei Gya, 20.02.2020) and was not evident during the field survey in late 1 3 August 2019. Thus, a slow water release must be assumed, also supported by the image from 11 September, when the lake level already began to sink.  1969-1993, 1993-2002 and 2002-2014 in the Trans-Himalayan region of Ladakh

Flood event in 2014
Field observations in September 2014 show that the moraine between the proglacial lake and the thermokarst hollow was not eroded and no fluvial sediments were accumulated. Fine-grained sediments with some small ripple marks were only found to a limited extent near the thermokarst hollow (Fig. 8). One individual plant (Thylacospermum caespitosum) was found in the channel, whereas sandy sediments were absent (Fig. 9). Furthermore, the satellite image taken at 11 a.m. on the day of the GLOF event displays a high water Landsat images were used, from 2014 on these were supplemented by 12 Sentinel and 7 PlanetScope data with higher spatial resolution (plus 9 images from various sensors); only on 6 satellite images, the lake surface was completely ice free, and on 14 images, some ice blocks could be detected on the lake. On all other images, the lake is totally or partly covered by ice discharge at around 5300 m a.s.l. (Fig. 8). These observations suggest that not a spillover but rather a subsurface (tunneling) drainage process occurred, which resulted from thawing of ice cores in the moraine. The observed rapid increase in the lake level can be attributed to increased snow and glacier melting or temporary blockage of the tunneling system. After the outburst, the lake reverted to its previous level as mapped for 2013 using RapidEye data. The Pléiades image from 6 August 2014 and the TanDEM-X indicate a reduction in Fig. 8 Proglacial lake of Gya-the lake levels before and after the GLOF as well as the lake level changes in 1998 and 2019 are shown. The contour lines derived from TanDEM-X indicate the depression and desiccated channel cutting the moraine. In front of the moraine, the large amount of melt water is detectable 1 3 vertical distance between the depression and the lake surface of about 2 m. Remote sensing imagery and the field survey in August 2019 showed a similar situation regarding lake level and extent (Fig. 10).
Downstream of the terminal moraine, drainage occurs in several river channels, but a sediment tail could be detected neither on remote sensing data nor during the field survey. The run-off was dammed by various blocky talus cones located between 0.4 and 2.2 km downstream of the terminal moraine, resulting in temporary waterlogging in the upper catchment (Figs. 8, 11). These waterlogged sections were identified in the field by fine-grained sediments and in some places by intact Kobresia turf. As water run-off was not detectable in some talus cones, subsurface flow through the cones can be assumed. In lower valley sections with slope gradients between 5° and 7° and broad valley bottoms, the flood dispersed over a width of up to 90 m, four times more than the average width of the channel. Intense fluvial incision and lateral erosion of unconsolidated material were detectable in the field, and the corresponding large sediments loads were accumulated in lower valley sections. These are also traceable on satellite images (Fig. 12). Due to the cumulative effect of temporary blockage in the upper valley section and lateral widening in some portions, the flood slowed down so that shepherds were able to warn villagers of the approaching hazard. Villagers reported the black colour and unusual smell of water, In the lower valley section above the village (Fig. 13), the channel narrows to less than 15 m. This constriction was aggravated by a concrete bridge with massive pillars, built in the 2000s, which was destroyed by the flood. A small footbridge and irrigation canal intakes were also completely damaged. Two sharp and narrow curves increased lateral erosion of the more than 5-m-thick fluvial terrace and inundation of further downstream, where 4000 m 2 of cultivated area, mostly located on the orographic right side, was eroded.
Further downstream, the water dispersed to the wider floodplain, where it destroyed two watermills and two houses and flooded two more houses in front of the Leh-Manali highway, which was only partly damaged. According to field observations, the flood reached a depth of more than 2 m in the village. After the flood, concrete walls (up to 2 m in height) were constructed along the undercut bank as flood protection measures. Some of them were partially damaged and dysfunctional in August 2019. Downstream, new gabions on both sides of the channel, two footbridges and one motorable concrete bridge were built, narrowing the floodplain to a dangerous extent. Interviews showed that villagers are aware of future GLOF risk and potential economic loss due to destroyed fields and damaged infrastructure.

Discussion
Due to the lack of standards for glacial lake mapping, the minimum size of assessed glacial lakes varies widely (   lakes larger than 9000 m 2 (with an estimated water volume of about 43,000 m 3 ) were considered. As many glacial lakes in Ladakh are ice-covered even in summer-evident in the time-series of the Gya lake-their detection by semi-automatic approaches is not possible; and time-consuming controls by visual interpretation are necessary. The integration of all lakes was not practicable considering the size of the mapped area and the spatial resolution of Landsat 5 imagery. Thus, the total number of glacial lakes in Ladakh is much higher than the 192 identified and used for this inventory.
The change detection analyses show that the number and size of glacial lakes have increased since 1969, which is similar to the observed trend in the entire Himalayan region (Bolch et al. 2019). The relative increase in the total glacial lake area, which amounts to 35% in the Trans-Himalayan region of Ladakh between 1993 and 2018, is lower than the observed trend for the eastern Himalaya and Spiti Lahaul (Gardelle et al. 2011;Nie et al. 2017;Sharma et al. 2018). Corresponding to the small size of the glaciers in the Trans-Himalayan region of Ladakh (Schmidt & Nüsser 2017), the glacial lakes are comparably small and only 6% are larger than 100,000 m 2 . Lake changes in Ladakh show a heterogeneous pattern varying from increasing to decreasing lake sizes and recurrent fluctuations. This is similar to observations made by Ashraf et al. (2017) for the adjacent Karakoram, by Shrestha et al. (2017) for the Koshi basin and by Nie et al. (2018) along the Himalayan arc.
The case of Gya underlines the potential hazard of these overlooked lakes. For instance, either due to the ice cover or due to its size being under the threshold of 0.01 km 2 and thereby assumed to not present relevant hazard potentials to downstream villages, the Gya lake was not classified as potentially hazardous by Worni et al. (2013). The omission of potentially hazardous lakes can be problematic, as exemplified with the Chorabari lake outburst event in 2013 (Bhambri et al. 2016). Various studies show that small glacial lakes might cause serious damages and can adversely affect local livelihoods, as witnessed in the case of the Domkhar flood in Ladakh (Narama et al. 2012;Ikeda et al. 2016) and several small floods in Hunza (Ashraf et al. 2012;Parveen et al. 2015). The examples show that even minor events may require further investigation as "emerging threats" (Hewitt and Liu 2010). This holds true for large areas of the Trans-Himalaya, where small proglacial lakes are located in many tributaries. Thus, the threshold for identifying potentially hazardous lakes needs to be adjusted according to local expansion of human habitat, infrastructure and agricultural areas (Wang et al. 2012). There is an urgent need to generate glacial lake inventories that include small lakes by using high spatial resolution satellite imagery (Bhambri et al. 2016). Moreover, the hydrological connectivity of several lakes may trigger the outburst of other lakes due to cascading effects (Erokhin et al. 2018;Kirschbaum et al. 2019;Petrakov et al. 2020). A sophisticated monitoring concept should also integrate satellite data with high temporal resolution to indicate the development of short-lived lakes on glaciers or on debris landforms with buried ice (Narama et al. 2012(Narama et al. , 2018 or fast glacial lake growth (Petrakov et al. 2020) as also evident in the case of Gya.
Beside an improved lake inventory, more bathymetrical surveys are required for a more precise estimation of lake volumes, as even lakes within similar geographical areas do not necessarily have equally predictable volumes (Cook and Quincey 2015). Nevertheless, various empirical formulae have been developed based on the relation between volume and area size using small numbers of measured lakes (Sakai et al. 2012 used 17 measurements;Veh et al. 2020 used 24 measurements;Wang et al. 2012 used 20 measurements;Yao et al. 2012 used 15 measurements). Due to limited bathymetrical surveys and differing lake types, the formulae differ significantly and volumes are over-or underestimated, with uncertainties of around ± 50% (Muñoz et al. 2020), ranging up to 400% (Cook and Quincey 2015). The comparison of Himalayan lake volumes listed in Veh et al. (2020; excluding the Chorabari lake) and calculated volumes based on the empirical formula by Huggel et al. (2002) indicates that in 8 cases the error is less than ± 25%, in 11 cases it ranges between ± 25-50%, and in 4 cases it is even larger than ± 50%. The estimated water discharge derived from lake level differences and empirical formulae differ by about 25 to 66% for the case of Gya. This suggests that the uncertainty in lake volume is primarily due to the absence of bathymetry data rather than errors in measuring lake area and formulae may need to be adapted to local conditions based on bathymetric surveys.
The change detection analyses show that in total 18 lakes released water of more than 200,000 m 3 between 1969 and 2018. It cannot be stated with confidence that each of these water releases caused a GLOF, as some lakes might have drained over longer periods as 1 3 shown for the Gya lake in August 2019. However, at least two cases match with GLOF events in Ladakh reported by Morup (2014). Many GLOFs might be unobserved or unreported due to the sparsely populated area, as in other Himalayan regions (Nie et al. 2018;Veh et al. 2018). On the other hand, some reported GLOFs (Morup 2014;Tabassum and Kanth 2013) could not be verified as no glacial lakes could be identified on satellite imagery within the relevant period. Therefore, in some cases it might be assumed that the reported floods have been caused by cloudbursts. This reflects the general problem of unknown hazards and understated risk that can result from inventories that rely on only a single method of classification or satellite imagery of lower spatial resolution.
As demonstrated in the case of the Gya GLOF, a false sense of security can result in maladaptive behaviour. This means that in spite of the additional effort required, establishing the robustness of such inventories through multiple independent methods may be necessary. Multi-criteria guidelines (Worni et al. 2013) may help to systematically categorize the individual lake hazard. However, the assessment of different variables by using remote sensing data can be difficult, as the freeboard height between lake level and moraine crest needs a DEM with a high spatial resolution; or even on satellite images with high spatial resolution the identification of buried ice in moraines as well as water discharge by tunneling processes cannot be detected. More field observations and measurements are therefore required. Tunneling processes can be assumed as a trigger of sudden water releases in various events (Erokhin et al. 2018;Ikeda et al. 2016;Narama et al. 2018). Furthermore, glacial lakes with an underground drainage are more susceptible to outburst floods than lakes with a stable surface drainage, because intra-moraine channels can be blocked by thermokarst processes (Narama et al. 2018), so that the blockage can lead to overflow and subsequent lake outburst through erosion and breach formation (Erokhin et al. 2018;Petrakov et al. 2012Petrakov et al. , 2020. The ice tunnel can be closed and the lake can be refilled again; thus, several floods can occur from the same glacial lake (Narama et al. 2018) or indicated by the lake level changes in Gya. The time series of Gya indicates a rapid increase in the lake level before the GLOF event in 2014 as well as in 2019. Whether these lake level rises are caused by a blockage of the tunnel or by an increase in melt water due to a massive heatwave as reported by the villagers of Gya in 2014 or high precipitation as observed in 2019 remains unclear. Both abrupt rise of temperature and torrential rainfall may act as triggering factors for GLOFs (Din et al. 2014;Erokhin et al. 2018;Petrakov et al. 2020). The widening of the tunnel system within the moraines and melting of buried ice may increase the future risk of GLOFs (Richardson and Reynolds 2000;Sakai 2012;Westoby et al. 2014;Worni et al. 2013).
As tunneling processes do not necessarily destroy the shape of moraines and drainage lines cross blocky material of scree and talus cones allowing sub-surface flow, the identification of former GLOF events based primarily on sediment tails (Veh et al. 2018(Veh et al. , 2019 or V-shaped trench (Komori et al. 2012;Nie et al. 2018) may be misleading. This means that the investigation of GLOF events caused by tunneling is important for further hazard assessments because it raises the question whether moraine dam failure, including overtopping, is the most common or just the most detectable form of GLOF. Due to the large number of unreported GLOF events and uncertainty about outburst mechanisms, this question deserves particular attention.
Besides the growing risk of GLOFs due to climate change, the vulnerability of local communities increases. Due to population growth and migration processes, habitations expand in most tributaries with dispersed settlement structures including new roads and bridges. As in other regions of the Himalaya (Hewitt and Mehta 2012), this results in potentially high exposure to cryosphere hazards along the river channels, even in areas without documented GLOFs. The integration of cryosphere change with its broader social context, as considered within the emerging field of socio-hydrology in climate-sensitive mountain regions (Nüsser 2017;Nüsser et al. 2019), offers a useful framework for such an analysis of the dynamic coupling of cryosphere hazards and societal response (Mukherji et al. 2019).

Conclusion
In the light of recurring GLOFs in Ladakh, the exposure of local populations to cryosphere hazards needs to be re-evaluated, with consequences for identification of vulnerability and the necessity for adaptation measures. This article offers a method for identifying exposure and the occurrence of GLOFs using analyses of satellite imagery. This is especially useful in cases where oral history is unreliable or where outbreak events occur far from habitation and such floods are simply classified as flash floods of uncertain origins. The reliability of this method could further be established using field research and ground truthing. The consilience and robustness of the methods described here can be thus established using two unrelated forms of detection.