Water quality and sediment transport issues in surface water

Proc. IAHS, 377, 67–76, 2018

[Link] Open Access
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.

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
Fatma Kotti1 , Laurent Dezileau3 , Gil Mahé2 , Hamadi Habaieb1 , Malik Bentkaya4 , Claudine Dieulin2 ,
and Oula Amrouni5
1 Institut National Agronomique de Tunis, 43, Av. Charles Nicolle 1082 Tunis-Mahrajène, Tunisie
2 UMR HydroSciences Montpellier/IRD, Montpellier, France
3 Géoscience Montpellier, Université de Montpellier, Montpellier, France
4 Centre National des Sciences et Technologies Nucléaires, Sidi Thabet, Tunisie
5 Institut National des Sciences et Technologies de la Mer, Salammbô, Tunisie

Correspondence: Fatma Kotti ([Link]@[Link])

Received: 15 June 2017 – Accepted: 4 October 2017 – Published: 16 April 2018

Résumé. 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 quan-
tification of sediment transport at the outlet. This study focuses on the largest watershed in Tunisia, the Wadi
Medjerda (23 600 km2 ). The main objective of this work is to assess the reduction of sediment transport follo-
wing 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 Sa-
lem, 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 sedimen-
tary 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 137 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.

Published by Copernicus Publications on behalf of the International Association of Hydrological Sciences.

68 F. Kotti et al.: Etude de l’impact des barrages sur la réduction des transports sédimentaires

1 Introduction guliers et peu turbulents, afin de récolter des informations

sur les dépôts réels sur plusieurs décennies. Nous avons uti-
Dans le contexte actuel des changements climatiques, le lisé des techniques d’analyse stratigraphique et de datation
suivi du transport sédimentaire dans les rivières est impor- des dépôts (137 Cs, 210 Pb, géochimie) pour évaluer la hauteur
tant pour mettre l’accent sur les problèmes environnemen- atteinte par les eaux lors de différents évènements d’inon-
taux souvent reliés à des changements observés dans le cycle dation, et proposer une chronologie des crues de la Med-
hydrologique : modification de la distribution des précipita- jerda. L’analyse pluridisciplinaire de ces archives sédimen-
tions, du ruissellement et de l’intensité des extrêmes hydro- taires permet de fournir de précieuses informations sur l’évo-
logiques (Bates et al., 2008 ; Hungtinton, 2006 ; Labat et lution climatique des dernières décennies, les événements ex-
al., 2004 ; Solomon, 2007). Le suivi du transport sédimen- trêmes et enfin sur les effets des actions anthropiques passées
taire permet également d’évaluer et de quantifier les trans- et récentes (Vilanova et al., 2006 ; Selby et Smith, 2007).
ferts de matières des continents vers les océans. Cependant,
le problème du nombre réduit de données relatives aux trans-
2 Matériel et méthode
ports solides a toujours freiné les investigations. Les me-
sures de transport solide effectuées en Tunisie et les études
2.1 Site d’étude du bassin versant de la Medjerda
ponctuelles détaillées restent insuffisantes pour traiter ce su-
jet (Ghorbel et Claude, 1982). Le défi consiste donc à utiliser Le fleuve Medjerda est le plus long de Tunisie, avec 484 km
des méthodes de mesures simples pour combler le manque de dont 350 km en Tunisie (Ben Mansoura et al., 2001). Il prend
connaissances. La présente étude s’intéresse au fleuve Med- sa source près de Souk-Ahras, dans le Constantinois algérien,
jerda en Tunisie, afin d’étudier en détail l’évolution des trans- puis coule vers l’est avant de se jeter dans la mer Méditerra-
ports solides à l’exutoire de son bassin versant. L’Oued Med- née dans le golfe de Tunis. La Medjerda draine 23 600 km2
jerda est devenu un fleuve presque totalement contrôlé depuis dont 7500 km2 en Algérie (Habaieb, 1992) (Fig. 1).
l’intensification des aménagements sur le bassin au début du Avant de se jeter dans la mer, la Medjerda traverse trois
XXème siècle, pour atténuer les effets des crues dans la basse types de plaines :
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 amont, la haute vallée (du côté de Ghardimaou) est
en particulier depuis l’ouverture du barrage Sidi Salem sur le constituée de bassins colmatés par des sédiments qua-
cours principal de la Medjerda, en 1981, responsable d’une ternaires. Le fleuve et ses affluents sont caractérisés par
modification majeure du régime hydrologique du fleuve en des méandres encaissés de 8 à 10 m. Ces régions ap-
aval (Zahar et al., 2008). L’alluvionnement des retenues des partiennent à des zones faiblement montagneuses. Les
aménagements constitue un problème récent et majeur (Ben fortes crues ne débordent qu’exceptionnellement en de-
Mamou et Louati, 2007). Or, l’histoire des apports sédimen- hors du lit du fleuve (Kotti et al., 2016) ;
taires dans cette zone d’étude n’est pas connue, ce qui en-
traîne deux difficultés : – en aval de la retenue de Sidi Salem, de Testour à Jdaida,
dite moyenne vallée de la Medjerda, le fleuve est en-
– comment évaluer l’apport solide à la mer si l’histoire caissé dans la plaine. L’encaissement repasse à plus de
des apports sédimentaires n’est pas connue ? 10 m en aval de Laroussia avant de s’atténuer rapide-
ment vers Jdaida. Dans tout ce secteur la plaine est inon-
– comment identifier les différents facteurs contrôlant dable surtout dans la partie centrale. En allant vers la
l’apport sédimentaire par la Medjerda ? mer les reliefs de ce bassin s’adoucissent ;
La paléohydrologie peut en partie répondre à ces questions. – dans la basse plaine, en aval de Jdaida, commence le
On tentera dans ce cas d’établir un historique de la sédimen- delta. La pente de la Medjerda est encore marquée et
tation dans le bief aval du fleuve Medjerda, avant la mer, et son tracé a fréquemment été modifié, soit naturelle-
de relier cet historique avec l’histoire des débits, l’évolution ment soit par des travaux (coupure de méandres, ou-
climatique et la construction des aménagements. La paléohy- verture de nouvel émissaire. . . ). Du point de vue géo-
drologie est la reconstruction de l’ampleur et de la fréquence logique, la Medjerda est occupée par des alluvions qua-
des inondations récentes et passées en utilisant des archives ternaires le plus souvent argilo-sableuses à argilo-sablo-
sédimentaires (Baker et al., 2002). Notre travail va se concen- limoneuses.
trer principalement sur la Medjerda, près de la station hydro-
logique de Jdaida située dans la basse vallée mais suffisam- La pluviométrie moyenne annuelle dépasse 1000 mm au nord
ment loin de la mer pour ne pas être impactée par la marée, du bassin versant de la Medjerda et diminue progressivement
zone d’étude pour laquelle il existe une quantité importante vers le sud. Il s’agit d’une pluviométrie souvent concentrée
de données historiques, particulièrement sur les débits des dans le temps, irrégulière et inégalement répartie sur les dif-
70 dernières années. Notre approche consiste à prélever des férentes saisons (Kotti et al., 2016). La base des données hy-
carottes de sédiment dans un méandre où les dépôts sont ré- drologique provient en grande partie de la Direction Générale

Proc. IAHS, 377, 67–76, 2018 [Link]/377/67/2018/

F. Kotti et al.: Etude de l’impact des barrages sur la réduction des transports sédimentaires 69

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

des Ressources en Eau. Par ailleurs le laboratoire HydroS- 2.3 Technique de prélèvements des sédiments,
ciences Montpellier (HSM) dispose d’une base de données méthodes de mesure et géochronologie
pluviométrique importante au pas d’espace du demi-degré
carré, qui propose des valeurs mensuelles de pluie par maille Le site du carottage a été choisi après plusieurs missions de
de 2500 km2 sur la période 1940–1999, complétée par des terrain (2013, 2014, 2015) qui avaient pour objectif d’iden-
données nationales (DGRE). tifier différents sites potentiellement intéressants pour la thé-
matique. Le site d’étude est localisé dans un méandre en aval
2.2 Régime hydrologique et stratégie de mobilisation du dernier barrage -Sidi Salem- mais en amont de la zone
des eaux de surface estuarienne, proche de la station hydrologique Jdaida, la der-
nière avant la mer. Trois carottes ont été échantillonnées sur
Le régime hydrologique est typiquement méditerranéen. Les cette même terrasse alluviale, C1 de 90 cm, C3 de 120 cm et
contrastes saisonniers sont très marqués puisque, pendant C5 de 300 cm de profondeur. L’approche paléohydrologique
l’été, le débit d’étiage peut se réduire à moins d’ 1 m3 s−1 tan- a été utilisée pour comparer les résultats sédimentologiques
dis que, pendant le reste de l’année, celui d’une crue de pé- et géochronologiques de cette archive sédimentaire aux in-
riodicité décennale atteint 1000–1200 m3 s−1 (Rodier et al., formations de la station hydrométrique de Jdaida située juste
1981). Le régime hydrologique de ce fleuve très hétérogène a en amont du site d’étude choisi.
été régulé par de nombreux ouvrages construits sur son cours L’identification et la mise en évidence des différents évé-
principal ou ses affluents. A partir de 1945, un programme nements de crues ont été réalisées grâce à l’inspection mi-
d’équipement de la Medjerda a été entrepris par la Direction nutieuse de chaque dépôt (décompte des lamines, identifica-
des Travaux Publics de Tunisie, destiné à mettre en valeur les tion de paléosols, détection des surfaces d’érosion, étude de
ressources encore inutilisées de cette vallée. Ce programme la variabilité granulométrique à partir d’un granulomètre la-
a permis la construction du barrage régulateur de Mellegue ser Beckman Coulter LS 13320 au laboratoire Géosciences
et la construction du barrage de Beni Metir (Strugo, 1955). Montpellier.
Depuis l’indépendance de la Tunisie, le programme de déve- Toute reconstitution paléo-environnementale, à partir de
loppement du bassin de la Medjerda a été poursuivi grâce l’étude des archives sédimentaires, nécessite l’établissement
à la réalisation de plusieurs ouvrages : Lakhmes (1966), d’une échelle chronologique la plus fine possible, afin de da-
Kasseb (1968), Bouhertma (1976), Siliana (1987), Sidi Sa- ter avec une bonne résolution les événements hydrologiques
lem (1981) et R’Mil (1999). En Algérie, un seul grand bar- passés. Les isotopes principalement utilisés pour dater les
rage a été identifié, Ain Dalia, mis en eau en 1988 et d’une événements du siècle dernier sont le 210 Pb et le 137 Cs (Dezi-
capacité de 82 Mm3 . Avec la mise en eau de ces barrages, le leau et al., 2014a). Ils permettent une datation relativement
régime hydrologique du fleuve Medjerda a été modifié (Fig. 1 précise, grâce à leurs propriétés chimiques (réactivité vis-
pour la localisation et tableau 1). à-vis des particules) et physiques (périodes de désintégra-

[Link]/377/67/2018/ Proc. IAHS, 377, 67–76, 2018

70 F. Kotti et al.: Etude de l’impact des barrages sur la réduction des transports sédimentaires

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

BARRAGES/ Superficie Capacité Envasement des Taux Taux d’envasement

Mise en eau du bassin totale initiale retenues Mesurée d’envasement des annuel des retenues
(km2 ) (Mm3 ) (Mm3 ) retenues (%) (%)
MELLEGUE (1954) 10 300 182 122 67 1,36
BENI METIR (1954) 103 62 1 ∗ ∗
KASSEB (1966) 170 81 5 6,08 0,19
SIDI SALEM (1981) 17 885 814 171 21 0,86
BOU HEURTMA (1976) 457 117 5 4,65 0,16
SILIANA (1987) 972 70 37 52,2 3,25
LAKHMESS (1968) 127 8,22 1 12,2 0,37
RMIL (1999) 221 4 ∗ ∗ ∗
AIN DALIA (1988) 193 82 (64) (22) 0,91 (ANBT)
∗ Données manquantes.

tion, rayonnement gamma). Le 137 Cs et le 210 Pb sont mesu- demi-vie (30,17 ans) proche de celle du 210 Pb, il permet de
rés par spectrométrie gamma (détecteur CANBERRA BEGe dater les séquences sédimentaires sur une même échelle tem-
3825) au laboratoire Géosciences Montpellier (France) et au porelle. Les profils d’activité du Césium 137 en fonction de
CNSTN à Sidi Thabet (Tunisie). la profondeur dans une coupe sédimentaire présentent des
pics dont le plus important est généralement attribué à l’an-
– Datation par le 210 Pbex (T1/2 = 22,3 ans)
née 1963, date du maximum de rejets atmosphériques par les
Le Plomb 210 est un isotope radioactif naturel issu de essais nucléaires. La présence de Césium 137 dans les dépôts
la chaîne de désintégration de l’Uranium 238. L’Uranium sédimentaires indiquera qu’ils ont été déposés postérieure-
238, contenu naturellement dans la croûte terrestre, subit ment aux années 1950. La mesure du 137 Cs sur le détecteur
la série de désintégration suivante : Uranium 238 – Ura- gamma CANBERRA se fait en intégrant le pic à 661 keV.
nium 234 – Thorium 230 – Radium 226 – Radon 222. Le L’erreur n’excède pas 6 %.
222Rn (T1/2 = 3,8 jours) ainsi produit, se trouve sous une
forme gazeuse. Lorsqu’il est produit à proximité de la sur-
3 Résultats
face, le gaz s’échappe dans l’atmosphère où il se transforme
rapidement en Polonium 218 qui se désintègre en quelques 3.1 Sédimentation des réservoirs
minutes en 210 Pb. Ce 210 Pb, dit atmosphérique (ou en excès)
retombe sous forme de particules humides ou sèches et s’ac- Le taux d’envasement de chaque barrage construit sur le bas-
cumule dans les sédiments. Son activité au sein de la colonne sin versant de la Medjerda met en évidence l’importance de
sédimentaire diminue avec le temps selon une loi de décrois- l’alluvionnement durant la durée de l’exploitation des rete-
sance radioactive. Le principe de la datation au 210 Pb s’ap- nues des barrages, dont Mellegue, Siliana et Sidi Salem qui
puie sur l’hypothèse selon laquelle le flux atmosphérique en sont sujettes à un alluvionnement plus important (tableau 1).
210 Pb est constant, ce qui implique que l’activité de l’isotope Il a été estimé que le volume initial de la retenue de Sidi
radioactif à la surface des sédiments reste toujours identique. Salem sera envasé jusqu’à 60 % 100 ans après la construc-
Cette hypothèse permet en mesurant l’activité du 210 Pb et tion du barrage, et l’envasement devrait être total pour les
en appliquant la loi de décroissance radioactive de détermi- retenues des deux barrages de Mellegue et Siliana dans 60 à
ner l’âge de nos sédiments. Au bout de cinq fois sa période, 65 ans (Ben Mammou et Louati, 2007). Des soutirages effi-
c’est à dire 110 ans, il a complètement disparu dans l’en- caces permettent dans ce cas de réduire le taux d’alluvionne-
vironnement (Dezileau et al., 2014b). La mesure du 210 Pb ment et d’augmenter la durée de vie de ces barrages. D’après
sur le détecteur gamma CANBERRA se fait en intégrant le Abid (1980), 59 à 64 % des apports solides dans la retenue
pic à 46,5 keV. Le 210 Pbex est estimé en soustrayant le pic du barrage de Mellegue, durant la période 1954–1980, ont
du 210 Pb au pic du 226 Ra à 186,2 keV. L’erreur n’excède pas été soutirés grâce aux manœuvres de dévasement pendant les
6 %. crues. La Fig. 2 montre le volume de vase évacué suite à des
événements exceptionnels entre 1954 et 2003 au niveau du
– Datation par le 137 Cs (T1/2 = 30,17 ans)
barrage Mellegue et entre 1982 et 2005 au niveau du barrage
Le Césium 137 est un isotope radioactif artificiel issu de Sidi Salem. En 1969, de très importants volumes de sédi-
réactions de fission nucléaire. Il est relâché dans l’atmo- ments évacués lors d’une crue se sont accumulés surtout dans
sphère lors de l’utilisation de bombes atomiques et d’ac- les zones où la pente diminue, ainsi que dans les secteurs
cidents dans des centrales nucléaires. Avec une période de évasés. En mars 1973, l’oued Medjerda a connu une autre

Proc. IAHS, 377, 67–76, 2018 [Link]/377/67/2018/

F. Kotti et al.: Etude de l’impact des barrages sur la réduction des transports sédimentaires 71

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.

crue plus intense provoquant de graves inondations avec une 3.2 Modélisation et reconstitution des transports solides
quantité exceptionnelle de matériaux transportée et ce mal-
gré les gros travaux d’endiguement effectués à Jdaida (entre L’objectif principal de cette section est de développer un mo-
1947 et 1957) (Ben Hassine and Rejeb, 2003). dèle simple qui montre l’atténuation du signal des apports so-
Depuis la mise en eau du barrage Sidi Salem en 1981 jus- lides à Jdaida avec les données disponibles sans tenir compte
qu’en 1990, le volume de vase évacué par ce barrage est des paramètres physiques du bassin. La méthodologie s’inté-
quasiment nul. On remarque le faible pourcentage de vase resse d’abord au premier barrage Mellegue construit puisque,
évacuée dû à l’arrêt des opérations de soutirage. Quelques on dispose à cette station de 1954 à 2005 d’un historique des
opérations de dévasement occasionnel ont été associées à apports solides annuels au barrage MES XMellegue et des vo-
des débits maximum exceptionnels : 595 m3 s−1 en 1990, lumes annuels de vase évacués MES lachées YMellegue . A la
671 m3 s−1 en 1995, 977 m3 s−1 en 2000 et 1020 m3 s−1 construction de chaque barrage correspond un stockage de
2003. Durant l’année 1995, les gestionnaires de Sidi Salem MES que l’on peut dans un premier temps estimer par un
ont libéré un volume total dévasé de 24,27 Mm3 . Les parti- pourcentage Yn d’une quantité de MES Xn . Le coefficient
cules solides se déposent en aval en fonction de la réduction annuel de réduction de transport solide pour le barrage de
de la vitesse du courant de turbidité, ainsi plusieurs facteurs Mellegue qui a été calculé est de l’ordre 51 %. Le fonction-
peuvent influencer les épaisseurs des dépôts accumulés tels nement sera étudié après pour tout le reste du bassin versant
que la vitesse du courant, la topographie, la présence d’obs- Medjerda, grâce à des valeurs calculées pondérées par la sur-
tacles perturbant les dépôts, etc. Ces facteurs créent une va- face de chaque barrage en appliquant le même coefficient de
riabilité latérale de dépôt importante, en aval de la retenue. réduction à Mellegue pour les autres barrages, ce qui consti-
Le lit de la Medjerda connaît de ce fait un engraissement tue bien sur une hypothèse réductrice, mais basée sur la seule
continu depuis 1981, date de la mise en service du barrage observation disponible.
Sidi Salem. De ce fait, et par manque de débits de crue im- Avant la mise en eau du barrage Sidi Salem en 1981, les
portant qui “ chassent ” les dépôts précédents, le débit occa- barrages Mellegue, Beni Metir, Bouherma, Kassab et La-
sionnant le débordement a diminué à moins de 200 m3 s−1 , khemes (tableau 2) déversaient vers la mer sans autre obs-
valeur nettement inférieure au débit de débordement avant la tacle. Depuis 1981, les quantités de MES qui ressortent de
construction du barrage qui était de 700 m3 s−1 . ces barrages arrivent au barrage d’ordre 2 Sidi Salem, si-
La Fig. 2a illustre parfaitement l’accroissement des vo- tué juste après la station Bousalem, sur le cours de la Med-
lumes solides évacués suite à la crue exceptionnelle de 2003. jerda (tableau 2). Depuis la mise en eau du barrage Siliana
Quatre fortes crues très rapprochées sont survenues dans le en 1987, les quantités qui ressortent du barrage Lakhemes
bassin de la Medjerda au début de l’année 2003, après quatre arrivent au barrage Siliana d’ordre 2. Cette étude permet de
années de sécheresse (Abdelhamid et Khalil, 2009). Elles ont classer l’ordre de grandeur des flux des matières solides ju-
cumulé un milliard de m3 d’apports à Sidi Salem en un mois. gés assez proche de la réalité. La quantité de MES qui arrive
La retenue a été assez vite remplie en début d’hiver. Ainsi, à Jdaida est calculée en additionnant les quantités de MES
après complet encaissement de la première crue dans le creux lâchées par chaque barrage avec les quantités de MES pro-
disponible, et déversement de la seconde crue, la troisième a duites dans le bassin intermédiaire, en appliquant à chaque
poussé le déversement sur la tour et réduit la tranche de la- fois la réduction de MES déjà calculée à Mellegue. Les ré-
minage résiduelle (Daoud et al., 2009). 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 sta-

[Link]/377/67/2018/ Proc. IAHS, 377, 67–76, 2018

72 F. Kotti et al.: Etude de l’impact des barrages sur la réduction des transports sédimentaires

L’archive sédimentaire montre d’importantes variations gra-

nulomé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 ex-
ceptionnels (crues intenses).
L’ensemble des travaux d’observation, des études granu-
lométriques et géochimiques de la carotte C5 permet d’iden-
tifier 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
Figure 3. Evolution des apports solides calculés à Jdaida et des de dépôts de crue est estimé à 18 et sont notés de FD1 à la
volumes calculés et mesurés de vase évacuée à Sidi Salem. 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 137 Cs (Fig. 4) indiquant que
tion de Jdaida depuis la mise en eau du barrage Sidi Salem.
tous les dépôts de crue de FD1 à FD18 sont tous postérieurs
Notons tout de même qu’il faut prendre les résultats obtenus
à 1950. D’autre part, la forte activité en cesium-137 mesu-
avec beaucoup de précaution, cependant cela produit un si-
rée au niveau de l’unité FD 7 (172 cm; 8,2 mBq g−1 ), nous
gnal d’atténuation des apports solides qui peut servir à inter-
permet de faire la supposition que ce dépôt de crue est asso-
préter les carottes prélevées dans la basse vallée. Nous avons
cié à la période de production maximale de Cesium-137 dans
eu recours à la modélisation à l’aide d’un modèle simpli-
l’atmosphère, c’est à dire au milieu des années 1960 (autour
fié pour pallier au manque des données réelles dans le but
de 1963, Fig. 4). Le 210 Pb est en excès dans l’unité FD1, ce
de l’estimation du transport solide dans la Medjerda, mais
qui indique que les dépôts de crue FD1- FD18 sont récents et
il faudrait mettre en place un programme de mesures de dé-
tous postérieurs aux années 1900. D’autre part, nous avons
bits solides au sortir des barrages pour valider ce modèle.
trouvé sur le terrasse alluviale entre 120 et 140 cm, des dé-
Cette approche doit être approfondie en prenant en compte
chets de pots de yaourth “ Sami ” fabriqué dans les années
les informations qui pourront être obtenues auprès des ges-
1970.
tionnaires de barrages, afin de mieux identifier la variabilité
Ces traceurs chronologiques, nous ont permis de mettre
locale de la production de MES et des taux de sédimentation
en place une chronologie approximative des dépôts de
dans les barrages. D’autres modèles seront peut-être mieux
crue (FD1 à FD18) que nous comparons à l’historique des
appropriés avec beaucoup plus de précision concernant entre
crues disponibles depuis 1930 à la station de Jdaida (Fig. 4).
autre la production de sédiments. Il est recommandable de
La comparaison de ces enregistrements permet d’attribuer
donner plus d’effort et d’investissements de la part des orga-
de manière plus fine des âges aux dépôts de crues de la
nismes concernés pour pouvoir donner une estimation exacte
séquence FD. L’événement de 1969, le premier événement
du transport solide.
puissant après 1963 enregistré par les archives instrumen-
tales (832 m3 s−1 ), a probablement déposé une unité de tex-
3.3 Approche paléohydrologique pour évaluer l’impact ture différente en silt limite sables fins (FD8 : 44 µm). Un
anthropique autre événement plus puissant encore, est enregistré en 1973
avec un débit estimé à 3500 m3 s−1 à Medjez El Bab. Les me-
3.3.1 Reconstitution de la variabilité passée des
sures hydrologiques réalisées pendant la crue montrent qu’un
transports sédimentaires du bassin versant de la
débit de 850 m3 s−1 est passé dans le lit du fleuve, mais que
Medjerda, Tunisie
le reste a débordé déclenchant des inondations exception-
L’utilisation de données granulométriques est indispensable nelles (Claude et al., 1977). Si nous supposons qu’il existe
à l’étude des archives sédimentaires. On se limite dans cette une corrélation entre l’intensité d’un événement et la taille
partie au résultat d’une seule carotte C5, la plus profonde et granulométrique qui caractérise son dépôt (Dezileau et al.,
la plus riche en couches de sédiments différentes à interpré- 2014b), alors nous pourrions associer l’événement de 1973 à
ter. Les deux autres carottes C1 et C3 ne sont pas assez pro- l’unité FD10 caractérisée par des sables fins (FD10 : 59 µm,
fondes pour atteindre le pic de Cesium, mais elles mettent en 20 cm d’épaisseur). Cela est cohérent par rapport aux data-
évidence une certaine variabilité latérale des dépôts, tout en tions au 137 Cs et la découverte dans ces dépôts de pots de
confirmant la nature essentiellement argileuse des premiers yaourt “ Sami ”. En nous fondant sur une relation possible
80 cm de dépôt. entre la granulométrie des sédiments et l’ampleur d’un évé-
La mise en place d’un log litho stratigraphique nous a per- nement, FD17 (FD 7 : 100 µm) pourrait être associé à 2003,
mis d’identifier visuellement les différents dépôts de crues. pour un débit 750 m3 s−1 à Jdaida et FD 12 avec du maté-

Proc. IAHS, 377, 67–76, 2018 [Link]/377/67/2018/

F. Kotti et al.: Etude de l’impact des barrages sur la réduction des transports sédimentaires 73

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

valeur moyenne de MES calculé (m3 )

Barrage Ordre de grandeur stockée dans le lâchée par le
des flux des barrage barrage
des matières solides
Mellegue Barrage d’ordre 1 2 136 732 3 269 028
Beni Metir Barrage d’ordre 1 21 359 32 677
Lakhemes Barrage d’ordre 1 28 770 42 291
Kassab Barrage d’ordre 1 38 070 56 020
Bouhertma Barrage d’ordre 1 131 856 144 417
Sidi Salem Barrage d’ordre 2 2 169 357 1 901 540
Siliana Barrage d’ordre 2 300 445 258 208
Passant à Jdaida avec les barrages 6 718 661
Passant à Jdaida sans les barrages 11 503 176

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

riel fin (D50 of 26 µm) à 1980 (Fig. 4). D’autres unités sous- fleuve Medjerda montre une réduction du flux sédimentaire
jacentes FD5 avec D50 de 212 µm, FD4 avec D50 de 212 µm à l’aval du barrage de Sidi Salem du fait des multiples amé-
et FD1 avec D50 de 224 µm ont été affectées à quatre événe- nagements.
ments de plus faible ampleur que celui de 1973 (538, 432 et
344 m3 s−1 ) qui ont eu lieu en 1959, 1957 et 1956 (ou 1952)
3.3.2 Etude de l’évolution des apports : impact
respectivement (Fig. 4). Les couches sédimentaires associées
anthropique
à ces crues sont caractérisées par des sables moyens (211 à
228 µm). Ceci est expliqué par le fait qu’avant 1969, les sec- La plaine de Medjerda est l’une des zones les plus peuplées
tions du cours d’eau étaient dégagées, la terrasse alluviale de de Tunisie, 288 834 habitants en 2004 (INS). Depuis l’indé-
la zone d’étude était en évolution. Cette approche nous a per- pendance en 1956, le pays a connu une croissance rapide de
mis de tracer l’histoire sédimentaire à travers l’étude de la la population, et les zones urbaines au nord ont augmenté de
carotte de sédiments prélevés. Les résultats indiquent un dé- 7,16 % entre 1950 et 2007 (Samaali, 2011). Compte tenu de
pôt de 75 cm entre 1981 et 2014, alors qu’il était de 230 cm l’augmentation des utilisations de l’eau dans les pôles éco-
entre 1956 (ou 1952) et 1981 (120 cm entre 1956/52 et 1963, nomiques urbains au nord-est et au centre-est, l’état a suivi
puis 105 cm entre 1963 et 1981). L’histoire sédimentaire du un plan de développement pour la mise en place d’un grand

[Link]/377/67/2018/ Proc. IAHS, 377, 67–76, 2018

74 F. Kotti et al.: Etude de l’impact des barrages sur la réduction des transports sédimentaires

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

Année 1950–1963 1963–1981 1981–2014

Nombre de grands barrages 2 5 9
Surface cumulée contrôlée km2 10 304 11 161 30 445
Taux de transport sédimentaire cm an−1 10 4,75 2,34
Nature de sédiment sable fin et silt et sable fin silt ou argile
sable grossier

nombre de retenues, dont 9 grands barrages construits dans pour permettre des études exhaustives sur le sujet. La modé-
le bassin versant de la Medjerda, afin de répondre à l’intérêt lisation du transport sédimentaire proposée dans cette étude
général. résulte donc d’une approche simplifiée destinée à montrer
Les résultats d’analyse granulométrique de la carotte C5 l’atténuation du signal des apports sédimentaires, calculé à
montrent principalement une succession de petites couches partir de données fragmentaires de transport solide et de taux
de matériel fin (silt ou argile) depuis la mise en eau du d’envasement disponibles sur quelques sites d’étude sur le
barrage Sidi Salem en 1981 et un taux de sédimentation bassin. Cette première approche indique une réduction crois-
faible (2,3 cm an−1 ), alors qu’il était beaucoup plus impor- sante des transports sédimentaires avec l’augmentation du
tant, environ 4.75 cm an−1 entre 1963 et 1981, et environ nombre de barrages, qui peut servir également à supporter
10 cm an−1 entre 1950 et 1963 (tableau 3). l’interprétation des archives sédimentaires prélevées dans les
Ceci illustre bien l’impact de la mise en eau des nombreux carottes de la basse vallée.
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. 5 Conclusion

La Medjerda est devenu un fleuve presque totalement

4 Discussion 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
Généralement, les barrages limitent les apports de sédiments beaucoup moins d’impact mécanique sur la région côtière.
à la mer et participent au déficit général de sédiments sur La mise en évidence de différents paléo-événements de crue
les côtes. Ils capturent jusqu’à 30 % des flux de sédiments qui ont été datés par des méthodes complémentaires (137 Cs,
produits sur les massifs (Vorosmarty et al., 2003). Depuis la 210 Pbex, géochimie, historiques des crues) montre que la ter-

construction des barrages sur le bassin versant de la Med- rasse FD a enregistré l’ensemble des événements de crue
jerda, dont le plus important, celui de Sidi Salem, le lit de entre 1952 (ou 1956) et 2014. Une succession visible de
la Medjerda a connu d’importants changements morpholo- couches sédimentaires correspondant aux dépôts des diffé-
giques dans sa basse vallée, et en particulier un important rentes crues qui se sont succédées sur le site d’étude a été
engraissement du lit mineur. Cet engraissement est dû en déterminée. Elle permet de reconstituer l’histoire des apports
grande partie à la réduction des débits de crue exceptionnels sédimentaires du fleuve Medjerda jusqu’à la mer sur une pé-
qui devraient “ chasser ” les dépôts anciens du lit mineur. riode de 60 ans, avant et après la construction de nombreux
Ceci a diminué la capacité de transit du fleuve ainsi que la barrages. Cette approche a permis d’évaluer le taux de ré-
vitesse d’écoulement entre la sortie du barrage et l’estuaire. duction et la variabilité interannuelle des transports des sédi-
Il en résulte l’apparition de débordements plus fréquents et la ments. Toutefois, il faudrait une étude multi-site pour valider
propagation des inondations pour des débits de plus en plus cette approche.
faibles. Par ailleurs, les crues ont également beaucoup moins D’après la carotte C5, le taux de sédimentation était d’en-
d’impact mécanique sur la région côtière. viron 10 cm an−1 entre 1952/56 et 1963, principalement
Les résultats de mesures de l’alluvionnement menées de- constitue de sables, puis il est passé à environ 4.75 cm an−1
puis une vingtaine d’année, montrent l’importance de la sé- entre 1963 et 1982, avec un équilibre plus net entre sables,
dimentation dans les retenues des barrages de Tunisie. Des silts et argiles. Après la mise en eau du barrage Sidi Salem
concentrations supérieures à 100 g L−1 ont été enregistrées en 1981 sur le cours d’eau principal, les apports consistaient
dans les eaux lors des crues de la Medjerda (Rodier et al., principalement en matériel fin, avec un taux de sédimenta-
1981). De ce fait, les retenues créées par les différents ou- tion faible (2,3 cm an−1 ), Il s’agit d’une diminution signifi-
vrages hydrauliques sont toutes confrontées à plus ou moins cative des apports sédimentaires à la mer, et d’une modifi-
long terme à l’alluvionnement. Les mesures de transport so- cation du type de sédiments, essentiellement fins avec une
lide effectuées en Tunisie et les études ponctuelles détaillées disparition du matériel grossier, dû aux barrages, et en par-
restent insuffisantes et n’ont pas été menées assez longtemps ticulier celui de Sidi Salem sur le fleuve Medjerda. Actuel-

Proc. IAHS, 377, 67–76, 2018 [Link]/377/67/2018/

F. Kotti et al.: Etude de l’impact des barrages sur la réduction des transports sédimentaires 75

lement la surface contrôlée par les barrages existants sur le Ben Mammou, A. and Louati, M. H. : Evolution temporelle de l’en-
bassin de la Medjerda est de 130 % de la surface du bassin vasement des retenues de barrages Tunisie. Revue des sciences
de la Medjerda, ce qui veut dire que l’eau de ruissellement de l’eau, Journal of water science, 20, 201–210, 2007.
est stockée 1,3 fois sur son parcours vers l’estuaire, avec les Ben Mansoura, A., Garchi, S., and Daly, H. : Analyzing forest users’
conséquences enregistrées dans la carotte sur la dynamique destructive behavior in Northern Tunisia, Land Use Policy, 18,
153–163, 2001.
sédimentaire pré-estuarienne. D’autres études (Benmoussa et
Benmoussa, T., Amrouni, O., Dezileau, L., Mahe, G., and Abdel-
al., 2017), illustrent l’impact de ce déficit d’apports sédimen- jaouad, S. : Impact of the Medjerda sedimentary fluxes on the
taires dans la morphodynamique côtière du golfe de Tunis. morphodynamic equilibrium of the northern coast of the Gulf
of Tunis (Medjerda-Raoued coast), Tunisia, PIAHS 377, IAHS
meeting, Port Elizabeth, South Africa, 2017.
Disponibilité des données. Des données sont disponibles sur Claude, J., Francillon, G., and Loyer, J. Y. : Les alluvions dépo-
demande. sées par l’Oued Medjerda lors des crues exceptionnelles de Mars
1973, Cah. ORSTOM, sér. Hydrol., vol. XIV, no. I1, 1977.
Daoud, A. and Jemmali, K. : Couplage d’une vacuateur vanne avec
Intérêts concurrents. Les auteurs déclarent qu’ils n’ont acun une tranche de laminage, cas du barrage Sidi Salem en Tunisie
conflit d’intérêts. Combining gated spillway and flood damping, Sidi Salem dam
in Tunisia. Colloque CFBR-SHF : “Dimensionnement et fonc-
tionnement des évacuateurs de crues”, 20–21 janvier 2009, Lyon
Déclaration de la numéro special. This article is part of the – Daoud Jemmali Goguel Leclerc, Couplage d’un évacuateur
special issue “Water quality and sediment transport issues in sur- vanné avec une tranche de laminage, cas du barrage de Sidi Sa-
face water”. It is a result of the IAHS Scientific Assembly 2017, lem en Tunisie, 2009.
Port Elizabeth, South Africa, 10–14 July 2017. Dezileau, L., Terrier, B., Berger, J. F., Blanchemanche,
P., Freydier, R., Bremond, L., Latapie, A., Paquier, A.,
Lang, M., and Delgado, J. L. : A multidating approach
Remerciements. Nous remercions l’Institut National applied to historical slackwater flood deposits of the
Agronomique de Tunisie et son Directeur Monsieur Mah- Gardon River, SE France, Geomorphology, 214, 56–68,
moud Elyes Hamza, pour le soutien logistique et financier, et pour [Link] 2014a.
l’octroi d’une bourse de stage à l’Université de Montpellier. Le Dezileau, L., Terrier, B., Berger, J. F., Blanchemanche, P., Latapie,
programme Campus France est particulièrement remercié pour Freydier, R., Paquier, A., Lang, M., and Delgado, J. L. : Paleohy-
l’attribution d’une bourse Eiffel de séjour long en France, ainsi que drological reconstruction of extreme floods of the Gardon River,
l’IRD pour une bourse de séjour scientifique long en France. Nous SE France, Houille Blanche, International Water Journal, 4, 44–
remercions aussi le Centre National des Sciences et Technologies 52, 2014b.
Nucléaires pour la disponibilité du matériel de carottage. Ghorbel, A. and Claude, J. : Mesure de l’envasement dans les rete-
nues de sept barrages en Tunisie : estimation des transports so-
Edited by : Christophe Cudennec lides, DRES-ORSTOM, No. : 049 34, Tunis, 1982.
Reviewed by : Didier Orange and one anonymous referee Habaieb, H. : Comparaison numérique des modèles de prévision des
crues, application à des bassins versants belge, français et tuni-
sien. Doctorat d’Etat en Sciences Agronomiques de l’Université
de Gand, Belgique, Mai 1992.
Références Huntington, T. : Evidence for intensification of the global water
cycle : Review and synthesis, J. Hydrol., 319, 83–95, 2006.
Abdelhamid, D. and Khalil, J. : Couplage d’un évacuateur vanné Kotti, F., Mahe, G., Habaieb, H., Dieulin, C., and Calvez, R. : Etude
avec une tranche de laminage, cas du barrage de Sidi Salem en des pluies et des débits sur le bassin versant de la Medjerda, Tu-
Tunisie, Colloque CFBR-SHF : “Dimensionnement et fonction- nisie, Bulletin de l’Institut Scientifique, 38, 19–28, 2016.
nement des évacuateurs de crues”, 20–21 janvier, Lyon, 2009. Labat, D., Godderis, Y., Probst, J., and Guyot, J. : Evidence for glo-
Abid, A. : Apports solides et soutirages aux barrages Nebeur sur bal runoff increase related to climate warming, Adv. Water Re-
l’oued Mellegue (1954–1980), Séminaire International d’Experts sour., 27, 631–642, 2004.
sur le Dévasement des retenues, Tunis, 1–4 juillet, 13 pp., 1980. Rodier, J., Colombani, J., and Kallel, R. : Monographie de la Med-
Baker, V. R., Webb, R. H., and House, P. K. : The scientific and so- jerda, ORSTOM, DGRE, 1981.
cietal value of paleoflood hydrology, in : Ancient Floods, Modern Samaali, H. : Etude de l’évolution de l’occupation et de l’utilisation
Hazards : Principles and Applications of Paleoflood Hydrology, du sol dans le delta de Medjerda par télédétection et systèmes et
edited by : House, P. K., Webb, R. H., Baker, V. R., Levish, D. systèmes d’informations géographiques, Thèse en vue de l’ob-
R., Water Science and Application Series, 5, 127–146, 2002. tention du Doctorat en Géographie. Université de Tunis Faculté
Bates, B. C., Kundzewicz, Z. W., Wu, S., and Palutikof, J. P. : Le des Sciences Humaines et Sociales Ecole Doctorale : Structures,
changement limatique et l’eau, document technique publié par le Systèmes et Modèles en Lettres et Sciences Humaines, 2011.
Groupe d’experts intergouvernemental sur l’évolution du climat, Selby, K. A. and Smith, D. E. : Late Devensian and Holocene sea-
Secrétariat du GIEC, Genève, 236 pp., 2008. level changes on the Isle of Skye, Scotland, UK, J. Quaternary
Ben Hassine, H. and Rejeb, M. M. : Crues et inondations dans la Sci., 22, 119–139, 2007.
moyenne Mejerda, Rapport de synthèse DGRE, 2003.

[Link]/377/67/2018/ Proc. IAHS, 377, 67–76, 2018

76 F. Kotti et al.: Etude de l’impact des barrages sur la réduction des transports sédimentaires

Solomon, S. : Climate change 2007 : the physical science basis, Vilanova, I., Prieto, A. R., and Espinosa, M. : Palaeoenvironmental
Cambridge University Press Cambridge, 2007. evolution and sea level fluctuations along the southeastern Pampa
Strugo, N. : La mise en valeur de la valeur de la vallée de la Med- grasslands coast of Argentina during the Holocene, J. Quaternary
jerda, BEST, 102, 49–95, 1955. Sci., 21, 227–242, 2006.
Vorosmarty, J. V., Meybeck, M., Fekete, B., Sharma, K., Green, P., Zahar, Y., Ghorbel, A., and Albergel, J. : Impacts of large dams on
and Syvitski, J. P. M. : Anthropogenic sediment retention : major downstream flow conditions of rivers : Aggradation and reduc-
global impact from registered river impoundments, Global Pla- tion of the Medjerda channel capacity downstream of the Sidi
net. Change, 39, 169–190, 2003. Salem dam (Tunisia), J. Hydrol., 351, 318–330, 2008.

Proc. IAHS, 377, 67–76, 2018 [Link]/377/67/2018/