Projected Effects of Climate Change on Future Hydrological Regimes in the Upper Yangtze River Basin , China

Climate change directly impacts the hydrological cycle via increasing temperatures and seasonal precipitation shifts, which are variable at local scales. ,e water resources of the Upper Yangtze River Basin (UYRB) account for almost 40% and 15% of all water resources used in the Yangtze Basin and China, respectively. Future climate change and the possible responses of surface runoff in this region are urgent issues for China’s water security and sustainable socioeconomic development. ,is study evaluated the potential impacts of future climate change on the hydrological regimes (high flow (Q5), low flow (Q95), and mean annual runoff (MAR)) of the UYRB using global climate models (GCMs) and a variable infiltration capacity (VIC) model. We used the eight bias-corrected GCM outputs from Phase 5 of the Coupled Model Intercomparison Project (CMIP5) to examine the effects of climate change under two future representative concentration pathways (RCP4.5 and RCP8.5).,e direct variance method was adopted to analyze the contributions of precipitation and temperature to future Q5, Q95, and MAR. ,e results showed that the equidistant cumulative distribution function (EDCDF) can considerably reduce biases in the temperature and precipitation fields of CMIP5 models and that the EDCDF captured the extreme values and spatial pattern of the climate fields. Relative to the baseline period (1961–1990), precipitation is projected to slightly increase in the future, while temperature is projected to considerably increase. Furthermore, Q5, Q95, and MAR are projected to decrease. ,e projected decreases in the median value ofQ95 were 21.08% to 24.88% and 16.05% to 26.70% under RCP4.5 and RCP8.5, respectively; these decreases were larger than those of MAR and Q5. Temperature increases accounted for more than 99% of the projected changes, whereas precipitation had limited projected effects onQ95 and MAR.,ese results indicate the drought risk over the UYRB will increase considerably in the future.


Introduction
e Yangtze is one of the longest rivers in the world, and the area surrounding this river is among the most developed, dynamic, densely populated, and highly concentrated industrial areas in China [1][2][3].Meanwhile, the Yangtze River Basin (YRB) is the main water resource of the South-to-North Water Transfer Project, and it is dominated by two types of monsoon current flow in a year: the Siberian northwest winter monsoon and the Asian southeast summer monsoon.ese two weather patterns are extraordinarily vulnerable to climate change due to their large seasonal and interannual variabilities in precipitation and temperature.In recent decades, the frequencies of floods and droughts in the YRB have been higher than elsewhere in China, which has led to much heavier socioeconomic losses [4][5][6][7][8].For example, the great flood of 1998 inundated 213,106 hectares (21 Mha) of land and destroyed five million houses in the YRB, causing an economic loss of over US$20 billion [5].Earth system models from Phase 5 of the Coupled Model Intercomparison Project (CMIP5) projected a warming trend of 0.47 °C/10 years and 0.73 °C/10 years over China in the 21st century under emission scenarios representative concentration pathways RCP4.5 and RCP8.5, respectively [9].ese projected climate changes may cause more natural droughts and floods, especially when coupled with local precipitation and temperature changes in this region.Additionally, population and economic growth will result in an increased demand for water for drinking and industrial activities.erefore, understanding the potential impacts of future climate change on the river discharge of the YRB has become an increasingly urgent issue of water security and is essential for assisting policymakers in adopting strategies to address extreme events.
Several previous studies have estimated the hydrological change response to climate change over the YRB.Xiao et al. [10] analyzed the intra-annual variations, tendencies, and mutational sites of runoff during the historical period in the YRB.Gu et al. [11] found that the runoff in the northern part of the Yangtze River will increase by approximately 10% from 2070 to 2099 due to climate change according to RegCM4.0 of the SRES A1B scenario.Qian et al. [12] considered terrestrial water storage changes within the YRB using the variable infiltration capacity (VIC) model under the SRES A2 and B2 climate scenarios.Koirala et al. [13] considered the runoff from 11 CMIP5 models together with a routing model and found that little change will occur in the discharge from the YRB due to the increased projected precipitation.Su et al. [9] used the VIC model and 20 CMIP5 global climate models (GCMs) and found that the total runoff is projected to increase by 10.7-21.4% during 2041-2070 under RCP8.5.Birkinshaw et al. [1] found that the projected changes in the basin's annual discharge range from −29.8 to +16.0% from 2041 to 2070 based on the results of 35 GCM models under RCP8.5.However, a large number of studies have focused on the monthly and annual mean discharges or the timings of the monthly mean discharge change during different climate scenarios rather than on the extreme discharge in the future.As reported by the IPCC [14], the increasing numbers of floods and droughts may both be associated with changes in the extremes of runoff and streamflow, which are more sensitive to climate changes within local regions, especially for large river source regions [13,15,16].e water resources of the UYRB account for almost 40% and 15% of the total water resources of the YRB and the whole of China.Possible changes in high flow and low flow over UYRB represent vital information for future water resource management and water security for the whole YRB.
erefore, in addition to considering the location, topography, and basin storage of the UYRB, projecting the climate change impacts on the spatiotemporal characteristics of hydrological regimes (mean runoff, high flow, and low flow) in the upstream region is more important for local water security.
Hydrological models driven by climate model projections are considered a reliable method for projecting future climate change and assessing its hydrological consequences [17].Statistical and dynamic downscaling methods are the main bridges spanning the resolution gap between GCMs and hydrological models.In both cases, statistical downscaling methods can identify potential sources of climate model bias and increase confidence in the simulated changes of a surface climate [18][19][20].e statistical downscaling method may provide a more appropriate approach than that of dynamical downscaling when station or extreme value points are required or when computational resources are limited [21].Within statistical downscaling methods, the quantile mapping bias-correction method is the most popular method.is method adjusts the distribution of the model output to match all statistical moments with respect to the observations.Li et al. [22] developed the EDCDF method using an adjusted bias to bias correct historical and projected model simulations, and the method proved to be efficient for reducing model biases by considering future uncertainty in climate fields [23,24].Aloysius et al. [25] and Yang et al. [26] proved that the EDCDF method is more efficient when model biases are reduced and can provide more reliable simulations and projections of the future extreme value of climate fields.
erefore, in this study, one of our main objectives is to validate the performance of the downscaled outputs of eight GCMs (precipitation and temperature) via the EDCDF method [26].For the same variables, we use the output of eight GCMs under the high and middle representative concentration pathways (RCP4.5 and RCP8.5) for the period from 2006 to 2099.In addition, coupled with the results of the eight bias-corrected CMIP5 models, the VIC model is applied to a 0.5 × 0.5 °grid cell to simulate the runoff of the whole basin from 1961 to 2099.ree hydrological regimes (mean annual runoff, high flow, and low flow) are selected to evaluate the impacts of climate change under the RCP4.5 and RCP8.5 scenarios.Finally, the contribution of precipitation and temperature to the hydrological regimes is quantified by the direct variance method.In Discussion, the potential effects of climate change on water resources and food security are discussed based on social statistical data.
is work provides a comprehensive study in the UYRB to explore the impact of climate changes on future hydrological regimes.e soil texture of the UYRB is derived from the 5-min Food and Agriculture Organization (FAO) dataset [27].Vegetation data and their characteristic parameters are based on the 1 km global vegetation dataset developed at the University of Maryland [28].All data are interpolated to a 0.5 × 0.5 °grid cell structure.

Datasets and Methods
In this study, the daily precipitation and temperature data from eight GCMs (bcc-csm1-1, CanESM2, CCSM4, CSIRO-Mk3-6-0, GISS-E2-R, MRI-CGCM3, MPI-ESM-LR, and NorESM1-M) participating in CMIP5 are used to project climate changes (Table 1). is study uses only the results from the most extreme RCP8.5 and moderate RCP4. 5 2 Advances in Meteorology emission scenarios.Based on the time series of observational data from available stations and the CMIP5 projection period, we interpolate the data for the historical period of 1961-2005 and the future projection period of 2006-2099.e daily outputs of the eight GCMs were bias corrected to the 0.5 × 0.5 °grid by using the EDCDF method and the observational data over the UYRB.

Methods.
To analyze the future changes in the variability of hydrological regimes with climate change, eight biascorrected GCMs outputs are used to force the VIC model to project the changes in extreme and average flows.We select 1961-1990 as the baseline period following the requirement of the World Meteorological Organization of using periods of at least 30 years for climate analysis.Furthermore, under the RCP4.5 and RCP8.5 scenarios, we divide the future period into three subperiods: 2010-2039 (near term), 2040-2069 (middle term), and 2070-2099 (long term).A multimodel ensemble (MME) is used to analyze the projected changes of hydrological regimes.ree hydrological regimes corresponding to the long-term availability of river discharge (MAR: mean annual runoff) and extreme flow (high flow Q 5 : exceeded 5% of the time; and low flow Q 95 : exceeded 95% of the time within a year) are selected.For each model, the MAR, Q 5 , and Q 95 of each year are calculated from the daily discharge simulation, resulting in 30 values for a 30-year period.

EDCDF Method.
e EDCDF method developed by Li et al. [22] applies a quantile-based mapping of the CDFs between both the historic and projection period and matches the climatic fields in the future projection period.For temperature, this method uses a four-parameter beta function to fit the temperature fields, and the projection period is bias corrected with the following equation: where x is a climatic variable, F −1 o−c is the quantile function corresponding to observation (o) for a historic training period (c), F m−c is the PDF of the model-simulated fields (m) for a historic training period, and p is the future projection climate fields.For precipitation, a two-parameter gamma distribution is used for the portion of a given time series with precipitation.e projected precipitation is bias corrected with the following equation: In the EDCDF method, the parametric distributions fit both the temperature and precipitation fields at each grid point and the distribution range parameters were taken as the extreme values from the data extended by half of one standard deviation of each grid cell.More details on this method are presented in Li et al. [22] and Yang et al. [26].

VIC Model.
e VIC hydrological model developed by the University of Washington [29] is widely used for hydrological simulations, and it has a strong physical basis and a more realistic representation of hydrological processes [9].e VIC model considers the dynamic variations in both the water and energy balances according to the subgrid heterogeneity represented by soil moisture storage, evaporation, and runoff generation [29][30][31].is model uses the variable infiltration curve to account for the spatial heterogeneity of runoff generation.Vegetation is described by the leaf area index (LAI), canopy resistance, and the relative proportion of roots in soil.Vertical water movement occurs within discrete soil layers through diffusion, and base flow is modeled from a lower soil moisture zone as a nonlinear recession.For the three soil layers, the top layer can quickly respond to rainfall via evapotranspiration.e second soil layer regulates infiltration and excess runoff and controls the quick flow component.e third layer controls the base flow generation.More details on the VIC-3L model are presented at the VIC website [29,30].

Statistical Method.
e multimodel ensemble (MME) average is calculated with equal weights and commonly used for reducing noise in the projections [25].e performances of the raw CMIP5 models and bias-corrected outputs against the observations are measured using the mean difference (bias) and probability of extreme values of climate fields.e bias, coefficient of determination (R 2 ), and Nash-Sutcliffe coefficient of efficiency (NSCE) [32] are used to evaluate VIC model performance.NSCE is defined as follows: where Q obs is the mean of observed discharges, Q sim is the modeled discharge, and Q obs (i) is the observed discharge at time i.

Direct Variance Method.
In this study, the direct variance method [33] is used to calculate the contribution ratio of precipitation and temperature to the runoff changes.e direct variance method for decomposition is as follows: where R o is the runoff of the baseline period (1961-1990), R pr is the simulated runoff by using RCP scenario precipitation and other input parameters for the baseline period data, and R tas is the same as R pr but for temperature.e contribution ratios of precipitation (C pr ) and temperature (C tas ) to the runoff changes are measured as follows: where the contribution ratio is between 0% and 100%.e sum of the contributions of the two climatic factors is 100%.

Performance of EDCDF Method.
We compare the spatial distributions of the bias of the monthly precipitation and temperature to verify the performance of the downscaled model outputs against the observational data (Figure 2).e results show that the bias-corrected model can dramatically reduce the bias of monthly precipitation and temperature from the original models.e MME of the original GCMs overestimates the precipitation varies from 0 to 144 mm per month in most regions but underestimates it by approximately 0.2 mm per month on some pixels of most of the east parts of the study area.e EDCDF shows remarkable skill in reducing the bias of precipitation, reducing it by −0.6 mm to −0.4 mm per month.On the contrary, the EDCDF can also obviously reduce the bias of the original model-simulated temperature.For temperature, the bias varies from −0.08 °C to 0.06 °C after downscaling, while the value for the original MME varies from −14 °C and 0 °C. is finding indicates that the EDCDF has an excellent ability to reduce the bias of the original models over the study area.
In addition, we evaluate the uncertainty of the EDCDF method's bias-corrected climate fields by using the probability plots of the extreme values of the monthly precipitation and temperature from the model simulations and observations at the 95 percent confidence level (Figures 3  and 4).In general, the EDCDF method has a strong ability to capture the precipitation and temperature distribution in the study area.However, a degree of uncertainty remains in the extreme precipitation values above the 0.9 probability level and under the 0.1 probability level compared to the observations.e EDCDF-downscaled models underestimate the extreme precipitation values in cases where the monthly precipitation is more than 140 mm per month and overestimate them in cases where the monthly precipitation is less than 10 mm per month.In the case of bias-corrected temperature, extreme values are more easily captured than 4 Advances in Meteorology they are for precipitation.Meanwhile, the bias-corrected monthly temperatures show significant agreement with the observations within the ranges of the 0.1 to 0.90 probability level.e bias-corrected models show a pattern of underestimation and a colder pattern when the monthly temperature is below 1 °C or above 19 °C.e aforementioned analyses show that bias correction is an important process during the statistical downscaling of the GCMs in regional or local climate studies.
e biascorrected data of this study can be used for further studies of climatic impact and adaptation.In the following analyses, projections from eight GCMs ensemble are built using the arithmetic mean for 30-year future time slices (2010-2039, 2040-2069, and 2070-2099) under the RCP4.5 and RCP8.5 emission scenarios.Furthermore, the changes in VIC model-simulated hydrological variables are calculated by comparing the MME outputs for the three given subperiods with that of the baseline period (1961-1990).

Calibration and Validation of the VIC Model.
e daily discharge observed at the Yichang station is used for model calibration and validation.According to the time-series discharge records of the Yichang hydrological station observations, the period of 1961-1970 is used for calibration, and the period of 1971-1990 is used for validation.e seven sensitive soil parameters (i, d 1 , d 2 , d 3 , D s , D smax , and W s ) are calibrated using the NSCE and bias of the agreement between the simulated and observed hydrographs [34].e specific meanings, ranges, and best fit values of the parameters are listed in Table 2. Finally, the validation process (1974-1987) is performed using the parameter set for the calibration period.
Figure 5 presents the calibration and validation results found using the daily discharge at the Yichang station.e NSCEs of the daily discharge used for calibration and validation are 0.708 and 0.797, respectively.e biases are −0.293 and −0.266 for the calibration and validation results, respectively.e R 2 values of the observation and simulation daily discharge are 0.889 and 0.934 for calibration and validation, respectively.e statistical analyses of the calibration and validation indicate that the model performed well in the daily discharge simulations.In the future period, we assume that the parameters used in the VIC model will not change with future climate changes, and the same parameters can be used for the hydrological simulations from 2006 to 2099 in the UYRB.

Projected Changes of Precipitation and Temperature.
Table 3 shows the change ratios of annual precipitation and temperature from the MME of the eight CMIP5 models in the near, middle, and long terms for RCP4.5 and RCP8.5.e UYRB exhibits large intermodel and interannual variability in both temperature and precipitation for the two RCPs.Under both RCPs, similar patterns of annual precipitation are projected, and a decreasing trend is projected in the near term, whereas increasing ones are projected over the middle and long terms.On the contrary, the two RCPs show different magnitudes of changes.In the near and middle terms, the two RCPs project similar magnitudes of change of the annual precipitation, showing decreases of less than 1% and increases of more than 3%.In the long term, both RCPs project a relatively greater increase against the baseline period, while RCP8.5 projects a greater increase (9.28%) than RCP4.5 (7.41%).In the future, the two RCPs are projected to show similar rising trends of the annual temperature, although the magnitudes of the changes are greater under RCP8.5 than under RCP4.5.e greatest warming, 4.88 °C, is projected in the long term under RCP8.5.
Table 4 presents the relative changes of monthly precipitation and temperature in summer and winter for the three future subperiods.
e results show that the precipitation and temperature are projected to have the same change trends under the two RCPs, but a decreasing tendency in precipitation is observed in the near term under RCP8.5.Under the RCP4.5 scenario, the monthly precipitation in summer (June, July, and August) increases by between 0.20% and 6.90%, the greatest increase is projected to occur in the long term, and the smallest increase is projected in the near term.
e RCP8.5 scenario shows variable trends in precipitation, with a decrease of −0.77% in the near term and increases of 3.23% and 8.02% in the middle and long terms, respectively.Both RCPs project increased precipitation in winter (December, January, and February) in the future.e change ratio of precipitation in the long term is larger than that of the middle and near terms.Similar to the RCP4.5 scenario, RCP8.5 also projects increased precipitation in winter, but it shows a wider potential range (2.47-20.29%)from the near term to the long term.
e projected monthly temperature under the two RCPs shows significant upward trends in summer and winter in the three future subperiods, but the ranges of the changes are narrower under RCP4.5.In addition, the magnitude of the temperature change is projected to undergo the greatest increase in summer in the long term.Under RCP8.5, the temperature in winter increases by 1.41 °C relative to the baseline value in the near term and by 5.14 °C in the long term.Note that the magnitude of the increase in temperature during winter is greater than that in summer in the near, middle, and long terms under both scenarios.

Projected Changes of Hydrological Regimes.
e relative mean changes of Q 5 , Q 95 , and MAR under two RCPs in the future are presented in Figure 6.Under both RCPs, decreases of   1961 1962 1964 1966 1968 1970 1972 1974 1976 1978 1980 1982 1984 1986 1988 1990 Observation Simulation Calibration NSCE = 0.708 Validation NSCE = 0.797     Advances in Meteorology smaller than that under RCP4.5 (−9.81%).e greatest decreases of MAR and Q 5 are projected to occur in the near term under RCP8.5, whereas the greatest decrease of Q 95 is projected to occur in the long term under the same RCP.e Q 5 presents considerably greater decreases under RCP4.5 than RCP8.5 in the long term.Q 95 decreases by a much greater magnitude and varies from −16.05% to 26.70% under RCP8.5, which is in contrast to the changes in MAR and Q 5 .Additionally, the relative changes of Q 95 in the near and middle terms under RCP4.5 are greater than those under RCP8.5.
To evaluate the relative changes of Q 5 and Q 95 over the whole UYRB, we statistically evaluate the magnitude of their changes in summer and winter for the three subperiods (Table 5).Under both RCPs, the relative changes of Q 5 are highly consistent for the three future subperiods and show decreases in the near and middle terms but increases in the long term.Q 5 is projected to undergo a greater decrease (−2.82%) and increase (1.79%) under RCP8.5 than RCP4.5 in the near and long terms.Notably, Q 5 decreases under both RCPs with similar temperature trends, although summer precipitation increases slightly (0.20%) under RCP4.5 and decreases (−0.77%) under RCP8.5 in the near term.Under both RCP4.5 and RCP8.5, similar decreasing trends of Q 95 in winter are projected for the three future subperiods.In the long term, RCP4.5 and RCP8.5 project decreases of Q 95 of −17.34% and −26.25%, respectively.e spatial patterns of the projected changes in Q 5 (high flow) in summer for the three future subperiods under RCP4.5 and RCP8.5 are shown in Figure 7. e results show that the spatial patterns of the relative changes in Q 5 are similar under the two RCPS during the three future subperiods.Under RCP8.5, the decreasing trend of Q 5 is observed in more than half of the whole basin (varying from 55.86% to 68.03%) for the near and middle terms.On the contrary, the area ratio of the decreasing trend in Q 5 reduces to 48.77% and 44.96% in the long term under RCP4.5 and RCP8.5, respectively.Specifically, the largest increase in Q 5 is projected to occur in the western upstream region of the basin, whereas the largest decrease is projected to occur in the northeastern parts of the basin.Although the spatial pattern of Q 5 is partially consistent for the two RCPs, the decreasing trends of Q 5 will occur in more regions for RCP8.5 than RCP4.5.Along the eastern border of the basin, Q 5 shows a positive tendency for RCP8.5 but a negative tendency for RCP4.5.
e spatial distribution of the changes in Q 95 (low flow) shows that RCP4.5 and RCP8.5 project similar spatial patterns of Q 95 changes during winter (Figure 8).Q 95 decreases over 75.20% of the whole basin during the near term under RCP4.5, which is a larger area than that of RCP8.5.e areas with decreasing changes tendency in Q 95 under RCP8.5 in the middle and long terms enlarge to 77.93% and 88.83% of the whole basin, respectively.In the near term, Q 95 increases in some parts of the upstream region and decreases in most other parts of the basin which are projected under the two RCPs.

Influence of Precipitation and
Temperature on Hydrological Regimes.Both precipitation and temperature have an influence on runoff.To evaluate the influence of climate fields on hydrological regimes, we analyze the contribution ratios of precipitation and temperature to the three hydrological regimes under the two RCPs (Figure 9).e results show that increases in temperature have a significant influence (more than 95%) on changes in Q 5 and MAR while precipitation plays a tiny role for the two RCPs in the near and middle terms.However, the contribution of precipitation to the Q 5 has increased by 20% for the long term.Both precipitation and temperature have almost equal contributions to Q 95 , but the ratio of temperature (more than 53%) is slightly higher than that of precipitation under both RCPs in the three future subperiods.
is finding demonstrates that the effect of temperature on hydrological regimes is more significant than that of precipitation in the study area, which is similar to the results found by Etter et al. [35] and Yu et al. [36].
We assume that the decreases in all hydrological regimes are driven by a remarkable increase in temperature, which leads to an increase in evaporation, and the slightly smaller increases in precipitation projected in the two scenarios.Furthermore, the decreases in Q 95 for the near, middle, and long terms are greater than those in MAR and Q 5 , which indicates that drought risk will increase in the future over the UYRB.
is result is confirmed by the research of Lu et al. [37], which states that the regional drought could become more severely prolonged and frequent in the future in this region according to the Soil Moisture Anomaly Percentage Index.Dai et al. [2] and Birkinshaw et al. [1] found that higher temperatures coupled with greater evapotranspiration will lead to a decrease in discharge and drought risks.

Discussion
We use the downscaled and bias-corrected outputs of eight GCMs to force the VIC model and evaluated the effects of climate change on hydrological regimes over the UYRB.Although this is a widely used and very helpful method for climate change studies, both GCMs and hydrological model simulations have many sources of uncertainty, e.g., in the predictors of GCMs and the parameters of the hydrological model.In this study, we use only those meteorological stations with observational data to bias correct the outputs of the GCMs without considering the uncertainties of the GCMs.In addition, we use the multimodel arithmetic average ensemble to analyze the climate changes and simulate river discharge.is widely used method has been considered appropriate for many studies [38,39].In this Advances in Meteorology method, the models with a poor simulation ability are averaged, and the simulation performance of each individual model is reduced.However, the optimal multimodel ensemble based on each model simulation performance requires further study.
e daily observed discharge values from the Yichang hydrological station are used to calibrate the parameters within the reasonable ranges of the VIC model, which mainly focuses on matching simulated and observed daily discharge.is process does not consider the uncertainties    Both precipitation and temperature influence the changes of Q 5 , Q 95 , and MAR.In this study, we use the direct variance method to quantify the contribution ratios of precipitation and temperature to the hydrological regimes.Under the projected increases in temperature, Q 5 , Q 95 , and MAR are projected to decrease, which will lead to increased drought risks in the future.Both climate change and human activities (e.g., human water consumption, irrigation, and reservoir) have effects on runoff, with population and economic changes showing an obviously increasing trend, especially after 2010 in the UYRB.Notably, the primary industry and secondary industry will increase and exceed the service industry, which requires more water resources.ese facts indicate that the water resource demands in this region will increase to accommodate the growing population and economy.erefore, the influences of human activities on water resources need further analysis.

Conclusions
We applied the EDCDF downscaling method to correct the bias of CMIP5 model outputs and evaluated their performances over the UYRB under two RCPs (RCP4.5 and RCP8.5).We used these tools together with the VIC model to project the potential changes of future MAR, Q 5 , and Q 95 .We considered eight GCM projections of the UYRB under RCP4.5 and RCP8.5 and examined the changes in precipitation and temperature between a historical period (1961-1990) and three future subperiods (2010-2039, 2040-2069, and 2070-2099).
e results showed that the bias-corrected models perform well in terms of the spatiotemporal variability of precipitation and temperature over the UYRB.In the future, the temperature is projected to increase, while precipitation presents a more variable shift.Under the effects of temperature and precipitation changes, MAR, Q 5 , and Q 95 present a decreasing trend in the future.
e values of both Q 5 and Q 95 decrease by a greater degree than MAR.In particular, Q 95 decreases in winter, with changes of −2.28% to −26.25% relative to the baseline value, as precipitation and temperature increase under the RCP8.5.Based on the results of the contribution ratios analysis, we found that the increase in temperature plays a decisive role in the reductions of Q 5 and annual average runoff while the increase in precipitation has limited effects under RCP4.5 and RCP8.5.e contribution of changes of temperature and precipitation to Q 95 is nearly equivalent.ese results indicate that droughts over the UYRB will increase considerably, with the MAR, Q 5 , and Q 95 showing unconspicuous decreasing trends.erefore, a sustainable water resource management strategy must be designed to address the demands of economic development, food security, and ecosystem conservation when considering future climate change.Policymakers and water management personnel need to pay more attention to the probabilities of natural hazards in the Yangtze River.

2. 1 .
Datasets.Daily observed precipitation and temperature data covering the period of 1961-2005 from 87 stations in the UYRB are provided by the National Climate Center of China Meteorological Administration.e daily discharge (1961-2005) at the Yichang hydrological station was obtained from the Yangtze River Hydrological Bureau of China.e locations of these stations and the digital elevation model (DEM) are shown in Figure 1.During the historical period, all of the stations are in operation and have the same period records of precipitation and temperature.

Figure 1 :
Figure 1: Elevation (in m) from DEM and locations of the meteorological stations and hydrological station in the Upper Yangtze River Basin (above the Yichang station).

Figure 5 :
Figure 5: Nash-Sutcliffe coefficient of efficiency (a) and deterministic coefficients (b, c) of the observed and simulated daily discharge during calibration (1961-1970) and validation (1971-1990) periods in the UYRB.

Figure 6 :
Figure 6: Change (%; vertical axis) of Q 5 (high flow), Q 95 (low flow), and MAR for RCP4.5 (blue) and RCP8.5 (pink).e box-whiskers show the 10 th , 25 th , 50 th , 75 th , and 90 th percentiles of the distributions of changes of the three hydrological indicators (high flow (Q 5 ), low flow (Q 95 ), and MAR) during the (a) near, (b) middle, and (c) long terms.e gray points denote the change values of every year from each model.

Figure 7 :
Figure 7: Relative changes (%; vertical axis) of Q 5 (high flow) in summer for RCP4.5 (a, c, e) and RCP8.5 (b, d, f ) in the (a, b) near, (c, d) middle, and (e, f ) long terms.Each figure is accompanied by a form that summarizes the percentage of the area that shows a changing magnitude.↓ denotes a decrease, and ↑ denotes an increase.

Figure 8 :
Figure 8: Relative changes (%; vertical axis) of Q 95 (low flow) in winter for RCP4.5 (a, c, e) and RCP8.5 (b, d, f ) in the (a, b) near, (c, d) middle, and (e, f ) long terms.Each figure is accompanied by a form that summarizes the percentage of the area that shows a changing magnitude.↓ denotes a decrease, and ↑ denotes an increase.

Figure 9 :
Figure 9: Contribution ratios of precipitation and temperature to Q 5 , Q 95 , and MAR for RCP4.5 (a) and RCP8.5 (b) in the near, middle, and long terms in UYRB.

Table 1 :
Information on eight global-coupled climate models.Model name Modeling center (or group) Resolution (lat × lon) bcc-csm1-1 Beijing Climate Center, China Meteorological Administration, China 2.8 × 2.8 °CanESM2 Canadian Centre for Climate Modeling and Analysis 2.8 × 2.8 °CCSM4 National Center for Atmospheric Research 1.25 × 0.9 °CSIRO-Mk3-6-0 Commonwealth Scientific and Industrial Research Organization in collaboration with Queensland Climate Change Centre of Excellence 1.875 × 1.875 °Advances in Meteorology

Table 2 :
Most sensitive parameters related to the UYRB, their description, the range used for the autocalibration, and the best fit value.
5, Q 95 , and MAR are projected in all three future subperiods, but the distribution of decreasing median values is generally wider under RCP8.5 than RCP4.5 in the middle and long terms.Furthermore, the distribution of Q 95 values in the near term is much larger under RCP4.5 than under RCP8.5.In the middle term, the decrease in MAR (−7.84%) under RCP8.5 is

Table 3 :
Change ratios of the annual precipitation and temperature in the near term, middle term, and long term from the MME in the UYRB.

Table 4 :
Change ratios of precipitation and temperature during summer and winter in the near term, middle term, and long term from the MME in the UYRB.

Table 5 :
Watershed average with the change (%; vertical axis) in Q 5 during summer and change in Q 95 during winter for RCP4.5 and RCP8.5 in the near term, middle term, and long term from the MME in the UYRB.