Spatiotemporal Variability Analysis of Glaciers in the Hindukush Region of Pakistan using Remote Sensing Data

Glaciers in northern Pakistan are a prime source of freshwater, providing headwater in the Indus river system and serving as a lifeline to millions of people in the region. These glaciers undergo continuous changes by melting due to global warming or accumulation due to snowfall/precipitation at higher altitudes. In this study, we used remote sensing data to quantify glacier changes in spatiotemporal variability in the past three decades. Five glaciers in the Gilgit region (near the junction of the Hindukush and Karakoram Mountains) with an extent of less than 5 square kilometers were selected, namely Phakor glacier, Karamber glacier, East Gammu glacier, Bhort glacier, and Bad-e-Swat glacier. The fluctuations in these glaciers were monitored using a digital elevation model (DEM) and a cloud-free continuous series of Landsat satellite pictures from the minimal snow cover season. The annual climatic trends were studied through spatially interpolated gridded climate data WοrldClim version-1 climate database for 1970 – 2000. We used it to study the variations of minimum and maximum temperature, solar radiation, and precipitation through the preparation of sub-sets from the original global grids. Because of its exact delineation in the Gilgit sub-basin, the characterized watersheds were visually compared to optical Landsat 8 OLI data for mountainous ridge matching, revealing that SRTM 30m (radar-based) demonstrated greater accuracy than other DEMs. The temporal assessment of Bhort, Bad-e-Sawat, East Gammu, Karamber, and other rivers was also carried out. It is observed that the glaciers in the Gilgit watershed are rather stable. The little variability of glaciers is due to their geographic condition, altitude, topography, and orientation. Validation of the mapped glacier classes has been performed to check the accuracy assessment through an error matrix method. The kappa coefficient from the error matrix has been calculated to be 84 %. The study makes a critical input to a greater understanding of watershed controlling and hydrological processes in the upper Indus catchment's Gilgit watershed.