Etude de l'impact des barrages sur la réduction des transports sédimentaires jusqu'à la mer par approche paléohydrologique dans la basse vallée de la Medjerda

The sedimentary contributions of the Medjerda to the coastal zone are poorly measured, and there is no chronicle of observations. In this context, the sediment monitoring appears indispensable for the quantification of sediment transport at the outlet. This study focuses on the largest watershed in Tunisia, the Wadi Medjerda (23 600 km²). The main objective of this work is to assess the reduction of sediment transport following anthropogenic intensification on the basin, especially since the construction of many large dams. In order to collect information on actual deposits over several years, the paleo-hydrological approach was applied through the study of sediment cores sampled in the low valley meanders on alluvial terraces, after the last dam (Sidi Salem, the largest water storage capacity over the basin), but before the estuary to avoid marine influence and near a hydrological station (Jdaida). The sedimentary deposits of the river provide key information on the past sedimentary inputs. A visible succession of sedimentary layers corresponding to the deposits of successive floods on the study site has been determined and the history of the sedimentary contributions of the Medjerda is reconstructed by this approach. The thickest layers of sedimentary deposits are related to exceptional events. They are mainly concentrated on the lower part of the core and are mainly composed of sands. The first 1.2 m of the core from the bottom upward relates to 10 years of river discharges, as can be determined from the ¹³⁷Cs datation. The next upward 1.05 m of core relates to the following 20 years of discharges, up to 1981, date of the construction of the Sidi Salem dam, and is composed of a mix of sand, silts and clays. The last 75 cm of core near the surface is only composed of clays with thin silt bands, and relates to a period of 32 years. We thus observe that there is no more sand deposits in the river bed since the construction of the Sidi Salem dam. The deficit of sediment supply to the sea is viewed as a major factor to be taken into account for better understanding of the dynamics of coastal areas in the context of global climate change.

1 Introduction

Dans le contexte actuel des changements climatiques, le suivi du transport sédimentaire dans les rivières est important pour mettre l'accent sur les problèmes environnementaux souvent reliés à des changements observés dans le cycle hydrologique  : modification de la distribution des précipitations, du ruissellement et de l'intensité des extrêmes hydrologiques (Bates et al., 2008 ; Hungtinton, 2006 ; Labat et al., 2004 ; Solomon, 2007). Le suivi du transport sédimentaire permet également d'évaluer et de quantifier les transferts de matières des continents vers les océans. Cependant, le problème du nombre réduit de données relatives aux transports solides a toujours freiné les investigations. Les mesures de transport solide effectuées en Tunisie et les études ponctuelles détaillées restent insuffisantes pour traiter ce sujet (Ghorbel et Claude, 1982). Le défi consiste donc à utiliser des méthodes de mesures simples pour combler le manque de connaissances. La présente étude s'intéresse au fleuve Medjerda en Tunisie, afin d'étudier en détail l'évolution des transports solides à l'exutoire de son bassin versant. L'Oued Medjerda est devenu un fleuve presque totalement contrôlé depuis l'intensification des aménagements sur le bassin au début du XXème siècle, pour atténuer les effets des crues dans la basse vallée très intensément exploitée. Le transport de matières solides a été considérablement réduit dans sa basse vallée, en particulier depuis l'ouverture du barrage Sidi Salem sur le cours principal de la Medjerda, en 1981, responsable d'une modification majeure du régime hydrologique du fleuve en aval (Zahar et al., 2008). L'alluvionnement des retenues des aménagements constitue un problème récent et majeur (Ben Mamou et Louati, 2007). Or, l'histoire des apports sédimentaires dans cette zone d'étude n'est pas connue, ce qui entraîne deux difficultés  :

• comment évaluer l'apport solide à la mer si l'histoire des apports sédimentaires n'est pas connue  ?

• comment identifier les différents facteurs contrôlant l'apport sédimentaire par la Medjerda  ?

La paléohydrologie peut en partie répondre à ces questions. On tentera dans ce cas d'établir un historique de la sédimentation dans le bief aval du fleuve Medjerda, avant la mer, et de relier cet historique avec l'histoire des débits, l'évolution climatique et la construction des aménagements. La paléohydrologie est la reconstruction de l'ampleur et de la fréquence des inondations récentes et passées en utilisant des archives sédimentaires (Baker et al., 2002). Notre travail va se concentrer principalement sur la Medjerda, près de la station hydrologique de Jdaida située dans la basse vallée mais suffisamment loin de la mer pour ne pas être impactée par la marée, zone d'étude pour laquelle il existe une quantité importante de données historiques, particulièrement sur les débits des 70 dernières années. Notre approche consiste à prélever des carottes de sédiment dans un méandre où les dépôts sont réguliers et peu turbulents, afin de récolter des informations sur les dépôts réels sur plusieurs décennies. Nous avons utilisé des techniques d'analyse stratigraphique et de datation des dépôts (¹³⁷Cs, ²¹⁰Pb, géochimie) pour évaluer la hauteur atteinte par les eaux lors de différents évènements d'inondation, et proposer une chronologie des crues de la Medjerda. L'analyse pluridisciplinaire de ces archives sédimentaires permet de fournir de précieuses informations sur l'évolution climatique des dernières décennies, les événements extrêmes et enfin sur les effets des actions anthropiques passées et récentes (Vilanova et al., 2006 ; Selby et Smith, 2007).

Figure 1 Localisation des grandes stations hydrométriques et des barrages sur le bassin versant de la Medjerda, et site du carottage.

2 Matériel et méthode

2.1 Site d'étude du bassin versant de la Medjerda

Le fleuve Medjerda est le plus long de Tunisie, avec 484 km dont 350 km en Tunisie (Ben Mansoura et al., 2001). Il prend sa source près de Souk-Ahras, dans le Constantinois algérien, puis coule vers l'est avant de se jeter dans la mer Méditerranée dans le golfe de Tunis. La Medjerda draine 23 600 km² dont 7500 km² en Algérie (Habaieb, 1992) (Fig. 1).

Avant de se jeter dans la mer, la Medjerda traverse trois types de plaines  :

• en amont, la haute vallée (du côté de Ghardimaou) est constituée de bassins colmatés par des sédiments quaternaires. Le fleuve et ses affluents sont caractérisés par des méandres encaissés de 8 à 10 m. Ces régions appartiennent à des zones faiblement montagneuses. Les fortes crues ne débordent qu'exceptionnellement en dehors du lit du fleuve (Kotti et al., 2016) ;

• en aval de la retenue de Sidi Salem, de Testour à Jdaida, dite moyenne vallée de la Medjerda, le fleuve est encaissé dans la plaine. L'encaissement repasse à plus de 10 m en aval de Laroussia avant de s'atténuer rapidement vers Jdaida. Dans tout ce secteur la plaine est inondable surtout dans la partie centrale. En allant vers la mer les reliefs de ce bassin s'adoucissent ;

• dans la basse plaine, en aval de Jdaida, commence le delta. La pente de la Medjerda est encore marquée et son tracé a fréquemment été modifié, soit naturellement soit par des travaux (coupure de méandres, ouverture de nouvel émissaire…). Du point de vue géologique, la Medjerda est occupée par des alluvions quaternaires le plus souvent argilo-sableuses à argilo-sablo-limoneuses.

La pluviométrie moyenne annuelle dépasse 1000 mm au nord du bassin versant de la Medjerda et diminue progressivement vers le sud. Il s'agit d'une pluviométrie souvent concentrée dans le temps, irrégulière et inégalement répartie sur les différentes saisons (Kotti et al., 2016). La base des données hydrologique provient en grande partie de la Direction Générale des Ressources en Eau. Par ailleurs le laboratoire HydroSciences Montpellier (HSM) dispose d'une base de données pluviométrique importante au pas d'espace du demi-degré carré, qui propose des valeurs mensuelles de pluie par maille de 2500 km² sur la période 1940–1999, complétée par des données nationales (DGRE).

2.2 Régime hydrologique et stratégie de mobilisation des eaux de surface

Le régime hydrologique est typiquement méditerranéen. Les contrastes saisonniers sont très marqués puisque, pendant l'été, le débit d'étiage peut se réduire à moins d' 1 m³ s⁻¹ tandis que, pendant le reste de l'année, celui d'une crue de périodicité décennale atteint 1000–1200 m³ s⁻¹ (Rodier et al., 1981). Le régime hydrologique de ce fleuve très hétérogène a été régulé par de nombreux ouvrages construits sur son cours principal ou ses affluents. A partir de 1945, un programme d'équipement de la Medjerda a été entrepris par la Direction des Travaux Publics de Tunisie, destiné à mettre en valeur les ressources encore inutilisées de cette vallée. Ce programme a permis la construction du barrage régulateur de Mellegue et la construction du barrage de Beni Metir (Strugo, 1955). Depuis l'indépendance de la Tunisie, le programme de développement du bassin de la Medjerda a été poursuivi grâce à la réalisation de plusieurs ouvrages  : Lakhmes (1966), Kasseb (1968), Bouhertma (1976), Siliana (1987), Sidi Salem (1981) et R'Mil (1999). En Algérie, un seul grand barrage a été identifié, Ain Dalia, mis en eau en 1988 et d'une capacité de 82 Mm³. Avec la mise en eau de ces barrages, le régime hydrologique du fleuve Medjerda a été modifié (Fig. 1 pour la localisation et tableau 1).

Tableau 1 Caractéristiques des aménagements construits sur le bassin versant de la Medjerda (DGBGTH, 2016).

^∗ Données manquantes.

2.3 Technique de prélèvements des sédiments, méthodes de mesure et géochronologie

Le site du carottage a été choisi après plusieurs missions de terrain (2013, 2014, 2015) qui avaient pour objectif d'identifier différents sites potentiellement intéressants pour la thématique. Le site d'étude est localisé dans un méandre en aval du dernier barrage -Sidi Salem- mais en amont de la zone estuarienne, proche de la station hydrologique Jdaida, la dernière avant la mer. Trois carottes ont été échantillonnées sur cette même terrasse alluviale, C1 de 90 cm, C3 de 120 cm et C5 de 300 cm de profondeur. L'approche paléohydrologique a été utilisée pour comparer les résultats sédimentologiques et géochronologiques de cette archive sédimentaire aux informations de la station hydrométrique de Jdaida située juste en amont du site d'étude choisi.

L'identification et la mise en évidence des différents événements de crues ont été réalisées grâce à l'inspection minutieuse de chaque dépôt (décompte des lamines, identification de paléosols, détection des surfaces d'érosion, étude de la variabilité granulométrique à partir d'un granulomètre laser Beckman Coulter LS 13320 au laboratoire Géosciences Montpellier.

Toute reconstitution paléo-environnementale, à partir de l'étude des archives sédimentaires, nécessite l'établissement d'une échelle chronologique la plus fine possible, afin de dater avec une bonne résolution les événements hydrologiques passés. Les isotopes principalement utilisés pour dater les événements du siècle dernier sont le ²¹⁰Pb et le ¹³⁷Cs (Dezileau et al., 2014a). Ils permettent une datation relativement précise, grâce à leurs propriétés chimiques (réactivité vis-à-vis des particules) et physiques (périodes de désintégration, rayonnement gamma). Le ¹³⁷Cs et le ²¹⁰Pb sont mesurés par spectrométrie gamma (détecteur CANBERRA BEGe 3825) au laboratoire Géosciences Montpellier (France) et au CNSTN à Sidi Thabet (Tunisie).

• Datation par le ²¹⁰Pbex (T1/2 = 22,3 ans)

Le Plomb 210 est un isotope radioactif naturel issu de la chaîne de désintégration de l'Uranium 238. L'Uranium 238, contenu naturellement dans la croûte terrestre, subit la série de désintégration suivante  : Uranium 238 – Uranium 234 – Thorium 230 – Radium 226 – Radon 222. Le 222Rn (T1/2 = 3,8 jours) ainsi produit, se trouve sous une forme gazeuse. Lorsqu'il est produit à proximité de la surface, le gaz s'échappe dans l'atmosphère où il se transforme rapidement en Polonium 218 qui se désintègre en quelques minutes en ²¹⁰Pb. Ce ²¹⁰Pb, dit atmosphérique (ou en excès) retombe sous forme de particules humides ou sèches et s'accumule dans les sédiments. Son activité au sein de la colonne sédimentaire diminue avec le temps selon une loi de décroissance radioactive. Le principe de la datation au ²¹⁰Pb s'appuie sur l'hypothèse selon laquelle le flux atmosphérique en ²¹⁰Pb est constant, ce qui implique que l'activité de l'isotope radioactif à la surface des sédiments reste toujours identique. Cette hypothèse permet en mesurant l'activité du ²¹⁰Pb et en appliquant la loi de décroissance radioactive de déterminer l'âge de nos sédiments. Au bout de cinq fois sa période, c'est à dire 110 ans, il a complètement disparu dans l'environnement (Dezileau et al., 2014b). La mesure du ²¹⁰Pb sur le détecteur gamma CANBERRA se fait en intégrant le pic à 46,5 keV. Le ²¹⁰Pbex est estimé en soustrayant le pic du ²¹⁰Pb au pic du ²²⁶Ra à 186,2 keV. L'erreur n'excède pas 6 %.

• Datation par le ¹³⁷Cs (T1/2 = 30,17 ans)

Le Césium 137 est un isotope radioactif artificiel issu de réactions de fission nucléaire. Il est relâché dans l'atmosphère lors de l'utilisation de bombes atomiques et d'accidents dans des centrales nucléaires. Avec une période de demi-vie (30,17 ans) proche de celle du ²¹⁰Pb, il permet de dater les séquences sédimentaires sur une même échelle temporelle. Les profils d'activité du Césium 137 en fonction de la profondeur dans une coupe sédimentaire présentent des pics dont le plus important est généralement attribué à l'année 1963, date du maximum de rejets atmosphériques par les essais nucléaires. La présence de Césium 137 dans les dépôts sédimentaires indiquera qu'ils ont été déposés postérieurement aux années 1950. La mesure du ¹³⁷Cs sur le détecteur gamma CANBERRA se fait en intégrant le pic à 661 keV. L'erreur n'excède pas 6 %.

Figure 2 (a) Le volume de vase évacué au niveau du barrage Mellegue entre 1954 et 2003. (b) Le volume de vase évacué au niveau du barrage Sidi Salem entre 1982 et 2005.

3 Résultats

3.1 Sédimentation des réservoirs

Le taux d'envasement de chaque barrage construit sur le bassin versant de la Medjerda met en évidence l'importance de l'alluvionnement durant la durée de l'exploitation des retenues des barrages, dont Mellegue, Siliana et Sidi Salem qui sont sujettes à un alluvionnement plus important (tableau 1). Il a été estimé que le volume initial de la retenue de Sidi Salem sera envasé jusqu'à 60 % 100 ans après la construction du barrage, et l'envasement devrait être total pour les retenues des deux barrages de Mellegue et Siliana dans 60 à 65 ans (Ben Mammou et Louati, 2007). Des soutirages efficaces permettent dans ce cas de réduire le taux d'alluvionnement et d'augmenter la durée de vie de ces barrages. D'après Abid (1980), 59 à 64 % des apports solides dans la retenue du barrage de Mellegue, durant la période 1954–1980, ont été soutirés grâce aux manœuvres de dévasement pendant les crues. La Fig. 2 montre le volume de vase évacué suite à des événements exceptionnels entre 1954 et 2003 au niveau du barrage Mellegue et entre 1982 et 2005 au niveau du barrage Sidi Salem. En 1969, de très importants volumes de sédiments évacués lors d'une crue se sont accumulés surtout dans les zones où la pente diminue, ainsi que dans les secteurs évasés. En mars 1973, l'oued Medjerda a connu une autre crue plus intense provoquant de graves inondations avec une quantité exceptionnelle de matériaux transportée et ce malgré les gros travaux d'endiguement effectués à Jdaida (entre 1947 et 1957) (Ben Hassine and Rejeb, 2003).

Depuis la mise en eau du barrage Sidi Salem en 1981 jusqu'en 1990, le volume de vase évacué par ce barrage est quasiment nul. On remarque le faible pourcentage de vase évacuée dû à l'arrêt des opérations de soutirage. Quelques opérations de dévasement occasionnel ont été associées à des débits maximum exceptionnels  : 595 m³ s⁻¹ en 1990, 671 m³ s⁻¹ en 1995, 977 m³ s⁻¹ en 2000 et 1020 m³ s⁻¹ 2003. Durant l'année 1995, les gestionnaires de Sidi Salem ont libéré un volume total dévasé de 24,27 Mm³. Les particules solides se déposent en aval en fonction de la réduction de la vitesse du courant de turbidité, ainsi plusieurs facteurs peuvent influencer les épaisseurs des dépôts accumulés tels que la vitesse du courant, la topographie, la présence d'obstacles perturbant les dépôts, etc. Ces facteurs créent une variabilité latérale de dépôt importante, en aval de la retenue. Le lit de la Medjerda connaît de ce fait un engraissement continu depuis 1981, date de la mise en service du barrage Sidi Salem. De ce fait, et par manque de débits de crue important qui “ chassent ” les dépôts précédents, le débit occasionnant le débordement a diminué à moins de 200 m³ s⁻¹, valeur nettement inférieure au débit de débordement avant la construction du barrage qui était de 700 m³ s⁻¹.

La Fig. 2a illustre parfaitement l'accroissement des volumes solides évacués suite à la crue exceptionnelle de 2003. Quatre fortes crues très rapprochées sont survenues dans le bassin de la Medjerda au début de l'année 2003, après quatre années de sécheresse (Abdelhamid et Khalil, 2009). Elles ont cumulé un milliard de m³ d'apports à Sidi Salem en un mois. La retenue a été assez vite remplie en début d'hiver. Ainsi, après complet encaissement de la première crue dans le creux disponible, et déversement de la seconde crue, la troisième a poussé le déversement sur la tour et réduit la tranche de laminage résiduelle (Daoud et al., 2009).

3.2 Modélisation et reconstitution des transports solides

L'objectif principal de cette section est de développer un modèle simple qui montre l'atténuation du signal des apports solides à Jdaida avec les données disponibles sans tenir compte des paramètres physiques du bassin. La méthodologie s'intéresse d'abord au premier barrage Mellegue construit puisque, on dispose à cette station de 1954 à 2005 d'un historique des apports solides annuels au barrage MES X_Mellegue et des volumes annuels de vase évacués MES lachées Y_Mellegue. A la construction de chaque barrage correspond un stockage de MES que l'on peut dans un premier temps estimer par un pourcentage Y_n d'une quantité de MES X_n. Le coefficient annuel de réduction de transport solide pour le barrage de Mellegue qui a été calculé est de l'ordre 51 %. Le fonctionnement sera étudié après pour tout le reste du bassin versant Medjerda, grâce à des valeurs calculées pondérées par la surface de chaque barrage en appliquant le même coefficient de réduction à Mellegue pour les autres barrages, ce qui constitue bien sur une hypothèse réductrice, mais basée sur la seule observation disponible.

Avant la mise en eau du barrage Sidi Salem en 1981, les barrages Mellegue, Beni Metir, Bouherma, Kassab et Lakhemes (tableau 2) déversaient vers la mer sans autre obstacle. Depuis 1981, les quantités de MES qui ressortent de ces barrages arrivent au barrage d'ordre 2 Sidi Salem, situé juste après la station Bousalem, sur le cours de la Medjerda (tableau 2). Depuis la mise en eau du barrage Siliana en 1987, les quantités qui ressortent du barrage Lakhemes arrivent au barrage Siliana d'ordre 2. Cette étude permet de classer l'ordre de grandeur des flux des matières solides jugés assez proche de la réalité. La quantité de MES qui arrive à Jdaida est calculée en additionnant les quantités de MES lâchées par chaque barrage avec les quantités de MES produites dans le bassin intermédiaire, en appliquant à chaque fois la réduction de MES déjà calculée à Mellegue. Les résultats de cette approch, illustrés par la Fig. 3 et le tableau 2, montrent la réduction de la quantité de MES calculé à la station de Jdaida depuis la mise en eau du barrage Sidi Salem. Notons tout de même qu'il faut prendre les résultats obtenus avec beaucoup de précaution, cependant cela produit un signal d'atténuation des apports solides qui peut servir à interpréter les carottes prélevées dans la basse vallée. Nous avons eu recours à la modélisation à l'aide d'un modèle simplifié pour pallier au manque des données réelles dans le but de l'estimation du transport solide dans la Medjerda, mais il faudrait mettre en place un programme de mesures de débits solides au sortir des barrages pour valider ce modèle. Cette approche doit être approfondie en prenant en compte les informations qui pourront être obtenues auprès des gestionnaires de barrages, afin de mieux identifier la variabilité locale de la production de MES et des taux de sédimentation dans les barrages. D'autres modèles seront peut-être mieux appropriés avec beaucoup plus de précision concernant entre autre la production de sédiments. Il est recommandable de donner plus d'effort et d'investissements de la part des organismes concernés pour pouvoir donner une estimation exacte du transport solide.

Figure 3 Evolution des apports solides calculés à Jdaida et des volumes calculés et mesurés de vase évacuée à Sidi Salem.

Tableau 2 Rapport sur valeur moyenne de MES calculé au niveau de chaque barrage et des valeurs moyennes de MES calculées à Jdaida.

3.3 Approche paléohydrologique pour évaluer l'impact anthropique

3.3.1 Reconstitution de la variabilité passée des transports sédimentaires du bassin versant de la Medjerda, Tunisie

L'utilisation de données granulométriques est indispensable à l'étude des archives sédimentaires. On se limite dans cette partie au résultat d'une seule carotte C5, la plus profonde et la plus riche en couches de sédiments différentes à interpréter. Les deux autres carottes C1 et C3 ne sont pas assez profondes pour atteindre le pic de Cesium, mais elles mettent en évidence une certaine variabilité latérale des dépôts, tout en confirmant la nature essentiellement argileuse des premiers 80 cm de dépôt.

La mise en place d'un log litho stratigraphique nous a permis d'identifier visuellement les différents dépôts de crues. L'archive sédimentaire montre d'importantes variations granulométriques en fonction de la profondeur, d'un matériel silteux-argileux à sableux (D50 compris entre 5 et 240 µm). Ces niveaux sableux sont d'un grain fin à moyen (de 63 à 449 µm) aux niveaux des profondeurs 125  cm, 150, 210, 250 et 290 cm, indiquant un régime hydraulique plus énergétique. Ces pics correspondent probablement à des événements exceptionnels (crues intenses).

Figure 4 Reconstitution historique des apports sédimentaire par les traceurs chronologiques (¹³⁷Cs) et la comparaison des enregistrements de débit maximum annuel au niveau de la station hydrométrique de Jdaida.

L'ensemble des travaux d'observation, des études granulométriques et géochimiques de la carotte C5 permet d'identifier de manière plus fine le nombre de dépôts de crues. Les variations granulométriques et géochimiques sont le reflet des modifications de la dynamique sédimentaire. Le nombre de dépôts de crue est estimé à 18 et sont notés de FD1 à la base de l'archive à FD18 au sommet. Dans la séquence des dépôts de crue de la carotte C5, l'unité de dépôt FD1 présente une activité non négligeable de ¹³⁷Cs (Fig. 4) indiquant que tous les dépôts de crue de FD1 à FD18 sont tous postérieurs à 1950. D'autre part, la forte activité en cesium-137 mesurée au niveau de l'unité FD 7 (172 cm; 8,2 mBq g⁻¹), nous permet de faire la supposition que ce dépôt de crue est associé à la période de production maximale de Cesium-137 dans l'atmosphère, c'est à dire au milieu des années 1960 (autour de 1963, Fig. 4). Le ²¹⁰Pb est en excès dans l'unité FD1, ce qui indique que les dépôts de crue FD1- FD18 sont récents et tous postérieurs aux années 1900. D'autre part, nous avons trouvé sur le terrasse alluviale entre 120 et 140 cm, des déchets de pots de yaourth “ Sami ” fabriqué dans les années 1970.

Ces traceurs chronologiques, nous ont permis de mettre en place une chronologie approximative des dépôts de crue (FD1 à FD18) que nous comparons à l'historique des crues disponibles depuis 1930 à la station de Jdaida (Fig. 4). La comparaison de ces enregistrements permet d'attribuer de manière plus fine des âges aux dépôts de crues de la séquence FD. L'événement de 1969, le premier événement puissant après 1963 enregistré par les archives instrumentales (832 m³ s⁻¹), a probablement déposé une unité de texture différente en silt limite sables fins (FD8  : 44 µm). Un autre événement plus puissant encore, est enregistré en 1973 avec un débit estimé à 3500 m³ s⁻¹ à Medjez El Bab. Les mesures hydrologiques réalisées pendant la crue montrent qu'un débit de 850 m³ s⁻¹ est passé dans le lit du fleuve, mais que le reste a débordé déclenchant des inondations exceptionnelles (Claude et al., 1977). Si nous supposons qu'il existe une corrélation entre l'intensité d'un événement et la taille granulométrique qui caractérise son dépôt (Dezileau et al., 2014b), alors nous pourrions associer l'événement de 1973 à l'unité FD10 caractérisée par des sables fins (FD10 : 59 µm, 20 cm d'épaisseur). Cela est cohérent par rapport aux datations au ¹³⁷Cs et la découverte dans ces dépôts de pots de yaourt “ Sami ”. En nous fondant sur une relation possible entre la granulométrie des sédiments et l'ampleur d'un événement, FD17 (FD 7 : 100 µm) pourrait être associé à 2003, pour un débit 750 m³ s⁻¹ à Jdaida et FD 12 avec du matériel fin (D50 of 26 µm) à 1980 (Fig. 4). D'autres unités sous-jacentes FD5 avec D50 de 212 µm, FD4 avec D50 de 212 µm et FD1 avec D50 de 224 µm ont été affectées à quatre événements de plus faible ampleur que celui de 1973 (538, 432 et 344 m³ s⁻¹) qui ont eu lieu en 1959, 1957 et 1956 (ou 1952) respectivement (Fig. 4). Les couches sédimentaires associées à ces crues sont caractérisées par des sables moyens (211 à 228 µm). Ceci est expliqué par le fait qu'avant 1969, les sections du cours d'eau étaient dégagées, la terrasse alluviale de la zone d'étude était en évolution. Cette approche nous a permis de tracer l'histoire sédimentaire à travers l'étude de la carotte de sédiments prélevés. Les résultats indiquent un dépôt de 75 cm entre 1981 et 2014, alors qu'il était de 230 cm entre 1956 (ou 1952) et 1981 (120 cm entre 1956/52 et 1963, puis 105 cm entre 1963 et 1981). L'histoire sédimentaire du fleuve Medjerda montre une réduction du flux sédimentaire à l'aval du barrage de Sidi Salem du fait des multiples aménagements.

Tableau 3 Evolution des sédiments transportés dans la basse vallée de Medjerda depuis 1950.

3.3.2 Etude de l'évolution des apports  : impact anthropique

La plaine de Medjerda est l'une des zones les plus peuplées de Tunisie, 288 834 habitants en 2004 (INS). Depuis l'indépendance en 1956, le pays a connu une croissance rapide de la population, et les zones urbaines au nord ont augmenté de 7,16 % entre 1950 et 2007 (Samaali, 2011). Compte tenu de l'augmentation des utilisations de l'eau dans les pôles économiques urbains au nord-est et au centre-est, l'état a suivi un plan de développement pour la mise en place d'un grand nombre de retenues, dont 9 grands barrages construits dans le bassin versant de la Medjerda, afin de répondre à l'intérêt général.

Les résultats d'analyse granulométrique de la carotte C5 montrent principalement une succession de petites couches de matériel fin (silt ou argile) depuis la mise en eau du barrage Sidi Salem en 1981 et un taux de sédimentation faible (2,3 cm an⁻¹), alors qu'il était beaucoup plus important, environ 4.75 cm an⁻¹ entre 1963 et 1981, et environ 10 cm an⁻¹ entre 1950 et 1963 (tableau 3).

Ceci illustre bien l'impact de la mise en eau des nombreux barrages sur le bassin, qui prive la basse vallée, l'estuaire et la zone littorale de la majeure partie des apports sédimentaires, en particulier grossiers.

4 Discussion

Généralement, les barrages limitent les apports de sédiments à la mer et participent au déficit général de sédiments sur les côtes. Ils capturent jusqu'à 30 % des flux de sédiments produits sur les massifs (Vorosmarty et al., 2003). Depuis la construction des barrages sur le bassin versant de la Medjerda, dont le plus important, celui de Sidi Salem, le lit de la Medjerda a connu d'importants changements morphologiques dans sa basse vallée, et en particulier un important engraissement du lit mineur. Cet engraissement est dû en grande partie à la réduction des débits de crue exceptionnels qui devraient “ chasser ” les dépôts anciens du lit mineur. Ceci a diminué la capacité de transit du fleuve ainsi que la vitesse d'écoulement entre la sortie du barrage et l'estuaire. Il en résulte l'apparition de débordements plus fréquents et la propagation des inondations pour des débits de plus en plus faibles. Par ailleurs, les crues ont également beaucoup moins d'impact mécanique sur la région côtière.

Les résultats de mesures de l'alluvionnement menées depuis une vingtaine d'année, montrent l'importance de la sédimentation dans les retenues des barrages de Tunisie. Des concentrations supérieures à 100 g L⁻¹ ont été enregistrées dans les eaux lors des crues de la Medjerda (Rodier et al., 1981). De ce fait, les retenues créées par les différents ouvrages hydrauliques sont toutes confrontées à plus ou moins long terme à l'alluvionnement. Les mesures de transport solide effectuées en Tunisie et les études ponctuelles détaillées restent insuffisantes et n'ont pas été menées assez longtemps pour permettre des études exhaustives sur le sujet. La modélisation du transport sédimentaire proposée dans cette étude résulte donc d'une approche simplifiée destinée à montrer l'atténuation du signal des apports sédimentaires, calculé à partir de données fragmentaires de transport solide et de taux d'envasement disponibles sur quelques sites d'étude sur le bassin. Cette première approche indique une réduction croissante des transports sédimentaires avec l'augmentation du nombre de barrages, qui peut servir également à supporter l'interprétation des archives sédimentaires prélevées dans les carottes de la basse vallée.

5 Conclusion

La Medjerda est devenu un fleuve presque totalement contrôlé, et les matières solides qu'il apporte ont été réduites de façon considérable. Ce qui indique aussi que les crues ont beaucoup moins d'impact mécanique sur la région côtière. La mise en évidence de différents paléo-événements de crue qui ont été datés par des méthodes complémentaires (¹³⁷Cs, ²¹⁰Pbex, géochimie, historiques des crues) montre que la terrasse FD a enregistré l'ensemble des événements de crue entre 1952 (ou 1956) et 2014. Une succession visible de couches sédimentaires correspondant aux dépôts des différentes crues qui se sont succédées sur le site d'étude a été déterminée. Elle permet de reconstituer l'histoire des apports sédimentaires du fleuve Medjerda jusqu'à la mer sur une période de 60 ans, avant et après la construction de nombreux barrages. Cette approche a permis d'évaluer le taux de réduction et la variabilité interannuelle des transports des sédiments. Toutefois, il faudrait une étude multi-site pour valider cette approche.

D'après la carotte C5, le taux de sédimentation était d'environ 10 cm an⁻¹ entre 1952/56 et 1963, principalement constitue de sables, puis il est passé à environ 4.75 cm an⁻¹ entre 1963 et 1982, avec un équilibre plus net entre sables, silts et argiles. Après la mise en eau du barrage Sidi Salem en 1981 sur le cours d'eau principal, les apports consistaient principalement en matériel fin, avec un taux de sédimentation faible (2,3 cm an⁻¹), Il s'agit d'une diminution significative des apports sédimentaires à la mer, et d'une modification du type de sédiments, essentiellement fins avec une disparition du matériel grossier, dû aux barrages, et en particulier celui de Sidi Salem sur le fleuve Medjerda. Actuellement la surface contrôlée par les barrages existants sur le bassin de la Medjerda est de 130 % de la surface du bassin de la Medjerda, ce qui veut dire que l'eau de ruissellement est stockée 1,3 fois sur son parcours vers l'estuaire, avec les conséquences enregistrées dans la carotte sur la dynamique sédimentaire pré-estuarienne. D'autres études (Benmoussa et al., 2017), illustrent l'impact de ce déficit d'apports sédimentaires dans la morphodynamique côtière du golfe de Tunis.