Next Article in Journal
Long Term Historic Changes in the Flow of Lesser Zab River, Iraq
Next Article in Special Issue
Improving Urban Runoff in Multi-Basin Hydrological Simulation by the HYPE Model Using EEA Urban Atlas: A Case Study in the Sege River Basin, Sweden
Previous Article in Journal
Hydrologic Trends in the Upper Nueces River Basin of Texas—Implications for Water Resource Management and Ecological Health
Previous Article in Special Issue
Quantifying Thermal Characteristics of Stormwater through Low Impact Development Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of HEC-HMS Model for Flow Simulation in the Lake Tana Basin: The Case of Gilgel Abay Catchment, Upper Blue Nile Basin, Ethiopia

1
Bahir Dar Institute of Technology, Bahir Dar University, P.O. Box 26 Bahir Dar, Ethiopia
2
Faculty of Agricultural and Environmental Sciences, University of Rostock, Satower 48, 18051 Rostock, Germany
*
Author to whom correspondence should be addressed.
Submission received: 16 January 2019 / Revised: 19 February 2019 / Accepted: 20 February 2019 / Published: 10 March 2019

Abstract

:
Understanding the complex relationships between rainfall and runoff processes is necessary for the proper estimation of the quantity of runoff generated in a watershed. The surface runoff was simulated using the Hydrologic Modelling System (HEC-HMS) for the Gilgel Abay Catchment (1609 km2), Upper Blue Nile Basin, Ethiopia. The catchment was delineated and its properties were extracted from a 30 m × 30 m Digital Elevation Model (DEM) of the Lake Tana Basin. The meteorological model was developed within HEC-HMS from rainfall data and the control specifications defined the period and time step of the simulation run. To account for the loss, runoff estimation, and flow routing, Soil Conservation Service Curve Number (SCS-CN), Soil Conservation Service Unit Hydrograph (SCS-UH) and Muskingum methods were used respectively. The rainfall-runoff simulation was conducted using six extreme daily time series events. Initial results showed that there is a clear difference between the observed and simulated peak flows and the total volume. Thereafter, a model calibration with an optimization method and sensitivity analysis was carried out. The result of the sensitivity analysis showed that the curve number is the sensitive parameter. In addition, the model validation results showed a reasonable difference in peak flow (Relative Error in peak, REP = 1.49%) and total volume (Relative Error in volume, REV = 2.38%). The comparison of the observed and simulated hydrographs and the model performance (NSE = 0.884) and their correlation (R2 = 0.925) showed that the model is appropriate for hydrological simulations in the Gilgel Abay Catchment.

1. Introduction

Adequate knowledge of rainfall-runoff processes is vital to estimate the amount of runoff produced within a given catchment. Knowing the amount of runoff within a given catchment is important for sustainable water resources project planning and management. The activities to estimate runoff volumes and flood peaks can be easily simplified by adopting a modelling concept and by understanding rainfall partitioning and the principal factors triggering runoff [1]. The type of the modelling approach normally depends on the purpose, data availability and ease of use [2].
Rainfall-runoff models are often used as a tool for a wide range of tasks, such as the modelling of flood events, the monitoring of water levels during different water conditions or the prediction of floods [3]. Generally, hydrological models can be classified as stochastic and deterministic. The stochastic models will produce outputs that have partial randomness but the deterministic models on the other hand do not give randomness. Cunderlik [4] further classified deterministic hydrologic models into three major categories. Firstly, the lumped model, which assesses the catchment response simply at the outlet without obviously counting for an individual sub-basins response. Secondly, the semi-distributed model, which is partly permitted to change in space with a division of the catchment into a number of sub-basins. The third type of model is the distributed model, which permits its parameters to change in place at a resolution normally chosen by the client.
Distributed hydrological models such as the European Hydrological System Model (MIKE-SHE) and Modular Modeling System (MMS), and semi-distributed models like the Hydrological Engineering Center Hydrological Modelling System, Soil and Water Assessment Tool (SWAT),Topography Based Hydrological Model (TOPMODEL), Hydrologiska Byråns Vattenbalansavdelning (HBV) and Hydrological Simulation Program-Fortran (HSPF) are developed for a runoff estimation based on the data availability and complexity of the hydrological systems. Flood modelling has been greatly improved in recent years with the advent of geographic information systems, ”radar-based” rainfall estimation using next-generation radar (NEXRAD), high-resolution digital elevation models, and distributed hydrological models [5]. More recently, flood modelling has been further improved with the advent of service-oriented architecture and numerical weather predictions [6]. In the case of flood predictions, rainfall-runoff models are very practical because they are even useful in the catchments with a limited amount of input data.
The Hydrologic Modelling System HEC-HMS, which is a hydrologic modelling software developed by the US Army Corps of Engineers Hydrologic Engineering Center (HEC) is an integrated modelling tool for all hydrologic processes of dendritic watershed systems. It consists of different component processes for rainfall loss, direct runoff, and routing. HEC-HMS has become very popular and been adopted in many hydrological studies because of its ability in the simulation of runoff both in short and longtime events, its simplicity to operate, and use of common methods [7]. Hydrographs developed by HEC-HMS either directly or in conjunction with other software’s are used for studies of urban drainage, water availability, future urbanization impact, flow forecasting, flood damage reduction, floodplain regulation, and systems operation [8].
Previous studies on HEC-HMS proved its ability to simulate and forecast streamflow based on different datasets and catchment types [9]. Most of these studies clearly indicated that the results of the model simulation were location specific, in that different combinations of a model set containing the loss methods, runoff transform methods, and baseflow separation techniques were found to respond variably [10]. In this regard, a study conducted by Al-abed et al. [11] on the Al-Zarqa Basin in Jordan using a monthly streamflow showed that the HEC-HMS model can give more acceptable results than other models. Other researchers who used the model in different areas also have gotten satisfactory results. Radmanesh et al. [12] calibrated and validated the HEC-HMS model in the Yellow River watershed in Southwestern Iran. Sardoii et al. [13] used the HEC-HMS and Geographic Information System GIS to simulate the rainfall-runoff process in the Amirkabir watershed. Fallah [14] applied the HEC-HMS model to simulate the river flow in the Mohammadabad watershed located in the North of Iran. Yusop et al. [15] used the HEC-HMS for predicting the runoff in the Oil Palm Catchment. Yener et al. [16] simulated the rainfall-runoff process by dividing the basin into three sub-basins: Kirazdere, Kazandere, and Serindere, whereby each sub-basin was modeled with its own parameters.
Although the HEC-HMS model has been tested and calibrated at a global scale, little effort has been made in the context of Ethiopian catchments [17]. The Upper Blue Nile (Abay) Basin, which covers an area of 199,000 km2, has extensive agricultural practices and unpredicted flood events. It has been seriously affected by land degradation and flood problems, especially in the lowlands [18]. The Gilgel Abay Catchment which is part of this basin faces the same problem. Therefore, an accurate estimation of the peak flow and total volume from the respective rainfall events is critically important to implement appropriate soil and water conservation, erosion control, and flood protection measures in time.
Despite the different modelling activities that are practiced in the basin, the HEC-HMS model was not tested, calibrated, and validated in the Gilgel Abay Catchment. Therefore, the objectives of this study were: To develop the rainfall-runoff model and assess the runoff potential of the catchment; and to calibrate and validate the model and fix the corresponding calibrated values for future hydrological investigations in similar catchments within the basin and the nearby basins.

2. Methodology

2.1. Study Area Description

The Gilgel Abay Catchment has an area of 1609 km2 at the gauging station near Piccolo, with elevations ranging between 1800 m and 3530 m above mean sea level (Figure 1). Geographically it is bounded by the northern latitudes from 10°56′ to 11°23′ and eastern longitudes from 36°48′ to 37°24′. Soils are characterized by clay, clay loam, and silt loam textures, with each texture sharing similar proportions of the catchment area [19]. The majority of the catchment is a basalt plateau with gentle slopes, while the southern part has a rugged topography. The climate is semi-humid with a mean annual rainfall of 1300 mm (1990–2012), more than 70% of which falls in the summer monsoon season (June to September).

2.2. Data

Input data as an important part of rainfall-runoff modelling can be categorized into hydro-meteorological (rainfall and stream flow) and physiographic (digital elevation model, land use/cover and soil type) databases.
The data required to build an HEC-HMS model are elevation, land use/cover, percent impervious area, soil, and hydrograph information. These datasets were used to determine stream/sub-basin characteristics and hydrologic parameter estimations.
The basic information on physiographic characteristics of the catchment (elevation and slope) was extracted from a 30 m × 30 m resolution digital elevation model of the Lake Tana Basin which was downloaded from the web for free. The Geographic Information System (ArcGIS 10.2) interface of the HEC-GeoHMS model was used to further process the DEM data. Terrain pre-processing and basin processing tools were used to generate the characteristic parameters of the basin and input files for HEC-HMS.
The calibration and validation of hydrological models are the fundamental operations that should be considered in many modelling studies which require a series of runoff data [20]. For this study the streamflow data (1990–2014), land use/cover, and soil data were collected from the Hydrology and GIS Departments of the Ministry of Water, Irrigation and Electricity of Ethiopia whereas precipitation data (1990–2014) were collected from the National Meteorological Agency of Ethiopia. Table 1 and Table 2 below indicate the events selected for the calibration and validation and physiographic characteristics of each sub-catchment respectively.

2.3. Models

2.3.1. HEC-HMS Model

The HEC-HMS model is physically based and conceptually semi-distributed model designed to simulate rainfall-runoff processes in a wide range of geographic areas, from large river basin water supplies and flood hydrology to small urban and natural watershed runoffs. It easily operates huge tasks in relation to hydrological studies, including losses, runoff transform, open channel routing, analysis of meteorological data, rainfall-runoff simulation, and parameter estimation [21,22]. Moreover, the HEC-HMS uses separate models that compute runoff volume, models of direct runoff, and models of baseflow. It has nine different loss methods, some of which are designed primarily for simulating events, while others are intended for continuous simulation. It also has seven different transformation methods. For example the Snyder Unit Hydrograph Yilma and Moges [17] and Clark Unit Hydrograph Banitt [23] methods have been applied successfully to simulate long-term stream flows elsewhere.
HEC-GeoHMS was used to create a basin model and a meteorological model, and to control specifications before running the model. The basin model and the basin features were taken as the background map file and imported to HEC-HMS 4.2. Since we do not have observation stations in each sub-basin, the precipitation values were estimated by the most commonly used Thiessen Polygon method and weights were worked out in HEC-GeoHMS software. The method was selected because of its easy application; accuracy depends on the sampling density. The boundaries are often oddly shaped, as transitions between polygons are often abrupt and continuous variables are often not well represented. During calibration, a daily simulation time step was used in estimating the parameters of the selected methods. Each model run combines a basin model, meteorological model, and control specifications with run options to obtain results. It contains four main components: (1) An analytical model to calculate overland flow runoff as well as channel routing, (2) an advanced graphical user interface illustrating hydrologic system components with interactive features, (3) a system for storing and managing data—specifically large, time variable data sets—and (4) a means for displaying and reporting model outputs [24].
The Soil Conservation Service Curve Number, SCS Unit hydrograph, and Muskingum routing methods were selected for each component of the runoff process as runoff depth, direct runoff, and channel routing respectively. These methods were chosen on the basis of applicability and limitations of each method, availability of data, suitability for the same hydrologic condition, stability, wide acceptability, and well established researcher recommendations. The baseflow was not considered because in short duration event modelling and for small sub-basins the contribution of base flow is insignificant for flood flows [25,26].
(i) Loss Model
The loss models in HEC-HMS normally calculate the runoff volume by computing the volume of water that is intercepted, infiltrated, stored, evaporated, or transpired and subtracting it from the precipitation. In this study, the Soil Conservation Service Curve Number loss method was selected to estimate direct runoff from a specific or design rainfall [27]. It has several advantages over other methods in that: It is a simple conceptual method for the estimation of the direct runoff amount from a storm rainfall event, and is well supported by empirical data; it relies only on the curve number, which is a function of the soil type and land use/cover that are the major runoff-producing watershed characteristics. However, there are several problems associated with the SCS-CN method. For example, it does not account for rainfall intensity and temporal variation of rainfall as well as for average ground slope [28]. Despite the above problems, the SCS-CN loss method was chosen for HEC-HMS analysis in this study because of the following reasons: It is commonly used in different environments and provides better results compared to initial and constant loss rate method [13]. Its calculation is made easier by the fact that only a few variables need to be estimated based on hydrologic soil group, land use and slope maps. And despite its simplicity, it yields results that are as good as those of complex models [29].
The SCS-CN model assumes that the accumulated rainfall-excess depends upon the cumulative precipitation, soil type, land use and the previous moisture conditions as estimated in the following relationship [21].
P e = ( P I a ) 2 P I a + S
where Pe is the accumulated precipitation excess at time t (mm); P is the accumulated rainfall depth at time t (mm); Ia is the initial abstraction (initial loss) (mm) = αS, α is 0.2 as a standard; and S is the potential maximum retention (mm), a measure of the ability of a watershed to abstract and retain storm precipitation.
In the curve number method, the runoff is directly proportional to the precipitation with an assumption that the runoff is produced after the initial abstraction of 20% of the potential maximum storage [30].
The maximum retention, S, and watershed characteristics are related through an intermediate dimensionless parameter, the curve number (CN) as:
S   =   25400 254 × C N C N
where CN is the SCS curve number used to represent the combined effects of the primary characteristics of the catchment area, including soil type, land use, and the previous moisture condition. The CN values range from 100 (water bodies) to approximately 30 for permeable soils with high infiltration rates [21].
(ii) The Transform Model
The transform prediction models in HEC-HMS simulate the process of the direct runoff of excess precipitation on the watershed, and they transform the precipitation excess in point runoff. In this study, the Soil Conservation Service Unit Hydrograph model was chosen to transform excess precipitation into runoff. It is a parametric model based on the average Unit Hydrograph (UH) derived from gauged rainfall and runoff data of a large number of small agricultural watersheds throughout the United States. The SCS proposed the Unit Hydrograph (UH) model, and it is included in the HEC-HMS program. The lag time (Tlag) is the only input for this method. It is the time from the center of mass of excess rainfall to the hydrograph peak and is calculated for each watershed based on the time of concentration Tc, as:
T l a g   =   0.6 T c
where Tlag and Tc are in minute.
The time of concentration can be estimated based on basin characteristics including topography and the length of the reach by Kirpich’s formula [31].
T c =   0.0078 × ( L 0.77 S 0.385 )
where L is the reach length in feet, and S is the slope in (ft/ft).
(iii) Routing Model
As the flood runoff travels through the channel reach, it becomes attenuated due to channel storage effects. The routing models available in HEC-HMS account for this attenuation. The Muskingum method, which was developed by McCarthy [32], is a popular lumped flow routing technique which was selected for this study.
The Muskingum routing method is a simple approximate method to calculate the outflow hydrograph at the downstream end of the channel reach from the inflow hydrograph at the upstream end. Among many models used for flood routing in rivers, it is a straightforward hydrological flood routing technique used in natural channels [33], and it has been extensively applied in river engineering practice since its introduction in the 1930s [34]. In this model calibration, two parameters are needed; travel time (K) of the flood wave through routing reach; and dimensionless weight (X) which corresponds to the attenuation of the flood wave as it moves through the reach. The routing parameters in the models are usually derived through calibration using measured discharge hydrographs [35].
S =   K   [ XI + ( 1 X )   Q ] ,
in which the prism storage in the reach is KQ, where K is a proportionality coefficient, and the volume of the wedge storage is equal to KX (I − Q), where X is a weighting factor having a range of 0 ≤ X ≤ 0.5.

2.3.2. HEC-GeoHMS Model

Background map files, basin model files, meteorological model files, and a grid cell parameter are created by HEC-GeoHMS which are input to a hydrological model HEC-HMS. The following steps adopted from Merwade [36] were used to extract the basin model: (i) Data collection such as DEM, soil and land use/cover; (ii) data assembly; (iii) terrain preprocessing. This latter part of the HEC-GeoHMS processing step helps to drive the drainage network from the input DEM processing, such as flow direction, flow accumulation, stream definition and segmentation, catchment polygon processing, drainage line processing and adjoint catchment processing; (iv) Hydrologic Modelling System (HMS) Project Setup: The input files for the HEC-HMS project were developed using the HMS project set up menu in HEC-GeoHMS, this helps to copy all the terrain preprocessing data to the HEC-HMS project; (v) basin processing; (vi) the extraction of basin and stream characteristics such as length, upstream and downstream elevations, and river slopes. It can help to extract the physical characteristics of sub-basins, such as longest flow path, basin centroid, centroid elevation, centroidal longest flow lengths, and basin slopes; (vii) the estimation of hydrologic parameters, such as the curve numbers, the percentage impervious area and time of concentration initial values are estimated using HEC-GeoHMS model processing; (viii) the creation of HMS model files, such as background shapefile, the basin model, meteorological model file, and a project file.
HEC-GeoHMS extension in ArcGIS was used to delineate the boundary of the area by considering the geographical reference point of the hydrological gauging station as the outlet point of the Gilgel Abay Catchment. Further HEC-GeoHMS processing of the DEM produced four sub-basins, five routing reaches, five junctions, and the major physiographic characteristics of the catchment as shown in Figure 1 and Table 2.

2.4. Model Calibration and Validation

A sensitivity analysis is usually undertaken in most modelling studies [15,20]. It is a necessary process to identify the key parameters and parameter precision required for calibration. The most fundamental sensitivity analysis technique utilizes partial differentiation, whereas the simplest method involves perturbing parameter values one at a time [37]. In this study, the sensitivity analysis was performed by changing the parameter values in the range of ±25% with 5% intervals. Zelelew and Melesse [10] and Majidi and Shahedi [38] conduct a sensitivity analysis by changing the parameter values in the range of ±30% with 5% intervals until the best fit between the observed and simulated parameters in the given hydrograph was obtained. Another study by Sardoii et al. [13] used a parameter change of ±15% with a 5% interval in order to analyze the sensitivity of the calibration parameters in the Amirkabir Dam Watershed in Iran. The sensitivity parameters were then selected based on their effect on peak discharge and total volume. Thereafter, the model calibration was done with the Univariate Gradient optimization package and Peak-Weighted Root Mean Square Error (PWRMS) objective function because of their simplicity and performance (Deng et al., 2010, as cited by Jin et al. [26]. In this study, ten flood events (six for calibration and four for validation) recorded in the period from 1990–2014 were selected for calibration and validation purposes (Table 1). Watershed parameters such as curve number, infiltration coefficients, time of concentration, and baseflow may need modification to produce the best fit between model and observations. For validation, the simulated data as predicted by the model must be computed and compared with the observed data and statistical tests of error functions must be carried on.
Statistical evaluation techniques like relative bias error functions proposed by Najim [39], Nash–Sutcliffe Efficiency (NSE) by Nash and Sutcliffe [40], which is equal to one minus the ratio of the mean square error and standard deviation of measured data time series [41], and coefficient of determination (R2) as described in Neter et al. [42] were applied to evaluate the performance of the model and the selected loss and transform methods.

3. Results and Discussion

3.1. Calibration

In analyzing a flood event, the most important aspect of the hydrograph is the peak flow, because the peak flow corresponds to the maximum downstream flooding. In contrast, peaks that are significantly less than the maximum may correspond to increased water levels, but not necessarily a flood event. The results of the hydrological model in this study showed a reasonable fit between the model and observations after optimization; the hydrograph shape and timing of peaks matched well, although the model tended to overestimate the runoff before optimization (Table 3). In the majority of events, the hydrograph shape was accurately reproduced in the model output. However, the volume of runoff was overestimated in events 3, 4, and 6 and underestimated for the remaining Events (1, 2, and 5).The calibration of the model improved the results greatly by decreasing and increasing the overestimated and underestimated volume of the runoff respectively.
The modelling results of peak discharge, total volume, and their relative errors with respect to the observed data, the Nash Sutcliffe Efficiency and coefficient of determination values during calibration are mentioned in Table 3 below. Figure 2 (Events 1, 3, and 6) and 3 (Events 1, 4, and 6) show the simulated and observed hydrographs and their correlation values respectively. Only sample calibrated events are presented in the Figure 2 and Figure 3.
As mentioned in Section 2.4, six extreme events were selected for calibration from a 25 year rainfall series data. The selection of rainfall-runoff events is a critical step for event hydrologic modelling and model calibration/validation [9]. The selection depends on many factors, such as rainfall characteristics, watershed properties and availability and completeness of rainfall and stream monitoring data. Flood analysis tends to use large events. However, considering that the scale of the events reflect the range of the applicability and performance of the model [43]. The events of this study were selected by comparing the runoffs produced for the respective rainfalls (a rainfall event with a higher amount of runoff is selected). Table 4 summarizes the initial and optimized parameters of the catchment and Table 5 shows the optimized parameter values of each event during calibration.
In the rainfall-runoff modelling of the study area, the parameters curve number scale factor and initial abstraction scale factor show a high variability and parameters of routing (K and X) are constant in the calibrated events. These results demonstrate the impact of basin topography, soil type, land use and other catchment characteristics.
According to the result of the flood volume relative error, REv, as can be seen in Table 3, the error is relatively small for all flood events with absolute values less than 10%. This may be due to the sub-humid attributes: The flood runoff is dominated by infiltration excess under the condition of short but intensive rainfall, or by the combination of infiltration and saturation excess under the condition of long-lasting rainfall of various intensities. From Figure 2 it is possible to conclude that runoff in all events is caused by the combination of infiltration and saturation excess due to strong intensity rainfall. The calculated values of the percent error both in total volume and peak flow between simulated and observed values in all simulations before optimization was very high, which falls in the range between the absolute value of 8.2 to 50.64% for the total volume and 8.9 to 96.6% for the peak flow, with mean values of 23.96% and 48.07%, respectively. Considering this result, a sensitivity analysis was done to identify the most sensitive parameter for the loss, transform, and routing methods. It was found that the curve number, initial abstraction, and lag time were more sensitive, less sensitive and insensitive parameters of the model, respectively. Nearly a similar result was obtained for the peak flow by Zelelew and Melesse [10] which ranges from 37.5 to 60.2%, with a mean value of 50.1% for a model set containing the SCS loss method with Clark’s unit hydrograph transform method combination. After optimization, the values in this study were reduced to 0.07 to 5.21% for the total volume and 0.28 to 4.51% for the peak flow, with mean values of 2.76% and 1.91%, respectively (Table 3). The result is very good according to Najim et al. [39] and Sabzevari et al. [44], who recommended that the acceptable ranges of relative percent errors between the observed and simulated values should be below ±20%. The study by Cheng et al. [45] also indicated that the runoff model is considered good if the percentage error of the runoff volume is less than 20% according to the criteria for flood forecasting in China. In this statistical evaluation criterion, the positive values indicate model underestimation bias, and the negative values indicate model overestimation bias [46].
Figure 3 and Table 3 also show a relatively close agreement between the observed and simulated peak flow values at the period of calibration (R2 = 0.842) after optimization. More specifically, Event 6 had the lowest REp (0.28) but an unexpected lower NSE (0.603) and R2 (0.689) values. These might be due to high peak discharges as a result of more extreme rainfall in 20–21 of Julyy 2008 and the modelling effect because the model simulation showed that the simulated peak is 24 h earlier than the observation. According to the classification range mentioned in Zou et al. [47] the mean correlation coefficient obtained in this study can be considered as strong (>0.8). Considering the Nash–Sutcliffe Efficiency (NSE) criteria, better results were obtained between the simulated and observed values, with a mean NSE value of 74.6% (Table 3). Therefore, the model performs well. The model simulation can be judged as satisfactory if Nash–Sutcliffe Efficiency is greater than 50%, good if it is greater than 65%, and very good if it is greater than 75% [48]. Overall, in this study the three statistical evaluation criteria with mean values of REP = 1.91%, REV = 2.76%, NSE = 0.746, R2 =0.842 showed good simulation between the estimated and observed values.

3.2. Validation

In this study, the validation of the simulated results (Table 6) of peak discharge, total volume, and their relative errors with respect to the observed data, the Nash–Sutcliffe Efficiency (NSE) and the coefficient of determination values are listed. Figure 4 (Events 8 and 9) and 5 (Events 7 and 9) showed the simulated and observed hydrographs and their correlation respectively. Only the sample validation events are presented in Figure 4 and Figure 5.
Four extreme events from 25 years rainfall series data, which are different from the events used for calibration, were selected for validation purpose (Table 1). The selection of the events was similar to the calibration events (see Section 3.1). The calibration model simulated the values of different runoff hydrographs, such as the peak flow and total volume by considering the mean of the optimized parameters shown in Table 5. The result showed that the simulated values are close to the observed ones for all the events (Table 6). The mean relative percent error between the observed and simulated values of peak flow and total volume is calculated to be 1.49% and 2.38%, respectively.
Figure 5 and Table 6 also show relatively close agreement between the observed and simulated peak flow values during validation (R2 = 0.925). Considering the Nash–Sutcliffe Efficiency (NSE) criteria, better results were obtained between the simulated and observed values, with a mean NSE value of 88.4%. According to Moriasi et al. [48], the model simulation can be judged as very good. Generally, in this study, the three statistical evaluation criteria with mean values of REP = 1.49%, REV = 2.38%, NSE = 0.884, and R2 =0925 show good simulation between the estimated and observed values. The validation results (Table 6) were better than the calibration (Table 3) and the statistical evaluation criteria also performed better. A similar result was obtained by Gharib et al. [49] when they used the Modclark model to evaluate the simulation of the rainfall-runoff in the Tangrah watershed, Iran.
The overall simulation results of events selected for validation show us that there is a possibility of getting a good performance of the HEC-HMS model in predicting the peak flow and total volume in Gilgel Abay Catchment by using SCS-CN and SCS UH methods. Studies by Yilma and Moges [17] and Majidi and Shahedi [38] also obtained credible simulation results of the peak flow and total volume in the HEC-HMS model containing a combination of different loss and transform methods.
From the results of the statistical evaluation criteria a good performance of the HEC-HMS model was obtained in simulating the runoff volume and peak flow with SCS loss and transform methods. However, in order to increase the performance of the model in simulating the runoff with these combinations, it is better to develop the local and regional curve number of the catchment instead of using the tabulated standardized CN values and to identify regional valid unit hydrographs. Zelelew [50] and Hawkins [51] also recommended the same.
Generally, the results of this study provide basic information on the extent of total volume and peak flow generated in the catchment from the respective rainfall events that in turn are useful for the planning, designing and management of different water resources activities. The results of these findings can also be used as a basic input data for further hydrological investigation for adjacent ungauged catchments by transferring the required calibrated parameters through regionalization and other techniques as explained by Jin et al. [52].

4. Conclusions

This study demonstrated that the HEC-HMS hydrological model is adaptable to tropical conditions. The model is based on the hydrological characteristics, topography, soil type and land use of the study area. Basin characteristics and initial values were analyzed using HEC-GeoHMS in ArcGIS in order to start the model calibration. The rainfall-runoff simulation was conducted using extreme rainfall events and initial results showed that there is a clear difference between observed and simulated peak flows and volumes in the study area. Therefore, the model calibration was conducted to optimize the parameters. To know the most influential parameter in the simulation a sensitivity analysis was carried out and the results showed that the curve number is the most sensitive parameter followed by the lag time. After optimization the peak flow and total volume of all events are very close to the observations with a very small error in peak and volume (<10%). Overall performance of the HEC-HMS model was very good in terms of relative error functions, Nash-Sutcliffe Efficiency, and coefficient of determination based on the selected loss, transform and flow routing methods. So, it is possible to suggest that the calibrated parameters can be further used to other catchments in the basin and nearby basins.
Note that the verification results are not worse than the results obtained after calibration. This result is surprising since the verification results are more or less worse in general. Therefore, we conclude that there are relatively unique input–output relations and that the runoff formation is dominated by the only mechanism.
Finally, the methodologies developed in this research can also be applied in other ungauged catchments and regions with similar characteristics. We suggest further studies in the study catchment to generate more detailed information for modelling work by reviewing the efficiency of meteorological and hydrometric measurements (poor in quantity and quality during this study) in order to establish the optimal number of stations and their adequate distribution in the catchment.

Author Contributions

Conceptualization, B.G.T.; data curation, B.G.T.; formal analysis, B.G.T.; investigation, B.G.T.; methodology, B.G.T.; supervision, M.A.B. and K.M.; validation, B.G.T.; writing–original draft, B.G.T.; writing–review & editing, B.G.T., M.A.B. and K.M.

Funding

This research received no external funding

Acknowledgments

The authors would like to acknowledge the Ethiopian Ministry of Water, Irrigation and Electricity for proving the flow, soil, and land use/cover data and the National Meteorology Agency of Ethiopia for providing the required rainfall data. We also want to acknowledge Mr. Aderaw Genetu Tassew for his professional assistance in language and grammar editing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, W.W.; Fu, B.J.; Meng, Q.H.; Zhang, Q.J.; Zhang, Y.H. Effects of land-use pattern change on rainfall-runoff and runoff-sediment relations: A case study in Zichang watershed of the Loess Plateau of China. J. Environ. Sci. 2004, 16, 436–442. [Google Scholar]
  2. Beven, K.J. Rainfall-Runoff Modelling: The Primer; John Wiley & Sons: Chichester, UK; Wiley-Blackwell: Hoboken, NJ, USA, 2012. [Google Scholar]
  3. Jia, Y.; Zhao, H.; Niu, C.; Jiang, Y.; Gan, H.; Xing, Z.; Zhao, Z.A. WebGIS-based system for rainfall-runoff prediction and real-time water resources assessment for Beijing. Comput. Geosci. 2009, 35, 1517–1528. [Google Scholar] [CrossRef]
  4. Cunderlik, M.J. Hydrologic Model Selection Fort the CFCAS Project: Assessment of Water Resources Risk and Vulnerability to Changing Climatic Conditions; Department of Civil and Environmental Engineering, The University of Western Ontario: London, ON, Canada, 2003. [Google Scholar]
  5. Bedient, P.B.; Holder, A.; Benavides, J.A.; Vieux, B.E. Radar-based flood warning system applied to Tropical Storm Allison. J. Hydrol. Eng. 2003, 8, 308–318. [Google Scholar] [CrossRef]
  6. Shi, H.; Li, T.; Liu, R.; Chen, J.; Li, J.; Zhang, A.; Wang, G. A service-oriented architecture for ensemble flood forecast from numerical weather prediction. J. Hydrol. 2015, 527, 933–942. [Google Scholar] [CrossRef]
  7. Halwatura, D.; Najim, M. Application of the HEC-HMS model for runoff simulation in a tropical catchment. Environ. Modell. Softw. 2013, 46, 155–162. [Google Scholar] [CrossRef]
  8. US Army Corps of Engineers. Hydrologic Modeling System, HEC-HMS. Quick Start Guide; US Army Corps of Engineers, Institute for Water Resources, Hydrologic Engineering Center: Davis, CA, USA, 2015.
  9. Chu, X.; Steinman, A. Event and continuous hydrologic modeling with HEC-HMS. J. Irrig. Drain. E-Asce. 2009, 135, 119–124. [Google Scholar] [CrossRef]
  10. Zelelew, D.; Melesse, A. Applicability of a Spatially Semi-Distributed Hydrological Model for Watershed Scale Runoff Estimation in Northwest Ethiopia. Water 2018, 10, 923. [Google Scholar] [CrossRef]
  11. Al-Abed, N.; Abdulla, F.; Khyarah, A.A. GIS-hydrological models for managing water resources in the Zarqa River basin. Environ. Geol. 2005, 47, 405–411. [Google Scholar] [CrossRef]
  12. Radmanesh, F.; Hemat, J.P.; Behnia, A.; Khond, A.; Mohamad, B.A. Calibration and assessment of HEC-HMS model in Roodzard watershed. In Proceedings of the 17th international conference of river engineering, Ahvaz, Iran, 2006. [Google Scholar]
  13. Sardoii, E.R.; Rostami, N.; Sigaroudi, S.K.; Taheri, S. Calibration of loss estimation methods in HEC-HMS for simulation of surface runoff (Case Study: Amirkabir Dam Watershed, Iran). Adv. Environ. Biol. 2012, 6, 343–348. [Google Scholar]
  14. Fallah, S. Simulation of maximum peak discharge in river tributaries using HEC-HMS model (Case study: Moha mmadabad watershed, Golestan province). Master’s Thesis, University of Gorgan, Gorgan, Iran, 2001. [Google Scholar]
  15. Yusop, Z.; Chan, C.; Katimon, A. Runoff characteristics and application of HEC-HMS for modelling stormflow hydrograph in an oil palm catchment. Water Sci. Technol. 2007, 56, 41–48. [Google Scholar] [CrossRef] [PubMed]
  16. Yener, M.K.; Sorman, A.U.; Sorman, A.A.; Sensoy, A.; Gezgin, T. Modeling studies with HEC-HMS and runoff scenarios in Yuvacik Basin, Turkiye. Int. Congr. River Basin Manag. 2007, 4, 621–634. [Google Scholar]
  17. Yilma, H.M.; Moges, S.A. Application of semi-distributed conceptual hydrological model for flow forecasting on upland catchments of Blue Nile River Basin, a case study of Gilgel Abbay catchment. Catchment Lake Res. 2007, 6, 1–200. [Google Scholar]
  18. Merrey, D.J.; Gebreselassie, T. Promoting Improved Rainwater and Land Management in the Blue Nile (Abay) Basin of Ethiopia: Annexes. ILRI: Nairobi, Kenya; Department for International Development: London, UK, 2011.
  19. Bitew, M.; Gebremichael, M. Assessment of satellite rainfall products for streamflow simulation in medium watersheds of the Ethiopian highlands. Hydrol. Earth Syst. Sci. 2011, 15, 1147–1155. [Google Scholar] [CrossRef]
  20. Azam, M.; San Kim, H.; Maeng, S.J. Development of flood alert application in Mushim stream watershed Korea. Int. J. Disast. Risk Re. 2017, 21, 11–26. [Google Scholar] [CrossRef]
  21. Feldman, A. Hydrologic Modeling System HEC-HMS Technical Reference Manua; US Army Corps of Engineers. Hydrologic Engineering Center: Davis, CA, USA, 2000.
  22. US Army Corps of Engineers. Hydrologic Modeling System (HEC-HMS) Application Guide Version 3.1.0; Institute for Water Resources: Davis, CA, USA, 2008.
  23. Banitt, A. Simulating a century of hydrographs e Mark Twain reservoir. In Proceedings of the 2nd Joint Federal Interagency Conference, Las Vegas, NV, USA, 27 June–1 Julyy 2010. [Google Scholar]
  24. Bajwa, H.; Tim, U. Toward immersive virtual environments for GIS-based Floodplain modeling and Visualization. In Proceedings of the 22nd ESRI User Conference, San Diego, CA, USA, 8–12 Julyy 2002. [Google Scholar]
  25. Kamali, B.; Mousavi, S. Automatic calibration of HEC-HMS model using multi-objective fuzzy optimal models. CEIJ 2014, 47, 1–12. [Google Scholar]
  26. Jin, H.; Liang, R.; Wang, Y.; Tumula, P. Flood-runoff in semi-arid and sub-humid regions, a case study: A simulation of jianghe watershed in northern China. Water 2015, 7, 5155–5172. [Google Scholar] [CrossRef]
  27. Environmental and Water Resources Instit. Curve Number Hydrology: State of the Practice; Hawkins, R.H., Ward, T.J., Woodward, D.E., Van Mullem, J.A., Eds.; American Society of Civil Engineers: Reston, VA, USA, 2009. [Google Scholar]
  28. Mockus, V. Estimation of Total (and Peak Rates of) Surface Runoff for Individual Storms. Exhibit A in Appendix B, Interim Survey Report, Grand (Neosho) River Watershed; USDA: Washington, DC, USA, 1949.
  29. Lastra, J.; Fernández, E.; Díez-Herrero, A.; Marquínez, J. Flood hazard delineation combining geomorphological and hydrological methods: An example in the Northern Iberian Peninsula. Nat. Hazards 2008, 45, 277–293. [Google Scholar] [CrossRef]
  30. Heshmatpoor, A. Identification runoff source area in tropical watershed. In Proceedings of the Postgraduate Qolloquium Semester, Kuala Lumpur, Malaysia, 26–29 October 2009. [Google Scholar]
  31. Kirpich, Z. Time of concentration of small agricultural watersheds. Civ. Eng. 1940, 10, 362. [Google Scholar]
  32. McCarthy, G.T. The unit hydrograph and flood routing. In Proceedings of the Conference of North Atlantic Division, Washington, WA, USA, 1938. [Google Scholar]
  33. Shaw, E.M. Hydrology in Practice; T.J. International Ltd.: Cornwall, UK, 1994. [Google Scholar]
  34. Tewolde, M.H.; Smithers, J. Flood routing in ungauged catchments using Muskingum methods. Water SA 2006, 32, 379–388. [Google Scholar] [CrossRef]
  35. Birkhead, A.; James, C. Muskingum river routing with dynamic bank storage. J. Hydrol. 2002, 264, 113–132. [Google Scholar] [CrossRef]
  36. Merwade, V. Watershed and Stream Network Delineation Using ArcHydro Tools; School of Civil Engineering, Purdue University: West Lafayette, IN, USA, 2012; pp. 1–7. [Google Scholar]
  37. Hamby, D. A review of techniques for parameter sensitivity analysis of environmental models. Environ. Monit. Assess. 1994, 32, 135–154. [Google Scholar] [CrossRef] [PubMed]
  38. Majidi, A.; Shahedi, K. Simulation of rainfall-runoff process using Green-Ampt method and HEC-HMS model (Case study: Abnama Watershed, Iran). J. Hydraul. Eng. 2012, 1, 5–9. [Google Scholar]
  39. Najim, M.M.M.; Babelb, M.S.; Loofb, R. AGNPS Model Assessment for a Mixed Forested Watershed in Thailand. ScienceAsia 2006, 32, 53–61. [Google Scholar]
  40. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  41. Gupta, H.V.; Kling, H.; Yilmaz, K.K.; Martinez, G.F. Decomposition of the mean squared error and NSE performance criteria: Implications for improving hydrological modelling. J. Hydrol. 2009, 377, 80–91. [Google Scholar] [CrossRef]
  42. Neter, J.; Wasserman, W.; Kutner, M.H. Applied Statistical Models; Richard D. Irwin, Inc.: Burr Ridge, IL, USA, 1990. [Google Scholar]
  43. Reed, S.; Koren, V.; Smith, M.; Zhang, Z.; Moreda, F.; Seo, D.J.; Participants, D.M.I.P. Overall distributed model intercomparison project results. J. Hydrol. 2004, 298, 27–60. [Google Scholar] [CrossRef]
  44. Sabzevari, T.; Ardakanian, R.; Shamsaee, A.; Talebi, A. Estimation of flood hydrograph in no statistical watersheds using HEC-HMS model and GIS (Case study: Kasilian watershed). J. Water Eng. 2009, 4, 1–11. [Google Scholar]
  45. Cheng, C.-T.; Ou, C.; Chau, K. Combining a fuzzy optimal model with a genetic algorithm to solve multi-objective rainfall–runoff model calibration. J. Hydrol. 2002, 268, 72–86. [Google Scholar] [CrossRef]
  46. Gupta, H.V.; Sorooshian, S.; Yapo, P.O. Status of automatic calibration for hydrologic models: Comparison with multilevel expert calibration. J. Hydrol. Eng. 1999, 4, 135–143. [Google Scholar] [CrossRef]
  47. Zou, K.H.; Tuncali, K.; Silverman, S.G. Correlation and simple linear regression. Radiology 2003, 227, 617–628. [Google Scholar] [CrossRef] [PubMed]
  48. Moriasi, D.N.; Arnold, J.G.; Van Liew, M.W.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  49. Gharib, M.; Motamedvaziri, B.; Ghermezcheshmeh, B.; Ahmadi, H. Evaluation Of Modclark Model For Simulating Rainfall-Runoff In Tangrah Watershed, Iran. Appl. Ecol. Env. Res. 2018, 16, 1053–1068. [Google Scholar] [CrossRef]
  50. Zelelew, D.G. Spatial mapping and testing the applicability of the curve number method for ungauged catchments in Northern Ethiopia. J. Soil Water Conserv. 2017, 5, 293–301. [Google Scholar] [CrossRef]
  51. Hawkins, R.H. Asymptotic determination of runoff curve numbers from data. J. Irrig. Drain. Eng. 1993, 119, 334–345. [Google Scholar] [CrossRef]
  52. Jin, X.; Xu, C.Y.; Zhang, Q.; Chen, Y.D. Regionalization study of a conceptual hydrological model in Dongjiang basin, south China. Quat. Int. 2009, 208, 129–137. [Google Scholar] [CrossRef]
Figure 1. Location of (a) Lake Tana Basin with hydro-meteorological stations and (b) Basin Model of Gilgel Abay Catchment.
Figure 1. Location of (a) Lake Tana Basin with hydro-meteorological stations and (b) Basin Model of Gilgel Abay Catchment.
Hydrology 06 00021 g001
Figure 2. Comparison of simulated and observed hydrographs during calibration. (a) Model output for event 3 (14 Julyy–02 August, 2006), (b) rainfall of contributing stations and simulated and observed hydrographs of Event 1 (02–26 Julyy, 1996), and (c) rainfall of contributing stations and simulated and observed hydrographs of Event 6 (02–23 Julyy, 2008).
Figure 2. Comparison of simulated and observed hydrographs during calibration. (a) Model output for event 3 (14 Julyy–02 August, 2006), (b) rainfall of contributing stations and simulated and observed hydrographs of Event 1 (02–26 Julyy, 1996), and (c) rainfall of contributing stations and simulated and observed hydrographs of Event 6 (02–23 Julyy, 2008).
Hydrology 06 00021 g002aHydrology 06 00021 g002b
Figure 3. Correlation between observed and simulated flow during calibration (a) Event 1, (b) Event 4, and (c) Event 6.
Figure 3. Correlation between observed and simulated flow during calibration (a) Event 1, (b) Event 4, and (c) Event 6.
Hydrology 06 00021 g003
Figure 4. Comparison of simulated and observed hydrographs during validation (a) Model output for event 8 (26 August–11 September, 2008) and (b) Rainfall of contributing stations and simulated and observed hydrographs of Event 9 (04–20 Julyy, 2012).
Figure 4. Comparison of simulated and observed hydrographs during validation (a) Model output for event 8 (26 August–11 September, 2008) and (b) Rainfall of contributing stations and simulated and observed hydrographs of Event 9 (04–20 Julyy, 2012).
Hydrology 06 00021 g004
Figure 5. Correlation between observed and simulated flow values after validation (a) Event 7 (02–15 August, 2008) and (b) Event 9 (04–20 Julyy, 2012).
Figure 5. Correlation between observed and simulated flow values after validation (a) Event 7 (02–15 August, 2008) and (b) Event 9 (04–20 Julyy, 2012).
Hydrology 06 00021 g005
Table 1. Rainfall events selected for calibration and validation.
Table 1. Rainfall events selected for calibration and validation.
EventsStart DateStart TimeEnd DateEnd TimeRemark
Event 102 Julyy 199600:0026 July 199600:00Calibration
Event 216 Aug199600:0031 August 199600:00Calibration
Event 314 July 200600:0002 August 200600:00Calibration
Event 405 August 200600:0017 August 200600:00Calibration
Event 501 September 200600:0019 August 200600:00Calibration
Event 602 July 200800:0023 July 200600:00Calibration
Event 702 August 200800:0015 August 200800:00Validation
Event 826 August 200800:0011 September 200800:00Validation
Event 904 July 201200:0020 July 201200:00Validation
Event 1008 August 201200:0023 August 201200:00Validation
Table 2. Physiographic characteristics of the catchment.
Table 2. Physiographic characteristics of the catchment.
Sub-Catchment CodeArea(km2)Perimeter (km)Slope (%)CN *Lag (h)
W340232.59124.80017.7385.003.26
W370480.07189.24023.3084.674.44
W480356.49155.56018.684.473.88
W540540.22175.52024.9784.872.94
*-Dimensionless.
Table 3. Simulated (before and after optimization) and observed peak discharges and total volume and after optimization error functions.
Table 3. Simulated (before and after optimization) and observed peak discharges and total volume and after optimization error functions.
EventsPeak Discharge (m3/s)Total Volume (mm)NSER2
SimulatedObservedREPSimulatedObservedREV
BOPAOPBOPAOP
Event 1341.8310.6313.81.02206.56276.02275.83−0.070.8120.857
Event 2365.2286.7277.7−3.24112.01146.00154.035.210.7070.856
Event 3505.8358.6355.5−0.87242.31154.67160.853.840.7690.803
Event 4511.9335.2330.1−1.54154.31137.3142.643.740.7940.971
Event 5427.2227.1217.3−4.51138.93160.3158.49−1.140.7880.876
Event 6429.2248.1248.80.28197.08131.34134.792.560.6030.689
Mean430.2294.4290.51.91175.2167.6171.12.760.7460.842
BOP = Before Optimization, AOP = After Optimization, REP = Relative Error peak, REV = Relative Error volume.
Table 4. Calculated and optimized parameters of the catchment.
Table 4. Calculated and optimized parameters of the catchment.
S.No.MethodsParameterValues
CalculatedOptimized
1Loss rate parameterCN scale factor1.18920.8977
Ia scale factor0.98151.0527
Impervious area (%)0.00.0
2Runoff TransformLag time217.738217.738
3Routing Method constantsK0.50.5
X0.10.1
Table 5. Optimized parameter values of each event during calibration.
Table 5. Optimized parameter values of each event during calibration.
ParametersIa Scale FactorCN Scale FactorTlagKX
Events
Event 10.836721.1721217.7380.50.1
Event 20.970460.8813217.7380.50.1
Event 31.5150.91532217.7380.50.1
Event 41.0150.81121217.7380.50.1
Event 50.988921.0074217.7380.50.1
Event 60.98980.59913217.7380.50.1
Mean1.05270.8977217.7380.50.1
Table 6. Simulated and observed peak discharge and total volume during validation.
Table 6. Simulated and observed peak discharge and total volume during validation.
EventsPeak Discharge (m3/s)Total Volume (mm)NSER2
SimulatedObservedREPSimulatedObservedREV
Event 7295.4298.00.87115.99119.052.570.9220.958
Event 8188.9192.41.82111.77112.540.6780.8480.908
Event 9195.5198.51.5179.5681.472.340.8470.892
Event 10242.1246.41.75121.31126.263.920.9190.943
Mean230.5233..81.49107.16109.832.380.8840.925

Share and Cite

MDPI and ACS Style

Tassew, B.G.; Belete, M.A.; Miegel, K. Application of HEC-HMS Model for Flow Simulation in the Lake Tana Basin: The Case of Gilgel Abay Catchment, Upper Blue Nile Basin, Ethiopia. Hydrology 2019, 6, 21. https://0-doi-org.brum.beds.ac.uk/10.3390/hydrology6010021

AMA Style

Tassew BG, Belete MA, Miegel K. Application of HEC-HMS Model for Flow Simulation in the Lake Tana Basin: The Case of Gilgel Abay Catchment, Upper Blue Nile Basin, Ethiopia. Hydrology. 2019; 6(1):21. https://0-doi-org.brum.beds.ac.uk/10.3390/hydrology6010021

Chicago/Turabian Style

Tassew, Bitew G., Mulugeta A. Belete, and K. Miegel. 2019. "Application of HEC-HMS Model for Flow Simulation in the Lake Tana Basin: The Case of Gilgel Abay Catchment, Upper Blue Nile Basin, Ethiopia" Hydrology 6, no. 1: 21. https://0-doi-org.brum.beds.ac.uk/10.3390/hydrology6010021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop