Multi-model assessment of glacio-hydrological changes in central Karakoram, Pakistan

The multi-model assessment of glacio-hydrological regimes can enhance our understanding of glacier response to climate change. This improved knowledge can uplift our computing abilities to estimate the contributing components of the river discharge. This study examined and compared the hydrological responses in the glacier-dominated Shigar River basin (SRB) under various climatic scenarios using a semi-distributed Modified Positive Degree Day Model (MPDDM) and a distributed Glacio-hydrological Degree-day Model (GDM). Both glacio-hydrological models were calibrated and validated against the observed hydro-meteorological data from 1988–1992 and 1993–1997. Temperature and precipitation data from Shigar and Skardu meteorological stations were used along with field estimated degree-day factor, temperature, and precipitation gradients. The results from both models indicate that the snow and ice melt are vital contributors to sustain river flow in the catchment. However, MPDDM estimated 68% of rain and baseflow contribution to annual river runoff despite low precipitation during the summer monsoon, while GDM estimated 14% rain and baseflow contribution. Likewise, MPDDM calculated 32%, and GDM generated 86% of the annual river runoff from snow and ice melt. MPDDM simulated river discharge with 0.86 and 0.78 NSE for calibration and validation, respectively. Similarly, GDM simulated river discharge with improved accuracy of 0.87 for calibration and 0.84 NSE for the validation period. The snow and ice melt is significant in sustaining river flow in the SRB, and substantial changes in melt characteristics of snow and ice are expected to have severe consequences on seasonal water availability. Based on the sensitivity analysis, both models’ outputs are highly sensitive to the variation in temperature. Furthermore, compared to MPDDM, GDM simulated considerable variation in the river discharge in climate scenarios, RCP4.5 and 8.5, mainly due to the higher sensitivity of GDM model outputs to temperature change. The integration of an updated melt module and two reservoir baseflow module in GDM is anticipated to advance the representation of hydrological components, unlike one reservoir baseflow module used separately in MPDDM. The restructured melt and baseflow modules in GDM have fundamentally enriched our perception of glacio-hydrological dynamics in the catchment.