The frequency and intensity of drought are increasing
dramatically with global warming. However, few studies have characterized
drought in terms of its impacts on ecosystem services, the mechanisms
through which ecosystems support life. As a result, little is known about
the implications of increased drought for resource management. This case
study characterizes drought by linking climate anomalies with changes in the
precipitation–runoff relationship (PRR) on the Loess Plateau of China, a
water-limited region where ongoing revegetation makes drought a major
concern. We analyzed drought events with drought durations
Drought is a complex and recurrent climate-related phenomenon that has far-reaching adverse impacts on agriculture, water resources, the environment and human life (Ghulam et al., 2007; Goddard et al., 2001). Few types of extreme events have been as economically and ecologically disruptive as droughts over the past half century (Dai, 2011). For example, multi-year droughts across the globe have triggered an increase in tree mortality that is linked to climate change (Allen et al., 2010). On regional scales, during 1998–2001, one-third of the territory of Iran was affected by the most severe drought in the history of the country. More than half of the country's population faced food shortages and a lack of drinking-water supplies (Raziei et al., 2009). The annual economic losses caused by droughts in the United States are estimated as being up to USD 6–8 billion (Wilhite, 2000). As a result of these impacts, increasing numbers of studies now focus on characterizing drought events, including drought identification and frequency analysis, and the resource management actions required by droughts (Michele et al., 2013).
Many drought indices have been developed to monitor the evolution of drought events on regional and global scales (Yan et al., 2016). Meteorological drought can be identified from anomalies in precipitation data, as assessed by the Palmer drought severity index (PDSI; Palmer, 1965), the standardized precipitation index (SPI; McKee et al., 1993), the precipitation anomaly (Moron, 1994) and the Bhalme–Mooley index (Bhalme and Mooley, 1980). Other researchers have recognized the necessity of developing drought indicators that reflect the causes and impacts of drought, and the resulting drought indices incorporate many different physical, biological and socioeconomic variables (Sheffield et al., 2004). For example, the soil moisture indicator (Xia et al., 2014), the crop drought indicator (Duff et al., 1997) and the crop water demand indicators are used to identify agricultural drought events, which are periods that feature dry soil conditions and result from below-average precipitation, intense but less frequent rain events or above-normal evaporation. All of these factors lead to reductions in crop production and plant growth. Hydrological droughts (Herbst et al., 1966; Mohan and Rangacharya, 1991) are usually identified using run-length theory and are characterized as periods when water demand exceeds supply because of long-term precipitation shortages. In addition, socioeconomic drought reflects the undesirable social and economic impacts induced by the other types of droughts mentioned above. All of these drought indicators employ single variables to quantify drought events. However, drought conditions are associated with multiple variables; thus, no single drought indicator provides a satisfactory characterization of the complex conditions associated with droughts or their broad impacts (Hao and Singh, 2015). At the same time, because ecosystem services are closely related to the living conditions experienced by human beings (Burkhard et al., 2014; Zheng et al., 2016), these drought indicators also cannot reflect the impacts of droughts on ecosystems and provide relatively little information to policy makers and resource managers.
The precipitation–runoff relationship (PRR) is an important issue in engineering hydrology, water resource planning and management, and watershed system evolution (Guo et al., 2016; Nourani et al., 2015). Current studies focus on changes in the PRR and the processes by which these changes occur. Charlier et al. (2015) applied the wavelet transform method to detect changes in the PRR. Sun et al. (2016) assessed hydrologic trends in urban catchments using a conceptual urban precipitation–runoff model, which is beneficial for storm water management and planning. There is no doubt that variations in the PRR reflect the actual integrated volume of precipitation and runoff. The ratio between the annual runoff and the annual precipitation, the so-called “runoff ratio” (Savenije, 1996; Feng et al., 2016b) is frequently employed to quantify the ecosystem service of water yield. This ecosystem service is of major concern in water-limited areas because it represents the water resource available to human beings. Studies have shown that protracted drought may affect runoff generation and cause changes in the PRR (Petrone et al., 2010; Saft et al., 2015). Therefore, the shift in the PRR caused by an extended drought will eventually have an adverse effect on the ecosystem service of water yield. Based on this perspective, we want to characterize drought events in terms of the changes in the PRR that they produce and investigate the characteristics and frequency of these dry events, thus contributing to further optimizations of ecosystem management.
Drought stress is the main environmental factor limiting terrestrial ecosystem productivity in arid and semi-arid areas (Boyer, 1982). The vegetation restoration programme in China represents the largest investment that has been made to restore the ecosystem in this developing country. Given the limited water resources on the Loess Plateau, the sustainability of vegetation restoration is a major concern of scientific research and policy makers there (Feng et al., 2016b). As a typical coupled human–natural system on the Loess Plateau of China, the ecosystem service of water yield is of concern for both the sustainability of the revegetation programme and the objective of raising the level of economic development (Wang et al., 2015; Feng et al., 2016a). In recent years, the climate of the Loess Plateau has become warmer and drier (Lü et al., 2014) as global temperatures have increased, a change that will eventually affect the balance between the water supply and demand of the region. Thus, analyzing drought characteristics in terms of the changes in the PRR that occur in response to multi-year dry periods is of great importance in estimating the effects of drought and the ecological reconstruction of the Loess Plateau as a whole.
The objectives of our study are to characterize droughts on the Loess Plateau. We first link climate conditions and the PRR to define drought. We then simulate the drought frequency and return period of drought events of different magnitudes. Finally, we discuss the policy implications of our study for ecosystem management.
Map of the study area.
Watershed characteristics and hydrological data for the study period 1961–1999.
ID, watershed identification number; PPT, annual precipitation.
The Loess Plateau (Fig. 1) is located between 100
In this study, we examine the 13 primary watersheds of the Loess
Plateau (Fig. 1), which make up approximately 35 % of the area of the Loess
Plateau. The river runoff from these 13 watersheds contributes 65 %
of the discharge into the middle reach of the Yellow River. The attributes
of each basin are shown in Table 1. Runoff and meteorological data for these
13 watersheds covering the period from 1961 to 1999 were obtained from the
Yellow River Conservancy Commission (
In this study, we define drought based on annual precipitation for two aspects. On the one hand, the amount of precipitation is the most important climatic control of drought conditions (Mishra and Singh, 2010). Moreover, because we are interested in determining whether the runoff response differs for multi-year droughts, we do not consider runoff in identifying drought events.
We first calculate the precipitation anomaly (PA) values in the studied
watersheds on the Loess Plateau. The time series of anomaly values are
divided by the mean annual precipitation and smoothed with a 3-year moving
average. Positive PA values indicate that the observed precipitation is
higher than the median. On the other hand, negative PA values indicate that
the observed precipitation is below the median and imply the possible
occurrence of a drought. Each drought event is characterized in terms of its
duration and severity. Studies have shown that the drought events with
shorter durations but greater intensities or lower intensities but greater
durations cause serious water-supply and other drought-related problems
(Shiau, 2006; Naresh et al., 2009). Therefore, the basic rules for
identifying drought events in this study are (1) a PA value for a single
year of
In this study, the cumulative PA values during each drought period are used
to measure drought severity (for convenience, drought severity is multiplied
by
We use the Kolmogorov–Smirnov (K–S; Massey, 1951) test to determine whether
annual precipitation and runoff data follow a roughly normal distribution. A
Box–Cox transformation is applied to those data that are not normally
distributed (Box and Cox, 1964). After identifying the major drought events,
we examine whether the change in the PRR is statistically significant
compared to the historic record using Student's
Drought can be characterized by multiple variables, such as duration,
severity and spatial extent (Steinemann and Cavalcanti, 2006; Hayes et al.,
2012), but how to determine the joint distribution between these variables
remains an important issue. Here, we use the copula function (Shiau, 2006).
We construct a joint distribution function using two main drought
characteristics, drought duration and severity. If the marginal distribution
functions of drought duration (
The maximum likelihood (ML), inference functions for margins (IFM) and
canonical maximum likelihood (CML) methods are commonly used in parameter
estimation (Mirabbasi et al., 2012; Lee et al., 2013). Here, we use ML and
IFM to estimate the parameters of the marginal distribution functions and
the joint copula functions of drought duration and severity, respectively.
We employ seven common distributions to describe the univariate probability
distributions as the candidate margins for drought duration and severity.
These distributions are the exponential, gamma, log-normal, extreme value,
generalized extreme value, Poisson and Weibull distributions. The K–S test
is used to establish the optimal marginal distribution function. The ML
method is used to estimate the distribution parameters of drought duration
(
Common two-dimensional copula function families.
Five commonly used two-dimensional functions are constructed using the
marginal distribution function of drought duration and drought severity
(Table 2). The goodness-of-fit test is performed by calculating the squared
Euclidean distance (SED) between the theoretical copula and the empirical
copula (Berg, 2009). The empirical copula and SED are respectively defined as
In addition to the SED method, the root mean square error (RMSE) and the Akaike information criterion (AIC) are adopted to further evaluate the fitted copula.
Of the 39 years that make up the study period, the precipitation anomaly is negative in 17 years. Based on the drought identification method developed in this study, seven dry periods are identified (Fig. 2), including both major dry periods and single-year dry periods. We find that the periods of 1970–1974, 1979–1983 and 1991–1999 all represent droughts with durations longer than 5 years. The corresponding drought severities are 0.51, 0.32 and 0.61, and the average precipitation values of these events are 231, 320 and 183 mm, respectively.
Time series with identified drought periods on the Loess Plateau during 1961–1999.
Parameters and goodness-of-fit values of the marginal distributions. The best fitted marginal distributions of duration and severity are marked with bold font.
Applying the K–S test to the precipitation–runoff data from 1961 to 1999
shows that the precipitation–runoff data in this time series approximate a
normal distribution, providing the premise for a linear relationship between
precipitation and runoff. The changes in the PRR during the three major
drought periods on the Loess Plateau show that the PRR may change
significantly during these drought periods (Fig. 3). No significant changes
can be identified during the drought periods of 1970–1974 (
The goodness-of-fit values of the copula functions. The optimal of copula function in each watershed is marked with bold font.
Table 3 lists the estimated parameters and the results of goodness-of-fit
tests for the marginal functions. We find that not all of the distributions
pass the K–S test at the 95 % (
Precipitation–runoff relationships during drought periods:
The joint probability distribution of drought duration and
drought severity and the joint return period:
Figure 4a and b show that the joint probability increases with increasing
drought duration and drought severity, regardless of the three-dimensional
joint cumulative probability or the contour lines. When the drought severity
ranges from 0.1 to 0.35 (which corresponds to precipitation values that
range from 306 to 423 mm), the contours are nearly vertical, and the
drought duration varies greatly (ranging from 1 to 9 years). When the
drought duration is less than 5 years, the cumulative probability of the
bivariate increases rapidly as drought severity increases. In contrast, when
the drought duration exceeds 5 years, the rate of increase in the joint
probability decreases as drought severity increases. When the drought
severity is not more than 0.4 (precipitation
Figure 4c and d show that, as drought duration and drought severity increase, the joint return period also shows an increasing trend. When the duration and severity of drought events reach a maximum during the study period, the joint return period of such drought events on the Loess Plateau is close to 22 years. For the three major drought events that occurred on the Loess Plateau between 1961 and 1999, the return period is about 5.74 years for drought durations of 5 years or drought severity of up to 0.51 (1970–1974). The return period reaches 4.32 years for the drought that occurred in 1979–1983, which displayed a duration of 5 years and a severity of 0.39. The return period of the most severe drought event in the study period, which had a duration of 9 years and a severity of 0.61, is 18.31 years. The return period of the drought period that occurred in 1970–1974 is approximately 5.74 years, given the corresponding drought characteristics. Therefore, the next drought event similar to the drought period that took place from 1970 to 1974 occurred around 1980. During 1979–1983, the drought duration reached 5 years, which is close to the estimated return period.
The annual precipitation–runoff scatter plot for each watershed.
Prolonged multi-year drought causes significant damages in natural
environments. Figure 5 demonstrates the range of changes in the PRR under
sustained precipitation decreases. According to the direction of change, the
dry period regression line is mainly located above or below the overall
regression line, and the PRR in the 13 studied watersheds exhibits no
significant change when the regression line of the dry period is above the
total regression line. Of the 15 cases in which dry events fell under the
overall regression line in the 13 watersheds from 1961 to 1999, significant
changes in the PRR (
Spatial distribution of watersheds with and without significant changes in the PRR during the drought periods.
Spatial distribution of the return period of drought
events in watersheds:
Trend in the annual runoff ratio from 1982 to 1999.
Anthropogenic factors of runoff ratio change in the 1980s vs. the 1990s (Terra, Dam, Tree and Pasture are the changes in the percentage area of terraces, check-dams, tree plantation and natural pastures, respectively. LAI is change in the GIMMS LAI for each catchment).
There were no significant changes in the PRR in 5 of the 13 watersheds during the study period. Figure 6 demonstrates that there is no geographical pattern in the spatial distribution of watersheds with and without significant change in the PRR. Compared with the annual average precipitation in separate watersheds during 1961–1999, the watersheds where no significant changes in the PRR occurred (Kuye, Dali, Qingjian, Yanhe and Jinghe) received greater amounts of precipitation. Thus, we conclude that the PRR responds strongly to protracted drought, and the persistence of drought conditions over many years is more likely to cause significant changes in the PRR in basins that receive less precipitation.
The return period varies with drought duration and severity. Therefore, the three major dry events that occurred on the Loess Plateau from 1961 to 1999 were selected to enable a study of the characteristics of the spatial distribution of the joint return period. Figure 7a, b and c show that, although the duration and severity of drought differ among these three major drought events, the spatial distribution of the drought return period is consistent. At the same time, the drought return period also differs between watersheds, and the spatially heterogeneous characteristics of the different catchments on the Loess Plateau can also be seen from the drought return period. In terms of spatial distribution, the return periods of drought events are longer in the southern and eastern watersheds of Jinghe, Beiluo, Xinshui and Fenhe, indicating that the frequency of drought events in these watersheds is relatively low. On the other hand, the watersheds to the north and west, such as Huangfu, Kuye, Tuwei and Weihe, display shorter drought return periods. Thus, these regions are more susceptible to drought under the same drought characteristics. We analyze the return periods corresponding to watersheds with significant changes in the relationship between precipitation and runoff. As shown in Fig. 7d, there are obvious differences in the return periods in these watersheds, which experience different drought characteristics. The return period in these watersheds is particularly important in ecosystem management.
Since the actual occurrence of drought in a region is complicated, and the
practical significance of various drought indicators is different, the
choice of drought definition is an important part of studying the process of
drought occurrence and development. Prolonged multi-year drought events
cause significant damages to both the natural environment and the
development of human societies (Belal et al., 2014). To ensure that the dry
periods would be sufficiently long and severe, Saft et al. (2015) used only
dry periods with drought duration
Ecosystem services represent a key concept for policy makers, and the variability of the relationship between precipitation and runoff is vital in the study of the ecosystem service of water yield. The PRR can be influenced by factors other than climate conditions, such as land use, the diversion of surface water, irrigation schemes, groundwater abstraction and the storage of water in catchments (Farley et al., 2005; Brown et al., 2005; Zheng et al., 2009). However, the catchments examined in our study lie in the part of the Loess Plateau with the greatest relief, and the Mu Us Desert is located in the northwest and the Weihe Plain is located in the southeast. The vegetation in the catchments is mostly rain-fed; thus, irrigation schemes can be neglected in the study area (Feng et al., 2016a). The thickness of the loess within the catchments is greater than 100 m (Derbyshire et al., 1998), and the groundwater is minimally impacted by the surface eco-hydrological processes; thus, groundwater recharge and groundwater discharge are not considered in the study area. In addition, Bouwer et al. (2006) concluded that increasing water consumption for irrigation and the degree of runoff variability caused by hydropower is 3 times higher than the variations in runoff under climate change in a densely populated region in the main agricultural irrigation area in India. The water consumed by the local communities on the Loess Plateau is fed to the residential areas, which are mainly located in the flat areas at the outlets of the catchments (the distribution of residential areas is shown in Fig. 1). Moreover, the population shows a tendency to move from the catchment area to the major cities, which are located along the mainstream of the river basin in the Weihe Plain (i.e., Baoji and Xi'an, these cities contain 57.35 % of the population within the study area) because of the accelerated urbanization that has taken place in this area since the 1980s (Hu et al., 2001). Thus, the water consumed by the local communities does not have a major effect on runoff in our study area. Finally, any diversions of surface water and water storage are found in the residential areas at the outlets of these catchments. Therefore, the impacts of these factors are not included in the PRR.
Instead, soil conservation measures, including the construction of terraces
and the construction of sediment-trapping dams, have been implemented on the
Loess Plateau since the 1950s (Wang et al., 2015), and these anthropogenic
factors may change water yields (Wang et al., 2009; Shi et al., 2013; Chang
et al., 2015). We used the partial correlation method to isolate the impacts
of anthropogenic activities from climatic factors. For the entire period of
1982–1999, the runoff ratio displayed a decreasing trend (Fig. 8). Terrace
construction played an important role in the reduction in the runoff ratio
from the 1980s to the 1990s (
The PRR responds to multi-year droughts in two different ways. The PRR either displays no significant change or a significant downward change. Saft et al. (2015) have shown that a significant upward trend in the PRR occurred during a period of drought in watersheds in southeastern Australia, but the probability of such effects is small. It is also not clear whether these observations reflect a real phenomenon or just sampling fluctuations. Analyzing the changes in the PRR during the three major dry periods on the Loess Plateau shows that the annual runoff tends to decrease gradually. Therefore, we believe that the occurrence of dry events will aggravate the reduction of runoff, which could lead to significant changes in the PRR. This study shows that sustained precipitation changes also have the capacity to transform the water balance of watersheds in water-limited areas.
The drought regressions fell under the total regression line 15 times in the
13 watersheds, indicating the occurrence of both significant and
non-significant changes during the study period. During drought periods, the
runoff decreased as expected when the precipitation decreased, which
decreases the runoff ratio and affects the ecosystem service of water yield.
Based on the analysis of the duration and severity of the 15 dry events
(Fig. 10a), we conclude that the occurrence of a drought event is more
likely to result in a significant decreasing trend in the PRR when the
drought duration is
In late 1999, China implemented a massive revegetation program in the name
of ecological restoration. This restoration involved returning croplands on
steep slopes (> 25
Drought events affect ecosystem productivity and reduce the carbon sink
capacity of terrestrial ecosystems (Ciais et al., 2005; Tian et al., 1998).
As an important ecosystem structural parameter, the leaf area index (LAI)
characterizes the physiologically functioning surface area through which
energy, mass (e.g., water and CO
The results of this study provide guidance for ecosystem management policy. In addition, the specific measures used in adapting to drought should also be improved. For example, drought-resistant crops, such as millet and maize, should be chosen. Farmers can also increase the soil water storage capacity by reclaiming level terraces and practicing contour strip farming (Panagos et al., 2015) and ridge-furrow cropping (Gan et al., 2013). Vegetation affects the conversion of surface energy, water, momentum and biochemicals through its physical and physiological processes, thus affecting atmospheric conditions (Bonan, 2008) and further changing regional precipitation and hydrological processes (Ellison et al., 2012). During afforestation, managers should select tree species that consume less water, reduce the density of planting, and consider the optimal distribution of woodlands, shrubs and grasslands to cope with the effects of warmer and drier conditions. For example, fish-scale pits (Wang et al., 2014) can be used to enhance the infiltration of atmospheric precipitation and soil moisture. Especially in areas with steep slopes, the use of rainwater harvesting measures, such as fish-scale pits, will boost the survival of trees and avoid drying of the soil.
Trend in the LAI from 1982 to 1999.
This article characterizes drought by linking climate anomalies with changes
in the PRR. We found that multi-year drought caused the PRR to display a
significant decreasing trend (
Our analysis revealed great spatial variability in drought across the Loess Plateau. We chose the Frank copula function as the optimal joint distribution function of drought duration and severity. The results demonstrated that the joint probability and return period gradually increased with increasing drought duration and severity. At the same time, the spatially heterogeneous characteristics of different watersheds on the Loess Plateau can be seen from the spatial distribution of the drought return period. By analyzing the annual LAI over the Loess Plateau during 1982–1999, we found that the LAI significantly decreases in drought periods compared to normal periods. Long-term drought is also a key factor in considering the influence of LAI trends in the future.
These results should lead to better water regulation and more effective strategies of ecosystem management. We can consider different plant species and the spatial variability in drought events to maximize the function of the ecosystem, based on the stability of the ecosystem structure.
Runoff and meteorological data used in this
research paper are available via the Yellow River Conservancy Commission
(
The authors declare that they have no conflict of interest.
We acknowledge the use of runoff data from the Yellow River Conservancy
Commission (