Next Article in Journal
Agricultural Diffuse Pollution and the EU Water Framework Directive: Problems and Progress in Governance
Previous Article in Journal
Bank Protection Structures along the Brahmaputra-Jamuna River, a Study of Flow Slides
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Runoff-Sediment Cointegration and Uncertainty Relations at Different Temporal Scales in the Middle Reaches of the Yellow River, China

State Key Laboratory of Hydraulic Engineering Simulation and Safety, Department of Hydraulic Engineering, School of Civil Engineering, Tianjin University, Tianjin 300350, China
*
Author to whom correspondence should be addressed.
Water 2020, 12(9), 2589; https://doi.org/10.3390/w12092589
Submission received: 28 July 2020 / Revised: 11 September 2020 / Accepted: 11 September 2020 / Published: 16 September 2020
(This article belongs to the Section Water Erosion and Sediment Transport)

Abstract

:
To understand the intricate runoff-sediment relationship in the middle Yellow River basin (MYRB), the Toudaoguai, Longmen, Tongguan and Huayuankou sites in the MYRB were selected to analyze the deterministic equilibrium and uncertainty relations of runoff-sediment based on 55-year hydrological data at multi-time scales. The Johansen test and wavelet neural network were used to verify the cointegration relationship among hydrological series. Runoff-sediment uncertain statistical relations and dynamics in the MYRB were also analyzed based on rating curves and hysteresis loops. The results showed that the logarithmic sequences of sediment load (SL), runoff and suspended sediment concentration (SSC) conformed to a linear cointegration relationship at the Toudaoguai station or in spring, winter or under small flow at other stations, but a nonlinear cointegration relationship was observed in other cases at other stations. Regarding runoff-sediment uncertain relationships, the rating curves, and hysteresis loops differed in stations (Toudaoguai and the other stations), as well as discharge (threshold: 1000 m3/s), season (ice-flood and rainy season) and saturation of flow at flood and monthly scales. At the annual scale, phased and unsynchronized characteristics of runoff and sediment load were evident with a decreasing trend. This study on the runoff-sediment relationship can rationally provide a theoretical basis for the management and development of the Yellow River and other similar rivers with sufficient sediment, especially for areas with serious soil erosion.

1. Introduction

The Yellow River in China is famous for its heavy sediment load (SL), as well as its different sources of runoff and sediment and their intricate spatio-temporal changes, especially in the middle Yellow River basin (MYRB), which is endowed with unconsolidated and erodible loess, contributing 90% of the sediment load [1]. The behavior of sediment in the MYRB is an obvious nonlinear dynamic system, which is controlled by the complex interactions of time-varying sediment availability, rainfall intensity and distribution [2,3], as well as the morphology (slope, cross section shape) [4] and composition of the riverbed. Therefore, this implies an inharmonious and imbalanced relationship between runoff and sediment and complicated management issues of the Yellow River, such as the evolution of the channels [5,6], ecology and landscape [7,8] and the lifespan of infrastructure [9].
To date, the research on the Yellow River has mainly focused on particular areas with heavy sediment load for the regulation of Yellow River, including reaches of Ningxia–Inner Mongolia [10,11] and the Hekouzhen–Longmen basin with abundant-coarse sediment [12,13], as well as main tributaries [6]. Due to the difference of grain size composition and geomorphologic agent, the regional characteristics of the Yellow River basin are obviously different. At present, there are few studies on runoff-sediment relationships across different characteristic regions. Due to the different sources of runoff and sediment, particle size composition and the effect of artificial regulation in the Yellow River basin, there is an obvious uncertain relationship between runoff and sediment, which has attracted the attention of numerous scholars [14,15,16,17]. In the study of the visualized runoff-sediment uncertain relationship during the process of sediment transport, a range of empirical models, such as sediment rating curves and hysteresis loops, are used to analyze and quantify suspended sediment concentration (SSC) and sediment track in rivers. Hysteresis patterns can provide useful insights into the sediment sources, sediment delivery and feedback mechanism [3,7,18]. Accurate estimates of sediment concentration by sediment rating curves are needed for erosion control structures, river morphological computations, management of ecosystems [2]. The regression parameters of the rating curve represent different physical meanings and explain the regional difference in sediment transport characteristics. For example, James et al. studied the correlations among SSC rating parameters, river basin morphology and climate [19]. The understanding of the uncertain relationship between runoff and sediment is of great significance to the construction of the water and sediment control system in the MYRB, the implementation of water resources management and configuration according to local conditions and the design of major water conservancy projects.
In previous research, the hydrological time series was generally assumed to be stable. However, the hydrological sequence generally includes long-term memory and is non-stationary [20]. Establishing a regression model based on non-stationary sequences for the runoff-sediment relationship is prone to the “pseudo-regression” problem [7]. Equilibrium relationships can be described as follows: the linear combination of non-stationary sequences is a stationary sequence [21]. It provides a scientific statistical method for non-stationary time series and interprets the long-term stable quantitative relationship among variables. The statistical relationship among variables with cointegration relationship no longer exists pseudo regression problem. Meanwhile, the error correction model (ECM), based on the cointegration relationship, can predict the hydrological variables effectively. Up to now, cointegration theory has been widely used in the field of econometrics, but it is rarely applied to watershed hydrology [22,23].
Focused on the above issues, this paper provides an in-depth study on the runoff-sediment stable equilibrium relation and uncertain relationship of Toudaoguai, Longmen, Tongguan and Huayuankou stations in the MYRB with different regional characteristics, followed by attempts to explore the hidden cointegration relationship, runoff-sediment dynamics, such as sediment origin and distribution, pollutant transport and the channel evolution and their statistical relation. This research rationally provides a theoretical basis for the management and development of the Yellow River and other similar rivers with sufficient sediment, especially for areas with serious soil erosion.

2. Study Area and Materials

The MYRB between the Toudaoguai and Huayuankou stations, located within 32°–42° N and 104°–113° E (Figure 1), significantly contributes to the total sediment of the Yellow River, due to the Loess Plateau. The area is characterized by loose soil, broken terrain, low vegetation coverage and frequent rainstorms [6,24]. Toudaoguai, Longmen, Tongguan and Huayuankou stations from upstream to downstream are distributed on the main stream of the MYRB, which locates in areas with different regional characteristics. The upstream reaches of the Toudaoguai are notable for their low sediment load and high water discharge. The reaches between the Toudaoguai and Longmen are mainly covered by coarse-grained and overburden sediment. From the Longmen to Sanmenxia, the river is characterized by high sediment load and fine sediment, and the reaches between the Sanmenxia and Huayuankou, including the Yiluo River and Qing River basins, have low sediment loads [25]. The length of this segment of the river is 1206.4 km (749.6 miles), which accounts for 22.1% of the total length of the river. The MYRB is 344,000 square kilometers (132,819 square miles), accounting for 45.7% of the total river basin [26]. Meanwhile, it accounts for 89% of the annual average sediment transport, and 90% of which comes from the rainy seasons. Tributaries with abundant sediment in the MYRB are plentiful, which carry an average of 56% of the river’s annual sediment load. In addition, its annual mean temperature ranges from about 6 to 11 °C [26], and the spatial and temporal distribution of precipitation within the basin is uneven. Annual average precipitation ranges from 320 mm (12.6 inches) in the north to 836 mm (32.9 inches) in the south, and the potential evapotranspiration ranges from 810 (31.8 inches) to 1260 mm (49.6 inches). To meet the needs of flood control, power generation, navigation and irrigation, a series of reservoirs including Longyangxia (1986), Liujiaxia (1968), Qingtongxia (1960), Sanmenxia (1960) and Xiaolangdi (1997) reservoirs from upstream to downstream were put into operation, which greatly change the runoff-sediment conditions.
The data used in this paper includes water discharge and SSC series of total 285 floods from 1960 to 2005 and corresponding daily series (including sediment transport rate) from 1960 to 2015 of the four stations. Based on the existing data, the homologous monthly and annual data of discharge, runoff, SSC and sediment load (suspended load and bed load) could be calculated. The data series were mostly collected from the Yellow River Water Conservancy Commission (YRCC) of the Ministry of Water Resources and the Yellow River Sediment Bulletin, which were all measured based on the Chinese national standard criterion.

3. Methods

3.1. Cointegration Theory

Non-stationary series are universal, which limits the use of the classical regression equation. Cointegration analysis is a statistical property that describes the long-term equilibrium relationship among time series and first proposed by Engle and Granger in 1978 to avoid the pseudo-regression problem among non-stationary sequences [21]. Its exact definition can be expressed as follows: if a time series becomes a stationary sequence after d differences, then the original sequence is called a d-order single integer sequence, which is recorded as I(d). If multiple time series X 1 , , X n are all of the same order I(d), and there is a vector ( α 1 , , α n ) that makes the weighted combination α 1 X 1 + + α n X n be a stationary sequence I(0), then this set of time series is cointegrated, which excludes the random trend caused by the unit root of the variable. The exploration of the cointegration relationship can reveal the long-term quantitative relationship among variables and avoid the pseudo-regression problem when it comes to statistical relationships among variables in a cointegration relationship. Meanwhile, the ECM based on cointegration relationship distinguishes between long-term equilibrium relations and short-term volatility relations, which can automatically correct the early errors and predict variables. Common cointegration test methods include the Engle-Granger (EG) two-step test (two variables) [21] and the Johansen test (more than two variables) [27]. The EG two-step test is a stationarity test based on regression residuals. The Johansen test was proposed by Johansen and Juselius to test the cointegration relationships of multi-variables with maximum likelihood estimation based on the VAR model. This study involves three variables, and there may be more than one cointegration relationship. Therefore, the Johansen test is used in this work to verify the cointegration relationship and determine the cointegration equation. Its specific principle is as follows [28]:
Presume that the vector y t is a p-order autoregressive model (VAR):
y t = A 1 y t 1 + + A p y t p + B x t + ε t
Then, convert this VAR into the following form of vector error correction model (VECM):
Δ y t = Π y t 1 + i = 1 p 1 Γ i Δ y t i + B x t + ε t
where Π = i = 1 p A i I , Γ i =   j = i + 1 p A j . If 0 < rank (Π) < n, there are two matrixs α and β whose order are all n × r . makes Π = α β . and Π = α β is I (0). When the rank of the matrix is r, then the matrix has r non-zero eigenroots, which is arranged by values as λ 1 , λ 2 , , λ r . The significance test for the number of eigenroots can be performed by the trace statistic λ t r a c e and the maximum eigenroot statistic λ m a x .
λ t r a c e = T i = r + 1 m l o g ( 1 λ i )
λ m a x = T l o g ( 1 λ r + 1 ) .
where λ i represents the estimate of eigenroots; T represents the sample size.
The relationship can be called “linear cointegration” if the time series satisfy a long-term linear equilibrium relationship [29,30]. According to the cointegration theory, an ECM can be established to describe the long-term static performance and short-term dynamic characteristics among the sequences. A first-order autoregressive distribution lag model ADL (1, 1) with only two variables was considered. The expression of the error correction model is as follows:
Δ y t = β 0 + β 1 Δ x t + α E c m ( t 1 ) + ε t
where Δ y t represents the short-term fluctuations of the explained variable; Δ x t represents the short-term fluctuations of the explanatory variables; E c m ( t 1 ) is the error correction term at the previous moment, reflecting the degree to which the previous term of the dependent variable deviates from its long-term equilibrium state in short-term fluctuations; α is the correction coefficient, representing the adjustment speed; ε t represents residuals; β 0 and β 1 are the regression coefficients.
Hurst discovered that hydrological time series has long-term memory characteristics in 1951 [31], showing that there is a nonlinear relationship between hydrological time series. In addition, many hydrological time series have fractal dimensions, so the long-term equilibrium relationship is nonlinear. For the vector time series X t = ( x 1 t , x 2 t , , x p t ) T , the sequence of { X t } components is called nonlinear cointegration if (1) x i t ( i = 1 , , p ) is a long-term memory mean (LMM) sequence, and (2) there is a nonlinear function f such that Z t = f ( x 1 t , x 2 t , , x p t ) is a short-term memory mean (SMM) sequence. The function f() is called a nonlinear cointegration function. The wavelet neural network (WNN), a typical nonlinear function, has been usually used to compute nonlinear cointegration function [32]. The WNN model can be expressed according to the following formula:
y t = n = 1 N ω n h ( ( m = 1 M v n m x m t b n ) / a n ) .
where M is the number of nodes of the input layer; N is the number of nodes of hidden layer; x m t is the input vector; the base function h ( ( x b n ) / a n ) (n = 1, 2,…N) is the wavelet function; a n is the expansion parameter and b n is the translation parameter; v n m is the connection weight from input layer to hidden layer; ω n is the connection weight from hidden layer to the output. y t is the output value.

3.2. Methodology for Runoff-Sediment Relationship

A sediment rating curve, reflecting the runoff-sediment relationship, describes the relation between discharge (Q) and SSC for a certain location of the river through a mathematical model [33,34]. Sediment rating curves based on the existence of cointegration relationship can avoid spurious regression in the study of runoff-sediment relationship and truly reflect the actual relationship among hydrological variables. The power function is usually applicable for runoff-sediment relation when it is unsaturated flow. Once the flow is saturated, the sediment rating curves (SRC) between maximum SSC and Q emerges with linear form (Equation (7)) [2].
S S C = c ( Q Q t ) + d
where c and d are fitting coefficients, and Q t represents the threshold discharge between unsaturated and saturated flows.
The hysteresis loop with clockwise, anticlockwise or figure-eight shape is also introduced to interpret runoff-sediment relationship. The empirically-based approach which compares the simultaneity or interval between the SSC peak and discharge maximum is widely applied [5,9].

4. Results and Discussion

4.1. Runoff-Sediment Cointegration Equilibrium Relation

To reduce the absolute value of the data, eliminate the heteroscedasticity and facilitate calculation, logarithmic data were used for the cointegration calculation. The annual data of SL, runoff and SSC from 1960 to 2000, as an example, were first used for the study of the cointegration relationship among variables. Evidently, there was a good positive correlation between discharge and SSC at the Toudaoguai station at an annual scale, and a facial linear relationship existed among the logarithmic hydrological sequences of SL, runoff and SSC. This is in agreement with Fan’s conclusion that there was a synchronous decrease of water and sediment at the Toudaoguai station, and SSC-Q fitted the power function when Q < 2000 m3/s [5,35]. However, at the other stations, obvious changes of SSC did not correspond with the discharge due to a sharp increase of sediment in the MYRB [36,37], and no obvious statistical model existed for the scattered data.

4.1.1. Linear Cointegration Relationship at Toudaoguai Station

In order to avoid problems such as pseudo-regression, it is necessary to verify the stationarity of the sequence before testing the cointegration relationship among the variables. The ADF (Augmented Dickey-Fuller) test was used to test the unit root of each time series at the Toudaoguai station. The test results are shown in Table 1.
It can be seen from the unit root test results in Table 1, although the ADF statistics of the sequences SL, runoff and SSC were all less than the 5% threshold, the P values of the trend items were also all less than 0.05, which illustrated that the above three series at the Toudaoguai station were all stationary with a decreasing tendency. As Zhao et al. and Wang et al. explained that the reduced rainfall, reservoir regulation, increased water use for industry and agriculture and the effect of water and soil conservation could induce the downward trend [24,38]. After a first-order difference, it was determined that the three sequences at the “none case” were all stationary according to the principle of minimum AIC (Akaike Information Criterion) and SC (Schwarz Criterion) values. This implied that the time series were all subordinated to I(1). Thus, it was necessary to conduct the Johansen test on a non-stationary time series SL, runoff and SSC, to investigate whether there was a long-term stable equilibrium relationship.
Before conducting cointegration tests on non-stationary hydrological sequences, a vector autoregressive model needed to be performed objectively to determine the maximum lag period of the relevant variables. The values of each measurement criterion with different lag orders are given in Table 2, indicating that the optimal lag order was an order of 1.
After determining the optimal lag order of the related variables, the Johansen cointegration test was used to verify the existence of long-term stable equilibrium relationships. The results of the cointegration test are shown in Table 3.
The P values obtained by the trace test and maximum eigenvalue test of three hypotheses of the number of cointegration relations were all less than 0.05. Thus, all three cases rejected the null hypothesis, which indicated that there were at least three cointegration equations among the variables. The following are the estimated cointegration relationships among the three variables.
{ S L = 0.533   r u n o f f + 1.429   S S C 1.865 S L = 1.916   S S C 0.807 S L = 2.125   r u n o f f 5.036 r u n o f f = 0.914   S S C + 1.985
The Johansen cointegration test results indicated three sequences, and all pairs of variables had long-term linear cointegration relationships; i.e., at the Toudaoguai station, runoff and SSC had long-term stable quantitative influence on SL. Furthermore, there was also a stable association between runoff and SSC. This was similar to the research results that there was a good positive correlation between discharge and SSC at the Toudaoguai station. Therefore, the linear cointegration relationship among the variables of the Toudaoguai station is also applicable at other time scales. The determination of the quantitative relationship of the three variables is helpful to the calculation and prediction of hydrological series at the Toudaoguai station. Meanwhile, the existence of the cointegration relationship among the variables also indicated that there was no pseudo-regression problem concerning the study of the uncertainty relationships between them. In addition, the existence of the linear cointegration relationship paves the way for the establishment of the error correction model for the more accurate prediction of variables. Annual SL, runoff and SSC data from 1960 to 2000 were used to perform the Johansen cointegration test and calibrate the error correction model (ECM). Annual data from 2001 to 2005 were used to verify the parameters of ECM. The estimation formulas of the error correction models are shown in equation 5. Moosa et al. pointed out that the coefficient of the difference term indicates the degree and direction of short-term fluctuations on the short-term changes of variables, and the coefficient of the correction term indicates the strength of adjusting the short-term fluctuations from the equilibrium state variable until the equilibrium state is reached [39].
{ E C M 1 : Δ S L = 1.101   Δ S S C + 0.917   Δ r u n o f f 0.838   E c m ( 1 ) + 0.0004 E C M 2 : Δ S L = 1.813   Δ S S C 0.876   E c m ( 1 ) + 0.0022 E C M 3 : Δ S L = 2.116   Δ r u n o f f 0.584   E c m ( 1 ) 0.0015 E C M 4 : Δ r u n o f f = 0.780   Δ S S C 0.595   E c m ( 1 ) + 0.0020
During the calibration period (1960–2000), the Nash coefficients between the calculated values of ECM and measured values were 0.96, 0.93, 0.88 and 0.81, respectively. The corresponding Nash coefficients between the calculated values of the ordinary least squares (OLS) and measured values were 0.95, 0.91, 0.84 and 0.75, respectively. The calibration results indicated that the calculation accuracy of ECM was significantly improved compared with OLS. Christoffersen et al. provided an explanation as to why the addition of the error correction term leads to the enhancement of forecasting [40]. Likewise, Jacobson et al. also showed that the performance of the ECM was better if long-run theory restrictions were imposed on the cointegrating relation [41]. Therefore, the use of the ECM is recommended to predict the hydrological variables, such as SL or runoff at the Toudaoguai station. More accurate prediction of runoff and sediment is of great significance to the flood control and disaster reduction, rational allocation of water resources, the solution of soil erosion and river deposition in the river basin. Meanwhile, the cointegration relationship among the three variables was the strongest with the highest Nash coefficient, and the relatively poor cointegration relationship between runoff and SSC was obvious. This was mainly dominated by the uncertain relationship between discharge and SSC, due to the incongruous relationship with reservoir regulation and the availability of sediment [35,42]. During the validation period (2001–2005), the comparative calculation results between the four ECMs and the OLS equations are shown in Table 4 and Table 5, the prediction of the four ECMs was better than those of the OLS except for a few years. Meanwhile, the error of ECM based on three variables was smaller than that of the other models on average in the prediction of SL, which is consistent with the results of the calibration period. Ultimately, the prediction errors were less than 20% except for ECM 2 in 2002, ECM1 and ECM3 in 2004, and the minimum prediction error was 0.36%, which verified the good prediction ability of ECM.

4.1.2. Nonlinear Cointegration Relationship at the Other Stations

Prior to examining nonlinear cointegration, the R/S statistic method was used to determine the long-term memory characteristics of the annual average SL, runoff and SSC of the Longmen, Tongguan and Huayuankou stations from 1987 to 2005. Hurst exponents are given in Table 6.
It can be seen from Table 6 that the Hurst index values were more than 0.5 for all time series at three stations; thus, they were all LMM time series. The empirical formula of the fractal dimension D is D = 2 − H; the fractal dimensions of the original sequences are shown in Table 7.
Table 7 clearly shows that the fractal dimensions of the three hydrological sequences at each site were different; thus, which meant none of the sequences conformed to a linear cointegration relationship [20]. Therefore, the wavelet neural network (WNN), served as a fitting method, was used to verify the existence of nonlinear cointegration among variables [22,23]. With the WNN method, WNNs with a single hidden layer were established, and the structure of WNN was 3–8–1, 2–8–1, 2–8–1, 2–8–1 for the input variables SL–runoff–SSC, SL–runoff, SL–SSC, and runoff–SSC, respectively. A short-term memory mean (SMM) zero time series served as the guidance output. The mallet wavelet was adopted for the wavelet base function. The training results showed that the function f(⋅) existed in all of the the above training combinations, which indicated that y t = f ( x 1 t , x 2 t , , x n t ) was a SMM time series. In other words, all of the above combinations showed nonlinear cointegration. The relevant parameters of the nonlinear cointegration function of different training combinations are shown in Table 8, Table 9, Table 10 and Table 11. Since the statistical relationships between hydrological series at the Longmen, Tongguan and Huayuankou stations were not significant and their nonlinear relationships were strong, the ECM based on nonlinear cointegration relationship was not studied in depth in this paper.

4.2. Runoff-Sediment Uncertainty Relation

4.2.1. Runoff-Sediment Relation at the Within-Flood Event Scale

As shown in Figure 2a, the relation between discharge and SSC at the Toudaoguai station could be described by a cubic polynomial function at the within-event scale. This was explained from a macro perspective by Fan et al. [39] that the decreasing velocity (caused by overbank discharge) promotes the formation of sandbars, which narrows the channel width and accelerates the further scouring of the channel bed. Dissimilarly, at the other three stations, wind-dominated geomorphologic agents (geomorphologic agents here refer to external agents, which are mainly generated by solar energy, gravitational energy and biological activitiesgravity, including hydraulic force, wind force, ice force, wave force etc.) in the spring and winter transport coarse sediment to slopes, gullies and river courses [43]. However, coarse particles cannot enter the slurry, because of the lack of sufficient fine sediment and low energy flow [44], resulting in a lower SSC with respect to Q < 1000 m3/s and emerging rating curves in the form of a power function. Simultaneously, the surface material composition of the MYRB has significant zonal characteristics: the composition is in the order of eolian sediment, sand loess covered by discontinuous aeolian sand, loess and clay loess from northwest to southeast [43,45], where the granularity composition gradually becomes finer. The high slurry viscosity of low energy flow for dominant fine sediment may form a laminar flow [46], which would lead to the synchronization of runoff and sediment. Therefore, the situation is similar to that of the Toudaoguai station at this time, and SL, runoff and SSC satisfy the linear cointegration relationship at the flood scale, and the same is true for other scales. When Q > 1000 m3/s, energy flow with abundant fine sediment will continue to carry the coarse sediment carried by the wind in winter and spring downward with a high amplitude of concentration, which promotes the formation of hyperconcentrated flow. The linear cointegration relation no longer applied. The SRCs between the maximum SSC values for a given discharge displayed a linear form once the flow was saturated and stable, implying a global maximum with an optimal particle size composition ratio [47] (Figure 2b–d). The utilization of optimal particle size ratio and discharge threshold (threshold: 1000 m3/s) is beneficial to the dredging of reservoirs. It was deduced that the sediment composition of the Longmen station was coarser than that of the Tongguan station, because hydraulic erosion plays a dominant role in the southeast, where the surface composition was dominated by fine particles from loess [48]. Therefore, the maximum SSC was larger than that of the Tongguan station, due to the appropriate particle size composition. In addition, the SSC values in the third stage were all larger than those in the previous periods at the three stations due to the large-scale utilization of clean water resources [49,50] and enhanced scouring capacity for the regulation in runoff of Longyangxia (relatively weak regulation of sediment because of location) [51]. Meanwhile, the optimized particle size ratio for the enhanced coarse sediment-carrying capacity and abundant fine sediment produced by “trapping coarse and discharging fines” of reservoirs and water-and-soil conservation measures also facilitated the increase in the SSC.
Sediment rating curves can not explain the internal runoff-sediment relationship from a microscopic point of view; thus, the hysteresis loop was introduced to determine the hysteresis behavior according to an empirically-based approach, which can provide significant information on suspended sediment sources. Clockwise loops, which occur in late autumn and winter frequently [7,8], can be attributed to the available and deposited sediment source nearby, indicating sufficient available sediment and a quick increase of SSC at the beginning [52]. Meanwhile, overbank floods may be another reason for the reduction of SSC with higher discharge. Sediment from upstream, tributaries and hillslopes (gravitational erosion with gradual penetration) is a major reason for anti-clockwise hysteresis, which may suggest the delay of sediment entrainment within a runoff event [7]; thus, slumps and slides usually occur at the falling stage of a flood, which usually occurs in early spring, summer and early autumn [7,8]. Simultaneously, the intensive sediment-carrying capability of hyperconcentrated flow results in sediment with minimal deposition, even with a lower energy flow in the falling stage of a flood, which also gives rise to an anticlockwise hysteresis loop [53]. Figure-eight shapes indicate more complex situations that mainly appear in early spring and summer [52]. To study the proportions of the three hysteresis types and the proportional changes of the same types at the four stations, 47, 62, 75 and 101 typical floods (an obvious flood process in which only one SSC peak occurs during a flow fluctuation) that occurred at the four stations were analyzed. It was observed that the clockwise, counter-clockwise, and figure-eight types accounted for 23.4%, 19.15% and 57.45% of floods at the Toudaoguai station, respectively. Correspondingly, the proportions of the three types at the Longmen, Tongguan and Huayuankou stations were 4.84%, 70.97% and 24.19%; 22.67%, 37.33% and 40%; and 19.80%, 45.54% and 34.65%, respectively. It was evident that counter-clockwise loops occurred at higher rates at the Longmen, Tongguan and Huayuankou stations, which was closely related to the massive sediment supply from abundant tributaries in the MYRB [24], the formation of the hyperconcentrated flow mentioned above and the sediment-delayed effect of soil and water conservation. The decrease of the proportion of anti-clockwise loops along the way was also related to the decrease of the formation frequency of hyperconcentrated flow with a predominance of fine particles and the decreased effect of soil and water conservation, which was verified by Xu et al. that the high frequency of the hyperconcentrated flow on the Loess Plateau occurred from northeast to southwest in space [54] and the sediment reduction of silt dam in the MYRB reached its peak in the 1970s, and has been declining ever since [55].

4.2.2. Runoff-Sediment Relation at the Monthly-Seasonal Scale

Apparently, sediment rating curves varied seasonally at the four stations for the entire period. In the ice-flood season, the runoff, which was of a small amount, mainly originated from melting ice, and the sediment condition was similar to that when Q < 1000 m3/s. Therefore, changes in SSC and Q were basically simultaneous, and the sediment rating curves derived from linear regression fitted the data well (Figure 3a,c,e,g). Specifically, the rating curve at the Toudaoguai station could be described by a cubic polynomial function with a high correlation coefficient (0.89) in rainy seasons (Figure 3b). Nevertheless, the correlations between SSC and Q were very low from the Longmen to Huayuankou stations due to the complex supply of sediment and the formation of hyperconcentrated flow with different combinations of coarse and fine sediment [43], shaping multiple SSC values for a given discharge (Figure 3d,f,h). Furthermore, the fluctuation frequency slowed down in the order of the Longmen, Tongguan and Huayuankou stations, which was determined by soil and water conservation measures, the operation of Xiaolangdi reservoir [56] and the relatively stable sediment composition [43,44] (dominated by fine particles) along the way.
At the monthly scale, the research period of the Toudaoguai station was divided into three stages with the same sudden change point of discharge and SSC caused by the operation of the Liujiaxia reservoir (1968) and the Longyangxia reservoir (1986) [5,35]. It was observed that hysteresis at the Toudaoguai station exhibited figure-eight patterns (Figure 4). Nevertheless, hysteresis loops with clockwise trends were found for the other stations (Figure 5). Analyses of the results showed that the sediment from the MYRB was much larger than that from upstream; in other words, the nearby suspended and deposited sediments were the main sediment source for the MYRB, which explained why the relations between Q and SSC showed hysteresis with clockwise trends at the other stations. Furthermore, the shapes of the C–Q hysteresis loops at the four stations evolved from linear loops at the Toudaoguai station to flattened loops with a sharp decrease of SSC and high discharge in September and October at the Longmen, Tongguan and Huayuankou stations, because the formation frequency of hyperconcentrated flow was greatly reduced due to the depletion of coarse sediment at the end of the rainy reason. Simultaneously, there were similarities; for example, the total suspended sediment load, affected by seasonal weather conditions, increased from May to August and fell dramatically from September to January [7].

4.2.3. Runoff-Sediment Relation at the Annual Scale

It can be seen in Figure 6 that the variations of annual mean discharge, SSC and fluctuation amplitude at the four stations all displayed marked decreasing trends, and the current runoff and sediment were the driest stages in the runoff-sediment cycle [57]. In particular, the changes of SSC were basically synchronous with discharge at the Toudaoguai station. Nevertheless, this was not the case at other stations because the abundant sediment supply with a fragile Loess Plateau triggered the weak correlation between SSC and Q [35,37]. Spatially, the runoff mainly comes from the upper reaches, leading to the same variation discipline of the runoff at the four stations under the operation of the Liujiaxia in 1986 and the joint control of the Liujiaxia and Longyang in 1968, while the variation rules of sediment load, mainly from the middle reaches, were different. Natural causes—such as precipitation, temperature, vegetation and human activities, such as soil and water conservation projects, the operation of dams and reservoirs, water consumption and frequent tillage—could principally account for the changes in runoff and sediment reported from many investigations [24,58].
Extreme precipitation or rainfall erosion has a tendency to increase from northwest to southeast [59,60]. Therefore, the phased and spatial characteristics in runoff and sediment load with a decreasing trend at the four stations are important to note. The years from 1960 to 1968 were all high-flow years with few exceptions, and the average annual runoff were 28.28 billion m3, 35.27 billion m3, 47.12 billion m3 and 52.84 billion m3, respectively. The years from 1969 to 1974 were basically consecutive normal flow years. Differently, dry years were obvious from 1990 to 2015 (except 2012) due to the reduction of precipitation and significant utilization of water resources [24,49]. The corresponding annual average runoff rates at the four stations were 16.16 billion m3, 19.05 billion m3, 23.84 billion m3 and 25.54 billion m3. The maximum runoff reached 4.30, 4.07, 4.68 and 6.04 times as much as the minimum values, respectively. The variation of sediment load at the Toudaoguai station was similar to the runoff, while sediment at the other three stations in the four stages of 1960–1980, 1981–1986, 1987–1998 and 1999–2015 showed significant stepwise changes. In the 1980s, the storage of dams and reservoirs dominated [61]. After 1999, the implementation of a large-scale “returning farmland to forest (grass)” policy became the main reason for the reduction of sediment [24]. The operation of the Longyangxia in 1986 reduced the quantity of runoff and increased the sediment carrying capacity, resulting in an increase in SSC from 1986 to 1999 [59]. The change in degree of runoff and sediment was unevenly distributed in space and the runoff and sediment reductions were not synchronized. Compared with 1960–1968, the annual average runoff of the four stations in 2000−2015 decreased by 41.77%, 47.36%, 50.76% and 51.84%, respectively. The corresponding reduction in sediment load were 78.05%, 86.97%, 82.05% and 91.58%, respectively. The decrement increased along the way. The maximum values were 18.73, 65.08, 44.58 and 159.07 times as much as the minimum value, respectively; thus, the change in annual sediment load was much greater than that in runoff. It is deduced that the impact of human activities especially for soil and water conservation measures on sediment load is more significant than on runoff [62]. It is of great significance for the establishment of regulation system of water and sediment in the MYRB, implement of soil and water conservation measures according to local conditions and design of major water conservancy projects to master the change rules of water and sediment and their relation at the space-time scale. Therefore, there was a significant decrease of correlation between Q and SSC at the Longmen, Tongguan and Huayuankou stations after 2000. From the analyzed data (Figure 7), the accumulated sediment load of Tongguan was the highest in all three stages, resulting in serious silt deposition. Abrupt changes occurred in 1997 at Huayuankou due to regulation of the Xiaolangdi reservoir. Due to the flood retardation and sediment release of Sanmenxia from 1965 to 1973 and the operation of storing clear water and releasing muddy water (with the inputs and outputs balanced) since 1974 [63], the slopes of the cumulative sediment load curves at the Tongguan and Huayuankou stations were basically parallel from 1974 to 1999. Compared with the runoff during 1960–1964, the runoff of Huayuankou from 1965 to 1973 and from 1974 to 1997 decreased by 19.49% and 31.04%, respectively, and the sediment load increased by 74.58% and 8.77%, respectively.

5. Conclusions

The current study presented definite equilibrium relationships among SL, SSC and runoff series of the Toudaoguai, Longmen, Tongguan, Huayuankou stations at different temporal scales. These logarithmic variables satisfied linear equilibrium relationships at the Toudaoguai station at multiple time scales, because the variation features of runoff and sediment were similar. At the other stations, the nonlinear cointegration was appropriate at the annual scale. The test of cointegration relationships reveals the quantitative relationship among SL, runoff and SSC and provides a theoretical equation for the calculation and prediction of hydrological variables. In addition, the existence of the cointegration relationship eliminates the pseudo-regression phenomenon in the subsequent study of the uncertainty relationship between variables. ECM, based on a linear cointegration equation, could improve the prediction accuracy compared with OLS. In general, ECM based on three variables was the best and worthy of recommendation in the forecasting of hydrological variables.
The different hydrological conditions induced by the variations of geomorphologic agent, available sediment, reservoir regulation and other natural and human factors result in uncertain relationships between runoff and sediment. The research on the uncertainty relationship between water and sediment of the Yellow River is of great significance to the construction of the water and sediment control system in the MYRB and the design of major water conservancy projects. The relation between discharge and SSC at the Toudaoguai station could be described by a cubic polynomial function at the within-event scale. Nevertheless, at other three stations, low energy flow with limited fine sediments resulted in a lower SSC with respect to Q < 1000 m3/s and emerged rating curves in the form of a power-law function. The SL, runoff and SSC satisfied the linear cointegration relationship at multiple time scales; When Q > 1000 m3/s, energy flow with different particle size ratio promotes the formation of hyperconcentrated flow. The SL, runoff and SSC satisfied the linear cointegration relationship at multiple time scales. Meanwhile, the SRCs between the maximum SSC values for a given discharge displayed a linear form once the flow was saturated and stable. Therefore, it is helpful for the calculation of the maximum SSC value for a particular discharge. Meanwhile, effective desilting can be achieved by controlling the volume of discharge (threshold: 1000 m3/s) and the proportion of sediment particle size. Different types of hysteresis loops occurred in different proportions at the four stations due to the different sources of runoff and sediment. At the monthly-seasonal scale, the sediment rating curve derived from a linear regression fitted the data well in the ice-flood season. The rating curve at the Toudaoguai station could also be described by a cubic polynomial function during the rainy season, but the correlations between SSC and Q were very low from Longmen to Huayuankou stations, due to a complex supply of sediment and the formation of hyperconcentrated flow. The hysteresis at the Toudaoguai station exhibited figure-eight patterns, and hysteresis loops with clockwise trends were found for the other stations at a monthly scale. At the annual scale, phased and spatial characteristics in the runoff and sediment load were evident. In addition, the degree of change in runoff and sediment was unevenly distributed in space and the runoff and sediment reduction were not synchronized, because of the impact of human activities on sediment load, which was more significant than on runoff.

Author Contributions

Conceptualization, X.W.; Methodology, X.W.; Software, D.L.; Formal Analysis, X.W. and D.L.; Data Curation, X.Q.; Writing—Original Draft Preparation, D.L. and X.Q.; Writing—Review & Editing, X.W. and X.Y.; Supervision, X.W.; Project Administration, P.Z.; Funding Acquisition, X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China (grant number 2018YFC1508403), the National Natural Science Foundation of China (grant number 51579173), the National Key Research and Development Program of China (grant number 2016YFC0401701) and the Science Fund for Creative Research Groups of the National Natural Science Foundation of China (grant number 51621092).

Acknowledgments

We wish to express our deep gratitude to the Yellow River Conservancy Commission for permission to access to the hydrological data. Meanwhile, we wish to acknowledge the assistance of all editors and reviewers and thank Sara J. Mason, MSc, from Liwen Bianji, (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sun, P.C.; Wu, Y.P.; Gao, J.N.; Yao, Y.Y.; Zhao, F.B.; Lei, X.H.; Qiu, L.J. Shifts of sediment transport regime caused by ecological restoration in the Middle Yellow River Basin. Sci. Total Environ. 2020, 698, 1–13. [Google Scholar] [CrossRef] [PubMed]
  2. Sadeghi, S.H.; Singh, V.P.; Kiani-Harchegani, M.; Asadi, H. Analysis of sediment rating loops and particle size distributions to characterize sediment source at mid-sized plot scale. Catena 2018, 167, 221–227. [Google Scholar] [CrossRef]
  3. Smith, H.G.; Dragovich, D. Interpreting sediment delivery processes using suspended sediment-discharge hysteresis patterns from nested upland catchments, south-eastern Australia. Hydrol. Process. 2009, 23, 2415–2426. [Google Scholar] [CrossRef]
  4. Lu, J.F. Effect of basin morphology on sediment yield in the middle reaches of the Yellow River. Geogr. Res. 2002, 21, 171–178. (In Chinese) [Google Scholar] [CrossRef]
  5. Fan, X.L.; Shi, C.X.; Shao, W.W.; Zhou, Y.Y. The suspended sediment dynamics in the Inner-Mongolia reaches of the upper Yellow River. Catena 2013, 109, 72–82. [Google Scholar] [CrossRef]
  6. Fang, H.Y.; Cai, Q.G.; Chen, H.; Li, Q.Y. Temporal changes in suspended sediment transport in a gullied loess basin: The lower Chabagou Creek on the Loess Plateau in China. Earth Surf. Process. Landf. 2008, 33, 1977–1992. [Google Scholar] [CrossRef]
  7. Sun, L.Y.; Yan, M.; Cai, Q.G.; Fang, H.Y. Suspended sediment dynamics at different time scales in the Loushui River, south-central China. Catena 2016, 136, 152–161. [Google Scholar] [CrossRef]
  8. Vercruysse, K.; Grabowski, R.C.; Rickson, R.J. Suspended sediment transport dynamics in rivers: Multi-scale drivers of temporal variation. Earth Sci. Rev. 2017, 166, 38–52. [Google Scholar] [CrossRef] [Green Version]
  9. Oeurng, C.; Sauvage, S.; Sanchez-Perez, J.M. Dynamics of suspended sediment transport and yield in a large agricultural catchment, southwest France. Earth Surf. Process. Landf. 2010, 35, 1289–1301. [Google Scholar] [CrossRef] [Green Version]
  10. An, C.H.; Lu, J.; Qian, Y.; Luo, Q.S.; Cui, Z.H. Spatial-temporal distribution characteristic and course of sedimentation in the Ningxia-Inner Mongolia reaches of the Yellow River. J. Hydraul. Eng. 2018, 49, 195–206. (In Chinese) [Google Scholar] [CrossRef]
  11. Pan, B.T.; Pang, H.L.; Zhang, D.; Guan, Q.Y.; Wang, L.; Li, F.Q.; Guan, W.Q.; Cai, A.; Sun, X.Z. Sediment grain-size characteristics and its source implication in the Ningxia-Inner Mongolia sections on the upper reaches of the Yellow River. Geomorphology 2015, 246, 255–262. [Google Scholar] [CrossRef]
  12. Ouyang, C.B.; Wang, W.L.; Tian, Y.; Tian, S.M. Evaluation on the variation of water-sediment and human activities in the He-Long Reach of the Yellow River over the past 60 years. J. Sediment Res. 2016, 4, 55–61. (In Chinese) [Google Scholar] [CrossRef]
  13. Liu, J.X.; Li, Z.G.; Zhang, X.P.; Li, R.; Liu, X.C.; Zhang, H.Y. Responses of vegetation cover to the Grain for Green Program and their driving forces in the He-Long region of the middle reaches of the Yellow River. J. Arid Land 2013, 5, 511–520. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, J.P.; Zhao, Y.; Ding, Z.H. Research on the Relationships between rainfall and meteorological yield in irrigation district. Water Resour. Manag. 2014, 28, 1689–1702. [Google Scholar] [CrossRef]
  15. Zhang, X.; Yu, G.Q.; Li, Z.B.; Li, P. Experimental study on slope runoff, erosion and sediment under different vegetation types. Water Resour. Manag. 2014, 28, 2415–2433. [Google Scholar] [CrossRef]
  16. Hu, J.F.; Zhao, G.J.; Mu, X.M.; Tian, P.; Gao, P.; Sun, W.Y. Quantifying the impacts of human activities on runoff and sediment load changes in a Loess Plateau catchment, China. J. Soils Sediments 2019, 19, 3866–3880. [Google Scholar] [CrossRef]
  17. Guo, W.X.; Li, Y.; Wang, H.X.; Cha, H.F. Temporal variations and influencing factors of river runoff and sediment regimes in the Yangtze River, China. Desalin. Water Treat. 2020, 174, 258–270. [Google Scholar] [CrossRef]
  18. Aich, V.; Zimmermann, A.; Elsenbeer, H. Quantification and interpretation of suspended-sediment discharge hysteresis patterns: How much data do we need? Catena 2014, 122, 120–129. [Google Scholar] [CrossRef]
  19. Syvitski, J.P.; Morehead, M.D.; Bahr, D.B.; Mulder, T. Estimating fluvial sediment transport: The rating parameters. Water Resour. Res. 2000, 36, 2747–2760. [Google Scholar] [CrossRef] [Green Version]
  20. Zhang, J.P.; Li, H.B.; Shi, X.X.; Hong, Y. Wavelet-nonlinear cointegration prediction of irrigation water in the irrigation district. Water Resour. Manag. 2019, 33, 2941–2954. [Google Scholar] [CrossRef]
  21. Engle, R.F.; Granger, C.W.J. Cointegration and error correction: Representation, estimation and testing. Econometrica 1987, 55, 251–276. [Google Scholar] [CrossRef]
  22. Zhang, J.P.; Zhao, Y.; Xiao, W.H. Multi-resolution cointegration prediction for runoff and sediment load. Water Resour. Manag. 2015, 29, 3601–3613. [Google Scholar] [CrossRef]
  23. Zhang, J.P.; Li, Y.Y.; Zhao, Y.; Hong, Y. Wavelet-cointegration prediction of irrigation water in the irrigation district. J. Hydrol. 2017, 544, 343–351. [Google Scholar] [CrossRef]
  24. Zhao, G.J.; Tian, P.; Mu, X.M.; Jiao, J.Y.; Wang, F.; Gao, P. Quantifying the impact of climate variability and human activities on streamflow in the middle reaches of the Yellow River basin, China. J. Hydrol. 2014, 519, 387–398. [Google Scholar] [CrossRef]
  25. Peng, J.; Chen, S.L.; Dong, P. Temporal variation of sediment load in the Yellow River basin, China, and its impacts on the lower reaches and the river delta. Catena 2010, 83, 135–147. [Google Scholar] [CrossRef]
  26. Wang, X.J.; Engel, B.; Yuan, X.M.; Yuan, P.X. Variation analysis of streamflows from 1956 to 2016 along the Yellow River, China. Water 2018, 10, 1123. [Google Scholar] [CrossRef] [Green Version]
  27. Johansen, S. Likelihood-Based Inference in Cointegrated Vector Autoregressive Models; Oxford University Press: Oxford, UK, 1995. [Google Scholar]
  28. Leybourne, S.; Kim, T.H.; Newbold, P. A more powerful modification of Johansen’s cointegration tests. Appl. Econ. 2008, 40, 725–729. [Google Scholar] [CrossRef]
  29. Haug, A.A.; Basher, S.A. Linear or nonlinear cointegration in the purchasing power parity relationship? Appl. Econ. 2011, 43, 185–196. [Google Scholar] [CrossRef]
  30. Zhang, L.Y.; Zhang, L.P.; Cao, F.L.; Song, X.Y. Annual runoff forecasting research based on the theory of cointegration and error correction model. Wuhan Univ. Eng. Sci. 2006, 39, 6–9. (In Chinese) [Google Scholar] [CrossRef]
  31. Hurst, H.E. Long-term storage capacity of reservoirs. Trans. Am. Soc. Civ. Eng. 1951, 116, 770–779. [Google Scholar]
  32. Liu, D.H.; Zhang, S.Y. Nonlinear error correction model and forecasting based on wavelet neural networks. Control Decis. 2006, 21, 1114–1118. (In Chinese) [Google Scholar] [CrossRef]
  33. Asselman, N.E.M. Fitting and interpretation of sediment rating curves. J. Hydrol. 2000, 234, 228–248. [Google Scholar] [CrossRef]
  34. Hassanzadeh, H.; Bajestan, M.S.; Paydar, G.R. Performance evaluation of correction coefficients to optimize sediment rating curves on the basis of the Karkheh dam reservoir hydrography, West Iran. Arab. J. Geosci. 2018, 11, 9. [Google Scholar] [CrossRef]
  35. Fan, X.; Shi, C.; Zhou, Y.; Shao, W. Sediment rating curves in the Ningxia-Inner Mongolia reaches of the upper Yellow River and their implications. Quat. Int. 2012, 282, 152–162. [Google Scholar] [CrossRef]
  36. Gao, P.; Mu, X.M.; Wang, F.; Li, R. Changes in streamflow and sediment discharge and the response to human activities in the middle reaches of the Yellow River. Hydrol. Earth Syst. Sci. 2011, 15, 1–10. [Google Scholar] [CrossRef] [Green Version]
  37. Hashimoto, H.; Takaoka, H.; Ikematsu, S. Hyper-concentrated flows in tributaries of the middle Yellow River. Monit. Simul. Prev. Remediat. Dense Debris Flows 2006, 90, 353–362. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, S.J.; Yan, M.; Yan, Y.X.; Shi, C.X.; He, L. Contributions of climate change and human activities to the changes in runoff increment in different sections of the Yellow River. Quat. Int. 2012, 282, 66–77. [Google Scholar] [CrossRef]
  39. Moosa, I.A.; Vaz, J.J. Cointegration, error correction and exchange rate forecasting. J. Int. Financ. Mark. Inst. Money 2016, 44, 21–34. [Google Scholar] [CrossRef]
  40. Christoffersen, P.F.; Diebold, F.X. Contegration and long-horizon forecasting. J. Bus. Econ. Stat. 1998, 16, 450–458. [Google Scholar]
  41. Jacobson, T.; Jansson, P.; Verdin, A.; Warne, A. Monetary policy analysis and inflation targeting in a small open economy: A VAR approach. J. Appl. Econom. 2001, 16, 487–520. [Google Scholar] [CrossRef]
  42. Hu, J.F.; Gao, P.; Mu, X.M.; Zhao, G.J.; Sun, W.Y.; Li, P.F.; Zhang, L.M. Runoff-sediment dynamics under different flood patterns in a Loess Plateau catchment, China. Catena 2019, 173, 234–245. [Google Scholar] [CrossRef]
  43. Xu, J.X. Optimal grain-size composition of hyperconcentrated flows in high-intensity coarse sediment producing area of the middle Yellow River basin and its implications in geomorphology. J. Sediment Res. 1999, 5, 3–5. [Google Scholar] [CrossRef]
  44. Xu, J.X. The optimal grain size composition of suspended sediment of hyperconcentrated flow in the middle Yellow River. Int. J. Sediment Res. 1997, 12, 170–176. [Google Scholar]
  45. Xu, J.X. Implication of relationships among suspended sediment size, water discharge and suspended sediment concentration: The Yellow River basin, China. Catena 2002, 49, 289–307. [Google Scholar] [CrossRef]
  46. Xu, J.X. Zonal distribution of river erosion and sediment yield in China. Chin. Sci. Bull. 1994, 39, 1356–1361. [Google Scholar]
  47. Fang, H.; Li, Q.; Cai, Q.; Liao, Y. Spatial scale dependence of sediment dynamics in a gullied rolling loess region on the Loess Plateau in China. Environ. Earth Sci. 2011, 64, 693–705. [Google Scholar] [CrossRef]
  48. Xu, J.X.; Cheng, D.S. Relation between the erosion and sedimemtation zones in the Yellow River, China. Geomorphology 2002, 48, 365–382. [Google Scholar] [CrossRef]
  49. Li, B.Q.; Liang, Z.M.; Zhang, J.Y.; Wang, G.Q.; Zhao, W.M.; Zhang, H.Y.; Wang, J.; Hu, Y.M. Attribution analysis of runoff decline in a semiarid region of the Loess Plateau, China. Theor. Appl. Climatol. 2018, 131, 845–855. [Google Scholar] [CrossRef]
  50. Miao, C.Y.; Shi, W.; Chen, X.H.; Yang, L. Spatio-temporal variability of streamflow in the Yellow River: Possible causes and implications. Hydrol. Sci. J. 2012, 57, 1355–1367. [Google Scholar] [CrossRef] [Green Version]
  51. Hou, S.H.; Wang, P.; Chu, W.B. Actions of the Longyangxia and Liujiaxia Reservoirs on the Runoff and Sediment of Yellow River; Yellow River Conservancy Press: Zhengzhou, China, 2005. [Google Scholar]
  52. Gentile, F.; Bisantino, T.; Corbino, R.; Milillo, F.; Romano, G.; Liuzzi, G.T. Monitoring and analysis of suspended sediment transport dynamics in the Carapelle torrent (Southern Italy). Catena 2010, 80, 1–8. [Google Scholar] [CrossRef]
  53. Fang, N.F.; Shi, Z.H.; Li, L.; Jiang, C. Rainfall, runoff, and suspended sediment delivery relationships in a small agricultural watershed of the Three Gorges area, China. Geomorphology 2011, 135, 158–166. [Google Scholar] [CrossRef]
  54. Qian, N.; Wan, Z.H. Mechanics of Sediment Transport; Science Press: Beijing, China, 1983. [Google Scholar]
  55. Xu, J.X. Erosion caused by hyperconcentrated flow on the Loess Plateau of China. Catena 1999, 36, 1–19. [Google Scholar]
  56. Wang, Y.J.; Wu, B.S.; Zhong, D.Y. Adjustment in the main-channel geometry of the lower Yellow River before and after the operation of the Xiaolangdi Reservoir from 1986 to 2015. J. Geogr. Sci. 2020, 30, 468–486. [Google Scholar] [CrossRef]
  57. Yao, W.Y.; Jiao, P. The change of water and sediment in the Yellow River and its research prospect. Soil Water Conserv. China 2016, 9, 55–63. (In Chinese) [Google Scholar] [CrossRef]
  58. Wang, S.; Fu, B.J.; Liang, W.; Liu, Y.; Wang, Y.F. Driving forces of changes in the water and sediment relationship in the Yellow River. Sci. Total Environ. 2017, 576, 453–461. [Google Scholar] [CrossRef]
  59. Zhang, Y.; Xia, J.; She, D.X. Spatiotemporal variation and statistical characteristic of extreme precipitation in the middle reaches of the Yellow River Basin during 1960–2013. Theor. Appl. Climatol. 2019, 135, 391–408. [Google Scholar] [CrossRef]
  60. Dang, S.Z.; Yao, M.F.; Liu, X.Y.; Dong, G.T. Variations and statistical probability characteristic analysis of extreme precipitation in the Hekouzhen-Longmen Region of the Yellow River, China. Asia Pac. J. Atmos. Sci. 2019, 55, 641–655. [Google Scholar] [CrossRef]
  61. Wenyi, Y.; Dachuan, R.; Jiangnan, C. Recent changes in runoff and sediment regimes and future projections in the Yellow River basin. Adv. Water Sci. 2013, 24, 607–616. [Google Scholar] [CrossRef]
  62. Zhao, G.J.; Mu, X.M.; Tian, P.; Wang, F.; Gao, P. The variation of streamflow and sediment flux in the middle reaches of Yellow River over the past 60 years and the influencing factors. Resour. Sci. 2012, 34, 1070–1078. (In Chinese) [Google Scholar]
  63. Mao, Z.P.; Peng, W.Q.; Zhou, H.D. Study on the influence of operation mode of Sanmenxia rservoir on wetland ecosystem in reservoir area. Water Resour. Dev. Res. 2005, 9, 12–17. (In Chinese) [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the middle reaches of the Yellow River.
Figure 1. Schematic diagram of the middle reaches of the Yellow River.
Water 12 02589 g001
Figure 2. Rating curves for SSC and discharge at the (a) Toudaoguai, (b) Longmen, (c) Tongguan and (d) Huayuankou station during from 1960 to 2005.
Figure 2. Rating curves for SSC and discharge at the (a) Toudaoguai, (b) Longmen, (c) Tongguan and (d) Huayuankou station during from 1960 to 2005.
Water 12 02589 g002
Figure 3. Seasonal rating curves for discharge and SSC in ice-flood seasons at (a) Toudaoguai, (c) Longmen, (e) Tongguan, (g) Huayuankou stations and in rainy seasons at (b) Toudaoguai, (d) Longmen, (f) Tongguan, (h) Huayuankou stations from 1960 to 2015.
Figure 3. Seasonal rating curves for discharge and SSC in ice-flood seasons at (a) Toudaoguai, (c) Longmen, (e) Tongguan, (g) Huayuankou stations and in rainy seasons at (b) Toudaoguai, (d) Longmen, (f) Tongguan, (h) Huayuankou stations from 1960 to 2015.
Water 12 02589 g003aWater 12 02589 g003b
Figure 4. Monthly average hysteresis loop of SSC and discharge during the period of (a) 1960−1968, (b) 1969−1986, (c) 1987−2015 at Toudaoguai station.
Figure 4. Monthly average hysteresis loop of SSC and discharge during the period of (a) 1960−1968, (b) 1969−1986, (c) 1987−2015 at Toudaoguai station.
Water 12 02589 g004
Figure 5. Monthly average hysteresis loop of SSC and discharge during the period of 1960−2015 at (a) Longmen, (b) Tonguan, (c) Huayuankou stations.
Figure 5. Monthly average hysteresis loop of SSC and discharge during the period of 1960−2015 at (a) Longmen, (b) Tonguan, (c) Huayuankou stations.
Water 12 02589 g005
Figure 6. Changes in annual mean water discharge and suspended sediment concentration at Toudaoguai station.
Figure 6. Changes in annual mean water discharge and suspended sediment concentration at Toudaoguai station.
Water 12 02589 g006
Figure 7. Cumulative runoff and sediment load (SL) from 1960 to 2015.
Figure 7. Cumulative runoff and sediment load (SL) from 1960 to 2015.
Water 12 02589 g007
Table 1. Unit root test results of sediment load (SL), runoff, and suspended sediment concentration (SSC) sequences.
Table 1. Unit root test results of sediment load (SL), runoff, and suspended sediment concentration (SSC) sequences.
SequencesTypeNull HypothesisADF Value%5 ThresholdP Value of ADF StatisticP Value of Trend ItemConclusion
SLTrend, interceptexist unit root−5.49−3.520.00030.0004Trendy, stationary
runoffTrend, interceptexist unit root−4.78−3.520.00200.0035Trendy, stationary
SSCTrend, interceptexist unit root−5.49−3.520.00030.0002Trendy, stationary
D(SL) 1noneexist unit root−8.39−1.950.0000stationary
D(runoff)noneexist unit root−8.23−1.950.0000stationary
D(SSC)noneexist unit root−8.00−1.950.0000stationary
1 D() represents the first order difference of the corresponding variable.
Table 2. P-step lag values of AIC 1, SC 2 and HQ 3 statistics.
Table 2. P-step lag values of AIC 1, SC 2 and HQ 3 statistics.
LagAICSCHQ
1−6.11−5.72−5.97
2−6.00−5.22−5.72
3−5.94−4.78−5.53
4−6.05−4.49−5.49
5−5.98−4.04−5.29
6−6.10−3.77−5.27
1 AIC represents the Akaike Information Criterion; 2 SC represents the Schwarz Criterion; 3 HQ represents the Hannan-quinn Criterion.
Table 3. The results of the Johansen cointegration test.
Table 3. The results of the Johansen cointegration test.
Trace TestMaximum Eigenvalue
Hypothesized No. of CE(s) (Cointegration Equation(s))Trace Statistic0.05 Critical Valuep ValueMax-Eigen Statistic0.05 Critical Valuep Value
None58.1229.800.000030.8421.130.0016
At most 127.2915.490.000616.8014.260.0194
At most 210.493.840.001210.493.840.0012
Table 4. Predicted results of sediment transport at Toudaoguai station.
Table 4. Predicted results of sediment transport at Toudaoguai station.
YearMeasured ValueOLS(Three Variables) 1ECM1(Three Variables) 2OLS(SL, SSC) 3ECM2(SL, SSC) 4OLS(SL, Runoff)ECM3(SL, Runoff)
Calculated ValueRe5Calculated ValueReCalculated ValueReCalculated ValueReCalculated ValueReCalculated ValueRe
20010.200.189.12%0.189.92%0.1714.66%0.1810.84%0.2313.57%0.200.36%
20020.270.3324.29%0.3114.40%0.3735.65%0.3633.45%0.256.72%0.257.46%
20030.280.293.98%0.280.61%0.3213.02%0.3110.24%0.2220.43%0.2318.64%
20040.240.2817.26%0.3023.62%0.2817.33%0.2816.90%0.2714.64%0.3024.69%
20050.400.4512.57%0.436.74%0.4818.21%0.4613.19%0.378.91%0.379.65%
1 OLS (three variables) represents the ordinary least squares (OLS) regression model among three variables; 2 ECM1 (three variables) represents the ECM1 among three variables; 3 OLS (SL, SSC) represents the OLS model between SL and SSC; 4 ECM2 (SL, SSC) represents the ECM2 between SL and SSC; 5 Re represents the relative error between variables. The other parameters have similar meanings.
Table 5. Predicted results of runoff at Toudaoguai station.
Table 5. Predicted results of runoff at Toudaoguai station.
YearMeasured ValueOLS(Runoff, SSC)ECM4(Runoff, SSC)
Calculated ValueReCalculated ValueRe
2001113.28100.6511.15%110.702.28%
2002122.75144.6617.85%144.5817.78%
2003115.57135.1416.93%128.3111.02%
2004127.61127.750.11%121.474.81%
2005150.21164.649.61%159.326.07%
Table 6. Hurst exponents of hydrological series SL, runoff, SSC.
Table 6. Hurst exponents of hydrological series SL, runoff, SSC.
SequenceSLRunoffSSC
Longmen0.660.710.63
Tongguan0.680.920.55
Huanyuankou0.790.910.89
Table 7. Fractal dimension of hydrological series SL, runoff, SSC.
Table 7. Fractal dimension of hydrological series SL, runoff, SSC.
SequenceSLRunoffSSC
Longmen1.341.291.37
Tongguan1.321.081.45
Huanyuankou1.211.091.11
Table 8. Related parameters of SL−runoff−SSC cointegration equation.
Table 8. Related parameters of SL−runoff−SSC cointegration equation.
LongmenTongguanHuayuankou
a n 1 b n 2 ω n 3 v n m 4 a n b n ω n v n m a n b n ω n v n m
0.864−0.7061.082−0.3752.333−1.4991.1180.2900.8150.2601.443−0.336−0.649−1.0370.7890.6210.8171.451
1.968−0.156−0.829−2.0871.244−1.5790.7800.401−0.921−0.368−1.209−0.5420.328−1.5730.6740.584−0.230−0.372
0.144−1.121−1.1940.9740.716−0.4030.733−0.447−1.149−0.343−2.0140.2501.483−0.3540.206−1.7590.2450.943
−0.626−0.5301.1382.009−1.166−2.0420.015−2.1760.292−2.253−0.493−0.262−0.942−0.335−0.504−0.224−0.7330.635
1.849−0.086−0.3730.5440.0740.642−1.064−1.222−1.240−0.8930.704−1.816−0.870−0.542−0.531−0.459−2.289−1.505
2.918−0.0110.231−0.782−1.112−0.2310.4600.2800.9711.397−0.8030.083−0.366−0.151−2.844−2.1800.5370.177
1.549−0.2690.5200.805−1.0260.815−1.2942.442−1.363−0.628−0.810−2.7310.417−1.048−1.315−1.127−1.4790.230
−0.9741.0241.101−1.680−1.136−1.3611.5570.7610.5640.357−1.6190.238−0.6590.5591.042−0.5252.5840.416
h ( ( x b n ) / a n ) (n = 1, 2,…N) is the wavelet function; 1 a n is the expansion parameter; 2 b n is the translation parameter; 3 ω n is the connection weight from hidden layer to the output; 4 v n m is the connection weight from input layer to hidden layer.
Table 9. Related parameters of SL−runoff cointegration equation.
Table 9. Related parameters of SL−runoff cointegration equation.
LongmenTongguanHuayuankou
a n b n ω n v n m a n b n ω n v n m a n b n ω n v n m
1.1050.6450.4271.1650.8340.5150.1490.382−0.903−1.3311.0370.024−0.8850.437−1.881
−0.5200.441−0.3471.1360.1170.178−0.376−0.832−0.1890.1741.134−0.573−0.7150.467−0.557
−0.1990.566−1.803−0.207−0.2730.095−1.1931.198−0.3250.854−1.2430.752−1.895−1.593−1.623
−1.7991.835−1.1210.756−2.1320.4150.6750.8912.1240.351−0.801−0.670−1.1171.222−1.236
−0.8600.5210.7841.6300.5731.035−0.287−0.074−0.606−0.7760.159−0.2980.6180.197−0.217
0.7170.3880.742−0.847−0.3775.8030.0040.267−0.325−1.547−0.380−0.750−1.139−1.693−0.373
−0.065−0.095−1.098−0.540−1.677−0.454−0.715−0.848−1.6142.2281.1590.1081.8831.297−0.327
0.839−1.3270.556−1.702−1.0240.953−1.569−0.533−0.6442.8321.4421.0560.529−1.6540.080
The parameters have the same meaning as Table 8.
Table 10. Related parameters of SL−SSC cointegration equation.
Table 10. Related parameters of SL−SSC cointegration equation.
LongmenTongguanHuayuankou
a n b n ω n v n m a n b n ω n v n m a n b n ω n v n m
−0.125−0.098−1.171−0.0481.304−0.9490.326−0.8753.2510.7360.8770.9422.041−0.9220.627
0.362−0.2500.5131.5922.143−0.462−0.445−1.033−0.2581.3161.095−0.8360.860−0.940−0.872
−0.2700.773−1.9441.065−0.8923.345−0.2691.797−0.590−4.2180.258−0.4380.7950.403−0.477
−0.378−0.4390.9420.218−1.966−0.6970.8631.900−0.271−0.3360.5432.995−2.532−2.595−0.096
0.028−0.7730.395−1.320−0.126−0.522−0.733−1.5141.362−0.208−0.686−0.6080.7920.255−0.518
1.3960.411−0.1360.8291.4350.8191.2512.6070.2080.701−1.1071.9340.633−0.831−0.369
0.0010.0021.0510.547−0.1980.3540.8511.7020.758−0.443−0.3561.048−0.4841.1631.099
−0.7831.2150.053−0.971−0.7030.8960.2700.224−0.9740.3520.198−0.5460.940−0.9580.707
The parameters have the same meaning as Table 8.
Table 11. Related parameters of runoff−SSC cointegration equation.
Table 11. Related parameters of runoff−SSC cointegration equation.
LongmenTongguanHuayuankou
a n b n ω n v n m a n b n ω n v n m a n b n ω n v n m
0.543−0.946−1.2980.6874.0251.0270.0670.492−1.0860.122−1.7890.697−0.0550.511−0.642
1.366−1.0121.0361.6133.357−1.0740.308−0.607−0.2780.676−0.1460.847−0.950−0.704−1.174
0.767−2.9631.494−2.682−1.0071.7590.6140.2672.8390.4250.2610.5181.1221.271−0.696
−0.3951.4721.3310.7533.346−0.395−1.461−0.532−1.0801.2970.5280.087−1.2040.120−1.374
0.1740.3201.2110.5781.204−1.163−0.3791.3140.488−3.2370.7990.832−0.547−2.0360.062
−0.716−0.818−2.101−1.2260.619−0.586−1.913−0.721−1.118−0.8640.235−0.885−1.124−1.0410.485
0.9041.3180.5500.5141.2081.010−0.428−0.2851.2640.710−0.9580.1130.250−2.124−0.289
−1.062−1.7940.9180.145−0.5011.0091.021−0.8480.308−0.010−0.5901.4450.7920.8020.124
The parameters have the same meaning as Table 8.

Share and Cite

MDPI and ACS Style

Wang, X.; Li, D.; Yuan, X.; Qi, X.; Zhang, P. Analysis of Runoff-Sediment Cointegration and Uncertainty Relations at Different Temporal Scales in the Middle Reaches of the Yellow River, China. Water 2020, 12, 2589. https://doi.org/10.3390/w12092589

AMA Style

Wang X, Li D, Yuan X, Qi X, Zhang P. Analysis of Runoff-Sediment Cointegration and Uncertainty Relations at Different Temporal Scales in the Middle Reaches of the Yellow River, China. Water. 2020; 12(9):2589. https://doi.org/10.3390/w12092589

Chicago/Turabian Style

Wang, Xiujie, Dandan Li, Ximin Yuan, Xiling Qi, and Pengfei Zhang. 2020. "Analysis of Runoff-Sediment Cointegration and Uncertainty Relations at Different Temporal Scales in the Middle Reaches of the Yellow River, China" Water 12, no. 9: 2589. https://doi.org/10.3390/w12092589

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop