304 Choukri et al.

: Analyse du fonctionnement hydro-sédimentaire d’un bassin versant

Analyse du fonctionnement hydro-sédimentaire du bassin versant

Tleta du Rif Occidental du Maroc à l’aide du modèle SWAT
F. CHOUKRI1, M. CHIKHAOUI1, M. NAIMI1, Y. PEPIN1,2, D. RACLOT1,2
(Reçu le 13/12/2018; Accepté le 06/03/2018)

Résumé
Le modèle SWAT (Soil and Water Assessment Tool) est un outil de simulation des flux dans le paysage, très utile pour orienter
la gestion et les aménagements des bassins versants. Dans cette optique, ce modèle a été testé pour reproduire, à un pas de temps
journalier, les flux d’eau et de sédiments arrivant à la retenue du barrage Ibn Batouta pour la période 1980-2010. Diverses don-
nées climatiques (précipitation et température), morphologiques (MNT) et thématiques (occupation du sol et pédologie) ont été
utilisées pour l’alimentation du modèle. Les mesures acquises au niveau de la retenue Ibn Batouta sur les périodes 1980-1982,
1983-1996 et 1997-2010 ont été utilisées respectivement pour les étapes d’initialisation, de calibration et de validation après
avoir exécuté une analyse de sensibilité pour déterminer les paramètres les plus sensibles. Les indices de performance R², NSE
et PBIAS ont été adoptés pour évaluer et contrôler la qualité des prédictions des débits mensuels liquides et solides. En effet
la performance de prédiction est très élevée pour les débits liquides mensuels avec R² (0,95; 0,90), NSE (0,94; 0,89) et PBIAS
(-6; 2) respectivement pour la période de calibration et la période de validation. La performance des prédictions des transports
solides mensuels est également très satisfaisante avec R² (0,83; 0,64), NSE (0,82; 0,61) et PBIAS (17; -6). Au final, le modèle
SWAT permet de reproduire correctement les flux mensuels d’eau et de sédiments à l’exutoire du bassin versant Tleta et son
utilisation pour tester des scénarios de gestion et d’aménagements, peut donc être envisagée.
Mots clefs: SWAT, Modélisation hydrologique, Érosion hydrique, Rif Occidental, Maroc.

Analysis of the hydro-sedimentary functioning of Tleta watershed of the Western Rif of Morocco
using the SWAT model
Abstract
SWAT Model (Soil and Water Assessment Tool) is a tool for simulating flows in the landscape, useful for guiding management
and ground works within catchment areas. Accordingly, this model was tested to reproduce, on a daily times pan, water and
sediment flows arriving to Ibn Batouta reservoir over the period 1980-2010. Various climatic (precipitation and temperature),
morphologic (MNT) and thematic data (land use and pedology) were used for feeding the model. The measurements acquired
at Ibn Batouta reservoir during the periods 1980-1982, 1983-1996 and 1997-2010 were used respectively for initialization,
calibration and validation after running the sensitivity analysis to determine the most significant parameters. The performances
indices of R², NSE and PBIAS were adopted to assess and control the quality of predictions of the monthly flows of water and
sediments. The performance of prediction of monthly water flows is very high with R² (0.95; 0.90), NSE (0.94; 0.89) and PBIAS
(-6; 2) respectively for both periods of calibration and validation. Meanwhile the performance of predictions of the monthly solid
transport is very satisfactory with R² (0.83; 0.64), NSE (0.82; 0.61) and PBIAS (17; -6) respectively. Finally, SWAT model allows
to properly reproducing monthly water and sediment flows within the Tleta catchment area and testing scenarios of management
and ground works to be considered.
Key words: SWAT, hydrological modeling, water erosion, Western Rif, Morocco.

Introduction aussi la durabilité des infrastructures en aval (barrages,

canaux, etc.). Partout dans le monde, l’érosion hydrique
Au Maroc, la région du Rif constitue une zone assurant présente une forme de dégradation dangereuse à l’égard de
une partie importante des ressources en eau du pays. la protection des ressources en eau et en sol. Selon Yang et
Malheureusement, l’érosion hydrique y est très impor- al. (2009), l’érosion du sol peut conduire à une réduction
tante du fait d’une combinaison de facteurs défavorables : de la fertilité des sols, la perte de nutriments et à la baisse
une topographie accidentée, un réseau hydrographique des rendements des cultures dans les terres agricoles.
très ramifié et très encaissé, une lithologie marneuse De même, Bakker et al. (2005 et 2004) ont trouvé qu’en
imperméable et donc très érodable, un régime de pluie moyenne l’érosion des sols réduit la productivité des
comprenant des averses avec de fortes intensités ou de cultures d’environ 4 % pour chaque 10 cm de sol perdu.
longue durée et une activité anthropique très marquée, par L’érosion excessive peut épuiser les réserves du sol, en
exemple en mettant à nu de très grandes étendues de sols azote et en phosphore, nécessaires pour la croissance des
à l’automne lors des labours. Cette forme de dégradation plantes, et l’érosion extrême peut dégrader le sol au point
des sols réduit fortement la performance agricole en amont qu’il sera incapable de soutenir la vie végétale (Neitsch
du fait de la perte de la partie la plus fertile du sol, mais et al., 2005).
1
Institut Agronomique et Vétérinaire Hassan II, Rabat, Maroc
2
IRD-UMR LISAH, Montpellier, France
©
Revue Marocaine des Sciences Agronomiques et Vétérinaires www.agrimaroc.org
 Rev. Mar. Sci. Agron. Vét. (2019) 7 (2): 304-312 305

Au Maroc, les pertes annuelles du sol à cause de l’érosion, pertinence pour le contexte nord-marocain du Rif doit
mesurées à l’exutoire de bassins versants de plusieurs être confirmée.
centaines de km² atteignent fréquemment 500 t km-² dans
Ce papier présente la mise en œuvre, le calage et la
le Moyen Atlas et parfois plus de 5 000 t km-² dans le Rif.
validation du modèle SWAT dans le bassin versant Tleta
La première conséquence est une chute progressive de
(180 km2) grâce à des données de suivi des flux d’eau et
la productivité des sols et une amplification de la dégra-
de sédiments acquises au niveau du barrage Ibn Batouta
dation des sols (15 millions d’hectares de terres érodées
durant la période 1980-2010. Ce travail de test de la perti-
sont concernés). Une seconde conséquence porte sur les
nence de l’outil de simulation SWAT doit être perçu comme
ressources en eau puisque l’envasement des retenues de
une première étape vers la possibilité de tester des scénarios
barrage est annuellement de l’ordre de 75 millions de m3,
d’aménagement et de conservation des eaux et des sols en
soit une réduction annuelle de 0,5 % de la capacité de
vue de lutter contre les changements globaux à venir.
stockage de l’eau par les barrages. Cela se traduit par des
pénuries d’eau ou une détérioration de la qualité de l’eau
qui entraînent une augmentation des frais de fonctionne- Méthodologie
ment et de gestion des barrages, des usines hydro-élec- Description du site
triques et des unités de traitement des eaux (http://www.
eauxetforets.gov.ma/fr/text.aspx?id=1025&uid=44). Situé au cœur du Tangérois, à mi-chemin entre deux
L’étude réalisée dans le Rif Occidental (Choukri et al., grandes villes marocaines de Tanger et de Tétouan (Nord
2016) sur la prédiction de l’évolution de l’érosivité des du Maroc), le bassin versant Tleta couvre une superficie
pluies selon des prévisions climatiques fournies par quatre de 180 Km² (Figure 1) à l’amont du barrage Ibn Batouta
MCGs et deux scénarios d’évolution du climat (A1B et construit en 1977 avec une capacité initiale de stockage de
A2) montre une faible évolution de ce facteur à l’échelle 43,5 Mm3. Les altitudes du bassin varient entre 23 et 672 m
annuelle, mais une forte évolution selon les saisons avec au point culminant avec des pentes oscillant entre 0 et
une diminution du facteur R en hiver et au printemps et 30 %. C’est une zone caractérisée par un climat atlantique
une augmentation marquée pendant l’été et l’automne. sous influence méditerranéenne, chaud et sec en été, froid
Les pluies exceptionnelles de la fin de l’été et du début et pluvieux en hiver, avec une période sèche entre mai et
de l’automne peuvent être classées comme dangereuses septembre et une période humide qui s’étale d’octobre à
pour les sols. Durant cette période, les sols de cultures mars. Le bassin versant est exposé à une double influence
annuelles sont nus et les précipitations, généralement sous marine (Mer Méditerranée et Océan Atlantique).C’est
forme d’orages violents, provoquent un ruissellement l’une des régions les plus arrosées du pays, connaissant
intense responsable d’une érosion des sols importante en ainsi une activité hydrologique très importante. Les pré-
nappe et en ravins (Sadiki, 2005). Il est donc important cipitations moyennes annuelles sont de 778  mm sur la
de réfléchir à la définition d’une stratégie de planification période 1979-2013 et la température moyenne est de 18°C
et de conservation des ressources en eau et en sol pour la sur la période 1976-2016, alors que l’évaporation au bac
région, qui permettent d’éviter une augmentation des flux Colorado est estimée à 1 690 mm sur la période 1983-2013.
érosifs et par conséquent une accélération de l’envasement Le bassin se caractérise par une dominance du substrat
des barrages à l’aval. marneux, ce qui rend le sous-sol très imperméable. La
lame d’eau écoulée annuellement est de 345 mm pour la
Actuellement, de nombreux modèles hydrologiques sont période 1979-2013, soit un coefficient de ruissellement
disponibles pour l’évaluation de l’impact de l’utilisation moyen interannuel de 0,43. Le bassin réagit rapidement
des terres et celui des changements climatiques sur la aux précipitations, le module de débit interannuel est de
quantité et la qualité des eaux mobilisées dans les bassins 1,56 m3s-1. Les caractères litho-stratigraphiques généraux
versants. Parmi ces modèles, le modèle agro-hydrologique de la chaîne rifaine sont bien représentés dans le bassin, à
SWAT (Soil and Water Assessment Tool)(Arnold and savoir le chevauchement de plusieurs nappes de charriage
Fohrer, 2005) est un modèle prédictif capable de simuler, (Flysh) constituant la ligne de crêtes au-dessus de l’unité
sur de longues périodes de temps, les effets de différents autochtone de Tanger (marneuse) formant les collines
scénarios de gestion du territoire sur les transferts hy- (Michard, 1976).Cette variété lithologique a permis le
driques et sur la mobilité des sédiments et des nutriments développement d’une mosaïque pédologique assez impor-
à l’échelle des grands bassins versants (Beaudin et al., tante: sols peu évolués, sols hydromorphes, lithosols, ver-
2006). Il est issu d’une série de modèles développés anté- tisols et sols brunifiés. Le couvert végétal est conditionné
rieurement (Arnold et al., 2012), notamment les modèles lui aussi par cette disposition morpho-pédologique, en
SWRRB (Williams et al., 1985), EPIC (Williams et al., effet, les faciès gréseux et flyshiques constituent le substrat
2008), CREAMS (Knisel, 1980) et GLEAMS (Leonard et d’une végétation de matorral parfois très dégradée, alors
al., 1987). SWAT a été appliqué sur des bassins agricoles que le faciès marneux est essentiellement cultivé et que
partout dans le monde pour analyser des problèmes de la forêt se limite à des surfaces restreintes à l’amont du
pollution, d’érosion ou d’aménagements (Chaponnière, bassin (Merzouk et al., 1996). Le couvert végétal a connu
2005). S’il a été testé et validé dans des contextes et à des une dégradation au cours de la période 1977-2009 avec
échelles spatio-temporelles très variées (Lévesque, 2007; une baisse de la surface du bassin classée en végétation
Ouessar et al., 2008; Laurent et al., 2007; Hunink et al., dense et moins dense de 46 à 34 %, au profit d’une pro-
2013; Bonumá et al., 2014; Zhang et al., 2008; Kim et al., gression des terrains de cultures (51,3 à 60 %) et des sols
2010; Mosbahi et al., 2011; Ficklin, 2010; Salou, 2009; nus (2,7 à 6 %).
Gatzke, 2010; Yang et al., 2008; Ficklin et al., 2009), sa
306 Choukri et al.: Analyse du fonctionnement hydro-sédimentaire d’un bassin versant

Données disponibles et préparation des entrées du du bassin versant Tleta et représentant de nombreux types
modèle de sols, associés au rapport pédologique de la région
Tanger-Tétouan (Inypsa, 1987), couvrant l’ensemble du
Données de contexte bassin versant Tleta et caractérisant l’ensemble des types
Plusieurs classes d’occupation du sol ont été reconnues de sols ont été mobilisés. Ce rapport pédologique nous
au niveau du bassin versant étudié (Figure 1a) grâce a permis de construire la base de données sol qui a été
à la classification supervisée d’une série des images: complétée par des échantillons prélevés sur le terrain.
LandSat MSS (90 m) de 1977 et LandSat ETM (30 m) Enfin, un modèle numérique du terrain (MNT) à 20 m de
de 1984, 1988, 2000, 2003, 2007 et 2009 téléchargées résolution issu d’images SPOT a été utilisé pour rendre
du site (http://earthexplorer.usgs.gov/), contrôlée par des compte de la morphologie de la zone d’étude.
vérités terrain. L’occupation du sol de 2009 a été utilisée
Données climatiques
comme référence pour le paramétrage du modèle. Le bas-
sin versant Telta est caractérisé par la prédominance des Le dispositif de mesures du climat a varié au cours du
terrains de culture, mais possède aussi de larges zones de temps. Il est composé de trois pluviomètres journaliers
matorral dense et clair et de sols nus. Les caractéristiques répartis à l’amont du bassin versant et d’une station clima-
physico-chimiques et hydrologiques des sols présentent, tique au niveau de la digue du barrage, relevée à 7 heures du
autre entrée importante du modèle SWAT, ont été définies matin et composée d’un pluviomètre journalier, d’un plu-
à l’échelle de l’unité pédologique. Pour cela, la carte viographe à auget basculeur de 0,1 mm, d’un bac d’évapo-
pédologique (1/50 000) (Figure 1b) couvrant l’ensemble ration type Colorado, de deux thermomètres (max. – min.),

a b

Figure 1. Localisation du bassin versant Tleta avec les sous-bassins et le modèle numérique de la zone. (a) Carte d’occupa-
tion du sol du bassin (2009). (b) Carte pédologique du bassin (CPCS, 1967)
 Rev. Mar. Sci. Agron. Vét. (2019) 7 (2): 304-312 307

d’un anémomètre et d’une girouette. Les températures démarre quand Qe devient positif et finit quand il devient
maximales et minimales sont lues sur les thermomètres nul, sauf si une autre crue vient la chevaucher, dans ce
(à 0,1°C), chaque température journalière est la moyenne cas nous arrêtons la première quand l’autre démarre. De
entre le maximum et le minimum du jour. Les données chacune des crues isolées, nous extrayons le volume total
climatiques exigées par SWAT sont la pluviométrie et la de la crue (V) et son débit maximal (Qmax).
température journalière (min. et max.). D’autres données
météorologiques optionnelles dépendent de la méthode Données de débits solides
de l’estimation de l’évapotranspiration potentielle (ETP).
Des campagnes bathymétriques sont réalisées environ
SWAT propose trois méthodes (Chaponnière, 2005) pour
tous les 5 ans, elles permettent d’établir les courbes Hau-
estimer l’ETP (Hargreaves, 1975 ; Priestley-Taylor, 1972
teur/Surface/Volume (H, V, S) utiles pour le calcul du
et Pennman-Monteith, 1948). Les données climatiques
bilan en eau du barrage, mais aussi pour le calcul des taux
(Pluie, Tmin., Tmax.) ont été prises sur les stations météo-
d’envasement de la retenue. Huit mesures de bathymétrie
rologiques gérées par l’Agence du Bassin Hydraulique du
sont disponibles pour le barrage Ibn Batouta (Tableau 1).
Loukkos (ABHL) ou par la Direction de la Météorologie
Nationale (DMN) (Tableau 1). Les flux des sédiments utilisés pour caler le modèle SWAT
sont issus des mesures bathymétriques et de l’analyse
Tableau 1. Données météorologiques utilisées
de la réponse hydrologique détaillée de chaque crue. Ils
Nom Latitude N Longitude W Mesures correspondent à la somme des débits solides piégés par la
retenue et des débits solides qui n’ont pas été retenus au
Béni Harchane 35,531 -5,720 Pjour sein de la retenue. Pour chaque crue, le volume de sédi-
Ibn Battouta 35,645 -5,733 Pjour ments piégés dans la retenue (Apiégé en m3) a été calculé
en déconvoluant le volume de sédiments accumulés dans
Kalaya 35,673 -5,747 Pjour la retenue entre deux mesures d’envasement (période de
Romane 35,704 -5,656 Pjour 5 ans) en fonction du volume global de la crue (V en m3)
Aéroport-Tanger 35,726 -5,906 Tjour
et du débit maximal instantané de la crue (Qmax en m3s-1)
à l’aide de l’équation de Williams (1975) (Equation 1).
Données de débits liquides Apiégé~ (V Qp)0,56 (eq. 1)
Un suivi du niveau d’eau dans la retenue du barrage Le flux de sédiment piégé par la retenue est ensuite déduit
accompagné du suivi des vidanges et des déversements en multipliant le volume de sédiments Apiégé par la densité
est réalisé au pas de temps horaire. La reconstitution des moyenne des sédiments de la retenue estimée à 1,2.
débits entrants dans la retenue est faite au pas de temps
horaire en utilisant le logiciel Hydraccess. Le calcul du On ajoute ensuite les flux de sédiments non piégés par le
débit entrant repose sur la conservation de l’eau à travers barrage. Pour évaluer ces flux de sédiments non piégés,
une équation de bilan (Albergel et al.,1999) qui fait inter- les flux de vidange ONEP sur la période 1980-2010 sont
venir les apports directs de la pluie sur la surface en eau supposés avoir des concentrations négligeables en sédi-
de la retenue, l’évaporation de la surface en eau, Le débit ments étant donné que le niveau de sédiments n’a atteint le
de déversement (m3s-1), le débit de vidange pour alimenter niveau des prises Jet Creux qu’en 2010, les flux de vidange
l’ONEP (m3s-1), le débit de vidange de fond (m3s-1) et la de fond sont supposés avoir une concentration moyenne
variation du volume stocké (m3s-1), ce dernier pouvant être de 0,5 g l-1 et les flux de déversement une concentration
positif ou négatif. Nous calculons ensuite, au pas de temps moyenne de 0,1 g  l-1. Les campagnes bathymétriques
horaire, le débit entrant Qe= dV/dt+ Qdev + Qvi où dV/dt effectuées sur Ibn Batouta (Tableau 3) font apparaître que
est la variation du volume stocké, Qdev est le déversement la retenue a perdu 33 % de sa capacité entre 1977 et 2013,
et Qvi le débit de vidange. Les débits journaliers puis men- soit une perte annuelle de 1 % (ou une durée de vie du
suels sont ensuite dérivés des valeurs horaires. Une crue barrage estimée à 100 ans).

Tableau 2. Base de données utilisée pour le paramétrage, la calibration et la validation du modèle SWAT
Échelle/
Type Source Description
Période
MNT SPOT 20 m Modèle numérique du terrain
Sol LCGS/Inypsa (1989) 1/50000 Classes du sol/ Caractéristiques pédologique (Figure 1b)
Classification d’une image Land-
Occupation du sol 1/50000 Classes d’occupation du sol (Figure 1a)
sat (2009, 30 m)
ABHL - Agence du Bassin Précipitations journalière (Station Ibn Batouta, Béni
Pluviométrie 1980-2010
Hydraulique de Loukkos Harchane, Romaine, Kalaya)
DMN - Direction de la Météoro-
Température 1980-2010 Températures journalières (min., max., moy.)
logie Nationale
Débits liquides ABHL 1980-2010 Débits liquides journaliers (Station Ibn Batouta)
Débits solides Débits reconstitués 1980-2010 Débits solides journaliers (Station Ibn Batouta)
308 Choukri et al.: Analyse du fonctionnement hydro-sédimentaire d’un bassin versant

Tableau 3. Mesures de bathymétrie disponibles au Les méthodes SCS (USDA-SCS, 1972) et Muskingum
niveau du barrage Ibn Batouta (Overton, 1966; Brakensiek, 1967) ont été retenues dans
cette étude pour estimer respectivement le ruissellement
Date Volume normal (Mm3) de surface et le transfert d’eau et de sédiments dans le ré-
01/05/1983 41 seau hydrographique. Le modèle SWAT utilise la méthode
MUSLE (équation 2) pour simuler l‘érosion au niveau des
01/04/1989 38,5 bassins versants (Neitsch et al., 2005).
01/05/1994 34,3
Sed = 11,8 * (Vtot*Qmax)0,56 * K * LS * C * P (équation 2)
15/05/2000 33,1
30/11/2003 32,5 Où Sed: flux de sédiments dans un jour donné (tonnes
31/10/2007 31,8 métriques), Vtot: volume de ruissellement de surface (en
m3), Qmax: débit de pointe (en m3 s-1), K: le facteur d’éro-
01/09/2013 29,12 dibilité du sol, C: le facteur de couverture et de gestion,
P est le facteur des pratiques antiérosives, LS: le facteur
Mise en œuvre du modèle SWAT topographique. Vtot est estimé, dans la modélisation SWAT
par l’équation empirique SCS, alors que Qmaxest estimé
Le modèle SWAT (Soil and Water Assessment Tool) est un par la méthode rationnelle modifiée (Neitsch et al., 2005).
modèle hydrologique semi distribué à base physique. Dans
son utilisation la plus courante, l’ensemble des processus MUSLE (Williams, 1975) est la version modifiée de
hydrologiques sont simulés sur un pas de temps journalier l’Universal Soil Loss Equation (USLE) développé par
à l’échelle de l’unité de réponse hydrologique (URHs), Wischmeier et Smith (1965, 1978). Elle permet le calcul de
l’unité de base pour la modélisation SWAT. Chaque unité l’érosion causée par les précipitations et le ruissellement,
présente une combinaison unique d’occupation de sol, de alors qu’USLE prédit l’érosion moyenne interannuelle en
type de sol à l›intérieur d›un sous bassin(Gassman et al., fonction de l’énergie de la pluie. Dans MUSLE, le facteur
2007).Le modèle SWAT décompose le fonctionnement d’énergie des précipitations est remplacé par un coefficient
du bassin en deux grandes phases : la phase terrestre où représentant le ruissellement. Il représente l’énergie utili-
il simule les déplacements de l’eau et des sédiments de sée dans l’arrachement et le transport des sédiments dans
la zone de production vers le réseau hydrographique et la les cours d’eau et pas seulement ceux liés au détachement
phase fluviale où il représente le routage de ces éléments comme pour l’équation classique d’USLE (Neitsch et al.,
tout au long des cours d›eau jusqu‘à l‘exutoire du bassin 2005). Les sorties en sédiments sont restituées en tonnes
(Neitsch et al., 2005). La simulation des apports en eaux par hectare à l’échelle de chaque URH, de chaque sous-
et en sédiments au niveau du bassin versant Tleta a été bassins, ou du bassin dans sa totalité.
faite à partir des données présentées dans le tableau 2 et
de la version ArcSWAT 2009. La structuration des don- Analyse de sensibilité, calage et validation du modèle
nées ainsi que la démarche suivie pour l’intégration de SWAT
ces informations dans le modèle a été réalisée en suivant La phase décalage consiste à ajuster certains paramètres
la procédure décrite dans la documentation du modèle. clés pour faire coïncider au mieux les données observées
Après la préparation des entrées de SWAT, nous sommes avec celles simulées par le modèle SWAT. Généralement,
passés à la phase de paramétrage qui englobe trois grandes une période de temps qui dépasse au moins trois ans est
phases principales: nécessaire pour cette étape de calage. Pour mieux satis-
• La prise en compte de la topographie à travers le MNT faire cette étape, une analyse de sensibilité permet préala-
(Modèle Numérique du Terrain)a permis de délimiter le blement de déterminer les paramètres clés qui méritent une
bassin et les sous bassins, de définir le réseau hydrogra- grande attention lors du calage. L’analyse de sensibilité
phique et de dériver la carte des pentes. Quatre classes de et l’estimation des paramètres optimaux du modèle ont
pente ont été considérées (0-5 %; 5-15 %; 15-30 % et > été réalisées en mixant des approches semi-automatiques
30 %) pour le bassin versant Tleta; regroupées dans l’outil libre SWAT-CUP, et une approche
manuelle («Manual calibration») pour finaliser l’ajus-
• La superposition des informations, relatives à la pente, tement des paramètres. Pour évaluer la performance du
l’occupation du sol et au type du sol a permis de discré- modèle à reproduire les observations, nous avons retenus
tiser le bassin d’étude en unités de réponse hydrologique plusieurs indicateurs statistiques recommandés par les
(URHs); utilisateurs de SWAT (Arnold et al., 2012): le coefficient
• La prise en compte des données de forçage clima- de détermination (R²),le coefficient de Nash-Sutcliffe
tique pour la période de référence de 31 ans (1980-2010) : (NSE) et l’indicateur de biais (%PBIAS). Pour donner
précipitations journalières et températures (min, max., un ordre de grandeur, Moriasi et al. (2007) considèrent
moyenne journalière) collectées au niveau des quatre comme satisfaisantes des valeurs de NSE dépassant 0,75
stations météorologiques (Tableau 3).En l’absence des (admissibles à partir de 0,65), des valeurs de NSE dépas-
données relatives à la vitesse du vent, l’humidité relative sant 50 %, et des valeurs de biais inférieurs à 25% pour
moyenne et la radiation solaire, la méthode Hargreaves les débits et inférieurs à 55% pour les sédiments. La pro-
(Hargreaves et Samani, 1982) basée uniquement sur les cédure traditionnelle du calage avec SWAT est d’ajuster
températures maximales et minimales journalières a été successivement les débits, les sédiments et finalement les
adoptée pour la simulation de l’ETP. nutriments.
 Rev. Mar. Sci. Agron. Vét. (2019) 7 (2): 304-312 309

Pour cette étude, la période de référence de 31 ans (1980- Résultats et discussion

2010) a été répartie en trois périodes: 1980-1982 pour
l’initialisation du modèle et pour prendre en compte des Analyse de sensibilité
conditions initiales du milieu, 1983-1996 pour le calage
et 1997-2010 pour la validation. L’analyse de sensibilité réalisée à l’aide de SWAT-CUP a
permis d’identifier un certain nombre de paramètres in-
Le calage des paramètres hydrologiques consiste à ajuster fluents. Pour la partie hydrologique, ces paramètres sont liés
au mieux les différents paramètres afin de reproduire les à la partition entre ruissellement et infiltration (CN2, curve
différents termes du bilan en eau (Rasolomanana, 2012). number), ou l’interaction entre le compartiment sous-ter-
Cette première étape du calage de SWAT est cruciale, car rain et l’écoulement dans l’oued (RCHRG_DP, GWQMN,
l‘eau est la composante la plus importante et la plus cri- REVAPMN, GW_DELAY), et un facteur de compensation
tique du fait qu’elle présente la force motrice qui contrôle pour le calcul de l’évapotranspiration (ESCO). Concernant
le fonctionnement des phénomènes qui se déroule au la dynamique d’exportation des sédiments, les paramètres
niveau de l’hydrosystème (érosion, envasement, transfert qui se sont avérés les plus sensibles sont les caractéris-
de polluants, etc.). Suivant une démarche classique (Car- tiques de longueur de pente (SLSUBBSN) et d’inclinaison
valho-Santos et al., 2015), le calage a consisté à reproduire (SLOPE), les coefficients de rugosité de Manning (OV_N)
au mieux les débits journaliers à l’exutoire du bassin et le facteur d’ajustement des débits de pointe pour le trans-
versant. Pour cette étude, la série de débits entrants dans port fluvial des sédiments (PRF).
le barrage Ibn Batouta pour la période 1983-1996 a été
considérée. Pour affiner le calage des processus hydrolo- Calibration et validation du modèle
giques, nous avons également cherché à bien reproduire
la contribution du ruissellement direct et la contribution La figure 2 montre la comparaison graphique des débits
du ruissellement retardé. Deux méthodes ont été utilisées mensuels observés et simulés par SWAT après calage.
pour séparer ces 2 contributions dans les données mesu- Les graphes de calibration et de validation présentent
rées : une méthode graphique sur quelques crues et une une bonne similarité entre les données de ruissellement
méthode de filtrage automatique sur l’ensemble de la série. simulées et observées pour la zone d’étude. Les critères
La séparation du ruissellement entre ces deux compo- statistiques d’estimation de la qualité du modèle confir-
santes par la méthode de filtrage automatique ne permet ment cette impression visuelle avec les valeurs suivantes:
pas de fournir un résultat fiable à l’échelle de chaque crue R² (0,95), NSE (0,94) et %PBIAIS (2) pour la période de
mais elle permet de fournir une indication assez fiable de calibration (1983-1996) contre R² (0,90), NSE (0,89) et
la contribution moyenne interannuelle de chaque compo- %PBIAS (-6) pour la période de validation (1997-2010).
sante. Si la fonction objective visée dans la phase de calage Les résultats obtenus montrent donc que SWAT est capable
reste une bonne reproduction des débits,la répartition des de reproduire correctement les données observées à l’exu-
différentes composantes du ruissellement sera également toire du bassin étudié au pas de temps mensuel.
analysée et discutée. Une fois les flux d’eau correctement Les performances mensuelles obtenues pour le calage des
reproduits, le calage des paramètres d’érosion les plus flux de sédiments transportés vers l’exutoire du bassin ver-
influents, a été réalisé en tentant de minimiser les écarts sant Tleta sont, selon les critères d’évaluation de Moriasi et
entre les flux de sédiments arrivant au barrage issus de la al. (2007), satisfaisantes pour la période de calibration (R²
mesure et du modèle sur la période 1983-1996. de 0.83, NSE de 0.82 et %PBIAIS de +17%) et acceptables

40
(a) nse = 0,94
Observé Simulé PBias = -6
30 R2 = 0,95
Débits (m3/s)

20

10

0
mai-83 mai-84 mai-85 mai-86 mai-87 mai-88 mai-89 mai-90 mai-91 mai-92 mai-93 mai-94 mai-95 mai-96

30
(b) Observé Simulé nse =0,90
PBias= 2
R2 = 0.89
Débits (m3/s)

20

10

0
janv.-97 janv.-98 janv.-99 janv.-00 janv.-01 janv.-02 janv.-03 janv.-04 janv.-05 janv.-06 janv.-07 janv.-08 janv.-09 janv.-10

Figure 2. Comparaison entre les débits mensuels observés et simulés par SWAT: 
(a) période de calibration et (b) période de validation.
310 Choukri et al.: Analyse du fonctionnement hydro-sédimentaire d’un bassin versant

pour la période de validation (R² de 0.64, NSE de 0.61 rect et de l’écoulement retardé étaient bien respectées dans
et %PBIAIS de -6%). Ces chiffres traduisent une légère les simulations. L’application de la méthode graphique
sous-estimation des sédiments observés pour la période de séparation des hydrogrammes a permis de quantifier
de calibration et une légère surestimation des sédiments à 20 % ou 25 % environ la contribution de l’écoulement
observés pour la période de validation (Figure 3). Le taux souterrain (incluant le débit de base et l’écoulement hypo-
d’érosion moyen simulé par SWAT au niveau de Tleta est dermique)au niveau du bassin versant Tleta, alors que le
de 29 t/ha/an pour la période de validation (1983-2010) méthode par filtrage automatique l’a évalué entre 15 à
alors que la dégradation spécifique sur ce bassin issue 25 % selon les paramètres de filtrage retenus. L’analyse du
des mesures bathymétriques corrigées des sédiments non bilan nous permet donc d’identifier que le modèle semble
piégés est de 28,8 t ha-1 an-1 sur la période (1983-2013, pas sous-estimer légèrement la part du ruissellement direct et
de bathymétrie disponible en 2010). Cette valeur globale surestimer légèrement la part du ruissellement retardé.
pour le bassin ne doit cependant pas masquer de fortes Les pistes de correction de cette légère anomalie reposent
disparités au sein du bassin versant. En effet de nombreux sur une meilleure estimation des paramètres régissant la
secteurs du bassin constituent aujourd’hui des badlands, partition entre ruissellement et infiltration ou l’interaction
zones sources de forts départs de sédiments. Les zones entre le compartiment sous-terrain et l’écoulement dans
cultivées ou en jachère sur fortes pentes peuvent aussi l’oued. D’autres sources d’incertitudes que le calage du
constituer des sources importantes de sédiments pour le modèle peuvent également expliquer ce léger décalage
barrage. Les sédiments sortants du barrage par vidange ou entre valeurs simulées et valeurs mesurées. Par exemple,
par déversement représentent moins de 5 % des sédiments la méthode de Hargreaves utilisée pour calculer l’éva-
entrant dans la retenue. Le taux d’arrachement du sol pour potranspiration potentielle est peu précise, puisqu’elle
le nord du Rif se situe généralement entre 1 000 et 3 000 implique seulement les paramètres min et max de la
m3km-2an-1, correspondant à des dégradations spécifiques température et que les valeurs utilisées pour ces deux
de 12 à 35 t ha-1an-1. paramètres ont été mesurées sur une station climatique
située à l’extérieur du bassin Tleta (aérodrome de Tanger),
Analyse du bilan hydrique faute de données disponibles au sein du bassin versant.
Les sorties du modèle SWAT entre 1983 et 2010 intégrant Tableau 4. Analyse des composantes hydrologiques
la période de calage et de validation permettent d’appré- issues du modèle SWAT calé et validé
hender les différents termes du bilan hydrique (Tableau
4). Sur cette période, les résultats de simulation montrent Quantités
Fractionnement du bilan hydrique
que sur une moyenne annuelle de 777 mm de précipitation (mm)
sur l’ensemble du bassin d’Ibn Batouta, environ 60  % Précipitation totale 777
(458,4  mm) retourne à l’atmosphère par évapotranspi- Ruissellement 160,3
ration. Le débit du cours d’eau tel que simulé par SWAT
totalise 283,4  mm qui se divise en ruissellement direct Mouvements latéraux d’eau dans le sol
9,4
(160,3 mm), flux latéraux dans le sol (9,4 mm) et contri- contribuant au débit du cours d’eau
bution des aquifères peu profonds (113,7 mm).Ce total est Contribution de l’eau souterraine (aquifère
113,7
légèrement supérieur à la valeur moyenne des flux d’eau peu profond) au débit du cours d’eau
arrivant au barrage annuellement (282 mm). Pour vérifier Évapotranspiration réelle 458,4
que les différents termes du bilan sont bien reproduits,
nous avons analysé les données mesurées entre 1983 et Évapotranspiration potentielle 1266,3
2010. Pour vérifier si les contributions de l’écoulement di-
12
(a) Observé Simulé nse = 0,82
10 PBias = 17
R2 = 0.83
Sédiments (105 t)

0
mai-83 mai-84 mai-85 mai-86 mai-87 mai-88 mai-89 mai-90 mai-91 mai-92 mai-93 mai-94 mai-95 mai-96

9
(b) nse = 0,61
8 Observé Simulé
PBias = -6
Sédiments (105t)

7
R2 = 0.64
6
5
4
3
2
1
0
janv.-97 janv.-98 janv.-99 janv.-00 janv.-01 janv.-02 janv.-03 janv.-04 janv.-05 janv.-06 janv.-07 janv.-08 janv.-09 janv.-10

Figure 3. Comparaison entre les flux mensuels de sédiments observés et simulés par SWAT :(a) période de calibration et (b)
période de validation
 Rev. Mar. Sci. Agron. Vét. (2019) 7 (2): 304-312 311

Conclusion Brakensiek, D. L. (1967). Kinematic Flood Routing. Tran-

sactions of the American Society of Agricultural Engineers
Le modèle SWAT est un outil très utile pour la gestion 10(3): 340-343.
orientée des bassins versants hydrologiques. La calibra-
tion et la validation du modèle sont des étapes clés dans la Carvalho-Santos, C., Nunes, J. P., Monteiro, A. T., Hein,
simulation et la réduction des incertitudes. L’étude réalisée L. and Honrado, J. P. (2015). Assessing the effects of land
cover and future climate conditions on the provision of
au niveau du bassin versant Tleta a montré que le modèle hydrological services in a medium-sized watershed of
SWAT peut simuler d’une façon satisfaisante les débits et Portugal.Hydrological processes.
les sédiments dans une zone sous influence atlantique et
méditerranéenne semi-aride. L’évaluation de la qualité de Chaponnière, A. (2005). Fonctionnement hydrologique
simulation a été réalisée par comparaison avec les mesures d’un bassin versant montagneux semi-aride: cas du bassin
de terrain suivant trois estimateurs statistiques. Les valeurs versant du Rehraya (Haut Atlas marocain).
trouvées de R² (0,95; 0,90), NSE (0,94; 0,89) et %PBIAS Choukri, F., Chikhaoui, M., Naimi, M., Raclot, D., Pépin,
(-6; 2) pour les débits respectivement pour la calibration et Y. and Lafia, K. (2016). Impact du changement climatique
la validation et les valeurs de R² (0,83; 0,64), NSE (0,82; sur l’évolution de l’érosivitédes pluies dans le Rif Occi-
0,61) et %PBIAS (17; -6) pour le transport solide confirme dental (Nord du Maroc). European Scientific Journal, ESJ
cette bonne qualité de simulation. Quelques pistes d’amé- 12(32).
lioration ont été évoquées, comme par exemple, renfor- Fadil, A., Rhinane, H., Kaoukaya, A., Kharchaf, Y. and
cer la connaissance des propriétés hydrodynamiques et Bachir, O. A. (2011).Hydrologic modeling of the Boure-
érosives des sols, renforcer les mesures météorologiques greg Watershed (Morocco) using GIS and SWAT model.
au sein même du bassin afin d’améliorer l’estimation de Journal of Geographic Information System 3(04): 279.
l’évapotranspiration potentielle (ETP), ou encore affiner Ficklin, D. L. (2010). Modeling the impacts of climate
le calage des paramètres régissant la partition infiltration/ change on hydrology and agricultural pollutant runoff in
ruissellement et les échanges nappe-rivière. Cependant, California’s Central Valley.University of California, Davis.
le modèle fournit des résultats largement satisfaisants et Ficklin, D. L., Luo, Y., Luedeling, E. and Zhang, M. (2009).
peut donc être utilisé dans le bassin versant de Tleta pour Climate change sensitivity assessment of a highly agricul-
évaluer l’impact combiné des scénarios d’occupation du tural watershed using SWAT.
sol et climatiques sur la quantité et la qualité des eaux
mobilisées par le barrage d’Ibn Batouta. Gassman, P. W., Reyes, M. R., Green, C. H. and Arnold, J.
G. (2007). The soil and water assessment tool: historical
development, applications, and future research directions.
Remerciements Center for Agricultural and Rural Development, Iowa State
Les auteurs remercient le programme JEAI Vecteur financé par University.
l’IRD et le projet ALMIRA (ANR-12-TMED-0003) financé par
l’ANR pour leur support financier et technique. Gatzke, S. E. (2010). Modeling the effects of agricultural
conservation practices and climate variability on sediment
yield and transport in the Colusa Basin, California.Univer-
Références sity of California, Davis.
Albergel J., Mansouri T., Pépin Y, Seguis L. (1999). Hargreaves, G. H. (1975). Moisture Availability and Crop
Reconstitution des crues sur les bassins versants des lacs Production, Trans. Am. Soc. Agric. Eng., 18(5): 980-984.
collinaires en Tunisie Centrale et au Cap-Bon. (Atelier sur
l’estimation des crues et des étiages. ENIT - INRS-eau Uni- Hargreaves, G. H. and Samani, Z. A. (1982). Estimating
versité du Québec - AUPELF réf. Tunis, 16-17 février 99. potential evapotranspiration. Journal of the Irrigation and
Drainage Division 108(3): 225-230.
Arnold, J., Moriasi, D., Gassman, P., Abbaspour, K., White,
M., Srinivasan, R., Santhi, C., Harmel, R., Van Griensven, Hunink, J., Niadas, I., Antonaropoulos, P., Droogers, P.
A. and Van Liew, M. (2012). SWAT: Model use, calibra- and De Vente, J. (2013). Targeting of intervention areas
tion, and validation. Transactions of the ASABE 55(4): to reduce reservoir sedimentation in the Tana catchment
1491-1508. (Kenya) using SWAT. Hydrological Sciences Journal
58(3): 600-614.
Arnold, J. G. and Fohrer, N. (2005). SWAT2000: current
capabilities and research opportunities in applied water- Kim, N. W., Lee, J. W., Lee, J. and Lee, J. E. (2010). SWAT
shed modelling. Hydrological processes 19(3): 563-572. application to estimate design runoff curve number for
South Korean conditions. Hydrological processes 24(15):
Bakker, M. M., Govers, G., Kosmas, C., Vanacker, V., 2156-2170.
Van Oost, K. and Rounsevell, M. (2005).Soil erosion as
a driver of land-use change. Agriculture, ecosystems & Knisel, W. G. (1980). CREAMS: a field scale model for
environment 105(3): 467-481. Chemicals, Runoff, and Erosion from Agricultural Mana-
gement Systems [USA]. United States. Dept. of Agricul-
Bakker, M. M., Govers, G. and Rounsevell, M. D. (2004). ture. Conservation research report (USA).
The crop productivity-erosion relationship: an analysis
based on experimental work. Catena 57(1): 55-76. Laurent, F., Ruelland, D. et Chapdelaine, M. (2007). Simu-
lation de l’effet de changements de pratiques agricoles
Bonumá, N. B., Rossi, C. G., Arnold, J. G., Reichert, J. M., sur la qualité des eaux avec le modèle SWAT. Revue des
Minella, J. P., Allen, P. M. and Volk, M. (2014). Simulating sciences de l’eau/Journal of Water Science 20(4): 395-408.
landscape sediment transport capacity by using a modified
SWAT model. Journal of Environmental Quality 43(1): Leonard, R., Knisel, W. and Still, D. (1987). GLEAMS:
55-66. Groundwater loading effects of agricultural management
systems. Transactions of the ASAE 30(5): 1403-1418.
312 Choukri et al.: Analyse du fonctionnement hydro-sédimentaire d’un bassin versant

Lévesque, É. (2007). Évaluation de la performance hydro- Yang, J., Reichert, P., Abbaspour, K., Xia, J. and Yang, H.
logique du modèle SWAT pour de petits bassins versants (2008).Comparing uncertainty analysis techniques for a
agricoles du Québec. Université Laval. SWAT application to the Chaohe Basin in China. Journal
of Hydrology 358(1): 1-23.
Merzouk, A., Fenjiro, I. et Laouina, A. (1996). Carto-
graphie de l’évolution des formes d’érosion dans le Rif Yang, Q., Meng, F.-R., Zhao, Z., Chow, T. L., Benoy, G.,
occidental (Maroc): étude multidate utilisant un SIG bassin Rees, H. W. and Bourque, C. P.-A. (2009). Assessing the
versant. Bulletin-Réseau Erosion (16): 444-456. impacts of flow diversion terraces on stream water and
sediment yields at a watershed level using SWAT model.
Michard, A. (1976). Eléments de géologie marocaine. Agriculture, ecosystems & environment 132(1): 23-31.
Éditions du Service géologique du Maroc.
Zhang, X., Srinivasan, R. and Van Liew, M. (2008).
Michaud, A., Deslandes, J. et Beaudin, I. (2006).Modé- Multi-site calibration of the SWAT model for hydrologic
lisation de l’hydrologie et des dynamiques de pollution modeling.Transactions of the ASABE 51(6): 2039-2049.
diffuse dans le bassin versant de la rivière aux Brochets à
l’aide du modèle SWAT–Rapport final. Sainte-Foy, Institut
de recherche et développement en agroenvironnement,[En
ligne]. http://www. irda. qc. ca/_documents/_Results/18.
pdf (5 octobre 2006).
Moriasi, D., Arnold, J., Van Liew, M., Bingner, R., Harmel,
R. and Veith, T. (2007). Model evaluation guidelines for
systematic quantification of accuracy in watershed simu-
lations. Trans. ASABE 50(3): 885-900.
Mosbahi, M., Benabdallah, S. and Boussema, M. (2011).
Hydrological Modeling in a Semi-Arid Catchment Using
SWAT Model. Journal of Environmental Science and
Engineering 5(12): 1695-1701.
Neitsch, S., Arnold, J., Kiniry, J., Williams, J. &King, K.
(2005). Soil and water assessment tool: theoretical docu-
mentation, version 2005. Texas, USA.
Ouessar, M., Bruggeman, A., Abdelli, F., Mohtar, R.,
Gabriels, D. and Cornelis, W. (2008).Modelling water-har-
vesting systems in the arid south of Tunisia using SWAT.
Hydrology and Earth System Sciences Discussions 5(4):
1863-1902.
Overton, D. (1966). Muskingum flood routing of upland
streamflow. Journal of Hydrology 4: 185-200.
Partie, I. (2006). Variabilité spatio-temporelle des exporta-
tions de sédiments et de phosphore dans le bassin versant
de la Rivière aux Brochets au sud-ouest du Québec.
Rasolomanana, S. D. (2012). Modélisation des petits bas-
sins versants en milieu rural: problématique des modèles
fortement paramétrés. Université Laval.
Salou, B. I. G. (2009). Estimation des apports au lac Saint-
Charles à l’aide du modèle hydrologique SWAT. Université
Laval.
Service, U. S. S. C. (1972). SCS national engineering
handbook, section 4: hydrology. The Service.
Williams, J., Arnold, J., Kiniry, J., Gassman, P. and Green,
C. (2008).History of model development at Temple, Texas.
Hydrological Sciences Journal 53(5): 948-960.
Williams, J. R. (1975). Sediment-yield prediction with
universal equation using runoff energy factor.
Williams, J. R., Nicks, A. and Arnold, J. G. (1985).Simula-
tor for water resources in rural basins. Journal of Hydraulic
Engineering 111(6): 970-986.
Wischmeier, W. H. and Smith, D. D. (1978).Predicting
rainfall erosion losses - A guide to conservation planning.
Predicting rainfall erosion losses - A guide to conservation
planning.