^{1}

^{2}

^{2}

The Sefrou watershed, located in the northern part of the Tabular Middle Atlas, is the primary watershed in the Sebou region. Due to increasing socio-economic development, the need for a stable water supply is critical. Therefore, this watershed’s physiography and morphology are studied to aim and facilitate further research in the Sebou region. We rely on the hydrological regime, frequency analysis of extreme flows, rainfall volume, and flood periods forecast to enable that characterization. Our findings determined that the Sebou region is of pluvial type with abundant rainfall, especially in winters due to floods. Moreover, substratum parameters induce a large volume of water resources from upstream to downstream. However, the weak structure of the hydrographic network and the water balance ha s highlighted a considerable loss of these resources, which handicaps the irrigation system and the supply of drinking water to the region’s cities. Finally, the annual modulus and the different coefficients show hydroclimatic fluctuations of semi-arid climatic conditions, with a temperate winter. Our results suggest a need to monitor water runoff and protect surface soils from flooding and rapid erosion (by planting trees). We also propose installing barriers and a small dam to recharge groundwater and artificially store surface water in this watershed.

Morocco is subject to various climatic conditions that vary from north to south due to its latitude and west to east due to its heterogeneous topography (Beaudet, 1969; Marchane et al., 2021). As a result, rainfall is highly irregular, especially in the Atlas Mountains. The precipitation and the consequent flows of the rivers of the Sebou river basin in the north of Morocco have shown significant climatic fluctuations with an irregular but cyclic trend of about 8 to 12 years (Bouabid & Chafai Elalaoui, 2010; Singla et al., 2010). Adaptation to this phenomenon requires implementing mechanisms and actions to recognize, preserve, and safeguard water resources.

The Sebou basin has been studied in its northern part (Ouergha watershed; Msatef et al., 2018) and also in the southeast (Haut Sebou watershed; Qadem, 2015), and its western part (Sefrou watershed) remains to be characterized.

In this paper, we focus on studying the Sefrou region, an essential component of the southeastern part of the Great Sebou Basin (primary rainfall receiver in northern Morocco; Agoumi & Debbarh, 2006), more specifically, the region’s primary watershed. Our choice is due to the ever-increasing socio-economic development that the region experiences. Such development implies a significant increase in water needs and supply for agriculture, industry, and drinking water in the future. Therefore, data collection and processing are necessary to accommodate the increase in water demand in the region.

The study we present here concentrates on the characteristics of the overall water dynamics of this watershed. Thus, our study will enable data-driven decision-making that will, in turn, facilitate the sustainable development of the Sefrou region. To achieve that goal, we propose a physical characterization analysis of the Sefrou watershed in addition to topographic, morphological, hydrographic, and climatic characterization studies of the studied region using Geographic Information Systems (GIS). The processed data is fundamental to the evaluation of water resources. Furthermore, it will enable the analysis of the watershed’s behavior and monitor its hydrological response.

The climate data used in this paper are spreading over a reasonable period (more than 60 years, for temperature and precipitation). Also, the physical data come from fieldwork over four years (2016 to 2020) and their processing in the laboratory by different methods (ArcGis, DTM, …). We carried out the hydroclimatic data from five stations of the Sefrou watershed: Ain Timedrine, Sefrou, Azzaba, ABHS-Fez, and Allal Elfassi Dam, data made available by the Sebou River Basin Agency (ABHS). They allow us to examine the watershed variables (precipitation, temperature, aridity) and make graphical representations to understand better this watershed’s functioning (maps and curves). We can see the geographical coordinates of these stations are as follows (

Our approach to the study and characterize the Sefrou watershed consists in analyzing the data from four perspectives:

- physical characteristics and geometry,

Station | Station number | X | Y | Z |
---|---|---|---|---|

Ain Timedrine | 576 | 578,700 | 350,350 | 650 |

Sefrou | 6562 | 549,370 | 359,200 | 930 |

Azzaba | 1496 | 569,650 | 359,570 | 478 |

ABHS-Fez | 3817 | 535,400 | 384,800 | 415 |

Allal Elfassi Dam | 876 | 566,650 | 370,540 | 430 |

- Physiographic characteristics,

- Climatic context,

- Hydrological network.

The Sefrou watershed is located between the parallels (33.41˚ to 34˚N) and meridians (4.43˚ to 4.56˚W). This area crosses two vital rivers, which take Aggay River (south of Sefrou) and then El Yhoudi River (north of the Sefrou), as seen in ^{2} with a perimeter of 127.42 km. The main watercourse comes from mountainous sources: the latter funnels water through the tabular Middle Atlas via Sefrou. The water flow will continue towards the El Yhoudi River’s upstream, ending up in the Sebou River, crossing different geographical regions and bioclimatic environments (Benaabidate, 2000).

The geographical and geological aspects show that this watershed is part of the Atlas mountain that extend over 2000 km between the Moroccan Atlantic margin in the west and Tunisia in the east. Major lineaments bound this system called southern and northern Atlas faults, with adjacent sub-atlasic depressions subdivided within Morocco into the Middle and High Atlas.

The geological data show that this watershed is positioned on the NE part of the Tabular Middle Atlas and the Saiss plain (Termier, 1936; Termier & Dubar, 1940; Martin, 1981; Benshili, 1989; Cirac, 1987; Charrière, 1984, 1989, 1992; Fedan, 1988; Ahmamou, 1987; Fedan & El Hassani, 2018). The stratigraphic series of this area are subdivided into four groups (

A Paleozoic basement, Ordovician to the Carboniferous, is particularly characterized by a complete Devonian sequence of conglomerates, limestones, reef limestones, and carbonate marls (Aboussalam et al., 2020).

The post-Paleozoic cover consists of an Upper Triassic clay-evaporite series (gypsiferous clay and mixed basalts); followed by Kandar Dolomitic Formation of Sebkha’s environment (Charrière, 1992) of lower Liasic age dolomitic deposits. Then brecciated dolomites, dolomitic sands, and bedded dolomites; The Liasic series continues with massive limestones with Oncolites, Gastropods, Terebratula, and imperforate foraminifera, and finally, oolite limestones with reefs that show frequent dolomitization (El Hammichi et al., 2008, 2016).

After major unconformity, above the Jurassic series, the northern outcrop of Sefrou (in the Sais plain and its borders) contains marl formations that fill the basin (Ahmamou, 1987).

This stratigraphic sequence ends with the Sais’s lacustrine limestones (middle to upper Pliocene), and finally, the Quaternary Formations (travertines, crusts, alluviums, glacis, and terraces (Gourari, 2001).

Vegetation activity and soil type are intimately linked, and their combined actions have a significant influence on surface runoff. The vegetation cover retains a variable proportion of atmospheric water, depending on its density, nature, and the amount of precipitation. As a result, the intercepted water is partially removed from the runoff. The forest plays an essential role in the retention of this water. An index (K) characterizes this retention (in %) and corresponds to the forest area’s ratio to the watershed area.

We calculate the forest cover index (K) from the obtained vegetation map of the Sefrou watershed (

We describe the physical characteristics of the Sefrou watershed and its geometry (area, morphology, shape, topography, relief) through the general data. Knowledge of the geological context and configuration of the watershed are crucial elements to understand surface and underground flow mechanisms, considering that the quality of the substrate influences not only the groundwater flow but also the surface flow. Moreover, it helps to determine the bedrock’s permeability and explain its influence on the differentiation of the hydrogeological context.

The resulting data from this step will enable determining its morphology, interactions with meteorological phenomena, and hydrological behavior. Furthermore, Hydroclimatic indicators (i.e., precipitation, temperature, vegetation, runoff) are necessary to measure the influence of relief on the watershed function; we will focus on:

- Hypsometric curve: To better understand the variations in altitude within the Sefrou watershed, and to translate the distribution of altitudes within the study area. We calculate the average altitude (H_{avg}) using the following equation:

H a v g = ∑ i = 0 N A i × h i A (1)

where: A_{i}: Area between two contour lines (km^{2}) per observation; h_{i}: Average altitude between two contour lines (m) per observation; A: Total area of the watershed (km^{2}), and N: Total number of observation.

The median altitude is the value found at 50% of the total watershed area.

Odd numbers:

H M e d = | N + 1 2 | t h observation (2)

Even numbers:

H M e d = | N 2 | t h + | N 2 + 1 | t h 2 observation (3)

- Concentration time is defined as the time needed for water to flow from the most remote point in a watershed to the watershed outlet (Haan et al., 1994).

The time concentration-time of water is not only a function of the shape of the basin but also of other parameters:

· lithology (permeable formations favor infiltration at the expense of runoff, which delays the concentration of water at the outlet),

· rainfall (the higher the rainfall intensity, the shorter the concentration-time, regardless of the soil type),

· And relief.

There are several formulas for its calculation. To achieve the latter, we use the formula of Giandotti (1937):

T c = 4 A + 1.5 × L 0.8 H a v g (4)

where: T_{c}: concentration-time (hours); A: area of the watershed (km^{2}); L: length of the main thalweg watercourse (km); H_{avg}: average altitude (m).

- The global slope index Ig (in m/km) allows the assessment of the basin’s relief importance. We define it as the ratio between the useful drop (D_{u}) and the length (L) of the equivalent rectangle, using the following equation:

I g = D u L e q = H 5 % − H 95 % L e q (5)

where: H_{5%}: altitude which corresponds to 5% of the cumulated surface; H_{95%}: altitude which corresponds to 95% of the cumulated surface; D_{u}: Difference in altitude H_{5%} - H_{95%}; and L_{eq}: length of the equivalent rectangle (km).

We approach the specific physiographic characteristics of the Sefrou watershed thanks to the different morphological indices (Gravelius compactness index, Horton compactness index).

Gravilius’s compactness index (K_{c}): According to Roche (1963), the index used by hydrologists to characterize the watershed shape is Gravilius’s compactness index (Gravelius, 1914). We calculate Gravilius’s compactness index by dividing the perimeter of the watershed by that of a circle of the same area.

K c = P 2 π A = 0.28 P A (6)

where: K_{c}: Gravilius’s compactness index; A: watershed surface (km^{2}); and P: Perimeter of the watershed (km).

Horton’s compactness index (Horton, 1932) provides information on the geometry of the watershed. We calculate it by dividing the average width of the watershed by the length of the mainstream, using the following equation:

K h = W a v g L 2 (7)

where: K_{h}: Horton’s compactness index; W_{avg}: average width of the watershed (km); and L: length of the mainstream (km).

We focus on studying the climatic context, necessary for the understanding of the supply mechanisms and circulation of surface water (rainfall, potential evapotranspiration, actual evapotranspiration, and amount assessment of the water falling). Furthermore, we develop hydrographic network characteristics (drainage density, torrentiality coefficient, network hierarchy, longitudinal profile of the watershed). Evapotranspiration is an essential component in water balance and irrigation scheduling models (Fisher & Pringle III, 2013).

Potential evapotranspiration (PET, expressed in mm) is a concept defined by Thornthwaite and Mather (1955); it expresses the quantity of water that would evaporate from the soil and be transpired by the vegetation if there were no water supply problems. It is therefore independent of the actual water availability of the soil. It depends only on the capacity of the ambient atmospheric environment to induce evapotranspiration and thus mainly involves energy indicators.

Actual evapotranspiration (AET, expressed in mm) is the amount of water evaporating from the soil/substrate and plants when the soil is at its natural moisture content. Actual evapotranspiration explains water and energy exchange between the soil, the land surface, and the atmosphere. Thus, actual evapotranspiration is a significant component of the hydrological cycle and an essential physical process in natural ecosystems.

Therefore, the AET is an index derived from the calculation of the soil water balance, which represents the water lost in the form of vapor through evaporation from the soil and plant cover and transpiration by plants (Thornthwaite & Mather, 1955). Therefore, its calculation requires combining each month’s temperature and solar radiation data to calculate PTEs with Turc’s formula and precipitation and maximum soil water reserve (MWR) values to estimate the available water stock (Piedallu et al., 2013).

Three methods are used to characterize the AET of this watershed:

Turc method (Turc, 1961): Turc’s method allows us to evaluate the average annual runoff of rivers as a function of precipitation (P in mm) and temperature (t in ˚C) using the following equation:

A E T = P 0.9 + P 2 L 2 with L = 300 + 25 t + 0.05 t 3 (8)

where: Actual evapotranspiration AET (mm); P: Average annual precipitation height (mm); and t: Average annual temperature (˚C).

Coutagne method (Coutagne, 1954), who suggest that the effective deficit of flow (AET) should be symbolically expressed as a function of precipitation (P) by a parabolic formula:

A E T = P − γ P 2 with γ = 1 0.8 + 0.14 × T (9)

where: P: is the average annual precipitation (m); and T: is the annual average temperature (˚C).

Thornthwaite balance method (Thornthwaite & Mather, 1955): The calculation of AET by this method depends on the Easily Usable Reserve (EUR) of the Soil. The maximum value for a soil to be saturated varies according to the region’s climate; for humid areas, Thornthwaite admitted that soil is saturated when it absorbs a water flow equivalent to 100 mm of precipitation. On the other hand, this value can be reduced to 50 mm in the case of arid and semi-arid areas because the soil is not constantly saturated. Thornthwaite (1948) drew a correlation between monthly average temperature with monthly potential evapotranspiration. The Thornthwaite method is based on the following equation (in Naoura, 2012):

P E T = 16 × ( 10 × T I ) α (10)

where: PET: Potential evapotranspiration; T: Average temperature of the month considered, expressed in ˚C; I: The sum of the monthly thermal indices i: I = ∑ i = ∑ i = 0 ( T 5 ) 1.514

In this equation α is a coefficient expressed as a function of the sum of the monthly thermal indices:

When

10 ˚ C ≤ I ≤ 80 ˚ C : α = ( 1.6 × I 100 ) + 0.5 (11)

Else:

α = ( 675 × 10 − 9 × I 3 ) − ( 771 × 10 − 7 × I 2 ) − ( 1792 × 10 − 5 × I ) + 0.49239 (12)

Evaluation of the waterfall amount

The amount assessment of the water falling on the Sefrou watershed is estimated by several methods (The arithmetic mean, Theissen polygon method, and Isohyet method) allowing the integration of the point data relative to the different stations, and help to estimate the spatial precipitations.

Arithmetic mean: The annual arithmetic mean can be estimated by the following equation:

P a v g = ∑ i = 1 n P i n (13)

where P_{avg}: Average precipitation over the watershed; P_{i}: Precipitation height at station I; and n: Number of stations in the watershed.

Theissen polygon method

Theissen polygon is a statistical method that averages the precipitation measured at the rainfall stations in the watershed region (Thiessen, 1911). The area of influence of each measuring station called Thiessen polygon is the weighting factor. These polygons are obtained graphically by plotting the perpendicular bisectors of the segments connecting two neighboring rain gauge stations on a geographical or topographical map.

The Thiessen polygon method is widely used because it is easy to apply and generally gives good results. It is particularly suitable when the rainfall network is not spatially homogeneous (irregularly distributed rain gauges). We perform the calculation using the following equation:

P a v g = ∑ i = 1 n A i × P i A t (14)

where P_{avg}: average precipitation over the watershed (mm); A_{t}: total area of the watershed (km^{2}); P_{i}: precipitation recorded at station I; and A_{i}: area of the polygon associated with the station (km^{2}).

Isohyet method

“The isohyetal method is used to estimate the mean precipitation across an area by drawing lines of equal precipitation. The method uses topographic and other data to yield reliable estimates” (Jain & Singh, 2005). When drawing the isohyet curves, we calculate the average rainfall as follows:

P a v g = ∑ i = 1 k A i × P i A t with P i = h i + h i + 1 2 (15)

where P_{avg}: average precipitation over the basin (mm); A_{t}: total area of the watershed (km^{2}); A_{i}: surface between two isohyets i and i + 1 (km^{2}); k: total number of isohyets; and P_{i}: average of the heights h of precipitation between two isohyets i and i + 1 (mm).

We study the hydrological network (drainage density, torrentiality coefficient, hierarchy of the watershed, longitudinal profile) and its characterization. It consists of collecting data and analyzing monthly and annual rate flows and rain flow relationships. It allows determining the flow regime in the watershed and frequency analysis to choose the return period of floods and the adjustment of flows by statistical laws (Average of monthly and annual flows). In this study, we use the following methods and equations:

Drainage density, first defined by Horton (1945), is defined as the total length of channel in a drainage basin divided by the watershed total area, according to the following equation:

D d = ∑ L A (16)

where: L: total length of the river system (km); and A: surface of the watershed (km^{2})

The torrentiality coefficient (or drainage texture) considers the frequency of elementary thalwegs (of low order, generally of order (1) by the drainage density. The torrentiality coefficient (for order 1 streams) estimates the adaptability of a watershed to surface runoff or infiltration.

It can be expressed by the Smith equation (Smith, 1950):

T = D d × F s (17)

where: D_{d}: drainage density, (D_{d} = 0.736 km^{–1}); F_{s}: designates the frequency of elementary thalwegs (F_{s} = N_{1}/A; with N_{1}: number of streams of order 1; F_{s} = 0.21).

To calculate the average annual flows, only the data of Ain Elouali station are available. The location of this station is in the North of the Sefrou watershed (_{sp}) measures of the annual average pluvial flow in a watershed is defined as the number of liters of water flowing on average each second per square kilometer of the basin. We use specific flow to express peak flows during floods, and in general, this peak decreases with increasing watershed size. For its calculation, we use the value of the flow Q (l/s or m^{3}/s) for the surface of the watershed (A in km^{2}), and the equation corresponds to:

Q s p = Q A (18)

In this methodology part, we emphasize the geographical and geological aspects, the occupation of soils (land use) as well as the characterization of the studied climatic zone according to the empirical climagram of Emberger (Sauvage, 1963).

We use the rainfall quotient (Emberger climate index, discussed in Calvet, 1979) to define the five different types of Mediterranean climates, from arid (desert) to humid (high mountain). Morocco is the only country in the Mediterranean region that has all five types. The Emberger rainfall quotient (Q2) is a synthetic expression of the Mediterranean climate, considering the annual average of rainfall and temperature (Sauvage, 1960). The latter take into consideration, on the one hand, the average of the minima of the coldest month (m) and, on the other hand, the average of the maxima of the warmest month (M). We use these two extreme thermal values to evaluate the average temperature (M + m)/2 and the average mean thermal amplitude (M − m). We calculate this quotient according to the following equation:

Q 2 = 2000 × P ( M + m ) ( M − m ) (19)

where Q2: Emberger’s rainfall quotient; P: average annual precipitation (P = 454.22 mm); M: average of maximum of hottest month (M = 307.79˚ Kelvin; M = 34.64 ˚C); m: average of the minima of the coldest month (m = 276.71˚ Kelvin; m = 3.56 ˚C).

Q2 = 50.1

Emberger bases his method on the surface limits of the various plant species, with Q2 values on the ordinate and minima average of the coldest month (˚C) on the abscissa of his climatogram. The graphical representation of the Sefrou watershed is as follows (

The Sefrou watershed relief characteristics study allows us to determine its morphology, interactions with meteorological phenomena, and hydrological behavior. The relief directly influences all the hydroclimatic indicators (i.e. precipitation, temperature, vegetation, runoff).

Hypsometric curve

To better understand the variations in altitude within the Sefrou watershed, we determined a hypsometric curve that allowed us to translate the distribution of altitudes within the study area. This curve also allows us to determine the characteristic altitudes (

From this curve, we conclude that the altitude varies enormously despite the relatively small surface area of the studied watershed, thus corresponding to a steep watershed.

Using the Equation (1) and the observation data (_{avg} = 928.36 m.

The median altitude is the value found at 50% of the total watershed area (Equations (2) and (3)).

On the hypsometric curve, the median altitude is H_{Med} = 905 m.

Concentration time

We calculate the Concentration time according to Equation (4) from the following data: A: area of the watershed (km^{2}) = 405 km^{2}; L: length of the main thalweg watercourse = 45.68 km; H_{avg}: average altitude (m) = 928.36 m.

The Sefrou watershed Concentration-time is: T_{c} = 6 h 11 min. This result classifies the Sefrou watershed as having a relatively short time of concentration.

The global slope index

We determine H_{5%} and H_{95%} using the hypsometric curve, thus D_{u} = 1080 m. In addition to that, from our initial parameters we know that: L_{eq} = 56.43 km. Therefore, the global slope index, according to Equation (5), will be:

I g = 1080 56.43 = 19.13 m / km

I g ≈ 0.019 %

The overall index of slope given by the formula is 19.13% or 0.019%. According to the classification of the relief given by the ORSTOM (Msatef et al., 2018), the value of the global index of slope of the Sefrou watershed shows that the relief is quite important (moderate relief).

The digitization of the Sefrou watershed covers 404.89 km^{2}, with a perimeter of 127.42 km. It shows that the shape of the hydrograph at the exit of the watershed depends on the latter’s form. This characterization is due to the different morphological indices that enable understanding and comparison with other watersheds. These indices are:

Gravilius’s compactness index (K_{c})

We use the Equation (6) to calculate K_{c} in the Sefrou watershed.

K_{c} = 1.786; this value allows us to conclude that the watershed is elongated with a potential linear erosion. For the same rainfall, this shape favors low flood peak flows because of the delay in water routing to the outlet.

Horton’s compactness index

In Sefrou’s watershed we found the that (Equation (7)): W_{avg} = 13.56 km et L = 45.76 km. Therefore, Horton’s compactness index is:

K h = 13.56 45.76 2 = 0.0647 km − 1

As the value of K_{h} is very low, this confirms that the Sefrou watershed is elongated.

In summary of the different physiographic characteristics, the Sefrou watershed (

The hydrographic network is relatively denser on the right bank than on the left one of the main river, and the two banks remain almost symmetrical according to their surfaces.

The Rainfall conditions of the watershed when examining the monthly rainfall averages shows that December is the wettest month of all stations, with a maximum recorded of 54.92 mm. Conversely, July is the driest month, with an average rainfall of 3.86 mm (

Shape characteristics | Perimeter (km) Perimeter (km) Gravelius Compactness Index (K_{c}) Horton compactness index (K_{h}) Length of equivalent rectangle (km) Width of equivalent rectangle (km) | 127.42 404.89 1.786 0.0647 57.089 7.094 |
---|---|---|

Characteristics of the relief | Maximum altitude (m) Minimum altitude (m) Average altitude (m) Median altitude (m) global slope index (%) Concentration time (T_{c}) | 1766 241 928.36 905 0.019 6 h 11 min |

Hydrological characteristics | Length of the hydrographic network (km) Drainage density (km -^{1} ) | 298.207 0.736 |

Annual precipitation decreases towards the SE of this watershed and directly influences runoff flow as follows: Sefrou: 517.2 mm/y; Anoceur station: 504 mm/y; Ain Timedrine: 407.4 mm/y; Allal Elfassi Dam: 335.3 mm/y; Azaba: 361 mm/y; and Fez: 396.2 mm/y.

The average precipitation for the study area is 420 mm/year, calculated from the results obtained in the six stations between 1957 and 2010. The orographic effect explains these variations due to the presence of sequences of high peaks in this region.

Potential evapotranspiration

The monthly PET values calculated for the watershed show significant fluctuations during the year. A minimum of 21.765 mm in January and a maximum of 130.65 mm in July and August explains this fluctuation. On the one hand, the low values recorded during the cold period are due to the low temperatures and the relatively high air humidity during the winter months. On the other hand, the high values of PET found in the hot period reflect conditions opposite to the cold period, i.e., high temperatures, low humidity, and low rainfall.

Furthermore, a third factor influencing the PET value is the low rate of transpiration caused by plants (in the Sefrou watershed, there is a relatively average rate of vegetation density).

Actual evapotranspiration

In Sefrou’s watershed, the actual evapotranspiration (AET) is very important because it explains the water and energy exchanges between the soil, the land surface, and the atmosphere. We have calculated it for the Sefrou watershed using the average of the three methods (Thornthwaite, 1948; Thornthwaite & Mather, 1955; Coutagne, 1954; Turc, 1955).

Turc method (Turc, 1961)

The value of the AET, calculated by the Turc method for the entire watershed is: L967.63˚C; Therefore the actual evapotranspiration is (from Equation (8)):

A E T = 454.22 ( 0.9 + 454.22 2 967.63 2 ) = 429.13 mm / year

Coutagne method (Coutagne, 1954)

The value of the AET calculated by the Coutagne method is (from Equation (9)):

γ = 0.29.

P = 0.45422 m.

AET = 0.39548 m/year = 395.48 mm/year.

Thornthwaite balance method (Thornthwaite et al., 1955)

To establish a balance that is undoubtedly correct, we use Equation (10), and we consider the two distinct periods:

· We have more precipitation during the Humid period, so when:

P + E U R ≥ P E T → A E T = P E T , we have a water excess ( P − A E T ), which will fill the EUR; the rest will feed the surface runoff and/or deep infiltration.

· During the Dry period, we have:

P + E U R ≤ P E T → A E T = P + E U R , andtwo cases arise:

- First, EUR is sufficient, then: AET = PET; in that case, PET is satisfied, and the EUR = EUR-complement to PET.

- Second, EUR is insufficient or null, then: P ≤ A E T ≤ P E T , this leads to the appearance of a deficit.

Analysis of the

- The monthly AET is maximum in the fall and spring months based on favorable weather conditions,

- The deficit reaches its maximum during the summer months due to the high PET values.

- The excess represents its maximum values for the winter months, which favors the watercourses feeding and the infiltration of the underground flows.

- The annual average AET value deduced from the Thornthwaite balance is 345.995 mm/year.

By using Equations (10, 11 and 12), In the case the Sefrou watershed, we found that I = ∑ i = 1 12 i = 79.50 ˚ C .

Since I is the range [10, 80], we will compute the coefficient α used the first formula:

α = ( 1.6 × I 100 ) + 0.5 = 1.772

The calculation of the AET by the Thornthwaite method is the following (

Average AET calculation

The average AET calculation using the three methods is (

The average amount of water lost over a year through evapotranspiration is 390.20 mm.

Arithmetic mean: P_{1}, P_{2}, P_{3}, P_{4}, P_{5}, and P_{6} are the annual average rainfall of the stations of Ain Timedrine, Sefrou, Azaba, Fez ABHS, Anoceur, and Allal Elfassi Dam, respectively over the whole series (1957-2010).

From paragraph (3.2.1.) we know that:

P (mm) | PET (mm) | P-PET | EUR | AET (mm) | Excess | Deficit | |
---|---|---|---|---|---|---|---|

Sep | 18.42 | 104.1 | −85.7 | 0 | 18.42 | - | 85.68 |

Oct | 41.1 | 56.215 | −15.15 | 0 | 41.1 | - | 15.12 |

Nov | 52.1 | 43.891 | 8.21 | 9.21 | 43.891 | - | - |

Dec | 54.92 | 27.788 | 27.14 | 27.13 | 27.788 | - | - |

Jan | 51.56 | 21.765 | 29.8 | 50 | 21.765 | 6.93 | - |

Feb | 51.68 | 27.833 | 23.84 | 50 | 27.833 | 23.84 | - |

Mar | 49.56 | 39.304 | 10.26 | 50 | 39.304 | 10.26 | - |

Apr | 51.04 | 46.634 | 4.41 | 50 | 46.634 | 4.41 | - |

May | 34.5 | 63.225 | −29.72 | 20.28 | 34.5 | - | - |

June | 16.36 | 98.04 | −81.68 | 0 | 36.64 | - | 61.4 |

July | 3.86 | 130.65 | −126.79 | 0 | 3.86 | - | 126.79 |

Aug | 4.26 | 129.68 | −125.42 | 0 | 4.26 | - | 125.42 |

Total | 429.36 | 789.12 | −360.8 | 256.62 | 345.99 | 45.44 | 414.41 |

Turc | Coutagne | Thornthwaite | Average | |
---|---|---|---|---|

ETR (mm) | 429.13 | 395.48 | 345.995 | 390.20 |

· P_{Sefrou} = 517.2 mm/y,

· P_{Anoceur station} = 504 mm/y,

· P_{Ain Timedrine} = 407.4 mm/y,

· P_{Allal ElFassi Dam} = 335.3 mm/y,

· P_{Azzaba} = 361 mm/y,

· P_{Fez} = 396.2 mm/y.

Accordingly to Equation (13) where n = 6, the annual arithmetic mean can be estimated as follow:

P a v g : = ∑ i = 1 n P i 6 = 420.18 mm

Theissen polygon method

We use the statistical Theissen polygon method (Equation (14)) to quantify the averages of the precipitation measured at the rainfall stations in the watershed region (

Station | A_{i}: Area (km^{2}) | P_{i} | P_{i} × A_{i} |
---|---|---|---|

Sefrou | 257.54 | 517.2 | 133,199.68 |

Allal Elfassi Dam | 72.78 | 335.3 | 24,401.79 |

Azzaba | 0.38 | 361 | 137.68 |

Anoceur | 74.19 | 504 | 37,392.77 |

Total (A_{t}) | 404.89 |

The amount of water falling on the Sefrou watershed P_{avg} = 481.93 mm.

Isohyet method

The isohyet map shows a progressive precipitation gradient between the different stations, thus facilitating the interpretation of precipitation. The precipitation gradient increases from north to south towards the higher altitudes (

The fallen water depth calculated by this method is P_{avg} = 460.87 mm.

Estimated mean of the water slide calculated by the three methods

We use the average value between the three methods to quantify the volume of water falling on the Sefrou watershed (

The annual rainfall volume (in m^{3}), for the Sefrou watershed, is estimated as follows:

V = P a v g × S

where V: volume of the falling water in the watershed (m^{3}/year); P_{avg}: average rainfall of the three methods (i.e.: 454.22 mm/year); and S: total surface of the watershed (i.e.: 404.89 km^{2}).

V = 454.22 × 404.89 = 183.909 × 10 6 m 3 / year

In Summary, the result obtained for the AET is 390.20 mm/year for the entire Sefrou watershed, while the annual average rainfall represents a value of 454.22 mm. Out of total annual precipitation of 183.909 million m^{3}/year for this watershed, the generated loss by this process is about 85.90% of the inflow to the water system, with a volume close to 157.988 million m^{3}/year. In comparison, the difference represents the volume of water drained and the amount of underground infiltration, which reaches 14.09% (or 25.921 million m^{3}/year), according to the following calculation:

The water retained volume by the soil is:

V R e t = ( 454.22 − 390.20 ) × 404.89 = 25.921 × 10 6 m 3 / year : ( 14.09 % ) ,

Isohyet | P_{i} (mm) | A_{i} (km^{2}) | P_{i}.A_{i} |
---|---|---|---|

P < 360 | 347 | 0.47 | 163.09 |

360 < P < 380 | 370 | 16.23 | 6005.1 |

380 < P < 400 | 390 | 35.49 | 13,841.1 |

400 < P < 420 | 410 | 35.33 | 14,485.3 |

420 < P < 440 | 430 | 36.79 | 15,819.7 |

440 < P < 460 | 450 | 41.69 | 18,760.5 |

460 < P < 480 | 470 | 39.29 | 18,466.3 |

480 < P < 500 | 490 | 133.37 | 65,351.3 |

P > 500 | 509 | 66.21 | 33,700.89 |

Total | 404.86 | 186,593.28 | |

P_{avg} (mm) | 460.87 |

Methods | Arithmetic average | Thiessen polygons | Isohyets | Average |
---|---|---|---|---|

P_{average} (mm) | 420.18 | 481.93 | 460.87 | 454.22 |

And the evaporated water volume is:

V e v a p = 390.20 × 404.89 = 157.988 × 10 6 m 3 / year : ( 85.90 % ) .

The hydrographical network refers to hierarchical and structured channels ensuring a watershed or region’s permanent or temporary drainage surface. This hierarchy is characterized by the increasing number of constituents from the first upstream ramification (order 1 of Strahler’s classification, Strahler, 1952) to the primary collector. This ranking increases with the watershed size, the number of tributaries, and drainage density (order 2, 3, 4, 5, or more). The density of a river’s tributaries increases under the following conditions: wet climate, steep slopes, and less porous rocks or surface formations.

Using the ArcHydro function of ArcGis, we were able to calculate the length of the main river and thalwegs of the Sefrou watershed. Thus, from upstream and high altitudes to the outlet, the main river extends over 45.68 km. The direction of flow of this watershed is from south to north.

We calculate the Sefrou watershed drainage density from Equation (16):

D d = 0.736 km − 1

The obtained value of D_{d} implies that the hydrographic network of the Sefrou watershed is dense.

The Sefrou watershed Torrentiality coefficient is (Equation (17)): T = 0.154.

According to Smith table (Smith, 1950), T indicates that it is a coarse texture. The value of T is relatively low since the torrentiality coefficient depends directly on the time of concentration (T_{c} = 6 h 11 min); this value is related to several natural indicators, such as climate, precipitation, vegetation cover, soil type, infiltration capacity, and terrain.

Using the profile along the main river (_{c} value (Concentration-Time). We note several breaks in slope indicating erosion at the break section, while the main tributaries occupy the right bank with appreciable density.

In Summary, to characterize the hydrographic network of the Sefrou watershed, we calculate the density of drainage, and the obtained value (D_{d}) implies that it is dense. The low permeability of some geological components (red Triassic

clays and Miocene marls in the northern part of the watershed), the vegetation cover’s weakness, and the mountainous relief in its southern part favor the acceleration of runoff and floods. The watershed’s morphology and the slope index’s calculation have shown several slope breaks, with intense erosion. These slopes are frequent in the southern (highest) part of the watershed and directly influence water runoff. Added to this is the high intensity of the hydrographic network of the main river right bank of the watershed (Wadi Aggay), which increases the runoff. Knowing that the Concentration-Time (T_{c}) influences the maximum flow observed, these slopes accelerate the surface flow.

The value obtained for the Torrentiality Coefficient (T) attests that this watershed has a coarse texture. The latter has a relatively low value and directly depends on the Concentration-Time (T_{c}). It’s linked to several natural indicators, such as climate, precipitation, vegetation cover, soil type, infiltration capacity, and relief.

The elongated form of the watershed, relatively strong stream slopes and low permeability outcrops are all indicators that influence the flow of very muddy surface water and the decrease in the infiltration/runoff ratio. Thus, the physical environment and climate may provide a suitable environment for the sudden flooding of rivers caused by heavy rainfall. So surface water’s response and Concentration-Time are shortened. This condition suggests carefully monitoring the protection of arable soils against flooding and erosion and installing barriers and a dam to recharge the area’s groundwater and artificially store surface water. Ouzerbane et al. (2019) described identical features in the Essaouira watershed (SW Morocco).

The curve obtained from the AHBS data (

Q s p = 1.121 ( l / s / km 2 )

For the Sefrou watershed, the highest interannual mean flow is 1000 l/s recorded in 1969, while the lowest is 280 l/s recorded in 2002.

The maximum monthly average flow appears in March (390 l/s), while the minimum occurs in August (325 l/s). The following hydrograph shows the variation of the monthly average flows at Ain Elouali station (

The graph lecture shows that water concentration requires seven months, i.e., from September to March, to reach its maximum (peak in March with 390 l/s). The time of flooding starts from the end of March until August (five months to reach the minimum quantity of this flow). The flood and receding times are approximate since we are working with an average monthly flow.

When comparing the hydrograph of average monthly flows (Q) and monthly precipitation (H) for the period 1968 to 2005 (

In Summary, the Sefrou watershed water balance assessment (

Inputs | Outputs | |||
---|---|---|---|---|

Precipitation | Evapotranspiration | Runoff | Infiltration | |

Volume m^{3}/year (10^{6}) | 183.96 | 161.28 | 11.21 | 11.47 |

Percentage | 100% | 87.67% | 6.09% | 6.23% |

Fluctuations and hydrological cycles characterize the rainfall values collected at the stations from 1957 to 2011. Moreover, the annual and monthly precipitation means are 454.22 mm and 35.45 mm, respectively. Furthermore, the yearly average temperature is 16.62˚C. Using the different methods, the combination of temperature and precipitation showed that the Sefrou watershed is globally under a semi-arid climate. The average value of actual evapotranspiration (AET) is 389.22 mm. Therefore, the study of the hydrological regime implies that an unstable regime characterizes the Sefrou watershed.

The analysis of the obtained data shows that the Rainfall-Flow Relationship is one of the essential components of the hydrological regime of the Sefrou watershed, fed mainly by rainfall. The lithological nature and morphology contribute to understanding its behavior related to precipitation.

The Sefrou watershed’s highest average annual flow is 1000 l/s recorded in 1969, while the lowest is 280 l/s recorded in 2002. Concerning the monthly average flow, its maximum occurs in March (390 l/s), while the minimum occurs in August (325 l/s). In the same way, the flooding starts from April until August (five months to reach the minimum quantity of this flow). Therefore, the flood and receding times are approximate since we work with a monthly average flow.

According to our results, the Sefrou watershed provides natural water resources. It ensures drinking water supply to the inhabitants of the northern Middle Atlas area and, more particularly, to the Sefrou habitats. Despite the usefulness of these essential capacities, they were causing risks represented by winter floods causing human and material damage. Therefore the creation of development projects in the region can solve the problem.

Finally, it seems that these characteristics are not specific to the Sefrou watershed, but somewhat identical in other watersheds, such as the adjacent watershed of Ouergha (Msatef et al., 2018), the Upper Sebou watershed (Qadem, 2015), and that of Wadi El Abid watershed in the High Atlas (Layati et al., 2021).

This study contributes to the physiographic and hydrological knowledge of the Sefrou watershed. It completes the global vision of the great Sebou basin in northern Morocco. This watershed shows that it is a semi-arid bioclimatic zone with temperate winter, abundant rainfall in winter, and very high flows. However, the summer period knows a drought which is sometimes strongly pronounced.

Rainfall measurements collected over the past 60 years show fluctuations between the different stations and hydrological cycles allowed to compute the annual average and monthly rainfall, 454.22 mm and 35.45 mm, respectively. Thus, the Sefrou watershed provides natural water resources and ensures drinking water supply to the inhabitants of the Middle Atlas regions, especially those of Sefrou city. However, regardless of the importance of these resources, the hierarchy of the hydrographic network of this basin (flow network) is poorly developed and branched, which leads to frequent winter floods, constituting a potential risk that can generate human and material damage.

The physiographic characteristics show an elongated relief in the Sefrou watershed by analyzing the water flow indicators. According to the hypsometric curve, the measurements between 500 and 900 m occupy most of the basin’s surface (moderate relief). However, based on the overall slope index and the specific gradient, the relief is relatively substantial and decreases towards the basin’s outlet (northern part of the study area).

The differences between the eastern and western sides of the drainage system are due to lithology and vegetation cover variations. In addition, the time of concentration is relatively short.

Furthermore, the high volume of these rains and the shortened concentration-time of these surface waters generally generate significant erosion, especially in the marl facies. These results suggest that decreasing the infiltration/runoff ratio is needed. For that, it is necessary to monitor water runoff, protect surface soils from flooding and rapid erosion (by planting trees), and install dikes and small dams to recharge the groundwater in this watershed.

The authors would like to thank the two anonymous reviewers who helped us structure our paper. We also thank: Prof. Dr. Ahmed El Hassani, Resident Member of the Hassan II Academy of Sciences and Technology, Morocco, for his support and verification of this paper’s scientific background and writing; Prof. Mohamed Chikhaoui, Professor of Agronomy at IAV Rabat, for his critical reading and advice on the formulation of this paper. We thank the Hassan II Academy of Sciences and Technology for its support through the project “Global warming Effects on the Desert landscape of Eastern and Southern Morocco” in which authors participate.

The authors declare no conflicts of interest regarding the publication of this paper.

El Hassani, F., Hattafi, Y., & Lahrach, A. (2021). Water Dynamics of the Sefrou Watershed, Northern Tabular Middle Atlas, Morocco. Journal of Geoscience and Environment Protection, 9, 102-130. https://doi.org/10.4236/gep.2021.99007