Simulation of the Potential Impacts of Projected Climate Change on Streamﬂow in the Vakhsh River Basin in Central Asia under CMIP5 RCP Scenarios

: Millions of people in Uzbekistan, Turkmenistan, Tajikistan, and Kyrgyzstan are dependent on the freshwater supply of the Vakhsh River system. Sustainable management of the water resources of the Vakhsh River Basin (VRB) requires comprehensive assessment regarding future climate change and its implications for streamﬂow. In this study, we assessed the potential impacts of projected climate change scenarios on the streamﬂow in the VRB for two future periods (2022–2060 and 2061–2099). The probable changes in the regional climate system were assessed using the outputs of ﬁve global climate models (GCMs) under two representative concentration pathways (RCPs), RCP4.5 and RCP8.5. The probable streamﬂow was simulated using a semi-distributed hydrological model, namely the Soil and Water Assessment Tool (SWAT). Evidence of a signiﬁcant increase in the annual average temperature by the end of the 21st century was found, ranging from 2.25 to 4.40 ◦ C under RCP4.5 and from 4.40 to 6.60 ◦ C under RCP8.5. The results of three GCMs indicated a decreasing tendency of annual average precipitation (from − 1.7% to − 16.0% under RCP4.5 and from − 3.4% to − 29.8% under RCP8.5). Under RCP8.5, two GCMs indicated an increase (from 2.3% to 5.3%) in the average annual precipitation by the end of 2099. The simulated results of the hydrological model reported an increasing tendency of average annual streamﬂow, from 17.5% to 52.3% under both RCPs, by the end of 2099. A shift in the peak ﬂow month was also found, i.e., from July to June, under both RCPs. It is expected that in the future, median and high ﬂows might increase, whereas low ﬂow might decrease by the end of 2099. It is concluded that the future seasonal streamﬂow in the VRB are highly uncertain due to the probable alterations in temperature and precipitation. The ﬁndings of the present study could be useful for understanding the future hydrological behavior of the Vakhsh River, for the planning of sustainable regional irrigation systems in the downstream countries, i.e., Uzbekistan and Turkmenistan, and for the construction of hydropower plants in the upstream countries.


Introduction
Water resources are essential for the existence of life on this planet. Sustainable access to these resources is necessary for our livelihood, agriculture, hydropower generation, economic growth, and ecosystems [1]. It is quite well documented that climate change is altering the distribution of freshwater resources, and is also influencing the eco-hydrological processes in different river systems of the world [2][3][4].
In Central Asia, desiccation of the Aral Sea from the 1960s to the late 1980s and shrinking of glaciers in the mountains of Pamir-Alay and Tian Shan are clear indications of climate change in this region. The Amu Darya River drains into the Aral Sea [5]; the Vakhsh River, which was the focus of this study, is the second largest headwater tributary of the Amu Darya River, originating in the high mountains of Pamir-Alay, and is among the most meltwater-dependent river systems in Central Asia [6]. The total annual flow of the Vakhsh River is 20.22 km 3 /year, and the area of irrigated land in this river system is about 172.2 thousand hectares [7]. The Vakhsh River is also a major source for the generation of green energy in Central Asia, with the Nurek Reservoir being the largest one in this region (water storage volume of 10.5 km 3 ).
The most threatening consequence of climate variation in mountainous regions is the rapid melting of snow/glaciers. Global warming is causing the premature melting of snow/glaciers in such regions and is also altering the patterns of precipitation, which will eventually influence the seasonality of the river flows [8][9][10][11]. Deng et al. [12] reported that the temperature lapse rate (in the spring and fall) and the precipitation (in the summer) are major factors influencing the discharge in the Kaidu River Basin in the mountainous region of Central Asia. Several previous hydroclimatic investigations in the South and Central Asian countries have observed an increasing tendency in river flows. The increasing tendencies of flows were mostly associated with the accelerated melting of glaciers [13]. For instance, Mahmood et al. [14] predicted a 15% increase in the future average annual flow of the Jhelum River, which is located in the western Himalayan mountain system. Garee et al. [15] projected that the future average annual flow of the Hunza River, in the Karakoram range, might increase up to 10%. Conversely, the annual flows of Brahmaputra River, one of the major rivers in Asia, are expected to decrease in the future [16]. White et al. [17] projected that the Amu Darya River flow, in Central Asia, is expected to decrease by between 10% and 20%. Kayumov et al. [18] studied the impacts of climate alteration on the glacier covers in Tajikistan, which is situated in Central Asia, and reported a decrease of the glaciated area of the country from 6.0% in the middle of the past century to about 4.8% in 2014. Although past trends of the streamflow of the largest river of Tajikistan (the Vakhsh River) under the context of regional climate change were documented by Kayumov et al. [18], clear manifestations of future climate change and its subsequent impacts on the streamflow in the basin remain unknown.
In the Tian Shan, Pamir, North Karakoram, Himalayan, and Tibetan Plateau mountainous system, in the middle and at the end of the 21st century, continuous warming was projected on annual and seasonal scales [19][20][21][22]. In Central Asia, warming (up to 7.0 • C according to Mannig et al. [23] and up to 6.5 • C according to Reyer et al. [24]) is projected until the end of this century. Bollasina et al. [25] reported an increasing trend in the concentration of aerosols in the atmosphere of Asia. In addition, Xin et al. [26] reported that the increasing concentration of atmospheric aerosols could cause a warm atmosphere in South and Central Asia.
The trends in precipitation in Asia are not consistent. In Central Asia, a persistent decrease in the annual precipitation was found at the end of the 20th century [27]. Unger-Shayesteh et al. [28] also reported a decreasing trend of future precipitation in eastern Pamir. Conversely, Meng et al. [29] indicated a rising trend of precipitation in some regions of Central Asia for the 21st century, for instance in the Tian Shan mountains. Malsy et al. [30] pointed out that over the Central Asian region, mean and maximum annual precipitation, based on the different datasets, will vary by 13% and up to 42%, respectively. Similarly, the study of future precipitation characteristics in river basins (Yellow and Xin Rivers) in China [31,32], in the Middle East [33], and in the westerly-dominated region of South Asia indicated a decreasing trend in summer precipitation [34]. However, an increasing trend of winter precipitation was reported by Luo et al. [35] for the Heihe River Basin. Babur et al. [36] projected possible annual variations in precipitation in the western part of the Himalayas and found a significant increase in precipitation (43% under representative concentration pathway (RCP)4.5 and 51% under RCP8.5). Omani et al. [13] projected changes in winter precipitation in the Pamir-Alay mountain range in Central Asia for the 2070-2099 period relative to the baseline  period for the 14 GCMs under RCP4.5 (from −4.48% to 35.91%) and RCP8.5 (from −8.62% to 61.29%). These contradictions in previous findings suggest that the ability to predict the increasing or decreasing tendency of winter precipitation may vary from model to model.
Hydrological models help to understand and investigate the dynamic behavior of river flows. Globally, many different hydrological models in different regions and at different scales, for instance, lumped, semi-distributed, and distributed, have been used to study the impacts of climate change on hydrological systems [37][38][39][40][41][42][43][44][45][46][47][48][49][50]. These hydrological models can predict hydrological variables and components of interest, particularly when the overall hydrological response of the watershed is of interest in the future. The Soil and Water Assessment Tool (SWAT) has been extensively used in different regions to analyze many hydrological problems, including the assessment of potential water flow variations in future climate scenarios [51][52][53][54].
The main scientific objectives of this study were as follows: (1) To simulate the streamflow of a mountainous catchment with minimum parameter uncertainty and high confidence in prediction using a semi-distributed hydrological model (SWAT); (2) to evaluate the changes in median extreme (low and high) flows, as well as the temporal shifts in peak flows; and (3) to assess the likely impact of climate variation on the monthly, seasonal, and interannual flows and snowmelt of the Vakhsh River Basin (VRB) in a multimodal experiment using the results of five global climate models (GCMs). Section 2 presents the materials and datasets used in this investigation. Section 3 describes the hydrological model, its calibration and validation, the uncertainty of the parameters, and the future scenarios. Section 4 presents the results and discussion. Finally, the conclusions are summarized in Section 5. The results of this study could be helpful for researchers and policy-makers alike.

Study Area
The Vakhsh River Basin, the northeast tributary of the Amu Darya River, is located in Central Asia and flows from north-central Kyrgyzstan to the southwest of Tajikistan. The basin is situated between 37.10 • and 39.74 • N and 68.31 • and 73.70 • E, and has a total length of 524 km with a drainage area of about 39,008 km 2 . The main tributaries of this river are Muksu and Obihingou. Elevation in the basin ranges from 302 to 7050 m above sea level ( Figure 1). The climate of the basin is continental under the domination of the westerlies, with very high local contrasts because of the orographic conditions. The temperature decreases with the increase in elevation, while precipitation has different patterns at different altitudes and aspects. The respective values of the precipitation and average maximum/minimum temperature over the Vakhsh River Basin from 2000 to 2013 were 120 mm, 3.9 • C, 1.5 km. The Vakhsh River is considered one of the most important rivers in terms of power production potential in Tajikistan. Several large-and medium-sized hydropower plants (HPPs) have been installed on this river.
Water 2020, 12, x FOR PEER REVIEW 4 of 34 rivers in terms of power production potential in Tajikistan. Several large-and medium-sized hydropower plants (HPPs) have been installed on this river.

Measured and GCM-Based Climate Datasets
The in situ gauge-based daily observations of the climatic parameters used in this study were obtained from the Agency of Hydrometeorology Committee on Environmental Protection under the Government of the Republic of Tajikistan. The daily observations of the river flow at the Darband hydrological station over the 2000-2013 period were obtained from the Ministry of Energy and Water Resources of the Republic of Tajikistan. The quality-ensured daily time series of maximum temperature (Tmax), minimum temperature (Tmin), average temperature (Tas), and precipitation (Pr) for ten meteorological stations (Khatlon, Yovon, Ganjina, Faizobod, Rasht, Chormagzak, Lakhsh, Bustonobod, Sanglokh, and Dehavz) were acquired for the 1966-2013 period. The gauge-based datasets were utilized for the calibration and validation of the SWAT model over the 2000-2013 period. The overall agreements between the gauge-based and GCM-based daily time series of the climatic parameters were assessed for the 1966-2004 period.
The simulated maximum/minimum temperature and precipitation data of the five GCMs under two RCPs (4.5 and 8.5) were downloaded from the www.2w2e.com website and downscaled using the Climate Change Toolkit [55] for the study area. Table 1 describes the outputs of the GCMs, which were based on the ISI-MIP (Inter-Sectoral Impact Model Intercomparison Project) [56], used in this study. The simulated results of Pr, Tmax, and Tmin of the five GCMs were extracted for the control period (1966-2004) and for two future time periods (2022-2060 and 2061-2099). GCM results were also obtained for medium (4.5 W/m 2 ) and high (8.5 W/m 2 ) concentration scenarios.

Measured and GCM-Based Climate Datasets
The in situ gauge-based daily observations of the climatic parameters used in this study were obtained from the Agency of Hydrometeorology Committee on Environmental Protection under the Government of the Republic of Tajikistan. The daily observations of the river flow at the Darband hydrological station over the 2000-2013 period were obtained from the Ministry of Energy and Water Resources of the Republic of Tajikistan. The quality-ensured daily time series of maximum temperature (Tmax), minimum temperature (Tmin), average temperature (Tas), and precipitation (Pr) for ten meteorological stations (Khatlon, Yovon, Ganjina, Faizobod, Rasht, Chormagzak, Lakhsh, Bustonobod, Sanglokh, and Dehavz) were acquired for the 1966-2013 period. The gauge-based datasets were utilized for the calibration and validation of the SWAT model over the 2000-2013 period. The overall agreements between the gauge-based and GCM-based daily time series of the climatic parameters were assessed for the 1966-2004 period.
The simulated maximum/minimum temperature and precipitation data of the five GCMs under two RCPs (4.5 and 8.5) were downloaded from the www.2w2e.com website and downscaled using the Climate Change Toolkit [55] for the study area. Table 1 describes the outputs of the GCMs, which were based on the ISI-MIP (Inter-Sectoral Impact Model Intercomparison Project) [56], used in this study. The simulated results of Pr, Tmax, and Tmin of the five GCMs were extracted for the control period  and for two future time periods (2022-2060 and 2061-2099). GCM results were also obtained for medium (4.5 W/m 2 ) and high (8.5 W/m 2 ) concentration scenarios. Table 1. Details of the general circulation models (GCMs) used for the projection of the precipitation and maximum/minimum temperature in the Vakhsh River Basin. The emission scenarios (RCPs) and the time periods considered in the calculation of the change factors and future climate change impacts on streamflow are also provided.

The Digital Elevation Model (DEM), Land Cover, and Soil Data
In this study, a digital elevation model (DEM) with a 90 m spatial resolution was obtained from the website of the Shuttle Radar Topographic Mission (SRTM) [57], which was utilized for the demarcation of the catchment boundary and the river network. The land use map of the study area at the 300 × 300 m grid scale was acquired from the Envisat Medium-Resolution Imaging Spectrometer (MERIS) [58]. Different land cover types such as wetlands (0.27%), water bodies (0.92%), urban areas (2.49%), forests (3.35%), agricultural areas (6.65%), bare land (9.80%), snow and ice (10.67%), bare rock gravel (24.30%), and grasslands (41.54%) were identified in the basin.
The soil map for the VRB was derived from version 1.2 of the Harmonized World Soil Database (HWSD), with a scale of 1:5,000,000. The HWSD is freely available on the website of the Food and Agriculture Organization (FAO) of the United Nations [59]. The soil layer obtained from the HWSD contained information on the chemical and physical properties of the soil, e.g., depth from the surface to the bottom of the soil, moist bulk density, clay content, silt content, and sand content. Therefore, a soil map was used directly in the hydrological model. Using the HWSD-based map, the silt, sand, and clay properties were computed on a percentage basis, which satisfied the model input data requirements, such as the saturated hydraulic conductivity and available water capacity. In this study, the prominent soil types were chromic cambisols (40.01%), ferralic arenosols (23.75%), calcic luvisols (22.68%), chromic luvisols (8.57%), haplic phaeozems (2.73%), dystric regosols (1.59%), and lithic leptosols (0.67%).

Soil and Water Assessment Tool (SWAT)
The Soil and Water Assessment Tool, developed by the Agricultural Research Service of the United States Department of Agriculture, is a physically based semi-distributed hydrology model. The simulation in this model can be performed in daily and subdaily time steps [60,61]. A schematic illustration of the SWAT model is presented in Figure 3. It is capable of simulating the impacts of climate change and land management measures on water availability, as well as the sediment transport in small-to-large ungauged river basins on the basis of soil types, land uses, and management circumstances [62]. Initially, this hydrological model was developed to quantify the impact of land-use Water 2020, 12, 1426 7 of 34 practices in large compound catchments. It is a flexible tool, and the model has been frequently applied in mountainous river basins [63][64][65]. It has been widely used for the assessment and adaptation of schemes to climate change, flood control, and hydropower management practices [66,67]. In the SWAT model, a catchment can be classified into different hydrological response units (HRUs) by overlaying soil characteristics, topography, land use, and management schemes. This model is able to calculate snowmelt proportions using a temperature index algorithm (TIA), and is a very efficient approach to analyzing the climate change impact on the streamflow in snow/glacier-dominated river basins. The SWAT model has the provision to divide each sub-basin into 10 elevation bands. The TIA approach was applied previously in the Yellow River basin [68], the Jhelum River Basin [36], and the Hunza River Basin [15] to simulate the streamflow in these mountainous watersheds. In this analysis, the TIA approach was utilized to estimate the snow/glacier meltwater. In the current research, each sub-basin in the Vakhsh River Basin was split into five elevation bands (EBs). The laps rate approach was used to estimate the precipitation, maximum temperature, and minimum temperature for each EB. The hydrological SWAT model is able to calculate evapotranspiration (ET) using the Penman-Monteith, Hargreaves, and Priestley-Taylor approaches, but we used the Hargreaves approach to evaluate the ET in the Vakhsh River Basin.
into five elevation bands (EBs). The laps rate approach was used to estimate the precipitation, maximum temperature, and minimum temperature for each EB. The hydrological SWAT model is able to calculate evapotranspiration (ET) using the Penman-Monteith, Hargreaves, and Priestley-Taylor approaches, but we used the Hargreaves approach to evaluate the ET in the Vakhsh River Basin.
For the current study, 19 sub-basins and 1048 HRUs in the Vakhsh River Basin were delineated using ArcSWAT 2012 (version 664). For model simulation, the daily measurement-based precipitation and maximum/minimum temperature data for a period of 14 years (2000-2013) were taken into consideration. The first three years (2000)(2001)(2002) were used as the warm-up period. The daily and monthly simulated streamflow was calibrated for the 2003-2008 time period and was then validated for the 2009-2013 period for the hydrological model, based on the flow data measured at the Darband hydrological station.

Sequential Uncertainty Fitting Version 2 (SUFI-2)
The semi-automatic calibration and validation approach was used in this study by applying the Sequential Uncertainty Fitting version 2 (SUFI-2) algorithm in the Soil and Water Assessment Tool-Calibration and Uncertainty Procedures (SWAT-CUP) tool package [69]. In the calibration, SUFI-2 expresses the errors and uncertainties during the hydrological modeling to a certain degree, as well as the ignorance and lack of understanding of the processes in natural systems. There is a close relationship between calibration and uncertainty. According to Abbaspour [70], without uncertainty, the calibration of the model is meaningless and deceptive. Moreover, any analysis with the calibrated model must include the uncertainty in the result based on the parameter uncertainties.  The semi-automatic calibration and validation approach was used in this study by applying the Sequential Uncertainty Fitting version 2 (SUFI-2) algorithm in the Soil and Water Assessment Tool-Calibration and Uncertainty Procedures (SWAT-CUP) tool package [69]. In the calibration, SUFI-2 expresses the errors and uncertainties during the hydrological modeling to a certain degree, as well as the ignorance and lack of understanding of the processes in natural systems. There is a close relationship between calibration and uncertainty. According to Abbaspour [70], without uncertainty, the calibration of the model is meaningless and deceptive. Moreover, any analysis with the calibrated model must include the uncertainty in the result based on the parameter uncertainties.

Elevation Bands
Taking into account the influence of the precipitation and temperature in mountainous regions, the division of the river basin into different elevation bands is regarded as a very important procedure. The SWAT model allows a basin to be divided into up to 10 elevation bands. In this study, the basin was divided into five elevation bands (EBs). Using the following mathematic equations, we adjusted the temperature and precipitation of each band.
where R band indicates the precipitation in an EB (mm); R day indicates the precipitation measured at the observed climate point (mm); EL band is the average elevation in the EB (m); EL c.s is the elevation at the climate point (m); plaps indicates the precipitation lapse rate (mm/km); days pcp,yr shows the mean annual number of precipitation days in the sub-basin; T max,band is the EB's maximum daily temperature (C); T min,band is the EB's minimum daily temperature (C); T mean,band shows the EB's mean daily temperature (C); T max is the maximum daily temperature measured at the climate station (C); T min is the minimum daily temperature measured at the climate station (C); T mean is the mean daily temperature measured at the climate station (C); tlaps shows the lapse rate of temperature (C/km); and 1000 is the conversion factor from meters to kilometers.

Calibration, Validation, and Sensitivity Analysis
In this study, the SWAT-CUP tool package was applied to calibrate the model parameters. This tool package is capable of offering two different sensitivity analyses, namely global and one-at-a-time sensitivity analyses. The one-at-a-time sensitivity analysis requires a separate sensitivity analysis of each parameter. In many cases, the sensitivity of one parameter depends on the other parameters. To facilitate work in a stationary computer and to avoid needing a more complicated sensitivity analysis to fix the range of values, in this calibration procedure, the global sensitivity analysis approach was used, which allowed us to analyze the sensitivity of all parameters in combination and to determine the key parameters in each hydrological sub-basin. Among the selected sensitive parameters, 14 were characteristically related to streamflow. By using the SUFI-2 model, the parameter sensitivity was evaluated based on the multiple regression system presented by Yang et al. [71] to estimate the effect of the selected value of the parameter (b i ) in the objective function (g). The following equation was used to calculate the parameter sensitivities.
The average changes of the parameter selected during the sensitivity analysis in the objective function could be computed using the above formula.

Assessment of the Hydrological Model's Performance
To evaluate the performance of the SWAT model, five indices were used: the coefficient of determination (R 2 ), the Nash-Sutcliffe efficiency (NSE) [72], the percent bias (PBIAS), the root-mean-square error standard deviation ratio (RSR), and the mean square error method (MSE). If the NSE and R 2 are both >0.5, then the SWAT model presents satisfactory outcomes; the range Water 2020, 12, 1426 9 of 34 of R 2 is −1 to 1 and that for NSE is 0 to 1. Less error and better performance of the simulated results are indicated by higher values of NSE and R 2 . The range of acceptable PBIAS values is −15 < PBIAS < +15; the best value is zero, and negative and positive values indicate model overestimation and underestimation, respectively. The RSR values range from zero to high positive values, and lower values of RSR show better model application [73]. The mathematical delimitations of the PBIAS, RSR, NSE, and R 2 objective functions are expressed below.
where n shows the entire number of sample pairs, Q obs i indicates the measured discharge station variable, Q obs indicates the average of the measured discharge station variables, Q sim i shows the simulated discharge variable, Q sim is the average of the simulated discharge variables, and i is the ith measured station or simulated data.

Additive and Multiplicative Methods of Change Factors
The five sets of GCM data  at the spatial resolution of 0.5 • were downloaded from the study domain www.2w2e.com and were downscaled using the Climate Change Toolkit (CCT); a detailed description of the CCT software is given in the supplementary material and can be downloaded along with the program from www.2w2e.com. The datasets from www.2w2e.com and CCT software have been used by many researchers, including Ashraf Vaghefi et al. [55] in America and Touseaf et al. [74] over Xijiang River Basin in China. In this study, under the two scenarios (RCPs 4.5 and RCP8.5), the meteorological variables, namely precipitation and maximum/minimum temperatures, were analyzed for two future time periods (2022-2060 and 2061-2099) with reference to the baseline condition, i.e., the 1966-2004 historical period. In order to obtain the future values of the current grid, the current variation was then coupled to the measured station variables. The mathematical formulation of additive change factor method (CFM) was used to evaluate the absolute change in the temperature of the GCM [75,76]. The observed values were multiplied by the change factor after calculating the arithmetic difference between the current and future GCM values. Characteristically, this method used for the precipitation of the GCM produces a reasonable evaluation of the value of a variable with respect to the relative change [75][76][77]. Figure 4 shows a schematic illustration of the method used for the evaluation of future scenarios using the multiplicative and additive change factors. The calculation of change factors may be applied to a single grid point of a GCM or to the average of several grid points. Equations (11) and (12) were used to estimate the monthly average values for GCM-based precipitation and maximum/minimum temperatures for the future scenarios and the baseline condition.
The GCMb i and GCMf i values show the GCM's future climate scenario and the GCM's baseline condition, GCMb and GCMf are the average values, and Nb and Nf present the number of variables at the temporal scale of the GCM's future climate scenario and the GCM's baseline condition, respectively.
Water 2020, 12, x FOR PEER REVIEW 10 of 34 The GCMb and GCMf values show the GCM's future climate scenario and the GCM's baseline condition, GCMb and GCMf are the average values, and Nb and Nf present the number of variables at the temporal scale of the GCM's future climate scenario and the GCM's baseline condition, respectively. To calculate the additive change factor (CF ) for the maximum/minimum temperatures and the multiplicative change factor (CF ) for the precipitation, Equations (13) and (14) were used.
The new locally scaled future additive (LSf , ) and the multiplicative (LSf , ) change factors from Equations (15) and (16) were obtained based on the values (CF ) and (CF ) [78].
In Equations (15) and (16), (LSf , ) and (LSf , ) represent the values of the variables of the future climate scenarios and were computed using the multiplicative and additive arithmetic formulation of the change factors. Lob presents the locally scaled observed values of the measured meteorological variable of a single measuring station or the mean measured time series for the ith time step. To calculate the additive change factor (CF add ) for the maximum/minimum temperatures and the multiplicative change factor (CF mul ) for the precipitation, Equations (13) and (14) were used.
The new locally scaled future additive (LSf add,i ) and the multiplicative (LSf mul,i ) change factors from Equations (15) and (16) were obtained based on the values (CF add ) and (CF mul ) [78].
In Equations (15) and (16), (LSf add,i ) and (LSf mul,i ) represent the values of the variables of the future climate scenarios and were computed using the multiplicative and additive arithmetic formulation of the change factors. Lob i presents the locally scaled observed values of the measured meteorological variable of a single measuring station or the mean measured time series for the ith time step.

Correlation between GCMs and Data Observed In Situ
With the aim of avoiding the systematic error between GCMs and data observed in situ (such as cumulative monthly precipitation and mean monthly temperature), the correlation of each downscaled GCM with the gauge-based data was assessed separately at the Rasht climate station in the Vakhsh River Basin for the 1966-2004 period.
The correlations between the gauge-based GCM data are represented in Figure 5. In each scatterplot of the GCMs, R 2 presents the coefficient of determination of the linear regression. From the scatterplot analysis, we found that the temperature correlation was considerably better than that of the precipitation ( Figure 5). This finding was consistent with those reported by Anjum et al. [52] for the Swat River Basin in Pakistan, by Bokke et al. [79] for the Nile River Basin in Ethiopia, and by Fang et al. [80] for the Kaidu River Basin in China. The correlation outcomes of all the applied GCMs (GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC, and NoerESM1-M) in the Vakhsh River Basin showed that the uncertainty of the precipitation was greater than that of the temperature. The GCM data had a comparatively stronger correlation with the observed data for the projected average monthly temperature. The average temperature showed almost identical accuracy to the maximum/minimum temperature. According to the execution of all models, the correlation coefficient between the results and the observations of the mean monthly temperature was higher than 0.99 for the 1966-2004 period. The correlation coefficient of the accumulated observed monthly precipitation and the GCMs data presented a value that was slightly higher than 0.95 for the 1966-2004 time period.

Changes in Annual Precipitation and Temperature
The mean annual temperature and precipitation projections compared to the baseline condition of two future time periods (2022-2060 and 2061-2099) under RCP4.5 and RCP8.5, coupled with the outputs of the five GCMs for the Vakhsh River Basin, are shown in Figure 6. The scatter plots in Figure 6 indicate the probable changes in the temperature of each GCM model and the GCM-based precipitation change rate in the Vakhah River Basin. All five selected GCMs (GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC, and NoerESM1-M) showed an increasing pattern of mean annual temperature under RCP4.5 and RCP8.5 in all future time periods. The three GCM simulations mostly showed a decreasing trend in mean annual precipitation over the basin, but two GCMs (MIROC and GFDL-ESM2M) displayed an increasing tendency in mean annual precipitation during the future time period 2022-2060 under RCP4.5. The pattern of annual average precipitation of the HadGEM2-ES model showed a significant increasing trend compared to the GFDL-ESM2M model, and the other three GCM models indicated a decreasing trend during the first future time period (2022-2060). Among the five GCMs, the MIROC model showed a significant increasing trend and the IPSL-CM5A-LR model showed a significant decreasing trend in annual precipitation over two future time periods, 2022-2060 and 2061-2099, under RCPs 4.5 and 8.5. The current findings indicate that the probable mean annual precipitation varied in the range of uncertainty. Similar results were reported by Anjum et al. [52] for the Swat River Basin in Pakistan and by Zhang et al. [32] for the Xian River Basin in China. Regarding the mean annual precipitation, the range of variation generally decreased, and the mean annual temperature increased over time. Furthermore, the multi-model ensemble (MME) predicted a decreasing trend in mean annual precipitation (from −4.46% to −7.42%) during the two future time periods under both RCPs in the Vakhsh River Basin.  the probable mean annual precipitation varied in the range of uncertainty. Similar results were reported by Anjum et al. [52] for the Swat River Basin in Pakistan and by Zhang et al. [32] for the Xian River Basin in China. Regarding the mean annual precipitation, the range of variation generally decreased, and the mean annual temperature increased over time. Furthermore, the multi-model ensemble (MME) predicted a decreasing trend in mean annual precipitation (from −4.46% to −7.42%) during the two future time periods under both RCPs in the Vakhsh River Basin. The projection of the climate variables under RCP4.5 displayed alterations in the mean annual temperature and precipitation during the future time periods 2022-2060 and 2061-2099, where the precipitation could vary from −22.1% to 6.5%, and the temperature could rise from 1.6 to 2.5 °C. Hence, for the second future period, the precipitation may vary from −16.0% to −1.7% and temperature may rise from 1.9 to 3.1 °C. Under RCP8.5, during the 2022-2060 future time period, the average annual precipitation is expected to decrease by 21%-0.05%, and the temperature may rise The projection of the climate variables under RCP4.5 displayed alterations in the mean annual temperature and precipitation during the future time periods 2022-2060 and 2061-2099, where the precipitation could vary from −22.1% to 6.5%, and the temperature could rise from 1.6 to 2.5 • C. Hence, for the second future period, the precipitation may vary from −16.0% to −1.7% and temperature may rise from 1.9 to 3.1 • C. Under RCP8.5, during the 2022-2060 future time period, the average annual precipitation is expected to decrease by 21-0.05%, and the temperature may rise from 1.8 to 3.7 • C; in the case of the 2061-2099 time period, the annual precipitation could vary from −29.8% to 5.3% and temperature may rise from 3.4 to 6.4 • C. Comprehensively, under RCP8.5, temperature and precipitation depicted a wide range of uncertainties, which is shown in Figure 6.  The projected maximum temperature (Tmax) and minimum temperature (Tmin) showed a coherent rising tendency in all seasons. The probable changes in Tmax and Tmin under RCP4.5 in both future time periods are expected to vary from 1.5 to 3.4 °C and from 1.5 to 3.1 °C in winter and from 1.8 to 3.6 °C and from 1.5 to 3.1 °C in spring, respectively. Under the same RCP scenario, the expected changes in Tmax and Tmin in both time periods in the summer season are 1.8-5.4 °C and 1-3.3 °C and in the fall season are 2.1-4.5 °C and 1.2-4 °C, respectively. Under RCP8.5, the respective probable increases in Tmax and Tmin in both time periods range from 2.7 to 6.8 °C and from 2.7 to 6.8 °C (winter) and from 2.6 to 7.1 °C and 2.1 to 6.1 °C (spring) compared with the range from 2.9 to 8.9 °C and from 1.1 to 6.8 °C (summer) and from 3.4 to 7.9 °C and from 1.3 to 7.1 °C (fall).

Changes in Future Seasonal Maximum/Minimum Temperature and Precipitation
The seasonal changes in the average maximum/minimum temperature and precipitation during the future time periods of 2022-2060 and 2061-2099, with reference to the baseline condition, are presented in Table 2 to 4.5% (HadGEM2-ES) in the respective future periods. In the summer season, under both RCPs, the decreasing tendency of mean precipitation was greater than in the other seasons; it decreased from 26.9% to 2.5%, and might vary from −57.1% to 26.4% at the end of this century. In the 2022-2060 time period under RCP4.5 and RCP8.5 in the summer, the average precipitation changed from −50.3% to 19.7% and −47.8% to 4.2%, respectively. Under both RCPs in the fall season, the mean precipitation for one GCM (MIROC) displayed a complete increase during the two future periods, and another GCM (NoerESM1-M) showed an increasing tendency (from 7.1% to 11%) in the first future period under RCP8.5 and in the second future period under RCP4.5; however, a decreasing trend (from 8.2% to 10.3%) was found in the 2022-2060 time period under RCP8.5 and in the 2022-2060 time period under RCP4.5. Furthermore, three GCMs (GFDL-ESM2M, HadGEM2-ES, and IPSL-CM5A-LR) generally indicated a decreasing trend in all of the future time periods. According to these three models, the respective mean precipitation in the fall season is expected to decrease from 8.6% to 23.2%, 10.4% to 25.4%, and 11.2% to 32%, respectively, and could vary from −41.2% to 5.1% over both time periods. The seasonal scale analysis of the mean precipitation changes highlighted large uncertainties among the five GCMs in the various future time periods under RCP4.5 and RCP8.5, as evidenced by the presence of both increasing and decreasing tendencies. However, the MME predicted an increase in winter precipitation in the 2022-2060 and 2061-2099 time periods (from 0.42% to 5.1% under RCPs 4.5 and 8.5) in the Vakhsh River Basin in Central Asia. The current findings in line with results of Omani et al. [13] and Ta et al. [81] where they used the CMIP5 MME for 16 GCMs and CMIP5 MME for 37 GCMs. Omani et al. [13] and Ta et al. [81] reported that winter precipitation is expected to increase in the Pamir-Alay mountain range in Central Asia. Previous studies have reported that both tropical and westerly circulation modes simultaneously impact the Central Asian climate, especially throughout the moist winter and spring seasons. In Central Asia, the weather types represent a regional expression of signs of tropical and extratropical circulation, which spread from the western and southern directions. Furthermore, previous studies indicate that hydroclimate conditions over Central Asia are strongly impacted by the location and strength of the westerly jet stream [82,83].
The projected maximum temperature (Tmax) and minimum temperature (Tmin) showed a coherent rising tendency in all seasons. The probable changes in Tmax and Tmin under RCP4.5 in both future time periods are expected to vary from 1.5 to 3.4 • C and from 1.5 to 3.1 • C in winter and from 1.8 to 3.6 • C and from 1.5 to 3.1 • C in spring, respectively. Under the same RCP scenario, the expected changes in Tmax and Tmin in both time periods in the summer season are 1.8-5.4 • C and 1-3.3 • C and in the fall season are 2.1-4.5 • C and 1.2-4 • C, respectively. Under RCP8.5, the respective probable increases in Tmax and Tmin in both time periods range from 2.7 to 6.8 • C and from 2.7 to 6.8 • C (winter) and from 2.6 to 7.1 • C and 2.1 to 6.1 • C (spring) compared with the range from 2.9 to 8.9 • C and from 1.1 to 6.8 • C (summer) and from 3.4 to 7.9 • C and from 1.3 to 7.1 • C (fall).

Hydrological Modeling Results
The simulation of the river discharge generated by applying the hydrological model is shown in Table 3. In this study, the SWAT-CUP tool was used to examine the sensitivity of the parameters, the uncertainty analysis, and the model calibration. The SWAT-CUP tool integrates procedures for several algorithms (Figure 3). In this study, the Sequential Uncertainty Fitting version 2 (SUFI-2) algorithm was selected to find parameters for baseflow simulation (GW QMN, GW_REVAP, REVAPMN, etc.), parameters for surface flow simulation (CN2, SOL_AWC, ESCO, etc.), and parameters for flow process line adjustment (ALPHA_BF, RCHRG_DP, etc.) for analysis based on experience and references [84,85]. Based on the statistical criteria of the t-stat and p-value, the sensitivity of each parameter and its relative significance were evaluated. The t-stat was computed from the coefficient of a variable in the multiple regression analysis divided by the standard error. A parameter was considered to be sensitive if the value of the coefficient was large in comparison to the standard error. The p-value was computed via comparison of the t-stat value with the Student's distribution table [86]. Table 3 shows the t-stat and p-values of the different effective parameters. In the hydrological modeling, the calibration was checked along with the sensitivity of the model parameters such as snow, groundwater recharge, soil, evaporation, and basin slope. In the first step, the model was calibrated without considering the different elevation bands for the 14 year period (2000-2013). The first three years of the period were selected as the warm-up period for the hydrological model. The period from 2003 to 2008 was considered the calibration period, and from 2009 to 2013 was the model validation period. The parameter ranges under different objective functions were adjusted to obtain better results; however, in this mountainous study area, the SWAT-CUP did not provide satisfactory results based solely on modifying the parameter ranges. Initially, the monthly discharge from 2003 to 2013 time period was calibrated, without applying the elevation band approach for the Darband discharge station. The uncertainty performance of the p-factor, r-factor, and the values of R 2 and NS were 0.48, 0.72, 0.41, and 0.04, respectively. It was essential to split the data into elevation bands in the mountainous watershed before fixing the range of parameters in the calibration procedure. We improved the simulation results after zoning (five elevations), as shown in Figures 8 and 9 and Table 4.
The model calibration and validation results, as well as a comparison of the observed and simulated flows at daily and monthly time intervals, are described in Table 4. The quantification of the goodness of calibration and the uncertainty performance were derived from the p-factor and r-factor. For discharge, the recommended values for the r-factor of <1.5 and p-factor of >0.7 or 0.75 are acceptable, depending on the sufficiency of the input calibrating data and the scale of the project [84]. In this study, the p-factor and the r-factor values were 0.87 and 1.15, and 0.86 and 0.83 for the calibration of the daily and monthly time steps, respectively. For the validation period, the values of the r-factor and the p-factor were 0.82 and 0.84, and 0.86 and 1.08 for the daily and monthly simulations, respectively. Statistical coefficients such as the p-factor and the r-factor showed that during the calibration and validation periods at the daily and monthly time scales, the semi-distributed SWAT model generated good statistical results. As expected, the model performed slightly better in the calibration period than in the period of validation. The calibration and validation results were slightly better for the monthly time scales than for the daily time scales. A comparison of the observed and simulated discharge and the 95% prediction uncertainty (95PPU) bandwidth indicated that the well-calibrated hydrological model could effectively simulate monthly and daily streamflow under different RCPs in the Vakhsh River Basin.
Water 2020, 12, x FOR PEER REVIEW 18 of 34 the monthly discharge from 2003 to 2013 time period was calibrated, without applying the elevation band approach for the Darband discharge station. The uncertainty performance of the p-factor, rfactor, and the values of R 2 and NS were 0.48, 0.72, 0.41, and 0.04, respectively. It was essential to split the data into elevation bands in the mountainous watershed before fixing the range of parameters in the calibration procedure. We improved the simulation results after zoning (five elevations), as shown in Figures 8 and 9 and Table 4. The model calibration and validation results, as well as a comparison of the observed and simulated flows at daily and monthly time intervals, are described in Table 4. The quantification of the goodness of calibration and the uncertainty performance were derived from the p-factor and rfactor. For discharge, the recommended values for the r-factor of <1.5 and p-factor of >0.7 or 0.75 are acceptable, depending on the sufficiency of the input calibrating data and the scale of the project [84]. In this study, the p-factor and the r-factor values were 0.87 and 1.15, and 0.86 and 0.83 for the calibration of the daily and monthly time steps, respectively. For the validation period, the values of the r-factor and the p-factor were 0.82 and 0.84, and 0.86 and 1.08 for the daily and monthly simulations, respectively. Figure    Statistical coefficients such as the p-factor and the r-factor showed that during the calibration and validation periods at the daily and monthly time scales, the semi-distributed SWAT model generated good statistical results. As expected, the model performed slightly better in the calibration period than in the period of validation. The calibration and validation results were slightly better for the monthly time scales than for the daily time scales. A comparison of the observed and simulated discharge and the 95% prediction uncertainty (95PPU) bandwidth indicated that the well-calibrated hydrological model could effectively simulate monthly and daily streamflow under different RCPs in the Vakhsh River Basin.  Table 4. The statistics of a comparison of the observed monthly and daily calibration and validation data for the SWAT model with the simulation band based on the p-factor and r-factor at the Darband discharge gauging station in the Vakhsh River Basin (VRB), used to estimate the model performance. NSE: Nash-Sutcliffe efficiency; R 2 : coefficient of determination; PBIAS: percentage bias; MSE: mean square error; RSR: root-mean-square error standard deviation ratio.

Factors and Statistical Indices
Monthly Daily

Future Changes in Extreme Discharges
Prediction of the streamflow as well as of extreme flows is of great importance for flow control in the Vakhsh River. The Vakhsh River serves as a water source for agriculture, municipal uses, and hydropower production in Central Asian countries. We used a SWAT model to predict streamflow and extreme flows in the Vakhsh River. Figure 10a Table 6. The high flow (Q5) in the two future time periods may vary from 30.9% to 53.3% and from 41.3% to 59.2% under RCP4.5, and under RCP8.5 high flow may vary from 40.1% to 49.5% and from 41.4% to 78%. For all five GCMs, the low flow (Q95) in the two future time periods may decrease from 51.4% to 57.2% and from 41% to 55.7% under both RCPs. However, the results showed that the changes in the median flow (Q50) were slightly different from the changes in extreme flows ( Table 6). The analysis of the median flow in the 2022-2060 time period presented a decreasing behavior for all five GCMs under both RCPs. The median flow might decrease from 10.8% to 22.1% and from 6% to 24.1% under RCP4.5 and RCP8.5, respectively, in the 2022-2060 time period. It is expected that the median flow might vary from −20.9% to 15.1% under RCP4.5 and might increase from 16.3% to 33.4% under RCP8.5 in the 2061-2099 time period. We revealed from the analysis of Q5, Q50, and Q95 compared to the low and median flows that the increase in the high flow could be greater. Similarly, the probable decrease in the low flow was greater than that in the median flow. We projected a significant increasing trend in high flow in Vakhsh River, which can be attributed to the rapidly melting snow and ice caused by constantly increasing air temperature. Increased high flow in the Vakhsh River implies a potential risk of flooding in the Vakhsh River Basin. Another phenomenon is that with increasing air temperatures, snowfall will decrease and rainfall will increase [13]. These variations may result in increased flooding during the early spring. Decreased low flow in the Vakhsh River  Table 6. The high flow (Q5) in the two future time periods may vary from 30.9% to 53.3% and from 41.3% to 59.2% under RCP4.5, and under RCP8.5 high flow may vary from 40.1% to 49.5% and from 41.4% to 78%. For all five GCMs, the low flow (Q95) in the two future time periods may decrease from 51.4% to 57.2% and from 41% to 55.7% under both RCPs. However, the results showed that the changes in the median flow (Q50) were slightly different from the changes in extreme flows ( Table 6). The analysis of the median flow in the 2022-2060 time period presented a decreasing behavior for all five GCMs under both RCPs. The median flow might decrease from 10.8% to 22.1% and from 6% to 24.1% under RCP4.5 and RCP8.5, respectively, in the 2022-2060 time period. It is expected that the median flow might vary from −20.9% to 15.1% under RCP4.5 and might increase from 16.3% to 33.4% under RCP8.5 in the 2061-2099 time period. We revealed from the analysis of Q5, Q50, and Q95 compared to the low and median flows that the increase in the high flow could be greater. Similarly, the probable decrease in the low flow was greater than that in the median flow. We projected a significant increasing trend in high flow in Vakhsh River, which can be attributed to the rapidly melting snow and ice caused by constantly increasing air temperature. Increased high flow in the Vakhsh River implies a potential risk of flooding in the Vakhsh River Basin. Another phenomenon is that with increasing air temperatures, snowfall will decrease and rainfall will increase [13]. These variations may result in increased flooding during the early spring. Decreased low flow in the Vakhsh River implies that the risk of drought may also increase during fall season due to less meltwater contributing to runoff in the basin.

Discharge and Snowmelt Changes under Future Scenarios and Shifts in the Peak Flows
Snow melting is essential in mountainous snow-covered basins. Figure 11 shows the changes in snowmelt for the 2020-2060 and 2061-2099 time periods under RCP4.5 and RCP8.5 for the five GCMs. In the two future time periods under both RCPs, an increasing snowmelt tendency was indicated. However, the significance of an increasing snowmelt trend varied from month to month. For instance, from October to June, the snowmelt presented a slightly increasing trend under both RCPs in both future time periods. In July and August, snowmelt may increase significantly compared with the other months in both time periods under both RCPs. In the Vakhsh River Basin, the months of March and April are warmer than those of May and June; however, snowmelt had a stronger increasing tendency in March and April than in May and June for all five GCMs under both RCPs in the two future time periods. Figure 11 shows that the Vakhsh River flow rate has great seasonal variability between winter and summer. It indicates that the Vakhsh is a snow-and glacier-fed river. Omani et al. [13] projected that under RCP8.5, at least 41% to 2% of the Vakshsh River Basin glaciers will remain. The multi-model ensemble under RCP4.5 and RCP8.5 scenarios projected that across the Vakhsh River Basin, by 2100, no glaciers shallower than 65 and 98 m.w.e will survive [13]. We conclude that increasing snowmelt in the future may be caused by the continuous increasing trend of air temperatures in the region, which will influence the water cycle and flow regime in the Vakhsh River.
In the two future time periods, by considering the average monthly discharge, the peak flow was estimated in the Vakhsh River under the RCP4.5 and RCP8.5 scenarios (Figure 12). The mean monthly discharge in the first time period showed an increasing trend in the months of April to June and a slight decreasing tendency in the months of August to May for all five GCMs under RCP4.5 and RCP8.5. Under both RCPs, in the 2061-2099 future time period, the mean monthly discharge increased for all five GCMs in the months February to July. The mean monthly discharge mostly showed a decreasing trend in the months of August to January under both RCPs in the second future time period. Under RCP4.5 in the first future time period, the mean monthly peak discharge was 1978 m 3 /s in In the two future time periods, by considering the average monthly discharge, the peak flow was estimated in the Vakhsh River under the RCP4.5 and RCP8.5 scenarios (Figure 12). The mean monthly discharge in the first time period showed an increasing trend in the months of April to June and a slight decreasing tendency in the months of August to May for all five GCMs under RCP4.5 and RCP8.5. Under both RCPs, in the 2061-2099 future time period, the mean monthly discharge increased for all five GCMs in the months February to July. The mean monthly discharge mostly showed a decreasing trend in the months of August to January under both RCPs in the second future time period. Under RCP4.5 in the first future time period, the mean monthly peak discharge was 1978     We revealed from this analysis that in the two future time periods (2022-2060 and 2061-2099) under RCP4.5 and RCP8.5, compared to the baseline condition, January was the month with the lowest mean monthly discharge. In general, discharge exhibited an increasing tendency, which suggested that for downstream regions, more water resources will be available in the future. Discharge increased in spring and summer seasons in all future time periods; similarly, a significant decreasing tendency was indicated in the months of August and September under both RCPs. Under both RCPs and in both future time periods, the monthly peak discharge could shift from the month of July to June for all GCMs (mean monthly peak discharge for the baseline). Moreover, the mean monthly peak discharge could shift to late spring (month of May) for two GCMs (IPSL-CM5A-LR and NoerESM1-M) in the 2061-2099 time period compared to the baseline mean monthly discharge in the months of July. In this scenario, the streamflow increases and the peak flow shifts earlier because of earlier snowmelt caused by global warming. We revealed from this analysis that in the two future time periods (2022-2060 and 2061-2099) under RCP4.5 and RCP8.5, compared to the baseline condition, January was the month with the lowest mean monthly discharge. In general, discharge exhibited an increasing tendency, which suggested that for downstream regions, more water resources will be available in the future. Discharge increased in spring and summer seasons in all future time periods; similarly, a significant decreasing tendency was indicated in the months of August and September under both RCPs. Under both RCPs and in both future time periods, the monthly peak discharge could shift from the month of July to June for all GCMs (mean monthly peak discharge for the baseline). Moreover, the mean monthly peak discharge could shift to late spring (month of May) for two GCMs (IPSL-CM5A-LR and NoerESM1-M) in the 2061-2099 time period compared to the baseline mean monthly discharge in the months of July. In this scenario, the streamflow increases and the peak flow shifts earlier because of earlier snowmelt caused by global warming.

Discussion
The Vakhsh River is the second largest tributary of the Amu Darya River after the Panj River. The climate of the basin is continental under the influence of the westerly wind, leading to significant seasonal changes in temperature and precipitation [87]. These changes are due to the mountainous topography of the basin, which has a very high local contrast [18,88]. The elevation ranges from 302 to 7050 meters above sea level, and mountains occupy nearly 90% of the basin's area; it is the most vulnerable territory in Central Asia to hydroclimatic changes [89]. Furthermore, the world's largest mountain systems neighbor Tajikistan-the Himalayas on the north-east and the Tian Shan on the

Discussion
The Vakhsh River is the second largest tributary of the Amu Darya River after the Panj River. The climate of the basin is continental under the influence of the westerly wind, leading to significant seasonal changes in temperature and precipitation [87]. These changes are due to the mountainous topography of the basin, which has a very high local contrast [18,88]. The elevation ranges from 302 to 7050 m above sea level, and mountains occupy nearly 90% of the basin's area; it is the most vulnerable territory in Central Asia to hydroclimatic changes [89]. Furthermore, the world's largest mountain systems neighbor Tajikistan-the Himalayas on the north-east and the Tian Shan on the south-east. In addition, westward Tajikistan is open to wet Caspian and Mediterranean winds in winter. In summer, the south-western parts of the basin are strongly affected by the dry heat waves from the deserts of Uzbekistan, Turkmenistan, and Afghanistan.
In the mountain region of the Central Asia, the large amount of precipitation leads to a substantial accumulation of snow during the winter and spring season. The stream flow increases rapidly with the beginning of the snow melting period in March and April. The Vakhsh River Basin receives its maximum precipitation in winter and spring, which results in a clear temporal separation in accumulation of the snow and the time of the peak flow. The Vakhsh River is characterized by a seasonal cycle of river flow with maximum flow in the summer season ( Figure 2). Summer runoff of the upstream catchments of the Central Asian rivers is controlled by the melting of snow and glaciers, which, in glacierized catchments, contributes up to 50% of the seasonal flow [90].
The downscaled precipitation and maximum/minimum temperature estimates from the five GCMs were in good agreement with the measurement-based data under the baseline condition . The product of multi GCMs of the Inter-Sectoral Impact Model Intercomparison Project (ISI-MIP) was used to assess the influences of the climate variation on the hydrological regimes, considered a suitable representative [2,91,92]. The meteorological data from ten climate stations and streamflow data from one hydrological station were used in this study. We analyzed the possible consequence of projected climate variation on interseasonal and annual water flow, as well as snowmelt in the VRB. We studied the responses of hydrological processes such as discharge and extreme and median flows to climate variation. The temporal impacts of the projected climate change were assessed by coupling a well-calibrated semi-distributed hydrological model (SWAT) with the results of the five GCMs under two greenhouse gas emission scenarios (RCP4.5 and RCP8.5). The effect of topography was corrected by applying the elevation band (temperature lapse rate and precipitation lapse rate) approach in this study, which improved the simulation results. We presented the results of the streamflow calibration and validation for the Vakhsh River at daily and monthly time scales. Additionally, by using the SWAT-CUP tool, we studied the sensitivity and prediction uncertainty of the model parameters, which was necessary to evaluate the strength of the calibrated model.
The results of the five GCMs presented continuous warming over the VRB at the annual and seasonal scales in the middle and at the end of this century. This finding is consistent with the results of previous studies in the Central and South Asian regions, including the Pamir-Alay [13,93], the Tian Shan-Pamir-North Karakoram [20], the Himalayas [21], and the Tibetan Plateau [22].
The probable change in water flow (presented in Table 6) indicated that the imbalances in the future water resources will be more severe. It is expected that the wet and dry seasons in the basin will become more severe than those in the baseline period. We showed that these phenomena are predicted to be stronger under RCP8.5 than under RCP4.5. A possible reason for these alterations is that in the future, more meltwater will be produced in early summer and more snow will be replaced by rain. The typical changes will accelerate the convergence of water flows and raise the flooding frequency and intensity. One of the possible causes of the rising temperature upstream of the Amu Darya River Basin could be the rising concentration of aerosols and greenhouse gases in the regional atmosphere. Folini et al. [94] reported that aerosol emissions in the 20th century might increase, in association with the enormous population and industrialization growth. Similarly, Bollasina et al. [25] confirmed that in Asia, the concentration of atmospheric aerosols has increased steadily. Xin et al. [26] reported that over China and Central Asia, a rising trend in the concentration of aerosols in the atmosphere could cause a substantial rise in temperature. The mean annual precipitation over the VRB is expected to rise in the 2022-2060 future time period under RCP4.5, as evidenced by the results of two GCMs (GFDL-ESM2M and MIROC), as well as the two GCMs (HadGEM2-ES and MIROC) that indicated a rising trend of precipitation in the second time period (2061-2099) under RCP8.5. The remaining GCM models that were analyzed in this study showed a decreasing trend of mean annual precipitation in the two future periods under both RCPs. The winter mean precipitation had a lower decreasing trend than that of the summer and fall seasons. The summer mean precipitation exhibited a greater decreasing tendency than the other seasons in the two future periods. Similarly, in the Yellow and Xin River Basins in China [32,68], in the Middle East [33], and in the westerly-dominated region of South Asia, a decreasing trend of summer precipitation was indicated [34]. The increasing or decreasing propensity of winter precipitation varies from model to model. Two GCMs showing a rising trend of winter precipitation and resembling this analysis were reported by Luo et al. [35] for the Heihe River Basin and by Omani et al. [13] for the Pamir-Alay mountains in Central Asia. However, the patterns of seasonal variations in precipitation for three GCMs presented in the Vakhsh River Basin are contrary to those in reported for the Hunza River Basin of the Karakoram mountains [15] and the Jhelum River Basin of the Himalayan mountains [36]. Pendergrass et al. [95] reported that the global winter precipitation increased over the second half of the 20th century, and they attributed this to the role of increasing moisture counteracted by weakening circulation. Li et al. [27] pointed out that in Central Asia, at the end of 20th century, there was a persistent decreasing trend of annual precipitation, and Meng et al. [29] confirmed that precipitation might increase in the middle of the 21st century in the south of the Tian Shan mountains. The multi-model ensemble result projected a decrease in average annual precipitation during the 2022-2060 and 2061-2099 time periods under both RCPs in the VRB. Hence, the MME projected an increase in winter precipitation in the 2022-2060 and 2061-2099 time period under RCPs 4.5 and 8.5. The interannual and seasonal scale analyses of the mean precipitation changes presented large uncertainties among the GCMs, as evidenced by both increasing and decreasing tendencies under RCP4.5 and RCP8.5 in various future time periods. In the Asian region, these contradictions in our findings may be associated with the rising concentration of anthropogenic absorbing aerosols and the westerlies system [25,94,96,97].
The analysis of the simulated Vakhsh River flows mostly described an increasing trend of the mean annual streamflow during the 2061-2099 time period under RCP4.5 and RCP8.5, as well as a decreasing trend during the 2022-2060 time period under RCP4.5 ( Figure 13). Mostly, the increase and a lesser decrease of annual future flow may be attributed to the similar projection of the different GCM models for the total annual precipitation. Similarly, a decreasing/increasing tendency for the different future periods of annual river flow was indicated in an arid alpine catchment in Karakoram [98]. A study of the extreme flows such as high flow (Q5) and low flow (Q95), as well as the median flow (Q50) revealed that the probable decrease in the low flow is larger than that in the median flow, and a rise in the high flow might be larger compared to the median and low flows. In addition, it is expected that the average monthly peak discharge may shift to earlier in the summer season, from July to June, for the two future time periods under both RCP4.5 and RCP8.5 for almost all GCMs, which is mainly due to the slight rise in precipitation in the spring and winter seasons, as well as because of an earlier snowmelt caused by global warming. Similarly, due to earlier snowmelt, Siegfried et al. [99] and Sorg et al. [100] projected the impacts of climate changes on flow seasonality and concluded that less water will be available in the summer months in the Syr Darya River Basin in Central Asia. Olsson et al. [101] confirmed, from trend analysis, this shift in seasonality of flow and predicted a possible decreasing trend of annual flow in the Zarafshan River Basin in Central Asia. A similar shift in the discharge peak (July to June) was pointed out by Liu et al. [98] for the Yarkant River Basin in Central Asia; in contrast, Babur et al. [36] reported that the discharge peak could be delayed (July to August) in the Jhelum River Basin. These discrepancies in detections might be related to various projected climate models for in the winter and summer seasonal precipitation in Central Asia. We found that the average monthly peak discharge in the Vakhsh River Basin indicated a significant decreasing tendency in August and September for the 2022-2060 and 2061-2099 time periods under RCP4.5 and RCP8.5 for all five of the GCMs. In this study, most of the GCM model outputs along with the multi-model ensemble results showed that the summer water flow in the Vakhsh River is expected to increase at the end of the 21st century under the two studied greenhouse gas emissions scenarios. Increases in summer water flow in the Vakhsh River can be ascribed to the rapid melting snow and ice caused by continuously increasing air temperature. Meanwhile, the projection of summer water availability is essential for the downstream countries, such as Uzbekistan and Turkmenistan in Central Asia, to adapt the national agricultural strategies to anomalous hydroclimatic conditions.

Conclusions
In this study, the semi-distributed hydrological SWAT model was applied in two future time periods (2022-2060 and 2061-2099) under the RCP4.5 and RCP8.5 scenarios in order to evaluate the effects of probable climate variation on the Vakhsh River flow. For the current basin, in order to predict climate variation, the results of the five GCMs were downscaled. The hydrological model was calibrated using a semi-automated SWAT-CUP tool and was validated for the 2003-2008 and 2009-2013 periods at daily and monthly time scales. We assessed the strength of the calibrated model by applying sensitivity and uncertainty analyses. Additionally, this mountainous catchment was divided into five elevation zones, which enhanced the accuracy of simulated outflows. The results of the values of the p-factor, the r-factor, and the statistical indices (NSE, R 2 , RSR, MSE, and PBIAS) for both the calibration and validation periods indicated that the hydrological model, which was built for the VRB, was quite well calibrated and reliable for use for the prediction of the possible influences of the climate variation on the flows. The main findings of this work were as follows.
(1) Based on the values of uncertainty and statistical evaluation indices of the simulated streamflow, it is concluded that under altering climatic conditions in the Vakhsh River Basin in Central Asia, the hydrological SWAT model is reliable to simulate the streamflow. . From the aspect of seasonal variation, under both RCPs, the MME indicated a decreasing trend of the fall and winter flows in the near future time period (from −24.06% to −24.36% and from −11.25% to −13.92%). However, in fall and winter seasons, flows are expected to increase in the end of the 21st century under both RCPs (3.75% and 12.75%). The MME revealed an increasing trend in summer and spring flows during the future time periods 2022-2060 and 2061-2099 under both RCPs 4.5 and 8.5 (from 1.72% to 2.65% and from 30.59% to 95.69%). In this study, uncertainty in flow simulation existed because we treated the ice melt of glaciers as snowmelt in the SWAT model. In the GCM models, uncertainties also exist, which are propagated into SWAT. The streamflow, snowmelt simulation, and results description could be influenced by GCM uncertainties. (5) Analysis of the flow duration curves revealed, for all GCMs, a possible increase in the high flows projected under both RCPs for the two future time periods. However, in the basin, the low flow was projected to decrease under RCP4.5 and RCP8.5 in the 2022-2060 and 2061-2099 time periods compared to the baseline condition . We found that the high flow was projected to increase more strongly in the future compared to the median and low flows. The possible decrease in the low flow was higher than the decrease in the median flow. (6) It is expected that snowmelt might increase continuously with increasing temperature, and the average monthly peak discharge in the Vakhsh River Basin might shift to earlier in the summer season, from July to June, while a significant decreasing tendency in the average monthly peak discharge was found in August and September for the two future time periods under both the RCP4.5 and RCP8.5 scenarios.
The findings of the current research might be useful for the formulation of new water resource management guidelines and adaptation strategies to climate change in current ungauged mountainous catchments. Additionally, these findings might be helpful for making new planning schemes for irrigation in downstream countries, such as Uzbekistan and Turkmenistan, as well as for the sustainable design and administration of water resources for constructing hydropower plants in upstream countries, such as Tajikistan and Kyrgyzstan.
Author Contributions: All authors were involved in the intellectual part of this paper. X.C. and T.L. designed and supervised this research. A.G. conducted the major data processing and modeling works and wrote the manuscript. N.G. helped in the data arrangement. M.N.A. and M.R. revised the manuscript structure. All authors read and approved the final manuscript.