Research Article Volume 6 Issue 1
1Department of Land and Water Conservation Engineering, Faculty of Agricultural Engineering and Technology, PMAS-Arid Agriculture University, Rawalpindi 46000, Pakistan
2Department of Farm Machinery and Precision Engineering, Faculty of Agricultural Engineering and Technology, PMAS Arid Agricultural University, Rawalpindi 46000, Pakistan
3Department of Energy system Engineering, Faculty of Agricultural Engineering and Technology, PMAS Arid Agricultural University, Rawalpindi 46000, Pakistan
4Department of Structure and Environmental Engineering, Faculty of Agricultural Engineering and Technology, PMAS Arid Agricultural University, Rawalpindi 46000, Pakistan
Correspondence: Muhammad Mashood Javaid, Department of Energy system Engineering, Faculty of Agricultural Engineering and Technology, PMAS Arid Agricultural University, Rawalpindi 46000, Pakistan
Received: January 10, 2022 | Published: January 17, 2022
Citation: Nadeem MU, Waheed Z, Ghaffar AM, et al. Application of HEC-HMS for flood forecasting in hazara catchment Pakistan, south Asia. Int J Hydro. 2022;6(1):7-12. DOI: 10.15406/ijh.2022.06.00296
Floods have become more severe and frequent as a result of climate change around the world, posing a hazard to public safety and economic development. This study investigates the use of distributed hydrological models in flash flood risk management in a small watershed in Hazara, Pakistan, with the goal of improving Pakistan's early warning lead time. First, the HEC-HMS model was built using geographic data and the river network's structure, then calibrated and verified using eight high rainfall events from 2013. demonstrating that the HEC-HMS model could simulate floods in the research area Second, given that rainfall and flood events have happened, this paper proposes an analysis approach for a flood forecasting and warning system, as well as criteria for sending urban-stream flash flood alerts based on rainfall, in order to provide sufficient lead time. The DEMs (digital elevation models) of the research regions were processed using HEC-Geo HMS, an ArcView GIS tool for catchment delineation, terrain pre-processing, and basin processing. The model was calibrated and verified using previously observed data. The proposed flood prediction and risk reduction methodology is nonstructural. The Hydrologic Modeling System (HEC-HMS), which provides a sufficient lead time forecast and computes the runoff/stage threshold conditions, is at the heart of the flood warning application. For flood risk assessment, data from the Pakistan Meteorological Department (PMD) is entered into a hydro-meteorological database and then into the HEC-HMS. A server-client application was utilised to visualise the real-time flood scenario and send out an early warning message. The outcomes of this study will be used to develop flood validation measures in the Hazara stream watershed to deal with potential flash floods.
Keywords: HEC-HMS, hydrological modeling, flash flood risk management, Hazara watershed Pakistan
Many of the ideological debates in pragmatic human ecology, critical realism political ecology, and the poststructural social nature thesis can be studied through flood danger research in the Hazara watershed. These discussions centre on the practical aspects of flood prevention, as well as the hazard-scape idea, which is founded on the geography concept of environment. Flash floods are typically caused by heavy rain, but they can also be caused by dam failures, levee failures, or even ice jams in rivers during the winter and spring months.1 As the world's population increases, so does the risk of flash flooding. Runoff from impervious surfaces, such as concrete or compacted bare soils, contributes a significant volume of runoff, which can quickly damage roads and structures.2 As a result, the intensity of the rainfall event and the spatial distribution of land cover have a major impact on the watershed's hydrologic response. As a result, flash floods should be examined at both spatial and temporal scales to identify the exact response time of the watershed (concentration time).
The primary purpose of a flood forecasting system is to reduce the loss of life, property, and business by providing accurate warnings with ample lead time to users and emergency managers specially in Pakistan.3 The lead time is determined by the hydrologic model's efficiency and forecast quality. The only way to combine meteorological forecasts with hydrologic models in flood forecasting systems is to use forecasted hydro-meteorological data.4 In meteorological forecasting, regional numerical weather prediction models are used. There are many different types of hydrologic models, and their use varies depending on the sort of findings needed and the availability of hydrological data. The hydrological processes can be characterised empirically, intellectually, or physically, and the deterministic models can be classified based on whether they provide a lumped or distributed description of the examined area.
The HEC-HMS Model was utilised for runoff prediction in this study since it is vital for early flood warning. As a substitute for HEC-1, the US Corps of Engineers' Hydrologic Engineering Center (HEC) took the initiative. HEC-HMS enhances HEC-1's capabilities by allowing for distributed modelling and continuous simulation. The application of HEC-HMS, as well as the precision of calibration and validation, are crucial for the development of EFWS. The HEC-HMS model can simulate precipitation-runoff and routing mechanisms in both natural and controlled environments5 HEC-HMS is a straightforward conceptual model that has been successfully developed by a number of academics throughout the world. It is superior for peak flow simulation than the Revitalized Flood Hydrograph model because of the semi-distributed modelling idea. Because the soil moisture accounting model, which is based on the natural hydrological cycle in natural catchments, is regarded as a suitable model for rural catchments. When the semi-distributed model HEC-HMS and the fully-distributed model Basin Pollution Calculation Center (BPCC) were used to simulate the rainfall runoff process in the mountain area, the findings revealed a little disagreement between the two models.
In Arc-GIS 10.7.1, the Digital Elevation Model (DEM) of the Hazara watershed was created using SRTM 90m (Shuttle Radar Topographic Mission). Aerial rainfall is calculated utilising the Thiessen polygon method (by multiplying gauge weights of metrological stations with situ-gauge rainfall). Arc-GIS was also used to demarcate the Hazara watershed. Diva-GIS is used to download Land Use/Land Cover (LU/LC) data. All maps and applicable formulas from prior studies are utilised to generate the parameters. Finally, the watershed is identified and then integrated into the HEC-HMS model for calibration and validation. Time of concentration, curve number, lag time, and area of watershed is calculated from suitable formulas and maps devolved in this study.
GIS (geographic information systems) has become an important part of hydrologic studies due to the spatial character of the causes and precipitation that affect hydrologic processes. In the parameterization of distributed hydrologic models, the usage of geographic information systems (GIS) is crucial. To avoid gross simplifications caused by parameter lumping at the river basin scale during representation, this is done. GIS applications are utilised to extract hydrologic information such as flow direction, flow accumulation, watershed boundaries, and stream networks from a DEM (digital elevation model). To assess the model's performance, this study used GIS and HEC-HMS.
Study area
Hazara stream sub-watershed is located in Pakistan city at Latitude: 35° 14' 60.00" N, Longitude: 73° 29' 59.99" Containing the area of 340km2. It is a region of Pakistan's Khyber Pakhtunkhwa province located in the northeast. It is situated to the east of the Indus River. The nine rainfall stations are present in Hazara watershed which is used to calculate aerial rainfall of stations Figure 1.
Thiessen polygon
In the Thiessen polygon, there is only a single point input feature. The associated point of a Thiessen polygon is closer to any point within it than any other point input feature. Perpendicular bisectors are generated for each triangle edge to construct the edges of the Thiessen polygons.6 Where the bisectors connect determines where the Thiessen polygon vertices are placed. The coordinate system is first set to WGS1984. The shapefile for Pakistan can be downloaded from the DIVA-GIS website. The vector data of metrological stations is acquired from the PMD website (Pakistan metrological department website), and the files are then imported into ArcGIS. Thiessen polygon is made in ArcGIS after that their weighted rainfall is measured which is known as Aerial Rainfall Figure 2, Table 1.
OBJECTID* |
Shape* |
Input_FID |
Shape |
Shape Area (m2) |
Area (m2) |
Parameter (m) |
weighted area (m2) |
1 |
Polygon |
4 |
2.4 |
0.2 |
2.53E+09 |
237700.1 |
0.331521 |
2 |
Polygon |
8 |
2.4 |
0.1 |
1.31E+09 |
249240.5 |
0.172167 |
3 |
Polygon |
7 |
0.6 |
0 |
1.69E+08 |
56305.4 |
0.022205 |
4 |
Polygon |
3 |
0.9 |
0 |
3.12E+08 |
93664.2 |
0.040924 |
5 |
Polygon |
1 |
2.2 |
0.1 |
1.43E+09 |
203963.3 |
0.187086 |
6 |
Polygon |
0 |
0.8 |
0 |
1.92E+08 |
74402.7 |
0.025234 |
7 |
Polygon |
6 |
1.2 |
0 |
4.4E+08 |
123391.2 |
0.057667 |
8 |
Polygon |
5 |
1.9 |
0.1 |
1.09E+09 |
180952 |
0.143663 |
9 |
Polygon |
2 |
0.7 |
0 |
1.49E+08 |
62096.4 |
0.019534 |
Table 1 Gauge weight calculation of study area
Digital elevation model
A digital elevation model (DEM) is a three-dimensional computer graphical representation of elevation data that is used to depict terrain on a planet, moon, or asteroid. The term "global DEM" refers to a discrete global grid. DEMs are the most common substrate for digitally produced relief maps, and they're often utilised in GIS systems.7 DSMs and DTMs that only provide height information without any additional surface characteristics are referred to as DEMs. In Arc-GIS 10.7, the Hazara watershed Digital Elevation Model (DEM) is based on SRTM 90m (Shuttle Radar Topographic Mission) Figure 3.
Delineation
Watershed delineation is one of the most common activities in hydrologic and environmental study. Watersheds can be estimated automatically using GIS technology employing digital elevation models (DEMs), which provide a fair representation of the topography. Since the early 1980s, automated watershed delineation techniques have been available and used in a number of GIS systems and specialised applications. These techniques work regardless of DEM quality and have become more popular as better resolution DEMs have become available, allowing for more detailed topographical depiction and analysis. The growth of desktop computing power, which allows complex procedures to be completed locally and swiftly enough, is another element contributing to the attraction of these techniques. These two considerations, together with a growing industry focus on drainage basin solutions, are putting more pressure on automated systems to deliver results quickly Figure 4.
Land use land cover (LU/LC)
The classification of raw satellite data into "land use and land cover" groupings yields Land Use and Land Cover (LULC) data. The way humans use the land and engage in socioeconomic activities is referred to as land use. Land use refers to how people use the landscape, such as whether it is developed, used for agricultural, or a combination of the two. Land cover is the physical material that covers the earth's surface. The amount of forest, wetlands, impermeable surfaces, agricultural, and other land and water types that cover a region is shown by land cover data. Land cover data can aid our understanding of the terrain and how humans interact with it Figure 5.8
HEC-HMS calibrations
HEC-HMS is a computer programme that simulates precipitation-runoff processes in dendritic drainage basins.9 It's meant to be used in a variety of settings and to solve as many problems as possible. Water supply and flood hydrology in large river basins, as well as local urban or natural watershed runoff, are all covered Figure 6.10,11 The program's hydrographs are used either directly or in conjunction with other software in studies of water availability, urban drainage, flow forecasting, future urbanisation impact, reservoir spillway design, flood damage reduction, floodplain regulation, and system operation. Two sub-basins were generated having different areas and different parameters calculated by suitable formulas and methods Figure 7,8.12–14,1
As stated in the emergence, each component of HEC-HMS replicates a different aspect of the precipitation-runoff process within a sub-basin of the basin.A set of parameters that indicate the component's unique properties, as well as quantitative properties that describe the physical processes, are essential for component representation (Hazara watershed). The standardized model parameters of each of the components represented in this model are shown in Table 1 and 2. Aside from the preset sub-areas, parameters were calibrated sequentially by adjusting their values until a good agreement between the observed and simulated hydrographs was attained Table 3.
Hydrological element |
Drainage (Km2) |
peak Discharge (m3/s) |
time of Peak |
Volume (mm3) |
Subbasin-1 |
239.5 |
369.4 |
14Jan2013, 19:00 |
45.2 |
Subbasin-2 |
240 |
293.8 |
14Jan2013, 19:00 |
32.66 |
Reach-1 |
239.5 |
340.7 |
14Jan2013, 20:00 |
45.2 |
Sink-1 |
479.5 |
624.9 |
14Jan2013, 19:00 |
38.93 |
Table 2 Final summary of HEC-HMS calibrations
Date |
Time |
inflow (m3/sec) |
outflow (m3/sec) |
13-Jan-13 |
0:00 |
0 |
0 |
13-Jan-13 |
1:00 |
0 |
0 |
13-Jan-13 |
2:00 |
0 |
0 |
13-Jan-13 |
3:00 |
0 |
0 |
13-Jan-13 |
4:00 |
0 |
0 |
13-Jan-13 |
5:00 |
0 |
0 |
13-Jan-13 |
6:00 |
0 |
0 |
13-Jan-13 |
7:00 |
0 |
0 |
13-Jan-13 |
8:00 |
0 |
0 |
13-Jan-13 |
9:00 |
0 |
0 |
13-Jan-13 |
10:00 |
0 |
0 |
13-Jan-13 |
11:00 |
0 |
0 |
13-Jan-13 |
12:00 |
0 |
0 |
13-Jan-13 |
13:00 |
0 |
0 |
13-Jan-13 |
14:00 |
0 |
0 |
13-Jan-13 |
15:00 |
0 |
0 |
13-Jan-13 |
16:00 |
0 |
0 |
13-Jan-13 |
17:00 |
0 |
0 |
13-Jan-13 |
18:00 |
1.2 |
0.3 |
13-Jan-13 |
19:00 |
26.2 |
6.7 |
13-Jan-13 |
20:00 |
64 |
30.4 |
13-Jan-13 |
21:00 |
71.1 |
57.9 |
13-Jan-13 |
22:00 |
40.5 |
61 |
13-Jan-13 |
23:00 |
18 |
40 |
14-Jan-13 |
0:00 |
9.3 |
21.1 |
14-Jan-13 |
1:00 |
7.6 |
11.6 |
14-Jan-13 |
2:00 |
13.3 |
9.8 |
14-Jan-13 |
3:00 |
13.9 |
12.6 |
14-Jan-13 |
4:00 |
9.7 |
12.6 |
14-Jan-13 |
5:00 |
11.9 |
10.9 |
14-Jan-13 |
6:00 |
14.6 |
12.3 |
14-Jan-13 |
7:00 |
15.9 |
14.4 |
14-Jan-13 |
8:00 |
27.6 |
18.3 |
14-Jan-13 |
9:00 |
48.9 |
30.4 |
14-Jan-13 |
10:00 |
64.6 |
48.3 |
14-Jan-13 |
11:00 |
81.6 |
64.7 |
14-Jan-13 |
12:00 |
79.3 |
77.2 |
14-Jan-13 |
13:00 |
103.1 |
84.3 |
14-Jan-13 |
14:00 |
119.4 |
102.5 |
14-Jan-13 |
15:00 |
117.6 |
115.1 |
14-Jan-13 |
16:00 |
156.5 |
126 |
14-Jan-13 |
17:00 |
247.3 |
170.4 |
14-Jan-13 |
18:00 |
348.3 |
252.9 |
14-Jan-13 |
19:00 |
369.4 |
331.2 |
14-Jan-13 |
20:00 |
283.4 |
340.7 |
14-Jan-13 |
21:00 |
142.5 |
264.1 |
14-Jan-13 |
22:00 |
60.2 |
151.6 |
14-Jan-13 |
23:00 |
29.3 |
74.2 |
15-Jan-13 |
0:00 |
9.6 |
35.1 |
15-Jan-13 |
1:00 |
3 |
14 |
15-Jan-13 |
2:00 |
0.8 |
5 |
15-Jan-13 |
3:00 |
7.6 |
3.3 |
15-Jan-13 |
4:00 |
15.7 |
8.5 |
15-Jan-13 |
5:00 |
12.6 |
13.3 |
15-Jan-13 |
6:00 |
7.1 |
11.5 |
15-Jan-13 |
7:00 |
11.3 |
9.1 |
15-Jan-13 |
8:00 |
14.8 |
11.6 |
15-Jan-13 |
9:00 |
10.8 |
13.2 |
15-Jan-13 |
10:00 |
6.7 |
10.4 |
15-Jan-13 |
11:00 |
5.5 |
7.3 |
15-Jan-13 |
12:00 |
3.4 |
5.4 |
15-Jan-13 |
13:00 |
10.9 |
5.6 |
15-Jan-13 |
14:00 |
32.5 |
14.6 |
15-Jan-13 |
15:00 |
42.4 |
30.7 |
15-Jan-13 |
16:00 |
72.2 |
46.6 |
15-Jan-13 |
17:00 |
64 |
64.4 |
15-Jan-13 |
18:00 |
26.9 |
55.5 |
15-Jan-13 |
19:00 |
14.2 |
30.5 |
15-Jan-13 |
20:00 |
7.6 |
16.4 |
15-Jan-13 |
21:00 |
2.5 |
8.4 |
15-Jan-13 |
22:00 |
0.8 |
3.5 |
15-Jan-13 |
23:00 |
0.2 |
1.3 |
16-Jan-13 |
0:00 |
0.1 |
0.4 |
16-Jan-13 |
1:00 |
0 |
0.1 |
16-Jan-13 |
2:00 |
0 |
0 |
16-Jan-13 |
3:00 |
0 |
0 |
16-Jan-13 |
4:00 |
0 |
0 |
16-Jan-13 |
5:00 |
0 |
0 |
16-Jan-13 |
6:00 |
0 |
0 |
16-Jan-13 |
7:00 |
0 |
0 |
16-Jan-13 |
8:00 |
0 |
0 |
16-Jan-13 |
9:00 |
0 |
0 |
16-Jan-13 |
10:00 |
2.5 |
0.6 |
16-Jan-13 |
11:00 |
2.7 |
2.1 |
16-Jan-13 |
12:00 |
7 |
3.5 |
16-Jan-13 |
13:00 |
6.9 |
6.2 |
16-Jan-13 |
14:00 |
2.3 |
5.7 |
16-Jan-13 |
15:00 |
0.8 |
2.7 |
16-Jan-13 |
16:00 |
2.2 |
1.6 |
16-Jan-13 |
17:00 |
2.2 |
2 |
16-Jan-13 |
18:00 |
0.7 |
1.8 |
16-Jan-13 |
19:00 |
0.2 |
0.9 |
16-Jan-13 |
20:00 |
0.1 |
0.3 |
16-Jan-13 |
21:00 |
0 |
0.1 |
16-Jan-13 |
22:00 |
0 |
0 |
16-Jan-13 |
23:00 |
0 |
0 |
17-Jan-13 |
0:00 |
0 |
0 |
Table 3 Calibrated & validated result of HEC-HMS run model
Final calibrations of is given below of sub-basin 1 and sub-basin 2, reach1 and reach 2 Figure 9.
Run model parameters
Sub-basin 1 area: 240 km2
Sub-basin 2 area: 239.5 km2
Muskingum; K (HR) = 1
Muskingum; X= 0.2
Lag time= 60
Curve Number=45
Impervious= 0% Figure 10–14.
The simulated and observed peak discharges occurred on the same day, with a maximum time difference of one hour, which is acceptable for flood predictions, as can be seen in the graphs above.
The model correctly predicted peak discharge based on known historical flood data, as seen in the findings above. The flood's size and timing were both fairly accurate. This demonstrates that the HEC-HMS is suitable for the catchments in question. We can draw from the findings that the complexity of the model structure has no influence on its applicability or efficiency. Despite its simplicity, the HEC-HMS framework is an effective flood forecasting tool. More testing of the HEC-HMS should be encouraged in order to determine its suitability for Pakistani catchments.
The author would like to thank Pakistan metrological department (PMD) and water and Power Development Authority (WAPDA) for providing the discharge and rainfall data used in this study.
The authors declares no conflict of interest.
©2022 Nadeem, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.