Bulletin de l’Institut Scientifique, Rabat, Section Sciences de la Terre, 2017, n° 39, 87-101.

e-ISSN : 2458-7184

Quantification de l’érosion hydrique et de la sédimentation à partir de

modèles empiriques dans le bassin versant de Tahaddart
(Rif nord occidental, Maroc)
Quantification of water erosion and sedimentation using empirical models
in the Tahaddart watershed (Northwestern Rif, Morocco)

Mounia TAHIRI1, Hassan TABYAOUI2, Fatima EL HAMMICHI2, Mohammed ACHAB3,

Abdelfatah TAHIRI3, Hassan EL HADI1

1. Applied Geology, Geomatic & Environnement Laboratory, Ben Msik Sciences Faculty, Hassan II University of
Casablanca, Morocco. ([email protected])
2. Natural Ressources and Environnement Laboratory, Polydisciplinary Faculty, Taza, B.P. 1223, Taza, Morocco.
3. Geo-Biodiversity and Natural Patrimony Laboratory (GeoBio), Scientific Institute, Geophysics, Natural Patrimony
Research Center (GEOPAC) Mohammed V University in Rabat, Morocco.
___________________________________________________________________________

Résumé. L'érosion des sols est un grave problème découlant de la dégradation des sols, de l'intensification de
l'agriculture et d'autres activités anthropiques. L'évaluation de l'érosion des sols est utile dans les travaux de
planification et de conservation dans un bassin versant. La modélisation fourni une approche quantitative et
cohérente pour estimer l'érosion des sols et le rendement des sédiments dans un large éventail de conditions.
Dans la présente étude, les modèles de quantification de l'érosion des sols dérivés de l'Équation universelle de
perte de sol (USLE) intégrée au SIG ont été utilisées pour estimer les pertes de sol et le module de sédimentation
dans le bassin versant de Tahaddart situé dans la partie Nord-ouest du Rif (Maroc). Ce bassin correspond à une
zone humide à influence atlantique et méditerranéenne ayant une aire d’environ 1141 km2. Il comprend quatre
sous-bassins allongés E-W. Les paramètres des modèles révisés (RUSLE) et modifiés (MUSLE) ont été estimés
en utilisant des données de télédétection et les zones de probabilité d'érosion ont été déterminées à l'aide de SIG.
L'érosivité estimée des précipitations (R), l’érodibilité (K) des sols, les facteurs topographiques (LS) et les
facteurs de gestion des cultures (C) varient de 39,21 Mjmm/ha/an, 0,22 à 0,31, 0,5 à 5 et 0,22 à 0,33
respectivement. Les résultats indiquent que pour le bassin versant de Tahaddart, la perte potentielle est de
202,337 t/ha/an, la perte réelle estimée du sol est d'environ 54528 t/an, la perte moyenne pour une crue est de
697894 t/ha, la dégradation spécifique est 34743,45 t/an. Les pertes en terre générées sur les versants l’emportent
sur les pertes dues au réseau hydrographique et la sédimentation domine dans les cours d’eau et participent à
l’envasement des barrages. Ces valeurs sont reparties de façon inégale d’un sous-bassin à l’autre. Les résultats
peuvent certainement aider à mettre en œuvre des pratiques de gestion et de conservation des sols pour réduire
l'érosion du sol par sous-bassin et dans le bassin versant de Tahaddart.
Mots-clés: Erosion hydrique, sédimentation, modèles empiriques, SIG, bassin versant Tahaddart, Maroc.

Abstract. Soil erosion is a serious problem arising from soil degradation, intensification of agriculture and other
human activities. Assessment of soil erosion is useful in planning and conservation studies in a watershed.
Modeling can provide a quantitative and coherent approach to estimate soil erosion and sediment yields under a
wide range of conditions. In this study, soil loss models derived from the Universal Soil Loss Equation (USLE)
integrated with GIS has been used to estimate soil losses and sedimentation models in the Tahaddart watershed
in the northwestern Rif, Morocco. It is a wetland with atlantic and mediterranean influence, with an area of about
1141 km2. It comprises four elongated E-W sub-basins. The parameters of the revised (RUSLE) and modified
(MUSLE) models were estimated using remote sensing data and erosion probability areas were determined using
GIS. The estimated rainfall erosivity (R), soil erodibility (K), topographic factors (LS) and crop management
factors (C) range from 39.21 Mjmm/ha/yr, 0.22 to 0.31, 0.5 to 5 and 0.22 to 0.33 respectively. The results
indicate that the potential loss for the Tahaddart watershed is 202.337 t/ha/yr, the estimated actual soil loss is
about 54528 t/yr, the average loss for a flood is 697894 t/ha, the specific degradation is 34743.45 t/year. Soil
losses on the slopes outweigh the losses due to the hydrographic network, and sedimentation predominates in the
rivers and contributes to the siltation of the dams. These values are distributed unequally from one sub-basin to
another. The results can certainly help in implementation of soil management and conservation practices to
reduce the soil erosion in each sub-basin and in the entire Tahaddart watershed.
Keywords: Water erosion, sedimentation, empirical models, GIS, Tahaddart watershed, Morocco.
 88 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental

Abridged English version (T) (Tab. 6 and 7), soil occupation (Tab. 8, Fig. 4),
erodibility K (Tab. 9, Fig. 5), physical characteristics
Introduction (LS) (Fig. 6) and conservation practices (P) at each
spatial unit of the watershed. The GIS allows to
The purpose of this study is to evaluate the annual evaluate the different equations and therefore to map
soil loss using the revised (RUSLE) and Modified the rate of erosion at the watershed scale. We used
(RUSLE) universal equations of soil loss, and the 30m x 30m surface unit (pixel size) as the
Sedimentation Models in the Tahaddart Watershed estimation base unit of each parameter.
(NW Rif, Morocco) as a river catchment area with an
outlet to the Atlantic. It synthesizes and completes Potential erosion
our detailed works in two of its four sub-basins: Oued Table 10 summarizes the different calculations of
Sania sub-basin (Tahiri et al. 2014c) and Oued erosion. It follows that the sub-basin El Hachef
Haricha sub-basin (Tahiri et al. 2015). The procedure shows the highest values of the potential erosion
followed is summarized in figure 1. (271.26 t/ha/yr). This area is characterized by slopes
Geology and geomorphology and steep slopes lengths that led to high values of
topographic factor LS whose contribution in the
The Tahaddart watershed, located between generation of losses by the model used is probably
Tangier and Assilah is geologically connected to the the heaviest. The potential loss to the Tahaddart
western Rif consisting of the overlap of several thrust watershed is 202.337 t/ha/yr.
sheets based on the Tangier unit (Suter 1980; Fig. 2).
On this geological complex structure, composed by Erosion on the slopes: RUSLE model
the juxtaposition of terrains with different The crossing of the RUSLE model factors maps
competence and permeability and affected by an (Roose & Georges 2004) allowed to quantify the
important brittle tectonics, a branched hydrographic actual erosion or average annual land losses in
network is established (Fig. 2). The main rivers (oued t/ha/yr. The results (Tab. 11) show that the Mharhar
El Hachef and oued Mharhar) extend from east to sub-basin has the highest value of A (RUSLE) at
west over 650 and 480 km respectively and average 64.68 t/ha/yr. The actual annual loss for this sub-
slopes of 1.20 to 4.10%. Such a situation may basin is 20632.92 t/yr. The Kharroub sub-basin has
encourage the development of stream flows and an A (RUSLE) value of 52.17 t/ha/yr, but given its
sediment deposits at the level of the river beds. Four larger surface area, the annual real loss is the highest
mobilization works control the flow at 37.5% of the loss with 22.260 t/yr. The average value of A
total area of the basin (Tab. 1). The Tahaddart basin (RUSLE) for the entire Tahaddart watershed is 47.79
comprises four sub-basins (SB) (Fig. 3 and Tab. 2): t/ha/yr. This allows deducing that the real annual
Oued Sania SB (with Ibn Battouta dam), Oued average loss per sheet erosion and rill in this
Mharhar SB, Oued Haricha SB (with 9 April 1947 watershed is 54528 t/yr.
dam) and Oued Hachef SB.
Soil Erosion by MUSLE model
Physiographic and Hydrological Parameters The crossing of the MUSLE model factors
(Williams 1975) and the volume of water runoff and
The physiographic parameters show that the instant maximum flows rates in each sub-basin were
Tahaddart basin has an elongated shape with a used to calculate the amount of sediment produced by
Gravilius compactness index of 1.55 over an area of sub-basin in t/yr (Tab. 12). High losses are observed
114100 ha. Its sub-basins are also elongated and in Mharhar and Kharoub (SE of Haricha) sub-basins.
show Kc values between 1.44 and 1.95 (Tab. 2). A These results are consistent with those developed by
Gradex method was used to adjust the series of the RUSLE model. This indicates that erosion is very
maximum rains and the corresponding flows (Table active at the downstream level of the Tahaddart
3). The flows selected from analysis and comparison watershed. Losses calculated for this watershed are
of the flows calculated by various methods (Mallet on average 697894 t/ha, of which 39% are involved
Gauthier, Fuller II, Fuller I, Hazen-Lazarevich) are in the silting up of dams Ibn Battouta and 9 April and
summarized in Table 4. Finally, calculation of the 61% discharged to the Atlantic.
concentration time from empirical formulas
(Formulas of Turazza, Ventura, Kripich, D. Specific degradation
Californian, H. Giandotti, Passini) shows in Table 5, Calculation of the specific degradation rates defined
that it is very long for the sub-basins of Kharoub and by Hession & Shanholz (1988) makes it possible to
Mharhar. This engenders a gathering at the outlet, at estimate the rate of erosion at the outlet of sub-basins
the lakes and sebkhas surrounding Oued Tahaddart. (Tab. 13). The average specific degradation or
quantity of solid transport in the Tahaddart basin is
Assessment and quantification of erosion 30.45 t/ha/yr. This value is close to the bathymetric
surveys conducted at the Ibn Batouta dam (about
The various indices have a spatial character and 36.52 t/yr for 1989-1994, compared with 26.21 t/y in
therefore vary according to the climatic aggression 1994-2000 and 15.91 from 2000 to 2003).
 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental 89

crossing maps of different factors in a GIS enabled

Land losses generated by hydrographic network assessment and quantification of erosion. Potential
The values of soil losses generated by the river erosion is estimated for the Tahaddart watershed at
network or linear erosion for each sub-basin are 202.337 t/ha/yr with more than 40% of the surface of
calculated based on the subtraction of the values sub-basin less than 50 t/ha/yr. The actual average loss
shown by the MUSLE model from those shown by per rill and sheet erosion in this watershed is 47.79
the RUSLE model (Tab. 14). Analysis of the results t/ha/yr, presenting an average level of risk. Total
shows that the losses due to the water system are losses at the outlet are on average 697894 t/ha (ie
about 16.6% for the Tahaddart watershed. Losses on 6.12 t/yr). These losses due to the water system, or
the slopes outweigh the losses due to the water linear erosion, are approximately 16%, confirming
system. that the loss on the slopes, or sheet erosion, outweigh
the losses due to the water system. The average
Empirical modeling: sedimentation (deposition) specific degradation is 30.45 t/ha/yr.
Empirical modeling by RUSLE and Soil losses from RUSLE and Sedimentation show
Sedimentation is used in each sub-basin. It helps that sheet erosion is dominant over the reliefs
clarify the balance of soil losses established in each surrounding the Tahaddart basin. The areas of low
individual parcel. Calculation results show that soil loss or deposition correspond to areas of low
regions with high losses dominate mainly on the sub- slope or plains (up to 50 t/ha/yr of sedimentation).
basin reliefs with 24% of the total area, losses of The results obtained in this study present only
regions between 0.0 to 50 t/ha/yr are concentrated at relative values that can help in the planning of soil
the outlet of sub-basins. Sedimentation dominates in conservation activities and a reduction in siltation of
the rivers of the major rivers. dams. These results provide important support to
decision-makers and developers to simulate scenarios
Conclusions for the evolution of the Tahaddart basin and, above
After careful analysis of the erosion factors all, to plan erosion control interventions.
involved in the loss prediction models for land, the
__________________________________________________________________________________________

INTRODUCTION dimensions (griffes, rigoles, ravines, etc.). Les causes

du phénomène érosif sont: la topographie par ses
Le phénomène d’érosion se manifeste par la deux composantes (longueur et degré de pente), les
dégradation de la terre par l’eau de ruissellement, la types de roches (généralement très friable : schistes et
pluie, le vent, la glace ou autres agents externes, y marnes), le relief ondulé et accidenté de la chaîne
compris les processus tels que l’ablation, le transport rifaine, l’état de la surface du sol (l’encroûtement de
en suspension, le charriage ou les mouvements de surface, la rugosité et la couverture végétale),
masse. A l’échelle nationale, les pertes en terres l’agressivité des pluies et l’opposition des deux
annuelles cumulées sont estimées à 100 millions de saisons, et finalement l’action anthropique. En plus
tonnes (Heusch 1970) ; 50 millions de mètre cube de l’érosion hydrique, sévit aussi une forte érosion
d’eau annuelle perdue par envasement au dépens de éolienne mais reste localisée : on note la formation de
la capacité de stockage (PNABV 1995). A l’échelle dunes sur les sables côtiers à l’Est de la ville
du domaine rifain, les études menées sur une d’Assilah au niveau des sols sableux cultivés. Les
superficie de 4 ha par INYPSA (1987) montrent que impacts négatifs des pertes des terres par l’érosion
la région de Tanger-Tétouan est exposée à un risque dans ces régions, s’aperçoivent sous de multiples
réel élevé ou très élevé de 42,4% et un risque aspects: la réduction de la superficie des sols
d’érosion moyen de 53,6% de la superficie totale de agricoles et la baisse des rendements des cultures,
cette zone. les provinces d’ (Al Hoceima- et de Taza) l’augmentation du taux d’envasement des retenues
présentent 22,9% de risque élevé ou très élevé, mais des barrages, l’aggravation des coulées boueuses
65,3% de risque moyen. Finalement, la région Nord mettant en péril les infrastructures avales, …
Orientale présente, dans l’ensemble, des risques réels La description, la cartographie et la quantification
d’érosion élevés ou très élevés plus modérés des phénomènes d’érosion permettent de bien cerner
représentant 13.7% et un risque moyen important ce problème à l’échelle du bassin et des sous-bassins.
d’érosion réelle (sur 78,2% de sa superficie). Différentes approches peuvent être considérées pour
Dans la région de Tanger, l’érosion hydrique est évaluer et quantifier les taux d’érosion. Certaines sont
le phénomène naturel le plus dominant. Les basées sur l’expertise (choix par l’expert des
différentes formes sont observées : l’érosion en nappe paramètres et de la pondération à appliquer), d’autres
se distingue nettement sur les terrains agricoles situés sur la modélisation (modèles empirique ou physique,
sur des pentes dépassant 10% et l’érosion linéaire globale ou spatialisée, déterministe ou statistique).
exprimée par tous les creusements linéaires qui Parmi ces modèles, certains modélisent le
entaillent la surface du sol suivant diverses formes et transport de la lame de ruissellement et la
sédimentation, tels que ANSWERS (Beasley 1977),
 90 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental

Figure 1. Organigramme de la méthodologie.

Figure 1. Organizational chart of the methodology.

CREAMS (Knisel 1980), SWRRB (Williams et al. Sous-bassin de l’Oued Sania (Tahiri et al. 2014c) et
1985), LISEM (Limburg Soil Erosion Model, De Sous-bassin de l’Oued Haricha (Tahiri et al. 2015).
Roo et al. 1996), KINEROS (KINematic Runoff and La démarche suivit est synthétisée sur la figure 1.
EROSion ; Woolhiser et al. 1970). D’autres
modélisent les processus physiques mis en jeu dans CADRE GEOLOGIQUE ET GEOMORPHOLOGIQUE
l’érosion hydrique tels que : GUEST (Rose et
al.1983, Hairshine & Rose 1992 a et b), WEPP Le bassin versant de Tahaddart, situé au sud de
(USDA Water Erosion Prediction Project, Lane & Tanger (Rif nord occidental, Maroc, Fig. 2), se
Nearing 1989, Nearing et al. 1989), EUROSEM rattache géologiquement au Rif occidental constitué
(Morgan et al. 1998) et SWAT (Arnold et al. 1998). par le chevauchement de plusieurs nappes de
Plusieurs équations ont été établies pour lier les charriage (nappe de Béni-Ider, nappe Numidienne et
facteurs d’érosion entre eux et quantifier les pertes de nappe de Melloussa) reposant sur l’unité externe de
sol (Zingg 1940, Musgrave 1949 in Ronard & Foster Tanger (Suter 1980, Fig. 3). Les terrains sont
1985, Wischmeier & Smith 1978). Plusieurs modèles représentés par des lithologies essentiellement
informatiques ont par la suite vu le jour, dont certains marno-argileuses et gréso-quartzitiques (et flysch)
sont basés sur l’équation universelle de perte de sol d’âges méso-cénozoïques (Fig. 3). Le relief est
(USLE) de Wischmeier et ses versions modifiées polymorphe et la topographie diffère de l’est à
(MUSLE) et révisées (RUSLE). l’ouest. Il est aplati à l’ouest et fortement accidenté à
L’objet de cette étude est d’évaluer les pertes l’Est. C’est un paysage de basses montagnes. Le
annuelles de sol en utilisant les équations universelles point haut de la zone d’étude culmine à 941 m. Il
de perte de sol révisée (RUSLE) et modifiée correspond au sommet de Jbel Sehâne à l’Est du
(MUSLE) et le module Sédimentation, dans le bassin bassin. Les pentes dépassent 35% au niveau des
versant de Tahaddart (Rif nord occidental, Maroc) en escarpements des formations gréseuses ou gréso-
tant que bassin versant de cours d'eau à exutoire vers marneuses se trouvant à l’Est.
l'Atlantique. Elle synthétise et complète nos travaux
de détail dans deux des ses quatre sous-bassins :
Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental 91

N
5 Km

SB Oued Mharhar

Barrage 9 Avril 1947

Figure 2. Carte de répartition des sous bassins (SB) de Tahaddart.

Figure 2. Repartition map of the Tahddart sub- (SB) basins.

Figure 3. Carte géologique simplifiée du Rif nord occidental; tirets: limite du bassin versant de Tahaddart.
Figure 3. Simplified geological map of north western Rif, dashes: Tahaddart watershed limit.
 92 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental

Tableau 1. Les ressources mobilisées dans le bassin versant de Tahaddart (Sources : DGH, 1993).
Table 1. Resources mobilized in Tahaddart watershed (Source: DGH 1993).

Ouvrage Oued Apport moyen Volume utile de Volume régularisé

(Mm3) retenue (Mm3) (Mm3/an)
Ibn Battouta El Kébir 56 10 26
9 avril Haricha 89 270 69
Prise de Bougdour El Kébir 49 - 9,5
Prise de Hachef Haricha 125 - 17,5
Nappe Charf El Akab 2

Tableau 2. Paramètres physiographique des sous-bassins de Tahaddart.

Table 2. Physiographic parameters of the Tahaddart sub-basins.

Indice de Cour d’eau

Pente
Sous- Surface Périmètre Gravelius Kc Altitude Altitude principale
moy.
bassin A (Km²) P (Km) max (m) min (m) Longueur
Valeur Forme SB (%) Nom
(Km)
El Hachef 219 76,48 1.44 Allongée 941 10 4.11 O. Haricha 24.63
Mharhar 319 67 1,95 Allongée 461 10 1,20 O. Mharhar 35,1
Kharoub 426,70 134,52 1,82 Allongée 915 10 1,49 O. Kharoub 70,31
O. Sania 177 69 1,45 Allongée 690 25 2,14 O. Sania 22,94

Tableau 3. Pluies journalières (Pj) de chaque période de retour avec le Gradex pour les différents sous-bassins.
Table 3. Daily rains (Pj) of each return period with the Gradex for the different sub-basins of Tahaddart.

Sous-bassin Station Gradex Pjmax Pjmax Pjmax Pjmax

24 h 10 ans 20 ans 50 ans 100 ans
SB Oued El Hachef 9 avril 36,19 151,31 176,45 209,001 233,39
SB Mharhar Ibn Battouta 27,92 129,18 148,38 173,24 191,86
SB Kharoub Jbel Hbib 33,22 163,67 187,59 218,55 241,76
SB Oued Sania Ibn Battouta 24,95 141,53 159,51 182,79 200,23

Tableau 4. Débit retenus des différentes méthodes Q (m3/s).

Table 4. Debit rates from the different methods Q (m3/s).

Sous-bassin T 10 (ans) T 20 (ans) T 50 (ans) T 100 (ans)

SB d’El Hachef 260 312 379 427

Q (m3/s) SB Mharhar 435 505 594 659

SB Kharoub 365,47 485 577 644
SB Oued Sania 250 300 340 380

contrôlent l’écoulement de la zone d’étude à raison

Sur ce bâti géologique complexe, composé par la de 37,5 % de la superficie total du bassin (Tab. 1).
juxtaposition de terrains à compétence et perméabilité Les eaux colléctées d’une superficie de 706 Km²
différentes et affecté par une importante tectonique (62,5% de la surface totale du bassin) ne sont que
cassante s'installe le réseau hydrographique du bassin partiellement cheminées vers les stations de
de Tahaddart (Tahiri 2014, Tahiri et al. 2014a, traitement et de production d’eau potable par le biais
2014b). L’Oued Tahaddart, le cours d’eau principal de deux prises de dérivation sans aucune capacité de
du bassin, a deux entités principales (Fig. 2) : 1- au stockage (Scapeni et al. 2008). Les eaux déversent
nord l’Oued Mharhar appelé vers l'est du bassin Oued dans l’océan en passant par la zone humide de
El Kébir puis Oued Sania (comportant le barrage Ibn Tahaddart. Le bassin de Tahaddart comprend alors
Battouta); 2- au sud l’Oued El Hachef a deux quatre Sous-Bassins (SB) (Fig. 2 et Tab. 2): le SB de
affluents principaux, l'Oued Haricha (avec le barrage l'Oued Sania (avec le barrage Ibn Battouta), le SB de
"9 Avril 1947") et Oued Kharroub au SE. Ces cours Oued Mharhar, le SB de Oued Haricha (avec le
d’eau prennent naissance à partir de sources barrage "9 Avril 1947") et le SB de Oued Hachef.
sporadiques (au contact des grés Numidien avec les Les variations annuelles du volume et de l'intensité
argilites et schistes). Quatre ouvrages de mobilisation des précipitations, l'humidité et la nature des sols et la
 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental 93

couverture végétale jouent un rôle clé sur l’érosion du engendre un rassemblement à l’exutoire, au niveau
sol (Tahiri et al. 2014c, Tahiri et al. 2015). des lacs et sebkhas qui entourent Oued Tahaddart.

PARAMETRES PHYSIOGRAPHIQUES DES Tableau 5. Temps de concentration (Tc) par différents

SOUS-BASSINS méthodes appliquées aux sous-bassins de Tahaddart
Table 5. Tc by different methods applied to the Tahaddart
Les paramètres physiographiques des sous- sub-basins.
bassins de Tahaddart sont résumés sur le tableau 2.
Ces paramètres sont importants pour le calcul des Sous-Bassin Valeur retenue (h)
pertes en sol. Il en ressort que les sous-bassins SB Oued El Hachef 2,70
diffèrent par leur surface et leur Indice de compacité
de Gravelius et ont tous des formes allongées. La SB Mharhar 13,67
valeur de Kc est comprise entre 1,44 et 1,95. SB Kharoub 7,70

PARAMETRES HYDROLOGIQUES DES SB Oued Sania 2,30

SOUS-BASSINS
DESCRIPTION DES FACTEURS D’EROSION
Ajustement des pluies journalières par la loi de
Gumbel Facteur d’érosivité pluviale (R)

L’application de la méthode Gradex a permis L’érosion en nappe dépend de la hauteur, de

d’ajuster les séries de pluies maximales et les débits l’intensité, de l’énergie cinétique et de la répartition
correspondants (Tab. 3). Les valeurs les plus faibles des pluies L’application de la formule de Masson
pour des Pj max de 10 à 100 ans, sont celle du sous- (1971), qui représente l’intensité moyenne sur 30min
bassin de Mharhar. Ceci est probablement lié à la des gouttes de pluie de chaque averse, à la station de
position de la station Ibn Battouta à l’amant du sous- Kayala a permis à Pepin & Hammouda (2012) de
bassin et non au centre, comme est le cas des autres déduire la valeur R = 54,75 MJ*mm/ha/h comme
sous-bassins. moyenne annuelle des érosivités sur 30 années (1980-
2010). Cependant, en absence de données sur les
Estimation du débit de pointe par des formules l’intensité moyenne des pluies chaque 30 min pour
régionales les autres stations, de nombreuses équations
permettent le calcul de R en fonction des
Les méthodes de calcul recensées en utilisant précipitations mensuelles, annuelles ou de l’indice de
différentes formules empiriques tiennent compte de Fournier (Kalman 1967, Roose 1977, Arnoldus 1980,
deux facteurs essentiels : la taille du bassin versant Renard et al. 1991). Ces équations de substitution ont
(ou sous-bassin) et le gradient topographique du été appliquées dans différents travaux impliquant une
bassin versant. Les différentes formules appliquées modélisation USLE (Morschel & Fox 2004, Sadiki et
(Mallet Gauthier, Fuller II, Fuller I, Hazen- al. 2004, El Gouarani et al. 2008) et ont montré des
Lazarevich) ont permis le calcul des débits de pointe résultats satisfaisants. L’application de l’équation
pour chaque période de retour. Les débits retenus à d'Arnoldus (1980) dans le bassin versant de
partir de l’analyse et la comparaison des débits Tahaddart, a permis de calculer le facteur
calculés par les différentes méthodes présentées d’agressivité climatique aux différentes stations
précédemment sont synthétisés dans le tableau 4. situées au centre et autour des sous-bassins (Tab. 6).

Calcul du temps de concentration Tableau 6. Les valeurs du Facteur R calculé par station par
la formule d’Arnoldus (1980).
Nous rappelons que le temps de concentration Tc Table 6. The R-Factor values calculated per station by the
des eaux sur un bassin versant se définit comme le Arnoldus formula (1980).
maximum de durée nécessaire à une goutte d'eau
Station Erosivité Période
pour parcourir le chemin hydrologique entre un point
(Mjmm/ha/an) (année)
du bassin et l'exutoire de ce dernier. Les valeurs Dar Chaoui 34,37 1941-2012
retenues pour chaque sous-bassin à partir de formules
empiriques (Formules de Turazza, Ventura, Kripich, 09-avril 33,55 1971-2012
D. Californienne, H. Giandotti, Passini) montrent Béni Harchane 35,35 1977-2012
(Tab. 5) montrent que le temps de concentration est Jbel Hbib 46,5 1970-2012
relativement long pour les sous-bassins d’Oued Sania
Kalaya 33,50 1969-2012
et El Hachef, dont l’exutoire est un barrage, ce qui va
engendrer un rassemblement lent des eaux vers Romane 33,97 1977-2012
l’exutoire. Le temps de concentration est très long Ibn Battouta 42,34 1983-2012
pour les sous-bassins de Kharoub et Mharhar. Ceci
 94 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental

La spatialisation des stations a permis de définir Tableau 8. Les moyennes retenues des valeurs du Facteur C
la valeur minimale, la valeur maximale, la moyenne pour chaque sous-bassin.
et l’écart type de l’indice R pour chaque sous-bassin. Table 8. Mean values of the C- Factor for each sub-basin.
Le tableau 7 montre les moyennes retenues pour
chaque sous-bassin. Sous-bassin Facteur C
SB El Hachef 0,21
Tableau 7. Les moyennes retenues des valeurs de l’indice R SB de Mharhar 0,31
pour chaque sous-bassin. SB de Kharoub 0,30
SB de Oued Sania 0,33
Table 7. Mean values of the R- Factor for each sub-basin.

Sous-bassin Facteur R (Mjmm/ha/an) Facteur d’érodibilité du sol (K)

SB de El Hachef 34,42
SB de Mharhar 43,54 L’érodibilité K est une mesure de la vulnérabilité
SB de Kharoub 39,93 des particules de sol au détachement et au transport
SB de Oued Sania 38,96 par la pluie et le ruissellement. Le facteur K peut être
déterminé suivant l’équation de Dumas (1965), la
La cartographie de l’indice R dans le nord du formule ajustée de Wischmeier et Smith (Römkens et
Maroc montre que la valeur moyenne calculée pour le al. in Renard et al. 1997) ou calculée à l’aide du
bassin versant de Tahaddart (R = 39,21 Mjmm/ha/an) triangle de texture (Brown, 2003) et du tableau de
est légèrement inférieure aux valeurs du Rif correspondance de Stone et Hilborn (2000). Les
occidental, central et son avant-pays. Dans le Rif valeurs de K pour chaque unité du sol ont été estimées
occidental l'indice R montre des valeurs variant entre en utilisant les analyses des différents types de sols
94,73 et 95,27 dans le bassin versant de Kalaya (Billaux & Bryssine (1967), Gharbaoui (1981), Alami
(Khali Issa et al. 2016), pour la partie centrale et son (in Osrirhi et al. 2007)). Les valeurs de K se situent
avant-pays. Dans le Rif central, R varie de 43 à 87,56 entre 0,01 (sols peu évolués d’apport colluvial) et 0,4
dans le bassin versant de Oued Lben (Rahhou 1999) (Lithosols). Les moyennes calculées pour chaque
et il est estimé de 51 à 194 dans la vallée de sous-bassin sont représentés sur la figure 5 et
l’Ouergha (Le Landais & Fabre 1975). Dans l’avant- résumées dans le tableau 9. L’indice K pour le bassin
pays oriental, cet indice est estimé à 56,3 dans le versant de Tahaddart varie alors de 0,22 à 0,31.
bassin versant de l’Oued Tleta (El Garouani et al.
2010), et varie de 40 à 50 MJmm/ha/an dans le bassin Tableau 9. Les moyennes retenues des valeurs du Facteur K
versant de l’Oued Boussouab (Sadiki et al. 2004). pour chaque sous-bassin.
Il faut noter que malgré les faibles valeurs de Table 9. Mean values of the K- Factor for each sub-basin.
l'indice d'agressivité des pluies, l'activité de l'érosion
est toujours intense comme en témoignent les Sous-bassin Facteur K
paysages très dégradés et les versants ravinés SB de El Hachef 0,31
observés sur le terrain. Cette situation paradoxale SB de Mharhar 0,28
SB de Kharoub 0,21
semble indiquer que les séries de pluies saturantes et
SB de Oued Sania 0,29
la faible épaisseur du sol sont à l'origine de ces
phénomènes de ravinement bien plus que l'agressivité
des pluies orageuses, très intenses mais limitées dans Facteur des pratiques culturales (P)
l'espace et dans le temps.
La représentation cartographique de la carte de
Facteur du Couvert végétal (C) pente de chaque sous-bassin montre que les pentes
inférieures à 25 degrés sont les plus dominants. Les
Le couvert végétal protège les sols et assure visites sur le terrain confirment l’absence de
l’amortissement des gouttes de pluie, le techniques de pratiques anti-érosives. Ceci nous
ralentissement du ruissellement et de l’infiltration. Il conduit à attribuer la valeur 1 pour ce facteur dans
pourrait varier de 0,001 pour les sols bien couverts chaque sous-bassin. Cette valeur 1 va neutraliser
jusqu’à 1 comme valeur maximale pour les sols l’effet de ce facteur.
complètement dénudés. Pour le bassin de Tahaddart,
ce facteur est déterminé à partir de la mise à jour sur Facteur LS
le terrain de la carte d’occupation du sol établie par
Kouraimi (1997) (Fig. 4). Les valeurs affectées aux La longueur, la forme et surtout l’inclinaison de
différents modes d’utilisation de sols se situent entre pente sont des paramètres qui influencent
0,0 (lacs et habitations rurales) et 0,7 (Céréaliculture considérablement l’érosion des sols (Roose 1994). La
et parcours). Les moyennes calculées pour chaque longueur de pente conditionne la vitesse de
sous-bassin sont résumées dans le tableau 8. L’indice ruissellement et le transport des particules croît en
C pour le bassin versant de Tahaddart varie de 0,22 et fonction de la longueur de la parcelle. De même, les
0,33. transports solides croissent de façon exponentielle
avec le pourcentage de pente (Zingg 1940, Hudson
 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental 95

Figure 4. Carte du facteur Couvert végétal du bassin versant de Tahaddart superposée à un MNT estompé à 315 sur 45°.
Figure 4. Map Vegetation factor of Tahaddart watershed superimposed on a DTM hillshade at N315 -45.

1973, Roose 1994). Les résultats du croisement des Le SIG s’apprête bien à l’évaluation de
couches de l’inclinaison (S) et la longueur (L) de la différentes équations et par conséquent à la
pente pour chaque sous-bassin montrent que les cartographie du taux d’érosion à l’échelle du bassin
valeurs de LS varient de 0,1 à 5. Ces valeurs versant. Nous avons utilisé l’unité de surface de 30 m
coïncident avec les situations les plus fréquentes de x 30 m (dimension du pixel) comme unité de base de
culture en Afrique de l’Ouest. Elles peuvent atteindre l’estimation de chaque paramètre. La superposition
20 en zones de montagne. Le relief à l’aval et au des différentes grilles de données dans le SIG nous
centre du bassin de Tahaddart (Fig. 6), présent une donne la valeur du taux d’érosion au niveau de
morphologie assez aplatie et montre le faible rôle chaque unité de surface.
que pourrait jouer cet indice dans la quantification de
l’érosion. L’érosion potentielle

EVALUATION ET QUANTIFICATION DE L’érosion potentielle désigne le pouvoir d’un sol

L’EROSION sous l’effet des dégradations des gouttelettes des
pluies, dans son environnement morphologique de
Les différents indices ont un caractère spatial et produire des sédiments. C’est une fonction
par conséquent varient en fonction des multiplicative de trois facteurs qui s’introduisent dans
caractéristiques physiques, d’occupation des sols, et la formule de Wishmeier ; elle est décrite comme
des pratiques de conservation au niveau de chaque suit :
unité spatiale du bassin versant.
 96 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental

Figure 5. Carte d’érodibilité du bassin versant de Tahaddart superposée à un MNT estompé à 315 sur 45°.
Figure 5. Map of erodibility of the Tahaddart basin superimposed on a DTM hillshade at N315 - 45°.

Ep = R.K.LS Tableau 10. Pertes potentielles pour chaque sous-bassin du

Avec bassin versant de Tahaddart.
Ep : Erosion potentielle (t/ha/an) Table 10. Potential losses for each sub-basin of the
R : Agressivité pluviale Tahaddart Watershed.
K : Erodibilité des sols
LS : Indice de pentes Sous-bassin Ep. (t/ha/an)
SB de El Hachef 271,25
Le croisement sur un SIG, des cartes synthétiques SB de Mharhar 210,97
SB de Kharoub 180,22
des différents facteurs de l’érosion potentielle a
SB de Oued Sania 146,9
permis de la quantifier. Le tableau 10 synthétise les
différents calculs de l’érosion. Il en ressort que le
L’érosion sur les versants : modèle de Wischmeier
sous-bassin d’El Hachef présente les valeurs les plus
élevées de l’érosion potentielle (271,26 t/ha/an).
Le modèle RUSLE permet d'estimer la perte
Cette zone est caractérisée par des pentes et
moyenne annuelle de sol. C’est une version révisée
longueurs de pentes élevées qui ont conduit à de
de l'équation universelle de pertes de sol (USLE) de
fortes valeurs du facteur topographique LS dont sa
Wischmeier & Smith (1978). Selon Roose & Georges
contribution dans la génération des pertes par le
(2004), l'érosion est une fonction multiplicative de
modèle utilisé est sans doute la plus pesante. La perte
l'érosivité des pluies (le facteur R) que multiplie la
potentielle pour le bassin versant de Tahaddart est de
résistance du milieu, laquelle comprend le facteur
202,337 t/ha/an.
 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental 97

Figure 6. Carte de pente (en pourcent) du bassin versant de Tahaddart superposée à un MNT estompé à 315 sur 45°.
Figure 6. Slope map (in percent) of the Tahaddart basin superimposed on a DTM hillshade at N315 - 45°.

d'érodibilité du sol (K), la topographie (LS), les annuelle par érosion en nappe et rigole dans ce bassin
pratiques antiérosives (P), le couvert végétal (C) et versant est de 54528,39 t/an.
les pratiques culturales (P). L’équation universelle
des pertes en sol est de la forme : Tableau 11. Erosion réelle pour chaque sous-bassin du
A = R*K*LS*C*P bassin versant de Tahaddart.
Table 11. Real erosion for each sub-basin of the Tahaddart
Le croisement des cartes synthétiques des watershed.
différents facteurs du modèle RUSLE a permis de
quantifier l’érosion réelle ou pertes annuelles Sous-bassin A (RUSLE) (t/ha/an)
moyennes en terre en t/ha/an. L’analyse des résultats SB de El Hachef 27,13
SB de Mharhar 64,68
pour les quatre sous-bassins (Tab. 11) montre que le SB de Kharoub 52,17
sous-bassin de Mharhar enregistre la plus forte valeur SB de Oued Sania 47,18
de A(RUSLE) avec 64,68 t/ha/an. La perte réelle
annuelle pour ce sous-bassin est de 20632,92 t/an. Le Ces résultats montrent que les valeurs enregistrées
sous-bassin de Kharroub, quand à lui, présente une dans le bassin de Tahaddart ainsi que les résultats des
valeur A (RUSLE) = 52,17 t/ha/an, mais vue sa surface études antérieures dans le Prérif (El Garouani et al.
qui est plus importante, la perte réelle annuelle est la 2010 , Heusch 1970, Ait Fora 1995, Siteri et al. 2003,
plus forte perte avec 22260,94 t/an. La valeur Laaroussi 2014) sont beaucoup plus élevés que les
moyenne de A(RUSLE) pour l’ensemble du bassin niveaux de tolérance (A(RUSLE) >11 t/ha/an). Cela
versant de Tahaddart est de 47,79 t/ha/an. Ceci montre la gravité de la dégradation spécifique des
permet de déduire que la perte moyenne réelle
 98 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental

sols due au type d’utilisation du sol, la lithologie et le Tableau 13. Dégradation spécifique pour chaque sous-
climat local qui favorisent l’érosion des sols. bassin du bassin versant de Tahaddart.
Table 13. Specific degradation for each sub-basin of the
L’érosion sur les versants : modèle Williams Tahaddart Watershed.

La formule de Williams (1975) (MUSLE) est une Sous-bassin DR (t/ha/an)

version modifiée du modèle USLE et permet SB de El Hachef 26,80
d’attribuer la part de l’érosion à chaque crue. La perte SB de Mharhar 23,85
SB de Kharoub 24,24
moyenne en sol pour une crue est une fonction SB de Oued Sania 46,92
multiplicative du volume de la crue (V en m3), le
débit de pointe de la crue (Qp en m3/s), l'érodibilité
Les pertes en terre générées par le réseau
du sol (K), l’indice de pente (S), la longueur de pente
hydrographique
(L), le couvert végétal (C) et les pratiques culturales
(P). La formule de Williams est estimée de la Les valeurs de la perte en terre générées par le
forme suivante : A = 11,8* (V*Qp) 0,56 réseau hydrographique ou érosion linéaire, pour
*K*L*S*C*P chaque sous-bassin, est calculée en se basant sur
l’hypothèse de soustraction des valeurs ressorties par
Le croisement des cartes synthétiques des le modèle MUSLE de celles ressorties par le modèle
différents facteurs du modèle MUSLE ainsi que les RUSLE (Tab. 14).
volumes d’eau ruisselé et les débits maximums
instantanés dans chaque sous-bassin ont permis de Tableau 14. Pertes en terre pour chaque sous-bassin du
calculer la quantité de sédiments produite par sous- bassin versant de Tahaddart.
bassin en t/an (Tab. 12). Table 14. Land losses for each sub-basin in the Tahaddart
watershed.
Tableau 12. Erosion réelle totale pour chaque sous-bassin
du bassin versant de Tahaddart.
Table 12. Total actual erosion for each sub-basin of the Sous-bassin A (RH) (t/an) P(A) (%)
SB de El Hachef 226 432,65 9,74 %
Tahaddart watershed.
SB de Mharhar 175 951,08 17,87%
Sous-bassin A (MUSLE) (t/an) SB de Kharoub 45 281,005 18,97 %
SB de El Hachef 106 106,25 SB de Oued Sania 441 106,33 19,80 %
SB de Mharhar 196 584,74
SB de Kharoub 176 963,20 L’analyse des résultats montre que les pertes dues
SB de Oued Sania 128 239,80 au réseau hydrographique sont d’environ 16,6% pour
le bassin versant de Tahaddart. Les pertes sur les
Les pertes élevées s’observent au niveau du sous- versants l’emportent sur les pertes dues au réseau
bassin de Mharhar et le sous bassin de Kharoub. Ces hydrographique.
résultats sont conformes à celles élaborées par le
modèle RUSLE. Ceci indique que l’érosion est très Modélisation empirique: sédimentation (déposition)
active au niveau avale du bassin de Tahaddart. Les
pertes calculées pour ce bassin versant sont en La modélisation empirique par RUSLE et
moyenne de 697894 t/ha dont 39 % participent à Sédimentation est utilisée dans chaque sous bassin.
l’envasement des barrages 9 Avril et Ibn Battouta et Elle permet d’apporter des précisions sur le bilan des
61 % déversent vers l’Atlantique. pertes en sol établies dans chaque parcelle
élémentaire. Les résultats des calculs des pertes en
Dégradation spécifique sol annuelles nettes montrent que :
- Les régions de fortes pertes dominent surtout sur les
Le calcul des taux de dégradation spécifique reliefs des sous-bassins. Elles correspondent à des
permet d’estimer le taux d’érosion à l’exutoire des zones de forte pente et occupent 24% de la surface
sous-bassins (Tab. 13). Ce coefficient est défini par totale.
Hession & Shanholz (1988) par la formule - Les régions de pertes entre 0,0 à 50 t/ha/an se
DR=10.D/L où D est la dénivelée du bassin versant concentrent à l’exutoire des sous-bassins. Ces
entre le point culminant et l’exutoire du bassin et L pertes participent à l’envasement des barrages Ibn
est la Longueur du talweg principal. La dégradation Battouta et 9 Avril à hauteur de 40% alors que
spécifique moyenne ou la quantité des transports 60% se déversent vers l’Oued Tahaddart.
solides dans le bassin de Tahaddart est de 30,45 - Les régions à pente faible correspondent à des
t/ha/an. Cette valeur est proche des levés endroits où dominent la sédimentation (ou
bathymétriques réalisés au niveau du barrage Ibn déposition) sur l’érosion. Ces endroits occupent le
Batouta (environ 36,52 t/an pour la période 1989 - centre des sous-bassins,
1994, contre 26,21 t/an en 1994-2000 et 15,91 de - La sédimentation domine dans les cours d’eau des
2000 à 2003). principaux Oueds.
 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental 99

CONCLUSION GENERALE REFERENCES

Le bassin versant de Tahaddart comprend quatre Ait Fora M. 1995. Modélisation spatiale de l’érosion
sous-bassins dont deux sont en amont de deux hydrique dans un bassin versant du Rif marocain:
barrages : barrage 9 avril et barrage Ibn Battouta. Le Validation de l’approche géomatique par la
bassin montre une forme allongée avec un indice de sédimentologie, les traceurs radioactifs et la
compacité de Gravilius est de 1,55 sur une superficie susceptibilité magnétique des sédiments. PhD,
Université de Sherbrooke, Québec, 200p.
de 114100 ha. Le réseau hydrographique est ramifié.
Arnold J.G., Srinivasan R., Muttiah R.S. et al. 1998. Large
Les principaux cours d’eau (Oued El Hachef et Oued Area Hydrologic Modeling and Assessment Part I:
Mharhar) s’allongent de l’est vers l’ouest sur des Model Development. Journal of American Water
distances respectives de 650 et 480 km et des pentes Resources Association, 34, 73-89.
moyennes de 1,20 à 4,10 %. Une telle situation Arnoldus H. M. J. 1980. An Approximation of the Rainfall
pouvant favoriser le développement des écoulements Factor in the Universal Soil Loss Equation. In: De
fluviatiles et les dépôts des sédiments au niveau des Boodt, M. and Gabriels, D., Eds., Assessment of
lits des Oueds. Le temps de concentration, déduit des Erosion. FAO Land and Water Development Division,
données pluviométriques est en moyen de 2,50 h pour Wiley & Sons, Chichester, 127-132.
Beasley D. B. 1977. ANSWERS. A Mathematical Model
les sous-bassins d’Oued Sania et El Hachef, dont
for Simulating the Effects of Land Use and
l’exutoire est un barrage, ce qui va engendrer un Management on Water Quality. Ph.D. Thesis, Purdue
rassemblement lent des eaux vers l’exutoire. Le University, West Lafayette, Indiana. 220p.
temps de concentration est très long pour les sous- Billaux P. & Bryssine G. 1967. Les sols du Maroc. In :
bassins de Kharoub et Mharhar. Ceci engendre un Congrès de pédologie méditerranéenne: Excursion au
rassemblement à l’exutoire, au niveau des lacs et Maroc. Cahiers de la Recherche Agronomique, 1, 59-
sebkhas qui entourent Oued Tahaddart. Après une 101.
analyse minutieuse des facteurs d’érosion entrant Brown R.B. 2003. Soil Texture, Soil and Water Science
dans les modèles de prédiction des pertes en terres, le Department, Florida Cooperative Extension Service,
Institute of Food and Agricultural Sciences. University
croisement des cartes synthétiques des différents
of Florida, Fact Sheet SL29, 8p.
facteurs dans un SIG a permis une évaluation et une Commission de Pédologie et de Cartographie des Sols
quantification de l’érosion. L’érosion potentielle est (CPCS). 1967. Classification des sols. Laboratoire de
évaluée pour le bassin versant de Tahaddart à Géologie et de Pédologie, Ecole Nationale Supérieure
202,337 t/ha/an avec plus de 40% de la surface des d’Agronomie, Grognon, France, 87p.
sous-bassins a moins de 50 t/ha/an. La perte moyenne De Roo A.P.J., Wesseling C.G., Ritsema C.J. 1996.
réelle par érosion en nappe et rigole dans ce bassin LISEM: a single event physically-based hydrologic and
versant est de 47,79 t/ha/an, présentant un niveau de soil erosion model for drainage basins: I. Theory, input
risque moyen. Les pertes totales à l’exutoire sont en and output. Hydrological Processes, 10, 8, 1107-1117.
DGH. 1993. Les ressources mobilisées dans le bassin de
moyenne de 697894 t/ha (soit 6,12 t/an). Ces pertes
Tahaddart. Rapport de la Direction Générale de
dues au réseau hydrographique, ou érosion linéaire, l’Hydraulique, Maroc. 56p.
sont d’environ 16% et confirment que les pertes sur Dumas J. 1965. Relation entre l'érodibilité des sols et leurs
les versants, ou érosion en nappe ou en rigole, caractéristiques analytiques O.R.S.T.O.M. Tunis, sér.
l’emportent sur les pertes dues au réseau pédologie, 3, 307-333.
hydrographique. La dégradation spécifique moyenne El Garouani A., Tribak A. & Abahrour M. 2010.
est de 30,45 t/ha/an. Les pertes en sols établis par Assessment the Effects of Land-Use Cover Changes on
RUSLE et Sédimentation montrent que l’érosion en Regional Soil Loss Susceptibility Using RUSLE Model
nappe est dominante sur les hauts versants du bassin and Remote Sensing Data. Red Books. International
Association of Hydrological Sciences Publication, 340,
versant de Tahaddart. Les régions de faible pertes en
343-349.
sol ou de déposition correspondent à des zones de El Gharbaoui A. 1981. La terre et l’homme dans la
faible pente ou à des plaines (jusqu’à 50 t/ha/an de péninsule tingitane: Etude sur l’homme et le milieu
sédimentation). naturel dans le Rif Occidental. Travaux de l’Institut
Les résultats obtenus dans cette étude ben qu’elle Scientifique, série Géologie Géographie Physique, 15,
ne présentent que des valeurs relatives, elles sont 1-439.
précieuses, Elles peuvent servir de base pour aider ESRI. 1999. ArcView GIS 3.2. Environmental Systems
dans la planification des activités de conservation des Research Institute, Inc., Redlands.
sols et une réduction de l’envasement des barrages. Hairshine P.B. & Rose C.W. 1992a. Modeling water
érosion due to overland flow using physical principles.
Ces résultats apportent une aide importante aux
1. Sheet flow. Water Resources Research, 28, 237-243.
décideurs et aménageurs pour simuler des scénarios Hairshine P.B. & Rose C.W. 1992b. Modeling water
d’évolution du bassin de Tahaddart et surtout érosion due to overland flow using physical principles.
planifier les interventions de lutte contre l’érosion. 2. Rill flow. Water Resources Research, 28, 245-250.
Hession W.C. & Shanholtz V.O. 1988. A geographic
REMERCIMENTS information system for targeting nonpoint-source
agricultural pollution. Journal of Soil Water Conserv.
Les auteurs remercient les évaluateurs anonymes pour leurs 43, 3, 264-266.
remarques et suggestions.
 100 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental

Heusch B. 1970. L’érosion dans le Prérif: Une étude Pepin Y. & Hammouda N. 2012. Variabilité de l’Intensité
quantitative de l’érosion hydraulique dans les collines des pluies dans la région de Tanger sur une période de
marneuses du Prérif occidental. Annales des recherches 30 ans. Revue marocaine des sciences agronomiques et
forestières, 12, 9-176. vétérinaires, 1, 23-27.
Hudson N. 1973. Soil Conservation. London, B.T. Batsford, PNABV. 1995. Rapport du Programme National
c1, 320p. d'Aménagement des Bassins Versants, Maroc.
Inypsa 1987. Etude des sols au 1/100.000 (Edition au Rahhou M. 1999. L’érosion dans le Prérif central, zone
1/50.000). Projet Intégré de développement Agricole de interfluviale Leben-Sebou-Ouergha, un prolongement
Tanger-Tétouan, Secteur de Tétouan. Inypsa -Maroc, de l’évolution naturel, une production sociale. Thèse
S.A. et Direction Provinciale de l’Agriculture de Université Mohammed V, Rabat, Morocco, 250p.
Tétouan. Renard K.G. & Freimund J.R. 1994. « Using monthly
Kalman R. 1967. Essai d’évaluation pour le pré-Rif du precipitation data to estimate the R factor in the revised
facteur couverture végétale de la formule de USLE », Journal of Hydrology, 157, 287-306.
Wischmeier de calcul de l’érosion. In: Rapport pour Renard K.G., Foster G.R., Weesies G.A. et al. 1991.
l’administration de la forêt et d’eau, Rabat, 1-12. RUSLE: Revised Universal Soil Loss Equation.
Khali Issa1 L., Ben Hamman Lech-Hab K., Raissouni A. et Journal of Soil and Water Conservation, 46, 30-33.
al. 2016. Quantitative Mapping of Soil Erosion Risk Renard K.G., Foster G.R., Weesies G.A. et al. 1997.
Using GIS/USLE Approach at the Kalaya Watershed Predicting Soil Erosion by Water: A Guide to
(North Western Morocco). Journal of Materials and Conservation Planning with the Revised Universal Soil
Environmental Sciences, 7, 8, 2778-2795. Loss Equation (RUSLE). US Department of
Knisel W.G. 1980. Erosion and Sediment Yield Models-- Agriculture, Washington DC.
An Overview, Proceedings of ASCE Symposium on Renard K.G., Foster G.R., Weesies G.A. et al. 1996.
Watershed Management, Boise, Idaho, 1, 141-150. Predicting soil erosion by water: A guide to
Kouraimi B. 1997. Cartographie et analyse de la conservation planning with the revised universal soil
production pastorale et ligneuse du matorral dans le loss equation (RUSLE). Agriculture Handbook Number
Rif Occidental : Cas du Bassin Versant de Tleta. 703, Agricultural Research Service, United States
Mémoire de 3ème cycle. Pastoralisme. Institut Department of Agriculture, Tucson, AZ, USA, 101-
Agronomique et Vétérinaire Hassan II, 50p. 142.
Laaroussi O. 2014. Etude de l’érosion et de la déposition Roose E. & Georges D. N. 2004. Recherches sur l‘érosion
des sols à l’échelle du bassin versant de l’Oued El hydrique en Afrique : revue et perspectives. Revue
Maleh (Fès, Maroc) par l’utilisation conjointe de la Sécheresse, 15, 1, 12-19.
télédétection, du SIG et des techniques de mesures sur Roose E. 1977. « Application of the Universal Soil Loss
le terrain. Doctorat, Faculty of Sciences and Equation of Wischmeier and Smith in West
Technology, Fès, Morocco, 141p. Africa », Soil Conservation Society of America,
Lane L. J. & Nearing M. A. eds. 1989. USDA Water Ankeny, Iowa, 50-71.
Erosion Prediction Project: Hillslope profile model Roose E. 1994. Introduction à la gestion conservatoire de
documentation. NSERL Report No. 2. West Lafayette, l’eau, de la biomasse et de la fertilité des sols (GCES).
Ind.: USDA‐ARS National Soil Erosion Research Bulletin pédologique de la FAO N°70, Rome. 420 p.
Laboratory. Rose C.W., Williams J.R., Sander G.C et al. 1983. A
Le Landais F. & Fabre G. 1975. Plan d’aménagement anti- mathematical model of Soil érosion and déposition
érosif du bassin versant de l’oued Ouergha en amont du process: I.Theory for a plane land element. Soil Science
barrage El Ouahda. Rapport inédit, Administration des Society of America, 47, 991-995.
Eaux et Forêts et de la Conservation des Sols. Sadiki A., Bouhlassa S., Auajjar J. et al. 2004. Utilisation
Masson J.M. 1971. L'Érosion des sols par l'eau en climat d’un SIG pour l’évaluation et la cartographie des
méditerranéen. Méthodes expérimentales pour l'étude risques d’érosion par l’Equation universelle des pertes
des quantités érodées à l'échelle du champ. Thèse de en sol dans le Rif oriental (Maroc): Cas du bassin
doctorat, Université des Sciences et techniques du versant de l’oued Boussouab. Bulletin de l’Institut
Languedoc-Montpellier 2, , 213p. Scientifique, Rabat, section Sciences de la Terre, 26,
Morgan R.P.C., Quinton, J.N., Smith, R.E. et al. 1998. The 69-79.
European Soil Erosion Model (EUROSEM): A Siteri H., Aït Brahim L., Sossey A.F. et al. 2003.
Dynamic Approach for Predicting Sediment Transport Intégration du modèle USLE dans un SIG pour la
from Fields and Small Catchments. Earth Surface détermination de la dégradation spécifique des sols
Processes and Landforms, 23, 527-544. dans le Rif septentrional (cas du bassin versant de
Morschel J. & Fox D. 2004. « Une méthode de l’oued Nakhla). Travaux de l’Institut Scientifique de
cartographie du risque érosif : application aux collines Rabat, Série Géologie et Géographie Physique, 21,
du Terrefort Lauragais », M@ppemonde, 76, 11 p. 241-276.
http://mappemonde.mgm.fr/num4/articles/art04404.pdf Stone R.P. & Hillborn D. 2000. Universal Soil Loss
Nearing M. A., Foster G. R., Lane L. J. et al. 1989. A Equation. Ministry of Agriculture and Food
process‐based soil erosion model for USDA Water (OMAFRA) Ontario, Canada. www.gov.on.ca/omafra/
Erosion Prediction Project technology. Trans. ASAE frech/engineer/facts/00-002.htm.
32, 5, 1587‐1593. Suter G. 1980. Carte structurale de la chaîne rifaine au
Osrirhi A., El Oumri M., Moussadek R. et al. .2007. 1/500 000. Notes et Mémoires du Service géologique du
Vocation agricole des terres de la zone d’Oujda, Maroc, 245b.
Rapport et cartes. Rapport Institut National de la Tahiri M. 2014. Contribution de l’outil de l’observation de
Recherche Agronomique. www.inra.org.ma/environ/ la Terre à l’analyse du milieu environnemental du
docs/cvat/oujda/oujda.pdf. bassin de Tahaddart (Rif nord occidental, Maroc).
 Tahiri et al.- Quantification de l’érosion hydrique et de la sédimentation, Tahaddart, Rif nord occidental 101

Doctorat, Hassan II University of Casablanca, Sub-Basin (Tahaddart Watershed, Western Rif,

Morocco, 256p. Morocco): Risk Assessment. Journal of Geoscience
Tahiri M., Achab M., Emran A. et al. 2014a. Lithology and Environment Protection, 4, 107-119.
Data Contribution in Hydrographic Network Williams J. R., Nicks A. D. & Arnold J. G. 1985. Simulator
Distribution Using Remote Sensing and GIS: Case of for Water Resources in Rural Basins. Journal of
the Tahaddart Basin, Northwestern Rif, Morocco. Hydrologic Engineering, ASCE 3, 6, 970-986.
International Journal of Advanced Research, 2, 380- Williams J. R. 1975. Sediment-Yield Prediction with
391. Universal Equation Using Runoff Energy Factor,
Tahiri M., Hakdaoui M., Emran, A. et al. 2014b. Present and Prospective Technology for Predicting
Contributions du SIG et de la Télédétection à la Sediment Yields and Sources, ARS-S-40, USDA,
corrélation entre le NDVI et la pluviométrie pour Agricultural Research Service, Washington D. C., 244-
l'étude de la variabilité climatique dans le bassin de 252.
Tahaddart (Rif, Maroc). European Journal of Scientific Wischmeier W.H. & Smith D.D. 1978. Predicting Rainfall
Research, 122, 253-274. Erosion Losses- A Guide to Conservation Planning.
Tahiri M., Tabyaoui H., El Hammichi F. et al. 2014c. USDA Handbook No. 537, U.S. Department of
Evaluation et Quantification de l’Erosion et la Agriculture in Cooperation with Purdue Agricultural
Sédimentation à Partir des Modèles RUSLE, MUSLE Experiment Station.
et Déposition Intégrés dans un SIG. Application au Woolhiser D. A., Hanson C. L. & Kuhlman R. 1970.
Sous-Bassin de l'Oued Sania (Bassin de Tahaddart, Rif Overland Flow on Rangeland Watersheds, New
nord occidental, Maroc). European Journal of Zealand. Journal of Hydrology, 9, 2, 336-356.
Scientific Research, 125, 157-178. Zingg A. W. 1940. Degree and Length of Land Slope as it
Tahiri M., Tabyaoui H., Tahiri A. et al. 2015. Modelling Affects Soil Loss in Runoff. Agricultural Engineering,
Soil Erosion and Sedimentation in the Oued Haricha 21, 2, 59-64.

Manuscrit reçu le 30/06/2017

Version révisée acceptée le 28/12/2017
Version finale reçue le 29/12/2017
Mise en ligne le 29/12/2017