Future Water Availability from Hindukush-Karakoram-Himalaya upper Indus Basin under Conflicting Climate Change Scenarios

Future of the crucial Himalayan water supplies has generally been assessed under the anthropogenic warming, typically consistent amid observations and climate model projections. However, conflicting mid-to-late melt-season cooling within the upper Indus basin (UIB) suggests that the future of its melt-dominated hydrological regime and the subsequent water availability under changing climate has yet been understood only indistinctly. Here, the future water availability from the UIB is presented under both observed and projected—though likely but contrasting—climate change scenarios. Continuation of prevailing climatic changes suggests decreased and delayed glacier melt but increased and early snowmelt, leading to reduction in the overall water availability and profound changes in the overall seasonality of the hydrological regime. Hence, initial increases in the water availability due to enhanced glacier melt under typically projected warmer climates, and then abrupt decrease upon vanishing of the glaciers, as reported earlier, is only true given the UIB starts following uniformly the global warming signal. Such discordant future water availability findings caution the impact assessment communities to consider the relevance of likely (near-future) climate change scenarios—consistent to prevalent climatic change patterns—in order to adequately support the water resource planning in Pakistan.


Introduction
Since the agrarian economy of Pakistan largely depends upon highly seasonal water supplies from the Hindukush-Karakoram-Himalayan (HKH) watersheds of the Indus basin, a pragmatic assessment of future water availability from these watersheds is utmost necessary for ensuring the sustainable socio-economic development in the country.The issue of future water availability from these HKH watersheds has further been accentuated in the wake of recent acceleration of climate change [1] and now established amplified warming over high altitudes [2,3]; thus substantial changes in the existing cryosphere and subsequently in the melt runoff regime are anticipated, threatening the overall hydrological equilibrium upstream and exacerbating the water management problems downstream.It is, therefore, fundamental to investigate the likely impacts of climate change on the future water availability from the HKH watersheds, particularly from the upper Indus basin (UIB-Figure 1) that contributes around half (2405 m 3 •s −1 ) [4] of the surface water availability in Pakistan, warranting a life to the immediate downstream Tarbela reservoir-a largest water storage structure in the country.
Only few studies have been performed in this regard that consistently suggest an increased future water availability from the UIB.For instance, the snow runoff model forced with SRES A2 climate change projections (2071-2100) from the PRECIS (Providing Regional Climates for Impact Studies) regional climate model (RCM) suggests a month-early snowmelt runoff peak, rise in glacier melt and rainfall runoffs and increase in the overall water availability by 7% [5].Similarly, the cryospheric-hydrological model driven with the climate change scenarios of representative concentration pathways (RCPs) 4.5 and 8.5 from the coupled model intercomparison project phase 5 (CMIP5) experiments has simulated an increased water availability at least until mid-21st century due to accelerated melting [6].Forcing the UBC (University of British Columbia) watershed model with the climate change projections from the Conformal-Cubic Atmospheric (CCAM) model and the ICTP regional climate model (RegCM) under RCP4.5 and RCP8.5 scenarios also suggests 50% and 80% increases in the water availability by mid to end of 21st century, respectively [7].Further, the UBC watershed model driven with the projections from PRECIS RCM under A2 and B2 climate change scenarios simulates likewise an increased water availability by around 60% and 35% till 2050s, respectively [8].
The above mentioned studies, and additionally those, focusing on the water availability from the sub-basins of the UIB [9][10][11][12][13][14], obtain the climate change information from the present-day global or regional climate models (GCMs/RCMs) that project a consistent warming for present and future climates under anthropogenic forcing [15,16].In contrast, the UIB has been reported to deciphering climatic changes that are discordant to the anthropogenic warming signals.For instance, observations from the few valley-bottom meteorological stations suggest that the UIB features mid-to-late melt (monsoon) season cooling throughout the second half of the 20th century [17][18][19].Such cooling phenomenon from the valley bottom stations has also been observed over recent three decades   [20][21][22].Lately, a comprehensive and systematic view of the prevailing climatic changes from the novel high-altitude observations (1995)(1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012), together with the updated valley-bottom stations record  confirms that the mid-to-late melt season cooling not only continues but is more profound during recent decades [4].On an extended retrospective timescale, the analysis of tree ring based paeleoclimatic temperature reconstructions also reveals that the Karakoram temperature changes consistently conflicted to the hemispheric temperature trends for past half millennium [23].
Interestingly, the observed hydrocryospheric changes in the UIB are quite consistent with the summer cooling phenomenon, and to somewhat, with the present state of increasing moisture accumulation as analyzed from the station record [4,18,21,22,24].For instance, analyses of daily to weekly MODIS snow cover products for early 21st century indicate an increasing extent of the snow cover over the UIB, over its sub-basins and over the abode Central Karakoram National Park (CKNP) [20,22,25,26].Similarly, analysis of the daily MODIS snow product for the 2001-2012 period reveals a falling tendency of the regional end-of-summer snow line altitudes for the UIB and for its sub-basins, indicating an overall positive change in the high-altitude water resources therein [25].Further, the glaciological studies and the analyses of remotely-sensed data reveal that the overall glacial extent within the UIB and its sub-regions is observed to be somewhat stagnant since 1980s [20,22,[27][28][29][30].In addition to stable glacial extent, balanced or slightly positive geodetic mass balance of the UIB glaciers has also been consistently observed, particularly since the beginning of 21st century [27,[31][32][33][34][35].Subsequently, the hydrometric observations indicate decreasing melt season discharges from the UIB and its sub-basins [4,17,19,36], which is also consistent with the observed cooling phenomenon [4].The continuation of discordant observed cooling for the melt-dominated hydrological regime of the UIB together with subsequent consistent responses from the abode cryosphere implies a broader spectrum of climatic uncertainty over the UIB, of which the climate change projections from the climate models merely serve a warming direction.In view of such observed hydro-cryo-climatic responses, it has been argued that the typical climate change scenarios-representing warming and glacier depletion-are irrelevant for assessing the near-future water availability from the UIB and that such scenarios are subject to a far future, if at all [25].Hence, aforementioned studies assess the future water availability from the UIB [5][6][7][8] or from its sub-basins [9][10][11][12][13][14] by addressing only partially the future climatic uncertainty.
Addressing a broader spectrum of climatic uncertainty over the UIB, this study assesses the future changes in the hydrological regime and the subsequent water availability from the UIB under two contrasting climate change scenarios: (1) one representing the continuation of observed mid-to-late (spring to early) melt season cooling (warming) and stable glacial extent till near-future, and; (2) one suggesting consistent warming as projected by seven fine scale (0.44 • ) CORDEX-SA experiments under the extreme emissions scenario of RCP8.5, together with intact-, 50%-and no-glacier cover.These climate change scenarios are forced to a less-data-extensive semi-distributed watershed model from the University of British Columbia (UBC).A semi-distributed model is advantageous in terms of its higher predictive skill and smaller parametric uncertainty than that of the fully distributed models in a data scarce region, as of the UIB [37].The UBC watershed model [38] has been setup for the UIB up to Besham Qila site with the best possible description of geophysical features, and has been driven for roughly two decades using the most representative high-altitude observations.Since the UBC model distinctly calculates snowmelt, glacial melt and rainfall runoff components, their changes under two likely but contrasting climate change scenarios have been assessed here.Thus, taking into account a wider spectrum of climatic uncertainty over the UIB, the study provides crucial information for the better management of future water resources, and supports in revisiting the near future plans for agriculture production and power generation in Pakistan.

Study Area
The UIB is located within the geographical range of 31-37 • N and 72-82 • E, covering an area of 162,393 km 2 up to the Besham Qila gauging site (Figure 1).The basin hosts three gigantic massifs-western Himalaya at lower latitudes of the basin, Karakoram Range in further north (Trans-Himalayas), and Hindu Kush Range in the west of Karakoram-Himalaya ranges.These ranges collectively host more than 13,000 glaciers covering 19,370 km 2 (roughly 12%) of the basin area where the Karakoram Range hosts most of the ice reserves [31,39].The mean annual snow cover varies between 3%-67% [25].The complex HKH terrain abode, sizeable horizontal and vertical extents and moisture inputs from two distinct large-scale circulation modes-westerly disturbances and summer monsoon-play a significant role in morphing the assorted hydrological regimes of the basin [4].Accordingly to defined hydrological regimes [40], a high flow regime of the UIB (April-September) can further be categorized into mountain nival and mountain glacial melt regimes, which observe their peak discharges during late-spring/early-summer and during late-summer, respectively.The snowmelt dominates in early part of the melt season while the glacier melt dominates for the rest of melt season.During the low flow period (October-March), the runoff comprises of mainly the groundwater component.The overall melt water contribution from the UIB exceeds 70%-80% of its total runoff yield [5,41].The nival regime markedly depends upon the accumulated precipitation during previous season while the concomitant temperatures largely drive the glacial regime [42].These regimes are mainly nourished by the solid moisture input from the western disturbances during winter and spring seasons [31,[43][44][45][46][47].During the summer monsoon period (July-September), the monsoonal incursions also bring considerable moisture at higher (lower) elevations in solid (liquid) form [24,43,44].Such moisture input is anomalously higher during the positive phase of the North Atlantic Oscillation (NAO) and during the warm phase of the El Niño-Southern Oscillation (ENSO) [48,49].

Data Used
The UIB hosts more than twelve high-altitude (above 2200 m asl) automated and six relatively low-altitude (below 2200 m asl) manual valley-bottom weather stations, featuring a medium-to-long term record lengths [4,25].In order to calibrate and validate the UBC watershed model, the meteorological observations from three high-altitude stations, namely Naltar, Hushe and Khunjrab at the elevation of 2898, 3075 and 4440 m asl, respectively, have been obtained from the snow and ice hydrology project (SIHP) of the water and power development authority (WAPDA), Pakistan (Table 1).The meteorological observations from these stations include daily values of temperature maximum (Tmax), temperature minimum (Tmin), and total precipitation (P) for the period 1995-2012.The considered meteorological stations: (1) feature relatively complete set of record for the investigated period; (2) cover possible far east-west extent of the Karakoram cryosphere; (3) spread well across most of the active hydrologic altitudinal range of the UIB (2500-5500 m asl) [41]; (4) and required a relatively less calibration effort.The same set of stations has been used for devising a hypothetical near-future (2013-2030) climate change scenario.While for the far-future climate change scenario (2087-2097), the outputs of seven fine-scale (0.44 • ) CORDEX-SA experiments that are obtained by dynamically downscaling six CMIP5 coupled climate models using three different RCMs (CCAM, RCA4 and REMO) for a high-end emissions scenario of RCP8.5 have been used (Table 2).For comparing simulated discharge with observations, daily stream flow record for the Besham Qila gauge was obtained from the surface water hydrology project (SWHP) of WAPDA for the period 1969-2012.
The watershed description parameters were estimated from the freely available databases.For instance, a lately available one arc-second (30 m) global digital elevation model (DEM) of the shuttle radar topographic mission (SRTM) was obtained from the USGS database.Similarly, the forest cover and forest canopy density datasets were obtained from the United Nations Environment Programme (UNEP) Environmental Data Explorer, as compiled from USGS/NASA-Distributed Active Archive Center [50].The glacier cover was extracted from the lately released Version 5.0 of the Randolf Glacier Inventory (RGI5) [39].

Hydrological Model
Originally developed for the mountainous watersheds around 3-4 decades ago and being continuously updated since then, the UBC is a conceptual, continuous, deterministic, and semi-distributed energy-balance model [38].In view of scarce data availability within the mountainous watersheds, as in case of the present study area, the UBC takes minimum set of point observations (Tmax, Tmin and P) together with geophysical watershed characteristics in order to reproduce the hydrological response at hourly and daily time scales.Since the vertical gradients of hydroclimatic parameters are key factors in determining the hydrological behavior of a watershed, model divides the watershed possibly into twelve elevation zones, within which all hydro-meteo-climatic conditions are assumed to be homogeneous.The UBC implements a simplified energy balance approach for the snowmelt.In this regard, the requisite parameters, such as, cloud cover, wind speed and albedo are estimated from the temperature record based upon the relationships established from the field investigations performed within the Karakoram-Himalaya region under the SIHP, WAPDA project [51].This specialty of the UBC model along with the afore-mentioned features makes its application highly suitable for the Himalayan environment.
The model melts the accumulated snowpack and the existing glaciers based on available energy, and provides discrete runoff contributions from snowmelt, glacier melt and rainfall.The generated runoff is subsequently allocated to fast, medium, slow and very slow runoff components by the soil moisture routine (Figure 2).For instance, after accounting for interception and sublimation losses, an amount of water reaching the surface is divided into fast (surface) runoff and infiltrated water, based upon the impermeable area that varies according to saturation conditions.Considering soil moisture deficit under evaporative demands, the infiltrated water is further divided into groundwater and interflow (medium runoff) components.The groundwater is allocated to upper (slow runoff) and deep (very slow runoff) zones based on a predefined percolation coefficient.The model applies a water-conserved linear reservoir approach for routing the runoff from each component towards the watershed outlet.The snow accumulation and melt routines are described in Appendix whereas for details about other components of the UBC model, earlier studies [51][52][53] are referred.

Model Setup
The UIB drainage area has been delineated up to the Besham Qila site by applying a GIS-based automated delineation procedure on a 30 m SRTM DEM.The automated delineation procedure erroneously includes the adjoining northeast watershed of Pangong Tso, which is in fact a closed basin as suggested by now substantial geological evidence [54,55].Excluding the Pangong Tso watershed, the UIB drainage area was estimated as 165,515 km 2 , which is quite in line with the SWHP, WAPDA estimates (i.e., 162,393 km 2 ).The details of actual delineation of the UIB considered here along with the implications of its overestimated drainage area on modeling runoff are given in [56].
In view of a large vertical extent of the UIB (594-8550 m asl), the basin was divided into twelve maximum possible elevation zones (bands) based on the area elevation curve (AEC) that was developed from the 30 m SRTM DEM. Figure 3 shows the demarcation of 12 elevation bands on AEC, their glacial coverage, and altitudes of meteorological stations used for model calibration.For each elevation band, the geophysical parameters, such as, total surface and glaciated areas along with their north-south orientations (0-1), forest cover and canopy density as well as mean elevations are estimated from the respective datasets (Table 3).
The UBC model allows adjustments to more than 90 parameters, of which usually 17 parameters of precipitation distribution and its adjustments, and of water allocation and its routing need to be calibrated while the rest take preset standard values.For the precipitation distribution within the UIB, a strong positive vertical gradient has been suggested by various studies [13,43,57,58].Instead, station observations suggest a negative precipitation gradient along an elevation increase.In fact deriving vertical precipitation lapse rate from a very sparse meteorological network within the complex terrain of HKH is extremely difficult, particularly when basin areas along the latitudinal and longitudinal transects come under the influence of either completely distinct or combination of active large scale precipitation regimes.Few field campaigns have suggested that the maximum precipitation within the UIB occurs roughly between 5000 and 6000 m asl [27,31,43,[58][59][60][61], where it becomes five to ten fold [27,61] and then decreases above such elevation range.Since, no long-term record is available above 4400 m asl, so far, observed precipitation datasets are neither fully representative of the basin discharge nor the glacial nourishment [10,42,62,63].Based on glacier mass balance proxy, a first-order approximation of the vertical precipitation lapse rate for the western Karakoram has been reported to be around 0.21% ± 12% at the height of 5500 m asl [63].Knowing the facts that large number of the UIB glaciers lying within the high-relief HKH ranges are mainly nourished by avalanche activities [31] and that the lapse rate is estimated only for a part of the UIB [63], there exists a huge uncertainty of precipitation distribution.Since the UBC allows the implementation of multiple elevation barriers in order to better describe the orographic precipitation variability, the height of such barriers and their vertical precipitation lapse rates were set by generalizing the available findings [27,61,63].The first barrier, E0LMID was set at 3000 m asl while the second barrier, E0LH I was set at 6000 m asl.The range between these two barriers covers most of the vertical extent of the glaciated area (Figure 3).The vertical precipitation lapse rate was set to 0.10% [27,61] for the elevation range below E0LMID height, and to 0.31% and −0.40% for up to and above E0LH I heights, respectively [63].
The debris-covered glaciers are very common within the UIB [64,65].Modelling the ablation from such debris-covered ice requires mapping of its extent, estimation of debris thickness and consideration of debris material composition [57]; former two are mainly considered in the modeling process.Presently, these estimates are not available for the whole UIB and estimating them from the remotely sensed datasets at such a large scale is nontrivial [66].Investigating glacier area of 3500 km 2 within the Hindukush-Karakoram region, a study reports roughly 20% of such area under the debris cover [28].Investigation over more than 700 glaciers within the CKNP suggests that around 11% of the total glacier area (around 4600 km 2 ) is under the debris cover [67].For the whole Hindukush-Karakoram, the supraglacial debris cover has generally been reported to be around 10% of the total glacier area [65].Most recent estimates suggest that around 2000 km 2 of the total glacial extent of the whole Indus basin is under the debris cover [68].Such estimates are in line with earlier studies [65,67] given the fact that debris-covered glaciers of the Indus basin are mainly present within the UIB.Therefore, the supraglacial debris cover from the studies [65,67,68] have been considered as reasonable estimates for the UIB.Such debris-covered glacial extent appears to be stable over time as either no or near-zero change has been reported for 93 Karakoram glaciers investigated over the 1977-2014 period [69].Further, the analysis over more than 700 CKNP glaciers since early 21st century also suggests insignificant changes in the supraglacial debris extent [20,22].
The supraglacial debris have a variable thickness, which substantially influences the rate of ablation from the underlying ice [64]; the ablation is typically highest at around 1 cm (lowest above 40 cm) of the debris thickness as suggested by various studies within the UIB (Rakhiot glacier [70], Rakhiot and Barpu glaciers [71], Baltoro glacier [72], Hinarche glacier [57], Dokriani Glacier [64]).An increase in the ablation rate for up to 1 cm of debris thickness is due to the conduction of absorbed radiative heating to underlying ice, while its decrease for further increase in debris thickness is due to the insulation effects; for the average debris thickness, the overall rate of debris-covered ice ablation typically remains around 0.5 ± 0.1 of the clean ice ablation rate [57,64,[70][71][72].Modeling the ablation from debris-free and debris-covered glaciers over wider area of around 4600 km 2 within the CKNP has also suggested, on average, similar rates of debris-cover ice ablation relative to clean ice ablation [67].Though the UBC model does not explicitly simulate the debris-covered ice ablation processes, however, their generalized implementation can be achieved through alternate mechanisms.For instance, the topographic parameters, such as orientation and slope affect the total ablation in a similar way the debris thickness does [72].However, quantitative estimates of the influences of slope and orientation on the ablation rates are not available for the study basin.On the other hand, reducing the extent of debris-covered part of the glaciers proportional to the factor of their ablation reduction (relative to clean ice glaciers) is another alternate to implement the debris-covered ice ablation in the models, in view of their typical surface melting property.Thus, the debris-covered glacial extent of the UIB that is around 1937 km 2 (10% of the total UIB glacier cover) has been reduced by a typical factor of 0.5 from the lower elevation bands, assuming the average thickness of the debris [57,64,[70][71][72].Such a generalized approach for implementing the debris-covered ice ablation is quite feasible for a semi-distributed watershed model.

Calibration and Validation
A semi-automated optimization procedure together with the manual tuning was adopted to obtain the optimized values of calibration parameters within their predefined ranges.For sound evaluation of calibration and validation of the model simulations, the use of several efficiency criteria has been recommended [73].In this regard, an objective function has been employed to assess the optimization (Equation ( 1)) together with four different shape and scale indices, such as, Nash-Sutcliffe Efficiency (NSE) [74], percentage volume error (%VE), root mean square error (RMSE) and weighted coefficient of determination (D).The NSE (Equation ( 2)) is a ratio of mean squared error (MSE) to the observed variance, subtracted from the unity.Ranging between −∞ (worst fit) and 1 (perfect fit), the NSE determines the scale and shape fit of the observed and simulated discharges.The %VE (Equation ( 3)) determines the scale fit, denoting a zero value as a perfect fit while positive and negative values as over-and under-estimations of the observed discharge, respectively.The RMSE (Equation ( 4)) with an optimal value of zero provides the residual between observed and simulated discharges.With an optimal value of 1, D (Equation 5) is a shape parameter, which determines the relationship between the simulated and the observed discharges.The simple D should not be used for error quantification since, relying only on the correlation, it can substantially overestimate the model efficiency [73].Here, a weighted form of D was used that takes into account the slope, b, and intercept, a, and performs much better than D alone [73].
where, Vsim and Vobs refer to simulated and observed flow volumes, while Qobs i and Qsim i are the observed and simulated discharges on day i, respectively.Qobs refers to the mean observed discharge for the simulation period containing n number of days.

Prevailing Hydro-Cryo-Climatic Changes
As mentioned afore in Section 1, there is now a substantial evidence that the UIB responses to the anthropogenic forcing are partly contrary to what has been observed over and across the neighboring regions and worldwide, and as projected by GCMs/RCMs under various emission scenarios.The rest observed changes that are consistent with a typical climate change signal mainly include a dominant spring warming and somewhat increase in winter precipitation [4,[17][18][19].Nevertheless, a summer cooling phenomenon is unique for the UIB with the present-day climate models being unable to represent it even at a qualitative scale [16].In Tabular Figure 4, yet a complete picture of prevailing climatic changes is presented from all operational high-and low-altitude stations within the UIB after [4].These changes exhibit significant cooling during the monsoon months, particularly in Tmax, where insignificant tendencies are also well agreed on a sign of change.For precipitation, most of the stations suggest an increase during the late summer monsoon (September) and winter (December to February) months.Hushe --0.5 --0.1 0.9 0.0 1.7 --0.6 --0.9 0.2 --2.0 --0.9 0.

Stations Jan Feb Mar Apr May Jun Jul Aug Sep Oct
Shigar --0.4 --0.8 --0.2 --0.8 --3.8 --1.5 --0.8 0.3 --0.1 --0.9 1.1 0.1 --0.2 --0.9 --0.9 --0.Bearing in mind that the above-mentioned cooling is observed from a sparse station network, which is further confined to the political boundary of Pakistan, such cooling may only sporadically exist across the basin.Against this background, trends in discharge-that exhibit integrated response to the prevailing climatic changes over a basin-as well as changes in the glacial mass balances can supplement the sparse observational network and can serve as proxy observations for the regions of no data record, such as, eastern Karakoram [4].A few studies analyzing these proxy observations from eastern Karakoram suggest against the observed cooling reports.For instance, arguing that a rise of late melt-season discharge from the eastern Karakoram over the 1973-2010 period is subject to enhanced glacial melt, the study [75] linked its speculation to a significant glacial mass loss of −0.10 m•yr −1 observed over a short period of 2003-2008 by [35].In contrast, the geodetic mass balance analysis over the same study zone but over a longer period (2000-2010) revealed a slightly positive mass balance for the eastern Karakoram glaciers [33].Furthermore, the areal extent of the eastern Karakoram glaciers has been reported to be growing since the end of 1980s [29].A significant falling trend of the end-of-summer regional snow line altitude for the eastern Karakoram over the period 2000-2012 also indicates a positive change in the high-altitude frozen water resources therein [25].For other main cryospheric regions of the UIB, the studies consistently report and support cooling phenomenon [4,17,19,21,27,28,30,35,36,75,76].Over the whole UIB, increasing snow cover in a recent decade [25,26] and conflict of tree-ring based paeleoclimatic Karakoram summer temperatures with the hemispheric warming for recent five centuries [23] have also been reported.
Though the aforementioned findings support a wide spread cooling over the UIB, the UIB discharge tendencies have been additionally assessed, in order to gain insight that whether such cooling phenomenon predominantly influences the overall hydrological regime, aside few exceptions might exist.For this, the long-term hydrograph along with the true (Sen's [77,78]) slopes of the long term trends for each day of a year (DOY), smoothed with 3-10 days moving average windows are plotted (Figure 5).It is noted that the UIB discharge tendencies are well explained by the prevailing climatic trends shown in Tabular Figure 4.For instance, cooling within the monsoon season (July-September) that mainly coincides with the glacier melt season, has resulted in a significant drop of flows during late-July to early-August period.Similarly, a significant rise (fall) in discharge for month of May (June) that coincides with the main snowmelt season, clearly indicate early shift of the snowmelt regime.These findings confirm the dominant influence of cooling on the glacial regimes of the UIB.
The observed cooling is attributed to the incursions of the summer monsoonal offshoots, which are seemed to be more frequent in near future, in view of the observed rising trends in the monsoon season cloud cover, [21] and precipitation [4] and due to projected extension in the monsoonal domain north and northwest ward over the UIB [47].Such obvious indications of prevailing hydroclimatic changes and their continuation emphasize the use of a plausible climate change scenario, at least for the near-future.Earlier, it has been argued that the climate change scenarios (temperature increase, glacier area decrease) adopted by previous studies are not relevant at least for the near-future and may only be realized in the far-future with no time stamp attached [25].Hence, in order to address larger spectrum of future climatic uncertainty, two plausible hypothetical climate change scenarios have been considered; one representing the continuation of prevailing contrasting climatic conditions while other consistent with the global warming patterns.

Near-Future Climate Change Scenario
Here, a climate change scenario has been empirically devised to represent the continuation of prevailing climatic conditions to the near-future.Contrary to typically devised hypothetical climate change scenarios that often perturb the climatic variables arbitrarily and uniformly through the annual course [44], here emphasis is given on a scenario that anticipates climatic changes on an intraseasonal scale, as also depicted from the Tabular Figure 4.
Thus, the monthly time series of considered high-altitude stations (1995-2012) have been extended for next 18-years (2013-2030) using a linear regression intercept and true slope obtained through Sen's slope method [77,78].The monthly change factors (Table 4) for the extended future period (2013-2030) were estimated relative to the observed period (1995-2012).A stochastic weather generator of LARS-WG [79] has been calibrated using the observed time series (1995-2012) from each station separately.Based on the calibrated observed statistics for the period (1995-2012), several synthetic ensemble members each of 30-year daily time series were stochastically generated.Then, the LARS-WG was validated by comparing the statistics of these stochastically generated synthetic ensemble members with those of the observed times series, which should not be statistically significantly different in an optimum case.The observed statistics for the period 1995-2012 were then perturbed by monthly change factors of the future period (2013-2030) in a validated LARS-WG and ten future ensemble members each of 30-year daily time series were generated.It is to note that the variability of future time series is different than that of observed time series, as it has been stochastically perturbed.
The future climatic changes are likely to change the extent of future cryosphere [10], which in turn can affect the future water availability [9,11].Thus, it is important to incorporate changes in the future cryosphere that are plausible under anticipated climate change scenario, while assessing the future water availability from a glacerized basin.Since the overall glacial extent of the UIB somewhat remained stagnant since 1980s [28][29][30]65] and particularly since early 21st century [20,22] under the ongoing regional warming, no change in the present-day glacial extent has been assumed for the near-future scenario (2013-2030), as also adopted by several studies [5,9,11].

Far-Future Climate Change Scenario
The climate change information have been obtained from seven fine scale (0.44 • ) CORDEX-SA simulations, which are performed by downscaling six CMIP5 GCMs employing three different RCMs (Table 2).In view of obtaining a far-future climate change scenario, projections for an extreme emissions scenario of RCP 8.5 are considered for the last decade of the 21st century.The RCP8.5 scenario, assuming ramp up of the radiative forcing up to 8.5 W•m −2 by the end of 21st century, is considered as an upper bracket of the climatic forcing [80,81].
It is pertinent to mention here that though CORDEX-SA experiments provide high resolution present and future climate simulations, however, their ability to reproduce the observed hydro-climatology over the steep HKH topography is by hardly any means better than their forcing CMIP5 GCMs, which being consistent to their older versions and downscaled experiments, substantially over-(under-) estimate precipitation (temperatures) and fail to reproduce summer cooling [7,15,16,[82][83][84][85][86][87].Such cold and wet biases in the CORDEX-SA and their CMIP5 forcing experiments are larger than their projected end of 21st century climatic changes under RCP8.5, indicating huge climatic uncertainty over the UIB [16].These projected warmer climates from the low fidelity CORDEX-SA experiments are considered here as a hypothetical warming scenario-assuming that the climate models have simulated these climates in a physically consistent way [88].Further, it is not realistic to attach a time stamp to this scenario, as it is not clear when the UIB will uniformly follow the anthropogenic warming.The CORDEX-SA experiments were adjusted for biases against the considered station datasets using a quantile-quantile mapping approach, prior to their use in the hydrological model.
For changes in the glacial extent, various studies have generally adopted an arbitrary reduction of 50% by mid-century under projected wetter and warmer climates [5,6,9,13].In contrast, employing a glacio-hydrological model over highly glacierzied (Shigar) sub-basin of the UIB, recent study has reported a negligible ice cover changes by the mid-century under warmer climates as projected under various RCP scenarios [14].Hence, in view of somewhat stagnant UIB glacial extent at present and conducive hydroclimatic conditions under ongoing regional and global warming, reduction in the future glacial extent under more wetter and warmer climates by the mid of 21st century needs to be rationalized for the UIB.Further, despite 50% reduction in the glacial extent by the mid-century, studies project consistent increase in the water availability due to enhanced glacier melt [5,6,13].In fact, glacier melt contribution and subsequent water availability will consistently rise under warmer climates until unless the glacial cover declines to a tipping point, as clearly shown by [14].Identifying timings of such tipping point glacial extent, however, is non-trivial, first in view of the huge multi-model spread of the projected warming [14,16], and second, the timings of occurrence of such projected warming is less likely to be realistic, and rather unknown, in view of conflicting cooling observed at present [4,17,20,21].Hence, for assessing the far-future water availability under the extreme-end forcing of RCP8.5 by the end of 21st century, three glacial extents are assumed regardless of their time of realization: (1) intact-glacier cover; (2) 50%-glacier cover, and; (3) no-glacier, at all.The 50% reduction in the glacier cover was obtained from the lowest elevation bands.

UBC Model
The calibration and validation results of the UBC watershed model are shown in Figure 6 and are quantitatively summarized in Table 5.The overall calibration results have suggested the estimates of NSE, D and EOPT as 0.9 with a negligible volume error (−0.01) and a RMSE of around 800 m 3 •s −1 .Within a 9-year calibration period (1994/1995-2002/2003), the results suggest maximum calibration efficiency for the year 2002-2003 where NSE and D are 0.96, and EOPT is 0.93.On the other hand, lowest performance was achieved for the flood year 1998-1999 (NSE of 0.82, D of 0.91 and EOPT of 0.60) followed by the flood year 1999-2000 (NSE of 0.81, D of 0.84 and EOPT of 0.67).These two consecutive years have experienced flood peaks of above than 13,000 m 3 •s −1 owing to intense July rainfall (high correlation with rainfall outflow), partially captured by the station observations but dampened during the calibration process in order to ensure homogeneous distribution across the large-extent study basin.The long-term flow duration curve and hydrograph shown in Figure 5 suggest that discharge peaks over and above 10,000 m 3 •s −1 are extremely rare for the UIB.It is pertinent to mention here that NSE is particularly sensitive to such peaks as D [73], thus EOPT too.Such situations may additionally happen within the large river basins, when the rainstorm only partially affects the basin in a way that it is not captured by the sparse observational network.Such phenomenon can further be illustrated from the validation statistics.The validation period (2003/2004-2011/2012) comprises of three extreme flood years in a series where the flood of year 2009-2010 was unprecedented.The 2010 flood event had partially influenced the study basin, yielding two peaked discharges (above 18,000 m 3 •s −1 ) at the basin outlet in August.The overall validation results suggest NSE, D, EOPT, %VE as 0.85, 0.83, 0.82 and −2.24, respectively.Excluding three exceptional years of partially affecting rain-fed flood events, it is noted that the overall validation results into higher NSE and D of 0.9 and 0.92, respectively, implying again that the model performance for these flood years is mainly restricted by an incomplete spatial rainfall distribution captured by the observational network.Notwithstanding, the overall performance of the UBC model for the period (1994/1995-2011/2012) suggests NSE, D, EOPT and %VE as 0.87, 0.86, 0.86, −1.15, respectively when including the flood years, and as 0.92, 0.86, 0.89, −2.82, respectively without considering the flood years.Such performance statistics suggest overall a successful application of the UBC model over the UIB.(d,e) for the validation period only; (f,g) for the whole period.

LARS-WG
The validation results of the LARS-WG are shown in Figure 7 and summerzied in Table 6.It is noted that LARS-WG performs well in reproducing the observed statistics despite the fact that the observed time series spans over a short period of 18-years (1995-2012).The synthetic temperature time series generated by LARS-WG features a mean annual bias ranging from 0.01 to −0.16 • C for the Tmin and −0.05 to 0.02 • C for the Tmax amid the considered meteorological stations.The absolute RMSE ranges from 0.146 to 0.340 • C for the Tmin and from 0.213 to 0.243 • C for the Tmax, where such RMSE range is below 0.02 for Hushe, 0.03 for Khunjrab and 0.04 for Naltar stations in the relative terms.As for the precipitation, the mean annual bias (−0.53 to 0.6) for the LARS-WG generated synthetic time series is negligible while the absolute RMSE is only 1.274 mm for Khunjrab, 3.642 mm for Hushe and 5.481 mm for Naltar stations.It is noted that the performance of LARS-WG is highest for Khunrab station followed by Hushe and Naltar stations, respectively.Similarly, the fractional positive or negative biases in the synthetic temperature time series suggest that the LARS-WG features higher skill for temperature variables as compared to precipitation, which is further higher for the Tmax than for the Tmin variables.Overall, the performance of the LARS-WG has been found quite satisfactory for the considered high-altitude stations within the study region.Such satisfactory performance of the LARS-WG has allowed its further application in devising a plausible near-future climate change scenario.The Figure 7 also shows the near-future climate changes scenario (Tmax, Tmin and P), which is devised by perturbing the statistics of synthetic validated LARS-WG time series with the change factors derived for the 2013-2030 future period relative to the 1995-2012 observed period (Table 4).

Future Water Availability
Results of both near-and far-future climate change scenario are given in Table 7, and are graphically shown in Figures 8 and 9, respectively.

Near-Future Climate Change Scenarios
Under the most plausible near-future scenario, a median ensemble change suggests that the mean UIB discharge will be reduced by roughly 8% (range 6%-10%) relative to the present-day mean discharge.Such reduction is a subsequent response from the individual discharge components, which exhibit distinct changes.For instance, the glacier melt contribution to the total flow will be reduced by 24% (range 21%-25%), indicating marked changes in the overall hydrological regime of the UIB.In contrast, snowmelt outflow has been projected to rise by 7% (range 2%-14%).Changes in other discharge components suggest 29% decrease in the rainfall outflow (range 0%-43%) while groundwater contribution will increase by roughly 2% (range 0%-3%).
Changes in the individual discharge components not only influence the amount of water availability but its timings too.The Figure 8 shows that increase in the snowmelt outflow components is accompanied by an early snowmelt, while reduction in the glacier melt contribution has resulted in its delayed peak by one month.The snowmelt peak has been shifted from June to May while the peak glacier melt is shifted from July to August.Subsequently, the mean annual hydrograph suggests an early onset due to early snowmelt while the peak annual discharge, coinciding with the glacier melt timings, is delayed by one month, shifting from July to August.

Far-Future Climate Change Scenario (Intact Glacier)
Given that the far-future climate change scenario is realized at some point in time and the present-day glacial extent remains unchanged by then (Figure 9), a substantial increase in the UIB discharge is projected; median of such an increase is 1.6 times (range 1 to 2.2 times) of the base period discharge.The maximum change is projected for the rainfall outflow component, which will increase by 3.4 times (range 2 to 5 times).This is followed by a threefold increase in the glacier melt contribution (range 2 to 4 times).The minimum rise is observed in the snowmelt contribution, which will be around 25% (range −2% to 74%) of the base period snowmelt contribution.The groundwater contribution will subsequently rise by 86% (range 56% to 100%).
Moreover, the timings of peak snowmelt runoff will advance by one month and will occur in May rather than in June.Similar advance is observed for the whole nival regime, which will span over March-July rather than April-August period.The spread of rainfall outflow is relatively large among the ensemble members.However, the median rainfall outflow will feature its peak in May instead of August, followed by in October.Interestingly, the glacial regime exhibits no shift but will expand from May-October to April-November, owing to an early onset and a delayed retreat while months of peak timings will remain unchanged.Similar is the case with the annual hydrograph, where the wet season will span over March-November period.Besides 50% reduction in the glacier cover, the median of the ensemble simulations suggests that the overall UIB discharge will increase by 24% (range −4% to 54%) relative to the base period.Interestingly, at the lower-end of the climatic uncertainty, the overall discharge will become slightly lower than that of the base period discharge.Similarly, though the glacier melt contribution to the overall discharge will rise by 50% and 100% for the median and maximum changes, respectively, however it will be slightly lower to that of the present-day for the lower-end of the climatic uncertainty.Such finding suggests that decreasing the glacier cover by more than 50% may yield the water availability lower than the base period, indicating a tipping point behavior of the glaciers under extreme forcing.The maximum change is again projected for the rainfall outflow component, which will increase by 2.8 times (range 1.5 to 4.5 times).The snowmelt contribution will rise by around 25% (range 1% to 62%) while the groundwater contribution will drop by 20% (range −5% to −25%) relative to the base period.
The timings of snowmelt outflows under 50%-glacier realization will be similar to those as projected under intact-glacier realization.Similarly, the median rainfall outflow will also feature its peak in May instead of August, followed by in October, where the peak in October is smaller in magnitude as compared to the intact-glacier scenario.Interestingly, the glacial melt regime agrees well with the glacier melt regime of the base period, exhibiting only minute extension in the late-glacier melt period.The groundwater contribution has been considerably dropped during the glacier melt period (June-September) while it will be maximum in May and June.The annual hydrograph will mostly remain close to the annual hydrograph of base period with a little extension owing to an early shift of the snowmelt regime, where the wet season will span over the March-November period.The maximum change is projected in the rainfall outflow, which will rise by 2.3 times (ranging from being double to 3.6 times) of its present-day contribution.After the rainfall outflow, the snowmelt outflow will rise by 19% as compared to its base period contribution.In contrast, a negative change of around 45% (range 29% to 48%) is projected for the groundwater component as there will be no groundwater recharge from the glacier melt throughout the glacier melt season.Besides positive changes in the rainfall and snowmelt outflows, the overall water availability from the UIB will be reduced by 40% (range 22% to 51%) relative to the present-day water availability, mainly due to absence of glacier and glacier melt.
The timings of the rainfall and snowmelt outflows under no-glacier realization will be similar to those as projected under intact-glacier and 50%-glacier realization.The peak groundwater period (June-October) for the intact-glacier realization will now confine to May-June only, as in case of 50%-glacier realization.The overall hydrological regime of the UIB will be shifted from glacial-nival to nival-pluvial regime.Hence, a peak discharge will be available in May, while wet (dry) season will span from March to July/August (September to February).

Discussion
The model performance over 9-year calibration and 9-year validation periods has been found very good, implying that model can reproduce the observed runoff with required accuracy using a minimum set of observations.However, limited skill of the UBC model in reproducing the rare extreme discharge peaks is mainly subject to extreme rainfall events partially affecting such a large study basin.This emphasizes either on the application of UBC model for the sub-basins of the UIB or employing a fully distributed hydrological model with a full set of available meteorological stations.Moreover, exceptions to the observed cooling at some places cannot be characterized in the UBC as it considers the hydroclimatic conditions within each elevation zone homogeneous, reinforcing the application of a fully distributed model.These models further allow spatially explicit distribution of snow pack and existing glacial extent and masses, which are primitive to reduce the typical uncertainties of compensation between the glacier and snowmelt.Further, the hydrological model calibration/validation against the multiple dataset of observed snow cover, glacial mass balances and discharges together improve the overall runoff estimation, and yields realistic snow and glacier melts components [89].This will be the focus of a follow up study.
The study aimed to address the climate change impacts on the hydrological regime of the UIB considering a broader climatic uncertainty; relevant climate change scenarios derived from the observational record, and those, projected by the low fidelity climate model experiments under extreme forcing scenarios.It is worth mentioning that the single scenario for the near-future may represent only an aspect of likely changes in the near future climate.Therefore, a wide range of plausible climate change scenarios that are devised empirically on similar lines-representative of the continuation of prevailing climatic trends-should be adopted in order to more adequately address the near-future climatic uncertainty, and subsequently, the water availability from the UIB to support the informed policy decisions.
The study has focused on changes in the future water availability and its temporal availability (seasonality) rather than changes in the mass balance of the glaciers.The results from the near-future scenario suggests the continuation of the observed hydroclimatic changes, and that, such changes will largely influence the melt dominated hydrological regime of the UIB.Under such scenario, the projected reduction of glacier melt by a quarter, and of rainfall runoff by half are due to the mid-to-late melt season cooling while a projected increase in the snowmelt up to 14% is a result of early-melt season warming.Subsequently, the overall water availability is projected to decline by up to 10%.Worth noting fact is that substantial reduction in the glacier melt, aside changes in precipitation that are also positive, indicates a positive change in the frozen water resources of the UIB.Such finding is largely consistent with the prevailing changes depicted from the long-term observed streamflow record (Figure 5) and with the reports of anomalously positive or stable glacial mass balances and extents [20,22,[27][28][29]31,[33][34][35], increasing snow covers and falling end-of-summer regional snow line altitudes [20,22,25,26], and subsequent reduction in the melt season discharges [4,19].On the other hand, a decrease in the glacier melt and in the overall discharge is in direct contrast to the reports from earlier studies, which addressing only partially the climatic uncertainty over the UIB assess its future water availability under warmer climates [5][6][7][8]90], and suggest rather enhanced glacial melt and rise in the water availability throughout the 21st century.Such contradictory findings of rise in the near-future water availability, lacking the relevance of near-future climate change in accordance with the prevailing hydroclimatic trends, are quite misleading for the immediate freshwater resources management in Pakistan.
Though abrupt drastic decline in the near-future runoff is less likely, considerable changes in its seasonality are rather obvious from the results.For instance, increasing snowmelt owing to early melt-season warming has resulted in a month early peak of the snowmelt runoff regime.In contrast, decreasing glacier melt owing to mid-to-late melt-season cooling has resulted in a delayed peak of the glacier melt runoff regime by one month.Subsequently, the overall hydrological regime exhibits similar changes as of the dominating glacier melt regime.Moreover, if the observed changes will continue roughly in the similar manner, the snow and glacier melt regimes seem to be shifting apart from each other, indicating more dependence on nival to pluvial rather than on glacial regime for the future water availability.These findings clearly suggest that the non-uniform climatic changes on an intra-seasonal scale have distinct impact on the discharge components.Hence, uniformly perturbing various input parameters throughout the annual course (which yield physically inconsistent changes) as typically adopted in the hypothetically devised climate change scenarios [44,91] will project superficial changes in the timings of the runoff components and in the overall seasonality of the runoff regime.
For the far-future scenario, an increase in the rainfall outflow under the intact-glacier realization is projected to be lower under the 50%-glacier realization, and further lower under the no-glacier realization.This is due to the fact that reduction in the glacier cover will expose more surface area to the soil layer infiltration and ground water percolation throughout the melt period and mainly during the glacier melt season.As a result, the direct rainfall outflow will be delayed due to infiltration and percolation to the groundwater and will become part of the groundwater contribution.Such relative reduction in the rainfall outflow during the main glacier melt season (July-September) can be seen from Figure 9. Similar case is observed for the snowmelt outflow increase, which being 24% and 19% for the 50%-and no-glacier realizations, respectively, will remain slightly lower than the increase for the intact-glacier realization.It is noted that increases in the rainfall outflows under no-glacier, 50%-glacier and intact-glacier realizations of the far-future climate change scenario are projected to be higher than the projected changes in the rest of hydrological components, which may partly be due to: (1) increase in the actual rainfall itself, and; (2) phase changes of the solid precipitation into liquid under the projected increase in the early-melt (spring) season temperatures (Figure 4), as the maximum increase in the rainfall outflow occurs during the main snowmelt period (Figure 9).
Besides 50% reduction in the glacier cover under the far-future scenario, the median change suggests that the glacier melt contribution and the overall water availability from the UIB will rise, consistent with earlier studies [5,6,13].However, for the lower-end of the climatic uncertainty under the far-future scenario, the glacier melt contribution and the overall water availability from the UIB is projected to become slightly lower than that of the base period, implying that a further decrease in the glacier extent (by more than 50%) may yield the water availability lower than the base period.Similar findings have been reported for the Shigar sub-basin of the UIB, for which a reduction in the glacial extent by half or more projects a reduction in the future water availability under SRES A2 scenario, relative to the base period [11].These contrasting findings of the future water availability increase or decrease with 50% reduction in the glacier extent mainly result from the future climatic uncertainty for the given forcing; the climate models project distinct climatic changes under the same forcing scenario owing to their structural uncertainty [16,47].
Summarizing, the far-future scenario suggests the overall dominance of the glacier melt regime owing to enhanced glacier melt under warmer climates, and subsequently, rise in the overall water availability until glaciers vanish (or reach a tipping point) where their contribution becomes so small that the overall water availability will reduce.These findings are largely in agreement with the earlier future water availability projections for the UIB under warmer climates as typically but inconsistently projected by the climate models under a variety of anthropogenic forcing [5][6][7][8], as explained in Section 1.The results of far-future water availability are further consistent with earlier studies focusing on the sub-basins of the UIB.For instance, a consistent increase in the mid-21st century water availability from the Shigar sub-basin has been reported under SRES A2 scenario due to enhanced glacier melt until the glacial extent reduces to 50% [11].Similar results of an increased future water availability due to enhanced glacier melt are reported for the Shigar sub-basin throughout the 21st century under various RCP scenarios as simulated by CMIP5 GCMs [14].The assessment of future water availability from the Hunza, Astore and Gilgit sub-basins of the UIB under SRES A2 scenario likewise suggests an increased water availability for the intact-glacier but drastic decrease for the no-glacier realizations by the end of 21st century.Further, a twofold increase in the future water availability from the Hunza sub-basin in response to the hypothetical warmer climates till the end of 21st century has also been reported [10].Nevertheless, projections of first increasing flows due to deglaciation and then abrupt drop on vanishing of glaciers and glacial melt (or reduction due to reaching tipping point) has most likely been delayed until the region comes uniformly under the influence of warming signal and starts following the same sign of change as projected by the present-day climate models under various emission scenarios.
Interestingly, prevailing observed changes and climate models' projections under anthropogenic forcing both consistently suggest a warming during the early-melt (spring) season in the future [16].Hence, all four realizations of the adopted climate change scenarios project similar changes for the nival regime of the UIB, suggesting a month early (in May instead of June) snowmelt outflow (peak) and its overall increase (7% to 25%) by near-to-far future, which is further consistent with the earlier reports.For instance, a month early snowmelt runoff peak has been reported for the UIB till the end of 21st century under SRES A2 emissions scenarios [5].Similar is the case with the Shigar sub-basin of the UIB, for which an earlier occurrence of the snowmelt outflow has been projected under various RCP scenarios throughout the 21st century [14].As for changes in the overall future water availability, a qualitative agreement is achieved amid the near-future scenario with intact-glacier and the no-glacier realization of the far-future scenario (reduction up to 10% and 50%, respectively), which emphasizes more on the importance of unique cooling phenomenon and its relevance than on changes in the cryosphere, while assessing the future water availability from the UIB.

Conclusions
The near-and far-future climate change scenarios that contradict on the direction of their projected climatic changes during the main glacier melt season suggest distinct changes in the overall hydrological regime of the UIB, and subsequently, in the future water availability from the basin.For instance, the near-future scenario suggests a substantial drop in the glacier melt and one month delay in its peak, and subsequently, reduction in the overall water availability.These changes, particularly in the view of increasing precipitation, indicate a positive change in the frozen water resources of the UIB.Further, strengthening of nival but suppression of the glacial melt regime under the near-future scenario implies both the regimes shifting apart, indicating substantial changes in the seasonality of the overall hydrological regime.In contrast to the near-future scenario, the far-future scenario project a strengthening of the glacial regime, no change in the glacier melt peak runoff timings and substantial increase in the overall water availability till the glacier cover either reaches a tipping point or fully vanishes.For the nival regime, however, both scenarios consistently suggest its strengthening that is accompanied by a month early occurrence of the snowmelt outflow.
In view of the huge adaptation costs involved, the contradictory future water availability projections for the UIB emphasize on considering the relevance of likely climate change scenarios for the foreseen future.Further, regardless of which scenario is realized, distinct changes in the seasonality of individual discharge components, and subsequently, of the overall hydrological regimes, and thus, of the timings of the future water availability are obvious, which will certainly impact the existing water resource management practices in Pakistan.Further, in view of sparse real and sporadic proxy observations and the low fidelity of present-day climate models, it is speculated that the challenges of future water availability assessment, such as, ascertaining prevailing changes from the observations, obtaining future climate projections and quantifying and reducing the uncertainties [92] seem to be far more severe for the UIB.Further, the non-climatic factors, such as, the socioeconomic developments, landuse/landcover changes, burgeoning population and inadequate management of existing water resources in the region pose multiple pressures on the future water resources, and may worsen the severity of climatic impacts on the future water availability.

Figure 3 .
Figure 3. (a) Area elevation curve (AEC), along with the demarcation of considered elevation bands, height of used meteorological stations, band-wise surface area and glacier cover; (b) Spatial description of twelve elevation bands identified based on AEC.Note: The elevations are given in the unit of meters while the areas are given in the unit of 1000 km 2 .

Figure 5 .
Figure 5.Long-term trends and variability of the UIB hydrological regime at the Besham Qila site.The blue bars indicate the times of a year when the UIB hydrological regime features significant negative trends while the red bars indicate the opposite.Inset shows the flow duration curve.

Figure 8 .
Figure 8. Results of near-future hypothetical climate change scenario (2013-2030), suggesting the continuation of prevailing climatic trends to near-future.Inset more clearly shows the snowmelt and rainfall outflow.

Table 1 .
High-altitude meteorological stations and outlet discharge gauge.

Table 2 .
The List of seven CORDEX South Asia experiments for RCP8.5 and their CMIP5 forcing.

Temperature Snowfall Rainfall 1. METEOROLOGICAL DATA INPUT DISTRIBUTION BY ELEVATION ZONES Infiltration Control Snow & Glacier Melt Evaporation Losses Soil Moisture Control Flash Behavior from High Intensity Rainfall Very Slow Runoff Slow Runoff Medium Runoff Fast Runoff Surface Runoff Interflow Upper Zone Groundwater Deep Zone Groundwater Lake or Reservoir Routing Control Generated Stream flow 5. MODIFICATION BY WATERSHED STORAGES 6. EVALUATION WITH RECORDED STREAM FLOW DATA 4. TIME DISTRIBUTION OF RUNOFF 3. RUNOFF COMPONENT ALLOCATION 2. WATERSHED MOISTURE BALANCE COMPUTATION Figure 2. The
schematic diagram of the UBC watershed model.

Table 3 .
Geophysical characteristics of the UIB watershed by elevation zone.Elevation is given in meters while areas are given in km 2 .

Table 4 .
Monthly change factors derived for the period 2013-2030 relative to 1995-2012 for devising a near-future scenario.

Table 5 .
Calibration and validation results of the UBC model applied over the UIB at Besham Qila.