Corps de l’article

1. Introduction

La gestion des ressources hydriques est une des priorités pour assurer la qualité de vie dans les sociétés actuelles. Au Maroc, cette gestion prend une importance particulière dans le développement du pays étant donné les distributions spatiale et temporelle inégales des précipitations (AGOUSSINE et BOUCHAOU, 2004; ALLOUZI et BOUAAM, 1994). Avec les périodes de sécheresse qui sont devenues de plus en plus fréquentes au cours des dernières décennies, les réservoirs collinaires ont pris une importance accrue dans la politique de gestion des ressources en eau (MAHFOUD, 1989). Au Maroc, plus de 104 lacs et 54 barrages collinaires ont été construits depuis 1980 au cours de différents projets de mobilisation des eaux superficielles. Ils ont été l’oeuvre de deux ministères : le ministère des Travaux publics chargé de l’équipement hydraulique du pays et le ministère de l’Intérieur chargé du programme de lutte contre les effets de la sécheresse.

Si ces ouvrages répondent parfaitement à leurs objectifs nationaux de mobilisation de l’eau, ils souffrent cependant, dans le nord du Maroc, d’un envasement rapide de leurs retenues dû à une forte érosion des sols sur leurs bassins versants, laquelle atteint des taux très élevés, compris entre 24 t•ha‑1•an‑1 sur parcelles et 140 t•ha‑1•an‑1 sur ravines dans les montagnes du Rif et du Pré-Rif marneux (HEUSH, 1970; MERZOUK et DAHMAN, 1998; TAYAA, 1997). En Afrique du Nord, la perte annuelle en volumes d’eau stockés par envasement des petites retenues est estimée comprise entre 4 et 5 % (ALBERGEL et al., 2003b). Elle est du même ordre de grandeur en Italie pour les petites retenues (GAZZALO et BASSI, 1969). Pour les grands barrages du Maroc, la perte annuelle en volumes d’eau stockés est estimée à 0,5 % (MARZOUKI, 1992).

Au Maroc, la recherche hydrologique à l’échelle des petits bassins versants n’a pas été suffisamment développée pour produire des modèles d’estimation fiables et précis des apports en eau et en sédiments dans les petites retenues (ALBERGEL et al., 2004). Des exemples connus de surdimensionnement et de sous-dimensionnement de la capacité de stockage des ouvrages font que l’estimation des apports hydro-sédimentaires dans les petites retenues est un enjeu économique important. La référence nationale concernant le calcul des débits de projet et le dimensionnement de ces petites retenues souffrent d’importantes lacunes. Les modèles empiriques, Curve Number (CN), modèles rationnels, hydrogramme unitaire, équations universelles de perte en sol MUSLE et RUSLE (WILLIAMS, 1975) donnent des résultats moyens sur des périodes de longue durée mais ne sont pas fiables pour évaluer l’impact d’une évolution du milieu naturel ou des pratiques culturales sur les régimes hydrologiques. Les modèles distribués à base physique sont de plus en plus nombreux et se veulent être de bons outils de simulation et de prévision (AMBROISE, 1999). Ils sont mis au point et testés sur des bassins versants connus pour lesquels on dispose d’importants jeux de données. Ils mettent en jeu de nombreuses variables à l’échelle de la surface concernée, variables qui ne sont pas faciles à déterminer. La calibration, à partir de jeux de données d’entrée et de sortie, est la technique la plus répandue pour s’affranchir des mesures que l’on ne sait pas faire. La validation, effectuée en utilisant d’autres jeux de données d’entrée et de sortie que ceux qui ont été utilisés pour le calage, permet d’évaluer l’adéquation du modèle à simuler le fonctionnement hydrologique du bassin versant.

Des travaux similaires ont été réalisés sur les petits bassins versants au sud de la Méditerranée, en Tunisie et en Algérie (ALBERGEL et al., 2003a; GOMER, 1994). Parmi ces modèles, on trouve le modèle « Mina », spécifiquement adapté aux petits bassins versants marneux, mais il nécessite une validation pour les autres régions d’Afrique du Nord que celle où il a été validé. En effet, le modèle « Mina » a été mis au point, calibré et validé, sur plusieurs petits bassins versants de 40 à 70 ha du grand bassin versant de l’Oued Mina au nord-ouest de l’Algérie (GOMER et al., 1992; GOMER, 1994). L’objectif principal de ce travail est donc de transposer ce modèle au bassin versant du barrage collinaire de l’Oued Saboun situé au nord-ouest du Maroc et de proposer, pour cette région du Rif, un outil de prévision des apports liquides et solides dans les petites retenues artificielles.

L’étude est basée sur un découpage du bassin versant en unités hydrologiques homogènes dont les valeurs des paramètres de l’infiltration sont déduites d’expérimentations in situ, à l’échelle du mètre carré sous simulations de pluies. La validation du modèle est présentée par confrontation des données simulées aux données observées à deux pas de temps : celui de la crue pour les apports liquides et celui de l’année pour les apports solides. Pour la spatialisation et la simulation hydrologique, ce travail utilise les outils des nouvelles technologies de géomatique (SIG) et d’acquisition des données hydrométéorologiques.

2. Matériels et méthodes

2.1 Le bassin versant de l’Oued Saboun

Le présent travail a été conduit sur le petit bassin versant marneux de l’Oued Saboun situé au nord-ouest du Maroc dans le grand bassin versant de l’Oued Tahadart qui draine les eaux de son bassin vers l’océan Atlantique (Figure 1). Un petit barrage en terre, de 15 m de hauteur, a été construit à l’exutoire du bassin versant de l’Oued Saboun. Mise en eau en 1991, sa retenue, d’une capacité de 1 100 000 m3, stocke les eaux de ruissellement afin d’irriguer un petit périmètre en aval de la zone de Fahs, près de la ville de Tanger (Tableau 1).

Figure 1

Situation du site d’étude.

Location of the study site.

Situation du site d’étude.

-> Voir la liste des figures

Tableau 1

Caractéristiques des retenues de Saboun et de Daimos et de leurs bassins versants.

Characteristics of the Saboun and Daimos dams and their catchment areas.

 

Saboun

Daimos

Type

Terre homogène

Terre homogène

Hauteur

15 m

5 m

Longueur de la digue

443 m

90 m

Apport moyen annuel

1 100 000 m3

250 000 m3

Capacité du réservoir

1 100 000 m3

30 000 m3

Surface du bassin versant

802 ha

228 ha

Altitude min du bassin versant

17 m

75 m

Altitude max du bassin versant

130 m

130 m

-> Voir la liste des tableaux

En amont de ce bassin versant, il existe également un lac collinaire créé en 1991, en même temps que le barrage Saboun, par la construction d’une digue en terre homogène. Ce lac présente une superficie de 1,3 ha à sa cote de déversement. Il est situé à l’exutoire d’un sous-bassin versant d’une superficie de 228 ha appelé Daimos, nom emprunté au village qui occupe la rive droite du sous-bassin. Il joue aussi un rôle de protection contre l’envasement du barrage Saboun. Son déversoir latéral est constitué de deux conduites en ciment de 80 et 40 cm de diamètre. En 1991, ces conduites étaient situées à environ 3,40 m au-dessus du fond de la retenue. En 2005, elles se situaient à 1,25 m du fond de la retenue, ce qui correspond à un envasement de 2,15 m de cette retenue sur la période 1991-2005.

Du point de vue climatique, la région de Tanger est caractérisée par un climat de type méditerranéen sub-humide (BEN ABID, 1982). Elle est considérée comme l’une des régions les plus pluvieuses du Maroc. Calculée sur une période de 30 ans, de 1969 à 1999 à la station de Kalaya, la plus proche du site de Saboun, la moyenne des précipitations annuelles est de 742 mm. La pluviométrie annuelle décennale sèche est de 490 mm, la pluviométrie annuelle décennale humide de 995 mm. Au cours des années sèches, les ressources en eau superficielles de cette région diminuent fortement et ne sont pas compensées par des ressources en eaux souterraines suffisantes, le déficit en eaux superficielles pouvant atteindre un taux très élevé, de l’ordre 84 % pour l’année 1998-99, année vicennale sèche (20 ans) et de 90 % pour l’année 1994-95 (ONEP, 1996). Le déficit en eaux souterraines de la région de Tanger s’explique par le fait qu’une grande partie du substratum de la région (Unité de Tanger) est formée d’argiles et de marnes d’âge Sénonien (TAUVIN, 1971). Les faciès argilo-marneux de l’Unité de Tanger empêchent toute percolation en profondeur, limitant ainsi la formation des nappes d’eaux souterraines. Les seules ressources en eau susceptibles d’être mobilisées en quantité importante sont donc les eaux superficielles dont les écoulements sont essentiellement générés par les fortes pluies.

2.2 Le modèle « Mina »

Le modèle « Mina » est un modèle déterministe semi-physique de type distribué, basé sur les concepts d’infiltration de Horton pour le calcul du ruissellement et d’énergie cinétique des pluies et de l’écoulement pour la production de l’érosion. Les mécanismes de transfert ou de routage sont assurés par l’équation de continuité mono-dimensionnelle de Saint-Venant (GOMER, 1994). C’est un modèle qui a été élaboré dans le but de permettre les quantifications de l’écoulement et de l’érosion des petits bassins versants marneux. Il a été choisi pour simuler les crues et l’érosion du bassin versant de Saboun parce que ses caractéristiques répondent parfaitement aux conditions de la zone d’étude où le processus d’infiltration dominant est a priori de type Hortonien (ALBERGEL et al., 2003a). Il présente en outre les avantages suivants :

  • c’est un modèle de type déterministe à base physique;

  • c’est un modèle distribué à mailles régulières;

  • c’est un modèle événementiel simulant l’écoulement crue par crue;

  • c’est un modèle adapté aux zones semi-arides et aux sols marneux; et

  • pour un modèle à base physique, c’est un modèle dont le nombre de paramètres est réduit.

Les principales données d’entrée du modèle « Mina » sont la pluie (les cumuls pluviométriques journaliers et les intensités horaires des pluies), la carte des unités hydrologiques, la carte des pentes et celle du réseau hydrographique, le taux d’infiltration et l’humidité du sol avant l’averse et l’énergie cinétique critique d’arrachement des particules du sol.

2.3 Constitution d’un échantillon de données hydrologiques

2.3.1 Équipement hydrométéorologique utilisé

Le barrage Saboun a été équipé d’une station hydrométéorologique automatique (MERZOUK et al., 1997) installée en novembre 1997 sur la digue de la retenue. Elle se compose :

  • d’une sonde piézo-résistive immergée permettant d’enregistrer automatiquement les variations de hauteur d’eau du lac à 1 cm près sur un pas de temps de 5 min;

  • d’un pluviographe à augets basculeurs permettant d’enregistrer le cumul des précipitations avec une précision de 0,5 mm;

  • d’un bac d’évaporation enterré, de type Colorado-ORSTOM, de 1 m2 d’ouverture et 60 cm de profondeur, la mesure étant effectuée manuellement chaque matin par l’observateur;

  • d’une batterie d’échelles limnimétriques implantées sur la digue;

  • d’une centrale d’enregistrement et de mémorisation des données; et

  • d’un déversoir calibré permettant de connaître les débits sortant de la retenue sur un pas de temps de 5 min (déversoir de forme rectangulaire de 8 m de largeur). Le déversement a lieu quand l’eau stockée dans la retenue atteint la cote 6,68 m (16,68 m NGM).

Les données limnimétriques et pluviométriques enregistrées sont transférées régulièrement au laboratoire pour établir une base de données hydrométéorologiques.

Le lac Daimos n’a pas été équipé d’une station hydrométrique. Cependant, nous avons considéré que la pluie tombée au niveau de la retenue du barrage Saboun est représentative de la pluie tombée sur l’ensemble du bassin. Par comparaison, le bassin versant de Kalaya (38 km2), situé environ 10 km à l’est, dispose de deux pluviomètres implantés à ses deux extrémités : l’un en aval à une altitude de 45 m (station hydrologique de Kalaya, pluviomètre nommé PKayala) et l’autre en amont à 300 m d’altitude (station de Romane, pluviomètre nommé PRomane). Sur la période 1974-2005, les totaux pluviométriques annuels de ces deux pluviomètres sont bien corrélés, avec un coefficient de corrélation de 0,96 : PKalaya = 0,96 x PRomane avec R2 = 0,95. À l’échelle journalière, le coefficient de corrélation entre pluies journalières aux deux stations est encore égal à 0,90 avec R2 = 0,78.

2.3.2 Reconstitution d’un échantillon de crues

La détermination du bilan hydrologique de la retenue du barrage Saboun est basée sur l’application de l’équation du bilan qui découle du principe de conservation des volumes d’eau exprimés en m3 (ALBERGEL et al., 1997), soit :

avec : ∆V la variation du volume, Vr le volume ruisselé, Vecs le volume des apports souterrains, Vp le volume des apports pluviométriques sur la retenue, Vev le volume évaporé, Vd le volume déversé, Vvi le volume vidangé, Vi le volume infiltré et Vu le volume utilisé.

Toutes les variables accessibles à l’observation (Vp, Vd, Vvi, Vev) ont été mesurées pour toute la période comprise entre novembre 1997 et octobre 1999. Le volume ruisselé Vr est la somme des débits entrant dans la retenue, il est obtenu par intégration de l’équation suivante : . Dans la dernière équation, les paramètres sont : Qr débit ruisselé, Qd le débit de déversé et Qp le débit dû à la pluie. Les paramètres non mesurés, le volume infiltré (Vi), le volume des eaux souterraines (Vecs) et le volume utilisé (Vu), sont globalement appréciés par bouclage du bilan.

Sur un petit pas de temps et durant les crues, il est possible de calculer des débits en dérivant l’équation (1) sur le pas de temps de mesure des volumes d’eau (5 min) et en faisant quelques approximations justifiées. Lors d’une crue de courte durée, les apports souterrains du bassin versant, l’évaporation du plan d’eau, l’infiltration et l’utilisation des eaux de la retenue peuvent être négligés. Ceci permet de reconstituer le ruissellement (ALBERGEL et al., 1997) en simplifiant les calculs du bilan à l’échelle de la crue.

Toutes les crues du bassin versant de Saboun ont été ainsi reconstituées entre novembre 1997 et août 1999 (ABDELLAOUI et al., 2002). Elles ont été représentées sur la figure 2 par leurs débits maximaux et leurs pluies journalières. Chaque crue est caractérisée par une chronique de débits sur un pas de temps de 5 min, un volume total et un débit maximal. La figure 3 présente la crue du 4/12/97, reconstituée à l’exutoire du bassin versant de l’Oued Saboun à partir des enregistrements de hauteurs d’eau du lac. Sur cette même figure, les débits déversés forment une crue très aplatie, fortement laminée par la retenue.

Figure 2

Crues et averses reconstituées à Saboun entre novembre 1997 et août 1999.

Stream flows and rainfall reconstituted at the Saboun dam between November 1997 and August 1999.

Crues et averses reconstituées à Saboun entre novembre 1997 et août 1999.

-> Voir la liste des figures

Figure 3

Crue du 04/12/1997, reconstituée au barrage de Saboun.

Flood of 04/12/1997, reconstituted at the Saboun dam.

Crue du 04/12/1997, reconstituée au barrage de Saboun.

-> Voir la liste des figures

2.4 Mesures de la sédimentation et calcul de l’envasement

La détermination du taux d’envasement de la retenue Saboun est basée sur des mesures ponctuelles de la bathymétrie du lac, complétées par une topographie fine des berges exondées (Figure 4). Le volume total des sédiments piégés entre deux campagnes de mesures bathymétriques est estimé par différence des volumes de la retenue et confrontation des modèles numériques de terrain (MNT) créés à partir des levés bathymétriques et topographiques. La différence des volumes d’eau stockés à la cote normale entre l’année 1991 et la dernière mesure bathymétrique représente la part des sédiments érodés sur le bassin versant et transportés par les eaux de ruissellement jusqu’au barrage de l’Oued Saboun depuis la construction du barrage en 1991 (ALBERGEL et al., 2003b).

Figure 4

Évolution de la bathymétrie entre 1999 et 2004 dans la retenue du barrage Saboun.

Evolution of bathymetry between 1999 and 2004 in the Saboun reservoir.

Évolution de la bathymétrie entre 1999 et 2004 dans la retenue du barrage Saboun.

-> Voir la liste des figures

Les résultats de ce contrôle bathymétrique, en matière de volume total envasé et de taux moyen d’envasement, sont présentés dans le tableau 2. Ils prennent en compte les mesures bathymétriques postérieures à la période considérée et permettent d’observer une certaine constance dans l’envasement de la retenue.

Tableau 2

Données disponibles de bathymétrie à Saboun.

Available bathymetric data for the Saboun reservoir.

 

Volume Initial

(m3)

Volume Restant

(m3)

Volume de vase

(m3)

Envasement

(Densité de Vase = 1,3) (t•an‑1)

Perte Annuelle

(%)

Période 1991‑1999

1 100 000

903 990

196 010

31 852

2,0

Période 1991‑2003

1 100 000

817 240

282 760

28 276

2,0

Période 1991‑2004

1 100 000

801 020

298 986

27 762

1,9

-> Voir la liste des tableaux

Le volume total des apports en sédiments arrivant au barrage est calculé en tenant compte du volume d’envasement et du volume des sédiments évacués au moment des déversements et des vidanges. Cette partie des sédiments évacués a été estimée en fonction des volumes totaux déversés et vidangés et de leurs concentrations moyennes en sédiments, concentrations déterminées par analyse de quelques échantillons d’eau prélevés au niveau du déversoir et près de la vidange de fond. Les concentrations en sédiments des déversements varient entre 0,01 et 0,12 g•L‑1, avec une moyenne de 0,06 g•L‑1. Cette faible concentration en sédiments est expliquée par l’orientation de l’axe d’écoulement principal des eaux entrant dans le barrage, l’écoulement principal étant parallèle à la digue. Par conséquent, le dépôt des sédiments se fait principalement dans la retenue, avant que les eaux de ruissellement n’arrivent au déversoir.

Au niveau du lac de Daimos, en amont du bassin versant, les concentrations en matières solides, mesurées grâce à des prélèvements effectués en 2005 au niveau du déversoir tout au long d’une crue, ont varié entre 0,01 g•L‑1 et 0,28 g•L‑1, avec une concentration moyenne de 0,10 g•L‑1. On peut donc en conclure que la concentration en sédiments des eaux qui sortent du lac de Daimos est nettement supérieure à celle des eaux qui sortent du barrage de Saboun, tout en restant relativement faible.

2.5 Cartographie des unités hydrologiques homogènes

Pour la mise en oeuvre du modèle « Mina », il nous a fallu définir des surfaces de sol à comportement hydrologique homogène, un fonctionnement hydrologique homogène étant défini par un ensemble de paramètres caractéristiques de la surface du sol et agissant sur son comportement hydrodynamique qui règle les conditions d’infiltration et de ruissellement de cette surface. En zone semi-aride méditerranéenne, ces paramètres sont les suivants (EL GAROUANI, 2001; GOMER, 1994; HOUSSA et al., 1998) :

  • la nature du sol, par ses caractéristiques texturale et structurale;

  • l’occupation du sol qui définit les variations saisonnières de son état de surface; et

  • l’état de surface du sol avec la présence ou l’absence de croûtes et d’organisations pelliculaires de surface.

Des campagnes d’échantillonnage des sols et de mesure de ces paramètres sur le terrain ont été réalisées entre 1997 et 1999. Une photographie aérienne (échelle 1:20 000 datant de 1996) a servi de support pour la localisation des parcelles et la spatialisation des paramètres. Les différentes couches d’information ainsi obtenues ont été superposées dans un système d’information géographique (SIG), ce qui nous a permis d’identifier six classes de sols caractérisant le fonctionnement hydrologique des sols du bassin versant de l’Oued Saboun. La répartition spatiale des ces unités hydrologiques a été considérée inchangée dans le temps, au moins pour les deux années d’études, car le parcellaire du bassin est complètement borné et la rotation des cultures comprend deux années de cultures céréalières et une année en jachère (GAY, 2004b; IDRISSI, 2004). Le tableau 3 donne les éléments des différentes couches d’informations retenues dans le SIG et la figure 5 donne le résultat de leur superposition.

Tableau 3

Caractéristiques des principales unités hydrologiques identifiées dans le bassin versant de Saboun.

Characteristics of the principal hydrological units identified in the catchment area of Saboun.

Type de sol

Occupation des sols

Texture

Structure

Croûte de surface

Classes hydrologiques

Peu évolué d’apport alluvial

Sol nu

Argileuse à charge grossière < 40 %

Sub-anguleux peu nette

Structurale

S1

Peu évolué d’apport alluvio-colluvial (terrasse)

Sol nu en jachère

Limoneux-argileuse

Sub-anguleux peu nette

Grossière

S2

Peu évolué d’apport colluvial à caractère vertique

Sol cultivé en jachère

Limoneux-argileuse

Fragmentaire polyédrique très nette

Battance

S3

Vertisol des bas fonds

Sol nu en jachère

Argileux-limoneuse

Fragmentaire polyédrique très nette

Pas de croûte

S4

Peu évolué d’apport colluvial à caractère vertique et calcimagnésique vertique

Sol labouré cultivé orge et fève

Limoneux-argileuse

Massive polyédrique nette

Pas de croûte

S5

Sols minéraux bruts

Affleurement rocheux et habitat

Cailloux et blocs

Grossière

Pas de croûte

S6

-> Voir la liste des tableaux

Figure 5

Carte des unités hydrologiques homogènes de bassin versant Saboun.

Map of the homogeneous hydrological units in the Saboun watershed.

Carte des unités hydrologiques homogènes de bassin versant Saboun.

-> Voir la liste des figures

2.6 Évaluation des paramètres du modèle « Mina »

2.6.1 Caractéristiques topographiques et hydrographiques

Les données topographiques, extraites de la carte au 1:50 000 d’El Manzla, ont été utilisées pour réaliser un Modèle Numérique de Terrain (MNT) formé de mailles carrées ayant 25 mètres de côté, une altitude étant affectée à chaque maille. Ces données ont été ensuite intégrées dans un Système d’Information Géographique (SIG) pour permettre la détermination des caractéristiques topographiques et hydrographiques du bassin versant, principalement la pente, le sens de l’écoulement superficiel et la taille des impluviums (FAIRFLEID et LEYMARIE, 1991). Le MNT utilisé est présenté sur la figure 6.

Figure 6

Modèle numérique d’altitude de bassin versant Saboun.

Digital elevation model of the Saboun watershed.

Modèle numérique d’altitude de bassin versant Saboun.

-> Voir la liste des figures

En commençant par le point le plus élevé et en suivant les mailles dont l’altitude est la plus faible, il est possible d’établir l’organisation générale du drainage des eaux de ruissellement en tout point du bassin versant et de simuler l’arborescence complète de la circulation des eaux superficielles (GARBRECHT et MARTZ, 1999). Chaque maille possède ainsi un seul sens d’écoulement et une direction d’écoulement vers l’aval qui est celle du gradient d’altitude. La simulation de l’écoulement sur chaque maille est approchée par l’équation de continuité monodimensionnelle de Saint-Venant (onde cinématique), le ruissellement étant calculé par différence entre la hauteur d’eau précipitée et la hauteur d’eau infiltrée.

2.6.2 Modélisation de l’infiltration

En tenant compte du phénomène de battance sur sols marneux, le modèle « Mina » exprime l’infiltration par l’infiltration relative (fr), rapport de la vitesse d’infiltration réelle (f) et de l’intensité de la pluie (Ip). L’infiltration réelle s’exprime en appliquant le modèle d’infiltration de HORTON (1933) en fonction de l’énergie cinétique de BRANDT (1989).

où fr est le taux d’infiltration relative, f la vitesse d’infiltration en mm•h‑1, f la vitesse d’infiltration finale en mm•h‑1, f0 la vitesse d’infiltration initiale en mm•h‑1, w le coefficient de courbature ou de forme, Ec l’énergie cinétique en j•m‑2.

Selon Gomer et al. (1992) et GOMER (1994), cette méthode permet de décrire de manière très exacte le comportement des sols développés sur un substrat marneux au moyen du modèle de Horton. Dans ces conditions, l’infiltration des sols dépend, d’une part, de l’intensité des précipitations et, d’autre part, de l’humidité initiale qui est représentée dans l’équation 2 par le coefficient de courbature ou de forme (w). La variation (f0‑ f) est considérée ici comme constante pour chaque unité hydrologique, puisque f0 et f sont des paramètres qui sont liés aux caractéristiques du sol, indépendantes de l’humidité initiale dans la modélisation de Horton.

Les données nécessaires à la résolution de cette équation ont été déterminées de façon expérimentale en utilisant les résultats d’expériences de simulation de pluies. Dix tests de simulation de pluies ont été réalisés à raison de un à deux tests par unité hydrologique homogène, en tenant compte de l’étendue spatiale de chacune des unités. L’intensité des précipitations appliquée fut maintenue constante, à une hauteur de 100 mm•h‑1 pendant une heure. Les données de l’infiltration pour chaque unité de surface ont servi ensuite à l’ajustement des paramètres du modèle de Horton en utilisant la méthode des moindres carrés. Ainsi, on a pu définir, pour chaque unité hydrologique, les vitesses d’infiltration réelle (initiales et finales) et les taux d’infiltration relative (initiale et finale) en divisant par l’intensité de la pluie (100 mm•h‑1), ainsi que les coefficients de forme (w). Les résultats de ces ajustements sont présentés dans le tableau 4 et la figure 7 donne un exemple graphique d’ajustement.

Tableau 4

Vitesses d’infiltration ajustées par le modèle de Horton, ainsi que l’humidité et le coefficient de forme (w) relatif pour chaque unité hydrologique.

Rates of infiltration, adjusted by the Horton model, as well as moisture and the relative form coefficient (w) for each hydrological unit.

Unité hydrologique

f initiale

(fin)

(mm•h‑1)

f finale

(ffin)

(mm•h‑1)

Coefficient de forme

w

Humidité initiale

vol %

Humidité finale

vol %

S1

102

18,7

0,002

31

34

S2

111

8,1

0,116

27

30

S3

103

31,4

0,013

25

33

S4

112

51,3

0,016

24

29

S5

101

5,2

0,011

16

39

S6

102

15,0

0,153

23

27

-> Voir la liste des tableaux

Figure 7

Ajustement des données d’infiltration sous pluie simulée par le modèle de Horton (unité S5).

Adjustment of infiltration data with simulated rain, Horton model (unit  S5).

Ajustement des données d’infiltration sous pluie simulée par le modèle de Horton (unité S5).

-> Voir la liste des figures

Les vitesses finales les plus faibles sont enregistrées pour les unités hydrologiques S2 et S5. L’unité S2 représente les surfaces des sols peu évolués d’apport alluvio-colluvial à texture limono-argileuse en présence d’une croûte dure de type grossière, ce qui réduit considérablement l’infiltration du sol. Son infiltration finale (f) est atteinte rapidement (w = 0,116), ce n’est pas le cas de l’unité S5 qui représente les sols cultivés à texture limono-argileuse dont le temps de passage de (f0) à (f) est très lent (w = 0,011). L’unité S5 regroupe les terrains sur sols à caractère vertique. Elle présente les mêmes caractéristiques pédologiques que l’unité S3, mais l’unité S3 diffère de l’unité S5 par l’existence d’un labour récent sur l’unité S3, ce qui augmente la porosité du sol et accroît très fortement le potentiel d’infiltration du sol.

Pour obtenir les valeurs de l’infiltrabilité initiale pour différents états d’humidité des sols, nous avons utilisé la relation établie par GOMER (1994) pour les sols algériens entre l’infiltrabilité initiale et l’indice d’antériorité des pluies (IK) déterminé selon la formule de KOLHER et LINSLEY (1951) :

où IKn est l’indice du jour n exprimé en mm, Pn-1 la pluie du jour n-1 exprimée en mm et α un coefficient de décroissance exponentielle fonction du type de sol et du climat exprimé en jour-1.

2.6.3 Arrachement des particules

L’arrachement des particules du sol peut être liée à l’action de deux forces de cisaillement : celle liée à l’énergie cinétique des gouttes de pluie lors de leur chute et celle liée à l’énergie cinétique du ruissellement lors de l’écoulement. L’érosion se réalise à condition que ces deux forces dépassent la résistance au cisaillement spécifique du type du sol, appelée aussi énergie cinétique critique de détachabilité des particules. Pour une maille i, de longueur l et de largeur b, le débit solide disponible est exprimé par l’équation suivante (SCHMIDT, 1992) :

avec :

  • IQi = ρflu Qi Vi

  • Irαi = sin(α) ρflu cos (α) ri li bi Vri (1-Vg)

  • Vri = 4,506 + 0,615 log (ri cos(α)).

où Qs est le débit solide (kg•s‑1), Qi le débit solide pour une maille i (kg•s‑1), li la longueur de la maille i (m), bi la largeur de la maille i (m), ri l’ intensité de précipitation (m•s‑1), IQi l’impulsion cinétique de l’écoulement (N•m2), Irαi l’impulsion cinétique efficace de la pluie (N•m2), mcr l’impulsion cinétique critique (N•m2), ρflu la densité du fluide, Vi la vitesse d’écoulement moyenne (m•s‑1), Vri la vitesse moyenne de chute (m•s‑1), α l’angle de la pente (degrés) et Vg le taux de couverture végétal au sol (%).

En pratique, ces seuils de détachabilité (impulsions cinétiques) sont difficiles à déterminer. Ils le sont par des essais de simulation de pluie sur des parcelles de grande taille, pour des longueurs de 22 m. Les résultats sont ensuite intégrés dans les équations ci-dessus pour le calcul de l’impulsion cinétique critique (mcr). En l’absence de ce type de simulateur sur le bassin versant de Saboun, les valeurs de l’impulsion cinétique critique de détachabilité pour chaque unité de surface du bassin versant Saboun ont été déduites des valeurs déterminées pour les mêmes types de sols sur le bassin versant de l’Oued Mina (Algérie) et figurant au tableau 5. Les valeurs de l’impulsion cinétique critique de détachabilité sont égales à : 0,0003 pour les sols peu évolués, 0,001 pour les sols très filtrants (vertisols et sols limono-argileux), 0,002 pour les sols argileux à charge grossière et 0,003 pour les bads-lands pierreux.

Tableau 5

Valeurs d’impulsion cinétique critique déterminées pour les sols du bassin versant d’Oued Mina et par analogie pour les surfaces du sol du bassin versant Saboun.

Critical kinetic impulse values for Wadi Mina catchment area soils and, by analogy, for the soil surface of the Saboun catchment.

Type de sol bassin versant

Oued Mina (GOMER, 1994)

Type de surface

bassin versant Saboun

Impulsion cinétique critique

N•m‑2 (d’après GOMER, 1994)

Vertisol très filtrant stable

S4

0,0010

Sol limono-argileux très filtrant

S3, S5

0,0010

Sol rouge très filtrant

-

0,0010

Sol argileux à charge grossière

S2

0,0020

Sol peu évolué

S1

0,0003

Sol argileux salin

-

0,0030

Sol salin très ruisselant

-

0,0030

Sol salin, ruissellant

-

0,0003

Sol argileux salin très ruisselant compact

-

0,0003

Bad-lands peu pierreuses

-

0,0001

Bad-lands pierreuses

S6

0,0030

-> Voir la liste des tableaux

2.6.4 Transfert des débits liquides et solides jusqu’à l’exutoire

Le routage et le transfert des eaux entre les mailles sont basés sur l’équation de continuité monodimensionnelle appelée aussi deuxième équation de Saint-Venant.

Pour distinguer entre un écoulement en nappe et un écoulement concentré en chenaux, nous avons utilisé la notion de surface limite d’apparition des rigoles. Sur les marnes tertiaires du bassin versant de l’Oued Mina, cette surface est évaluée à 1 000 m2 au moyen d’un Système d’Information Géographique par confrontation de la carte des aires contributives à la carte du réseau de drainage (GOMER et al., 1992). Dans le bassin versant de l’Oued Saboun, sur les marnes du Miocène, elle a été évaluée par la même méthode à 1 200 m2 qui correspond à une longueur de la pente comprise entre 30 et 40 m à partir de la ligne de partage des eaux.

Le transfert des particules solides, par analogie avec l’équation de continuité pour le transfert de la phase liquide, se fait de la même manière en appliquant l’équation de continuité adaptée à la phase solide. Les dépôts et les reprises sont gérés par le calcul d’un bilan entre la sédimentation et l’érosion, fonction du coefficient de dépôt Cdep qui varie entre 0 (dépôt) et 1 (érosion). Ce coefficient de dépôt est calculé comme suit :

où h est la hauteur de la lame d’eau, ω la vitesse de sédimentation, ∆x la distance et v la vitesse de flux liquide.

3. Résultats et discussions

Le modèle « Mina » a été appliqué pour toutes les averses enregistrées à la station Saboun au cours des deux années hydrologiques 1997-98 et 1998-99. Nous avons choisi ces deux années hydrologiques parce qu’elles permettent de bien encadrer le régime hydrologique du bassin versant de l’Oued Saboun. En prenant les deux stations de Tanger DPE et Kalaya comme stations de référence, avec une hauteur pluviométrique annuelle de 797 mm, l’année 97-98 présente une période de retour de trois ans humides, alors que l’année 98-99, pour une hauteur pluviométrique annuelle de 444 mm, présente une période de retour de vingt ans secs.

La simulation des débits liquide et solide a été réalisée au pas de temps de 5 min, identique au pas de temps d’enregistrement des hauteurs d’eau du lac par le limnigraphe. Quant aux caractéristiques topographiques et pédologiques, elles ont été définies pour chaque maille carrée de 25 m de côté par l’utilisation d’un SIG intégrant le MNT et la carte pédologique du bassin versant.

Le modèle est appliqué de façon distincte sur les parties amont et aval du bassin versant :

  • la partie amont correspond au bassin versant de Daimos d’une superficie de 228 ha;

  • la partie aval correspond au bassin versant de Saboun sans le bassin versant de Daimos, sa superficie est de 574 ha.

Cette désagrégation du bassin versant de l’Oued Saboun en deux parties imbriquées est appliquée de manière à introduire une fonction de stockage de l’eau et des sédiments dans le lac de Daimos. Ne disposant pas de mesures hydrométriques à Daimos, nous sommes partis sur l’hypothèse que le lac était vide en début de simulation (1er octobre 1997), qu’il se remplit à chaque crue et déverse à partir d’un volume de 30 000 m3. Le niveau de remplissage de la retenue de Daimos avant l’averse tient compte de son niveau de remplissage à la fin de l’averse précédente et de l’évaporation de l’eau stockée depuis l’averse précédente. Nous avons calculé le volume journalier évaporé sur le lac de Daimos en multipliant l’évaporation journalière mesurée au bac Colorado enterré près du barrage de Saboun par la surface du plan d’eau au barrage Daimos le même jour. Pour les lacs collinaires de la zone semi-aride, l’évaporation mesurée au bac Colorado enterré est en effet peu différente de l’évaporation du lac (GAY, 2004a). L’infiltration dans la retenue de Daimos a été négligée.

La figure 8 compare les crues entrant dans les lacs de Saboun et Daimos à la date du 03/02/1998 et la figure 9 montre la reconstitution, par le modèle « Mina », des débits entrant dans le barrage de Saboun avant et après calibration du paramètre d’infiltrabilité initiale pour l’événement du 18/12/1997.

Figure 8

Crue du 03/02/1998 (en violet, le débits simulés sortant du lac Daimos ; en rouge les débits reconstitués entrant au barrage Saboun)

Flood of 03/02/1998 (in purple, simulated flows leaving lake Daimos; in red, reconstituted flows entering the Saboun reservoir).

Crue du 03/02/1998 (en violet, le débits simulés sortant du lac Daimos ; en rouge les débits reconstitués entrant au barrage Saboun)

-> Voir la liste des figures

Figure 9

Comparaison des débits observés et simulés entrant au barrage Saboun pour la crue du 18/12/1997, avant et après calibrage.

Comparison of observed and simulated flows entering Saboun reservoir for the flood of 18/12/1997, before and after calibration.

Comparaison des débits observés et simulés entrant au barrage Saboun pour la crue du 18/12/1997, avant et après calibrage.

-> Voir la liste des figures

La figure 10 présente les écarts entre les volumes observés et les volumes simulés pour l’ensemble des crues en fonction de l’indice d’antériorité des pluies (IK) avant calibration et après calibration du paramètre d’infiltrabilité. Le tableau 6 présente les résultats de cette calibration en matière de volume d’eau entrant au barrage et d’envasement de la retenue en fonction des volumes de sédiments déversés.

Figure 10

Comparaison entre les volumes des crues simulées et observées en fonction de l’indice IK avant et après correction d’infiltrabilité.

Comparison between volumes of simulated and observed flows according to index IK, before and after correction of the permeability.

Comparaison entre les volumes des crues simulées et observées en fonction de l’indice IK avant et après correction d’infiltrabilité.

-> Voir la liste des figures

Tableau 6

Apports en eau et en sédiment simulés crue par crue à l’exutoire du barrage Saboun entre 1997 et 1999.

Water and sediment fluxes, simulated for individual flood events, at the Saboun dam outlet between 1997 and 1999.

Date

Pr

(mm)

I5 max

mm•h‑1

IK

Volumes ruisselés mesurés

(m3)

Volumes déversés mesurés

(m3)

Volumes ruisselés simulés

(m3)

Sédiments simulés

(kg)

Sédiments déversés estimés

(kg)

Volume d’envasement de la retenue, estimé

(m3)

Masse d’envasement de la retenue estimé

(kg)

07/11/97

15,5

21

0

9 382

0

3 974

968 418

0

745

968 418

18/11/97

18,5

66

0

4 149

0

4 866

1 358 708

0

1 045

1 358 708

21/11/97

37,5

48

20

65 333

0

3 482

699 049

0

538

699 049

24/11/97

16

30

13

11 770

0

2 683

654 732

0

504

654 732

25/11/97

31,5

24

15

18 152

0

3 812

1 294 247

0

996

1 294 247

4‑5/12/97

47

36

34

239 694

146 830

114 278

4 344 404

8 810

3 322

4 318 600

16/12/97

39

48

2

16 504

1 720

9 111

2 746 957

103

2 113

2 746 900

17‑18/12/97

82,5

78

36

181 453

82 615

121 307

3 111 806

4 957

2 382

3 096 600

19/12/97

5

24

43

12 848

64 860

3 326

445 431

3 892

334

434 200

22/12/97

24

36

23

68 450

32 354

56 372

1 836 060

1 941

1 408

1 830 400

12/01/98

24,5

36

36

43 340

17 194

35 804

1 313 880

1 032

1 008

1 310 400

01/02/98

11,5

20

24

22 061

16 084

7 083

126 571

965

95

123 500

3‑4/2/98

70

49

40

251 731

234 486

97 760

2 839 833

14 069

2 152

2 797 600

06/02/98

12

14

20

68 732

55 805

42 061

1 229 083

3 348

938

1 219 400

26/09/98

60

84

6

13 850

0

204 570

11 171 031

0

8 593

11 171 031

29‑31/12/98

65

54

0

11 106

0

7 882

860 497

0

662

860 497

12/01/99

36,5

30

14

3 824

0

11 410

3 498 475

0

2 691

3 498 475

17/01/99

17,5

18

5

4 109

0

9 456

315 626

0

243

315 626

21/01/99

23

12

0

4 480

0

4 442

569 262

0

438

569 262

23/01/99

20,5

12

5

9 188

0

3 844

498 370

0

383

498 370

24‑25/02/99

25

30

37

11 180

0

7 277

725 631

0

558

725 631

28/02/99

10

12

4

1 202

0

1 642

177 561

0

137

177 561

24/03/99

53

78

5

24 218

0

6 626

642 575

0

494

642 575

01/05/99

14

12

1

3 895

0

2 637

263 470

0

203

263 470

Pr : précipitation ; I5 : intensité maximale de pluie en 5 min ; IK : indice d’humidité de Kolhler.

-> Voir la liste des tableaux

La confrontation des résultats simulés et observés montre que le modèle sous-estime l’écoulement pour les averses à fort indice d’humidité des sols et surestime l’écoulement pour les averses isolées à indice d’humidité faible, proche de 0. Ceci est expliqué par la complexité des processus qui génèrent le ruissellement sur les sols du bassin versant de Saboun. Pour les sols à texture limono-argileuse du bassin versant de Saboun, comme les unités S2, S3 et S5, le modèle de Horton peut être appliqué sans inconvénient car les simulations de pluies ont révélé que ce type de sol possédait les taux les plus faibles d’infiltration avec des coefficients de ruissellement atteignant 70 %. En effet, l’aptitude au ruissellement des sols à texture limono-argileuse est liée principalement à la dégradation structurale de la surface du sol sous l’action de la pluie et à la formation de la croûte de battance. Cependant, le caractère vertique de certains sols (unités S3 et S5) accroît leur capacité à l’infiltration lorsque les sols sont secs, en raison de l’existence de fentes de retrait plus ou moins larges et profondes selon l’humidité du sol. Pour les vertisols, stables et filtrants (unité S4), ils se comportent en début de saison des pluies comme des réservoirs ayant une certaine capacité de stockage, leur capacité au ruissellement dépend donc de leur teneur en eau et non de la formation de la battance. Pour une averse qui survient dans des conditions antérieures d’humidité faible sur un sol argileux, les fentes et fissures entraînent une importante infiltration de la pluie dans le sol avant le ruissellement (GERMANN, 1990; MRABET, 1989), favorisant beaucoup plus des écoulements hypodermiques que des écoulements de surface directs sur les versants. Ce phénomène a été également observé par Heusch (1970) sur le bassin versant de Mda (région de Ouezzane) sur le même type de sol (sol marneux) du Rif occidental. Il a constaté que plus de 70 % des eaux écoulées sur les versants, rejoignent les talwegs par des écoulements hypodermiques liés principalement à la fissuration. Ces écoulements sont plus importants en année sèche qu’en année humide.

Le caractère plus ou moins argileux et vertique des sols du bassin versant de l’Oued Saboun provoque donc une sous-estimation de l’infiltration pour les averses isolées dont l’indice des pluies antérieures est faible (IK proche de 0). Par contre, pour une averse qui survient dans des conditions d’humidité des sols importante (IK supérieur à 15 mm), le sol très humide perd la cohésion de ses particules. Il s’ensuit un effondrement de la structure macroscopique et la disparition d’une partie des pores assurant l’infiltration et la rétention d’eau. L’espace poral ou intergranulaire n’est plus interconnecté, ce qui entraîne une réduction de l’infiltration et sa surestimation par le modèle de Horton.

En matière d’érosion, les résultats de la simulation du transport solide ont été comparés aux données de l’envasement issues de la bathymétrie. L’envasement calculé par le modèle est estimé à 22 850 tonnes pour l’année 1997-1998 et 18 700 tonnes pour l’année 1998-1999, alors que la moyenne interannuelle estimée par bathymétrie sur la période 1999-2004 est de 27 750 tonnes. Si on tient compte des précipitations des années 2001-2002 et 2003-2004, exceptionnellement humides et dépassant 1 000 mm, l’érosion semble être correctement simulée par le modèle « Mina ». Les résultats de cette comparaison montrent donc que la transposition sur le bassin versant de l’Oued Saboun des impulsions cinétiques critiques déterminées sur les sols du bassin versant de l’Oued Mina pour les mêmes types de sols expriment bien les processus érosifs sur le bassin versant de l’Oued Saboun.

À partir de ces résultats de simulation, nous avons cherché à reconstituer le remplissage et l’envasement du barrage à l’échelle mensuelle sur la période comprise entre septembre 1997 et août 1999. Les apports mensuels en eau au barrage sont calculés en sommant les volumes des crues; l’envasement mensuel est calculé en sommant les volumes du transport solide simulé dont on soustrait la part des sédiments sortant par le déversoir. Nous divisons le tonnage simulé par la densité moyenne mesurée de sédiment (dmoy = 1,3) pour obtenir le volume sédimenté dans le barrage. La part des sédiments déversés est obtenue par le produit du volume déversé et de la concentration moyenne en matières solides des déversements. Les sorties en eau du barrage sont, d’une part, l’évaporation et, d’autre part, l’utilisation de l’eau du lac pour l’alimentation du périmètre irrigué. L’évaporation est calculée à l’échelle journalière à partir des mesures au bac Colorado enterré près de la digue du barrage et l’utilisation de l’eau pour l’irrigation est calculée à partir de la consommation électrique de la station de pompage. Les volumes mensuels utilisés, prélevés dans le lac de Saboun, sont consignés au tableau 7.

Tableau 7

Utilisation mensuelle de eau du barrage Saboun en m3.

Monthly use of water from the Saboun reservoir in m3.

 

Sept

Oct

Nov

Déc.

Janv.

Fév.

Mars

Avril

Mai

Juin

Juillet

Août

Année pluvieuse

1997/1998

0

0

0

0

0

0

0

0

7 300

18 100

8 200

15 000

Année sèche

1998/1999

14 300

7 500

8 200

6 700

5 000

3 600

9 800

14 600

9 900

11 000

35 000

65 000

-> Voir la liste des tableaux

Les résultats de cette reconstitution sont présentés à la figure 11. Ils montrent une bonne corrélation entre les données de remplissage mesurées et simulées au barrage de Saboun. En ce qui concerne l’envasement, le contrôle bathymétrique effectué dans le lac en octobre 1998 correspond assez bien à l’envasement déterminé par simulation à la même date. En outre, la simulation montre que les apports en sédiments ont eu lieu principalement au cours de la période hivernale de décembre 1997 à février 1998.

Figure 11

Variation de l’état de remplissage du barrage Saboun de novembre 1997 à octobre 1999 : (comparaison entre les données simulées en rouge et les données mesurées en bleu).

Variation of the state of filling of the Saboun’s dam from November 1997 to October 1999: (comparison of simulated data, in red and measured data, in blue).

Variation de l’état de remplissage du barrage Saboun de novembre 1997 à octobre 1999 : (comparaison entre les données simulées en rouge et les données mesurées en bleu).

-> Voir la liste des figures

4. Conclusion

Cette étude avait pour objectif de simuler les apports en eaux et en sédiments au barrage collinaire de l’Oued Saboun en utilisant le modèle « Mina », modèle hydrologique distribué à base physique basé sur l’équation d’infiltration de Horton pour le calcul de l’écoulement et de l’énergie cinétique critique de détachabilité des particules du sol pour le calcul de l’érosion.

La confrontation des résultats aux données observées a montré, pour un calage des paramètres d’infiltration basé sur des expériences de simulation de pluies effectuées au printemps (mars et avril 1999), une nette surestimation des volumes écoulés simulés en début de saison des pluies, pour de faibles indices d’antériorité des pluies correspondant à des sols secs. Cet écart peut être atténué par la prise en compte de la présence de fissures dans les sols marneux en augmentant le taux d’infiltration initial dans l’équation d’infiltration de Horton. Inversement, la confrontation des résultats aux données observées a montré une nette sous-estimation des volumes écoulés simulés pour les forts indices d’antériorité des pluies correspondant à des sols très humectés. Cet écart peut être atténué par une réduction du taux d’infiltration final dans l’équation de Horton. Ces améliorations ont permis une meilleure adaptation du modèle « Mina » au comportement hydrodynamique des sols du bassin versant de l’Oued Saboun.

Pour améliorer la qualité des ajustements de l’équation de Horton aux sols argileux et limono-argileux vertiques, il semble donc nécessaire de développer des expériences de simulation de pluies à différentes périodes de l’année en tenant compte des caractéristiques hydriques initiales des sols qui déterminent leurs capacités à l’infiltration. La prise en compte d’une évolution saisonnière de cette capacité d’infiltration, en relation avec l’humidité du sol, semble indispensable pour rendre compte de manière satisfaisante du fonctionnement hydrologique du bassin versant de l’Oued Saboun et probablement de tous les bassins versants du sud de la Méditerranée qui possèdent les mêmes types de sols.