Correction: Xu, M., et al. Modeling Glacier Mass Balance and Runo ﬀ in the Koxkar River Basin on the South Slope of the Tianshan Mountains, China, from 1959 to 2009. Water 2017, 9 , 100 Correction: Xu, M. et al. Modeling Glacier Mass Balance and Runoff in the Koxkar River Basin on the South Slope of the Tianshan Mountains, China, from 1959 to 2009

The authors wish to make the following corrections to this paper [...]


Introduction
Glaciers concentrated in mountain regions are important water resources in arid and semi-arid areas in west China [1][2][3].Rising temperatures have resulted in the loss of more than half of the glaciers in Asia since 1900 [4].In the northwest arid region of China, the importance of meltwater is particularly prominent.Meltwater accounts for 22% of the runoff from mountains.Melted water from glaciers from the Tianshan Mountains supplied 38.5% of the total runoff in the Tarim River basin [5].Since the 1980s, runoff has significantly increased, flowing out from the Tianshan Mountains.The glacial retreat has major impact on water resources in arid regions [6].
Hydrology in the arid region has become an important topic [7].The hydrological model can simulate hydrological processes; it has practical significance and scientific value for exploring new ideas about management of water resources in the river basin [8].However, due to the extremely high altitude and difficult environment of the Tianshan Mountains, hydro-meteorological data are rare in Water 2017, 9, 100 2 of 19 this areas, especially in the glaciated watersheds.Additionally, there have been few applications of hydrological models to the south slope of Mount Tumur, Tianshan Mountains.There are more than 2.52 million people and 1.33 million acres of irrigated agricultural land in the lower reaches of the river basin which depend on the runoff coming from high mountains.Therefore, it is urgent to investigate water resources from glaciers in this area.It will also help understand the changes in water resources under the climate change.
The HBV (Hydrologiska Byrans Vattenavdelning) model was developed by the Swedish Meteorological and Hydrological Institute (SMHI) for cold regions in 1973.The HBV model, with its strong practical aspects and improved simulation results, has been applied successfully to forecast runoff for more than 200 river basins in Sweden and has been successfully applied in more than 40 countries of the world to study hydrological processes under climate change, changes of permafrost, and changes of land use [9][10][11][12][13][14].The HBV model is a semi-distributed conceptual model; it produces better results of simulation when less data input is required [15].It has also been successfully used in different watersheds in the cold regions of west China.Kang et al. [16] established a conceptual hydrological model for upstream of Hei river in Qilian Mountains of west China based on the principle of the HBV model.The model is accordance with the characteristics of formation of mountain runoff in study area, where the basin was divided into two basic landscape zones, alpine tundra-snow and mountain vegetation [16].In addition, the model has been applied to glaciated watersheds in the Tibetan Plateau, where it also performed well; the HBV light model includes a glacier module which is more applicable in mountain regions where glaciers are distributed [17].
In this study, the HBV light model and the principle of water balance were employed to simulate the change in depth of runoff, evaporation, and glacier mass balance (GMB) in the Koxkar River Basin (KRB) between 1959 and 2009.The work aims to reconstruct climate change in the past 50 years in the region in order to understand changes in glacier mass balance and runoff in response to climate change.This work may provide knowledge for water resource management in the arid regions of Central Asia.

Study Area
The KRB is located in the south of Mount Tumur of Western Tianshan Mountains; it is one of the source regions of the Akesu River (Figure 1).The area of KRB is 116.5 km 2 , the area of glacier is 89.6 km 2 (76.9%), and its elevation ranges from 2960 to 6342 m.a.s.l.The length of Koxkar glacier is 25.1 km, the ice volume is 15.79 km 3 .Annual precipitation is about 533 mm in the KRB, which is supplied mainly by moisture from the Atlantic and Arctic.The precipitation which occurs from June to August accounts for 50%, 70% of precipitation occurs from May to September, and about 30% occurs during the cold season.Annual precipitation near the snow line is between 750 and 850 mm [18][19][20].There is one automatic weather station (AWS) and one hydrological station located in the terminal of glacier (Base Camp) (2974 m.a.s.l.) (Figure 1), and four other AWSs located in the glacier with elevations of 3200, 3400, 3700, and 4200 m.a.s.l., respectively (Figure 1).

Geographical Data
The HBV light model requires geographical data including a digital elevation model (DEM), aspect, and land cover.The 30-m-resolution Global DEM (http://datamirror.csdb.cn)and a digital vector map of modern glaciers were obtained from the Chinese Glacier Inventory [20].The basin was divided by altitude zones, aspect zones, glacier zones and nonglacier zones (Figure 2).The land cover data were taken from the Science Data Center for Cold and Arid Regions (http://westdc.westgis.ac.cn/).The types of land cover include bare land, shrub, grassland, snow cover, glacier, and water (Table 1).

Geographical Data
The HBV light model requires geographical data including a digital elevation model (DEM), aspect, and land cover.The 30-m-resolution Global DEM (http://datamirror.csdb.cn)and a digital vector map of modern glaciers were obtained from the Chinese Glacier Inventory [20].The basin was divided by altitude zones, aspect zones, glacier zones and nonglacier zones (Figure 2).The land cover data were taken from the Science Data Center for Cold and Arid Regions (http://westdc.westgis.ac.cn/).The types of land cover include bare land, shrub, grassland, snow cover, glacier, and water (

Observations and Meteorological Data
The HBV light model requires daily air temperature, daily precipitation, gradients of temperature and precipitation data for input.Five AWSs have been established in KRB since 2007, they have operated until present (Figure 1).Precipitation was monitored by both manual rain gauge and automatic snow/rain gauge (T-200B weighting rain gauge).This study used air temperature and

Observations and Meteorological Data
The HBV light model requires daily air temperature, daily precipitation, gradients of temperature and precipitation data for input.Five AWSs have been established in KRB since 2007, they have operated until present (Figure 1).Precipitation was monitored by both manual rain gauge and automatic snow/rain gauge (T-200B weighting rain gauge).This study used air temperature and precipitation from the AWS at the Base Camp (Camp-AWS) as input data, they are of higher quality.Based on the altitudes of the five AWSs, we calculated the gradients of temperature and precipitation: Water 2017, 9, 100 5 of 19 they are 0.66 • C (100 m) −1 and 4.6% (100 m) −1 , respectively.Since meteorological observations started in 2007 in the KRB, the long-term air temperature and precipitation  data were reconstructed using four national meteorological stations near to KRB (NMS´s) (Figure 1, Table 2) (for the method see Section 4).

Hydrological and Glacial Data
The hydrological station was established at the base camp (BC) in 2007 (Figure 1).Water stage was measured by an automatic HOBO Water Level Logger (Onset, Bourne, MA, USA) with an interval of 30 min; velocity of flow was measured manually by a hydrodynamometer.Through several measurements at different water stages, the relationship between water stage and discharge was established.Then, the discharge was calculated at the BC from 2007 to 2011.

Hydrologiska Byrans Vattenavdelning (HBV) Model
The HBV model is a semi-distributed model, which means that a catchment can be separated into different elevation and vegetation zones as well as into different subcatchments [21,22].The HBV light is a conceptual mode which is based on the HBV model [23].It was coupled with a more detailed snow and glacier melt subroutine, employing the degree-day approach [24].The subroutine has a physical basis.For this study, the HBV light model, which takes into account various aspect classes in each elevation belt and couples Monte Carlo and genetic algorithms to calibrate the parameters automatically [25,26], was applied to model runoff depth, glacier runoff, and GMB in the KRB.The model is composed mainly of the input data, which were processed as semi-distributed by four modules: the glacier ablation module, soil module, response module, and flow concentration module (Figure 3) [22].
precipitation from the AWS at the Base Camp (Camp-AWS) as input data, they are of higher quality.Based on the altitudes of the five AWSs, we calculated the gradients of temperature and precipitation: they are 0.66 °C (100 m) −1 and 4.6% (100 m) −1 , respectively.Since meteorological observations started in 2007 in the KRB, the long-term air temperature and precipitation  data were reconstructed using four national meteorological stations near to KRB (NMS´s) (Figure 1, Table 2) (for the method see Section 4).

Hydrological and Glacial Data
The hydrological station was established at the base camp (BC) in 2007 (Figure 1).Water stage was measured by an automatic HOBO Water Level Logger (Onset, Bourne, MA, USA) with an interval of 30 min; velocity of flow was measured manually by a hydrodynamometer.Through several measurements at different water stages, the relationship between water stage and discharge was established.Then, the discharge was calculated at the BC from 2007 to 2011.

Hydrologiska Byrans Vattenavdelning (HBV) Model
The HBV model is a semi-distributed model, which means that a catchment can be separated into different elevation and vegetation zones as well as into different subcatchments [21,22].The HBV light is a conceptual mode which is based on the HBV model [23].It was coupled with a more detailed snow and glacier melt subroutine, employing the degree-day approach [24].The subroutine has a physical basis.For this study, the HBV light model, which takes into account various aspect classes in each elevation belt and couples Monte Carlo and genetic algorithms to calibrate the parameters automatically [25,26], was applied to model runoff depth, glacier runoff, and GMB in the KRB.The model is composed mainly of the input data, which were processed as semi-distributed by four modules: the glacier ablation module, soil module, response module, and flow concentration module (Figure 3) [22].The main formulas of each module are described as follows.The model simulates daily discharge using daily precipitation and temperature as input data.Precipitation is simulated to be either snow or rain depending on whether the temperature is above or below a threshold temperature, TT ( • C).All precipitation simulated to be snow, i.e., falling when the temperature is below TT, is multiplied by a snowfall correction factor, SFCF.In the snow and glacier routine, the amount of melting water is calculated by the degree-day method [27], which uses positive degree-day (while daily mean temperature is above 0 • C) multiplied by a factor.This model also considers the influence of different aspects and different melt intensity of glacier and snow [25] as shown in Equation (2).Meltwater and precipitation are retained within the snowpack until the amount exceeds a certain fraction, CWH, of the water equivalent of the snow.Refreezing of meltwater is also considered in this model [28][29][30].Liquid water within the snowpack refreezes according to Equation (3).Rainfall and snowmelt (P) are divided into water filling the soil box and groundwater recharge depending on the relation between water content of the soil box (SM, in mm) and its largest value (FC, in mm) as shown in Equation ( 4).Actual evaporation from the soil box equals the potential evaporation if SM/FC is above LP (the threshold of Soil water content and Field water holding capacity), whereas a linear reduction is used when SM/FC is below LP, as shown in Equation ( 5) [23].Groundwater recharge is added to the upper groundwater box (SUZ, in mm).PERC (mm•day −1 ) defines the maximum percolation rate from the upper to the lower groundwater box (SLZ, in mm).Runoff from the groundwater boxes is computed as the sum of two or three linear outflow equations depending on whether or not SUZ is above a threshold value, UZL (mm), as shown in Equation ( 6).This runoff is finally transformed by a triangular weighting function defined by the parameter MAXBAS, as shown in Equation ( 7) to give the simulated runoff (mm•day −1 ).If different elevation zones are used, the changes in precipitation and temperature with elevation are calculated using the two parameters P CALT (%/100 m) and T CALT ( • C/100 m), as shown in Equations ( 8) and ( 9).The long-term mean of the potential evaporation, E pot , M, for a certain day of the year can be corrected to its value at day t, E pot (t), by using the deviations of the temperature, T(t), from its long-term mean, TM, and a correction factor, CET ( • C −1 ), as shown in Equation (10) [31,32].

Reconstruction of Temperature and Precipitation Data
Understanding the climate, discharge, and GMB in past decades is important for predicting future trends.In this work, we used the HBV hydrological model to simulate runoff and GMB from 1959 to 2009.Long-term observed data are unavailable for the basin, so it was necessary to reconstruct the historical runoff and GMB.The daily air temperature and precipitation in the KRB were reconstructed from 1959 to 2009 based on in-situ meteorological observations in the KRB between 2007 and 2009 and four NWSs around the KRB for the period 1959-2009.

Reconstruction of Temperature
There are many methods to interpolate temperature, e.g., inverse distance weight, trend surface, kriging interpolation, and cokring interpolation; these methods are usually used to interpolate the spatial distribution of temperature.However, they require a high number of data points for calculation.In this study, we had only four NMSs for interpolation (Table 2).Thus, the previous temperature data were interpolated by the statistical method of multiple regression analysis, which depended only on the temperature recorded at the four stations.Two equations were formulated because the four stations had records of various lengths.
T 0 is the daily air temperature at Camp-AWS, T 1 , T 2 , T 3 , and T 4 are the daily air temperatures of Akesu, Baichen, Aheqi, and Zhaosu, respectively.The result calculated by Equation (13) (R 2 = 0.942) was better than that calculated by Equation ( 14) (R 2 = 0.918), so we used Equation (13) to reconstruct the daily temperature.The reconstructed daily temperatures fit well to the observation ones from 2007 to 2009 (R 2 = 0.91; p < 0.01; Figure 4).The temperature increased 0.32 • C (10a −1 ; p < 0.01) in the past semi-century in the KRB (Figure 5), which is a little higher than the average change in the southern slope of the Tianshan Mountains.(0.30 • C, 10a −1 ) [34].

Reconstruction of Precipitation
The precipitation, which has large temporal variation, was influenced by the terrain in the Tianshan Mountains.The reconstruction of precipitation was conducted as follows.We used stepwise multiple regression analysis to eliminate stations that recorded precipitation farther from the Camp-AWS precipitation in a step-by-step manner.Stations were eliminated in the order of Akesu, Baichen, Aheqi, and Zhaosu, respectively.The correlation between Camp-AWS and Akesu was not significant due to Akesu being located in the lowland city.Therefore, we used data from other three stations located in the mountain regions to interpolate the Camp-AWS precipitation.
The gradient of precipitation and inverse distance weighting (IDW) methods were combined to reconstruct Camp-AWS precipitation from 1959 to 2009.The gradients of precipitation were 2.05 mm/100 m, 3.54 mm/100 m, and 1.6 mm/100 m between Camp-AWS and Baichen, Aheqi, and Zhaosu, respectively, from 2007 to 2009.Based on the three precipitation gradients, we calculated the precipitations at the same elevation as Camp-AWS from 1959 to 2009.Then, we used IDW to calculate the precipitation at Camp-AWS from 1959 to 2009, as shown in Equation (15).

Reconstruction of Precipitation
The precipitation, which has large temporal variation, was influenced by the terrain in the Tianshan Mountains.The reconstruction of precipitation was conducted as follows.We used stepwise multiple regression analysis to eliminate stations that recorded precipitation farther from the Camp-AWS precipitation in a step-by-step manner.Stations were eliminated in the order of Akesu, Baichen, Aheqi, and Zhaosu, respectively.The correlation between Camp-AWS and Akesu was not significant due to Akesu being located in the lowland city.Therefore, we used data from other three stations located in the mountain regions to interpolate the Camp-AWS precipitation.
The gradient of precipitation and inverse distance weighting (IDW) methods were combined to reconstruct Camp-AWS precipitation from 1959 to 2009.The gradients of precipitation were 2.05 mm/100 m, 3.54 mm/100 m, and 1.6 mm/100 m between Camp-AWS and Baichen, Aheqi, and Zhaosu, respectively, from 2007 to 2009.Based on the three precipitation gradients, we calculated the precipitations at the same elevation as Camp-AWS from 1959 to 2009.Then, we used IDW to calculate the precipitation at Camp-AWS from 1959 to 2009, as shown in Equation (15).

Reconstruction of Precipitation
The precipitation, which has large temporal variation, was influenced by the terrain in the Tianshan Mountains.The reconstruction of precipitation was conducted as follows.We used stepwise multiple regression analysis to eliminate stations that recorded precipitation farther from the Camp-AWS precipitation in a step-by-step manner.Stations were eliminated in the order of Akesu, Baichen, Aheqi, and Zhaosu, respectively.The correlation between Camp-AWS and Akesu was not significant due to Akesu being located in the lowland city.Therefore, we used data from other three stations located in the mountain regions to interpolate the Camp-AWS precipitation.
The gradient of precipitation and inverse distance weighting (IDW) methods were combined to reconstruct Camp-AWS precipitation from 1959 to 2009.The gradients of precipitation were 2.05 mm/100 m, 3.54 mm/100 m, and 1.6 mm/100 m between Camp-AWS and Baichen, Aheqi, and Zhaosu, respectively, from 2007 to 2009.Based on the three precipitation gradients, we calculated the precipitations at the same elevation as Camp-AWS from 1959 to 2009.Then, we used IDW to calculate the precipitation at Camp-AWS from 1959 to 2009, as shown in Equation (15).

Model Calibration and Validation
Except for meteorological data, there are 13 parameters in the HBV model that need to be calibrated (Table 3).To obtain a proper set of parameters, we gave every parameter a reasonable range according to the literature.A set of parameters representing the characteristics of the KRB was then automatically obtained by observation and by employing the Monte Carlo method [35][36][37]  In Equation ( 15), P0 is the daily precipitation of Camp-AWS and P1, P2, and P3 are the daily precipitations of the NMSs at Baichen, Aheqi, and Zhaosu, respectively.P0 annual, P1 annual, P2 annual, and P3 annual are the average annual precipitation of Camp-AWS, Baichen, Aheqi, and Zhaosu, respectively, between 2007 and 2009.S1, S2, and S3 are the distances between Camp-AWS and Baichen, Aheqi, and Zhaosu, respectively.
Comparison of reconstructed with observed monthly precipitation is shown in Figure 6 (R 2 reached 0.91; p < 0.01).Because precipitation is very important for the simulation of runoff and GMB, we needed to obtain a more vigorous validation for the reconstructed result.Thus, the coefficient of efficiency (Reff); Equation ( 11) was used to evaluate the reconstructed daily result.Reff was 0.88 from 2007 to 2009, which is fairly satisfactory considering the great varying of precipitation in the mountainous regions (Figure 6).By applying the above method, the annual precipitation at Camp-AWS from 1959 to 2009 was reconstructed (Figure 7).The precipitation increased by 4.8 mm•a −1 (p < 0.05), consistent with other studies in the regions [34].

Model Calibration and Validation
Except for meteorological data, there are 13 parameters in the HBV model that need to be calibrated (Table 3).To obtain a proper set of parameters, we gave every parameter a reasonable range according to the literature.A set of parameters representing the characteristics of the KRB was then automatically obtained by observation and by employing the Monte Carlo method [35][36][37].We

Model Calibration and Validation
Except for meteorological data, there are 13 parameters in the HBV model that need to be calibrated (Table 3).To obtain a proper set of parameters, we gave every parameter a reasonable range according to the literature.A set of parameters representing the characteristics of the KRB was then automatically obtained by observation and by employing the Monte Carlo method [35][36][37].We selected hydrological data from 2007/2008 and 2008/2009 for calibration.In addition, data for 2009/2010 and 2010/2011 were used for validation.For the period of calibration, R eff and R 2 are 0.77 and 0.76 for 2007/2008, and 0.87 and 0.84 for 2008/2009, respectively.For the period of validation, R eff and R 2 are 0.89 and 0.89 for 2009/2010 and 0.80 and 0.81 for 2010/2011, respectively (Figure 8).The results are acceptable on the daily scale, especially considering the uncertainty in both observation and simulation data in such a mountainous area.The modeling result reflects the actual runoff well (Figure 9).All R eff values of the daily runoff are greater than 0.76.R eff is even better (0.95) for the monthly result (Figure 10), and there is an only 6% difference between the observed and simulated annual runoff.Considering the dominant role of glacier runoff in the KRB, the reason for the above differences may be that increasing the temporal resolution led to decreases in the accuracy of the degree-day model [38].
Water 2017, 9, 100 10 of 19   8).The results are acceptable on the daily scale, especially considering the uncertainty in both observation and simulation data in such a mountainous area.The modeling result reflects the actual runoff well (Figure 9).All Reff values of the daily runoff are greater than 0.76.Reff is even better (0.95) for the monthly result (Figure 10), and there is an only 6% difference between the observed and simulated annual runoff.Considering the dominant role of glacier runoff in the KRB, the reason for the above differences may be that increasing the temporal resolution led to decreases in the accuracy of the degree-day model [38].

Simulation of Runoff and GMB for 1959-2009
We assumed that the set of parameters did not change from 1959 to 2009, and the depth of runoff and evaporation were simulated (Figure 6a,b) using the reconstructed data.The runoff depth increased 5.5 mm•a −1 (p < 0.01), and it increased more significantly after the 1990s than before the 1990s.The simulated result is consistent with the results for other watersheds in which runoff is observed in the Tianshan Mountains.[39].Evaporation increased by 3.5 mm•a −1 in the KRB.On the basis of the water balance principle shown in Equation ( 16), the change in water volume in the basin (∆S) equals the total precipitation (P) minus the runoff (R), soil evaporation and sublimation (E s ), and evaporation at the glacier surface (E g ).
∆S is the sum of the soil moisture change (∆S S ) and the GMB, as shown in Equation ( 17).The long-term GMB can be calculated according to Equation (18).Because there is no permafrost in the KRB, it can be assumed that ∆S S did not change.The contribution of the change in ice content in permafrost to the water balance is less than 1% on a yearly scale [40].Precipitation was calculated by considering the precipitation gradient.E S was simulated, and E g was regarded as constant at 137.2 mm (Figure 11c) [41].Because there is no long-term observation of GMB in the KRB, we chose the two nearest glaciers (No. 1, in the headquarter of the Urumuqi river, and Tuyuksu) in the Tianshan Mountains, which have been observed in situ, for comparison from 1959 to 2009.As can be seen from Table 4, the value of GMB in Koxkar Glacier is between the values for the No. 1 and Tuyuksu glaciers, thus the simulated result is reasonable and credible.Glacial ablation in the KRB increased at a rate of 4.2 mm•a −1 (p < 0.01) from 1959 to 2009, and the average elevation of glacier decreased by 18.53 m (16.68 m w.Equation), or at a rate of 370 mm•a −1 (Figure 11c).

Simulation of Runoff and GMB for 1959-2009
We assumed that the set of parameters did not change from 1959 to 2009, and the depth of runoff and evaporation were simulated (Figure 6a,b) using the reconstructed data.The runoff depth increased 5.5 mm•a −1 (p < 0.01), and it increased more significantly after the 1990s than before the 1990s.The simulated result is consistent with the results for other watersheds in which runoff is observed in the Tianshan Mountains.[39].Evaporation increased by 3.5 mm•a −1 in the KRB.On the basis of the water balance principle shown in Equation ( 16), the change in water volume in the basin (ΔS) equals the total precipitation (P) minus the runoff (R), soil evaporation and sublimation (Es), and evaporation at the glacier surface (Eg).
ΔS is the sum of the soil moisture change (ΔSS) and the GMB, as shown in Equation ( 17).The long-term GMB can be calculated according to Equation (18).Because there is no permafrost in the KRB, it can be assumed that ΔSS did not change.The contribution of the change in ice content in permafrost to the water balance is less than 1% on a yearly scale [40].Precipitation was calculated by considering the precipitation gradient.ES was simulated, and Eg was regarded as constant at 137.2 mm (Figure 11c) [41].Because there is no long-term observation of GMB in the KRB, we chose the two nearest glaciers (No. 1, in the headquarter of the Urumuqi river, and Tuyuksu) in the Tianshan Mountains, which have been observed in situ, for comparison from 1959 to 2009.As can be seen from Table 4, the value of GMB in Koxkar Glacier is between the values for the No.

Sensitivity Analysis of Runoff
In this study, we used the HBV model to analyze the sensitivity of runoff to changes of temperature and precipitation in the KRB.Climate change in the area is manifested mainly as changes in temperature and precipitation, and it has had a significant impact on the hydrological processes of the watershed, which is partly covered by glacier.Since the required input of the HBV model is daily data and the time interval of predicted meteorological data is monthly, in order to quantitatively study the sensitivity of runoff to changes in temperatures and precipitation, we assumed that the set of parameters and the glacier did not change in the future.Thus, the temperature was increased by 1, 2, and 4 °C, the precipitation was increased by 10%, 20%, and 30%, and the runoff was simulated by the HBV model.
When the temperature increased by 1 °C to 4 °C (T + 1 °C, T + 2 °C, T + 4°C), the annual runoff increased by 27.4%, 56.2%, and 127%, respectively (Figure 12, Table 5).As can be seen from Figure 12, in the case of rising temperature, increases in the amount of runoff occur mainly from May to October.The month with the largest runoff is August, which shows that the runoff was supplied mainly by meltwater from the glacier.The distribution of internal annual runoff changes when temperature is higher (Figure 12).Significantly increased runoff normally appears in May, and significantly decreased runoff normally appears in October.However, when temperature is higher, the significantly increased and decreased runoff periods shift to earlier and later, respectively.

Sensitivity Analysis of Runoff
In this study, we used the HBV model to analyze the sensitivity of runoff to changes of temperature and precipitation in the KRB.Climate change in the area is manifested mainly as changes in temperature and precipitation, and it has had a significant impact on the hydrological processes of the watershed, which is partly covered by glacier.Since the required input of the HBV model is daily data and the time interval of predicted meteorological data is monthly, in order to quantitatively study the sensitivity of runoff to changes in temperatures and precipitation, we assumed that the set of parameters and the glacier did not change in the future.Thus, the temperature was increased by 1, 2, and 4 • C, the precipitation was increased by 10%, 20%, and 30%, and the runoff was simulated by the HBV model.
When the temperature increased by 1 , the annual runoff increased by 27.4%, 56.2%, and 127%, respectively (Figure 12, Table 5).As can be seen from Figure 12, in the case of rising temperature, increases in the amount of runoff occur mainly from May to October.The month with the largest runoff is August, which shows that the runoff was supplied mainly by meltwater from the glacier.The distribution of internal annual runoff changes when temperature is higher (Figure 12).Significantly increased runoff normally appears in May, and significantly decreased runoff normally appears in October.However, when temperature is higher, the significantly increased and decreased runoff periods shift to earlier and later, respectively.When precipitation increases by 10%, 20%, and 30% (P + 10%, P + 20%, P + 30%), monthly runoff increases by 6.1%, 11.3%, and 19.8%, respectively (Figure 13, Table 6).Increased amounts of runoff occur mainly in the ablation season (from May to October).However, when precipitation increases, the runoff does not increase very significantly.This is mainly because the runoff in the Tianshan Mountains is limited by temperature conditions.When only precipitation increases, runoff does not have a significant increase.Under the precipitation increases, increased runoff is concentrated mainly in June, July, and August.Within this period of time, radiation is sufficient and increased precipitation on the glacier leads to melt, which leads to increased runoff in the watershed.Increased precipitation causes a reduction in runoff from October to April (Figure 13).The increase of air temperature and precipitation has a great influence on the seasonal distribution of runoff.As the temperature increases, the increased degree of ablation decreases.The coefficient variations (Cv) of runoff increase at first and then decrease (Figure 14).However, the differences in CV are not obvious when temperature increases 1 °C and 2 °C, indicating that the fluctuation of runoff is not significant in this situation.When the temperature increases 4 °C, Cv decreases, and the fluctuation of the runoff decreases.With the increase of precipitation, Cv becomes greater, and the fluctuation of runoff increases, but fluctuation is lower than that caused by the changes of temperature.When precipitation increases by 10%, 20%, and 30% (P + 10%, P + 20%, P + 30%), monthly runoff increases by 6.1%, 11.3%, and 19.8%, respectively (Figure 13, Table 6).Increased amounts of runoff occur mainly in the ablation season (from May to October).However, when precipitation increases, the runoff does not increase very significantly.This is mainly because the runoff in the Tianshan Mountains is limited by temperature conditions.When only precipitation increases, runoff does not have a significant increase.Under the precipitation increases, increased runoff is concentrated mainly in June, July, and August.Within this period of time, radiation is sufficient and increased precipitation on the glacier leads to melt, which leads to increased runoff in the watershed.Increased precipitation causes a reduction in runoff from October to April (Figure 13).The increase of air temperature and precipitation has a great influence on the seasonal distribution of runoff.As the temperature increases, the increased degree of ablation decreases.The coefficient variations (Cv) of runoff increase at first and then decrease (Figure 14).However, the differences in Cv are not obvious when temperature increases   When precipitation increases by 10%, 20%, and 30% (P + 10%, P + 20%, P + 30%), monthly runoff increases by 6.1%, 11.3%, and 19.8%, respectively (Figure 13, Table 6).Increased amounts of runoff occur mainly in the ablation season (from May to October).However, when precipitation increases, the runoff does not increase very significantly.This is mainly because the runoff in the Tianshan Mountains is limited by temperature conditions.When only precipitation increases, runoff does not have a significant increase.Under the precipitation increases, increased runoff is concentrated mainly in June, July, and August.Within this period of time, radiation is sufficient and increased precipitation on the glacier leads to melt, which leads to increased runoff in the watershed.Increased precipitation causes a reduction in runoff from October to April (Figure 13).The increase of air temperature and precipitation has a great influence on the seasonal distribution of runoff.As the temperature increases, the increased degree of ablation decreases.The coefficient variations (Cv) of runoff increase at first and then decrease (Figure 14).However, the differences in CV are not obvious when temperature increases 1 °C and 2 °C, indicating that the fluctuation of runoff is not significant in this situation.When the temperature increases 4 °C, Cv decreases, and the fluctuation of the runoff decreases.With the increase of precipitation, Cv becomes greater, and the fluctuation of runoff increases, but fluctuation is lower than that caused by the changes of temperature.results is consistent with increased glacier ablation in the Tianshan Mountains and is reasonable and credible at the annual scale.
Influence of debris on glacial runoff.Rock glaciers are expected to become increasingly important as current warming continues; they are an important part of the cryospheric water store.Research suggests that they could be important for local water management [53,54].Glacier melt can be significantly altered by debris cover, but knowledge about the fraction and thickness of debris cover on the Tian Shan glaciers is still sparse [55].The structure of glacier is complex, the storage and release of water within the large complex glacier will have a great impact on the runoff process.The debris can temporarily store large amounts of meltwater, it regulates runoff at the end of the glacier in monthly and seasonal scales.The accumulation and release of this water is a very complex process in glaciers, it is not linearly applied to a regulation factor that can be accurately described.Therefore, the process of simulated results cannot be accurately fitted to the observed runoff and this is an important factor in errors of modeling.

Conclusions
In this study, variations of runoff and glacial characteristics over the past 50 years were simulated and the sensitivities of runoff under different temperature and precipitation variations in the KRB were analyzed.Temperature and precipitation were found to increase by 0.32 • C (10a) −1 and 4.8 mm•a −1 , respectively, from 1959 to 2009 at Camp-AWS in the KRB.Runoff depth increased at 5.6 mm•a −1 .The simulated GMB of the KG was approximately 370 mm•a −1 .The glaciers thinned 18.53 m from 1959 to 2009.When the temperature was increased by 1, 2, and 4 • C, the annual runoff increased by 27.4%, 56.2%, and 127%, respectively.When the precipitation increased by 10%, 20%, and 30%, the monthly runoff increased by 6.1%, 11.3%, and 19.8%, respectively.The runoff is more sensitive to variation of temperature than variation of precipitation.Increases in precipitation and temperature will lead to changes in the internal annual distribution of runoff.Temperatures, precipitation, and other climatic and topographical factors are not isolated, changes of these factors are closely related and affect the process of runoff.The drainage system and internal complex storage of large-scale glacier surfaces lead to the differences between simulated runoff and observed runoff.

Figure 1 .
Figure 1.The locations of the Koxkar River Basin KRB (a), sites of observation (four meteorological stations around the KRB are shown at top right (b) and hydrological station (c).

Figure 1 .
Figure 1.The locations of the Koxkar River Basin KRB (a), sites of observation (four meteorological stations around the KRB are shown at top right (b) and hydrological station (c).

Figure 2 .
Figure 2. (a) Elevation and (b) aspect in the Koxkar River Basin.

Figure 2 .
Figure 2. (a) Elevation and (b) aspect in the Koxkar River Basin.

Figure 4 .Figure 5 .
Figure 4. Comparison of simulated and observed daily temperature at Camp-AWS from 2007 to 2009.

Figure 4 .
Figure 4. Comparison of simulated and observed daily temperature at Camp-AWS from 2007 to 2009.

Figure 4 .Figure 5 .
Figure 4. Comparison of simulated and observed daily temperature at Camp-AWS from 2007 to 2009.

Figure 6 .
Figure 6.Comparison of simulated and observed monthly precipitation (P) at Camp-AWS from 2007 to 2009.
. We selected hydrological data from 2007/2008 and 2008/2009 for calibration.In addition, data for 2009/2010 and 2010/2011 were used for validation.

Figure 6 .
Figure 6.Comparison of simulated and observed monthly precipitation (P) at Camp-AWS from 2007 to 2009.

Figure 6 .
Figure 6.Comparison of simulated and observed monthly precipitation (P) at Camp-AWS from 2007 to 2009.
selected hydrological data from 2007/2008 and 2008/2009 for calibration.In addition, data for 2009/2010 and 2010/2011 were used for validation.

For
the period of calibration, Reff and R 2 are 0.77 and 0.76 for 2007/2008, and 0.87 and 0.84 for 2008/2009, respectively.For the period of validation, Reff and R 2 are 0.89 and 0.89 for 2009/2010 and 0.80 and 0.81 for 2010/2011, respectively (Figure

Figure 9 .
Figure 9.Comparison of monthly runoff between observation and simulation.

Figure 9 .
Figure 9.Comparison of monthly runoff between observation and simulation.

Figure 9 . 19 Figure 8 .
Figure 9.Comparison of monthly runoff between observation and simulation.

Figure 9 .
Figure 9.Comparison of monthly runoff between observation and simulation.

Figure 10 .
Figure 10.Comparison of monthly simulated (R sim ) and observed runoff depth (R obs ) from 2007 to 2011.

Table 4 .
Comparison of simulated and observed glacier mass balance accumulation between Koxkar Glacier and adjacent glaciers.
1 and Tuyuksu glaciers, thus the simulated result is reasonable and credible.Glacial ablation in the KRB increased at a rate of 4.2 mm•a −1 (p < 0.01) from 1959 to 2009, and the average elevation of glacier decreased by 18.53 m (16.68 m w.Equation), or at a rate of 370 mm•a −1 (Figure 11c).

Figure 12 .
Figure 12.Sensitivity of runoff under various temperature variations relative to the period 2007-2011.

Figure 13 .
Figure 13.Sensitivity of runoff under different precipitation variations relative to the period 2007-2011.

Figure 12 .
Figure 12.Sensitivity of runoff under various temperature variations relative to the period 2007-2011.

1 •
C and 2 • C, indicating that the fluctuation of runoff is not significant in this situation.When the temperature increases 4 • C, Cv decreases, and the fluctuation of the runoff decreases.With the increase of precipitation, Cv becomes greater, and the fluctuation of runoff increases, but fluctuation is lower than that caused by the changes of temperature.

Figure 12 .
Figure 12.Sensitivity of runoff under various temperature variations relative to the period 2007-2011.

Figure 13 .
Figure 13.Sensitivity of runoff under different precipitation variations relative to the period 2007-2011.

Figure 13 .
Figure 13.Sensitivity of runoff under different precipitation variations relative to the period 2007-2011.

Table 1 .
Types of vegetation in the Koxkar River Basin.

Table 1 .
Types of vegetation in the Koxkar River Basin.

Table 2 .
Details of the Camp-AWS and the four national meteorology stations (NMS) around the Koxkar River Basin.AWS: automatic weather station.

Table 2 .
Details of the Camp-AWS and the four national meteorology stations (NMS) around the Koxkar River Basin.AWS: automatic weather station.

Table 3 .
Definition of and set of optimized parameters of HBV for the Koxkar River Basin.

Table 3 .
Definition of and set of optimized parameters of HBV for the Koxkar River Basin.

Table 5 .
Change of runoff under various changes of temperature.

Table 5 .
Change of runoff under various changes of temperature.

Table 6 .
Change of runoff under various changes of precipitation (P).

Table 5 .
Change of runoff under various changes of temperature.

Table 6 .
Change of runoff under various changes of precipitation (P).

Table 6 .
Change of runoff under various changes of precipitation (P).