Al-Bayan Center for Planning and Studies

SIMULATION RIVER FLOOD WAVE IN CITY DUE TO A HYPOTHETICAL MOSUL BREAK

Thair M. Al-Taiee _ and Anass M. M. Rasheed About

Al-Bayan Center for Planning and Studies is an independent, nonprofit think tank based in , . Its primary mission is to offer an authentic perspective on public policy issues related to Iraq and the neighboring region. Al-Bayan pursues its vision by conducting autonomous analysis, as well as proposing workable solutions for complex issues that concern academia and policymakers.

The opinions expressed in this paper do not necessarily represent the views of the center

www.bayancenter.org [email protected] PO Box 2268 Jadiriya, Baghdad SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK

SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK

1 2 Thair M. Al-Taiee and Anass M. M. Rasheed

ABSTRACT Mosul dam is one of the biggest hydraulic structures in Iraq. It is located on Tigris River north Iraq and was selected as a case study to predict flood disasters caused by a hypothetical Mosul dam failure due to its foundation defect in which the dam had been survive since 1986 (initial time of dam operation). It is an earth-fill zoned dam with 105 meter height and 11.1 billion cubic meter of volume capacity. The objectives of the present research were to predict the characteristics of the flood wave due to a hypothetical Mosul dam break and to estimate damaged areas downstream, specially in Mosul city (the second largest city in Iraq), applying the simplified dam break model (SMPDBK) and utilizing the geographic information system. Many scenarios of dam operation were utilized to simulate dam failure at various degrees of severity. In case of the dam failure of high severity (330 m.a.s.l.) Mosul city locating 60 km from Mosul dam will be inundated with the highest flooding level of 235.2 meter above sea level, the maximum flood discharge will be (207632) 3m /sec with an average velocity of 3.5 m/sec, while the highest water level will reach (25.3) meter above the natural Tigris River bed within (9) hours from the beginning of the dam failure. The flooded area between the dam site and Mosul city was 252km2, while the maximum percentage of the inundated area within Mosul city will be about (54%) of the whole area of the city. The output of the channel routing process between Mosul dam site and Mosul city will be the attenuation of the peak discharge during the dam failure by a percentages ranging from (43 to 38)%. Finally, the main constructions which will be inundated in Mosul city were recognized. Relevant organizations in public sectors can utilize this information and prepare strategic plans to reduce the impact of life loss and assets of government and local people in the downstream areas of the dam.

1. Assist. Prof., Research Center for and Water Resources, Mosul University, Iraq E-mail: [email protected] 2. Lecturer, Research Center for Dams and Water Resources, Mosul University, Iraq 3 Al-Bayan Center for Planning and Studies

Keywords: Inundation, Mosul dam, Flood, Dam Failure.

INTRODUCTION Dam is a ubiquitously disputed subject ranging from its benefits, the effects of dam construction towards the ecology system, the collapse of local communities, and epidemic diseases. In addition, the failure of large dam that might immediately occur could lead to tremendous loss of lives and assets. Risk assessment studies consider the failure of embankment dams often make use of breach parameter prediction methods that have been developed from analysis of historic dam failures. Similarly, prediction of peak breach outflow can also be using relations developed from case study data. Dam failure can be caused by overtopping a dam due to un-sufficient spillways capacity during large inflow to the reservoir, by seepage or piping through the dam or along internal conduits, slope embankment slides, earthquake damage and liquidation of earthen dams from earthquakes or landslide generated waves within reservoir. When a dam is breached catastrophic flash flooding occurs as the impounded water escape through the gap into downstream channel. The outflow flood hydrograph from a dam failure depends on many factors. The primary factors are the physical characteristics of the dam, the volume of reservoir and mode of failure. The parameters which control magnitude of the peak discharge and the slope of hydrograph include; the breach dimensions; the manner and length of time for breach to develop; the depth and volume of water stored in reservoir and inflow to reservoir at time of failure. A study conducted by Middlebrooks (1952) into causes of 220 earth dam failure during the period 1850-1950 summarizes observed caused agent and their frequency of occurrence prior to 1980s, little analytical work had been done on numerically describing the mechanics of failure. The extent of an affected downstream area by a flood wave resulting from a theoretical dam breach is a function of both the height of the flood wave and the downstream distance and width of the river at a particular location. An associated and important factor is the flood wave travel time (Attalla, 2002). The flood wave should be routed downstream the point where no the effect of the failure will be longer have negative consequences. There are numerous tools available today for analyzing dams’ failures and their

4 SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK resulting outflow hydrographs. Some of the best known and most widely used are the National weather service (NWS) Dam Break Flood forecasting model (DAMBRK); the U.S. Army corps of Engineer Hydrologic Engineering Center Flood Hydrograph package, HEC-1 (Hydrologic Engineering Center, 1981); and the NWS simplified Dam-Break Flood forecasting model, SMPDBRK (Wetmore and Fread 1983). One of these models is DAMBRK which is the most widely used. Wurbs (1987) discussed and compared the state of the art models available at that time and recommended the use of DAMBRK or SMPDBK depending on the level of accuracy required and resources and input data available. Dam break flood routing models e.g. DAMBRK (Fread, 1984) and FLDWAV (Fread, 1993) simulate the outflow from a reservoir and through the downstream valley resulting from a developing breach in a dam. These models focus their computational effort on the routing of the breach outflow hydrograph. The development of the breach is idealized as a parametric process, defined by the shape of the breach, its final size and the time required for its development (often called failure time). Thompson (1987) also compared DAMBRK and SMPDBK and recommended the use of SMPDBK as a screening tool and DAMBRK for more detailed analysis. Babb and Mermel (1988) summarized over 600 dam incidents throughout the world but high quality detailed information was lacking in most cases. During the 1980s several authors compiled data bases of well documented case studies for breach parameters or peak breach outflows (SCS, 1981; Singh and Snorrason 1982; MacDonald and Langridge_Monopolis, 1984; Costa, 1985; Froehlich, 1987, 1995b; and Singh and Scarlators, 1988).

STATEMENT OF THE PROBLEM As being Mosul dam is one of the most significant strategic projects in Iraq since any defect or failure that may happen in the dam God forbid will be a catastrophic effects on the population living downstream of this dam specially on Mosul city located on Tigris River 60 km downstream the dam site and having 1.7 million of population, Figure (1).

5 Al-Bayan Center for Planning and Studies

Figure (1) Location Map of the study area

The basic problem in this dam is in the continuous deterioration in the dam foundations which contain gypsum and anhydrite formations that dissolve under the effect of storing water in the dam reservoir head resulting in some phenomena such as fractures and infiltration that need a continuous treatment to reinforce this foundation. Finishing the building of Mosul dam and being into action in 1986, the problem of gypsum dissolving under the dam foundations came out and since the treatment has begun through grouting and cements injections to secure the dam safely which has been going on until now. As a result of this problem number of worrying geological phenomena such as infiltration in the left side, slides beside the right levee of the dam head, cavities in the foundation depth have been occurred and a fractures are appeared very close to the left side of the dam at distance 60 meter from the spillway with the increasing of infiltration quantities from the springs in this area (WWW.The Mosul dams press statement.mht, 2007). Therefore it is an important and mandatory to predict the flood caused by dam failure due to these reasons. This task can be accomplished by developing a model of dam collapse such as SMPDBK model and considering volume of water storage, water level, velocity and time taken for flood wave from dam to downstream area. A map

6 SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK demonstrating the flooded area and estimating the loss from such incident can then be produced. This information can be useful in planning, disasters announcement measures and evacuation of people living downstream area.

THE AIM OF THE RESEARCH For all the above reasons, the aim of this research work is to simulate the flood wave in Tigris River downstream Mosul dam as a result of a hypothetical failure describing the inundation area focusing on Mosul city predicting the maximum depth covering the city by preparing maps delineating these flooded areas along the Tigris River and indicating the travel time of flood wave on every significant habitation along the river reach. A computer model called simplified dam break model (SMPDBK) will be used for such problem. This model is adopted due to its simplicity and it allows for the investigation of partial failure occurring over a finite interval of time. The basic equations utilized by the model is the Saint Venant equations of unsteady flow which are the continuity and momentum equations which in-turn based on conservation of mass and momentum equations. Such research work will be useful in accomplishing the evacuation plans that should be taken into consideration by the responsible authority in case if there is a probability of dam collapse to eliminate the hazards and the damages that may occur in the life and the assets downstream the dam.

Description of Mosul Dam Mosul dam which has been built since 1985 is 100 meter height and of zoned earth filling type having a mud central core. It is a multipurpose project, for flood control and to store water with a volume of 11.11 billion cubic meters at the maximum operation level of the dam for the generating electric power and exploiting the reservoir behind for irrigation and tourist purposes. It is situated on the Tigris River at a distance 60 km north of Mosul city with a 1.7 million population. Figure (2) shows a dam cross-section with the operation levels.

7

Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt 287

Al-Bayana distance Center for 60 Planning km north and of Studies Mosul city with a 1.7 million population. Figure (2) shows a dam cross-section with the operation levels.

341m Crest level 338 m W. L. (PMF) 335 m (1000 Year flood)

330 m Max. operation level Earth fill shell (11.1 * 109 m3 storage) 8.16 * 109 m3 usable active storage for irrigation and 300 m Min. operation level power. (2.95 * 106 m3 storage)

245 m Bed level

Clay core

Figure (2): Schematic diagram of Mosul dam cross section. Figure (2): Schematic diagram of Mosul dam cross section.

Predicting Breach Parameters Predicting Breach Parameters To carry out a dam flood routing simulation, breach parameters were estimated and To carryprovided out a as dam input flood to the routing dam break simulation, and flood breach routing parameters simulation model.were estimated Several and providedmethods as areinput available to the fordam estimating break and breach flood parameters. routing All simulation of these methods model. are Several methodsbased areon regression available analysis for estimating of data collected breach from parameters. actual dam failures. All of An these equat ions methods are predictedbased on by regression Von Thum analysis and Gillette of data (1990), collected Froehlich from (1995a) actual and dam Bureau failures. of An reclamation (1988) and MacDonald and Langridge Monopolis (1984) and Singh and equationsSnorrason predicted (1984) were by selected Von Thum from the and literature Gillette and (1990), applied to Froehlich estimate the (1995a) Breach and Bureauwidth of and reclamation failure time (1988) or time andfor the MacDonald breach to be and developed Langridge completely Monopolis (Table 1).(1984) m Bed level 245 and SinghThese equations and Snorrason were depended (1984) on were the analysis selected of fromlarge casethe studyliterature actual and dams applied had to been failed due to many reasons, the equations used a basic dam dataClay concerning core estimateheight the of waterBreach in thewidth reservoir, and failurestorage volumetime or in time the reservoirs for the breachand height to ofbe breach. developed completelyAll of those (Table selected 1). These equations equations was examined were and depended calibrated on in the many analysis failed dams of large case confirmingstudy actual their dams feasibility had been to apply failed except due MacDonald to many reasons, and Langridge the equations Monopolis used a basicequations dam whichdata concerninggave large ove heightrestimated of watervalues fromin the the reservoir, others; therefore storage it was volume not in depended in the present work. the reservoirs and height of breach. All of those selected equations was examined and Table calibrated 2 shows in the many present failed recommended dams confirming values of the breach their width feasibility and time toof failure apply except MacDonaldfor five scenarios and Langridge of water levelMonopolis of reservoirs equations that was which fed to gavethe SMPBRK large overestimated model as input data for the simulation process. The water level in the reservoir represents the valuesnormal from or the sunny others; day loading therefore condition it was (i.e. not the dependeddam during innormal the prpoolesent failure); work. this is Tablebecause 2 shows the presentthe present hypothetical recommended dam failure values was assumed of the due breach to piping width or seepagand timee of failure in the dam foundation as the real condition now in Mosul dam. failure for five scenarios of water level of reservoirs that was fed to the SMPBRK model as input data for the simulation process. The water level in the reservoir represents the normal or sunny day loading condition (i.e. the dam during normal pool failure); this is because the present hypothetical dam failure was assumed due to piping or seepage failure in the dam foundation as the real condition now in Mosul dam.

8

288 Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK 288 Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt

Table (1) Breach width and failure time Table (1) Breach widthTable and (1) failure Breach timewidth and failure time Water level in Water height Breach width Failure time Author Waterreservoir level in Waterin reservoir height Breach(m) width Failure(hr) time Author reservoir(m.a.s.l) in reservoir(m) (m) (hr) (m.a.s.l)290 (m)40 155 1.05 300290 5040 189155 1.051.25 Von Thum 310300 6050 205189 1.251.45 (1990)Von Thum (1990) 320310 7060 230205 1.451.65 330320 8070 255230 1.651.85 290330 4080 388255 1.853.55 300290 5040 431388 3.554.23 Froehlich 300 50 431 4.23 (1995a)Froehlich 310 60 533 6.00 (1995a) 320310 7060 612533 6.007.57 330320 8070 694612 7.579.30 290330 4080 120694 9.301.32 Bureau of 300290 5040 150120 1.321.65 ReclamationBureau of 310300 6050 180150 1.651.98 (1988)Reclamation 320310 7060 210180 1.982.31 (1988) 330320 8070 240210 2.312.64 290330 4080 175240 2.640.25 Singh and 300290 5040 210175 0.250.44 SnorrasoSingh and 310300 6050 245210 0.440.62 nSnorraso (1984) 320310 7060 280245 0.620.82 n (1984) 330320 8070 315280 0.821.00 330 80 315 1.00

Table (2) RecommendedTable (2) Recommended breach width breach and width failure and time failure time Table (2) Recommended breach width and failure time Water Storag Surface Water Recommended Recommended Scenario levelWater in eStorag areaSurface of heightWater RecommendedBreach width RecommendedFailure time numberScenario thelevel dam in volumee reservoirarea of (m)height 9 6 2 Breach(m) width Failure(hr) time number (m.a.s.l.)the dam (10volume) reservoir(10 ) m (m) (109) 6 2 (m) (hr) 1 (m.a.s.l.)290 40 1.8m3 (10120) m 209 1.53 21 300290 5040 2.51.8m3 170120 243209 1.531.85 32 310300 6050 4.852.5 230170 291243 1.852.44 43 320310 7060 4.857.5 300230 333291 2.442.99 54 330320 8070 11.17.5 385300 376333 2.993.56 5 330 80 11.1 385 376 3.56

Concerning breach shape, for earth dams it was described as varying from triangular to trapezoidalConcerning as breach the breach shape, processed. for earth dams The greatit was majority described in asthe varying dam breach from triangulardescribed toas trapezoidal as the breach processed. The great majority in the dam breach described as trapezoidal in the literature. Froehlich (1995a) suggested a value for the side slope9 of breachtrapezoidal assuming in the to literature. be 1.4 for Froehlich overtopping (1995a) failure suggested or 0.9 for a valueother failurefor the modesside slope noting of thatbreach the assumingaverage side to be slope 1.4 factorfor overtopping for the 63 failurecase studies or 0.9 of for failed other dams failure was modes nearly noting 1. In thethat pre thesent average simulation side slope process factor using for the SMPDBK 63 case studiesmodel, ofseepage failed damsor piping was nearlyin the 1.dam In bodythe pre wassent proposed simulation to be process the reason using of SMPDBK dam failure. model, Concerning seepage breach or piping shape, in althoughthe dam body was proposed to be the reason of dam failure. Concerning breach shape, although Al-Bayan Center for Planning and Studies

Concerning breach shape, for earth dams it was described as varying from triangular to trapezoidal as the breach processed. The great majority in the dam breach described as trapezoidal in the literature. Froehlich (1995a) suggested a value for the side slope of breach assuming to be 1.4 for overtopping failure or 0.9 for other failure modes noting that the average side slope factor for the 63 case studies of failed dams was nearly 1. In the present simulation process using SMPDBK model, seepage or piping in the dam body was proposed to be the reason of dam failure. Concerning breach shape, although the model assumes rectangular shaped breach, a trapezoidal breach was analyzed by specifying a rectangular breach width that is equal to the average width of the trapezoidal breach. The river channel from tail water cross section downstream to each routing point is approximated as prismatic channel. Cross sections at each roting point are obtained from a topographic maps or satellite images for the area. Concerning manning roughness coefficient which used in the model, it was recommended to use n = 0.03 for the river channel. This value was calculated and calibrated by applying manning equation for a real cross-section near Mosul hydrological station for many measured discharges had been recorded in the river. For the river bank 0.04 was recommended for the manning roughness coefficient, while 0.1 was suggested within the Mosul city due to the existed building and many other obstructions.

The Simulation Process Several simulation runs were made using the simplified dam break model (SMPDBK) which was constructed by Fread (1983). The model requires input data related to reservoir bathymetry, dam geometry, embankment material and topographic maps for the downstream river reach with an estimation for the manning roughness coefficient along the river reach including the river channel, the flood plain along the river and also the obstruction existed within Mosul city such as buildings, etc. The SMPDBK model was fed with the required input data for the five scenarios of dam failure as mentioned in Table 2. The input data

10 SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK includes the calculated recommended breach width and failure time, surface area of reservoir, storage volume, maximum total discharge from the dam, bottom outlet elevation of the breach and manning roughness coefficient of the river reach. As a result of model operation, the output results of the model for the five scenarios of Mosul dam failure (i.e. 330, 320, 310, 300, 290) m.a.s.l in which include the maximum flow, maximum elevation, maximum depth and time to reach maximum depth for the flood wave along the distance downstream Mosul dam site were drawn along the distance downstream the dam site as shown in Figures (3,4,5,6). It is clear that the percentages of the maximum discharge attenuation between the first cross-section of the rive reach (0 km) downstream dam site and Mosul city cross-section for the five scenarios was (45%, 43%, 40%, 39% and 38%) respectively. The time to reach the maximum depth in Mosul city for the five scenarios of dam failure beginning from the end of failed breach development was (8.74, 7.49, 6.58, 5.99 and 5.52) hours respectively. 290 Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt

Reservoir water surface level (m.a.s.l) 4.E+05 EL.330 4.E+05 EL.320

3.E+05 EL.310

3.E+05 EL.300 EL.290 2.E+05

2.E+05

1.E+05

5.E+04

0.E+00 0 10 20 30 40 50 60 70 80 90 100 Distance from the dam to downstream (Km)

Figure (3) Flood discharges from the failed dam along the Tigris River reach

Reservoir water surface

level (m.a.s.l) 300 EL.330

290 EL.320 280 EL.310 270 EL300 260 EL.290 250

240 230

220 210

200 0 10 20 30 40 50 60 70 80 90 100 11 Distance from the dam to dow nstream (Km)

Figure (4) Water elevation of the flood wave along the Tigris River reach

Reservoir water surface 40 level (m.a.s.l) EL.330 EL.320

EL.310 30 EL.300 EL.290

20

10 0 10 20 30 40 50 60 70 80 90 100 Distance from the dam to dow nstream (Km )

Figure (5) Water depths of the flood wave along the Tigris River reach

290 Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt

Reservoir water surface level (m.a.s.l) 4.E+05 EL.330 4.E+05 EL.320

3.E+05 EL.310

3.E+05 EL.300 EL.290 2.E+05

2.E+05

1.E+05

5.E+04

0.E+00 Al-Bayan Center for Planning0 and10 Studies 20 30 40 50 60 70 80 90 100 Distance from the dam to downstream (Km)

Figure (3) Flood discharges from the failed dam along the Tigris River reach

Reservoir water surface

level (m.a.s.l) 300 EL.330

290 EL.320 280 EL.310 270 EL300 260 EL.290 250

240 230

220 210

200 0 10 20 30 40 50 60 70 80 90 100 Distance from the dam to dow nstream (Km)

Figure (4) Water elevation of the flood wave along the Tigris River reach

Reservoir water surface 40 level (m.a.s.l) EL.330 EL.320

EL.310 30 EL.300 EL.290

20

10 0 10 20 30 40 50 60 70 80 90 100 Distance from the dam to dow nstream (Km )

Figure (5) Water depths of the flood wave along the Tigris River reach Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt 291

Reservoir water surface level (m.a.s.l) 12 EL.330 11 EL.320 10 EL.310

9 EL.300 8 EL.290 7 6

5 4 3 2 0 10 20 30 40 50 60 70 80 90 100

Distance from the dam to dow nstream (Km)

Figure (6) Time to reach maximum depth of the flood wave along the Tigris River reach Figure (6) Time to reach maximum depth of the flood wave along the Tigris River reach

12 The maximum flood wave maps were determined for the Tigris River reach between dam site and Mosul city. Figure (7) shows the inundation area for the whole river reach for the fourth scenario of dam failure (320 ma.s.l.). While Figures (8,9,10,11 and 12) show the inundation maps within Mosul city for the five scenarios of dam failure (sunny day dam failure), (290, 300, 310, 320, 330) meter above sea level. A comparison between two statuses for Mosul city before and after dam failure for the fourth scenario (320 m.a.s.l.) was shown in figure 13. During the evaluation and analysis of the inundation maps of the flood wave within Mosul city for the five scenarios of dam failure, it was observed that the maximum water elevation within the city will be about 235.2 meter above see level for the last scenario of dam failure (i.e. 330 m.a.s.l.) and the flood wave will cover about 53.6% of the Mosul city in which the areas of the eastern bank of the city will be subjected to flood slightly more than the western bank as shown in Table 3 due to the topographical reasons and the large extension of the building in the eastern bank of Mosul city. The flooded area from the dam site and Mosul city for the five dam failure scenarios ranges between (152 and 252) km2 respectively (Table 4) in which includes a very good agricultural lands cultivated by the villagers living there. The time to reach the maximum depth in Mosul city for the five scenarios of dam failure beginning from the end of failed breach development was (8.74, 7.49, 6.58, 5.99 and 5.52) hours respectively. SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK

The maximum flood wave maps were determined for the Tigris River reach between dam site and Mosul city. Figure (7) shows the inundation area for the whole river reach for the fourth scenario of dam failure (320 ma.s.l.). While Figures (8,9,10,11 and 12) show the inundation maps within Mosul city for the five scenarios of dam failure (sunny day dam failure), (290, 300, 310, 320, 330) meter above sea level. A comparison between two statuses for Mosul city before and after dam failure for the fourth scenario (320 m.a.s.l.) was shown in figure 13. During the evaluation and analysis of the inundation maps of the flood wave within Mosul city for the five scenarios of dam failure, it was observed that the maximum water elevation within the city will be about 235.2 meter above see level for the last scenario of dam failure (i.e. 330 m.a.s.l.) and the flood wave will cover about 53.6% of the Mosul city in which the areas of the eastern bank of the city will be subjected to flood slightly more than the western bank as shown in Table 3 due tothe topographical reasons and the large extension of the building in the eastern bank of Mosul city. The flooded area from the dam site and Mosul city for the five dam failure scenarios ranges between (152 and 252) km2 respectively (Table 4) in which includes a very good agricultural lands cultivated by the villagers living there. The time to reach the maximum depth in Mosul city for the five scenarios of dam failure beginning from the end of failed breach development was (8.74, 7.49, 6.58, 5.99 and 5.52) hours respectively.

13 Al-Bayan Center for Planning and Studies

292 Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt

Mosul Dam

Tigris River

Flood Elevations (m) Mosul City 290 281 272 263 254 245 236 227 218 209

Figure (7): The whole case study flood area for the worst scenario of dam failure (330 m.a.s.l.)

Flood depth (m) 36 30

24

18

12

6

0

Figure (8) Flood depth map at Mosul city for the fifth scenario of dam failure (330 m.a.s.l.)

14 SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK

Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt 293

Flood depth (m) 36 30

24

18

12

6

0

Figure (9) Flood depth map at Mosul city for the fourth scenario of dam failure (320 m.a.s.l.)

Flood depth (m) 36

30

24

18

12

6

0

Figure (10) Flood depth map at Mosul city for the third scenario of dam failure (310 m.a.s.l.)

15 Al-Bayan Center for Planning and Studies

294 Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt

Flood depth (m) 36 30 24 18 12 6 0

Figure (11) Flood depth map at Mosul city for the second scenario of dam failure (300 m.a.s.l.)

Flood depth (m)

30 24 18 12 6 0

Figure (12) Flood depth map at Mosul city for the first scenario of dam failure (290 m.a.s.l.)

Table (3) Percentages of flooded areas in the western and eastern banks of Mosul city

% of Flooded % of flooded % of flooded Water Flooded flooded area in area from the area from the level in area in the area from the total area of total area of the dam western the whole eastern the eastern the western (m.a.s.l.) bank (km2) area of bank bank bank Mosul city (km2) 16 290 32.41 21.5 35.61 33.08 34.50 300 39.29 23.33 43.18 35.90 40.10 310 44.13 25.03 48.49 38.50 44.30 320 50.68 26.75 55.69 41.16 49.60 330 53.69 30.02 58.99 46.18 53.60

294 Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt

Flood depth (m) 36 30 24 18 12 6 0

Figure (11) Flood depth map at Mosul city for the second scenario of dam failure (300 m.a.s.l.)

Flood depth (m)

30 24 18 12 6 0

SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK Figure (12) Flood depth map at Mosul city for the first scenario of dam failure (290 m.a.s.l.)

Table (3) Percentages of flooded areas in the western and eastern banks of TableMosul (3) city Percentages of flooded areas in the western and eastern banks of Mosul city

% of Flooded % of flooded % of flooded Water Flooded flooded area in area from the area from the level in area in the area from the total area of total area of the dam western the whole eastern the eastern the western (m.a.s.l.) bank (km2) area of bank bank bank Mosul city (km2) 290 32.41 21.5 35.61 33.08 34.50 300 39.29 23.33 43.18 35.90 40.10 310 44.13 25.03 48.49 38.50 44.30 320 50.68 26.75 55.69 41.16 49.60

330Thirteenth International53.69 Water30.02 Technology Conference,58.99 IWTC 1346.18 2009, Hurghada,53.60 Egypt 295

TableTable (4) (4) FloodedFlooded areasareas between between dam dam site siteand andMosul Mosul city for city the for five the dam five failure dam failure scenarios scenarios

Water level in the dam (m.a.sl.) 290 330 310 320 330 Flooding area (km2) 152 174 196 226 252

A A selected selected cross cross-section-section of of the inundation inundation area area within within Mosul Mosul city city (67.4 (67.4 km km downstream dam site) for the five scenarios of dam failure was drawn in Figure 13(a downstream dam site) for the five scenarios of dam failure was drawn in Figure and b) in which the maximum width of the inundation area was about 6.4 km. While a 13(a and b) in which the maximum width of the inundation area was about 6.4 discharge-water level curve for the flood wave at Mosul city for the five scenarios of damkm. failure While was a discharge-water shown in Figure level 13(c curve). This for relationship the flood willwave be at usefulMosul in city determining for the thefive required scenarios values of ofdam flowing failure discharge was shown during in Figure the flood 13(c). wave This forrelationship any selected will waterbe level.useful in determining the required values of flowing discharge during the flood wave for any selected water level. AsAs a a result result of of thethe inundationinundation of Mosul Mosul city city due due to to the the failure failure of ofMosul Mosul dam, dam, many many commercial,commercial, educational, educational, industrial, industrial, demographical demographical and and medical medical constructions constructions in- addition to major governmental and non governmental building will be submerged by in- addition to major governmental and non governmental building will be the expected flood wave. Table (5) shows some of the main constructions which will submerged by the expected flood wave. Table (5) shows some of the main be inundated in Mosul city only according to the worst case (fifth scenario) of Mosul damconstructions failure. which will be inundated in Mosul city only according to the worst case (fifth scenario) of Mosul dam failure. Evacuation and emergency plans Evacuation and emergency plans

Evacuation plans should be prepared and implemented by the local jurisdiction17 controlling inundation areas. The assistance of local civil defense personnel should be requested in preparation of the evacuation plan. State and local law enforcement agencies usually will be responsible for the execution of much of the plan. The plans may include delineation of the area to be evacuated; routes to be used; traffic control measures; methods for providing emergency transportation; special procedures for the evacuation and care of people from institutions such as hospitals, procedures for securing the perimeter and for interior security of the area; procedures for the lifting of the evacuation order and reentry to the area; the details indicating which organization are responsible for specific functions, equipment and personnel resources required. Finally a copy of the plan should be furnished informed to the dam agency. Al-Bayan Center for Planning and Studies

Evacuation plans should be prepared and implemented by the local jurisdiction controlling inundation areas. The assistance of local civil defense personnel should be requested in preparation of the evacuation plan. State and local law enforcement agencies usually will be responsible for the execution of much of the plan. The plans may include delineation of the area to be evacuated; routes to be used; traffic control measures; methods for providing emergency transportation; special procedures for the evacuation and care of people from institutions such as hospitals, procedures for securing the perimeter and for interior security of the area; procedures for the lifting of the evacuation order and reentry to the area; the details indicating which organization are responsible for specific functions, 296 Thirteenth International Water Technology Conference, IWTC 13 2009, Hurghada, Egypt equipment and personnel resources required. Finally a copy of the plan should be furnished informed to the dam agency.

Mosul dam Fig. (13a) Tigris river Flow direction

Mosul boundary

Selected Flood wave cross-section

Reservoir water surface level (m.a.s.l.) Fig. (13b)

Fig. (13b)

Fig. (13c)

Fig. (13c)

Figure (13) Inundation cross-section and rating curve for the flood wave along Tigris river reach at the middle of Mosul city

Table (5) Numbers of the expected inundated constructions in Mosul city according to the fifth scenario of dam failure

18 SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK

CONCLUSIONS The simplified dam break model SYMDBK was selected due to its simplicity and efficiency to simulate and predict the flood wave of a hypothetical Mosul dam failure at different failure scenarios. It was utilized to explore various aspects of the dam break hazard evaluation and to compute characteristics of downstream flood wave of Mosul dam failure. River channels geometry represented by cross- sections which are obtained from topographic maps and satellite images to describe the width variations of river channel and valley was fed to the SYMDBK using the GIS facilities. The main findings demonstrate that an area of 252 km2 along Tigris River between dam site and Mosul city which was classified as very good agricultural lands and was utilized by the villagers there will be flooded due to Mosul dam failure in the worst failure scenario, the maximum flood depth at Mosul city will be 25 meter and the flood elevation will be 235.2 meter above see level and the maximum flood discharge will be (207632) m3/sec with an average flow velocity 3.5 m/sec. The flood wave will cover about 54% of the Mosul city in which the areas of the eastern bank of the city will be subjected to flood slightly more than the western bank. Finally, the time to reach the maximum depth in Mosul city for the five scenarios of dam failure beginning from the end of failed breach development was (8.74, 7.49, 6.58, 5.99 and 5.52) hours respectively.

19 Al-Bayan Center for Planning and Studies

REFERENCES:

Attalla, T.A. (2002). A review on dams and breach parameters estimation. Master of Science Hydrosystem Engineering. Dept. of Civil and Environment. Virgina Polytechnic Institute and State University. Blacksburg. V.A. Babb, A.O. and T.W. (1968). Cataloge of dam disasters. Failure and Accidents. Bureau of Reclamation, Washington, D.C. Bureau of reclamation (1988). Downstream hazard classification guidelines. ACER. Tech. Memorandum No. 11, U.S. Dept. of the Interior. Denver, 57. Costa, J.E. (1985). Floods from dam failure. U.S. geological survey, open File, Rep. No. 85-560, Denver, 54. Fread, D.L. (1985). DAMBRK: The NWS dam-break flood forecasting model. National weather service. Office of Hydrology. Silver Spring Md. Fread, D.L. (1993). NWS FLDWAV model; The replacement of Dam-BRK for dam break flood prediction. Dam safety 93 proc.th 10 Annual ASDSO Conf. Association State Dam Safety Officials. Lexington, Ky, 177-184. Froehlich, D.C. (1987). Embankment dam break parameters. Hydraulic Eng. Proc. ASCE National Conf. on Hydraulic Eng., New York, 570-575. Froehlich, D.C. (1995a). Embankment dam breach parameters revisited. Water Resources Eng. Proc. 1995 ASCE Conf. on Water Resources Eng., New York. 887- 891. Froehlich, D.C. (1995b).Peak outflow from breached embankment dam. J. Water Resources Plan Manage. Div. Am. Soc. Civil Eng. 121(1),90-97. HEC-1 (1981). Flood Hydrograph package, USERs Manual, Hydrologic Eng. Center, U.S. Army Corps of Engineer, Davis, CA. MacDonald, T.C. and Langridge-Monopolis, J. (1984). Breaching characteristics of dam failure. J. Hydraulic Eng. 110(5), 567-586. Middlebrooks, T.A. Earth dam practice in the united stats, ASCE Centennial Transaction paper 2620, 1952. pp 697-722. Singh, K.P. and Snorrason, A. (1984). Sensitivity of outflow peaks and flood

20 SIMULATION TIGRIS RIVER FLOOD WAVE IN MOSUL CITY DUE TO A HYPOTHETICAL MOSUL DAM BREAK stages to the selection of dam breach parameters and simulation models, J. Hydrology, 68, 295- 310. Singh, V.P. and Scarlatos, P.D. (1985). Breach erosion of earthfill dam flood routing Beed Model, Research Report, Army office, Battelle, Research Triangle Park, North Carolina, 131p. Soil Conservation Service (SCS). (1981). Simplified dam breach routing procedure; Tech Release No. 66 (Rev.1), 39. Von Thum, J.L. and Gillette, D.R. (1990). Guidance on breach parameters. Internal Memorandum, U.S. Dept. of the Interior, Bureau of Reclamation Denver, 17. Wetmore, J.N. and D.L. Fread (1983). The NWS simplified dam break model exution, Brief national weather service office of hydrology. Silver Spring, Maryland. Wurbs, R.A. (1987). Dam breach flood wave models. Journal of Hydraulic Engineering, Vol. 113. WWW. The Mosul dam Press. (2007).

21