Abstract

Based on monthly streamflow and precipitation data from 1960 to 2010 in the Jinjiang River Basin of China, Standardized Precipitation Index (SPI) and Standardized Streamflow Index (SSI) were used to represent meteorological and hydrological drought, respectively. The response of hydrological drought to meteorological drought under the influence of Shanmei reservoir was investigated. The results indicate that SPI and SSI have a decreasing trend during recent several decades. Monthly scales of SSI series have a significant decreasing trend from November to the following February and a significant increasing trend from May to July at Shilong hydrological station. There are three significant periodic variations with a cycle of 6-7 years, 11-12 years, and 20-21 years for annual scales of SSI series. SPI series have the same periodic variations before the 1980s, but they have not been synchronous with SSI since the 1980s at Shilong due to influences of Shanmei reservoir, especially at the periodic variations of 20-21 years. The variation of the lag time of hydrological drought in response to meteorological drought is significant at the seasonal scale. The lag time of hydrological drought to meteorological drought extends one month on average in spring, summer, and autumn but about three months in winter.

1. Introduction

Global warming has increased the frequency and intensity of extreme climate events such as drought [1]. The global annual economic loss owing to the drought ranges from 6 to 8 billion dollars and is higher than other meteorological disasters [24]. China is one of the countries that experiences the most frequent droughts and the most severe damage. Overall, the evolution of drought in China has exhibited a correlation with global drought [5, 6].

Drought typically is categorized as meteorological drought, agricultural drought, hydrological drought, and socioeconomic drought [7, 8]. Indices of meteorological drought are typically based on precipitation, and hydrological drought is based on streamflow. A meteorological drought may develop quickly, while a hydrological drought lags behind meteorological drought. The lag time depends upon watershed characteristics such as evapotranspiration capacity, the properties of the underlying surface (e.g., geology, soil, and vegetation conditions), and management factors (e.g., construction of reservoirs and dams) of the drainage basin [9, 10]. In recent decades, the construction of large reservoirs has shown a significant impact on surface hydrological processes and runoff mechanisms, and it also can interfere with the natural evolution of hydrological drought in the drainage basins [11, 12]. The evaluation of the response relationship between meteorological drought and hydrological drought impacted with large reservoirs is important for the realization of early warning and drought prevention as well as the mitigation of hydrological drought.

A drought index is a measurement that reflects the cause and extent of drought. In principle, a good drought index should be well defined, physically meaningful, based on conveniently collected data, and computationally stable. The Standardized Precipitation Index (SPI) [13] can meet these analysis requirements under different water resource conditions [13, 14]. Previous studies applying this method have been conducted in different regions (e.g., USA, China, Ethiopia, Korea, and Iran) [1421] and confirmed that this method is of practical significance. Having similar characteristics as SPI, Standardized Streamflow Index (SSI) [22], which uses streamflow as the characteristic index and evolved from SPI, has been widely applied to the hydrological drought studies [10, 19, 20].

Currently, two primary methods are used to investigate the response of hydrological drought to meteorological drought: the first is based on the simulation of hydrological models to investigate the relationship of the temporal-spatial dynamic changes of two types of drought [19, 20]. Because the definition of large reservoir parameters is ambiguous in hydrological models, if the simulation results are not properly verified, large errors may occur in the response evaluation results [19, 20]. The other method is the statistical analysis. Because the correlation analysis of the drought characteristics is based on measured hydrological and meteorological data, this method possesses high operability and has been widely applied [2326]. Wen et al. [23] used long term streamflow data to analyze the Murrumbidgee River Basin in Australia and showed that the regulation of upstream reservoirs alleviates the extent of hydrological drought in the downstream irrigation regions. In a study of the Diyala River Basin in Central Asia, Al-Faraj and Scholz [24] showed that the construction of upstream reservoirs significantly changed the downstream hydrological conditions and exacerbated the extent of downstream hydrological drought. Furthermore, in an investigation of the change in the upstream and downstream hydrological conditions of a drainage basin, Huang et al. [25] and Zhang et al. [26] showed that the correlation between a hydrological drought and a meteorological drought is closely related to the construction of hydraulic projects. These studies focused on revealing how the construction of large reservoirs has changed the hydrological conditions and the variation of hydrological drought. However, the changes in response of hydrological drought to meteorological drought under the effect of large reservoirs remain to be further investigated.

Therefore, taking the Jinjiang River Basin in the southeastern coast of China and selecting the Shilong hydrological station which is affected by a large reservoir and Anxi hydrological station which is not affected by a reservoir, a comparative analysis was made to explore the differences in drought response in this study. Based on the monthly streamflow data from 1960 to 2010 and the corresponding precipitation data of three meteorological stations, the changes in response of hydrological drought (SSI) to meteorological drought (SPI) under the impact of large reservoirs (Shanmei Reservoir) in the drainage basin were investigated through a comparative analysis using the nonparametric Mann-Kendall test [27, 28], a wavelet analysis [29], and a Pearson correlation analysis.

2. Overview of the Study Area

Jinjiang River Basin is located in Quanzhou City, Fujian Province, which is the starting point of China’s 21st century “Silk Road on the Sea.” The drainage area is 5629 km2, the total river length is 302 km, and the length of the mainstream is 182 km. Jinjiang River Basin is the third largest river in Fujian Province. Dongxi and Xixi Rivers are tributaries in the drainage basin, and these two tributaries converge downstream at Shilong hydrological station. The drainage area controlled by Shilong hydrological station is 5024 km2 and accounts for 89.60% of the total Jinjiang River Basin. The drainage basin upstream of this hydrological station is selected as the study area (Figure 1). Climate in the drainage basin is characterized as South Asia humid subtropical monsoon climate with a mean annual rainfall of 1868 mm (1960–2010) and a mean temperature of around 20°C. The area is characterized by significant monsoon climate, with clear wet and dry seasons [30]. Precipitation from March to September in the wet season accounts for 83.8% of total annual precipitation, while precipitation in the five months of the dry season, from October to February, accounts for only 16.20%. This area is one of the poorest water resource areas in Fujian Province with annual mean surface water resources of 9.63 × 109 m3, which is only 8% of the Fujian Province’s average (http://www.fjwater.gov.cn/). This area frequently experiences hydrological droughts [31]. To meet the regional water usage requirement, several large and small reservoirs have been constructed, among which Shanmei Reservoir located in the middle tributary of Dongxi Basin has a water area of 1,023 km2 and a total storage capacity of 6.55 × 108 m3. The reservoir is with multiple-annual regulating functions such as irrigation, power generation, and water supply. It is the only large multifunction reservoir in Jinjiang River Basin. The dam filling was completed in 1972, while the protection and expansion were conducted in 1979. At the end of 1982, the construction of the reservoir was fully completed. In addition, when the Longmentan Water Diversion Project was completed in 1989, the drainage area upstream of the diversion dam was 360 km2, and the annual water diversion was 4.06 × 108 m3. Using a transbasin diversion, this project introduced or delivered the Dazhangxi mainstream streamflow of the Min River Basin to the lateral lake, Huyangxi, of Dongxi, which is a tributary of Jinjiang River Basin. The streamflow from Min River Basin then converged to the Shanmei Reservoir and was transported to the downstream Shilong hydrological station by regulation of the reservoir.

3. Materials and Methods

3.1. Materials

Hydrological and meteorological data were obtained from the Fujian Provincial Hydrology and Water Resources Survey Bureau and the Fujian Provincial Meteorological Bureau. The data included the monthly streamflow of two hydrological stations (Shilong and Anxi) between 1960 and 2010. Shilong station is located at the confluence of the Dongxi and Xixi tributaries. The streamflow was affected by Longmentan Water Diversion Project in the Dongxi tributary and the Shanmei Reservoir. Considering the relatively small study area and the even distribution of meteorological stations, the monthly precipitation was obtained from the calculation of the arithmetic mean of the three meteorological stations (Nan’an, Anxi, and Yongchun) from 1960 to 2010. According to the four-season classification method in meteorology for China, March–May is the spring, June–August is the summer, September–November is the autumn, and December–February is the winter.

3.2. Methods
3.2.1. Standardized Precipitation/Streamflow (SPI/SSI)

Standardized Precipitation Index (SPI) is used to represent meteorological drought, while Standardized Streamflow Index (SSI) is for hydrological drought. The 51-year monthly precipitation and streamflow data from 1960 to 2010 were used to calculate SPI and SSI at different time scales. Therefore, the suitable probability distribution types for the variables were first determined, and then SPI and SSI were obtained after normal standardization. Calculated SPI and SSI can be used to monitor hydrological and meteorological drought characteristics under multiple temporal and spatial scales and thus provide a reliable method for obtaining the response patterns of hydrological drought to meteorological drought [1620]. Shukla determined that compared to other probability density functions the streamflow series conforms more to a Log-Normal function [21]. Angelidis et al. [15] investigated and compared the fitting effect of Gamma, Normal, and Log-Normal functions on precipitation series at different time scales. Their results showed that the Log-Normal function had good applicability on the precipitation series at time scales greater than a month. In order to confirm the applicability of this conclusion in Jinjiang River Basin, the Kolmogorov-Smirnov (K-S) test method [32, 33] is performed to examine the goodness of fit of Log-Normal distribution for the precipitation and streamflow at 4 time scales (1-month, 3-month, 6-month, and 12-month), as listed in Table 1. The Log-Normal fits better at different time scales for streamflow overall and at longer time scales (6-month and 12-month) for precipitation, which is consistent with the conclusions previously obtained for the Xijiang River Basin (southeastern China) [20]. Therefore, Log-Normal function was selected as the probability distribution type for variables. The detailed calculation procedures can be referred to in previous studies [15, 18].

3.2.2. Mann-Kendall Algorithm

Mann-Kendall (M-K) nonparametric trend test method is recommended by the WMO (World Meteorological Organization), as a nonparametric approach for calculating the trends of hydrometeorological series [27, 28]. M-K test was used to analyze the trend of SSI and SPI series in this study. The detailed calculation method is based on the test statistic , which is expressed as follows [34]:where , are series of data values, and the sign function is defined as

The variance of is calculated as follows:

” represents the length of data series. Mann and Kendall have suggested that when , the statistic is approximately normally distributed with zero mean [27, 28]. Then, the standardized test statistic is calculated by

The standardized M-K nonparametric statistic follows the standard normal distribution with mean of zero and variance of one. If , a significant trend is accepted to exist in the time series, where is the value of the standard normal distribution with an exceeding probability of . The significance level is set as 0.05 and 0.01 in this study, and their corresponding value is 1.96 and 2.38, respectively. When it means a significant increasing trend and means the increasing trend is extremely significant, while and indicates a significant and extremely significant decreasing trend, respectively [34].

3.2.3. Continuous Wavelet Transform (CWT)

Continuous wavelet transform (CWT) [29] is a useful mathematical method to study the information in time-frequency domains from the perspective of multi-time scales, which has been widely used to identify periodic oscillations of signals [6, 35]. And the basis of a Morlet wavelet (Φ) is given aswhere is the Morlet wavelet function, is defined as imaginary symbol of a complex, is the dimensionless frequency, and is the dimensionless time parameter. Given a time series and , then the CWT can be computed by where is wavelet variance; is the wavelet transform coefficient, which can represent the characteristics of the signal change at different times scales; and are called a scaling parameter and a translation parameter, which measures the degree of compression or scale and determines the time location of the wavelet, respectively; and is the complex conjugate of . In the derivation process of wavelet filtering, it is assumed that the signals are infinite; however, the signals that will be processed are finite in practical progress; therefore the border of composition signals will have great error, which is called “boundary effect” [36]. Thus when the CWT was used to eliminate the “boundary effect” on the two boundaries of the variable (i.e., the time series), symmetry-extending was applied to the time series. According to analysis of the wavelet coefficients at time-frequency domains from the perspective of multi-time scales, we can obtain the differences of periodic components between two time series.

3.2.4. Pearson Correlation Coefficients (PCC)

Pearson correlation coefficient (PCC) was used to quantitatively measure the correlation of SSI and SPI at different time scales to compare the evolution characteristics and the relationship between meteorological drought and hydrological drought. The PCC is defined as [37]where and represent the two time series of values, with . PCC values range from −1 to 0 indicates that sets are negatively correlated, while values from 0 to 1 mean they are positively correlated. In order to analyze how hydrological drought responds to meteorological drought, we used PCC method to analyze the time lag between hydrological drought and meteorological drought by using the monthly SSI (SSI-1) series and SPI at different time scales (SPI-1 to SPI-12) for each month, by finding the differences between the drought durations. For example, SSI-1 and SPI-3 have higher PCC values; this indicates that the time response of hydrological drought to meteorological drought is about two months. The higher PCC suggests that this response time is more sensitive and vice versa [38].

4. Analysis of Results

4.1. Response of Hydrological Drought to the Evolution of Meteorological Drought
4.1.1. Response of Trend Change

The M-K trend test approach was utilized to analyze the trends of monthly SSI values and total series at Shilong, Anxi hydrological station, and SPI series in study area. The M-K statistic values of SSI and SPI in the total series on a monthly scale are all greater than zero (Table 2). Shilong station passed the 0.01 significance test, Anxi station passed the 0.05 significance test, and SPI total series passed the 0.05-significant-level test. M-K values of SSI and SPI exhibited a trend of pronounced increase, indicating that the Jinjiang River Basin did not exhibit a significant drought trend over 51 years. This conclusion is consistent with findings by Y. Chen and X.-W. Chen [39], who found a trend of general development towards flood in Fujian Province as observed in 1960–2006 using SPI. According to Yang et al. [40] and Wu [41], a rapid temperature increase occurred in Fujian in recent decades, combining with a remarkable increase in precipitation, which is more pronounced in the southeastern coast of Fujian. These results suggest a good consistency between the nonsignificant trends of droughts in the Jinjiang River Basin and warming-wetting climate in the southeastern coast of Fujian.

The results of monthly M-K trend test indicate that the SSI and SPI sequence values have statistically significant decreasing in May–July, while the other months exhibit an increasing trend. In general, the dry season (October to February of the following year) exhibits an apparent trend of increasing, which is more pronounced in Shilong station than Anxi station, particularly in winter (December to February of the following year). Comparing SPI value’s changes in the drainage basin, the evolution of SSI values at the two hydrological stations did not coincide with or correspond to SPI values. This result is particularly evident at Shilong station. Overall, the monthly SSI values in the Jinjiang River Basin have an apparent intra-annual variation characteristic.

4.1.2. Decadal Change Response

The cumulative distribution functions (CDFs) of SSI-1 at Shilong and Anxi stations for different decades (1960–1969, 1970–1979, 1980–1989, 1990–1999, and 2000–2010) were plotted, represented by 1960s, 1970s, 1980s, 1990s, and 2000–2010, respectively. The decadal evolution of hydrological drought frequency in the drainage basin was compared between two hydrological stations (Figure 2). All the CDFs passed K-S tests at chosen confidence level of 99%. The frequency variation of hydrological drought for 1960s, 1970s, and 1980s at Anxi station is consistent with Shilong station for the drainage basin by the distance between CDFs. In contrast, the variation between Shilong and Anxi station shows inconsistency after 1980s. In particular, the frequency of hydrological drought in the 1990s significantly decreased.

4.1.3. Response of Periodic Change

Morlet wavelet was used to conduct a periodic analysis of SPI series in the drainage basin on the annual scale. The real part of the wavelet transform coefficient, the modulus square time-frequency distribution, and the real part of the wavelet transform coefficient of the main period were obtained, as shown in Figure 3. For convenient analysis, the positive of the real part of the wavelet coefficient was characterized as a solid line to indicate that the SPI value at this stage is high (e.g., prone to wetness). In contrast, the negative of the real part of the wavelet coefficient was characterized as a dashed line to indicate that the SPI value at this stage is low (e.g., prone to drought).

From Figure 3(a) we can see that the annual scale of SPI time series includes many strong and weak periodic oscillations at different time scales that reflect the characteristics of alternating variation between wet and dry. The modulus square time-frequency distribution of wavelet coefficient (Figure 3(b)) shows that three apparent periodic oscillations exist, 6-7 years, 11-12 years, and 20-21 years for the annual SPI time series, among which the oscillation of the 20-21-year time scale is more significant than others but is weaker after 1980s. These results suggest good consistency with Minjiang River Basin, which is similar to the geographical features of the Jinjiang River Basin by using the EMD method [39].

Figure 3(c) shows SPI main periodic oscillations; the periodic oscillation of 20-21 years has been gradually weakened but the oscillation in 11-12 years has been gradually strengthened, suggesting that under global warming the occurrence or shifting of drought-flood in the drainage basin became more frequent. As a result, measures are necessary to prevent the risks of drought-flood and secondary disasters.

The calculation of the real part of wavelet transform coefficients for Shilong and Anxi stations at the same time scales shows that three apparent periodic oscillations existed around 6-7 years, 11-12 years, and 20-21 years for SSI series of the two hydrological stations. Figures 4(a), 4(b), and 4(c) show the comparison of the main period wavelet coefficient of SSI at two hydrological stations with SPI. The wavelet coefficient of SSI lags behind SPI, indicating the time lag between hydrological drought and meteorological drought. For the periodic oscillation around 6-7 years, the wavelet coefficient of the two hydrological stations keeps synchronous fluctuation with SPI in the drainage basin. The nonsynchronization was evident for the periodic oscillation around 11-12 years, particularly after 1980s. In contrast, the SSI of Anxi station remained almost synchronous with SPI for the periodic oscillation around 20-21 years, but the nonsynchronization for Shilong station was evident at the end of 1980s. To better explain the changes before and after 1980s about the main periodic oscillations, we analyzed the correlations between SPI and SSI wavelet transform coefficients for Shilong and Anxi station in different periods. From Table 3 we can clearly see that the correlations between SPI and SSI of main period wavelet coefficient at Shilong station are lower than at Anxi station, especially during 1980–2010. The higher correlation of SSI and SPI supports the greater synchronization between them.

4.2. Response Time of Hydrological Drought to Meteorological Drought at the Shilong Station

Previous studies have shown that SPI exhibits multiple time scale characteristics [13, 15, 17, 18]. An appropriate choice of SPI time scale can effectively indicate dry and wet conditions and also can indicate the response time of hydrological drought to meteorological drought [1720]. The above analysis shows that the lack of synchronization between hydrological drought and meteorological drought at Shilong station became more evident after the 1980s. The dam and impoundment expansion project for Shanmei Reservoir was completed in 1982, and that year was taken as the breakpoint between the period with no reservoir effect (1960–1982) and the period with a reservoir effect (1983–2010). The relationships between monthly scales SSI and SPI on different time scales (from SPI-1 to SPI-12) for the two periods (1960–1982 and 1983–2010) were analyzed at different months (from January to December). The results illustrate the degree of consistency of SSI and SPI at the time scale: greater PCC values are more consistent and indicate that this response time is more sensitive and vice versa. Overall, Figures 5(a) and 5(b) show that the response time of hydrological drought to meteorological drought exhibited seasonal characteristics. However, compared with the period with no reservoir effect, correlations between monthly scale SSI and short time scales SPI (e.g., SPI-1 and SPI-2) at each month become weaker, which indicates that the response time of hydrological drought to meteorological drought in the period with a reservoir effect at Shilong station becomes longer, in general. The response times in spring, summer, autumn, and winter have changed from 1–7, 1-2, 1–3, and 3–5 months to 2–7, 1–3, 1–4, and 6–9 months, respectively.

The years with significant drought were further selected to verify the changes of the correlation between SSI and SPI for Shilong station under the regulation of the Shanmei Reservoir. Because this region is one of the regions with the most scarce water resources in Fujian Province, the drought in this drainage basin is close to the drought in the entire province [31]. Due to limited data in the basin, the years with severe drought throughout the entire province were selected for this analysis. According to the Fujian Provincial Bureau of Statistics [30] and related materials (http://www.doc88.com/p-652201703763.html) over 340000 ha of land was affected by drought in 1950–1995; the recorded years with severe drought were 1963, 1971, 1977, 1983, 1986, 1991, 2003, and 2009. The monthly SPI and SSI of the above 8 years were extracted for analysis (due to the lag effect of hydrological drought to meteorological drought, advanced and delayed 6 months were extended for the severe drought year). As seen in Figure 6, the fluctuation of SSI corresponded to SPI in 1963, 1971, and 1977 before the 1980s. It exhibits an excellent correlation with an average Pearson coefficient of 0.82 (Table 3). However, for the years after the 1980s (i.e., 1983, 1987, 1991, 2003, and 2009), the average lag time significantly increased. It is particularly obvious for the years 1987, 1991, and 2003. In addition, the correlation of SSI and SPI was poor compared with that before the 1980s, with an average Pearson coefficient of 0.66 (Table 4). It indicates that the function of regulation by Shanmei Reservoir affects the correlation between SSI and SPI, resulting in a change in the response time of hydrological drought to meteorological drought at Shilong hydrological station.

4.3. Discussion

Although the evolution of hydrological drought at two hydrological stations showed a similar decreasing trend (decreased droughts) and the hydrological drought was generally correlated with meteorological drought, M-K test results for SSI and SPI on a monthly scale showed that the evolution of hydrological drought in the drainage basin exhibited a nonsynchronized change with meteorological drought. This change is more evident at the Shilong station than the Anxi station; the M-K values in Shilong hydrological station were significantly higher than at Anxi station for November to February but lower for June to July. The frequency variation of hydrological drought of each decade at Anxi station is almost synchronized with meteorological drought for the drainage basin, while Shilong station exhibited a different and nonsynchronized pattern. The periodic oscillation of Shilong station was almost consistent with the corresponding period of Anxi station as well as the drainage basin before the end of 1980s. However, the difference became significant on the 20-21-year cycle scale.

The SPI is based on precipitation deficits and therefore is purely dependent on climate anomalies [13, 42]. The evolution and periodic oscillation of hydrological drought at two hydrological stations (Shilong and Anxi station) in Jinjiang River Basin are consistent with meteorological drought in general, which indicates that climate change plays a role in the evolution of hydrological drought. Under a similar climate condition, the inconsistency of SSI series and SPI series at the corresponding time scale was more significant for Shilong station than Anxi station, which is strongly related to human activities (e.g., dam construction and operation). Overall, the human activities in the Jinjiang River Basin primarily included a change in land use and the construction of hydraulic projects. Lin et al. [43] investigated the land use data of the Jinjiang River Basin in 1985 and 2006 and applied the Soil and Water Assessment Tool (SWAT) model to simulate the response of flood-drought streamflow to land use changes. Their results showed that due to a decrease in woodlands and an increase in agricultural land, under the same precipitation conditions, no significant difference of the streamflow rate change was observed between the two stations. Therefore, the evolution of hydrological drought at Shilong station was caused by other human activities, such as the construction of large reservoirs. The completion of the Shanmei Reservoir with multiple-annual regulating functions, including the completion of dam protection and impoundment expansion projects in 1979–1982, has produced a significant impact on the interannual and particularly the intra-annual variations in streamflow of the downstream Shilong station. Meanwhile, the Longmentan Water Diversion Project, completed in 1989, introduced the streamflow from the Dazhangxi mainstream of the Min River to Huyangxi of upstream Dongxi of the Jinjiang River by transbasin diversion. This part of streamflow entered into the Shanmei Reservoir for unified regulation. The direct result of these projects is the significant effect on the interannual increase and intra-annual stability of streamflow for Shilong station. Figure 7 shows the change trends of the coefficient of variation (Cv) (Figure 7(a)) and concentration (Cr) (Figure 7(b)) of the intra-annual streamflow for Shilong and Anxi stations, respectively. The intra-annual streamflow change stabilized, and the decreasing trend was less affected by the reservoir than Anxi station.

5. Conclusions

In general, the evolution of hydrological drought was associated with meteorological drought at the two hydrological stations in the Jinjiang River Basin, and they both exhibited a decreasing trend. However, the hydrological drought decreased during November to February of the following year and increased during May to July, showing a nonsynchronized change with that of meteorological drought. This result was more evident at Shilong station than Anxi station. The primary reason is that an important human activity in the drainage basin, the multiple-annual regulation of Shanmei Reservoir, has significantly affected the runoff at Shilong station and correspondingly affected the evolution of hydrological drought.

Due to the impact of Shanmei Reservoir, the response time of hydrological drought to meteorological drought for Shilong station exhibited a changed seasonal variation (i.e., spring, summer, and autumn have extended 1 month on average, and winter has extended 3-4 months on average).

To a large extent, the regulation function of large reservoirs usually can change the temporal and spatial evolution of runoff and hydrological drought in the drainage basin. Therefore, further investigation of the response pattern of hydrological drought to meteorological drought under the impact of large reservoirs is crucial to the realization of early warning, drought prevention, and mitigation.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

The study was supported by the Science and Technology Major Program of Fujian Province (Grant no. 2015Y4002) and the National Natural Science Foundation of China (Grant nos. 41301031 and 41501106). The hydrological data is provided by Fujian Provincial Hydrology and Water Resources Survey Bureau, and the meteorological data is provided by the Fujian Provincial Meteorological Bureau. The authors wish to thank Greg Mierle for kind proofreading.