Next Article in Journal
A Numerical Study on Travel Time Based Hydraulic Tomography Using the SIRT Algorithm with Cimmino Iteration
Previous Article in Journal
Impact of Unsteady Flow Events on Bedload Transport: A Review of Laboratory Experiments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatio-Temporal Analysis of Drought Indicated by SPEI over Northeastern China

1
School of Hydraulic Engineering, Dalian University of Technology, Dalian 116024, China
2
China Yangtze Power Company Limited, Yichang 443133, China
3
State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116024, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(5), 908; https://doi.org/10.3390/w11050908
Submission received: 13 March 2019 / Revised: 25 April 2019 / Accepted: 26 April 2019 / Published: 30 April 2019
(This article belongs to the Section Hydrology)

Abstract

:
Drought is a natural extreme climate event which occurs in most parts of the world. Northeastern China is one of the major agricultural production areas in China and also a typical vulnerable climate zone. To understand the spatio-temporal characteristics of drought over northeastern China, we first assessed the trends of precipitation and temperature. Drought events were then characterized by Standardized Precipitation Evapotranspiration Index over various temporal scales. The Trend Free Prewhitening Mann–Kendall test and distinct empirical orthogonal function, were used to investigate the trends and spatio-temporal patterns of droughts. The results indicate precipitation increasing trends are mostly detected in Heilongjiang and Jinling provinces, however, the majority of the trends are insignificant. Temperature increasing trends are detected over the entire northeastern China and most of them are significant. Decreasing drought trends are observed in Heilongjiang province and some bordering area in Jilin province, whereas increasing trends are noticed in Liaoning province and some bordering area in Jilin province. Two main sub-regions of drought variability—the Liaohe River Plain and the Second Songhua River basin (LS region), and the Songnen Plain and the Lesser Hinggan Mountains (SL region) are identified, and the detected droughts for the two sub-regions correspond well with recorded drought loss. The results will be beneficial for regional water resource management and planning, agriculture production, and ecosystem protection in northeastern China.

1. Introduction

In response to the exacerbating climate change, particularly the successive and rapid global warming [1], a wealth of regional studies suggest an increasing trend in the intensity and frequency of droughts in many parts of the world [2,3,4,5,6]. Specifically, drought is a naturally recurring extreme climate event over land characterized by below-normal precipitation over a period of months to a few decades [7]. It occurs with spatio-temporal variation in frequency, severity, duration, and can have devastating impacts on regional water resources, agriculture, industry, and other social-ecological systems, with profound economic and societal impacts in an increasingly globalized and uncertain world [8,9]. In contrast to the permanent aridity in arid areas, drought is a temporary dry period occurs in most parts of the world, even in wet and humid regions.
Due to the wide range of sectors affected by droughts, their varied spatio-temporal variations, the interdependence across climatic, hydrologic, geologic, geomorphic, ecological, and societal variables, and the demand placed on water supply by different users, it is very difficult to adopt a definition that fully describes the drought phenomena and the respective impacts [10,11,12]. Generally, based on the degree of water deficit, droughts are often classified into four types in practice: meteorological (lack of precipitation for a period of time), agricultural (lack of soil moisture for plant growth for a period of time without the effect of surface water), hydrological (lack of surface water for specific water consumption from a water resource management system), and socioeconomic (system failure in supplying water demand for a good water economy) [13,14,15]. Among all types of droughts, meteorological drought normally triggers the other three drought types. A nice summary of drought classification and definition can be found in Reference [16].
Our particular interest is the spatio-temporal analysis of drought over northeastern China, a region known as China’s granary, which is situated on one of the world’s few fertile black earth belts and has the highest endowment of cropland per capita in China. Being a typical vulnerable climate zone in middle- and high-latitude, droughts in northeastern China are affected by both global climate change and site characteristics. The climate in northeastern China has recently become warmer and drier, evidenced by the increasing frequency and intensity of droughts since the 1990s, which has become a serious threat to food production and economic development in the region [17]. Yang et al. [18] analyzed the affected area and frequency of drought during 1950–2010 and found the increasing frequency and intensity of droughts in northeastern China. Shen et al. [19] used the Standardized Precipitation Evapotranspiration Index (SPEI) to analyze droughts in northeastern China. Ma et al. [20] analyzed the projections of the future drought under a moderate scenario for greenhouse gas emissions (SRES A1B) and concluded that the crop growing season in northeastern China will continue to be drier during 2011–2100. Liu et al. [21] investigated the spatial and temporal characteristics of precipitation in Jilin province and found that almost 90% of the stations showed decreasing trends in annual precipitation and evaporation. Yu et al. [22] analyzed drought risk in agricultural production season based on Standardized Precipitation Index (SPI) index, which showed that drought occurred more frequently and severely in 1996–2009 than in 1965–1983, especially in the western and northern regions where precipitation is scarce over northeastern China. Yang et al. [23] selected SPI and agricultural drought indicators to explore the causes of drought, which showed that spring drought was the most serious problem in northeastern China, while drought in winter chiefly occurred in Heilongjiang province and the northwest of Liaoning province. Yue et al. [24] analyzed the relationship between the reconnaissance drought index (RDI) and main large-scale atmospheric circulation patterns, which showed that droughts in northeastern China were less frequent, but had greater severity and duration, mainly in the western and southeastern part.
Although many previous studies have investigated the duration, tendency, and frequency of droughts over northeastern China, their work mostly focuses on the temporal characteristics of droughts; even for few studies on spatial patterns of droughts, the results are quite contradicted and different. In addition, considering that agriculture is vital in northeastern China, drought analysis during crop growing season should be highlighted, which has not been the usual case used in previous studies, where monthly and annual were often studied. The drought characteristics in northeastern China remains uncertain and doubtful, which poses considerable challenges to water security and sustainable development, especially under changing climatic conditions.
Our primary objective is to provide a comprehensive spatio-temporal analysis of droughts over northeastern China using the most up-to-date data. Drought events are characterized by SPEI over various temporal scales including 1, 3, 6, and 12 consecutive months, and 5-month together from May to September, especially for crop growing season. First, the Trend Free Prewhitening Mann–Kendall (TFPW-MK) test are used to investigate the trends in the SPEI series. Major patterns of long-term change of drought events are then characterized using the Distinct Empirical Orthogonal Function (DEOF). Finally, the detected droughts are compared with the recorded historical drought in identified sub-regions.
To achieve this objective, the paper is structured as follows: Section 2 introduces the study area and the data used. Section 3 describes the methodology for analyzing drought characteristics. In Section 4, the results and discussion are revealed. Section 5 provides the summary and conclusions of the study.

2. Study Area and Data

Northeastern China (38°41′ N–53°33′ N, 118°50′ E–135°4′ E) encompasses the three provinces of Liaoning, Jilin, and Heilongjiang, with a total area of 790,000 km2 and a population of about 110 million. The northeast is China’s most important source of soybeans, corn, and japonica rice, accounting for 41% of the country’s soybean production, 34% of its corn output, and 30–50% percent of its japonica rice crop. The landscape is mostly characterized by plains and mountains with elevations from 3 to 1017 m. The high latitude makes northeast the coolest region in China with long durations of severely cold winters. The annual precipitation decreases from southeast (900 mm) to northwest (400 mm) as the distance from the sea increases, and the annual temperature varies from −4.2 to 10.9 °C, decreasing from south to north. The precipitation is generally concentrated in the summer and autumn, and it coincides with the crop growing season, in which the main crops are sowed at the beginning of May and harvested at the end of September.
The monthly ground-based meteorological datasets were obtained from China Meteorological Administration (CMA; http://www.cma.gov.cn/), including precipitation, average temperature, maximum and minimum temperature, relative humidity, sunshine hours, and wind speed. To analyze the most up-to-date drought characteristics over northeast China, the selected period ranges from 1960 to the newest 2017 altogether 696 months, before 1960 most of the meteorological stations had changed their observational instruments. There are 91 meteorological stations in northeastern China, however, 8 of 91 meteorological stations have insufficient record length, and 83 meteorological stations were finally selected across this region, with 22 stations in Liaoning province, 31 in Jilin province, and 30 in Heilongjiang province. Northeastern China and meteorological station locations are shown in Figure 1. Meteorological data from the CMA has been carefully assessed to ensure consistency, continuity, and quality by taking into account instrument location and relocation of the station when necessary, to mitigate the effects of changes in vegetation and urbanization [25]. It should be noted that there are only two meteorological stations in the north-west part of Heilongjiang province, therefore, the results and analyses regarding this area can be doubtful as interpolated values bring huge uncertainties.

3. Methodology

3.1. Standardized Precipitation Evapotranspiration Index (SPEI)

Drought indices are quantitative measures of drought level by assimilating single or multiple atmosphere variables (e.g., precipitation and temperature) into a single numerical value, which prove to be more useful than the original variables [26]. Based on the analysis of the advantages and disadvantages of Palmer Drought Severity Index (PDSI) [27] and Standardized Precipitation Index (SPI) [28], Vicente-Serrano et al. [29] proposed a new drought index, Standardized Precipitation Evapotranspiration Index (SPEI). This index takes both precipitation and temperature into account, it combines the response of drought to evapotranspiration which PDSI pays attention to and the advantages of simple calculation and multi-time scale of SPI. SPEI has been successfully used and validated to analyze the spatial and temporal characteristics of drought in many parts of China [30,31,32].
The SPEI calculation in this paper adopted the SPEIcalc, an SPEI calculator from Consejo Superior de Investigaciones Científicas (CSIC; http://spei.csic.es/index.html) [29]. SPEI is very easy to calculate, which is based on the monthly difference between precipitation and potential evapotranspiration (PET), at different time scales of interest. The SPEI at different time scales can represent different climatic water balances. Vicente-Serrano et al. [29] applied the Thornthwaite procedure [33] to estimate the PET, which only requires monthly mean temperature and latitude of the site. With a value for PET, the difference between the precipitation P and PET for the month i can be calculated as follows:
D i = P i P E T i
The calculated Di values are aggregated at different time scales. The difference in a given month j and year i depends on the chosen timescale k. For example, the accumulated difference for one month in a particular year i with a 12-month timescale can be calculated as follows [29]:
X i , j k = { l = 13 k + j 12 D i 1 , l + l = 1 j D i , l ,   j < k     l = j k j D i , l ,    j k
where Di,l is the PPET difference in the l month of year i.
Based on the behavior at the most extreme values, the log-logistic distribution adapts very well to standardize the D series to obtain the SPEI [29]. The probability density function of a three-parameter log-logistic distributed variable is expressed as:
f ( x ) = β α ( x γ α ) β 1 [ 1 + ( x γ α ) β ] 2
where α, β, and γ are scale, shape, and origin parameters, respectively, for D values in the range (γ > D < ∞).
According to Ahmad et al. [34], the L-moment procedure is the most robust and easy approach to estimate the parameters of the Log-logistic distribution. When L-moments are calculated, the parameters of the Pearson III distribution can be obtained following Singh et al. [35]:
β = 2 w 1 w 0 6 w 1 w 0 6 w 2
α = ( w 0 2 w 1 ) β Γ ( 1 + 1 / β ) Γ ( 1 1 / β )
γ = w 0 α Γ ( 1 + 1 β ) Γ ( 1 1 β )
where the k-th probability-weighted moment can be estimated as
w k = 1 n i = 1 n x i ( 1 i 0.35 n ) k , k = 0 , 1 , 2 ,
where xi is an ordered random sample x1x2xn, and n is the sample size.
The probability distribution function of the D series, according to the log-logistic distribution, is given by
F ( x ) = [ 1 + ( α x γ ) β ] 1
The SPEI can be easily obtained as the standardized values of F(x) [36].
SPEI = W C 0 + C 1 W + C 2 W 2 1 + d 1 W + d 2 W 2 + d 3 W 3
where
W = 2 ln ( P ) , P 0.5
and P is the probability of exceeding a determined D value, P = 1 − F(x). If P > 0.5, then P is replaced by 1 − P and the sign of the resultant SPEI is reversed. The constants are: C0 = 2.515517, C1 = 0.802853, C2 = 0.010328, d1 = 1.432788, d2 = 0.189269, and d3 = 0.001308.
The SPEI of one month usually represents meteorological drought, while the timescale of 3–6 months is considered as agricultural drought index. Longer scales, such as 6 months and 12 months, are used to indicate hydrological drought and to monitor surface water resources [37,38]. In addition to the consecutive 1, 3, 6, and 12 consecutive months, and 5-month together from May to September is also selected to analyze the drought characteristic for crop growing season. The SPEI class is shown in Table 1.

3.2. Trend Free Prewhitening Mann–Kendall Nonparametric Test

The Mann–Kendall (MK) nonparametric test, proposed by Mann [39] and then modified by Kendall [40], is based on the correlation between the ranks of a time series and their chronological order. This method is often used to analyze the trend change of hydro-meteorological time series such as precipitation, runoff, temperature, drought index, etc. [41,42,43,44]. The test has the advantage of being rarely disturbed by outliers, it is easy to calculate, and it does not assume any special form for the probability distribution function of the data [45].
However, some hydrological time series may usually have serial correlation, which will increase the probability that the MK test detects a significant trend, leading to a disproportionate rejection of the null hypothesis of no trend, whereas the null hypothesis is actually true [46]. On this basis, Serinaldi et al. [47] compared the results of standard and improved versions of the MK test, and emphasized the effect of serial correlation on trend detection. Clarke [48] also pointed out the problem of circular reasoning in hydrological statistics and the failure to account for spatial correlations between variables. Poppick et al. [49] took the temperature trend as an example to compare parametric statistical methods and nonparametric methods. Therefore, Trend Free Prewhitening MK (TFPW-MK) [46], which detrends and prewhitens the time series before MK test, is adopted in this paper. The specific steps are as follows:
Step 1. The slope b of a trend in time series X t is estimated by the Theil-Sen approach (TSA) [50]. If the slope is almost equal to zero, then it is unnecessary to continue to conduct trend analysis. If it differs from zero, then it is assumed to be linear, and the time series are detrended by:
X t = X t T t = X t b t
Step 2. The lag-1 serial correlation coefficient r1 of the detrended series X t is computed using Equation (11) and then the Auto Regression Model (1) is removed from the X t by Equation (12).
r k = 1 n k t = 1 n k [ X t E ( X t ) ] [ X t + k E ( X t ) ] 1 n t = 1 n [ X t E ( X t ) ] 2
Y t = X t r 1 X t 1
This pre-whitening procedure after detrending the series is referred to as the TFPW procedure. The residual series after applying the TFPW procedure should be an independent series.
Step 3. The identified trend T t and the residual Y t are blended by
Y t = Y t + T t
Step 4. The MK test is applied to the blended series to assess the significance of the trend.
We can get the statistic Z through the calculation of TFPW-MK test. The time series shows an increasing trend when Z is greater than 0, while the trend is decreasing when Z is less than 0, and the greater the absolute value of Z, the more obvious the trend of increase or decrease. Thus, in the bilateral trend test, if |Z| ≥ Z1−α/2 at a given α confidence level, the original hypothesis is unacceptable. That is to say, at the 1 − α confidence level, the time series has a distinct upward or downward trend. The Z values of ±1.96, and ±1.64 are equal to the confidence levels of 0.05, and 0.1, respectively. Based on these confidence levels, the detected trend could be classified into six zones according to Z value [42]: (1) Z < −1.96, indicating significant decreasing trend; (2) Z ∈ [−1.96 to −1.64), indicating weak decreasing trend; (3) Z ∈ [−1.654 to 0), indicating no significant decreasing trend; (4) Z ∈ (0 to 1.64], indicating no significant increasing trend; (5) Z ∈ (1.64 to 1.96], indicating weak increasing trend; and (6) Z > 1.96, indicating the significant increasing trend.

3.3. Distinct Empirical Orthogonal Function

The analysis presented in this paper adopt the Distinct Empirical Orthogonal Function (DEOF) [51], which has proved to obtain better results in studies of spatial variability patterns identification. The DEOF method is based on normal EOF method, and is then compared against a fitted isotropic diffusion process [52]. In DEOF, the EOF-modes are considered as a representation of a continuous spectrum of spatial patterns resulting from a stochastic process. In the spatial structure, some spatial structures will be more dominant than others. The most distinguished EOF-modes (DEOFs) from the isotropic diffusion null hypothesis can be found by rotating the leading EOF-modes. The time series of the DEOFs are referred to as the distinguished principle components (DPCs). The details about the method can be found in http://users.monash.edu.au/~dietmard/deof-analysis.html.

4. Results and Discussion

4.1. Variation Characteristics of Precipitation and Temperature

To better understand the drought patterns over northeastern China, the variation characteristics of the two elements in SPEI, precipitation, temperature are analyzed first. TFPW-MK trend test was applied to the 83 annual precipitation and temperature series from 1960 to 2017 at each meteorological station, and then, the Z values were interpolated over northeast China, as shown in Figure 2. It can be seen that more than half of northeastern China has increasing trends in precipitation and the regions with increasing trends are mostly located in Heilongjiang and Jilin provinces. However, most increasing trends are insignificant. Most of the Liaoning province and regions around the north of Lesser Hinggan Mountains in Heilongjiang province present decreasing trends; however, the decreasing trends are insignificant, either.
For temperature, it is more straightforward that the entire northeastern China shows increasing trends and the majority of them are significant. The strong increasing trends of temperature is consistent with [8], which reports each of the past three decades has been warmer than all preceding decades since 1850 and the period from 1983–2012 was “likely” the warmest 30-year period of the past 1400 years in the Northern Hemisphere.
The areal average precipitation and temperature anomaly curves relative to the average of entire study period at annual and crop growing seasonal scales are presented in Figure 3 and Figure 4. Northeastern China experienced dry periods in two separate decades, at both annual and crop growing seasonal scales, during the 1970s and years after 2000. The annual scale deviations are −23.9 mm and −5.3 mm for the two periods, respectively. The dry period presence is much stronger for the crop growing season, where the seasonal deviations are −49.0 mm and −42.5 mm, respectively. These analyses indicate that the amplitude of variation for precipitation at crop growing season is larger, which means even a mild water shortage year could still have a severe water shortage during crop growing shortage, consistent with Yu et al. [22].
It is apparent from Figure 4 that the warming trends are quite obvious at both annual and crop growing seasonal scales, with almost strictly monotone increasing at decade level. The climate (or temperature) in northeastern China can be divided into a warm and a cold period at 1990. For annual temperature, the increases for the 1990s and the years after 2000 are almost identical; while the increase in years after 2000 is stronger than in the 1990s for temperature during crop growing season.

4.2. Results and Uncertainty of SPEI

SPEI was calculated for 83 meteorological stations to analyze the regional conditions of drought at various time scales. The areal average SPEI series of northeastern China at different time scales are shown in Figure 5. At areal average, SPEI1 and SPEI3 have increasing trend and SPEI6May–Sep has decreasing trend. While for SPEI6 and SPEI12, the changing trends are not very obvious.
Since all the analyses in this paper are based on the SPEI results, it is very necessary to study the uncertainty regarding SPEI calculation. The accuracy and reliability of SPEI calculation basically depend on the observed data, especially on the length of the sample. Similar to Hu et al. [53] about SPI sampling uncertainty, the uncertainty of SPEI caused by sampling uncertainty is calculated as follows: (1) Resampling technology is used to extract a certain number of samples from the original sample (the precipitation and temperature data), which allows repeated sampling. (2) Calculate the SPEI according to the extracted samples. (3) Repeat the above 1000 times, and thus get 1000 samples of SPEI. (4) Calculate the sample variance of the above 1000 samples of SPEI, and get the variance of the statistic [53]. The expected estimation of SPEI is the average value of the 1000 sets of SPEIj, j = 1, 2, …, 999, 1000, and the lower bound (LB = SPEI51) and upper bound (UB = SPEI950) of the 90% confidence interval (CI) estimation are calculated, accordingly.
Taking SPEI6May–Sep as an example, Figure 6 shows the results of uncertainty of SPEI6May–Sep due to sampling in three typical meteorological stations. It can be seen that the 90% CI is very narrow compared to the expected value, which means sampling uncertainty is small in our paper. Therefore, the SPEI results are directly used for further analyses.

4.3. Spatial Patterns of Temporal Trends in Drought

TFPW-MK trend test was then applied to the SPEI series at each meteorological station and the Z values were interpolated over northeastern China. Figure 7 depicts the spatial structure of long-term SPEI trends. Decreasing drought trends (i.e., increase of SPEI) are observed in Heilongjiang province and some bordering area in Jilin province, whereas increasing trends are noticed in Liaoning province and some bordering area in Jilin province. The situation of long-term drought trends is quite different at different timescales. As can be seen from Figure 7, the decreasing drought trends are observed in majority of northeastern China at 1- and 3-month timescales, while with the increasing of timescale, the increasing trends of drought have also been noticed in more and more regions of northeastern China. Specifically, the increasing trends of drought are dominant at 12-month timescale. Overall, the significant increase in drought severity mostly occurs in the coastal region of Liaodong Gulf, southeastern Liaodong Peninsula, and the downstream region of Hunhe River Basin, while the significant decrease is mostly observed in regions such as the northwest corner of Heilongjiang province, the southeast of Lesser Hinggan Mountains and the northeast of Changbai Mountains.
Although the changing trends are insignificant over northeastern China at crop growing seasonal scale, the vast majority of the area experiences insignificant or weak increasing trends, except the very northwest corner of Heilongjiang province, the northeast of Changbai Mountains and the central part of Songnen Plain have observed insignificant decreasing drought trends. This is not surprising given the fact that the situation of precipitation and temperature are more adverse during the crop growing season. From the analyses above, we could conclude that the drought situation may be even worse during crop growing season as the vast majority of northeastern China has an increasing drought trend.

4.4. Spatial Variability of Drought Using DEOF

Table 2 displays percentages of variance explained by the first two DEOFs relative to the considered timescales. The first two DEOFs are able to explain about 54.1 to 58.6% of the total explained temporal variances of drought depending on the timescale, with the maximum and minimum variances observed for SPEI1 and SPEI12, respectively.
The patterns of the leading two DEOFs for SPEI time series were computed for 83 meteorological stations and interpolated over northeastern China, shown in Figure 8. The congruence coefficient proposed in Richman and Lamb [54] was used to check whether the patterns represented potential physical modes of SPEI variability, or were only an artifact of the DEOF analysis. The congruence coefficient ranges from 0.81 to 0.91 within the regions enclosed by the contour in Figure 8, which indicates the identified variability modes are acceptable as described in in Richman and Lamb [54]. For all SPEI series, the loading patterns of the first DEOF (DEOF1) with a maximum explained the temporal variance of 45.2% and 40.2% are relative to SPEI1 and SPEI12, respectively. The DEOF1 values are all positive and the spatial distribution of the DEOF1 indicates the high positive loading in the Liaohe River Plain and the Second Songhua River basin (LS region). This shows that the LS region has a similar temporal distribution of drought as that detected under various time series. This region is the main channel for tropical cyclone and other southern cyclone systems to penetrate northeastern China. These cyclones usually make landfall from the west of Liaoning province, bringing abundant water vapour and thus making this region wet. Droughts often occur in the years without cyclone. Besides, LS region has many similar hilly and mountainous terrains, which is another reason for the similar drought characteristics in this region. The unique natural weather and land surface conditions are the main reasons for the similarity of drought characteristics in this region.
The DEOF2s explain a maximum variance of around 14.1% in the case of SPEI6May–Sep and 13.4% for SPEI1, respectively. The spatial distribution of DEOF2 indicates the high positive spatial behavior of drought in the Songnen Plain and the Lesser Hinggan Mountains (SL region). The SL region belongs to the cold temperate zone where the climate characteristics are significantly affected by the Lesser Hinggan Mountains, while Liaohe River Plain and area surrounding its east are near the ocean which mainly receive the influence of seawater vapor. In addition, this region is also affected by the air currents of the Sea of Japan and the Sea of Okhotsk. Although two DEOFs are unable to extract the full temporal variation of drought across the region, their loadings seem to divide the whole region into two sub-regions that is LS and SL regions, characterized by different drought variabilities.
The two drought sub-regions DEOF1 and DEOF2 coincide basically with the drought-prone Heilongjiang and northwest Liaoning described by Yu et al. [22] and Yang et al. [23], which also show the importance of dividing the drought sub-regions for research. Therefore, by studying the changes of drought in DEOF1 and DEOF2 over time, we can understand the main drought evolutions in northeastern China.
The DPC scores are also investigated to characterize the temporal variability of drought severity. Figure 9 presents DPC1 and DPC2 scores for the value of SPEI6May–Sep, since crop growing season is our primary concern. Figure 9a reveals that the LS region (DPC1) experienced severe drought in years 1982 and 1999 to 2001, which is consistent with Yue et al. [24] describing that drought hot spots mainly occur in the western part of the study area in 2000–2004. These detected droughts correspond well with the China’s Drought Statistics Yearbook, which recorded that northeastern China suffered from severe drought in 1982. Especially in the Liaohe River Plain and the Second Songhua River basin, most of the moisture was merely 5–10% for 50-cm deep surface soil. It also recorded that 3-year consecutive drought has occurred in Liaoning province during 1999–2001, with 90% rivers were witnessed cutoff. For the SL region (DPC2), the only year extreme drought was observed in 1995 with DPC score less than −2.5. In 1995, nine counties and cities in Heilongjiang province were recorded suffering from extreme drought, and more than half of the agricultural areas were in a drought state. The rest 1964, 2001, 2004, 2007, and 2010 years detected as severe dry were all recorded severe drought in the LS region.

5. Summary and Conclusions

This paper provides a comprehensive analysis of spatio-temporal patterns over northeastern China during the period 1960–2017. The widely used SPEI drought index is selected to evaluate drought characteristics at various timescales, especially the crop growing season from May to September. The results are based on long-term meteorological datasets and widely used methods, and are consistent with previous studies and observed drought statistics. Therefore, the results are credible to some degree and can be used to assist with drought mitigation. The main conclusions are summarized as follows.
(1)
The annual precipitation changing trends increase from south to north, with decreasing trends in nearly the entire Liaoning province, and increasing trends in Heilongjiang and Jilin provinces. However, the majority of the trends are insignificant. The warming trends are more obvious and straightforward, as increasing trends are detected over entire northeastern China.
(2)
TFPW-MK test detects the changing trends of SPEI at various timescales. Overall, significant increasing drought trends are observed in the coastal region of Liaodong Gulf, southeastern Liaodong Peninsula, and the downstream region of Hunhe River Basin, while significant decreasing trends are noticed in the northwest corner of Heilongjiang province, the southeast of Lesser Hinggan Mountains and the northeast of Changbai Mountains. Further, drought increasing trends are more dominant at crop growing seasonal scale, thus the drought situation may be even worse during crop growing season than the rest of the year.
(3)
DEOF is used to identify two main sub-regions of drought variability—the Liaohe River Plain and the Second Songhua River basin (LS region), and the Songnen Plain and the Lesser Hinggan Mountains (SL region). Based on the crop growing seasonal DPC scores, the LS region experienced severe droughts in the years 1982 and 1999 to 2001. In the SL region, severe to extreme droughts were observed in 1964, 1995, 2001, 2004, 2007, and 2010.
The detected drought trend in this paper is very important as it provides reference for the future study of drought prediction in northeastern China. Special attention should be paid to the obvious growth of temperature, which can produce positive interaction and feedback that intensify drought. The identification of these sub-regions with similar drought variability and characteristics can be useful for drought risk management at a regional scale in northeastern China. Especially for management departments, drought resistance strategy and irrigation schemes can be formulated purposefully according to drought sub-regions.
Like any natural phenomena, drought is very difficult to characterize. Although various drought indices have been proposed to monitor drought, this paper only adopted the widely used SPEI index because the present analyses are already very complicated with the temporal trends and spatial variability of drought over various temporal scales. The SPEI index is merely an index belongs to meteorological drought, so the comparable studies of multi drought indices in the future would provide more comprehensive analysis of drought. It also should be noted that there is only two meteorological station in the north-west part of Heilongjiang province, therefore, the results and analyses regarding this area can be very doubtful as interpolated values will bring huge uncertainties. If gridded data can be obtained for drought analysis, the reliability of the results will be greatly improved. Another factor contributing to the uncertainty is that, although our results are consistent with the measured data, there is no quantitative relationship between the degree of loss and SPEI.

Author Contributions

L.Y. designed the study and wrote the manuscript; K.S. performed the experiments and analyzed the data; H.Z. revised the manuscript according to reviewers’ comments; Z.X. assisted with the results and conclusions sections; J.H. assisted with the data collection and preprocessing; C.Z. provided a feedback on the structure of the manuscript and reviewed the manuscript.

Funding

This study was jointly supported by the Supported by the National Key Research and Development Program of China (No. 2017YFC0406001), and supported by the National Natural Science Foundation of China (No. 91647201, No.91747102, No. 51579027).

Acknowledgments

Special thanks are given to the two anonymous reviewers for their constructive remarks.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Field, C.B.; Barros, V.R.; Mastrandrea, M.D.; Mach, K.J.; Abdrabo, M.K.; Adger, N.; Anokhin, Y.A.; Anisimov, O.A.; Arent, D.J.; Barnett, J.; et al. Summary for Policymakers. In Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2014; pp. 1–32. [Google Scholar]
  2. Barriopedro, D.; Gouveia, C.M.; Trigo, R.M.; Wang, L. The 2009/10 Drought in China: Possible Causes and Impacts on Vegetation. J. Hydrometeorol. 2012, 13, 1251–1267. [Google Scholar] [CrossRef] [Green Version]
  3. Zhang, W.; Jin, F.F.; Turner, A. Increasing Autumn Drought over Southern China Associated with Enso Regime Shift. Geophys. Res. Lett. 2014, 41, 4020–4026. [Google Scholar] [CrossRef]
  4. Cai, Q.; Liu, Y.; Liu, H.; Ren, J. Reconstruction of Drought Variability in North China and Its Association with Sea Surface Temperature in the Joining Area of Asia and Indian–Pacific Ocean. Palaeogeogr. Palaeoclim. Palaeoecol. 2015, 417, 554–560. [Google Scholar] [CrossRef]
  5. Spinoni, J.; Naumann, G.; Vogt, J.; Barbosa, P. European Drought Climatologies and Trends Based on a Multi-Indicator Approach. Glob. Planet. Chang. 2015, 127, 50–57. [Google Scholar] [CrossRef]
  6. Wang, Z.; Li, J.; Lai, C.; Zeng, Z.; Zhong, R.; Chen, X.; Zhou, X.; Wang, M. Does Drought in China Show a Significant Decreasing Trend from 1961 to 2009? Sci. Total Environ. 2017, 579, 314–324. [Google Scholar] [CrossRef]
  7. Dai, A. Drought under Global Warming: A Review. Wiley Interdiscip. Rev. Clim. Chang. 2011, 2, 45–65. [Google Scholar] [CrossRef]
  8. Stocker, T.F.; Qin, D.; Plattner, G.K.; Tignor, M.M.B.; Allen, S.K.; Boschung, J.; Nauels, A.; Xia, Y.; Bex, V.; Midgley, P.M. Ipcc, 2013: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2013; p. 1535. [Google Scholar]
  9. Sternberg, T. Regional Drought Has a Global Impact. Nature 2011, 472, 169. [Google Scholar] [CrossRef]
  10. Asong, Z.E.; Wheater, H.S.; Bonsal, B.; Razavi, S.; Kurkute, S. Historical Drought Patterns over Canada and Their Teleconnections with Large-Scale Climate Signals. Hydrol. Earth Syst. Sci. 2018, 22, 3105–3124. [Google Scholar] [CrossRef]
  11. Panagoulia, D.; Dimou, G. Definitions and Effects of Droughts. In Proceedings of the Conference on Mediterranean Water Policy: Building on existing experience, Mediterranean Water Network, Valencia, Spain, April 1998. [Google Scholar]
  12. Awass, A.A. Hydrological Drought Analysis—Occurrence, Severity, Risks: The Case of Wabi Shebele River Basin. Ph.D. Thesis, Universität Siegen, Siegen, Germany, January 2009. [Google Scholar]
  13. Majid, M.; Amirataee, B. Comprehensive Stochastic Assessment of Meteorological Drought Indices. Int. J. Climatol. 2017, 37, 998–1013. [Google Scholar]
  14. Babak, A.; Montaseri, M. The Performance of Spi and Pnpi in Analyzing the Spatial and Temporal Trend of Dry and Wet Periods over Iran. Nat. Hazards 2017, 86, 89–106. [Google Scholar]
  15. Wilhite, D.A.; Glantz, M.H. Understanding: The Drought Phenomenon: The Role of Definitions. Water Int. 1985, 10, 111–120. [Google Scholar] [CrossRef]
  16. Mishra, A.K.; Singh, V.P. A Review of Drought Concepts. J. Hydrol. 2010, 391, 202–216. [Google Scholar] [CrossRef]
  17. Zhang, S.J.; Zhang, Y.S.; Chen, P.S.; Cai, F.; Jin-Wen, W.U.; Rui-Peng, J.I. Change characteristics of humid index and dry/wet climateboundary in Northeast China. Agric. Res. Arid Areas 2011, 29, 226–232. (In Chinese) [Google Scholar]
  18. Guiyu, Y.; Dongmei, H.; Yiming, C. Studies on drought and flood transformation in northeast China from 1950 to 2010. China Water Resour. 2014, 5, 45–48. (In Chinese) [Google Scholar]
  19. Guoqiang, S.; Haifeng, Z.; Zhenfeng, L. Spatiotemporal analysis of meteorological drought (1961—2014) in NortheastChina using a standardized precipitation evapotranspiration index. Acta Ecol. Sin. 2017, 37, 5882–5893. (In Chinese) [Google Scholar]
  20. Jiangyong, M.; Jie, P.; Yinlong, X.; Jiang, J. Drought Trend in Northeast China in the Future under SRES A1B Scenario. Arid Zone Res. 2013, 30, 329–335. (In Chinese) [Google Scholar]
  21. Liu, B.; Xiao, C.; Liang, X. Evaluation of Spatial and Temporal Characteristics of Precipitation Variations in Jilin Province, Northeast China. Theor. Appl. Climatol. 2015, 122, 129–142. [Google Scholar] [CrossRef]
  22. Yu, X.; He, X.; Zheng, H.; Guo, R.; Ren, Z.; Zhang, D.; Lin, J. Spatial and temporal analysis of drought risk during the crop-growing season over northeast China. Nat. Hazards 2014, 71, 275–289. [Google Scholar] [CrossRef]
  23. Yang, G.; Shao, W.; Wang, H.; Han, D. Drought Evolution Characteristics and Attribution Analysis in Northeast China. Procedia Eng. 2016, 154, 749–756. [Google Scholar] [CrossRef]
  24. Yue, Y.; Shen, S.H.; Wang, Q. Trend and Variability in Droughts in Northeast China Based on the Reconnaissance Drought Index. Water 2018, 10, 318. [Google Scholar] [CrossRef]
  25. Li, Q.; Dong, W.; Li, W.; Gao, X.; Jones, P.; Kennedy, J.; Parker, D. Assessment of the Uncertainties in Temperature Change in China during the Last Century. Chin. Sci. Bull. 2010, 55, 1974–1982. [Google Scholar] [CrossRef]
  26. Zargar, A.; Sadiq, R.; Naser, B.; Khan, F.I. A Review of Drought Indices. Environ. Rev. 2011, 19, 333–349. [Google Scholar] [CrossRef]
  27. Palmer, W.C. Meteorological Drought. In Research Paper No. 45; Us Dept. of Commerce: Washington, DC, USA, 1965; p. 58. [Google Scholar]
  28. McKee, T.B.; Doesken, N.J.; Kleist, J. The Relationship of Drought Frequency and Duration to Time Scales. In Proceedings of the 8th Conference on Applied Climatology, Anaheim, CA, USA, 17–22 January 1993. [Google Scholar]
  29. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A Multiscalar Drought Index Sensitive to Global Warming: The Standardized Precipitation Evapotranspiration Index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef] [Green Version]
  30. Yu, M.; Li, Q.; Hayes, M.J.; Svoboda, M.D.; Heim, R.R. Are Droughts Becoming More Frequent or Severe in China Based on the Standardized Precipitation Evapotranspiration Index: 1951–2010? Int. J. Climatol. 2014, 34, 545–558. [Google Scholar] [CrossRef]
  31. Xu, K.; Yang, D.; Yang, H.; Li, Z.; Qin, Y.; Shen, Y. Spatio-Temporal Variation of Drought in China During 1961–2012: A Climatic Perspective. J. Hydrol. 2015, 526, 253–264. [Google Scholar] [CrossRef]
  32. Zuo, D.; Cai, S.; Xu, Z.; Li, F.; Sun, W.; Yang, X.; Kan, G.; Liu, P. Spatiotemporal Patterns of Drought at Various Time Scales in Shandong Province of Eastern China. Theor. Appl. Climatol. 2018, 131, 271–284. [Google Scholar] [CrossRef]
  33. Thornthwaite, C.W. An Approach toward a Rational Classification of Climate. Geogr. Rev. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  34. Ahmad, M.I.; Sinclair, C.D.; Werritty, A. Log-logistic flood frequency analysis. J. Hydrol. 1998, 98, 205–224. [Google Scholar] [CrossRef]
  35. Singh, V.P.; Guo, H.; Yu, F.X. Parameter estimation for 3-parameter log-logistic distribution (LLD3) by Pome. Stoch. Hydrol. Hydraul. 1993, 7, 163–177. [Google Scholar] [CrossRef]
  36. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions; Dover Publications: New York, NY, USA, 1965. [Google Scholar]
  37. Hayes, M.; Svoboda, M.; Wall, N.; Widhalm, M. The Lincoln Declaration on Drought Indices: Universal Meteorological Drought Index Recommended. Bull. Am. Meteorol. Soc. 2011, 92, 485–488. [Google Scholar] [CrossRef] [Green Version]
  38. Beguería, S.; Vicente-Serrano, S.M.; Reig, F.; Latorre, B. Standardized Precipitation Evapotranspiration Index (Spei) Revisited: Parameter Fitting, Evapotranspiration Models, Tools, Datasets and Drought Monitoring. Int. J. Climatol. 2014, 34, 3001–3023. [Google Scholar] [CrossRef]
  39. Mann, H.B. Nonparametric Tests against Trend. Econom. J. Econom. Soc. 1945, 13, 245–259. [Google Scholar] [CrossRef]
  40. Kendall, M.G. Rank Correlation Methods; Charles Griffin: London, UK, 1955. [Google Scholar]
  41. Dawood, M. Spatio-Statistical Analysis of Temperature Fluctuation Using Mann–Kendall and Sen’s Slope Approach. Clim. Dyn. 2017, 48, 783–797. [Google Scholar]
  42. Wang, W.; Zhu, Y.; Xu, R.; Liu, J. Drought Severity Change in China During 1961–2012 Indicated by Spi and Spei. Nat. Hazards 2015, 75, 2437–2451. [Google Scholar] [CrossRef]
  43. Da Silva, R.M.; Santos, C.A.; Moreira, M.; Corte-Real, J.; Silva, V.C.; Medeiros, I.C. Rainfall and River Flow Trends Using Mann–Kendall and Sen’s Slope Estimator Statistical Tests in the Cobres River Basin. Nat. Hazards 2015, 77, 1205–1221. [Google Scholar] [CrossRef]
  44. Ay, M.; Kisi, O. Investigation of Trend Analysis of Monthly Total Precipitation by an Innovative Method. Theor. Appl. Climatol. 2015, 120, 617–629. [Google Scholar] [CrossRef]
  45. Serrano, A.; Mateos, V.L.; Garcia, J.A. Trend Analysis of Monthly Precipitation over the Iberian Peninsula for the Period 1921–1995. Phys. Chem. Earth Part B Hydrol. Oceans Atmos. 1999, 24, 85–90. [Google Scholar] [CrossRef]
  46. Yue, S.; Pilon, P.; Phinney, B.; Cavadias, G. The influence of autocorrelation on the ability to detect trend in hydrological series. Hydrol. Process. 2002, 16, 1807–1829. [Google Scholar] [CrossRef]
  47. Serinaldi, F.; Kilsby, C.G.; Lombardo, F. Untenable nonstationarity: An assessment of the fitness for purpose of trend tests in hydrology. Adv. Water Resour. 2018, 111, 132–155. [Google Scholar] [CrossRef]
  48. Clarke, R.T. On the (mis)use of statistical methods in hydroclimatological research. Hydrol. Sci. J. 2010, 55, 139–144. [Google Scholar] [CrossRef]
  49. Poppick, A.; Moyer, E.J.; Stein, M.L. Estimating trends in the global mean temperature record. Adv. Stat. Climatol. Meteorol. Oceanogr. 2017, 3, 33–53. [Google Scholar] [CrossRef] [Green Version]
  50. Hirsch, R.M.; Slack, J.R.; Smith, R.A. Techniques of trend analysis for monthly water quality data. Water Resour. Res. 1982, 18, 107–121. [Google Scholar] [CrossRef]
  51. Dommenget, D. Evaluating EOF modes against a stochastic null hypothesis. Clim. Dyn. 2007, 28, 517–531. [Google Scholar] [CrossRef]
  52. Dommenget, D. An objective analysis of the observed spatial structure of the tropical Indian Ocean SST variability. Clim. Dyn. 2011, 36, 2129–2145. [Google Scholar] [CrossRef]
  53. Hu, Y.M.; Liang, Z.M.; Liu, Y.W.; Wang, J.; Yao, L.; Ning, Y. Uncertainty analysis of SPI calculation and drought assessment based on the application of Bootstrap. Int. J. Climatol. 2015, 35, 1847–1857. [Google Scholar] [CrossRef]
  54. Richman, M.B. Review article: Rotation of principal components. J. Climatol. 1986, 6, 293–335. [Google Scholar] [CrossRef]
Figure 1. Locations of northeastern China as well as meteorological stations.
Figure 1. Locations of northeastern China as well as meteorological stations.
Water 11 00908 g001
Figure 2. Trends of precipitation and temperature over northeastern China. (The Z values of ±1.64, ±1.96 and ±2.58 are equal to the confidence levels of 0.1, 0.05 and 0.01, respectively.).
Figure 2. Trends of precipitation and temperature over northeastern China. (The Z values of ±1.64, ±1.96 and ±2.58 are equal to the confidence levels of 0.1, 0.05 and 0.01, respectively.).
Water 11 00908 g002
Figure 3. Areal average precipitation anomaly curves. (a) Annual. (b) Crop growing season. (The dotted lines represent the decadal averages).
Figure 3. Areal average precipitation anomaly curves. (a) Annual. (b) Crop growing season. (The dotted lines represent the decadal averages).
Water 11 00908 g003
Figure 4. Areal average temperature anomaly curves. (a) Annual. (b) Crop growing season. (The dotted lines represent the decadal averages).
Figure 4. Areal average temperature anomaly curves. (a) Annual. (b) Crop growing season. (The dotted lines represent the decadal averages).
Water 11 00908 g004
Figure 5. Areal average SPEI series over different time scales.
Figure 5. Areal average SPEI series over different time scales.
Water 11 00908 g005
Figure 6. Estimations of the SPEI6May–Sep considering the sampling uncertainty.
Figure 6. Estimations of the SPEI6May–Sep considering the sampling uncertainty.
Water 11 00908 g006
Figure 7. Spatial structure of long-term SPEI trends at various timescales. (The Z values of ±1.64, ±1.96 and ±2.58 are equal to the confidence levels of 0.1, 0.05 and 0.01, respectively.).
Figure 7. Spatial structure of long-term SPEI trends at various timescales. (The Z values of ±1.64, ±1.96 and ±2.58 are equal to the confidence levels of 0.1, 0.05 and 0.01, respectively.).
Water 11 00908 g007
Figure 8. Patterns of the first two DEOFs for SPEI at various timescales.
Figure 8. Patterns of the first two DEOFs for SPEI at various timescales.
Water 11 00908 g008aWater 11 00908 g008b
Figure 9. DPC scores of the first two DEOFs for SPEI at crop growing seasonal scale.
Figure 9. DPC scores of the first two DEOFs for SPEI at crop growing seasonal scale.
Water 11 00908 g009
Table 1. Categories of dryness/wetness degree according to the SPEI values.
Table 1. Categories of dryness/wetness degree according to the SPEI values.
CategoriesSPEI Values
Extremely wet[2.0, +∞)
Severely wet[1.5, 2.0)
Moderately wet[1.0, 1.5)
Slightly wet[0.5, 1.0)
Near normal(−0.5, 0.5)
Slightly dry(−1.0, −0.5]
Moderately dry(−1.5, −1.0]
Severely dry(−2.0, −1.5]
Extremely dry(−∞, −2.0]
Table 2. Percentage of variance explained by the first two DEOFs of the SPEI at various timescales.
Table 2. Percentage of variance explained by the first two DEOFs of the SPEI at various timescales.
CategoriesExplained Variance
DEOF1 (%)DEOF2 (%)
SPEI145.213.4
SPEI344.613.6
SPEI643.313.7
SPEI1240.213.9
SPEI6May–Sep42.514.1

Share and Cite

MDPI and ACS Style

Ye, L.; Shi, K.; Zhang, H.; Xin, Z.; Hu, J.; Zhang, C. Spatio-Temporal Analysis of Drought Indicated by SPEI over Northeastern China. Water 2019, 11, 908. https://doi.org/10.3390/w11050908

AMA Style

Ye L, Shi K, Zhang H, Xin Z, Hu J, Zhang C. Spatio-Temporal Analysis of Drought Indicated by SPEI over Northeastern China. Water. 2019; 11(5):908. https://doi.org/10.3390/w11050908

Chicago/Turabian Style

Ye, Lei, Ke Shi, Hairong Zhang, Zhuohang Xin, Jing Hu, and Chi Zhang. 2019. "Spatio-Temporal Analysis of Drought Indicated by SPEI over Northeastern China" Water 11, no. 5: 908. https://doi.org/10.3390/w11050908

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop