Abstract

Quantifying drought and wetness fluctuations is of great significance to the regional ecological environment and water resource security, especially in the fragile Qinghai-Tibetan Plateau (QTP). In this paper, the standardized precipitation evapotranspiration index (SPEI) was calculated based on the observed data and China Meteorological Forcing Dataset (CMFD) in the QTP for the period of 1979–2015, and the drought and wetness evolution based on the SPEI series and respective contribution of temperature and precipitation were also analyzed. Results indicated that meteorological stations are mainly concentrated in the eastern part of the plateau, which cannot reflect the drought and wetness trend of the whole QTP. The linear trend and Mann–Kendall test revealed that SPEI series calculated based on CMFD data at 1-, 3-, 6-, 9-, 12-, and 24-month time scales all showed significant upward trend , indicating that the QTP as a whole tended to be wetter. Spatially, the regions with significant drying and increased drought probability were mainly concentrated in the Qaidam Basin and the southern part of the QTP, and the mean contribution rates of temperature and precipitation variability to SPEI trend in these regions were 60% and −11%, respectively. The regions with significant wetting and decreased drought probability were mainly concentrated in the northeast, central, and western parts of the plateau, and the mean contribution rates of temperature and precipitation variability to SPEI trend were −9% and 61% in these regions. From the statistics in different climatic regions, most of the arid and humid regions in the QTP tended to be drier, while the semiarid regions tended to be wetter.

1. Introduction

Drought and wetness fluctuations, which are closely related to precipitation and evapotranspiration variation [1], have profound influences on agriculture and ecosystems. Over the past few decades, the global climate has changed remarkably [2]. Ground observation, remote sensing, and reanalysis data revealed that temperature is rising in most regions of the world [3, 4]. Although the spatial difference of precipitation change was relatively large [5], the frequency of extreme precipitation had increased extensively in most parts of the world [6, 7]. With the changes of temperature and precipitation, the state of drought and wetness was also changing. Globally, climate change is leading dry regions to become drier and wet regions to become wetter [8, 9], and the frequency of extreme drought and flood events would further increase in the 21st century [1012]. These extreme events can cause billions of dollars in losses and affect more people than any other climate-related disaster [13]. Therefore, it is of great significance to comprehensively understand the drought and wetness change trend.

Many meteorological drought indexes had been developed for quantifying drought and wetness fluctuations, including the Palmer drought severity index (PSDI) [14], standardized precipitation index (SPI) [15], China-Z index (CZI) [16], standardized precipitation evapotranspiration index (SPEI) [17], and soil wetness deficit index (SWDI) [18]. Among these indexes, SPI and SPEI are easy to obtain and widely used for drought and wetness changes and their impact on agriculture and ecological environment in many countries and regions of the world [1925]. The main advantage of the SPI is that it is calculated at multi-time scales and can monitor different drought types [26, 27]. However, only precipitation is used in the calculation of the SPI [17]. The index cannot reflect the influence of temperature change on the drought and wetness conditions. Compared with SPI, SPEI is calculated based on the process of surface water balance and also introduces the multi-time scale character used in SPI calculation [17, 28]. Chen et al. [29] compared the monitoring results of SPI, PDSI, and SPEI in China and found that the meteorological drought monitored by SPEI was better than that monitored by SPI and PDSI. The results obtained by Zarei et al. [30] also showed that the accuracy of SPEI in monitoring drought events is higher than that of SPI and reconnaissance drought index (RDI) in Iran. Furthermore, previous studies had also shown that SPEI variation is highly correlated with vegetation index, crop growth state, and soil water storage change [3133], so it has a wide application potential in water resource management, agricultural production, and other fields in the future.

The QTP is recognized as the water tower of Asia and can influence the Asian monsoon through its mechanical and thermal forcing [3437]. As a sensitive and fragile area, the evolution of drought and wetness in the QTP has attracted extensive attention. Shi and Liu [38] analyzed the characteristics of drought events over the continent in the Eastern Hemisphere based on the SPEI data and found that the QTP had been one of the regions with high variance. Wang et al. [39] used the high-accuracy self-calibrating Palmer drought severity index (sc_PDSI) to investigate the drought variation in China and found that northern Tibetan Plateau and southern Tibetan Plateau experienced a wetting trend between 1961 and 2009. Feng et al. [40] investigated spatial-temporal patterns of meteorological drought in the QTP and its surrounding areas and revealed that the QTP became wetter in spring during 1970–2017. Chen et al. [41] divided China into six climatic regions and analyzed the characteristics of drought and wetness trend based on SPEI, PDSI, and sc_PDSI and proposed that the QTP became wetter during 1961–2012. It should be noted that most of these conclusions are based on the observed data. However, the meteorological stations are mainly concentrated in the eastern part of the plateau, while the drought and wetness conditions in the west with few stations had received little attention. Additionally, the precipitation and temperature variability is the main factor that directly affects the change trends of drought/wetness, and previous studies had shown that the temperature and precipitation in the QTP both showed increasing trend [4244]. However, it is still not clear that which is the main factor that dominates the drought and wetness trend in the plateau.

In this study, both the observed data from 50 meteorological stations and CMFD high-resolution meteorological forcing data were used to analyze drought and wetness condition evolutions over the QTP. The main objectives of this study were to (1) assess the spatiotemporal variation characteristics of drought and wetness in the QTP during 1979 to 2015 and (2) quantify the respective contributions of temperature and precipitation variability to drought and wetness trend in the QTP. The rest of the article is arranged as follows. The basic conditions and meteorological data in the QTP and the methods used in this paper are introduced in Section 2. In Section 3, the spatiotemporal evolution of drought and wetness is analyzed, and the respective contributions of temperature and precipitation variability to drought and wetness trend are calculated. The discussion and conclusions are given in Sections 4 and 5, respectively.

2. The Study Area, Data, and Methods

2.1. Study Area and Data Sources

The QTP lies in the southwest of China, including most of the Qinghai Province and Tibet Autonomous Region and a few parts of Xinjiang Uygur Autonomous Region, Gansu, Sichuan, and Yunnan (Figure 1). The Altun Mountains, Qilian Mountains, and Qaidam Basin are located in the north of the plateau, and the Himalaya Mountains are located in the southwest. The QTP covers a total area of 257 × 104 km2, of which 122.2 × 104 km2 is permafrost area [45]. With an average elevation of more than 4500 meters, the QTP is the highest plateau in the world, so it is also known as the “Third Pole.” In addition, the QTP is also the headstream of many major rivers (including the Yangtze River, Yellow River, Lancang River, Nujiang River, and so on) in the world, and China’s largest lake (Qinghai Lake) is also located in the northeast of the plateau, so it is also known as the “Water Tower of Asia.” The climate on the plateau is alternately controlled by monsoon in the summer and westerlies in the winter [46]. Precipitation is mainly concentrated in the warm season, and the light precipitation is dominant in winter [47]. Spatially, the central and southern parts of the plateau are mainly controlled by the Indian monsoon, and the annual maximum precipitation can reach more than 1000 mm, while the western and northern parts are controlled by the westerlies, and the annual precipitation is relatively low, and the annual precipitation in the Qaidam Basin in the north is even less than 100 mm [48]. The temperature is relatively high due to low latitude and altitude in the southeast of the plateau, while in the west, the temperature is relatively low due to high altitude [49].

Observed monthly temperature and precipitation were used for SPEI calculations in this study. The observed data (1979–2015) for 50 sites in the QTP (Figure 1) were derived from the National Meteorological Information Center (https://www.data.cma.cn/). These data were obtained from the daily value statistics of meteorological elements, and the production process strictly followed the standards and methods formulated by the National Meteorological Administration of China. The accuracy of the data had been verified and also had been widely used in the research of China’s climate change and hydrological process [50, 51].

Considering the sparsity of observation sites in the QTP, the CMFD data offered by the National Tibetan Plateau Data Center (https://www.data.tpdc.ac.cn/) were also used for SPEI calculations. The CMFD data constitute a near-surface meteorological forcing dataset and cover the China domain at a 3-hourly time step and a spatial resolution of 0.1°, spanning from 1979 to 2015. It was made through fusion of ground-based observations of approximately 700 stations in China with several gridded datasets including GLDAS, MERRA, GEWEX-SRB, and TRMM 3B42 v7 [52]. The meteorological elements include near-surface temperature, pressure, specific humidity, wind speed, ground downward short-wave radiation, downward long-wave radiation, and precipitation rate. The CMFD data had been widely used in the study of the climate and hydrology of the QTP [53, 54].

2.2. SPEI for Determining Drought and Wetness Conditions

SPEI is an improved drought index of SPI designed by Vicente-Serrano et al. [17]. It is calculated based on the process of surface water balance and also introduces the multi-time scale character used in SPI calculation, which fully combines the advantages of SPI and PDSI. Because of its advantages, it had been widely used in drought monitoring in many studies [40, 41]. Based on SPEI calculations, drought and wetness can be classified into seven categories as listed in Table 1. Calculation of the index is performed using the computer program SPEI Calculator, which is developed and maintained by the Institutional Repository of the Spanish National Research Council. Software and documentation are available online for downloading (https://www.digital.csic.es/). The default method for calculating the potential evapotranspiration in the software is the Thornthwaite method. Although the potential evapotranspiration estimated by this method in winter and spring in China has certain differences from other methods, the overall estimation of the potential evapotranspiration form Thornthwaite and other methods is very comparable [40, 55, 56]. By providing a text file containing temperature and precipitation data, geographic latitude, time, and other basic information, the program can automatically calculate SPEI series at different time scales. In this study, SPEI values were calculated both based on the observed and CMFD data in the QTP for the period of 1979–2015.

According to the classification criteria in Table 1, the frequency of each drought and wetness category during a period can be obtained, and the probability of each category is calculated as follows [57]:where pi and fi are the probability (%) and frequency of ith drought and wetness category during a period, respectively, and F is the total number of data points.

2.3. Statistical Methods
2.3.1. Methods for Detecting Drought and Wetness Trends

Drought and wetness trends are indicated by the linear trends of SPEI series in this study. The Mann–Kendall test is used to determine the significance of the trends. The Mann–Kendall (MK) test does not require samples to follow a certain distribution and is also not affected by some outliers [58, 59]. In previous studies, the MK method had been widely used in the trend analysis of time series in meteorological, hydrological, and ecological fields [19, 60]. The method is based on the Z value to judge the change trend of time series and is computed as follows [40, 61].

Calculate the test statistical variable S based on the SPEI series :where sgn(x) represents the sign function and can be expressed as follows:

The variance of S can be expressed as

The standard test statistical variable Z is computed as

When the Z value is positive, it represents an increasing trend, while negative Z value means a decreasing trend. When |Z| is greater than 1.28, 1.64, and 2.32, it means that the corresponding trend of time series has passed the significance test at , , and levels, respectively.

2.3.2. Quantification of the Respective Contributions of Temperature and Precipitation

A simple method adapted based on the algorithm proposed by Wu and Chen [62] was used for quantifying the respective contributions of the temperature and precipitation variability to the drought and wetness trends, and the calculation flow of this method is as follows:(1)The monthly temperature and precipitation data from 1979 to 2015 were divided into 12 groups by month, respectively, and the linear detrending method is used for these groups. The linear detrending method is described as follows:where Yi is the detrended monthly precipitation or temperature in the ith year , yi is the original monthly precipitation or temperature in the ith year, and ϒ is the precipitation or temperature trend during 1979–2015.(2)According to the results of the previous step, the SPEI series under 3 climate scenarios can be obtained. The first SPEI time series (SPEIo) was calculated based on the original meteorological data, the second (SPEIDp) was calculated using the original temperature records and detrended precipitation data, and the third (SPEIDt) was calculated using the detrended temperature data and observed precipitation records.(3)The factors affecting the variability of SPEIo can be divided into three categories, including temperature, precipitation, and other factors. The other factors mean the factors that can affect drought and wetness conditions except temperature and precipitation. When the slope of SPEIo series (S_SPEIo) is greater than 0, it means that this grid point tends to be wetter during 1979–2015, and the contribution rates (Cr) of the three factors are calculated as follows:where S_SPEIo, S_SPEIDp, and S_SPEIDt represent the linear slopes of SPEIo, SPEIDp, and SPEIDt time series and Cr_T, Cr_P, and Cr_O are the Cr values of temperature, precipitation, and other factors. If Cr > 0, it means that this factor has a positive contribution to the wetting trend, and the greater the absolute value is, the greater the contribution is. If Cr < 0, it means that this factor has a negative contribution to the wetting trend. If Cr = 0, it means there is no contribution. When S_SPEIo<0, it means that this grid point tends to be drier, and the contribution rates are calculated as follows:

If Cr > 0, it means that this factor has a positive contribution to the drying trend. If Cr < 0, it means that this factor has a negative contribution to the drying trend. If Cr = 0, it means there is no contribution. It should be noted that only the contributions of temperature and precipitation variability to the drought and wetness trends were analyzed in this study.

3. Results

3.1. Spatiotemporal Variation of SPEI Series in the QTP

Based on the temperature and precipitation data recorded at 50 meteorological stations and 25711 CMFD grid points from 1979 to 2015, the SPEI values at 1-, 3-, 6-, 9-, 12-, and 24-month time scales were calculated. The variation trend of SPEI time series at each station and grid point is given by using the linear slope and MK test (Figure 2). As can be seen from Figures 2(a)2(f), the variation trend of SPEI time series at different time scales is consistent and the SPEI time series of more than 84% stations showed a downward trend. Since the meteorological stations are mainly concentrated in the eastern part of the QTP, it is believed that the eastern part of the plateau tended to be drier during 1979 and 2015. Figures 2(g)2(l) show the results obtained based on CMFD data. According to statistics, the grid points with a decreasing trend of SPEI are less than 38%, mainly concentrated in the southeast of QTP and Qaidam Basin, which is consistent with the conclusion obtained based on meteorological stations. However, in the central and western parts of the plateau, SPEI time series showed upward trend during 1979 and 2015, indicating that these areas tended to be wetter.

In order to characterize the drought and wetness trend in the whole QTP, the SPEI series at 1-, 3-, 6-, 9-, 12-, and 24-month time scales were averaged over the 25711 CMFD pixels, and the temporal evolution of these SPEI series is displayed in Figure 3. It can be seen from Figure 3 that the fluctuation of SPEI series calculated by small time scale is relatively large, while the fluctuation is relatively small by large time series. The SPEI values at all time scales were mainly negative before 2001 and positive after 2001. The results of the linear slope and MK test showed that all SPEI time series showed an obvious upward trend , which means that the QTP tended to be wetter as a whole during 1979 and 2015.

3.2. Probability of Drought and Wetness Occurrences

The probabilities of occurrence for each drought and wetness category at all CMFD grids in the QTP were calculated. In order to analyze the dynamic change of occurrence probability of each drought and wetness category, the whole period was divided into two subperiods: the first is from 1979 to 2000 (Period I) and the second is from 2001 to 2015 (Period II). As shown in Figure 4, the difference between probabilities of drought and wetness occurrences based on SPEI at different time scales were small in the whole plateau. During 1979–2000, the probability of mild drought was the highest while the probability of slightly wet was the highest during 2001–2015. The probabilities of extreme drought and moderate drought did not change much between the two subperiods. Compared with Period I, the probability of extreme drought increased slightly in Period II, and the probability of moderate drought decreased slightly. However, the probabilities of mild drought, slightly wet, moderately wet, and extremely wet were different between the two subperiods. In Period II, the probability of mild drought decreased, while the probabilities of slightly wet, moderately wet, and extremely wet increased.

The spatial differences of drought and wetness probability variations in the QTP were also analyzed. Considering that the probabilities calculated by different time scales were close, only SPEI calculated at 12-month time scales was used for statistics. Figure 5 shows the probability of each drought and wetness category of all CMFD pixels in 2001–2015 minus that in 1979–2000. If the value obtained by subtraction is positive, it means that the probability in 2001–2015 has increased compared with that in 1979–2000; otherwise, the probability has decreased. As shown in Figure 5(a), the probability of extreme drought in most areas of the southeastern, southwestern, and northwestern parts of the QTP increased during 2001–2015 compared with that during 1979–2000. However, except for the Qaidam Basin, Southern Xinjiang, and the southeast of Tibet, the increased probabilities of extreme drought in other regions were less than 10%. In the central and northeastern parts of the plateau, the probability of extreme drought decreased, with a decrease of less than 10% in almost all regions. From Figures 5(b) and 5(c), it can be seen that the probabilities of moderate drought and mild drought in the Qaidam Basin and the southeast of plateau increased obviously during 2001–2015, while the probability mainly decreased in other regions, but the decreased value of the probability of mild drought was higher than that of moderate drought. During 2001–2015, the probabilities of slightly wet and moderately wet showed the opposite trend. The probabilities of slightly wet and moderately wet decreased in the Qaidam Basin and the southeast of the plateau but increased in other regions (Figures 5(d) and 5(e)). It can be seen from Figure 5(f) that the probability of extremely wet increased in most other regions except the Qaidam Basin. According to this result, the probability of drought had increased and the probability of wetness had decreased obviously in the regions that tended to be drier. However, the results were opposite in the regions that tended to be wetter.

3.3. The Respective Contributions of Temperature and Precipitation Variability to the Drought and Wetness Trend

Using the method introduced in Section 2.3.2., the Cr values of the temperature and precipitation variability to the drought and wetness trends for 25711 grids in the QTP were calculated. Generally, the increase of the temperature and precipitation results in positive and negative contribution rates for the SPEI trend, respectively, in the drought-prone regions. However, in the regions that tend to be wetter, the increase of the temperature and precipitation results in negative contribution and positive contribution rate for the SPEI trend, respectively. Figures 6(a) and 6(b) show the Cr values of the precipitation and temperature variability in drought-prone regions. It can be seen that the contribution rates of temperature in these regions range from −33% to 100%, with an average value of 60%. The contribution rates of precipitation are between −100% and 50%, and the average is −11%. The increase of temperature was the main reason that the Qaidam Basin and southeast and southwest of Qinghai-Tibetan Plateau tended to be drier during 1979–2015. Spatially, the Cr values of temperature to SPEI trend in the Qaidam Basin were positive, while the contributions of precipitation were negative. However, in some areas of the southeast and southwest of the QTP, the contribution rates of temperature and precipitation to the drought trend were both positive. Figures 6(c) and 6(d) show the Cr values of the precipitation and temperature variability in the regions that tend to be wetter. The Cr values of temperature in these regions ranged from −50% to 74%, with an average value of −9%. The contribution rates of precipitation were between −46% and 100%, and the average is 61%. According to this result, the variation of precipitation was the main reason that the northeast, central, and western parts of the QTP tended to be wetter, which was different from the regions that tend to be drier. Spatially, the Cr values of temperature and precipitation are both positive in the northeast and western parts of the QTP. However, in most of the central parts, Cr values of precipitation are positive, but Cr values of temperature are negative.

4. Discussion

4.1. Comparison with Previous Studies

In previous studies, Li et al. [63] applied SPEI to characterize the drought conditions in the southeast part of the QTP during 1982–2012, and the results explicitly showed a drying trend. Chen et al. [41] calculated six drought indexes with meteorological grid data at 0.5 degree resolution, analyzed the change trend of drought and wetness in different climatic regions of China, and concluded that the QTP had an obvious trend of wetting. The above conclusions are highly consistent with ours, indicating that the conclusions obtained in this study are credible. In the western part, Li et al. [64] concluded that dryness conditions were aggravating during the period 1971–2014 across Tibet. However, we found that Tibet was showing a wetting trend in the period of 1979–2015. It should be noted that Li et al. [64] assessed the drought condition over Tibet only based on meteorological data from 38 stations, and these stations are mainly concentrated in the south of Tibet. In our study, both stations and CMFD pixels showed that drought conditions were aggravating in the same regions. However, in the north and central Tibet, most areas showed obvious wetting trends, which means that only the observed data cannot truly reflect the overall drought and wetness trend in Tibet.

4.2. Drought and Wetness Variability in Different Climate Regions

Previous studies have shown that most wet regions were becoming wetter and dry regions were becoming drier under global warming [9, 65], and the trend of drought and wetness in different climatic zones of the QTP is not clear. The China Meteorological Administrations takes the cumulative temperature with the daily average temperature of no less than 10°C and the multi-year average value of the temperature in the coldest month as the heat index and the dryness as the moisture index and divides China into 32 secondary climate regions [66], including 15 climate regions in the QTP (Figure 7). These 15 climate regions are Southern Xinjiang (IIID1), Northern Tibet (HD2), Southern Tibet (HC3), Central Tibet (HC2), Qinba (IVA2), Sichuan (VA3), Changdu (HB2), Bomi and Western Sichuan (HA1), Dawang and Zayu (HVA1), Qilian and Qinghai Lake (HC1), Southern Qinghai (HB1), Qaidam (HD1), Central Inner Mongolia (IIC2), Inner Mongolia and Gansu (IID1), and North Yunnan (VA5). Among these climatic zones, Qinba, Sichuan, Bomi and Western Sichuan, Dawang and Zayu, and North Yunnan belong to the humid regions, Changdu and Southern Qinghai belong to the semihumid regions, Southern Xinjiang, Northern Tibet, Qaidam, and Inner Mongolia and Gansu belong to the arid regions, and Southern Tibet, Central Tibet, Qilian and Qinghai Lake, and Central Inner Mongolia belong to the semiarid regions.

Previous studies had shown that the SPEI at 12-month time scale (SPEI-12) is suitable for describing the long-term drought and wetness trend [67]. In this study, we also selected SPEI-12 to analyze the characteristics of drought and wetness changes in different climatic regions of the QTP. Following the variation trend of SPEI-12 from each CMFD grid, the percentages of areas tending to be drier and wetter significantly in different climate zone were also calculated (Figure 8). As shown in Figure 8, climate regions tending to be drier significantly include Southern Xinjiang, Qinba, Sichuan, Changdu, Bomi and Western Sichuan, Qaidam, Inner Mongolia and Gansu, and North Yunnan. Climate regions tending significantly to be wetter include Northern Tibet, Southern Tibet, Central Tibet, Dawang and Zayu, Qilian and Qinghai Lake, Southern Qinghai, and Central Inner Mongolia. According to this result, most of the arid and humid regions in the QTP tended to be drier, while the semiarid regions tended to be wetter.

5. Conclusion

This study investigated the drought and wetness trend and the respective contributions of temperature and precipitation variability by using SPEI calculated based on the observed and CMFD dataset in the QTP during 1979–2015. The main results are summarized as follows:(1)The SPEI series based on the observed and CMFD data showed that the southeast of the QTP and the Qaidam Basin in the northeast tended to be drier during 1979–2015, but most areas in the middle and west of the plateau tended to be wetter. The SPEI sequence obtained by averaging all CMFD grid points revealed that the whole QTP tends to be wetter during the study period.(2)On the whole, the probability of drought was decreasing and the probability of wetness was increasing in the QTP during 1979–2015. In regions that tended to be drier, the probabilities of mild drought, moderate drought, and extreme drought were also increasing, while in areas that tend to be wetter, the probabilities of all l grades of drought were decreasing.(3)In the drier regions, temperature is the dominant factor controlling the change trend of SPEI, the average contribution rate of temperature is 60%, and the contribution rate of precipitation is only −11%. In the wetter regions, the contribution rate of precipitation is much higher than that of temperature. The average contribution rate of precipitation is 61% and that of temperature is −9%.(4)The statistical results in different climatic regions of the QTO showed that most of the arid and humid regions in the QTP tended to be drier, while the semiarid regions tended to be wetter.

Nomenclature

pi :Probability of ith drought and wetness category
fi :Frequency of ith drought and wetness category
p:Significance level of the Mann–Kendall test
SPEIo :SPEI series calculated using the original meteorological data
SPEIDp :SPEI series calculated using the original temperature and detrended precipitation
SPEIDt :SPEI series calculated using the detrended temperature and original precipitation
S_SPEIo :Linear slope of SPEIo series
S_SPEIDp :Linear slope of SPEIDp series
S_SPEIDt :Linear slope of SPEIDt series
Cr_T :Contribution rate of temperature to SPEI variation
Cr_P :Contribution rate of precipitation to SPEI variation
Cr_O :Contribution rate of other factors to SPEI variation.

Data Availability

The CMFD dataset can be downloaded directly from the National Tibetan Plateau Data Center (https://data.tpdc.ac.cn/). Climatic zone data of the QTP were adapted from the climatic regionalization data of China and can be downloaded from the Resource and Environment Data Science Cloud Platform and Data Center (https://www.resdc.cn). The observed data were provided by the National Meteorological Information Center (https://data.cma.cn/).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (42171467, 42001060, and 41705139) and the Basic Research Project of Qinghai Province (2021-ZJ-947Q). The authors also thank the National Tibetan Plateau Data Center, Resource and Environment Science and Data Center in Beijing, and China Meteorological Administration for providing the meteorological data for this study.