A deadly cascade: A catastrophic landslide in Uttarakhand state in India on February 2021 damaged two hydropower plants, and more than 200 people were killed or are missing. Shugar et al. describe the cascade of events that led to this disaster. A massive rock and ice avalanche roared down a Himalayan valley, turning into a deadly debris flow upstream from the first of the two hydropower plants. The sequence of events highlights the increasing risk in the Himalayas caused by increased warming and development. Science , abh4455, this issue p. 300

Knowledge of the ice thickness distribution of glaciers and ice caps is an important prerequisite for many glaciological and hydrological investigations. A wealth of approaches has recently been presented for inferring ice thickness from characteristics of the surface. With the Ice Thickness Models Intercomparison eXperiment (ITMIX) we performed the first coordinated assessment quantifying individual model performance. A set of 17 different models showed that individual ice thickness estimates can differ considerably – locally by a spread comparable to the observed thickness. Averaging the results of multiple models, however, significantly improved the results: on average over the 21 considered test cases, comparison against direct ice thickness measurements revealed deviations on the order of 10 ± 24 % of the mean ice thickness (1σ estimate). Models relying on multiple data sets – such as surface ice velocity fields, surface mass balance, or rates of ice thickness change – showed high sensitivity to input data quality. Together with the requirement of being able to handle large regions in an automated fashion, the capacity of better accounting for uncertainties in the input data will be a key for an improved next generation of ice thickness estimation approaches.

Ice volume estimates are crucial for assessing water reserves stored in glaciers. Due to its large glacier coverage, such estimates are of particular interest for the Himalayan–Karakoram (HK) region. In this study, different existing methodologies are used to estimate the ice reserves: three area–volume relations, one slope-dependent volume estimation method, and two ice-thickness distribution models are applied to a recent, detailed, and complete glacier inventory of the HK region, spanning over the period 2000–2010 and revealing an ice coverage of 40 775 km². An uncertainty and sensitivity assessment is performed to investigate the influence of the observed glacier area and important model parameters on the resulting total ice volume. Results of the two ice-thickness distribution models are validated with local ice-thickness measurements at six glaciers. The resulting ice volumes for the entire HK region range from 2955 to 4737 km³, depending on the approach. This range is lower than most previous estimates. Results from the ice thickness distribution models and the slope-dependent thickness estimations agree well with measured local ice thicknesses. However, total volume estimates from area-related relations are larger than those from other approaches. The study provides evidence on the significant effect of the selected method on results and underlines the importance of a careful and critical evaluation.

On 3 October 2023, a multihazard cascade in the Sikkim Himalaya, India, was triggered by 14.7 million cubic meters of frozen lateral moraine collapsing into South Lhonak Lake. The impact generated an ~20-meter tsunami-like impact wave, which breached the moraine and drained ~50 million cubic meters of the lake’s water. The ensuing glacial lake outburst flood (GLOF) eroded ~270 million cubic meters of sediment, which overwhelmed infrastructure, including hydropower installations along the Teesta River. The physical scale and human and economic impacts of this event prompt urgent reflection on the role of climate change and human activities in exacerbating such disasters. Insights into multihazard evolution are pivotal for informing policy development, enhancing early warning systems (EWS), and spurring paradigm shifts in GLOF risk management strategies in the Himalaya and other mountain environments.

In the course of glacier retreat, new glacier lakes can develop. As such lakes can be a source of natural hazards, strategies for predicting future glacier lake formation are important for an early planning of safety measures. In this article, a multi-level strategy for the identification of overdeepened parts of the glacier beds and, hence, sites with potential future lake formation, is presented. At the first two of the four levels of this strategy, glacier bed overdeepenings are estimated qualitatively and over large regions based on a digital elevation model (DEM) and digital glacier outlines. On level 3, more detailed and laborious models are applied for modeling the glacier bed topography over smaller regions; and on level 4, special situations must be investigated in-situ with detailed measurements such as geophysical soundings. The approaches of the strategy are validated using historical data from Trift Glacier, where a lake formed over the past decade. Scenarios of future glacier lakes are shown for the two test regions Aletsch and Bernina in the Swiss Alps. In the Bernina region, potential future lake outbursts are modeled, using a GIS-based hydrological flow routing model. As shown by a corresponding test, the ASTER GDEM and the SRTM DEM are both suitable to be used within the proposed strategy. Application of this strategy in other mountain regions of the world is therefore possible as well.

The Cordillera Blanca in Peru has been the scene of rapid deglaciation for many decades. One of numerous lakes formed in the front of the retreating glaciers is the moraine-dammed Lake Palcacocha, which drained suddenly due to an unknown cause in 1941. The resulting Glacial Lake Outburst Flood (GLOF) led to dam failure and complete drainage of Lake Jircacocha downstream, and to major destruction and thousands of fatalities in the city of Huaráz at a distance of 23 km. We chose an integrated approach to revisit the 1941 event in terms of topographic reconstruction and numerical back-calculation with the GIS-based open-source mass flow/process chain simulation framework r.avaflow, which builds on an enhanced version of the Pudasaini (2012) two-phase flow model. Thereby we consider four scenarios: (A) and (AX) breach of the moraine dam of Lake Palcacocha due to retrogressive erosion, assuming two different fluid characteristics; (B) failure of the moraine dam caused by the impact of a landslide on the lake; and (C) geomechanical failure and collapse of the moraine dam. The simulations largely yield empirically adequate results with physically plausible parameters, taking the documentation of the 1941 event and previous calculations of future scenarios as reference. Most simulation scenarios indicate travel times between 36 and 70 min to reach Huaráz, accompanied with peak discharges above 10 000 m³ s⁻¹. The results of the scenarios indicate that the most likely initiation mechanism would be retrogressive erosion, possibly triggered by a minor impact wave and/or facilitated by a weak stability condition of the moraine dam. However, the involvement of Lake Jircacocha disguises part of the signal of process initiation farther downstream. Predictive simulations of possible future events have to be based on a larger set of back-calculated GLOF process chains, taking into account the expected parameter uncertainties and appropriate strategies to deal with critical threshold effects.