Applsci 12 00647 v2
Applsci 12 00647 v2
Applsci 12 00647 v2
sciences
Article
Prediction of River Sediment Transport Based on Wavelet
Transform and Neural Network Model
Zongyu Li 1 , Zhilin Sun 1 , Jing Liu 1 , Haiyang Dong 2 , Wenhua Xiong 3 , Lixia Sun 2 and Hanyu Zhou 4, *
Abstract: The sedimentation problem is one of the critical issues affecting the long-term use of rivers,
and the study of sediment variation in rivers is closely related to water resource, river ecosystem
and estuarine delta siltation. Traditional research on sediment variation in rivers is mostly based
on field measurements and experimental simulations, which requires a large amount of human
and material resources, many influencing factors and other restrictions. With the development of
computer technology, intelligent approaches have been applied to hydrological models to establish
small information in river areas. In this paper, considering the influence of multiple factors on
sediment transport, the validity of predicting sediment transport combined with wavelet transforms
and neural network was analyzed. The rainfall and runoff cycles are extracted and decomposed
into time series sub-signals by wavelet transforms; then, the data post-processing is used as the
neural network training set to predict the sediment model. The results show that wavelet coupled
Citation: Li, Z.; Sun, Z.; Liu, J.; Dong, neural network model effectively improves the accuracy of the predicted sediment model, which can
H.; Xiong, W.; Sun, L.; Zhou, H. provide a reference basis for river sediment prediction.
Prediction of River Sediment
Transport Based on Wavelet Keywords: wavelet transform; sediment prediction; rainfall and runoff; estuaries
Transform and Neural Network
Model. Appl. Sci. 2022, 12, 647.
https://doi.org/10.3390/
app12020647 1. Introduction
Academic Editor: Kelin Hu Due to the tremendous variations in midstream and downstream areas of rivers and
in the vicinity of estuaries, under the influence of climate change or strong anthropogenic
Received: 19 November 2021
factors at present, high efficiency and accurate tools are urgently needed for describing and
Accepted: 5 January 2022
predicting runoff and sediment movement conditions of rivers [1–3]. The hydro-sediment
Published: 10 January 2022
variability of rivers has brought significant impacts on river evolution, hydropower plant
Publisher’s Note: MDPI stays neutral construction, river ecology, channel development and estuarine delta siltation [4]. Sediment
with regard to jurisdictional claims in load is transported by water flow to tributaries, reservoirs and affected by the processes of
published maps and institutional affil- sediment trapping, and eventually outflow to the ocean. These processes can be difficult
iations. to measure. In the midstream and downstream areas of the rivers, sediment-carrying
capacity of flow is weakened by the interception of dams and other hydraulic facilities,
which has brought about fundamental alterations in the runoff processes and sediment
movement conditions of the rivers as well as in the ecological. Whereas, the transport
Copyright: © 2022 by the authors.
of sediment through an estuary, which has emerged as one of the potential threats to the
Licensee MDPI, Basel, Switzerland.
This article is an open access article
coastal areas, as the change of river sediment supply. The changes in sediment transport
distributed under the terms and
not only seriously affect the construction of water resources, but also have a profound
conditions of the Creative Commons impact on the production and life of people in riverine areas. Among them, modeling
Attribution (CC BY) license (https:// of runoff-sediment processes is highly variable and nonlinear in nature. Bonaldo et al.
creativecommons.org/licenses/by/ noted that sediment discharge from natural rivers is one of the most powerful drivers of
4.0/). coastal morphology evolution and proposed the use of integrated information resources
ical and water resources systems with multi-timescale variability. This paper analyzes the
precipitation and runoff changes of the Lancang River system based on wavelet transform,
and uses the first main cycle of rainfall and runoff as a training set to provide relevant
flow path and hydrological information for predicted sediment in the basin, which helps
to understand the evolution of the river to improve water conservancy construction and
ecological civilization in inland and coastal areas, and provides a reference basis for water
resources utilization.
2. Methods
2.1. Data
Mengsheng hydrological station, located in Lincang City, Yunnan Province (99◦ 25E, 23◦ 22N),
is the main stream hydrological station of Nanbi River in the Lancang Water System, with a
catchment area of 1766 km2 , inter-annual average precipitation of 1687.9 mm, inter-annual
average runoff of 161,766,000 m3 and inter-annual average flow of 51.3 m3 /s, in which
water resources of 334,269,000 m3 are available in Lincang city area.
The hydrological data contains a total of 54 years of observed data for rainfall, runoff
and sediment from 1964 to 2017 at Mengsheng hydrological station in the Lancang River
system. Before studying the relationship between rainfall, runoff and sediment, the hydro-
logical data will be preprocessed [6]:
x (i ) − min x(i)
y (i ) = , i = 1, 2, · · · , a (1)
max x(i) − min x(i)
where max x(i) , min x(i) are the maximum and minimum values of data x (i ), respectively; a
is the number of the data.
t2
ψ(t) = e− 2 eiwt , w ∈ R (3)
where ψ(t) is called the fundamental wavelet and can compose a transform function by
scaling and shifting:
t−n
1
ψm,n (t) = |m|− 2 ψ (4)
m
where ψm,n (t) is sub-wavelet function; m is a non-zero real number; n is a real number.
The discrete wavelet transform is applied to the hydrological series of rainfall and
runoff by the convolution of a wavelet function with awaiting signal:
N
k∆t − n
1
W f (m, n) = |m|− 2 ∆t ∑ f (k∆t)ψ
m
, k = 1, 2, · · · , N (5)
k =1
where W f (m, n) is the wavelet transforms coefficients; N is the number of discrete points.
The real part and mode square of wavelet transform coefficients are used for hydro-
logical period determination. Additionally, the period of the hydrologic time series was
reflected by the wavelet variance with different scales:
N
1
Var (m) =
N ∑ W 2f (m, n) (6)
k =1
Appl. Sci. 2022, 12, 647 4 of 11
W f = A + D1 + D2 + . . . D N (8)
where a means the number of samples and b represents the number of input variables; θ
is harmonic parameter to make the model work best, which is searched to minimize the
cost function. The most appropriate parameter θ was defined when the cost function P(θ )
was minimized. x (i) and y(i) represent respectively the value of all input variables and
measured value for samples i.
Appl. Sci. 2022, 12, 647 5 of 11
Figure 2. (a) Real part of rainfall wavelet transform coefficients; (b) Square of rainfall wavelet
transform coefficient mode; (c) Real part of runoff wavelet transform coefficients; (d) Square of
rainfall wavelet transform coefficient mode.
It is clear from Figure 2 that there is an obvious interannual variation in the rainfall
and runoff series, and both of them are nearly the same, indicating that runoff is mainly
controlled by rainfall in the region. Additionally, it can be found that there are four main
cycle variations about rainfall and runoff, which are 3–6 years, 8–12 years, 15–20 years and
25–30 years, respectively. The annual rainfall and runoff period, mainly on the 3–6 years,
was more significant from 1972 to 1980, accompanied by relatively ambiguous abundance
and depletion variations. The variation of annual runoff series with 8–12 years is mainly
shown in 1965–1975 and 1998–2008, showing regional oscillations of alternating abundance
Appl. Sci. 2022, 12, 647 6 of 11
and depletion. In main cycle variation of 15–20 years, annual rainfall mainly occurred in
1965–1975 and annual runoff mostly in 1965–1980. In main cycle variation of 15–20 years,
annual rainfall mainly occurred in 1965–1975 and annual runoff mostly in 1965–1980. The
variation of hydrological time series with a period of 25–30 years is relatively stable in
the whole–time domain and combined with the results of wavelet mode squared contours
and variograms, the first main cycle of rainfall and runoff has a consistent performance of
28 years.
In this paper, Figures 5 and 6 show the discrete wavelet transform is applied to the
monthly average runoff and rainfall information of 28-year main cycle, and the high frequency
and low frequency parts of the information are extracted respectively, then the runoff and
rainfall wavelet coefficients are used as the input layer to make predictions. In essence, the
main cycle runoff and rainfall information is used as a training set for the prediction of
wavelet coefficients in addition to the prediction of the original data. The advantage is that the
magnitude and structure of the weights of the high frequency and low frequency components
of the prediction model reflect the relationship between the constituent components of the
hydrological time series. Meanwhile, the periodic runoff and rainfall data are used as the
training set to improve the accuracy of predicting sediment load.
The sequence after discrete wavelet decomposition for data of rainfall and runoff
is shown from Figures 7–10, in which the detail signals (D1, D2 . . . D4, D5 and D6) are
high-frequency components and approximate signal (A) can be included in the trend term.
As seen in Figures 9 and 10, the fluctuations of D3, D4, D5 and D6 are slight, while the D1
Appl. Sci. 2022, 12, 647 7 of 11
and D2 as the random item of the original sequence with high volatility. The detail signals
(D1, D2 . . . D4, D5 and D6) and approximate signal (A) can be used to expressed for time
series data characteristics.
Figure
Figure 11.11. Observed
Observed and
and predicted
predicted values.
values.
ForFor predicting
predicting model
model effects,
effects, the the
RootRoot
MeanMean Squared
Squared ErrorError (RMSE),
(RMSE), the fitthe fit rating
rating (S)
and the Correlation Coefficient (r) are commonly used to determine the accuracy of the of
(S) and the Correlation Coefficient (r) are commonly used to determine the accuracy
the model’s
model’s predictions.
predictions. The RMSE
The RMSE represents
represents the distribution
the distribution of data of data around
points points around
the best-the
best-fit line. For a perfect match between the observed and predicted
fit line. For a perfect match between the observed and predicted values, the rootvalues, the root mean
mean
square error is zero. The closer the observed and predicted values are the
square error is zero. The closer the observed and predicted values are the smaller thesmaller the RMSE
value [27]. The RMSE can be calculated by using the following relationship:
RMSE value [27]. The RMSE can be calculated by using the following relationship:
s
𝛴 𝑄 Σ−( Q𝑄mi − Qci )2
RMSE = =
RMSE 𝑀 M
Appl. Sci. 2022, 12, 647 9 of 11
where Qmi is the measured sediment on month i and Qci is the predicted sediment on
month i; M is the total number of values.
The general form of the fit rating (S) proposed by Willmott [28] is expressing as:
Σ ( Q mi − Q ci )2
S = 1− 2
∑ Q mi − Q mi + Q ci − Q ci
where Qmi is the measured average sediment on month i; Qci is the predicted average
sediment on month i. If the S is close to 1, it indicates that the error between the predicted
and measured results is small.
The Correlation Coefficient (r) demonstrates the closeness between the observed and
predicted values:
2
∑ Q mi − Q mi Q ci − Q ci
r= q 2 q 2
∑ Q mi − Q mi ∑ Q ci − Q ci
Based on the prediction results, the Root Mean Squared Error (RMSE), the fit rating
(S) and the Correlation Coefficient (r) are 0.05, 0.98 and 0.99, respectively, indicating
that the prediction model can predict the overall trend and volatility of sediment load
more accurately.
4. Conclusions
Predicting sediment transport is essential for effective river basin planning and man-
agement. Among them, the simulation of river sediment transport is highly variable and
nonlinear in nature, hence the difficulty of the runoff–sediment production hydrologi-
cal process remains challenging in terms of sediment prediction. Additionally, sediment
load data are highly dependent on runoff, while runoff is directly related to rainfall in-
tensity. In this study, based on wavelet transform analysis and neural network principle,
the rainfall-runoff-sediment model was established, so as to improve the accuracy of
sediment prediction.
Using the wavelet transform, this article obtains the study of interannual cycle vari-
ation of rainfall and runoff in the study area, both of which have the first main cycle of
28 years, and the main period data were decomposed to obtain the detail signals (D1, D2 . . .
D4, D5 and D6) of high-frequency components and approximate signal (A) on time series.
Then, the monthly average wavelet decomposition of rainfall and runoff was coupled into
the neural network training set as the input layer and the monthly average sediment is
used as the output quantity for predicting the monthly average sediment quantity. The
Root Mean Squared Error (RMSE), the fit rating (S) and the Correlation Coefficient (r) were
used to compare the observed and predicted monthly average sediment transport values;
the results show that the RMSE, the S and the r were 0.05, 0.98 and 0.99, respectively. The
model establishes the relationship between rainfall, runoff and sediment in river basins,
which effectively improves the accuracy of the predicted sediment model and can enhance
the recognition of river systems, thus providing a reference for the establishment of small
information sediment prediction in river areas, and provide a reference basis for sediment
prediction and improve the utilization rate of hydropower energy.
Author Contributions: Conceptualization, Z.L.; Data curation, W.X.; Funding acquisition, Z.S.
and H.Z.; Methodology, Z.L. and H.D.; Project administration, Z.S., L.S. and H.Z.; Visualiza-
tion, J.L.; Writing-original draft, J.L. All authors have read and agreed to the published version
of the manuscript.
Funding: The article was supported by was supported by the National Natural Science Foundation
of China (Grant No. 91647209) and the Public Welfare Project of Zhoushan City (No. 2022C31046).
Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.
Appl. Sci. 2022, 12, 647 10 of 11
Data Availability Statement: The data presented in the present study are available on request from
the corresponding author.
Acknowledgments: Great thanks were owed to the Hydrology and Water Resources Bureau in the
Lincang city where data was monitored and those two anonymous referees for their invaluable
suggestions to improve this article.
Conflicts of Interest: The author declares no conflict of interest.
References
1. Han, Z.; Long, D.; Huang, Q.; Li, X.; Zhao, F.; Wang, J. Improving Reservoir Outflow Estimation for Ungauged Basins Using
Satellite Observations and a Hydrological Model. Water Resour. Res. 2020, 56, e2020WR027590. [CrossRef]
2. Bonaldo, D.; Benetazzo, A.; Sclavo, M.; Carniel, S. Modelling wave-driven sediment transport in a changing climate: A case study
for northern Adriatic Sea (Italy). Reg. Environ. Chang. 2015, 15, 45–55. [CrossRef]
3. Nezhad, A.; Akbari, G. A new approach for prediction of estuary flow-sediment variation at the downstream end of a reach.
In Proceedings of the 5th IAHR Symposium on River, Coastal and Estuarine Morphodynamics, Enschede, NL, USA, 17–21
September 2007; CRC Press: Boca Raton, FL, USA, 2019.
4. Chai, Y.; Zhu, B.; Yue, Y.; Yang, Y.; Li, S.; Ren, J.; Xiong, H.; Cui, X.; Yan, X.; Li, Y. Reasons for the homogenization of the seasonal
discharges in the Yangtze River. Hydrol. Res. 2020, 51, 470–483. [CrossRef]
5. Bonaldo, D.; Antonioli, F.; Archetti, R.; Bezzi, A.; Correggiari, A.; Davolio, S.; De Falco, G.; Fantini, M.; Fontolan, G.; Furlani,
S.; et al. Integrating multidisciplinary instruments for assessing coastal vulnerability to erosion and sea level rise: Lessons and
challenges from the Adriatic Sea, Italy. J. Coast. Conserv. 2019, 23, 19–37. [CrossRef]
6. Bajirao, T.S.; Kumar, P.; Kumar, M.; Elbeltagi, A.; Kuriqi, A. Superiority of Hybrid Soft Computing Models in Daily Suspended
Sediment Estimation in Highly Dynamic Rivers. Sustainability 2021, 13, 542. [CrossRef]
7. Idrees, M.B.; Lee, J.; Kim, D.; Kim, T. Complementary Modeling Approach for Estimating Sedimentation and Hydraulic Flushing
Parameters Using Artificial Neural Networks and RESCON2 Model. KSCE J. Civ. Eng. 2021, 25, 3766–3778. [CrossRef]
8. Kamruzzaman, M.; Shahriar, M.; Beecham, S. Assessment of Short Term Rainfall and Stream Flows in South Australia. Water
2014, 6, 3528–3544. [CrossRef]
9. Reisenbüchler, M.; Bui, M.D.; Rutschmann, P. Reservoir Sediment Management Using Artificial Neural Networks: A Case Study
of the Lower Section of the Alpine Saalach River. Water 2021, 13, 818. [CrossRef]
10. Khan, M.A.; Stamm, J.; Haider, S. Assessment of Soft Computing Techniques for the Prediction of Suspended Sediment Loads in
Rivers. Appl. Sci. 2021, 11, 8290. [CrossRef]
11. Ibrahim, K.S.M.H.; Huang, Y.F.; Ahmed, A.N.; Koo, C.H.; El-Shafie, A. A review of the hybrid artificial intelligence and
optimization modelling of hydrological streamflow forecasting. Alex. Eng. J. 2022, 61, 279–303. [CrossRef]
12. Campozano, L.; Mendoza, D.; Mosquera, G.; Palacio Baus, K.; Célleri, R.; Crespo, P. Wavelet analyses of neural networks based
river discharge decomposition. Hydrol. Process. 2020, 34, 2302–2312. [CrossRef]
13. French, M.N.; Krajewski, W.F.; Cuykendall, R.R. Rainfall forecasting in space and time using a neural network. J. Hydrol. 1992,
137, 1–31. [CrossRef]
14. Jain, S.; Das, A.; Srivastava, D. Application of ANN for Reservoir Inflow Prediction and Operation. J. Water Resour. Plan. Manag.
1999, 125. [CrossRef]
15. Fan, H.; Jiang, M.; Xu, L.; Zhu, H.; Cheng, J.; Jiang, J. Comparison of Long Short Term Memory Networks and the Hydrological
Model in Runoff Simulation. Water 2020, 12, 175. [CrossRef]
16. Ali, S.; Shahbaz, M. Streamflow forecasting by modeling the rainfall–streamflow relationship using artificial neural networks.
Modeling Earth Syst. Environ. 2020, 6, 1645–1656. [CrossRef]
17. Javan, K.; Lialestani, M.R.F.H.; Nejadhossein, M. A comparison of ANN and HSPF models for runoff simulation in Gharehsoo
River watershed, Iran. Modeling Earth Syst. Environ. 2015, 1, 56390. [CrossRef]
18. Wang, L.; Zhang, S.; Duan, Q. The Analysis of Wavelet Neural Network In Oilfield Production Prediction. Appl. Mech. Mater.
2013, 368–370, 1804–1807. [CrossRef]
19. Yang, Y.; Jia, X.; Wendroth, O.; Liu, B. Estimating Saturated Hydraulic Conductivity along a South-North Transect in the Loess
Plateau of China. Soil Sci. Soc. Am. J. 2018, 82, 1033–1045. [CrossRef]
20. Quilty, J.; Adamowski, J. A stochastic wavelet-based data-driven framework for forecasting uncertain multiscale hydrological
and water resources processes. Environ. Model. Softw. 2020, 130, 104718. [CrossRef]
21. Shi, W.; Yu, X.; Liao, W.; Wang, Y.; Jia, B. Spatial and temporal variability of daily precipitation concentration in the Lancang River
basin, China. J. Hydrol. 2013, 495, 197–207. [CrossRef]
22. Malekpour Heydari, S.; Aris, T.N.M.; Yaakob, R.; Hamdan, H. Data-Driven Forecasting and Modeling of Runoff Flow to Reduce
Flood Risk Using a Novel Hybrid Wavelet-Neural Network Based on Feature Extraction. Sustainability 2021, 13, 11537. [CrossRef]
23. Ian, G.; Bengio, Y.; Courville, A. Deep Learning; MIT Press: Cambridge, MA, USA, 2016.
24. Hecht-Nielsen, R. Theory of the backpropagation neural network. In Proceedings of the IEEE TAB Neural Network Committee,
San Diego, CA, USA, 24–27 July 1988.
25. Li, H.; Song, W. Characteristics of Climate Change in the Lancang-Mekong Sub-Region. Climate 2020, 8, 115. [CrossRef]
Appl. Sci. 2022, 12, 647 11 of 11
26. Irannezhad, M.; Liu, J.; Chen, D. Influential Climate Teleconnections for Spatiotemporal Precipitation Variability in the Lancang-
Mekong River Basin From 1952 to 2015. J. Geophys. Res. Atmos. 2020, 125, e2020JD033331. [CrossRef]
27. Kumar, M.; Kumari, A.; Kushwaha, D.P.; Kumar, P.; Malik, A.; Ali, R.; Kuriqi, A. Estimation of Daily Stage–Discharge Relationship
by Using Data-Driven Techniques of a Perennial River, India. Sustainability 2020, 12, 7877. [CrossRef]
28. Willmott, C.J. On the Validation of Models. Phys. Geogr. 1981, 2, 184–194. [CrossRef]