Advances in Meteorology
Volume 2020 (2020), Article ID 2409068, 26 pages
https://doi.org/10.1155/2020/2409068

Research Article

Meteorological Drought, Hydrological Drought, and NDVI in the Heihe River Basin, Northwest China: Evolution and Propagation

Fanglei Zhong,¹ Qingping Cheng,1,2 and Ping Wang³

¹Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou 730000, China
²University of the Chinese Academy of Sciences, Beijing 100039, China
³School of Tourism and Geographical Sciences, Yunnan Normal University, Kunming 650500, China

Correspondence should be addressed to Qingping Cheng; qpchengtyli@foxmail.com

Received 24 October 2019; Revised 10 February 2020; Accepted 22 February 2020; Published 9 June 2020

Academic Editor: Budong Qian

Copyright © 2020 Fanglei Zhong et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Understanding the evolution and propagation of different drought types is crucial to reduce drought hazards in arid and semiarid regions. Here, Standardized Precipitation Index (SPI), Streamflow Drought Index (SDI), and Vegetation Condition Index (VCI) were used to investigate the spatiotemporal variation of different drought types and correlations between Pre (Pre-R)/post (Pos-R)-reservoir. Results showed that the average peak/intensity/duration/severity of meteorological droughts (MD) were greater in the Pre-R than in the Pos-R period in the upstream Heihe River Basin (UHRB), while there was little change between the Pre-R and Pos-R periods in the midstream Heihe River Basin (MHRB). The average peak/intensity/duration/severity of hydrological drought (HD) decreased in the mainstream for Yingluoxia (Ylx) but increased for Zhengyixia (Zyx) station in the Pos-R period. Propagation time decreased by 3 months (negative effect) in Ylx and increased by 8 months (positive effect) in Zyx compared with the Pre-R period. In the Pos-R period, propagation time increased (1–3 months) for tributaries (positive effect). Propagation times for the mainstream and tributaries varied for different seasons and time periods. Pearson’s correlation coefficient values were lower at short timescales (1–3 months) but higher at long timescales for the Pos-R period in Ylx and Zyx for SDI-1 with different timescales of SPI. The SDI and SPI had no lag in the UHRB and MHRB. However, VCI with SPI had a significant lag correlation at short timescales in the UHRB (lag 6 months) and MHRB (lag 4 months), and the VCI with SDI had a significant lag correlation for 1 month in the MHRB. The propagation time from MD to HD has been reduced for Pos-R in the UHRB. There was a positive effect (prolonged MD propagation HD time) in Pos-R but still faces serious drought stress in the MHRB.

1. Introduction

Drought is recognized as the world’s most costly and pressing natural hazard, giving rise to significant losses in the fields of the economy, ecology, and environment (e.g., crop losses, degradation and desertification, urban water supply shortages, and forest fires) [1–4]. The damage from droughts is expected to increase in severity [5]. Compared with other natural hazards, the spatial extent of drought is very large and its time of influence is commonly much longer. Under global warming, more frequent and severe droughts have been projected in the 21st century [6–8], particularly in the mid-latitudes [9, 10]. Severe and extreme droughts have become more frequent in the eastern and central portions of northern China since the late 1990s [11, 12].

In general, drought can be classified into four types: meteorological, hydrological, agricultural, and socioeconomic [13]. Meteorological droughts (MD) are the water shortages caused by an imbalance in precipitation and evaporation; precipitation is commonly used for MD analysis [14, 15]. Hydrological droughts (HD) are related to a period with inadequate surface and subsurface water resources for established water uses in a given water resource management system. Streamflow data are widely applied for HD analysis [14, 16]. Because the influence of drought on terrestrial ecosystems is becoming increasingly acute, the vegetation response to drought is a crucial topic in the domain of climate research [17–19]. The Normalized Difference Vegetation Index (NDVI) is a good measure to estimate green biomass, leaf area index, and patterns of productivity and has been widely used to assess vegetation degradation, ecosystem features, and the physiological drought conditions of vegetation [20, 21].

Drought indices are the best ways to monitor drought and drought events at present. Lloyd-Hughes [22] counted over 100 drought indicators that have been developed for different types of drought [23]. A commonly used index is the Standardized Precipitation Index (SPI) [24]. Its calculation is simple because it only needs precipitation as an input, and it can be calculated over different timescales, which allows SPI to monitor both short-term droughts, such as agricultural droughts, and long-term droughts, such as HD [3, 25–27]. For HD index, the Streamflow Drought Index (SDI) proposed by Nalbantis and Tsakiris [28] overcomes the problems of predicting drought onset and duration using cumulative streamflow volumes and the areal extent of drought using a spatially integrated streamflow at a basin outlet; it is also widely used in different regions of the world [28–33]. The Vegetation Condition Index (VCI) is a normalization of the NDVI, which filters out the contribution of local geographic resources to the spatial variability of the NDVI [17, 34]. It reflects agricultural (vegetation) drought accurately in many parts of the world and has been applied in daily drought monitoring by the National Oceanic and Atmospheric Administration of the United States (NOAA) and the National Satellite Meteorological Centre of China [35].

Many researchers in related fields have focused on drought evolution, drought prediction, and developing reliable drought indices using observed and simulated data [4, 36–39]. Recently, studying the propagation time of VCI and HD based on MD using Pearson’s or Spearman’s correlation coefficients has become a research focus [8, 10, 18, 19, 21, 23, 40–46]. These two kinds of droughts can reflect the different stages of drought development. Commonly, MD develops and ends relatively quickly, while HD is the result of MD [4]. The corresponding propagation time depends on local landscape conditions [46]. Therefore, investigating the propagation time from MD to HD [4] and from MD to VCI along with its potential influencing factors is necessary for establishing an effective monitoring and warning system for HD and VCI based on MD.

However, little research has been done into the propagation of MD to HD and VCI and HD propagation to VCI, in particular in subarid and arid watersheds, which is one of the aims of this paper. The investigation of VCI and HD propagation time from MD at different timescales under the influence of global warming and human activities (e.g., construction of reservoirs or dams, land use, and land cover change) will help to elucidate the impacts of droughts on terrestrial ecosystems and inform planning and optimizing water resource allocations during droughts. It is important to identify the drought propagation process and mechanisms to help establish a drought early warning system. The primary objectives of this study are as follows: (1) to inspect the spatiotemporal characteristics of MD and HD duration, intensity, severity, and peak, while detecting the spatiotemporal characteristics of VCI; (2) to investigate VCI and HD propagation times based on MD at different timescales and VCI propagation times based on HD at different timescales; and (3) to explore the potential factors influencing drought characteristics and propagation times.

2. Materials and Methods

2.1. Study Area

The Heihe River Basin (HRB) is the second largest inland river basin in the arid zone of northwest China. The HRB flows through Qinghai Province, Gansu Province, and the Inner Mongolia Autonomous Region. The HRB is located along the land Silk Roads (“One Belt, One Road”) between 97.1°E–102.0°E and 37.7°N–42.7°N with a total area of approximately 1.43 × 10⁵ km² [47, 48]. The HRB includes three sections from south to north: the upstream Heihe River Basin (UHRB) from the Qilian Mountains to Yingluoxia station (outlet of the mountains) belongs to the subhumid and semiarid temperate continental monsoon climate zone; the midstream Heihe River Basin (MHRB) running from Yingluoxia station to Zhengyixia station belongs to the arid climate zone; and downstream HRB terminates in the Juyan Lakes (east and west branches, respectively) [49]. From upstream to downstream, the elevation decreases, the water availability decreases, and the landscape changes from glaciers and alpine biomes in the UHRB to steppes and agricultural ecosystems in the MHRB to riparian ecosystems and vast areas of desert in the downstream area [47]. In this study, the UHRB and MHRB were selected as the study area. Figure 1 shows the study area and the hydrometeorological stations. Analyses were also conducted at the sub-basin level in the UHRB and MHRB. In the MHRB, the annual average precipitation, potential evapotranspiration, and streamflow were 226 mm, 924 mm, and 88 mm, respectively. Over 90% of the population, grain production, and major industries are concentrated in the MHRB. Approximately 84% of the total available water was consumed for irrigation. On average, a decline in groundwater of about 1.86 m is attributable to water consumption for irrigation of the farmland area in Heihe River from 1981 to 2010 [50] with demand constantly increasing [51, 52]. There are a series of reservoirs (one large reservoir, 9 medium size reservoirs, and 89 small reservoirs) and a relatively complete irrigation system consisting of more than 893 main canals and branch canals [53, 54]. Land use and land cover have also significantly changed in the HRB. Urban, cultivated, and forest land have increased, while water area and grassland have become increasingly degraded.

(a)

(b)

(c)

Figure 1: The upstream and midstream in the Heihe River Basin with meteorological and hydrological stations and digital elevation map (a), land use and land cover (LULC) in 2011 (b), and vegetation types (c).

2.2. Data

Daily precipitation datasets from 1961 to 2012 were collected by five meteorological stations (Table 1) and derived from the China Meteorological Data Sharing Service System of National Meteorological Information Centre V3.0 (http://www.nmic.gov.cn/). Rigorous quality control was conducted by the China National Meteorological Information Centre before the data were released. Monthly observed streamflow data, which were also of high quality, for eight hydrological stations (Table 2) were collected from the Hydrological Bureau of Gansu Province. Yingluoxia (Ylx) and Zhengyixia (Zyx) stations were located in the main stream of the HRB. The monthly NDVI (1 km spatial resolution) dataset was the satellite remote sensing data of SPOT/VEGETATION NDVI and MODIS from 1999 to 2012. The data were obtained from the resource and environmental science data center of the Chinese Academy of Sciences (http://www.resdc.cn), and LULC data in the UHRB (2000 and 2010) and MHRB (1975, 1985, 1995, 2000, and 2010) were obtained from the Cold and Arid Regions Science Data Centre (http://westdc.westgis.ac.cn/). The historical disaster records were derived from the China Meteorological Disaster Yearbook (Gansu volume). Zhang et al. [49] used Pettitt tests to determine abrupt change points in the annual streamflow series from 1945 to 2012 and from 1957 to 2012 at the Ylx and Zyx stations to ascertain the impact of the reservoir activities in the HRB. The results indicated only one significant abrupt change point in 1979 at Ylx station and in 1984 at Zyx station [49]. The impact of human activities is more evident after the Chinese Government initiated the Ecological Water Diversion Project (EWDP) in mainstream Heihe River for 2000 [54]. Therefore, in this paper, we focused on the mainstream Ylx and Zyx stations, between the Pre-R and Pos-R periods, before and after EWDP. The abrupt change points for the tributaries for six stations are shown in Table 2 (Pettitt tests).

Table 1: Information on the meteorological stations.

Sub-basin
Station
Longitude (°E)
Latitude (°N)
Elevation (m)
Data period

UHRB
Qilian (Ql)
100.24°
38.19°
2787.4
1961–2012

Yeniugou (Yng)
99.58°
38.41°
3286.0
1961–2012

MHRB
Zhangye (Zy)
100.46°
38.91°
1482.7
1961–2012

Shandan (Sd)
101.08°
38.77°
1764.6
1961–2012

Gaotai (Gt)
99.79°
39.36°
1332.2
1961–2012

Table 2: Information on the hydrological stations.

Sub-basin
Station
River
Longitude (°E)
Latitude (°N)
Catchment area (km²)
Elevation (m)
Data period
Abrupt year

UHRB
Yingluoxia (Ylx)
Heihe River
100.18
38.82
10009
1700
1945–2012
1979

Shuangshusi (Sss)
Hongshui River
100.83
38.32
578
2490
1957–2012
1974

Wafangcheng (Wfc)
Dazhuma River
100.48
38.43
334
2440
1957–2012
1985

Liyuanbao (Lyb)
Liyuan River
100.00
38.97
2641
1760
1956–2012
1976

MHRB
Zhengyixia (Zyx)
Heihe River
99.47
39.82
35634
1280
1957–2012
1984

Lijiaqiao (Ljq)
Maying River
101.13
38.52
1205
2150
1957–2012
1986

Shunhua (Sh)
Haichaoba River
100.72
38.59
191.8
1948
1957–2012
1995

Suyoukou (Syk)
Suyou River
100.58
38.79
473.4
1565
1957–2012
1971

Abrupt change year was detected by the Pettitt test.

2.3. Methods

2.3.1. Standardized Precipitation Index (SPI)

The SPI index was developed by McKee et al. [24]. Within a certain geographic area, the precipitation usually fluctuates regularly. If the precipitation is less than the average annual precipitation, a drought may occur. If precipitation exceeds the annual average, flooding may occur. To calculate the SPI, a frequency distribution function is first constructed from a series of long-term precipitation observations. In this paper, a three-parameter log-logistic distribution was used because the frequencies of precipitation accumulated at different timescales are well modeled using these statistical distributions. Because SPI is based on the cumulative probability of a given timescale, here the total amounts of precipitation in the current month and previous i months (i = 1, 2, 3, ......) were used to calculate SPI on a timescale of i + 1 month. In our study, we used the SPI (1, 3, 6, 9, ......, 48 months) to analyze the characteristics of drought change in the HRB. For drought classifications, please see Table 3.

Table 3: Classification of SPI, SDI, and VCI.

Value
Drought grades

SPI/SDI ≤ −2.0, VCI < 15
Extreme drought

−2.0 < SPI/SDI ≤ −1.5, 15 ≤ VCI < 30
Severe drought

−1.5 < SPI/SDI ≤ −1.0, 30 ≤ VCI < 45
Moderate drought

−1.0 < SPI/SDI ≤ 0, 45 ≤ VCI < 60
Mild drought

2.3.2. Streamflow Drought Index (SDI)

The SDI was developed by Nalbantis and Tsakiris [28] to characterize hydrological drought based on developing the concepts of the SPI. It is worth noting that the procedure should be performed separately for each month and can be calculated for any timescale (e.g., 1, 3, 6, 9, and 12 months). The total streamflow in a given month j and year i depends on the chosen timescale k and can be calculated as [55]

where Vi-1,l and Vi,l denote streamflow volumes in years i − 1 and i, respectively. Due to the three-parameter log-logistic distribution used to calculate the SPI, we also chose this probability distribution to standardize the observed streamflow data used to get the SDI to facilitate comparison in the HRB. It is emphasized here that positive SDI values reflect wet conditions, while negative values indicate a hydrological drought.

2.3.3. Vegetation Condition Index (VCI)

Kogan and Sullivan [56] proposed the VCI index in 1993. It is here assumed that the change of the VCI index is only affected by weather factors; in this way, the ecosystem factors and extreme weather factors are separated, the impact of extreme weather on vegetation is evaluated, and agricultural drought is monitored. A smaller VCI index correlates to a greater drought. The VCI index is calculated as follows:

where NDVI, NDVImax, and NDVImin are monthly NDVI, multiyear maximum NDVI, and multiyear minimum NDVI (in this study, 14 years), respectively, for each grid cell. The VCI changes from 0 to 100 corresponding to changes in vegetation condition from extremely unfavorable to optimal. In the case of an extremely dry month, the vegetation condition is poor, and the VCI is close to or equal to zero. A VCI of 50 reflects fair vegetation conditions. At optimal vegetation conditions, the VCI is close to 100 [34]. According to t-test, the trend was divided into extremely significant increase (slope > 0, p < 0.01), significant increase (slope > 0, 0.01 ≤ p < 0.05), insignificant change (p ≥ 0.05), extremely significant decrease (slope < 0, p < 0.01), and significant decrease (slope < 0, 0.01 ≤ p < 0.05).

2.4. Drought Variables and Run Theory

According to McKee et al. [24] and Spinoni et al. [57], the derived drought variables based on run theory [58] follow certain definitions (Figure 2). When the monthly SPI/SDI values were below −1, the corresponding month was potentially identified as a drought event (e.g., the four drought events for D₀, D₁, d₀, and d₂ shown in Figure 2). Thresholds of −1 (moderate drought), −1.5 (severe drought), and −2 (extreme drought) were used to identify drought events and drought characteristics. Drought duration was defined as the number of months from the first month in which the indicator goes lower than −1 to the last month with a negative value before the indicator returns to positive values. If the duration of a drought event only contained one month where SPI/SDI < −1, this month was regarded as a single drought event (e.g., D₁ with a magnitude S₁). Drought severity was defined as the sum of the monthly absolute values of the index when the index was ≤−1.0 (−1.5, −2). A drought event may contain a few consecutive months with negative SPI/SDI (e.g., duration D₀ and its magnitude S₀). If the duration of a drought event contained two branches, such as d₀ and d₂, and the interruption period d₁ between d₀ and d₂ was less than 6 months in which 0 < SPI/SDI < 1 [59, 60], these months were still regarded as a single drought event (D₂ = d₀ + d₁ + d₂). The corresponding magnitude was then defined as S₂ = s₁ + s₂ [59, 60]. Intensity stands for the number of months in which the drought indicator was lower than −1 (−1.5, −2). Drought peak referred to the month in the “drought event” with the lowest value of the indicator [61].

Figure 2: Definition of drought characteristics including peak, intensity, duration, and magnitude using the SPI/SDI and run theory.

2.5. Correlation Analysis

Pearson’s correlation coefficient (PC) was calculated between SDI-1 (one month SDI) and SPI (1, 3, 6, 9, ....., 48 months), monthly (April–October) VCI, and different timescales of SPI and SDI (1, 3, 6, 9, and 12 months) series indicates the characteristics of interannual propagation at 95% and 99% confidence levels. The PC ranges from −1 to 1; a PC > 0 indicates a positive correlation and a PC < 0 indicates a negative correlation. In this study, the PC method was also used to detect the intra-annual propagation time of the SDI-1 to the different timescales of SPI and monthly VCI to the different timescales of SDI and SPI conditions. For instance, if the maximum correlation in the March SDI and the March series of SPI at the timescales of 1, 3, …, 48 months is recorded at the 3-month timescale, this means that the SDI in March is mostly determined by the cumulative water balance occurring from January to March [19, 21, 23, 41]. Meanwhile, to determine whether there is a lag between the SPI (accumulation periods of 1–48 months) and SDI-1, cross-correlations were calculated for SDI-1 series which were lagged by 0 to 6 months after the SPI series. In this case, the SPI accumulation period with the strongest correlation with SDI-1 was denoted as the lagged SPI-n (1, 3, 6, 9, ....., 48 months) [23]. The same interpretation applies to the lag correlation of the VCI with different timescales SPI and VCI with different timescales SDI.

3. Results

3.1. MD and HD for the Mainstream Characteristics between Different Timescales and Thresholds

The distribution of drought characteristics based on thresholds of −1, −1.5, and −2 for different timescales at different periods, including Pre-R (1965–abrupt change year), Pos-R (abrupt change year +1–2012), and total period (1965–2012), was evaluated in the UHRB and MHRB using boxplots. The drought characteristics of SPI-1 and SDI-1 in the MHRB (Table S1) were almost consistent compared with historical records and, for different combinations of timescales, thresholds and periods (Table S2) are shown in the supporting information. The drought peak and intensity (Figures S1-S2), duration, and severity for different timescales at different thresholds are shown in Figure 3 for MD and in Figure 4 for HD. More drought events were identified at shorter accumulation periods (1–12 months for MD and 1–6 months for HD) based on thresholds equal to −1.0, −1.5, and −2.0 in the UHRB (Ylx) and MHRB (Zyx). For MD, the average intensity, duration, and severity of fluctuations increased with cumulative timescale increases between Pre-R and Pos-R periods. The changes were more obvious in Pos-R compared with Pre-R periods for all thresholds in the UHRB. The average peak, intensity, duration, and severity of the MD did not change much for Pre-R, whereas Pos-R fluctuations increased with cumulative timescale increases using thresholds of −1.0, −1.5, and −2.0. The average peak, intensity, duration, and severity of value changed little in the MHRB. With regard to HD, the duration and severity were larger for HD than for MD because a short-term MD may not develop into an HD [42] and as the number of events decreases, duration lengthens, and severity worsens. The average intensity (Figure S2), duration, and severity of the HD were significantly different in Ylx and Zyx for different periods based on the thresholds −1.0, −1.5, and −2.0. In Ylx, the characteristics of HD were similar to MD in that, compared to the Pre-R period, the average intensity, duration, and severity of the HD decreased significantly, especially for SDI-3 and SDI-12. However, in Zyx, the HD were quite inconsistent with MD, and the average intensity, duration, and severity of the HD increased significantly during the Pos-R period in Zyx, especially for SDI-1, SDI-6, and SDI-12.

(a)

(b)

(c)

(d)

Figure 3: Boxplots showing MD duration and severity in the UHRB (a, c) and in the MHRB (b, d) based on SPI using thresholds of −1, −1.5, and −2 for different timescales at different periods (the dot and midline in the boxplots are the mean and median; the upper and lower boundaries of box are the interquartile range).

(a)

(b)

(c)

(d)

Figure 4: Boxplots showing HD duration and severity for Ylx (a, c) and for Zyx (b, d) based on SDI using thresholds of −1, −1.5, and −2 for different timescales at different periods (the dot and midline in the boxplots are the mean and median; the upper and lower boundaries of box are the interquartile range).

3.2. Comparison of MD Characteristics for SPI-1 and HD for SDI-1

The characteristics of MD and HD based on SPI-1 and SDI-1 using thresholds of −1, −1.5, and −2 in the UHRB and MHRB are shown in Tables 4-5 and Figure 5. More MD events occurred at Yng station (names and abbreviations are shown in Table 1) in the UHRB (Figure 5) and at Sd station in the MHRB. More HD events occurred at Zyx (Figure 5). It can be seen from Tables 4-5 that the average, maximum, and minimum peak/intensity/duration/severity in Zyx, Yng, and Ql stations were identified as longer and more severe based on thresholds of −1.0, −1.5, and −2.0 for MD and HD. HD characteristics based on these thresholds for the average peak, maximum peak, and minimum peak in Sh, respectively, were more obvious, and longer average, maximum, and minimum intensities were found in Ylx, Ljq, and Sh (Table 5). Longer, more severe HD occurred in Ylx, Ljq, Sh, and Wfc.

Table 4: The distribution peak, intensity, duration, and severity for MD characteristics based on SPI-1 using thresholds of −1, −1.5, and −2 in the UHRB and MHRB.

Threshold
Peak
Intensity
Duration
Severity

Mean
Max
Min
Mean
Max
Min
Mean
Max
Min
Mean
Max
Min

Ql
−1.0
1.7
3.6
1.1
1.7
4.0
1.0
4.6
15.0
1.0
3.1
10.2
1.1

−1.5
2.0
3.6
1.2
1.3
3.0
1.0
4.6
15.0
1.0
3.6
10.2
1.5

−2.0
2.5
3.6
2.0
1.0
1.0
1.0
4.3
15.0
1.0
4.1
10.2
2.0

Yg
−1.0
1.6
3.1
1.7
1.7
5.0
1.0
4.4
23.0
1.0
3.8
12.9
1.0

−1.5
2.2
3.1
2.5
1.6
3.0
1.0
5.5
23.0
1.0
5.4
12.9
1.8

−2.0
2.5
3.1
2.0
1.2
2.0
1.0
6.9
23.0
1.0
5.1
12.5
2.0

Gt
−1.0
1.4
2.9
1.0
1.5
4.0
1.0
4.1
14.0
1.0
3.1
8.1
1.0

−1.5
1.8
2.9
1.5
1.1
2.0
1.0
4.5
14.0
1.0
3.6
8.1
1.6

−2.0
2.7
2.9
2.5
1.0
1.0
1.0
2.0
3.0
1.0
3.5
4.0
2.9

Zy
−1.0
1.6
2.7
1.0
1.6
6.0
1.0
4.7
30.0
1.0
3.5
18.8
1.0

−1.5
2.0
2.7
1.5
1.5
4.0
1.0
5.5
30.0
1.0
4.4
18.8
1.9

−2.0
2.3
2.7
2.0
1.1
2.0
1.0
4.0
15.0
1.0
3.9
10.0
2.0

Sd
−1.0
1.6
3.0
1.0
1.6
5.0
1.0
4.2
17.0
1.0
3.2
9.4
1.0

−1.5
1.9
3.0
1.5
1.2
3.0
1.0
3.5
14.0
1.0
3.4
9.4
1.5

−2.0
2.3
3.0
2.0
1.1
2.0
1.0
3.1
11.0
1.0
3.7
9.4
2.0

Table 5: The distribution of peak, intensity, duration, and severity for HD characteristics based on SDI-1 using thresholds of −1, −1.5, and −2 in the UHRB and MHRB.

Threshold
Peak
Intensity
Duration
Severity

Mean
Max
Min
Mean
Max
Min
Mean
Max
Min
Mean
Max
Min

Ylx
−1.0
1.6
2.7
1.1
4.3
18.0
1.0
11.1
31.0
1.0
8.6
30.8
1.1

−1.5
1.9
2.7
1.5
4.0
18.0
1.0
12.3
31.0
1.0
10.9
30.8
1.5

−2.0
2.3
2.7
2.1
1.8
3.0
1.0
19.5
31.0
10.0
18.0
24.6
8.9

Sss
−1.0
1.8
3.3
1.0
3.7
14.0
1.0
12.4
36.0
1.0
8.1
24.7
1.0

−1.5
2.0
3.3
1.5
2.1
5.0
1.0
12.4
30.0
1.0
8.5
22.2
1.5

−2.0
2.6
3.3
2.0
1.3
2.0
1.0
9.5
20.0
2.0
7.7
12.2
2.2

Wfc
−1.0
1.7
3.0
1.0
3.4
17.0
1.0
8.9
51.0
1.0
7.0
37.0
1.2

−1.5
2.0
3.0
1.6
2.7
8.0
1.0
12.2
51.0
1.0
9.8
37.0
1.6

−2.0
2.5
3.0
2.0
2.0
4.0
1.0
19.3
54.0
6.0
15.5
40.4
6.6

Lyb
−1.0
1.8
3.0
1.0
4.1
16.0
1.0
10.2
30.0
1.0
8.7
35.0
1.0

−1.5
2.1
3.0
1.6
2.8
8.0
1.0
10.4
30.0
1.0
9.9
35.0
1.7

−2.0
2.4
3.0
2.0
1.4
4.0
1.0
8.3
28.0
1.0
9.4
32.0
2.0

Zxy
−1.0
1.8
5.0
1.0
3.0
21.0
1.0
7.0
40.0
1.0
6.2
41.3
1.1

−1.5
2.2
5.0
1.5
2.6
12.0
1.0
6.9
38.0
1.0
7.6
38.9
1.5

−2.0
2.7
5.0
2.0
1.6
6.0
1.0
9.5
37.0
1.0
10.2
37.3
2.1

Ljq
−1.0
1.8
4.0
1.0
3.4
17.0
1.0
7.0
32.0
1.0
6.8
33.7
1.1

−1.5
2.4
4.0
1.5
3.3
12.0
1.0
10.3
32.0
1.0
10.6
33.7
1.5

−2.0
2.9
4.0
2.1
3.0
12.0
1.0
13.6
32.0
3.0
14.0
29.8
3.8

Sh
−1.0
2.1
4.3
1.0
4.3
18.0
1.0
12.9
43.0
1.0
9.7
36.3
1.7

−1.5
2.5
4.3
1.6
3.0
9.0
1.0
8.9
33.0
1.0
9.2
35.3
1.7

−2.0
2.8
4.3
2.0
1.8
5.0
1.0
8.4
32.0
1.0
10.3
33.7
2.3

Syk
−1.0
1.7
3.8
1.0
2.4
8.0
1.0
7.6
36.0
1.0
5.2
18.7
1.0

−1.5
2.0
3.8
1.5
1.6
3.0
1.0
7.6
36.0
1.0
5.9
18.7
1.8

−2.0
2.9
3.8
2.4
2.0
3.0
1.0
10.8
35.0
2.0
9.2
18.7
3.2

(a)

(b)

Figure 5: Spatial distribution of MD (SPI-1) and HD (SDI-1) events in the UHRB and MHRB.

3.3. Variation in VCI

Vegetation condition in the extremely dry months of April, May, and October (Figure 6) was poor, especially for alpine vegetation, cultivated vegetation, grassland, and marshy grassland. From June to September, the vegetation was in optimal condition for coniferous forests, shrubs, and marshy grassland in the UHRB and for coniferous forests and cultivated vegetation in the MHRB. It is worth noting that the vegetation condition was poor for alpine vegetation over the whole growing season.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 6: The spatial variation in the average VCI for the growing season (April–October) in the UHRB and MHRB.

We further analyzed the trend and significance level of the growing season for VCI. Figure 7 shows that, over the entire growing season, alpine vegetation (UHRB) and grassland (MHRB) showed a significant decreasing trend, while cultivated vegetation (MHRB), coniferous forest, shrubs, and grassland (UHRB) showed a significant increasing trend.

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 7: The spatial distribution of VCI trends for the growing season in the UHRB and MHRB.

3.4. Correlation of VCI, HD, and MD

3.4.1. Interannual Propagation Time between VCI, HD, and MD at Different Timescales in Different Periods

Figures 8-9 describe the relationship between SDI-1 and VCI-1 with all SPI accumulation periods (1–48 months) and VCI with all SDI accumulation periods (1–12 months) for different periods (Pre-R and Pos-R, before and after EWDP, and total period) in the mainstream for Ylx (SPI calculated using the arithmetic mean of precipitation from Yng and Ql stations), Zyx (SPI calculated using the arithmetic mean of precipitation from Gt, Sd, and Zy stations), and tributaries (SPI calculation is based on the area flowing through or the adjacent area because all meteorological stations were not available, e.g., Sh-Sd, Sss-Ql, Lyb-Yng, Ljq-Sd, Wfc-Ql, and Syk-Zy). Figure 8 indicates that the propagation times were 1, 1, 3, 1, 1, and 1 in the Pre-R period, respectively; 1, 3, 9, 3, 6, and 6 in the Pos-R period, respectively; and 1, 3, 6, 1, 6, and 6 in the total period for six tributaries, respectively. Figure 9(a) shows that the SDI-1 and SPI had the highest PC (PC = 0.27, ) in the Pre-R period for the 9-month period, and, in the Pos-R period, the highest PC (PC = 0.46, ) was observed for the 6-month period. Figure 9(a) also shows the highest PC for SDI-1 and SPI (PC = 0.38, ) was for the 6-month period before the EWDP period, and, after the EWDP period, the highest PC (PC = 0.48, ) was observed for the 6-month period in Ylx. Figure 9(b) reveals the highest PC (PC = 0.19, ) was found for the 1-month period Pre-R and the 9-month period (PC = 0.24, ) Pos-R. Figure 8(b) also demonstrates the highest PC (PC = 0.17, ) was found for the 3 months before the EWDP period and over 36 months (PC = 0.35, ) after the EWDP period in Zyx. These results demonstrate that the changes in the streamflow were more sensitive to precipitation in the Pre-R period for short timescales (1–3 months), and Pos-R propagation time increased (2–6 months) for tributaries, while the propagation time decreased by 3 months in Ylx and increased by 8 months in Zyx compared to the Pre-R period. Thus, reservoir operation increased the propagation time of the streamflow variables to the SPI during the Pos-R period in Zyx, whereas the propagation time of the streamflow variables to the SPI was decreased during the Pos-R period in Ylx. There was no change in Ylx before and after the EWDP, while after the EWDP, the propagation time of the HD to the MD increased on a 33-month timescale compared with before the EWDP in Zyx. Figure 9(c) shows that the VCI and SPI strongest negative PC was found for the 6-month (PC = −0.09) period in the UHRB, and the strongest positive PC was found for the 36-month (PC = 0.13) period in the MHRB. Figure 9(d) shows that the VCI and SDI highest PC was found for the 6-month period (PC = 0.23, ) in the UHRB and the 3-month period (PC = 0.43, ) in the MHRB. The results indicated that correlation between VCI and MD was not significant in the UHRB and MHRB. However, the correlation between VCI and HD was significant in the UHRB (Ylx) and MHRB (Zyx), while the VCI propagation time to HD was found at short timescales.

(a)

(b)

(c)

Figure 8: Pearson’s correlation coefficien between the SDI-1 series and the SPI series at different timescales for six tributaries.

(a)

(b)

(c)

(d)

Figure 9: Pearson’s correlation coefficien between SDI-1 and SPI (a, b) for different timescales at different periods, and VCI (c, d) with different timescales of SPI and SDI in the UHRB (Ylx) and in the MHRB (Zyx).

3.4.2. Intra-Annual and Lag Correlations between VCI, HD, and MD at Multiple Scales in Different Periods

The SDI-1 with SPI at different timescales was investigated in relation to the propagation from MD to HD. Figure 10 shows the propagation times from MD to HD had significant seasonal characteristics for tributaries in different periods. In spring (March–May), the drought propagation time was relatively short (1–12 months) in Sss, Ljq, and Wfc, whereas the propagation time was relatively long (1–48 months) in Sh, Lyb, and Syk. In summer (June–August), the drought propagation time was relatively long, except for Wfc. Drought propagation times in autumn (September–November) were relatively short (1–6 months) in Sss, Lyb, and Syk and relatively long (1–48 months) in Sh, Ljq, and Wfc. Drought propagation times in winter (December–February) were similar to those in summer and were relatively long (1–48 months), except for Lyb in the Pre-R period. Drought propagation times in autumn and winter were relatively long, whereas in spring and summer, they were relatively short in the Pos-R period. The total period was similar to the Pos-R period. To summarize, in the Pos-R period, the drought propagation times in winter did not change and remained long, whereas drought propagation times were increased in autumn. Propagation time was longer for Ss, Lyb, and Syk and shorter for Sss, Ljq, and Wfc in spring; and propagation times in summer were shorter, except for Wfc and Ss. It can be clearly seen from Figure 11 that the intra-annual correlations between the monthly SDI-1 and SPI at multiple scales and in different periods have strong seasonal characteristics in Ylx and Zyx. The propagation times in spring, summer, autumn, and winter were approximately 3, 9, 6 and 12 months in the Ylx for the Pre-R period (Figure 11(a)), respectively. The propagation times in the spring, summer, autumn and winter were approximately 6, 9, 1, and 6 months in Ylx for the Pos-R period (Figure 11(b)), respectively. The propagation times in spring, summer, autumn, and winter were 6 months in Ylx for the total period (Figure 11(c)). The propagation times in spring, summer, autumn, and winter were approximately 6, 1, 1, and 15 months in Zyx for the Pre-R period (Figure 11(d)), respectively. The propagation times in spring, summer, autumn, and winter were 9, 12, 36, and 9 months in Zyx for the Pos-R period (Figure 11(e)), respectively. The propagation times in spring, summer, autumn, and winter were approximately 9, 9, 3, and 6 months in Zyx for the total period (Figure 11(f)), respectively. These results show that the propagation time in the Pos-R period increased in spring, decreased in autumn and winter, and did not change in summer in Ylx compared to the Pre-R period, whereas in Zyx it increased in spring, summer, and autumn, and it decreased in winter. Compared with the Pre-R period, the PC was lower on a 1-month timescale in Ylx (Figure 11(b)) and on a 1–3-month timescale in Zyx (Figure 11(e)). However, the PC between SDI-1 and different timescale of SPI in the Pos-R period was higher than that in the Pre-R period for long-term timescales (>12 months) in Ylx (except for January and June) and Zyx (except for June, July, November, and December). It can be seen from Figures S3(a)-S3(b) that, after the EWDP, the PC of SDI-1 and the timescale of SPI were higher in spring and autumn and lower in summer (except for August) and winter than before the EWDP in Ylx. However, Figures S3(c)-S3(d) show that, after the EWDP, the PC of SDI and SPI was higher in spring, summer, autumn, and winter, specifically in spring, in Zyx. In different periods (Pre-R and Pos-R and before and after EWDP) and different timescales, the PC of the Pos-R period and after the EWDP was higher in Ylx (except for 1, 36, and 48 months) and Zyx, especially in Zyx than in the Pre-R period and before the EWDP. Figure 12(a) shows that the higher PC between VCI and SPI in April at 6 months was 0.63 (), whereas the lower PC in May at 1 month and 48 months was −0.06. The highest PC between VCI and SDI in May at 12 months was 0.78 (), whereas the lowest PC in March at 3 months was −0.69 () in the UHRB (Figure 12(b)). Figure 12(c) shows that the highest PC between VCI and SPI was 0.84 () in July at 9 months, whereas the lowest PC was −0.19 in October at 1 month. The highest PC between VCI and SDI was observed at 0.71 in May at 6 months, and the lowest PC was −0.23 in March at 6 months in the MHRB (Figure 12(d)).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 10: Intra-annual strongest PC (a–c) and the drought propagation time in months (e–f) for SPI and SDI-1 at six tributaries.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 11: Heat map showing Intra-annual PC of SPI accumulation periods of 1–48 months with SDI-1 in the UHRB (Ylx) and MHRB (Zyx) (note: the x-axis represents the correlation between the SDI-1 and SPI at different timescales and the y-axis represents the 12 months of the year; the up arrow represents the correlation coefficients that increased in the postreservoir period (Pos-R) in the Ylx and Zyx compared to the period of prereservoir (Pre-R), or vice versa); (a) (1965–1979, PC ≥ 0.52, p < 0.05), (b) (1980–2012, PC ≥ 0.35, p < 0.05), and (c) (1965–2012, PC ≥ 0.29, p < 0.05) in the UHRB (Ylx), and (d) (1965–1984, PC ≥ 0.45 or ≤− 0.47, p < 0.05), (e) (1985–2012, PC ≥ 0.38 or ≤−0.48, p < 0.05), and (f) (1965–2012, PC ≥ 0.29 or ≤ −0.33, p < 0.05) in the MHRB (Zyx).

(a)

(b)

(c)

(d)

Figure 12: Heat map showing growing season Pearson’s correlation coefficients between VCI and different timescale of SPI, VCI, and different timescale of SDI for the period of 1999–2012 in the UHRB (Ylx). (a) (PC ≥ 0.53, p < 0.05), (b) (PC ≥ 0.55 or ≤ −0.69, p < 0.01) and MHRB (Zyx), (c) (PC ≥ 0.55, p < 0.05), and (d) (PC ≥ 0.54, p < 0.05).

We further analyzed whether there was a lag between the different timescales of SPI with SDI-1. The PC between SDI-1 lagged by 0–6 months, with SPI accumulation periods of 1–48 months for different timescales at different periods (Pre-R, Pos-R, and total), and the strongest correlation was detected at a lag of 0 months (i.e., no lag) in Ylx and Zyx (Figure 13). We also found no lag before and after the EWDP (Figure S5), whereas there was a higher PC at 3–12 months before the EWDP and after the EWDP at 6–15 months in Ylx. The highest PC was at 1–12 months before the EWDP and after EWDP at 9–48 months in Zyx. We also analyzed the lag correlation of the SPI with VCI-1 and SDI (accumulation periods of 1–12 months) with VCI-1 (Figure 14), and PC was calculated for the VCI-1 series, which were lagged by 0 to 6 months after the SPI and SDI series. The results demonstrated that the highest PC (PC = 0.21, ) was found at a 6-month lag between VCI and SPI, while the strongest correlation was at a lag of 0 months (i.e., no lag) between VCI and SDI in the UHRB (Ylx) during 1999–2012. The highest PC (PC = 0.26, ) was observed at the 4-month lag scale between VCI and SPI, and the strongest PC (PC = 0.47, ) was found at a lag of 1 month between VCI and SDI during 1999–2012 in the MHRB (Zyx). As mentioned above, it can be concluded that there was no lag between SDI and SPI; however, VCI with SPI had significant lag correlations in the short term in the UHRB and MHRB. Additionally, the VCI with SDI had a significant 1-month lag correlation in the short term in the MHRB (Zyx).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 13: Heat maps showing PC between SDI-1 lagged by 0–6 months and SPI accumulation periods of 1–48 months for different timescales at different periods in the UHRB (Ylx) and MHRB (Zyx). (a) (1965–1979, PC ≥ 0.15 or ≤ ‒0.16, p < 0.05), (b) (1980–2012, PC ≥ 0.14, p<0.01), and (c) (1965–2012, PC ≥ 0.11, p <0.01) in the UHRB (Ylx) and (d) (1965–1984, PC ≥ 0.16 or ≤‒0.11, p < 0.05), (e) (1985–2012, PC ≥ 0.12, p < 0.01), and (f) (1965–2012, PC ≥ 0.10, p < 0.05) in the MHRB (Zyx).

(a)

(b)

(c)

(d)

Figure 14: Heat maps showing PC between VCI lagged by 0–6 months and SPI (1-48)/SDI (1-12) for the period of 1999–2012 in the UHRM (Ylx). (a) PC ≥ 0.16, p < 0.05), (b) PC ≥ 0.18 or ≤ ‒0.15, p < 0.05) and in the MHRB (Zyx), (c) PC ≥ 0.15, p < 0.05), and (d) PC ≥ 0.16 or ≤ ‒0.18, p < 0.05).

4. Discussion

Our results indicated that MD were more serious in Pos-R than in Pre-R in the UHRB, there was little change in the MHRB, and MD were more serious in the UHRB than in the MHRB for duration and severity based on thresholds of −1.0, −1.5, and −2.0 for the total period. These results are inconsistent with [44], probably because the precipitation data types and SPI probability distribution function were different. In our study, precipitation data types were based on meteorological stations, and the SPI probability distribution function was a three-parameter log-logistic distribution. However, Ma et al. [44] used 0.5° resolution precipitation data and their SPI probability distribution function was an empirical Gringorten plotting position. The Pos-R period had a significant positive impact on the evolution of different timescales based on thresholds of −1.0, −1.5, and −2.0 in Ylx, which was consistent with Wu et al. [60], who showed that Pos-R period in the HD had reduced duration and severity compared with Pre-R period. However, there were more serious duration and severity, specifically at 6 and 12 months, than in the Pre-R period in Zyx. For the total period (Table S2), HD events were more frequent at Zyx than at Ylx for different timescales based on thresholds of −1.0, −1.5, and −2.0. These results were consistent with Ma et al. [44]. Figure 1(c) and Table 6 show that the main vegetation type in the UHRB was grassland, and there was an increasing trend to cultivated land for 2000–2010. The UHRB had fewer reservoirs and less arable land than the MHRB, and the drought duration and severity were more likely to be related to climate change. For example, Ma et al. [44], Qiu et al. [52], Gao and Zhang [62], Yang et al. [63], Cong et al. [64], and Cheng et al. [65] showed that an increase in streamflow was caused by climate change (increased precipitation) factors rather than human activities in the UHRB. At Zyx, the duration and severity of drought were greater than those in the Pre-R period. More extreme MD events occurred in 1985, 1992, 2000, and 2005 (i.e., Pos-R) [66], and the decrease in precipitation may have led to a decrease in streamflow and thus more serious HD events compared with the Pre-R period (Table S2). The area of cultivated land increased, and unused land and grassland decreased from 1975 to 2010. The total irrigation volume in the area increased by 0.57 × 10⁸ m³/year in 1990–2010 in the MHRB [54]. Along with the rapid increase in the number of pumping wells, groundwater irrigation has increased 2.11 × 10⁸ m³/year [54]. It is clear that human activities have changed the land surface and altered hydrological processes, including evapotranspiration, infiltration, surface runoff, and storage of water and; in this way, they affect the development of drought and thereby exacerbate drought [67]. Zhu et al. [68] also pointed out that human activities, such as the expansion of irrigation, rapid population growth, and socioeconomic development, have deeply modified hydrological processes in the MHRB. Therefore, the abovementioned climate change and human activities may affect changes in drought characteristics between the different periods in the UHRB and MHRB.

Table 6: Changes in the area of the primary types of land use and transition matrix and the primary types of land use changes in the UHRB and MHRB (10⁴ hm²).

Period
Land use types
Cultivated
Forest
Grassland
Water area
Built-up area
Unused land

UHRB
2000–2011
Cultivated
0.10
0.02
0.16
0.01
0.02
0.00

Forest
0.05
1.84
16.78
0.48
0.00
1.81

Grassland
0.29
0.65
42.13
1.27
0.03
6.13

Water area
0.04
0.06
0.96
0.91
0.01
0.63

Built land
0.03
0.00
0.04
0.00
0.02
0.00

Unused land
0.00
0.08
7.83
0.57
0.00
16.87

MRHB
1975–2010
Cultivated
14.41
0.04
0.26
0.05
0.28
0.15

Forest
0.42
1.20
0.21
0.01
0.01
0.12

Grassland
3.97
0.26
12.86
0.16
0.06
1.70

Water area
0.53
0.02
0.34
1.86
0.01
0.19

Built land
0.06
0.00
0.00
0.00
1.54
0.00

Unused land
1.84
0.38
1.46
0.14
0.22
62.76

MRHB
1975–1987

14.40
1.52
15.33
2.22
1.60
65.52

1987–1992
15.25
1.21
14.49
2.04
1.34
66.88

1992–2001
16.48
1.23
14.29
2.11
1.41
67.04

2000–2010
17.88
1.25
13.03
1.80
1.71
63.87

Previous research has examined the propagation time and lags of HD to MD in different regions of the world [4, 23, 43, 44]. In this study, we found that the PC at Ylx was larger than that at Zyx. The correlation with climate change and human activities changed from positive to negative at Zyx. Ma et al. [44] found that climate change was inclined to increase streamflow and propagation time, contributing from −57% to 63% in the UHRB, whereas human activities played a dominant role in water consumption with a contribution rate greater than −89% to further alter HD characteristics and propagation time in the MHRB. In this study, we found the MD propagation to HD was different in different catchments and had a significant seasonal discrepancy at different timescales for different periods. The reservoirs increased the propagation time of the hydrological variables to MD during the Pos-R period at Zyx, and the results were in line with Wu et al. [60]. The PC values between SDI and SPI in the Pos-R period were lower on a short-term scale (1–3 months) than those in the Pre-R period, whereas they were higher for the long-term timescale (>12 months) at Ylx and Zyx, which was also consistent with Wu et al. [60]. The propagation time in the Pos-R period increased in spring, decreased in autumn and winter, and did not change in summer at Ylx compared with the Pre-R period, while it increased in spring, summer, and autumn and decreased in winter at Zyx. The possible causes of the uneven spatial variations of precipitation are temperature and melting of snow and glaciers in the Qilian Mountains when the temperature is rising, which may lead to the differences in HD variations [44, 62]. The propagation time for tributaries would also vary in different periods in relation to natural and social environmental factors [4, 23, 69]. For instance, Deng et al. [70] demonstrated that the 33 tributaries in the MHRB no longer joined the mainstream after 1980s, and they gradually disappeared and formed independent irrigation oases. These may diverge in terms of different climate and human activity impacts. The lack of lag in the SDI and SPI was in line with Barker et al. [23] at Ylx and Zyx.

The high correlation between monthly NDVI and the MD index at different timescales is an indicator that can describe the impact of drought on vegetation, and the month of highest correlation means the greatest sensitivity of vegetation to drought [10, 41]. Our results exhibited a negative correlation on short timescales (≤9 months) between VCI and SPI, consistent with Li et al. [19]. Meanwhile increasing drought timescales (≥12 months) suggested positive correlations in the UHRB and MHRB (Figure 8(c)). Niu et al. [45] also found that variations in the Enhanced Vegetation Index (EVI) generally followed the evolution of soil moisture in the top 1 m of the soil profile. There were both positive and negative time lags. For example, during a 2001 drought event, it was observed that the Soil Moisture Anomaly Index reached a negative peak value. However, the EVI did not show a consistent response, indicating it was not affected by natural drought events because it was maintained by irrigation activities in cropland in the MHRB. The strongest PC values between VCI and SPI were for 6-month and 36-month timescales for the growing season in the UHRB and MHRB, respectively. These results were inconsistent with Vicente-Serrano et al. [41], who found that, in arid and humid regions, vegetation responded to MD at short timescales as well as in semiarid and subhumid regions at long timescales. The main reason may be the relationship with different vegetation types, anthropogenic activity, and climate change, although the UHRB is a semiarid region. In Figure 1(c), there is little coniferous forest, and the region is mainly grassland, shrub, and alpine vegetation. However, cultivated vegetation dominates in the MHRB, and there are 24 irrigation districts (Figure S5) and 71 reservoirs (four under construction) with thousands of canals and over 6000 pumping wells for irrigation [71]. The total irrigation volume in the area has increased in the MHRB [52]. Thus, crop improvements and water-saving irrigation technologies would increase water use efficiency and drought tolerance capacity [70]. Our findings also differed from Zhang et al. [18], who found significant annual high PC values in arid and semiarid regions, and the corresponding drought timescales were 3–6 months (January–December). This may be related to the different timescales (April–October). We further found a significant positive and high PC value (PC ≥ 0.53, ) between VCI and SPI on an intra-annual scale, implying that water availability is the critical factor for vegetation’s various spatiotemporal activities [52] in the UHRB and MHRB, especially in the MHRB (PC ≥ 0.70, ). Our results also revealed that, with an increasing drought timescale, the PC values between VCI and SPI were higher in the MHRB, specifically in August (the PC values fluctuated from 0.43 to 0.82 for 1–48 months) and September (the PC values fluctuated from 0.07 to 0.83 for 1–48 months). These results suggested that croplands can suffer from extreme and prolonged drought conditions [10] in the MHRB. Furthermore, different vegetation types have different propagation times for MD; the grassland and cultivated vegetation at a 3-month timescale and a 12-month timescale for the shrubland, coniferous forest, and broadleaf forest [8, 21]. Figures 14 and 6 demonstrate that alpine vegetation and grassland were poor (especially in April and May) and there was an extremely significant decrease for the growing season, whereas cultivated vegetation was good, and there was an extremely significant increase for the growing season in the UHRB and MHRB. Additionally, VCI with SPI had a significant lag correlation in the short term in the UHRB (6-month lag) and MHRB (4-month lag), and the VCI with SDI had a significant 1-month lag correlation in the MHRB, similar to Lin et al. [35]. Additionally, the correlation between a single timescale of the HD and VCI with the different timescales of the MD did not fully reflect their relationships. The reason is that the streamflow and vegetation propagation times with precipitation differ at different timescales [10, 18, 19, 21, 72, 73].

Our study had some limitations. The quantification of drought is a difficult task and drought is intrinsically multiscalar. There is no unique physical variable we can measure to quantify drought intensity. Vegetation propagation time to droughts is still an open scientific problem due to the complexity of droughts and limited knowledge of physiological processes. However, understanding the relationship between these mechanisms and the characteristics of droughts is crucial for improving our knowledge of vegetation vulnerability to climate fluctuations and climate change [41]. Meanwhile, with global warming, the propagation time of VCI and HD, VCI and MD, and HD and MD to climate change also has changed [5], which poses new challenges that need further research. Furthermore, increasing human activity has dramatically changed natural droughts; to manage drought effectively, we need to acknowledge that human influence is as integral to drought as natural climate variability. The complex interactions between natural and human processes need to be considered [67]. In addition, there are still some issues worthy of in-depth discussion. Often, people are considered as passive recipients at the end of the propagation cascade from MD via soil moisture drought to HD [74] but, in fact, people actively influence drought propagation [67]. Anthropogenic changes to the land surface alter hydrological processes, including evapotranspiration, infiltration, surface streamflow, and storage of water [67, 75]. Understanding how irrigation [76] (the effect of long-term irrigation in the HRB on the evolution of drought [77, 78]), land use/land cover change [79], and construction of reservoirs or dams [42, 63, 64, 80, 81] impact MD, HD, and VCI characteristics and MD propagation to HD and VCI need further research. Trnka et al. [82] presented a set of 60 priority questions and considered interdisciplinary characteristics to optimize future drought research, primarily covering the following drought-related topics: monitoring, impacts, drought forecasts, climatology, adaptation, and planning. Therefore, in future drought research, it is necessary to comprehensively consider priority questions and interdisciplinary characteristics, especially for irrigated agricultural areas in arid regions.

5. Conclusions

Drought affects land surface dynamics, and the strongest impact of drought on vegetation occurs in arid and semiarid areas. Therefore, understanding the evolutionary characteristics and propagation of different drought types may provide useful information about appropriate adaptation and mitigation strategies against the effects of drought on agricultural production and vegetation production and even alleviate the losses caused by drought. Our results highlight the following points:

(1)More drought events occurred at shorter accumulation periods based on thresholds = −1.0, −1.5, and −2.0 for MD and HD. Compared to the Pre-R period, the average peak, intensity, duration, and severity of HD decreased at Ylx, whereas the effects were more serious at Zyx. The PC values were lower at short-term scales (1–3 months) than at long-term scales for the Pos-R period at Ylx and Zyx. The PC values for MD and HD after the EWDP were higher than those before the EWDP in the UHRB and MHRB, especially in the MHRB. Human activities at different timescales (Pre-R period and Pos-R period and before EWDP and after EWDP) may affect the correlation between drought and the timescale of MD lag to HD.(2)The propagation time decreased 3 months at Ylx and increased 8 months at Zyx compared to the Pre-R period and between VCI and SPI at 6 months in the UHRB and 9 months in the MHRB. However, VCI with SDI at 6 months at Ylx and 3 months at Zyx did significantly affect the growing season. On an intra-annual scale, for Ylx, Zyx, and six tributaries, propagation times were different in the Pre-R and Pos-R periods, and there were seasonal differences between Pos-R and Pre-R periods.(3)The SDI and SPI showed no lag at Ylx and Zyx. However, VCI with MD had a significant lag correlation at the short-term scale in the UHRB (6-month lag) and MHRB (4-month lag), while the VCI with HD has significant 1-month lag correlation at Zyx.(4)Overall, the average peak/intensity/duration/severity of MD and HD are weakening and the propagation time from MD to HD is also reduced for Pos-R period compared with Pre-R period in the UHRB. Constructing reservoirs has had a positive effect, prolonging MD propagation to HD times, but the drought characteristics (average peak/intensity/duration/severity) have increased from the perspective of climatology. Therefore, positive drought prevention measures are necessary to consider the characteristics of propagation between MD and HD as well as seasonal differences in the MHRB.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest in this paper.

Acknowledgments

This research was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant nos. XDA19070502, XDA20100104, and XDA19040500), the National Natural Science Foundation of China (Grant nos. 41571516 and 41690144), and the Fundamental Research Funds for the Central Universities (Grant no. 2019jbkyjd013).

Supplementary Materials

Figure S1: boxplots showing MD peak and intensity in the UHRB ((a) and (c)) and in the MHRB ((b) and (d)) based on SPI using thresholds of −1, −1.5, and −2 for different timescales at different periods (the dot and midline in the boxplots are the mean and median; the upper and lower boundaries of box are the interquartile ranges). Figure S2: boxplots showing HD peak and intensity at Ylx ((a) and (c)) and at Zyx ((b) and (d)) based on SDI using thresholds of −1, −1.5, and −2 for different timescales at different periods (the dot and midline in the boxplots are the mean and median; the upper and lower boundaries of box are the interquartile ranges). Figure S3: the growing season spatial distribution of VCI during 2000–2001 in the UHRB and MHRB. Figure S4: heat map showing interannual Pearson’s correlation coefficient of SPI accumulation periods of 1–48 months with SDI-1 in the UHRB (Ylx) and MHRB (Zyx) (note: the x-axis represents the correlation between the SDI-1 and SPI at different timescales and the y-axis represents the 12 months of the year; the up arrow represents that the correlation coefficients are increased in the period after EWDP in the UHRB (Ylx) and MHRB (Zyx) compared to the period before EWDP, or vice versa); (a) (1965–1999, PC ≥ 0.35, ), (b) (2000–2012, PC ≥ 0.57, ) notes in the Ylx and (c) (1965–1999, PC ≥ 0.35 or ≤ −0.38, ), (d) (2000–2012, PC ≥ 0.56, ) notes in the Zyx. Figure S5: heat maps showing PC between VCI lagged by 0–6 months and SPI (1–48)/SDI (1–12) before (1965–1999) and after (2000–2012) EWDP; (a) (1965–1999, PC ≥ 0.10, ), (b) (2000–2012, PC ≥ 0.16, ) in the UHRM (Ylx) and (c) (1965–1999, PC ≥ 0.11 or PC ≤ −0.11, ), (d) (2000–2012, PC ≥ 0.17, ) in the MHRB (Zyx). Figure S6: reservoirs under construction and construction (a) irrigation system and area (b) in the MHRB. Table S1: identification of MD (threshold = −1, SPI-1) and HD (Zyx) characteristics (threshold = −1, SDI-1) based on run theory with comparison with historical drought records in MHRB. Table S2: MD and HD drought characteristics for different timescales for different threshold for prereservoir period, postreservoir period, and total period in UHRB (Ylx) and MHRB (Zyx). (Supplementary Materials)

References

S. E. Nicholson, C. J. Tucker, and M. B. Ba, “Desertification, drought, and surface vegetation: an example from the west african sahel,” Bulletin of the American Meteorological Society, vol. 79, no. 5, pp. 815–829, 1998.

A. T. DeGaetano, “A temporal comparison of drought impacts and responses in the New York city metropolitan area,” Climatic Change, vol. 42, no. 3, pp. 539–560, 1999.

A. K. Mishra and V. P. Singh, “A review of drought concepts,” Journal of Hydrology, vol. 391, no. 1-2, pp. 202–216, 2010.

S. Huang, P. Li, Q. Huang, G. Leng, B. Hou, and L. Ma, “The propagation from meteorological to hydrological drought and its potential influence factors,” Journal of Hydrology, vol. 547, pp. 184–195, 2017.

A. Dai, “Drought under global warming: a review,” Wiley Interdisciplinary Reviews: Climate Change, vol. 2, no. 1, pp. 45–65, 2011.

L. Wang and W. Chen, “A CMIP5 multimodel projection of future temperature, precipitation, and climatological drought in China,” International Journal of Climatology, vol. 34, no. 6, pp. 2059–2078, 2014.

G. Leng, Q. Tang, and S. Rayburg, “Climate change impacts on meteorological, agricultural and hydrological droughts in China,” Global and Planetary Change, vol. 126, pp. 23–34, 2015.

Z. Li, T. Zhou, X. Zhao et al., “Assessments of drought impacts on vegetation in China with the optimal time scales of the climatic drought index,” International Journal of Environmental Research and Public Health, vol. 12, no. 7, pp. 7615–7634, 2015.

J. Sheffield and E. F. Wood, “Projected changes in drought occurrence under future global warming from multi-model, multi-scenario, IPCC AR4 simulations,” Climate Dynamics, vol. 31, no. 1, pp. 79–105, 2008.

H.-J. Xu, X.-P. Wang, C.-Y. Zhao, and X.-M. Yang, “Diverse responses of vegetation growth to meteorological drought across climate zones and land biomes in northern China from 1981 to 2014,” Agricultural and Forest Meteorology, vol. 262, pp. 1–13, 2018.

M. Yu, Q. Li, M. J. Hayes, M. D. Svoboda, and R. R. Heim, “Are droughts becoming more frequent or severe in China based on the Standardized Precipitation Evapotranspiration Index: 1951–2010?” International Journal of Climatology, vol. 34, no. 3, pp. 545–558, 2014.

L. Zhang, J. Xiao, Y. Zhou, Y. Zheng, J. Li, and H. Xiao, “Drought events and their effects on vegetation productivity in China,” Ecosphere, vol. 7, no. 12, Article ID e01591, 2016.

R. R. Heim, “Review of twentieth-century drought indices used in the United States,” Bulletin of the American Meteorological Society, vol. 83, pp. 1149–1166, 2002.

T. J. Chang, “Investigation of precipitation droughts by use of kriging method,” Journal of Irrigation and Drainage Engineering, vol. 117, no. 6, pp. 935–943, 1991.

E. A. B. Eltahir, “Drought frequency analysis of annual rainfall series in central and western Sudan,” Hydrological Sciences Journal, vol. 37, no. 3, pp. 185–199, 1992.

J. A. Dracup, K. S. Lee, and E. G. Paulson Jr., “On the statistical characteristics of drought events,” Water Resources Research, vol. 16, no. 2, pp. 289–296, 1980.

S. M. Quiring and S. Ganesh, “Evaluating the utility of the Vegetation Condition Index (VCI) for monitoring meteorological drought in Texas,” Agricultural and Forest Meteorology, vol. 150, no. 3, pp. 330–339, 2010.

Q. Zhang, D. Kong, V. P. Singh, and P. Shi, “Response of vegetation to different time-scales drought across China: spatiotemporal patterns, causes and implications,” Global and Planetary Change, vol. 152, pp. 1–11, 2017.

C. Li, W. Leal Filho, J. Yin et al., “Assessing vegetation response to multi-time-scale drought across inner Mongolia plateau,” Journal of Cleaner Production, vol. 179, pp. 210–216, 2018.

Y. Zhang, J. Gao, L. Liu, Z. Wang, M. Ding, and X. Yang, “NDVI-based vegetation changes and their responses to climate change from 1982 to 2011: a case study in the Koshi River Basin in the middle Himalayas,” Global and Planetary Change, vol. 108, pp. 139–148, 2013.

A. Zhao, A. Zhang, S. Cao, X. Liu, J. Liu, and D. Cheng, “Responses of vegetation productivity to multi-scale drought in Loess Plateau, China,” CATENA, vol. 163, pp. 165–171, 2018.

B. Lloyd-Hughes, “The impracticality of a universal drought definition,” Theoretical and Applied Climatology, vol. 117, no. 3-4, pp. 607–611, 2014.

L. J. Barker, J. Hannaford, A. Chiverton, and C. Svensson, “From meteorological to hydrological drought using standardised indicators,” Hydrology and Earth System Sciences, vol. 20, no. 6, pp. 2483–2505, 2016.

B. T. Mckee, J. Nolan, and J. Kleist, “The relationship of drought frequency and duration to time scales,” in Proceedings of the 8th Conference on Applied Climatology, vol. 17–22, pp. 179–184, Anaheim, CA, USA, January 1993.

D. Wang, M. Hejazi, X. Cai, and A. J. Valocchi, “Climate change impact on meteorological, agricultural, and hydrological drought in central Illinois,” Water Resources Research, vol. 47, no. 9, 2011.

A. AghaKouchak and N. Nakhjiri, “A near real-time satellite-based global drought climate data record,” Environmental Research Letters, vol. 7, no. 4, p. 044037, 2012.

J. H. Stagge, L. M. Tallaksen, L. Gudmundsson, A. F. Van Loon, and K. Stahl, “Candidate distributions for climatological drought indices (SPI and SPEI),” International Journal of Climatology, vol. 35, no. 13, pp. 4027–4040, 2015.

I. Nalbantis and G. Tsakiris, “Assessment of hydrological drought revisited,” Water Resources Management, vol. 23, no. 5, pp. 881–897, 2009.

S. Li, L. Xiong, L. Dong, and J. Zhang, “Effects of the three gorges reservoir on the hydrological droughts at the downstream yichang station during 2003–2011,” Hydrological Processes, vol. 27, no. 26, pp. 3981–3993, 2013.

H. Tabari, J. Nikbakht, and P. Hosseinzadeh Talaee, “Hydrological drought assessment in northwestern Iran based on streamflow drought index (SDI),” Water Resources Management, vol. 27, no. 1, pp. 137–151, 2013.

E. Rimkus, E. Stonevičius, V. Korneev, J. Kažys, G. Valiuškevičius, and A. Pakhomau, “Dynamics of meteorological and hydrological droughts in the Neman river basin,” Environmental Research Letters, vol. 8, no. 4, p. 045014, 2013.

T. Fischer, M. Gemmer, B. Su, and T. Scholten, “Hydrological long-term dry and wet periods in the Xijiang River basin, South China,” Hydrology and Earth System Sciences, vol. 17, no. 1, pp. 135–148, 2013.

X. Hong, S. Guo, Y. Zhou, and L. Xiong, “Uncertainties in assessing hydrological drought using streamflow drought index for the upper Yangtze River basin,” Stochastic Environmental Research and Risk Assessment, vol. 29, no. 4, pp. 1235–1247, 2015.

A.-A. Belal, H. R. El-Ramady, E. S. Mohamed, and A. M. Saleh, “Drought risk assessment using remote sensing and GIS techniques,” Arabian Journal of Geosciences, vol. 7, no. 1, pp. 35–53, 2014.

Q. Lin, Z. Wu, V. P. Singh, S. H. R. Sadeghi, H. He, and G. Lu, “Correlation between hydrological drought, climatic factors, reservoir operation, and vegetation cover in the Xijiang Basin, South China,” Journal of Hydrology, vol. 549, pp. 512–524, 2017.

I. Cordery and M. McCall, “A model for forecasting drought from teleconnections,” Water Resources Research, vol. 36, no. 3, pp. 763–768, 2000.

S. Huang, J. Chang, Q. Huang, and Y. Chen, “Spatio-temporal changes and frequency analysis of drought in the wei river basin, China,” Water Resources Management, vol. 28, no. 10, pp. 3095–3110, 2014.

A. Mondal and P. P. Mujumdar, “Return levels of hydrologic droughts under climate change,” Advances in Water Resources, vol. 75, pp. 67–79, 2015.

J. Niu, J. Chen, and L. Sun, “Exploration of drought evolution using numerical simulations over the Xijiang (West River) basin in South China,” Journal of Hydrology, vol. 526, pp. 68–77, 2015.

J. Lorenzo-Lacruz, S. Vicente-Serrano, J. González-Hidalgo, J. López-Moreno, and N. Cortesi, “Hydrological drought response to meteorological drought in the Iberian Peninsula,” Climate Research, vol. 58, no. 2, pp. 117–131, 2013.

S. M. Vicente-Serrano, C. Gouveia, J. J. Camarero et al., “Response of vegetation to drought time-scales across global land biomes,” Proceedings of the National Academy of Sciences, vol. 110, no. 1, pp. 52–57, 2013.

A. F. Van Loon, M. H. J. Van Huijgevoort, and H. A. J. Van Lanen, “Evaluation of drought propagation in an ensemble mean of large-scale hydrological models,” Hydrology and Earth System Sciences, vol. 16, no. 11, pp. 4057–4078, 2012.

M. Peña-Gallardo, S. M. Vicente-Serrano, J. Hannaford et al., “Complex influences of meteorological drought time-scales on hydrological droughts in natural basins of the contiguous Unites States,” Journal of Hydrology, vol. 568, pp. 611–625, 2019.

F. Ma, L. Luo, A. Ye, and Q. Duan, “Drought characteristics and propagation in the semiarid Heihe river basin in northwestern China,” Journal of Hydrometeorology, vol. 20, no. 1, pp. 59–77, 2019.

J. Niu, S. Kang, X. Zhang, and J. Fu, “Vulnerability analysis based on drought and vegetation dynamics,” Ecological Indicators, vol. 105, pp. 329–336, 2019.

R. P. Pandey and K. S. Ramasastri, “Relationship between the common climatic parameters and average drought frequency,” Hydrological Processes, vol. 15, no. 6, pp. 1019–1032, 2001.

X. Li, G. Cheng, Y. Ge et al., “Hydrological cycle in the Heihe river basin and its implication for water resource management in endorheic basins,” Journal of Geophysical Research: Atmospheres, vol. 123, no. 2, pp. 890–914, 2018.

N. Zhao, T. Yue, C. Chen, M. Zhao, and Z. Fan, “An improved statistical downscaling scheme of Tropical Rainfall Measuring Mission precipitation in the Heihe River basin, China,” International Journal of Climatology, vol. 38, no. 8, pp. 3309–3322, 2018.

A. Zhang, C. Zheng, S. Wang, and Y. Yao, “Analysis of streamflow variations in the Heihe River Basin, northwest China: trends, abrupt changes, driving factors and ecological influences,” Journal of Hydrology: Regional Studies, vol. 3, pp. 106–124, 2015.

J. Niu, X.-G. Zhu, M. A. J. Parry et al., “Environmental burdens of groundwater extraction for irrigation over an inland river basin in Northwest China,” Journal of Cleaner Production, vol. 222, pp. 182–192, 2019.

Y. Ge, X. Li, C. Huang, and Z. Nan, “A decision support system for irrigation water allocation along the middle reaches of the Heihe river basin, northwest China,” Environmental Modelling & Software, vol. 47, pp. 182–192, 2013.

L. Qiu, D. Peng, Z. Xu, and W. Liu, “Identification of the impacts of climate changes and human activities on runoff in the upper and middle reaches of the Heihe River basin, China,” Journal of Water and Climate Change, vol. 7, no. 1, pp. 251–262, 2015.

X. Deng and C. Zhao, “Identification of water scarcity and providing solutions for adapting to climate changes in the Heihe river basin of China,” Advances in Meteorology, vol. 2015, Article ID 279173, 13 pages, 2015.

M. Zhang, S. Wang, B. Fu, G. Gao, and Q. Shen, “Ecological effects and potential risks of the water diversion project in the Heihe River Basin,” Science of the Total Environment, vol. 619-620, pp. 794–803, 2018.

A. A. Paulo, L. S. Pereira, and P. G. Matias, “Analysis of local and regional droughts in southern Portugal using the theory of runs and the standardised precipitation index,” in Tools for Drought Mitigation in Mediterranean Regions, G. Rossi, A. Cancelliere, L. S. Pereira, T. Oweis, M. Shatanawi, and A. Zairi, Eds., Springer, Dordrecht, Netherlands, 2003.

F. Kogan and J. Sullivan, “Development of global drought-watch system using NOAA/AVHRR data,” Advances in Space Research, vol. 13, no. 5, pp. 219–222, 1993.

J. Spinoni, T. Antofie, P. Barbosa et al., “An overview of drought events in the Carpathian Region in 1961–2010,” Advances in Science and Research, vol. 10, no. 1, pp. 21–32, 2013.

V. M. Yevjevich, An Objective Approach to Definitions and Investigations of Continental Hydrologic Droughts, Colorado State University, Fort Collins, CO, USA, 1967.

Y. L. Zhou, P. Zhou, J. L. Jin, and J. Li, “Establishment of hydrological drought index on sources of regional water supply and its application to drought frequency analysis for Kunming,” Journal of Hydraulic Engineering, vol. 45, pp. 1038–1047, 2014.

J. Wu, Z. Liu, H. Yao et al., “Impacts of reservoir operations on multi-scale correlations between hydrological drought and meteorological drought,” Journal of Hydrology, vol. 563, pp. 726–736, 2018.

S. Mitra and P. Srivastava, “Spatiotemporal variability of meteorological droughts in southeastern USA,” Natural Hazards, vol. 86, no. 3, pp. 1007–1038, 2017.

L. Gao and Y. Zhang, “Spatio-temporal variation of hydrological drought under climate change during the period 1960–2013 in the Hexi Corridor, China,” Journal of Arid Land, vol. 8, no. 2, pp. 157–171, 2016.

L. Yang, Q. Feng, Z. Yin et al., “Identifying separate impacts of climate and land use/cover change on hydrological processes in upper stream of Heihe River, Northwest China,” Hydrological Processes, vol. 31, no. 5, pp. 1100–1112, 2017.

Z. Cong, M. Shahid, D. Zhang, H. Lei, and D. Yang, “Attribution of runoff change in the alpine basin: a case study of the Heihe Upstream Basin, China,” Hydrological Sciences Journal, vol. 62, no. 6, pp. 1013–1028, 2017.

Q. Cheng, X. Zuo, F. Zhong, L. Gao, and S. Xiao, “Runoff variation characteristics, association with large-scale circulation and dominant causes in the Heihe River Basin, Northwest China,” Science of the Total Environment, vol. 688, pp. 361–379, 2019.

K. G. Wen, China Meteorological Disaster (Gansu Volume), China Meteorological, Beijing, China, 2006, In Chinese.

A. F. Van Loon, T. Gleeson, J. Clark et al., “Drought in the anthropocene,” Nature Geoscience, vol. 9, no. 2, pp. 89–91, 2016.

G. Zhu, Y. Su, C. Huang, Q. Feng, and Z. Liu, “Hydrogeochemical processes in the groundwater environment of Heihe River Basin, northwest China,” Environmental Earth Sciences, vol. 60, no. 1, pp. 139–153, 2010.

J. I. López-Moreno, S. M. Vicente-Serrano, J. Zabalza et al., “Hydrological response to climate variability at different time scales: a study in the Ebro basin,” Journal of Hydrology, vol. 477, pp. 175–188, 2013.

X.-P. Deng, L. Shan, H. Zhang, and N. C. Turner, “Improving agricultural water use efficiency in arid and semiarid areas of China,” Agricultural Water Management, vol. 80, no. 1–3, pp. 23–40, 2006.

X. Xu, Y. Jiang, M. Liu, Q. Huang, and G. Huang, “Modeling and assessing agro-hydrological processes and irrigation water saving in the middle Heihe River basin,” Agricultural Water Management, vol. 211, pp. 152–164, 2019.

J. Niu, J. Chen, K. Wang, and B. Sivakumar, “Multi-scale streamflow variability responses to precipitation over the headwater catchments in southern China,” Journal of Hydrology, vol. 551, pp. 14–28, 2017.

H. Liu, M. Zhang, Z. Lin, and X. Xu, “Spatial heterogeneity of the relationship between vegetation dynamics and climate change and their driving forces at multiple time scales in Southwest China,” Agricultural and Forest Meteorology, vol. 256-257, pp. 10–21, 2018.

R. Orth and G. Destouni, “Drought reduces blue-water fluxes more strongly than green-water fluxes in Europe,” Nature Communications, vol. 9, no. 1, p. 3602, 2018.

N. S. Diffenbaugh, D. L. Swain, and D. Touma, “Anthropogenic warming has increased drought risk in California,” Proceedings of the National Academy of Sciences, vol. 112, no. 13, pp. 3931–3936, 2015.

P. Droogers, W. G. M. Bastiaanssen, M. Beyazgül, Y. Kayam, G. W. Kite, and H. Murray-Rust, “Distributed agro-hydrological modeling of an irrigation system in western Turkey,” Agricultural Water Management, vol. 43, no. 2, pp. 183–202, 2000.

Y. Chen, J. Niu, S. Kang, and X. Zhang, “Effects of irrigation on water and energy balances in the Heihe River basin using VIC model under different irrigation scenarios,” Science of The Total Environment, vol. 645, pp. 1183–1193, 2018.

J. Niu, Q. Liu, S. Kang, and X. Zhang, “The response of crop water productivity to climatic variation in the upper-middle reaches of the Heihe River basin, Northwest China,” Journal of Hydrology, vol. 563, pp. 909–926, 2018.

B. I. Cook, R. L. Miller, and R. Seager, “Amplification of the North American “Dust Bowl” drought through human-induced land degradation,” Proceedings of the National Academy of Sciences, vol. 106, no. 13, pp. 4997–5001, 2009.

Q. Huang, Z. Sun, C. Opp, T. Lotz, J. Jiang, and X. Lai, “Hydrological drought at dongting lake: its detection, characterization, and challenges associated with three gorges dam in central yangtze, China,” Water Resources Management, vol. 28, no. 15, pp. 5377–5388, 2014.

S. Rangecroft, A. F. Van Loon, H. Maureira, K. Verbist, and D. M. Hannah, “Multi-method assessment of reservoir effects on hydrological droughts in an arid region,” Earth System Dynamics Discussion, pp. 1–32, 2016.

M. Trnka, M. Hayes, F. Jurečka et al., “Priority questions in multidisciplinary drought research,” Climate Research, vol. 75, no. 3, pp. 241–260, 2018.