Climate Change Impact on the Hydrology of a Typical Watershed in the Tianshan Mountains

1State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Urumqi 830011, China 2University of Chinese Academy of Sciences, Beijing 100049, China 3Department of Geography, Ghent University, 9000 Ghent, Belgium 4National Institute of Water and Atmospheric Research, Christchurch 8011, New Zealand 5Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China


Introduction
The Tianshan Mountains, regarded as the "water tower of Central Asia" [1], are located in the innermost center of the Eurasian continent.The long distance to the surrounding oceans causes a dry climate, especially for the surrounding basins.Rivers starting in the mountainous regions provide agricultural and domestic water for the surrounding basins and oases.With their distinctive topographic and landscape features, the Tianshan Mountains show a unique energy balance and hydrological cycle and are expected to be sensitive to climate change [2,3].
Many reports show a widespread climatic and hydrologic change in the Tianshan Mountains during the past few decades [4].For example, temperature demonstrated a significant rising trend (significant level is smaller than 0.001) at a rate of 0.33∼0.34∘ C/decade during 1960∼2010, which is higher than China (0.25 ∘ C/decade) and the entire globe (0.13 ∘ C/decade) [5,6]; precipitation increased substantially in most regions especially for the middle and high latitudes; glacier area decreased by 11.5% and the thickness of snowpack has also decreased [5,7].Pan evaporation and wind speed have also changed [3].The annual runoff increased as well, for example, for the Urumqi River, the Kaidu River, and the Aksu River [8,9].
Future changes in the streamflow and watershed hydrology have become increasingly important to water resource management in the Tianshan Mountains.However, only a limited number of studies currently focus on impact of future climate change on hydrology; for example, Sorg et al. [1] indicated that the total runoff is likely to remain stable or even increase slightly in the near future but it will decrease 2 Advances in Meteorology at the end of the 21st century for Central Asia.There are also researches demonstrating that the annual runoff will decrease to some extent in the first half of the 21st century [2,10].Previous studies seldom address implications of climate change on the hydrological cycle and hydrological components (e.g., ET, surface flow, and groundwater).To complement these studies, this paper aims at understanding the future hydrological system and assessing the responses of the hydrologic system to climate change.
In the present study, two sets of climatic data, that is, RCM outputs and SCC data, are used to force SWAT [11] and are applied to the Kaidu River Basin, a typical watershed on the south slope of the Tianshan Mountains, to assess future changes of the hydrologic cycle and the hydrological effects of changes in climate variables.Questions that are addressed include the followings.(1) How will the future climate and hydrology change in this region?(2) What is the effect of climate change on the hydrologic cycling?(3) What is the difference between simulations with RCM outputs and SCC? Understanding these issues will enable assessing the future hydrological change and its unique hydrometeorological processes better.The remaining is constructed as follows: Section 2 introduces the study area and data; Section 3 describes the hydrological model, bias correction methods, and analysis procedures; Section 4 provides the results and discussion, followed by conclusions in Section 5.

Study Area and Data
The Kaidu River Basin, with a drainage area of 18,634 km 2 above the Dashankou hydrological station, is considered to be a very typical watershed in the Tianshan Mountains (Figure 1).Its altitude ranges from 1,342 m to 4,796 m above sea level (asl) with an average elevation of 2,995 m.The climate here is temperate continental with alpine climate characteristics and obvious seasonal variation.This river provides water resources for agricultural activity and the ecological environment of the oasis, with an area of over 70,000 km 2 and a population of over 1.15 million, which is mainly stressed by water scarcity [12].Therefore, projecting the impact of future climate change on water resources is urgent for the sustainable development of this region and it also provides information on the implications of climate change on the water tower in Central Asia.
The daily observed meteorological data (including precipitation, maximum/minimum temperature, wind speed, and relative humidity of two meteorological stations Bayanbulak and Baluntai) from 1970 to 2005 are from the China Meteorological Data Sharing Service System (http://www .cma.gov.cn/2011qxfw/2011qsjgx/).The annual mean temperature at the Bayanbulak meteorological station amounts to −4.1 ∘ C and the mean annual precipitation is 278 mm (Figure 1).Generally, precipitation falls as rain from May to September each year and as snow from October to April.
The observed streamflow data at the Dashankou hydrologic station are from the Xinjiang Tarim River Basin Management Bureau.The average flow at the Dashankou hydrological station amounts to approximately 120 m 3 /s (equivalent to 202 mm runoff/year), ranging from 15 m 3 /s to 973 m 3 /s.

Methodology
3.1.Regional Climate Model and RCP Scenarios.The outputs of a regional climate model (RegCM4.0)[13] forced by a global climate model (Beijing Climate Center Climate System Model, BCC CSM1.1) [14,15] at a horizontal resolution of 50 km are used as future climate data.Firstly, the RCM model was validated with the observational data set over China for the period from 1976 to 2005 and then it was used to predict the future climate change under the new emission scenarios of the Representative Concentration Pathways, RCP4.5 (lower emission scenario) and RCP8.5 (higher emission scenario).RCP4.5 is a stabilization scenario with the total radiative forcing rising until 2070 and then remaining at a stable centration of 4.5 W/m 2 .In contrast, RCP8.5 is a continuously rising radiative forcing pathway (at a target of 8.5 W/m 2 in 2100) with a further enhanced residual circulation and significant CH 4 increases [16,17].The RCM validation shows reasonable simulations of temperature and precipitation were obtained over China and, compared to the BCC CSM1.1 model, marked improvement of the RCM was achieved in reproducing present day precipitation and temperature (for more details refer to [18]).

Bias Correction Methods.
Five precipitation and three temperature correction methods were selected to bias-correct the raw RCM outputs.Precipitation correction methods include linear scaling, local intensity scaling, power transformation, distribution mapping, and quantile mapping.Temperature correction methods include linear scaling, variance scaling, and distribution mapping.They are combined into Advances in Meteorology 3 15 schemes to evaluate their performances in simulating streamflow.It turns out that the precipitation correction methods have more significant influence than the temperature correction methods on streamflow simulation, and the power transformation and quantile mapping perform best in terms of frequency based statistics.Thereafter, the quantile mapping method (for precipitation) and the distribution mapping (for temperature) are selected to correct the raw RCM outputs for the future period (for more details see [19]).

Hydrologic Model and Uncertainty Analysis Method.
SWAT has been extensively used for the comprehensive modeling of the impact of management practices and climate change.SWAT simulates the hydrologic and sedimentary processes, plant growth, river routing, and in-stream water quality process, among which the surface runoff is calculated from daily rainfall and snowmelt with a modified Soil Conservation Service (SCS) curve number method [20]; water routing is simulated using variable storage or the Muskingum river routing method [11].
The SWAT model input includes the digital elevation model (DEM), soil textural and physicochemical properties, and land use data.The meteorological variables including daily precipitation, max./min.temperature, relative humidity, solar radiation, and wind speed were used to force the SWAT model.SWAT uses elevation bands to represent the topographic effects on precipitation and temperature.Within each elevation band, the precipitation and temperature are estimated based on their lapse rates.For more details, refer to the SWAT manual (http://www.brc.tamus.edu/).
The SWAT model (forced by the observed meteorological data) was calibrated against the observed streamflow.The calibration period is from 1986 to 1989 and the validation period from 1990 to 2005 [21].The calibrated optimal parameters are then kept fixed in the following simulations.The evaluation indices for the hydrological model include NS, PBIAS, and the determination coefficient where  obs  and  sim  are the th observed and simulated flows,  mean is the mean of the observed data, and  is the number of observations.Normally NS > 0.50, |PBIAS| < 25%, and  2 > 0.6 are taken as the criteria for satisfactory modeling of the river discharge and the model performance can be evaluated as excellent if NS > 0.75 and |PBIAS| < 10% [22].GLUE (generalized likelihood uncertainty estimation) [23] is an uncertainty analysis technique, in which the parameter uncertainty accounts for all sources of uncertainty, such as input uncertainty, structure uncertainty, parameter uncertainty, and response uncertainty [24].In GLUE, the parameter uncertainty is described as a set of discrete "behavioral" parameter sets with corresponding "likelihood weights." The procedure of a GLUE analysis consists of three steps.Firstly, after the definition of the "generalized likelihood measure, " (), a large number of parameter sets are randomly sampled from the prior distribution and each parameter set is assessed as either "behavioral" or "nonbehavioral" by comparing its value of () to the threshold value.Secondly, each behavioral parameter set is given a "likelihood weight" and we gave them equal weights in this study.Finally, prediction uncertainty is represented by 5% and 95% quantiles of the cumulative distribution of the behavioral parameter sets.
Two indices are used to quantify the quality of the uncertainty performance.Those indices are the percentage of measurements bracketed by the 95% prediction uncertainty band (-factor) and width of band (-factor, calculated by the average width of the band divided by the standard deviation of the corresponding measured variable).

SCC Data Description and Analysis
Procedures.In the following section, temperature and precipitation are denoted as  and  and the absolute and relative changes are represented by Δ and .For example, Δ refers to an absolute temperature change and  a relative precipitation change.The hydrological processes analyzed in this study include streamflow, surface runoff, subsurface runoff, and evapotranspiration, which are denoted as ,   ,   , and ET, and their relative changes are described as ,   ,   , and ET, respectively.
By investigating the transient evolution of climate change in the corrected RCM outputs on decadal scales, five periods (each spanning 20 years) are defined: 1986∼2005 (control period) and 2020∼2039, 2040∼2059, 2060∼2079, and 2080∼ 2099.Due to the intra-annual characteristics of the hydrometeorology in the Kaidu River Basin (Figure 1), the wet season (from April to September) and dry season (from October to March next year) are defined based on the intra-annual distribution of  and ; for example,  and  in the wet season account for 88% and 73% of their annual amounts.The climatic and hydrological changes are classified into three categories, that is, a significant change, small change, and insignificant change, to clearly demonstrate the changing magnitude according to the values of relative change for precipitation and hydrological components and absolute change of temperature.These categories are presented in Table 1.2) indicate excellent performances for both calibration and validation periods, with "NS"s and " 2 "s over 0.80, which is highly acceptable according to Moriasi et al. [22].Concerning the monthly streamflow, the "NS" is 0.89 during 1986∼2005 and it indicates that the SWAT model captured the natural monthly streamflow variability adequately.The performances of bias-corrected RCM outputs (compared to observed meteorological data) are listed in Table 2.The "NS"s are −0.57(0.57) and 0.77 (0.95) for daily (monthly) precipitation and temperature for 1990∼2005, respectively.And the statistics of the streamflow simulated with the biascorrected RCM outputs shows acceptable results with "NS"s equal to 0.46 and 0.62 and PBIAS within 10% for daily and monthly streamflows.

RCM Projected Hydrometeorologic Changes
4.2.1.Changes in Temperature and Precipitation.Temperature is highly likely to increase in the future, with a basin warming of 1.0∼2.2∘ C and 1.6∼4.6 ∘ C under RCP4.5 and RCP8.5 in the 21st century (Table 3).Temperature increases continuously under both scenarios but the magnitude is larger under RCP8.5 (Figure 3).
Precipitation shows an overall increasing trend in the 21st century with an annual increase of 2%∼16% and 7%∼ 24% under RCP4.5 and RCP8.5, which confirms the previous arguments of Sorg et al. [1].However, precipitation change varies substantially among seasons (Figure 3).Normally, a small increase in the wet season (−2%∼16%) and a significant increase during dry season (18%∼78%) are projected.Note that the relative increase (not the absolute increment) of precipitation for the dry season is much bigger than for the wet season, which is in line with the climate changes in other regions, for example, the semiarid Colorado River Basin [25] and the wet Ganges-Brahmaputra-Meghna basin [26].

Changes in the Hydrological
Cycle.The changes in precipitation and temperature cause changes in potential streamflow.The average annual streamflow rises by −1%∼ 18% and 4∼20% under the RCP4.5 and RCP8.5 in the 21st century, based on the average annual streamflow of 194 mm for the control period (1986∼2005) (Table 3).Note that the streamflow stopped increasing in 2080∼2099 (end of 21st century) under RCP8.5 despite the rise in precipitation, which confirms the finding of Sorg et al. [1] and may aggravate water scarcity in this region.Figure 3 also shows the projected changes in surface runoff (  ), subsurface runoff (  ), and evapotranspiration (ET) under RCP4.5 and RCP8.5.Overall, changes of hydrologic components are bigger for RCP8.5 than for RCP4.5.The annual change of   is insignificant (<5%) but with obvious seasonal variability; for example, changes of   range from −22% to 2% for the wet season and in 4%∼78% for the dry ET (mm) ET (%) season under RCP8.5.The annual   changes by −0.7%∼17% and 4%∼18% for RCP4.5 and RCP8.5, which is consistent with the changes of .ET increases continuously in the 21st century with average increases of 2%∼10% and 7%∼24% under RCP4.5 and RCP8.5.(2) the streamflow is also influenced by temperature dominated snowmelt (snowfall accounts for about 17% of watershed precipitation), which reduces the dependence of streamflow on precipitation and therefore results in a low streamflow elasticity [27].The response of  to Δ depends on the magnitude of Δ. decreases slightly when 0 < Δ ≤ 2.0 ∘ C while it decreases dramatically when Δ > 2.0 ∘ C for both wet and dry seasons.For example, when Δ = 2.0 ∘ C, a 40% precipitation increase results in an average value of  being 240 mm (23% increase compared to the average streamflow of 194 mm) but when Δ = 4.0 ∘ C, the same precipitation increase only generates an average  of 180 mm (about 7% decrease) (Figure 4).

Response of Hydrological
The responses of   ,   , and ET to climate change are also demonstrated in Figure 4.For   , the responses of   to Δ are quite different for the wet and dry seasons: the higher Δ, the lower   for the wet season but the higher   for the dry season.Since   in the dry season only accounts for 13% of the annual   , the response of annual   is consistent with that of the wet season.For   , the responses of   to Δ and  are similar to the responses of  due to the dominant role of groundwater recharge in water yield in the Kaidu River Basin.For ET, it is mainly influenced by Δ with temperature sensitivity (ET/Δ) being 7.3%/ ∘ C. To verify this result, we firstly investigated basin-scale energy and water budget using the Budyko method [28,29].It is shown that ET is mainly energy limited rather than water limited (average ET/ = 0.67 and PET/ = 0.88).Secondly, the high determinate coefficient  2 = 0.75 (significant level is smaller than 0.01) between the mean annual  and ET also indicates that ET has a strong correlation with .This is consistent with previous studies, which have shown that a significant variation in ET is expected to follow changes in air temperature [30,31].
In addition, simulations with RCM outputs are shown in Figure 4 to analyze the differences between simulations of these two data sets.Two typical periods of RCM simulations are selected, that is, 2020∼2039 under RCP4.5 and 2080∼2099 under RCP8.5, to represent mild and intense climate change scenarios (shown as blue and red stars in Figure 4).It is indicated that the simulations of hydrological components with RCM outputs for 2020∼2039 under RCP4.5 (mild climate change) are close to these of the nearby contour lines (simulations with SCC), which suggests that similar results of ,   ,   , and ET are obtained for RCM outputs and for SCC under mild climate change scenarios.However, for 2080∼2099 under RCP8.5 with  = 24% and Δ = 4.6 ∘ C, the simulated values of ,   ,   , and ET deviate from the simulations of SCC.There are two possible reasons: (1) changes of other meteorological inputs, that is, solar radiation, wind speed, and humidity, are slightly smaller for ET (mm) The simulations with RCM outputs for 2020∼2039 under RCP4.5 and for 2080∼2099 under RCP8.5 (their corresponding meteorological changes are listed in Table 3) are indicated using blue and red stars with labels.2020∼2039 under RCP4.5 than those for 2080∼2099 under RCP8.5 (−0.8%, 2.6%, and 0.9% compared to −2.2%, 4.1%, and 1.4%); (2) for 2080∼2099 under RCP8.5, precipitation increases by 24% with great seasonal variation, which may alter the hydrological regime; for example, precipitation increases by 139% for March, April, and May while it decreases by −0.1% for June, July, and August.Since changes of solar radiation, wind speed, and humidity are within ±5%, the second reason, that is, the shift of the precipitation temporal distribution, contributes a lot to the deviation of simulations with RCM outputs from simulations with SCC.
Furthermore, the exceedance probability curves of the annual runoff in response to climate change are demonstrated in Figure 5.The exceedance probability curves are almost parallel when Δ ranges in 0∼6 ∘ C.However, the responses of  to  are not the same for each exceedance probability: high sensitivity of  with probabilities less than 0.1 and low sensitivity of  with probabilities larger than 0.9.A comparison of the simulations with RCM outputs (four future periods under RCP4.5 and RCP8.5) and with SCC (the same changes in  or  with the corresponding RCM outputs) indicates that differences between simulation with RCM outputs and SCC are becoming greater as climate change gets more intense; for example, the simulation with RCM for 2080∼2099 under RCP8.5 overestimates the corresponding simulations with SCC (Figure 5(h)), which collaborates the conclusion that under intense climate change scenarios the simulated hydrology with RCM deviates from that simulated with SCC.
The contributions of hydrologic components to water yield are displayed by the De Finetti diagram in Figure 6.For the control period, the averages of   ,   , and ET are 0.22, 0.28, and 0.50.For SCC, as Δ increases from 0 to 6 ∘ C, the contribution of ET increases rapidly from 0.49 to 0.73 and the contributions of   and   decrease from 0.22 to 0.11 and from 0.29 to 0.16.Δ has a more significant influence on the proportion change than .As  changes from −20% to 60%, ET decreases from 0.71 to 0.58 and   and   increase from 0.13 to 0.15 and from 0.16 to 0.27.For simulations with RCM outputs, proportions of hydrological components do not change significantly under RCP4.5, while the proportion of ET shows a significant increase under RCP8.5.

Sources of Uncertainty and Other Considerations.
There are uncertainties in estimating climate change impact on hydrology.As indicated by previous studies [32], the sources  3).
the land surface under the future climate changes.Effects of land cover change on streamflows and other components of the hydrological cycle are not considered.Though uncertainty in hydrologic modeling was quantified with the GLUE method, it only accounts for part of the total uncertainty in climate change impact studies [34].Uncertainties associated with the climate model and downscaling were not considered here although two emission scenarios were included.Any uncertainty associated with them may cause the results to deviate from reality.However, we are dedicated to pursuing a thorough investigation of the response of the hydrological cycle to future climate change for this region and we believe this study is an important first step in achieving this goal.

Conclusions
This study assessed the implications of climate change on hydrology in a typical watershed in the Tianshan Mountains with two sets of climatic data, that is, RCM outputs and SCC, loosely coupled to a hydrological model (SWAT).Major conclusions can be summarized as follows: (1) The hydrological model shows excellent performance with "NS"s over 0.8 for the daily streamflow for both calibration and validation periods.And the selected bias correction methods were effective in downscaling RCM outputs, with "NS"s being 0.57 and 0.95 regarding monthly precipitation and temperature.
( It is valuable to quantify the future responses of hydrology to climate change in the Tianshan Mountains.This study will provide useful information for water resource management and will serve as a basis for further climate change impact studies.

Figure 1 :
Figure 1: Location (top left), topography, and hydrometeorologic stations (bottom) of the Kaidu River Basin and the daily average temperature () and precipitation () at the Bayanbulak station and streamflow () at the Dashankou station (top right).

Table 1 : 4 . 1 .
Classification of magnitude for climatic and hydrological changes. and Δ represent relative change and absolute change.Validation of the Hydrological Model and the Bias Correction Methods.Performance of the hydrological model forced by observed meteorological data and the 95% prediction uncertainty bands are shown in Figure 2. The simulated streamflow agrees quite well with the observation for both calibration period (1986∼1989) and validation period (1990∼  2002).For the uncertainty analysis, NS is used as () and 0.70 as threshold value with 10,000 initial parameter sets; 288 sets were selected as behavioral points.The results show that most of the observations are bracketed by the 95% prediction uncertainty band (-factor being 87% and 80% for calibration and validation periods and -factor being 1.18 and 1.19, resp.).The lower -factor for the validation period can be partly attributed to operation of hydropower station since 1991 (Figure2) which leads to great fluctuation in winter streamflow.Statistics of model efficiency (Table

Figure 3 :
Figure 3: Summary of future climate inputs ( and ) and simulated hydrologic components (,   ,   , and ET) under RCP4.5 and RCP8.5, compared to their values in the control period (1986∼2005).All these hydrometeorologic factors are presented in terms of wet season, dry season, and annual values.

Figure 4 :
Figure 4: Response surfaces of streamflow (), surface runoff (  ), subsurface runoff (  ), and evapotranspiration (ET) to climate change.The simulations with RCM outputs for 2020∼2039 under RCP4.5 and for 2080∼2099 under RCP8.5 (their corresponding meteorological changes are listed in Table3) are indicated using blue and red stars with labels.

Table 2 :
Statistics of bias-corrected RCM outputs and SWAT simulated streamflows forced by the observed climate variables and biascorrected RCM outputs.
[19]as correction methods used are quantile mapping for precipitation and distribution mapping for temperature[19].b" Daily" or "monthly" in the brackets means the time step used to calculate the statistics.
)  increases almost linearly with  while the response of  to  depends on the magnitude of Δ and  decreases significantly when Δ is greater than 2 ∘ C. (4) Similar responses of ,   ,   , and ET to  and  are obtained for the RCM outputs and for SCC under mild climate change scenarios.However, for intense climate change scenarios, simulations of ,   ,   , and ET with RCM outputs (e.g., for 2080∼2099 under RCP8.5)deviate from simulations with SCC.(5) Δ has more significant influence on the proportion change of each hydrologic component than  does.As Δ increases from 0 to 6 ∘ C, the contribution of ET increases rapidly from 0.49 to 0.73 and   and   decrease by 0.11 and 0.13.As  changes from −20% to 60%, ET,   , and   change by −0.13, 0.02, and 0.11 as a result.