Characterizing drought in terms of changes in the precipitation–runoff relationship: a case study of the Loess Plateau, China

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 ≥ 5 years and mean annual precipitation anomaly (PA) values ≤−5 % during drought periods. The results show that continuous precipitation shifts are able to change the water balance of watersheds in water-limited areas, and multiyear drought events cause the PRR to change with a significantly decreasing trend (p < 0.05) compared to other historical records. For the Loess Plateau as a whole, the average runoff ratio decreased from 10 to 6.8 % during 1991– 1999. The joint probability and return period gradually increase with increasing of drought duration and severity. The ecosystem service of water yield is easily affected by drought events with durations equal to or greater than 6 years and drought severity values equal to or greater than 0.55 (precipitation ≤ 212 mm). At the same time, multi-year drought events also lead to significant changes in the leaf area index (LAI). Such studies are essential for ecosystem management in water-limited areas.


Introduction
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 20 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). At regional scales, during 1998-2001, one-third of the territory of Iran was affected by the most severe drought in the history of the 2 country. More than half of the country's population faced food shortages and a lake 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 US$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). 5 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 10 Bhalme-Mooley index (Bhalme and Mooley, 2009). Others 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 15 periods that feature dry soils 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 where water demand exceeds supply because of long-term precipitation shortages. In addition, socioeconomic drought reflects 20 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 25 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 30 resource planning and management, and watershed system evolution Nourani et al., 3 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 PPR. Sun et al. (2016) assessed hydrologic trends in urban catchments using a conceptual urban precipitation-runoff model, which is beneficial for the stormwater management and planning. There is no doubt that variations in the PPR reflect the actual integrated volume of precipitation and runoff. The ratio between 5 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 PPR (Petrone et al., 2010;Saft et al., 2015). Therefore, the shift in the PPR caused by an extended 10 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. 15 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 20 system on the Loess Plateau of China, the ecosystem service of water yield is of concern for both the sustainability of revegetation programme and the objective of raising the level of economic development Feng et al., 2016a). In recent years, the climate of the Loess Plateau climate 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, analysing 25 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 30 conditions and the PRR to define drought. We then simulate the drought frequency and return period of 4 drought events of different magnitudes. Finally, we discuss the policy implications of our study for ecosystem management.

Study area
The Loess Plateau (Fig. 1) is located between 100°54É-114°33É and 33°43Ń-41°16Ń and covers a 5 total area of 624,000 km 2 . It has a semi-arid and semi-humid climate. From 1960 to 2000, the annual precipitation ranged from 200 mm in the northwest to 700 mm in the southeast. It is a key area in current and past efforts to conserve soil and water in the Yellow River basin. The ecology of the Loess Plateau is sensitive to climate change. 10 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 thirteen 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 15 (http://www.yellowriver.gov.cn/) and the National Meteorological Information Centre (NMIC; http://cdc.nmic.cn/home.do) of the Chinese Meteorological Administration (CMA), respectively. The AVHRR GIMMS LAI3g datasets covering the period from 1982 to 1999 were also used in this study.
These datasets were generated from AVHRR GIMMS NDVI3g data using an Artificial Neural Network (ANN) based model (http://sites.bu.edu/cliveg/). Some studies have shown that the Grain for Green 20 Program(GGP), which began in 1999, resulted in a reduction of runoff from the Loess Plateau (Zhang et al., 2008;Feng et al., 2016b), so the study period is selected from 1961 to 1999.

Drought identification
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,25 2010). Moreover, because we are interested in determining whether the runoff response differs for multiyear 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 30 5 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 5 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 ≤ -10% or (2) mean PA values of less than 0 for more than three consecutive years. Note that the PA value of the starting year of each drought period is negative. 10 In this study, the cumulative PA values during each drought period are used to measure drought severity (for convenience, drought severity is multiplied by -1 to obtain a positive value). Based on the rules mentioned above, we identified all of the drought events that occurred in each watershed. To reflect the response of the PRR to drought events over the years, we must ensure that the dry periods are sufficiently long and severe. In the subsequent analysis, we consider only drought events with durations ≥ 5years and 15 mean annual PA values ≤ -5% during the drought period. Finally, the dry events are classified into major dry period and single-year dry period.
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 20 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 t-test (p ≤ 0.05). The historical records refer to scatterplots of annual precipitation-runoff during the period of 19611999, except for particular major drought periods. For example, when the drought that occurred in 19701974 is considered, the corresponding historical record includes a 25 precipitation-runoff scatter plot that includes data from 19611969 and 19751999. For the drought that occurred in 19911999, the corresponding historical record refers to a precipitation-runoff scatterplots containing data from 19611990.

Drought frequency analysis
Drought can be characterized by multiple variables, such as duration, severity, and spatial extent 30 (Steinemann and Cavalcanti, 2006;Hayes et al., 2012), but how to determine the joint distribution 6 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(d) and drought severity(s) are F D (d) and F S (s) respectively, a copula C exits that combines these two marginal distributions to give the joint distribution function, F D,S (d , s): If the marginal distributions F D (d) and F S (s) are continuous, f D (d) and f S (s) are the density functions corresponding to F D (d) and F S (s),respectively, and the joint probability density function becomes: where c is the density function of C, which is defined as:

Parameter estimation
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., 15 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 20 optimal marginal distribution function. The ML method is used to estimate the distribution parameters of drought duration ( 1 , 1 ) and severity ( 2 , 2 ), and the parameters of copula function are then calculated using the following formula: where lnL C is the log-likelihood function of the copulas.

The determination of the optimal joint distribution function
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 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. 10 N is the drought series length and n is the number of drought events. The return period of a single variable can be obtained from the definition of the copula function as:

Drought return period
The joint distribution function of drought duration and severity is: The joint return period of the two characteristic variables is calculated as:

Drought events on the Loess Plateau
Of the 39 years that make up the study period, the precipitation anomaly is negative in 17 years. Based on 20 the drought identification method developed in this study, 7 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 five years. The corresponding drought severities are 0.51, 0.32 and 0.61, and the average precipitation values of these events are 231 mm, 320 mm and 183 mm, respectively. 25 8 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 5 (p=0.692) and 1979-1983 (p=0.514), although the regression lines for these dry periods deviate from the overall regression line. However, a significant decrease in the PRR can be identified for 1991-1999 (p=0.000). In this period, the dry period regression lines are lower than nearly all of the other points, indicating unprecedentedly low runoff generation rates for a given amount of precipitation. The runoff ratio is approximately 10 percent from 1961 to 1990. In the dry period of 1991-1999, the average runoff 10 ratio decreased to 6.8 percent. Table 3 lists the estimated parameters and the results of goodness-of-fit tests for the marginal functions.

Drought frequency on the Loess Plateau
We find that not all of the distributions pass the K-S test at the 95% (=0.05) significance level. Further considering the RMSE, the marginal distributions that provide the best fits to the drought duration and 15 severity values are Weibull and gamma distributions, respectively. The results for these distributions are shown in bold font and underlined in Table 3. The goodness-of-fit for the joint function was evaluated by calculating the SED, RMSE and AIC. These metrics indicate that the Frank copula function is best suited to fitting the duration and severity of the drought events in the study area, except for the Jialu, Dali and Beiluo watersheds. Note that, for the three watersheds mentioned above, the values of the SED, 20 RMSE and AIC are also relatively small. Finally, we choose the Frank copula function to represent the joint distribution function of drought duration and severity. Fig. 4a and 4b 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 25 drought severity ranges from 0.1 to 0.35 (which corresponds to precipitation values that range from 306 mm 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 five 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 30 severity is not more than 0.4 (precipitation ≥282 mm), the joint probability of drought increases rapidly 9 with increasing drought duration, and the rate of increase of the joint probability decreases as drought duration increases when drought severity is greater than 0.4 (precipitation≤282 mm). The density of the contours shows that, on the Loess Plateau from 1961 to 1999, drought events with drought duration ranging from 2 to 5 years and drought severity ranging from 0.2 to 0.5 (for which the corresponding precipitation values are 235 mm -470 mm) occurred easily. The joint probabilities of drought 5 characteristics are key to drought management. The probability that both drought duration and severity simultaneously exceed certain thresholds is useful information for environmental and government agencies responsible for water system management under drought conditions. Fig. 4c and 4d show that, as drought duration and drought severity increase, the joint return period also 10 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)(1971)(1972)(1973)(1974).
The return period reaches 4.32 years for the drought that occurred in 1979-1983, which displayed a 15 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. In 1979-1983, the drought duration reached 5 years, which is close to the 20 estimated return period.

Variability of the PRR during the dry periods in each watershed
Prolonged multi-year drought causes significant damages in natural environments. Fig. 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 25 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 (p < 0.05) occurred in 9 cases, accounting for approximately 60% of the total cases. In these cases, the dry period regression line lies lower than nearly all the other parts of the historical record, indicating 30 unprecedentedly low runoff generation rates for the given precipitation values. Thus in a sequence of 10 years with decreased precipitation, we can conclude that lower runoff not only relate to the lower precipitation, but also less runoff than expected caused by the multi-year drought period.
There were no significant changes in the PRR in 5 of the 13 watersheds during the study period. Fig. 6 demonstrates that there is no geographical pattern in the spatial distribution of watersheds with and 5 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.

Spatial variability of drought frequency on the Loess Plateau
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 study of the characteristics of the spatial distribution of the joint return period. Figs. 7a,7b and 7c show that, although the duration and severity of drought differ among these three major drought events, the spatial distribution of the drought 15 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, 20 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 analyse 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 25 in these watersheds is particularly important in ecosystem management. 11 4. Discussion

Reliability of the identification of drought events
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 5 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 ≥7 years and severity ＜-5%. In this study, the drought duration was limited to not less than 5 years. Linking climate and the change in the PRR to identify drought events on the Loess Plateau from 1961 to 1999 shows that 1962, 1965, 1986-1987 and 1989 are single-year dry periods, 10 whereas1970 -1974, 1979-1983 and 1991-1999 are major dry periods. The years when major dry events occurred in the Yellow River basin from 1961from to 1999from were 1965from ,1972from ,1980from ,1995from and 1997from (Fu et al., 2008; http://www.mwr.gov.cn). For example, a severe drought in northern China in 1997 damaged 1.94 million hectares of crops in the Yellow River basin and resulted in 226 days of zero flow from Henan province to Shandong province. The total length of the river with zero flow was approximately 687 km. 15 Thus, the dry events identified by the method described in this study are consistent with historical data, further illustrating the reasonableness of the method presented here.

The influence of multi-year drought events on the ecosystem service of water yield
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 20 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, 25 irrigation schemes can be neglected in the study area (Feng et al., 2016). 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 12 water consumption for irrigation and the degree of runoff variability caused hydropower is three 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 5 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, 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 10 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 , 15 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 (p=0.048, Fig. 9). The effects of other anthropogenic activities, including dam construction, tree 20 plantations and pasture, did not cause the observed change in the runoff ratio in this period. Terrace construction contributed 25% of the reduction in the runoff ratio in the 1990s. Thus, drought events are the major factor driving the reduction in runoff in the study area.
The PRR responds to multi-year droughts in two different ways. The PRR either displays no significant 25 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. Analysing the changes of 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 30 occurrence of dry events will aggravate the reduction of runoff, which could lead to significant changes in 13 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 5 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 ≥ 6 years and the drought severity is ≤ 0.55 (i.e., the precipitation ≤ 212 mm). The occurrence of less runoff than expected 10 causes problems for the ecosystem service of water yield. The shift in the PRR induces a contradiction between the expected runoff and the actual amount of water in a watershed if the runoff is not predicted accurately when the PRR changes. This result suggests that we must consider the impacts of prolonged drought events on ecosystems and optimize the modelling techniques used to assess the PRR to cope with the effects of long-term droughts that occur in response to changed climatic conditions. 15

Policy implications for ecosystem management
In late 1999, China implemented a massive revegetation programme in the name of ecological restoration. This restoration involved returning croplands on steep slopes ( > 25°) to woodlands, shrublands or grasslands. Based on the spatial distribution of drought events on the Loess plateau, the 13 watersheds are further divided into four regions. The drought return periods in the northern and western 20 regions are lower than those in the central and eastern regions, and we can see the distribution of the effects of the project in the four different regions (Fig. 10b). The different spatial distributions of the return period further reflect regional differences.
Drought events affect ecosystem productivity and reduce the carbon sink capacity of terrestrial 25 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 CO2) and momentum are exchanged between the vegetated land surface and the planetary boundary layer (Myneni., 2002). During 19821999, the LAI displays an obvious increasing trend in the 1980s over the Loess Plateau (Fig. 11). To illustrate the impacts of drought on the 30 14 ecosystems, we choose the major drought period of 19911999 and the normal period of 19841990.
We find that the LAI decreases significantly (p=0.032, Student's t-test) in 19911999 compared to that of 19841990. Since different vegetation types respond differently to drought events over different time scales, the regional drought return period and vegetation types should be taken into account in future policies, and appropriate policies should be formulated according to the specific regional conditions to 5 control the adverse effects of drought on ecosystem productivity. We must pay more attention to these drought characteristics, which can induce significant changes in the PRR, and the return period (Fig. 7d).
This practice will help avoid the waste of capital investment and will effectively improve the implementation of the GGP. 10 The results of this study provides 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 maizes, 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 15 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 20 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.

Conclusions
This article characterizes drought by linking climate anomalies with changes in the PRR. We found that 25 multi-year drought caused the PRR to display a significant decreasing trend (p < 0.05) compared to the historical records. For the Loess Plateau as a whole, the average runoff ratio decreased to 6.8 percent of the average annual precipitation in 19911999 compared to 10 percent of the annual average precipitation in . In 8 of the 13 studied watersheds, significant changes in the PRR occurred during 1961-1999. When we compared the annual average precipitation in separate watersheds and 30 analysed their drought characteristics, we concluded that this situation is likely to occur when the drought duration is not less than 6 years and the drought severity is equal to or greater than 0.55 (i.e., the annual precipitation is ≤212 mm).
Our analysis revealed great spatial variability in drought across the Loess Plateau. We chose the Frank 5 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 analysing the annual LAI over the Loess Plateau during 1982-1999, we found that the LAI significantly 10 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 15 maximize the function of the ecosystem, based on the stability of the ecosystem structure.