The monthly standardized precipitation evapotranspiration index (SPEI) can
be used to monitor and assess drought characteristics with 1-month or
longer drought duration. Based on data from 1961 to 2018 at 427
meteorological stations across mainland China, we developed a daily SPEI
dataset to overcome the shortcoming of the coarse temporal scale of monthly
SPEI. Our dataset not only can be used to identify the start and end dates
of drought events, but also can be used to investigate the meteorological,
agricultural, hydrological, and socioeconomic droughts with a different timescales. In the present study, the SPEI data with 3-month (about 90 d)
timescale were taken as a demonstration example to analyze spatial distribution
and temporal changes in drought conditions for mainland China. The SPEI
data with a 3-month (about 90 d) timescale showed no obvious intensifying
trends in terms of severity, duration, and frequency of drought events from
1961 to 2018. Our drought dataset serves as a unique resource with daily
resolution to a variety of research communities including meteorology,
geography, and natural hazard studies. The daily SPEI dataset developed is
free, open, and publicly available from this study. The dataset
with daily SPEI is publicly available via the figshare portal (Wang et al.,
2020c), with
Drought is one of the most destructive natural hazards worldwide. It can lead to adverse effects to the ecological system, industrial production, agricultural practices, drinking water availability, hydrological processes, and water quality (Bussi and Whitehead, 2020; Lai et al., 2019; Vicente-Serrano et al., 2012; Wang et al., 2014, 2017). Drought caused ca. USD 221 billion in loss from 1960 to 2016 reported by the International Disaster Database (EM-DAT), and the drought events in South Asia influenced over 60 million residents from 1998 to 2001 (Agrawala et al., 2001). Unfortunately, drought is expected to increase in frequency and intensity due to the future warming air temperature (Trenberth et al., 2014; Zambrano et al., 2018). The exacerbated drought conditions have promoted some national legislation (such as drought preparedness and planning) to carry out risk management and adaptive strategy for drought disasters (Garrick et al., 2017).
The various drought types result in the difficulty of drought characterization and assessment. Drought definition is not unique. Some studies proposed defining drought according to the water deficit (Wilhite and Glantz, 1985), while others defined drought based on the period of abnormal arid conditions (Eslamian et al., 2017). The most common drought can be classified into four types including (1) meteorological, (2) agricultural, (3) hydrological, and (4) socioeconomic (Mishra and Singh, 2010). Meteorological drought results from precipitation deficit or evaporation increases (McKee et al., 1993). Meteorological drought can propagate into agricultural drought with lower soil moisture availability, and it also can lead to hydrological drought with lower streamflow and socioeconomic drought with lower water availability (Barella-Ortiz and Quintana-Seguí, 2019; Gevaert et al., 2018). In general, drought indices are normally used to monitor and assess the condition or spatial–temporal characteristics of drought.
Many drought indices have been developed for drought characterization and assessment, such as the Palmer drought severity index (PDSI) (Dai et al., 2004), standardized precipitation index (SPI) (McKee et al., 1993), vegetation water supply index (VWSI) (Carlson et al., 1994), vegetation health index (VHI) (Kogan, 2002), vegetation temperature condition index (VTCI) (Wan et al., 2004), and other drought indices (Men-xin and Hou-quan, 2016; Wang et al., 2015, 2017). PDSI and SPI are the most popular drought studies worldwide (Dai et al., 2004; McKee et al., 1993); however, they have some limitations. PDSI is only suitable to the agricultural drought through characterization of the soil water deficit, and it cannot identify the meteorological, hydrological, and socioeconomic droughts (Feng and Su, 2019). In addition, PDSI limits the spatial comparability of drought due to the fact that it is heavily dependent on data calibration (Sheffield et al., 2009; Yu et al., 2014). Although the SPI can be used to monitor and assess different drought types on multiple spatial scales at a monthly time step, it only considers the precipitation factor and neglects effects of evaporation stemming from temperature and other meteorological factors (Wang et al., 2014, 2017; Yang et al., 2018). To solve the above problems, the standardized precipitation evapotranspiration index (SPEI), which considers the advantages of both PDSI and SPI, was developed to monitor and assess droughts (Vicente-Serrano et al., 2010). It not only accounts for the effect of evaporation on drought, but also has the capability of spatial comparability and characterization of different drought types with multiple timescales (Feng and Su, 2019; Wang et al., 2015). SPEI can be used to delineate spatial–temporal evolution of drought, drought characteristics, and impacts of drought at the regional and global scales (Mallya et al., 2016; Wang et al., 2014).
However, the commonly used SPEI fails to identify droughts with a duration less than 1 month (Van der Schrier et al., 2011; Vicente-Serrano et al., 2010). With future climate change, flash droughts have been recently categorized as a type of extreme climate event. Flash droughts occur along with sudden onset, rapid aggravation, and sudden end of drought and could lead to severe consequences (Pendergrass et al., 2020). It is imperative to characterize flash droughts with short-term duration (e.g., several days). To use the sub-month-resolution drought index, we have developed the daily SPEI for the first time, and our daily SPEI has been used to assess the drought and its impacts in previous studies (Wang et al., 2015, 2017). The new SPEI can not only identify drought with a duration of 1 month and more than 1 month, but also monitor drought with a duration of several days. In addition, our new daily SPEI has filled the gap in the capability to monitor the onset and duration of droughts. Our daily SPEI has similar principles as the commonly used monthly SPEI in terms of time accumulation effects (Vicente-Serrano et al., 2010; Wang et al., 2015; Yu et al., 2014). The daily SPEI data with different timescales can also meet the requirements of characterizing and assessing different drought types (meteorological drought, agricultural drought, and hydrological drought) at multiple timescales (Wang et al., 2014).
The SPEI can be calculated by the difference between daily precipitation and daily potential evapotranspiration (PET) (Vicente-Serrano et al., 2012). Precipitation can generally be directly obtained by the meteorological observation stations (Wang et al., 2015). But PET can only be estimated by drivers of meteorological data or remote sensing data (Wang et al., 2018, 2017). Although there are at least 50 methods to calculate the PET, the methods estimate the inconsistent and different values due to diverse assumptions, data inputs, and climatic regions (Grismer et al., 2002; Lu et al., 2005). PET plays an important role in understanding fluxes of the heat and mass of the atmospheric system at local and global scales (Thomas, 2000). Thus, it is necessary to choose a suitable method to estimate PET. The choice of candidate probability distributions for SPEI calculation is also very important (Vicente-Serrano et al., 2010, 2012). The chosen distribution for SPEI generally needs a location parameter because climatic water balance may have negative values (when PET is greater than precipitation in certain periods) (Wang et al., 2015, 2017). Distributions for SPEI normalization have a generalized logistic distribution, Pearson type III distribution, normal distribution, and generalized extreme value (GEV) distribution (Stagge et al., 2015). The four candidate SPEI distributions have the best goodness of fit for the accumulated climatic water balance (Stagge et al., 2015; Wang et al., 2015; Wang et al., 2017). However, the GEV distribution has the best performance among all four probability distributions across all of continental Europe, because of the lower rejection frequencies of GEV using several tests (Kolmogorov–Smirnov (K–S), Anderson–Darling (A–D), and Shapiro–Wilk (S–W)) (Stagge et al., 2015). Therefore, we choose the GEV distribution fitting the accumulated climatic water balance to calculate SPEI. The SPEI is suited to investigate the effects of climate change and global warming on drought severity. SPEI has been widely used in diverse studies on drought variability and impact and drought monitoring systems (Boroneant et al., 2011; Fuchs et al., 2012; Potop et al., 2014; Sohn et al., 2013).
The aim of this study, therefore, is to produce a long (1961–2018) daily drought index dataset for all of mainland China. Specifically, we used the new daily SPEI algorithm to produce the multi-timescale drought dataset at a daily time resolution. Meteorological data with 427 stations including multiple factors (daily precipitation, daily average air temperature, daily minimum air temperature, daily maximum air temperature, and sunshine) are used. The developed drought dataset at the national scale has the potential to be used to monitor and assess droughts and their impacts for the sectors including the agricultural sector, forest sector, hydrological sector, ecological sector, environmental sector, and so on.
Daily meteorological data from 1960 to 2018 were collected from the National
Meteorological Science Data Sharing Service Platform (
The location of meteorological stations across mainland China.
The daily SPEI can be calculated by the difference between daily
precipitation and daily potential evapotranspiration. Because air
temperature and solar radiation explained at least 80 % of
evapotranspiration variability (Martí et al., 2015; Priestley and
Taylor, 1972), the Hargreaves model based on temperature and solar radiation
can be used to estimate the daily potential evapotranspiration
(Hargreaves and Samani, 1982; Mendicino and Senatore, 2013; Wang et al.,
2015). The daily potential evapotranspiration can be obtained by the
following formula:
SPEI calculation depends on the accumulating deficit or surplus (
The daily SPEI dataset was calculated in five accumulating periods (30, 90, 180, 360, 720 d) based on the water balance (difference between precipitation and PET). The classifications for the SPEI drought classes are presented in Table 1.
Categorization of drought and wet grade according to the SPEI (Wang et al., 2014).
We used the method described by Yevjevich (1967) to define the
drought characteristics (severity, duration, and intensity). A drought event
can first be determined by drought start and end dates, and its duration
and severity were then assigned. Thus, we accounted for the continuity of
drought propagation. The continuous days with SPEI values less than the
threshold (such as
Schematic diagram of drought and wet events (the red shaded area denotes the drought events; the blue shaded area denotes the wet events).
The SPEI data based on 90 d (3-month) timescales can be used to identify soil moisture or agriculture droughts (Wang et al., 2014). Due to its important applications, we selected the SPEI data with the 90 d timescales as the example data for analysis in the present study. To investigate the spatial–temporal characteristics of the example data, we defined three variables including annual total drought severity (ATDS), annual total drought duration (ATDD), and annual total drought frequency (ATDF). The three variables were obtained by summing the severity, duration, and frequency of all the drought events in each year at 427 stations.
We also used the non-parametric Mann–Kendall (MK) test to detect monotonic
trends (Kendall, 1948; Mann, 1945), because the MK test does not
require data normality (Mann, 1945; Wang et al., 2020a, b). We computed slopes for ATDS, ATDD, and ATDF using the Sen method
(Sen, 1968). These statistical methods are commonly used in analyses of
water resources, climate, and ecology data. For the MK test, the global
trend for the entire series is significant when
The ATDS can be used to identify hot spots with more severe drought conditions. Figure 3 shows the calculated ATDS values across mainland China. We categorized ATDS values into two main groups, with higher ATDS values indicating more severe drought conditions. The distribution of ATDS values shows that, in general, the northeastern parts of China had more severe drought conditions than the southern parts. However, our results also indicate that the humid climate zone in the south also experienced severe drought conditions, though not as much as for northern parts of China (Fig. 3).
The spatial distribution of ATDS across mainland China.
Figure 4 shows that ATDD values ranged from 100 to 110 d for most stations across mainland China. This indicates that there was nearly one-third of a year when most stations were experiencing drought conditions. More stations with ATDD values ranging from 100 to 110 were found compared with stations with ATDD values of 120–130 (Fig. 4). For drought years, the duration days of drought events are expected to be longer. The ATDD had similar spatial distribution characteristics to those of the ATDS, indicating that droughts also occurred in the humid climate zone.
The spatial distribution of ATDD across mainland China.
Figure 5 shows the spatial distribution of ATDF values across mainland China. In general, most stations had four to six annual drought events. There were fewer stations with six to eight annual drought events compared with stations with two to four annual drought events. We also detected that drought events could occur in both arid and humid regions based on spatial distributions of ATDF values (Fig. 5). Since the ATDF indicated only the annual average drought events, we could expect that for the more severe drought years the ATDF would have greater values for different stations.
The spatial distribution of ATDF across mainland China.
The changing trends of ATDS can be used to detect whether drought severity
is weakening or intensifying with time. Figure 6 shows the spatial
distribution of changing trends of ATDS from 1961 to 2018 across
mainland China. In general, there were more stations with weakening trends
in drought severity than those with intensifying trends across all stations
(Fig. 6). It seems that both weakening and intensifying absolute values
were largest in the northeast, northwest, and central China compared with
other parts. However, after scrutiny, we found that drought severity tended
to weaken in the northeast, northwest, and central China, with more stations
having significant weakening trends according to statistical tests (
The spatial distribution of the changing trends of ATDS (the red
and green triangles indicate increasing and decreasing trends,
respectively. “***” denotes
The changing trends of ATDD can be used to detect whether drought duration
is getting shorter or longer. Figure 7 shows the spatial distribution of
changing trends for the ATDD across all stations. In general, stations in
the southeast demonstrated downward trends with shortening drought duration,
while stations in the northwest had upward trends for the ATDD with
increasing drought duration (Fig. 7). Note that the increasing or
decreasing trends for ATDD were significant (
The spatial distribution of the changing trends of ATDD (the red
and green triangles indicate increasing and decreasing trends,
respectively. “***” denotes
The changing trends of ATDF can be used to detect whether the frequency of
drought events is increasing or decreasing with time. Figure 8 shows the
spatial distribution of changing trends of ATDF across all stations. Most
stations demonstrated no significant trend in the frequency of drought
events, except for dozens of stations in western China having significant
upward trends (
The spatial distribution of the changing trends of ATDF (the red
and green triangles indicate increasing and decreasing trends,
respectively. “***” denotes
The reason for selecting a 90 d (3-month) timescale to assess spatial and temporal characteristics of drought conditions across mainland China is because the SPEI with the 90 d (3-month) timescale can indicate the agricultural drought (or soil moisture) (Van der Schrier et al., 2011; Wang et al., 2014, 2017), and its results are comparable with the PDSI (Dai et al., 2004; Van der Schrier et al., 2011) and other drought indices including the surface water supply index (SWSI) and moisture adequacy index (MAI) (Doesken et al., 1991; McGUIRE and Palmer, 1957). The commonly used monthly SPEI has been used to assess drought characteristics and their impacts worldwide from the regional scale to the global scale (Stagge et al., 2015; Vicente-Serrano et al., 2010; Wang et al., 2014). The SPEI with different timescales is relevant for meteorological drought (1-month timescale), agricultural drought (3–6-month timescale, about 90–180 d), hydrological drought (12-month timescale, about 360 d), and socioeconomic drought (24-month timescale, about 720 d), respectively (Homdee et al., 2016; Potop et al., 2014; Tirivarombo et al., 2018; Vicente-Serrano et al., 2010).
Our new SPEI dataset with multiple timescales was developed and compiled
using the daily SPEI algorithm in the previous study (Wang et al.,
2015). The daily SPEI has been used in drought characterization and
assessment and was validated by drought characterization and assessment
(Jevšenak, 2019; Jia et al., 2018; Salvador et al., 2019; Wang et
al., 2015, 2017). The global SPEI database with monthly
temporal resolution and 0.5
Although the daily SPEI has better performance in drought characterization and assessment (Jevšenak, 2019; Wang et al., 2017), the uncertainty of daily SPEI still needs to be evaluated in future works. Our daily SPEI dataset used the simple Hargreaves model based on temperature and solar radiation to estimate daily potential evapotranspiration (Hargreaves and Samani, 1982; Wang et al., 2017). We will further investigate effects of various evapotranspiration models (such as the CRAE model, Penman algorithm, Thornthwaite algorithm, Makkink algorithm, and Priestley–Taylor algorithm) on the calculation of SPEI (Makkink, 1957; Morton, 1983; Penman, 1948; Priestley and Taylor, 1972; Thornthwaite, 1944). We only chose SPEI based on the 90 d (3-month) timescale as an example to analyze drought characteristics, and the results demonstrated that there were no obvious intensifying trends for drought across mainland China, which is consistent with other studies (Han et al., 2020). Meanwhile, our newly developed daily SPEI will be further validated in other regions of the world. In addition, SPEI values at different timescales can be used as a proxy for other types of droughts, but it lacks the complete picture (no soil moisture condition, streamflow, etc.) (Zargar et al., 2011).
Our long-term daily SPEI dataset has contributed significantly to our understanding of drought evolution, especially flash drought. The dataset can be used to monitor and assess different drought types (meteorological drought, agricultural drought, and hydrological drought) through different timescale data. It can also identify the start and end dates for drought. The dataset is valuable to meteorological research and natural hazard communities for various purposes such as assessment of extreme climate or drought effect evaluation.
All daily SPEI datasets including data and their description at 427 observed
meteorological stations are provided open access via
figshare (Wang et al., 2020c), available at
In the present study, we have produced a daily SPEI dataset from 1960 to 2018 at 427 meteorological stations across mainland China. Our open-access dataset is an important contribution to drought assessment, and it can overcome the disadvantages of the commonly used monthly SPEI database. Our daily dataset can help monitor and assess the spatial and temporal characteristics of droughts. It can be used to assess the impacts of droughts on ecological systems, hydrological processes, and other natural resources. Our daily multi-timescale SPEI dataset can be widely used in studies on meteorological drought (1-month timescale), agricultural drought (3–6-month timescale), hydrological drought (12-month timescale), and socioeconomic drought (24-month timescale). The dataset will reduce the time spent on research and avoid the duplication of efforts, which will be highly attractive to meteorological, geographical, and natural hazard researchers and researchers from other areas.
QW led the study, developed the method, and wrote the manuscript with input from all the authors. JQ and XZ discussed the results and revised the manuscript. All the authors contributed to the final paper. QW, JY, RZ, XW, and XZ collected and analyzed data over time, providing statistics and material (figures and tables) for the paper.
The authors declare that they have no conflict of interest.
This research received financial support from the National
Natural Science Foundation of China (41601562), the Strategic Priority
Research Program of the Chinese Academy of Sciences (XDA13020506), and the China
Scholarship Council (201806655014). The authors sincerely thank James Howard Stagge for his
help on the codes and calculation of SPEI. Special thanks go to the
meteorological data provider from the China Meteorological Administration
(
This research received financial support from the National Natural Science Foundation of China (grant no. 41601562), the Strategic Priority Research Program of the Chinese Academy of Sciences (grant no. XDA13020506), and the China Scholarship Council (grant no. 201806655014).
This paper was edited by Ge Peng and reviewed by five anonymous referees.