^{1}

^{*}

^{1}

^{1}

^{1}

^{2}

^{3}

The objective of this study is to model the hydrology in the Sidi Jabeur basin, located in Bouregreg watershed at the north-central of Morocco, using the spatially distributed model (ATHYS) in order to understand and determine the different watershed hydrological processes. The study requires the collection of a series of data as inputs models namely rainfall data, water quantity, soil occupation, digital terrain model and requires also a calibration in order to evaluate the model in validation phase. The simulation results are obtained from the validation phase aim to replicate the operation of Sidi Jabeur watershed, 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 in improving the efficiency and presents a high advantage in anticipation of runoff volume.

The water is the most important natural resource especially in the arid or semi-arid zones that face high population growth, scarcity of freshwater, irregularity of rainfall, excessive land use change and increasing vulnerability to risks such of drought, desertification and pollution. Thus, the availability and the sustainable use of this resource become the core of the local and national strategies and politics in these regions [

Managing water resources are mostly required at watershed scale [

The objective of this study is to model the dynamics of Sidi Jabeur watershed functioning with the reconstitution of river flow at outlet (Station Sidi Jabeur). The simulation of water intake in the hydrological station has been done with the use of a spatially distributed model ATHYS. However, the objective is to test the adaptability and ensures a good performance of the model. In this work, the time step of modeling is daily, this choice due to compare modeled flows with occasional discharge measurements on the hydrometric Sidi Jabeur station.

Sidi Jabeur basin is located at the south of the Bouregreg watershed north-central of Morocco (^{2} with an elevation ranging from 1602 m to 330 m (Sidi Jabeur outlet).

Between 1977-2005, the minimum annual rainfall is 288 mm/year in Ouljet Haboub station, while the maximum is 300 mm/year in Sidi Jabeur station. Furthermore, the annual and monthly average is respectively 294 mm and 28 mm (

ATHYS (A Telier 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 calculate the hydrograph produced by each watershed outlet, this calculation is made from the result obtained under the production function.

The spatially distributed model ATHYS are used to simulate the flow rates daily time in Sidi Jabeur station from 1996 to 2004. The calibration and validation period is respectively 1996-1999 and 2000-2004. To estimate the parameters of production and transfer models, a manual procedure by the sensitivity analysis had been used. From input data (digital terrain model, precipitation), the hydrological model ATHYS is used to define the daily runoff flow for each cells of Sidi Jabeur watershed, and also to obtain the drainage network, slope classes… etc. From the effective rainfall, SCS function defines the runoff depth. Moreover, unit hydrograph transfer function calculates the flows for each sub basins outlet, identified as nodes calculation (Ouljet Haboub and Sidi Jabeur). 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 447 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 spatially on the catchment area. The method of Thiessen polygons was used to define the influence’s area for each station (Ta- ble 1).

The following table shows the breakdown of the areas of Sidi Jabeur watershed based in the production classes (

The area of high runoff concerns only 5.15% of the whole watershed area. This area corresponds to the reliefs with bare soil (and/or) with high slopes.

To evaluate our model three statistic coefficients were adopted: Coefficients of Determination (R^{2}), Nash- Sutcliffe efficiency coefficient (NSE) [

Station | Thiessen coefficient | Area (Km^{2}) | Extent % |
---|---|---|---|

Sidi Jabeur | 0.367 | 1067 | 36.70 |

Ouljet Haboub | 0.633 | 1843 | 63.30 |

Production class | Area (km^{2}) | Extent % |
---|---|---|

1) Low runoff | 990 | 34.02 |

2) Middle runoff | 1770 | 60.82 |

3) Rapid runoff | 150 | 5.15 |

Total | 2910 | 100 |

In this study, the functions used by ATHYS model are SCS and unit hydrograph for the production and transfer function. The sensitivity analysis of model has allowed to show the best values of the parameter S (SCS) 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 (

The statistic evaluators showed a good correlation between the monthly observed and simulated river discharge with R^{2} of 0.92, NSE of 0.83, and RSR of 0.41 for the calibration period. The validation period revealed good values for R^{2} 0.90, NSE 0.81 and RSR 0.39, this 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 approach we adopted was to model the dynamics of Sidi Jabeur watershed functioning with the reconstitution of river flow at outlet (Station Sidi Jabeur) using the spatially distributed model (ATHYS). The model was successfully calibrated in the Sidi Jabeur watershed and produced good simulation results. The evaluation of the

Parameter | Prod class 1 | Prod class 2 | Prod class 3 | |
---|---|---|---|---|

Production function | STO [mm] | 39 | 43 | 51 |

INF [mm∙h^{−1}] | 10.3 | 9.3 | 9.3 | |

F [m^{−1}] | 41 | 39.4 | 38.4 | |

Ds [j^{−1}] | 0.1 | 0.08 | 0.08 | |

Transfer Function | V0 [m∙s^{−1}] | 1.35 | ||

Alfa [ad] | 0.6 | |||

K0 [ad] | 0.2 | |||

K1 [mn] | 0.9 |

model performance was carried out successfully with the recommended statistical coefficients, the comparison of observed and simulated flow stream revealed a Nash-Sutcliffe coefficient and R² superior to 0.81 for both calibration and validation periods. This study 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.

MouradKhattati,MostaphaSerroukh,IsmailRafik,HakimMesmoudi,BrirhetHassane,YassineBouslihim, (2016) Hydrological Modelling of Sidi Jabeur Watershed (Morocco) Using Spatially Distributed Model ATHYS. Journal of Geoscience and Environment Protection,04,77-83. doi: 10.4236/gep.2016.41009