Next Article in Journal
Climatological Drought Monitoring in Switzerland Using EUMETSAT SAF Satellite Data
Previous Article in Journal
Adaptive Modulation and Coding for Underwater Acoustic Communications Based on Data-Driven Learning Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Glacial Lake Outburst Flood Hazard and Risk Assessment of Gangabal Lake in the Upper Jhelum Basin of Kashmir Himalaya Using Geospatial Technology and Hydrodynamic Modeling

1
Department of Geography and Disaster Management, School of Earth and Environmental Sciences, University of Kashmir, Srinagar 190006, India
2
Department of Water Resources Development & Management, Indian Institute of Technology Roorkee, Roorkee 247667, India
3
National Institute of Hydrology, Roorkee 247667, India
4
Department of Ecology, Environment and Remote Sensing, Government of Jammu and Kashmir, Srinagar 190018, India
5
Geological Survey of India, Jammu and Kashmir, Srinagar 190008, India
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(23), 5957; https://0-doi-org.brum.beds.ac.uk/10.3390/rs14235957
Submission received: 5 October 2022 / Revised: 10 November 2022 / Accepted: 15 November 2022 / Published: 24 November 2022

Abstract

:
Climate warming-induced glacier recession has resulted in the development and rapid expansion of glacial lakes in the Himalayan region. The increased melting has enhanced the susceptibility for Glacial Lake Outburst Floods (GLOFs) in the region. The catastrophic failure of potentially dangerous glacial lakes could be detrimental to human life and infrastructure in the adjacent low-lying areas. This study attempts to assess the GLOF hazard of Gangabal lake, located in the Upper Jhelum basin of Kashmir Himalaya, using the combined approaches of remote sensing, GIS, and dam break modeling. The parameters, such as area change, ice thickness, mass balance, and surface velocity of the Harmukh glacier, which feeds Gangabal lake, were also assessed using multitemporal satellite data, GlabTop-2, and the Cosi–Corr model. In the worst-case scenario, 100% volume (73 × 106 m3) of water was considered to be released from the lake with a breach formation time (bf) of 40 min, breach width (bw) of 60 m, and producing peak discharge of 16,601.03 m3/s. Our results reveal that the lake area has increased from 1.42 km2 in 1972 to 1.46 km2 in 1981, 1.58 km2 in 1992, 1.61 km2 in 2001, 1.64 km2 in 2010, and 1.66 km2 in 2020. The lake area experienced 17 ± 2% growth from 1972 to 2020 at an annual rate of 0.005 km2. The feeding glacier (Harmukh) contrarily indicated a significant area loss of 0.7 ± 0.03 km2 from 1990 (3.36 km2) to 2020 (2.9 km2). The glacier has a maximum, minimum, and average depth of 85, 7.3, and 23.46 m, respectively. In contrast, the average velocity was estimated to be 3.2 m/yr with a maximum of 7 m/yr. The results obtained from DEM differencing show an average ice thickness loss of 11.04 ± 4.8 m for Harmukh glacier at the rate of 0.92 ± 0.40 m/yr between 2000 and 2012. Assessment of GLOF propagation in the worst-case scenario (scenario-1) revealed that the maximum flood depth varies between 3.87 and 68 m, the maximum flow velocity between 4 and 75 m/s, and the maximum water surface elevation varies between 1548 and 3536 m. The resultant flood wave in the worst-case scenario will reach the nearest location (Naranaag temple) within 90 min after breach initiation with a maximum discharge of 12,896.52 m3 s−1 and maximum flood depth and velocity of 10.54 m and 10.05 m/s, respectively. After evaluation of GLOF impacts on surrounding areas, the area under each inundated landuse class was estimated through the LULC map generated for both scenarios 1 and 2. In scenario 1, the total potentially inundated area was estimated as 5.3 km2, which is somewhat larger than 3.46 km2 in scenario 2. We suggest a location-specific comprehensive investigation of Gangbal lake and Harmukh glacier by applying the advanced hazard and risk assessment models/methods for better predicting a probable future GLOF event.

1. Introduction

Glaciers are highly sensitive to climatic fluctuations [1,2,3,4,5] and are therefore considered the key indicators of climate change [3,6,7,8]. The Himalayan region is one of the few regions where climate change effects are most prominent [9,10]. Most studies have reported that Himalayan glaciers are in a state of continuous recession [11,12,13,14,15], except the glaciers in the Karakoram region, which are stable or advancing [12]. Retreating glaciers at a faster rate are believed to have triggered many catastrophic events in recent decades; for example, GLOFs are becoming a serious and emerging threat in the high mountainous regions of the world in general and the Himalayan region in particular [4,13,14,16]. GLOF is a type of cryospheric natural hazard that can cause catastrophic fatalities and damage to neighboring areas, and the Himalayan region is one of the well-known GLOF-vulnerable regions in the world [17].
The Himalayan region, being geologically young and fragile, has witnessed significant warming in the present century compared to the past century. Warming has been reported throughout the Himalayan region at rates ranging from 0.15 °C to 0.60 °C/10 yr, which is very high when compared to the global mean warming rate of 0.74 °C/100 yr [18,19,20]. As a result of this warming trend, Himalayan glaciers are losing mass at a faster pace [21,22,23], resulting in the evolution and expansion of various glacial lake types [24,25,26]. Among these, potentially dangerous glacial lakes pose a major threat to the communities and infrastructure in the neighboring regions.
The study of Himalayan glaciers has substantially advanced in the last decade, and the spatiotemporal frequency and degree of changes in glacial lakes have been reported as well [27]. According to climate model projections based on high-emission scenarios, it has been reported that approximately 65% of the glaciers’ ice mass in High Mountain Asia (HMA) will vanish by the end of the 21st century [28]. As a result, the size and number of proglacial lakes (glacier-fed) and supraglacial may increase, which coalesce to form large glacier lakes. Further, it is essential to mention that most of these lakes are formed behind unconsolidated (loose) and inherently unstable moraine dams.
Moraine dam failures induced by the high rate of lake area expansion, seepage, dead-ice melting, degradation of permafrost, seismic activity, cloudburst, a flood wave generated by a rock avalanche/fall or due to the snow/ice avalanche into the lake, and an abrupt rise in the water level may cause these lakes to burst [29,30,31]. Such outburst events often have very high destructive potential, as an enormous amount of water is discharged in a short span of time with a capacity to erode the unconsolidated/loose material, potentially resulting in a strong and violent flow that can travel a long distance in the downstream regions.
Earlier, glacial lake studies in the former state of Jammu and Kashmir were carried out at different levels, for example, at the international level (glacial lake inventory of HMA) by Ives et al. [32], Wang et al. [33], and Chen et al. [34]; at the national level by the Central Water Commission [35], National Remote Sensing Center Hyderabad [36], Ives et al. [32], Fujita et al. [37], and Worni et al. [38]; and at the basin level by Dar et al. [39] and Ahmed et al. [40,41]. The Potentially Dangerous Glacial Lakes (PDGLs) in the area have been identified by seven previous studies. According to Kumar and Goel [42], Gangabal lake and Konsar-Nag lake are susceptible lakes. In recent investigations on glacier lakes in the Indian Himalayas, Mal et al. [43] categorized 25 glacial lakes as PDGLs. Among these, four lakes are located in Jammu and Kashmir; however, only one of them, namely Sona-Sar lake, is found in the area under study. Sona-Sar lake was also confirmed as a high-risk Lake (HRL) by Dubey and Goel [44] and a high-hazard glacial lake (HHGL) by Ahmed et al. [41]. Dubey and Goel [44] have reported four glacial lakes (Hirubagwan, Krishna-Sar, Gad-Sar, and Har-Nag) as very high-risk lakes (VHRL) in the upper Jhelum basin. Dimri et al. [45] have identified 40 lakes as potentially dangerous and two as critical lakes in Jammu and Kashmir. GLOFs pose a severe threat to populations living in mountainous regions. To demarcate vulnerable areas for GLOF hazard management and mitigation, it is necessary to identify potentially dangerous lakes and simulate outburst floods.
Thus, keeping in view the hazard level of the dangerous lakes assessed by various regional studies, in-depth research evaluating the GLOF hazard and risk assessment of the individual lakes is quite rare. In the present study, we particularly focus on Gangabal lake (potentially dangerous), the largest lake in volume and area, located in the Sind catchment of the UJB of the northwestern Himalayas. The specific objectives of the current study include: (1) analyze the spatiotemporal dynamics of the lake (Gangabal) and its feeding glacier (Harmukh) by utilizing multitemporal satellite data between 1972 and 2020; (2) study the surface velocity, geodetic mass changes, and ice thickness of the Harmukh glacier; (3) use coupled HEC-RAS 1D and 2D unsteady flow analysis to model a worst-case GLOF scenario and locate exposed parts in the closest settlement, which is situated 9.2 km from the outlet of the lake; and (4) analyze the impact of the probable GLOF event on downstream landuse and landcover.

2. Review of Literature

GLOFs have become a severe threat in recent decades as a result of increased human settlements, anthropogenic activities, and unplanned infrastructural development in the fragile high alpine environment, which were not previously inhabited [40,46]. The Himalayan region is home to around 200 potentially hazardous glacial lakes and has witnessed 40 GLOF events in the past forty years [32,47]. The major GLOF events which are well documented in the literature include: Chog Kumdan event in 1929 [48], Gelhaipuco 1964 [49,50], Nare 1977 [51,52], Nagma Pokhari 1980 [53], Cirenamco 1981 [29], Ghulkin GLOF 1981 [54], Zhangzangbo Boqu 1981 [55,56,57,58], Dig Tsho 1985 [55], Luggey Tsho 1994 [59], Tampokhari 1998 [51], and Sabai Tsho 1998 [60]. The Chorabari GLOF event in 2013 [4] caused heavy loss to property and life in the downstream region.
GLOFs are frequently associated with glacier mass loss and are thought to be one of the most essential geomorphological and hazardous consequences of melting glaciers [61]. Considering the ongoing glacier mass-loss induced by climate change, the study of past and probable future GLOFs is a very pertinent topic [62,63]. GLOFs and their related processes have major ramifications in a variety of domains, for example, risk management, disaster risk reduction, designing hydropower plants, and geomorphology and sediment yield [64,65]. GLOFs in the Himalayas have previously resulted in catastrophic deaths and destruction in neighboring regions [66]. With the continuous glacier mass loss in the region [67], GLOF risk has been exacerbated and needs much more attention owing to the potentially catastrophic damage [17]. It is necessary to construct a database of GLOF events from the existing literature in order to expose the process and mechanism of GLOFs for mitigation, hazard assessment, and risk management of glacial lakes. We compiled an inventory of historical GLOFs in the Himalayan region (see Figure 1). A total of 158 GLOF events were documented in the Indian Himalayan region, and the analysis revealed that most of the GLOF events had been seen in the Ladakh region of the northwestern Himalayas (Figure 1). However, few of them are repeated events.
Many glacial lakes surrounded by unconsolidated moraines pose a severe threat to the people living in the world’s mountainous regions. The dams could fail unexpectedly, causing devastating floods with peak discharges far beyond average. Glacier-dammed and moraine-dammed lakes are the two main sources of GLOFs in the Himalayan region [68]. The lake outbursts occur when glacier ice initiates to float on the dammed lake or when the stored water erodes the dam laterally, subglacially, or superficially [69]. Moraine-dammed lakes are susceptible to outbursts due to a variety of conditioning factors and triggering mechanisms.
The majority of these glacial lakes are inaccessible due to their remote locations, tough terrain, and harsh weather conditions. Therefore, continuous monitoring and assessing these glacial lakes at regular intervals through field-based investigations and observations are impossible. Satellite remote sensing and GIS serve the purpose of monitoring and assessing these glacial lakes at different spatial and temporal scales. Furthermore, to prevent and monitor the GLOF hazard and to evaluate the damage and destruction that may occur in the future, the integration of various geospatial techniques and modeling approaches could play a key role in identifying hazardous/dangerous lakes and monitoring GLOF events in real-time. Such robust techniques have been effectively used worldwide to monitor and assess glaciers, glacial lakes, and associated hazards. The GIS and RS-based tools and techniques reasonably allow rapid investigation of huge glaciated terrains and are therefore considered suitable for glacial lake and GLOF studies [40].
Geospatial techniques coupled with hydraulic and hydrologic simulations can be used to obtain information about the status of glaciers and glacier lakes at high spatial and temporal resolutions and map flood extent [70]. Satellite images have already been utilized in several global studies to monitor glaciers and glacier lakes, as well as hydraulic and hydrologic models such as HEC-RAS, HEC-HMS, HEC-GeoRAS, HEC-Geo-HMS, Arc Hydro, etc., for GLOF mapping. The HEC-RAS model helps us analyze the possible GLOF impacts, especially in generating flood hydrographs and several other flood wave parameters such as flood water depth, water surface elevation, and flow velocity. The key objective of this study was to use RS, GIS, and hydrodynamic modeling (HEC-RAS 5.0.7) to analyze the hazard and risk associated with the potentially dangerous glacial lake (Gangabal) located in the Sind catchment of the Jhelum basin. GLOF modeling has the ability to provide crucial insights for predicting future events, which can aid in damage assessment, mitigation, and planning.

3. Study Area

This study focused on the largest glacial lake, Gangabal, located in the Sindh catchment of the Jhelum basin in the Western Himalayas, with an area of 1.65 km2, per recent estimates. The lake is positioned between 34°25′55.49″N longitude and 74°55′27.12″E lattitude at an altitude of 3576 m asl [71]. The lake is directly fed by the meltwater of the Harmukh glacier, one of the largest glaciers in the Jhelum basin [14]. Water from the lake is directly supplied to another lake in the downstream region, Nundkol, located at a distance of 0.9 km from the outlet. It then merges with the Sind tributary through Wangath Nalah (Figure 2). The distance between the lake and the feeding glacier is 0.78 km, with a maximum slope of 32 degrees and an elevation gain of 393 m from the upper point of the lake to the glacier front. The lake has a maximum length of 2.5 km and maximum width of 1 km. It is located in the Ganderbal district near Ganga Gali. The lake has an average depth of 42.5 m, obtained through the improved area volume scale technique. The water level of the Gangabal lakes is rising at a negligible rate of 0.009 m/yr, as observed by Srivastava et al. [72]. Gangabal lake has glacial lake id no 0143J1501489 as per the Glacial Lake Atlas of Indus River Basin prepared by NRSC in 2020 [73]. Long-term climate data (1981–2020) of the meteorological observatory nearest to Gangabal lake, i.e., Shalimar, located in the Srinagar district, shows that Tmax has increased by 1.79 °C from 19.0 °C (1981) to 20.79 °C (2020), with a slope of 0.04 and R2 value of 0.19 (Figure 3). Tmin, on the other hand, does not depict any remarkable impact, decreasing by 0.40 °C in the last 40 years. The mean temperature has shown a significant increased mean of 0.7 °C during the previous 40 years, with an R2 value of 0.13. A decreasing trend was also observed for the annual precipitation, as depicted in Figure 3.
Gangabal lake has been classified as a potentially dangerous lake by various researchers in the past, for example, a high-risk lake by Gupta et al. [36], a potentially dangerous lake by Guru et al. [74], medium-risk lake by Ahmed et al. [41], and important GLOF lakes by Kumar and Goel [42] and Kumar et al. [15]. The GLOF susceptibility of Gangabal lake is provided in Section 5.6. Considering the GLOF potential of the lake, the present study prioritized it for hazard and risk assessment using the integrated approaches of geospatial technology and hydrodynamic modeling.

4. Materials and Methods

4.1. Data

Area change analysis was performed for Gangabal lake and its feeding glacier using multitemporal satellite data. In this study, we used high-resolution Resourcesat-2 LISS-IV data scenes (5.8 m), Sentinel-2A (12 m), Landsat TM (30 m), Landsat ETM (30 m), Landsat ETM+ (30 m), and Landsat OLI (30 m) to map the glacial lake outlines from 1972 to 2020 and glacier outlines from 1990 to 2020. All the procured scenes of Resourcesat-2 were downloaded from the NRSC Hyderabad web portal https://www.nrsc.gov.in/, accessed on 25 October 2021. Sentinel-2A data scenes were acquired from the “Copernicus internet hub of European Space Agency” (ESA), and Landsat imageries were downloaded from the USGS Earth Explorer data portal https://earthexplorer.usgs.gov/, accessed on 12 July 2021. Cloud cover seems to be the major hindrance during the selection of scenes; therefore, the imagery from the adjacent month or year was preferred.
Advanced Land Observing Satellite (ALOS) Phased Array Synthetic-Aperture Radar (PALSAR) radiometrically terrain corrected (RTC) Digital Elevation Model (DEM) was obtained from the web portal, i.e., https://asf.alaska.edu/ (accessed on 12 July 2021) and was utilized to obtain the topographic information (elevation, slope, aspect, hill shade, etc.) for the GLOF modeling. ALOS PALSAR has a spatial resolution of 12.5 m. It has been widely used for hydrodynamic modeling to generate geometry data such as river center lines, bank lines, flow lines, cross sections, and breach invert levels. Google Earth images were also used in this study for validation and cross-checking. SRTM-DEM for 2002 and TAN-DEM X for 2012, with a spatial resolution of 90 m, were used for the mass balance calculation of the glacier. Recent Landsat-8 panchromatic (scene-8) scenes were utilized to estimate glacier velocity. RGI glacier extents that possess the required parameters for the depth estimation were used as input in the GlabTop-2 model to derive the depth values of the glacier. The details of satellite data and other ancillary data sources used in the study, including their source, number of bands, spatial and temporal resolution, and purpose, are given in Table 1. The freely available CubeSat data downloaded from the Planet web portal https://www.planet.com/ accessed on 2 August 2021 were used to calculate various glacier and glacial lake parameters.
A step-by-step methodology was employed to understand and monitor the glacier and glacial lake-associated hazards in the study region [44,66,75,76,77,78]. The methodology for assessing the status of the Harmukh glacier and the hazard assessment of Gangabal lake and analyzing the glaciological characteristics of its feeding glacier (Harmukh) are summarized in the following section. The overall methodology of this work is illustrated in Figure 4.

4.2. Methods

4.2.1. Characteristics of the Adjacent Glacier (Harmukh Glacier)

The glacier mapping and change detection analysis was carried out using multi-temporal satellite data. In the past studies, glacier mapping has been done using different methods and techniques. In this study, we have used NDGI technique combined with visual interpretation to identify and delineate glaciers as suggested by [3,79], which is considered the most accurate method for glacier mapping [80].

4.2.2. Glacier Depth

Glacier ice thickness estimation is generally obtained through mathematical equations, for example, volume–area scaling [81,82,83,84]. In this study, the depth of the Harmukh glacier was calculated using the GlabTop-2 model. GlabTop-2 is a slope-dependent, gridded ice thickness estimation model that is based on the flow mechanics of the average basal shear stress (τ), shape factor (f), and surface slope [85,86]. The formula is as follows (Haeberli and Hoelzle [87]):
I c e   t h i c k n e s s   h = τ f ρ g   sin α
where in Equation (1) h represents the ice thickness, τ indicates the basal shear stress, f denotes the shape factor, and ρ is the ice density.
Glacier polygons are transformed into a raster matching the DEM cells. Cells are differentiated as inner glacier cells (light yellow), glacier marginalized cells (orange), cells adjacent to the glacier (dark cyan), and nonglacier cells (white). Dark green cells depict randomly selected cells (r) for which local ice thickness is measured; the black square represents the buffer of variable size, which is enlarged (black-dashed square), until an elevation extent of hmin is attained inside the buffer. The value for parameter f is considered 0.8, which generally works well for Himalayan glaciers [86,88]. The hmin is considered as 20, and fractions of cells are taken as 30% for random cell pick-up, g is the acceleration due to gravity (9.8 m s−2), α is the slope, and ρ is the ice density, which is 900 kg m−3 [89]. A similar projection system is necessary for all the datasets used as an input file in the Python-based GlabTop-2 model to run the code smoothly. The working principle of the GalbTop-2 model is illustrated in Figure 5.

4.2.3. Glacier Surface Velocity

Glaciers move downslope under the influence of gravity and their underlying mass, which is highly dependent upon the topographical and geomorphological characteristics [90]. Therefore, glacier velocity (movement) plays a significant role in mass distribution received at the accumulation and ablation zone of the respective glaciers. The surface displacement/velocity of the Harmukh glacier was derived through the “Co-registration of the Optically Sensed Images and Correlation” (Cosi–Corr) model, a plug-in for the Environment for Visualizing Images (ENVI) interface. Cosi–Corr is freely available software from the California Institute of Technology and provides an important tool to accurately coregister and orthorectify and correlate optical remotely sensed scenes/images to obtain deformations on the ground surface from multitemporal satellite images. The Cosi–Corr model was initially developed and released by the tectonics observatory at Caltech in 2008 [91] to extract the crustal deformations caused by earthquake movements.
Local displacements among the satellite images (pre- and post-event) with different spatial resolutions can be obtained through the Cosi–Corr software [92] and have been extensively used for velocity estimation in various glaciated regions of the world [93,94,95,96,97]. The method coregisters optical satellite images and correlates them to determine the resulting displacement. The present study used two Landsat 8 scenes (pan bands) having 15 m spatial resolution for two consecutive years, i.e., 2020 and 2021, to obtain the glacier velocity. The surface displacement images were thoroughly checked and corrected for possible errors before the glacier surface velocity estimation. Three bands—east–west displacement, north–south displacement, and signal-to-noise ratio—are provided by the Cosi–Corr generated output. The correlation of two Landsat scenes (pre-event and post-event) was carried out using 32 × 32 initial and 16 × 16 final window size with 2 XY steps. The mask threshold was set at 0.90, and robust iteration 2 was used, as suggested by [98]. The results of the correlation yield a displacement image of east–west and north–south directions together with a signal-to-noise ratio image (SNR). The image with north–south and east–west directions was used to calculate the resultant surface movement with an SNR value greater than 0.9. The vector file showing direction was also generated from the discarded image with an SNR value of less than 0.9. The errors related to the calculated surface velocity of the glacier were assessed over the stable patches surrounding the glacier body by presuming less or no movement over the stable areas of a glacier [98,99]. The mean and standard deviation of displacements were summed across stable areas of a glacier and then divided by the interval of time between the two correlated satellite images. The final accuracy of the results generated from the model generally depends upon the exact coregistration, orthorectification, less cloud cover, and shadow of the satellite scenes.

4.2.4. Glacier Mass Balance

Understanding glacier mass balance is important to decode the response of glaciers to climate change and how glacier changes affect water supplies and sea level rise [99]. The glacier mass balance can be estimated through several techniques, such as Accumulation Area Ratio (AAR), Improved Accumulation Area Ratio (IAAR), Geodetic, energy balance, and glaciological methods. Among all these techniques, a glaciological method for mass balance estimation is the most recommended owing to its accuracy, reliability, and insight into the field conditions [100]. However, the glaciological method was not possible for the present glacier due to its remote location and financial and instrumental constraints. Therefore, in this study we used the geodetic balance assessment technique as it is considered a much simpler procedure and quick method to analyze glacier mass changes through time [101,102,103,104]. The geodetic technique quantifies the ice thickness changes by differencing the DEMs of two-time points. The derived thickness changes can be subsequently converted to mass loss by employing the glacier area and density assumptions [104,105]. The geodetic mass changes of Harmukh glacier were estimated using SRTM-DEM for the years 2000 and 2015. The inaccurate or no coregistration of the DEMs may result in enormous mass change values that lead to uncertainties in the final results. Therefore, coregistration of the DEMs was performed in the ERDAS imagine software before the DEM differencing to lessen the uncertainties. Furthermore, the interquartile range method recommended by Huber et al. [106] was used to find and eliminate the outliers in the off-glacier and on-glacier elevations. The widely applied deterministic interpolation technique, i.e., the Inverse Distance Weighted (IDW) method, was utilized for the voided areas derived from the DEM differencing results.
The errors in the glacier mass changes resulted due to the geodetic mass balance method (∇M) associated with the glacier area (A) and the density of ice (ρ), which is assumed to be 850 ± 60 kg m−3, together with the surface elevation change h, as stated in [102]. Uncertainty associated with the elevation changes ( σ Δ H ) was estimated by Equation (3).
Δ M M = Δ d h d h 2 + Δ ρ i c e ρ i c e 2 + Δ A A 2
σ Δ H = σ d e m 2 + σ V o i d   f i l l 2 + σ T D X d a t e 2  
where σ d e m 2 in the equation represents the uncertainty in the DEM and σ V o i d   f i l l 2 denotes the uncertainty that resulted during the void-filling process. The uncertainty related to the date is represented in the equation by σ T D X d a t e 2 , which was calculated as ±two times the elevation change per year.

4.2.5. Glacial Lake Mapping

Several automated and semiautomated methods and techniques are available for extracting and mapping glacial lakes using geospatial technology, e.g., supervised and unsupervised classification, band rationing technique, NDWI, etc., as used by previous studies in different parts of the Himalayas [13,107]. NDWI, combined with visual interpretation, seems to be the most appropriate technique for identifying and mapping glacial lakes. We used the NDWI technique coupled with manual delineation, as suggested by Huggel et al. [108]. Shadow, cloud cover, and water-rich soils are the main challenges faced during the creation of NDWI maps. Meanwhile, DEM-derived hill shade, aspect, and slope values were removed from the final NDWI map to overcome such problems. High-resolution Google Earth imageries were used for smooth visual interpretation. The uncertainty in the present study was computed for the glacier and lake according to their areal extents following similar methods applied by [109,110]. The resolution of the satellite imagery determines the accuracy of the glacier and lake delineation. High-resolution images cause fewer uncertainties in comparison to low-resolution images.

4.2.6. Glacial Lake Depth, Volume, and Peak Discharge Estimation

The surface area of a glacier lake becomes a key indicator for measuring lake volume and hazard because the depth of a glacial lake is difficult to determine using remote sensing-based methods. In such cases, empirical equations developed in the various glacial lake-bearing river catchments are frequently used. Estimating parameters such as lake depth, volume, and potential peak discharge is essential before assessing and simulating a GLOF event. There are various equations available in the literature to estimate such parameters. Here, we used the widely employed empirical formula given by Huggel et al. [108] to calculate the lake depth and were subsequently compared with other formulae given by Fujita et al. [37] and Patel et al. [111]. In contrast, volume was calculated by several formulas generated for different regions of the Himalayas, such as for the Koshi basin of Nepal Himalayas by Liu et al. [112], the Sikkim Himalayas by Sharma et al. 2018, the Himalayan region by Fujita et al. [37], and the Chandra basin in the western Himalayas by Patel et al. [111]. Additionally, empirical equations developed for regions outside the Himalayas were also included, such as for the Swiss Alps by Huggel et al. [108] and the Cordillera of western Canada by Evans [113]. The results obtained by employing all these equations showed considerable variation. Thus, an average of all these equations was used in the present study, as suggested by Emmer and Vilmek [114] and Chowdhury et al. [115]. Similarly, average peak discharge was calculated using three empirical equations given by Evans [113], Popov [116], and Huggel et al. [108]. The details of the equations used to calculate the depth, volume, and peak discharge are given in Table 2.

4.2.7. GLOF Hazard Assessment Using HEC-RAS 5.0.7

Hydraulic modeling is an essential part of developing a reliable flood forecasting system. The results of hydraulic model simulations may be used to create inundation maps that can be utilized by community officials or the general public to assess their flood risk. GLOF hazard and risk modeling has been carried out across the world using a variety of methods. For instance, few studies have employed geometric models such as modified single flow [118,119], Monte Carlo Least Cost Path [120,121], and Random Walk Process [122] to analyze and estimate probable GLOF flood inundation. Such types of models need minimal input data and may be applied to a wide range of lakes, but they can only provide a first-order evaluation and cannot provide accurate flood maps [123]. However, more complex modeling software such as HEC-RAS 1D, HEC-RAS-2D, FLO-2D, RAMS, SOBEK, DAM BREAK, and BASEMENT has also been used in different parts of the world [38,50,101,124,125,126,127,128]. These models can be used to produce flood hydrographs and other important flood wave data, such as water depth, flow velocity, and water surface height, which can be used to develop risk management and mitigation plans.
The HEC-RAS software is an integrated and open-source software developed in a multitasking environment [54]. It has been widely used to study glacial lake hazards [129]. It allows for performing complex one-dimensional (steady flow) and two-dimensional (unsteady flow) calculations, water temperature, water quality modeling, and sediment transport computations. Graphical User Interface (GUI), data management and storage system, separate analysis components, graphics, mapping, and reporting tools are essential tools embedded within the system. It is user-friendly software designed to perform hydraulic calculations for natural and artificial channels, floodplain areas, over banks, and levee-protected areas. The model is also utilized to simulate flood scenarios resulting from glacial lake breach outburst and is based on one-dimensional St. Venant equations. In the present study, we used the 2D hydrodynamic model, i.e., Hydrologic Engineering Center’s River Analysis System (HEC-RAS) version 5.0.7, to simulate the GLOF event in the Sind catchment of the Jhelum basin. In this study, two hydrodynamic modeling GLOF scenarios, i.e., 100% and 75%, were performed for Gangabal lake based on the volume calculated through the empirical equations (see Section 4.1) with varied breach formation time and breach width. In the first stage, 1D HEC-RAS was used to create a breach hydrograph at the lake site. The essential breach parameters that are used as input to calculate the breach hydrograph include breach formation time, breach width, breach height, and volume of the lake. Breach width and breach formation time were derived through Froehlich’s empirical equations:
Bw = 0.1803Ko (Vw)0.32(hb)0.19
Tf = 0.00254(Vw)0.53(hb)0.9
where Bw and Vw are the breach width (m) and volume of the lake (m3), whereas hb and Tf are the breach height (m) and breach formation time (h). Equation (1), mentioned in the above section, was used to compute the volume of the lake. Bw is the function Vw and hb. These equations have been widely used in past studies to approximate such parameters [104,130]. The volume of the lake was estimated through the improved volume area equation given by Qi et al. 2022. The breach heights (hb) were estimated through the ALOS-DEM [104]. In the present study, the first potential GLOF scenario was considered a worst-case event wherein 100% volume (73 × 106 m3) of water was considered to be released from the lake with breach formation time (bf) of 40 min and breach width (bw) 80 m. In the second scenario, 75% of the lake volume was considered for the potential GLOF routing with a breach formation time (bf) of 30 min and breach width (bw) of 40 m (Table 3). Subsequently, generated breach hydrographs were routed through a river starting at the lake site to the last selected location, namely Drag Tanga village, located at a distance of 29.4 km.
Unsteady flow analysis was performed for Gangabal lake to observe the peak flow propagation, extent of flood inundation, water surface elevation, and flood water depth. Flow and geometry data are the two important datasets required to run the unsteady flow analysis in the HEC-RAS environment. The complete details of the software are given in the web link https://www.hec.usace.army.mil/. We used an implicit finite volume approach to simulate the flood as an unsteady clear water flow. In actuality, sediment may be entertained in the water flow, altering its rheology and reducing the traveled distance. As a result, assuming a clean water flow reduces the validity of our findings, which is a typical issue in GLOF research throughout the world. However, since modeling erosion, deposition, and transport requires comprehensive field-based data, which is not accessible, this simplification must be accepted [131]. The drainage basin and other elevation parameters were derived from the DEM. Geometry data such as centerline (river), bank lines, flow paths, storage area, and cross sections were created in the HEC-GeoRAS interface and were subsequently imported to the HEC-RAS model. The cross-section profile, XYZ perspective plot, and river bed profile downstream from Gangabal lake are depicted in Figure 6.
The unsteady flow analysis requires upstream and downstream boundary conditions to run the model. The inflow hydrograph is used as an upstream boundary condition, whereas normal depth derived from the river bed profile is taken as a downstream boundary condition. The downstream GLOF routing and the assessment of the flood hydraulics via the channel were done by using the outflow hydrograph (model output) as a boundary condition. Subsequently, 2D hydrodynamic modeling was carried out across the river channel. The 2D flood routing mainly requires the DEM as primary input data. The other data sets include cross sections (500 m on both sides of the channel), flow lines, meanings coefficient, normal depth, etc. Chow [132] suggested a Manning coefficient ranging between 0.03 and 0.07 for the hilly areas having steep slopes with less or no vegetation cover, gravels, boulders, cobbles, and bushes present along the river banks. Considering this, we used a Manning coefficient range of 0.04 to 0.06 for the GLOF routing in the current study.
The computational time was set as 24 h so that the flow wave (modeled) has enough time to propagate downstream to the last selected station. The model provides output in the form of flood inundations, peak flood, flood heights, maximum flood depths, and velocity and routed flood hydrographs throughout the flood channel. In the final stage, the potential impacts of the GLOF event were evaluated at various important locations in the downstream region.

4.2.8. Downstream Impacts of GLOFs

Landuse and landcover maps of the downstream area of the lake were generated to analyze the impact of the GLOF event on the various landuse classes, e.g., agriculture, forest, settlements, etc. Sentinel-2A data from the recent year with a fine spatial resolution of 10 m were utilized for this purpose. A maximum likelihood supervised classification algorithm was used to classify the satellite imagery. Supervised classification has been widely used for the classification of satellite imageries across the world [133,134]. Supervised classification is a three-step process that includes training samples, allocation of classes, and testing. ArcGIS 10.3 software was used to carry out this classification, and at the final stage, the results derived from the classification algorithm were tested for accuracy using the kappa coefficient.

5. Results and Discussion

5.1. Glacier and Lake Changes

Multitemporal satellite data for the years 1972, 1981, 1992, 2001, 2010, and 2020 were used to analyze the changes in the area of Gangabal lake. The data analysis shows that Gangabal lake has progressively increased in the area during the last 48 years. The lake area increased from 1.42 km2 in 1972 to 1.46 km2 in 1981, 1.58 km2 in 1992, 1.61 km2 in 2001, 1.64 km2 in 2010, and 1.66 km2 in 2020. Overall, the lake area experienced a 17% growth in its area between 1972 and 2020 at an annual rate of 0.005 km2 per year. Previously, Dar et al. [39] have shown a 0.47 km2 (36%) increase in the area of Gangabal lake between 1980 and 2016. The variation in the results may have been caused by digitization errors and the satellite imagery used. However, in this study, the post-correction (manual corrections) of the vector polygon layer derived from NDWI was performed by more than one interpreter to lessen the uncertainty in the lake area. Furthermore, the lake area for the years 2010 and 2020 was validated with the area of the lake polygon digitized on the high-resolution Google Earth imagery. The lake expansion is depicted in Figure 7 on different Landsat images.
The spatiotemporal analysis of the glacier, namely the Harmukh glacier, associated with Gangabal lake, was also carried out for a period of 30 years during 1990–2020. The data analysis indicated that the area of the glacier in 1990 was 3.6 km2, which decreased to 2.9 in 2020. The glacier area for the years 1992, 2001, 2010, and 2020 was estimated at 3.6, 3.3, 3.1, and 2.9 km2. Deglaciation of 0.7 km2 (19 ± 5%) since 1990 at the rate of 0.023 km2 per year was observed (Figure 8). The glacier also revealed an upward shift in ELA from 35 to 98 m during the last three decades. The glacier snout has also shown a retreating trend, with a 234 m retreat between 1990 and 2020 at 7.8 m per year. In a recent study, an area loss of 17% (0.6 km2) between 1980 and 2018 was also reported by Mir et al. [135].
An area loss of 0.67 km2 between 1969 and 2020 was also reported in past studies by Ahmad et al. [82], with a length change retreat of 5.12 m per year. Murtaza et al. [136] observed total deglaciation of 1.23 km2 (16 ± 0.24%) between 1992 and 2018 in the Harmukh range of the northwestern Himalayas. The Harmukh glacier, the largest glacier in the entire range, lost an area of 0.34 km2 (10 ± 0.88%), with a snout retreat of 187 m observed during the 28-year study period [136]. The other individual glacier studies across the Kashmir Himalayan basin have shown a retreating pattern; for example, Nehnar glacier has shown a retreat of 1.41 km2 (50.35%) during 55 years [84], and the Kolahoi glacier lost an area of 23% since 1962 [137], Bodpathri glacier lost 0.5 km2 (19.5%) during the past 56 years Ahmad et al. [82]. Additionally, glacier area loss, snout retreat, and length changes have been reported in various glaciers of the central and eastern Himalayas.

5.2. Geodetic Mass Balance

Glaciological (direct) and geodetic (indirect) techniques are two important methods to estimate the mass balance of a glacier. In the present study, the geodetic method was employed to determine the mass balance of the Harmukh glacier, as direct field-based measurements are not feasible owing to difficult and unapproachable terrain. The results derived from the DEM differencing show an average ice thickness loss of 11.04 ± 4.8 m for the Harmukh glacier at the rate of 0.92 ± 0.40 m/yr between 2000 and 2012. The ice thickness changes over the glacier surface range from 0.31 to 1.02 m/yr. The spatial distribution of elevation changes (ice thickness) is illustrated in Figure 9. The findings reveal negative mass changes over the ablation parts of the glacier and positive changes over the accumulation zone.
The findings on average ice thickness loss of the Harmukh glacier are in close agreement with a study carried out by Abdullah et al. [138] over the Upper Indus basin, which reports an average ice thickness loss of 1.12 ± 0.40 m/yr for the Jhelum basin glaciers. However, the average ice thickness changes derived from the DEM differencing technique vary greatly across the different river catchments of the Himalayan region. For example, an average ice thickness loss of 0.46 ± 0.26 m/yr in the Ladakh range, 1.17 ± 0.4 m/yr in the Zanskar range, 1.12 ± 0.40 m/yr in the Greater Himalayan range, 1.69 ± 0.60 m/yr were reported by Abdullah et al. [138]. Another study by Majeed et al. [104] showed an average ice thickness loss of 0.38 m/yr for the Pangong group of glaciers in Trans-Himalayan Ladakh. The accuracy of the glacier thickness changes derived from the DEM differencing technique largely depends upon the accuracy (coregistration) of the utilized DEMs and their spatial resolution. We recommend that with the use of high-resolution DEMs, better results can be inferred through the geodetic mass balance estimation technique.

5.3. Glacier Depth

The glacier outlines of “Randolph Glacier Inventory (RGI) version 6.0″ were extracted and used to calculate the ice thickness of the Harmukh glacier. The RGI is an open-source database that consists of important parameters (flowlines, elevation, slope, etc.) that are prerequisites for the GlabTop-2 model to compute the glacier ice thickness along with the DEM. The minimum and maximum ice thickness for the Harmukh glacier were estimated to be 7.3 and 85 m, respectively, whereas the mean thickness was calculated as 23.46 m. The spatial distribution of ice thickness was also analyzed by reclassifying the ice thickness map into four classes at equal intervals to estimate the area in each thickness class. As mentioned in Table 4, the maximum area of 47.18% (1.36 km2) of the Harmukh glacier area comes under ice thickness class 1 (0–30 m), 41.85% area is covered with ice thickness class 31–60 m, whereas only 10.95% area of the glacier is covered with the ice thickness greater than 60 m. The maximum ice thickness was observed in the central parts (central flow lines) and near the ELA of the glacier. In contrast, the lowest thickness (pixel value) was observed at the glacier boundaries.
The model has been widely used and tested for different glaciers worldwide to calculate the ice thickness and bed over-deepening sites [88,97,102,103,139,140,141,142,143,144]. The model calibration and validation performed on the Chota Shigri and Bara Shigri glaciers located in Central Himalaya by Sattar et al. [97] revealed a reasonable agreement between GPR data and estimated ice thickness.

5.4. Glacier Velocity

A remote sensing-based technique was used in the present study to estimate glacier surface velocity/displacement due to the lack of field-based velocity measurement, which is laborious, time-consuming, and requires financial support. Among various tools and techniques (for example, Im-GRAFT, CIA, and IMCOR) that are used for the estimation of the surface velocity of the glacier, the Cosi–Corr model is appropriate and reliable for the estimation of glacier surface displacement (direction and magnitude) by employing optical images [145], and it has various processing tools that provide additional information [91]. The model yielded a mean velocity of 3.2 m/year using the pre-event (2019) and post-event (2020) Landsat-8 pan bands with a spatial resolution of 15 m. The minimum and maximum surface velocity vary between 0 and 7 m/yr. The velocity was observed to be higher in the upper parts of the ablation zones, varying between 4 and 6 m/yr and then gradually decreasing toward the glacier’s snout, as depicted in Figure 10. The same phenomena have been observed in the Dhauliganga basin of central Himalaya by Sattar et al. [97]. The Harmukh glacier is a clean surface glacier; therefore, surface velocity may be greater than other debris-covered glaciers of the Jhelum basin. As illustrated in Figure 10, a significant variation in velocity and flow was observed across the glacier (accumulation zone, ablation zone, terminus). Several studies have carried the glacier surface velocity using satellite images at the catchment level [33,93,94,146] and on individual glaciers across the Himalayan region [147,148,149]. Garg et al. [150] have studied the glacier surface velocity of the Sakchum, Chota Shigri, and Bara Shigri glaciers, revealing an average velocity of 10.2, 20.9, and 25.3 m/yr, respectively, using ASTER images. The field-based measurements of the surface velocity of the Batal glacier, which is in a similar size class to the Harmukh glacier, showed an average velocity of 6.2 m/yr [151]. The surface flow of the glaciers mainly depends upon the ice thickness, surface gradient, and morphology of the ice surface. The surface velocity was higher near the crevasses (slope break), as illustrated in the red and yellow squares shown on the velocity map and validated on the Google Earth imagery. At the same time, the homogenous flow was observed on the gentle slopes.

5.5. Depth, Volume, and Peak Discharge Estimation

5.5.1. Depth

The depth of the lake was calculated using three empirical equations given in Table 5. The results show a huge variation; therefore, an average of all these equations was taken as the mean depth of Gangabal lake, as suggested by Emmer and Vilmek [114] and Chowdhury et al. [115]. The results further suggest that the depth derived from equations given by Patel et al. [111] and Fujita et al. [37] shows higher values as compared to Huggel et al. [108].

5.5.2. Volume

The volume of the proglacial lake (Gangabal) was derived from five empirical equations developed in different regions of the world. The results show that the volume obtained from the equations given by Evans [113], Huggel et al. [108], and Liu et al. [112] are almost in a similar range. However, volume calculated through the equations of Fujita et al. [37] and Patel et al. [111] shows overestimated values of 102.4 × 106 and 117 × 106 m3. On the other hand, the volume obtained by Sharma et al. 2018 shows a value of 1.54 × 106 m3, which indicates underestimated values (Table 5). Thus, in this study, an average of all these equations was taken from the volume of the proglacial lake. The volume of Gangabal lake calculated through an average of all these equations was estimated as 73.76 × 106 m3, which is approximately equal to the volume estimates derived from the recently developed equation given by Qi et al. [117]. Therefore, it can be concluded that an empirical equation for volume estimation given by Qi et al. [117] is suitable for estimating the volume of those lakes that are not accessible for the bathymetry survey.

5.5.3. Peak Discharge

To ascertain the peak discharge of the glacial lake, three equations given in Table 5 were employed. Average volume, i.e., 73.54 × 106 m3 calculated from the equations given in Table 5, was used as input data to obtain the peak discharges from these equations individually. However, like depth and volume, peak discharge values derived from the empirical equations show a huge variation. Therefore an average of all these equations was taken as a maximum peak discharge of the lake. Allen et al. [75] have suggested that the equation of Fujita et al. [37] is more reliable than the equations given by Huggel et al. [108] and Popov [116] because, in the case of the latter, the peak discharge is overestimated. Although we did not consider the equation of Fujita et al. [37] for the peak flow, instead, we used the average of all equations as suggested by Emmer and Vilmek [114] and Chowdhury et al. [115].
Table 5. Results (depth, volume, and peak discharge) are derived from different empirical equations.
Table 5. Results (depth, volume, and peak discharge) are derived from different empirical equations.
S. No.SourceDepthVolume 106 m3Peak Discharge (Qmax m3)
1Evans [113]-74.14 10,657.59
2Huggel et al. [108]42.4970.44-
3Fujita et al. [37]62.33102.477,409.92
4Patel et al. [111]71.05117-
5Sharma et al. 2018-1.54-
6Liu et al. [112]-77.69-
7Popov [116]--53,883.24
8Qi et al. [117]-73.12-
9Average58.6273.7647,316.91

5.6. GLOF Susceptibility

As a result of the current glacier recession, proglacial lakes are forming and expanding all over the world. Between 1980 and 2018, the number and area of proglacial lakes in the world rose by almost 50% [152]. With the ongoing climate warming, hazards related to glaciers and glacial lakes are expected to increase in the future. GLOFs pose a severe threat to the people living in the mountain regions of the world. Therefore, it has received serious attention in recent years, primarily because of catastrophic damage and fatalities associated with them. The Himalayan region has witnessed several glacier-associated flood events [153] which have flooded the downstream valleys. The valley of Kashmir, where the present lake is located, is no exception. It has also witnessed many flood events dating back to July 1959, when a massive glacial flood hit the valley because of a persistent rain for four consecutive days [154], and a recent being the September 2014 flood event [155]. Because of previous GLOF episodes in Jammu and Kashmir (including Ladakh) and the flood events, it is necessary to monitor glacial lakes associated with glaciers in order to identify lakes that may cause downstream area destruction. According to Mal et al. [43], Jammu and Kashmir have the highest aggregated level of GLOF hazard, followed by Arunachal Pradesh (AP), Sikkim, Himachal Pradesh (HP), and Uttarakhand. A recent study by Dubey and Goyal [44] also categorized JK at the second position in terms of the GLOF hazard level, whereas Sikkim is on the top and followed by HP and AP. Gupta et al. 2021 reported four lakes as potentially dangerous in the Upper Jhelum basin. The present lake has been classified as a potentially dangerous lake by various researchers in the past, for example, a high-risk lake by Gupta et al. [36], potentially dangerous lake by Guru et al. [74], medium-risk lake by Ahmed et al. [41], and important GLOF lakes by Kumar and Goel [42] and Kumar et al. [15].
The lake has an area of 1.66 km2 according to recent estimates (2020), which is greater than the threshold of 0.1 km2 reported in past studies for categorizing the potentially dangerous glacial lake [31,156,157,158,159]. The lake has expanded by 17% during the period 1972 to 2020. The lake is fed by the Harmukh glacier, which has retreated more than 15% during the past 48 years. With the increasing temperature and glacier recession, the area of the lake is anticipated to increase in the coming years. The water supply from Gangabal lake directly flows into Nund-Kol lake, located in the downstream region at a distance of 880 m from the lake outlet. In case of a lake breach outburst event or overtopping failure, the presence of a lake in the downstream course of Gangabal lake may witness a cascading effect, resulting in a catastrophic disaster, as the breaching of one lake can cause the bursting of another. The average slope between the lake and the glacier was estimated at 20.30 degrees (Figure 11a). A steep slope between a lake and its feeding glacier may cause mass movement and ice calving, leading to dam collapse [13,76,78]. The lake is also susceptible to potential rock/ice avalanches and other mass movements owning to its steep slope between the lake and glacier snout (Figure 11b). Rock and ice avalanche zones are visible on high-resolution Google Earth imagery. Glacial lakes with water volumes more than 10 × 106 m3 are regarded as dangerous as they may cause a devastating flash flood downstream. The volume of Gangabal lake estimated using the improved volume estimation equation reveals 73 × 106 m3 of water in the lake. Overall, the lake has been categorized as a high susceptibility lake. A detailed analysis of the GLOF susceptibility of Gangabal lake is available in Ahmed et al. [41].

5.7. Hazard Assessment of Gangabal Lake

The occurrence of the GLOF event is governed mainly by two factors (1) dam stability and (2) the possibility of a triggering factor. As far as the dam condition of Gangabal lake is concerned, we have considered it partly unstable in this study and, therefore, susceptible to dam breach/failure. The lake is also susceptible to potential ice/rock avalanches owing to the steep slope between the glacier and the glacial lake, as depicted in Figure 11a,b. Typically, triggering factors such as ice/snow avalanches, rock falls, landslides, or glacier-calving processes produce impulse waves that might start an overtopping failure, eventually resulting in a GLOF event [97]. The other factors that may result in the GLOF event are cloudbursts, a sudden expansion of the lake under warming temperatures, and incessant precipitation, which may cause the lake’s water level to rise rapidly, leading to the dam overtopping or a piping failure. The mean water level of Gangabal lake using high-resolution GE imagery and the DEM was measured as 3571 m; the dam is almost flat with its elevation rising from 3576 to 3583.
Using Froehlich’s equation, a series of breach hydrographs were calculated for the different breach scenarios with varied breach widths, viz, 20, 40, 60, and 80 m (Figure 12). All these scenarios produced different peak flood discharge estimates mainly due to the different breach parameters. The peak flood at the time of the lake outburst varied from 10,620.03 m3/s (20 m) to 17,787.79 m3/s (80 m). In the first scenario (worst-case), 100% volume (73 × 106 m3) of the lake was considered to be released for the GLOF routing with a breach formation time (bf) of 40 min and breach width (bw) of 80 m. The resultant breach hydrograph produced a peak flood of 17,787.79 m3/s that was achieved within 10 min of initiation of the breach event. The modeled breach hydrographs at various breach formation times are depicted in Figure 12. In the second potential GLOF scenario, 75% of the lake volume (52 × 106 m3) was modeled with a 40 m breach width and 30 min of breach formation time (Figure 13). In scenario 2, the peak discharge of 9155.36 m3 was attained after 12 min after the breach initiation.
Two-dimensional hydraulic routing of generated breach hydrographs using HEC-RAS 5.0.7 was performed for both scenarios (a) scenario 1 (bf = 40 min and bw = 80 m) and (b) scenario 2 (bf = 30 min and bw = 40 m) to assess the potential impacts of GLOFs from Gangabal lake on the downstream valley of the lake up to a distance of 29.3 km. The breach flood, i.e., GLOFs generated by the model, was routed through the river channel, including various locations (villages, bridges, roads, agricultural fields, etc.).
The assessment of the GLOF hydraulics in worst case scenario (scenario 1) revealed that maximum flood depth varies between 3.87 and 68 m, maximum flow velocity between 4 and 75 m/s, and maximum water surface elevation varies between 1548 and 3536 m (Figure 14a–c). The resultant flood potentially inundated areas at nine settlement locations, including 98–104 buildings. It is pertinent to mention that most buildings and other structures were inundated at the end of the channel, which is within the elevation zone of 1726 to 1869 m asl. However, depth, velocity, and peak flow were at their highest in the upstream regions due to the steep slope of the channel walls. The spatial distribution of flood inundation is shown in Figure 15, Figure 16 and Figure 17.
The peak flood and other important breach flood wave parameters such as flood depth, flood velocity, and WSE derived from the model output at important cross-sections (locations) were found to be maximum at Naranaag temple (XCS-45500), situated at a distance of 11.67 km from the lake site and attained at 1:00 h after the breach initiation. The peak flood was 13,995.36 m3 s−1 for the same cross-section (Naranaag temple) with a maximum flood depth of 14.43 m and maximum flow velocity of 8.83 m/s. The details of the derived parameters at each important location (cross-section) are provided in Table 6, whereas GLOF hydrographs are provided in Figure 17.
At 1 h 10 min, an outburst flood with a maximum discharge of 12,896.52 m3 s−1 and maximum flood depth and velocity of 10.54 m and 10.05 m/s would hit the Khanan village region, located 16.63 km from the lake site. The unnamed settlement located 17.67 km from the lake outlet, including 13 buildings and an agricultural area (see Figure 16a,b) was potentially inundated with a peak flood of 11,708.66 m3 s−1, which would be attained at 1:10 with a maximum flood depth of 13.23 m and maximum flood velocity of 8.46 m/s.
The flood propagation parameters were also analyzed for the most important site at which the routed channel meets with the main river, i.e., Sind, located 24.76 km from the outlet of the selected lake. The maximum peak discharge at this location was modeled as 10,888.83 m3 s−1 with maximum flood depth and velocity of 9.21 m and 5.86 m/s, together with a maximum WSE of 1762 m. At this confluence, the maximum inundation was observed which includes 20–25 buildings, two bridges, one recreational park, and a massive agricultural area, all of which came under the flood wave.
The analysis further revealed that flood depths at various cross-sections shows sudden fluctuations: the flood height increases at 13.37 km from the lake outlet, then suddenly decreases to 9.2 m at a distance of 14.67 km. Again, the flood heights gradually rise to 17.6 m near the Baba Nagre Khanam bridge and then gradually decrease to 3.8 m at the Sind river confluence. A similar phenomenon has been observed in the GLOF-routed channel of Thulagi lake [46,160] and Lower Barun lake [160]. This phenomenon resulted due to the narrow and shallow depth of the routed channel. At the narrow channel locations, the depth, velocity, and peak flood were observed as high, while it decreases at the locations where the width of the river channel is broader, which results in an increase in water spread and, thereby, a decrease in the magnitude of flood wave parameters. Further downstream, the peak discharge of 10,848 m3 s−1 was estimated at bridge site 5 with 8.98 m flood depth and 6.86 m/s flood velocity, respectively. At the last selected location, i.e., Drug Tanga village, the flood discharge was estimated as 10,223.98 m3 s−1, whereas flood height and flood velocity were observed as 9.26 m and 5.37 m/s, respectively, shown in Figure 17.
In the second scenario, the maximum flow depth and velocity were modeled as 48 m and 52 m/s, respectively, whereas the maximum WSE rises to 3546 m at the lake site, which decreases to 1689 m at the Drag Tanga village. The peak flood of 5231 m3 s−1 was estimated at the Naranaag temple cross-section with a peak depth and velocity of 9 m and 5.2 m/s, respectively. The total number of inundated buildings in the second scenario decreased to 75–80. For example, in the case of scenario 1, a total number of nine buildings were inundated near the unnamed settlement, as depicted in Figure 18a. However, in the second scenario, only four buildings were potentially inundated due to the resultant GLOF (Figure 18b). In the first scenario, the total inundated area was estimated as 5.31 km2, which also decreased to 3.46 km2 in scenario 2. The inundation map derived from scenario 2 is depicted in Figure 18. Flow depth, velocity, and peak discharge decrease as the routed flow reaches the downstream areas having lower elevations. The analysis revealed that the breach hydrograph is the primary factor that governs the depth, velocity, WSE, and peak flow in the downstream channel. The maximum flow depth and velocity decrease as the Manning coefficient is increased. Flow velocity was observed higher at the places with steep slopes because it is dependent upon the slope of the channel. Depth was found highest in the upstream locations compared to the downstream areas; in the former case, the width of the channel is less, while it is more in the latter case. The channel widens as it proceeds toward the low-elevation regions, thus providing a sufficient water spread area.
The maximum area is inundated in the downstream region due to the gentle slope and greater channel width, but values of flood wave parameters remain minimum in such areas. The investigation showed that the flow hydraulics of a routed GLOF is determined by the initial breach hydrograph. The initial breach hydrograph significantly impacts the hydraulic characteristics of a GLOF at any point downstream. The flow hydraulics were also examined in relation to the channel’s top width. The top width of the channel increases in a similar manner to the flow velocity, decreasing both inundation depth and flow velocity [161].

5.8. Downstream Impact Assessment

The threat of GLOFs is increasing not only in the hilly terrains but also in the downstream valley regions. Therefore, a thorough evaluation of the GLOF hazard, risk, and susceptibility using reliable scientific methodologies is the current need (Khanal et al. 2015). In this study, after evaluating the GLOF impacts on the downstream course of the river and its adjoining areas, the area under each inundated landuse class was estimated through the LULC map generated for both scenarios (scenario 1 and 2). The LULC map prepared for scenarios 1 and 2 is shown in Figure 19a,b. The analysis reveals that in the worst-case scenario, 5.3 km2 of the geographical area would be potentially inundated due to the outburst flood from the lake outlet to the last selected location (Drag Tanga village).
The inundated landuse classes in scenario 1 in terms of area and percentage include water area 0.8 km2 (15.09%), forestland 0.7 km2 (13.2%), grassland 0.4 km2 (7.55%), cropland 0.6 km2 (11.32%), scrub/shrub 1.1 km2 (20.75%), built-up area 0.3 km2 (5.66%), and barren 1.4 km2 (26.42%). In the second scenario, wherein 75% of the lake water volume was considered for the GLOF routing, the inundated area decreased by 1.84 km2 which is 34% of the total inundated area in scenario 1 (Table 7). The total area potentially inundated in scenario 2 was estimated as 3.46 km2, which includes a water area of 0.58 km2 (16.76%), forestland of 0.51 km2 (14.73%), grassland of 0.43 km2 (12.42%), cropland 0.35 km2 (10.11%), scrub/shrub 0.62 km2 (17.91%), built-up area 0.19 km2 (5.49%), and barren 0.78 km2 (22.54%). The investigation shows that the exposure level was greater in the low-lying areas compared to the areas situated at higher elevations. Validation and calibration are crucial parameters in flood inundation models to evaluate and match inundation extents and depths. However, owing to the lack of field data and observed flow data, it is often difficult to validate the hydraulic properties derived from the flooding models. The derived hydraulic calculations cannot be validated in the present study because no past GLOF event has been witnessed or reported for Gangabal lake, although the model reliability has been verified with various studies that have documented GLOF events in various parts of the world, such as Peru, China, India, Nepal, and Tibet.
Glacial lakes are susceptible to breach outbursts under the conditions of rapid warming and incessant rainfall (extreme weather events). GLOFs may also occur due to the sudden ice/rock avalanche, landslides, cloud bursts, glacier surges, ice core melting, and glacier snout calving [33]. The temperatures in the entire Indian Himalayan region are significantly increasing, as evident from recent studies [45,162,163]. Averaged temperature mean, minimum, and maximum over India are showing a warming trend with an increase of 0.15, 0.15, and 0.13 °C per decade, respectively, between 1986 and 2015 [164]. The changing precipitation regimes and increasing temperatures over the Indian Himalayan region are major concerns for the region’s snow cover and glacier health [10]. The temperatures in the Jhelum basin, where Gangabal lake and its feeding glacier are located, have shown a significant increasing trend and decreasing precipitation trend [84,165,166], while glaciers, on the other hand, are rapidly retreating [41,84,166], and glacial lakes are expanding, and their frequency is increasing [39,40,41]. The warming temperatures and, thereby, rapid glacier melting in the area may result in sudden glacial lake outburst events, enhancing the GLOF susceptibility of the region.

5.9. Limitations

The HEC-RAS model was successfully applied in this study and in various past studies to obtain the various hydraulic properties (flood depth, flood velocity, water surface profiles) of the outburst flood and aid in the development of inundation maps. Evaluating the accuracy of the outputs derived from the physically based hydraulic models is difficult since they are sensitive to numerous input parameters [125]. The topographic and geomorphological characteristics of the routed channel also affect the nature of GLOFs. The peak flood discharge and flood velocity through a particular channel are directly affected by Manning’s roughness coefficient, which substantially affects GLOFs. The increase or decrease in Manning’s “n” values plays a significant role in the peak discharge variations. GLOF modeling through hydrological models is nothing more than an approximation of an event. Despite that, these models can be utilized to analyze a physical phenomenon and its corresponding effects on the development of water resource structures and flood management. Most assumptions in GLOF modeling are based on breach parameters (breach depth and breach width), which impact peak flood discharges and the time of flood arrival. The results contain errors and uncertainties since there is a paucity of field-based bathymetry data and a number of other parameters related to the dam break and the routing of the breach hydrograph.
The use of a DEM with coarse spatial resolution also overestimates or underestimates the flood inundation and the magnitude of the hydraulic parameters. The DEM simply defines the elevation of a point in a given area at a particular spatial resolution. Based on the variation observed in the model output, we suggest considering uncertainties associated with the DEM when evaluating the GLOF flow routing and its potential socioeconomic impacts. The accuracy of the water depth is also influenced by the DEM data used to extract the cross-sections; if the DEM data are very accurate, the error will be minimal. The surface area of the glacial lake, considered a critical input parameter to generate a breach hydrograph, was obtained through the empirical equation. The empirical formulae used to calculate the volume and depth of Gangabal lake are based on the regression function of its area. The depth and volume estimates derived from the empirical equations for Gangabal lake may provide rough estimates and overestimates the values due to its larger area. However, when field-based data are not available or possible to obtain, empirical equations can be used to gather first-hand information. The model’s main limitation lies in its taking into account the clear water while simulating unsteady flow. In actuality, sediment may be entrained in the flow, altering its rheology and reducing the distance traveled. As a result, assuming a clean water flow reduces the validity of our findings, which is a typical issue in GLOF research worldwide. The mechanism of glacial lake bursting and the genesis of glacial lake breaches are not well understood. Despite the various drawbacks mentioned above and assumptions, hydrodynamic modeling is useful because it provides some valuable estimates of the glacial lake breach outburst flood, enables an accurate estimation of the design flood, and offers some guidance for the prevention and mitigation of the potential flood.

6. Conclusions

Generally, glaciers across the Himalayas are retreating, resulting in the formation and growth of numerous glacial lakes under the direct influence of climate change. The glacial lakes are expanding at a remarkable pace. In this study, a detailed study of the potentially dangerous glacial lake, Gangabal, and its feeding glacier, Harmukh, was conducted using the integrated approaches of remote sensing, GIS, and hydraulic modeling using HEC-RAS. Our data analysis reveals that the glacial lake has gradually expanded from 1.42 km2 in 1972 to 1.46 km2 in 1981, 1.58 km2 in 1992, 1.61 km2 in 2001, 1.64 km2 in 2010, and 1.66 km2 in 2020. Overall, the lake area experienced a 17% growth in its area from 1972 to 2020 at the rate of 0.005 km2 per year. On the other hand, the glacier has shown a substantial decrease in its area, from 3.6 km2 in 1990 to 3.3 km2 in 2000, 3.1 km2 in 2010, and 2.9 km2 in 2020. The region observed a deglaciation of 0.7 km2 (19 ± 5%) since 1990 at the rate of 0.023 km2 per year. The results derived from the DEM differencing show an average ice thickness loss of 11.04 ± 4.8 m for the Harmukh glacier at the rate of 0.92 ± 0.40 m/yr during the years 2000 and 2012. The minimum and maximum ice thickness for the Harmukh glacier derived from the GlabTop-2 model was estimated to be 7.3 and 85 m, respectively, whereas the mean thickness was calculated as 23.46 m. The Cosi–Corr model yielded a mean velocity of 3.2 m/yr using the pre-event (2019) and post-event (2020) pan bands with a spatial resolution of 15 m from Landsat-8. The minimum and maximum surface velocity varied between 0 and 7 m/yr. Assessment of GLOF propagation in the worst-case scenario (scenario 1) revealed that the maximum flood depth varies between 3.87 and 68 m, the maximum flow velocity ranges between 4 and 75 m/s, and the maximum water surface elevation varies between 1548 and 3536 m. The resultant flood wave in the worst-case scenario would reach the nearest location (Naranaag temple) within 1 h 10 min after the initiation of the breach with a maximum flood discharge of 12,896.52 m3 s−1 and maximum flood depth and velocity of 10.54 m and 10.05 m/s, respectively. After the GLOF impact assessment on downstream exposed areas, the area under each inundated landuse class was estimated through the LULC map generated for both scenarios 1 and 2. In the first scenario, the total potentially inundated area was estimated as 5.3 km2, which was reduced to 3.46 km2 in scenario 2. Owing to the tremendous hazard potential and incidence of GLOFs in the Himalayan region, it is crucial to consider the probable impacts of GLOFs while planning, designing, and building key infrastructure along the flow channels. This is especially important when constructing dams, as they are located on the flow path of a river channel. We attempted a preliminary assessment of both the Harmukh glacier and Gangabal glacial lake by applying an integrated research approach. We suggest utilizing high-resolution satellite data, climate data, field investigations, and robust hydrological models to comprehensively investigate benchmark glaciers and dangerous glacial lakes in the region.

Author Contributions

R.A.: original design of the article, data analysis, writing of original draft, and validation of findings; M.R., G.F.W. and S.T.A.: investigation, writing, formal analysis, validation, review, and editing; R.A.M., G.M. and A.F.R.: investigation, conceptualization, methodology, and review—editing; P.A., S.K.J., G.M. and M.F.: supervision, writing—editing; P.A., S.K.J., G.M. and M.F.: final review. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data will be available from the corresponding author on request.

Acknowledgments

The authors of this original research work are very much indebted to the United States Geological Survey (USGS) for freely providing the satellite data to carry out this study. The first authors especially acknowledge the support of the University Grants Commission in the form of grants made available under the Maulana Azad National Fellowship scheme for pursuing the Ph.D. We would like to thank all the three anonymous reviewers whose constructive comments helped us to improve our manuscript and special thanks to Academic Editor for providing valuable suggestions for improving this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kääb, A.; Frauenfelder, R.; Roer, I. On the response of rockglacier creep to surface temperature increase. Glob. Planet. Chang. 2007, 56, 172–187. [Google Scholar] [CrossRef]
  2. Lama, S.; Devkota, B. Vulnerability of mountain communities to climate change and adaptation strategies. J. Agric. Environ. 2009, 10, 76–83. [Google Scholar] [CrossRef] [Green Version]
  3. Bolch, T.; Peters, J.; Yegorov, A.; Pradhan, B.; Buchroithner, M.; Blagoveshchensky, V. Identification of potentially dangerous glacial lakes in the northern Tian Shan. In Terrigenous Mass Movements; Springer: Berlin/Heidelberg, Germany, 2012; pp. 369–398. [Google Scholar]
  4. Allen, S.K.; Linsbauer, A.; Randhawa, S.S.; Huggel, C.; Rana, P.; Kumari, A. Glacial lake outburst flood risk in Himachal Pradesh, India: An integrative and anticipatory approach considering current and future threats. Nat. Hazards 2016, 84, 1741–1763. [Google Scholar] [CrossRef] [Green Version]
  5. Kaushik, S.; Singh, T.; Bhardwaj, A.; Joshi, P.K. Long-term spatiotemporal variability in the surface velocity of Eastern Himalayan glaciers, India. Earth Surf. Process. Landforms 2022. [Google Scholar] [CrossRef]
  6. IPCC. Climate Change 2007: The Physical Science Basis, Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2007; p. 996. [Google Scholar]
  7. Kääb, A.; Berthier, E.; Nuth, C.; Gardelle, J.; Arnaud, Y. Contrasting patterns of early twenty-first-century glacier mass change in the Himalayas. Nature 2012, 488, 495–498. [Google Scholar] [CrossRef]
  8. Mir, R.A.; Jain, S.K.; Saraf, A.K.; Goswami, A. Detection of changes in glacier mass balance using satellite and meteorological data in Tirungkhad basin located in Western Himalaya. J. Indian Soc. Remote Sens. 2014, 42, 91–105. [Google Scholar] [CrossRef]
  9. Shrestha, A.B.; Mukherji, A.; Mishra, A.; Wester, P. The Hindu Kush Himalaya Assessment; Springer Nature: Berlin, Germany, 2019. [Google Scholar]
  10. Negi, H.S.; Kumar, A.; Kanda, N.; Thakur, N.K.; Singh, K.K. Status of glaciers and climate change of East Karakoram in early twenty-first century. Sci. Total Environ. 2021, 753, 141914. [Google Scholar] [CrossRef]
  11. Mehta, M.; Dobhal, D.P.; Bisht, M.P.S. Change of Tipra glacier in the Garhwal Himalaya, India, between 1962 and 2008. Prog. Phys. Geogr. 2011, 35, 721–738. [Google Scholar] [CrossRef]
  12. Bahuguna, I.M.; Rathore, B.P.; Brahmbhatt, R.; Sharma, M.; Dhar, S.; Randhawa, S.S.; Kumar, K.; Romshoo, S.; Shah, R.D.; Ganjoo, R.K.; et al. Are the Himalayan glaciers retreating? Curr. Sci. 2014, 106, 1008–1013. [Google Scholar]
  13. Mir, R.A.; Jain, S.K.; Lohani, A.K.; Saraf, A.K. Glacier recession and glacial lake outburst flood studies in Zanskar basin, western Himalaya. J. Hydrol. 2018, 564, 376–396. [Google Scholar] [CrossRef]
  14. Ahmed, R.; Wani, G.F.; Ahmad, S.T.; Sahana, M.; Singh, H.; Ahmed, P. A Review of Glacial Lake Expansion and Associated Glacial Lake Outburst Floods in the Himalayan Region. Earth Syst. Environ. 2021, 5, 695–708. [Google Scholar] [CrossRef]
  15. Kumar, V.; Shukla, T.; Mehta, M.; Dobhal, D.P.; Bisht, M.P.S.; Nautiyal, S. Glacier changes and associated climate drivers for the last three decades, Nanda Devi region, Central Himalaya, India. Quat. Int. 2021, 575, 213–226. [Google Scholar] [CrossRef]
  16. Nie, Y.; Sheng, Y.; Liu, Q.; Liu, L.; Liu, S.; Zhang, Y.; Song, C. A regional-scale assessment of Himalayan glacial lake changes using satellite observations from 1990 to 2015. Remote Sens. Environ. 2017, 189, 1–13. [Google Scholar] [CrossRef] [Green Version]
  17. Nie, Y.; Liu, Q.; Wang, J.; Zhang, Y.; Sheng, Y.; Liu, S. An inventory of historical glacial lake outburst floods in the Himalayas based on remote sensing observations and geomorphological analysis. Geomorphology 2018, 308, 91–106. [Google Scholar] [CrossRef]
  18. Bhutiyani, M.R.; Kale, V.S.; Pawar, N.J. Long-term trends in maximum, minimum and mean annual air temperatures across the Northwestern Himalaya during the twentieth century. Clim. Chang. 2007, 85, 159–177. [Google Scholar] [CrossRef]
  19. Shrestha, U.B.; Gautam, S.; Bawa, K.S. Widespread climate change in the Himalayas and associated changes in local ecosystems. PLoS ONE 2012, 7, e36741. [Google Scholar] [CrossRef] [Green Version]
  20. Negi, H.S.; Shekhar, M.S.; Gusain, H.S.; Ganju, A. Winter climate and snow cover variability over north-west Himalaya. In Science and Geopolitics of the White World; Springer: Cham, Switzerland, 2018; pp. 127–142. [Google Scholar]
  21. Sabin, T.P.; Krishnan, R.; Vellore, R.; Priya, P.; Borgaonkar, H.P.; Singh, B.B.; Sagar, A. Climate change over the Himalayas. In Assessment of Climate Change over the Indian Region; Springer: Singapore, 2020; pp. 207–222. [Google Scholar]
  22. Nie, Y.; Pritchard, H.D.; Liu, Q.; Hennig, T.; Wang, W.; Wang, X.; Liu, S.; Nepal, S.; Samyn, D.; Hewitt, K.; et al. Glacial change and hydrological implications in the Himalaya and Karakoram. Nat. Rev. Earth Environ. 2021, 2, 91–106. [Google Scholar] [CrossRef]
  23. Sharma, R.K.; Kumar, R.; Pradhan, P.; Sharma, A. Climate-Induced Glacier Retreats and Associated Hazards: Need for Robust Glaciers and Glacial Lake Management Policy in Sikkim Himalaya, India. In Climate Change; Rani, S., Kumar, R., Eds.; Springer: Cham, Switzerland, 2022; pp. 161–182. [Google Scholar]
  24. Shugar, D.H.; Burr, A.; Haritashya, U.K.; Kargel, J.S.; Watson, C.S.; Kennedy, M.C.; Bevington, A.R.; Betts, R.A.; Harrison, S.; Strattman, K. Rapid worldwide growth of glacial lakes since 1990. Nat. Clim. Chang. 2020, 10, 939–945. [Google Scholar] [CrossRef]
  25. Pandey, P.; Ali, S.N.; Champati Ray, P.K. Glacier-glacial lake interactions and glacial lake development in the central Himalaya, India (1994–2017). J. Earth Sci. 2021, 32, 1563–1574. [Google Scholar] [CrossRef]
  26. Li, D.; Shangguan, D.; Wang, X.; Ding, Y.; Su, P.; Liu, R.; Wang, M. Expansion and hazard risk assessment of glacial lake Jialong Co in the central Himalayas by using an unmanned surface vessel and remote sensing. Sci. Total Environ. 2021, 784, 147249. [Google Scholar] [CrossRef] [PubMed]
  27. Haritashya, U.K.; Kargel, J.S.; Shugar, D.H.; Leonard, G.J.; Strattman, K.; Watson, C.S.; Shean, D.; Harrison, S.; Mandli, K.T.; Regmi, D. Evolution and controls of large glacial lakes in the Nepal Himalaya. Remote Sens. 2018, 10, 798. [Google Scholar] [CrossRef]
  28. Kraaijenbrink, P.D.; Bierkens, M.F.P.; Lutz, A.F.; Immerzeel, W.W. Impact of a global temperature rise of 1.5 degrees Celsius on Asia’s glaciers. Nature 2017, 549, 257–260. [Google Scholar] [CrossRef]
  29. Vuichard, D.; Zimmermann, M. The 1985 catastrophic drainage of a moraine-dammed lake, Khumbu Himal, Nepal: Cause and consequences. Mt. Res. Dev. 1987, 7, 91–110. [Google Scholar] [CrossRef]
  30. Fujita, K.; Suzuki, R.; Nuimura, T.; Sakai, A. Performance of ASTER and SRTM DEMs, and their potential for assessing glacial lakes in the Lunana region, Bhutan Himalaya. J. Glaciol. 2008, 54, 220–228. [Google Scholar] [CrossRef] [Green Version]
  31. Jha, L.K.; Khare, D. Detection and delineation of glacial lakes and identification of potentially dangerous lakes of Dhauliganga basin in the Himalaya by remote sensing techniques. Nat. Hazards 2017, 85, 301–327. [Google Scholar] [CrossRef]
  32. Ives, J.D.; Shrestha, R.B.; Mool, P.K. Formation of Glacial Lakes in the Hindu Kush-Himalayas and GLOF Risk Assessment; ICIMOD: Kathmandu, Nepal, 2010; pp. 10–11. [Google Scholar]
  33. Wang, Y.; Li, J.; Wu, L.; Guo, L.; Hu, J.; Zhang, X. Estimating the Changes in Glaciers and Glacial Lakes in the Xixabangma Massif, Central Himalayas, between 1974 and 2018 from Multisource Remote Sensing Data. Remote Sens. 2021, 13, 3903. [Google Scholar] [CrossRef]
  34. Chen, F.; Zhang, M.; Guo, H.; Allen, S.; Kargel, J.S.; Haritashya, U.K.; Watson, C.S. Annual 30 m dataset for glacial lakes in High Mountain Asia from 2008 to 2017. Earth Syst. Sci. Data 2021, 13, 741–766. [Google Scholar] [CrossRef]
  35. CWC. June 2016 Report on “Monitoring of Glacial Lakes/Water Bodies in the Himalayan Region of Indian River Basins for June 2016; CWC: Roma, Italy, 2016. [Google Scholar]
  36. Gupta, A.; Maheshwari, R.; Guru, N.; Sweta, R.B.S.; Raju, P.V.; Rao, V.V. Quantitative prioritization of potentially critical glacial Lakes in the Indus River basin using satellite derived parameters. Geocarto Int. 2021, 1–23. [Google Scholar] [CrossRef]
  37. Fujita, K.; Sakai, A.; Takenaka, S.; Nuimura, T.; Surazakov, A.B.; Sawagaki, T.; Yamanokuchi, T. Potential flood volume of Himalayan glacial lakes. Nat. Hazards Earth Syst. Sci. 2013, 13, 1827–1839. [Google Scholar] [CrossRef] [Green Version]
  38. Worni, R.; Huggel, C.; Stoffel, M. Glacial lakes in the Indian Himalayas—From an area-wide glacial lake inventory to on-site and modeling based risk assessment of critical glacial lakes. Sci. Total Environ. 2013, 468, S71–S84. [Google Scholar] [CrossRef]
  39. Dar, T.; Rai, N.; Kumar, S.; Bhat, M.A. Climate change impact on cryosphere and streamflow in the Upper Jhelum River Basin (UJRB) of north-western Himalayas. Environ. Monit. Assess. 2022, 194, 140. [Google Scholar] [CrossRef]
  40. Ahmed, R.; Wani, G.F.; Ahmad, S.T.; Mir, R.A.; Almazroui, M.; Jain, S.K.; Ahmed, P. Spatiotemporal dynamics of glacial lakes (1990–2018) in the Kashmir Himalayas, India using Remote Sensing and GIS. Discov. Water 2021, 1, 7. [Google Scholar] [CrossRef]
  41. Ahmed, R.; Ahmad, S.T.; Wani, G.F.; Mir, R.A.; Ahmed, P.; Jain, S.K. High resolution inventory and hazard assessment of potentially dangerous glacial Lakes in Upper Jhelum Basin, Kashmir Himalaya, India. Geocarto Int. 2022, 1–32. [Google Scholar] [CrossRef]
  42. Kumar, S.; Goel, N.K. Monitoring of Two Typical Glacier Lakes in Indus Basin. Int. J. Eng. Res. Technol. 2015, 3, 1–7. [Google Scholar]
  43. Mal, S.; Allen, S.K.; Frey, H.; Huggel, C.; Dimri, A.P. Sectorwise Assessment of Glacial Lake Outburst Flood Danger in the Indian Himalayan Region. Mount. Res. Dev. 2021, 41, R1. [Google Scholar] [CrossRef]
  44. Dubey, S.; Goyal, M.K. Glacial lake outburst flood hazard, downstream impact, and risk over the Indian Himalayas. Water Resour. Res. 2020, 56, e2019WR026533. [Google Scholar] [CrossRef]
  45. Dimri, A.P.; Allen, S.; Huggel, C.; Mal, S.; Ballesteros-Canovas, J.A.; Rohrer, M.; Shukla, A.; Tiwari, P.; Maharana, P.; Bolch, T.; et al. Climate change, cryosphere and impacts in the Indian Himalayan Region. Curr. Sci. 2021, 120, 774–790. [Google Scholar] [CrossRef]
  46. Khanal, N.R.; Hu, J.M.; Mool, P. Glacial lake outburst flood risk in the Poiqu/Bhote Koshi/Sun Koshi River basin in the Central Himalayas. Mt. Res. Dev. 2015, 35, 351–364. [Google Scholar] [CrossRef]
  47. Veh, G.; Korup, O.; von Specht, S.; Roessner, S.; Walz, A. Unchanged frequency of moraine-dammed glacial lake outburst floods in the Himalaya. Nat. Clim. Chang. 2019, 9, 379–383. [Google Scholar] [CrossRef]
  48. Gunn, J.P. The Shyok flood, 1929. Himal. J. 1930, 2, 35–47. [Google Scholar]
  49. Richardson, S.D.; Quincey, D.J. Glacier outburst floods from Ghulkin Glacier, upper Hunza Valley, Pakistan. In Proceedings of the EGU General Assembly Conference Abstracts, Vienna, Austria, 19–24 April 2009; p. 12871. [Google Scholar]
  50. Bajracharya, S.R.; Mool, P.K.; Shrestha, B.R. Impact of Climate Change on Himalayan Glaciers and Glacial Lakes: Case Studies on GLOF and Associated Hazards in Nepal and Bhutan; International Centre for Integrated Mountain Development—ICIMOD: Patan, Nepal, 2007. [Google Scholar]
  51. Yamada, T.; Sharma, C.K. Glacier Lakes and Outburst Floods in the Nepal Himalaya; IAHS Publications-Publications of the International Association of Hydrological Sciences: Wallingford, UK, 1993; Volume 218, pp. 319–330. [Google Scholar]
  52. Ghimire, M. Review of studies on glacier lake outburst floods and associated vulnerability in the Himalayas. Himal. Rev. 2004, 35–25, 49–64. [Google Scholar]
  53. Govindha Raj, K.B. Remote sensing based hazard assessment of glacial lakes: A case study in Zanskar basin, Jammu and Kashmir, India. Geomat. Nat. Hazards Risk 2010, 1, 339–347. [Google Scholar] [CrossRef]
  54. Hussain, A.; Nasab, N.; Bano, D.; Karim, D.; Anwar, W.; Hussain, K.; Uddin, N. Glacier lake outburst flood modeling of Khurdopin glacier lake using HEC-RAS and GIS. In Селевые Пoтoки: Катастрoфы, Риск, Прoгнoз, Защита; Universal Publishing: New York, NY, USA, 2020. [Google Scholar]
  55. Xu, D. Characteristics of debris flow caused by outburst of glacial lake in Boqu river, Xizang, China, 1981. GeoJournal 1988, 17, 569–580. [Google Scholar] [CrossRef] [Green Version]
  56. Shrestha, A.B.; Eriksson, M.; Mool, P.; Ghimire, P.; Mishra, B.; Khanal, N.R. Glacial lake outburst flood risk assessment of Sun Koshi basin, Nepal. Geomat. Nat. Hazards Risk 2010, 1, 157–169. [Google Scholar] [CrossRef]
  57. Chen, N.S.; Hu, G.S.; Deng, W.; Khanal, N.; Zhu, Y.H.; Han, D. On the water hazards in the trans-boundary Kosi River basin. Nat. Hazards Earth Syst. Sci. 2013, 13, 795–808. [Google Scholar] [CrossRef] [Green Version]
  58. Wang, W.; Gao, Y.; Anacona, P.I.; Lei, Y.; Xiang, Y.; Zhang, G.; Li, S.; Lu, A. Integrated hazard assessment of Cirenmaco glacial lake in Zhangzangbo valley, Central Himalayas. Geomorphology 2018, 306, 292–305. [Google Scholar] [CrossRef]
  59. Watanbe, T.; Rothacher, D. The 1994 Lugge Tsho glacial lake outburst flood, Bhutan Himalaya. Mt. Res. Dev. 1996, 16, 77–81. [Google Scholar] [CrossRef]
  60. Komori, J.; Koike, T.; Yamanokuchi, T.; Tshering, P. Glacial lake outburst events in the Bhutan Himalayas. Glob. Environ. Res. 2012, 16, 59–70. [Google Scholar]
  61. Clague, J.J.; Huggel, C.; KOrup, O.; MCguire, B. Climate change and hazardous processes in high mountains. Rev. De La Asoc. Geológica Argent. 2012, 69, 328–338. [Google Scholar]
  62. Huss, M.; Bookhagen, B.; Huggel, C.; Jacobsen, D.; Bradley, R.S.; Clague, J.J.; Vuille, M.; Buytaert, W.; Cayan, D.R.; Greenwood, G.; et al. Toward mountains without permanent snow and ice. Earth’s Future 2017, 5, 418–435. [Google Scholar] [CrossRef]
  63. Clague, J.J.; O’Connor, J.E. Glacier-related outburst floods. In Snow and Ice-Related Hazards, Risks, and Disasters; Elsevier: Amsterdam, The Netherlands, 2021; pp. 467–499. [Google Scholar]
  64. Hewitt, K. Environmental disasters in social context: Toward a preventive and precautionary approach. Nat. Hazards 2012, 66, 3–14. [Google Scholar] [CrossRef]
  65. Schwanghart, W.; Worni, R.; Huggel, C.; Stofel, M.; Korup, O. Uncertainty in the Himalayan energy–water nexus: Estimating regional exposure to glacial lake outburst foods. Environ. Res. Lett. 2016, 11, 74005. [Google Scholar] [CrossRef] [Green Version]
  66. ICIMOD. Glacial Lakes and Glacial Lake Outburst Floods in Nepal; ICIMOD: Kathmandu, Nepal, 2011. [Google Scholar]
  67. Song, C.; Sheng, Y.; Wang, J.; Ke, L.; Madson, A.; Nie, Y. Heterogeneous glacial lake changes and links of lake expansions to the rapid thinning of adjacent glacier termini in the Himalayas. Geomorphology 2017, 280, 30–38. [Google Scholar] [CrossRef] [Green Version]
  68. Richardson, S.D.; Reynolds, J.M. An overview of glacial hazards in the Himalayas. Quat. Int. 2000, 65, 31–47. [Google Scholar] [CrossRef]
  69. Iturrizaga, L. Glacier Lake outburst floods. In Encyclopedia of Snow, Ice and Glaciers; Singh, V.P., Singh, P., Haritashya, U.K., Eds.; Springer Science: Dordrecht, The Netherlands, 2011; pp. 381–399. [Google Scholar]
  70. Gilany, N.; Iqbal, J. Geospatial analysis and simulation of glacial lake outburst flood hazard in Shyok Basin of Pakistan. Environ. Earth Sci. 2020, 79, 139. [Google Scholar] [CrossRef]
  71. Kumar, D.; Singh, A.K.; Singh, D.S. Spatio-temporal fluctuations over Chorabari glacier, Garhwal Himalaya, India between 1976 and 2017. Quat. Int. 2021, 575, 178–189. [Google Scholar] [CrossRef]
  72. Srivastava, P.; Bhambri, R.; Kawishwar, P.; Dobhal, D.P. Water level changes of high altitude lakes in Himalaya–Karakoram from ICESat altimetry. J. Earth Syst. Sci. 2013, 122, 1533–1543. [Google Scholar] [CrossRef] [Green Version]
  73. Rao, B.S.; Gupta, A.; Guru, N.; Maheshwari, S.R.; Raju, P.V.; Rao, V.V. Glacial Lake Atlas of Indus River Basin; National Remote Sensing Centre, ISRO: Hyderabad, India, 2020; pp. 1–273. [Google Scholar]
  74. Guru, N.; Gupta, A.; Maheshwari, R.; Rao, B.S.; Raju, P.V.; Rao, V.V. Identification of potential lakes susceptible to GLOF in Jhelum Subbasin of Indus Basin. Hydraulics 2019, 1452–1462. [Google Scholar]
  75. Allen, S.K.; Zhang, G.; Wang, W.; Yao, T.; Bolch, T. Potentially dangerous glacial lakes across the Tibetan Plateau revealed using a large-scale automated assessment approach. Sci. Bull. 2019, 64, 435–445. [Google Scholar] [CrossRef]
  76. Che, T.; Xiao, L.; Liou, Y.A. Changes in glaciers and glacial lakes and the identification of dangerous glacial lakes in the Pumqu River Basin, Xizang (Tibet). Adv. Meteorol. 2014, 2014, 903709. [Google Scholar] [CrossRef] [Green Version]
  77. Ahmed, R.; Ahmad, S.T.; Wani, G.F.; Mir, R.A.; Almazroui, M.; Bansal, J.K.; Ahmed, P. Glacial Lake Changes and the Identification of Potentially Dangerous Glacial Lakes (PDGLs) under Warming Climate in the Dibang River Basin, Eastern Himalaya, India; Geocarto International: Hong Kong, China, 2022; pp. 1–24. [Google Scholar]
  78. Jain, S.K.; Lohani, A.K.; Singh, R.D.; Chaudhary, A.; Thakural, L.N. Glacial lakes and glacial lake outburst flood in a Himalayan basin using remote sensing and GIS. Nat. Hazards 2012, 62, 887–899. [Google Scholar] [CrossRef]
  79. Paul, F.; Bolch, T.; Kääb, A.; Nagler, T.; Nuth, C.; Scharrer, K.; Shepherd, A.; Strozzi, T.; Ticconi, F.; Bhambri, R.; et al. The glaciers climate change initiative: Methods for creating glacier area, elevation change and velocity products. Remote Sens. Environ. 2015, 162, 408–426. [Google Scholar] [CrossRef] [Green Version]
  80. Paul, F.; Kääb, A.; Maisch, M.; Kellenberger, T.; Haeberli, W. The new remote-sensing-derived Swiss glacier inventory: I. Methods. Ann. Glaciol. 2002, 34, 355–361. [Google Scholar] [CrossRef] [Green Version]
  81. Frey, H.; Machguth, H.; Huss, M.; Huggel, C.; Bajracharya, S.; Bolch, T.; Kulkarni, A.; Linsbauer, A.; Salzmann, N.; Stoffel, M. Estimating the volume of glaciers in the Himalayan–Karakoram region using different methods. Cryosphere 2014, 8, 2313–2333. [Google Scholar] [CrossRef] [Green Version]
  82. Ahmad, S.T.; Ahmed, R.; Wani, G.F.; Sharma, P.; Ahmed, P.; Mir, R.A.; Alam, J.B. Assessing the Status of Glaciers in Upper Jhelum Basin of Kashmir Himalayas Using Multi-temporal Satellite Data. Earth Syst. Environ. 2022, 6, 375–389. [Google Scholar] [CrossRef]
  83. Bahr, D.B.; Pfeffer, W.T.; Kaser, G. A review of volume-area scaling of glaciers. Rev. Geophys. 2015, 53, 95–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Shafiq, M.U.; Islam, Z.U.; Bhat, I.A.; Ahmed, P. Spatio-temporal behaviour of Nehnar Glacier from 1962 to 2017, Jhelum basin, Kashmir Himalayas, India. Phys. Geogr. 2020, 41, 517–536. [Google Scholar] [CrossRef]
  85. Nye, J.F. The mechanics of glacier flow. J. Glaciol. 1952, 2, 82–93. [Google Scholar] [CrossRef] [Green Version]
  86. Paterson, W.S.B. The sliding velocity of Athabasca glacier, Canada. J. Glaciol. 1970, 9, 55–63. [Google Scholar] [CrossRef] [Green Version]
  87. Haeberli, W.; Hölzle, M. Application of inventory data for estimating characteristics of and regional climate-change effects on mountain glaciers: A pilot study with the European Alps. Ann. Glaciol. 1995, 21, 206–212. [Google Scholar] [CrossRef] [Green Version]
  88. Linsbauer, A.; Paul, F.; Haeberli, W. Modeling glacier thickness distribution and bed topography over entire mountain ranges with GlabTop: Application of a fast and robust approach. J. Geophys. Res. Earth Surf. 2012, 117. [Google Scholar] [CrossRef] [Green Version]
  89. Farinotti, D.; Huss, M.; Bauder, A.; Funk, M.; Truffer, M. A method to estimate the ice volume and ice-thickness distribution of alpine glaciers. J. Glaciol. 2009, 55, 422–430. [Google Scholar] [CrossRef] [Green Version]
  90. Benn, D.I.; Evans, D.J. Glaciers & Glaciation; Routledge: Oxfordshire, UK, 2014. [Google Scholar]
  91. Ayoub, F.; Leprince, S.; Keene, L. User’s Guide to COSI-CORR Co-Registration of Optically Sensed Images and Correlation; California Institute of Technology: Pasadena, CA, USA, 2009; Volume 38, p. 49s. [Google Scholar]
  92. Leprince, S.; Barbot, S.; Ayoub, F.; Avouac, J.P. Automatic and precise orthorectification, coregistration, and subpixel correlation of satellite images, application to ground deformation measurements. IEEE Trans. Geosci. Remote Sens. 2007, 45, 1529–1558. [Google Scholar] [CrossRef] [Green Version]
  93. Romshoo, S.A.; Murtaza, K.O.; Shah, W.; Ramzan, T.; Ameen, U.; Bhat, M.H. Anthropogenic climate change drives melting of glaciers in the Himalaya. Environ. Sci. Pollut. Res. 2022, 29, 52732–52751. [Google Scholar] [CrossRef] [PubMed]
  94. Das, S.; Sharma, M.C. Glacier surface velocities in the Jankar Chhu Watershed, western Himalaya, India: Study using Landsat time series data (1992–2020). Remote Sens. Appl. Soc. Environ. 2021, 24, 100615. [Google Scholar] [CrossRef]
  95. Khadka, M.; Steiner, J.; Hassan, J.; Maussion, F.; Fujita, K.; Gurung, T.; Zhao, C. Performance Evaluation of Existing Algorithms and Datasets for Glacier Surface Velocity Estimation in Himalayan Glaciers. In AGU Fall Meeting Abstracts; AGU: New Orleans, LA, USA, 2021; Volume 2021, p. C13B-04. [Google Scholar]
  96. Pandey, P.; Banerjee, D.; Ray, P.K.C. A satellite-based comprehensive observation of glaciological characteristics of Shunkalpa (Ralam) Glacier, Central Himalaya, India. J. Earth Syst. Sci. 2021, 130, 139. [Google Scholar] [CrossRef]
  97. Sattar, A.; Goswami, A.; Kulkarni, A.V.; Das, P. Glacier-surface velocity derived ice volume and retreat assessment in the dhauliganga basin, central himalaya–A remote sensing and modeling based approach. Front. Earth Sci. 2019, 7, 105. [Google Scholar] [CrossRef]
  98. Scherler, D.; Strecker, M.R. Large surface velocity fluctuations of Biafo Glacier, central Karakoram, at high spatial and temporal resolution from optical satellite images. J. Glaciol. 2012, 58, 569–580. [Google Scholar] [CrossRef] [Green Version]
  99. Ruiz, L.; Berthier, E.; Viale, M.; Pitte, P.; Masiokas, M.H. Recent geodetic mass balance of Monte Tronador glaciers, northern Patagonian Andes. Cryosphere 2017, 11, 619–634. [Google Scholar] [CrossRef] [Green Version]
  100. Huss, M.; Bauder, A.; Funk, M.; Hock, R. Determination of the seasonal mass balance of four Alpine glaciers since 1865. J. Geophys. Res. Earth Surf. 2008, 113. [Google Scholar] [CrossRef]
  101. Wang, P.; Li, Z.; Li, H.; Wang, W.; Yao, H. Comparison of glaciological and geodetic mass balance at Urumqi Glacier No. 1, Tian Shan, central Asia. Glob. Planet. Chang. 2014, 114, 14–22. [Google Scholar] [CrossRef]
  102. Ramsankaran, R.A.A.J.; Pandit, A.; Parla, A. Decadal estimates of surface mass balance for glaciers in Chandra Basin, Western Himalayas, India—A geodetic approach. In Climate Change Signals and Response; Springer: Singapore, 2019; pp. 109–125. [Google Scholar]
  103. Rashid, I.; Majeed, U. Retreat and geodetic mass changes of Zemu Glacier, Sikkim Himalaya, India, between 1931 and 2018. Reg. Environ. Chang. 2020, 20, 125. [Google Scholar] [CrossRef]
  104. Majeed, U.; Rashid, I.; Najar, N.A.; Gul, N. Spatiotemporal dynamics and geodetic mass changes of glaciers with varying debris cover in the pangong region of Trans-Himalayan Ladakh, India between 1990 and 2019. Front. Earth Sci. 2021, 9, 748107. [Google Scholar] [CrossRef]
  105. Muhammad, S.; Tian, L.; Khan, A. Early twenty-first century glacier mass losses in the Indus Basin constrained by density assumptions. J. Hydrol. 2019, 574, 467–475. [Google Scholar] [CrossRef]
  106. Huber, J.; McNabb, R.; Zemp, M. Elevation changes of west-central Greenland glaciers from 1985 to 2012 from remote sensing. Front. Earth Sci. 2020, 8, 35. [Google Scholar] [CrossRef] [Green Version]
  107. Khadka, N.; Zhang, G.; Thakuri, S. Glacial lakes in the Nepal Himalaya: Inventory and decadal dynamics (1977–2017). Remote Sens. 2018, 10, 1913. [Google Scholar] [CrossRef] [Green Version]
  108. Huggel, C.; Kääb, A.; Haeberli, W.; Teysseire, P.; Paul, F. Remote sensing based assessment of hazards from glacier lake outbursts: A case study in the Swiss Alps. Can. Geotech. J. 2002, 39, 316–330. [Google Scholar] [CrossRef]
  109. Shukla, A.; Garg, P.K.; Srivastava, S. Evolution of glacial and high-altitude lakes in the Sikkim, Eastern Himalaya over the past four decades (1975–2017). Front. Environ. Sci. 2018, 6, 81. [Google Scholar] [CrossRef] [Green Version]
  110. Begam, S.; Sen, D. Mapping of moraine dammed glacial lakes and assessment of their areal changes in the central and eastern Himalayas using satellite data. J. Mt. Sci. 2019, 16, 77–94. [Google Scholar] [CrossRef]
  111. Patel, L.K.; Sharma, P.; Laluraj, C.M.; Thamban, M.; Singh, A.; Ravindra, R. A geospatial analysis of Samudra Tapu and Gepang Gath glacial lakes in the Chandra Basin, Western Himalaya. Nat. Hazards 2017, 86, 1275–1290. [Google Scholar] [CrossRef]
  112. Liu, M.; Chen, N.; Zhang, Y.; Deng, M. Glacial lake inventory and lake outburst flood/debris flow hazard assessment after the Gorkha earthquake in the Bhote Koshi Basin. Water 2020, 12, 464. [Google Scholar] [CrossRef] [Green Version]
  113. Evans, S.G. The maximum discharge of outburst floods caused by the breaching of man-made and natural dams. Can. Geotech. J. 1986, 23, 385–387. [Google Scholar] [CrossRef]
  114. Emmer, A.; Vilímek, V.; Zapata, M.L. Hazard mitigation of glacial lake outburst floods in the Cordillera Blanca (Peru): The effectiveness of remedial works. J. Flood Risk Manag. 2018, 11, S489–S501. [Google Scholar] [CrossRef]
  115. Chowdhury, A.; Kroczek, T.; De, S.K.; Vilímek, V.; Sharma, M.C.; Debnath, M. Glacial Lake Evolution (1962–2018) and Outburst Susceptibility of Gurudongmar Lake Complex in the Tista Basin, Sikkim Himalaya (India). Water 2021, 13, 3565. [Google Scholar] [CrossRef]
  116. Popov, N. Assessment of glacial debris flow hazard in the north Tien-Shan. In Soviet-China-Japan Symposium and Field Workshop on Natural Disasters; 1991; pp. 384–391.
  117. Qi, M.; Liu, S.; Wu, K.; Zhu, Y.; Xie, F.; Jin, H.; Gao, Y.; Yao, X. Improving the accuracy of glacial lake volume estimation: A case study in the Poiqu basin, central Himalayas. J. Hydrol. 2022, 610, 127973. [Google Scholar] [CrossRef]
  118. Allen, S.K.; Schneider, D.; Owens, I.F. First approaches towards modelling glacial hazards in the Mount Cook region of New Zealand’s Southern Alps. Nat. Hazards Earth Syst. Sci. 2009, 9, 481–499. [Google Scholar] [CrossRef] [Green Version]
  119. Prakash, C.; Nagarajan, R. Glacial Lake changes and outburst flood hazard in Chandra basin, North-Western Indian Himalayas. Geomat. Nat. Hazards Risk 2018, 9, 337–355. [Google Scholar] [CrossRef]
  120. Watson, C.S.; Carrivick, J.; Quincey, D. An improved method to represent DEM uncertainty in glacial lake outburst flood propagation using stochastic simulations. J. Hydrol. 2015, 529, 1373–1389. [Google Scholar] [CrossRef] [Green Version]
  121. Rounce, D.R.; Watson, C.S.; McKinney, D.C. Identification of hazard and risk for glacial lakes in the Nepal Himalaya using satellite imagery from 2000–2015. Remote Sens. 2017, 9, 654. [Google Scholar] [CrossRef] [Green Version]
  122. Mergili, M.; Schneider, J.F. Regional-scale analysis of lake outburst hazards in the southwestern Pamir, Tajikistan, based on remote sensing and GIS. Nat. Hazards Earth Syst. Sci. 2011, 11, 1447–1462. [Google Scholar] [CrossRef] [Green Version]
  123. Kougkoulos, I.; Cook, S.J.; Edwards, L.A.; Clarke, L.J.; Symeonakis, E.; Dortch, J.M.; Nesbitt, K. Modelling glacial lake outburst flood impacts in the Bolivian Andes. Nat. Hazards 2018, 94, 1415–1438. [Google Scholar] [CrossRef] [Green Version]
  124. Klimeš, J.; Benešová, M.; Vilímek, V.; Bouška, P.; Cochachin Rapre, A. The reconstruction of a glacial lake outburst flood using HEC-RAS and its significance for future hazard assessments: An example from Lake 513 in the Cordillera Blanca, Peru. Nat. Hazards 2014, 71, 1617–1638. [Google Scholar] [CrossRef]
  125. Anacona, P.I.; Mackintosh, A.; Norton, K. Reconstruction of a glacial lake outburst flood (GLOF) in the Engaño Valley, Chilean Patagonia: Lessons for GLOF risk management. Sci. Total Environ. 2015, 527, 1–11. [Google Scholar] [CrossRef]
  126. Petrakov, D.; Shpuntova, A.; Aleinikov, A.; Kääb, A.; Kutuzov, S.; Lavrentiev, I.; Stoffel, M.; Tutubalina, O.; Usubaliev, R. Accelerated glacier shrinkage in the Ak-Shyirak massif, Inner Tien Shan, during 2003–2013. Sci. Total Environ. 2016, 562, 364–378. [Google Scholar] [CrossRef] [Green Version]
  127. Westoby, M.J.; Glasser, N.F.; Brasington, J.; Hambrey, M.J.; Quincey, D.J.; Reynolds, J.M. Modelling outburst floods from moraine-dammed glacial lakes. Earth-Sci. Rev. 2014, 134, 137–159. [Google Scholar] [CrossRef] [Green Version]
  128. Somos-Valenzuela, M.A.; Chisolm, R.E.; Rivas, D.S.; Portocarrero, C.; McKinney, D.C. Modeling a glacial lake outburst flood process chain: The case of Lake Palcacocha and Huaraz, Peru. Hydrol. Earth Syst. Sci. 2016, 20, 2519–2543. [Google Scholar] [CrossRef]
  129. Klimeš, J.; Novotný, J.; Novotná, I.; de Urries, B.J.; Vilímek, V.; Emmer, A.; Strozzi, T.; Kusák, M.; Rapre, A.C.; Hartvich, F.; et al. Landslides in moraines as triggers of glacial lake outburst floods: Example from Palcacocha Lake (Cordillera Blanca, Peru). Landslides 2016, 13, 1461–1477. [Google Scholar] [CrossRef]
  130. Sattar, A.; Haritashya, U.K.; Kargel, J.S.; Leonard, G.J.; Shugar, D.H.; Chase, D.V. Modeling lake outburst and downstream hazard assessment of the Lower Barun Glacial Lake, Nepal Himalaya. J. Hydrol. 2021, 598, 126208. [Google Scholar] [CrossRef]
  131. Hagg, W.; Ram, S.; Klaus, A.; Aschauer, S.; Babernits, S.; Brand, D.; Guggemoos, P.; Pappas, T. Hazard Assessment for a Glacier Lake Outburst Flood in the Mo Chu River Basin, Bhutan. Appl. Sci. 2021, 11, 9463. [Google Scholar] [CrossRef]
  132. Chow, V.T. Open-Channel Hydraulics; McGraw-Hill: New York, NY, USA, 1959; p. 680. [Google Scholar]
  133. Alam, A.; Bhat, M.S.; Maheen, M. Using Landsat satellite data for assessing the land use and land cover change in Kashmir valley. GeoJournal 2020, 85, 1529–1543. [Google Scholar] [CrossRef] [Green Version]
  134. Imdad, K.; Rihan, M.; Sahana, M.; Parween, S.; Ahmed, R.; Costache, R.; Chaudhary, A.; Tripathi, R. Wetland health, water quality, and resident perceptions of declining ecosystem services: A case study of Mount Abu, Rajasthan, India. Environ. Sci. Pollut. Res. 2022, 1–27. [Google Scholar] [CrossRef]
  135. Mir, R.A.; Jain, S.K.; Lohani, A.K.; Ahmed, R.; Ahmad, S.T.; Wani, G.F. Application of Mike 11 for One-Dimensional GLOF Modeling of a Rapidly Expanding Dalung Proglacial Lake, Indus River Basin, Western Himalaya. In Challenges of Disasters in Asia; Springer: Singapore, 2022; pp. 147–161. [Google Scholar]
  136. Murtaza, K.O.; Dar, R.A.; Paul, O.J.; Bhat, N.A.; Romshoo, S.A. Glacial geomorphology and recent glacial recession of the Harmukh Range, NW Himalaya. Quat. Int. 2021, 575, 236–248. [Google Scholar] [CrossRef]
  137. Rashid, I.; Majeed, U.; Aneaus, S.; Pelto, M. Linking the recent glacier retreat and depleting streamflow patterns with land system changes in Kashmir Himalaya, India. Water 2020, 12, 1168. [Google Scholar] [CrossRef] [Green Version]
  138. Romshoo, S.A.; Murtaza, K.O.; Abdullah, T. Towards understanding various influences on mass balance of the Hoksar Glacier in the Upper Indus Basin using observations. Sci. Rep. 2022, 12, 15669. [Google Scholar] [CrossRef]
  139. Pieczonka, T.; Bolch, T.; Kröhnert, M.; Peters, J.; Liu, S. Glacier branch lines and glacier ice thickness estimation for debris-covered glaciers in the Central Tien Shan. J. Glaciol. 2018, 64, 835–849. [Google Scholar] [CrossRef] [Green Version]
  140. Magnin, F.; Haeberli, W.; Linsbauer, A.; Deline, P.; Ravanel, L. Estimating glacier-bed overdeepenings as possible sites of future lakes in the de-glaciating Mont Blanc massif (Western European Alps). Geomorphology 2020, 350, 106913. [Google Scholar] [CrossRef]
  141. Cathala, M.; Magnin, F.; Linsbauer, A.; Haeberli, W. Modelling and characterizing glacier-bed overdeepenings as sites for potential future lakes in the deglaciating French Alps. Geomorphol.-Relief Process. Environ. 2021, 27, 19–36. [Google Scholar] [CrossRef]
  142. Bello, C.; Suarez, W.; Cruz, R.; Rueda, N.F.B. Ice thickness measurements and volume estimates for Znosko glacier (Antarctica) using Glabtop model. Earth Space Sci. Open Arch. ESSOAr 2020. [Google Scholar]
  143. Gantayat, P.; Kulkarni, A.V.; Srinivasan, J. Estimation of ice thickness using surface velocities and slope: Case study at Gangotri Glacier, India. J. Glaciol. 2014, 60, 277–282. [Google Scholar] [CrossRef] [Green Version]
  144. Majeed, U.; Rashid, I.; Sattar, A.; Allen, S.; Stoffel, M.; Nüsser, M.; Schmidt, S. Recession of Gya Glacier and the 2014 glacial lake outburst flood in the Trans-Himalayan region of Ladakh, India. Sci. Total Environ. 2021, 756, 144008. [Google Scholar] [CrossRef]
  145. Jawak, S.D.; Kumar, S.; Luis, A.J.; Pandit, P.H.; Wankhede, S.F.; Anirudh, T.S. Seasonal comparison of velocity of the eastern tributary Glaciers, Amery Ice Shelf, Antarctica, using SAR offset tracking. ISPRS Annals of the Photogrammetry. Remote Sens. Spat. Inf. Sci. 2019, 4, 595–600. [Google Scholar]
  146. Nela, B.R.; Singh, G.; Kulkarni, A.V. Ice thickness distribution of Himalayan glaciers inferred from DInSAR-based glacier surface velocity. Environm.l Monit. Assess. 2023, 195, 1–20. [Google Scholar] [CrossRef]
  147. Tiwari, R.K.; Gupta, R.P.; Arora, M.K. Estimation of surface ice velocity of Chhota-Shigri glacier using sub-pixel ASTER image correlation. Curr. Sci. 2014, 106, 853–859. [Google Scholar]
  148. Sahu, R.; Gupta, R.D. Spatiotemporal variation in surface velocity in Chandra basin glacier between 1999 and 2017 using Landsat-7 and Landsat-8 imagery. Geocarto Int. 2021, 36, 1591–1611. [Google Scholar] [CrossRef]
  149. Singh, D.K.; Thakur, P.K.; Naithani, B.P.; Dhote, P.R. Spatio-temporal analysis of glacier surface velocity in dhauliganga basin using geo-spatial techniques. Environ. Earth Sci. 2021, 80, 11. [Google Scholar] [CrossRef]
  150. Shukla, A.; Garg, P.K. Spatio-temporal trends in the surface ice velocities of the central Himalayan glaciers, India. Glob. Planet. Chang. 2020, 190, 103187. [Google Scholar] [CrossRef]
  151. Patel, L.K.; Sharma, P.; Singh, A.T.; Pratap, B.; Oulkar, S.; Thamban, M. Spatial surface velocity pattern in the glaciers of Chandra Basin, Western Himalaya. Geocarto Int. 2022, 37, 5327–5344. [Google Scholar] [CrossRef]
  152. Zhang, T.; Wang, W.; Shen, Z.; Zhan, N.; Wang, Z.; An, B. Understanding the 2004 glacier detachment in the Amney Machen Mountains, northeastern Tibetan Plateau, via multi-phase modeling. Landslides 2022, 1–16. [Google Scholar] [CrossRef]
  153. Mool, P.; Bajracharya, S.R.; Joshi, S.P. Inventory of Glaciers, Glacier Lakes and Glacial Lake Outburst Floods, Monitoring and Early Warning System in the Hindu Kush-Himalayan Region; International Centre for Integrated Mountain Development: Kathmandu, Nepal, 2001. [Google Scholar]
  154. Bhatt, C.M.; Rao, G.S.; Farooq, M.; Manjusree, P.; Shukla, A.; Sharma, S.V.S.P.; Kulkarni, S.S.; Begum, A.; Bhanumurthy, V.; Diwakar, P.G.; et al. Satellite-based assessment of the catastrophic Jhelum floods of September 2014, Jammu & Kashmir, India. Geomat. Nat. Hazards Risk 2017, 8, 309–327. [Google Scholar]
  155. Malik, I.H.; Najmul Islam Hashmi, S. Ethnographic account of flooding in North-Western Himalayas: A study of Kashmir Valley. GeoJournal 2020, 87, 1265–1283. [Google Scholar] [CrossRef]
  156. Campbell, J.G.; Pradesh, H. Inventory of Glaciers, Glacial Lakes and the Identification of Potential Glacial Lake Outburst Floods (GLOFs) Affected by Global Warming in the Mountains of India, Pakistan and China/Tibet Autonomous Region; International Centre for Integrated Mountain Development: Lalitpur, Nepal, 2005. [Google Scholar]
  157. Aggarwal, A.; Jain, S.K.; Lohani, A.K.; Jain, N. Glacial lake outburst flood risk assessment using combined approaches of remote sensing, GIS and dam break modelling. Geomat. Nat. Hazards Risk 2016, 7, 18–36. [Google Scholar] [CrossRef]
  158. Washakh, R.M.A.; Chen, N.; Wang, T.; Almas, S.; Ahmad, S.R.; Rahman, M. GLOF risk assessment model in the Himalayas: A case study of a hydropower project in the Upper Arun River. Water 2019, 11, 1839. [Google Scholar] [CrossRef] [Green Version]
  159. Iribarren Anacona, P.; Norton, K.P.; Mackintosh, A. Moraine-dammed lake failures in Patagonia and assessment of outburst susceptibility in the Baker Basin. Nat. Hazards Earth Syst. Sci. 2014, 14, 3243–3259. [Google Scholar] [CrossRef] [Green Version]
  160. Maskey, S.; Kayastha, R.B.; Kayastha, R. Glacial lakes outburst floods (GLOFs) modelling of Thulagi and lower Barun glacial lakes of Nepalese Himalaya. Prog. Disaster Sci. 2020, 7, 100106. [Google Scholar] [CrossRef]
  161. Sattar, A.; Goswami, A.; Kulkarni, A.V.; Emmer, A. Lake evolution, hydrodynamic outburst flood modeling and sensitivity analysis in the Central Himalaya: A case study. Water 2020, 12, 237. [Google Scholar] [CrossRef] [Green Version]
  162. Dimri, A.P.; Palazzi, E.; Daloz, A.S. Elevation dependent precipitation and temperature changes over Indian Himalayan region. Clim. Dyn. 2022, 59, 1–21. [Google Scholar] [CrossRef]
  163. Kuniyal, J.C.; Kanwar, N.; Bhoj, A.S.; Rautela, K.S.; Joshi, P.; Kumar, K.; Sofi, M.S.; Bhat, S.U.; Rashid, I.; Lodhi, M.S.; et al. Climate change impacts on glacier-fed and non-glacier-fed ecosystems of the Indian Himalayan Region: People’s perception and adaptive strategies. Curr. Sci. 2021, 120, 888. [Google Scholar] [CrossRef]
  164. Sanjay, J.; Revadekar, J.V.; Ramarao, M.V.S.; Borgaonkar, H.; Sengupta, S.; Kothawale, D.R.; Ratnam, J.V. Temperature changes in India. In Assessment of Climate Change over the Indian Region; Springer: Singapore, 2020; pp. 21–45. [Google Scholar]
  165. Parvaze, S.; Ahmad, L.; Parvaze, S.; Kanth, R.H. Climate change projection in Kashmir Valley (J and K). Curr. World Environ. 2017, 12, 107. [Google Scholar] [CrossRef]
  166. Romshoo, S.A.; Fayaz, M.; Meraj, G.; Bahuguna, I.M. Satellite-observed glacier recession in the Kashmir Himalaya, India, from 1980 to 2018. Environ. Monit. Assess. 2020, 192, 597. [Google Scholar] [CrossRef]
Figure 1. Distribution of GLOF events reported in the India Himalayan region.
Figure 1. Distribution of GLOF events reported in the India Himalayan region.
Remotesensing 14 05957 g001
Figure 2. Location map of the study area.
Figure 2. Location map of the study area.
Remotesensing 14 05957 g002
Figure 3. Annual temperature (Tmean, Tmin, and Tmax) and precipitation at Shalimar station located in Srinagar (SKUAST-Met station).
Figure 3. Annual temperature (Tmean, Tmin, and Tmax) and precipitation at Shalimar station located in Srinagar (SKUAST-Met station).
Remotesensing 14 05957 g003
Figure 4. Overall methodology of this study.
Figure 4. Overall methodology of this study.
Remotesensing 14 05957 g004
Figure 5. Schematic representation of the working principle of the GlabTop-2 model. Glacier polygons are transformed into a raster matching the DEM cells.
Figure 5. Schematic representation of the working principle of the GlabTop-2 model. Glacier polygons are transformed into a raster matching the DEM cells.
Remotesensing 14 05957 g005
Figure 6. (a) Cross sections at various locations across the river channel; (b) XYZ perspective plot of the river channel; and (c) river bed profile depicting distance with respect to the elevation of the river channel. * XS interpolations.
Figure 6. (a) Cross sections at various locations across the river channel; (b) XYZ perspective plot of the river channel; and (c) river bed profile depicting distance with respect to the elevation of the river channel. * XS interpolations.
Remotesensing 14 05957 g006
Figure 7. Gangabal lake expansion from 1972 to 2020.
Figure 7. Gangabal lake expansion from 1972 to 2020.
Remotesensing 14 05957 g007
Figure 8. Glacier area changes from 1990 to 2020.
Figure 8. Glacier area changes from 1990 to 2020.
Remotesensing 14 05957 g008
Figure 9. Spatial distribution of elevation changes over Harmukh glacier from 2000 to 2012.
Figure 9. Spatial distribution of elevation changes over Harmukh glacier from 2000 to 2012.
Remotesensing 14 05957 g009
Figure 10. Glacier surface velocity/displacement with vectors showing flow direction. The red and yellow squares depict the maximum velocity near the crevasses.
Figure 10. Glacier surface velocity/displacement with vectors showing flow direction. The red and yellow squares depict the maximum velocity near the crevasses.
Remotesensing 14 05957 g010
Figure 11. (a) Slope map showing potential zones of mass movement between the glacier and lake. (b) Potential zones of mass movement on Google Earth imagery.
Figure 11. (a) Slope map showing potential zones of mass movement between the glacier and lake. (b) Potential zones of mass movement on Google Earth imagery.
Remotesensing 14 05957 g011
Figure 12. Moraine breach hydrographs showing peak discharge with varying beach widths (20, 40, 60, 80).
Figure 12. Moraine breach hydrographs showing peak discharge with varying beach widths (20, 40, 60, 80).
Remotesensing 14 05957 g012
Figure 13. Breach hydrograph in scenario 2 releasing 75% volume of the water from Gangabal lake.
Figure 13. Breach hydrograph in scenario 2 releasing 75% volume of the water from Gangabal lake.
Remotesensing 14 05957 g013
Figure 14. Plots showing (a) max. flood depth (m), (b) max. flood velocity (m/s), (c) WSE (m).
Figure 14. Plots showing (a) max. flood depth (m), (b) max. flood velocity (m/s), (c) WSE (m).
Remotesensing 14 05957 g014
Figure 15. Spatial distribution of flood inundation from the lake outlet to the last selected location: (a) zoomed view of inundation at Khanam bridge; (b) routed channel confluence with Sind river.
Figure 15. Spatial distribution of flood inundation from the lake outlet to the last selected location: (a) zoomed view of inundation at Khanam bridge; (b) routed channel confluence with Sind river.
Remotesensing 14 05957 g015
Figure 16. (a) GLOF propagation at various locations downstream; (b) zoomed view of the inundated settlement; (c) flood depth at Settlement 4.
Figure 16. (a) GLOF propagation at various locations downstream; (b) zoomed view of the inundated settlement; (c) flood depth at Settlement 4.
Remotesensing 14 05957 g016
Figure 17. Routed hydrographs at six important locations along the flow channel: (a) Naranag temple; (b) Khanan bridge; (c) settlement 1; (d) confluence; (e) bridge; and (f) Drag Tanga village.
Figure 17. Routed hydrographs at six important locations along the flow channel: (a) Naranag temple; (b) Khanan bridge; (c) settlement 1; (d) confluence; (e) bridge; and (f) Drag Tanga village.
Remotesensing 14 05957 g017
Figure 18. GLOF inundation modeled in scenario 2: (a) inundated settlement and bridge; (b) decrease in the inundated area at unnamed settlement in scenario 2; (c) flood velocity at Naranaag temple.
Figure 18. GLOF inundation modeled in scenario 2: (a) inundated settlement and bridge; (b) decrease in the inundated area at unnamed settlement in scenario 2; (c) flood velocity at Naranaag temple.
Remotesensing 14 05957 g018
Figure 19. LULC maps generated for GLOFs: (a) scenario 1 and (b) scenario 2.
Figure 19. LULC maps generated for GLOFs: (a) scenario 1 and (b) scenario 2.
Remotesensing 14 05957 g019
Table 1. Data sets used in the study.
Table 1. Data sets used in the study.
S. No.DataYearNo. of BandsSpatial ResolutionSourcePurpose
1Planet CubeSat202083Planet
https://www.planet.com/
(accessed on 2 August 2021)
Identification of lake parameters
2Resourcesat-2202035.8NRSC Bhuvan portal
https://bhuvan.nrsc.gov.in/home/index.php (accessed on 25 October 2021)
Glacier and lake mapping
3Sentinel-2A20201310USGS Earth Explorer
http://www.earthexplorer.usgs.gov/
(accessed on 12 July 2021)
Glacier and lake mapping
4Landsat 5 TM1972630USGS Earth Explorer
http://www.earthexplorer.usgs.gov/
(accessed on 12 July 2021)
Glacier and lake mapping
Landsat 7 ETM19809
Landsat 7 ETM+1990, 20009
Landsat 8 OLI20109
5Landsat 8 Pan bands2019,
2020
115USGS Earth Explorer
http://www.earthexplorer.usgs.gov/
(accessed on 12 July 2021)
Glacier velocity estimation
6Goole Earth images1980 to 2020-0.5–30Google Earth
https://earth.google.com/web/
(accessed on 02 November 2021)
Cross check
7ALOS-PALSAR DEM2015112.5USGS Earth Explorer
http://www.earthexplorer.usgs.gov/
(accessed on 02 November 2021)
Hydrodynamic modeling
8SRTM-DEM2000190USGS Earth Explorer
http://www.earthexplorer.usgs.gov/
(accessed on 02 November 2021)
Geodetic mass balance
9TAN-DEM X2012190GEOC German Observation Centre
https://download.geoservice.dlr.de/TDM90/
(accessed on 02 November 2021)
Geodetic mass balance
Table 2. The empirical equations derived from the literature: (a) depth, (b) volume, and (c) peak discharge.
Table 2. The empirical equations derived from the literature: (a) depth, (b) volume, and (c) peak discharge.
S. No.Empirical EquationSource
Depth
1D = 0.104 × A0.42[108]
2Dm = 55 × A0.25[37]
3Dm = 4 × 10−5 × A + 5.0564[111]
Volume
1V = 0.035 × A1.5[113]
2V = 0.104 × A1.42[108]
3Dm = 55 × A0.25
V = A × Dm
[37]
4Dm = 4 × 10−5 × A + 5.0564
V = A × Dm
[111]
5V = 0.0522 × A1.1766Sharma et al. 2018
6V = 0.0578 × A1.4683[112]
7V = 40.67 × A1.184 3.218 × Ratio (mxw/mxl)[117]
Peak Discharge
1Qmax = 0.72 × V0.53[113]
2Qmax = 0.0048 × V0.896[116]
3Qmax = 0.00077 × V1.017[108]
Table 3. Breach parameters calculated from Froehlich’s equation.
Table 3. Breach parameters calculated from Froehlich’s equation.
GLOF ScenariosReleased Volume (106 m3)Percentage Volume of the Lake Released (%)Breach Width
(m)
Breach Depth (m)Breach Formation
Time (min)
Peak Discharge (m3)
Scenario 173.310080704017,787.89
Scenario 252.4754050309155.36
Table 4. Thickness classwise area in Harmukh glacier.
Table 4. Thickness classwise area in Harmukh glacier.
ValueThickness Classes (m)Pixel CountArea (km2)%
10–3090631.3647.18
231–6080501.2141.85
3>6021060.3210.95
Total19,2162.9100
Table 6. Flood wave parameters at various important locations.
Table 6. Flood wave parameters at various important locations.
Distance from Lake (km)Location NameLat.Long.Cross SectionVolume (cu. m)Channel Elevation (m)Water Surface Elevation (m)Flow Depth (m)Flow Velocity (m/s)Peak Discharge (Qmax)Peak Discharge Arrival Time
11.67Narang Temple34.353374.97745,50089,647.662229.282243.7114.438.8313,995.3601:00
11.9Narang Kangan Temple34.351974.975545,230.589,671.382220.222234.213.989.1913,513.1501:00
13.72Bridge 134.33974.963543,50090,308.612143.742160.9617.228.1913,351.2701:00
14.67Bridge 234.331374.960242,266.591,021.122093.892103.519.6211.1712,987.0301:00
16.63Khanan village34.319474.944840,25091,389.12012.22022.7410.5410.0512,896.5201:10
17.25Dam34.315874.939839,50091,589.311988.582002.313.729.5112,205.8401:10
17.67Bridge and Settlement34.31374.935939,25091,712.171982.081995.3113.238.4612,103.5101:10
18.69Babanagre Khanan Bridge34.306374.927637,94092,335.891932.621949.9917.378.9911,834.6401:10
19.54Settlement34.302974.922637,25092,504.231909.121923.9614.8412.4611,708.6601:10
21.78Settlement34.288774.904434,763.593,416.971833.581850.4916.9111.0311,485.6901:20
23.23Settlement34.280874.893533,25094,336.231788.051799.5211.476.8711,320.7501:20
23.94Road and bridge34.275974.888632,50093,936.861767.611768.40.790.5711,270.0101:20
24.23Road and bridge34.274874.886132,27393,951.161763.491767.343.850.9411,127.3701:20
24.76Tail race park34.271274.882231,77394,122.861753.241762.459.215.8610,888.8301:20
26.11Bridge34.276474.868930,50096,439.31734.21742.468.263.810,848.6801:30
26.62Drag Tanga village34.283274.858428,00097,429.231699.91710.7110.814.2410,223.9801:40
Table 7. The inundated area under different LULC classes in scenario 1 and scenario 2.
Table 7. The inundated area under different LULC classes in scenario 1 and scenario 2.
S. NoScenario 1Scenario 2
Lulc ClassArea km2Percentage (%)Lulc ClassArea km2Percentage (%)
1River course0.815.09River course0.5816.76301
2Forestland0.713.21Forestland0.5114.73988
3Grass0.47.55Grass0.4312.42775
4Cropland0.611.32Cropland0.3510.11561
5Scrub/shrub1.120.75Scrub/shrub0.6217.91908
6Built-up area0.35.66Built-up area0.195.491329
7Barren 1.426.42Barren 0.7822.54335
Total5.3100Total3.46100
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ahmed, R.; Rawat, M.; Wani, G.F.; Ahmad, S.T.; Ahmed, P.; Jain, S.K.; Meraj, G.; Mir, R.A.; Rather, A.F.; Farooq, M. Glacial Lake Outburst Flood Hazard and Risk Assessment of Gangabal Lake in the Upper Jhelum Basin of Kashmir Himalaya Using Geospatial Technology and Hydrodynamic Modeling. Remote Sens. 2022, 14, 5957. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14235957

AMA Style

Ahmed R, Rawat M, Wani GF, Ahmad ST, Ahmed P, Jain SK, Meraj G, Mir RA, Rather AF, Farooq M. Glacial Lake Outburst Flood Hazard and Risk Assessment of Gangabal Lake in the Upper Jhelum Basin of Kashmir Himalaya Using Geospatial Technology and Hydrodynamic Modeling. Remote Sensing. 2022; 14(23):5957. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14235957

Chicago/Turabian Style

Ahmed, Rayees, Manish Rawat, Gowhar Farooq Wani, Syed Towseef Ahmad, Pervez Ahmed, Sanjay Kumar Jain, Gowhar Meraj, Riyaz Ahmad Mir, Abid Farooq Rather, and Majid Farooq. 2022. "Glacial Lake Outburst Flood Hazard and Risk Assessment of Gangabal Lake in the Upper Jhelum Basin of Kashmir Himalaya Using Geospatial Technology and Hydrodynamic Modeling" Remote Sensing 14, no. 23: 5957. https://0-doi-org.brum.beds.ac.uk/10.3390/rs14235957

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

Article Metrics

Back to TopTop