Journal of
Engineering and Technology Research

  • Abbreviation: J. Eng. Technol. Res.
  • Language: English
  • ISSN: 2006-9790
  • DOI: 10.5897/JETR
  • Start Year: 2009
  • Published Articles: 195

Full Length Research Paper

PSO-ANN’s based suspended sediment concentration in Ksob basin, Algeria

Baazi Houria*, Kalla Mahdi and Tebbi Fatima Zohra
Natural Hazards and Territory Planning Laboratory (LRNAT), Hadj Lakhdar University, Batna (UHLB), Algeria.
Email: [email protected]

  •  Received: 19 October 2014
  •  Accepted: 24 November 2014
  •  Published: 28 December 2014


Suspended sediment concentration estimation has a major influence on river basin planning and management. Prediction of such parameter with artificial neural network (ANN) has shown its performance, because of this, Back Propagation Neural Network model trained with particle swarm optimization (PSO) is used to forecast the daily sediment concentration for Ksob river, Tebessa using 22 years data set from Morsott gauging station; the recorded daily suspended sediment concentrations and correspondent daily discharges were used to train the ANN model. PSO is used to allow ANN architecture to be easily optimized. Simulation of both ANN and PSO-ANN models has shown more accurate results compared with the traditional sediment rating curve.


Key words:  Ksob basin, Algeria, sediment rating curve, neural networks, particle swarm optimization.


The estimation of sediment yield is of crucial importance, especially to ensure the efficient operation of hydraulic structures management and also for environmental issues. For many years, suspended sediment concentrations have been collected in gauging stations and to estimate sediment concentrations during periods when samples are not collected, relations between sediment concentrations or load and stream flows are developed and known as sediment rating curve (SRC) (Asselman, 2000; Cohn et al., 1989; Horowitz, 2002). This technique has been use for a long time despite errors involved in using it (Walling, 1977). 
Artificial neural network (ANN) is an alternative and complementary set of techniques to traditional models (Abrahart et al., 2005; Cigizoglu, 2008). They may be treated as an universal approximators (ASCE, 2000)  and most of their attractive features are self-learning, self-adapting, good robustness and capability of dealing with non-linear problems (Dhar, 2010), without the physics being explicitly provided to them (Nagesh, 2004).
Among many successful bio-inspired Swarm Intelligence, the well known approach particle swarm optimization (PSO) which was inspired from the social behaviors of birds flocking or fish schooling and was developed by Kennedy and Eberhart (1995). It is similar to Genetic Algorithms (GA) in the sense that is a population-based search method.
Combining ANN and Swarm Intelligence has received much attention in several engineering issues in order to overcome single model deficiencies (Chen, 2014). In the present study, an hybrid model of a ANN combined with particle swarm optimizer (PSO) in order to improve ANN training for suspended sediment modelling of Morsott gauging station data in the Ksob basin of Tebessa town, Algeria.
Study area
The Ksob catchment (1304 km2) is located in the extreme northeast of Algeria (Figure 1). It is a sub –basin of Wadi Medjerda, and a part of the Sahara Atlas. It is wedged between  of  35° 21` to 35° 43 N and longitude 7° 43` to 8° 21` E. The region is characterized by semi-arid climate with very hot and dry summer and very cold and wet winter. Inter annual average rainfall is about 339 mm. These rains are often in stormy form. The annual average temperature is around 16°C with a maximum in July (25.87°C) and a minimum in January (6.52°C).
Vegetation consists of forests, maquis and reforestation, the remaining area is divided between agricultural lands and wastelands. The development of agriculture is heavily compromised by soil characteristics, runoff, and intensity of human action. The basin consists mainly of old and recent alluvium, clay, sand stone and limestone gravel. These formations gave good soil permeability.


Sediment rating curve
A sediment rating curve (SRC) is a relationship established between sediment  concentration,  C,  and  water  discharge, Q, or between sediment load Qs, and discharge Q. It is a ‘black box’ type model not directly related to any physical parameters (Tebbi et al., 2012)and having a standard form of:
Coefficients a, a’ and b, b’ are empirically determined where a represents the sediment concentration for a discharge of 1.0 m3s-1, a' represents the sediment load for a discharge of 1.0 m3s-1 and b and b’ reflect load response to changes in discharge and typically fall in the range of 1 to 2 or possibly higher (Julien, 2010).
Large number of studies  were conducted to estimate suspended sediment yields using SRC method in Algerian rivers (Achite and Ouillon, 2007; Khanchoul et al., 2010; Khanchoul and Jansson, 2008; Terfous et al., 2001; Touaibia et al., 2001). Most of these studies preferred the form of Equation (2) of the method (Qs=a’Qb’).
Artificial neural networks
Artificial neural networks (ANN) consist of layers of interconnected artificial neurons (Figure 2). A neuron (also called a “node” or “unit”) is the basic unit of an artificial neural network, simulating a biological neuron. It performs a weighted sum of inputs and passes this to an activation function to produce the output of the neuron. A multi-layer perceptron is feedforward neural network that consists of neurons arranged in a distinct layered topology. The input layer nodes serve to introduce the values of the input variables. The hidden and output layer neurons are each connected to all of the units in the preceding layer.
For most applications fully-connected networks are better. There is usually some weight associated with every connection. Input layer represents the non processed information that is introduced into the network. The activity of each hidden unit is determined by the activities of the input units and the weights on the connections between the input units and the hidden units. The behavior of the output units depends on the activity of the hidden units and the weights between the hidden units and output units (Sivanandam and Deepa, 2006). To adjust the weights, the most common learning algorithm is called Back Propagation.
ANN are being used alone or hybridized  increasingly to predict and forecast water resources variables and  have seen significant activity in various hydrology related areas such as rainfall-runoff (Kim, 2005), stream flow (Dolling and Varas, 2002), ground water (Daliakopoulos et al., 2005; Sreekanth et al., 2009), precipitation forecasting and water quality issues (Eslamian and Lavaei, 2009; Shamim et al., 2004), and for designing optimal strategies in reservoir operation (Chandramouli and Deka, 2005; Chaves and Chang, 2008; Hasebe and Nagayama, 2002). Also, ANN have been widely used for sediment transport  modeling (Abrahart et al., 2008; Adib and Tagavifar, 2010; Chutachindakate, 2009; Cigizoglu and Alp, 2006; Cigizoglu and Kisi, 2006; Firat and Güngör, 2010; Kisi, 2007; Liu et al., 2013; Rajaee et al., 2011; Wang et al., 2008).
In Algeria, Boukhrissa et al. (2013) investigated the ability of ANN  models to improve the accuracy of streamflow-suspended sediment relationships in daily and annual suspended sediment estimation for the El Kebir catchment and shown that the ANN models have the highest efficiency to reproduce the daily sediment load and the global annual sediment yields.
Kisi (2012) using daily stream-flow and suspended sediment concentration data from two stations on the Eel River in California, found that Least Square Support Vector Machine performs better than ANN models. Kisi et al. (2012a) applied Genetic Programming technique for estimating the daily suspended sediment load in two stations in Cumberland River in U.S.; results were comparecd with those of the Adaptive Neuro-Fuzzy Inference System, Artificial Neural Networks and Support Vector Machine, Genetic Programming give beter accuracy. Most of actual studies attempt to improve ANN models accuracy by hybridization with several evolutionary  techniques.
Ramezani et al. (2014) used an ANN approach to predict sediment for Maroon River in Iran by optimizing the ANN connection weights with social based algorithm (SBA). Kisi et al. (2012b) compared neural networks with Artificial Bee Colony algorithm model with those of the Neural Differential Evolution, Adaptive Neuro-Fuzzy, Neural Networks and Rating Curve models using data from two stations, Rio Valenciano Station and Quebrada Blanca Station results showed that the ANN-ABC was able to produce better results than other models.
Kisi (2005) found that neuro-fuzzy gives better estimates than neural networks. Using 11-year data (1994 to 2004) of  Doiraj River located in Iran, Kalteh (2013)  applied both ANN and SVM models for predicting suspended sediment load and found that ANN models based on  a quasi-Newton method named Broyden, Fletcher, Goldfarb and Shanno (BFGS) algorithm has a better performance comparing with other models.
Particle swarm optimization
The basic particle swarm optimization (PSO) algorithm consists of three steps, namely, generating particle’s positions and velocities, velocity update, and finally, position update.  First, the positions, Xid(k) and velocities, Vid(k) of the initial swarm of particles are randomly generated.
The positions and velocities are given in a vector format for the ith particle at time k. The second step is to update velocities of all particles at time k +1 using the particles objective or fitness values, which are, function of the particles current positions in the design space at time k.
The fitness function value of a particle determines which particle has the best global value in the current swarm, , and determines the best position of each particle over time, pid, that is, in current and all previous moves.
The velocity update formula uses these two pieces of information for each particle in the swarm along with the effect of current motion, vid(k), to provide a search direction, Vid(k+1), for the next iteration.
Where i=1,2,…n . n is the number of particles. d=1,2,3,…m, m is the number of input variables to be optimized.W=weights trading off the impact of the local best and global best solutions’ on the particle’s total velocity. C1,C2 Position update is the last step in an iteration. φ1, φsamples a uniform random distribution [0,1].
The three steps of velocity update, position update, and fitness calculations are repeated until a desired convergence criterion is met.


Figure 3 is constructed using 294 data of both mean daily suspended sediment concentrations and mean daily discharges of Morsott gauging station in the period of 02/1975 to 10/2002. Goodness of sediment rating curve fit using Equation (1) is checked using coefficient of determination R². R² obtained for the best fit by a single curve is 0.527 and corresponding root mean square error is 6.0359 g l-1.
ANN model
To be able to compare SRC method and other methods we used daily discharges (m3/s) as  the  unique  inputs  in ANN_BP model. Matlab’s neural network fitting tool creates a two-layer feed-forward with sigmoid hidden neurons and linear output neurons and it is trained with Levenberg-Marquardt backpropagation algorithm. Also training, validation and test sets are sampled respectively taking 70, 15 and 15%. Number of neurons in hidden layer is taken as 3. Figure 4 compares outputs of the ANN_BP model and observed suspended sediment concentrations.
PSO- ANNmodel
For combining particle swarm optimization algorithm and BP_ANN algorithm using a Matlab® code, the first step is initialization of random particles that represents a set of weights, and then simulating neural network, then evaluate fitness value of initialized particles. Determination of optimal solution individual Pbest and global Gbest according to the fitness value of each particle, and updating the velocity and position of particles according to Formulas (3) and (2). Best particle of current particles is stored. Repeat iteration until maximum number of iteration is met (Figure 5).
Estimating suspended sediment concentration by SRC, ANN and PSO_ANN’s respectively enable us to estimate suspended sediment concentrations in the Ksob basin. Figure 6 shows linear regression between observed and simulated concentration using different transfer functions and different training algorithms of PSO_ANN model for test period. We can observe that R² is improved about.
Table 1 indicates models performances using R square coefficient and root mean square error and shows that using hybrid ANN_BP model with principal transfer functions in hidden layer and training   algorithms  at  test period for same sample improves RMSE and R²; results show that PSO-ANN’s improves ANN’s model at least RMSE 21% and R² 18% and both ANN and PSO-ANN model are better than SRC model.


Sediment rating curve, neural networks back propagation and hybrid PSO based neural networks models were employed in this study in order to model suspended sediment concentrations in Ksob river using Morsott gauging station data. Compared with the ANN model, PSO-ANN’ appears to be more suitable for suspended sediment concentration modeling. Due to the lack of data, multiple linear regression and recurrent ANN’s models were not considered, the proposed model based on PSO can be further investigated with several improved PSO models also other models in future studies.


The authors would like to acknowledge the Algerian Agency of Water Resources (ANRH, Algiers) for providing us the required data.


The authors have not declared any conflict of interest.


Abrahart RJ, Kneale PE, See LM (2005). Neural networks for hydrological modelling. Taylor & Francis, LONDON, UK.
Abrahart RJ, See LM, Heppenstall AJ, White SM (2008). Neural Network Estimation of Suspended Sediment: Potential Pitfalls and Future Directions. In: Abrahart RJ, See LM, Solomatine DP (eds) Practical Hydroinformatics, vol. 68. Water Science and Technology Library. Springer Berlin Heidelberg, pp. 139-161.
Achite M, Ouillon S (2007). Suspended sediment transport in a semiarid watershed, Wadi Abd, Algeria (1973-1995). J. Hydrol. 343:187-202.
Adib A, Tagavifar A (2010). Evaluation and Comparison Different Methods of Preparation of Sediment Rating Curve in Telezang Station of the Dez River. Austr. J. Basic Appl. Sci. 4(5):717-723
ASCE TC (2000) Artificial Neural Networks in Hydrology II. Hydrologic Applications. J. Hydrol. Eng. 5:124-137.
Asselman NEM (2000). Fitting and interpretation of sediment rating curves. J. Hydrol. 234:228-248.
Boukhrissa ZA, Khanchoul K, Le Bissonnais Y, Tourki M (2013). Prediction of sediment load by sediment rating curve and neural network (ANN) in El Kebir catchment, Algeria. J. Earth Syst. Sci. 122:1303-1312.
Chandramouli V, Deka P (2005) Neural Network Based Decision Support Model for Optimal Reservoir Operation. Water Resour. Manag. 19:447-464
Chaves P, Chang FJ (2008). Intelligent reservoir operation system based on evolving artificial neural networks. Elsevier, Kidlington, ROYAUME-UNI. 31(6):926-936.
Chen ZY (2014). A Hybrid Algorithm by Combining Swarm Intelligence Methods and Neural Network for Gold Price Prediction. In: Wang LL, June J, Lee CH, Okuhara K, Yang HC (eds) Multidisciplinary Social Networks Research, 473. Communications in Computer and Information Science. Springer Berlin Heidelberg, pp. 404-416.
Chutachindakate C (2009). Integrated sediment approach and impacts of climate change on reservoir sedimentation. Kyoto University.
Cigizoglu HK (2008). Artificial Neural Networks In Water Resources. In. pp.115-148.
Cigizoglu HK, Alp M (2006). Generalized regression neural network in modelling river sediment yield Adv. Eng. Softw. 37:63-68. 
Cigizoglu HK, Kisi O (2006). Methods to improve the neural network performance in suspended sediment estimation. J. Hydrol. 317:221-238
Cohn TA, Delong LL, Gilroy EJ, Hirsch RM, Wells DK (1989). Estimating constituent loads. Water Resour. Res. 25:937-942.
Daliakopoulos IN, Coulibaly P, Tsanis IK (2005). Groundwater level forecasting using artificial neural networks. J. Hydrol. 309(1-4):229-240.
Dhar V (2010). Comparative performance of some popular artifcial neural network algorithms on benchmark and function approximation problems PRAMANA. J. Phys. 74(2):307-324.
Dolling OR, Varas EA (2002). Artificial neural networks for stream flow prediction. J. Hydraul. Res. 40(5):547-554.
Eslamian S, Lavaei N (2009). Modelling nitrate pollution of groundwater using artificial neural network and genetic algorithm in an arid zone. Int. J. Water. 5(2):194-203.
Firat M, Güngör M (2010). Monthly total sediment forecasting using adaptive neuro fuzzy inference system. Stoch. Environ. Res. Risk Assess. 24:259-270.
Hasebe M, Nagayama Y (2002). Reservoir operation using the neural network and fuzzy systems for dam control and operation support Adv. Eng. Softw. 33:245-260 doi:
Horowitz AJ (2002).The use of rating (Transport) curves to predict suspended sediment concentration: Am matter of temporal resolution. Paper presented at the Turbidity and Other Sediment Surrogates Workshop, Reno, NV.
Julien PY (2010). Erosion and Sedimentation (2nd Edition). Cambridge University Press.
Kalteh AM (2013). Monthly river flow forecasting using artificial neural network and support vector regression models coupled with wavelet transform. Comput. Geosci. 54:1-8
Kennedy J, Eberhart R (1995). Particle swarm optimization. Neural Networks, 1995.Proceedings IEEE Int. Conf. 4:942-1948. 68.
Khanchoul K, El Abidine Boukhrissa Z, Acidi A, Altschul R (2010). Estimation of suspended sediment transport in the Kebir drainage basin, Algeria Quaternary International In Press,
Khanchoul K, Jansson MB (2008). Sediment Rating Curves Developed On Stage And Seasonal Means In Discharge Classes For The Mellah Wadi, Algeria Geografiska Annaler: Series A, Phys. Geogr. 90:227-236. doi:10.1111/j.1468-0459.2008.341.x.
Kim D-IJaY-O (2005). Rainfall-runoff models using artificial neural networks for ensemble stream flow prediction Hydrol. Processes 19:3819-3835
Kisi O (2005) Suspended sediment estimation using neuro-fuzzy and neural network approaches. Hydrol. Sci. J. 50(4):683-696
Kisi O (2012). Modeling discharge-suspended sediment relationship using least square support vector machine J. Hydrol. 456-457:110-120.
Kisi Ö (2007). Development of Streamflow-Suspended Sediment Rating Curve Using a Range Dependent Neural Network. Int. J. Sci. Technol. 2(1):49-61.
Kisi O, Dailr AH, Cimen M, Shiri J (2012a). Suspended sediment modeling using genetic programming and soft computing techniques. J. Hydrol. 450-451:48-58.
Kisi O, Ozkan C, Akay B (2012b). Modeling discharge–sediment relationship using neural networks with artificial bee colony algorithm. J. Hydrol. 428–429:94-103.
Liu QJ, Shi ZH, Fang NF, Zhu HD, Ai L (2013). Modeling the daily suspended sediment concentration in a hyperconcentrated river on the Loess Plateau, China, using the Wavelet–ANN approach. Geomorphology 186:181-19.
Nagesh K (2004). ANN Applications in Hydrology - Merits and Demerits, Integrated Water Resources Planning and Management. Jain Brothers, New Delhi.
Rajaee T, Nourani V, Zounemat-Kermani M, Kisi O (2011). River Suspended Sediment Load Prediction: Application of ANN and Wavelet Conjunction Model. J. Hydrol. Eng. 16(8):613-627. doi:doi:10.1061/(ASCE)HE.1943-5584.0000347.
Ramezani F, Nikoo M, Nikoo M (2014). Artificial neural network weights optimization based on social-based algorithm to realize sediment over the river. Soft Comput. pp.1-13 doi:10.1007/s00500-014-1258-0.
Shamim MA, Ghumman AR, Ghani U (2004). Forecasting Groundwater Contamination Using Artificial Neural Networks. Paper presented at the International Conference on Water Resources and Arid Environment, King Saud University, Riyadh, Saudi Arabia.
Sivanandam SN, Deepa SN (2006). Introduction to Neural Networks Using Matlab 6.0. Tata McGraw-Hill.
Sreekanth PD, Geethanjali N, Sreedevi PD, Ahmed S, Kumar NR, Jayanthi PDK (2009). Forecasting groundwater level using artificial neural networks. Curr. Sci. 96(7):933-939.
Tebbi FZ, Dridi H, Morris GL (2012). Optimization of cumulative trapped sediment curve for an arid zone reservoir: Foum El Kherza (Biskra, Algeria). Hydrol. Sci. J. 57:1368-1377.
Terfous A, Megnounif A, Bouanani A (2001). Etude du transport solide en suspension dans l'Oued Mouilah (Nord Ouest Algérien). Rev. Sci. Eau. 14:173-185.
Touaibia B, Aidaoui A, Gomer D, Achite M (2001). Temporal quantification and variability of sediment discharge in a semiarid area in northern Algeria. Hydrol. Sci. J. 46:41-53.
Walling DE (1977). Limitations of the rating curve technique for estimating suspended sediment loads, with particular reference to British rivers. Paper presented at the Erosion and Solid Matter Transport in Inland Waters, Paris.
Wang Y-M, Traore S, Kerh T (2008). Using artificial neural networks for modeling suspended sediment concentration. Paper presented at the Proceedings of the 10th WSEAS International Conference on Mathematical Methods and Computational Techniques in Electrical Engineering, Sofia, Bulgaria.