Reconstruction and Characterisation of Past and the Most Recent Slope Failure Events at the 2021 Rock-Ice Avalanche Site in Chamoli, Indian Himalaya

Frequent ice avalanche events are being reported across the globe in recent years. On the 7 February 2021, a flash flood triggered by a rock-ice avalanche with an unusually long runout distance, caused significant damage of life and property in the Tapovan region of the Indian Himalaya. Using multi-temporal satellite datasets, digital terrain models (DTMs) and simulations, here we report the pre-event and during-event flow characteristics of two large-scale avalanches within a 5-year interval at the slope failure site. Prior to both the events, we observed short-term and longterm changes in surface velocity (SV) with maximum SVs increasing up to over 5 times the normal values. We further simulated the events to understand their mechanical characteristics leading to long runouts. In addition to its massive volume, the extraordinary magnitude of the 2021 event can partly be attributed to the possible remobilisation and entrainment of the colluvial deposits from previous ice and snow avalanches. The anomalous SVs should be explored further for their suitability as a possible remotely observable precursor of ice avalanches from hanging glaciers. This sequence of events highlights that there is a need to take into account the antecedent conditions, while making a holistic assessment of the hazard.


Introduction
Slope failure events in high-mountains globally hamper infrastructure and economic development, and cause hundreds of fatalities each year [1][2][3]. Contemporary research [4][5][6][7] linked the increasing mass movement hazards in high mountains to the increase in temperatures and erratic precipitation conditions. However, many of these studies also attribute the cascading effects of rapidly changing periglacial landscape and local topography for these mountain disasters [3][4][5][6]. Therefore, it is relevant to have a basic understanding of various slope failure processes in high mountains [8]. Mass movements in high mountains can be characterised based on three main parameters [9]: (1) constituents involved (regolith, rock, glacial ice, and snow), (2) propagation velocity (slow, intermediate, and fast), and (3) characteristics of the movement (debris cloud, slurry, and coherent mass). An understanding of these characterising parameters further allows for a distinct classification of high-mountain slope failure events within categories, such as mudflow, debris flow, avalanche, landslide, rockfall, and debris fall. Falls are abrupt movements of rocks or debris after detaching from steep slopes along structural discontinuities, such as joints, faults, bedding planes, and other fractures, and movement can occur in form of free-fall, bouncing, and rolling. A mudflow is a very rapid movement of saturated, fine-grained deposits with water contents equal to or higher than the liquid limit, changing the phase of the deposits from solid to plastic phase [10]. Debris flows are also very rapid movements temporal satellite datasets and robust rock-ice avalanche modelling approaches. We focus on quantifying and discussing the mechanical characteristics (momentum, pressure, and shear stress) of this mass movement leading to its translation into an enormous debris flow. We selected two massive avalanche events, temporally spaced by a~4.5-year gap, which ran through the same valley stretch, but were constituently different. The first event of September 2016, hereafter referred to as "the first slope failure (SF1)", was a frontal ice-block failure of a hanging glacier. The second event of the 7 February 2021, hereafter referred to as "the second slope failure (SF2)", was a rock-ice avalanche and was over~2.5 times more voluminous than the 2016 ice avalanche. The successive occurrence of these two major events in the same topographical setting provides us an excellent opportunity to understand the difference in their mass movement behaviour owing to varying degree of ice content within them. The manifestation of these two massive slope failure events, SF1 and SF2, originating from nearly the same release height and hitting the same valley within a period of 5 years, is unique and offers us an unprecedented natural testbed to understand how frequent mass movements in the coming years can vary in terms of their runout and destructive nature. A thermomechanical simulation of SF1 offers us the ability to modify the valley topography to include deposits from SF1 in analysing their possible influence on the debris flow runout of SF2 [25]. In the following sections, we provide a brief overview of the avalanche sites, and separately present our findings on the estimated surface velocity (SV) and the modelled accounts of the two slope failures. Considering that SF1 was a frontal ice-block failure, while SF2 was a landslide causing ice-bedrock failure, in this paper, we used SV as an all-inclusive term to account for all the movements (sliding and deformational) occurring at the avalanche sites. Thus, this paper investigates three main research questions: (1) are there any quantifiable pre-event changes in SVs? (2) how do SV changes vary between two constituently different slope failures?, and (3) what are the differences in mechanical characteristics of two large slope failures originating from the same release height and hitting the same valley?

Study Area
The study area is situated in Central Himalaya within a transitional zone between the wetter Southern Himalayan ranges dominated by Indian Summer Monsoons (ISMs) and the cold arid Tibetan Highlands in the north [26]. Although the predominant precipitation source is ISM during the warmer months, the region also receives considerable precipitation from mid-latitude westerlies in the form of winter snowfall [27] between December and March. The average annual precipitation is~1100 mm with sub-zero winter temperatures and summer temperature oscillating between 5 • C and 29 • C [28]. These weather systems display a significant control over regional cryosphere, ecosystems, and landscape evolution, ultimately influencing the socioeconomics of the region [26]. The rugged topography of this region lies between~2000 and~7800 m asl [29]. Topographically, the region is a hotbed for high-velocity mass movement events as it is characterised by high peaks, narrow river valleys, and deep gorges [26]. The study area partly falls within the Alaknanda River catchment and is abundant in several mountain streams and tributaries fed by meltwaters [26].
Geologically, the study area is located between the Vaikrita Thrust to the south and the Malari Fault to the North [26], and consists of the Vaikrita Group of high-grade metamorphic rocks with both acidic and basic intrusions, which are well-exposed along the Rishiganga River valley [30]. These metamorphic rocks are dominated by schist, amphibolite, and gneisses [31], and such a stratified geology is considered prone to episodes of glaciation, deglaciation, and neotectonics [32,33]. Thus, the geology of the study area is structurally fragile and exposed to the contrasting regional climatic conditions [26]. The elevation range of the detached SF1 ice slab was~5200-5400 m asl while that of the avalanched SF2 rock-ice body was~4900-5600 m asl. Interestingly, the equilibrium line altitude (ELA), i.e., the dividing altitude between the accumulation and ablation zones of a glacier, reported for glaciers in this region is at~5400 m asl [28]. At ELA, the annual snow Remote Sens. 2022, 14, 949 4 of 32 accumulation equals the annual snow ablation over a period of several years, making this altitude very sensitive to climate fluctuations [28]. The fact that both these events occurred at the slopes close to the reported ELA, indicates their particular sensitivity towards both winter accumulation and summer ablation. Further details on the geomorphic and dimensional characteristics of these slope failures are provided in the Results and Discussion section for better continuity and understanding of the events. The avalanched debris hit the valley bottom of Ronti Gad stream, which merges in the Rishiganga River before Raini Village. Rishiganga contributes to the Dhauliganga River after crossing Raini Village. SF2 followed this route for~32 km as a flash flood [29], causing a widespread destruction of life and property and visibly changing the geomorphology of the valley. A short YouTube video [34] captured the enormity and the destructive nature of the flash flood triggered by the SF2 event up to~8 km downstream from the deposition site. Figure 1 provides a detailed overview of the study area.
Remote Sens. 2022, 14, x FOR PEER REVIEW 4 of 33 lanched SF2 rock-ice body was ~4900-5600 m asl. Interestingly, the equilibrium line altitude (ELA), i.e., the dividing altitude between the accumulation and ablation zones of a glacier, reported for glaciers in this region is at ~5400 m asl [28]. At ELA, the annual snow accumulation equals the annual snow ablation over a period of several years, making this altitude very sensitive to climate fluctuations [28]. The fact that both these events occurred at the slopes close to the reported ELA, indicates their particular sensitivity towards both winter accumulation and summer ablation. Further details on the geomorphic and dimensional characteristics of these slope failures are provided in the Results and Discussion section for better continuity and understanding of the events. The avalanched debris hit the valley bottom of Ronti Gad stream, which merges in the Rishiganga River before Raini Village. Rishiganga contributes to the Dhauliganga River after crossing Raini Village. SF2 followed this route for ~32 km as a flash flood [29], causing a widespread destruction of life and property and visibly changing the geomorphology of the valley. A short YouTube video [34] captured the enormity and the destructive nature of the flash flood triggered by the SF2 event up to ~8 km downstream from the deposition site. Figure 1 provides a detailed overview of the study area.

Materials and Methods
In addition to the high-resolution images obtained from Google Earth for visual analysis, we mainly used two remote sensing datasets for performing quantitative analyses: (1) multi-temporal satellite images for SV estimation and (2) high-resolution multi-temporal digital terrain models (DTMs) for Rapid Mass Movement Simulation (RAMMS) modelling. The details of these datasets are provided in Table 1. The image IDs of the satellite images used for visual analysis are provided in the respective figure captions. The overall workflow, designed to investigate the pre-event and during-event flow characteristics, is depicted in Figure 2.

Surface Velocity (SV) Change Estimation
We used the Co-Registration of Optically Sensed Images and Correlation (COSI-Corr) tool [41], which works as a plug-in module to the Environment for Visualizing Images (ENVI) software for estimating horizontal displacement fields. COSI-Corr implements a phase correlation method based on the Fourier Shift Theorem and the retrieved relative displacement between an image pair is a resultant of an iterative, unbiased process, estimating the phase difference of their Fourier transform [41]. The standard outputs are in the form of East-West (EW) displacement, North-South (NS) displacement, and signal-to-noise ratio (SNR), and the absolute displacement is calculated as the square root of the sum of squares of the EW and NS displacement components. This displacement can be divided by various units of time, i.e., day in the present study, to obtain the SVs. The typical correlation accuracy achieved at sub-pixel resolution by COSI-Corr is as high as 1/20-1/10 of the spatial resolution of an image pixel [42]. ysis, we mainly used two remote sensing datasets for performing quantitative analyses: (1) multi-temporal satellite images for SV estimation and (2) high-resolution multi-temporal digital terrain models (DTMs) for Rapid Mass Movement Simulation (RAMMS) modelling. The details of these datasets are provided in Table 1. The image IDs of the satellite images used for visual analysis are provided in the respective figure captions. The overall workflow, designed to investigate the pre-event and during-event flow characteristics, is depicted in Figure 2.  The panchromatic band (Band 8) of Landsat 8 Operational Land Imager (OLI) and blue band (B2) of Sentinel-2 MultiSpectral Instrument (MSI) are the preferred bands for obtaining glacier and mass movement SVs [43,44]. Landsat 8 images were downloaded from the United States Geological Survey's (USGS's) EarthExplorer portal [45] and Sentinel-2 images were downloaded from the Copernicus Open Access Hub [46]. A 15 m/pixel resolution Landsat 8 Panchromatic band was selected for the SV estimation of SF1, as appropriate Sentinel-2 data was not available during September 2015. We found cloud-free Landsat 8 images covering 1-17 September 2015 and 3-19 of September 2016, with identical surface conditions for comparing the same duration SVs in the two consecutive years prior to SF1. We did not go back farther than one year before SF1 as it was primarily a frontal block failure and the critical ice load conditions for such events vary significantly from one year to another in this region [47]. The selected image pairs clearly highlighted the effect of this aspect on the estimated SVs. In fact, SF1 is a recurring frontal ice-block failure and we identified a previous similar event in early 2000 ( Figure 3). In contrast, SF2 being primarily a landslide causing an ice-bedrock failure, we opted to perform yearly SV observations for it, starting from 2016 with good temporal availability of Sentinel-2 images. We again ensured cloud-free and similar surface conditions by selecting all the images from October, and this approach significantly reduced uncertainties. For SV observations~40 days prior to SF2, we again opted for Band 2 of Sentinel-2. The presence of seasonal snow was unavoidable in this case. In order to effectively deal with mountain shadow, which was covering the entire SF2 hanging glacier uniformly during these~40 days, we first clipped and then contraststretched the SF2 site in all the scenes by comparative histogram matching. Afterwards, we performed COSI-Corr processing only for this contrast-stretched part and the obtained results were considerably good with acceptable uncertainties ( Table 2).
vations ~40 days prior to SF2, we again opted for Band 2 of Sentinel-2. The presence of seasonal snow was unavoidable in this case. In order to effectively deal with mountain shadow, which was covering the entire SF2 hanging glacier uniformly during these ~40 days, we first clipped and then contrast-stretched the SF2 site in all the scenes by comparative histogram matching. Afterwards, we performed COSI-Corr processing only for this contrast-stretched part and the obtained results were considerably good with acceptable uncertainties (Table 2).    Lacroix et al. [43] suggested to keep the window size for COSI-Corr processing as a compromise between the robustness of large windows and the ability of smaller windows to detect motions of rather small objects. Thus, we tried several window sizes and the most effective window size for our study site with the least uncertainties for both Landsat 8 and Sentinel-2 was found to be 32 × 32 window. We used this window for further processing all the data presented in this study. In the post-processing stage, we applied two of the filters effectively used by previous studies (e.g., [48]) to make the results more robust.
Using the SNR filter, pixels with low SNR values (SNR < 0.90) were masked as poorly correlated pixels. Using the magnitude filter, we masked the abrupt values, assuming that the movement may not change abruptly but gradually [42], thus obtaining spatially smooth SV maps. However, we intentionally did not apply the magnitude filter for observing yearly movements as we wanted to identify the localised zones showing anomalous movements on a yearly scale. This approach perfectly served the purpose, as we could see the gradual spread of the stress field along the indentation for SF2 hanging glacier with each passing year until 2019.

Surface Velocity Uncertainty Estimation
The uncertainty introduced due to possible orthorectification errors may result in some minor sub-pixel-scale horizontal shifts in the Landsat 8 and Sentinel-2 images. However, in case of imageries for the same path and row used within this study, the orthorectification uncertainties affecting SV estimations are negligible [49]. Eventually, for estimating the cumulative uncertainties in our SV estimates, we adopted the approach suggested by previous studies (e.g., [48,50]). We sampled all the pixels representing the stable terrain devoid of ice and vegetation. For these stable pixels, the overall mean and standard deviation of displacements were summed up and divided by the temporal separation of the correlated images to obtain the overall uncertainty per correlated scene-pair. The estimated uncertainties are provided in Table 2 and were of the same order of magnitude as reported by the previous studies (e.g., [48,50]), except for the last pair of 26 January 2021-5 February 2021, owing to extreme snow cover in 5 February 2021 image.

Rapid Mass Movement Simulation (RAMMS)
The RAMMS tool for mass movement modelling, developed at the Swiss Institute for Snow and Avalanche Research, was successfully implemented for high mountain regions, such as the Alps (e.g., [51][52][53]), Andes (e.g., [54]), Himalaya (e.g., [55]), Southern Alps (e.g., [56]), Tatras (e.g., [57]), and others for realistic thermomechanical simulations. The RAMMS tool [13,25,58] has a user-friendly interface that allows for the integration of state-of-the-art numerical solution methods with input parameters derived from DTMs and satellite images or maps. An efficient second-order numerical solution of the depthaveraged avalanche dynamics equations forms the core of this programme [25]. More specifically, the model relies on the two-parameter Voellmy model, which is actually field calibrated using data from a real scale avalanche test site in Vallée de la Sionne, canton Valais, Switzerland [25]. The usual input parameters consist of topography, release area, release height, forest extent, and calibrated frictional coefficients, while the modelled output includes all types of mechanical parameters, such as avalanche velocity, momentum, flux, pressure, runout, and depositional height. Generating input data for this programme is straightforward as all the input parameters can easily be derived from the DTM and satellite images in the form of geographical information system (GIS) compatible shape files. The software manual provides sufficient details on how the tool works and how various frictional coefficients and input parameters can be manipulated to obtain the most realistic simulations. Detailed information on RAMMS and its input and output parameters can be downloaded from the software manual webpage [59]. More importantly, not long ago, a rock/ice avalanche model capable of modelling both extreme cases, i.e., avalanches containing mainly ice (e.g., SF1 in the present study) as well as avalanches that mainly consist of rock debris (e.g., SF2 in the present study), was embedded in the RAMMS software and was already successfully used by several case studies worldwide [60].
The hazard potential of many of the recent rock/ice avalanches, including SF2 in Chamoli, has been complex to comprehend, due to the extreme flow transitions among the constituents, i.e., ice, snow, and debris, and the temperature of the material. The proportion of the constituents (ice, snow, and debris) usually decides whether the dry fluidised (powder avalanches) flow regime would prevail, or a fluid-lubricated (debris flows) flow regime. The thermomechanical properties of the material entrained by the avalanche often determine its flow regime and modelling rock/ice avalanche speed and runout, and therefore requires a thermomechanical model accounting for the melting and lubricating properties of the entrained material, which is usually a mixture of snow, ice, water, and debris [13,58]. RAMMS is capable of modelling the thermomechanical regime of various types of avalanches, and therefore we used it as a tool to model during-event flow characteristics of SF1 and SF2.
High-resolution DTMs are extremely useful for geomorphology and mountain hazard research [61,62] and are most preferred for RAMMS simulations, as they inform the model about all the obstructions and natural dams to employ the most appropriate frictional parameters for a particular event. Fortunately, we have the availability of both pre-event and post-event high-resolution (2 m/pixel) DTMs covering SF1 and SF2 (Table 1). These DTMs were used in a study by Shugar et al. [21] and were derived by Bhushan and Shean [36] and Shean et al. [40] from the stereopairs of WorldView-1, WorldView-2, WorldView-3, and GeoEye-1 satellite images. Table 1 provides all the related references detailing the generation and characteristics of these DTMs. These DTMs allowed us to precisely demarcate and estimate the release area and volumes of SF1 and SF2 through DTM differencing as input to RAMMS. These areal and volumetric estimates are detailed in the Results and Discussion Sections. Another model input has to be the avalanche return period, which is pre-defined in the model as 10 years, 30 years, 100 years, and 300 years, representing various known or hypothetical recurrences. For SF1, we examined historical data and found it to be a recurring event, once the hanging glacier reaches its critical geometry. We could confirm that nearly 16 years ago, an event of same areal scale occurred at the SF1 site ( Figure 3). Therefore, an average return period of 10 years was selected for SF1. This also translated in the simulation results in which the simulated SF1 avalanche followed the exact depositional boundaries that were observed in the after-event satellite images. We further ran the model using 30 years as the return period and the results were not equally comparable to the after-event satellite images that we derived using 10 years as the return period. On the other hand, SF2 predominantly consisted of rock debris [21] and was more suited for simulation within debris flow module of RAMMS, without a need to define the return period. Considering that both SF1 and SF2 were over 10 × 10 6 m 3 , the avalanche volume category within RAMMS was set to "Large". The model also takes into account the forest cover present along the avalanche path as an important input parameter, since vegetation provides resistance to the avalanche front. We digitised the forest mask using satellite images acquired prior to both the events. SF1 was mostly uninterrupted by any forest cover, while SF2, being a larger event, encountered dense vegetation along its path in lower reaches of the valley. Considering that SF1 was majorly an ice avalanche, for the avalanche core at rest, we used a density of 850 kg m −3 as recommended by Margreth et al. [53]. For SF2, considering the suggested rock:ice proportion [21] of 80:20 and the predominance of metamorphic rocks at the avalanche site [23], we inputted a density of~2050 kg m −3 [63]. As the final vital input to RAMMS, the Coulomb friction (µ) and the turbulent friction (ε) are defined. The friction values µ and ε strongly depend on the return period and avalanche volume, and, in our study, owing to high-resolution DTMs and temporal image analysis, we were certain of these input parameters. The friction values, µ and ε, depend on the configurational energy content of the avalanche core [13], and in the avalanche module used for SF1 simulation in RAMMS, if a calculation is obtained with constant friction values, no terrain undulations and forest areas are considered. Therefore, for SF1, we opted for the variable friction values, which are automatically estimated by RAMMS, based on DTM-derived parameters (slope angle, altitude, and curvature), and we further provided forest information and the global parameters of the return period and avalanche volume. In case of SF2, the debris flow module of RAMMS needed to be initialised through a carefully selected combination of µ and ε, calibrated using the documented [21][22][23][24] timeand flow-related details of the event. This calibration was achieved through the observed runout extent on the post-disaster satellite imagery and the reported debris flow zones, flow heights, and velocities [21][22][23][24]. Thus, the model was precisely calibrated through extensive iterations to define the optimal friction values of µ and ε as 0.13 and 200 ms −2 , respectively. These values correspond well with the suggested range of values in the RAMMS manual.
A general topic of previous investigations for SF2 has been to explain the possible source of entrainment material, which led to such an enormous debris flow. A recent study [29], based on field observations, suggests that the vast deposits of the 2016 event might have been a possible source of material, which significantly enlarged the volume and extent of SF2. To investigate this aspect, we further used the erosion module in RAMMS while simulating SF2. We provided variable depths of sediment erosions between 2-8 m, based on the simulations by Jiang et al. [22], in the upper reaches of Ronti Gad, in order to model the increase in volume and velocity of SF2 as it travelled along the channel. This erosion corresponds to the material that aids the entrainment of the avalanche front. Supplementary Materials Files (Video S1: Simulated SF1 event, Video S2: Simulated SF2 event) of this paper provide a 3D visualisation of the RAMMS modelling results.

Results and Discussion
We present and discuss the main results of this research under the following subheadings, focussing on the pre-event and during-event flow characteristics of SF1, the long-term and short-term pre-event flow characteristics of SF2, and the simulated during-event runout characteristics of SF2.

SF1: Key Observations and Pre-Event and During-Event Flow Characteristics
SF1 site (Figures 1 and 3) appears to be a hotspot for ice avalanching from the same initiation zone. The maximum cross-sectional lengths and area of the broken part of the hanging glacier during SF1 were~0.57 km (across the flow),~0.92 km (along the flow), and 0.28 km 2 , respectively. SF1 was a frontal block failure and the multi-temporal DTMs of 2021 and 2015 provided by Shugar et al. [21] helped us to estimate a volume of~10 × 10 6 m 3 , corresponding predominantly to the glacier ice and possibly some seasonal snow. Using multi-temporal satellite images, we observed at least two major and comparable slope failure events at the SF1 site during the past two decades. These events were spaced 16.5 years apart and they produced exactly the same indentation (Figure 3), marked by the dashed red curve in Figure 1c. For the earlier event, the 3D perspective views in Figure 3a depict the avalanched zone, and by using all the available Landsat images, we could narrowdown the timing of this event between January and April 2000. Similarly, the 3D perspective views in Figure 3b highlight the repetition of another avalanche (called SF1 in this paper) at the same indentation site. Using freely accessible Sentinel-2 and PlanetScope [64] images, we could narrow-down the timing of SF1 between 19 and 24 September 2016.
To understand the enormity of SF2, there is a need to consider these events from 2000 and 2016 at the SF1 site. These events were large-volume (~10 × 10 6 m 3 ) ice avalanche events leading to large deposits of ice and debris (orange rectangles in Figure 3) in the avalanche deposition zone in a narrow mountain stream named Ronti Gad. This deposition, following the 2000 and 2016 events, covered a~3.5 km length of the stream (Figures 3 and 4). The thick ice-debris deposits from these events tend to persist for multiple years after the event. For example, some of the possible remnants of the deposits from the 2000 event could be observed even in 2005, more than 5 years after the event (Figure 5b).
Another observation made on multiple years images ( Figure 5) is the significantly seasonally active nature of the avalanche slope (Figure 5h). In addition to the ice avalanches, these slopes are always susceptible to snow avalanches during the winter months ( Figure 5). A~500 m-long stretch of the valley, directly beneath the avalanche slopes, is mostly covered by compacted snow, ice, and debris ( Figure 5). These deposits also often hinder the full flow of the narrow stream, as is evident in Figure 4b, building up the temporary zones of increased hydraulic gradient and developing the dammed water pockets along the stream. We identified one such seasonally active avalanche slope below SF1 and SF2 sites ( Figures 5 and 6). While in the 3 November 2013 image of Figure 5 we observe a clear stream devoid of any avalanche deposit from the previous event of 2000 or from subsequent snow avalanches, in the 3 June 2014 and the 7 May 2015 images, we see thick fresh avalanche deposits covering the stream valley. Figure 5h highlights the seasonally active avalanche slope below the SF1 and SF2 sites, which accumulates significant snow mass during winter, often leading to slope failures relatively smaller than SF1 and SF2, but still large enough to form a deposition zone of~500 m to 1 km in length. This large volume of snow, ice, and debris can accentuate any subsequent slope failure events, such as SF2, which occurred in the winter of 2021. We observed the available images of these slopes, days before SF2. Figure 6 highlights the significant snow accumulation on the seasonally active slope on 2 February 2021, while on 5 February 2021 Sentinel-2 image, it could be seen that a fresh cumulative snowfall of~112 mm [22] significantly added to the snow cover. Figure 6 additionally shows the developed indentation from where the hanging glacier broke on 7 February 2021, causing SF2. Terrain profiles for SF1 and SF2 slopes (Figure 7), representing the terrain before either of the reported events occurred, estimate the maximum slope reaching as high as~37 • . The reported average slope for the entire SF2 mass [24] was even higher at~40 • . These slope values are entirely within the most avalanche-prone slope range reported from this region [18].   Another observation made on multiple years images ( Figure 5) is the significantly seasonally active nature of the avalanche slope (Figure 5h). In addition to the ice avalanches, these slopes are always susceptible to snow avalanches during the winter months ( Figure 5). A ~500 m-long stretch of the valley, directly beneath the avalanche slopes, is mostly covered by compacted snow, ice, and debris ( Figure 5). These deposits also often hinder the full flow of the narrow stream, as is evident in Figure 4b, building up the temporary zones of increased hydraulic gradient and developing the dammed water pockets along the stream. We identified one such seasonally active avalanche slope below SF1 and SF2 sites (Figures 5 and 6). While in the 3 November 2013 image of Figure 5 we observe a clear stream devoid of any avalanche deposit from the previous event of 2000 or from subsequent snow avalanches, in the 3 June 2014 and the 7 May 2015 images, we see thick fresh avalanche deposits covering the stream valley. Figure 5h highlights the seasonally active avalanche slope below the SF1 and SF2 sites, which accumulates significant snow mass during winter, often leading to slope failures relatively smaller than SF1 and SF2, but still large enough to form a deposition zone of ~500 m to 1 km in length. This large volume of snow, ice, and debris can accentuate any subsequent slope failure events, such as SF2, which occurred in the winter of 2021. We observed the available images of these  [41] is an efficient tool to estimate surface deformations and mass movements through the correlations between temporal images, even in the absence of groundcontrol points (GCPs). COSI-Corr was conclusively established in estimating the motion in mass movements [43,[65][66][67][68], glaciers [44], and sand dunes [69], with the relative accuracy of co-registration and correlation up to~1/20th of a pixel [41]. Thus, there were two prime objectives for us to assess the flow characteristics of SF1. First, to use COSI-Corr for quantifying the SV changes leading to SF1, and, second, to perform the simulation and flow analyses of the broken ice volume during the SF1 event using RAMMS model. We used the cloudless and snow-free Landsat 8 satellite Panchromatic band (15 m/pixel) to estimate the SV during 1-17 September 2015 (a year before the 2016 event) and for a similar period during 3-19 September 2016 (Figure 8a,b), immediately before SF1 and at the end of the yearly ablation season. We also estimated the yearly SVs during 17 September 2015 to 19 September 2016, highlighting the higher SVs in the avalanched zone of the hanging glacier (Figure 8a,b). Compared to 1-17 September 2015, the maximum SV in the avalanched part of the hanging glacier is~5 times higher during 3-19 September 2016 (Figure 8b) and the formed stress is strikingly evident in the region that leads to SF1 (middle and rightmost images of Figure 8a). The average velocity at the yearly scale (~12.3 ± 1.5 cm/day) for the avalanched part of the hanging glacier is also~2 times higher than the average short-term velocities (~6.8 ± 1.7 cm/day) (Figure 8b). The maximum SVs show a distinctly increasing trend before the event, while the overall average SVs do not follow a pattern. sonally active slope on 2 February 2021, while on 5 February 2021 Sentinel-2 image, it could be seen that a fresh cumulative snowfall of ~112 mm [22] significantly added to the snow cover. Figure 6 additionally shows the developed indentation from where the hanging glacier broke on 7 February 2021, causing SF2. Terrain profiles for SF1 and SF2 slopes (Figure 7), representing the terrain before either of the reported events occurred, estimate the maximum slope reaching as high as ~37°. The reported average slope for the entire SF2 mass [24] was even higher at ~40°. These slope values are entirely within the most avalanche-prone slope range reported from this region [18]. COSI-Corr [41] is an efficient tool to estimate surface deformations and mass movements through the correlations between temporal images, even in the absence of groundcontrol points (GCPs). COSI-Corr was conclusively established in estimating the motion in mass movements [43,[65][66][67][68], glaciers [44], and sand dunes [69], with the relative accuracy of co-registration and correlation up to ~1/20th of a pixel [41]. Thus, there were two prime objectives for us to assess the flow characteristics of SF1. First, to use COSI-Corr for quantifying the SV changes leading to SF1, and, second, to perform the simulation and flow analyses of the broken ice volume during the SF1 event using RAMMS model. We used the cloudless and snow-free Landsat 8 satellite Panchromatic band (15 m/pixel) to The visualisation of the surface deformations caused by the load effects, through an estimation of the surface displacements, has proven effective in observing the development and trajectory of mass movements [67,70]. As early as 1973, SVs were measured for the first time on an unstable cold-based hanging glacier in order to predict its collapse [71]. The observed velocities were shown to increase as a power-law function of time, up to infinity, at the theoretical time of failure [71]. This is also the characteristic signature of the stress buildup as a critical phenomenon in the case of various other naturally occurring ruptures, such as earthquakes [72], landslides [73], and rockfalls [74]. Velocity measurements performed on cold-based Weisshorn and Mönch hanging glaciers in Switzerland further confirmed this behaviour [75]. The validity of an SV monitoring approach was demonstrated in 2014 by the successful prediction of a hanging glacier "break-off" from the south face of the Grandes Jorasses (Monte Bianco, Italy), 10 days prior to the event [76]. However, owing to logistic issues, such efforts were understandably limited in their spatiotemporal domain, and all relied on field data collection for SV estimations. Our observations in the cases of SF1 and SF2, highlight that remote sensing-derived SVs have large spatial coverage and sufficient continuity and can represent a cost-effective and safer way to monitor hanging glaciers, and possibly even predict large and dangerous ice-rock avalanches in certain cases. However, any possible future predictability certainly requires more research in this direction to identify the types of slope failures which show statistically significant trends in SVs.  (Figure 8a,b), immediately before SF1 and at the end of the yearly ablation season. We also estimated the yearly SVs during 17 September 2015 to 19 September 2016, highlighting the higher SVs in the avalanched zone of the hanging glacier (Figure 8a,b). Compared to 1-17 September 2015, the maximum SV in the avalanched part of the hanging glacier is ~5 times higher during 3-19 September 2016 ( Figure 8b) and the formed stress is strikingly evident in the region that leads to SF1 (middle and rightmost images of Figure 8a). The average velocity at the yearly scale (~12.3 ± 1.5 cm/day) for the avalanched part of the hanging glacier is also ~2 times higher than the average short-term velocities (~6.8 ± 1.7 cm/day) (Figure 8b). The maximum SVs show a distinctly increasing trend before the event, while the overall average SVs do not follow a pattern. The visualisation of the surface deformations caused by the load effects, through an estimation of the surface displacements, has proven effective in observing the development and trajectory of mass movements [67,70]. As early as 1973, SVs were measured for the first time on an unstable cold-based hanging glacier in order to predict its collapse [71]. The observed velocities were shown to increase as a power-law function of time, up to infinity, at the theoretical time of failure [71]. This is also the characteristic signature of the stress build-up as a critical phenomenon in the case of various other naturally occur- The RAMMS simulation further revealed the flow characteristic of this entire SF1 event (Figures 8 and 9), which lasted for~3 min as per the modelled results. Based on the mechanical characteristics, the event can be divided into three prominent phases ( Figure 9): (1) break-off and avalanching in transect AB, (2) impact and shearing in transect BB', and (3) deposition in transect B'A'. Overall, the maximum velocity, maximum flow height, and maximum pressure reached were~90 m/s,~100 m, and~6000 kPa, respectively. The avalanching following the frontal break-off continued for~2300 m before it struck the valley. The flow velocity, momentum, and pressure reached their peaks during this phase. Immediately after the impact, the maximum shear stress reached its peak (~1200 kPa), and this zone of~200 m signifies the shearing and fragmentation of the ice core ( Figure 9). Finally, a depositional analysis in RAMMS simulates the deposition zone within a~3000 m length along the Ronti Gad stream (Figure 8c). The maximum deposition height is as much as~48.10 m in the middle regions of the deposition zone (Figure 8c), highlighting the magnitude of the SF1 event. This also points towards the possibility that SF1 in 2016 and the previous event in 2000 might have preconditioned the valley significantly in terms of sedimentation availability and topographic smoothness, leading to the enormity of debris flow during SF2 in 2021. We further generated the deposition outline for SF1 and it remarkably matched the ice-debris extent observed in 9 October 2016 post-SF1 Sentinel-2 image. In the Himalayan context, the modelled flow characteristics and mechanical parameters are comparable to those reported for the Gyari glacier avalanche in the Shyok Basin of the Indus River, which killed 139 people [55]. The volume of the Gyari avalanche was 4.31 × 10 6 m 3 , as opposed to~10 × 10 6 m 3 of SF1, but the maximum height of the deposits reached as high as 40 m [55]. The difference between the simulated maximum pressure for Gyari (~2500 kPa) and SF1 (~6000 kPa) was also proportional to the mass difference between the two events [55]. ring ruptures, such as earthquakes [72], landslides [73], and rockfalls [74]. Velocity measurements performed on cold-based Weisshorn and Mönch hanging glaciers in Switzerland further confirmed this behaviour [75]. The validity of an SV monitoring approach was demonstrated in 2014 by the successful prediction of a hanging glacier "break-off" from the south face of the Grandes Jorasses (Monte Bianco, Italy), 10 days prior to the event [76]. However, owing to logistic issues, such efforts were understandably limited in their spatiotemporal domain, and all relied on field data collection for SV estimations. Our observations in the cases of SF1 and SF2, highlight that remote sensing-derived SVs have large spatial coverage and sufficient continuity and can represent a cost-effective and safer way to monitor hanging glaciers, and possibly even predict large and dangerous ice-rock avalanches in certain cases. However, any possible future predictability certainly requires more research in this direction to identify the types of slope failures which show statistically significant trends in SVs. The RAMMS simulation further revealed the flow characteristic of this entire SF1 event (Figures 8 and 9), which lasted for ~3 min as per the modelled results. Based on the mechanical characteristics, the event can be divided into three prominent phases ( Figure  9): (1) break-off and avalanching in transect AB, (2) impact and shearing in transect BB', and (3) deposition in transect B'A'. Overall, the maximum velocity, maximum flow height, and maximum pressure reached were ~90 m/s, ~100 m, and ~6000 kPa, respectively. The avalanching following the frontal break-off continued for ~2300 m before it struck the valley. The flow velocity, momentum, and pressure reached their peaks during this phase.

SF2: Long-Term and Short-Term Surface Velocity Evolution
We further monitored the temporal evolution in the SVs of SF2 hanging glacier. For the long-term yearly SV estimation, we used available best quality multi-temporal Sentinel-2 images from the October months of all the years, to ensure similar shadow and surface conditions. Cloud-free images from the September months were not available for some of the years. 10 m/pixel bands of Sentinel-2 images have been reported to be extremely useful in detecting precursory motions in the mass movements using COSI-Corr [43]. We generated a time-series of the yearly SVs (Figure 10a) to observe the long-term changes and the zones of the hanging glacier showing anomalous SVs. During October 2016-October 2017, we can observe anomalous ice movement fields developing in SF2 break-off zone of the hanging glacier, below its bergschrund (a large crevasse often marking the head of a mountain glacier), with the maximum SV reaching as high as~39.2 ± 1.2 cm/day, while the average SV for the entire hanging glacier is~2.9 ± 1.2 cm/day (Figure 10a). The average SV kept on fluctuating through 2017-2019;~5.9 ± 1.2 cm/day during 2017-2018 and 3.1 ± 1.6 cm/day during 2018-2019 (Figure 10a). However, the most interesting observation is the distinct stress field growing across the entire break-off zone, below the bergschrund of the hanging glacier, with the maximum SVs reaching up to~51.3 ± 1.2 cm/day during 2017-2018 (Figure 10a).
Immediately after the impact, the maximum shear stress reached its peak (~1200 kPa), and this zone of ~200 m signifies the shearing and fragmentation of the ice core ( Figure 9). Finally, a depositional analysis in RAMMS simulates the deposition zone within a ~3000 m length along the Ronti Gad stream (Figure 8c). The maximum deposition height is as much as ~48.10 m in the middle regions of the deposition zone (Figure 8c), highlighting the magnitude of the SF1 event. This also points towards the possibility that SF1 in 2016 and the previous event in 2000 might have preconditioned the valley significantly in terms of sedimentation availability and topographic smoothness, leading to the enormity of debris flow during SF2 in 2021. We further generated the deposition outline for SF1 and it remarkably matched the ice-debris extent observed in 9 October 2016 post-SF1 Sentinel-2 image. In the Himalayan context, the modelled flow characteristics and mechanical parameters are comparable to those reported for the Gyari glacier avalanche in the Shyok Basin of the Indus River, which killed 139 people [55]. The volume of the Gyari avalanche was 4.31 × 10 6 m 3 , as opposed to ~10 × 10 6 m 3 of SF1, but the maximum height of the deposits reached as high as 40 m [55]. The difference between the simulated maximum pressure for Gyari (~2500 kPa) and SF1 (~6000 kPa) was also proportional to the mass difference between the two events [55].  ing the head of a mountain glacier), with the maximum SV reaching as high as ~39.2 ± 1.2 cm/day, while the average SV for the entire hanging glacier is ~2.9 ± 1.2 cm/day ( Figure  10a). The average SV kept on fluctuating through 2017-2019; ~5.9 ± 1.2 cm/day during 2017-2018 and ~3.1 ± 1.6 cm/day during 2018-2019 (Figure 10a). However, the most interesting observation is the distinct stress field growing across the entire break-off zone, below the bergschrund of the hanging glacier, with the maximum SVs reaching up to ~51.3 ± 1.2 cm/day during 2017-2018 (Figure 10a).  These estimations are in accordance with the conclusions provided by Shugar et al. [21], where they report the destabilisation in SF2 mass starting from 2016, i.e., after SF1, and displaying maximum displacements during the summer months of 2017 and 2018. These SV values indicate that SF2 mass was already unstable post-SF1, and experienced significant displacement by 2020. The indentation that was almost invisible in 2016, grew to be as large as~60-80 m by October 2020 (Figure 11), highlighting the major displacements during this period. The heavy precipitation [26,29] in the early weeks of February 2021 further added to this large, destabilised volume (~27 × 10 6 m 3 ) of rock, ice, and snow, leading to SF2. In fact, as per the Supplementary Data provided by Jiang et al. [22], the cumulative snowfall only on 5 February 2021 was~112 mm, considerably higher than any of the days of that winter season. SF1 in itself was a fairly large event that exerted a huge pressure of~6000 kPa at its impact point, and the localised tremor could have initiated the destabilisation in the bedrock of SF2 hanging glacier below its bergschrund. A valid question here might be that why a similar intensity 2000 event at SF1 site did not trigger the destabilisation in the bedrock of SF2 hanging glacier. The possible answer lies in the timing of both these events at SF1 site. The 2000 event occurred in the winter months, while the 2016 SF1 occurred in September, i.e., at the end of summer, when the bedrock is more prone to detachment. The powerful impact of SF1 might have been sufficient to cause such detachment in the bedrock of SF2 hanging glacier, which continues to grow during the following years ( Figure 11). However, it will need a separate detailed investigation to conclusively prove if impacts such as SF1 can trigger further glacier detachments. Here, we also report the short-term SV changes during the~40 days prior to SF2 (Figure 10b). While the average SV during 27 December 2020-16 January 2021 was~11.9 ± 2.5 cm/day, it nearly doubles to~23.2 ± 3.2 cm/day during 16-26 January 2021, and reaches over eight times this figure to~105.2 ± 31.7 cm/day during 26 January-5 February 2021 (Figure 10b). However, an even more dramatic increase is recorded in the maximum SVs, showing a~10-fold rise during these~40 days, with the maximum SVs reaching~481.8 ± 31.7 cm/day by 5 February 2021, 2 days prior to SF2 (Figure 10b). However, these estimates should be interpreted with care due to the proportionally high uncertainties associated with the SV measurements for the 26 January-5 February 2021 period (Figure 10b), owing to significantly high snow cover and poor illumination conditions in 5 February 2021 Sentinel-2 image (Figures 6b and 11f). Such conditions are reported to significantly affect the quality and accuracy of SV estimations [44], and therefore, we cannot conclusively say if the derived SVs for the 26 January-5 February 2021 period refer to actual displacements or are majorly a result of image artifact. Nevertheless, these short-term SVs provide a sequential picture of the pre-event flow characteristics, during a period of approximately one month, leading to SF2. These estimations are in accordance with the conclusions provided by Shugar et al. [21], where they report the destabilisation in SF2 mass starting from 2016, i.e., after SF1, and displaying maximum displacements during the summer months of 2017 and 2018. These SV values indicate that SF2 mass was already unstable post-SF1, and experienced significant displacement by 2020. The indentation that was almost invisible in 2016, grew to be as large as ~60-80 m by October 2020 (Figure 11), highlighting the major displacements during this period. The heavy precipitation [26,29] in the early weeks of February 2021 further added to this large, destabilised volume (~27 × 10 6 m 3 ) of rock, ice, and snow, leading to SF2. In fact, as per the supplementary data provided by Jiang et al. [22], the cumulative snowfall only on 5 February 2021 was ~112 mm, considerably higher than any of the days of that winter season. SF1 in itself was a fairly large event that exerted a huge pressure of ~6000 kPa at its impact point, and the localised tremor could have initiated the destabilisation in the bedrock of SF2 hanging glacier below its bergschrund. A valid question here might be that why a similar intensity 2000 event at SF1 site did not trigger the destabilisation in the bedrock of SF2 hanging glacier. The possible answer lies in the timing of both these events at SF1 site. The 2000 event occurred in the winter months, while the 2016 SF1 occurred in September, i.e., at the end of summer, when the bedrock is more prone to detachment. The powerful impact of SF1 might have been sufficient to cause such detachment in the bedrock of SF2 hanging glacier, which continues to grow during the following years (Figure 11). However, it will need a separate detailed investigation to conclusively prove if

SF2: Key Observations and During-Event Flow Characteristics
The total area of the avalanched part of SF2 hanging glacier was~0.27 km 2 with the maximum length of~850 m. The destabilisation in the bedrock started after SF1 in 2016 at an elevation of~5600 m asl (Figure 11), and on 7 February 2021, this sliding mass of rock and ice crashed into the Ronti Gad valley at an elevation of~3800 m asl. We simulated this SF2 event using high-resolution DTMs and RAMMS model. We contrast-stretched a high-resolution image of 11 February 2021 (Figure 12e), and the scar certainly suggested a significant removal of bedrock, as also modelled by Shugar et al. [21]. This event started as a landslide and turned into a rock-ice avalanche. Real-time images provided by PlanetScope [64] suggest that before Raini village, SF2 turned into a more fluidised debris flow, causing a flash flood in Rishiganga River. This flash flood took nearly 27 min to cover the~6 km of the distance between Raini village and Tapovan town (Figure 13), clocking a disastrous velocity of~4 m/s. This velocity is comparable to those of several reported flash floods in mountain catchments and is capable of carrying a significant amount of debris and sediment to a large distance, destroying any obstructions in the path [77,78]. The large debris and sedimentation spread along the river channel after SF2 is visible on the 10 February 2021 Sentinel-2 image of Figure 14d.  The RAMMS simulation provided us an informative account of the mechanical characteristics of SF2 (Figure 15), which further helped us in deciphering the various zones of this entire event ( Figure 16). The reported accounts [21][22][23][24]29] based on field validation, air surveys, and the visual interpretation of satellite images suggest that the total distance between SF2 site and the end of the rock-ice avalanche deposition zone was~13 km (~2 km of mountain slope and~11 km of river valley). From this point, the deposits started becoming more fluidised and turned into a debris flood [21]. The same published sources [21][22][23][24]29], using seismic modelling and videographic evidence provided by the local population, also put a constraint on the timeline of this event. The initial slope failure occurred at 10:21 Indian Standard Time (IST) (4:51 Greenwich Mean Time (GMT)) and it took slightly less than a minute to hit the valley [21]. PlanetScope image taken at 5:01:51 GMT (Figure 13), shows the SF2 front at the exact bend before Raini village where the rock-ice avalanche started turning into a debris flow. This indicates that SF2 rock-ice avalanche front took 10 min to travel to the transition point. However, this transition bend was also one of the sharpest and narrowest along the valley, stopping the runout and depositing the majority of the avalanched and entrained mass. This further consolidation of rock debris and ice along Rishiganga valley could have taken a few more minutes, as a long trail of avalanche can still be seen in Figure 13a. Thus, we calibrated our RAMMS model for the aforementioned distance and time to reach the optimal friction values of µ and ε. The final model outcome presented us with a time of~13.5 min for the entire SF2 avalanche event, from slope failure until complete deposition. This is in accordance with the estimated timeline above and the simulated timeline provided by Jiang et al. [22]. This consistency in the spatiotemporal premise for the final outcomes gave us the needed confidence in the simulated mechanical parameters for SF2 (Figures 15 and 16). The improving satellite remote sensing capabilities, such as the temporal resolutions, help to study natural hazards more efficiently, as in this case too, PlanetScope was able to capture the first image within~10 min of the event, giving us a reliable timestamp for calibrating the model.  The RAMMS simulation provided us an informative account of the mechanical characteristics of SF2 (Figure 15), which further helped us in deciphering the various zones of this entire event ( Figure 16). The reported accounts [21][22][23][24]29] based on field validation, air surveys, and the visual interpretation of satellite images suggest that the total distance between SF2 site and the end of the rock-ice avalanche deposition zone was ~13 km (~2 km of mountain slope and ~11 km of river valley). From this point, the deposits started becoming more fluidised and turned into a debris flood [21]. The same published sources [21][22][23][24]29], using seismic modelling and videographic evidence provided by the local population, also put a constraint on the timeline of this event. The initial slope failure occurred at 10:21 Indian Standard Time (IST) (4:51 Greenwich Mean Time (GMT)) and it took slightly less than a minute to hit the valley [21]. PlanetScope image taken at 5:01:51 GMT (Figure 13), shows the SF2 front at the exact bend before Raini village where the rock-ice avalanche started turning into a debris flow. This indicates that SF2 rock-ice avalanche front took 10 min to travel to the transition point. However, this transition bend was also one of the sharpest and narrowest along the valley, stopping the runout and depositing Figure 14. Temporal changes (a-d) in the SF1 deposition zone, as seen on Sentinel-2 images. Yellow outlines in the 18 October 2020 image highlight the regions along the channel, which significantly changed after SF1 and accumulated significant ice-debris deposits. Red arrows mark the avalanche slope. Sentinel-2 image IDs: S2A_OPER_MSI_L1C_TL_SGS__20160919T104034_A006494_T44RKU, S2A_OPER_MSI_L1C_TL_SGS__20161009T053631_A006780_T44RLU, S2A_OPER_MSI_L1C_TL_ VGS1_20201018T072817_A027801_T44RLU, and S2B_MSIL1C_20210210T051939_N0209_R062_T44RLU. The maximum velocity reached by SF2 was~65 m/s, while most of the flow height remains below 120 m ( Figure 16). These values also correspond well with other studies; Shugar et al. [21] reported flow height of up to 120 m and the maximum velocity of 60 m/s for this event, while Jiang et al. [22] modelled the maximum flow height of 100 with the maximum velocity of 65 m/s. The maximum pressure reaches~6000 kPa, while the momentum reaches~4000 kgms −1 in the GG' transect ( Figure 16). Based on the modelled mechanical characteristics, SF2 as a rock-ice avalanche can be divided into five prominent phases before it turns into a more fluidised debris flow in Dhauliganga River ( Figure 16 sistency in the spatiotemporal premise for the final outcomes gave us the needed confidence in the simulated mechanical parameters for SF2 (Figures 15 and 16). The improving satellite remote sensing capabilities, such as the temporal resolutions, help to study natural hazards more efficiently, as in this case too, PlanetScope was able to capture the first image within ~10 min of the event, giving us a reliable timestamp for calibrating the model.  phases before it turns into a more fluidised debris flow in Dhauliganga River ( Figure 16): (1) slope failure and sliding as slope runout (initial ~1800 m of transect GG', i.e., GA); (2) impact and shearing (~1800-~2400 m of transect GG', i.e., AB); (3) entrainment from valley sediments (~2400-~5000 m of transect GG', i.e., BC); (4) liquefaction and rock-ice avalanche front propagation as valley runout (~5000-~8900 m of transect GG', i.e., CD); and (5) deposition, detrainment, and transition into a more fluidised debris flood (~8900-12,500 m of transect GG', i.e., DG'). Video S2) confirms that while the avalanche was airborne briefly, it soon became grounded and that could explain the obliteration of this spur. We perceive that the presence of this spur offered a topographic barrier to the lowermost part of the eastern glacier, causing a slight pause in SF2 propagation (hence the dip in the mechanical parameters), containing it within a channel, and becoming obliterated in turn. This topographic barrier and friction also caused some deformation and possible breakage of the SF2 mass, as is evident from the spikes in the shear stress profile in Figure 16.  In the first phase of the landslide, the detached mass travelled~1800 m before it struck the valley. Our model results highlight an interesting and previously unreported aspect of this phase. At nearly 750 m below the terminus of the SF2 hanging glacier, we observed a dip in all the mechanical parameters, and, remarkably, we could also observe a dip in the difference DTM, highlighting the removal of significant rock volume by the SF2 rock-ice mass before it struck the valley (Figure 16). SF2 majorly consisted of rock volume and started as a landslide, and, therefore, when the large rock body hit the spur, it eroded it. The red ellipse in Figure 12 highlights this spur in the path of SF2 and the middle panel of Figure 16 shows the difference DTM highlighting the obliterated spur during SF2 propagation. This spur had a height range of~4-74 m with an average height of~25 m; The total area covered by it was~31,078 m 2 and the total volume of this spur removed by SF2 was~795,810 m 3 . This observation signified two aspects of this phase: (1) the model worked significantly well in predicting a hindrance in the flow, leading to the possible removal of this spur owing to the mechanical characteristics (high momentum and power) of the landslide, and (2) a significant rock volume of~0.8 × 10 6 m 3 was further added to the avalanched volume of~27 × 10 6 m 3 even before it reached the valley. Shugar et al. [21] postulated that the reason behind the intact lowermost part of the neighbouring larger eastern glacier post-SF2, might be the possibility that the avalanche became airborne for a while. A 3D visualization of our modelling result (Supplementary Materials Video S2) confirms that while the avalanche was airborne briefly, it soon became grounded and that could explain the obliteration of this spur. We perceive that the presence of this spur offered a topographic barrier to the lowermost part of the eastern glacier, causing a slight pause in SF2 propagation (hence the dip in the mechanical parameters), containing it within a channel, and becoming obliterated in turn. This topographic barrier and friction also caused some deformation and possible breakage of the SF2 mass, as is evident from the spikes in the shear stress profile in Figure 16.
In the impact and shearing phase, the SF2 mass hit the valley floor and all the mechanical parameters reached closer to their peak values during this phase. Immediately after the impact, the maximum shear stress reaches its peak (~1200 kPa) and this zone of~600 m signifies the shearing and fragmentation of the rock and ice core ( Figure 16). While the ice-debris from SF1 in 2016 mostly remains within~3 km of the deposition zone, covering the main channel and surroundings, the SF2 event, starting from slightly steeper slopes, releases over~2.5 times more of the rock-ice volume on the valley floor with previously deposited avalanche debris, significant seasonal snow (Figure 6b), and possibly dammed water ponds (Figures 4 and 5). This impact could have led to considerable ice and snow melt through the substantial change in momentum (Figure 16), leading to the transmission of energy by the falling and sliding rock and ice along the stretch. This immense rate of change in momentum and power, once the avalanche hit the valley, also caused the large burst whose bang could reportedly be heard from as far as~15 km from the avalanche site, and was later identified in the seismic readings from a station located over 100 km away [21]. For a similar event at Kolka in the North Ossetia-Alania region of the Caucasus Mountains in 2002, where a sliding hanging glacier of~20 × 10 6 m 3 failure volume hit and eroded a small glacier, a broad estimate, considering that the substantial amount of energy was dissipated during the fall and impact, suggested that~0.44 × 10 6 m 3 of ice could be instantly melted by a falling volume of 8 × 10 6 m 3 [79]. Considering that SF2 had a total volume of over~27 × 10 6 m 3 and it travelled an almost similar drop height as Kolka, it could have melted at least 3 times of~0.44 × 10 6 m 3 of ice instantly, i.e.,~1.32 × 10 6 m 3 of ice, corresponding to the~1/4th of the total SF2 ice volume, due to the direct impact and frictional sliding. This can explain why in contrast to most previously documented rock avalanches, very little debris was found at the base of the failed slope [21]. The Kolka event too, similar to SF2, did not allow the fallen ice and debris to deposit in the vicinity of the fall, but led to an unprecedented runout distance of~20 km, followed by flash mudflows. In such voluminous avalanches, granular friction can effectively be reduced by ice melt and associated fluidisation effects, leading to much larger travel distances than usual [79].
The next phase of SF2 consisted of entrainment from the valley sediments, mostly within the~2400-~5000 m of transect GG', i.e., in the BC region ( Figure 16). Point C in Figure 16 marks the sharp turn in the channel where previous ice avalanches of 2000 and 2016 from SF1 site had stopped. While the majority of the ice-debris deposit from SF1 appears to have melted away before the SF2 event, we cannot ignore the enormity of SF1 deposits (Figures 4, 5 and 8), when the maximum deposition height reached as high as 48 m. The SF1 event had not just altered the valley profile, but had also enriched it with significant volume of sediments and percolated meltwater in subsequent years. A mystery in explaining the long runout of SF2 has been that how only 20% of ice volume could displace nearly the entire 80% of rock volume, so far down the valley [29]. We propose that the long runout of SF2 was facilitated by the considerable entrainment from the valley sediments of past ice and snow avalanches (Figures 4 and 5) during this phase of SF2 in the BC transect ( Figure 16). To evaluate this hypothesis, we modelled two scenarios: (1) with the inclusion of differential erosion areas and erosion depth based on the modelled results of Jian et al. [22], and (2) without specifying the erosion areas and depth. By default, the model includes a critical shear stress that permits erosion only when the shear stress in any given cell exceeds the critical shear stress value for the onset of erosion. However, with the available information, i.e., in this case, the extent of SF1 deposits and modelled erosion by Jiang et al. [22], the inclusion of defined erosion zones can improve the model accuracy. RAMMS also provides the option to input erosion area parameters, and these were appropriately adjusted as per the instructions in RAMMS manual. Scenario 2 of not specifying the erosion areas and depth also meant that we had to consider a smaller turbulent friction (ε) parameter. While, for Scenario 1, after the model calibration, we opted for a value of ε = 200 m/s 2 , representing a transition from a solid-dominated granular flow to a muddier flow owing to sediment entrainment, for Scenario 2, we chose a value of ε = 100 m/s 2 , representing a relatively drier solid-dominated flow. As we expected, compared to Scenario 1, the rock-ice avalanche in Scenario 2 took nearly twice the time with an unrealistically low avalanche front momentum (only the~1/3rd of Scenario 1), for reaching point G' in Figure 16. This supports our assumption that the role of past avalanche deposits in entraining SF2 was considerably significant.
The subsequent phase of SF2 between~5000 m and~8900 m of transect GG' (i.e., CD) does not only show a relatively constricted valley profile, but is also marked by dense vegetation along the channel, which offered considerable resistance to SF2 propagation. The rapidly changing shear stress, flow pressure, and maximum flow height values in Figure 16 highlight this aspect. This friction led to further liquefaction of the avalanche core and the front took the shape of a well-defined rock-ice avalanche between this~5000-~8900 m of transect GG', i.e., in section CD. Point D represents another sharp turn in the valley and allowed for the temporary accumulation of material, leading to rise in the flow heights ( Figure 16).
The last phase of SF2 is marked by the deposition of transported material betweeñ 8900 and~12,500 m of transect GG', i.e., section DG'. Transect HH' in Figure 16 represents the place where Ronti Gad meets Rishiganga River and, from this point onwards, a vast volume of sediments and rocks were deposited, temporarily damming the river and forming a~700 m long lake [21]. A deposition analysis in RAMMS resulted in the total deposition volume of~15 × 10 6 m 3 . The maximum flow height reached up to~120 m. While Shugar et al. [21] also reported the maximum deposition height in the same range, based on DTM differencing, they reported the total deposition volume of~8 × 10 6 m 3 , between the Ronti Gad-Rishiganga River confluence and point G' in Figure 16. There are two reasons why our deposition volume is larger: (1) we observed and quantified the deposition on a larger stretch of transect GG', i.e., DG', and (2) our deposition estimates are in real-time, representing the terrain just after the rock-ice avalanche stopped at point G'. Figure 17 clarifies this point further. The 2021 post-SF2 DTM was generated using stereopairs acquired on 10 and 11 February and, by then, a significant volume of deposits fluidised and escaped as the debris flood. After 24 h of the event, the sediment plume could be seen on satellite images~150 km downstream of the source, and~80 times of the normal turbidity could be observed~500 km downstream within 8 days of the event [21]. This means that the deposited avalanche material continuously kept on fluidising, contributing to a constant, albeit gradually diminishing, debris flow. To evaluate our modelled deposition results, we hypothesised that while DTM difference elevations would show lower deposition heights, considering that they were generated using 10/11 February stereopairs, the profile pattern should match our modelled deposition heights. When we performed the profile analysis along four transects in the deposition region (Figure 16), we found that there was a remarkable similarity between the profile patterns between the modelled curve representing the deposition height immediately after the SF2 event on 7 February and the DTM difference curve representing the terrain on 11 February (Figure 17). This further proved the authenticity of our simulated results.
The SF2 event could be a result of the combined contribution of both the gradually destabilising slope and a deformational stress building up due to the winter snow accumulation load. A recent paper [47] suggested a possible polythermal regime for glaciers north of the study area, promoting both sliding and deformational forms of motion. Both pressure melting and creep-rate enhancement driven by stress concentration act as prime factors in producing noticeable surface motion under polythermal regime [80]. An approximate range for the existence of cold, polythermal, and temperate ice based on the relation between the mean annual air temperature (MAAT) and critical slope for failure on ramp-type glaciers [81] suggests that, for the average slope of SF2 site and assuming the polythermal nature of ice, a slope failure is possible within an MAAT range of −2.5 • C to −5 • C. The MAAT for SF1 and SF2 corresponds to~−5 • C, based on the published data for Chamoli district [82], extrapolated to the SF2 site using a constant environmental lapse rate of −6.5 • C [83]. This means that such slope failures can become increasingly possible, irrespective of the seasons, with changing climatic conditions. In their review on avalanching glacier instabilities, Faillettaz et al. [84] identified three different types of instabilities depending on the thermal properties of the ice/bedrock interface. If cold, the glacier experiences a critical acceleration following the power-law up to the break-off, and a prediction of the final break-off is possible. The type of motion observed in case of SF1 is closer to a possible cold regime, indicating a possibility of predicting its future break-offs. In case of a temperate regime, water plays a key role in the development of instability, and critical conditions promoting the final instability can still be identified [84]. However, in the case of a partly temperate or polythermal regime, predicting the actual break-off becomes difficult. While SF1 and SF2 are different in their mechanical characteristics, in both these cases, we observe certain abrupt but characteristically different changes in SVs, which raises another research question concerning whether their thermal regimes are different despite being in proximity, or if they experience a transition/evolution in their thermal regimes. The temporal evolution observed in the case of SF2 might represent the instability in a glacier experiencing a rapid transition from a cold to a temperate glacier bed due to percolation through the detached headwall. This possibility makes our SV estimates relevant to possibly open-up a new research frontier on understanding the transition in thermal regimes of hanging glaciers in these mountains, and how this transition might make them more prone to slope failures.

Conclusions
This research started with targeting three key research questions (Figure 2) mentioned in the Introduction Section, and we successfully derived several useful conclusions for them. The change in satellite-derived SVs should further be explored as an important precursor of ice-rock avalanching, and its temporal observation in high mountains can prove to be extremely helpful for predicting several SF1-type future slope failure events. Our SV estimates prove that continuously improving temporal and spatial resolution of satellite datasets can be effectively employed in performing such observations at unprecedented spatiotemporal scales. We also observed how SV changes vary between two constituently different slope failures, i.e., SF1 and SF2. In the case of the landslide causing ice-bedrock failure, i.e., SF2, SV anomalies developed and disappeared gradually and could be observed even at yearly scales. However, that makes predicting the actual time of failure difficult. The observed increase in SVs days-to-weeks prior to SF2 should be interpreted with caution, as relatively higher uncertainties are associated with them, due to high snow cover and poor illumination in satellite images. More of similar case studies are needed to infer if other SF2-type events also experience an immediate surge in SVs prior to their occurrence. On the other hand, for the frontal block failure, i.e., SF1, the SV anomalies are more localised and pronounced closer to the time of slope failure, as the glacier attains its critical geometry ( Figure 18). We further performed calibrated RAMMS simulations to infer the differences in the mechanical characteristics of the two large slope failures, i.e., SF1 and SF2, originating from the same release height and hitting the same valley floor. The simulated mechanical parameters were able to describe the during-event flow characteristics of both the events, and the changes in these parameters further helped in identifying the various zones of these mass movements.
Although the past events at SF1 site had not inflicted any direct damage to life and property, it is difficult to predict their future impacts in combination with other glacial hazards, such as SF2. While SF2 is an ice-bedrock failure, SF1 shows the characteristics of a frontal block failure as it tends to periodically replicate itself. More importantly, we can observe new ice accumulation at the avalanched slope of SF1 site since its last failure in 2016 ( Figure 18). The hanging glacier is recovering its previous critical geometry before SF1, and~0.052 km 2 (~20% of the total lost area in SF1) was recovered by October 2020 (Figure 18). Based on the available high-resolution images on Google Earth, we can observe that the full recovery usually takes 12-15 years, as happened after the 2000 event at SF1 site. Our results provide the first-hand account of noticeable SV changes days before SF1, and indicate that temporal SV monitoring of hanging glaciers might help to predict similar future events However, it is more important to first develop statistically robust "anomaly" characterisation routines for a variety of ice avalanche events to identify the cases with the most potential to be predicted through remote sensing. This chain of events clearly highlights the complexity of high mountain hazards, where considering the hazardous nature of an event in isolation is no longer sufficient. Instead, there is an increasing need to take into account the antecedent conditions, while making a holistic assessment of any high mountain hazard. similar future events However, it is more important to first develop statistically robust "anomaly" characterisation routines for a variety of ice avalanche events to identify the cases with the most potential to be predicted through remote sensing. This chain of events clearly highlights the complexity of high mountain hazards, where considering the hazardous nature of an event in isolation is no longer sufficient. Instead, there is an increasing need to take into account the antecedent conditions, while making a holistic assessment of any high mountain hazard.