This study aims to compare the performance of two hydrological models, conceptual (HEC-HMS) and spatial (ATHYS) on the Aguibat Ezziar watershed. The comparative analysis is based on the performances of simulation in terms of Nash-Sutcliffe and RSR. The study requires the collection of a series of data as inputs models namely rainfall data, water quantity, soil occupation, DTM and requires also a calibration in order to evaluate these models in validation phase. The simulation results were obtained from the validation phase aiming to replicate the operation of watershed Aguibat Ezziar, and present a suitable adjustment perspective of the observed hydrograph. These results show that the objective is achieved and a model distributed like ATHYS plays an effective role to improve the efficiency and present a high advantage in anticipation of runoff volume comparing with other models.
Recent years have shown that the water resource mobilization reached a limit, that why it is essential now that this vital resource is managed in an efficient and sustainable manner. Hydrological models are tools to reproduce the operation of water systems, computing a period of rate return data for the design of a hydraulic structure, determination of a flood of a project, real-time flow forecasting on streams [
The Aguibat Ezziar watershed covers an area of 3640 km² and an elongated shape. It is located at the northwest of the Bouregreg watershed north central Morocco (
Aguibat Ezziar watershed is situated between an elevation ranging from 1615 m to 80 m at downstream of Bouregreg watershed. Between 1977-2004, the minimum annual rainfall is 336 mm/year in Lala Chafia station, while the maximum is 434 mm/year in Tsalat station. Furthermore, the annual and monthly average is respectively 394 mm and 39 mm (
ATHYS (ATelier Hydrologique Spatialisé) is rainfall-runoff transformation model, developed by the Research Institute for Development (IRD) of Montpellier. The basic principles of ATHYS are a hydrological environment for distributed modeling, including a series of models, DEM processing, hydrological and rainfall data and geographical display, spatial data interpolation. The flow generated for a rain event (rainfall-runoff transformation) is calculated for each cell, the amount of rain is estimated by the production model which will contribute to runoff. The transfer model calculated the hydrograph produced by each watershed outlet, this calculation is made
from the result obtained under the production function. The HEC-HMS is a reliable model developed by the US Army Corps of Engineers that could be used for many hydrological simulations [
Both models are used to simulate the flow rates daily time in Aguibat Ezziar station from 1996 to 2003. The calibration and validation period is respectively 1996-1999 and 2000-2003. To estimate the parameters of production and transfer models, a manual procedure by the sensitivity analysis had been used. From input data (DTM, precipitation), the hydrological model ATHYS is used to define the daily runoff flow for each cells of Aguibat Ezziar watershed, and also to obtain the drainage network, slope classes…etc. From the effective rainfall, SCS and Top Model functions define the runoff depth. Moreover, Unit hydrograph transfer function calculates the flows for each sub basins outlet, identified as nodes calculation (Tsalat, Lala Chafia and Aguibat Ezziar). The DEM used is a set of elevation points (three coordinates, x, y and z) located every 90 m. The topography is reconstructed by the assembly of these points (
In this study, three different production classes have been identified from the topography (slopes), as a file reclassified between 459 classes of natural slopes. Low and middle runoff zone correspond to the plains with slopes from 0 to 200 m/km, while the high runoff zone corresponds to the high slopes mainly on reliefs.
Climatologically information is accurately known in each station. This information must be spatializing on the catchment area. The method of Thiessen polygons was used to define the influence’s area for each station (
The area of high runoff concerns only 7.31% of the whole watershed area. This area corresponds to the reliefs with bare soil (and/or) with high slopes.
Data preparation for HEC-HMS model was performed manually for the three stations. The SCS Curve Number function is used to model the losses, and is based on knowledge of the type and land use. Runoff modeling is made by the SCS unit hydrograph model that integrates the coefficient Lag Time which can be theoretically determined by the formula “0.6 × concentration times”. To evaluate ours models two statistic coefficients were adopted: Nash-Sutcliffe efficiency coefficient (NSE, Equation (1)) [
Station | Thiessen coefficient | Area (Km²) | Extent % |
---|---|---|---|
Tsalat | 0.658 | 2395 | 65.90 |
Lala Chafia | 0.296 | 1077 | 29.60 |
Aguibat Ezziar | 0.045 | 164 | 4.50 |
Production class | Area (km2) | Extent % |
---|---|---|
1―Low runoff | 1331 | 36.57 |
2―Middle runoff | 2043 | 56.13 |
3―Rapid runoff | 266 | 7.31 |
Total | 3640 | 100 |
Yiobs: Observed discharges; Yisim: Simulated discharges; Yimean: Mean of observed discharges.
In this study, the functions used by ATHYS model are SCS and unit hydrograph for the production and transfer function. The sensitivity analysis has allowed to show the best values of the parameter S (SCS) or STO (Top Model) affecting mainly the peak flow. The decline of the hydrograph explained by the emptying time is generally influenced by the DS parameter. Preferentially, we fix the influenced parameters on the general curve of the hydrograph, namely, the diffusion coefficient (K0), the parameters K1 and Alfa, and the transfer speed of V0, which figures on the transfer function (
Parameter | Production class 1 | Production class 2 | Production class 3 | |
---|---|---|---|---|
Production function | STO [mm] | 36 | 45 | 49 |
INF [mm∙h−1] | 10.1 | 9.5 | 9.5 | |
F [m−1] | 40 | 39.6 | 38.5 | |
Ds [j−1] | 0.1 | 0.09 | 0.09 | |
CN | A.EZ [ | L.CHA [ | TSA [ | |
Impervious [%] | A.EZZIAR [0.33] | L.CHAFIA [0.42] | TSALAT [0.57] | |
Transfer Function | V0 [m∙s−1] | 1.25 | ||
Alfa [ad] | 0.5 | |||
K0 [ad] | 0.1 | |||
K1 [mn] | 0.9 | |||
Lag time [min] | A.EZZIAR [ | L.CHAFIA [ | TSALAT [ |
The statistic evaluators showed a good correlation between the monthly observed and simulated river discharge with NSE ATHYS of 0.74, NSE HEC-HMS of 0.73, and RSR ATHYS of 0.39 and RSR HEC-HMS of 0.36 for the calibration period. The validation period revealed good values for NSE 0.73 and RSR 0.32, these model performance for both calibration and validation periods is evaluated as “very good performance rating” which is defined by the flowing ranges: 0 to 0.5 for RSR, 0.75 to 1 for NSE (
The average difference between the observed and simulated monthly total flow is 4% which confirms a good model’s calibration for the daily time step.
Comparing ATHYS and HEC-HMS models on a long event for a long term (24 h) shows that the behavior of the model HEC-HMS is correct with reproductive hydrograph peaks in calibration, in other hand, ATHYS performance at Aguibat Ezziar station is very good for the stability of the model in continuous mode.
The comparison of models performance (ATHYS and HEC-HMS) on Aguibat Ezziar watershed shows that a distributed model like ATHYS plays an effective role to improve the efficiency and present a high advantage in anticipation of runoff volume. While, The HEC-HMS model is very correct with the estimate peaks and so it can be a more effective model for flood forecasting. This study had demonstrated the utility of distributed model and also showed the ability of ATHYS model to be used to simulate the water quantity in semi-arid regions [
The authors thank the Ministry of Energy, Mines, Water and Environment (S.Water-Dr.Brirhet) for providing the data on the studied boreholes.
MouradKhattatı,MostaphaSerroukh,IsmailRafık,HakimMesmoudı,HassaneBrırhet,YassineBouslıhım,FatimaHara, (2016) Hydrological Modeling of Aguibat Ezziar Watershed (Morocco), Comparative Study of Two Different Hydrological Models. Journal of Geographic Information System,08,50-56. doi: 10.4236/jgis.2016.81005