Study on the Spatial and Temporal Characteristics of Mesoscale Drought in China under Future Climate Change Scenarios

: In this study, precipitation, and temperature data from HadGEM2-ES under Represen-tative Concentration Pathways (RCPs) 4.5 and 8.5 were used to evaluate drought in China in the 21st century. The K-means clustering algorithm was used to analyze the regional characteristics of the dry hazard index (DHI) in China, and the impact of climate change on the variation trend and periodicity of regional drought in China was explored. The results show that the temperature and potential evapotranspiration (PET) of all clusters have an increasing trend under the two RCPs, and the precipitation of most clusters shows a signiﬁcantly increasing trend. The drought index calculated by the standardized precipitation-evapotranspiration index (SPEI) is higher than those calculated by the standardized precipitation index (SPI) and standardized effective precipitation evapotranspiration index (SP*ETI). The variation trends of drought intensity and frequency in China are not signiﬁcant in the 21st century; however, the local variation trends are signiﬁcant. The droughts in most parts of the Xinjiang Province, northern Tibet and western Qinghai Province show signiﬁcantly increasing trends. According to the DHI analyses and the variations in the drought area ratio, with increases in greenhouse gas concentrations, the droughts in central and western China will become more severe, and drought will spread to the eastern areas of China. In the case that both precipitation and temperature may increase in the future, the increase in evapotranspiration caused by temperature rise will greatly affect drought dynamics. The main drought periodicity in China in the 21st century is 1~3.6 years. Drought is affected by climate change but not signiﬁcantly.


Introduction
Drought is a long-term sustainable natural phenomenon with a wide spatial distribution [1], changing the flow characteristics of the underlying surface and affecting hydrological process. The latest studies show that with the further development of global warming, drought will be more severe in the future [2][3][4][5][6]. Therefore, there is an urgent need to build a reliable information system to predict the occurrence and distribution of drought, which may help to prevent and mitigate the associated disasters, and play an important role in hydrological forecast [7]. Generally, precipitation, temperature, and evaporation data can be used to calculate drought indices and predict drought trends.
China is a country with frequent and severe drought disasters [8]. In recent decades, China has suffered several long-term extreme droughts, which caused huge economic losses and ecological damage [9][10][11]. From 2009 to 2010, severe drought in southwestern China affected 61.3 million people, with the area of crop failures reaching 1.1 million hectares. Moreover, global warming is likely to worsen China's drought situation. Therefore, the analysis and prediction of drought development trends are conducive to forecast the changes of various elements in the hydrological process and reduce the impact of drought on the overall development of the country.
Research on the change in drought characteristics in China has mainly analyzed the temporal and spatial changes in the drought intensity, frequency, and area [5,6,[44][45][46][47][48]. Wang et al. [44] analyzed the spatial distribution characteristics of the drought frequency and intensity in Southwest China. Huang et al. [46] analyzed the spatial distribution characteristics of the drought frequency and duration in South China. Yao et al. [6] partitioned drought regions in China and analyzed the changes in the drought frequency, intensity, and duration in each region. Although the drought partition method, based on individual drought characteristics, was applied in the above research and drought analysis of each partition, a single drought index could only correspond to a situation and it was still difficult to comprehensively reflect the drought changes.
Above all, the impact of regional climate change on drought in China has been extensively studied, and many valuable results have been obtained, although there are still some deficiencies in the following aspects. (1) In the study of regional drought in China, for the observed data for China over the past 50 years, the drought evaluation effect of SPEI is better than that of SPI. However, for the drought forecast in the next 100 years, the advantages and disadvantages of various indices have not been analyzed. (2) The partition rules of drought areas have not been clear and have been mostly based on an individual method or experience. (3) Studies of the temporal variation in drought have mostly focused on the variation trend, but there is little research on the comprehensive considerations of the drought trend, frequency, area, and periodicity.
To analyze the impact of climate change on regional drought in China, the following three aspects of the research are undertaken. First, the applicability of RCMs is analyzed. Second, a partition method based on distribution features of DHI is proposed to comprehensively assess drought characteristics in China, with a comprehensive consideration of drought trends, frequency, area, and periodicity. Finally, the constructed predicting model of mesoscale drought in China is used to analyze the spatial and temporal characteristics of the study area, based on DHI partitions.

Sources of Observed Data and Study Area
The observed data are based on the daily value dataset of surface climate data in China (Version 3.0) provided by the China meteorological data sharing service network (https://www.resdc.cn/data.aspx?DATAID=263, accessed on 1 June 2020). The location of each station is shown in Figure 1

Selection of RCMs for Future Projected Precipitation Data
In this paper, an RCM based on the GCM (global climate model) model is adopted; in the GCM, there are several prediction models [49,50], so the forecast of climate change is uncertain [51]. There are many models to predict future climate change, so it is necessary to determine the proper RCMs for describing the regional climate change in China. According to the research results on the applicability of RCM [49,50,52], 25 models were selected as the preliminary RCMs. Considering tremendous computation, it is difficult to evaluate the applicability of each RCM. Based on the evaluation results of 25 preliminary RCMs in previous studies and the performance of the RCMs (see Section 3.2) tested by the comparison of the measured and simulated values, the RCMs which can describe the regional climate change in China is finally selected.
Sperber et al. [50] used the GCM to study summer precipitation in Asia and found that the simulation results of models CNRM-CM5, NCAR-CCSM4, NorESM1-M, GFDL-CM3, and GFDL-cm2.1 had little deviation from the actual rainfall. McSweeney et al. [52] used the HadGEM2-ES, CCSM4, CNRM-CM5, MPI-ESM-LR, MPI-ESM-MR, GFDL-ESM2G, GFDL-ESM2M, and GFDL-CM3 models to establish a satisfactory model of special climate conditions in South Asia. Jiang et al. [53] found that the CSIRO-Mk3. 6.0, BNU-ESM, MIROC5, NorESM1-M, HadGEM2-AO, and GFDL-ESM2M models were effective in simulating the average wind speed. Based on the spatial distribution characteristics of annual average precipitation in Central Asia simulated by Wu et al. [54], the low values of CanESM2, CNRM-CM5, and MRI-CGCM3 are relatively close to the actual low values. According to the applicability of the RCMs and the evaluation of the accuracy of various RCMs in the above literature, the RCMs to be optimized are shown in Table 1. The data in Table 1 was obtained from the website as follows (https://climate-scenarios.canada.ca/?page=gridded-data, accessed on 1 June 2020). The IDW method is used in all the models to obtain the 0.25 • × 0.25 • grid data.

Drought Indices
For most parts of China, annual precipitation is distributed unevenly in time, with more monthly monsoon rainfall and less rainfall in other periods, so a time scale of 12 months is used. The three indices calculated on a time scale of 12 months can be represented as SPI12, SPEI12 and SP*ETI12. Many studies use a time scale of 12 months to reflect the characteristics of long-term drought [5,17,18,55].
SPEI and SP*ETI can be estimated using the following equations: First, precipitation is selected as the statistic to calculate the SPEI and SP*ETI.
where i represents the series number of the ith time segment; D(i) is water shortage content on ith; P is the precipitation; and PET is the evapotranspiration, which can be calculated by the Thornthwaite method [56].
EP is the effective precipitation and can be calculated by the Soil Conservation Service (SCS) formula [57].
D is aggregated at a different time scale (n) as: The three-parameter log-logistic distribution [24] is selected as the optimum fitting distribution. The probability density function f (x) and cumulative distribution function F(x) are listed below: where α, β, and γ are scale, shape, and location parameters, respectively.
Finally, the SPEI is calculated by transforming F to the standard normal distribution as shown in Equation (5).
The detailed calculation for Equation (6) refers to the research of Vicente-Serrano et al. [24], and the selection of the fitting distribution refers to the research of McKee et al. [23].
The next process of drought analysis is to determine the key parameters. The Lmoment approach is adopted, in which Bayesian information criteria (BIC) are used to determine the optical fitting distribution.

Regionalization of Droughts
Generally, the regional characteristics of droughts are described to divide the drought area. The severity of a drought depends on the duration, intensity, frequency, and extent of specific drought episodes [43]. In this study, the DHI [43] is used to comprehensively reflect the characteristics of drought. Cluster analysis is used to analyze the distribution characteristics of the DHI in China; then, the regionalization of droughts is conducted. The calculation of DHI requires four sub-indicators [29,43], which can be shown as follows: (1) The number of drought events observed in a particular period (FRQ); (2) The number of drought events with durations greater than 24 months (FRQ24); (3) The maximum observed severity across the observed episodes (Smax); (4) The maximum duration in months across the drought episodes (Dmax).
where w 1 , w 2 , w 3 , and w 4 are the relevant weights that are estimated using an analytical hierarchical process (AHP) [29,58], and the four relevant weights in this research are 0.25, and score represents the score of the rating district of the study area. The detailed calculation methods are as follows. First, the start time, end time of drought incidents, and drought intensity are selected according to drought indices. Second, FRQ (the time length is a year), FRQ (24), Smax and Dmax are calculated. Third, the range of the subindicator weight series data of each grid in the study area is analyzed. The entire interval is divided equally into four segments, and the score is given according to the position of the data to be evaluated in the subindicator interval. Last, the detailed evaluation method provided previous research [29,43].

Cluster Analysis
The goal of cluster analysis is to classify the collected data on the basis of similarity; cluster analysis is widely used in many fields. There are many methods of cluster analysis for drought regionalization, such as the region of influence method [59], entropy-based method [60], hierarchical clustering, fuzzy clustering and K-means clustering [61,62]. After studying and comparing these methods, it is found that machine learning algorithms are realized simply and easily [61,63]. Because of the simplicity and efficiency of the K-means clustering algorithm, it has become the most widely used clustering algorithm. However, it can only be applied to continuous data, and the number of categories must be manually specified before a clustering operation. Sometimes, the local optimum cannot be avoided and local optimum is not as good as the global optimum, and the convergence speed is relatively slow in large-scale data. Therefore, the K-means clustering algorithm [64] is used to regionalize the drought area. After determining the method of drought regionalization, appropriate cluster validity indices are used to determine the number of clusters, i.e., the number of partitions. There are many clustering validity indices, such as the Dunn index [65], average silhouette width index [66] and Davies and Bouldin index (DBI) [67]. In this study, the average silhouette width index (ASWI) and the DBI are used to determine the number of clusters. Details of these indices can be found in Goyal and Gupta [63].
The objective function that the k-means algorithm attempts to minimize can be shown as Equation (8).
where, Y i − C j 2 is the squared Euclidean distance between the ith data point and jth cluster center in multidimensional space of data attributes. M is the total number of data points and K is the number of clusters. Initially, the coordinates of K-cluster centers are guessed randomly, and each point is assigned to the nearest cluster center. Cluster centers are updated by averaging the coordinates in particular clusters and reassignment is made to obtain new clusters. This process is repeated until convergence is reached. Since the final result is very sensitive to the initially used cluster center, the algorithm is repeated again and again with different initial random cluster centers for obtaining the best results.

Trend and Periodicity Analysis
The modified Mann-Kendall test is adopted in this study [68][69][70]. In this study, Fourier transform is used to analyze the periodic variation in drought; the specific theory can be seen in the research of Moreira et al. [71].
Above all, the detailed methodology of this study can be summarized as follows. First, the SPI, SPEI, and SP*ETI are selected to evaluate the drought characteristics in China. Second, based on the above four indices, the drought hazard index is used to obtain characteristic clusters in the study area, and the applicability of RCMs is analyzed. Finally, the comprehensive impact of climate change on spatial and temporal characteristics of drought in China is constructed.

Identification of Drought Clusters
Drought characteristics, such as the average severity, duration, and interval are obtained using the run approach method. According to the run approach method, the duration of drought is defined as the time period when the drought indices (SPI12, SPET12 and SP*ET12) are less than −0.5. The cumulative sum of the difference between the absolute value of SPI and the critical value of 0.5 in the drought duration is defined as the drought intensity. The values of w 1 , w 2 , w 3 , and w 4 are all set to 0.25 according to published data [29,43]. Based on the drought index, and DHI is calculated, as shown in Figure 2. K-means clustering is utilized to identify clusters in the study area. Lower values of DBI and higher values of ASWI correspond to well-separated clusters in which the DHI values are close. In contrast, the DHI values in clusters vary widely. According to historical DHI values, the average DHI is obtained ( Figure 2). The variations in DBI and ASWI corresponding to the average DHI with the partition number of drought clusters is shown in Figure 3. There are two optimal drought clusters in Figure 3, which clearly cannot describe the hydrological and meteorological characteristics of the regional climate in China. The spatial distribution differences of regional climate data in China are relatively large. Considering the differences in climate change in China, Chen and Sun optimally divided China into six drought clusters [5]. From Figure 3, the number of suboptimal drought clusters is 8 and 9. Since the clusters obtained by the two partition methods are similar and the ninth cluster is relatively small (in the partition of 9 clusters), eight clusters are used in this paper ( Figure 4).
higher values of ASWI correspond to well-separated clusters in which the DHI values are close. In contrast, the DHI values in clusters vary widely. According to historical DHI values, the average DHI is obtained ( Figure 2). The variations in DBI and ASWI corresponding to the average DHI with the partition number of drought clusters is shown in Figure 3. There are two optimal drought clusters in Figure 3, which clearly cannot describe the hydrological and meteorological characteristics of the regional climate in China. The spatial distribution differences of regional climate data in China are relatively large. Considering the differences in climate change in China, Chen and Sun optimally divided China into six drought clusters [5]. From Figure 3, the number of suboptimal drought clusters is 8 and 9. Since the clusters obtained by the two partition methods are similar and the ninth cluster is relatively small (in the partition of 9 clusters), eight clusters are used in this paper ( Figure 4).

Performance Evaluation of RCMs for Drought Simulation
Based on the observed and simulated values of rainfall and temperature, the performance of various RCMs was examined using the Taylor diagram method ( Figure 5). The values of SC and RRMSE of model 6 are better than those of model 7, so model 6 (i.e., HadGEM2-ES) is selected preferentially for further analysis. To test the performance of the HadGEM2-ES model, the measured and simulated values of annual average precipitation ( Figure 6) and temperature ( Figure 7)   the hydrological and meteorological characteristics of the regional climate in Chi spatial distribution differences of regional climate data in China are relatively larg sidering the differences in climate change in China, Chen and Sun optimally d China into six drought clusters [5]. From Figure 3, the number of suboptimal d clusters is 8 and 9. Since the clusters obtained by the two partition methods are and the ninth cluster is relatively small (in the partition of 9 clusters), eight clust used in this paper ( Figure 4).

Performance Evaluation of RCMs for Drought Simulation
Based on the observed and simulated values of rainfall and temperature, the mance of various RCMs was examined using the Taylor  values, the average DHI is obtained ( Figure 2). The variations in DBI and ASWI corre sponding to the average DHI with the partition number of drought clusters is shown i Figure 3. There are two optimal drought clusters in Figure 3, which clearly cannot describ the hydrological and meteorological characteristics of the regional climate in China. Th spatial distribution differences of regional climate data in China are relatively large. Con sidering the differences in climate change in China, Chen and Sun optimally divided China into six drought clusters [5]. From Figure 3, the number of suboptimal drough clusters is 8 and 9. Since the clusters obtained by the two partition methods are simila and the ninth cluster is relatively small (in the partition of 9 clusters), eight clusters ar used in this paper ( Figure 4).

Performance Evaluation of RCMs for Drought Simulation
Based on the observed and simulated values of rainfall and temperature, the perfor mance of various RCMs was examined using the Taylor diagram method (

Performance Evaluation of RCMs for Drought Simulation
Based on the observed and simulated values of rainfall and temperature, the performance of various RCMs was examined using the Taylor diagram method ( Figure 5). The values of SC and RRMSE of model 6 are better than those of model 7, so model 6 (i.e., HadGEM2-ES) is selected preferentially for further analysis. To test the performance of the HadGEM2-ES model, the measured and simulated values of annual average precipitation ( Figure 6) and temperature ( Figure 7) are obtained by the IDW method and compared. The linear deviation correction method [72] is used for data processing. The variation processes of monthly precipitation and temperature in the clusters from January 1959 to December 2005 are shown in Figures S1 and S2, respectively. The calculated and observed values of precipitation and temperature in the clusters of the HadGEM2-ES model fit well.
To analyze the impact of simulated precipitation and temperature of HadGEM2-ES on reflecting drought, with the drought clusters in Figure

Trend Analyses of Drought Severity and Frequency
The time span from 1960 to 2100 is divided into three periods: the historical period (Hist, 1960(Hist, -2005, the near-future period (NF, 2021-2050), and the far future period (FF, 2071-2100), unit of frequency is times/year. Mann-Kendall and Sen's slope tests were used to evaluate drought severity and the variation trend of drought frequency over the above three periods, as shown in Figure Table 3.

Trend Analyses of Drought Severity and Frequency
The time span from 1960 to 2100 is divided into three periods: the historical period (Hist, 1960(Hist, -2005, the near-future period (NF, 2021-2050), and the far future period (FF, 2071-2100), unit of frequency is times/year. Mann-Kendall and Sen's slope tests were used to evaluate drought severity and the variation trend of drought frequency over the above three periods, as shown in  Table 3.

Trend Analyses of Drought Severity and Frequency
The time span from 1960 to 2100 is divided into three periods: the historical period (Hist, 1960(Hist, -2005, the near-future period (NF, 2021-2050), and the far future period (FF, 2071-2100), unit of frequency is times/year. Mann-Kendall and Sen's slope tests were used to evaluate drought severity and the variation trend of drought frequency over the above three periods, as shown in  Table 3.    (Table 2). Except for Cluster 4 and Cluster 5, precipitation in other areas shows an obvious increasing trend. In addition, under RCP 4.5, the Sen's slope of PET is higher than that of precipitation only in Cluster 1, Cluster 4 and Cluster 5. Under RCP 8.5, the Sen's slope of PET is commonly higher than that of precipitation except in Cluster 2 and Cluster 3.

Trend Analyses of Drought Severity and Frequency
The time span from 1960 to 2100 is divided into three periods: the historical period (Hist, 1960(Hist, -2005, the near-future period (NF, 2021-2050), and the far future period (FF, 2071-2100), unit of frequency is times/year. Mann-Kendall and Sen's slope tests were used to evaluate drought severity and the variation trend of drought frequency over the above three periods, as shown in Figure 8 Table 3. The variation trend of drought severity and the frequency of the clusters in China over the 21st century can be summed as follows. The drought severity and frequency of most areas of Cluster 1 show an increasing trend over the historical period. According to the SPI and SP*ETI, drought severity shows an increasing trend under RCP 4.5. Drought severity will decrease over NF and increase over FF under RCP 4.5, based on the SPEI analysis results. Under RCP 8.5, both the SPI and SPEI drought indicators show that the drought severity of Cluster 1 shows a decreasing trend.
For Cluster 2, based on the analysis results of the SPEI and SP*ETI, the drought severity decreased under RCP 4.5. Over NF and FF, the drought severity and frequency of Cluster 1 and Cluster 2 increase significantly in Northern Xinjiang Province, Northwestern Gansu Province, and western Inner Mongolia.
For Cluster 3, drought severity and frequency show a decreasing trend over Hist. However, drought frequency increased over NF and FF based on the analysis results of the SPEI and SP*ETI. Under RCP 4.5, drought severity does not show obvious variation trends over NF and FF.
The drought severity and frequency of Cluster 4 and Cluster 5 in northeastern China show an increasing trend over Hist. According to the SPI and SPEI, drought severity over NF and FF shows an increasing trend.
The drought severity and frequency values of Cluster 6, Cluster 7, and Cluster 8, located in central and southern China, show a decreasing trend over Hist. According to the SPI drought indicator, the drought severity and frequency of the 3 clusters will have an increasing trend over NF and FF. However, the drought severity and frequency of the three clusters will decrease over NF and FF based on the analysis results of the SP*ETI and SPI (Table 3).

Trend Analyses of Drought Severity and Frequency
The time span from 1960 to 2100 is divided into three periods: the historical period (Hist, 1960(Hist, -2005, the near-future period (NF, 2021-2050), and the far future period (FF, 2071-2100), unit of frequency is times/year. Mann-Kendall and Sen's slope tests were used to evaluate drought severity and the variation trend of drought frequency over the above three periods, as shown in Figure 8. All drought indicators over Hist show a similar pattern. The average values of the Sen's slope of clusters corresponding to Figure 8a-c are shown in Table 3           Generally, the variation trend of drought severity and frequency in China in the 21st century is not obvious, but the local variation trends are much more significant.

The Variation in the Drought Area Ratio
The variations in drought area with time under different conditions are shown in Figure 9. The calculated results of the modified Mann-Kendall test are shown in Table 4. Compared with the varying curves with the drought indicators of the SPI, SPEI and SP*ETI, the drought area ratio calculated with the SPI is typically lower than that calculated with the SPEI and SP*ETI, indicating that the drought hazard is more serious when the effective precipitation and potential evaporation are considered. The drought area ratio based on the evaluation results of the SPEI is usually greater than that of the SPI and SP*ETI, which indicates that potential evaporation related to temperature has a significant impact on drought severity. As shown in Figure 9b,c, the drought areas of all clusters over NF show an increasing trend, which can be proven by the positive values of the modified Mann-Kendall test in Table 4. According to the negative values of the Mann-Kendall test for Cluster 1, Cluster 2, and Cluster 6 shown in Table 4, the drought areas in these clusters over FF show a decreasing trend.
Generally, drought areas will first increase and then decrease in most parts of China over the period from 2006 to 2100; the most remarkable area is Cluster 1.

Analyses of the DHI Variations
To study the temporal and spatial characteristics of drought from 1960 to 2100, the DHI distribution over Hist  and TF (2006-2100) is shown in Figure 10. As shown in Figure 10, the DHI over TF is greater than that over Hist, which indicates that drought in most of China will have an increasing trend. The DHI distribution under RCP

Analyses of the DHI Variations
To study the temporal and spatial characteristics of drought from 1960 to 2100, the DHI distribution over Hist  and TF (2006-2100) is shown in Figure 10. As shown in Figure 10, the DHI over TF is greater than that over Hist, which indicates that drought in most of China will have an increasing trend. The DHI distribution under RCP 4.5 indicates that drought in the central and western areas of China, especially in northern Xinjiang Province, northwestern Gansu Province, and western Inner Mongolia, shows an increasing trend in the future. The DHI distribution under RCP 8.5 indicates that drought will spread to East China as the drought severity increases in the central and western areas of China. The DHI value calculated by the SPEI is greater than those calculated by the SPI and SP*ETI, which suggests that a predicted drought is more severe when considering potential evaporation ( Figure 10 and Table 5). This finding shows that a higher emission concentration will aggravate drought, and the effect on temperature will be more obvious. Considering the potential evapotranspiration related to temperature, with a rise in temperatures, the predicted drought would be more severe, and the drought hazard would be greater. of China. The DHI value calculated by the SPEI is greater than those calculated by the SPI and SP*ETI, which suggests that a predicted drought is more severe when considering potential evaporation ( Figure 10 and Table 5). This finding shows that a higher emission concentration will aggravate drought, and the effect on temperature will be more obvious.
Considering the potential evapotranspiration related to temperature, with a rise in temperatures, the predicted drought would be more severe, and the drought hazard would be greater.   Figure 10. The distribution of the drought hazard index for different clusters.

Periodicity Analyses
Periodograms of the SPI, SP*ETI and SPEI series were calculated to examine the hidden cycles of drought in different clusters for each grid point across China. To analyze the change in periodicity over time, three time periods, including historic , future period 1 (2021-2050) and future period 2 (2071-2100), were used. As shown in Figure 11, the percentage of the total area with a periodicity of 1~3.6 years is more than 35%, and the percentage of the total area with other grouped periodicity values is mostly less than 20%. In addition, the area percentages of some regions (Cluster 1-Cluster 6) over different periods are similar, which indicates that the periodicity of these clusters does not change with time. In Cluster 7 and Cluster 8, the percentage of the cluster area with periodicities of 1-3.6 years, 3.6-5 years, and 5-10 years increases to a certain extent. These areas are close to the Pacific Ocean and may be influenced by the El Niño phenomenon. The grouped periodicities of the clusters under RCP 4.5 and RCP 8.5 are similar, indicating that carbon emission concentrations have a certain influence on drought periodicity. Generally, climate variation has a certain influence on drought periodicity, while the influence of climate variation on other areas is not substantial.

Spatial and Temporal Characteristics of Drought in China
Studies have shown that the intensity, frequency and periodicity of drought in areas of China have changed significantly due to changes in climate factors, such as precipita-

Spatial and Temporal Characteristics of Drought in China
Studies have shown that the intensity, frequency and periodicity of drought in areas of China have changed significantly due to changes in climate factors, such as precipitation, temperature and wind speed [1, 5,6,48,73]. Using various CMIP5 models to predict future precipitation, it is concluded that the annual precipitation in most areas of China will increase in the 21st century [6,16,74,75]. Some scholars believe that an increase in precipitation will lead to a decrease in drought. For example, according to the research of Lin et al. [76], the increase in precipitation would decrease the drought in Northwest China. The research results of Zarch et al. [16] indicate that changes in climate will alter the areal extents of aridity zones in the future. Moreover, based on the research of predecessors and the work in this study, the predicted precipitation under RCP 8.5 is greater than that under other RCPs [77,78]. Temperature anomalies are also considered to play an important role in drought occurrences, greatly amplifying the evaporative demand and thereby increasing the overall drought intensity and impact [4,[19][20][21]. Although precipitation will increase in the future, increasing temperatures will also promote evaporation and lead to more severe droughts. For example, the droughts in most parts of the Xinjiang Province, northern Tibet, and western Qinghai Province show significantly increasing trends.

Impact of Selection of the Drought Index on Drought Characteristics
It is necessary to be cautious in the selection of drought indicators. For example, the research of Oguntunde et al. [79] suggests that the drought evaluation results of SPI and SPEI are similar in the short term, but will diverge in the long term. According to the results of Ma et al. [73], Sobol's global sensitivity principle reveals that the frequency of drought events induced by PET is greater than the precipitation deficit in the eastern monsoon, but the frequency of drought events dominated by precipitation is higher than that of PET in the northwest arid region. SPEI evaluates a drought as much more severe than the same drought evaluated by SPI and SP*ETI. The research results of this study show that the temperature in China will increase in the future, and evaporation will continue to be aggravated. Therefore, drought caused by increasing temperatures must be considered, and the drought indicators combining precipitation and evaporation should be optimized.
Some limitations of the current study should be acknowledged. In the calculation process of SPEI, potential evaporation based on temperature is used instead of evaporation, which increases the prediction of actual evaporation and aggravates the drought. The selection of the calculation method of PET might affect the evaluation results of drought. According to the research of Zhou et al. [80], PET can be estimated by a set of models. The set consists of four temperature-decisive models (Hamon, Hargreaves-Samani, Oudin, Thornthwaite), two radiation-decisive models (Energy-Only and Priestley-Taylor) and two synthesis models (Penman and Penman-Monteith). The research results show that both P-M and Thornthwaite [56] are suitable for assessing drought conditions on the national scale, but the drought assessed by the Thornthwaite method is worse than that by the Penman-Monteith method in 1990-2017 when warming is faster than previously in the arid region. It is impossible to accurately obtain the actual regional evaporation and analyze whether the calculated value is close to the actual value. Thus, we can only determine that the calculated evaporation by a certain method is relatively close to the average value. The Penman-Monteith method accounts for changes in energy, humidity, and wind speed. Due to the limitations of the data availability for the CMIP5 simulations, the Thornthwaite equation instead of the Penman-Monteith method is used to calculate evapotranspiration in this study. Some early studies [21] have indicated that the role of temperature is exaggerated for droughts over North China when using the Thornthwaite equation, compared to the Penman-Monteith method; however, there is almost no difference for South China. Therefore, we believe that using the Thornthwaite equation in this study provides robust results, although the increase in droughts in response to anthropogenic warming may be somewhat overestimated for North China.

The Impact of Carbon Emission on Drought
Previous Intergovernmental Panel on Climate Change (IPCC) assessment reports have performed assessments of rising temperature and possible risks in the climate system under various emission scenarios. Some research results and this study found that the precipitation of RCP8.5 is greater than that of the other RCPs. In this study, the analysis of GCM simulated data shows that the high carbon emission leads to the more obvious increasing trend of temperature and PET and to the increase in precipitation of most clusters ( Table 2, Figures S6-S8). Some previous studies suggest that high carbon emissions would increase the severity of regional droughts. According to previous research, with the comparison of two scenarios, the severity of different drought hazards, including the drought frequency, duration, and trends, is higher under RCP 8.5 than RCP 4.5. The research of [47], indicates that the occurrence probability of hot drought events (SPEI less than −1.0) will increase to nearly 100% by the year 2050, even though the future annual precipitation is projected to increase under the RCP 4.5 scenario, and these conditions would become even worse under the RCP 8.5 scenario. In addition, for the Pearl River Basin, drought events lasting 1-2 months would be decreased by 7.7%, those lasting 3-4 months would be increased by 4.3%, and those lasting more than five months would be increased by 3.4% under RCP 4.5, respectively.
In the present study, carbon emissions are shown to affect the drought intensity and frequency in China, which varies greatly in different regions. For example, the drought intensity and frequency in eastern Xinjiang Province, the western part of Inner Mongolia, and northern Qinghai Province show significantly increasing trends. However, under RCP 8.5, the drought intensity and frequency of most regions of Cluster 4 and Cluster 5 show significantly decreasing trends in the short term. The drought frequency shows a significantly decreasing trend, while the drought intensity shows a significantly increasing trend in the long term. Therefore, in the analysis of carbon emission, attention should be paid to its spatial and temporal characteristics. The cluster with high carbon emissions may not aggravate drought conditions, and some areas may even show decreasing trends of drought.

Periodicity Analysis
Few studies have examined the drought periodicity in China, and most of them have focused on drought frequency analyses. The research of Ma et al. [73] suggests that the frequency of drought events in China is consistent with the global change trend. The drought frequency gradually declines from 1961 to 1980 and increases from 1981 to 2017. The analysis of periodicity shows that the main periodicity of drought in areas in China is 1-3.6 years, and the percentage of coastal areas with this periodicity will increase. Wang et al. [81] found that sea-level cycles of 3 and 5 years in the Pearl River Estuary are related, and Feng and Li [82] suggested that enhanced rainfall is associated with the decaying phase of EP El Niño events in southern China, whereas reduced rainfall occurs when CP El Niño events arise in the spring following maturity.

Uncertainties Analysis
It is well known that there are many uncertainties in the assessment of climate change impacts. The uncertainties in GCM projections result from errors in the model structure, scenarios, and initial conditions. Uncertainty due to the GCM structure would make different models produce different climate projections for the same emission scenarios; model predictions can also vary significantly over time.
The Third national assessment report of China on climate change indicates that one of the reasons for uncertainty is the lack of information. Due to few observation stations in western China (especially in the Qinghai-Tibet Plateau Region) and the lack of hydrology meteorological observation information, climate evaluation results of GCM simulations lack accuracy. As shown in Figures 6 and 7, the fewer data that are available, the more uncertainty there is in the simulation results. A similar problem also arises in the research of [83]. The largest errors appear in western China, based on the validation with the China-CMIP5 dataset during the period from 1976 to 2005.
Because it is difficult to estimate carbon emissions in the future, many scenarios have been set up in previous studies to consider possible circumstances. Studies, such as [84] assume an instantaneous removal of all anthropogenic CO 2 ; [85] prescribes a 1% per year decrease in atmospheric CO 2 concentrations from quadrupled preindustrial CO 2 levels, while [86] prescribes a decline in CO 2 concentrations that reflects atmospheric CO 2 in the RCPs.
The uncertainty of temperature and precipitation data leads to the uncertainty of drought assessment in GCM projections. Selection of the drought index, grid-scale, and data correction method also affect the accuracy of a drought assessment. In the future, a better understanding of the hydrological cycle may lead to better and less uncertain projections, thus resulting in better future climatic predictions that may reduce uncertainty in drought analysis.

Conclusions
Based on precipitation and temperature data of HadGEM2-ES under RCP 4.5 and RCP 8.5, the variations in drought during the 21st century in China are predicted. The impact of climate change on the variation trend and periodicity of regional drought in each cluster of China is explored. The main conclusions in this study are drawn as follows: (1) Three drought indices (SPI, SP*ETI, and SPEI) in the historical period are close, diverging in the future period. (2) Droughts in China over different periods have regional characteristics. All drought indices show similar variation trends. Droughts in most parts of Xinjiang Province, northern Tibet, and western Qinghai Province show a significantly increasing trend. (3) The effect of climate change on the drought cycle is not significant. With the increase in greenhouse gas concentrations, drought in central and western China will become more severe, and it will spread to the eastern areas in China. A higher emission concentration will lead to an increase in the potential evaporation rate, and the corresponding evapotranspiration will lead to a higher drought hazard. (4) in the study of drought changes in the future, more consideration should be given to comprehensively reflecting the regional drought characteristics by using various drought indices related to temperature and evapotranspiration.
The limitations of this study are as follows: (1) potential evapotranspiration based on temperature is used instead of the actual evapotranspiration data; (2) selection of the PET calculation formula and uncertainty regarding the drought assessment affect the GCM projections. Further efforts are needed to be made to explore and study methods that can more accurately calculate actual evapotranspiration and simulate future climate change to assess droughts more accurately.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/w13192761/s1, Figure S1: Comparison of cluster averages of observed precipitation and precipitation ofHadGEM2-ES for various clusters. The blue line indicates observed data and the red line denotes HadGEM2-ES data, Figure S2: Comparison of cluster averages of observed temperature and temperature of HadGEM2-ES for various clusters. The blue line indicates observed data and the red line denotes HadGEM2-ES data, Figure S3: Comparison of cluster averages of observed SPI12 and SPI12 of HadGEM2-ES for various clusters. The blue line indicates observed data and the red line denotes HadGEM2-ES data, Figure S4