Supra-glacial debris cover is key to glacier ablation through increasing (thin debris layer) or decreasing (thick debris layer) melt rates, thereby regulating the mass balance of a glacier and its meltwater runoff. The thickening or lateral expansion of supra-glacial debris cover correlates with a reduction of glacier ablation and, consequently, runoff generation, which is also considered to be an influential factor on the rheology and dynamics of a glacierized system. Studies on supra-glacial debris cover have recently attracted wide attention especially for glaciers in the Himalayas and Karakoram, where the glaciers have heterogeneously responded to climate change. In this study, we used 32 images from the Landsat Thematic Mapper, Enhanced Thematic Mapper Plus, and Operational Land Imager archive, going back to 1990, which are available on the Google Earth Engine cloud-computing platform, to map the supra-glacial debris cover in the Hunza Valley, Karakoram, Pakistan, based on a band ratio segmentation method (normalized difference snow index [NDSI] < 0.4), Otsu thresholding, and machine learning algorithms. Compared with manual digitization, the random forest (RF) model was found to have the greatest accuracy in identifying supra-glacial debris, with a Kappa coefficient of 0.94 ± 0.01 and an overall accuracy of 95.5 ± 0.9%. Overall, the supra-glacial debris cover in the study area showed an increasing trend, and the total area expanded by 8.1–21.3% for various glaciers from 1990 to 2019. The other two methods (Otsu thresholding and NDSI < 0.4) generally overestimated the supra-glacial debris covered area, by 36.3 and 18.8%, respectively, compared to that of the RF model. The supra-glacial debris cover has migrated upward on the glaciers, with intensive variation near the equilibrium-line altitude zone (4,500–5,500 m a.s.l.). The increase in ice or snow avalanche activity at high altitudes may be responsible for this upward expansion of supra-glacial debris cover in the Hunza Valley, which is attributed to the combined effect of temperature decrease and precipitation increase in the study area.

The terrestrial water storage anomaly (TWSA) gap between the Gravity Recovery and Climate Experiment (GRACE) and its follow-on mission (GRACE-FO) is now a significant issue for scientific research in high-resolution time-variable gravity fields. This paper proposes the use of singular spectrum analysis (SSA) to predict the TWSA derived from GRACE. We designed a case study in six regions in China (North China Plain (NCP), Southwest China (SWC), Three-River Headwaters Region (TRHR), Tianshan Mountains Region (TSMR), Heihe River Basin (HRB), and Lishui and Wenzhou area (LSWZ)) using GRACE RL06 data from January 2003 to August 2016 for inversion, which were compared with Center for Space Research (CSR), Helmholtz-Centre Potsdam-German Research Centre for Geosciences (GFZ), Jet Propulsion Laboratory (JPL)’s Mascon (Mass Concentration) RL05, and JPL’s Mascon RL06. We evaluated the accuracy of SSA prediction on different temporal scales based on the correlation coefficient (R), Nash–Sutcliffe efficiency (NSE), and root mean square error (RMSE), which were compared with that of an auto-regressive and moving average (ARMA) model. The TWSA from September 2016 to May 2019 were predicted using SSA, which was verified using Mascon RL06, the Global Land Data Assimilation System model, and GRACE-FO results. The results show that: (1) TWSA derived from GRACE agreed well with Mascon in most regions, with the highest consistency with Mascon RL06 and (2) prediction accuracy of GRACE in TRHR and SWC was higher. SSA reconstruction improved R, NSE, and RMSE compared with those of ARMA. The R values for predicting TWS in the six regions using the SSA method were 0.34–0.98, which was better than those for ARMA (0.26–0.97), and the RMSE values were 0.03–5.55 cm, which were better than the 2.29–5.11 cm RMSE for ARMA as a whole. (3) The SSA method produced better predictions for obvious periodic and trending characteristics in the TWSA in most regions, whereas the detailed signal could not be effectively predicted. (4) The predicted TWSA from September 2016 to May 2019 were basically consistent with Global Land Data Assimilation System (GLDAS) results, and the predicted TWSA during June 2018 to May 2019 agreed well with GRACE-FO results. The research method in this paper provides a reference for bridging the gap in the TWSA between GRACE and GRACE-FO.

Multi-temporal glacier inventories provide key information about the glaciers, their characteristics, and changes and are inevitable for glacier modelling and investigating geodetic mass changes. However, to date, no consistent multi-temporal glacier inventory for the whole of the Karakoram exists, negatively affecting the monitoring of spatio-temporal variations in glaciers' geometric parameters and their related applications. We used a semi-automatic method combining automatic segmentation and manual correction and produced a multi-temporal Karakoram glacier inventory (KGI) compiled from Landsat TM/ETM+/OLI (Thematic Mapper, Enhanced Thematic Mapper Plus, and Operational Land Imager) images for the 1990s, 2000s, 2010s, and 2020s. Our assessments using independent multiple digitisation of 37 glaciers show that the KGI is sufficiently accurate, with an overall uncertainty of ±3.68 %. We also performed uncertainty evaluation for the contiguous glacier polygons using a buffer of half a pixel, which resulted in an average mapping uncertainty of ±5.21 %. We calculated more than 20 attributes for each glacier, including coordinates, area, supraglacial debris area, date information, and topographic parameters derived from the ASTER GDEM (Advanced Spaceborne Thermal Emission and Reflection Radiometer global digital elevation model). According to KGI-2020s, approximately 10 500 alpine glaciers (>0.01 km² each) cover an area of 22 510±828 km² of which 10.18±0.38 % (2290±84 km²) is covered by supraglacial debris. Over the past 3 decades, the glaciers experienced a loss of clean ice and/or snow area but a gain in supraglacial debris. Supraglacial debris cover has increased by 17.63±1.44 % (343.30±27.95 km²), while non-debris-covered glaciers decreased by 1.56±0.24 % (319.85±49.92 km²). The total glacier area was relatively stable and showed only a slight insignificant increase of 23.45±28.85 km² (0.10±0.13 %). The glacier area has declined by 3.27±0.24 % in the eastern Karakoram, while the glacier area slightly increased in central (0.65±0.10 %) and western Karakoram (1.26±0.11 %). Supraglacial debris has increased over the whole of Karakoram, especially in areas above 4200 m a.s.l. (above sea level), showing an upward shift. The glacier area changes were characterised by strong spatial heterogeneity, influenced by surging and advancing glaciers. However, due to global warming, the glaciers are on average retreating. This is in particular true for small and debris-free glaciers. The multi-temporal KGI data are available at the National Cryosphere Desert Data Center of China: https://doi.org/10.12072/ncdc.glacier.db2386.2022 (F. Xie et al., 2022).

The China–Pakistan international Karakoram Highway passes through the core area of the “Karakoram Anomaly,” whose glaciers have maintained or increased their mass during a period when most glaciers worldwide have receded. We synthesized the literature and used remote-sensing techniques to review the types, distribution, characteristics, causes and frequency of major glacial hazards along the Karakoram Highway. We found that the glacier-related hazards could be divided into direct and indirect hazards, including glacier surges, glacial lake outburst floods, and glacial floods, which are concentrated in East Pamir and the Hunza River Basin. In the past 100 years, hazards from glaciers surges and glacial floods only occurred once and twice, respectively, which appear suddenly, with the hazard-causing process being short-lived and occurring mainly in the summer. Glacial lake outburst floods mainly occur in the spring and summer in the Hunza River Basin. Among these, ice-dammed lakes have the highest frequency of flooding, their formation and outbursts being closely related to the sudden advancement of surge-type glaciers. Under the background of global climate warming, we speculate that the glacier surge cycle may shorten and the frequency of the formation and outbursts in the glacial lakes may increase. In the future, we should combine models and new field observations to simulate, and deepen our understanding of the physical mechanisms of different glacier-related hazards. In particular, on-site monitoring should be carried out, to include the evolution of glaciers subglacial hydrological systems, the thermal state at the base of the glaciers, and the opening and closing of drainage channels at the base of the ice dams.