AmericanJournalofWaterResources SimulationofFloods
AmericanJournalofWaterResources SimulationofFloods
net/publication/329389505
CITATIONS READS
21 3,400
3 authors, including:
All content following this page was uploaded by Mohammed Sharif on 04 December 2018.
Received August 17, 2018; Revised October 13, 2018; Accepted October 25, 2018
Abstract Flooding in Delhi Segment of River Yamuna has a devastating effect on the life of the inhabitants,
particularly those residing close to the banks. The River Yamuna experienced floods 33 times during the last century,
often endangering the life of people as well as the important infrastructural facilities that exist along the banks of the
river. The present paper describes the development and application of a hydrodynamic model based on HEC-RAS –
modelling system developed by the Hydrologic Engineering Centre at United States Army Corps of Engineers - for
the simulation of floods in the Delhi segment of River Yamuna. The HEC-RAS model was first calibrated and
validated and then applied for the simulation of historical floods of 2010 and 2013. With the HEC-RAS model, the
vulnerability assessment of different bridges and barrages in the Delhi segment has been carried out. It can be
concluded that the results presented herein could provide valuable aid to policy makers in formulating mitigation
strategies to counteract the adverse impacts of flooding in the Yamuna River basin.
Keywords: Yamuna, HEC-RAS, flood, mitigation, Delhi
Cite This Article: Azhar Husain, Mohammed Sharif, and Mohammed Lateef Ahmad, “Simulation of Floods
in Delhi Segment of River Yamuna Using HEC-RAS.” American Journal of Water Resources, vol. 6, no. 4
(2018): 162-168. doi: 10.12691/ajwr-6-4-3.
Snead [5] presented the application of two unsteady flow hydraulic flood routing model, with time-varying roughness
hydraulic models, Mike 11 and HEC-RAS for flood updating to simulate flows through natural channels.
routing and visualization for the Mill Creek Watershed The authors solved Saint Venant’s equation using the
located in Cincinnati, Ohio. Zadeh et al. [6] estimated the quasi-steady dynamic wave and full dynamic wave theory.
extent of the flood zone and economic damages over an A case study of unsteady flood modelling through
8.2 km reach of the perennial LaeenSoo River in the HEC-RAS was carried out for the Karad - Kurundwad
northern Khorasan Province, Iran using HEC-GEORAS, a reach of the Krishna River. Sankhua et al. [22] focused
combination of HEC-RAS and Arc View GIS software. on the concepts of hydraulic flood routing model, with
An application of Mike 11 to the simulation of flow in time-varying roughness updating to simulate flows
Euphrates River in Iraq has been presented by Kamel [7]. through natural channels, based on the quasi-steady
The results of the study demonstrated that the model dynamic wave and full dynamic wave theory, emphasizing
provides a good simulation of the observed flows. Pareta the solving of the intricate Saint Venant’s equation.
[8] developed a model for simulating flood inundation Several studies have been carried out in Yamuna River
boundaries. With the model developed by Pareta [8], it basin as well. For example, Sharma et al. [23] focused on
was possible to precisely overview the extent of flooding the monitoring of the diffuse pollution characteristics of
in the area of interest. Sharif et al. [9] investigated trends the agricultural land confining the River Yamuna in Delhi.
in magnitude and timings of flows in the Upper Indus Vijay et al. [24] presented hydrodynamic simulation of the
basin – a climatically sensitive basin in the Himalayan river Yamuna under different designated flood flows to
region. Trends in several hydro-meteorological variables delineate the land availability under existing and modified
were investigated by Khattak et al. [10] in the upper Indus riverbed geometry including channel dredging and
River basin (UIRB). riverbed dressing. Kumar et al. [25] describes flood
An analysis of flow dynamic process and the quantification management strategies for the national capital region of
of the peak discharges by using HEC-RAS has been Delhi, India. Husain. A [26] describe the method flood
presented by Villazón et al. [11] for the Piraí River. The peak estimation at Hathinkund and Okhla barrage at
authors simulated unsteady flows using derived hydrographs Yamuna River for the national capital region of Delhi,
from measured hourly water levels as boundary conditions. India. An application of POT approach for estimation of
Al-Fahdawi [12] applied a numerical model to hydrodynamic floods of varying return periods has been described by
modelling of the Euphrates River in its reach between Kumar et al. [27]. The major objective of the present
Haditha Dam and Hit city. The author used HEC-RAS research is to a develop a hydrodynamic model for the
model to determine different parameters from a given simulation of floods in the Delhi segment of River
flood caused by a hypothetical foundation failure of Yamuna. The intent is to assess the vulnerability of
Haditha Dam. Yarrakula et al. [13] developed a flood hydraulic structures in the Delhi segment of River
forecasting model for the Subrnarekha River. The digital Yamuna to floods associated with different return-levels.
elevation Model (DEM) of the study area was prepared
using the high resolution CARTOSAT-1 imageries. Cecile
Ashwanden et al. [14] compared water surface profiles for 2. Study Area
Tar River generated from steady and unsteady flow
hydraulic models. Brych et al. [15] used HEC-RAS to River Yamuna is one of the major rivers of India and
prepare flood inundation maps for the Orlice valley in the the largest tributary of the Ganga. The river originates
Czech Republic. Extreme floods were considered to from the Yamunotri glacier of the lower Himalayas at an
calibrate the hydraulic model of the Orlice river system. elevation of 6320 m above mean sea level (Figure 1). The
Khattak et al. [16] describes the application of HEC-RAS river has a total length of 1376 km and drains an area of
model for the development of floodplain maps for the part 345848 km2. In the upper reaches of River Yamuna, there
of Kabul River that lies in Pakistan. are several hill streams that join together to form the main
Several flood modelling studies have been carried out stream. There are four main rivers that join Yamuna in the
in Indian basins over the last decade. Vijay et al. [17] high Himalayan ranges. These include Rishi Ganga, which
describes a hydrodynamic model called River CAD for the joins on the right bank of the Yamuna, whereas the
assessment of extent of inundation, and evaluation Hanuman Ganga joins on the left bank. In the lower
of different development scenarios in the basin of interest. Himalayan ranges, the Yamuna River receives Kamal,
Mazumder [18] describes a methodology for the Tons, Giri and Beta on its right bank and Aglag and Asan
determination of waterway for a bridge using detailed on left banks. The Chambal, Betwa, Sindh and Ken are the
topographic, hydrologic, and morphological investigations. important tributaries joining Yamuna on the right bank in
Pramanik et al., [19] used extracted river cross-sections to the plains. The Hindon river joins River Yamuna on its
simulate the magnitude of flood in the deltaic reaches of left bank. Among all these tributaries, Tons at hills and
Brahmani river basin located in the eastern India. Timbadiya Chambal at plains are the most important tributaries in
et al. [20] developed an integrated hydrodynamic model of terms of their discharges. The Tons is the principal source
the lower Tapi River, India. Firstly, the one-dimensional of water in mountainous range and generally carries more
hydrodynamic model was calibrated for Manning’s roughness water than the main stream. In plains, during non-
of river channel and, subsequently one-dimensional and monsoon period, River Chambal contributes about 5-10
two-dimensional integrated hydrodynamic model was times more water to the Yamuna than its own flow.
used to ascertain the sensitivity of Manning’s ‘n’ on However, since the year 2003, there is a significant
coastal flood plain depth of the lower Tapi River. reduction in the water quantity that River Chambal
Doiphode and Ravindra [21] focused on the concepts of discharges into the Yamuna River.
American Journal of Water Resources 164
There are five barrages that control the flow of Yamuna found within the highly braided river stretch. These
waters during its journey from the Himalayas to the barrages are operated only during the monsoon period,
National Capital Territory of Delhi (NCT). These barrages when excess of water leads to flooding in the adjoining
are located at Dak Patthar (about 160 km from origin in areas. The total discharge at Tajewala (now Hathnikund)
Uttaranchal), Hathnikund (172 km distance from origin, often exceeds 500,000 cusec during the monsoon peak
just at foothills in Haryana), Wazirabad (in NCT of Delhi, flows, and the Hathnikund Barrage itself has been
396 km distance from origin), Okhla (in NCT – Delhi, 418 designed for a 100-yr peak discharge of 776,900 cusec
km distance from origin, and at Mathura (Near Gokul (22,000 cumec). On the other hand, the flow during the
village in Uttar Pradesh about 570 km distance from summer decreases to only a few thousand cusec [28].
origin). These barrages are the major water abstraction
locations on the river. The water is contributed into the
Yamuna River, not only through its tributaries, but also by 3. Data
the canals and drains originating from various urban centers.
The three major cities - Delhi (the capital of India), In the present work, simulation of water surface profile
Agra (the city known for Taj Mahal) and Mathura (the at different locations of interest has been carried out using
birthplace of Lord Krishna) are situated on the banks of HEC-RAS. Two main types of data, namely the geometric
River Yamuna. The Delhi stretch of the river Yamuna data, and hydraulic data, are required for the hydraulic
from Wazirabad barrage to Okhla barrage is around 23 km analysis of the stream channel geometry. These data along
long, and is located between 28° 28′ 05″–28° 54′ 36″ N with the flow data have been used in this research. The
and 77° 09′–77° 24′ E at the downstream. The Yamuna basic geometric data include establishing the connectivity
enters in the National Capital Region (NCR) of Delhi, of the river system (river network); cross sections; reach
approximately 1.65 Km north of Palla Village. It runs for lengths, and hydraulic structures such as bridge piers.
about 45 kms in the southeast direction before leaving The River Schematic in HEC-RAS allows to define
NCR of Delhi at a point to the east of Jaitpur downstream river network and reference cross-sections and control
of Okhla Barrage. The total area of the river in the Delhi structures in the river as well as to graphically obtain an
stretch is about 9,700 ha of which 1,600 ha of land is overview of model information in the current simulation.
under water (river extent) and the remaining 8,100 ha is The river system schematic should be developed before
dry land (floodplains). Three barrages namely Wazirabad any other data can be entered. The river centerline of
barrage, Yamuna/ITO barrage and Okhla barrage are the study area contains x- and y-coordinate data in a 2-D
165 American Journal of Water Resources
plane so as to spatially connect the unsteady flow models data. For the manual entry of the data, the details of
to the corresponding terrain models. The model automatically 77 cross-sections along the study reach of 23 km between
looks for a companion file to the image. This file contains Wazirabad and Okhla Barrage was obtained from
information about the image, including the coordinate the Irrigation and Flood Control Department of
system and the extents of the image. Figure 2 shows the the government of the National Capital Territory of
River system schematic for the study reach. Delhi. A typical cross section at Old Railway Bridge is
shown in Figure 3. Using the data from the DEM,
3.1. Geometric Data the geometric data for the study reach can be created
using HEC- GeoRAS extension of ArcGIS. A typical
The geometric data can be created either by using cross section at Old Railway Bridge is shown in
Geo-RAS or through manual entry of the cross-section Figure 3.
210 Legend
204
202
200
198
0 200 400 600 800 1000 1200
Station (m)
Figure 3. Cross sections of Yamuna River at section u/s of Old Railway Bridge
American Journal of Water Resources 166
215 Legend
WS Flood 2010
210 Ground
205
200
Wazirabad Barrag
Nizamuddin Bri...
195
DND Bridge
ITO Barrage
ISBT Bridge
190
0 5 10 15 20 25
Main Channel Distance (km)
Figure 4. Simulated and observed water surface profile for 2010 flood
167 American Journal of Water Resources
Table 1. Output of HEC-RAS simulations for the 2010 flood Figure 5 shows the simulated and observed water levels
Water Surface Elevation (m) for the 2013 flood. The simulated water surface profile is
shown in red, whereas the dots indicate the observed water
River Station (Number) Simulated Observed Observed-Simulated
levels on 19 June 2013. The numerical values of simulated
Wazirabad Barrage (775) 208.98 208.35 -0.63 water levels are shown in Table 2. It can be seen from
Old Railway Bridge (535) 206.96 207.0 0.04 Table 2 that the difference between the simulated and
ITO Barrage (395) 206.16 205.5 -0.66 observed water levels is negligible at the Old Railway
Bridge, whereas there was some overestimation of water
Okhla Barrage (2) 201.38 200.45 -0.93
levels at the other three stations. However, the amount of
overestimation is not significantly high.
It can be clearly seen from Table 1 that the difference
between the observed and simulated water levels is Table 2. Observed and HEC-RAS simulated water levels under 2013
negligible at Old Railway Bridge. At the other three flood
locations, namely Wazirabad Barrage, ITO Barrage and Water Surface Elevation (m)
Okhla Barrage the simulated water levels were found to be River Station (Number) Simulated Observed Observed-simulated
slightly higher than the observed levels. Thus, it may be
Wazirabad Barrage (775) 209.21 208.45 -0.76
concluded that the model tends to overestimate the water
levels at some locations, although the magnitude of Old Railway Bridge (535) 207.04 207.18 0.14
overestimation is not significantly high. ITO Barrage (395) 206.23 205.5 -0.73
Case 2: 2013 Flood Okhla Barrage (2) 201.52 200.6 -0.92
The HEC-RAS model was used to simulate 2013 flood.
215 Legend
WS Flood 2013
210 Ground
205
200
Wazirabad Barrag
Nizamuddin Bri...
195
DND Bridge
ITO Barrage
ISBT Bridge
190
0 5 10 15 20 25
Main Channel Distance (km)
by increasing the nodes in the network currently considered [13] Yarrakula, K., Deb, D., and Samanta, B. (2010). Hydrodynamic
in the modelling process. For greater reliability in the modeling of Subernarekha River and its floodplain using remote
sensing and GIS techniques. Journal of Scientific &Industrial
assessment of vulnerability of bridges/barrages under Research vol.69:529-536.
inputs to HEC-RAS from a larger number of RCMs, [14] Cecile, A., Keren, J.C., Sean, R. (2010). Inundation mapping using
rather than a single model should be utilized. Depending hydraulic model’s case studies of steady and unsteady models on
upon the availability of the data, the future studies would the Tar river, NC’, Advanced Hydrologic Prediction Service
(AHPS) web pages published by North Carolina University,
employ a larger ensemble of RCMs and development PP. 1-10.
scenarios for understanding the range of climate change [15] Brych, k., Dittrt F., Elias V., (2002). Development of flood
impacts on the hydrologic regime in the Yamuna River boundary maps of urban areas using HEC – RAS’ , Proceedings of
basin. Information related to streamflow projections could the fourth International FRIEND conference, held at Cape town,
be incorporated into planning activities. The model developed South Africa. PP.237-242.
[16] Khattak, M. S., Anwar, F, Saeed, T. U., Sharif, M. (2016).
herein could be applied to simulate future floods in the Floodplain Mapping Using HEC-RAS and ArcGIS: A Case
basin, and therefore aid in the development of strategies Study of Kabul River, Arb J. of Sc. & Engg. (Springer), 41(4),
for mitigating the adverse impacts of floods in the basin. 1375-1390.
[17] Vijay, R., Sargoankar, A. and Gupta, A. (2007). Hydrodynamic
simulation of river Yamuna for riverbed assessment: A case study
References [18]
of Delhi region, Environ Monit Assess, 130: 381-387.
Mazumder, S. K. (2009). Determination of waterway under a
bridge in Himalayan region - some case studies, Journal of the
[1] Ahmad, S. (1999). Comparison of one-dimensional and two-dimensional Indian Roads Congress, paper no. 551.
hydrodynamic modeling approaches for Red River Basin. Natural [19] Pramanik, N., Panda, R. K. and Sen, D. (2010). One dimensional
Resources Institute Facility for Intelligent Decision Support, hydrodynamic modeling of river flow using DEM extracted river
University of Manitoba Winnipeg, Manitoba, Canada. cross-sections, Water Resour Manage, 24:835-852.
[2] Tate, E. C., and Maidment, D. (1999). Floodplain mapping using [20] Timbadiya, P. V., Patel, P .L, and Porey, P. D. (2012). Calibration
HEC-RAS and ArcView GIS, CRWR Online Report 99-1, Center of channel and flood plain roughness using 1d/2d integrated
for Research in Water Resources, Bureau of Engineering Research, hydrodynamic model: a study of surat city on lower tapi river.
University of Texas, Austin, TX 78712-4497. Proceedings of the 10th Intl. Conf.on Hydroscience and
[3] Husain. A (2017). Flood modelling by using HEC-RAS. Engineering, Nov. 4-7, Orlando, Florida, U.S.A.
International Journal of Engineering Trends and Technology [21] Doiphode, S. L., Oak, R. A., (2012). Dynamic Flood routing
(IJETT) – Volume 50 Number 1 August 2017. model and unsteady flow modeling .A case study of Upper
[4] Yang, X. and Rystedt, B. (2002). Predicting flood inundation Krishna River.
and risk using gis and hydrodynamic, Model: A case study at [22] Sankhua, R. N., Sathe, B. K and Srivastava, A. K. (2012). A case
Eskilstuna, Sweden, MFDM-04, Indian Cartographer, 183-191. study on dynamic wave routing and unsteady flood modelling of
[5] Snead, D. B. (2000). Development and application of unsteady part of krishna basin with HEC-RAS. India water week water,
flood models using geographic information systems. thesis, Faculty energy and food security.
of the Civil Engineering Department, University of Texas, Austin. [23] Sharma, D., Gupta, R., Singh, R. K. and Kansal, A. (2012).
[6] Zadeh, A. H., Ahmadi, M. Z., Sharifi, M. B. and Masoudian, M. Characteristics of the event mean concentration (EMCs) from
(2005). Flood zoning using the HEC-RAS hydraulic model in a rainfall runoff on mixed agricultural land use in the shoreline zone
GIS environment, Journal of Environmental Hydrology, vol. 13, of the Yamuna River in Delhi, India, Appl Water Sci, 2: 55-62.
paper 2. [24] Vijay, R., Sargoankar, A. and Gupta, A.(2009). A hydrodynamic
[7] KameL, A. H. (2008). Application of a hydrodynamic mike 11 approach to address Yamuna riverbed development in Delhi. J.
model for the Euphrates River in Iraq, Slovak Journal of Civil Civ. Eng. 36: 1155-1163.
Engineering. [25] Kumar, M, Sharif, M, and Ahmed, S., (2017). Flood Risk
[8] Pareta, K. (2012). Flood monitoring simulation model: a case Management Strategies for National Capital Territory of Delhi,
study of Red river delta, Vietnam, International Journal of India. ISH Journal of Hydaulic Engineering.
Research in Engineering. IT and Social Sciences, Vol 2, No. 6: [26] Husain. A (2015).Flood Peak Estimation at Hathnikund and
2250-0588. Okhla Barrage. Journal of Civil Engineering and Environmental
[9] Sharif, M, Archer, D, Fowler, H, and Forsythe, N. (2013). Trends Technology: volume 2, number 13, 2015 pp 40-45
in magnitude and timings of flow in the Upper Indus River Basin, [27] Kumar, M, Sharif, M, and Ahmed, S. (2018). Flood Estimation in
Journal of Hydrology and Earth System Sciences. River Yamuna (India) Using Peak-Over-Threshold Method. ISH
[10] Khattak, M. S, Babel, M. S., and Sharif, M. (2011). Hydro- Journal of Hydaulic Engineering.
meteorological trends in upper Indus river basin. J. Climate [28] Gopal, B. and Calhan, M. (2007). River Yamuna from Source to
Research, 46: 103-119. Delhi: Human Impacts and Approaches to Conserv ation. Martin,
[11] Villazón, M. F., Aranibar, C. R. and Willems, P. (2009). P., Gopal, B. and Southey, C. (Editors).
Hydrodynamics unsteady flow model applied to Pirai River for [29] Prasuhn, A.L. Fundamentals of Hydraulic Engineering. Oxford
flood events estimation in Santa Cruz-Bolivia, 1ST International University Press, New York, NY, 1992.
Congress of Hydro climatology, Cochabamba-Bolivia. [30] Hydrologic Engineering Center (RD-42), (1995). "Flow Transient
[12] Al-Fahdawi, S. A. (2009). Numerical modeling of flood wave in Bridge Back water Analysis", U.S. Army Croup of Engineers,
behavior with meandering effects (Euphrates River, Haditha-Hit), Davis, California.
thesis, Building and Construction Engineering Department, [31] HEC-RAS, River Analysis System, (2010). "User’s Manual",
University of Technology, Iraq. Version 4.1, U.S. Army Corps of Engineers Hydrologic
Engineering Center, Davis, California, November.