Global warming is one of the main challenges of recent times. The glaciers are melting faster than expected which has resulted in global mean sea level rise and increased the risk of floods. The development of modern remote sensing technology has made it possible to obtain images more frequently than ever before. On the other side, the availability of high-performance computing hardware and processing techniques have made it possible to provide a cost-effective solution to monitor the temporal changes of glaciers at a large scale. In this study, supervised machine learning methods are investigated for automatic classification of glacier covers from multi-temporal Sentinel-2 imagery using texture, topographic, and spectral data. Three most commonly used supervised machine learning techniques were investigated: support vector machine (SVM), artificial neural network (ANN) and random forest (RF). The proposed method was employed on the data obtained from Passu watershed in Hunza Basin located along the Hunza river in Pakistan. Three main classes were considered: glaciers, debris-covered glaciers and non-glaciated areas. The data was split into training (70%) and testing datasets (30%). Finally, an area-based accuracy assessment was performed by comparing the results obtained for each classifier with the reference data. Experiments showed that the results produced for all classifiers were highly accurate and visually more consistent with the depiction of glacier cover types. For all experiments, random forest performed the best (Kappa= 0.95, f-measure =95.06%) on all three classes compared to ANN (Kappa= 0.92, f-measure =92.05) and SVM (Kappa= 0.89, f-measure =91.86% on average). The high classification accuracy obtained to distinguished debris-covered glaciers using our approach will be useful to determine the actual available water resources which can be further helpful for hazard and water resource management.

Gravity Recovery and Climate Experiment (GRACE) and satellite altimetry are suitable for the precise measurement of terrestrial water storage (TWS) and lake water level variations from space. In this study, two GRACE solutions, namely, spherical harmonics (SH) and mascon (MSC), are utilized with the Global Land Data Assimilation System (GLDAS) model to estimate the spatial and temporal variations of TWS in the Upper Indus Basin (UIB) for the study period of January 2003 to December 2016. The TWS estimated by SH, MSC, and the GLDAS model are consistent and generally show negative trends of -4.47 ± 0.38 mm/year, -4.81 ± 0.49 mm/year, and -3.77 ± 0.46 mm/year, respectively. Moreover, we use the GLDAS model data to understand the roles of variations in land surface state variables (snow water equivalent (SWE), soil moisture, and canopy water storage) in enhancing or dissipating the TWS in the region. Results indicate that SWE, which has a significant contribution to GRACE TWS variability, is an important parameter. Spearman's rank correlations are calculated to demonstrate the relationship of the GLDAS land surface state variables and the GRACE signals. A highly positive correlation between SWE with TWS is estimated by SH and MSC as 0.691 and 0.649, respectively, indicating that the TWS signal is mainly reliant on snow water in the study region. The ground water storages estimated by SH and MSC solutions are nearly stable with slight increasing trends of 0.63 ± 0.48 mm/year and 0.29 ± 0.51 mm/year, respectively. We also take advantage of the potential of satellite altimetry in measuring lake water level variations, and our result indicates that Cryosat-2 SARin mode altimetry data can be used in estimating small water bodies accurately in the high mountainous region of the UIB. Moreover, the climate indices data of El-Niño Southern Oscillation and Pacific Decadal Oscillation are analyzed to determine the influence of pacific climatic variability on TWS.

Mountains regions like Gilgit-Baltistan (GB) province of Pakistan are solely dependent on seasonal snow and glacier melt. In Indus basin which forms in GB, there is a need to manage water in a sustainable way for the livelihood and economic activities of the downstream population. It is important to monitor water resources that include glaciers, snow-covered area, lakes, etc., besides traditional hydrological (point-based measurements by using the gauging station) and remote sensing-based studies (traditional satellite-based observations provide terrestrial water storage (TWS) change within few centimeters from the earth’s surface); the TWS anomalies (TWSA) for the GB region are not investigated. In this study, the TWSA in GB region is considered for the period of 13 years (from January 2003 to December 2016). Gravity Recovery and Climate Experiment (GRACE) level 2 monthly data from three processing centers, namely Centre for Space Research (CSR), German Research Center for Geosciences (GFZ), and Jet Propulsion Laboratory (JPL), System Global Land Data Assimilation System (GLDAS)-driven Noah model, and in situ precipitation data from weather stations, were used for the study investigation. GRACE can help to forecast the possible trends of increasing or decreasing TWS with high accuracy as compared to the past studies, which do not use satellite gravity data. Our results indicate that TWS shows a decreasing trend estimated by GRACE (CSR, GFZ, and JPL) and GLDAS-Noah model, but the trend is not significant statistically. The annual amplitude of GLDAS-Noah is greater than GRACE signal. Mean monthly analysis of TWSA indicates that TWS reaches its maximum in April, while it reaches its minimum in October. Furthermore, Spearman’s rank correlation is determined between GRACE estimated TWS with precipitation, soil moisture (SM) and snow water equivalent (SWE). We also assess the factors, SM and SWE which are the most efficient parameters producing GRACE TWS signal in the study area. In future, our results with the support of more in situ data can be helpful for conservation of natural resources and to manage flood hazards, droughts, and water distribution for the mountain regions.