Hydrological Modeling of Upper OumErRabia Basin (Morocco), Comparative Study of the Event-Based and Continuous-Process HEC-HMS Model Methods

Human population growth and land-use changes raise demand and competition for water resources. The Upper OumErRabia River Basin is experiencing high rangeland and matorral conversion to irrigated agricultural land expansion. Given Morocco’s per capita water availability, River-basin hydrologic modelling could potentially bring together agricultural, water resources and conservation objectives. However, not everywhere have hydrological models considered events and continuous assessment of climatic data. In this study, HEC-HMS modelling approach is used to explore the event-based and continuous-process simulation of land-use and land cover change (LULCC) impact on water balance. The use of HEC-GeoHMS facilitated the digital data processing for coupling with the model. The basin’s physical characteristics and the hydro-climatic data helped to generate a geospatial database for HEC-HMS model. We analyzed baseline and future scenario changes for the 1980-2016 period using the SCS Curve-Number and the Soil Moisture Accounting (SMA) loss methods. SMA was coupled with the Hargreaves evapotranspiration method. Model calibration focused on reproducing observed basin runoff hydrograph. To evaluate the model performance for both calibration and validation, the Coefficient of determination (R2), Nash-Sutcliffe efficiency (NSE), Root Mean Square Error (RSR) and Percent Bias (PBIAS) criteria were exploited. The average calibration NSE values were 0.740 and 0.585 for event-based (daily) and continuous-process (annual) respectively. The R2, RSR and PBIAS values were 0.624, 0.634 and +16.7 respectively. This How to cite this paper: Msaddek, M., Kimbowa, G. and El Garouani, A. (2020) Hydrological Modeling of Upper OumErRabia Basin (Morocco), Comparative Study of the Event-Based and Continuous-Process HEC-HMS Model Methods. Computational Water, Energy, and Environmental Engineering, 9, 159-184. https://doi.org/10.4236/cweee.2020.94011 Received: September 3, 2020 Accepted: October 18, 2020 Published: October 21, 2020 Copyright © 2020 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access


Introduction
Morocco like the most countries, experiences anthropogenic pressure and land-use dynamics particularly in and around mountainous basins [1] [2]. The interactions among land-uses, soil types, vegetation cover, climate variability and subsequent impact on hydrologic behaviors of watersheds have been studied by various researchers [3]- [9]. For the recent years, more people have settled around the OumErRabia (OER) due to favorable climate thus increasing the overall basin's population. This has consequently raised demand and competition for land, water and food resources [10] and thus leading to irrigated agriculture expansion. In the Upper OumErRabia River Basin (UOERRB), land-use and land cover change (LULCC) trends mainly involve the conversion of rangeland and matorral to irrigated agriculture and forestry. Land use and climate related studies done in Morocco for example [11]- [16] indicate a need for comprehensive mechanisms to conserve river waters and strategic planning to enable water availability for various competing uses. However, most recent land-use studies conducted in the OumErRabia (OER) basin are limited to water quality variation [17] and soil fertility impact [18].
The unsustainable land-use practice involving overgrazing, intense cultivation, over-exploitation of forests resources especially by the poor communities is pointedly shifting the hydrologic characteristics of the OUERRB. This has resulted in amplified runoff, erosion, sedimentation and degradation threatening the basin's water availability and reliability. The growing water demands, recurring and prolonged droughts associated with increasing climate variability and changing LULC are pointedly responsible for several water allocation challenges in the UOERRB [12]. Consequently, this has caused tension among water resources planning, developing and management authorities [19] towards allocation of the already scarce water resources among the increasing users [20] [21]. However, [22] noted limited or hardly any studies conducted to back up these insights.
The Morocco government and local authorities realize the disastrous impact of recurrent droughts on valued water resources. However, much of the government's effort has focused on sensitization of people about sustainable agri-cultural practices, water and related resources exploitation. The bearing of anthropogenic actions and climate variability studies is left less vibrant yet existing government policies are hardly implemented. According to an assessment study by [23], the distressing reduction in the basin's surface and ground water resources has been majorly linked to the regional climate variability in the previous decades. There is increased establishment of private wells within the watershed further adversely impacting the region's ground water resources. There is need for sustainable land-use and management plans and strategies especially at the moment when the OumErRabia Hydraulic Basin Agency locally known as l'Agence du Bassin Hydraulique de l'Oum-Er-Rbia (ABHOER) is extending irrigated area to fulfil the country's food demands and agricultural industries' raw materials needs. The ABHOER recognizes anticipated depletion of major available water resources owing to frequent droughts. However, the quantification of LULCC and impact-assessment of these processes on upstream river basin hydrology is still inadequate. The idea of hydrologic modelling [24] [25] [26] has been in use for quite time [24] [27] and is widely applied for instance in management of water resources [28], impact-assessment of LULCC and climate change on hydrology, water resource evaluation planning and allocation [16] [18] [29] [30] [31] [32] [33]. The HEC-HMS was selected for this study because it is a well-documented semi distributed hydrologic model flexible for both temporal and spatial scales and can be easily set-up and used with medium expertise [34]. Besides, the model usually requires less input calibration parameters.
In Morocco, HEC-HMS modeling approach has been used in recent studies such as [35]- [40]. Despite several national and regional studies evaluating hydrologic responses to LULCC and climate variability across the OER basin [9] [41], there is hardly any study that exclusively focuses on the UOERRB watershed accounting for local trends. Therefore, analyzing the impact of LULCC on hydrologic behavior of a watershed is vital for sustainable water resources management. Thus the objectives of the study were to: 1) explore and understand impact of land-use and land cover change (LULCC) on the water balance of the UOERRB using HEC-HMS-modelling approach; and 2) compare the performance of the event-based and continuous-process methods of the HEC-HMS model. We focus on the role of forests and agroforestry and the impacts of forest transition and forest degradation on the water availability in a case study area UOERRB.

Study Area Description
The OER River draining the middle Atlas and the north-western High Atlas [11] is a major source of hydroelectric power and irrigation for a significant population [42] and thus a key ecosystem to the Moroccan economy. The 550 km long main stream of the basin, covering 35,000 km 2 of the surface, is an agriculturally rich region containing 50 percent of Morocco's public irrigated agriculture [23] [43]. Nevertheless, the study focused on the 1049 km 2 UOERRB: a hydrological unit in the Middle Atlas with an equivalent length and perimeter of 119 km and 255 km respectively. The basin is located between latitude 33˚ and 33˚05' North and longitude 5˚01' and 5˚08' West, with basin outlet located at the Tarhat weather station. This basin contributes to the large watershed drained by the OER River whose source is at an altitude of 1800 m, 47 km of the Khénifra city ( Figure 1). Geologically, the basin major soil substrates include limestone and dolomites, with the eastern boundary of Khénifra constituting Triassic formations. The Khénifra essentially comprises sandstones chists attributed to the ordovician with the limestone and dolomites located on flushing levels on the slopes. In the valleys, are the triassic red clays, constituting places of intense ravines mostly in the Middle Atlas margin slopes [44] [45]. At the Atlas mountain side, the basin experiences a Mediterranean climate characterized with high temperatures in the summer and warm autumn. There is an increasing rainfall distribution from downstream zones to upstream Atlas mountainous areas with significant spatial variation within the basin. The UOERRB's average annual temperature and rainfall are 16.1˚C and 600 mm respectively. The basin's mean minimum and maximum temperatures are 7˚C and 22˚C respectively. The minimum and maximum precipitation is 295 and 1300 mm respectively, with October to May being the wettest months contributing about 85% of the total annual rainfall. The basin's vegetation cover is characterized by the presence of matorrals with a predominance of olive, orange orchards and forests. The forest cover is mostly concentrated at the basin downstream. The basin's main economic activity is agriculture in the plains and irrigated fields. The basin's topography varies from 2400 m in the upstream south to 864 m in the downstream north.

HEC-HMS Model Inputs
Daily mean rainfall and hourly discharge data for the two gauging stations of Tarhat and Tamchachte ( Figure 2) were collected from the ABHOER. The watershed has limited number of gauging and weather stations. Weather data was required for SMA continuous simulation of the watershed stream flow. However, due to lack of observed measurements, the weather data used for this study was estimated by interpolation on measured climate data from the 4 weather stations across the study area obtained from Global Weather Data for SWAT. The data was partitioned into sets for testing, calibration and for validation. Table 1 presents the model input data sources and processing.

Data Pre-Processing
The 2016 LULC map was derived from Lands at satellite imagery data (https://earthexplorer.usgs.gov/) with 2.5 m resolution. The image was classified into seven classes by using Arc-GIS 10.3 and photo-interpretation method [46]. The major LULC in UOERRB (Table 2 and Figure 3) is forest, matorral, arable    The HSG map comprising of numerous soil textures was generated by digitizing from the pedological map of Khenifra region ( Figure 4). From the two gauging stations of Tarhat and Tamchachte (Table 3)    Tarhat respectively.

Methods
The study is based on analysis of climate and stream flow data used for HEC-HMS model calibration and validation. Field visits were done to ascertain the current land use in the basin for the model setup. The major inputs for the HEC-HMS model included the DEM, soil data, land use and land management.
The model was first calibrated and validated before simulating hydrologic processes scenarios. Multiple parameters were individually and manually adjusted (maximum-minimum) and the model was executed with identified reasonable ranges of the most sensitive parameters. Figure 6 illustrates detailed re-  HEC-HMS offers the opportunity to represent all natural or artificial entities installed in a basin [50]. These elements influence the rainfall-flow transformation from sub-basins, outfalls to rivers and diversion water channels. The meteorological model's objective is to distribute the rainfall seized over whole study area.
The HEC-HMS has seven different precipitation methods [52] for describing meteorology. Storm-applies to the daily precipitation depth one of project downpour distributions defined by the SCS; and HMR 52 storm-uses specified storm area and area duration precipitation curves to compute total precipitation depth. The choice of best method depends on: modeler's objectives (determination of the project flow, hydrologic study of the watershed among others) and data availability that is the quantity (number of stations, measurement period, etc.); the quality (rain gauges, measure, etc.) of precipitation data; and the type of modeling envisaged. However, almost all methods have the specificity of homogeneous spatial and temporal rainfall depth distribution [53]. For this study, the Specified Hyetograph was adopted.

Application of the HEC-HMS
We analyzed LULC change and climate variability relationship with hydrologic processes in the UOERRB. The HEC-HMS model was selected for the study as it is: widely used to assess hydrology and water availability in agricultural catchments around the world [54] [55] [56] [57]; open source and user-friendly in relations to handling input data [58]. HEC-HMS is a semi-distributed conceptual hydrological model potentially capable of spatio-temporal simulations of rainfall-runoff relations within a basin [59] [60]. The model can be used for computation of other downstream processes for instance channel routing and reservoir routing. Globally, HEC-HMS has fruitfully been applied in catchment modeling of several river basins [51] [61] [62]. In Morocco, there are several studies on the HEC-HMS application for instance [40] [63] [64] [65] [66]. Most of the studies assessed the impact of climate and land-use on hydrology at watershed level. Reference [22] established water availability of the OER basin using different models. However, failure by most findings to account for LULCC and climatic variability effects sparked this research.

LULC Change Scenarios
We considered LULC maps of 2002 and 2016 for impact analysis. Scenario 1 is about a negative effect of deforestation and urbanization on peak flows. The CN values were estimated by increasing the surface of the urbanized/built up area, and assuming low forest density vegetation cover. For Scenario 2, we assumed that the pressure exerted on the forest cover with parallel actions to reforest the bare soils would result in clear dense forests. The percent conversion of arable land and rangeland to forest cover was assumed to be 80 and 20 respectively.

Event-Based and Continuous-Process Calibration and Validation of the HEC-HMS Model
The calibration was done for the event-based and the continuous-process me-Computational Water, Energy, and Environmental Engineering thods and some of the results have been selectively presented. The event-based model calibration results for the 10 events between 1991-1996 and 2001-2006 are presented in Table 4. Out of the 10 events, the resultant hydrographs of six sampled events are presented in Figure 7 and Figure 8.  During model calibration, ranges of parameters were adopted. Table 5 summarizes the range of event-based parameters used during event-based model calibration.
Event-based validation was based on precipitation and flow data of 2006-2011 and a sample of the resultant hydrographs is presented in Figure 9. Table 6 summarizes the volume and peak flow for the 3 selected validation events. The 1991-2001 continuous-process calibration hydrograph is shown in Figure 10.
The model validation results with the 2001-2011 precipitation and stream flow data are presented in Table 7. A comparison of observed and simulated output from the SMA continuous validation modelling of daily flow series for the UOERRB is presented in Figure 10.

Model Performance Evaluation
The calibration and validation performance results of the HEC-HMS at the Tarhat gauging station, the UOERRB outlet indicate good predictability of both the event-based and continuous-process model. Table 8 summarizes the performance of the 10 calibration events considered for the study. The average NSE, R 2 , RSR and PBIAS were 0.697, 0.474, 0.828, and −8.128 respectively.
The performance rating of the three validation events is depicted in Table 9 with average NSE, R 2 , RSR and PBIAS values of 0.581, 0.692, 0.596, and 3.213 respectively.
For the continuous-process, the NSE, R 2 , RSR and PBIAS for bothcalibration and validation models are presented in Table 10.

Simulation of Future Land-Use and Land Cover Scenarios
A simple future LULC distribution pattern in Scenario 1, used as a direct input to the calibrated event No. 6 generated various CNs that resulted in numerous volumes and peak flows (Table 11). Similarly, Scenario 2 simulation (increase in forest cover) resulted in volume and peak flow values presented in Table 12.

Typical Peak Streamflow Trends and Duration
Vegetation plays an attenuating and important role during flood periods. An increase in vegetation results in a corresponding delay in the surface runoff thus the flood point is attenuated. From Table 13 and Figure 11, a decrease in vegetation cover reflected in CN increase directly resulted in peak flow increase.  Climate variability was assessed as the difference between the flow rates of

Discussion
The basin characteristics such as the elongated shape delays the flow of water to the outlet during a flood period and thus justify the results of the HEC-HMS model. We observed strong annual streamflow variation and the average inter-annual monthly flow was 21 m 3 /s. Variation in basin precipitation was found influential to streamflow trends thus the positive correlation between precipitation and peak flows. The changes in precipitation are directly reflected in surface runoff pattern (magnitude and frequency) of peak flows for the study period. From the LULC maps, forest cover accounts for 32.6 percent of total watershed area in 2002 compared to the 39.9 percent in 2016. The expansion rate of forest land attributed to reforestation directly affects the magnitudes of runoff volume and peak flows as presented in scenarios. The information on LULC change serves as a vital tool in decisions making and policy formulation by the local authority (ABHOER). From the 2016 LULC map, we categorize the UOERRB as an agricultural basin and propose irrigated land expansion to be sustainably balanced with re/afforestation to offset the likely tremendous impact on the basin's available water resources.

Land-Use and Land Cover Change Impact Analysis
Considering the 10 calibration events, the performance evaluation ratings were generally good. However, events 3 and 4 produced relatively poor results. This is reflected in better average validation performance rating values compared to the calibrated ones. Out of the 10 calibration events, only one had an NSE value less than 0.5 indicating a very good performance accordingly. The calibration results indicate an adequately done process where average and maximum NSE values are 0.697 and 0.996 respectively. In general, most of the RSR and PBIAS values are considered good according to performance ratings. Despite achieving relatively less acceptable RSR values, the overall performance is considered satisfactory. Generally considering average performance values, the calibration model was accurate.
The ratings indicate a strong correlation between simulated and observed runoff volumes and peak discharges for all chosen events. From the event-based validated model (2006-2011), the simulated daily hydrograph matched with observed stream flow with an underestimation of peak flows. The event-based validation performance results obtained were satisfactory and acceptable to simulate the basin runoff for future LULC pattern projection.

Continuous-Process Simulation Analysis and Scenario
The 1991-2001 continuous process calibration model comparisons between the observed and simulated flow rates hydrographs at the Tarhat station are presented in Figure 10. The calibration results indicate a good generated hydrograph thus an adequately done process. The performance rating values indicate a strong relationship between observed and simulated values achieved during calibration. However, there was a multitude of spikes during the calibration period. The simulated model timing of the peak matches well, but the peak flow is underestimated by 18.5%. Besides, there is under prediction of low flows for the simulated model. Considering the continuous-process simulation, the model fairly reproduced the observed hydrographs for 1991-2001. However, there are some scatters in peak flows. One can conclude that the calibrated HEC-HMS model simulates the runoff and it can be used to evaluate the effects of LULC change on runoff generation in the UOERRB.
Similarly, comparison of the simulated and observed runoff hydrographs with the continuous-process method based on the same performance criteria, indicate good model results. Relative to calibration, the event-based model validation performance during validation is generally reduced. However, the results remain adequately acceptable with average NSE, R 2 , RSR and PBIAS values of 0.581, 0.692, 0.596, and 3.213 respectively. For 2001-2011 continuous-process validation, the simulated flow values fairly correlate with the observed values. However, out of the 4-performance rating, the model satisfactorily fulfilled only the RSR and PBIAS. Besides, the NSE and R 2 of 0.341 and 0.361 respectively are close to 0.400 considered satisfactory. Therefore, the prediction of the validation model can be considered acceptable.

Sensitivity Analysis of SCS-CN and SMA Parameters
From the model results, imperviousness; simple canopy-max storage; SMA's soil storage, tension storage, and maximum infiltration were the most sensitive parameters for estimating runoff by continuous-process simulation. During the non-winter months, the basin follows a wetting and drying sequence. However, this study assumed only dry periods during SMA continuous-process modelling. The calibrated GW-1 storage depth varied between 1.45 and 4.5 mm, while the GW-2 storage depth values ranged between 22.5 to 36 mm. The GW-1 storage coefficient varied between 18 and 85 hours. The GW-2 storage coefficient values ranged from 36 to 1485 hours. This variability clearly indicates a non-uniformity behavior between the interflow and groundwater flow throughout the water years [67]. Overall, the model is reliable to reproduce stream flows in the basin and the abnormalities in the results can be explained by the impact of precipitation uncertainties.
Considering the UOERRB, the measured precipitation of Tarhat and Tamchachte stations are considered non-representative due to high spatial variability in the basin rainfall. The HEC-HMS results obtained are similar to those of [40] [64] [65] [68] [68]. The mean NSE values for calibration and validation were 0.65 and 0.62 [64] and 0.74 and 0.73 [40]. The average NSE values were 0.99 for studies by [65] [69]. However, all these studies were carried out from other basins rather than UOERRB and most of the research studies were event-based.
In SCS-CN event-based model, initial soil moisture is normally classified into dry, moderate and wet conditions that may lead to increased simulation error. However, the initial soil moisture content was not measured in field. The study was based on AMC-II conditions and [70] notes difficult in correct selection of CN values from available handbook tables. Development of SMA continuous-process model requires several parameters for quality model outputs. However, most of these parameters were not readily available and were thus deduced from published research. The Lands at images used in this study were taken during different weather seasons and thus separate analysis of these images could yield precise representative LULC and better SMA validation model results.

Conclusion
In this study, the HEC-HMS modelling approach is used to explore the event-based (SCS-CN) and the continuous-process (SMA) simulation of land-use and land cover change (LULCC) impact on water balance. The HEC-HMS performance results were good with average calibration NSE values of 0.740 and 0.585 for event-based and continuous-process respectively. We argue that the basin's hydrologic response to positive and negative LULCC scenarios is significant. The HEC-HMS model simulation of the daily stream flow at basin outlet was considered satisfactorily good and acceptable. However, there is a slight under and over prediction of the high flows. Like with hydrologic research, the results of the study should be carefully interpreted with some measure of uncertainty. The results of the model indicated that an increase in forest cover and irrigated agricultural land expansion negatively impacted the peak flow and river stream discharge volumes. This hydrological modeling could also be used to facilitate subsequent hydrological studies within the region. The outcomes of this study can be vital in the future flood prediction and warning. For further development of this study, we recommend: detailed LULC survey within the basin to verify the HEC-HMS; adopting field-based data approaches to establish the SMA parameters; localized PET measurements to give more representative results; and evaluating significance of snow effect on stream flows.