Simulation of a semiarid stream flow using the 1D model (Rubarbe) : case of Ksob Wadi in M. Hasbaia, H. Adoui, André Paquier

To cite this version:

M. Hasbaia, H. Adoui, André Paquier. Simulation of a semiarid stream flow using the 1D model (Rubarbe) : case of Ksob Wadi in Algeria. Online Journal of Space Communication, Ohio University, 2015, 25, pp.120-126. ￿10.1016/j.proenv.2015.04.017￿. ￿hal-01281106￿

HAL Id: hal-01281106 https://hal.archives-ouvertes.fr/hal-01281106 Submitted on 1 Mar 2016

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Available online at www.sciencedirect.com ScienceDirect

Procedia Environmental Sciences 25 ( 2015 ) 120 – 126

7th Groundwater Symposium of the International Association for Hydro-Environment Engineering and Research (IAHR) Simulation of semiarid stream flow using the 1D model (Rubarbe), case of Ksob wadi in Algeria

M. Hasbaiaa*, H. Adouia and A. Paquierb aUniversity of Msila, PO Box 166 Ichebilia 28000 Msila, Algeria bIrstea, UR HHLY,5 rue la Doua,CS 70077, 69626 Villeurbanne Lyon, France

Abstract

This study aims to show the ability of 1D model to simulate the wadi flow during the flood. The wadi is a semiarid stream characterized by the non uniform geometry and an intermittent flow. These complexities can be modelled by a robust 1D model such as Rubarbe. In this article we study the 1994 flood of Ksob wadi that caused a large inundation of M’sila town in the center of Algeria. The simulated reach located downstream the Ksob dam, crosses M’sila over a length of 6800 m with a bed slope of 0.6 %. The recorded flow peak of the 1994 flood is about 1430 m3/s and flood duration is 24 days. The flood propagation is simulated through the studied reach using the Rubarbe code (1D model) developed in Irstea (former Cemagref). The results show that the wadi overflows in some sections at the upstream and the center of the reach. The major part of the overflow is observed just around the peak of the flood, but the duration is variable from a few minutes to several hours. The overflow hazard limits are plotted in each section from the intersections of the highest water elevation recorded during the flood with the ground level of the two banks of the same section, taking into account the floodplain slope. According to our investigations, the calculated flooded area is very close to the real one, especially in the central part of the reach.

© 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license © 2015 The Authors. Published by Elsevier B.V. (http://creativecommons.org/licenses/by-nc-nd/4.0/). PeerPeer-review-review under under responsibility responsibility of the of Scientific the Scientific Committee Committee of the IAHR of theGroundwater IAHR Groundwater Symposium 2014 Symposium 2014.

Keywords: 1D flow model, Rubarbe, Ksob wadi, flood, inundation.

* Corresponding author. E-mail address: [email protected]

1878-0296 © 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/). Peer-review under responsibility of the Scientific Committee of the IAHR Groundwater Symposium 2014 doi: 10.1016/j.proenv.2015.04.017 M. Hasbaia et al. / Procedia Environmental Sciences 25 ( 2015 ) 120 – 126 121

1. Introduction

From the beginning of human life on the Earth, floods were one of the natural disasters that threaten safety and properties. Through the history, flood control was a big stake for all the civilizations. Nowadays, all the world countries - among which Algeria - are trying to use all the possible means in order to forecast and fight against these disasters. According the International Institute of Stockholm (SIWI), floods are causing about 20000 casualties per year [1].Mankind participated to the aggravation of flooding process by degrading the natural environment along most of the streams (impermeabilization of the catchment). The land impermeabilization by buildings and infrastructures is well known to increase the runoff coefficients and to accelerate the water flows and thus to aggravate the discharges received downstream. Then, the reduced capacity of sewage networks is considered as one of the causes of the flooding processes. Among the Mediterranean regions, Algeria is one that is faced to processes of flash floods above all in the semi-arid and arid areas. These floods that suddenly occur are often difficult to forecast. Flow is rising rapidly, area discharge is high, rainfall rate is high and the concerned basins have a rather small area. Several disasters caused by these types of floods were registered in Algeria (M’sila in 1994, the September 2000 flood in Boussaâda area, Ghardaïa and Bechar in October 2008, etc.). In this paper, the study is focussed on one of these events, the case of Ksob wadi. The overflow of the 1994 year flood is investigated starting from the definition of the flow hydrograph. Then, the flow hydrograph is propagated along the reach of Ksob wadi (6800 m long) that crosses the city using the RubarBE code (1-D model) that was developed at Irstea (France).

2. Presentation of the study area

The basin of an area of 26 000 km² is the fifth large basin of Algeria (Figure 1), It is located at about 150 km southwards from the Mediterranean coast (Bejaïa bay). The elevations of the summits of decrease from East to West and oscillate between 1900 and 1000 m while in the southern area, a few summits located in the Saharian Atlas reach 1200 m. The location of the Hodna basin between two series of mountains at North and South structures the basin around a nearly flat closed basin at an elevation of 400 m and that receives the surface flow of the region. At the centre of the basin, the chott El Hodnahas an area of 1 150 km².The sub basin of Ksob wadi is located in the northern part of Hodnaand covers an area of 1456 km2. It is bounded to the north and northwest by the mountain range , south and southwest by the mountains of Hodna and east by the high plains of Setif. The catchment of Ksob wadi consists of two separate physical units with mountainous borders, between which lies a vast plain called the high plains of BordjBouArreridj. It is characterized by a dense hydrographic network.This basin is characterized by a semi-arid climate, interannual average precipitation of 225 mm in 1993/1994 with runoff velocity1.52 km / h and the average temperature is 17 ° C.

122 M. Hasbaia et al. / Procedia Environmental Sciences 25 ( 2015 ) 120 – 126

Fig. 1. Location of k’sob basin

3. Hydrodynamic modelling

The studied reach is a section of the K’sob wadi situated downstream of Ksob dam, through the city of M'sila 6800 m long, equipped with a hydrometric station (Medjaz). The upstream end is near of Moulions Erriad « Hodna »and the downstream end Sidi Amara. The reach contains several meanders, the width varying from 34 to 230 m with low slope of 0, 6%. Material from the bottom is characterized by irregular diameter that varies between 0.001 to 42 mm. According to the hydraulic services, since the completion of the dam K'sob. Figure 2 represents the wadi location on a photo from Google Earth and location of the cross section used in RubarBE calculation. Figure 3 shows the longitudinal profile and an example of cross section from the dataused in modelling.

K’sobwadi

Fig. 2. Location of k’sobwadi

M. Hasbaia et al. / Procedia Environmental Sciences 25 ( 2015 ) 120 – 126 123

-2444

Fig. 3. Geometry of Ksob wadi.

4. Description of the code RubarBE

The details of the 1-D code are provided in [2], [3], [4] and the key features of the model are summed up in [5], [6]. The Barred Saint Venante quations (Eqs. 1 & 2) are solved by a Godunov type finite volume scheme issued of research studies [7] about the 1-D modelling of shock problems (snow avalanches, dam break waves, landslides in a reservoir, etc). ∂A ∂Q + = q (1) ∂t ∂x ∂Q ∂ ⎛ Q² ⎞ ∂z f QQ + ⎜β + P⎟ = − − gAgA + B (2) ∂ ∂xt ⎝ A ⎠ ∂x D² where A is the wetted section, Q the flow discharge, Zf the bed elevation, x the distance along the channel, D the conveyance, B the lateral pressure, P the pressure, β the velocity coefficient, g the gravity, t the time. Then, the code was modified to take into account sediment transport and change of topography during the floods. Consequently, it uses several sediment transport capacity equations and several head loss formulas.

5. Method and results

For simulated flood of 23/09/1994, the ANBT (National Agency for Dams and Transfers) provided flow discharges evacuated by the dam during the flood, this flood hydrograph (downstream from the dam) is shown in Figure 4. For the friction coefficient, Strickler formula (k = 21/d501/6) is used in which the d50 is determined by particle size analysis of each survey section based on engineering study carried out by TESCO (Hungarian engineering office) of which the results are summed up in Table 1

124 M. Hasbaia et al. / Procedia Environmental Sciences 25 ( 2015 ) 120 – 126

Fig. 4. Hydrographof the 1994 flood (measured flow at dam site)

Table 1. The median diameter.

Cross 1R 2L 7R 8L 14L 15R 18L 19R 25R 26L 28L 38R 44R 50L section

d50 [mm] 18.00 0,01 0,12 0,10 0,11 0,09 0,55 0,01 0,09 6.00 0,14 0,125 1,75 0,14

The flow propagation along the reach described here above using RubarBE code provides the water elevation , the velocity and the flow discharge in every section and every time step during the whole flood duration. Using these latter results, the overflow depth Hd is calculated by Equation 3.

Hd= Zw – min(ZL, ZR) (3)

where Zw is the water elevation, ZL and ZR the left and right bank elevations respectively. Using Surfer software, the overflow depth of the studied flood is drawn in every section (figure 5).

Fig. 5. Overflow depth along Boussaâda wadi (100 year flood). M. Hasbaia et al. / Procedia Environmental Sciences 25 ( 2015 ) 120 – 126 125

Figure 5 shows that K’sob wadi overflows in a few sections upstream and nearly all the sections in the downstream half. The overflow duration varies between one and 30 hours along the reach. The hazard map (Figure 6) represents the hydraulic constraint that concerns every section because of the hydrologic processes in the basin and hydraulic behaviour in the stream. This map synthetizes the knowledge of hazard estimated for a reference event from the available information. The hazard limits are obtained linking the crossings of the higher water elevations with the bed elevation in every section on the two banks of the concerned wadi and respecting the slope of the flood plain. The intersection the higher water line recorded during the flood with the natural terrain on a map topographic a scale 1/25000 allowed us to define the overflow boundaries reached during this flood. Then, overflow boundaries are displayed on a photo of Google Earth (figure 6).

Fig. 6. Hazard mapping of Ksob wadi (1994 year flood).

To check the quality of our results on field we conducted a survey to identify the actual limits recorded in 1994. The comparison of real and simulated limits shows that our results can be accepted in a forecasting framework. 126 M. Hasbaia et al. / Procedia Environmental Sciences 25 ( 2015 ) 120 – 126

The extrapolation of these results to other settings (layout ...) requires an extra effort to view data and modeling.

6. Conclusion

This study used the 1-D code RubarBE to simulate the 1994 year flood along a 6800 m long reach of K’sob wadi that crosses the city. Data of flow discharge provided by ANBT are inputs to the model. The calculation permit to draw the flood limits (on a topographic map) with suitable accuracy as checked on the filed by ground surveys. The results of this study can be improved coupling a 1-D calculation in the main channel and a 2-D calculation in the floodplain.

References

[1] Yahiaoui, A., Inondations torrentielles cartographie des zones vulnérables en Algérie du nord (cas de l’oued Mekrra, wilaya de Sidi BelAbbés), PhDthesis, Ecole Nationale Polytechnique, 2012. [2] Balayn, P., Contribution à la modélisation numérique de l’évolution morphologique des cours d’eau aménagés lors de crues, PhD thesis, Université Claude Bernard, Lyon, France, 2001. [3] El Kadi K., Evolution d’un lit de rivière en fonction des apports. PhD thesis, Université Claude Bernard, Lyon, France, 2006. [4] Paquier, A., 2013. Logiciels Rubar3- RubarBE - Notice d’utilisation. Irstea, Villeurbanne, France [5] El Kadi K, Paquier A, and Gay A., One-dimensional numerical modelling of dam-break waves over movable beds: application to experimental and field cases, Environmental Fluid Mechanics, 8 (2008):169-198. [6] El Kadi K, and Paquier A., One-dimensional numerical mode ling of sediment transport and bed deformation in open channels, Water Resources Research 45(2009):W05404. [7] Vila, J.P., Sur la théorie et l'approximation numérique de problèmes hyperboliques non linéaires. Applications aux équations de SAINT-VENANT et à la modélisation des avalanches de neige dense, PhD thesis, Université Paris VI, 1986.