A Multi-GCM Assessment of the Climate Change Impact on the Hydrology and Hydropower Potential of a Semi-Arid Basin ( A Case Study of the Dez Dam Basin , Iran )

In this paper, the impact of climate change on the climate and discharge of the Dez Dam Basin and the hydropower potential of two hydropower plants (Bakhtiari and Dez) is investigated based on the downscaled outputs of six GCMs (General Circulation Models) and three SRES (Special Report on Emission Scenarios) scenarios for the early, mid and late 21st century. Projections of all the scenarios and GCMs revealed a significant rise in temperature (up to 4.9 ◦C) and slight to moderate variation in precipitation (up to 18%). Outputs of the HBV hydrologic model, enforced by projected datasets, show a reduction of the annual flow by 33% under the climate change condition. Further, analyzing the induced changes in the inflow and hydropower generation potential of the Bakhtiari and Dez dams showed that both inflow and hydropower generation is significantly affected by climate change. For the Bakhtiari dam, this indicates a consistent reduction of inflow (up to 27%) and electricity generation (up to 32%). While, in the Dez dam case, the inflow is projected to decrease (up to 22%) and the corresponding hydropower is expected to slightly increase (up to 3%). This contrasting result for the Dez dam is assessed based on its reservoir and hydropower plant capacity, as well as other factors such as the timely releases to meet different demands and flow regime changes under climate change. The results show that the Bakhtiari reservoir and power plant will not meet the design-capacity outputs under the climate change condition as its large capacity cannot be fully utilized; while there is room for the further development of the Dez power plant. Comparing the results of the applied GCMs showed high discrepancies among the outputs of different models.


Introduction
Anthropogenic global warming and its consequences, especially in the arid and semi-arid regions, received particular attention in recent years as many scholars documented the occurrence and dominance of droughts, a rise in the temperature, and an increase in the atmospheric water demand, accompanied by a reduction in the precipitation and runoff [1][2][3][4].Additionally, in Iran, many studies support the fact that during recent decades the climate has experienced variations, mostly toward hot and dry conditions [5][6][7][8][9][10][11][12][13][14].
Anticipated climatic changes can alter the hydrological regimes, such as the amount of discharge or the timing of the surface flow on both the regional and local (catchment) scales [15][16][17], with socio-economic and environmental consequences.Projections of the impact of climate change on the hydrological conditions on a global scale show that for low and mid-latitude regions, a reduction of freshwater could exacerbate their water-management problems; while the higher latitudes will experience higher amounts of surface flow [18,19].
Several studies worldwide predicted that the combination of the temperature increase and precipitation variation during the present century will result in a significant change of the runoff.According to the study by Tong et al. on the variability of the discharge, under the B1 scenario and the urbanization in the Las Vegas Wash watershed in the USA, climate change is the main deriving factor of future changes in the watershed, causing a wintertime discharge decrease and a summertime discharge increase [20].The study on the future water availability in Bangladesh revealed that climate change has a significant impact on runoff and evapotranspiration because the region will face a higher irrigation demand, a decline of groundwater, and a variability of rainfall and runoff (both increasing and decreasing) [21].
The projection of future climatic conditions in Senegal by Tall et al. showed an increase in the temperature, evaporation, and precipitation by the mid-21st century and a decrease of these parameters by the late-21st century.They also reported that dependent on the applied scenarios, the runoff will change (both increasingly and decreasingly).This could lead to an arid climate dominance in the region [22].He et al. assessed the hydrologic sensitivity to climate change in the upper San Joaquin River basin in California by employing projected temperature and precipitation datasets, the reported temperature rise (between 1.5 and 4.5 • C) and precipitation variations (between 80 and 120%) [23].They showed that climate change can lead to annual streamflow variations between −41 and 16%.They also detected an earlier shift of most of the streamflow by 15 to 46 days as a result of the temperature increase, which is the cause of the higher seasonal variability of streamflow.Modeling the future climatic and hydrologic response to climate change in Spain revealed a 1.5-3.3• C temperature increase, a 6-32% precipitation decrease, and a 2-54% runoff decrease [24].Xu and Luo, by employing seven GCMs (General Circulation Models) and the A1B scenario in two semi-arid and humid regions in China, reported dramatic changes in the temperature (up to +8.6 • C), precipitation (up to 139%), and seasonal discharge (up to 304%), as well as an increase in the extreme flows and seasonal shifts of discharge [25].Investigating the scenarios of climate change impact on the river flow in Western Kenya using different GCMs and SRES (Special Report on Emission Scenarios) showed that climate change has the potential to significantly alter the river flow [26].
Additionally, assessments of climate change impact on the river flow in Iran show considerable variations of surface water resources across different parts of the country.Simulation of streamflow in the North of Iran, through forcing the hydrologic model with climatic projections based on different SRES scenarios, revealed increases and decreases of the discharge for the wet and dry seasons, respectively, with an overall increase in the annual discharge [27].According to another study on the future changes of the climatic condition and runoff across Iran, the future temperature and precipitation are projected to vary by ±6 • C and ±60%, respectively [28].Additionally, these changes will result in higher rates of annual evaporation and a runoff reduction, as well as seasonal variations of runoff (an increase in winter and a decrease in spring).Rafiei Emam et al. showed that the hydrologic response of the Raza-Ghahavand region (a semi-arid region in Iran) to climate change is a decrease in precipitation and an increase in temperature, which leads to less groundwater recharge and a lower soil water content [29].A study on the impact of climate change on the surface water supply in the Zayandeh-Rud River Basin in Central Iran showed an increase of 0.4-0.76• C in annual temperature, a decrease of 14-38% in the precipitation and a decrease of 8-43% in the runoff [30].Modeling the future climate and water resources under the A1B, A2, and B2 scenarios in the Karkheh River Basin (located in the West and Southwest of Iran) revealed a temperature increase and reduction of the water yield in the basin.This reduction is considerable from April to September as a result of the temperature increase and precipitation decrease [31].
It is widely documented that based on the future simulations, the modeled wintertime discharge increase is usually followed by a springtime discharge decrease, which is due to the temperature increase which causes the snow to melt sooner and for there to be less snow [15][16][17]32,33].
Additionally, a global scale study on the river flow variations caused by the A2 and B1 emission scenarios suggested an increase in the seasonality of the river flow, i.e., an increase in the high flows and a decrease in the low flows, for about one-third of the Earth [34].Elsner et al., using the A1B and B1 scenarios, showed that climate change causes significant shifts and changes in the amount and timing of the runoff in the Pacific Northwest [35].Pervez and Henebry reported that a combination of the land use and climatic change under different SRES scenarios can cause some changes and shifts in the amount of discharge and flow timing in the Brahmaputra River Basin [36].Boyer et al. studied the projected variability of the wintertime and springtime discharges at a regional scale and found that under climate change, the wintertime discharge increases, while the springtime discharge decreases [16].
One of the vulnerable industries to climate change is hydropower generation, which completely relies on the amount of precipitation, snow cover, snowmelt, streamflow, and the timing of the flow, all of which show high inter-annual variability [37][38][39].Based on the IPCC report, as a result of climate change, the hydropower generation potential is projected to drop by up to 6% [40].Therefore, it is essential to adapt the water resources management to the future climatic changes [41].In Iran, climatic changes and their repercussion on hydropower generation have been studied by Jahandideh et al. [42] and Jamali et al. [43], focusing on the Karun and Karkheh river basins, respectively.They reported a reduction of the hydropower generation potential in both of the studied basins.

Problem Statement and Objectives of the Study
The Dez Dam Basin is located in Iran's Southwest and the discharge from this basin is planned to provide water for different sectors, such as the agriculture, industry, fishery and hydropower generation sectors through existing and currently-in-constructed dams.In the planning and management of water resources, the base period recorded dataset is normally considered.However, the records from recent years and future projections of climate models suggest that global warming and climate change could alter the climate indicators and hydrologic conditions.Therefore, for basins like the Dez Dam Basin, which is not completely developed yet, it is crucial to assess the variations in the amount of discharge on a basin scale with respect to the climate change scenarios; this is essential for climate change adaptation.In this study, we investigated the variations of hydroclimatic conditions induced by climate change in the Dez Dam Basin and its consequences with regards to the hydropower generation potential through two large dams in the basin.For this aim, first, the precipitation and temperature values were projected based on the three SRES scenarios.Then, using the HBV hydrologic model, the discharge was simulated under climate change conditions.Next, the hydropower generation potential was calculated for the two hydropower plants of the Bakhtiari and Dez dams.Since many scholars identified GCMs as one of the significant sources of uncertainty in hydro-climatic studies [25,33,[44][45][46][47][48][49], in this study, an ensemble of six GCMs and three emission scenarios were employed to consider the different climatic conditions within a range, offered by the outputs of different GCMs and scenarios.This assessment provides a useful means to the modify water resources management strategies, considering the repercussions of climate change on the surface water resources and the hydropower plants of the Dez Dam Basin.

Study Area and Data
The Dez Dam Basin (hereafter referred to as DDB) is located in the Southwest of Iran between 31 • 35 51"-34 • 7 46" N and 48 • 9 15"-50 • 18 37" E and is the upstream tributary of the great Karun catchment.Figure 1 shows the location of the Dez Dam Basin in Iran, as well as the streams and hydrometric stations in this basin, based on the data layers acquired from the Iran Water Resources Management Company.Figure 1 also illustrates the four sub-basins of the DDB (delineated by the ArcSWAT tool in ArcGIS), namely, Tireh, Marbereh, Sazar, and Bakhtiari.The two main rivers of Sazar and Bakhtiari drain the basin and join at the point known as Tange Panj to form the Dez River.
Figure 1 shows the location of the Dez dam, which is currently operational, and the Bakhtiari dam, which is under construction.As depicted in Figure 1, the Bakhtiari dam is located at the end point of the Bakhtiari sub-basin.Additionally, the inflow to the multi-purpose Dez dam comes from the Sazar and Bakhtiari rivers.The Bakhtiari dam is mainly designated for hydropower generation.While the storage of the Dez reservoir provides water for the domestic, agricultural, and industrial sectors, it is also used to generate electricity.Table 1 presents a description of the studied sub-basins and the hydro-meteorological indicators of the study area.The presented coordinates and elevations in Table 1 belong to the centroid of each sub-basin.
Water 2018, 10, x FOR PEER REVIEW 4 of 22 storage of the Dez reservoir provides water for the domestic, agricultural, and industrial sectors, it is also used to generate electricity.Table 1 presents a description of the studied sub-basins and the hydro-meteorological indicators of the study area.The presented coordinates and elevations in Table 1 belong to the centroid of each sub-basin.The catchment is located within a mountainous area known as the Zagros Mountains with limited shallow aquifers in some parts of the sub-basins 1 to 3. The basin's climate is characterized as semi-arid to Mediterranean with warm summers and cold winter and less than 800 mm of precipitation per year.

GCM-Scenario Ensemble
To project the future T and P, six different GCMs were applied, including CCSM3, ECHAM5-OM, GFDL-CM2.1,HadCM3, INM-CM3.0, and IPSL-CM4.Each of the applied GCMs couple different components of the Earth system on different grid resolutions.Applying a multi-model ensemble of GCMs enables us to consider a wide range of predictions.Table 2 provides a description of the applied GCMs.The three SRES emission scenarios of A1B, A2, and B1 were applied to project the future T and P for the 2011-2030 (the 2020s), 2046-2065 (the 2050s), and 2081-2100 (the 2080s) time horizons.Therefore, in this study, an ensemble of 54 combinations of GCMs, scenarios, and time horizons (hereafter referred to as GSTs) were used to assess the climatic change and its impact on the discharge and hydropower potential of the DDB.
The related assumptions regarding each SRES scenario and the corresponding CO2  The catchment is located within a mountainous area known as the Zagros Mountains with limited shallow aquifers in some parts of the sub-basins 1 to 3. The basin's climate is characterized as semi-arid to Mediterranean with warm summers and cold winter and less than 800 mm of precipitation per year.

GCM-Scenario Ensemble
To project the future T and P, six different GCMs were applied, including CCSM3, ECHAM5-OM, GFDL-CM2.1,HadCM3, INM-CM3.0, and IPSL-CM4.Each of the applied GCMs couple different components of the Earth system on different grid resolutions.Applying a multi-model ensemble of GCMs enables us to consider a wide range of predictions.Table 2 provides a description of the applied GCMs.The three SRES emission scenarios of A1B, A2, and B1 were applied to project the future T and P for the 2011-2030 (the 2020s), 2046-2065 (the 2050s), and 2081-2100 (the 2080s) time horizons.Therefore, in this study, an ensemble of 54 combinations of GCMs, scenarios, and time horizons (hereafter referred to as GSTs) were used to assess the climatic change and its impact on the discharge and hydropower potential of the DDB.The related assumptions regarding each SRES scenario and the corresponding CO 2 concentrations can be found in the IPCC's fourth assessment technical report (AR4) [50].
To downscale the future T and P in the study area, under three SERS scenarios (A1B, B1, and A2) during the 2020s, the 2050s, and the 2080s, the stochastic weather generator of LARS-WG was applied.This generated the future time series based on the probability distribution of the base period data and the correlations between the observations.The detailed description of LARS-WG is provided by Semenov [51] and Semenov and Stratonovitch [52].

Hydrological Modeling
The HBV-light (HBV-light-GUI, V. 4.0.0.6) semi-distributed hydrological model is used to simulate the streamflow in the DDB.The model considers snow routine, soil moisture routine, response function, and the flood routing of the basin.This model simulates Q with T, P, and E as the input data.The model provides the option to link sub-basins, therefore, for SUB-1, SUB-2, and SUB-3, the semi-distributed mode was applied.A detailed description of the HBV model structure and routines are provided by Seibert and Vis [53].Prior to the application of the model, it is necessary to calibrate its parameters by the trial and error procedure, as recommended by Bergström [54].After model calibration, validation, and after ensuring its efficiency, it was run with a projected climate series to simulate the future runoff in the basin.

Model Calibration and Validation
There are several statistics to identify the efficiency of a model.In this research, we examined the goodness of fit with different criteria (Equations ( 1)-( 3)) to diagnose the efficiency of the model in the simulation of the streamflow in the calibration and validation periods.In details, the Nash-Sutcliff measure (R eff ), the coefficient of determination (R 2 ), and the mean annual difference (M diff ) are utilized to evaluate the model performance regarding efficiency, the timing of the flow, and the average error, respectively.
where, Q sim and Q obs are the simulated and observed discharge data and n is the number of data.

Modeling the Two Reservoirs System
The two reservoirs systems of the Bakhtiari and Dez dams are modeled based on the water balance equation (Equation (4)) by considering all the restrictions on the operation of the reservoirs and the power plant, for power generation.
where, S t+1 is the reservoir volume at the end of the t period (beginning of the t + 1 period); S t is the reservoir volume at the beginning of the t period; Q t is the inflow to the reservoir during the t period; RE t is the outflow from the reservoir to generate energy during the t period; RD t is the outflow from the reservoir to meet the downstream demand during the t period; SPL t is the spill from the reservoir during the t period; EVAP t is the evaporation from the reservoir during the t period; and ADD t is the volume of the added flow from the upstream reservoir during the t period.
Additionally, based on the minimum operational reservoir volume (S min ) and the maximum reservoir volume (S max ), the following conditions are considered in the model:

•
If S min < S t < S max , the outflow from the reservoir and the hydropower generation is equal to the water and energy demands of that specific month and neither deficit nor spill will occur.

•
If S max < S t , considering the upper limit for the reservoir volume (S t ≤ S max ), the reservoir volume (S t ) is equal to S max and the excessive amount of water (S max − S t ) will spill.In this condition, there is no water or energy deficit and the secondary energy could be produced.

•
If S t < S min , considering the lower limit for the reservoir operation (S min ≤ S t ), the reservoir volume will be substituted with the minimum operational reservoir storage, i.e., S min .Therefore, in that month, the deficit is equal to (S t − S min ).In this condition, some or all the demands may not be met.If S t − S min ≥ 0, the water is released based on the priorities to meet high prioritized demands.Additionally, the energy generation will be affected in accordance with the reduction of the amount of water flowing through a turbine.
The power generation is associated with the installed capacity, efficiency, and plant factor as well as to the inflow and hydraulic head.Equation (5) shows the relationship that is used to calculate the energy generation.
where, P t is the power generated during the t period (W); γ is water specific weight (N/m 3 ); Q t is the inflow to the turbine during the t period (m 3 /s); H t is the net hydraulic head on the turbine during the t period (m); and e t is the power plant efficiency during the t period.
The following constraints are applied for the calculation of the hydropower generation based on the limits of the hydropower plant (Equations ( 6)-( 8)).
where, PCC is the installed capacity (MW); Q min is the minimum turbine inflow (m 3 /s); Q max is the maximum turbine inflow (m 3 /s); H min is the minimum required head to operate the power plant (m); H max is the maximum head to operate the power plant (m).

Projected Impact of Climate Change on Temperature and Precipitation Rates
Figure 2 illustrates the average of the projected temperatures, obtained by different GSTs and the average temperature of the base period (1986-2010).As depicted in Figure 2, all the GCMs and scenarios suggest a small temperature change during the 2020s.However, for SUB-4, a temperature rate increase up to 1.2 • C is projected for the first time horizon.While, based on all the scenarios, a significant rise in the temperature is projected for the time horizons of the 2050s and the 2080s.Additionally, the average of the projected temperatures in the 2080s is higher than in the 2050s.Generally, the B1 scenario suggests a smaller temperature increase during the 2050s and the 2080s; while, scenarios A2 and A1B revealed a higher temperature during these time horizons.Additionally, regarding the differences between the GCMs, it seems that INCM3 mostly shows the lowest rates of the temperature increase, especially during the 2050s and 2080s; while GFDL, in most cases, revealed the highest temperature increase during the 2050s, based on the scenarios A1B, A2, and B1, respectively.Additionally, during the 2080s, the highest temperatures based on all the scenarios are suggested by ECHAM5.As it is shown in Figure 2, the amplitude of the multi-GCM projections of temperature suggests greater uncertainties in the 2080s, compared to the 2050s and the 2020s.Conversely, the values of the projected temperatures by different GCMs are approximately similar for the first time horizon (the 2020s).
Water 2018, 10, x FOR PEER REVIEW 7 of 22 Additionally, the average of the projected temperatures in the 2080s is higher than in the 2050s.Generally, the B1 scenario suggests a smaller temperature increase during the 2050s and the 2080s; while, scenarios A2 and A1B revealed a higher temperature during these time horizons.Additionally, regarding the differences between the GCMs, it seems that INCM3 mostly shows the lowest rates of the temperature increase, especially during the 2050s and 2080s; while GFDL, in most cases, revealed the highest temperature increase during the 2050s, based on the scenarios A1B, A2, and B1, respectively.Additionally, during the 2080s, the highest temperatures based on all the scenarios are suggested by ECHAM5.As it is shown in Figure 2, the amplitude of the multi-GCM projections of temperature suggests greater uncertainties in the 2080s, compared to the 2050s and the 2020s.
Conversely, the values of the projected temperatures by different GCMs are approximately similar for the first time horizon (the 2020s).The maximum temperature increase is consistently identified by all models for the 2080s time horizon under the A2 scenario, except for SUB-1.In SUB-1, the greatest rise in temperature is revealed by projections of ECHAM5 under the A1B scenario.The greatest temperature increase in the 2080s time horizon compared to the base period ranges from 4 °C (SUB-1) to 4.9 °C (SUB-4).
In a similar way, Figure 3 illustrates the average of the projected annual precipitation for different GSTs, as well as the average of the precipitation for the base period (1986-2010).As it is obvious, most of the GSTs suggest increasing the precipitation amounts compared to the base period.However, the projections obtained by employing IPSL mostly suggest a reduction of the amount of precipitation in the future time horizons, except for the outputs of the IPSL-A1B-2050s for SUB-1 and The maximum temperature increase is consistently identified by all models for the 2080s time horizon under the A2 scenario, except for SUB-1.In SUB-1, the greatest rise in temperature is revealed by projections of ECHAM5 under the A1B scenario.The greatest temperature increase in the 2080s time horizon compared to the base period ranges from 4 • C (SUB-1) to 4.9 • C (SUB-4).
In a similar way, Figure 3 illustrates the average of the projected annual precipitation for different GSTs, as well as the average of the precipitation for the base period (1986-2010).As it is obvious, most of the GSTs suggest increasing the precipitation amounts compared to the base period.However, the projections obtained by employing IPSL mostly suggest a reduction of the amount of precipitation in the future time horizons, except for the outputs of the IPSL-A1B-2050s for SUB-1 and SUB-3; these show a slight increase in the amount of precipitation.Additionally, projections for SUB-2 show both a rise and decline of the precipitation, compared to the base period, which is mostly inclining towards a decrease of the precipitation rate in the 2020s.Based on the results, the greatest precipitation changes mostly occur in the 2080s time horizon.Regarding the differences between GCMs, it is clear that the IPSL performs differently.Additionally, GFDL shows greater deviations in the 2050s from the base period.In addition, some inconsistent results were found among the sub-basins.For example, the projections mostly suggest a precipitation decrease in the first time horizon (the 2020s) in SUB-2, while for other sub-basins the precipitation is mostly projected to be increased.
In total, the projections of most of the GSTs indicate a warmer future and a moderate increase in the precipitation amount.However, the results of the precipitation are more anomalous, showing both an increase and decrease of the amount of precipitation.Additionally, for the late 21st century, the GCM-scenario ensemble suggests a wider range of projections for the temperature and precipitation, implying that the uncertainties of the projections increase in a wider time span.The ensemble revealed similar results for the 2020s.

Hydrological Modeling of the Dez Dam Basin
To simulate the streamflow in the DDB, the semi-distributed HBV model is applied.First, the model calibration was completed using at least 12 years of observed discharge data.After the calibration, the model performance was evaluated using at least 5 years of observed discharge data, independent of the calibration period.It is taken into consideration that the calibration and validation periods cover both dry and wet conditions to ensure that the model is capable of working with different conditions.Regarding the differences between GCMs, it is clear that the IPSL performs differently.Additionally, GFDL shows greater deviations in the 2050s from the base period.In addition, some inconsistent results were found among the sub-basins.For example, the projections mostly suggest a precipitation decrease in the first time horizon (the 2020s) in SUB-2, while for other sub-basins the precipitation is mostly projected to be increased.
In total, the projections of most of the GSTs indicate a warmer future and a moderate increase in the precipitation amount.However, the results of the precipitation are more anomalous, showing both an increase and decrease of the amount of precipitation.Additionally, for the late 21st century, the GCM-scenario ensemble suggests a wider range of projections for the temperature and precipitation, implying that the uncertainties of the projections increase in a wider time span.The ensemble revealed similar results for the 2020s.

Hydrological Modeling of the Dez Dam Basin
To simulate the streamflow in the DDB, the semi-distributed HBV model is applied.First, the model calibration was completed using at least 12 years of observed discharge data.After the calibration, the model performance was evaluated using at least 5 years of observed discharge data, independent of the calibration period.It is taken into consideration that the calibration and validation periods cover both dry and wet conditions to ensure that the model is capable of working with different conditions.
Table 3 represents the information regarding the length of the calibration and validation periods, as well as the model performance at each sub-basin.Overall, the statistics show an acceptable performance of the calibrated model.As it is obvious for the validation period, negative values of Mdiff (%) are obtained for SUB-1, 2, and 4, which are mainly due to the reduction of the annual flow during the validation period which coincides with the final years of the observation period.Figure 4a-d shows the observed and simulated streamflow at the outlet of the four sub-basins during the calibration and validation periods.

Hydrological Simulation under Climate Change Scenarios
After ensuring the efficiency of the calibrated model, the streamflows were simulated under different climate change scenarios with the projected T and P time series as inputs.Figures 5-8 illustrate the annual pattern of the simulated streamflows under the considered GSTs, as well as the observed discharge during the base period (1989-2009), for SUB-1 to SUB-4, respectively.As it is shown, the amplitude of the future simulations shows a significant variability implying

Hydrological Simulation under Climate Change Scenarios
After ensuring the efficiency of the calibrated model, the streamflows were simulated under different climate change scenarios with the projected T and P time series as inputs.Figures 5-8 illustrate the annual pattern of the simulated streamflows under the considered GSTs, as well as the observed discharge during the base period (1989-2009), for SUB-1 to SUB-4, respectively.As it is shown, the amplitude of the future simulations shows a significant variability implying inconsistencies between the simulated streamflows by different GCMs, mainly for the 2050s and 2080s.Obviously, the differences between the simulations are smaller during the 2020s, while the highest variations are associated with the late 21st century.These findings are in agreement with Vidal and Wade [55] who found a highly increasing variance in the late 21st century due to the spread of the outputs of different GCMs.In addition, the greatest inconsistencies between the simulated streamflows correspond to the higher amounts of discharge, which mainly occur during April and March.Based on the simulations presented in Figures 5-8, not only the amount of discharge, but also the flow regime experiences some changes under the studied climate change conditions.As for SUB-1 (Figure 5), the advancing shift of the peak flows compared to the base period is obvious for all time horizons, while for SUB-3 (Figure 7), the advancing shift occurs during the 2080s, and for SUB-4, (Figure 8) it was found that the peak of the annual hydrograph occurs earlier during both the 2050s and 2080s time horizons.Similar results indicating the advancing shifts of the peak flows under climate change conditions are widely documented by different researchers.Regarding the timing of the peak flow, Nohara et al. reported that under climate change conditions, the peak flows shift earlier [18].This is because of the temperature increase which causes the snowmelt to start sooner.Similarly, the hydrologic regime change and advancing shifts of the peak flows is reported by Gan et al. in the Naryn river basin (Central Asia) [17].Based on the simulations presented in Figures 5-8, not only the amount of discharge, but also the flow regime experiences some changes under the studied climate change conditions.As for SUB-1 (Figure 5), the advancing shift of the peak flows compared to the base period is obvious for all time horizons, while for SUB-3 (Figure 7), the advancing shift occurs during the 2080s, and for SUB-4, (Figure 8) it was found that the peak of the annual hydrograph occurs earlier during both the 2050s and 2080s time horizons.Similar results indicating the advancing shifts of the peak flows under climate change conditions are widely documented by different researchers.Regarding the timing of the peak flow, Nohara et al. reported that under climate change conditions, the peak flows shift earlier [18].This is because of the temperature increase which causes the snowmelt to start sooner.Similarly, the hydrologic regime change and advancing shifts of the peak flows is reported by Gan et al. in the Naryn river basin (Central Asia) [17].Table 4 summarizes the percentage of the variations of discharge for each sub-basin based on the applied GSTs.
As it is obvious in Table 4, the sub-basins of the DDB respond differently to climate change.For example, the surface flow of the SUB-3 (in most cases) seems to be less affected by climate change compared to other sub-basins.Table 4 summarizes the percentage of the variations of discharge for each sub-basin based on the applied GSTs.
As it is obvious in Table 4, the sub-basins of the DDB respond differently to climate change.For example, the surface flow of the SUB-3 (in most cases) seems to be less affected by climate change compared to other sub-basins.As shown in Figures 5-8 and Table 4, the results vary based on the GCMs used.For instance, the amount of annual flow of SUB-3 for the 2080s is projected to increase based on INCM3 (A2 and B1).However, IPSL suggests a reduction of flow in SUB-3 for the same time horizon and scenarios.Generally, among the six applied GCMs, IPSL revealed higher reductions in the flow.Therefore, not only were the results for the studied sub-basins different, but the projected response of a certain subbasin by different GCMs could also be dissimilar.This result is in good agreement with other studies as many researchers have so far identified the GCMs as the main source of uncertainty and reported that the GCM selection can cause significant deviations in the results of the climate change impact assessments [25,44,46].In this regard, Graham et al., by examining different RCMs, GCMs and hydrological models, concluded that the choice of GCMs is more determinative than other factors [45].Likewise, Habets et al. applied several climate models, hydrological models, downscaling methods, and emission scenarios to simulate the future water resources and reported that the uncertainties caused by the climate models are 3-4 times greater than other factors [47].As mentioned by Turco et al., significant differences and deviations between the results of different GCMs/RCMs imply large uncertainties regarding the use of a certain combination of GCMs and RCMs [48].Additionally, Fang et al. documented that most of the streamflow uncertainties correspond to the structural uncertainty of GCMs [33].Vidal and Wade mentioned that the highly increasing variance of the late 21st century is mostly due to the spread of the results obtained from different GCMs [55].As shown in Figures 5-8 and Table 4, the results vary based on the GCMs used.For instance, the amount of annual flow of SUB-3 for the 2080s is projected to increase based on INCM3 (A2 and B1).However, IPSL suggests a reduction of flow in SUB-3 for the same time horizon and scenarios.Generally, among the six applied GCMs, IPSL revealed higher reductions in the flow.Therefore, not only were the results for the studied sub-basins different, but the projected response of a certain sub-basin by different GCMs could also be dissimilar.This result is in good agreement with other studies as many researchers have so far identified the GCMs as the main source of uncertainty and reported that the GCM selection can cause significant deviations in the results of the climate change impact assessments [25,44,46].In this regard, Graham et al., by examining different RCMs, GCMs and hydrological models, concluded that the choice of GCMs is more determinative than other factors [45].Likewise, Habets et al. applied several climate models, hydrological models, downscaling methods, and emission scenarios to simulate the future water resources and reported that the uncertainties caused by the climate models are 3-4 times greater than other factors [47].As mentioned by Turco et al., significant differences and deviations between the results of different GCMs/RCMs imply large uncertainties regarding the use of a certain combination of GCMs and RCMs [48].Additionally, Fang et al. documented that most of the streamflow uncertainties correspond to the structural uncertainty of GCMs [33].Vidal and Wade mentioned that the highly increasing variance of the late 21st century is mostly due to the spread of the results obtained from different GCMs [55].

Variation of the Inflow to Dez and Bakhtiari Reservoirs and Hydropower Generation under Climate Change
The number of annual inflows to the Bakhtiari and Dez reservoirs under climate change conditions and the percentage of changes of the simulated inflows compared to the base period are obtained based on the two reservoir system of the Bakhtiari and Dez dams.The percentage of changes in the inflow to each of these reservoirs compared to the base period are shown for each GST in Table 5.The long-term averages of the discharges at the outlet of Bakhtiari and Sazar are 5106.88and 2938.21Mm 3 , respectively.Considering that the Bakhtiari Dam is located at the end-point of the Bakhtiari basin and that the Dez dam is located downstream of Tange Panj, the inflow to the Bakhtiari reservoir equals the Bakhtiari basin's discharge and the inflow to the Dez reservoir is the total discharge from the upstream basins, which equals 8045.1 Mm 3 .Based on the results, the percentages of changes of the simulated inflow to the Bakhtiari reservoir under climate change varies between −1.7% (HadCM3-B1-2050) and −26.8% (IPSL-A1B-2080).Additionally, the inflow to the Dez reservoir for future time horizons is simulated to decrease by up to 21.8% (IPSL-A1B-2080) and increase by up to 3.6% (HadCM3-B1-2080) compared to the base period.The hydropower generation potential for the base period and for the future time horizons were also calculated using the inflows obtained from the previous step.Studying the monthly average of hydropower generation for the Bakhtiari and Dez reservoirs under climate change conditions versus the base period (Figures 9 and 10) shows that the ranges of simulations by the use of different GCMs and scenarios are greater in the late 21st century, while the amount of produced energy for different GCMs in the 2020s time horizon indicates a higher agreement of the results.and scenarios are greater in the late 21st century, while the amount of produced energy for different GCMs in the 2020s time horizon indicates a higher agreement of the results.The percentages of the deviations of the projected hydropower generation potential from the base period for both reservoirs are presented in Table 5 for all the GSTs.Based on the results of the different GSTs, the potential of the hydropower generation of the Bakhtiari power plant was simulated to decrease between 0.9% (HadCM3-B1-2050) and 31.8%(IPSL-A1B-2080) compared to the base period.This is in agreement with the reduction of discharges to the Bakhtiari reservoir.On the contrary, the hydropower generation potential of the Dez power plant is simulated to be slightly higher than the base period, for all the GSTs because the highest increase reaches up to 2.8% (CCSM-A2-2080).Meanwhile, its inflow is mostly projected to decrease for future time horizons.The percentages of the deviations of the projected hydropower generation potential from the base period for both reservoirs are presented in Table 5 for all the GSTs.Based on the results of the different GSTs, the potential of the hydropower generation of the Bakhtiari power plant was simulated to decrease between 0.9% (HadCM3-B1-2050) and 31.8%(IPSL-A1B-2080) compared to the base period.This is in agreement with the reduction of discharges to the Bakhtiari reservoir.On the contrary, the hydropower generation potential of the Dez power plant is simulated to be slightly higher than the base period, for all the GSTs because the highest increase reaches up to 2.8% (CCSM-A2-2080).Meanwhile, its inflow is mostly projected to decrease for future time horizons.As mentioned above, in the Dez reservoir case, there are some inconsistencies between the changes in its inflow and hydropower generation caused by climate change.While the results for the Bakhtiari reservoir shows more consistency.The difference between the responses of the two reservoirs in terms of their hydropower generation potential could be attributed to the different operational rules due to different purposes, as well as due to other factors such as the size of the reservoirs and the installed capacity of their hydropower plants.To study the possible reasons, the time series of inflow, outflow, spill, and the water level in the reservoirs for the base period and for the future time horizons were assessed.The time series pertaining to the base period and the future projections under climate change (for one of the GSTs as an example) are presented in Figures 11 and  12.As mentioned above, in the Dez reservoir case, there are some inconsistencies between the changes in its inflow and hydropower generation caused by climate change.While the results for the Bakhtiari reservoir shows more consistency.The difference between the responses of the two reservoirs in terms of their hydropower generation potential could be attributed to the different operational rules due to different purposes, as well as due to other factors such as the size of the reservoirs and the installed capacity of their hydropower plants.To study the possible reasons, the time series of inflow, outflow, spill, and the water level in the reservoirs for the base period and for the future time horizons were assessed.The time series pertaining to the base period and the future projections under climate change (for one of the GSTs as an example) are presented in Figures 11 and 12.
Regarding the different responses of the two studied reservoirs, the following facts should be taken into consideration:

•
It should be noted that the Dez dam is a multi-purpose dam which provides water for different purposes during specific times to meet specified demands.Therefore, the releases from its reservoir are planned and only the part of the release or spill, which is not greater than the penstock or turbine capacity, contributes to the power generation.However, the release from Bakhtiari is only for hydropower generation purposes.

•
Having considered that the whole capacity of the Dez hydropower plant is relatively small, compared to its inflow and releases, a significant proportion of the releases (or spills) does not contribute to power generation.Meanwhile, the Bakhtiari reservoir, with a capacity of 5.16 Bm 3 and an average inflow of 5.11 Bm 3 , can save most of the inflows with negligible spills, both in the base period and in the future time horizons (Figures 11a and 12a).

•
Additionally, the large capacity of the hydropower plant does not pose any limitation on the energy production.Therefore, in the case of Bakhtiari, there is a direct relationship between the changes in the rates of inflow and energy generation.

•
A comparison between the simulated inflow of the Dez reservoir during the future time horizons and the base period suggests that under the climate change conditions, a fewer number of floods and fewer inflows and peak flows would lead to fewer losses through spill (Figures 11b and 12b), which means that more water could be saved in the reservoir to be used to generate electricity.
Considering the reasons explained above, the future changes induced by climate change in the hydropower generation potential of the Dez power plant are not consistent with the changes in the inflow of the Dez reservoir.Additionally, the slight increase in the hydropower generation potential could be attributed to the changing of the regime of discharges with smaller peaks, leading to fewer spills and more water remaining to produce electricity.Regarding the different responses of the two studied reservoirs, the following facts should be taken into consideration:

•
It should be noted that the Dez dam is a multi-purpose dam which provides water for different purposes during specific times to meet specified demands.Therefore, the releases from its reservoir are planned and only the part of the release or spill, which is not greater than the penstock or turbine capacity, contributes to the power generation.However, the release from Bakhtiari is only for hydropower generation purposes.

•
Additionally, the small capacity of Dez reservoir (2.7 Bm 3 ), compared to the discharge from its draining catchment (8 Bm 3 ), causes considerable spills (Figures 11b and 12b).Considering the reasons explained above, the future changes induced by climate change in the hydropower generation potential of the Dez power plant are not consistent with the changes in the inflow of the Dez reservoir.Additionally, the slight increase in the hydropower generation potential could be attributed to the changing of the regime of discharges with smaller peaks, leading to fewer spills and more water remaining to produce electricity.

Conclusions
The impact of climate change on the climate, discharge, and the hydropower generation potential in the Dez Dam Basin was studied based on the downscaled outputs from six GCMs and three SRES scenarios for the three time horizons.The study revealed that the basin experiences a significant temperature rise in the mid and late 21st century (up to 4 °C).These changes are accompanied by variations in precipitation, which mostly lean towards a slight increase in the amount of precipitation.In total, the projections of all the scenarios and GCMs indicate a warmer future and small-to-moderate increase in the amount of precipitation.However, the obtained results for precipitation are more anomalous, showing both an increase and a decrease in the amount of precipitation.To simulate the future discharge at the outlet of each sub-basin under climate change conditions, the calibrated HBV hydrologic model was enforced with the projected temperature and precipitation time series.The results mostly suggest a reduction of the annual discharge in the study area.The most significant reduction of the annual flow, compared to the base period, reaches up to 33% in the two sub-basins located upstream of the basin, namely Tireh and Marbereh.Meanwhile, the greatest reductions of the annual flow of the Sazar and Bakhtiari sub-basins reach up to 17 and 27%, respectively.Considering the climate change impacts on the temperature, precipitation, and discharge in different sub-basins of the DDB, it was found that the responses of the four sub-basins

Conclusions
The impact of climate change on the climate, discharge, and the hydropower generation potential in the Dez Dam Basin was studied based on the downscaled outputs from six GCMs and three SRES scenarios for the three time horizons.The study revealed that the basin experiences a significant temperature rise in the mid and late 21st century (up to 4 • C).These changes are accompanied by variations in precipitation, which mostly lean towards a slight increase in the amount of precipitation.In total, the projections of all the scenarios and GCMs indicate a warmer future and small-to-moderate increase in the amount of precipitation.However, the obtained results for precipitation are more anomalous, showing both an increase and a decrease in the amount of precipitation.To simulate the future discharge at the outlet of each sub-basin under climate change conditions, the calibrated HBV hydrologic model was enforced with the projected temperature and precipitation time series.The results mostly suggest a reduction of the annual discharge in the study area.The most significant reduction of the annual flow, compared to the base period, reaches up to 33% in the two sub-basins located upstream of the basin, namely Tireh and Marbereh.Meanwhile, the greatest reductions of the annual flow of the Sazar and Bakhtiari sub-basins reach up to 17 and 27%, respectively.Considering the climate change impacts on the temperature, precipitation, and discharge in different sub-basins of the DDB, it was found that the responses of the four sub-basins are different in many cases, highlighting the noteworthiness of analyzing the impacts of climate change on local scales.
Moreover, in this study, the impact of climate change on the hydropower generation potential of the two hydropower plants in the DDB was investigated.Based on the results, climate change has the potential to significantly alter the hydropower generation potential in this basin.The results showed a reduction in the inflow and electricity generation for the Bakhtiari reservoir.Meanwhile, for the Dez reservoir, the reduction in the inflow was accompanied by a slight increase in the generated electricity.This contrasting result obtained for the Dez reservoir was assessed and has been attributed to the small size of the reservoir (2.7 Bm 3 ) compared to the basin's discharge (8 Bm 3 ), the low capacity of hydropower plant, the different purposes and timely releases from its reservoir, as well as the flow regime changes in the future which cause less spills due to the lower peaks of floods.Overall, the results showed a reduction of the electricity generation at the Bakhtiari power plant.Therefore, its capacity seems to be high considering the future climatic conditions, while, based on the findings, there is room for the further development of the Dez power plant in order to increase its capacity for the production of more electricity.
This study used multiple GCMs for the projection of temperature, precipitation, discharge, and hydropower generation based on the three future time horizons scenarios.The results showed considerable discrepancies in the projected variables when obtained from different GCMs, which indicate the important role of the GCMs in future climatic impact assessments.Therefore, it is highly recommended that for future projections, different GCMs be employed to cover a range of likely projections.

Figure 1 .
Figure 1.The map of the Dez River Basin; (a) the sub-basins, main streams, and location of the hydrometric stations and reservoirs (b) its location in Iran.

Figure 1 .
Figure 1.The map of the Dez River Basin; (a) the sub-basins, main streams, and location of the hydrometric stations and reservoirs (b) its location in Iran.

Figure 4 .
Figure 4.The observed and simulated streamflow during the calibration and validation periods in the sub-basins: SUB-1 (a); SUB-2 (b); SUB-3 (c); and SUB-4 (d).The striped line shown on the graphs separates the calibration and validation periods.

Figure 4 .
Figure 4.The observed and simulated streamflow during the calibration and validation periods in the sub-basins: SUB-1 (a); SUB-2 (b); SUB-3 (c); and SUB-4 (d).The striped line shown on the graphs separates the calibration and validation periods.

Water 2018 ,
10, x FOR PEER REVIEW 10 of 22 streamflows correspond to the higher amounts of discharge, which mainly occur during April and March.
There are similar reports by other researchers about the different responses to climate change of nearby sub-basins within a region.Nazif and Karamouz, by studying the variability of streamflows in Central Iran under climate change, showed that the streamflows are significantly altered [56].However, they noticed that the responses of three adjacent basins to climate change were different.Additionally, Ashraf Vaghefi et al. employed the CGCM model, forced by three SRES scenarios, to project the future water resource availability in the Karkheh River Basin, located in West and Southwest of Iran and found that the freshwater availability increases in the northern parts of the basin but it decreases in the southern regions [57].Additionally, Musau et al.reported that the sensitivity of the four adjacent watersheds, considering their hydrologic response to climate change, were different[26].This highly spatial variability of the responses of the adjacent sub-basins underlines the importance of studying the impact of climate change on the hydrological conditions in local scales in order to get a better perspective of the behavior of each basin, rather than a holistic view.
There are similar reports by other researchers about the different responses to climate change of nearby sub-basins within a region.Nazif and Karamouz, by studying the variability of streamflows in Central Iran under climate change, showed that the streamflows are significantly altered [56].However, they noticed that the responses of three adjacent basins to climate change were different.Additionally, Ashraf Vaghefi et al. employed the CGCM model, forced by three SRES scenarios, to project the future water resource availability in the Karkheh River Basin, located in West and Southwest of Iran and found that the freshwater availability increases in the northern parts of the basin but it decreases in the southern regions[57].Additionally, Musau et al. reported that the sensitivity of the four adjacent watersheds, considering their hydrologic response to climate change, were different[26].This highly spatial variability of the responses of the adjacent sub-basins underlines the importance of studying the impact of climate change on the hydrological conditions in local scales in order to get a better perspective of the behavior of each basin, rather than a holistic view.

22 Figure 11 .
Figure 11.The monthly inflows, outflows, spills, and reservoir storages for the base period; (a): inflow, outflow, and spill of the Bakhtiari reservoir; (b): inflow, outflow, and spill of the Dez reservoir.

Figure 11 .
Figure 11.The monthly inflows, outflows, spills, and reservoir storages for the base period; (a): inflow, outflow, and spill of the Bakhtiari reservoir; (b): inflow, outflow, and spill of the Dez reservoir.

Figure 12 .
Figure 12.The monthly inflows, outflows, spills, and reservoir storages for the 2080s future time horizon based on the projections of the IPSL under the B1 scenario; (a): inflow, outflow, and spill of the Bakhtiari reservoir; (b): inflow, outflow, and spill of the Dez reservoir.

Figure
Figure The monthly inflows, outflows, spills, and reservoir storages for the 2080s future time horizon based on the projections of the IPSL under the B1 scenario; (a): inflow, outflow, and spill of the Bakhtiari reservoir; (b): inflow, outflow, and spill of the Dez reservoir.

Table 1 .
The description of the sub-basins of the Dez River basin.

Table 1 .
The description of the sub-basins of the Dez River basin.

Table 2 .
The description of the GCMs (General Ciculation Models) and the IPCC-AR4 SRES (Special Report on Emission Scenarios) emission scenarios.

Table 3 .
The model performance statistics for the calibration and validation periods.

Table 3 .
The model performance statistics for the calibration and validation periods.

Table 5 .
The percentage of variation of the inflows and the hydropower generation of the Bakhtiari and Dez reservoirs under climate change conditions projected by the six GCMs and three SRES scenarios for the 2020s, 2050s, and 2080s time horizons compared to the base period.