Impact of Climate Change on Drought in the Upstream Yangtze River Region

Based on Coupled Model Intercomparison Project Phase 5 (CMIP5) dataset and a variable infiltration capacity (VIC) hydrological model, this study assesses the possible influence of climate change in the upstream region of the Yangtze River on droughts in the future 30 years. Long-term daily soil moisture content were simulated by VIC model at a 50 km × 50 km resolution from 1951 to 2013. Regional historical drought events were then recognized based on soil moisture anomaly percentage index and validated with field data. Five relatively independent representative global circulation models were selected and the outputs of them were downscaled temporally and spatially as the inputs of VIC model for daily soil moisture content simulations both in the period of 1971–2000 for the present-day climate and in the period of 2021–2050 for the future. The results show that the projected annual mean temperature is likely to increase from 1.4 ◦C to 1.8 ◦C. The projected change in mean annual precipitation could be increased slightly by 0.6% to 1.3%, but the trend of precipitation change in summer and autumn might be opposite of that. Comparing the drought characteristics values recognized in 1971–2000, seven to eight additional regional drought events are likely to happen in 2021–2050. Drought duration and drought intensity are also likely to extend for 18 d to 25 d and increase by 1.2% to 6.2%, respectively. But, drought area could decrease slightly by 1.3% to 2.7% on average. These changes in drought characteristics values suggest that regional drought could become more severely prolonged and frequent in future.


Introduction
Drought is one of the major natural disasters of the world.According to the statistics from 1995 to 2014, there were up to 327 global large-scale droughts, with an average of about 5.4 million people affected annually, and about USD 430 million in economic losses.Compared with the values from 1961 to 1980, these are 2.6, 2.6, and 9.6 times greater, respectively [1].Future global climate change may lead to frequent occurrences of extreme climate and weather events [2].Climate change is projected to alter the frequency and magnitude of droughts and, in some seasons and regions, drought intensity may be further strengthened [3].
Drought occurrence and the severity of droughts are often determined by the drought index.The condition of the soil moisture content over a large area is the key to drought development, and studies have shown that a drought index based on soil moisture content can effectively reflect the changes in drought development [4][5][6].With the advance in terrestrial hydrological models and the use of long-term hydrological measurement data as constraints, the simulation of the soil moisture content has become one of the effective ways in establishing extensive long-term soil moisture content databases [7,8].A variable infiltration capacity (VIC) model [9,10] is a gridded large-scale semi-distributed land surface hydrological model widely employed in drought research.For example, Andreadis et al. [11] reconstructed the drought history over the North American continent, and Wu et al. [12] reconstructed and analyzed the temporal and spatial variations and characteristics of historical droughts from 1951 to 2009 in China.Since this method considers the integrated effects of infiltration, evaporation, vegetation, and soil properties on droughts and is capable of producing a drought index on a daily basis through the continuous simulation of daily soil moisture content, it can accurately reflect the progress and extent of drought development through relatively strict physical constraints [13].
With the data from future climate change scenarios as input, future soil moisture content data can be obtained by using hydrological models and the variation in future droughts can be studied.The climate models are the key tools for climate change projections, and the simulation results from more than 40 climate models were used in the recent IPCC AR5.However, by far the largest contribution to uncertainty in future projections stems from the fact that climate models are imperfect and therefore their projections are uncertain [14].In some applications, the use of multiple models has been shown to provide predictions superior to those of a single model [15,16].Nevertheless, studies demonstrated that in many models, the calculations of temperature and precipitation are not independent; that is, they share the same formula, parameterization processes, or methods in some cases [17][18][19].The fact that parameterization processes in many models are identical or belong to the same category highlights the concern about the dependence and relevance among the models, leading to the problem of the simulation results of some scenarios being similar and the credibility of the conclusions being affected by the number of consistent prediction results [20].Therefore, when multiple models are applied, the effects of model similarity on the results should be considered.
The upstream Yangtze River region is located in China's eastern and western political, economic, and cultural transition and convergence zone.Observational data from this region in recent years indicate a significant decrease in the autumn precipitation and runoff in Yichang at the watershed outlet and frequently recurring droughts [21].To date, there are few studies on droughts in the upstream region of the Yangtze River.Most of these employed the regional precipitation and runoff data to consider meteorological drought [22] and hydrological drought [23] problems and seldom considered the soil moisture variations.Xu [24] combined the results, such as the runoff depth in the drainage basins, actual evaporation, and surface soil moisture content simulated by a distributed hydrological model, with the Palmer drought severity index (PDSI) drought mode to construct a gridded biweekly drought assessment and forecast model, the application of which achieved a good result.
By simulating soil moisture content, this study assesses the impact of climate change on regional droughts over the upstream Yangtze River region for the future 30 years.The main study includes four steps: (1) with a large-scale hydrological simulation system in a VIC hydrological model based on the precipitation and temperature data measured at selected sites, a long-term regional soil moisture content database was simulated and established; (2) based on the simulated soil moisture content, the drought index was formulated and the historical drought events were recognized and verified with the measured data; (3) considering similarities in precipitation simulation using the climate models, several independent and representative models were selected and downscaled temporally and spatially to build future regional climate change scenarios; and (4) using the regional climate change scenarios, the VIC model was run to simulate future soil moisture content and to quantitatively assess the impact of climate change on the temporal and spatial distribution of drought characteristics, such as the number of occurrences, duration, areal extent, and intensity of drought events in the upstream Yangtze River.

Study Area
In this study, the upstream Yangtze River region refers to the area from the river mouth to Yichang in the Hubei Province (Figure 1), an area of about one million km 2 .In this study, based on the classification and division as a Class II water resource, the study area is divided into six subregions: Region I Jinsha River (above Shigu), Region II Jinsha River (below Shigu), Region III Mintuo River, Region IV Jialing River, Region V Wu River, and Region VI Yibin to Yichang.

Soil Moisture Content Simulation
This study employed a VIC model to simulate and establish the daily soil moisture content.The fundamental principles of the VIC model are referred to in the literature (http://www.hydro.washington.edu/Lettenmaier/Models/VIC/index.shtml).Nine typical hydrographs covering the upstream Yangtze River region were selected to evaluate the applicability of the VIC model.The data series lasted for 7 to 21 years.The basic information on the catchment and control stations are shown in Table 1.At a 50 km × 50 km resolution, the upstream Yangtze River region was divided into 502 grid cells.For each grid, the meteorological data required for the model input were obtained by inverse distance interpolation, and vegetation and soil parameters were set according to the global land vegetation cover database at a 1 km × 1 km resolution [25] and the global soil database at a 10 km × 10 km resolution [26], respectively.
Based on the method that integrated the Rosenbrock algorithm [27] and manual intervention, the hydrological parameters were optimized and determined.Figure 2 shows the observed and simulated

Soil Moisture Content Simulation
This study employed a VIC model to simulate and establish the daily soil moisture content.The fundamental principles of the VIC model are referred to in the literature (http://www.hydro.washington.edu/Lettenmaier/Models/VIC/index.shtml).Nine typical hydrographs covering the upstream Yangtze River region were selected to evaluate the applicability of the VIC model.The data series lasted for 7 to 21 years.The basic information on the catchment and control stations are shown in Table 1.At a 50 km × 50 km resolution, the upstream Yangtze River region was divided into 502 grid cells.For each grid, the meteorological data required for the model input were obtained by inverse distance interpolation, and vegetation and soil parameters were set according to the global land vegetation cover database at a 1 km × 1 km resolution [25] and the global soil database at a 10 km × 10 km resolution [26], respectively.
Based on the method that integrated the Rosenbrock algorithm [27] and manual intervention, the hydrological parameters were optimized and determined.Figure 2 shows the observed and simulated monthly hydrograph at the Yichang station covering 1980-2000.As observed in Figure 2, the VIC model can simulate the hydrograph relatively well.From the calibration results through the nine typical hydrographs (Table 1), the relative error for the simulated flow is between −16.2% and 21.7%, the daily flow efficiency index is between 0.60 and 0.94, and the monthly average flow efficiency index is between 0.79 and 0.97 with an average of 0.87.In general, the VIC model is capable of better simulating the watershed runoff process.The soil moisture anomaly percentage index (SMAPI) is able to describe large-scale drought characteristics [12].SMAPI is calculated as the difference between the actual soil moisture content and the corresponding long-term moisture content at a specific time divided by the corresponding long-term content.Its expression is as follows: The parameters θ and θ denote the present soil moisture content and the appropriate soil moisture content for the climate under concern at that period, respectively; and, therefore, θ can be viewed as the mathematically expected value from the time series of the soil moisture content.
The soil moisture classifications based on SMAPI are shown in Table 2.The soil moisture anomaly percentage index (SMAPI) is able to describe large-scale drought characteristics [12].SMAPI is calculated as the difference between the actual soil moisture content and the corresponding long-term moisture content at a specific time divided by the corresponding long-term content.Its expression is as follows: The parameters θ and θ denote the present soil moisture content and the appropriate soil moisture content for the climate under concern at that period, respectively; and, therefore, θ can be viewed as the mathematically expected value from the time series of the soil moisture content.
The soil moisture classifications based on SMAPI are shown in Table 2.The drought events are described by their duration, intensity, severity, and the affected areal extent.For model grids, the drought duration is defined as the number of successive days when SMAPI values are 5% below the drought threshold, the intensity is the average SMAPI value during the drought, and a drought duration of 60 days is given as the threshold for one drought event to be recognized [12].
In regional terms, the areal extent is an important characteristic in large-scale drought event recognition, and it is defined as the percentage of grids where SMAPI values are 5% below the drought threshold over the region.The drought intensity is the average SMAPI of those grids where SMAPI values are 5% below the threshold.Regional drought duration is different from that of a grid and is defined as the number of successive days when the daily drought areal extent is greater than 30%.The severity is defined the same way as that for a grid.A regional drought duration of 60 days is defined as the threshold for one regional drought event to be recognized.

Regional Climate Change Scenarios
This study applied three emissions scenarios from the representative concentration pathways (RCPs) 2011 scenario series: the low emission scenario (RCP2.6),medium emission scenario (RCP4.5),and high emission scenario (RCP8.5).The monthly mean temperature and precipitation outputs from 47 global circulation models (GCMs; see Table 3) within Coupled Model Intercomparison Project Phase 5 (CMIP5) are used in this study.We selected the baseline period of 1971-2000 as the present-day climate and the period of 2021-2050 (2030s) as the future.
The multi-model ensemble approach is used in this study in an effort to reduce the uncertainty from the internal variability and inter-model difference.Previous work [14] suggests that the behavior of the CMIP3 multi-model ensemble maybe similar to a set of probably only five to ten independent models.Considering model independence in multi-model ensemble projection, we classified CMIP5 models based on their simulation values-the amount by which models differ in some outputs in the same set of conditions reflects the level of their independence [19].
Water 2016, 8, 576 In this study, we examine the present-day climate simulations in precipitation for model independence.The main calculating procedures are as follows [28]: (1) Normalizing the differences between simulated (f ) and observational (o) data based on the observed standard deviation (σ n ), which are expressed for model m, grid-point n: For model m, the errors in (2) form spatial patterns expressed as vector: where N is the number of grid-points in the regional domain.
(2) Calculating the multi-model error pattern e: where M is the number of models.(3) Removing the effect of e and calculate the error vector d m : where (•) * indicates statistical standardization, and r is the correlation between the m th model's error filed and the e.After the removal of e, the correlation between model m error pattern and e becomes zero.(4) Calculating the correlation coefficient matrix r i,j : Following the above procedures, the average correlation coefficient matrix r i,j for precipitation in the four seasons was obtained.Based on a weighted pair-wise average distance algorithm developed for the Interactive Data Language (IDL), the hierarchical clustering analysis was carried out, and the results are shown in Figure 3.
According to the results from the model hierarchical clustering analysis (Figure 3), the 47 CMIP5 models were divided into five groups.In each group, we selected the model with the best precipitation simulation results in the present-day period and that with the projected data in the future period available under each of the RCP scenarios.In addition, if the precipitation simulation results from the models do not differ significantly, the results of the temperature simulations and the previous application of models in the study area were considered as model selection criteria.
Based on the above selection method and criteria, BCC-CSM1.1(m)(China), GISS-E2-R (USA), HadGEM2-ES (UK), IPSL-CM5A-LR (France), and MPI-ESM-MR (Germany) were chosen as the five representative models for analysis.It is assumed that these relatively independent representative models reflect the simulation quality of the study area precipitation in different model groups.
In order to solve the scale compatibility problem between the output of GCMs and the input of the VIC model, the output of the five GCMs were statistically downscaled both temporally and spatially.The detailed procedures are as follows.According to the results from the model hierarchical clustering analysis (Figure 3), the 47 CMIP5 models were divided into five groups.In each group, we selected the model with the best precipitation simulation results in the present-day period and that with the projected data in the future period available under each of the RCP scenarios.In addition, if the precipitation simulation results from the models do not differ significantly, the results of the temperature simulations and the previous application of models in the study area were considered as model selection criteria.
Based on the above selection method and criteria, BCC-CSM1.1(m)(China), GISS-E2-R (USA), HadGEM2-ES (UK), IPSL-CM5A-LR (France), and MPI-ESM-MR (Germany) were chosen as the five representative models for analysis.It is assumed that these relatively independent representative models reflect the simulation quality of the study area precipitation in different model groups.
In order to solve the scale compatibility problem between the output of GCMs and the input of the VIC model, the output of the five GCMs were statistically downscaled both temporally and spatially.The detailed procedures are as follows.
The climate model output data was subjected to inverse distance interpolation so that they could fit into 50 km × 50 km grids.Temperature and precipitation data were then corrected by using a linear-scaling approach at a monthly scale [29].
Firstly, the mean monthly values of the observed ( , obs i x ) and present-day simulated ( GCMs_c,i x ) data at month i (i = 1, 2, …, 12) were calculated in the baseline period, respectively.Then, using Equation ( 7), the ratio ( GCMs  ) between the observed and simulated values was calculated to be the revised coefficient.Finally, with Equation ( 8), the future simulated monthly value G C M s_f ,i x for the future period of 2021-2050 was multiplied by the correction coefficient to obtain the corrected simulated value G C M s_f ,i x  .
GCMs_f , GCMs, GCMs_f , Temperature correction is similar to the correction for precipitation.The revised coefficient was the ratio between the observed and simulated monthly values over the years in the baseline period.
The corrected simulated values were the product of the simulated monthly values and the correction coefficient.The climate model output data was subjected to inverse distance interpolation so that they could fit into 50 km × 50 km grids.Temperature and precipitation data were then corrected by using a linear-scaling approach at a monthly scale [29].
Firstly, the mean monthly values of the observed (x obs,i ) and present-day simulated (x GCMs_c,i ) data at month i (i = 1, 2, . . ., 12) were calculated in the baseline period, respectively.Then, using Equation ( 7), the ratio (λ GCMs ) between the observed and simulated values was calculated to be the revised coefficient.Finally, with Equation ( 8), the future simulated monthly value x GCMs_f,i for the future period of 2021-2050 was multiplied by the correction coefficient to obtain the corrected simulated value x GCMs_f,i . x Temperature correction is similar to the correction for precipitation.The revised coefficient was the ratio between the observed and simulated monthly values over the years in the baseline period.The corrected simulated values were the product of the simulated monthly values and the correction coefficient.
Temporal disaggregation is adopted from the earlier work of Wood et al. [30] and Raff et al. [31].We assumed that, a month into the future, the daily variations of the climatic parameters are the same as the time distribution in a month randomly selected from the past, allowing temperature and precipitation data series in the past to be divided monthly into four classes: "hot and wet", "cold and wet", "hot and dry" and "cold and dry".Based on the future climate change forecast data, the time solution was solved for the class that was randomly sampled from the corresponding climate class (with a random simulation of 10 times).

Past Drought Event Recognition
From 1953-2013, 12 regional drought events were recognized in the upstream Yangtze River region.Table 4 shows the drought events in the subregions in the study area, including the number of total events, the beginning and ending dates T 1 /T 2 , duration D (d), and areal extent A (%) for the top three most severe drought events.The longest duration was 124 d (about 4.1 months), while the shortest was 111 d (about 3.7 months).The areal extent was between about 42% and 51%.Table 4.The list of the total regional drought events in the upstream Yangtze River region identified by the SMAPI-based procedure during the period 1953 to 2013, with the top three most severe drought events.The regional drought events are sorted according to their severity.The beginning and ending dates (T 1 /T 2 ) of a drought event are given, as well as the duration (D), and the drought areal extent (A).Since droughts are processes that require a relatively long time to develop, they differ from other natural hazards, such as floods or earthquakes, in that the latter have well-defined and obvious beginning and ending times.Therefore, this study only recognized typical regional drought events and compared these events with actual historical drought information and relevant records in order to verify whether the recognition results were reasonable.The actual information was adopted from the study by Wang et al. [32], which summarized, organized, and analyzed all droughts in China in 2011.From 29 June 2011 to 1 November 2011, the second ranking event of Region VI and the third recognition ranking event of Region III were compared with actual drought information for verification.

Regions Total Drought Events
Figure 4 shows the variation of the geographic distribution of SMAPI.To facilitate the comparison with the data from Wang et al. [32], the boundaries of the provincial administrative districts in southwestern China were added in Figure 4.For the second ranking event of Region VI, the drought started on 29 June.As shown in Figure 4, in July, more than 70% area of Region VI suffered from mild drought.In early August, the drought area extent was decreased, but the drought intensity grew in the southern part of Region VI.From the middle of August to early September, the drought intensity and area grew gradually and severe drought occurred in some parts of Region VI.In the middle of September, although the drought area was decreased, severe droughts still occurred in the southern part of Region VI.In September, the affected population reached 925,000 and the direct economic losses equaled 1.03 billion Yuan in Chongqing.On 1 November, the regional drought event ended.For the third ranking event of Region III, the drought started on 8 August.From the middle of August to the beginning of September, the drought area expanded and the subregional areas that suffered from moderate drought grew.From the middle of September, the drought intensity weakened and the drought area diminished gradually.On 1 November, the regional drought event ended.In general, the beginning and ending dates, the variations in drought intensity, and the areal extent recognized in this study agree well with the actual information.Due to factors such as the drought definition, data coverage duration, and survey units, there is little information on verifying drought event recognition over the upstream Yangtze River region.However, through the above comparison and analysis for a typical drought event, the drought event recognition method based on SMAPI can reconstruct the situation during historic droughts relatively well, and the recognized grid results and regional drought events can be applied in analyzing temporal and spatial variation trends of drought characteristics.

Projected Changes in Temperature
All five representative models project increases in temperature by the 2030s and show a good agreement (Table 5).Changes in mean annual temperature are 1.4 5, respectively.For the seasonal mean temperature, the increases in autumn and winter are more significant and about 0.1 • C higher than the increase in annual average, while the increases in the spring and summer are less.In terms of projected changes in temperature in each subregion (Figure 5), the degree of temperature rise varies.Regions I, II, and III experience a considerable temperature increase in spring and winter, while Regions IV, V, and VI have a smaller increase.In particular, the increase in Region I is the most prominent.In summer and autumn, the changes in all subregions are similar, with a slightly greater increase in Region IV.In terms of projected changes in temperature in each subregion (Figure 5), the degree of temperature rise varies.Regions I, II, and III experience a considerable temperature increase in spring and winter, while regions IV, V, and VI have a smaller increase.In particular, the increase in Region I is the most prominent.In summer and autumn, the changes in all subregions are similar, with a slightly greater increase in Region IV.
For the projected precipitation change in each subregion (Figure 6), the percentage increase in spring varies, but is always less than 10%.Moreover, the precipitation in Region I increases for all four seasons.In this region, there is a greater increase of approximately 10% in winter but a smaller increase of approximately 5% in summer.The precipitation in Region II also increases slightly.
For Regions III, IV, V, and VI, summer and autumn precipitation decrease.In particular, Regions IV, V, and VI experience a more significant decrease in autumn precipitation.Under the RCP8.5 scenario, the decrease is more than 5%.In winter, Regions III and IV are more likely to receive more precipitation, while Regions V and VI receive less.For the projected precipitation change in each subregion (Figure 6), the percentage increase in spring varies, but is always less than 10%.Moreover, the precipitation in Region I increases for all four seasons.In this region, there is a greater increase of approximately 10% in winter but a smaller increase of approximately 5% in summer.The precipitation in Region II also increases slightly.For Regions III, IV, V, and VI, summer and autumn precipitation decrease.In particular, Regions IV, V, and VI experience a more significant decrease in autumn precipitation.Under the RCP8.5 scenario, the decrease is more than 5%.In winter, Regions III and IV are more likely to receive more precipitation, while Regions V and VI receive less.7 shows the projected changes in number of drought events by the 2030s.Under RCP2.6, RCP4.5 and RCP8.5, the number of events increases by 8 (5 to 11) in RCP2.6, 8 (6 to 9) in RCP4.5, and 7 events (5 to 9) in RCP8.5.The increase in each subregion varies.Under the RCP2.6 scenario, Region IV experiences a larger increase with 6 events (1 to 12), whereas, Region V has the least growth with 3 events (1 to 7).Under the RCP4.5 scenario, Region I experiences a smaller increase of 2 events (−5 to 6); while the other subregions have approximately 5 events.Similarly, under the RCP8.5 scenario, Region I has a slight increase of only 3 events (−2 to 7), but the number increases by at least 5 events for the other subregions.In particular, the increase in Region II is more than 8 events (4 to 14).

Projected Changes in Drought Duration
As shown in Table 8, the regional drought duration is likely to increase by the 2030s.From grid drought duration, under the RCP2.6 scenario (the first row in Figure 7), the duration increase varies for regional droughts, and the areas with longer lasting droughts occupy 80% of the drainage area under study.About 52% of the area experiences an increase within 15 d.Greater increases are found in downstream Region II and Regions IV and V, and upstream Region VI.At some localities, the increase exceeds 30 d.The areas with reduced duration are concentrated in Region IV with a drop within 15 d.Compared with the results under the RCP2.6 scenario, the results under the RCP4.5 scenario indicate fewer areas with extended drought durations.The areas with longer drought duration occupy about 74% of the drainage area.Most areas experience a duration increase of 15 d, which accounts for 44% of the entire drainage area.Downstream Region II, part of Regions IV, Region V, and upstream Region VI have larger increases, with an increase over 30 d at some localities.The areas with reduced duration are concentrated in Region I, upstream Region II, midstream and upstream Region IV, and midstream and downstream Region VI, with an increase often within 15 d.under the RCP8.5 scenario, the area with extended drought duration accounts for about 81% of the drainage area.Most of the regions have a duration increase within 15 d, occupying about 53% of the drainage area.Region I, downstream Region II, Regions V and upstream Region VI have a greater increase, some of which experience an increase of over 30 d.The areas with shorter droughts are clustered in upstream Region I and midstream Region IV, with a drop within 15 d.From grid drought duration, under the RCP2.6 scenario (the first row in Figure 7), the duration increase varies for regional droughts, and the areas with longer lasting droughts occupy 80% of the drainage area under study.About 52% of the area experiences an increase within 15 d.Greater increases are found in downstream Region II and Regions IV and V, and upstream Region VI.At some localities, the increase exceeds

Projected Changes in Drought Areal Extent
The projected changes in drought areal extent (Table 9) is likely to decrease by 1.5% (−11.1% to 3.9%) in RCP2.6, 1.3% (−10.7% to 9.2%) in RCP4.5, and 2.7% (−10.0% to 7.1%) in RCP8.5.Region IV shows a decreasing trend.In this region, the drought area decreases by 5.8% (−15.3% to 1.4%) in RCP2.6, 1.0% (−9.3% to 10.5%) in RCP4.5, and 2.1% (−10.5% to 5.7%) in RCP8.5.However, other subregions demonstrate an increasing trend, with an increase within 10%.Among them, Region III shows the greatest increase of 6.9% (−7.2% to 31.2%) under the RCP2.6 scenario, whereas, Region V has the greatest increase of 7.1% (0.7% to 11.8%) under the RCP4.5 scenario.Under the RCP8.5 scenario, Regions III and VI experience the most significant increase of 9.7% (−10.4% to 17.0% and 0.6% to 33.2%, respectively).The intensity of the regional drought is likely to increase over the upstream Yangtze River region (Table 10).The intensity increases by 2.6% (−11.2% to 1.2%) in RCP2.6, 6.2% (0.6% to 12.5%) in RCP4.5, and 1.2% (−3.9% to 5.0%) in RCP8.5, respectively.In each subregion, strengthened droughts are also noted.Among them, Region IV has the greatest intensity increase, which is 13.4% (1.4% to 43.0%) in RCP2.6, 19.4% (8.3% to 44.6%) in RCP4.5, and 14.0% (−2.0% to 43.1%) in RCP8.5.In all scenarios, the increase in Region IV exceeds 10%.Region II experiences the smallest increase.Under the RCP8.5 scenario, the rise in Region II is just 0.3% (−3.7% to 12.4%).For the variations of grid drought intensity (Figure 8) under the RCP2.6 scenario, the area with stronger drought occupies 82% of the drainage area, but the intensity increase differs.Region II demonstrates a greater increase, with an average over 30%.The areas with less intensive droughts are mostly located in upstream Region I and midstream and downstream of Region II, with a decrease within 10%.Compared with that under the RCP2.6 scenario, the output under the RCP4.5 scenario shows more areas with stronger droughts, accounting for 89% of the drainage area under study.The increase in Regions III, IV, V, and VI often exceeds 20%.The areas with reduced drought intensity are concentrated midstream and downstream of Region II, with a drop within 10% under the RCP8.5 scenario; the areas with increased intensity account for 83% of the entire drainage area, which agrees with that under the RCP2.6 scenario.The difference between the results in RCP8.5 and RCP2.6 show that the intensity increases in upstream Region I but decreases in some areas in Region V under the RCP8.5 scenario.The areas with a greater increase are mainly located in Region I and upstream Region III, while those with decreased intensity are concentrated in the midstream and downstream areas of Region II, with a reduction falling within 10%.For the variations of grid drought intensity (Figure 8) under the RCP2.6 scenario, the area with stronger drought occupies 82% of the drainage area, but the intensity increase differs.Region II demonstrates a greater increase, with an average over 30%.The areas with less intensive droughts are mostly located in upstream Region I and midstream and downstream of Region II, with a decrease within 10%.Compared with that under the RCP2.6 scenario, the output under the RCP4.5 scenario shows more areas with stronger droughts, accounting for 89% of the drainage area under study.The increase in Regions III, IV, V, and VI often exceeds 20%.The areas with reduced drought intensity are concentrated midstream and downstream of Region II, with a drop within 10% under the RCP8.5 scenario; the areas with increased intensity account for 83% of the entire drainage area, which agrees with that under the RCP2.6 scenario.The difference between the results in RCP8.5 and RCP2.6 show that the intensity increases in upstream Region I but decreases in some areas in Region V under the RCP8.5 scenario.The areas with a greater increase are mainly located in Region I and upstream Region III, while those with decreased intensity are concentrated in the midstream and downstream areas of Region II, with a reduction falling within 10%.

Discussion
Although the soil moisture content would increase with an increase of 0.6%-1.3% in precipitation, that could be not enough to offset the impacts derived from temperature increases (1.4-1.8 °C), which would lead to higher evaporation rates.For a better understanding of VIC model sensitivities to changes in temperature and precipitation, a series of hypothetical climate scenarios are constructed through altering the baseline temperature and precipitation data, e.g., temperature changes of ∆T = 0.0 °C, 1.4 °C, 1.8 °C with precipitation changes of ∆P = 0.0%, 0.6%, 1.3%.Comparing with the baseline soil moisture content, if the precipitation remained unchanged and the temperature rose from 1 °C to 2 °C, then the soil moisture content changes would decrease from −0.8% to −1.5%.Meanwhile, a precipitation increase of 0.6%-1.3%combined with unchanged temperatures are likely to result in soil moisture content merely changing from 0.1% to 0.2%.The results indicate that the drought is more sensitive to the variation of temperature than to the precipitation under the future climate scenarios.
Changes in drought frequency, duration, areal extent, and intensity suggest that regional drought would become more severely prolonged and frequent in the future.The trends of projected changes of drought characteristics agree well with the findings by Leng et al. [33].However, there is still a large uncertainty in the quantitative drought assessment on the impact of climate change.Some previous studies indicate that the largest uncertainty sources may arise from GCM structure [23,34,35].In this study, the similarity in precipitation simulation between climate models was considered for GCM selections, which could reduce the uncertainty from GCM structure to a certain

Discussion
Although the soil moisture content would increase with an increase of 0.6%-1.3% in precipitation, that could be not enough to offset the impacts derived from temperature increases (1.4-1.8 • C), which would lead to higher evaporation rates.For a better understanding of VIC model sensitivities to changes in temperature and precipitation, a series of hypothetical climate scenarios are constructed through altering the baseline temperature and precipitation data, e.g., temperature changes of ∆T = 0.0 • C, 1.4 • C, 1.8 • C with precipitation changes of ∆P = 0.0%, 0.6%, 1.3%.Comparing with the baseline soil moisture content, if the precipitation remained unchanged and the temperature rose from 1 • C to 2 • C, then the soil moisture content changes would decrease from −0.8% to −1.5%.Meanwhile, a precipitation increase of 0.6%-1.3%combined with unchanged temperatures are likely to result in soil moisture content merely changing from 0.1% to 0.2%.The results indicate that the drought is more sensitive to the variation of temperature than to the precipitation under the future climate scenarios.
Changes in drought frequency, duration, areal extent, and intensity suggest that regional drought would become more severely prolonged and frequent in the future.The trends of projected changes of drought characteristics agree well with the findings by Leng et al. [33].However, there is still a large uncertainty in the quantitative drought assessment on the impact of climate change.Some previous studies indicate that the largest uncertainty sources may arise from GCM structure [23,34,35].In this study, the similarity in precipitation simulation between climate models was considered for GCM selections, which could reduce the uncertainty from GCM structure to a certain extent by using multi-model ensembles.However, it should be noted that no general all-purpose metric satisfactorily utilizes the various outputs of GCMs to the present.In addition, based on current study, the correlations among drought characteristics are needed for future study.

Conclusions
This study evaluated the effects of climate change on the temporal and spatial distribution of drought characteristics, such as the number of drought events, drought duration, drought area, and drought intensity in the upstream Yangtze River region for the next 30 years.The main results are summarized as follows: (1) The simulated soil moisture content based on variable infiltration capacity (VIC) model was applied to construct a recognition system for historical drought events.This can better reconstruct the occurrence, development, and termination processes for typical historical drought events in order to analyze the temporal and spatial variation patterns of drought events.(2) Compared with the baseline period, the projected changes in temperature are likely to increase from 1.4 • C to 1.8 • C by the 2030s, with a more significant increase in autumn and winter.
The mean annual precipitation is likely to increase slightly by 0.6% to 1.3%.However, the seasonal mean precipitation would only increase under the RCP4.5 scenario by merely 0.8% in autumn.
In spring, the mean seasonal precipitation increase differs in the subregions.In Regions III, IV, V, and VI, the precipitation is likely to decrease in summer and autumn.In Regions V and VI, the precipitation in winter could also be decreased.(3) The projected changes in drought frequency, duration, areal extent, intensity suggest that regional drought would become more severely prolonged and frequent by the 2030s.Seven to eight additional regional drought events are likely to happen and the drought duration could extend for 18 d to 25 d.The areas with shorter droughts often have a duration drop within 15 d.The drought area decreases by 1.3% to 2.7% on average.The regional drought intensity increases by 1.2% to 6.2%, with the greatest increase of over 10% in Region IV.The drought intensity increases for more than 80% of the study area, and the increase is over 30% in some localities.
a Class II water resource, the study area is divided into six subregions: Region I Jinsha River (above Shigu), Region II Jinsha River (below Shigu), Region III Mintuo River, Region IV Jialing River, Region V Wu River, and Region VI Yibin to Yichang.

Figure 1 .
Figure 1.Location of the study area and locations of the hydrological and meteorological stations.

Figure 1 .
Figure 1.Location of the study area and locations of the hydrological and meteorological stations.

Figure 3 .
Figure 3. Hierarchical cluster of the 47 Coupled Model Intercomparison Project Phase 5 (CMIP5) models for the precipitation based on the model error correlation for the upstream Yangtze River region.

Figure 3 .
Figure 3. Hierarchical cluster of the 47 Coupled Model Intercomparison Project Phase 5 (CMIP5) models for the precipitation based on the model error correlation for the upstream Yangtze River region.

Figure 4 .
Figure 4.The geographic distribution of the Soil Moisture Anomaly Percentage Index (SMAPI) values from 29 June 2011 to 1 November 2011.

Figure 5 .
Figure 5. Projected changes in mean seasonal temperature by the 2030s relative to the baseline period.

Figure 5 .
Figure 5. Projected changes in mean seasonal temperature by the 2030s relative to the baseline period.

Figure 6 .
Figure 6.Projected changes in mean seasonal precipitation by the 2030s relative to the baseline period.

Figure 6 .
Figure 6.Projected changes in mean seasonal precipitation by the 2030s relative to the baseline period.

3. 3 .
Analysis on Drought Variation Characteristics under Future Climate Change 3.3.1.Projected Changes in Numbers of Drought Events Table The projected duration increases by 21 d (3 d to 58 d) in RCP2.6, 18 d (12 d to 26 d) in RCP4.5, and 25 d (13 d to 63 d) in RCP8.5.Each subregion increases by a different degree.Among them, Region V has the most considerable increase of 24 d (−5 d to 47 d) in RCP2.6, 17 d (5 d to 40 d) in RCP4.5, and 21 d (10 d to 33 d) in RCP8.5.The increase in Region I is the smallest, being 1 d (−25 d to 26 d) in RCP2.6, 1 d (−28 d to 30 d) in RCP4.5, and 2 d (−33 d to 28 d) in RCP8.5.Although a few models project a decrease in drought duration in some subregions, the projection results from multiple-model mean still indicate a significant increase in drought duration.
30 d.The areas with reduced duration are concentrated in Region IV with a drop within 15 d.Compared with the results under the RCP2.6 scenario, the results under the RCP4.5 scenario indicate fewer areas with extended drought durations.The areas with longer drought duration occupy about 74% of the drainage area.Most areas experience a duration increase of 15 d, which accounts for 44% of the entire drainage area.Downstream Region II, part of Regions IV, Region V, and upstream Region VI have larger increases, with an increase over 30 d at some localities.The areas with reduced duration are concentrated in Region I, upstream Region II, midstream and upstream Region IV, and midstream and downstream Region VI, with an increase often within 15 d.under the RCP8.5 scenario, the area with extended drought duration accounts for about 81% of the drainage area.Most of the regions have a duration increase within 15 d, occupying about 53% of the drainage area.Region I, downstream Region II, Regions V and upstream Region VI have a greater increase, some of which experience an increase of over 30 d.The areas with shorter droughts are clustered in upstream Region I and midstream Region IV, with a drop within 15 d.

Figure 7 .
Figure 7.The geographic distribution of projected changes in drought duration by the 2030s relative to the baseline period.

Figure 7 .
Figure 7.The geographic distribution of projected changes in drought duration by the 2030s relative to the baseline period.

Figure 8 .
Figure 8.The geographic distribution of projected changes in drought intensity by the 2030s relative to the baseline period.

Figure 8 .
Figure 8.The geographic distribution of projected changes in drought intensity by the 2030s relative to the baseline period.

Table 1 .
Description of calibration catchment characteristics and results of calibration.

Table 1 .
Description of calibration catchment characteristics and results of calibration.

Table 2 .
Soil moisture classifications based on the SMAPI (Soil Moisture Anomaly Percentage Index).

Table 2 .
Soil moisture classifications based on the SMAPI (Soil Moisture Anomaly Percentage Index).

Table 3 .
Description of 47 global circulation models from the CMIP5 models used in this study.
T1Japan Agency for Marine-Earth Science and Technology, Atmosphere and Ocean Research Institute (The University of Tokyo), and National Institute for Environmental Studies, JapanU1Atmosphere and Ocean Research Institute (The University of Tokyo), National Institute for Environmental Studies, and Japan Agency for Marine-Earth Science and Technology, Japan

Table 5 .
Projected changes ( • C) in mean annual and seasonal temperature by the 2030s relative to the baseline period.

Table 6 .
Projected changes (%) in mean annual and seasonal precipitation by the 2030s relative to the baseline period.

Table 7 .
Projected changes in total drought events by the 2030s relative to the baseline period.

Table 8 .
Projected changes (d) in drought duration by the 2030s relative to the baseline period.

Table 9 .
Projected changes (%) in drought area by the 2030s relative to the baseline period.

Table 10 .
Projected changes (%) in drought intensity by the 2030s relative to the baseline period.