Assessment of Climate Change Impact on Reservoir Inflows Using Multi Climate-models under Rcps—the Case of Mangla Dam in Pakistan

Assessment of climate change on reservoir inflow is important for water and power stressed countries. Projected climate is subject to uncertainties related to climate change scenarios and Global Circulation Models (GCMs). This paper discusses the consequences of climate change on discharge. Historical climatic and gauging data were collected from different stations within a watershed. Bias correction was performed on GCMs temperature and precipitation data. After successful development of the hydrological modeling system (SWAT) for the basin, streamflow was simulated for three future periods (2011–2040, 2041–2070, and 2071–2100) and compared with the baseline data (1981–2010) to explore the changes in different flow indicators such as mean flow, low flow, median flow, high flow, flow duration curves, temporal shift in peaks, and temporal shifts in center-of-volume dates. From the results obtained, an overall increase in mean annual flow was projected in the basin under both RCP 4.5 and RCP 8.5 scenarios. Winter and spring showed a noticeable increase in streamflow, while summer and autumn showed a decrease in streamflow. High flows were predicted to increase, but median flow was projected to decrease in the future under both scenarios. Flow duration curves showed that the probability of occurrence of high flow is likely to be more in the future. It was also noted that peaks were predicted to shift from May to July in the future, and the center-of-volume date of the annual flow may vary from −11 to 23 days in the basin, under both RCP 4.5 and RCP 8.5. As a whole, the Mangla basin will face more floods and less droughts in the future due to the projected increase in high and low flows, decrease in median flows and greater temporal and magnitudinal variations in peak flows. These outcomes suggest that it is important to consider the influence of climate change on water resources to frame appropriate guidelines for planning and management.


Introduction
Global mean temperature has increased by 0.6 • C [1] over the course of the 20th century.Climate models estimate [2][3][4][5][6][7][8] that the global average temperature is likely to increase 4.0 • C by the conclusion of the 21st century [9].Reliable prediction of climate is pre-requisite to comprehend its impacts on hydrology and water resources [10].
Various authors used SRES scenarios for climate change impact studies [11][12][13][14][15], nowadays those scenarios have become outdated.Most of the research to date in the Jhelum and Upper Indus Basin have utilized a few GCMs under SRES scenarios for climate change impact studies [16][17][18][19][20][21][22].The SRES scenarios exaggerate resource accessibility and are unlikely in upcoming production outputs from fossil fuels [23].Representative concentration pathways (RCPs) are new scenarios and these overcome the shortcoming of the SRES scenarios.The RCPs are not linked with exclusive socioeconomic assumptions or emissions scenarios.However, these are based on the groupings of economic, technological, demographic, policy, and future institutional challenges of mitigation and adaptation.Another benefit of RCPs is their better resolution that helps in performing regional and local comparative studies [23].The uncertainties in future climate originate due to internal climate variability, model uncertainty, and scenario uncertainty [24].Majone et al. (2016) reported a rise in the average temperature of the Noce basin, which is located in Italy, from 2 to 4 • C depending on the climate model.The study indicated a rise in an annual average precipitation from 2% to 6% with more changes in winter and autumn.The water yield showed an increase under SRES scenarios [25].
Few studies have been done in Jhelum Basin or its tributaries regarding climate change [18,[26][27][28][29].Only two studies projected future climate using GCMs under SRES scenarios by three or fewer GCMs [27,28].Akhtar et al. (2008) observed trends of rainfall in three decades  for Upper Indus Basin (UIB) near the Jhelum Basin.It was stated that the values of decadal escalation in rainfall at the major station of UIB, i.e., Skardu, Shahpur, and Dir climate stations were 22, 103 and 120 mm, respectively [30].Furthermore, the mean annual temperature and mean annual precipitation will rise in UIB in the 21st century from 0.3 to 4.8 • C, 19% to 113%, respectively [30].In that study, a few GCMs were used under SRES scenarios.A gap in the literature was found, and to fill that gap this study is much needed.
This study aims to answer the following major questions: What will be a probable climate of Mangla watershed?How will inflow into Mangla Reservoir change because of climate change using GCMs under RCPs?The selection of GCMs intensely affects the projected climate [31].In this paper, seven GCMs: BCC-CSM 1.1-m, CCSM4, CSIRO BOM ACCESS1-0, GFDL-CM3, MIROC5, MRI-CGCM3, and UKMO-HadGEM2 under two RCPs 4.5 and 8.5 were selected to simulate future discharge using hydrological model SWAT.The first novelty of this study is the use of different GCMs to cover a wide range of possibilities in future climate.The second novelty of this research is the use of new scenarios (RCPs) to cover the uncertainty regarding emission scenarios.This study will be useful to planners and decision makers when planning and applying suitable water management practices for water resources to adapt to the impacts of climate change.
A principal amount of water enters in the Mangla reservoir from March to August.Generally, in May, the maximum quantity of water enters the reservoir.During October to February very little flow comes into the reservoir, mostly less than 400 m 3 /s.Due to snowmelt, the flow of water starts increasing in March and increases to an absolute peak in the middle of May.More than 75% of the flow comes into the reservoir from March to August.In the remaining six months, less than 25% of the flows reach the reservoir.
The climate of the Jhelum River Basin is principally related to the altitude variations in the basin.The altitude changes from the north to the south, so climate also follows the same trend.Temperature varies from the subtropical in the southern part of the watershed where the elevation is less than Mean monthly temperature and precipitation in Mangla watershed and sub-basins are specified in Figure 2. In general, precipitation distribution in the watershed is bi-modular.First, a larger peak comes in March in the form of snowfall and a second lower peak comes in July in the form of monsoon rainfall.During the observed period , the highest amount of precipitation happens in the northern hilly part of the Kunhar sub-basin.There is a noteworthy spatial variation in precipitation over the Mangla watershed.The average annual precipitation in the northern and the southern parts are 1893 mm, and 846 mm, respectively.About half of the annual precipitation for the northern area occurs from December to March in the form of snow.Besides rain and snowfall, permanent glaciers are the sources of stream flow.The temperature in the watershed varies extremely.On one hand, in the northern part of the basin, the temperature regularly drops below 0 • C from December to March.On the other hand, in the southern part of the basin, temperature can reach 50  The temperature data from observed climatological gauges were used to estimate the lapse rate for the Mangla basin.With the increase in elevation, temperature follows a negative trend.The Lapse rate for the basin was calculated based on the relationship between the average daily temperature and the altitude of the climate stations.Figure 3 displays a very good correlation between temperature and elevation.The Lapse rate for the study area is −6.7 °C/km.From Figure 4 it is clear that temperature in the observed period has increased.The rate of increase in Tmax and Tmin is 0.339 °C and 0.165 °C, respectively, per decade.The rate of increase in Tmax is more than the rate of increase in Tmin.

Meteorological Data
The observed daily meteorological data (1979-2010) were collected from the Pakistan Meteorological Department (PMD), Surface Water Hydrology Project (SWHP) of Water and Power Development Authority (WAPDA) Pakistan, Indian Meteorological Department (IMD), and National Centers for Environmental Prediction's Climate Forecast System Reanalysis (CFSR).The SWHP is mainly for measurement of discharge, but also observes certain climatic variables, namely precipitation and temperature.The data collected includes daily, Tmax and Tmin, precipitation, solar radiation, wind speed and relative humidity.An inventory of meteorological stations is presented in Table 1.The location of the climatological stations is displayed in Figure 1.If only observed weather station data is utilized, then amount precipitation in three sub-basins is less than total flows.Practically, this is impossible.This means that installed weather stations are not sufficient to represent the climate of mountainous watershed.
CFSR data can be used in the data-scarce region [33,34].The observed climatic data were missing in some stations and do not cover the entire basin, that is why CFSR data were used in the Mangla Basin to overcome this limitation.Data of six out of 26 stations used in this study is taken from CFSR database and in the other five stations missing data for the significant period were filled in with CFSR data.Daily precipitation, Tmax, Tmin, solar radiation, wind speed and relative humidity, having a resolution of 0.50 • × 0.50 • is available from 1979 to 2011.
The temperature data from observed climatological gauges were used to estimate the lapse rate for the Mangla basin.With the increase in elevation, temperature follows a negative trend.The Lapse rate for the basin was calculated based on the relationship between the average daily temperature and the altitude of the climate stations.Figure 3 displays a very good correlation between temperature and elevation.The Lapse rate for the study area is −6.7 • C/km.The temperature data from observed climatological gauges were used to estimate the lapse rate for the Mangla basin.With the increase in elevation, temperature follows a negative trend.The Lapse rate for the basin was calculated based on the relationship between the average daily temperature and the altitude of the climate stations.Figure 3 displays a very good correlation between temperature and elevation.The Lapse rate for the study area is −6.7 °C/km.From Figure 4 it is clear that temperature in the observed period has increased.The rate of increase in Tmax and Tmin is 0.339 °C and 0.165 °C, respectively, per decade.The rate of increase in Tmax is more than the rate of increase in Tmin.From Figure 4 it is clear that temperature in the observed period has increased.The rate of increase in Tmax and Tmin is 0.339

Discharge Data
The river network of the Jhelum Basin and its tributaries are shown in Figure 1.Mangla Dam is at the outlet point of the entire basin.Therefore, the watershed area of the Jhelum River Basin, which contributes to the runoff into the Mangla dam was considered for the hydrological analysis.The daily discharge data from eight flow gauging stations was collected from the WAPDA, Pakistan.The data availability period along with other characteristics is presented in Table 2.The average monthly flow (1979-2010) from the different sub-basins and also from the Jhelum River Basin at Mangla is presented in Table 3.The mean monthly runoff of the Jhelum River at Mangla Dam varies between 309 and 1929 m 3 /s.The minimum flow occurs in November, while the maximum flow happens in May.Variations in the temperature and precipitation pattern lead to prominent changes in the stream flows.The mean annual discharge at the Domel gauging station, Muzaffarabad and Azad Patten is 326 m 3 /s, 326 m 3 /s and 833 m 3 /s, respectively (Table 3).The contribution of Kanshi and kahan tributary is negligible.Temperature (°C)

Discharge Data
The river network of the Jhelum Basin and its tributaries are shown in Figure 1.Mangla Dam is at the outlet point of the entire basin.Therefore, the watershed area of the Jhelum River Basin, which contributes to the runoff into the Mangla dam was considered for the hydrological analysis.The daily discharge data from eight flow gauging stations was collected from the WAPDA, Pakistan.The data availability period along with other characteristics is presented in Table 2.The average monthly flow (1979-2010) from the different sub-basins and also from the Jhelum River Basin at Mangla is presented in Table 3.The mean monthly runoff of the Jhelum River at Mangla Dam varies between 309 and 1929 m 3 /s.The minimum flow occurs in November, while the maximum flow happens in May.Variations in the temperature and precipitation pattern lead to prominent changes in the stream flows.The mean annual discharge at the Domel gauging station, Muzaffarabad and Azad Patten is 326 m 3 /s, 326 m 3 /s and 833 m 3 /s, respectively (Table 3).The contribution of Kanshi and kahan tributary is negligible.

DEM
The Digital Elevation Model (DEM) was downloaded from NASA Shuttle Radar Topographic Mission (SRTM).The SRTM DEM has a resolution of 30 × 30 m at the equator and was delivered in mosaicked 1-arc second product pans for easy download and use.The basin and sub-basins defined mechanically from the DEM in the SWAT model (Table S5).Neelam and Kunhar sub-basins are the main source of snow-melt runoff in the Mangla Dam, having more than 70% area above 3000 m MSL.DEM statistics showed that maximum, minimum, mean elevation and the standard deviation is 6276 m, 261 m, 3041 m and 1581 m, respectively.

Soil Data
Soil data was downloaded from the Digital Soil Map of the World (DMSW).DSMW vector data from FAO Soils Portal for South Asian countries was used and projected to WGS-1984 UTM Zone-43N coordinate system.Soil data was classified into 11 classes (Figure 5).S5).Neelam and Kunhar sub-basins are the main source of snow-melt runoff in the Mangla Dam, having more than 70% area above 3000 m MSL.DEM statistics showed that maximum, minimum, mean elevation and the standard deviation is 6276 m, 261 m, 3041 m and 1581 m, respectively.

Soil Data
Soil data was downloaded from the Digital Soil Map of the World (DMSW).DSMW vector data from FAO Soils Portal for South Asian countries was used and projected to WGS-1984 UTM Zone-43N coordinate system.Soil data was classified into 11 classes (Figure 5).The Dominant soil group is Gleyic Solonchak which covers 48.61% of the area of the basin.The remaining sets comprise Calcaric Phaeozems, Mollic Planosols, Haplic Chernozems, Haplic Solonetz, Calcic Chernozems, Gelic Regosols, Gleyic Solonetz, Luvic Chernozems, Lithic Leptosols and Dystric Cambisols, which occupy 22.95%, 20.24%, 1.62%, 1.51%, 1.15%, 1.09%, 1.09%, 0.71%, 0.70% and 0.35% of the area, respectively.The sand, silt, clay and rock amounts of each soil class, as well as their soil parameters, were determined by the Food and Agriculture Organization of the United Nations (FAO) digital soil images of the globe (Table 4).Soil properties such as soil bulk density, texture, soil electric conductivity, soil composition, and soil available water capacity of clay, silt and sand can be obtained from the dataset [32].The soil data set can be found in polygon or in a grid format.The Dominant soil group is Gleyic Solonchak which covers 48.61% of the area of the basin.The remaining sets comprise Calcaric Phaeozems, Mollic Planosols, Haplic Chernozems, Haplic Solonetz, Calcic Chernozems, Gelic Regosols, Gleyic Solonetz, Luvic Chernozems, Lithic Leptosols and Dystric Cambisols, which occupy 22.95%, 20.24%, 1.62%, 1.51%, 1.15%, 1.09%, 1.09%, 0.71%, 0.70% and 0.35% of the area, respectively.The sand, silt, clay and rock amounts of each soil class, as well as their soil parameters, were determined by the Food and Agriculture Organization of the United Nations (FAO) digital soil images of the globe (Table 4).Soil properties such as soil bulk density, texture, soil electric conductivity, soil composition, and soil available water capacity of clay, silt and sand can be obtained from the dataset [32].The soil data set can be found in polygon or in a grid format.

Future Climate Data
The future climate data under RCPs were downloaded.Seven GCMs under two scenarios from CMIP5 were considered for this study.These GCMs cover diverse resolutions, varying from 0.94 • × 1.25 • to 2.8 • × 2.8 • , come from different climate centers all around the world and are updated beyond the year 2000 [35].The data for these GCMs, for selected RCPs, were downloaded for Tmax, Tmin, and precipitation.The Tmax, Tmin and precipitations used are from seven GCMs under two RCPs.The forcing intensities of these two RCPs are 4.5 W/m 2 and 8.5 W/m 2 , respectively, and approximately conform to the medium and high condition.The GCMs used for climate projection in the study area are presented in Table 5.These GCMs cover the period from 1979 to 2100, which is divided into base period (1979-2010) and future three-time horizons (the 2020s: 2011-2040, 2050s: 2041-2070 and 2080s: 2071-2100).S1 and S2).The combination of GCMs can be used even though the selected GCMs may not necessarily be the best models for the area [36].Linear scaling (LS) aims to perfectly match the monthly average of corrected values with observed ones [37][38][39].The monthly corrected values are constructed upon the differences between observed and raw GCMs data.The temperature is typically corrected with an additive and precipitation is typically corrected with a multiplier on a monthly basis.Monthly differences of the climate data, are obtained using observed period (1981-2010) of raw GCMs and observed data.Following Equations ( 1) and ( 2) are applied to correct GCMs future precipitation and temperature data.
T f uture,daily = T GCM,daily, f uture + T (observed,monthly) − T (GCM, observed,monthly) , P f uture,daily = P GCM,daily, f uture × ( P (observed,monthly) Then climate data was analyzed relative to baseline climate on the annual and seasonal basis for future horizons.

The SWAT Model Description
A wide range of different available hydrological models was assessed and explored in order to select the most appropriate model to conduct this research.The selection considers not only the research objectives that determine the complexity and structure of the model, but also the data requirements because a model's development can be hindered by data availability for the model's calibration and validation.For this research, the selection criteria for a hydrological model needed to include the hydrology of the area of the catchment, and the spatial effects of climate change on runoff.SWAT can simulate discharge, sediment yield, water quality and land management practices.For this study, SWAT is used for simulation of discharge only.
SWAT is used to simulate the flow of very small to very large watershed not only in the US but also in the whole world [40][41][42][43].SWAT can simulate the flow process in a broad range of watersheds [42,44,45].The SWAT hydrological model was established and has been used by Hydro-Quebec for 20 years.It is presently used for forecasting of inflows on all ranges of the watershed.The aforementioned studies showed that SWAT can work efficiently for the simulation of hydrological studies.The following conclusions have been drawn, which are also on the basis of previous studies: i The efficiency of the SWAT model is very high for the hydrological studies for the large catchment.ii Satisfactory simulation is obtained for daily, monthly, seasonally and annual runoffs.iii The performance of the snow-melting process of SWAT is satisfactory.iv Projection of streamflows under climate change is possible.v SWAT is in the public domain.
The Hydrologic Response Unit (HRU) is based on unique land use, soil, and the slope is the smallest unit that SWAT generates.Each HRU simulates discharge separately and is then routed to obtain the total discharge from the watershed.The study area was divided into 26 sub-basins and 375 HRUs were created (Table S6; Figure S3).Bearing in mind the variation in elevation in the Mangla watershed, each sub-basin was divided into elevation bands (Table S7).This would imply a better analysis of snowfall and snowmelt in the basin.SWAT allows the splitting of each sub-basin into a maximum of 10 elevation bands.With the elevation bands, precipitation, Tmax and Tmin are calculated separately for each band.For the ET Hargreaves evapotranspiration method is used because future data of solar radiation, wind speed, and relative humidity was not available.
The hydrology constituents of the model described as follows [46]: The hydrologic replication in SWAT is made on the water balance equation:

Model Calibration and Validation
The ArcSWAT for ArcGIS 10.2.1 interface for SWAT 2012 is used in this study to set up the model to simulate the discharge at the outlet of the Jhelum River Basin at Mangla Dam.This approach has been used successfully [47][48][49][50].Calibration of a model is a procedure in which the model parameters are adjusted in such a manner that the simulated flow captures the discrepancies of the observed flow [51].Calibration is achieved by SWAT calibration and uncertainty programs SWAT-CUP using 5000 iterations and manual calibration to achieve better agreement between simulated and observed values.Measured discharge data of the Jhelum River basin are collected from 1981 to 2010 at all hydrological stations, data from 1986 to 1995 is used for calibration and 1996 to 2005 data is used for validation including two years as a warmup period (1979-1980).

Performance Evaluation
Three performance evaluation parameters were used to check the performance of SWAT to project flow, namely: the coefficient of determination (R 2 ), Nash-Sutcliffe efficiency (NSE), and Percent Bias (PBIAS) [52].R 2 describes the degree of collinearity between simulated and measured data, i.e., the proportion of the variance.R 2 ranges from 0 to 1, with higher values indicating less error variance, and, typically, values greater than 0.5 are considered acceptable.NSE displays how fine the observed plot fits the simulated plot.NSE ranges from 0 to 1, with higher values indicating a less error, and, typically, values greater than 0.5 are considered acceptable.PBIAS measures the average tendency of the simulated data to be larger or smaller than their observed counterparts, in other words, it characterizes the percent mean deviation between observed and simulated flows.PBIAS can be positive or negative, positive means underestimation and negative means overestimation, typically, values of −15% < PBIAS < +15% are considered acceptable [53].Moreover, simulated and observed data are compared graphically to discover how fine simulated flow captures the low and high observed flows.

Impact of Climate Change on Discharge
Hydrological simulations were applied for each of the climate sequences [35].The impact of climate change on the average annual and seasonal discharge was analyzed relative to base period flows.Different indicators such as low flow, median flow, mean flow, high flow, flow duration curves, temporal shift in peaks, and temporal shifts in center-of-volume dates were calculated for the three periods and the results were compared to the baseline period's data so as to explore the impact of climate change on the streamflow in the basin.When analyzing streamflow to construct an installation such as a reservoir streamflow occurrence in the future and magnitude of the streamflow is required.Flow duration curves can give streamflow occurrence and magnitude.These curves present the percentage of times that the flow in a stream is likely to exceed or be equal to a specified value of the flow.These curves can be applied in different kinds of studies such as hydropower management, water resource management, and low and high flow studies [54].The following equation is used to construct the flow duration curves: ), P or the probability of flow is equal to or exceeds a specified value (% of time), M is the rank of events, and N is the number of events in a specified period of time.In the present study, the daily time series were used to construct the flow duration curves for the base period (1981-2010) and for the three future periods: the 2020s, 2050s, and 2080s.

Annual
Figure 6 displays the yearly deltas of temperature and precipitation at each progressive time horizon, compared to the observed period.The scatter plot shows the changes in average annual temperature as a function of the changes in average annual precipitation projected by each GCM.Performance of each GCM is discussed in detail in Table S2.
The annual tendency of the GCMs is bi-vocal: five of GCMs (CSIRO BOM ACCESS1-0, GFDL-CM3, MIROC5, MRI-CGCM3, and UKMO-HadGEM2) projected a rise in annual temperature and precipitation while two GCMs (BCC-CSM 1.1-m and CCSM4) projected a rise in temperature but a decline in precipitation.This means that in the future a wide range of uncertainties is possible in precipitation.The magnitude of temperature increase is accentuated with time.While the range of precipitation variations increased over the horizons.
Under RCP 4.5, change in average annual temperature and precipitation is projected in three horizons.The 2020s, 2050s and 2080s may vary from 0.92 to 1.47 • C, −8.65% to 43%, 1.80 to 3.09 • C, −13% to 37% and 2.09 to 4.64 • C, −7.18% to 43%, respectively, by using seven GCMs.Under RCP 8.5, projected change in average annual temperature and precipitation in horizon the 2020s, 2050s and 2080s may vary from 0.84 to 1.52 • C, −12% to 34%, 2.33 to 3.61 • C, −16% to 51% and 4.17 to 7.81 • C, −21% to 51%, respectively, by using seven GCMs (Figure 6).RCP 8.5 covered a wide range of uncertainties while RCP 4.5 covered a short range of possibilities in projected temperature and precipitation.

Annual
Figure 6 displays the yearly deltas of temperature and precipitation at each progressive time horizon, compared to the observed period.The scatter plot shows the changes in average annual temperature as a function of the changes in average annual precipitation projected by each GCM.Performance of each GCM is discussed in detail in Table S2.6 displays the seasonal deltas of temperature and precipitation at each progressive horizon, compared to the observed base period using seven GCMs under two RCPs.The change in average seasonal precipitation is less in winter and spring while change is precipitation is more for summer and autumn using seven GCMs.Monsoon rainfall is projected to be more intense while, in winter less snowfall is expected due to a projected increase in temperature.
There is a continuous rise in projected Tmax and Tmin in all seasons.The large increase in temperature is expected in the winter.For horizon 2080s, the autumn temperatures could increase by 8.9 • C under GFDL RCP 8.5.Average seasonal change in Tmax, Tmin, and precipitation is projected to vary from 0.2 to 5.3 • C, 0.6 to 2.0 • C and −45.1% to 128.4%, respectively, by using seven GCMs under RCP 4.5.While under RCP 8.5 Tmax, Tmin, and precipitation are projected to vary from 0 to 8.9 • C, 0.4 to 8.2 • C and −37.3% to 180.3%, respectively.

Model Calibration and Validation
For sensitivity analysis, 27 parameters were considered out of which eight were found to be relatively sensitive.Past studies were the point of focus for calibration to get sensitive parameters [44,[55][56][57][58] (Table 7).The t-test offers a measure of sensitivity, the largest absolute value represents higher sensitivity and p-value determined the significance of sensitivity.A value close to zero has more significance.Sensitivity rank is shown in Table 7. Rank of parameters increase from top to bottom as p-value is close to 0 and t-value is very high at the top of Table 7. SFTMP is the least important parameter because of the very high p-value and very low t-value.SMTMP is the most important influencing parameter on flow with very low p-value and very high t-value.Among the sensitive parameters, snow parameters (SMFMX, SMFMN, SFTMP, SMTMP, and TIMP) were found to be more sensitive for Neelam and Kunhar sub-basins as these basins are mainly snow-fed, while for Poonch, Kanshi, Lower Jhelum and Kahan sub-basins other parameters (ALPHA_BF, GWQMN and SOL_AWC) were found to be more sensitive.A brief description of each parameter is mentioned in the SWAT user's manual [59].The optimum value of each parameter is given in Table 8.The calibrated model parameters were applied to simulate the discharge.Table 9 shows the model's evaluation parameters (R 2 , NSE, and PBIAS) calculated using the observed and simulated discharge for calibration (1986)(1987)(1988)(1989)(1990)(1991)(1992)(1993)(1994)(1995) and validation (1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005) at different measuring stations.In the case of calibration, the values of R 2 , NSE, and PBIAS ranged from 0.66% to 0.88%, 0.60% to 0.85% and −10.81% to 13.81% and in the case of validation the values of R 2 , NSE, and PBIAS, vary from 0.64% to 0.83%, 0.55% to 0.81% and −14.09% to 14.52%.Simulated and observed discharge at Azad Pattan has R 2 , NSE, and PBIAS 0.88%, 0.85% and −1.70%, respectively (Table 9).Comparison between observed and simulated flows, including precipitation for the calibration and validation are presented in Figure 7.At different measuring stations, shapes of observed flow were well captured by the shapes of the simulated flow (Figures S5 and S6).Nevertheless, some high peaks and low flows were not matched well.At Azad Pattan, a few peak and low flows were underestimated by the model while on other occasions a few were overestimated.This underestimation/overestimation might be due to the scarcity of precipitation gauges available in the basin.

Annual and Seasonal Variations
The impact of climate change on the average annual and seasonal flow is illustrated in Table 10.Flows are expressed as a percentage change relative to baseline period, for each climate projection and horizon.
Simulated flows using the baseline data  show that inflows in winter are 367 m 3 /s, in spring are 1502 m 3 /s, in summer are 1639 m 3 /s, and in autumn are 509 m 3 /s.The changes in average seasonal flows are very high related to annual variations in discharge (Table 10).Variations in winter and spring discharge are mostly positive, even with the decrease in precipitation.However, the changes in flows are mostly negative for summer and autumn, when using the seven GCMs. Te largest increase in inflows is projected for winter.For the 2080s horizon, the winter flows could increase by more than 150% under MRI-CGCM3 RCP 8.5.The changes in average seasonal flow under

Annual and Seasonal Variations
The impact of climate change on the average annual and seasonal flow is illustrated in Table 10.Flows are expressed as a percentage change relative to baseline period, for each climate projection and horizon.Simulated flows using the baseline data  show that inflows in winter are 367 m 3 /s, in spring are 1502 m 3 /s, in summer are 1639 m 3 /s, and in autumn are 509 m 3 /s.The changes in average seasonal flows are very high related to annual variations in discharge (Table 10).Variations in winter and spring discharge are mostly positive, even with the decrease in precipitation.However, the changes in flows are mostly negative for summer and autumn, when using the seven GCMs. Te largest increase in inflows is projected for winter.For the 2080s horizon, the winter flows could increase by more than 150% under MRI-CGCM3 RCP 8.5.The changes in average seasonal flow under RCPs 4.5 and 8.5 are projected to be −29.1% to 130.7% and −49.4% to 171.0%, respectively, when using the seven GCMs (Table 10).Under RCP 4.5, changes in average seasonal flow are projected to vary in winter, spring, summer, and autumn from −17.2% to 107.7%, −22.7% to 38.4%, −29.1% to 97.0%, and −11.8% to 130.7%, respectively, using the seven GCMs (Table 10).Under RCP 8.5, changes in average seasonal flow are projected to vary in winter, spring, summer, and autumn from 5.5% to 171%, −19.6% to 40.3%, −49.4% to 146.4%, and −15.3% to 141.1%, respectively, using the seven GCMs (Table 10).
Changes in average annual flow are projected to vary from −22.9% to 61.1% and −27.2% to 97.7% under RCP 4.5 and RCP 8.5, respectively, using the seven GCMs (Table 10).A more than 20% increase in discharge is expected when using MRI-CGCM3 and GFDL-CM3, while a less than 20% change in discharge is projected when using MIROC5, UKMO-HadGEM, CISRO BOM ACESS1-0, BCC-CSM 1.1-m, and CCSM4 under RCP 4.5.A more than 30% increase in discharge is expected when using MRI-CGCM3, MIROC5, and GFDL-CM3, while a less than 30% change in discharge is projected using UKMO-HadGEM, CISRO BOM ACESS1-0, BCC-CSM 1.1-m, and CCSM4 under RCP 8.5.CCSM4, under both RCPs, showed a decrease in average annual runoff; BCC-CSM 1.1-m, under RCP 4.5, also showed a decrease in mean annual discharge due to a decrease in future precipitation.Under RCP 4.5, changes in average annual flow projected for the horizons of the 2020s, 2050s, and 2080s may vary from −8.7% to 61.1%, −22.9% to 56.3%, and −12.3% to 59.7%, respectively, when using the seven GCMs (Table 10).Under RCP 8.5, changes in average annual flow projected for the horizons of the 2020s, 2050s, and 2080s may vary from −10.6% to 70.4%, −20.4% to 97.7%, and −27.2% to 74.2%, respectively, when using the seven GCMs.RCP 8.5 covered a wide range of uncertainties while RCP 4.5 covered a short range of possibilities in projected annual flows.The annual tendency of the GCMs is bi-vocal: six GCMs (CSIRO BOM ACCESS1-0, GFDL-CM3, MIROC5, MRI-CGCM3, BCC-CSM 1.1-m, and UKMO-HadGEM2) projected a rise in annual flows while one GCM (CCSM4) projected a decrease in projected flows.This means that there is a wide range of possibilities in projected flows when using these GCMs under the two RCPs.
For discharge anomalies, the reference value is the average discharge over the observed period of 1979-2010.A positive anomaly value indicates that the observed discharge is less than the average discharge from 1979 to 2010, while a negative anomaly indicates that the observed discharge is more than the average discharge from 1979 to 2010.Annual discharge anomalies vary from −81% to 320% and −78% to 322%, respectively, when using the seven GCMs under RCPs 4.5 and 8.5 (Figures 8 and 9).Maximum and minimum anomalies of annual discharge were respectively covered by MRI-CGCM 3.1 and CCSM4.Upon using the other five GCMs, annual anomalies under RCPs 4.5 and 8.5 are mostly between −50% and +50%.

Changes in Flow Duration Curve as well as Low, Medium, and High Flows
Figures 10-12 display the flow duration curves in three future periods (2020s, 2050s and 2080s) under RCP 4.5 and RCP 8.5 scenarios.This comparison shows that the probability of occurrence of low and high flow could be higher in the future in the Mangla basin under both scenarios.Tables 11-13 show the projected changes in high flow (Q5), median flow (Q50), and low flow (Q95) in the 2020s, 2050s, and 2080s with respect to the baseline under both scenarios.Under both scenarios, Q5 and Q95 were projected to increase in all three future periods using five or more GCMs in the Mangla basin.However, Q50 was predicted to decrease by more than five GCMs.This decrease in median flow is most likely due to the decrease in median annual precipitation.The change in high flow and flow duration curve show that the frequency of floods and their magnitudes will increase in the future which will create a lot of management problems in the basin.Flooding can not only cause economic   11-13 show the projected changes in high flow (Q5), median flow (Q50), and low flow (Q95) in the 2020s, 2050s, and 2080s with respect to the baseline under both scenarios.Under both scenarios, Q5 and Q95 were projected to increase in all three future periods using five or more GCMs in the Mangla basin.However, Q50 was predicted to decrease by more than five GCMs.This decrease in median flow is most likely due to the decrease in median annual precipitation.The change in high flow and flow duration curve show that the frequency of floods and their magnitudes will increase in the future which will create a lot of management problems in the basin.Flooding can not only cause economic  11-13 show the projected changes in high flow (Q 5 ), median flow (Q 50 ), and low flow (Q 95 ) in the 2020s, 2050s, and 2080s with respect to the baseline under both scenarios.Under both scenarios, Q 5 and Q 95 were projected to increase in all three future periods using five or more GCMs in the Mangla basin.However, Q 50 was predicted to decrease by more than five GCMs.This decrease in median flow is most likely due to the decrease in median annual precipitation.The change in high flow and flow duration curve show that the frequency of floods and their magnitudes will increase in the future which will create a lot of management problems in the basin.Flooding can not only cause economic losses but also loss of life.Nonetheless, with proper utilization and management of the increased flow, Pakistan can actually increase hydropower and food production in the basin.

Temporal Shifts in Peak Flows
In Figure 13, the mean monthly discharge of the baseline period (1981-2010) is plotted against the mean monthly discharge in the future periods (2020s, 2050s, and 2080s) to explore temporal shifts and magnitudes of peak flows at Mangla dam.At this site, a definite advance/delay and increase in peak flows were projected in all three future periods under both scenarios.The peak flows were projected to shift, with increase under both scenarios.This shows that the basin will not only face an increase in frequency and magnitude of floods (mentioned in previous sections) but will also face the shift of these floods from May to July/August.The projected monthly discharge under RCP 4.5 may vary from 200 to 2400 m 3 /s and December and August are respectively the months of minimum and maximum flow using six GCM except MRI-CGCM3.The projected monthly discharge under RCP 8.5 may vary from 243 to 2820 m 3 /s and December and August are the months of minimum and maximum flows using six GCM except MRI-CGCM3.A few GCMs show a shift in discharge, which is due to the increase in temperature and early snowmelt.MRI-CGCM3, MIROC5, GFDL CM3, and UKMO-HadGEM projected increases in flow in summer due to an increase in the intensity of monsoon events.In Figure 13, the mean monthly discharge of the baseline period (1981-2010) is plotted against the mean monthly discharge in the future periods (2020s, 2050s, and 2080s) to explore temporal shifts and magnitudes of peak flows at Mangla dam.At this site, a definite advance/delay and increase in peak flows were projected in all three future periods under both scenarios.The peak flows were projected to shift, with increase under both scenarios.This shows that the basin will not only face an increase in frequency and magnitude of floods (mentioned in previous sections) but will also face the shift of these floods from May to July/August.The projected monthly discharge under RCP 4.5 may vary from 200 to 2400 m 3 /s and December and August are respectively the months of minimum and maximum flow using six GCM except MRI-CGCM3.The projected monthly discharge under RCP 8.5 may vary from 243 to 2820 m 3 /s and December and August are the months of minimum and maximum flows using six GCM except MRI-CGCM3.A few GCMs show a shift in discharge, which is due to the increase in temperature and early snowmelt.MRI-CGCM3, MIROC5, GFDL CM3, and UKMO-HadGEM projected increases in flow in summer due to an increase in the intensity of monsoon events.To determine the impact of climate change on the timing of streamflows, an indicator such as center-of-volume date (CVD)-a date at which half of the total volume of streamflow passes through at a gauging station for a specific time period-was used in the present study.CVD was calculated according to the equation described in [44].Table 14 shows the changes in CVD under RCP 4.5 and RCP 8.5, with respect to the baseline period, in the three future periods at Mangla dam.Table 14  To determine the impact of climate change on the timing of streamflows, an indicator such as center-of-volume date (CVD)-a date at which half of the total volume of streamflow passes through at a gauging station for a specific time period-was used in the present study.CVD was calculated according to the equation described in [44].Table 14 shows the changes in CVD under RCP 4.5 and RCP 8.5, with respect to the baseline period, in the three future periods at Mangla dam.Table 14 shows that about half of the flow, on an average, of each year, in the Mangla basin, passed through by 13 June for the baseline period .However, this is predicted to shift in the future.The positive values show delay in flow and negative values show advance in flows.The change in CVD were projected under RCP 4.5 and RCP 8.5 in all three future periods, with −11-23 days.On one hand, by using the best selected GCM (Table S2) CVD is expected to advance the total time by 4-9 days, on the other hand using the least appropriate GCM CVD is expected to cause a delay of 7-20 days.

Conclusions
Pakistan is one of the most water-stressed countries in the world and its water resources are greatly vulnerable to changing climate.In the present study, the possible impacts of climate change on the water resources of the Mangla basin were assessed under RCP 4.5 and RCP 8.5 using seven GCMs.The Jhelum River originates from the northwestern side of Pirpanjal and is one of the main tributaries of the Indus River.The SWAT hydrological model was used to simulate streamflow in the basin for the future.The model was calibrated and validated for the periods of 1985-1995 and 1996-2005 respectively, at different gauging stations.Three indicators (Coefficient of determination, Nash Efficiency, and percentage deviation), and graphical representations of differences between observed and simulated data were used to check the performance of the model.Bias-corrected temperature and precipitation data for the period of 1979-2100 under RCP 4.5 and RCP 8.5 scenarios of seven GCMs were fed into SWAT model to simulate the streamflow for the future.In this study, the simulated streamflow data was divided into three future periods (2020s, 2050s, and 2080s) and was compared with the baseline period .Different indicators like changes in mean flow, low flow, median flow, high flow, flow duration curves, temporal shift in peaks, and temporal shifts in center-of-volume dates were used to investigate the changes in streamflow under RCP 4.5 and RCP 8.5.The main conclusions of the study are the following:

Figure 1 .Figure 1 .
Figure 1.The Mangla River basin showing climatic stations in red dots and flow gauging stations in green dots and sub-basins.

Figure 2 .
Figure 2. Mean monthly temperature and precipitation in Mangla watershed and sub-basins.

Figure 3 .
Figure 3. Relationship between temperature and elevations of the Jhelum Basin.

Figure 2 .
Figure 2. Mean monthly temperature and precipitation in Mangla watershed and sub-basins.

Figure 2 .
Figure 2. Mean monthly temperature and precipitation in Mangla watershed and sub-basins.

Figure 3 .
Figure 3. Relationship between temperature and elevations of the Jhelum Basin.

Figure 3 .
Figure 3. Relationship between temperature and elevations of the Jhelum Basin.
Model (DEM) was downloaded from NASA Shuttle Radar Topographic Mission (SRTM).The SRTM DEM has a resolution of 30 × 30 m at the equator and was delivered in mosaicked 1-arc second product pans for easy download and use.The basin and sub-basins defined mechanically from the DEM in the SWAT model (Table

Figure 5 .
Figure 5. Soil map of the study area.

Figure 5 .
Figure 5. Soil map of the study area.

Figure 6 . 5 Figure 6 .
Figure 6.Annual deltas of climate change projections.The annual tendency of the GCMs is bi-vocal: five of GCMs (CSIRO BOM ACCESS1-0, GFDL-CM3, MIROC5, MRI-CGCM3, and UKMO-HadGEM2) projected a rise in annual temperature and precipitation while two GCMs (BCC-CSM 1.1-m and CCSM4) projected a rise in temperature but a

Figure 8 .
Figure 8. Annual discharge anomalies under RCP 4.5 using the seven GCMs.

Figure 9 .
Figure 9. Annual discharge anomalies under RCP 8.5 using the seven GCMs.

Figure 9 .
Figure 9. Annual discharge anomalies under RCP 8.5 using the seven GCMs. 5.3.2.Changes in Flow Duration Curve as well as Low, Medium, and High Flows Figures 10-12 display the flow duration curves in three future periods (2020s, 2050s and 2080s) under RCP 4.5 and RCP 8.5 scenarios.This comparison shows that the probability of occurrence of low and high flow could be higher in the future in the Mangla basin under both scenarios.Tables11-13show the projected changes in high flow (Q5), median flow (Q50), and low flow (Q95) in the 2020s, 2050s, and 2080s with respect to the baseline under both scenarios.Under both scenarios, Q5 and Q95 were projected to increase in all three future periods using five or more GCMs in the Mangla basin.However, Q50 was predicted to decrease by more than five GCMs.This decrease in median flow is most likely due to the decrease in median annual precipitation.The change in high flow and flow duration curve show that the frequency of floods and their magnitudes will increase in the future which will create a lot of management problems in the basin.Flooding can not only cause economic

Figure 9 .
Figure 9. Annual discharge anomalies under RCP 8.5 using the seven GCMs.
loss of life.Nonetheless, with proper utilization and management of the increased flow, Pakistan can actually increase hydropower and food production in the basin.

Table 1 .
Inventory of climate stations.
• C and 0.165 • C, respectively, per decade.The rate of increase in Tmax is more than the rate of increase in Tmin.

Table 2 .
Hydrological stations in the Jhelum Basin.

Table 2 .
Hydrological stations in the Jhelum Basin.

Table 4 .
Soil constituents and parameters in Mangla River Basin.

Table 5 .
Global Circulation Models (GCMs) used for climate projection in the study area under representative concentration pathway (RCP) 4.5 and 8.5.

Table 6 .
Projected change in mean annual and seasonal Tmax, Tmin, and PPT in three-time slices.

Table 7 .
Sensitivity analysis of parameters using SWAT-CUP.

Table 8 .
Model parameters used to calibrate discharge.

Table 9 .
Statistical analysis for calibration and validation.

Table 9 .
Statistical analysis for calibration and validation.

Table 10 .
Change in discharge (%) using the seven GCMs under the two RCPs.

Table 11 .
Percent future changes in high flows with respect to the baseline flow (1981-2010) under RCP 4.5 and RCP 8.5 scenarios in the Mangla River Basin.

Table 12 .
Percent future changes in median flows with respect to the baseline flow (1981-2010) under RCP 4.5 and RCP 8.5 scenarios in the Mangla River Basin.

Table 13 .
Percent future changes in low flows with respect to the baseline flow (1981-2010) under RCP 4.5 and RCP 8.5 scenarios in the Mangla River Basin.

Table 14 .
Future changes in center-of-volume dates (CVD) with respect to the baseline period (1981-2010) at under both scenarios, RCP 4.5 and RCP 8.5, in the Mangla basin.