Tendency of Runoff and Sediment Variety and Multiple Time Scale Wavelet Analysis in Hongze Lake during 1975–2015

Hongze Lake plays a key role in flood and waterlogging prevention, analyzing the variation process and characteristics of multi-time scales will have a great practical significance to water resources management and regulation in the Huaihe River basin of China. This research proposed a combinatorial mutation test method to study the interannual variation trends and change points of runoff and sediment flowing into and out of Hongze Lake during the period 1975–2015. It is concluded that the annual variation trend of the inflow and outflow runoff time series is consistent, with no obvious decreasing trend and change point, while the inflow and outflow sediment time series showed a decreasing trend, and the change point was 1991. Then, the runoff and sediment time series were analyzed by the wavelet method. The results showed that the time series has multi-time scale characteristics. The annual inflow runoff and sediment would enter into the dry period in a short time after 2015, and both would be in the valley floor stage. Among the influencing factors, the variation of rainfall in the basin was the main factor affecting the runoff variation. Changes in heavy rainfalls pattern, the construction of hydraulic engineering projects, and land use/cover change (LUCC) are the main reasons for the significant decrease and mutation variation of inflow sediment.


Introduction
Runoff plays a vital role in the hydrologic cycle. It is of great significance to explore the changing trend of water flow in the field of water resources [1]. Hydrologic events such as variations in runoff and sediment are partly results of variations in the hydrologic cycle [2], these variations can be divided into two categories: trend and shift [3]. Hydrological regime-changing has important significance in planning and designing of hydraulic structures along the lake, understanding these variations plays an important role in lake engineering management [4].
Flood and drought have always been the two most important factors restricting social stability and economic development in the Huaihe River basin. At present, most studies on the hydrological cycle variations of the Huaihe River basin focused on the mainstream and branches water system, Zhou et al. [5] studied the dynamic characteristics of flood time series in the Huaihe River basin based on chaos theory. Lin et al. [6] simulated the heavy precipitation and flood disaster in the Huaihe River basin by using the framework of atmospheric-hydrological simulation. Zhu et al. [7] analyzed the impacts of climate variability and anthropogenic activities on runoff in the Huaihe River basin based on the 1964-2010 long series of hydrological time series data. Ye et al. [8] and Yu et al. [9] explained the spatial characteristics and variation trends of precipitation in the Huaihe River basin.
Hongze Lake as one of the four largest freshwater in China, located in the junction of middle and lower Huaihe River basin, its storage, and drainage capacity is crucial to flood control in the middle and lower Huaihe River. Affected by climate variability and anthropogenic activities, the runoff and sediment flowing into and out of Hongze Lake affected the flood regulation and storage capacity. Many scholars have focused on the variation of rainfall and flood disasters in the Hongze Lake area. Jiang et al. [10], based on the historical water level and discharge data of Hongze lake, analyzed the long and short time series of flood peak discharge into the lake, and elaborated on the characteristics of flood disaster. Yin et al. [11] analyzed the relationship between the maximum water level variation and climate variability and anthropogenic activities from 1736 to 2005 in Hongze Lake. Duan et al. [12] studied the influence of anthropogenic activities (artificial sand mining) on the hydrological regime of Hongze Lake. There is short of research on the variation trend and causes of runoff and sediment in the long time series of Hongze Lake.
In the context of global warming, the frequency and intensity of extreme precipitation events are increasing [13], and this directly leads to the constant variation of the hydrological time series frequency domain characteristics. However, the frequency domain feature is difficult to identify by traditional hydrological analysis methods. The variation of runoff and sediment characteristics directly affects lake regulation, flood control, and navigation. Therefore, it has prominent theoretical and practical significance to study the variations trends and causes of runoff and sediment in Hongze Lake. Also this is beneficial to flood and waterlogging disaster prevention and water resources management.
Because time series analysis has fewer assumptions and uncertainties, the historical data time series analysis is more helpful to understand the variation behavior of runoff compared with the simulation model [14]. There are many kinds of discriminant methods for trend variation, among which the statistical test can analyze the trend variation significance level and has a rigorous statistical basis. The Mann-Kendall (M-K) [15,16] trend test can deal with non-normal distribution and explain the existence of missing values and non-stationarity in data, these advantages make it a widely used tool for analyzing hydrological variables variation trends [17,18]. However, the M-K test has poor processing ability for nonlinear data, so other analytical methods are needed. Decomposing the original time series will beneficial to observe the variation trend of hydrological data in different periods, so the wavelet analysis method was chosen for its ability to process nonlinear data and decompose the original time series [19]. Many scholars have studied the characteristics of runoff and sediment in different rivers and lakes by using wavelet analysis [20][21][22], however there is a lack of research on the multi-time-scale variation characteristics of runoff and sediment time series in Hongze Lake by the wavelet analysis method.
There are many methods to study the change point of hydrological time series, including the Mann-Whitney test [23], the Kruskal-Wallis test [24], the Pettitt test [25], and the Student's t-test [26]. Among these tests, Pettitt test is widely used in hydroclimatology because of its high accuracy [27]. In the previous studies [28][29][30], the purpose of this method is to identify the change points of the variable and analyze the mutation causes. However, due to the difference of presuppositions and applicable conditions, the results are often different, so the analytical method of hydrological time series change point needs to be improved.
Also, the research on the influence of climate variability and anthropogenic activities on water circulation and water resources is a hot topic in climatology and hydrology. In recent years, extreme weather events such as a rainstorm, flood, strong storm surge, drought, high temperature and heat wave have become more and more intense [31,32]. Climate variability, mainly characterized by global warming, affects the regional land water circulation by changing the atmospheric temperature and precipitation, and drives the changes of hydrological factors such as runoff and sediment [33,34]. Anthropogenic activities mainly include the construction of hydraulic engineering, water and soil conservation measures, unreasonable farming, etc., of which the construction of hydraulic engineering is a typical representative, sediment is trapped in the reservoir, resulting in reduced sediment transport downstream of the river [35]. As an artificial "disturbance system", LUCC is not only an important part of regional environmental evolution, but also one of the main causes of physical process changes on the surface [36,37]. Ervinia et al. [38] found that the land use change enhanced the impact of climate change on runoff change. Dao et al. [39] found that climate change from 1978 to 2000 had a stronger impact on the hydrological process of the Vietnam river basin than the land use change. Apollonio et al. [40] showed the impact of land use on the hydrological response is closely related to the watershed scale. Bussi et al. [41] simulated the effects of future climate and land use changes on the transport of suspended sediment in the Thames River.
The purpose of this research is to improve the accuracy of the trend detection technology of time series, then clarify the variation trend, change points and periodicity of the runoff and sediment time series of the Hongze Lake, and further analyze the influencing factors from climate variability and anthropogenic activities. In this paper, the hydrologic stations of the surrounding water system of Hongze Lake are selected to explore the variation trend and change points of runoff and sediment flowing into and out of the Hongze Lake basin from 1975 to 2015. A combinatorial mutation test method based on the M-K test and Pettitt test is proposed to identify the catastrophe points. Based on statistical analysis, wavelet analysis is used to analyze the spatiotemporal variation characteristics of the hydrological time series, then predict its variation trend in a short time in the future. By analyzing the causes of the change point, revealing the impacts of climate variability and anthropogenic activities on runoff and sediment variations in Hongze Lake. It is of great significance for the management and economic development to understand the characteristics, variation trends, and influencing factors of runoff and sediment in Hongze Lake scientifically.

Research Areas and Data
This paper mainly takes Hongze Lake as the research area, which is located at the junction of the middle and lower Huaihe River basin of China. The hydrological data are the annual runoff and sediment time series of ten hydrological stations flowing into and out of Hongze Lake from 1975 to 2015. The surrounding water system of Hongze Lake [42] is shown in Figure 1.
Water 2020, 12, x FOR PEER REVIEW 3 of 18 sediment transport downstream of the river [35]. As an artificial "disturbance system", LUCC is not only an important part of regional environmental evolution, but also one of the main causes of physical process changes on the surface [36,37]. Ervinia et al. [38] found that the land use change enhanced the impact of climate change on runoff change. Dao et al. [39] found that climate change from 1978 to 2000 had a stronger impact on the hydrological process of the Vietnam river basin than the land use change. Apollonio et al. [40] showed the impact of land use on the hydrological response is closely related to the watershed scale. Bussi et al. [41] simulated the effects of future climate and land use changes on the transport of suspended sediment in the Thames River. The purpose of this research is to improve the accuracy of the trend detection technology of time series, then clarify the variation trend, change points and periodicity of the runoff and sediment time series of the Hongze Lake, and further analyze the influencing factors from climate variability and anthropogenic activities. In this paper, the hydrologic stations of the surrounding water system of Hongze Lake are selected to explore the variation trend and change points of runoff and sediment flowing into and out of the Hongze Lake basin from 1975 to 2015. A combinatorial mutation test method based on the M-K test and Pettitt test is proposed to identify the catastrophe points. Based on statistical analysis, wavelet analysis is used to analyze the spatiotemporal variation characteristics of the hydrological time series, then predict its variation trend in a short time in the future. By analyzing the causes of the change point, revealing the impacts of climate variability and anthropogenic activities on runoff and sediment variations in Hongze Lake. It is of great significance for the management and economic development to understand the characteristics, variation trends, and influencing factors of runoff and sediment in Hongze Lake scientifically.

Research Areas and Data
This paper mainly takes Hongze Lake as the research area, which is located at the junction of the middle and lower Huaihe River basin of China. The hydrological data are the annual runoff and sediment time series of ten hydrological stations flowing into and out of Hongze Lake from 1975 to 2015. The surrounding water system of Hongze Lake [42] is shown in Figure 1.

M-K Test
In addition to trend analysis, the M-K test theory can also be used for the change point test, which is effective for verifying the variation from a relatively stable state to another state. For the time series X containing n samples, we construct an order column in the equation [15,16]: where Under the assumption of random independence of time series, statistics are defined as the equation: where, S k is the cumulative value of the number of values at time i is greater than that at time j, UF 1 = 0; E(S k ), Var(S k ) are the mean and variance of cumulative S k respectively, when x 1 , x 2 , . . . , x n are independent of each other, E(S k ) and Var(S k ) have the same distribution, then: UF k is the standard normal distribution, which is a statistical time series calculated by time series X. At a given level of significance α, according to the normal distribution table, when UF k > U α , there are obvious trend variations in the time series.
Inverting the time series X and repeating the process, setting UB k = −UF k (k = n, n − 1, . . . ,1), UB 1 = 0. To calculate the UF k and UB k , if UF k is greater than 0, the time series shows an upward trend; if the UF k less than 0, it shows a downward trend. At the significance level α=0.05, the mutation probability increases when UF k exceeds the critical value ±1.96. Within the confidence interval, if there is an intersection point between curve UF k and UB k , it is a possible catastrophe point. However, the M-K test does not applicable to the time series with multiple or multi-scale mutations time series, that is, when there are multiple intersection points in the confidence interval, there may be pseudo catastrophe points, and the clutter points need to be removed [43].

Pettitt Test
The basic principle of the Pettitt test is to assume that the sample is divided into two-time sub-series arbitrarily and the mean is the same. That means the random variables time series is divided into two parts arbitrarily, x 1 , x 2 , . . . , x to and x to+1 , x to+2 , . . . , x T . If the distributions of the random variables in the two parts are F 1(x) and F 2(x) , and F 1(x) F 2(x) , the catastrophe point is supposed to happen at t o . The statistics are defined as the equation [25]: When the time series follow a continuous distribution, statistic U t,T is obtained by the equation: Water 2020, 12, 999

of 18
For t = 2, 3, . . . , T, The point satisfying Equation (8) is considered to be the most likely position for the mutation.
The corresponding significance probability can be obtained by the equation: The Pettitt test can only identify one change point in a long time series, but there may be more than one change point in a hydrological time series. In the research, all the change points obtained using the method proposed by Li et al. [44]. The first-order change point was identified by the Pettitt test. Then, the hydrological long time series is divided into two parts based on the first-order change point, and the new change points are tested until all the change points are found.

Combinatorial Test
To obtain a more accurate change point, the research proposed a combinatorial mutation test method based on the M-K test and Pettitt test. The test process is as follows: ( This method uses the Pettitt test to offset the defect that the M-K test cannot identify the true and false change points, compared with the single Pettitt test, the combinatorial test can identify the change points and regions accurately, ensuring that the change points are significant, so that all significant points can be identified more scientifically and accurately.

Wavelet Function
The basic idea of wavelet analysis is to use a cluster of wavelet functions to represent or approximate a certain signal or function. The wavelet function is the key to wavelet analysis. Wavelet function is ϕ(t)∈L 2 (R) and it satisfies Equation (10).
where ϕ(t) is the base wavelet function, which can form a cluster function system by scaling and shifting on the time axis: where ϕ a,b (t) is the daughter wavelet; a is the scale factor, which can reflect the period length of the wavelet; b is the translation factor, which can reflect the shift in time.
At present, many wavelet functions can be used [44], for example, Mexican hat wavelet, Haar wavelet, Morlet wavelet, and Meyer wavelet. Among, the Morlet wavelet is widely used to identify periodic oscillations of the real life signals, which can detect the time-dependent amplitude and phase Water 2020, 12, 999 6 of 18 for different frequencies [45,46], it is a very accurate method for obtaining the periodicity [47]. Given the hydrological characteristics of the runoff and sediment time series, a continuous wavelet transform should be used to detect the existence of oscillations and their period, so the research selected Morlet wavelet function to analyze the multi-time scale variations of the runoff and sediment time series in Hongze Lake.

Wavelet Transform
Let L 2 (R) represents the measurable square-integrable function space defined on the real axis, if the function f (t)∈L 2 (R) satisfies Equation (12).
Then, for a given energy-limited signal f (t)∈L 2 (R), the continuous wavelet transform is defined as the equation: where

Wavelet Variance
For discrete-time series, the wavelet variance can be calculated according to the equation: where W(a,x j ) is the wavelet transform coefficient with scale factor a and translation factor x j . For complex wavelet function W 2 (a,x j ), it is the square of the wavelet transform coefficient module; n is the total number of wavelet transform coefficients obtained under scale a. According to Equation (14), the wavelet variance Var(a) is equal to the mean value of the square of the corresponding wavelet transform coefficient W 2 (a,x j ) of each sample at the scale a. The wavelet variance graph is the process of wavelet variance changing with time scale a, which can reflect the energy distribution of signal fluctuation in different time scales a. Therefore, the wavelet variance diagram can be used to determine the main period of the studied signal.

Runoff and Sediment Trend Analysis
From the analysis of hydrological time series, the average annual runoff and sediment flowing into and out of Hongze Lake fluctuated in waves from 1975 to 2015, as shown in Figure 2. There is little difference between the annual inflow and outflow runoff, while the annual inflow sediment is greater than the outflow except for some special years, the lake keeps in deposition state. The variation of runoff and sediment is closely related, which means that the sediment increases with the runoff, and when the runoff decreases, the sediment will decrease at the same time.
From the analysis of hydrological time series, the average annual runoff and sediment flowing into and out of Hongze Lake fluctuated in waves from 1975 to 2015, as shown in Figure 2. There is little difference between the annual inflow and outflow runoff, while the annual inflow sediment is greater than the outflow except for some special years, the lake keeps in deposition state. The variation of runoff and sediment is closely related, which means that the sediment increases with the runoff, and when the runoff decreases, the sediment will decrease at the same time. The M-K test is conducted on the runoff and sediment time series of Hongze Lake from 1975 to 2015, the statistics are shown in Table 1. Taking 1975 as the starting point of the calculation, the calculated M-K statistics positive series of the annual mean inflow and outflow runoff are shown in Figure 3. The M-K test is conducted on the runoff and sediment time series of Hongze Lake from 1975 to 2015, the statistics are shown in Table 1. Taking 1975 as the starting point of the calculation, the calculated M-K statistics positive series of the annual mean inflow and outflow runoff are shown in Figure 3. The inflow and outflow runoff M-K statistics are -0.98 and −0.60 during 1975-2015, and the decreasing trend is not obvious. While the M-K statistics are close to the critical value of −1.96 during 1994-2005, and the decreasing trend is relatively significant. However, during the whole study period, the M-K statistics are all lower than the significant level, with no significant variation trend, and the variation of runoff and sediment flowing into and out of Hongze Lake in normal fluctuation. The

Analysis of Wavelet Transform
Because the inflow variation law of runoff and sediment is similar to the outflow, so only the inflow runoff and sediment time series analyzed by wavelet theory. The real part contour map of wavelet coefficients can reflect the periodic variation law of the hydrological time series on various time scales [48], as shown in Figure 6a,b. The solid line shows that the real part of the wavelet

Analysis of Wavelet Transform
Because the inflow variation law of runoff and sediment is similar to the outflow, so only the inflow runoff and sediment time series analyzed by wavelet theory. The real part contour map of wavelet coefficients can reflect the periodic variation law of the hydrological time series on various time scales [48], as shown in Figure 6a,b. The solid line shows that the real part of the wavelet coefficient is positive, which represent the flood season of Hongze Lake, that's the orange area in the figure. The dotted line shows that the real part of the wavelet coefficient is negative, which represents the dry season of Hongze Lake, that's the light blue area in the figure. The thick solid line shows the zero isolines, that is, the borderline between rich period and dry period.
research domain, and with a global character.
The modulus square value of the wavelet coefficient is equivalent to the wavelet energy spectrum, and the oscillation energy of different periods can be analyzed, the larger the value is, the stronger the energy of the corresponding period is. Orange color is used to represent the area with a large modulus square value of the wavelet coefficient, which can be seen in the strong and weak distribution of annual inflow runoff in each time scale, as shown in Figure 6c. Among them, the energy of 26-32 time scale is strongest, the periodic distribution is obvious and almost occupying the whole study domain, the oscillation center is around 2004. The energy of 9-14 time scale is also strong, though the period is significant but locality, it mainly occurred before 2003, and the center of oscillation is around 1996. The 2-5 time scale mainly occurred in 1997-2006, but the energy is weak and locality, the center of oscillation is around 2003.
In the research time domain, the strong and weak distribution of annual inflow sediment in each time scale is shown in Figure 6d. The energy of 26-32 time scale is the strongest, and the periodic distribution is obvious, almost occupying the whole study time domain, and the oscillation center is around 2004 or 2012. The energy of 10-14 time scale is also strong, and the period is significant, and it also has the global character, the oscillation center is around 1991. The time scale of 5-8 mainly occurred before 1987, the energy performance is relatively significant but locality. The time scale of 3-4 also has periodic distribution, but the energy performance is very weak and locality.  There are four types of periodic variation rules in the evolution process of inflow runoff, namely 23-32a, 14-20a, 9-13a and 5-8a, as shown in Figure 6a. Among them, there are quasi-second oscillations between high and low runoff in time scale of 23-32a, which is relatively stable in the whole research time domain. There are quasi-third oscillations between high and low runoff in time scale of 14-20a, which mainly occurred after 1986. There are quasi-sixth oscillations in the time scale of 9-13a, which are also global. The periodic variations in the time scale below 8a are relatively chaotic, it shows that the inflow runoff changed frequently and had poor regularity in a small scale cycle. There are three types of periodic variation rules in the evolution process of inflow sediment, namely 23-32a, 9-15a and 5-8a as shown in Figure 6b. There are quasi-second oscillations between high and low runoff in time scale of 23-32a, and quasi-fifth oscillations in time scale of 9-15a. It can be seen that the periodic variation of the above two time scales are relatively stable in the whole research domain, and with a global character.

Analysis of wavelet variance and periodic characteristic
The modulus square value of the wavelet coefficient is equivalent to the wavelet energy spectrum, and the oscillation energy of different periods can be analyzed, the larger the value is, the stronger the energy of the corresponding period is. Orange color is used to represent the area with a large modulus square value of the wavelet coefficient, which can be seen in the strong and weak distribution of annual inflow runoff in each time scale, as shown in Figure 6c. Among them, the energy of 26-32a time scale is strongest, the periodic distribution is obvious and almost occupying the whole study domain, the oscillation center is around 2004. The energy of 9-14a time scale is also strong, though the period is significant but locality, it mainly occurred before 2003, and the center of oscillation is around 1996. The 2-5a time scale mainly occurred in 1997-2006, but the energy is weak and locality, the center of oscillation is around 2003.
In the research time domain, the strong and weak distribution of annual inflow sediment in each time scale is shown in Figure 6d. The energy of 26-32a time scale is the strongest, and the periodic distribution is obvious, almost occupying the whole study time domain, and the oscillation center is around 2004 or 2012. The energy of 10-14a time scale is also strong, and the period is significant, and it also has the global character, the oscillation center is around 1991. The time scale of 5-8a mainly occurred before 1987, the energy performance is relatively significant but locality. The time scale of 3-4a also has periodic distribution, but the energy performance is very weak and locality.

Analysis of Wavelet Variance and Periodic Characteristic
The wavelet variance diagram can reflect the distribution of wave energy with time scale a of the studied hydrological time series, it can be used to determine the main period in the evolution process of hydrological series [49].
There are four obvious peaks in the annual runoff time series, and they are 30a, 11a, 6a, and 4a, as shown in Figure 7a. The maximum peak value corresponds to the time scale of 30a, which indicates that the cycle oscillation of 30a is the strongest, it is the first main period of runoff variation in Hongze Lake during 1975-2015. The time scale of 11a corresponds to the second peak value, which is the second dominant period. The time series of 6a and 4a also correspond to two peak values, but the variance value is relatively small. This shows that the fluctuation of the above four periods control the variation characteristics of the annual runoff inflow in the whole study time domain.
Water 2020, 12, x FOR PEER REVIEW 11 of 18 The wavelet variance diagram can reflect the distribution of wave energy with time scale of the studied hydrological time series, it can be used to determine the main period in the evolution process of hydrological series [49].
There are four obvious peaks in the annual runoff time series, and they are 30 , 11 , 6 , and 4 , as shown in Figure 7a. The maximum peak value corresponds to the time scale of 30 , which indicates that the cycle oscillation of 30 is the strongest, it is the first main period of runoff variation in Hongze Lake during 1975-2015. The time scale of 11 corresponds to the second peak value, which is the second dominant period. The time series of 6 and 4 also correspond to two peak values, but the variance value is relatively small. This shows that the fluctuation of the above four periods control the variation characteristics of the annual runoff inflow in the whole study time domain.
There are four obvious peaks in the annual sediment time series from 1975-2015, they are 30 , 12 , 7 and 3 , as shown in Figure 7b. Among them, the maximum peak value corresponds to the time scale of 30 , which indicates that the cycle oscillation of 30 is the strongest, it is the first dominant period of sediment variation in Hongze Lake, which is the same as runoff series variation rules. The time scale of 11 and 7 have a similar peak value. The fluctuation of the above four periods controlled the variation characteristics of the annual sediment inflow in the whole study time domain.
According to the results of the wavelet variance test, drawing the real part variation process of the four main period wavelet coefficients of the runoff time series evolution, as shown in Figure 8. The average period and variation characteristic in the evolution process of annual inflow runoff under different time scales can be analyzed, the analysis results show in Table 3. The distribution characteristics of the runoff time series from 1975 to 2015 in the whole study time domain are uneven, and with significant localization characteristics. The variation process of high and low runoff under different characteristic time scales is different and closely related to the time scale. So, the inflow runoff time series exist two main periods, which are time scale of 11 and 30 . The real part variation process of four main period wavelet coefficients of sediment time series evolution is drawn, as shown in Figure 9. Moreover, the analysis results of different characteristic time scales are shown in Table 4. So, the inflow sediment time series exist two main periods, which are the time scale of 12 and 30 . From the analysis, the annual inflow runoff and sediment time series exist the same first main period, it can be predicted that the annual inflow runoff and sediment maintain in a low-value period, and will reach the valley floor value after 2015.  There are four obvious peaks in the annual sediment time series from 1975-2015, they are 30a, 12a, 7a and 3a, as shown in Figure 7b. Among them, the maximum peak value corresponds to the time scale of 30a, which indicates that the cycle oscillation of 30a is the strongest, it is the first dominant period of sediment variation in Hongze Lake, which is the same as runoff series variation rules. The time scale of 11a and 7a have a similar peak value. The fluctuation of the above four periods controlled the variation characteristics of the annual sediment inflow in the whole study time domain.
According to the results of the wavelet variance test, drawing the real part variation process of the four main period wavelet coefficients of the runoff time series evolution, as shown in Figure 8. The average period and variation characteristic in the evolution process of annual inflow runoff under different time scales can be analyzed, the analysis results show in Table 3. The distribution characteristics of the runoff time series from 1975 to 2015 in the whole study time domain are uneven, and with significant localization characteristics. The variation process of high and low runoff under different characteristic time scales is different and closely related to the time scale. So, the inflow runoff time series exist two main periods, which are time scale of 11a and 30a. The real part variation process of four main period wavelet coefficients of sediment time series evolution is drawn, as shown in Figure 9. Moreover, the analysis results of different characteristic time scales are shown in Table 4. So, the inflow sediment time series exist two main periods, which are the time scale of 12a and 30a. From the analysis, the annual inflow runoff and sediment time series exist the same first main period, it can be predicted that the annual inflow runoff and sediment maintain in a low-value period, and will reach the valley floor value after 2015.
From the analysis, the annual inflow runoff and sediment time series exist the same first main period, it can be predicted that the annual inflow runoff and sediment maintain in a low-value period, and will reach the valley floor value after 2015.          Enter into a relatively rich water period Figure 8c  4 The distribution characteristics are not obvious Figure 8d Note: is the time scale.  Maintain in a low water period and will reach the valley floor value   Enter into a relatively rich water period Figure 8c 4a The distribution characteristics are not obvious Figure 8d Note: a is the time scale. Maintain in a low sediment period and will reach the valley floor value  Enter into a relatively rich sediment period Figure 9b 7a The distribution characteristics are not obvious Figure 9c  3a The distribution characteristics are not obvious Figure 9d 3.3. Cause Analysis of Runoff and Sediment Variation

Cause Analysis of Runoff Trend Variation
The main factors for the variation of runoff in the river basin are rainfall and water consumption by human activities. Considering that the variation of inflow runoff and outflow runoff are relatively consistent, the research only analysis the influence of rainfall on the variation of inflow runoff. The variation of rainfall and runoff are relatively synchronous, and the runoff varies with the variation of rainfall, as shown in Figure 10a. In 1975In , 1983In , 1991In , 2003In , and 2007, great floods occurred in the Huaihe River basin, and the runoff inflow into the lake increased significantly. Meanwhile, in 1978Meanwhile, in , 1986Meanwhile, in , 1992Meanwhile, in , 1994Meanwhile, in , 2000Meanwhile, in , and 2001, the rainfall in the Huaihe River basin decreased, the drought led to a decrease in the runoff inflow into the lake. From 1992 to 2003, the inflow runoff decreased significantly, the M-K statistic is close to the critical value, the results of the M-K test agree with the actual variation.
Water 2020, 12, x FOR PEER REVIEW 13 of 18

Cause Analysis of Runoff Trend Variation
The main factors for the variation of runoff in the river basin are rainfall and water consumption by human activities. Considering that the variation of inflow runoff and outflow runoff are relatively consistent, the research only analysis the influence of rainfall on the variation of inflow runoff. The variation of rainfall and runoff are relatively synchronous, and the runoff varies with the variation of rainfall, as shown in Figure 10a. In 1975In , 1983In , 1991In , 2003In , and 2007, great floods occurred in the Huaihe River basin, and the runoff inflow into the lake increased significantly. Meanwhile, in 1978Meanwhile, in , 1986Meanwhile, in , 1992Meanwhile, in , 1994Meanwhile, in , 2000Meanwhile, in , and 2001, the rainfall in the Huaihe River basin decreased, the drought led to a decrease in the runoff inflow into the lake. From 1992 to 2003, the inflow runoff decreased significantly, the M-K statistic is close to the critical value, the results of the M-K test agree with the actual variation.
Through the M-K trend analysis of the multi-year rainfall time series, the M-K statistic of rainfall from 1975 to 2015 is −0.05, and the rainfall shows no obvious decreasing trend. The variation trend of inflow runoff is consistent with the rainfall, as shown in Figure 10b. Sun et al. [50] have pointed out that in the last 50 years, the annual fluctuation of precipitation in the Huaihe River basin was relatively strong, and the rainfall intensity showed an insignificant decreasing trend, which was consistent with the research results.
Before 1993, the variation trend of rainfall was consistent with the inflow runoff, and rainfall played a leading role in the runoff variation. After 1993, rainfall tends to increase while inflow runoff tends to decrease. So, 1975-1993 was selected as the base period of the whole study period. According to the regression analysis [51], the relationship between inflow runoff Y and rainfall X in the base period is: Y = 1.0197X − 577.18 (R 2 = 0.8171). According to the rainfall-runoff relationship, the average annual inflow runoff is 32.75 billion m 3 during 1993-2015, the runoff decreased by 1.09 billion m 3 compared with 1975-1993. The measured annual average inflow runoff is 27.10 billion m 3 during 1993-2015, the runoff decreases by 3.92 billion m 3 compared with 1975-1993. Thus, among the major factors influencing the decrease of runoff from 1993 to 2015, the decrease in rainfall accounted for 28%, and the increase in water use for human activities accounted for 72%.
Under investigation, the water consumption increased from 39.86 billion m 3 to 54.02 billion m 3 during 1993-2015. The M-K trend was conducted on the water consumption time series, the results showed that the statistic Z = 4.72 is far more than 95% significant level, and the water consumption increased significantly. The economy of the Huaihe River basin increased by 1.35 times and the population increased by 8 million from 1994 to 2004. With the growth of population, the rapid development of domestic water, construction and production water, agricultural irrigation water, aquaculture and other activities, the degree of exploitation and utilization of water resources has been increasing and exceeded 60% [52]. Anthropogenic activities have become the main reason for the decrease of inflow runoff after 1993.  Through the M-K trend analysis of the multi-year rainfall time series, the M-K statistic of rainfall from 1975 to 2015 is −0.05, and the rainfall shows no obvious decreasing trend. The variation trend of inflow runoff is consistent with the rainfall, as shown in Figure 10b. Sun et al. [50] have pointed out that in the last 50 years, the annual fluctuation of precipitation in the Huaihe River basin was relatively strong, and the rainfall intensity showed an insignificant decreasing trend, which was consistent with the research results.
Before 1993, the variation trend of rainfall was consistent with the inflow runoff, and rainfall played a leading role in the runoff variation. After 1993, rainfall tends to increase while inflow runoff tends to decrease. So, 1975So, -1993 was selected as the base period of the whole study period. According to the regression analysis [51], the relationship between inflow runoff Y and rainfall X in the base period is: Y = 1.0197X − 577.18 (R 2 = 0.8171). According to the rainfall-runoff relationship, the average annual inflow runoff is 32.75 billion m 3 during 1993-2015, the runoff decreased by 1.09 billion m 3 compared with 1975-1993. The measured annual average inflow runoff is 27.10 billion m 3 during 1993-2015, the runoff decreases by 3.92 billion m 3 compared with 1975-1993. Thus, among the major factors influencing the decrease of runoff from 1993 to 2015, the decrease in rainfall accounted for 28%, and the increase in water use for human activities accounted for 72%.
Under investigation, the water consumption increased from 39.86 billion m 3 to 54.02 billion m 3 during 1993-2015. The M-K trend was conducted on the water consumption time series, the results showed that the statistic Z = 4.72 is far more than 95% significant level, and the water consumption increased significantly. The economy of the Huaihe River basin increased by 1.35 times and the population increased by 8 million from 1994 to 2004. With the growth of population, the rapid development of domestic water, construction and production water, agricultural irrigation water, aquaculture and other activities, the degree of exploitation and utilization of water resources has been increasing and exceeded 60% [52]. Anthropogenic activities have become the main reason for the decrease of inflow runoff after 1993.

Trend Variation of Sediment and Cause Analysis of Change Point
Sediment transport is mainly related to climate variability and anthropogenic activities in the basin. In the trend analysis, the variation trends of sediment flowing into and out of Hongze Lake are related to runoff. The M-K test statistics of rainfall and runoff in the study period did not exceed the significance level, while the sediment exceeded the significance level and showed a more obvious decreasing trend. It meant rainfall and runoff had no obvious effects on the decreasing trend of sediment, anthropogenic activities may had a great influence on sediment transport.
Since 1949, large-scale reservoir construction has been carried out in the Huaihe River basin to develop hydropower resources rationally and reduce flood disasters. At present, 18 large reservoirs have been built successively in the Huaihe River basin, with the maximum capacity of flood storage reaching 8.89 billion m 3 . Because the period from 1975-2015 is not the prosperous period of reservoir construction in the Huai River basin, therefore, the research period is divided by the completion time of large and medium-sized reservoirs and storage capacity changing time in the Huaihe River basin. For example, the Banqiao reservoir began to be rebuilt in 1978, and the Huaihe River basin was flooded in 1991. Therefore, the period is divided into 1975-1978, 1979-1991, and 1992-2015. The reservoir capacity increased from 14.75 billion m 3 to 15.854 billion m 3 , the average inflow sediment decreased from 7.98 million t to 5.22 million t, and the average outflow sediment decreased from 5.24 million t to 2.27 million t. As shown in Figure 11, the inflow sediment and outflow sediment decreased with the increase of the total reservoir capacity, and the silt trapping of the reservoir has a significant effect on the inflow and outflow sediment.
For example, the Banqiao reservoir began to be rebuilt in 1978, and the Huaihe River basin was flooded in 1991. Therefore, the period is divided into 1975-1978, 1979-1991, and 1992-2015. The reservoir capacity increased from 14.75 billion m 3 to 15.854 billion m 3 , the average inflow sediment decreased from 7.98 million t to 5.22 million t, and the average outflow sediment decreased from 5.24 million t to 2.27 million t. As shown in Figure 11, the inflow sediment and outflow sediment decreased with the increase of the total reservoir capacity, and the silt trapping of the reservoir has a significant effect on the inflow and outflow sediment. Figure 11. The relationship between sediment flowing into and out of Hongze Lake and the total volume of the upstream reservoir.
According to the above study, the change point of sediment flowing into and out of Hongze Lake occurred in 1991. although there was no significant change in rainfall from 1975-2015, but heavy rainfalls affect the variation of runoff and sediment in a short time. For example, in 1991, the average rainfall of the Huaihe River was 1082 mm, 20% more than the average year. The rainfall is concentrated from June to August, then the runoff and sediment flowing into Hongze Lake increased significantly, and the largest value occurred in July. Meanwhile, after 1991, the rainfall decreased and the M-K test statistic of rainfall is −0.37 < 0 during 1991-2015, which showed that the rainfall decreased after 1991. Rainfall reinforced the decreasing trend of runoff and sediment.
In terms of anthropogenic activities, LUCC can result in alterations in rainfall-sediment yield relationship in the Huaihe River basin [53]. According to statistics, anthropogenic activities strongly disturbed the surface of the Huaihe River basin, causing the changes of land use in the region. Before the 1980s, the Huaihe River basin was dominated by cultivated land and forest. With the increase of the regional population, more natural wetlands and forests were transformed into construction land and farmland, water and soil erosion in the basin was becoming more and more serious, the ecological environment of the river basin was seriously damaged. In the 1980s, the government started soil and water conservation work in the basin, focusing on reducing sloping farmland, increasing vegetation, and closing mountains for forest and grass, thus controlling soil erosion to a Figure 11. The relationship between sediment flowing into and out of Hongze Lake and the total volume of the upstream reservoir.
According to the above study, the change point of sediment flowing into and out of Hongze Lake occurred in 1991. although there was no significant change in rainfall from 1975-2015, but heavy rainfalls affect the variation of runoff and sediment in a short time. For example, in 1991, the average rainfall of the Huaihe River was 1082 mm, 20% more than the average year. The rainfall is concentrated from June to August, then the runoff and sediment flowing into Hongze Lake increased significantly, and the largest value occurred in July. Meanwhile, after 1991, the rainfall decreased and the M-K test statistic of rainfall is −0.37 < 0 during 1991-2015, which showed that the rainfall decreased after 1991. Rainfall reinforced the decreasing trend of runoff and sediment.
In terms of anthropogenic activities, LUCC can result in alterations in rainfall-sediment yield relationship in the Huaihe River basin [53]. According to statistics, anthropogenic activities strongly disturbed the surface of the Huaihe River basin, causing the changes of land use in the region. Before the 1980s, the Huaihe River basin was dominated by cultivated land and forest. With the increase of the regional population, more natural wetlands and forests were transformed into construction land and farmland, water and soil erosion in the basin was becoming more and more serious, the ecological environment of the river basin was seriously damaged. In the 1980s, the government started soil and water conservation work in the basin, focusing on reducing sloping farmland, increasing vegetation, and closing mountains for forest and grass, thus controlling soil erosion to a certain extent. These anthropogenic activities have changed the land use and cover of the study area, thus causing the decrease of runoff and sediment to some extent. Besides, hydraulic engineering projects played an important role to decrease the sediment, due to the series flood disaster in the Huaihe River basin, China has designed 19 projects to control the Huaihe River, by 1998 (when the ShiMantan reservoir was completed), the upstream reservoir of Hongze Lake had increased by 0.77 billion m 3 , part of the sediment was intercepted by the reservoir. At the same time, the implementation of the projects to control the Huaihe River is beneficial to soil and water conservation, since 1991, the upstream of Hongze Lake has been treated with 2184 km 2 of soil erosion.
So, the climate variability and anthropogenic activities have a significant impact on the sediment flowing into and out of Hongze Lake, while changes in heavy rainfall patterns, the construction of hydraulic engineering projects, and LUC are the main factors.

Conclusions
This paper based on the M-K test and Pettitt test, a combinatorial mutation test method with high precision was proposed to analyze the trend of runoff and sediment time series of Hongze Lake from 1975 to 2015. The results show that the runoff flowing into and out of Honze Lake were within the confidence interval of significance level a = 0.05, but didn't exceed the significance level, and there were no obvious trend variation and change points. While the sediment flowing into and out of Honze Lake showed a decreasing trend, and the change point is 1991.
The Morlet wavelet analysis method was used to analyze the periodicity of annual runoff and sediment time series flowing into Hongze Lake from 1975 to 2015. There existed four scales of the periodic variation in the time series, and in different time scales. The rich and dry variation trends of inflow runoff and sediment time series were different and closely related to the time scales. The annual inflow runoff and sediment time series had a similar variation period, with the same first main period and the characteristic scale was 30a, which played a leading role in the rich and dry variation trend of annual inflow runoff and sediment in Hongze Lake. Under this characteristic scale, it can be predicted that the annual inflow runoff and sediment would enter into the low-value period respectively in a short time after 2015, and both would be in the valley floor stage.
Finally, analyzing the variation trend and mutation cause of runoff and sediment in Hongze Lake from the climate variability and anthropogenic activities. Changes in heavy rainfalls pattern, the construction of hydraulic engineering projects and LUCC in the Huaihe River basin, leading to the decreasing trend of sediment flowing into and out of Hongze Lake. The results show that the combinatorial mutation test method has higher accuracy, the time series of annual runoff and sediment flowing into and out of Hongze Lake had the characteristics of multi-time scale variation, the variation trends of runoff were mainly affected by rainfall, and the sediment variation trends and change point were mainly affected by changes in heavy rainfall patterns, the construction of hydraulic engineering projects, and LUCC. The analysis methods and results provide the theoretical basis for flood prevention and control in the Huaihe River basin and further study on runoff and sediment in Hongze Lake.