In contrast to widespread glacier retreat evidenced globally, glaciers in the Karakoram region have exhibited positive mass balances and general glacier stability over the past decade. Snow and glacier meltwater from the Karakoram and the western Himalayas, which supplies the Indus River Basin, provide an essential source of water to more than 215 million people, either directly, as potable water, or indirectly, through hydroelectric generation and irrigation for crops. This study focuses on water resources in the Upper Indus Basin (UIB) which combines the ranges of the Hindukush, Karakoram and Himalaya (HKH). Specifically, we focus on the Gilgit River Basin (GRB) to inform more sustainable water use policy at the sub-basin scale. We employ two degree-day approaches, the Spatial Processes in Hydrology (SPHY) and Snowmelt Runoff Model (SRM), to simulate runoff in the GRB during 2001-2012. The performance of SRM was poor during July and August, the period when glacier melt contribution typically dominates runoff. Consequently, SPHY outperformed SRM, likely attributable to SPHY's ability to discriminate between glacier, snow, and rainfall contributions to runoff during the ablation period. The average simulated runoff revealed the prevalent snowmelt contribution as 62%, followed by the glacier melt 28% and rainfall 10% in GRB. We also assessed the potential impact of climate change on future water resources, based on two Representative Concentration Pathways (RCP) (RCP 4.5 and RCP 8.5). We estimate that summer flows are projected to increase by between 5.6% and 19.8% due to increased temperatures of between 0.7 and 2.6 degrees C over the period 2039-2070. If realized, increased summer flows in the region could prove beneficial for a range of sectors, but only over the short to medium term and if not associated with extreme events. Long-term projections indicate declining water resources in the region in terms of snow and glacier melt.

In many high altitude river basins, the hydro-climatic regimes and the spatial and temporal distribution of precipitation are little known, complicating efforts to quantify current and future water availability. Scarce, or non-existent, gauged observations at high altitudes coupled with complex weather systems and orographic effects further prevent a realistic and comprehensive assessment of precipitation. Quantifying the contribution from seasonal snow and glacier melt to the river runoff for a high altitude, melt dependent region is especially difficult. Global scale precipitation products, in combination with precipitation-runoff modelling may provide insights to the hydro-climatic regimes for such data scarce regions. In this study two global precipitation products; the high resolution (0.1 degrees x 0.1 degrees), newly developed ERA5-Land, and a coarser resolution (0.55 degrees x 0.55 degrees) JRA-55, are used to simulate snow/glacier melts and runoff for the Gilgit Basin, a sub-basin of the Indus. A hydrological precipitation-runoff model, the Distance Distribution Dynamics (DDD), requires minimum input data and was developed for snow dominated catchments. The mean of total annual precipitation from 1995 to 2010 data was estimated at 888 mm and 951 mm by ERA5-Land and JRA-55, respectively. The daily runoff simulation obtained a Kling Gupta efficiency (KGE) of 0.78 and 0.72 with ERA5-Land and JRA-55 based simulations, respectively. The simulated snow cover area (SCA) was validated using MODIS SCA and the results are quite promising on daily, monthly and annual scales. Our result showed an overall contribution to the river flow as about 26% from rainfall, 37-38% from snow melt, 31% from glacier melt and 5% from soil moisture. These melt simulations are in good agreement with the overall hydro-climatic regimes and seasonality of the area. The proxy energy balance approach in the DDD model, used to estimate snow melt and evapotranspiration, showed robust behaviour and potential for being employed in data poor basins. (c) 2021 Published by Elsevier B.V.

We analyse an ensemble of statistically downscaled Global Climate Models (GCMs) to investigate future water availability in the Upper Indus Basin (UIB) of Pakistan for the time horizons when the global and/or regional warming levels cross Paris Agreement (PA) targets. The GCMs data is obtained from the 5th Phase of Coupled Model Inter-Comparison Project under two Representative Concentration Pathways (RCP4.5 and RCP8.5). Based on the five best performing GCMs, we note that global 1.5 degrees C and 2.0 degrees C warming thresholds are projected in 2026 and 2047 under RCP4.5 and 2022 and 3036 under RCP8.5 respectively while these thresholds are reached much earlier over Pakistan i.e. 2016 and 2030 under RCP4.5 and 2012 and 2025 under RCP8.5 respectively. Interestingly, the GCMs with the earliest emergence at the global scale are not necessarily the ones with the earliest emergence over Pakistan, highlighting spatial non-linearity in GCMs response. The emergence of 2.0 degrees C warming at global scale across 5 GCMs ranges from 2031 (CCSM4) to 2049 (NorESM) under RCP8.5. Precipitation generally exhibits a progressive increasing trend with stronger changes at higher warming or radiative forcing levels. Hydrological simulations representing the historical, 1.5 degrees C and 2.0 degrees C global and region warming time horizons indicate a robust but seasonally varying increase in the inflows. The highest inflows in the baseline and future are witnessed in July. However, the highest future increase in inflows is projected in October under RCP4.5 (37.99% and 65.11% at 1.5 degrees C and 2.0 degrees C) and in April under RCP8.5 (37% and 62.05% at 1.5 degrees C and 2.0 degrees C). These hydrological changes are driven by increases in the snow and glacial melt contribution, which are more pronounced at 2.0 degrees C warming level. These findings should help for effective water management in Pakistan over the coming decades. (c) 2021 Elsevier B.V. All rights reserved.

Investigating the trends in the major climatic variables over the Upper Indus Basin (UIB) region is difficult for many reasons, including highly complex terrain with heterogeneous spatial precipitation patterns and a scarcity of gauge stations. The Weather Research and Forecasting (WRF) model was applied to simulate the spatiotemporal variability of precipitation and temperature over the Indus Basin from 2000 through 2015 with boundary conditions derived from the Climate Forecast System Reanalysis (CFSR) data. The WRF model was configured with three nested domains (d01-d03) with horizontal resolutions increasing inward from 36 km to 12 km to 4 km horizontal resolution, respectively. These simulations were a continuous run with a spin-up year (i.e., 2000) to equilibrate the soil moisture, snow cover, and temperature at the beginning of the simulation. The simulations were then compared with TRMM and station data for the same time period using root mean squared error (RMSE), percentage bias (PBIAS), mean bias error (MBE), and the Pearson correlation coefficient. The results showed that the precipitation and temperature simulations were largely improved from d01 to d03. However, WRF tended to overestimate precipitation and underestimate temperature in all domains. This study presents high-resolution climatological datasets, which could be useful for the study of climate change and hydrological processes in this data-sparse region.

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 & x2212;4.47 & x00B1; 0.38 mm/year, & x2212;4.81 & x00B1; 0.49 mm/year, and & x2212;3.77 & x00B1; 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 & x2019;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 & x00B1; 0.48 mm/year and 0.29 & x00B1; 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 & x00F1;o Southern Oscillation and Pacific Decadal Oscillation are analyzed to determine the influence of pacific climatic variability on TWS.