Drought Propagation Patterns under Naturalized Condition Using Daily Hydrometeorological Data

Drought propagation pattern forms a basis for establishing drought monitoring and early warning. Due to its regional disparity, it is necessary and significant to investigate the pattern of drought propagation in a specific region. With the objective of improving understanding of drought propagation pattern in the Luanhe River basin, we first simulated soil moisture and streamflow in naturalized situation on daily time scale by using the Soil and Water Assessment Tool (SWAT) model. ,e threshold level method was utilized in identifying drought events and drought characteristics. Compared with meteorological drought, the number of drought events was less and duration was longer for agricultural and hydrological droughts. ,e results showed that there were 3 types of drought propagation pattern: from meteorological drought to agricultural/hydrological drought (M-A/H), agricultural/hydrological drought without meteorological drought (NM-A/H), and meteorological drought only (M). To explain the drought propagation pattern, possible driven factors were determined, and the relations between agricultural/hydrological drought and the driven factors were built using multiple regression models with the coefficients of determination of 0.4 and 0.656, respectively. ,ese results could provide valuable information for drought early warning and forecast.


Introduction
Drought is a complicated, recurrent natural hazard with characteristics of covering extensive areas and lasting for months to years [1,2], which have significant negative impacts on economy and ecology [3][4][5][6].Drought is commonly defined as a persistent occurrence of below-normal natural water availability [7], which is applicable to each part of the terrestrial hydrological cycle.In general, drought is distinguished into different types related to components of hydrological cycle: meteorological drought (a precipitation deficit), agricultural or soil moisture drought (a belownormal storage in the unsaturated zone), hydrological drought (water availability lower than the normal in aquifers and/or streams), and socioeconomic drought (occurring when the demand of various commodities for water exceeds the supply) [8].
Drought starts with periods of precipitation deficit and propagates through the terrestrial hydrological system, which in turn results in drought in soil moisture and groundwater or streamflow [7,9].e developing process is called drought propagation.e term "drought propagation" was introduced as a theoretical framework by Changnon [10], nevertheless, created as a term by Eltahir and Yeh [11].Precipitation deficit with a prolonged period decreases streamflow, subsurface water and groundwater storage, and generated hydrological drought [12,13].Based on Van Loon [14], the features of drought propagation include lag, lengthening, attenuation, and pooling.Hisdal and Tallaksen [3] showed streamflow droughts are less frequent and more persistent than precipitation droughts.Van Lanen [15] demonstrated propagation of recharge drought led to slighter hydrological drought in different climate regions.
In previous studies, it was demonstrated that drought propagation was affected by climate seasonality [13,[16][17][18].Sung and Chung [13] pointed that streams with lower than the average discharge during high-stream seasons might affect drought development.Huang et al. [16] examined the propagation time from meteorological to hydrological drought using the cross wavelet analysis, and found the propagation time has seasonality with short time in spring and summer and long duration in autumn and winter.Based on drought propagation processes in different seasons, Van Loon and Van Lanen [18] divided hydrological drought into six types, which resulted from interaction of precipitation and temperature in various seasons.Van Loon et al. [17] revealed that how seasonality of climate influenced drought propagation.
e key steps of drought propagation studies are to identify drought event and define the characteristics of drought events.Drought frequency (Number of drought events) and severity (drought duration and deficit volume) are vital features to define characteristics of drought events, which can reflect translation of drought signals in a hydrological system [7,19,20].Tallaksen et al. [20] showed deficit volume was a robust feature to weigh on the severity of drought event over the catchment area.Many studies have concentrated on defining drought and drought characteristics using various indices, which include Standardized Precipitation Index (SPI), Standardized Streamflow Index (SSI), Standardized Precipitation Evapotranspiration Index (SPEI), and Standardized Runoff Index (SRI) [16,21,22].
Apart from indices, the above-mentioned drought features also can be obtained by the threshold level method.e threshold level method can truncate continuous time series with a defined threshold value, and it includes fixed and variable threshold.Variable threshold, changing over the year, reflects the seasonality, which has been widely applied to the study of drought propagation [7,[23][24][25][26][27].Compared with drought indices, threshold level method can calculate deficit volume, which is an important feature in water management.It does not require a preferred knowledge of probability distributions [21].On the other hand, all drought categories of hydrometeorological variables can make comparison by threshold level method, which is necessary when studying drought propagation.It has been testified that this method has a potential capacity to analyze dailybased drought [19,28,29].
In order to analyze drought propagation including agricultural drought and groundwater drought, hydrological models are usually applied, such as the conceptual, semidistributed rainfall-runoff model-HBV [14,19], combined SWAP and MODFLOW models [30], and spatially distributed physically based model SIMGRO [15].ese models were employed to simulate the soil moisture and groundwater time series, which cannot be observed directly.In regions where human activities affect runoff processes, hydrological models are essential to distinguish the naturalized and human-induced drought.
In this study, the primary objectives are (1) to investigate drought propagation patterns on daily scale with naturalized situations (without human effect) and ( 2) to explain the driven factors that affect drought propagation in the Luanhe River basin.

Study Area
e Luanhe River basin is located in the northeast part of China, which is between 115 °30′E∼119 °15′E and 39 °10′N∼42 °30′N (Figure 1).e Luanhe River originates from Mongolia Plateau, passes through the Yanshan Mountains, and eventually flows into the Bohai Bay with a total drainage area of 44600 km 2 .
ere are three geographic types in the basin, which are plain, mountains, and plateaus.
e altitude ranges from 2 to 2205 m with an average of 766 m.
e basin belongs to temperate semiarid continental monsoon climate with cold winter and hot summer.Annual average temperature is −0.3∼11 °C.Annual potential evaporation is 950∼1150 mm.Average annual precipitation is 535 mm, with 70%∼80% of annual rainfall in the flood season (June-September). is attributes to unstable characteristics of the duration, intensity, and location of the subtropical high over the northern Pacific in summer [31].
e Luanhe River basin has a strong seasonality in meteorological forcing, with relatively low precipitation and low potential evapotranspiration in winter and spring and relatively high precipitation and high potential evapotranspiration in summer.
is results in a strong seasonality in discharge.Highest river flows occur in summer and the lowflow season is winter.
Luanhe River basin suffered from severe droughts in 1972, 1980-1984, and 1997-2005.e annual runoff in this basin showed significant decreasing trend and resulted in a long hydrological drought during last decade, which affected the water supply to Tianjin and Tangshan cities. erefore, it is of great significance to select Luanhe River basin as the study area to research drought propagation.

Data
For the study area, daily rainfall data were available from 13 rain gauges and  [34,35] is one of the mostly widely applied nonparametric tests for trend detection in hydrologic time series.Before employing the MK trend test, the trendfree prewhitening (TFPW) procedure proposed by Yue et al. [36] is usually used to efficiently eliminate the influence of serial correlation on the MK test.
For identifying the true slope of MK trend analysis, the Sen's slope estimator developed by Sen [37] is commonly adopted [38,39].e Sen's produce applied following the MK test measures the magnitude of any significant trend found in the MK test.
In this study, the MK trend test with TFPW procedure and Sen's method were used to determine the temporal trends in precipitation.Details of the three methods can be found in Kisi and Ay [40]; Yue et al. [36], and Da Silva et al. [41] respectively.

Identification of Change Points.
ere are many methods to determine the change point of hydrological series, such as moving F test method [42], R/S analysis method [43], and Bayesian model [44].
e nonparametric Mann-Kendall-Sneyers test [34,35,45] was applied in this study to determine the occurrence of a change point.
is test is widely used to detect abrupt change in hydrological data, because it has the advantage of not assuming any distribution form for the data and has similar power to its parametric counterparts.e normalized variable statistic UF k of the forward sequence and UB k of the backward sequence are calculated and used to plot the forward and backward curves.If the intersection of the two curves occurs within the given confidence interval, then it indicates a change point.
In addition, the nonparametric Pettitt test developed by Pettitt [46] was further employed in this study to examine the change points obtained by the Mann-Kendall-Sneyers test.e Pettitt test detects a significant change in the mean value of observed series when the exact time of the change is unknown [46].A version of the Mann-Whitney statistic U t,N is used to test whether two divided segments of a sequence of random variables are from the same population.
Readers may refer to Chen et al. [47] and Xie et al. [48] for a detailed discussion on the Mann-Kendall-Sneyers test and the Pettitt test, respectively.

Description of SWAT Model.
In this paper, the Soil and Water Assessment Tool (SWAT) model was applied to simulate daily streamflow and soil moisture time series in Luanhe River basin.It is a semidistributed, process-based hydrological model.In SWAT model, the watershed is subdivided into Hydrologic Response Units (HRUs), in which runoff generation is simulated by CN method or the Greeen-Ampt infiltration method.
e SWAT model has been widely applied to runoff simulation on different time scales.
SWAT-CUP, an automatic parameter estimation method, is used to estimate the parameters in the SWAT model.In SWAT-CUP, algorithm of SUFI-2 is selected to perform calibration and sensitivity analysis.
e model performance is assessed by Nash-Sutcliffe model efficiency coefficient (Ens) [49] and the coefficient of determination (R 2 ), calculated as where S i and O i are simulated and observed discharge, respectively.O is the average of the observed discharge.In general, monthly Ens is deemed satisfactory at >0.5 and daily and monthly R 2 are satisfactory at >0.5 [50].
Firstly, a sensitivity analysis of threshold value to drought event identification was carried out to select a proper threshold for the considered study area.e next step is to calculate threshold value according to the selected percentile of the monthly duration curves.Because monthly threshold series is discrete (a "staircase" pattern), a centered moving average of 30 days is applied to smoothen the threshold values.Subsequently, drought event and its drought characteristics (drought duration; onset and end date of drought; drought deficit volume) can be calculated based on the selected threshold value.In order to assure comparability of drought characteristics, it is crucial to apply the same threshold level (the percentile) in drought propagation studies [14].
A drought event is defined to begin when a variable declines below the predefined threshold and continue until the variable exceeds the threshold [58].
To identify whether the variable x lies in a drought situation on day t, a binary variable δ is represented as the following equation: where δ refers to the binary variable and τ is the predefined threshold value.e duration (Δ i ) of a drought event i can be calculated as 4 Advances in Meteorology where T is the total length of the variable x and Δt is the time interval.
where d(t) is the deviation from threshold value τ on the day t. Di where D i refers to the deficit volume of drought event i.
where d i, max is the maximum deviation of drought event i. e deficit volume is viewed as the most appropriate characteristic to measure severity of drought.
Note that deficit volume (D i ) can only be calculated for flow variables, for example rainfall and streamflow with dimension (mmT −1 ).Soil moisture is a state variable, similar to groundwater head, thus a measure of storage.e storage coefficient must be known to convert soil moisture into a flow variable.It requires detailed information about the properties of soil layers, which was not available.erefore, deficit volume cannot be calculated for soil moisture and instead maximum deviation (d i, max ) is commonly used to measure severity of drought.
Let DS i represents the average severity of drought event i, calculated as where D i and Δ i are the deficit volume and duration of drought event i, respectively.Two restrictions are used for identifying minor droughts according to Hisdal et al. [59]: (1) drought events with duration shorter than 5 days or (2) drought events with a deficit volume less than 0.5% of the maximum deficit volume in all the drought events.To reduce the effect of too much minor droughts, the minor droughts which satisfy the two restrictions were not considered in this study.

Drought Propagation.
Drought propagation is a dynamic process which focuses on the propagation from meteorological droughts to agricultural and/or hydrological droughts.Based on the identified drought events for all drought types in the above section, drought propagation patterns will be determined through the analysis of drought characteristics to answer 2 questions: (1) Does meteorological drought definitely result in agricultural and/or hydrological droughts?(2) Are agricultural and/or hydrological droughts necessarily caused by meteorological drought?en, the driven factors for drought propagation patterns will be found, and drought propagation equation will be established to explain the drought propagation patterns in the study area.
In this study, 3 different drought propagation patterns were identified.
e first one is that only meteorological drought occurs without other drought types following (M).
e second one is the case in which a meteorological drought occurs followed by an agricultural or hydrological drought.In its definition, the maximum lag time of onset from meteorological to agricultural/hydrological drought is 45 days.
is drought propagation pattern is termed as agricultural or hydrological drought originated from meteorological drought (M-A/M-H).
If the lag time between agricultural/hydrological drought and its previous meteorological drought is longer than 45 days, it is considered that the agricultural/hydrological drought is not caused by meteorological drought, namely, the third propagation pattern (NM-A/NM-H).

Temporal Changes in Precipitation.
In this work, a proper modelling of the hydrological response of the Luanhe River basin is a critical issue which was carried out by SWAT model.e model was calibrated in naturalized conditions and then applied to simulate discharge and soil moisture data for the period 1980∼2012, with the purpose of removing human influences.Indeed, this approach relies on the assumption that precipitation pattern remains unchanged over the period 1963∼2012.Hence, several statistical methods were employed to verify the stationary precipitation pattern.
Time series of annual precipitation during 1963∼2012 over the Luanhe River basin was obtained by using iessen polygon technique.
e statistic Z values of MK test with TFPW procedure and Sen's slope estimator (Q med ) were calculated for the annual precipitation series, with the values of −0.842 and −0.906, respectively.e results indicated a nonsignificant decreasing trend in the annual precipitation at the significance level α � 5% (|Z 1−α/2 | �1.96).
e Mann-Kendall-Sneyers test and the Pettitt test were applied at the confidence level of 95% to detect the change point of the annual precipitation series.From the result of Mann-Kendall-Sneyers test (shown in Figure 2 Based on the above results, the precipitation over the Luanhe River basin is free of significant temporal change, implying stationary precipitation pattern during the period 1963∼2012.[50], if the value of Ens is greater than 0.50, the SWAT model performs well.

Calibration and
Figure 3 depicts the simulated and observed runoff values on daily scale at Luanxian hydrological station, and it can be seen that the peak of simulated runoff is lower than the measured values, especially in 1964.
For the underestimation of peak flow in model calibration and validation, there might be two major reasons as follows.Firstly, when calibrating the model parameters, more attention was paid on the performance of simulated streamflow processes during dry periods.Secondly, due to the little accumulation of snow during winter, snowmelt made no contribution to the streamflow during spring in the Luanhe River basin.ere is no significant snow influence in the basin.us, rainfall and its characteristics are the dominant factor driving peak stream flows.However, on account of the limited observed materials available, precipitation data from only 13 rainfall gauge stations were used as the precipitation inputs for SWAT modelling.Accordingly, these limited data could not completely represent the precipitation field of the entire basin.Moreover, with the unevenly distributed precipitation both in space and time, precipitation fields in the Luanhe River basin are spatially variable.As for a rainstorm-generated stormflow event, the model is most unlikely to accurately capture the real center of rainfall, whereupon leading to some deviations existing between the actual and modeled precipitation processes.Consequently, the simulated peak flows showed poor performance.
Figure 4 shows the scatter plot of observed and simulated daily discharge.e low flow is very close to 1 : 1-line, which agrees well with the observed in undisturbed period.Compared with high flow, the model has a good performance to simulate low flow.Hence, this model can be used to simulate discharge without human influences during 1980-2012.

Sensitivity Analysis of reshold Value.
To evaluate the impact on drought characteristics and select a proper threshold for the study area, the 60th, 70th, 80th, 90th, and 95th percentile of the monthly duration curves were considered to calculate the threshold values, using the observed precipitation data in Luanxian station.According to the different threshold values, 5 sets of meteorological drought events and related characteristics (duration (Δ i ) and deficit volume (D)) were identified.We analyzed and compared the five groups of results for both drought duration and deficit volume.Table 2 summarizes the total number, as well as the average value and standard deviation of Δ i and D, for the 5 sequences of drought events.
In general, as the percentile decreased, the threshold value increased, leading to the number of drought events rising, as well as the average values of both Δ i and D. For the too-low threshold value (e.g., the one corresponding to the 95th percentile), most of the defined drought events were essentially extreme events, unexpectedly having too short drought duration and undersize deficit volume.However, for the too-high threshold value (e.g., the one related to the 6 Advances in Meteorology 60th percentile), extreme drought events were overestimated, and their duration and severity both oversize.In addition, there seemed to be too many small droughts, which is inconsistent with the fact.As for the 70th, 80th, and 90th percentile, the related 3 sets of drought events had similar average values of drought duration and deficit volume, while the standard deviation descended with the percentile increasing.Several previous studies investigated the features (including drought duration, drought severity, and drought peak) of meteorological droughts in the Luanhe River basin (e.g., Ren et al., Zhou et al., Wang et al., and Wang et al.).By comparing the previous researchers' study results with ours, the drought features defined by the threshold value based on 80th percentile show performance more superior.Consequently, the 80th percentile seems to be applicable to define a proper threshold for the study area.e 80th percentile, locating in the range of 70th to 95th percentile, is often used for perennial rivers [60,61].
According to Van Loon [14], the same type of threshold is applied to precipitation, runoff, and soil moisture data for the sake of comparison.erefore, the threshold value was defined as the 80th percentile of the monthly duration curves.

Characteristics of Different Drought Types.
According to the simulation results from SWAT model, several characteristics of different drought types were determined.Table 3 shows the characteristics of different drought types.e number of meteorological drought events is obviously the largest, which is 1.4 and 1.5 times larger than hydrological drought and agricultural drought, respectively.However, meteorological droughts have considerably shorter duration, which is about half of agricultural drought on average.From meteorological drought to hydrological drought, the mean deficit volume attenuates considerably.
e longest duration of hydrological drought is more than 200 days, which is comparable to agricultural drought.80% of drought events are less than 70 days for all drought types.Drought events with less than 30 days account for 74%, 67%, and 66% for meteorological, agricultural, and hydrological drought, respectively.e largest deficit volume of meteorological and hydrological drought is 44 mm and 10 mm, respectively.80%  Advances in Meteorology of meteorological drought events are less than 8.2 mm in deficit volume, and 80% of hydrological drought events are less than 1.67 mm.Deficit volume of drought events is markedly larger for precipitation, and the largest deficit volume is 4 times higher than hydrological drought.e reason is that rainfall is larger and more variable, leading to higher deviation from the threshold [18].e upper Luanhe River basin is located on the Bashang Plateau, the middle reaches flow through Yanshan Mountains, while the downstream area has relatively flat topography.Due to the complex terrain and uneven distributions of precipitation, the streamflow regimen of the basin is dominated by precipitation.In addition, temperature is playing an increasingly important role under the impact of global warming.Accordingly, not only the agricultural drought caused by soil water deficiency would develop into the hydrological drought caused by lack of streamflow, but also the meteorological drought, which results from precipitation shortage and/or increased evaporation, would directly lead to the hydrological drought. is is the main reason behind our results that the identified agricultural drought events are less than hydrological drought events following Table 3.

Drought Propagation Types.
In general, meteorological drought develops into agricultural drought and then hydrological drought.Based on the previous drought identification, each drought event was determined, and there were 3 different drought propagation patterns (Table 4).83 out of 131 meteorological droughts developed into agricultural or hydrological drought, which accounts for 63.4% of the total meteorological drought events (M-A/M-H).But not all meteorological droughts develop into agricultural and hydrological drought, which means only meteorological drought occurs, and no other drought types follow (M).48 out of 131 meteorological droughts belonged to this propagation pattern.
e third propagation type is that agricultural or hydrological drought is not caused by meteorological drought (NM-A/NM-H).
Table 5 shows the characteristics of each drought type for different drought propagation patterns.For pattern M, the average duration of meteorological droughts is 25.0 days, and the average deficit of meteorological droughts is 1.9 mm.For those meteorological droughts developing into agricultural or hydrological droughts (M-A/M-H), the average meteorological drought duration are 23.6 days and 22.2 days, and the average meteorological drought deficit are 7.9 mm and 7.6 mm, respectively, which are 4.16 times and 4.00 times larger than those for pattern M. erefore, severe meteorological droughts are prone to develop into agricultural and hydrological droughts.
As to agricultural drought, M-A propagation pattern has a larger average duration value of 46.7 days than NM-A pattern, which has the similar results with maximum deviation.For hydrological drought, the average duration and deficit are 42.9 days and 1.4 mm, respectively, for M-H pattern, which are larger than NM-H pattern. is demonstrates that meteorological droughts can lead to relatively

Drought propagation types
Number of events M-A/M-H 83 M 48 NM-A/NM-H Advances in Meteorology severe hydrological droughts.Although there are no meteorological droughts for NM-A and NM-H patterns, which means the daily precipitation is higher than the threshold value, the precipitation is below the multiyear average before the occurrence of agricultural and hydrological droughts for several days.

Lag Time from Meteorological to Agricultural/ Hydrological Drought.
e onset and end time are two important characteristics for drought events, and the lag time from meteorological to agricultural/hydrological drought can express drought development.erefore, we selected drought events from M-A/M-H drought propagation type which include meteorological, agricultural, and hydrological droughts occurring concurrently, to calculate the average lag time of onset and end time from meteorological to agricultural/hydrological drought.e results are listed in Table 6.
54% of meteorological drought events lead to agricultural droughts, and then hydrological droughts, and there are 39% drought events for which hydrological droughts started before agricultural droughts.e average lag of onset from meteorological to agricultural and hydrological drought is 21.3 days and 26.8 days, respectively.e average lag of end time from meteorological to agricultural and hydrological drought is 36.6 days and 59.2 days, respectively, in which 48% agricultural droughts ended after hydrological droughts.

From Meteorological Drought to Hydrological Drought.
According to the NM-A/NM-H drought propagation patterns, meteorological drought is not the only driven factors for agricultural and hydrological droughts.ere must be some other factors affecting the occurrence of agricultural and hydrological drought, i.e., precipitation before the onset of agricultural/hydrological drought, evapotranspiration, and so on [62].
Based on monthly precipitation and runoff data, Wu et al. [63] built a nonlinear relationship between hydrological drought duration and meteorological drought duration and hydrological drought severity and meteorological drought severity, respectively.But they did not consider other factors.In this paper, the following factors are chosen as the driven factors for hydrological drought severity (DS H ): (1) simulated daily evapotranspiration during hydrological drought (E), (2) meteorological drought severity (DS M ), (3) simulated mean soil moisture deficit (SD), and (4) the summation of the difference between the daily precipitation and the long-term average daily precipitation 10 (PD 10 ), 20 (PD 20 ), 30 (PD 30 ), and 40 (PD 40 ) days before hydrological drought.en the relations between hydrological drought severity and each driven factor was analyzed using Pearson correlation analysis, and the results are shown in Table 7.
It can be seen that SD, DS M and E are positively related to DS H significantly, and PD 10 is negatively correlated with DS H significantly, but PD 20 , PD 30 and PD 40 are not significant with the significance level of 0.01.With the increase of SD, DS M and E, together with the decrease of PD 10 , hydrological drought may occur.
en the relationship between DS H and each significant driven factor is described by an exponential function, shown in Figure 5.
To further illustrate the drought propagation pattern of M-H and NM-H, the relationship between drought severity (DS H ) and all the driven factors is built using multiple regression analysis, which can be expressed as DS H � 0.011 * EXP(1.130DSM + 0.540E + 0.012SD − 0.035PD 10 ) − 0.008.
e coefficient of determination (R 2 ) is 0.656.e calculated DS H and observed DS H are plotted in Figure 6. e plots are scattered systematically along the 1 : 1 line.erefore, the equation can reflect the comprehensive impacts of the selected driven factors on hydrological drought.e conditions of these 3 propagation types are listed in Table 8.

From Meteorological Drought to Agricultural Drought.
Similar to the propagation from meteorological to hydrological drought, some contributing factors affecting agricultural drought severity (SD) are chosen: (1) DS M , (2) evapotranspiration during agricultural drought (E A ), (3) PD 10 , PD 20 , PD 30 , and PD 40 .e Pearson correlation coefficients between SD and each contributing factor are listed in Table 9.It can be seen that DS M , EA, and PD 10 passed the significance test with the significance level of 0.01.DS M and EA have a positive impact on agricultural drought, e coefficient of determination of Equation ( 9) is 0.4.e calculated SD and the observed values are plotted in Figure 7, which scattered near the 1:1 line. is model could explain the drought propagation type of M, M-A, and NM-A.e conditions of the 3 drought propagation types from meteorological to agricultural drought are listed in Table 10.

Discussion
is study uses daily rainfall, soil moisture, and runoff data to identify drought propagation patterns.Compared with monthly data, daily data perform better in identifying the drought characteristics (the precise drought onset, duration, end time, and number of drought events).In some previous studies in which monthly data were used, the relations between meteorological drought index and hydrological drought index with l-month (l �1, 2, 3. ..) lag time were built, and l month with the highest correlation coefficient was considered as the lag time, which was invariant for a given watershed [16,64,65].However, because of different drought characteristics for each drought event, the lag time must be changed from event to event.Drought lag time on the basis of daily data could avoid such a problem.To remove the impacts of human activities, we applied the SWAT model to simulate streamflow from 1963 to 2012 under the 1970 land use condition.is concept stemmed from the observation-modelling framework presented by Van Loon and Van Lanen [18]; who applied a conceptual hydrological model of HBV in a case study area in Spain.Because there are no observed soil moisture values in the entire watershed, another advantage by using SWAT is to provide simulated daily soil moisture to identify agricultural droughts.Groundwater system was not considered for hydrological drought in this study, since the data are not available.However, streamflow drought could represent the hydrological drought in many previous studies [66][67][68].us, rainfall, soil moisture, and streamflow data could be used to analyze drought propagation patterns.e deficit volume, which was viewed as the best drought characteristic to study drought propagation [58], was not used to calculate for soil moisture storage [14].ereby, deficit value was only used to compare meteorological drought and hydrological drought.According to the obtained drought propagation patterns, meteorological drought unnecessarily results in agricultural and hydrological droughts, and vice versa.In order to explain the 3 drought propagation patterns, we tried to find the possible driven factors of the formation of agricultural and hydrological droughts.en the relations between hydrological drought and the driven factors and between agricultural drought and the driven factors were established with the coefficients of determination 0.656 and 0.4, respectively.Although the accuracy is not high enough to make hydrological or agricultural drought forecast, it indeed interprets drought propagation patterns in the Luanhe river basin.
80th percentile is widely used as the threshold value in drought identification, which was also applicable in the Luanhe River basin [69].Different threshold value may have impact on drought characteristics and then on drought propagation patterns.And also droughts with duration longer than 5 days were counted and those shorter than 5 days were neglected, which also influenced propagation patterns.
erefore, the uncertainty of the selection of

Advances in Meteorology
threshold values and drought events should be further studied.Human activities are not considered here, which is another issue to be studied in the future.

Conclusions
is paper contributed to identification of drought propagation patterns in a watershed, and in particular, the explanation through the daily observed and simulated data.
e conclusions are as follows: (1) To remove the impact of human activities and obtain the soil moisture data, SWAT model was employed to get the hydrometeorological data under natural condition.e model performed very well in the Luanhe river basin, and the simulated results could be used for study of drought propagation.
(2) 3 types of drought propagation patterns were obtained, which are M-A/H, M-NA/NH, and NM-A/H.Meteorological drought unnecessarily resulted in agricultural/hydrological drought, and agricultural/ hydrological drought was not necessarily caused by meteorological drought.(3) Possible driven factors of agricultural and hydrological drought were found, and the relations between agricultural/hydrological drought and the driven factors were built by multiple regression method.ese relations could explain the 3 drought propagation patterns.
(a)), no intersection of the curves within the confidence intervals indicates the absence of a significant change point for the precipitation series.Together with the Pettitt test result (shown in Figure 2(b)), it is demonstrated abrupt changes in the annual precipitation did not occur from 1963 to 2012.

Figure 2 :
Figure 2: Results of (a) Mann-Kendall-Sneyers test and (b) Pettitt test for the annual precipitation during 1963∼2012.e fine dotted lines represent the critical value corresponding to the 95 % confidence level.

Figure 5 :
Figure 5: e exponential function between DS H and SD, DS M , E PD 10.

Figure 7 :
Figure 7: e scatters of simulated and observed SD.

Table 1 .
Validation of SWAT Model.e period before 1979 is undisturbed by human activities, especially without check dams and reservoir in Luanhe River basin.In the calibration period of 1963-1975, Ens and R 2 are greater than 0.65 and 0.65, respectively.Especially in monthly simulation, they are both 0.89.In validation period (1976-1979), both Ens and R 2 of monthly simulation exceed 0.90, and the values of daily simulation is greater than 0.74.According to Moriasi et al.
erefore years of 1963-1979 are considered as undisturbed Advances in Meteorology period for model calibration and validation based on daily discharges.Runoff data from Luanxian hydrological station are selected to calibrate and validate the SWAT model in the undisturbed period (1963-1979).eresults are listed in

Table 1 :
Results of calibration and validation of the River basin.

Table 2 :
Summary of drought characteristics of meteorological drought events based on the 5 different threshold values.

Table 3 :
Drought characteristics during the period of 1963-2012.

Table 5 :
Drought characteristics for each drought propagation pattern.

Table 6 :
Lag time from meteorological to agricultural/hydrological drought.

Table 7 :
Pearson correlation analysis between hydrological drought severity and each driven factor.

Table 9 :
Pearson correlation coefficients between SD and the driven factors.Figure 6: Scatter plots of the calculated DS H and observed DS H.

Table 10 :
e conditions of propagation from meteorological to agricultural drought.