Skip to navigation – Site map

HomeSectionsPaleoenvironment, Geoarchaeology,...2026Synthèse des dynamiques anciennes...

2026

Synthèse des dynamiques anciennes et analyse spatio-temporelle des modifications récentes, actuelles et futures de la ligne de rivage dans le delta de la Majerda, Tunisie

Synthesis of Past Dynamics and Spatio-Temporal Analysis of Recent, Current and Future Shoreline Changes in the Majerda Delta, Tunisia
Luc Simiroré Diatta, Katia Schörle, Semah Bettaieb, Luc Padhiou Mendy, Ameur Oueslati and Luc Lapierre

Abstracts

Since ancient times, the Majerda delta has been marked by strong morphological dynamics. Analysis of existing literature, based on stratigraphic, geoarchaeological and sedimentological data, has made it possible to reconstruct the evolution of coastal dynamics and reposition ancient coastlines since the 12th century BC. Prospective mapping of the shoreline and areas potentially subject to submersion by 2050, 2100 and 2124, based on recent rates of change and sea level rise scenarios, indicates that the Majerda delta will undergo considerable landscape changes, which will impact on socio-economic activities and ecosystems. The superimposition of past, recent and future coastline positions reveals a trend, already underway, to the sea encroaching on low-lying land, with a coastline likely to return, in places, to positions close to those of the period between the 6th and 10th centuries AD, if we consider the boundaries of potentially submersible areas, estimated at 4,380 ha ± 100 ha, 7,414 ha ± 1545 ha and 14,374 ha ± 2934 ha, for 2050, 2100 and 2124 respectively.

Top of page

Full text

L. Diatta remercie vivement Christophe Morhange, Edward Anthony et un relecteur anonyme pour la correction constructive d’une première version du manuscrit.

1. Introduction

1Les deltas sont des espaces façonnés grâce à l’interaction de conditions hydrodynamiques et climatiques favorables, ainsi qu’à des interventions anthropiques (Boudet et al., 2017). Offrant des terres arables, des ressources en eau et divers services écosystémiques (Anthony et al., 2025), ils attirent depuis longtemps des populations, accueillant entre 350 et 500 millions d’habitants à l’échelle mondiale (Anthony et al., 2024). Cependant, dans le contexte actuel du changement climatique (GIEC, 2022, p. 2913) et de l’élévation du niveau marin relatif, combinés aux pressions anthropiques croissantes et à la diminution des apports sédimentaires (Xiqing et al., 2005 ; Estournel et al., 2023), ces territoires sont de plus en plus vulnérables à l’érosion côtière et à la submersion marine.

2Le delta de la Majerda (fig. 1), formé par l’un des cours d’eau les plus importants du Maghreb et le seul à être pérenne mais aussi le plus grand de Tunisie, n’échappe pas à cette vulnérabilité (Hzami et al., 2024). L’analyse de ses transformations récentes apparaît donc essentielle dans la perspective d’une gestion durable. Ce delta ayant conservé les traces des changements survenus, présente un intérêt scientifique majeur (Reyniers, 1951 ; Lézine, 1966 ; Delile et al., 2015a ; Pleuger et al., 2019a ; Paskoff et al., 1995). Il bénéficie ainsi d’une riche base d’informations géographiques sur les tendances d’évolution de son environnement côtier, grâce non seulement à une documentation photographique aérienne, mais également des archives archéologiques (Reyniers, 1951 ; Lézine, 1966 ; Ben Jerbania, 2017), stratigraphiques, sédimentologiques (Delile et al., 2015a ; Pleuger et al., 2019a) et cartographiques (Jauzein, 1971 ; Paskoff et al., 1995). Cette base de données permet une analyse des dynamiques côtières sur le temps long. Par ailleurs, la disponibilité des données géospatiales multi-sources et multi-spectrales offre la possibilité de quantifier les dynamiques récentes et actuelles, constituant ainsi une base robuste pour l’élaboration de cartographies prospectives des évolutions futures à différentes échelles temporelles. L’élaboration de cartes prospectives est, en effet, d’une importance capitale pour la gestion des territoires côtiers, car favorisant une meilleure appréhension des dynamiques côtières (Hequette et al., 2022) et permettant d’anticiper sur les changements morphologiques futures, ainsi que les réponses éventuelles à apporter afin de mieux définir les aptitudes à l’aménagement et d’éviter l’aggravation des risques (Akacha, 2021).

3Cet article a pour objectif de reconstituer les positions anciennes, récentes et actuelles des lignes de rivage dans le delta de la Majerda depuis le XIIe siècle avant J.‑C. Il propose également une cartographie prospective des dynamiques futures du trait de côte, fondée sur les évolutions observées, afin de fournir une cartographie intégrée des dynamiques anciennes, récentes, actuelles et futures de ce système deltaïque.

2. Présentation de la zone d’étude

4Le delta de la Majerda se situe dans la partie occidentale du golfe de Tunis, au nord-est de la Tunisie (fig. 1). Il couvre une superficie d’environ 1 310 km2 et est encadré par une topographie de petites montagnes et de collines peu hautes, mais assez accidentées et bien marquées dans le paysage. Cette zone deltaïque présente une morphologie basse et d’importants bas-fonds occupés par des zones humides lacustres et marécageuses à eau douce ou de type chotts et sebkhas à eaux et sol salés. Exutoire des flux hydrosédimentaires des oueds qui prennent source dans les reliefs locaux et surtout ceux de l’oued Majerda, elle occupe l’emplacement d’une ancienne baie, le sinus Uticensis ou baie d’Utique (Paskoff, 1994 ; Chelbi et al., 1995 ; Slim et al., 2004). Le colmatage de cette baie est survenu en raison des apports terrigènes acheminés par l’oued Majerda, entrainant un recul de la mer à une distance actuelle d’environ 10 km (Pleuger, 2020). Sur ses 460 km de long environ, depuis sa source en Algérie où elle a un caractère torrentiel (Delile et al., 2015a) jusqu’à sa zone deltaïque, en passant par son cours inférieur, la Majerda traverse un bassin versant de 23 700 km2. Ce dernier a une géologie variée et ses sols (fig. 1.A) ont été exploités depuis l’époque phénicienne. L’extension des cultures et les défrichements qu’elle a nécessité ont exposé ces sols ainsi que de vastes étendues de versants à lithologies variées et souvent contrastés, à l’érosion hydrique et au transport de grands volumes de sédiments, dont est fait l’essentiel des dépôts deltaïques (Paskoff, 1994 ; Moldenhauer et al., 2008).

Figure 1

Figure 1

A : Esquisse géologique du delta de la Majerda. B : Géographie du delta. C : Rose des vents actifs au niveau du golfe de Tunis

Crédits : A : Samaali, 2011 ; B : L.S. Diatta ; C : Khsiba, 2018

5Au cours du temps, et jusqu’à la période actuelle, la Majerda a connu plusieurs défluviations (Paskoff, 1994) (fig. 1.B), transportant d’importantes quantités d’alluvions qui se sont accumulées sur de longues périodes. Son débit moyen varie selon les saisons. Lors d’une crue de périodicité décennale, il atteint 1 000 à 1 200 m3/s, et jusqu’à 3 500 m3/s (crue exceptionnelle de 1973), tandis qu’en été, le débit d’étiage peut se réduire à moins 1 m3/s (Paskoff, 1994).

6Cependant, depuis les années 1950, avec l’anthropisation croissante du bassin versant, notamment l’implantation de barrages, les apports sédimentaires vers la mer ont fortement diminué (Sliti, 1990 ; Kotti et al., 2018). Cette situation a entraîné une nouvelle dynamique du rivage, marquée par la prédominance de l’érosion marine, qui constitue le facteur dominant de l’évolution. Ceci est d’autant plus manifeste que la côte est principalement constituée de plages sableuses, dépendantes des apports sédimentaires. De fait, le long d’un linéaire littoral d’environ 40 km, la côte deltaïque présente des plages à sédiments fins (D50 ~ 200 µm) sur l’estran et du sable sous-marin sur l’avant-côte (Louati et al. 2014 ; Hzami, 2021). Ces plages, doublées par une importante flèche littorale, ne sont pas, à la différence de celles de la partie orientale du golfe de Tunis ou celles situées plus au nord, exposées aux vents dominants dans la région, et sont baignées par des eaux moins profondes, ce qui amortit l’effet des vagues. Néanmoins, cette côte est sous l’influence d’une dérive littorale monodirectionnelle qui la rend particulièrement vulnérable à certains aménagements, à l’instar des ouvrages susceptibles d’interrompre le transit sédimentaire. Cette dérive, orientée vers le sud, est générée par les vagues de directions nord-nord-est à nord-est, qui constituent les principaux agents hydrodynamiques des plages du golfe de Tunis (El Arrim, 1996). Elle prend naissance au niveau de la pointe de Rass Ettarf et sa face méridionale rocheuse caractérisée par des côtes rocheuses basses et des falaises.

3. Données et méthodes

7La méthodologie adoptée dans le cadre de cette étude repose sur une approche combinant l’analyse des dynamiques anciennes et l’évaluation des évolutions récentes et futures du delta de la Majerda. Elle s’appuie, d’une part, sur une synthèse des anciennes dynamiques, réalisée à partir de la bibliographie, proposant des résultats hypothétiques sur les processus de comblement ainsi que les positions spatio-temporelles des lignes de rivage dans le compartiment nord de l’ancienne baie d’Utique. Ce compartiment a en effet été davantage étudié que celui du sud, du fait qu’il est considéré comme le secteur probable de localisation de l’ancien port d’Utique (Jauzein, 1971 ; Paskoff, 1994 ; Chelbi et al., 1995 ; Slim et al., 2004 ; Delile et al., 2015a). D’autre part, la méthodologie repose sur une analyse diachronique de données géospatiales récentes et actuelles, mobilisées afin de quantifier les évolutions du trait de côte et d’élaborer des cartes prospectives à court, moyen et long terme.

3.1. Synthèse des anciennes dynamiques

8Plusieurs hypothèses, basées sur différentes approches, ont été avancées concernant le déplacement du lit de la Majerda et les processus d’enclavement de la ville d’Utique. Nous nous intéressons dans cet article à celles issues de l’analyse de données sédimentologiques, combinées aux données SIG (Delile et al., 2015a ; 2015b), géoarchéologiques et stratigraphiques (Pleuger et al., 2019a ; 2019b). Ces travaux permettent d’apporter des précisions sur la chronologie d’un tel déplacement et sont principalement orientés sur les processus de comblement du compartiment nord de l’ancienne baie d’Utique (voir fig. 1.B).

9Delile et al. (2015a ; 2015b) axent leur étude sur l’analyse de la progression du delta, précisément dans le couloir (voir fig. 1.B) séparant le compartiment sud et le compartiment nord de l’ancienne baie d’Utique. Ils apportent une démarche nouvelle, en prenant en compte non seulement les données géoarchéologiques et sédimentologiques, mais aussi en associant une base de système d’information géographique (SIG), constituée de cartes topographiques et géologiques (extraction des courbes de niveau et de points altimétriques, calcul des pentes), ainsi que des photographies aériennes, pour proposer des hypothèses sur les processus hydro-morpho-sédimentaires de la dynamique du rivage en direction de la mer.

10Plus récemment, Pleuger et al. (2019a ; 2019b) combinant des approches géomorphologiques, géographiques et géoarchéologiques, se basent sur l’analyse de données sédimentologiques à partir de carottes d’une profondeur comprise entre 15 et 20 m, pour reconstituer les paléo‑environnements et l’histoire de la progradation du compartiment nord du delta depuis l’Holocène.

3.2. Processus de cartographie des dynamiques futures

11Deux approches d’analyse prévisionnelle ont été appliquées pour l’élaboration des cartes prospectives d’évolution futures du trait de côte dans le delta de la Majerda. La première se base sur une projection aux échéances 2050, 2100 et 2124 des statistiques quantitatives des dynamiques historiques des 52 dernières années (1972-1985-1995-2005-2024), et en intégrant la valeur d’un recul maximal (Rmax) lors d’un événement extrême (Paquier et al., 2021). La seconde se fonde sur l’identification des terres submersibles par croisement d’une valeur d’élévation statique du niveau de la mer avec un Modèle Numérique de Surface (MNS).

3.2.1 Projection en fonction des taux d’évolution historiques et intégration d’un Rmax

12La quantification des évolutions historiques prises en compte pour la projection repose sur une mosaïque d’images satellitaires photographiques Corona de 1972, trois images Landsat 5 Thematic Mapper (TM) de 1985, 1995 et 2005, une image Landsat 8 Operational Land Imager (OLI) de 2015 et une image Pléiades à très haute résolution de 2024 (tab. 1).

Tableau 1

Tableau 1

Caractéristiques des données géospatiales utilisées.

13Le choix des années considérées s’explique tout d’abord par la disponibilité des données, mais également par la volonté de disposer d’informations couvrant des périodes intra-acquisitions relativement importantes. En outre, le choix est guidé par la volonté de mettre en évidence les effets à long terme d’événements environnementaux ayant caractérisés le delta, à l’image de la crue de la Majerda de 1973 (Claude et al., 1976), de même que la multiplication des aménagements humains qui ont influé sur la dynamique du trait de côte. La qualité des données a également été un élément déterminant dans le choix, car influençant fortement sur la précision du suivi des dynamiques par rapport aux caractéristiques du milieu et à la vitesse de variation de la ligne de rivage.

14Ces données ont nécessité des prétraitements, dans le but d’améliorer leur qualité avant leur exploitation. Il s’agit notamment de l’ortho-rectification de l’image Pléiades, du recalage de la photographie satellitaire Corona, et du ré-échantillonnage de l’ensemble des données multi‑spectrales.

15Après prétraitement des données, les lignes de références représentant le trait de côte (limite de la végétation, limite entre le sable sec et le sable humide, limite supérieure des aménagements, limite du sommet de falaise), ont été extraites sur chaque type de côte (côte de marais, côte sableuse, côte aménagée, côte de falaise). Les traits obtenus sont ensuite utilisés pour l’analyse de la cinématique du trait de côte, en calculant les taux d’évolution à partir des indices LRR (Linear Regression Rate) et EPR (End Point Rate) du module DSAS (Digital Shoreline Analysis System) développé par Thieler et al. (1999). Cette analyse a révélé une évolution variable du trait de côte dans le delta de la Majerda, mettant en évidence des secteurs en érosion (orange et rouge), en accrétion (en vert) et d’autres stables (en jaune) (fig. 2).

16L’analyse diachronique a montré que la phase de progradation était plus intense entre 1972 et 1985 (tab. 2), principalement autour de l’ancienne embouchure de la Majerda (Foum El Oued), due aux apports sédimentaires de ce dernier lors de la crue de 1973 (Claude et al., 1976). En revanche, la période plus récente (1985 à 2024) se caractérise par l’intensification d’une évolution régressive, avec une dynamique des plages perturbée par endroit par les interventions humaines le long du rivage (construction de ports, de digues et d’épis) (fig. 2).

Figure 2

Figure 2

Évolution globale du linéaire côtier du delta de la Majerda, à partir des positions du trait de côte en 1972, 1985, 1995, 2005, 2015 et 2024. Zoom sur quelques secteurs

Crédit : L. S. Diatta

Tableau 2

Tableau 2

Statistiques de la mobilité historique du trait de côte en mètre/an dans le delta de l’Oued Majerda entre 1972-1985-1995-2005-2015-2024. L’expression « Pas d’épis » renvoie aux périodes durant lesquelles les épis du port n’étaient pas encore prolongés, et ne séparaient donc pas la flèche de Foum El Oued en deux segments : nord et sud.

17Ces statistiques récentes (1985-2024) ont été considérées pour les projections futures car traduisant l’évolution globale du littoral deltaïque après la multiplication des aménagements réalisés dans le bassin versant de la Majerda, ainsi que sur le front de mer. Ces rythmes d’évolution récents sont plus appropriés pour évaluer les positions futures du trait de côte. Ainsi, pour l’élaboration des cartes prospectives à court, moyen et long terme, ces taux d’évolution diachroniques calculés ont été projetés sur 26 ans, 76 ans et 100 ans à partir de la position du trait de côte de 2024 afin d’estimer sa position en 2050, 2100 et 2124.

18La méthode de projection s’appuie sur un traitement simple de deux fichiers générés par le module DSAS lors de la quantification des dynamiques et contenant les statistiques d’évolution historiques du trait de côte. Il s’agit, d’une part, du fichier regroupant les transects qui sont perpendiculaires aux traits de côtes, suivant un intervalle de 50 m, et dans lequel sont enregistrées les statistiques d’évolution. D’autre part, le second fichier concerne l’ensemble des points d’intersections entre les transects et les différents traits de côte. Ces fichiers ont été importés dans le système d’information géographique Qgis 3.34.9 où la chaine de traitement qui va suivre a été effectuée. Les valeurs de l’indice WLR (Weighted Linear Rate) du module DSAS ont été pris en compte pour cette projection, car permettant d’inclure l’incertitude liée à la position et à l’extraction de chaque trait de côte.

19La chaine de traitement suivante est effectuée sous SIG (Qgis 3.34.9) pour la projection du trait de côte :

  • Calcul de l’orientation en degrés des transects (l’azimut).

20Cette opération est réalisée à partir du fichier contenant les transects dans la calculatrice de table, en ajoutant l’expression suivante :

Round (degrees(azimuth(start_point($geometry), end_point($geometry))),1).

21Avec « Round » qui est une fonction mathématique permettant d’arrondir une valeur au nombre de décimales indiqué (1 pour notre cas d’étude) ; « degrees » qui est une fonction mathématique qui convertit des angles en radians vers des degrés ; « azimut » qui est une fonction mathématique qui renvoie à l’azimut par rapport au nord sous forme d’angle en radians mesuré dans le sens des aiguilles d’une montre à partir de la verticale entre un point a (point de départ) et un point b (point d’arrivée) ; et « $geometry » qui est une fonction géométrique qui renvoie à la géométrie d’une entité (@Qgis 3.34.9).

  • Jointure des fichiers.

22Après le calcul de l’azimut des transects, le fichier contenant les transects a été joint à celui des points d’intersections entre les transects et les traits de côte, permettant ainsi d’attribuer à chaque point une orientation en degré de même que les valeurs des indices statistiques calculées.

  • Projection futures et restitution cartographique des valeurs projetées.

23La projection des positions futures du trait de côte s’est faite en prenant comme référence le trait de côte de 2024. La formule suivante a été appliquée pour déterminer les positions futures :

Dr = Tx * Na (Paquier et al., 2021)

24Avec Dr la distance de recul entre le trait de côte de référence de 2024 et celui de l’horizon temporel considéré, Tx le taux de recul du trait de côte et Na le nombre d’années entre la date du trait de côte de référence (2024) et l’échéance temporelle considérée.

25La restitution cartographique des valeurs projetées s’est faite à partir de l’outil « géométrie par expression » du même logiciel Qgis 3.34.9, en ajoutant l’expression suivante :

Project($geometry, "Proj2050", radians(( "azimut" + 180.0 ) % 360.0)).

26Avec « Project » qui est une fonction géométrique qui permet de projeter à partir d’un point, en utilisant une distance, un relèvement (azimut) et une élévation en radians ; « $geometry » qui est une fonction géométrique qui renvoie à la géométrie d’une entité ; « Proj2050 » qui correspond à la couche des valeurs projetées que l’on souhaite cartographier ; « Radians » qui convertit des degrés en radians, et « Azimut » qui renvoie à la couche des valeurs d’orientation en degré (azimut).

27En outre, la projection des tendances d’évolution futures du trait de côte a nécessité l’identification d’un recul majeur événementiel. L’intégration de ce paramètre dans la projection se justifie par l’augmentation de la fréquence des événements extrêmes (GIEC, 2023), pouvant induire d’importantes modifications du littoral avec des phases de recul exceptionnelles du trait de côte. Ils constituent une source d’informations prioritaires pour l’élaboration des projections futures du trait de côte (BRGM-CEREMA, 2022). Le recul événementiel correspond ainsi au recul maximum (Rmax) observé lors de la manifestation d’un phénomène extrême à l’instar des impacts de tempêtes (Winter et al., 2017).

28Ainsi, même en ayant une ampleur limitée sur la façade orientale du littoral tunisien, ces tempêtes accélèrent l’érosion et engendrent des dommages parfois irréparables, principalement au niveau des secteurs aménagés (Brahim et al., 2023). La prise en compte du taux de recul à la suite d’un événement exceptionnel est d’une importance capitale, car permettant ainsi de mieux anticiper les évolutions à moyen et long terme et de limiter les risques de sous-évaluation du périmètre pouvant être touché lors de ces événements (BRGM-CEREMA, 2022). Pour ce faire, au niveau du delta de la Majerda, la projection du recul événementiel est déterminée en fonction des données bibliographiques de la tempête du 20 au 21 janvier 1981 (Sliti, 1990). Il convient de noter qu’une tempête plus récente a affecté la zone d’étude le 19 janvier 2026 (Institut National de Météorologie INM, @Le courrier de l’Atlas), sans avoir, à notre connaissance, fait l’objet de mesures. Ainsi, même si aucune quantification par mesure directe du recul du rivage n’a été faite au lendemain des tempêtes de janvier 1981 et de 2026, une estimation peut se faire en considérant quelques travaux ayant abordés les impacts de celle 1981.

  • 1 Laboratoire Central d'Hydraulique de France.

29Les travaux du LCHF1 (1982) permettent une quantification relative du taux de recul et affirment que la manifestation de cette tempête a engendré un recul de « quelques dizaines de mètres en quelques heures ». Oueslati (2017) apporte des données plus précises en affirmant que cette tempête, qui a frappé le golf de Tunis, a entrainé l’érosion de la quasi-totalité des plages dont certaines étaient larges de plus de 15 m à la veille de l’événement.

30Cette bibliographie, constituant une source d’information du recul du rivage lors de la tempête de janvier 1981, permet ainsi d’estimer que le trait de côte aurait reculé, par endroits, d’environ 10 à 15 m après cet événement. Néanmoins, par précaution et par soucis d’un excès d’extrapolation excessive, la projection du recul événementiel (Rmax) se basera sur une valeur de 10 m, en l’absence de mesures précises au niveau de la zone d’étude. Il convient toutefois de préciser de que cette valeur pourrait être dépassée.

31La formule appliquée pour la projection future du trait de côte, avec la prise en compte du recul événementiel est la suivante :

Dr-max = Tx * Na + Rmax (Paquier et al., 2021).

32Avec Dr-max, la distance potentielle de recul lors d’un événement extrême, Tx le taux de recul historique, Na le nombre d’années entre la date du trait de côte de référence (2024) et l’échéance considérée (2050, 2100, 2124) et Rmax la valeur maximale de recul événementiel considérée.

3.2.2 Identification des terres submersibles

33Le rapport spécial du GIEC (2019) met en évidence une élévation globale de + 15 cm au cours du XXe siècle, avec 3,6 mm/an sur la période de 2006 à 2015 (GIEC, 2019, p. 323). En se basant sur un scénario climatique d’émissions élevées (Representative Concentration Pathway, RCP 8.5), le rapport de synthèse du GIEC (2023) confirme la poursuite, voire l’accélération, de cette tendance au-delà de 2100. Il projette ainsi une élévation du niveau marin d’environ 0,29 m en 2050, de 1,1 m en 2100 et de 1,88 m en 2150 (GIEC, 2023). Cette évolution progressive du niveau marin est susceptible d’entrainer une augmentation de la fréquence et de l’ampleur des événements extrêmes liés au niveau marin (GIEC, 2023), ainsi qu’une accentuation du recul du trait de côte et la submersion des terres basses.

34En tenant compte des scénarios d’élévation globale du niveau de la mer proposés par le GIEC (2023, p. 80) et les estimations régionales pour la Tunisie (CITCCNUCC, 2001 ; Hinkel et al., 2015), sur lesquelles sont basées les précédentes études portant sur le littoral tunisien (APAL-PNUD, 2015 ; Akacha, 2021 ; Hzami, 2021) des élévations statiques du plan d’eau de 0,5 m, 1 m et 2 m ont été retenues pour déterminer les zones potentiellement submersibles, respectivement aux horizons 2050, 2100 et 2124. Le choix de ces valeurs se justifie, d’une part, par leur cohérence avec les hypothèses adoptées dans les travaux antérieurs (CITCCNUCC, 2001 ; Dasgupta et al., 2009 ; Hinkel et al., 2015 ; Akacha, 2021 ; Thiéblemont et al., 2021 ; GIEC, 2021). Il est, d’autre part, guidé par la volonté de prendre en compte un scénario de forte élévation du niveau de la mer avec une augmentation qui pourrait atteindre 0,5 m d’ici 2050 (Sweet et al., 2017), 1,7 m (Thiéblemont et al., 2021) et 2 m (Dasgupta et al., 2009 ; Bamber et al., 2019) d’ici 2100, dans un scénario à fortes émissions, correspondant à une situation où les défis d’adaptation se manifesteraient beaucoup plus tôt (GIEC, 2023) avec une augmentation constante des niveaux extrêmes de la mer (Fox-Kemper et al., 2021). En effet, le long des côtes de la mer Méditerranée par exemple, Vecchio et al. (2024) soulignent que les projections du GIEC sous-estiment l’élévation future du niveau de la mer, dans la mesure où les effets de la tectonique et d’autres facteurs locaux n’ont pas été correctement pris en compte (Vecchio et al., 2024).

  • 2 DSM-OPT : Digital Surface Models from OPTical stereoscopic very-high resolution imagery.
  • 3 ForM@Ter : pôle de données et de services Terre solide.
  • 4 THEIA : pôle Thématique national des surfaces continentales.
  • 5 DINAMIS : Dispositif Institutionnel National d'Accès Mutualisé en Imagerie Satellitaire.

35La détermination des terres submersibles dans le delta de la Majerda repose ainsi sur une approche par croisement du niveau d’eau considéré avec la topographie du terrain. Un Modèle numérique de surface (MNS) à très haute résolution, créé à partir d’images stéréoscopiques Pléiades de 2024 a été utilisé. Il est obtenu via le service de calcul en ligne DSM-OPT2, opéré par Form@ter3 en collaboration avec THEIA4 et le dispositif DINAMIS5. Les MNS Pléiades étant nativement calculés par rapport à l’ellipsoïde terrestre, nous avons procédé au préalable à une correction du MNS acquis, en utilisant le géoïde terrestre de 2008. Cette correction passe par une soustraction par rapport au géoïde terrestre, afin de ramener le MNS au niveau moyen de la mer. À partir du MNS corrigé, on applique un rehaussement dont la valeur ici est de 1,5 m afin de corriger l’erreur d’altitude résiduelle sur le calcul du MNS en ajustant son altitude 0 m sur le niveau de la plage, le plus près de la mer.

3.3 Évaluation des incertitudes liées à l’ajustement des taux d’évolution considérés pour les projections, ainsi que pour les terres submersibles

36L’analyse diachronique de la dynamique du trait de côte (52 ans) dans le delta de la Majerda a révélé différentes tendances d’évolution en fonction des périodes et des sous-cellules hydrosédimentaires. Certains secteurs présentent globalement la même tendance d’évolution sur toute la période d’observation, tandis que d’autres ont enregistré au cours des périodes récentes une variation de leur dynamique. Pour le premier cas, les taux calculés sur le long terme ont été considérés (1985‑2024). En revanche, pour le second cas, l’observation d’une nouvelle tendance d’évolution a conduit à privilégier les taux d’évolution sur le court terme (2021, 2022 et 2024). Ce choix d’ajustement est opéré dans le but de prendre en compte les taux d’évolution du trait de côte les plus représentatifs des tendances d’évolution actuelles qui caractérisent chaque secteur. Toutefois, la méthode de projection présente des incertitudes variables dont l’identification et l’estimation permettent ainsi une meilleure interprétation des tendances hypothétiques d’évolution future de ce littoral deltaïque. Le tableau 3 présente une synthèse des incertitudes liées aux différentes étapes de cartographie des projections futures.

Tableau 3

Tableau 3

Identification et qualification des sources d’incertitudes lors des différentes étapes de la cartographie des positions futures du trait de côte.

37Certaines de ces incertitudes, notamment celles issues des taux d’évolution bruts calculés sous DSAS, sont quantifiables. Dans le cadre de cette étude, les principaux facteurs d’imprécision pris en compte et pouvant être estimés concernent : 1) l’erreur de géoréférencement liée au calage des images, 2) l’erreur liée à l’incertitude de la précision d’extraction du trait de côte, qui dépend à la fois de la résolution spatiale (taille du pixel) et de l’interprétation de l’information par l’algorithme, et 3) l’erreur induite par l’oscillation de la marée (tab. 4). À partir de ces erreurs, il est possible de déterminer l’erreur globale moyenne, calculée à partir des équations suivantes :

Avec Ept : erreur de position globale du trait de côte de l’année considérée ; Eg : erreur de géoréférencement de l’image ; Ee : erreur d’extraction du trait de côte ; Eo : erreur d’oscillation de marée ; E : erreur globale moyenne en m/an.

Tableau 4

Tableau 4

Marges d’erreurs estimées pour l’analyse de la cinématique du trait de côte au niveau du delta de la Majerda.
La marge d’erreur liée à la quantification des dynamiques du rivage de toute la période étudiée (1985‑2024) est de ± 0,39 m/an.

38D’autres erreurs, liées à la pérennité des ouvrages et à la complexité de l’évolution des bandes sableuses, sont difficilement quantifiables. Néanmoins, l’évaluation qualitative à partir de l’élaboration d’un organigramme synthétisant les critères d’analyse permet d’attribuer à chaque secteur un niveau d’incertitude (fig. 3). Les critères retenus se résument à la 1) révision des taux d’évolution par rapport à la tendance actuelle, 2) la présence d’une embouchure dont l’évolution est difficile à prévoir compte tenu de la mobilité résultant d’interactions complexes des processus fluviaux, marins et anthropiques, de la présence d’ouvrages dont l’influence future est difficile à prévoir (maintien ou démantèlement), de même que la mobilité complexe de cette côte deltaïque. L’erreur calculé des taux d’évolution bruts issus des statistiques du DSAS a également été prise en compte pour l’évaluation qualitative des incertitudes. En effet, chaque taux d’évolution modélisé est associé à une erreur standard, renseigné par l’indice WSE (Weighted Standard Error) du module DSAS. Pour ce faire, la moyenne de ces erreurs a été calculée pour chaque secteur de la zone d’étude et est dénommée M_WSE. Afin de définir les classes d’incertitudes liées aux taux d’évolution, les moyennes de chaque secteur ont été comparées à la moyenne des M_WSE calculée pour l’ensemble des secteurs ainsi que l’écart type associé, en utilisant une méthode de discrétisation classique basée sur la catégorisation suivante (Paquier et al., 2021) :

  • Si le M_WSE du secteur considéré est inférieur à la moyenne des M_WSE – 1 écart‑type, il est considéré que l’incertitude chiffrée est faible, le score est donc de 0 ;

  • Si le M_WSE du secteur considéré est compris entre la valeur moyenne des M_WSE – 1 écart-type et la moyenne des M_WSE, alors le score est de 1.

  • Si le M_WSE du secteur considéré est compris entre la valeur moyenne des M_WSE et la moyenne des M_WSE + 1 écart-type, alors le score est de 2 ;

  • Si le M_WSE du secteur considéré est supérieur à la moyenne des M_WSE + 1 écart type, le score est alors de 3.

39Cette méthode de discrétisation permet d’accroitre graduellement le score des secteurs dont l’incertitude liée aux taux d’évolution du trait de côte augmente par rapport aux autres secteurs, tout en ajoutant les scores de l’ensemble des critères considérés. Les scores variant ainsi de 1 (faible incertitude) à 6 (forte incertitude). Quand le score est compris entre 0 et 3, l’incertitude est considérée comme faible, de 3 à 4, elle est considérée comme moyenne et de 5 à 6, elle est considérée comme forte.

Figure 3

Figure 3

Organigramme de la méthode d’évaluation des incertitudes des taux d’évolution du trait de côte utilisés pour les projections futures

Figure adaptée de Paquier et al., 2021

40L’évaluation qualitative des incertitudes révèle que les taux d’évolutions projetées, présentent des marges d’erreur qu’il convient d’intégrer lors de l’interprétation. Malgré ces incertitudes, la régression du trait de côte dans le delta devrait se poursuivre, conformément aux tendances d’évolution diachroniques observées.

41En ce qui concerne l’estimation des terres submersibles en fonction d’une élévation du plan d’eau, les incertitudes associées aux variations du niveau marin sont difficiles à quantifier. La plupart des travaux abordant cette problématique, y compris ceux consacrés au littoral tunisien ne précisent pas les marges d’incertitude associées à leurs projections. Par ailleurs, le GIEC (2021, p. 324) indique que les écarts régionaux d’élévation du niveau marin peuvent varier d’environ ± 30 %, sous l’effet combiné de la dilatation thermique, de la dynamique océanique et de la fonte des glaces terrestres (GIEC, 2021, p. 324). Sur cette base, les marges d’incertitude des surfaces potentiellement submersibles et estimées dans la présente étude ont donc été définies en considérant 30 % de la valeur d’élévation du plan d’eau. Ce qui correspond à 0,15 m, 0,3 m et 0,6 m, respectivement pour les horizons 2050, 2100 et 2124.

4. Résultats et discussions

4.1. Hypothèses des dynamiques anciennes du delta de la Majerda : une importante mobilité des lignes de rivage

42La méthodologie utilisée par Delile et al. (2015a ; 2015b) a permis de reconstruire la morphologie du delta, qui, combinée à l’analyse des photographies aériennes (photographie de 1954), permettent de mieux identifier les processus hydro-morpho-sédimentaires responsables de la topographie du delta. Delile et al. (2015a) déduisent avec l’appui des données archéologiques, que le passage de la Majerda dans le compartiment nord s’est effectué au plus tard avant la fin de la période phénicienne, et peut-être avant la fondation de la ville d’Utique, au XIIe siècle av. J.‑C. Ils précisent que ce passage de la Majerda n’a pas dû constituer une menace d’alluvionnement de l’ancien port d’Utique, en raison de la manifestation de la dérive littorale dominante, qui contraignait les paléo‑lits à infléchir leur cours vers le sud et transportant ainsi les apports sédimentaires vers le sud. Cette situation est encore bien perceptible dans le delta avec le déplacement des sédiments du nord de la flèche de Foum El oued vers le sud, ainsi qu’au niveau de l’actuelle embouchure de la Majerda où les apports sédimentaires sont transportés puis s’accumulent au sud. Par ailleurs, les résultats de Delile et al. (2015a) révèlent l’existence d’un double système de cônes alluviaux sur l’ensemble du delta (fig. 4), dont la première composante, de type torrentiel, est présente sur l’ensemble de la frange deltaïque. La seconde, particulière aux sites de défluviation de l’oued, forme les lobes alluviaux coalescents. Les auteurs avancent l’idée que ce serait la coalescence de ces deux systèmes morphologiques à proximité d’Utique qui serait un des moteurs des processus d’ensablement, plus complexe, du port de l’ancienne ville d’Utique.

Figure 4

Figure 4

Cônes alluviaux dans le delta de la Majerda

Source : Delile et al., 2015a

43À partir de ces éléments géomorphologiques et de l’analyse des données archéologiques, ces auteurs proposent une hypothèse de comblement et des modifications de la ligne de rivage dans le compartiment nord du delta, depuis le XIIe siècle av. J.‑C. Delile et al. (2015b) affirment qu’un corridor maritime se situait à l’ouest du compartiment nord du delta. Celui-ci s’est rétréci à l’époque romaine et a complètement disparu au début du Ve siècle et plus tard au VIe siècle après J.‑C, isolant définitivement le port d’Utique avec le développement d’une tourbière au nord du promontoire (fig. 5). Cet enclavement était étroitement lié aux cônes d’alluvions qui s’épaississaient et comblaient progressivement la baie d’Utique, faisant avancer le trait de côte (fig. 5).

Figure 5

Figure 5

Carte des anciennes positions hypothétiques du trait de côte et des processus de comblement du compartiment nord du delta de la Majerda

Sources : Delile et al. 2015a ; 2015b

44Pleuger et al. (2019a ; 2019b) ont apporté plus de clarté et ont reconstruit un schéma spatio-temporel de l’évolution du paysage du compartiment nord de l’ancienne baie d’Utique (fig. 6). Ces auteurs confirment alors l’hypothèse de Delile et al. (2015a ; 2015b) faisant état de l’existence d’une bathymétrie qui dépasse, par endroits, –15 m à une certaine distance du rivage de la baie et affirment qu’au IVe siècle après J.-C, il était encore possible d’amarrer le long de la façade nord du promontoire d’Utique. Ils soulignent que l’abandon progressif de la ville antique résulte d’une augmentation des apports sédimentaires fluviaux associée à une crise érosive en amont, et non d’un changement de cours du fleuve survenu vers 2 600 av. J.‑C. Ces auteurs soutiennent qu’entre 2 600 et le Ier siècle av. J.-C. le comblement progressif du corridor entre le promontoire d’Utique et celui de Kalâat Landlouss est entamé, associé aux premiers sédiments pro-deltaïques (Pleuger et al., 2019a).

Figure 6

Figure 6

Carte du processus de comblement hypothétique du compartiment nord du delta de la Majerda

Source : Pleuger et al., 2019a

45Pleuger et al. (2019a) identifient le passage du lit de la Majerda dans le compartiment nord entre le Ier siècle av. J.-C. et le IVe siècle apr. J.-C., alors que Delile et al. (2015a) situent ce passage au XIIe av. J.-C. Cette divergence chronologique pourrait s’expliquer par les différences méthodologiques des deux études : Delile et al. (2015a) fondant leur interprétation sur l’étude des cônes alluviaux deltaïques et des données archéologiques, tandis que Pleuger et al. (2019a) s’appuient sur l’analyse de carottes sédimentaires profondes de 15 à 20 m. Néanmoins, on peut retenir de ces travaux que la défluviation de l’oued vers le compartiment nord s’est suivie d’une augmentation significative des apports sédimentaires entre le IVe et le VIIe siècle apr. J.‑C. (Pleuger et al. 2019a ; 2019b) entrainant un déplacement du trait de côte (fig. 7). Pleuger et al. (2019a ; 2019b) affirment que vers le Xe siècle de notre ère, la ville d’Utique est complètement isolée de la mer, avec un scellement complet vers le XVe siècle et un retrait de la côte à l’avantage des sédiments du fleuve vers le XVIIIe siècle apr. J.‑C. Ce n’est qu’après la grande crue de 1973 que la Majerda occupe son lit actuel, empruntant un canal creusé après la seconde Guerre Mondiale, débouchant entre Kalâat Landlouss et Raoued.

Figure 7

Figure 7

Carte des anciennes positions hypothétiques du trait de côte du compartiment nord du delta de la Majerda

Sources : Pleuger et al., 2019a ; 2019b

46Ces hypothèses complémentaires mettent en évidence une grande mobilité des anciennes lignes de rivage du delta de la Majerda, comme c’est le cas dans plusieurs deltas de la Méditerranée où certains ports ont subi un envasement rapide en raison des apports massifs de sédiments générés par les nouvelles pratiques agricoles à partir du Néolithique (Anthony et al., 2014). Ainsi, ces évolutions anciennes, combinées aux dynamiques récentes issues de l’analyse de données géospatiales, constituent une base d’informations géographiques sur les transformations survenus dans ce delta. Cette base est enrichie par l’élaboration de cartographies prospectives des dynamiques futures, permettant une analyse globale du fonctionnement de cet espace deltaïque, dans l’optique d’une gestion côtière efficace qui lui est mieux adaptée.

4.2. Essai cartographique des projections futures du trait de côte aux horizons 2050, 2100 et 2124 : une tendance à l’intensification du retrait du rivage par son érosion

47La projection du trait de côte, basée sur les dynamiques récentes et actuelles, permet d’identifier les zones potentiellement érodables. Elle repose sur l’hypothèse d’une continuation des dynamiques observées, en l’absence de nouvelles interventions ou de nouveaux ouvrages de protection. En effet, face à l’érosion côtière, la construction d’ouvrages de protection est une pratique courante le long des littoraux. Sur le linéaire côtier étudié, ces aménagements se limitent encore à quelques épis/digues au niveau des ports et de certains secteurs touristiques à l’instar de la bande de Sidi Ali El Mekki. Cependant, bien qu’ils constituent une solution technique pour la protection des installations humaines, ces aménagements contribuent à l’intensification du phénomène érosif le long du littoral tunisien (Oueslati, 2016), en favorisant le recul du trait de côte vers l’aval des ouvrages, en fonction de la direction de la dérive littorale dominante. Par ailleurs, leur efficacité nécessite un suivi et un entretien continu, impliquant parfois des moyens élevés et une forte responsabilité pour les gestionnaires. Ces derniers, peuvent dans le cas d’une mauvaise gestion ou par le manque de moyens financiers, mais également suite à la réévaluation des priorités socio-économiques, arrêter d’entretenir l’infrastructure, ce qui occasionnera sa détérioration, voire sa destruction ou son démantèlement. C’est ce qui justifie l’importance d’évaluer différentes hypothèses d’évolutions du trait de côte afin d’anticiper les scénarios possibles et de prévoir des mesures de prévention plus adaptées à une gestion côtière efficace. Pour cette étude, deux hypothèses ont été considérées : une première qui se base sur un maintien et une résistance des ouvrages de protection existants, tandis que la seconde considère la destruction de certains ouvrages, notamment des épis (à l’exception de ceux des ports) et l’ouvrage longitudinal (digue) situé en arrière-plan du port de Kalâat Landlouss.

4.2.1 Hypothèse de maintien des ouvrages existants

48La figure 8 montre la projection du trait de côte du delta de la Majerda aux horizons 2050, 2100 et 2124, en fonction des dynamiques passées. Les valeurs enregistrées (1985‑2024) ont permis de projeter un recul moyen du trait de côte à 104 m en 2050, 221 m à l’horizon 2100 et 309 m en 2124, avec respectivement des marges d’incertitude de ± 10 m, ± 30 m et ± 39 m. Les valeurs projetées et présentées dans la suite de cet article, doivent être analysées en ayant à l’esprit ces marges d’incertitude.

Figure 8

Figure 8

Carte prospective de l’évolution future du trait de côte dans le delta de la Majerda à l’horizon 2050 (trait jaune), 2100 (trait rose) et 2124 (trait rouge). Zoom sur quelques segments caractéristiques des tendances d’évolution : 1. Bande de Sidi Ali El Mekki ; 2. La brèche de la sebkha de Ghar El Melh ; 3. Ancienne embouchure de la Majerda et bande nord de la flèche de Foum El Oued ; 4. Bande sud de la flèche de Foum El Oued ; 5. Embouchure actuelle de la Majerda ; 6. Bande allant de Raoued à Gammarth

Crédit : L. S. Diatta

49Les secteurs concernés par ce recul futur sont ceux qui enregistrent déjà une régression du rivage à savoir la brèche de la lagune de Ghar El Melh, dont l’ouverture actuelle est d’environ 91 m. Cette brèche risque de s’élargir sous l’action du recul du trait de côte, projeté à 134 m en 2050, 404 m en 2100 et 565 m en 2124 (fig. 8.2). La manifestation d’un tel recul entraînerait la disparition d’une partie de la bande sableuse de la lagune de Ghar El Melh et de surcroît la disparition du paysage de la plage de Coco Beach, dont certains équipements se situent directement en bordure du rivage et sur l’eau sous forme de cabanes sur pilotis (fig. 9).

50Conjointement, la partie sud de cette bande devrait connaître une tendance à la progradation, liée au transit sédimentaire nord-sud transportant les sédiments érodés du nord du golfe vers le sud. Elle enregistrerait des taux de +28 m en 2050, +80 m en 2100 et +105 m en 2124. Par ailleurs, notons que de cette lecture des résultats, qui suggère une progradation, la réalité peut être différence et traduirait une érosion en raison de la vulnérabilité de ce linéaire face aux agents hydrodynamiques. Une situation d’érosion de la partie nord et de progradation de la partie sud de la bande sableuse pourrait conduire probablement à un déplacement progressif des installations touristiques vers le sud, qui parait protégé face à l’érosion côtière, mais non à l’abri de la submersion marine et des conséquences de l’élévation du niveau marin relatif.

Figure 9 - Illustration des enjeux de la bande sableuse de Ghar El Melh

Figure 9 - Illustration des enjeux de la bande sableuse de Ghar El Melh

1) Vue d’ensemble de la partie nord de la bande. 2) Installations touristiques en bordure du rivage. 3) Installations touristiques sur pilotis. 4) Vue d’ensemble de la partie sud de la bande de Ghar El Melh

Tous droits réservés

51Le secteur de l’ancienne embouchure de la Majerda, ainsi que la flèche de Foum El Oued vont connaitre une forte érosion, conduisant à la disparition d’une grande partie de la flèche de Foum El Oued dès 2050 (voir fig. 8.3.4). Précisons que la projection du trait de côte sur ce secteur a été réalisée en tenant compte des dynamiques observées après la construction (entre 2019 et 2020) des épis du port de Kalâat Landlouss, lesquels ont fortement influencé l’évolution récente de cette flèche. L’analyse diachronique des positions du trait de côte depuis 1972, confirme les travaux antérieurs (Paskoff, 1994 ; Oueslati et al., 2017) et montre que la flèche s’est progressivement allongée vers le sud, jusqu’à son rattachement avec la plage de Kalâat Landlouss. Cet allongement s’accompagne d’un amincissement de la flèche, dû en partie à l’absence d’apports sédimentaires depuis le nord. Cependant, depuis la mise en place des épis du port (début 2020), on constate une dynamique plus complexe avec une légère progradation sur quelques mètres au nord de l’épi, tandis que la partie sud de l’ouvrage connait une érosion plus intense.

52L’analyse des modifications récentes entre 2021, 2022 (image Landsat 8 OLI) et 2024 (image Pléiades) révèle une généralisation du recul du trait de côte sur l’ensemble de la flèche de Foum El Oued, de part et d’autre du port, à l’exception de quelques mètres au nord de l’épi septentrional, constituant la zone de séquestration des sédiments. Cette évolution négative de la quasi-totalité de la flèche s’explique, d’une part, par la baisse des apports sédimentaires, dont une partie est séquestrée par les épis mis en place au niveau de la bande de Sidi Ali El Mekki (voir fig. 8.1). D’autres part, la tendance d’évolution négative de la flèche est liée à l’intensification des agents hydrodynamiques, qui entraînent l’ouverture périodique de brèches au nord du port de Kalâat Landlouss (fig. 10) ainsi que l’élargissement de la brèche entre le port et la bande sud.

Figure 10

Figure 10

Illustration du recul actuel du segment de côte situé immédiatement au nord du port de Kalâat Landlouss : ouverture de brèches périodiques

Crédit : tous droits réservés

53Notons au passage, que ces évolutions morphologiques observées sur la bande sableuse de Kalâat Landlouss présentent des similitudes avec celles des côtes sénégalaises, à l’instar de la langue de Barbarie dans la Grande côte (Diatta, 2004 ; Sy, 2013), et les flèches sableuses de la Casamance (Thior, 2020 ; Diatta et al., 2022). Ces bandes sableuses, également allongées vers le sud, sont soumises à des processus naturels et anthropiques, entraînant l’ouverture, l’élargissement, voire la multiplication de brèches, et de ce fait un recul considérable de la ligne de rivage, comme observé sur la bande sableuse de Cabrousse, en Casamance (Diatta et al., 2022).

54Ainsi, compte tenu des évolutions observées, la dynamique récente de la flèche de Foum El Oued, enregistrant des rythmes d’évolution actuels depuis la construction des épis du port, a été considérée pour la projection du trait de côte aux horizons 2050, 2100 et 2124. Il convient de préciser que des ajustements ont été apportés au niveau de l’ancienne embouchure de la Majerda, en considérant les rythmes d’évolution observés entre 1985 et 2024. Cet ajustement se justifie par l’observation d’une dynamique régressive continue sur l’ensemble de la période d’étude (entre 1985 et 2024), y compris après la mise en place des épis du port.

55Les valeurs enregistrées permettent de projeter un recul considérable de la flèche à l’horizon 2050, susceptible d’aboutir à une connexion de la lagune de Kalâat Landlouss, formée suite à la mobilité de la flèche, avec la mer. Des phases de sédimentation ponctuelles seront probablement notées et vont favoriser une régénération de la bande au nord du port, du fait que les sédiments en provenance de la dérive littorale seront séquestrés par l’épi du port. Cependant, cette flèche restera instable et vulnérable à l’érosion, notamment en période de fortes tempêtes. En 2100, certaines zones agricoles à proximité de la lagune seront érodées suivant un recul moyen de 867 m par rapport au trait de côte de 2024 et de 357 m par rapport au trait de côte projeté en 2050.

56À l’échéance de 2124, soit les 100 prochaines années, le recul de cette flèche serait plus inquiétant avec une régression moyenne de 1 123 m par rapport à la situation actuelle, 630 m par rapport au trait de côte projeté en 2050 et 244 m par rapport au trait de côte projeté en 2100. À noter que la projection au niveau de la lagune de Kalâat Landlouss (projection de 2050 à 2100 puis de 2100 à 2124) présente une importante incertitude. En effet, la projection des positions futures du trait de côte s’est faite par rapport à la dynamique de la flèche de Foum El Oued qui ne présente pas les mêmes caractéristiques physiques ni morphologiques que le rivage de la lagune de Kalâat Landlouss. Ce dernier est plus résistant car étant constitué de sédiments plus compacts à la différence de ceux qui caractérisent la flèche sableuse, qui sont facilement érodables et sous l’influence directe du transit sédimentaire. De fait, même si le recul futur du trait de côte sur ce secteur sera important, les valeurs projetées peuvent être relativement moins élevées. Il faut dans ce cas interpréter la localisation projetée du trait de côte le long de cette lagune, comme une tendance d’évolution et considérer la surface entre le trait de côte de 2024 et les traits projetés comme correspondant à un espace au sein duquel l’érosion côtière pourrait se manifester aux horizons temporels considérés. Néanmoins, cette projection à partir des dynamiques récentes de la flèche de Foum El Oued a l’avantage de fournir un aperçu des tendances d’évolution basé sur la situation la plus extrême, permettant ainsi d’identifier les zones de ce linéaire côtier potentiellement exposées à l’érosion côtière.

57Par ailleurs, la partie sud de la plage de Kalâat Landlouss devrait quant à elle continuer à bénéficier d’une partie des apports sédimentaires de la dérive littorale en provenance du nord du golfe, auxquels s’ajouteront les sédiments érodés au niveau de la flèche de Foum El Oued et plus tard ceux de la lagune de Kalâat Landlouss qui seront transportés par la dérive littorale. Ce secteur va ainsi connaître une progradation de 50 m en 2050, de 141 m en 2100 et de 196 m en 2124 (voir fig. 8.4).

58Au niveau de l’actuelle embouchure de la Majerda, le recul du rivage qui a débuté depuis la fin des années 1970 et le début des années 1980, avec le grignotage des sédiments déposés lors de la crue de 1973, va se poursuivre pour atteindre 152 m en 2050, 426 m en 2100 et 561 m en 2124 (voir fig. 8.5). Ce recul entraînerait en 2100 le grignotage des parcelles agricoles et en 2124 les dangers porteront sur la vulnérabilité de quelques habitations éparses localisées en arrière-plan. En revanche, le sud immédiat de cette embouchure continuerait à enregistrer une accumulation de sédiment, comme c’est le cas actuellement d’après l’analyse de l’évolution entre 1985 et 2024. Ainsi, le trait de côte va enregistrer une progradation de 67 m à l’horizon 2050, de 131 m en 2100 et de 194 m en 2124.

59La bande allant de Raoued à Gammarth, dont les premiers signes de recul ont été notés dans les années 1980 (Sliti, 1990), devrait connaitre un recul plus modéré que les autres secteurs. Le recul moyen y atteindrait 10 m, 27 m et de 45 m, respectivement pour les échéances 2050, 2100 et 2124 (voir fig. 8.6). Néanmoins, un tel recul pourrait entraîner la disparition des plages de ce secteur et plus tard l’annexion de certaines installations touristiques ainsi que les habitations dès 2100.

4.2.2 Hypothèse de démantèlement/effacement de certains ouvrages

60La figure 11 montre la projection future du trait de côte aux horizons 2050, 2100 et 2124 avec l’hypothèse d’une destruction de certains ouvrages de protection, à l’exception des épis des ports de Ghar El Melh et de Kalâat Landlouss, maintenus en raison de leur importance socio‑économique. Cette hypothèse correspond à un scénario sécuritaire qui se base sur une approche qui maximise le recul du trait de côte sur les secteurs disposant d’ouvrages et affectés par l’érosion en aval et parfois en amont à une certaine distance de ces aménagements. Elle a été appliquée en reconduisant les taux de recul observés soit sur les secteurs adjacents (zone d’épis), soit sur la bande située en aval (flèche de Foum El Oued). En effet, comme mentionné plus haut, les ouvrages de protection peuvent se détériorer et ne plus jouer de façon efficace leur rôle de protection, dans le cas d’une mauvaise gestion ou par le manque de moyens financiers, mais également suite à la réévaluation des priorités socio-économiques. L’évaluation de ce scénario de destruction/démantèlement de certains ouvrages existants, permet d’anticiper les éventuels risques liés à une défaillance de ces aménagements. Les traits de côte correspondant à cette hypothèse sont matérialisés par des pointillés. La destruction de ces ouvrages entraînerait un recul plus important du rivage, à l’exception des zones caractérisées par une évolution positive « naturelle », favorisée par les apports sédimentaires de la dérive littorale nord-sud. Ces zones vont poursuivre leur progradation comme l’indique la projection de l’évolution future du trait de côte en fonction de l’hypothèse de maintien des ouvrages existants. Quant aux autres secteurs disposant d’ouvrages, à savoir la bande de Sidi Ali El Mekki ainsi que la flèche de Foum El Oued, la détérioration des ouvrages engendrerait une régression plus intense sur chacun de leurs profils côtiers.

Figure 11

Figure 11

Carte prospective de l’évolution future du trait de côte dans le delta de la Majerda dans l’hypothèse d’un effacement des ouvrages (épis) à l’horizon 2050 (trait jaune), 2100 (trait rose) et 2124 (trait rouge). Zoom sur les secteurs disposant d’ouvrages : 1) Bande de Sidi Ali El Mekki ; 2) Ancienne embouchure de la Majerda et flèche de Foum El Oued

Crédit : L. S Diatta

61Sur la bande de Sidi Ali El Mekki, l’effacement des deux épis qui se trouvent au nord du port va se traduire par la manifestation de l’érosion côtière sur l’ensemble du linéaire de la bande littorale (fig. 11.1). Une telle situation de recul du rivage entrainerait un grignotage progressif de la plage, soit en 2050 une réduction considérable de sa largeur qui mesure de nos jours entre 250 et 300 m environ. Dans les 100 prochaines années, cette bande va probablement disparaître et les conséquences directes porteront sur la destruction des biens, des installations touristiques et habitations qui seront annexés par la mer, mais aussi sur l’écosystème lagunaire.

62Concernant la flèche de Foum El Oued, on considère que la tendance actuelle de recul généralisée du trait de côte va se poursuivre et que l’ouvrage longitudinal (digue) qui sépare en partie le port et la lagune ne jouera plus son rôle (fig. 11.2). L’inefficacité de cet ouvrage provoquerait un recul de part et d’autre du port, affectant une surface plus étendue et de surcroit plus de parcelles agricoles et d’installations humaines.

63Ainsi, à l’horizon 2050, le trait de côte suivrait à ce niveau la limite côté terre de la lagune, de part et d’autre du port, constituant le point de sapement des vagues. En 2100, les terres humides et les parcelles agricoles en arrière-plan de la lagune seront annexées par la mer et en 2124, les dégâts porteront sur les installations humaines qui se situent à environ 1,3 km de la lagune de Kalâat Landlouss et dont certains risqueront d’être submergées.

64Globalement, l’analyse des dynamiques futures du trait de côte dans le delta de la Majerda, à partir des taux d’évolution passés, montre que, quelle que soit l’hypothèse qui se présentera, la mobilité des lignes de rivage va se traduire par la genèse d’un nouveau paysage et une forte érosion sur certaines zones. Cette situation est susceptible d’imposer à moyen et long termes une relocalisation des populations, suite à la destruction de leurs biens, des installations touristiques et autres activités socio‑économiques.

4.3 Cartographie des terres submersibles dans le delta de la Majerda

65L’analyse des dynamiques du trait de côte du delta de la Majerda au cours des dernières décennies révèle une fragilité croissante de son système côtier, avec une tendance à la multiplication des secteurs exposés à l’érosion côtière. De fait, face à la perspective d’une variation positive du niveau de la mer dans les prochaines années (GIEC, 2023, p. 80), cette zone deltaïque présente une grande vulnérabilité.

66La question de la vulnérabilité du front deltaïque à une élévation du niveau de la mer a ainsi été soulevée par Oueslati et al. (2015). Dans la même perspective, la fig. 12 présente une cartographie des terres potentiellement submersibles aux échéances 2050, 2100 et 2124. Ainsi en 2050, suivant l’hypothèse d’une élévation statique de 0,5 m, la bande côtière deltaïque présente une forte vulnérabilité sur une distance moyenne d’environ 3 km par rapport au trait de côte actuel. Les terres potentiellement submersibles à l’horizon 2050 sont estimées à 4 380 ha ± 100 ha. Les espaces à proximité des zones humides ainsi que les embouchures, sont les plus concernés. Une grande partie de ces territoires, correspond déjà à des secteurs fréquemment inondés avec le niveau marin actuel (Oueslati et al., 2015).

Figure 12

Figure 12

Cartographie globale des terres submersibles dans le delta de la Majerda aux échéances 2050, 2100 et 2124

Crédit : L. S. Diatta

67Autour de la lagune de Sidi Ali El Mekki, les terres submersibles en 2050 se situent sur la bande séparant cette lagune de celle de Ghar El Melh. Par contre, la bande est séparant la lagune de la mer, y compris le linéaire côtier de Rass Ettarf, est faiblement concernée par le phénomène de submersion à l’horizon 2050. Cette situation s’explique par la séquestration des sédiments charriés par la dérive littorale dominante depuis la construction du port de Ghar El Melh (1974) ainsi que les deux épis au nord du port (1993). En effet, l’accumulation des matériaux s’est traduite par un élargissement de la bande littorale, soit une progradation estimée à 70 m entre 1978 et 1991 (El Arrim, 1996), favorisant le colmatage des anciennes brèches ouvertes sur cette bande (Oueslati et al., 2015). Il s’en est suivi alors une occupation humaine, matérialisée par la création de petits polders artificiels et surélevés pour des besoins agricoles (Oueslati et al., 2015), ainsi que par le développement d’habitations.

68À l’inverse, la lagune de Ghar El Melh dont l’actuelle ouverture sur la mer résulte de la construction du port et dont la largeur risque d’augmenter au vu des dynamiques récentes du trait de côte, présente une part importante des terres qui seront potentiellement submersibles dès 2050. Cette submersion concernerait toute la longueur du lido séparant la lagune de la mer, ainsi que la flèche de Foum El Oued et l’ancienne embouchure de la Majerda. Une partie des terres basses se trouvant à proximité de la limite côté terre de la lagune de Ghar El Melh, ainsi que la lagune de Kalâat Landlouss seront submergées dans les 26 prochaines années. La submersion de ces zones lagunaires et de surcroit leur connexion avec la mer, impliquerait des conséquences négatives avec la perturbation de l’écosystème lagunaire, susceptible d’entrainer le passage d’un système lagunaire à celui marin, affectant non seulement l’activité agricole mais également la biodiversité actuelle du milieu.

69Entre l’actuelle embouchure de la Majerda et la ville de Raoued, la submersion intéresserait des surfaces plus étendues, incluant la Tunis Bay Golf Club, ainsi que la vaste étendue de terres basses de Garâat Ben Ammar (fig. 13).

Figure 13

Figure 13

Panorama des terres basses entre l’actuelle embouchure de la Majerda et la ville de Raoued

Crédit : L. S. Diatta

70Au sud du delta, les zones submersibles à l’horizon 2050 sont principalement localisées autour de la sebkha d’Ariana, dont une partie des terres est de nos jours périodiquement inondée. Cette sebkha qui constituait à l’origine une baie peu profonde, s’est comblée grâce aux apports d’alluvions (Babbou et al., 2006) et a connu une évolution remarquable du bâti depuis les années 1980 (Samaali, 2011). Les conséquences à ce niveau porteront sur l’extension des installations humaines occupant de part et d’autre les contours de la sebkha et dont certaines constructions se trouvent à proximité immédiate du territoire d’influence de la sebkha, sur des terres à très basses altitudes constamment gorgées d’eau (Chouari, 2016) (fig. 14).

Figure 14 - Exemple de terres submersibles de basses altitudes à Raoued

Figure 14 - Exemple de terres submersibles de basses altitudes à Raoued

A) terres basses urbanisées. B) Lotissements dans la sebkha

Crédit : Chouari, 2016

71À l’horizon 2100, les terres potentiellement submersibles couvriraient une superficie de 7 414 ha ± 1 545 ha et sont tributaires à un prolongement des terres submersibles en 2050. Cette estimation confirme les travaux de Oueslati et al. (2015) qui identifiaient, pour un scénario d’élévation de 0,5 m à l’horizon 2100, 6 550 ha submersibles dans le delta de la Majerda.

72Ainsi, la bande sableuse séparant la lagune de Sidi Ali El Mekki de la mer, qui présente une certaine protection à l’échéance 2050, va connaitre dans les prochaines 76 ans une très forte vulnérabilité face à la submersion marine. Cette situation engendrerait probablement une ouverture totale soit une connexion des deux lagunes à la mer. Le trait de côte suivrait la limite actuelle des habitations du village d’Aousaja, tandis que le nord des deux lagunes, limité par la ville de Ghar El Melh, serait protégé de la submersion, du fait des altitudes plus hautes au contact du pied méridional de Jbel Ennadhour.

73Entre l’ancienne embouchure et le sud de l’actuelle embouchure de la Majerda, la superficie des zones submersibles va croitre pour atteindre plus de parcelles agricoles et des installations touristiques à l’instar de la Tunis Bay Golf Club, située à proximité des terres basses de Garâat Ben Ammar.

74La sebkha d’Ariana connaitrait une évolution moins prononcée de ses terres submersibles, qui seront légèrement supérieures à celles de l’échéance 2050. Ceci est due aux reliefs environnants, à l’image des Jbels Ennahli et Sidi Amor au sud-ouest et des collines de Gammarth, de la Marsa et de Sidi Bou Saïd du côté oriental, constituant une barrière contre la submersion marine. Toutefois, la frange littorale entre Raoued et le nord de Gammarth, qui sépare la sebkha de la mer, serait concernée par la submersion, impactant fortement les activités socio-économiques de la ville de Raoued. Au sud-ouest de la sebkha, les dégâts porteront également sur les habitations longeant la limite de la zone d’influence de la sebkha.

75À l’horizon 2124, la superficie des terres susceptibles d’être submergées atteindraient 14 374 ha ± 2 934 ha du territoire deltaïque. Les lagunes de Ghar El Melh et de Sidi Ali El Mekki seraient totalement connectées à la mer. Le trait de côte suivrait alors le linéaire du village de Ghar El Melh, la limite sud du village d’Aousaja et la zone comblée au nord du promontoire d’Utique, identifiée par Delile et al. (2015a) et Pleuger et al. (2019a) comme la localisation de l’ancien port d’Utique.

76À l’est du delta, la limite des zones submersibles correspondra aux contours du promontoire de Kalâat Landlouss. Le trait de côte retrouverait approximativement l’une de ses anciennes positions, celle qu’il avait occupée entre le VIe et le Xe siècle apr. J.‑C (Pleuger et al., 2019b ; Delile et al., 2015a). Cette reprise de territoires par la mer engendrerait des conséquences néfastes au vue de l’intérêt économiques actuel du delta de la Majerda, qui constitue un espace d’importance capitale pour l’activité agricole (Samaali, 2011) et le tourisme.

77Au niveau de l’actuelle embouchure, les terres submersibles s’étendraient jusqu’au sud du promontoire de Kalâat Landlouss. Elles concerneraient la totalité de la Tunis Bay Golf Club, ainsi que certains villages longeant la partie est de la route de Sidi Amor‑Boukhitiwa.

78Autour de la sebkha d’Ariana, les zones submersibles vont moyennement évoluées vers le sud ainsi que vers la zone de Gammarth, du fait de la présence d’une topographie plus élevée. Néanmoins, certaines terres basses à l’instar de la zone de Soukra, seront concernées par cette submersion. Vers le nord de la sebkha, les zones submersibles pourraient se connecter à celles situées au sud de l’actuelle embouchure de la Majerda.

79Ainsi, quelle que soit l’échéance considérée, le delta de la Majerda présente une sensibilité importante à la variation relative du niveau de la mer. Dès 2050, l’élévation de cette dernière occasionnerait des modifications remarquables du paysage côtier.

4.4. Essai cartographique des dynamiques anciennes, récentes et futures du littoral

80Le littoral deltaïque de la Majerda se particularise par des dynamiques variables, marquées ces dernières années par une intensification du phénomène érosif. Les impacts sont perceptibles avec l’anthropisation croissante du delta, notamment depuis le début du XXe siècle, contribuant à une augmentation de la vulnérabilité côtière. Ainsi, même étant conscient de cette vulnérabilité des installations humaines en bordure de mer, le littoral reste toujours un espace fortement attractif, où la population et les surfaces artificialisées ne cessent de croître. Cette anthropisation du système côtier devrait s’accentuer au cours des prochaines années, avec une occupation humaine qui prend de plus en plus d’ampleur (Samaali, 2011), qu’il s’agisse de l’habitat, des activités touristiques en bordure de mer ou encore pour des besoins agricoles. La multiplication des enjeux justifie la nécessité d’élaborer des cartographies prospectives des dynamiques futures, basées sur les évolutions passées, afin d’identifier les territoires susceptibles d’être exposés à court, moyen et long terme en fonction des positions futures du trait de côte. Ces cartographies prospectives constituent ainsi une base d’informations géographique, visant à mieux penser et orienter la gestion et l’aménagement du littoral.

81La démarche adoptée repose sur l’intégration dans un système d’information géographique, des dynamiques historiques et récentes du trait de côte. Les positions anciennes sont issues des travaux de Delile et al. (2015a ; 2015b) et de Pleuger et al. (2019a ; 2019b), ayant permis de reconstituer l’évolution du littoral depuis le 12e siècle. Les positions récentes ont, quant à elles, été obtenues grâce à l’analyse de données géospatiales, permettant de cartographier les positions diachroniques du trait de côte pour les années 1972, 1985, 1995, 2005, 2015 et 2024. À partir d’une quantification des taux de variations de ces positions diachroniques du trait de côte, auxquels d’autres positions intermédiaires (2021 et 2022) ont été ajouté sur certains tronçons, nous avons proposé des hypothèses d’évolution futures du trait de côte dans le delta de la Majerda. Ces hypothèses traduisent les tendances d’évolution probables dans les prochains 26, 76 et 100 ans. La figure 14 synthétise les positions anciennes, récentes et futures du trait de côte dans ce delta.

Figure 14

Figure 14

Reconstruction des positions anciennes, actuelles et futures du trait de côte dans le delta de la Majerda

Crédit : L. S. Diatta

82L’analyse de la figure 14 révèle une progradation ancienne du delta de la Majerda, illustrée par les anciens traits de côte hypothétiques (traits en vert). Cette progradation a été favorisée par la présence de caractéristiques différentes de celles actuelles, notamment l’importance des apports sédimentaires. À cela s’ajoute le fait que le niveau de la mer était plus bas que celui actuel qui se traduit une élévation continue (GIEC, 2021).

83Les travaux de Sliti (1990) confirment cette progradation du rivage deltaïque jusqu’aux années 1950, période pendant laquelle on note une forte sensibilité du système côtier à l’érosion côtière. L’analyse des données géospatiales mobilisées confirme également cette progradation globale du delta jusqu’à la fin des années 1970 et début des années 1980, bien que les taux observés durant cette période soient moins remarquables, comparés à ceux de la période antique.

84Ainsi, avec l’anthropisation croissante du bassin de la Majerda, combiné à l’élévation actuelle du niveau de la mer, la dynamique se caractérise par une tendance au recul du trait de côte, même si l’on note par endroits des phases de progradation. Ces taux d’évolution des dernières années (contraste de gris), ont été considérés pour la projection des positions futures du trait de côte aux horizons 2050, 2100 et 2124 (traits en rouge).

85La comparaison des positions anciennes (traits verts), actuelles (traits gris) et futures (traits rouges) du trait de côte suggère une recolonisation progressive des terres par la mer. Cette recolonisation des terres déjà entamée, risque de s’intensifier dans les prochaines décennies, si l’on considère les variations futures du niveau de la mer.

86La figure 14 montre qu’en 2124, le trait de côte se trouverait localement à une position proche de celle occupée entre le VIe et Xe siècle apr. J.-C (Delile et al., 2015a ; Pleuger et al., 2019a), si l’on considère la limite des terres submersibles. Il suivrait la limite actuelle du linéaire côtier du promontoire de Rass Ettarf, la limite nord de la lagune de Sidi Ali El Mekki et celle de Ghar El Melh, annexant quelques terres basses de la ville de Ghar El Melh. Vers l’intérieur du delta, le trait de côte suivrait la limite sud du village d’Aousaja, puis le nord de l’ancien couloir du promontoire d’Utique ou se localise probablement l’ancien port d’Utique (Delile et al., 2015b ; Pleuger et al., 2019b). De là, il formerait une courbe concave en direction de Ghar El Melh, suivant approximativement la limite des anciens secteurs comblés par l’oued Majerda (Delile et al., 2015a ; Pleuger et al., 2019a). Il contournerait ensuite le secteur sud du promontoire de Kalâat Landlouss, avant de suivre au sud de ce dernier, une position parallèle au rivage actuel, suivant la limite côté terres des zones humides qui caractérisent le delta de la Majerda, à savoir Garâat Ben Ammar au centre et la sebkha Ariana au sud.

Conclusion

87Le delta de la Majerda, à l’instar de nombreux deltas méditerranéens, a connu une évolution remarquable depuis l’Antiquité. Bénéficiant d’une base d’information géographique sur les tendances d’évolution de son environnement côtier, l’exploitation de la bibliographie existante a permis de reconstituer les dynamiques côtières anciennes et de repositionner les anciennes lignes de rivage. Cette reconstitution met en évidence une tendance générale à la progradation deltaïque, avec la genèse d’un nouveau paysage qui a favorisé l’installation humaine et le développement des activités socio‑économiques.

88Cependant, au XXe siècle, avec l’intensification des aménagements humains dans le bassin versant de la Majerda, combinée aux changements climatiques actuels, le processus de sédimentation est fortement réduit, engendrant de ce fait le grignotage progressif de la frange littorale par la mer. L’analyse diachronique de données géospatiales multi-sources et multi-spectrales montre, ces dernières décennies, une intensification du phénomène érosif. Dans ce contexte de dynamique régressive, qui impacte les activités humaines et les écosystèmes, il nous est paru pertinent et utile d’élaborer des cartographies prospectives représentant les positions futures du trait de côte ainsi que les zones potentiellement submersibles aux horizons 2050, 2100 et 2124. Ces projections, fondées sur les taux d’évolution récents et les scénarios l’élévation future du niveau de la mer, sans oublier l’histoire plus ancienne, montrent que, quelle que soit l’échéance considérée, le delta de la Majerda connaitra des modifications considérables de son paysage, qui affecteront indéniablement les activités socioéconomiques et l’environnement. Il convient toutefois de souligner que des événements naturels exceptionnels, à l’instar de la récente tempête du 19 janvier 2026, peuvent accélérer cette évolution régressive du trait de côte. La superposition des positions anciennes, récentes et futures du trait de côte, met ainsi en évidence une tendance, déjà entamée, à la recolonisation marine des terres basses. Ainsi, la cartographie prospective proposée constitue de ce fait un outil essentiel pour la gestion de ce littoral deltaïque. Elle permet d’anticiper sur les changements morphologiques futures dans la frange marine, ainsi que sur les réponses éventuelles à apporter afin d’éviter l’aggravation des risques dans l’optique d’une meilleure définition des aptitudes à l’aménagement et à une gestion efficace de ce territoire deltaïque.

Top of page

Bibliography

Akacha, Z., (2021), L'île de Jerba (Tunisie orientale) : indicateurs et facteurs de l'évolution récente de la position de rivage et vulnérabilité à une élévation du niveau marin. Thèse de doctorat en géographie, Université de Tunis, Tunisie, 279 p.

Anthony, E., Syvitski, J., Cohen, K. M., Saito, Y., Zainescu, F., Vespremeanu-Stroe A., Nicholls R. J. , Marriner N., Amorosi A., Maselli V., Minderhoud P. S. J., Tamura T., Day J., Woodroffe C. D., Preoteasa L., Tatui F., Sabatier F., Morhange C., Besset M., , P. Kemp P., Chen Z., (2025), A 7000-Year Record of Human Influence on Global River Deltas: Geomorphology, Stratigraphy, Sediment Dynamics, and the Anthropocene. Earth Science Reviews 271, 35 p. https://doi.org/10.1016/j.earscirev.2025.105302

Anthony, E., Syvitski, J., Zainescu, F., Nicholls, R. J., Cohen, K. M., Marriner N., Saito, Y., Day, J., Minderhoud, P. S. J., Amorosi, A., Chen Z., Morhange, C., Tamura, T., Vespremeanu-Stroe, A., Besset, M., Sbatier, F., Kaniewski, D., Maselli, V., (2024), Delta Sustainability from the Holocene to the Anthropocene and Envisioning the Future. Nature Sustainability 7, p. 1235‑1246. https://doi.org/10.1038/s41893-024-01426-3

Anthony, E.J., Marriner, N., Morhange, C., (2014), Human Influence and the Changing Geomorphology of Mediterranean Deltas and Coasts over the Last 6000 years: from Progradation to Destruction Phase ? Earth‑Science Reviews 139, p. 336–361. https://doi.org/10.1016/j.earscirev.2014.10.003

Babbou, A. C., Darragi, F., Addad, A., (2006), Étude minéralogique et géochimique d’une carotte de la sebkha d’Ariana. Revue Méditerranéenne de l’environnement, p. 160-169.

Bamber, J.L., Oppenheimer, M., Kopp, R.E., Aspinall, W.P., Cooke, R.M., (2019), Ice Sheet Contributions to Future Sea-Level Rise from Structured Expert Judgment. PNAS 116, p. 11195‑11200. https://doi.org/10.1073/pnas.1817205116

Ben Jerbania, I., (2017), La production des amphores ovoïdes de type « Africaines » à Utique. Antiquités africaines 53, p. 175‑192. https://doi.org/10.4000/antafr.303

Brahim, F., Bada, D., (2023), Aménagement et érosion accélérée des plages sur la côte orientale de la Tunisie : exemples au Cap Bon et au Sahel, Dynamiques Environnementales 52, p. 1‑43. https://doi.org/10.4000/12arv

Boudet, L., Sabatier, F., Radakovitch, O., (2017), Modeling of Sediment Transport Pattern in the Mouth of the Rhone Delta : Role of Storm and Flood Events, Estuarine. Coastal and Shelf Science 198, p. 568-582. https://doi.org/10.1016/j.ecss.2016.10.004

BRGM/Cerema, (2022), Recommandations pour l’élaboration de la carte locale d’exposition au recul du trait de côte, Co-édition BRGM et Cerema, 95 p.

Chelbi, F., Paskoff, R., Trousset, P., (1995), La baie d’Utique et son évolution depuis l’Antiquité : une réévaluation géoarchéologique. Antiquité africaines 31, p. 7‑51

Chouari, W. (2016) Spatialisation et Analyse du Risque d’Inondation dans le Bassin versant de la Sebkha de l’Ariana (Tunisie Nord-orientale). Bulletin de la Société de Géographie d’Égypte, 89, p. 123‑154. https://doi.org/10.21608/bsge.2016.90344

CITCCNUCC, (2001), Communication initiale de la Tunisie à la convention cadre des Nations Unies sur les changements climatiques. Tunis, ministère de l’Environnement et de l’Aménagement du territoire, Tunisie, 211 p.

Claude, J., Francillon, G., Loyer, J.Y., (1976), Les alluvions déposées par l'oued Medjerda lors des crues exceptionnelles de mars 1973, ministère de l’Agriculture. Convention DRES-ORSTOM action de type A, 179 p. https://horizon.documentation.ird.fr/exl-doc/pleins_textes/divers14-05/08532.pdf

Dasgupta, S., Laplante, B., Meisner, C., Wheeler, D., Yan, J., (2009), The Impact of Sea Level Rise on Developing Countries: A Comparative Analysis. Climatic Change 93, p. 379‑388.

Delile, H., Goiron, J.-P., Bravard, J.-P., (2015a), Contribution du SIG à la reconstitution des processus de progradation de la Medjerda : quels enjeux pour le port d’Utique ? In Actes de la 1re Rencontre de l’ATEQ pour l’Étude du Quaternaire en Tunisie, Tunis, p. 133-153. https://shs.hal.science/halshs-01099895v1/document

Delile, H., Abichou, A., Gadhoum, A., Goiran, J-P., Pleuger, E., Monchambert, J-Y., Wilson, A., Fentrass, E., Quinn, J., Ben Jerbania, I., Ghozzi, F., (2015b), The Geoarchaeology of Utica, Tunisia: The Paleogeography of the Mejerda Delta and Hypotheses Concerning the Location of the Ancient Harbor. Geoarchaeology 30, p. 291‑306. https://doi.org/10.1002/gea.21514

Diatta, L. S., Fall, L. C. A., Sané, Y., (2022), Impacts de la dynamique du littoral entre Cabrousse et Boudiédiète (commune de Diembéring), Basse-Casamance, Revue Espaces Société Marocaine 58, p. 173‑198. https://doi.org/10.34874/IMIST.PRSM/EGSM/31136

Diatta, I., (2004), L’ouverture d’une brèche à travers la Langue de Barbarie (Saint Louis du Sénégal). Les autorités publiques et les conséquences de la rupture. Mémoire de maitrise, Université Gaston Berger, 120 p.

El Arrim, A., (1996), Étude d’impact de la dynamique sédimentaire et des aménagements sur la stabilité du littoral du Golfe de Tunis. Thèse de doctorat en géologie, Université de Tunis, Tunisie, 208 p.

Estournel, C., Mikolajczak, G., Ulses, C., Bourrin, F., Canals, M., Charmasson, S., Doxaran, D., Duhaut, T., Durrieu de Madron, X., Marsaleix, P., Palanques, A., Puig, P., Radakovitch, O., Sanchez-Vidal, A., Verney, R., (2023), Sediment dynamics in the Gulf of Lion (NW Mediterranean Sea) during two autumn-winter periods with contrasting meteorological conditions, Progress in Oceanography 210, 38 p. https://doi.org/10.1016/j.pocean.2022.102942

Fox-Kemper, B., Hewitt, H.T., Xiao, C., Aðalgeirsdóttir, G., Drijfhout, S.S., Edwards, T.L., et al., (2021), Ocean, Cryosphere and Sea Level Change. In: Climate Change 2021: The Physical Science Basis. IPCC, 2021, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, p. 1211‑1362. https://doi.org/10. 1017/9781009157896.011

Hinkel, J., Athanasios, J., Albrecht, C., Lincke, D., Wolff, C., (2015), Assessment of Costs of Sea-Level Rise in the Republic of Tunisia Including Costs and Benefits of Adaptation. PAC/RAC, Technical report, p. 1-62.

Hzami, A., Amrouni, O., Heggy, E., Mahé, G., Missaoui, H., (2024), On the Growing Socioeconomic Vulnerability of Southern Mediterranean Coastal Lagoons. Proceedings of IAHS 385, p. 377‑385. https://doi.org/10.5194/piahs-385-377-2024

Hzami, A., Heggy, E., Amrouni, O., Mahé, G., Maanan, M., Abdeljaouad, S., (2021), Alarming Coastal Vulnerability of the Deltaic and Sandy Beaches of North Africa, Scientific Reports 11, 15 p. https://doi.org/10.1038/s41598-020-77926-x

IPCC., (2021), Climate Change : The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. [Masson-Delmotte, V., P. Zhai, A. Pirani, S.L. Connors, C. Péan, S. Berger, N. Caud, Y. Chen, L. Goldfarb, M.I. Gomis, M. Huang, K. Leitzell, E. Lonnoy, J.B.R. Matthews, T.K. Maycock, T. Waterfield, O. Yelekçi, R. Yu, B. Zhou (eds.)]. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, https://doi.org/10.1017/9781009157896

IPCC, (2019), Special Report on the Ocean and Cryosphere in a Changing Climate [H.-O. Pörtner, D.C. Roberts, V. Masson-Delmotte, P. Zhai, M. Tignor, E. Poloczanska, K. Mintenbeck, A. Alegría, M. Nicolai, A. Okem, J. Petzold, B. Rama, N.M. Weyer (eds.)]. Cambridge University Press, Cambridge, UK and New York, NY, USA, p. 321–445. https://doi.org/10.1017/9781009157964.006

IPCC, (2013), Climate Change : The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [Stocker, T.F., D. Qin, G.-K. Plattner, M. Tignor, S.K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex, P.M. Midgley (eds.)]. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 1535 p.

Jauzein, A., (1971), Évolution récente du delta de la Medjerda : les agents de la morphogenèse. Travaux du Laboratoire de Géologie de l'Ecole Normale supérieure, vol. 5, p. 128-151.

Kotti, F., Dezileau, L., Mahé, G., Habaieb, H., Bentkaya, M., Dieulin, C., Amrouni, O., (2018), Etude de l’impacts 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. Journal des sciences de la Terre africaines, vol. 142, p. 226-233. https://doi.org/10.5194/piahs-377-67-2018

LCHF, (1982), Protection des plages sud de Tunis. 2e partie, Rapport du Laboratoire central d'hydraulique de France, ministère de l’Équipement, Direction des services aériens et maritimes, 37 p.

Louati, M., Saïdi, H., Zargouni, F., (2014), Shoreline change assessment using remote sensing and GIS techniques : a case study of the Medjerda delta coast, Tunisia. Arabian Journal of Geosciences 8, p. 4239-4255. http://dx.doi.org/10.1007/s12517-014-1472-1

Moldenhauer, K. M., Zielhofer, C., Faust, D., (2008), Heavy metals as indicators for Holocene sediment provenance in a semi-arid Mediterranean catchment in northern Tunisia, Quaternary International 189, p. 129‑134. https://doi.org/10.1016/j.quaint.2007.09.006

Oueslati, A., (2017), Urbanisation, aménagements touristiques et érosion des plages dans la berge sud de la Méditerranée : exemples au Maroc, en Tunisie et en Égypte. Actes du 3e Colloque La Méditerranée à l’invité de la géographie, Sousse, Tunisie, p. 50-94. https://www.bio.unifi.it/upload/sub/progetti/wadi/Oueslati_2017.pdf

Oueslati, A., Ben Charrada, R., Amira, S., (2017), Sur l’évolution récente de la flèche littorale de l’ancienne embouchure de Oued Majerda et ses conséquences (Tunisie nord-orientale). Revue Tunisienne de Géographie 46‑47, p. 181‑195.

Oueslati, A., Labidi, O., Elamri, Th., (2015), Atlas de la vulnérabilité du littoral tunisien à l’élévation du niveau marin, Publ. APAL, 67 p.

Oueslati, A., El Aroui, O., Sahtout, N., (2015), Sur la grande vulnérabilité du lido du complexe lagunaire de Ghar El Melh et de ses terres humides (Tunisie septentrionale), Méditerranée 125, p. 65‑73. https://doi.org/10.4000/mediterranee.7944

Paskoff, R., (1994), Le delta de la Medjerda (Tunisie) depuis l’Antiquité, Études rurales 133‑134, p. 15-29. https://www.persee.fr/doc/rural_0014-2182_1994_num_133_1_3451

Paquier, A. E., Mugica, J., Maspataud, A., (2021), Estimation de l’évolution potentille du trait de côte aux horizons 2040 et 2100 sur la Plaine orientale (Haute-Corse), Rapport final BRGM/FR-70722-FR, 213 p.

Pleuger, E., Goiron, J. Ph., Mazzini, I., Delile, H., Abichou, A., Gadhoum, A., Djerbi, H., Piotrowska, N., Wilson, A., Fentress, E., Ben Jerbania, I., Fagel, N., (2019a), Palaeogeographical and palaeoenvironmental reconstruction of the Medjerda delta (Tunisia) during the Holocene, Quaternary Science Reviews, vol. 220, p. 263-278. https://doi.org/10.1016/j.quascirev.2019.07.017

Pleuger, E., Goiron, J. Ph., Delile, H., Gadhoum, A., Abichou, A., Wilson, A., Fentress, E., Ben Jerbania, I., Ghozzi, F., Fagel, N., (2019b), Exploration of the maritime façade of Utica : The potential location of the Phoenician and Roman harbours. Quaternary International, vol. 511, p. 140-152. https://doi.org/10.1016/j.quaint.2019.04.007  

Reynier, F., (1950), Note sur le comblement de la basse Medjerda, Tunis, 18 p.

Samaali, H., (2011), Étude de l’évolution de l’occupation et de l’utilisation du sol dans le delta de Mejerda par télédétection et SIG. Thèse de doctorat en géographie, Université de Tunis, Tunisie, 376 p. NNT : fftel-00612952v2

Sammari, C. et Gana, S., (1995), Revue de l’hydrodynamique au large des côtes tunisiennes, Bulletin de l’Institut National Scientifique et Technique d’Océanographie et de Pêche, Salammbô 22, 25 p. file:///C:/Users/Saint-Luc/Downloads/hayadi,+1995-VOL+22-ART2.pdf

Slim, H., Trousset, P., Paskoff, R., Oueslati, A., Bonifay, M., Lenne, J., (2004), Le littoral de la Tunisie. Étude géoarchéologique et historique. Édit. CNRS, Paris, Antiquités africaines, Vol. 1, p. 3-308. https://www.persee.fr/doc/etaf_0768-2352_2004_mon_1_1

Sliti, M., (1990), Fonctionnement des brises lames dans le système marin littoral du golfe de Tunis. Thèse de doctorat en océanologie, Université de Bordeaux 1, Bordeaux, France, 257 p.

Sy, B. A., Bilbao I. A., Sy A. A., Sanchez-Pérez, I., Valido, S. R, (2013), Résultats du suivi 2010-2012 de l’évolution de la brèche ouverte su la langue de Barbarie au Sénégal et de ses conséquences. Physico-Géo vol. 7, p. 223-242. https://doi.org/10.4000/physio-geo.3569

Sweet, W. V., Kopp, R. E., Weaver, C. P., Obeysekera, J., Horton, R. M., Thieler E. R., Zervas, C., (2017), Global and regional sea level rise scenarios for the United States. NOAA Technical report NOS CO-OPS 083, 75 p.

Thior, M., (2020), Dynamique du littoral de la Casamance : caractéristiques morphodynamique, changements environnementaux et impacts socioéconomique. Thèse de doctorat en géographie, Université Assane SECK de Ziguinchor, Sénégal, 389 p.

Thieler, E. R., Hammar-Klose, E. S., (1999), National Assessment of Coastal Vulnerability to Sea-Level Rise : Preliminary Results for the U.S. Atlantic Coast, U.S. Geological Survey, Open-File Report 99-593, 1 sheet. https://pubs.usgs.gov/dds/dds68/reports/eastrep.pdf

Thiéblemont, R., Le Cozannet, G., Rohmer, J., Toimil, A., Álvarez-Cuesta, M., Losada, I.J., (2021), Deep Uncertainties in Shoreline Change Projections: an Extraprobabilistic Approach Applied to Sandy Beaches. Natural Hazards and Earth System Science 21, p. 2257‑2276. https://doi.org/10.5194/nhess-21-2257-2021

Vecchio, A., Anzidei, M., Serpelloni, E., (2024), Sea Level Rise Projections up to 2150 in the Northern Mediterranean Coasts, Environmental research19. 15 p. https://doi.org/10.1088/1748-9326/ad127e

Winter, R. C., Ruessink, G., (2017), Sensitivity Analysis of Climate Change Impacts on Dune Erosion: Case Study for the Dutch Holland Coast, Climatic Change 141, p. 685‑701. https://doi.org/10.1007/s10584-017-1922-3

Xiqing C., Erfeng Z., Hongqiang M., Zong Y., (2005), A Preliminary Analysis of Human Impacts on Sediment Discharges from the Yangtze, China, into the Sea, Journal of Coastal Research 21, p. 515‑521. https://doi.org/10.2112/03-0034.1

Top of page

Notes

1 Laboratoire Central d'Hydraulique de France.

2 DSM-OPT : Digital Surface Models from OPTical stereoscopic very-high resolution imagery.

3 ForM@Ter : pôle de données et de services Terre solide.

4 THEIA : pôle Thématique national des surfaces continentales.

5 DINAMIS : Dispositif Institutionnel National d'Accès Mutualisé en Imagerie Satellitaire.

Top of page

List of illustrations

Title Figure 1
Caption A : Esquisse géologique du delta de la Majerda. B : Géographie du delta. C : Rose des vents actifs au niveau du golfe de Tunis
Credits Crédits : A : Samaali, 2011 ; B : L.S. Diatta ; C : Khsiba, 2018
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-1.jpg
File image/jpeg, 390k
Title Tableau 1
Caption Caractéristiques des données géospatiales utilisées.
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-2.png
File image/png, 37k
Title Figure 2
Caption Évolution globale du linéaire côtier du delta de la Majerda, à partir des positions du trait de côte en 1972, 1985, 1995, 2005, 2015 et 2024. Zoom sur quelques secteurs
Credits Crédit : L. S. Diatta
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-3.jpg
File image/jpeg, 308k
Title Tableau 2
Caption Statistiques de la mobilité historique du trait de côte en mètre/an dans le delta de l’Oued Majerda entre 1972-1985-1995-2005-2015-2024. L’expression « Pas d’épis » renvoie aux périodes durant lesquelles les épis du port n’étaient pas encore prolongés, et ne séparaient donc pas la flèche de Foum El Oued en deux segments : nord et sud.
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-4.png
File image/png, 60k
Title Tableau 3
Caption Identification et qualification des sources d’incertitudes lors des différentes étapes de la cartographie des positions futures du trait de côte.
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-5.png
File image/png, 82k
Caption Avec Ept : erreur de position globale du trait de côte de l’année considérée ; Eg : erreur de géoréférencement de l’image ; Ee : erreur d’extraction du trait de côte ; Eo : erreur d’oscillation de marée ; E : erreur globale moyenne en m/an.
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-6.jpg
File image/jpeg, 8.1k
Title Tableau 4
Caption Marges d’erreurs estimées pour l’analyse de la cinématique du trait de côte au niveau du delta de la Majerda.La marge d’erreur liée à la quantification des dynamiques du rivage de toute la période étudiée (1985‑2024) est de ± 0,39 m/an.
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-7.png
File image/png, 17k
Title Figure 3
Caption Organigramme de la méthode d’évaluation des incertitudes des taux d’évolution du trait de côte utilisés pour les projections futures
Credits Figure adaptée de Paquier et al., 2021
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-8.jpg
File image/jpeg, 228k
Title Figure 4
Caption Cônes alluviaux dans le delta de la Majerda
Credits Source : Delile et al., 2015a
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-9.jpg
File image/jpeg, 584k
Title Figure 5
Caption Carte des anciennes positions hypothétiques du trait de côte et des processus de comblement du compartiment nord du delta de la Majerda
Credits Sources : Delile et al. 2015a ; 2015b
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-10.jpg
File image/jpeg, 561k
Title Figure 6
Caption Carte du processus de comblement hypothétique du compartiment nord du delta de la Majerda
Credits Source : Pleuger et al., 2019a
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-11.jpg
File image/jpeg, 573k
Title Figure 7
Caption Carte des anciennes positions hypothétiques du trait de côte du compartiment nord du delta de la Majerda
Credits Sources : Pleuger et al., 2019a ; 2019b
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-12.jpg
File image/jpeg, 690k
Title Figure 8
Caption Carte prospective de l’évolution future du trait de côte dans le delta de la Majerda à l’horizon 2050 (trait jaune), 2100 (trait rose) et 2124 (trait rouge). Zoom sur quelques segments caractéristiques des tendances d’évolution : 1. Bande de Sidi Ali El Mekki ; 2. La brèche de la sebkha de Ghar El Melh ; 3. Ancienne embouchure de la Majerda et bande nord de la flèche de Foum El Oued ; 4. Bande sud de la flèche de Foum El Oued ; 5. Embouchure actuelle de la Majerda ; 6. Bande allant de Raoued à Gammarth
Credits Crédit : L. S. Diatta
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-13.jpg
File image/jpeg, 502k
Title Figure 9 - Illustration des enjeux de la bande sableuse de Ghar El Melh
Caption 1) Vue d’ensemble de la partie nord de la bande. 2) Installations touristiques en bordure du rivage. 3) Installations touristiques sur pilotis. 4) Vue d’ensemble de la partie sud de la bande de Ghar El Melh
Credits Tous droits réservés
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-14.jpg
File image/jpeg, 478k
Title Figure 10
Caption Illustration du recul actuel du segment de côte situé immédiatement au nord du port de Kalâat Landlouss : ouverture de brèches périodiques
Credits Crédit : tous droits réservés
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-15.jpg
File image/jpeg, 252k
Title Figure 11
Caption Carte prospective de l’évolution future du trait de côte dans le delta de la Majerda dans l’hypothèse d’un effacement des ouvrages (épis) à l’horizon 2050 (trait jaune), 2100 (trait rose) et 2124 (trait rouge). Zoom sur les secteurs disposant d’ouvrages : 1) Bande de Sidi Ali El Mekki ; 2) Ancienne embouchure de la Majerda et flèche de Foum El Oued
Credits Crédit : L. S Diatta
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-16.jpg
File image/jpeg, 556k
Title Figure 12
Caption Cartographie globale des terres submersibles dans le delta de la Majerda aux échéances 2050, 2100 et 2124
Credits Crédit : L. S. Diatta
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-17.jpg
File image/jpeg, 377k
Title Figure 13
Caption Panorama des terres basses entre l’actuelle embouchure de la Majerda et la ville de Raoued
Credits Crédit : L. S. Diatta
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-18.jpg
File image/jpeg, 338k
Title Figure 14 - Exemple de terres submersibles de basses altitudes à Raoued
Caption A) terres basses urbanisées. B) Lotissements dans la sebkha
Credits Crédit : Chouari, 2016
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-19.jpg
File image/jpeg, 189k
Title Figure 14
Caption Reconstruction des positions anciennes, actuelles et futures du trait de côte dans le delta de la Majerda
Credits Crédit : L. S. Diatta
URL http://journals.openedition.org/mediterranee/docannexe/image/16587/img-20.jpg
File image/jpeg, 638k
Top of page

References

Electronic reference

Luc Simiroré Diatta, Katia Schörle, Semah Bettaieb, Luc Padhiou Mendy, Ameur Oueslati and Luc Lapierre, “Synthèse des dynamiques anciennes et analyse spatio-temporelle des modifications récentes, actuelles et futures de la ligne de rivage dans le delta de la Majerda, Tunisie”Méditerranée [Online], Paleoenvironment, Geoarchaeology, Historical Geography, Online since 29 April 2026, connection on 07 May 2026. URL: http://journals.openedition.org/mediterranee/16587; DOI: https://doi.org/10.4000/165j9

Top of page

About the authors

Luc Simiroré Diatta

Aix Marseille Université, CNRS, IRD, INRAE, CEREGE, Aix-en-Provence, France

Katia Schörle

Aix Marseille Université, CNRS, CCJ, Aix-en-Provence, France

Semah Bettaieb

Faculté des Lettres et sciences humaines de Sfax, Sfax, Université de Sfax, Tunisie

Luc Padhiou Mendy

Aix Marseille Université, CNRS, IRD, INRAE, CEREGE, Aix-en-Provence, France

Ameur Oueslati

Université de Tunis, Tunisie

By this author

Luc Lapierre

Française de Photogrammétrie et de Télédétection, Paris, AOrOc, France

Top of page

Copyright

CC-BY-SA-4.0

The text only may be used under licence CC BY-SA 4.0. All other elements (illustrations, imported files) may be subject to specific use terms.

Top of page
Search OpenEdition Search

You will be redirected to OpenEdition Search