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.

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).

Water resources are highly dependent on climatic variations. The quantification of climate change impacts on surface water availability is critical for agriculture production and flood management. The current study focuses on the projected streamflow variations in the transboundary Mangla Dam watershed. Precipitation and temperature changes combined with future water assessment in the watershed are projected by applying multiple downscaling techniques for three periods (2021–2039, 2040–2069, and 2070–2099). Streamflows are simulated by using the Soil and Water Assessment Tool (SWAT) for the outputs of five global circulation models (GCMs) and their ensembles under two representative concentration pathways (RCPs). Spatial and temporal changes in defined future flow indexes, such as base streamflow, average flow, and high streamflow have been investigated in this study. Results depicted an overall increase in average annual flows under RCP 4.5 and RCP 8.5 up until 2099. The maximum values of low flow, median flow, and high flows under RCP 4.5 were found to be 55.96 m³/s, 856.94 m³/s, and 7506.2 m³/s and under RCP 8.5, 63.29 m³/s, 945.26 m³/s, 7569.8 m³/s, respectively, for these ensembles GCMs till 2099. Under RCP 4.5, the maximum increases in maximum temperature (Tmax), minimum temperature (Tmin), precipitation (Pr), and average annual streamflow were estimated as 5.3 °C, 2.0 °C, 128.4%, and 155.52%, respectively, up until 2099. In the case of RCP 8.5, the maximum increase in these hydro-metrological variables was up to 8.9 °C, 8.2 °C, 180.3%, and 181.56%, respectively, up until 2099. The increases in Tmax, Tmin, and Pr using ensemble GCMs under RCP 4.5 were found to be 1.95 °C, 1.68 °C and 93.28% (2021–2039), 1.84 °C, 1.34 °C, and 75.88%(2040–2069), 1.57 °C, 1.27 °C and 72.7% (2070–2099), respectively. Under RCP 8.5, the projected increases in Tmax, Tmin, and Pr using ensemble GCMs were found as 2.26 °C, 2.23 °C and 78.65% (2021–2039), 2.73 °C, 2.53 °C, and 83.79% (2040–2069), 2.80 °C, 2.63 °C and 67.89% (2070–2099), respectively. Three seasons (spring, winter, and autumn) showed a remarkable increase in streamflow, while the summer season showed a decrease in inflows. Based on modeling results, it is expected that the Mangla Watershed will experience more frequent extreme flow events in the future, due to climate change. These results indicate that the study of climate change’s impact on the water resources under a suitable downscaling technique is imperative for proper planning and management of the water resources.

Previous studies have shown that glacier changes were heterogeneous in the western Karakoram, with the coexistence of retreating, advancing, and surging glaciers. However, it remains unclear that the mechanisms driving these changes. Based on the Shuttle Radar Topography Mission (SRTM) DEM and TerraSAR-X/TanDEM-X images (2014), this study presents glacier surface height changes in the Hunza Basin of the western Karakoram, employing the method of differential synthetic aperture radar interferometry (DInSAR). A slight negative glacier mass balance was observed in the Hunza Basin during 2000–2014. Surge-type glaciers would not have an obvious effect on overall mass balance in regional assessments over long-time scales. Further, glacier surface velocities in the Hunza Basin were estimated from Landsat images for the period of 1990–2018 by utilizing published data sets and Landsat images. Compared to the annual glacier surface velocities, 22 surge events were observed in seven surge-type glaciers in the Hunza Basin. Glacier flow can be attributed to thermally and hydrologically control, and the geomorphological characteristics of different individuals. This study gives us a new insight into the situation of the “Karakoram anomaly” under the background of glacier mass loss in the high mountains of Asia (HMA).