Introduction

1Les côtes méditerranéennes sont aujourd’hui parmi les plus fortement aménagées au monde, avec un taux d’urbanisation d’environ 64 % en 2015 (PNUD, 2001 ; Lafitte et al., 2012 ; Cramer et al., 2018). L’enchevêtrement d’impacts anthropiques de diverses formes (touristique, portuaire, hydraulique, etc.) et de conditions naturelles complexes (géologique, morpho-sédimentaire, etc.) n’a fait qu’aggraver la vulnérabilité du littoral face au risque d’érosion du rivage (Anthony, 2014 ; Anthony et al., 2022 ; Besset, 2017 ; Hzami et al., 2021). Appartenant à la rive sud de la Méditerranée, la côte orientale de la Tunisie, qui s’étire entre Kelibia et Chebba (figure 1), n’a pas échappé de cette situation désastreuse, si bien qu’aujourd’hui 70 % de ce littoral souffre d’érosion grave rendant son budget sédimentaire de plus en plus déficitaire. L’étude de la cinématique de ce littoral, surtout depuis les années 1980, révèle qu’il est le plus souvent confronté à des risques d’érosion atteignant dans bien des cas des rythmes effrayants de plus de 8 m/an (cas de la plage au Sud d’El Kasba à Hammamet (Bada, 2020)). En fait, l’importance de la place qu’occupent les compartiments sédimentaires caractérisés par des évolutions négatives s’explique d’abord par des conditions naturelles. Ceci apparait surtout à travers la faiblesse de la topographie, la grande extension des terres humides et la médiocrité du réseau hydrographique (Oueslati, 1993, 2004 ; APAL, 1998, 2003 ; Brahim, 2001 ; Bada, 2007, 2017a). Cette vulnérabilité, ainsi que le risque d’érosion, vont s’accentuer davantage avec la multiplication hâtive des aménagements côtiers entrepris depuis quelques décennies. La construction des barrages et l’endiguement des cours d’eau, qui ont réduit encore les apports solides vers la mer, mais aussi le prélèvement de sables sur certains estrans, ont provoqué un appauvrissement remarquable du budget sédimentaire de la plage. Il en est aussi de la construction de ports et d’ouvrages de défense contre la mer qui ont largement perturbé la circulation des sédiments le long du rivage (Paskoff, 1993 ; Oueslati, 1999, 2004 ; Brahim, 2001, 2017 ; Fathallah et al., 2010 ; Labaïed & Oueslati, 1995 ; Bada, 2017a, 2020). Dans cette conjoncture menaçante va s’ajouter un autre indicateur de vulnérabilité. Il s’agit de l’élévation future du niveau moyen de la mer (Oueslati, 2004, 2015 ; Touzi & Ben Zakour, 2015 ; Bada, 2017b).

2Ce travail vise ainsi à déterminer le poids des indicateurs de la vulnérabilité morpho-sédimentaire et anthropique du littoral Kelibia-Chebba conjuguée avec la nouvelle conjoncture d’élévation du niveau moyen de la mer en effectuant une analyse multivariée d’interdépendance ou exploratoire multidimensionnelle d’analyse factorielle en composantes principales. L’objectif est de comprendre la structure de ces divers indicateurs ou variables de vulnérabilité et de déterminer, par classification des cellules, leurs contributions à la dynamique du littoral. Les indicateurs de vulnérabilité les plus contributifs au risque d’érosion marine vont servir par la suite à la production des cartes de synthèse.

Présentation de la zone d’étude : un cadre naturel déjà vulnérable, aggravé par une mauvaise politique d’aménagement

3La façade orientale du littoral tunisien, au Nord du golfe de Gabès, est limitée par deux promontoires bien marqués dans le paysage, celui de Kelibia au Nord et de Chebba au Sud. De tracé irrégulier, cette côte, qui s’étend sur environ 275 km, se distingue par une morphologie très variée (figure 2). Les plages sableuses y occupent toutefois la place dominante (196 km, soit 71 % de l’ensemble du trait de côte), alternant avec une série de promontoires rocheux (Kelibia, Hammamet, Hergla, Monastir, Eddimess, Mahdia, Salakta et Chebba), les plus étendues se trouvant dans le golfe de Hammamet et sur la côte d’El Ghedhabna. Cette morphologie de côtes basses sableuses s’explique par la faiblesse de la topographie où le rivage est souvent bordé par des plaines alluviales très basses : il est ainsi rare qu’elles dépassent 5 m d’altitude notamment à l’approche de la mer en favorisant ainsi l’extension des terres humides de type sebkhas et chotts. Ceci constitue en fait un facteur de fragilité au risque d’érosion surtout dans la nouvelle conjoncture de l’élévation du niveau marin. Par ailleurs, la vulnérabilité de la côte Kelibia-Chebba trouve son origine également à travers la modestie et l’intermittence du réseau hydrographique. C’est ainsi que 70 % du rivage ne bénéficie pas des apports sédimentaires des oueds qui sont soit absents soit étalent leurs dépôts dans des dépressions ou des terres humides littorales (notamment dans la moitié méridionale du golfe de Hammamet, entre Bouficha et Hergla, mais aussi dans le Sahel méridional, entre Rass Edimess et Chebba).

4En mer, l’effet morpho-dynamique de la marée, de type semi diurne et dont l’amplitude moyenne est faible (0,4 m), est généralement masqué par l’action de la houle dont la hauteur moyenne est d’environ 1,5 m. En fait, cette dernière, malgré l’abaissement de la bathymétrie et la remarquable extension de la plate-forme continentale, notamment à la hauteur du promontoire de Monastir (l’isobathe 5 m se situe entre 400 m et 18 km du rivage), pouvait connaitre des amplitudes exceptionnelles dépassant 4 m notamment pendant les fortes tempêtes (Oueslati, 2016 ; Brahim, 2014 ; Bada, 2017a). La prédominance des secteurs Nord et Est de la houle ne peuvent conduire qu’à des courants côtiers adjacents dans une direction principale Nord-Sud le long du rivage.

5Les vents, quant à eux, principalement de secteurs N.O et N.E, sont régulièrement violents. En effet, les données de l’Office National de la Météorologie (entre 1982 et 2011), témoignent de la faiblesse des fréquences annuelles des vents calmes à moyens notamment à Monastir (11,5 %). Ceux qui sont efficaces, pouvant déplacer les grains de sables (≥ à 5 m/s), présentent des pourcentages assez élevés atteignant par exemple 74,5 % à Nabeul et 88,5 % à Monastir. Dans cette dernière, les vents ont atteint des vitesses moyennes records dépassant 90 km/h en octobre 1984.

6La vulnérabilité de la côte Kelibia-Chebba est accrue par la multiplication des aménagements, souvent implantés de manière inadaptée au bon fonctionnement naturel des différentes unités du littoral. Très faiblement anthropisée jusqu’au milieu du vingtième siècle, cette côte est devenue parmi les plus fortement aménagées de la rive Sud de la Méditerranée (Oueslati, 2010). Ainsi, en 1962, les espaces bâtis, souvent dispersés, ne couvraient que 1273 ha. En 2010, ces espaces occupaient 14 594 ha formant souvent de grands pôles urbains de très forte densité, comme par exemple la zone de Nabeul-Hammamet et de Sousse-Monastir (Bada, 2017a). C’est notamment l’extension urbaine, l’implantation des unités touristiques et la multiplication des aménagements portuaires (22 ports) et aussi hydrauliques (33 barrages collinaires) qui caractérisent la forte anthropisation de la zone d’étude et qui ne peuvent pas être sans impact sur l’équilibre du bilan sédimentaire (Pasckoff, 1985 ; Oueslati, 1993, 2004, 2010 ; Ben Haj 2004 ; Brahim 2001, 2017 ; Fathallah et al., 2010 ; Souayed et al., 2011 ; Bada 2017b, 2020).

Données et méthodes

7Le transport sédimentaire littoral est animé par la houle et opère généralement dans le cadre de cellules sédimentaires aux limites naturelles comme les falaises rocheuses, ou artificielles, comme les digues de port (Carter, 1988 ; BRGM, 2015 ; OEC, 2019 ; Bertoni et al., 2021). Le concept de cellule littorale a été généralement utilisé dans le cadre d’études du budget sédimentaire dans lequel les processus morphodynamiques peuvent être non négligeables. L’accent y est mis sur la définition de chaque cellule et les gains ou pertes de sédiments qui déterminent l’accrétion ou l’érosion, et donc potentiellement la vulnérabilité (Bowen et al., 1966 ; Bray et al., 1995 ; Schwartz, 2005 ; Patsch, 2006 ; van Rijn, 2011 ; Anthony et al., 2019). Ainsi, pour identifier et préciser les limites des cellules littorales de la côte Kelibia-Chebba, cinq critères ont été retenus, basés principalement sur le bilan sédimentaire et le sens de transport des sédiments, à savoir :

• La quantification de l’évolution du trait de côte (m/an) par le script DSAS de l’U.S. Geology Survey (Digital Shoreline Analysis System). Il a pour objectif de dégager les secteurs qui avancent, dits d’accumulation ou sites puits, et ceux qui reculent, dits d’érosion ou sources. La reconstitution des traits de côte a été effectuée par fusion de diverses sources (photos aériennes aux échelles 1/5000 (1962) et 1/10 000 (1996 et 2002) et images satellitaires de 2004, 2008, 2010 de QuickBird et 2021 et 2023 de Maxar (GEO1) avec respectivement une résolution de 2,4 et 0,4 m au sol). La procédure d’orthorectification des photos aériennes, combinée au processus de vectorisation des lignes de rivage par l’outil ArcGis, a comporté des incertitudes géométriques estimées au total à ± 6 mètres ;

• La quantification du volume du transport littoral et la détermination de sa direction par l’équation du Coastal Engineering Research Center (CERC, 2020), dérivée de Komar et Inman (1970). Son objectif est d'obtenir des estimations recevables sur le transport sédimentaire ou la dérive littorale potentielle afin de prévoir le processus de dépôt et/ou d’érosion dans chaque secteur de différentes cellules sédimentaires ;

• L’analyse des accumulations sableuses liées aux ouvrages de défense côtière et aux ports qui occupent généralement tous les promontoires du littoral Kelibia-Chebba constituant ainsi des limites artificielles des unités sédimentaires. C’est ainsi que sous l’effet de leurs jetées qui constituent des vrais barrages au transport des sédiments de la dérive littorale et par conséquent des changements dans la position du trait de côte, deux bilans sédimentaires ont été dégagés. Un premier bilan tend à devenir déficitaire correspondant à un secteur source ou d’érosion et un deuxième bilan, positif, représente un secteur d’accumulation ou puits alimenté par un secteur d’érosion ou de transit ;

• L’analyse de l’évolution (surtout de la direction) des flèches sableuses et des orientations des embouchures des cours d’eau exoréiques notamment sur la partie d’avant plage. Elles constituent de bons indicateurs sur le sens du transit sédimentaire et par conséquent sur la succession de différents secteurs formant une cellule littorale ;

• L’analyse de la variation longitudinale de la granulométrie des plages, notamment la couche active du bas de plage pour déceler s’il existe une variation systématique de la granulométrie des sédiments. En fait, ce phénomène est lié en particulier à la variation du niveau d'énergie des houles et des courants qu’elles induisent. Au sein d’une cellule littorale, la granulométrie diminue assez régulièrement en partant du secteur d’érosion, souvent en falaise où la fraction grossière domine, vers un autre d’accumulation où la granulométrie fine l’emporte.

8L’application de cette approche multicritère a permis d’identifier 13 cellules sédimentaires de différentes longueurs et caractéristiques (Bada, 2017a). Ces cellules sont compartimentées en 35 secteurs de différents états de fonctionnement (tableau 1).

Tableau 1 - Les cellules sédimentaires de la côte Kelibia-Chebba

Id	Cellule sédimentaire	Secteur	Longueur (km)	Total (km)	Pourcentage par rapport au total
1	Kelibia-Mâamoura	Erosion	17,5	51,6	20,3
		Accumulation	34,1
2	Mâamoura-Port de Beni Khiyar	Erosion	0,5	1,9	0,7
		Transit	1,2
		Accumulation	0,2
3	Port de Beni Khiyar-Rive gauche de l’Oued Belgaied	Erosion	11,7	20,3	8
		Transit	6,5
		Accumulation	2,1
4	Rive droite de l’Oued Belgaied-Port de plaisance de Yasmine El Hammamet	Erosion	3,2	7,1	2,8
		Transit	1,2
		Accumulation	2,7
5	Port de plaisance de Yasmine El Hammamet- Port de Hergla	Erosion	0,8	38,9	15,3
		Transit	12,9
		Accumulation	23,5
		Erosion	1,7
6	Port de Hergla-Port de Plaisance de Kantaoui	Erosion	7,3	19,1	7,5
		Accumulation	11,8
7	Port de plaisance de Kantaoui-Port de Sousse	Erosion	4,9	8,8	3,5
		Accumulation	4,0
8	Port de Sousse-Bassin de tranquillisation à Sidi Abdelhamid	Erosion	3,0	4,0	1,6
		Accumulation	1,0
9	Bassin de tranquillisation à Sidi Abdelhamid-Marina de Monastir	Erosion	3,6	14,5	5,7
		Accumulation	5,9
		Erosion	5,0
10	Marina de Monastir-Flèche littorale d’El Enf	Erosion	5,8	9,0	3,5
		Accumulation	3,2
11	Khnis-Mahdia	Erosion	12	49,1	19,3
		Transit	25
		Accumulation	3,5
		Erosion	8,6
12	Mahdia-Salakta	Erosion	3,5	11,5	4,5
		Transit	7,1
		Accumulation	0,9
13	Salakta-El Ghedhabna	Erosion	8,0	18,2	7,2
		Accumulation	10,2
		Total	254	254	100

9Dans ce travail, une liste de 14 variables ou indicateurs (naturels et anthropiques) a été établie. Ces variables sont susceptibles de renseigner sur la vulnérabilité des cellules sédimentaires en rapport avec les actions d’aménagement et des modes d’occupation combinés au facteur de l’élévation future du niveau moyen de la mer. La liste des indicateurs a été réduite par la suite au strict minimum en évitant les variables redondantes ou totalement corrélées par des combinaisons arithmétiques (Bada, 2017a). Notons ici que l’indicateur du rythme d’évolution sédimentaire faisant référence au bilan sédimentaire volumique a été préféré aux autres (linéaire et surfacique) qui sont restreints à une seule dimension spatiale linéaire (Sabatier, 2001).

10Ainsi, dans l'objectif de mettre en valeur les relations potentiellement structurantes entre ces indicateurs et les différentes cellules littorales de la côte Kelibia-Chebba, nous n’avons retenu, pour les 13 individus correspondant aux cellules littorales, que 10 indicateurs qui ont été jugés a priori les plus pertinents et les plus significatifs (figure 3).

11La quantification des poids des indicateurs de vulnérabilité de la côte Kelibia-Chebba face au risque d’érosion va reposer, dans ce travail, sur un traitement statistique exploratoire dans un “hyper-espace” multidimensionnel par la méthode d’Analyse factorielle en Composantes Principales (ACP). Les indicateurs retenus notamment de type qualitatifs (ports, morphologie et géologie, ouvrages de défense côtière) ont été initialement convertis en des données numériques par processus de codage. Pour atténuer l’effet d’échelle, les 10 indicateurs ont été, dans un second temps, centrés et réduits sous le logiciel Statistica.

Résultats

Extraction des axes factoriels optimaux

12L’application du processus d’Analyse factorielle en Composantes Principales, avec la recommandation de Kaiser (1958) proposant une valeur propre supérieure à 1, nous a permis de retenir les 4 premiers facteurs ou axes expliquant 85,6 % de l’inertie totale (tableau 2). Ainsi, nous allons passer d’une dimension de 10 à une autre de 4 ce qui va réduire considérablement la dimension initiale.

Tableau 2 - Valeurs propres et statistiques associées des facteurs liés à la vulnérabilité des cellules sédimentaires de la côte entre Kelibia et Chebba au risque d’érosion

Factoriel	Valeur propre	% de la variance	% cumulé de la variance	Histogramme
1	3,68	36,78	36,78
2	2,18	21,82	58,59
3	1,64	16,42	75,01
4	1,06	10,57	85,58
5	0,72	7,20	92,79
6	0,49	4,87	97,66
7	0,13	1,30	98,95
8	0,09	0,86	99,81
9	0,01	0,13	99,94
10	0,01	0,06	100,00

Examen de la structure des relations points-indicateurs ou variables de vulnérabilité 

13Afin d'obtenir une structure factorielle clairement marquée par de fortes corrélations avec les indicateurs de vulnérabilité ou variables, nous avons appliqué en premier temps une rotation Varimax brut pour maximiser les variances des poids factoriels dans l'espace original (Forina et al., 1988). La nouvelle structure factorielle est donnée dans le tableau 3 avec des poids marqués α ≥ à 0,70.

Tableau 3 - Poids factoriels des variables avec une rotation Varimax brut

ID	Variable	Fact 1	Fact 2	Fact 3	Fact 4
1	Rythme d’évolution sédimentaire volumique (m3/an)	0,04	0,89	0,01	-0,12
2	Ouvrages de défense	0,56	0,10	-0,40	0,74
3	Infrastructures hydrauliques	0,14	0,18	0,95	0,09
4	Ports	-0,08	-0,31	0,14	0,84
5	Espaces bâtis (ha)	0,81	0,14	0,05	0,35
6	Densité de la population côtière (h/km²)	0,84	0,32	-0,14	0,19
7	Espaces touristiques (ha)	-0,92	0,47	0,14	-0,14
8	Fréquentation des plages (estivant/m²)	0,87	0,06	0,02	-0,14
9	Morphologie et géologie	-0,02	-0,79	-0,34	0,22
10	Elévation du niveau marin	-0,22	-0,70	0,58	-0,16

14Ainsi, il ressort de ce tableau que :

• Le premier axe ou composante principale est constitué par 4 variables (5, 6, 7 et 8). Les variables 5, 6 et 8 sont loin du centre, mais proches les unes par rapport aux autres (figure 4). Cette configuration permet de conclure qu’il existe entre ces trois variables des corrélations fortement positives avec le 1^er facteur, c’est-à-dire qu’à de fortes valeurs du facteur correspondent de fortes valeurs de ces variables et vice-versa ;

• Le second axe est lié au milieu physique et à son évolution dans une conjoncture future d’élévation du niveau de la mer. En effet, il est composé par les variables 1, 9 et 10 ;

• Le troisième axe où la corrélation de la 3^ème variable est la plus forte (r = 0,95) et est liée aux infrastructures hydrauliques (3) ;

• Enfin, une amélioration nette se voit sur le quatrième axe après rotation où nous avons gagné une deuxième variable concernant les ouvrages de défense (2) qui vient de s’ajouter à celle des aménagements portuaires (4). 

15Les contributions de différentes variables (en %) sont très fluctuantes d’un axe à un autre (tableau 4). Ainsi, sur :

• Le premier axe factoriel, les variables liées à la population et l’occupation de son espace ont contribué largement dans la construction de cet axe avec des pourcentages souvent supérieurs au seuil moyen calculé à 10 %. Ces contributions varient entre 20 et 25 %. Pour les autres variables de ce même axe, les contributions sont faibles, voire même nulles ;

• Sur le deuxième axe, les contributions qui sont liées au milieu physique et son évolution sont les plus importantes. Elles vont de 22 à 36 % ;

• Quant au 3^ème axe, la contribution est la plus forte, mais elle n’intéresse que la seule variable liée aux infrastructures hydrauliques avec 57 %. Ceci doit être vu avec prudence car cette variable est très différente des autres et peut gêner l’étude du reste des variables ;

• Enfin, les contributions des variables du 4^ème axe, qui sont liées aux ouvrages de défense et aux infrastructures portuaires, sont aussi importantes avec respectivement des variances de 29 et 50 %.

Tableau 4 - Contributions des variables par axe factoriel après rotation Varimax brut (en %)

ID	Variable	Fact 1	Fact 2	Fact 3	Fact 4
1	Rythme d’évolution sédimentaire volumique (m3/an)	0,04	35,68	0,01	1,09
2	Ouvrages de défense	9,18	0,48	10,37	28,81
3	Infrastructures hydrauliques	0,57	1,52	57,18	0,54
4	Ports	0,20	4,36	1,18	49,96
5	Espaces bâtis (ha)	19,61	0,89	0,17	8,81
6	Densité de la population côtière (h/km²)	21,11	4,71	1,34	2,68
7	Espaces touristiques (ha)	25,24	2,33	1,31	1,48
8	Fréquentation des plages (estivant/m²)	22,56	0,16	0,03	1,48
9	Morphologie et géologie	0,01	28,12	7,17	3,42
10	Elévation du niveau marin	1,47	21,76	21,24	1,73

Examen de la structure des relations points-individus

16Nous avons réalisé une deuxième analyse factorielle en composantes principales qui nous a permis d’obtenir un second point de vue sur les observations ou les individus correspondant aux différentes cellules sédimentaires. Pour les résultats relatifs aux individus, nous remarquons que leurs coordonnées ou scores sur les axes factoriels marquent une opposition très nette notamment sur les deux premiers axes (figure 5). Ces deux axes résument près de 59 % de l’information totale du tableau initial.

17Cette opposition se voit surtout entre deux groupes de cellules sédimentaires. Un premier groupe de cellules qui se caractérise, notamment, par des activités touristiques dominantes à savoir : les cellules situées entre le port de Beni Khiyar et la rive gauche de l’Oued Belgaied, entre la rive droite de l’Oued Belgaied et le port de Yasmine El Hammamet, entre le port de plaisance de Kantaoui et le port de Sousse et entre le bassin de tranquillisation à Sidi Abdelhamid et la Marina de Monastir. Les scores factoriels de ces cellules sur le premier facteur sont négatifs et varient entre -2,21 et -2,59 (tableau 5). Pour le deuxième groupe de cellules, les actions de boisements sur les dunes littorales (figures 6 et 7) et les plages larges, très souvent relayées par des sebkhas qui priment. Nous désignons les cellules Port de Yasmine El Hammamet-Port de Hergla, Salakta-El Ghedhabna et Kelibia-Mâamoura dont les scores sont positifs et fluctuent sur le premier axe entre 2,41 et 2,74.

18Alors que sur le deuxième facteur, qui résume 21,8 % de l’inertie totale, l’opposition indique plutôt les variantes du premier axe. Ainsi sur cet axe, on remarque, par exemple, que les cellules qui se situent entre le port de plaisance de Kantaoui et le port de Sousse et entre le bassin de tranquillisation à Sidi Abdelhamid et la Marina de Monastir s’opposent parfaitement aux cellules qui s’étirent entre Beni Khiyar et le port de Yasmine El Hammamet.

Tableau 5 - Coordonnées factorielles des individus basées sur les corrélations

ID	Individu	Fact 1		Fact 2	Fact 3	Fact 4
1	Kelibia-Mâamoura	2,41		0,99	-2,55	-0,82
2	Mâamoura-Port Beni Khiyar	0,25		-1,27	0,63	1,48
3	Port de Beni Khiyar-Rive gauche de l'Oued Belgaied	-2,59		-1,45	-1,31	-1,80
4	Rive droite de l'Oued Belgaied-Yasmina Hammamet	-2,45		-2,16	-1,75	1,29
5	Port Yasmine El Hammamet-Port de Hergla	2,53		1,75	-1,64	0,87
6	Port de Hergla-Port de Plaisance de Kantaoui	-0,64		-0,12	0,32	1,42
7	Port de plaisance de Kantaoui- Port de Sousse	-2,21		1,33	0,56	-0,67
8	Port de Sousse-Bassin de tranquillisation à Sidi Abdelhamid	1,71		0,30	1,66	0,12
9	Bassin de tranquillisation à Sidi Abdelhamid-Marina de Monastir	-2,33		2,41	0,73	0,30
10	Marina de Monastir-Flèche littorale d'El Enf	-0,32		0,83	0,86	-1,30
11	Khnis-Mahdia	-0,31		1,00	0,52	0,38
12	Mahdia-Salakta	1,22		-1,94	1,43	-0,70
13	Salakta-El Ghedhabna	2,74		-1,68	0,54	-0,59

19Comme pour les variables, nous avons calculé les contributions de différentes cellules sédimentaires par axe factoriel après rotation Varimax brut. Les résultats, fournis dans le tableau 6, montrent que les contributions significatives qui dépassent le seuil moyen de 7,7 % sont importantes notamment sur le premier axe factoriel. Sur cet axe, 54 % des individus ont des contributions qui varient, en termes de rapport, entre 1,3 et 2. Pour les autres axes, le nombre d’individus dont les contributions sont supérieures au seuil moyen a légèrement glissé à la baisse (entre 5 et 6). Toutefois, il faut signaler que leurs fluctuations autour du seuil moyen sont plus prononcées sur le 3^ème axe où nous avons calculé des rapports entre 1,0 et 3,9, soient des pourcentages qui fluctuent entre 8,0 et 30,4 %.

Tableau 6 - Contributions des individus par axe factoriel après rotation Varimax brut (en %)

ID	Individu	Fact 1	Fact 2	Fact 3	Fact 4
1	Kelibia-Mâamoura	12,13	3,44	30,36	4,94
2	Mâamoura-Port Beni Khiyar	0,13	5,72	1,84	16,05
3	Port de Beni Khiyar-Rive gauche de l'Oued Belgaied	14,04	7,40	8,00	23,47
4	Rive droite de l'Oued Belgaied-Yasmina Hammamet	12,54	16,43	14,28	12,02
5	Port Yasmine El Hammamet-Port de Hergla	13,38	10,81	12,59	5,47
6	Port de Hergla-Port de Plaisance de Kantaoui	0,86	0,05	0,49	14,68
7	Port de plaisance de Kantaoui- Port de Sousse	10,21	6,26	1,46	3,26
8	Port de Sousse-Bassin de tranquillisation à Sidi Abdelhamid	6,09	0,31	12,85	0,11
9	Bassin de tranquillisation à Sidi Abdelhamid-Marina de Monastir	11,37	20,47	2,50	0,68
10	Marina de Monastir-Flèche littorale d'El Enf	0,22	2,41	3,46	12,22
11	Khnis-Mahdia	0,20	3,56	1,25	1,06
12	Mahdia-Salakta	3,13	13,25	9,57	3,52
13	Salakta-El Ghedhabna	15,69	9,90	1,34	2,51

Classification des indicateurs-individus sur les axes de l’ACP

20Dans un second temps, les coordonnées des variables et des individus sur les facteurs de l’ACP ont fait l’objet d’une classification grâce à une représentation simultanée dans l’espace des facteurs. Le graphique de Biplot ou cercle de corrélation (Gower et al., 1996) est réservé à ce type de représentation (figure 8). Le graphique de Biplot met en lumière une typologie de trois partitions (figure 9) dont des variantes ou sous-partitions ont été retenues aussi en se basant sur le poids du deuxième axe factoriel résumant 21,8 % du tableau initial. Les résultats de cette classification sont fournis dans le tableau 7.

21Dans cette partie, nous analyserons chaque partition, ainsi que les variantes qui en résultent, en mettant l’accent, notamment, sur les variables les plus contributives à la vulnérabilité de différentes cellules sédimentaires face au risque d’érosion.

Tableau 7 - Partitions et sous-partitions des cellules sédimentaires littorales de la côte Kelibia-Chebba selon les coordonnées des variables et des individus sur les facteurs de l’ACP

ID	Individu ou cellule sédimentaire	Partition
		Principale	Sous-partition
1	Kelibia-Mâamoura	1	1,1
5	Port Yasmine El Hammamet-Port de Hergla	1	1,1
13	Salakta-El Ghedhabna	1	1,2
3	Port de Beni Khiyar-Rive gauche de l'Oued Belgaied	2	2,1
4	Rive droite de l'Oued Belgaied-Yasmina Hammamet	2	2,1
7	Port de plaisance de Kantaoui- Port de Sousse	2	2,2
9	Bassin de tranquillisation à Sidi Abdelhamid-Marina de Monastir	2	2,2
2	Mâamoura-Port Beni Khiyar	3	3,1
12	Mahdia-Salakta	3	3,1
6	Port de Hergla-Port de Plaisance de Kantaoui	3	3,2
8	Port de Sousse-Bassin de tranquillisation à Sidi Abdelhamid	3	3,2
10	Marina de Monastir-Flèche littorale d'El Enf	3	3,2
11	Khnis-Mahdia	3	3,2

La première partition des cellules littorales : une partition de cellules potentiellement très vulnérable au risque d’érosion

22Cette partition comprend les trois plus longues cellules sédimentaires de la côte Kelibia-Chebba. Il s’agit des cellules situées entre Kelibia et Mâamoura, entre la marina de Yasmine El Hammamet et le port de Hergla et enfin, celle située entre Salakta et El Ghedhabna. Deux sous-partitions ont y été toutefois distinguées selon leurs cordonnées sur les facteurs de l’ACP.

• La première sous-partition renferme les deux cellules Kelibia-Mâamoura et Port de plaisance de Yasmine El Hammamet-Port de pêche de Hergla. Dans ces deux cellules, au moins 60 % des indicateurs (6/10) participent fortement à la vulnérabilité et la fragilisation des côtes. Les contributions figurent notamment sur le premier axe factoriel qui tient à lui seul 40 % des variables dont les pourcentages varient entre 20 et 25 %. En fait, les quatre indicateurs ou variables à savoir « Espaces touristiques (ha) », « Fréquentation des plages (estivant/m²) », « Densité de la population côtière (h/km²) » et « Espaces bâtis » sont les plus contributifs du premier axe factoriel avec une variance totale de 88,5 %. Ainsi, d’après le tableau 8, il ressort que l’indicateur « Espaces touristiques » a contrôlé plus de 25 %. Cette importante contribution se voit dans la cellule sédimentaire qui s’étire entre le port de plaisance de Yasmine El Hammamet et le port de pêche de Hergla et précisément dans le secteur d’érosion entre les hôtels Vincci Lella Beya et El Maheri où 41 unités hôtelières ont été construites depuis 1995. La comparaison des photographies aériennes de 1962 et 1996 révèle que ce secteur était en équilibre dominé par des terrains souvent marécageux et sa plage était large de quelques dizaines de mètres.

23En plus des espaces touristiques et la fréquentation des plages, la densité de la population est également un autre indicateur explicatif de la vulnérabilité des côtes dans ces deux cellules littorales. En fait, cette variable caractérise le premier axe factoriel par un poids qui vaut près de 21 %. Le poids de la densité de la population se voit très nettement surtout dans les villes côtières de la cellule sédimentaire de Kelibia-Mâamoura où les taux d’accroissement globaux de la densité de la population entre 1984 et 2022 ont dépassé 150 % notamment à Beni Khiyar (321,8 %). Une telle évolution des densités de la population est souvent accompagnée d’une extension des espaces bâtis. Ces derniers ont contribué sur le premier axe factoriel de l’ACP pour 19,6 %. Les villes de la cellule sédimentaire de Kelibia-Mâamoura sont les lieux favoris de cette extension. Ainsi, entre 1962 et 2022, la demande en espace bâti a été évaluée en moyenne de 11,8 ha/an. Les villes de Korba (figure 10) et de Menzel Temime en sont des exemples très illustratifs puisque nous avons constaté une extension des zones urbaines de plus en plus rapide notamment en direction des plages.

24Quant au deuxième axe factoriel de l’ACP, la variable érosion marine exprimée par le rythme d’évolution sédimentaire volumique (m³/an) est la seule dont la contribution à la vulnérabilité est supérieure au seuil moyen avec près de 36 %. Cette importante contribution se manifeste notamment dans le secteur d’érosion situé au Sud de la cellule sédimentaire qui s’étend entre le port de plaisance de Yasmine El Hammamet et le port de pêche de Hergla où nous avons calculé, entre 2002 et 2008, une vitesse de recul exceptionnelle atteignant par endroit 6 m/an (figure 11). Un tel rythme a engendré une déperdition considérable des sédiments, évaluée, par la formule CERC (Coastal Engineering Research Center, 1984), à 23 120 m³. Le poids important de la variable érosion marine se manifeste également entre le port de Kelibia et l’hôtel Mamounia où le rythme de retrait a été compris, pour la période 1962-2010, entre 0,8 et 1,4 m/an, soit une perte sédimentaire volumique estimée à 102 927 m³.

25Pour le troisième axe de l’ACP, la variable « Infrastructure hydraulique » est de loin la plus contributive à la vulnérabilité du littoral avec près de 57 % (5,7 fois le seuil moyen). Ceci est associé notamment à la cellule sédimentaire située entre Kelibia et Mâamoura où divers travaux de CES (Travaux de Conservation des Eaux et des Sols), un grand barrage (Lebna), 13 barrages collinaires et 6 lacs (Direction Générale des Barrages et des Grands Travaux Hydrauliques / DGBGTH, 2014) dont la plupart ont été construits entre 1995 et 2001, ont fortement contribué à la réduction, voire même à l’arrêt total du transport solide vers la plage. Le barrage collinaire d’El Harouni, situé en amont d’une côte en érosion qui s’étend au Sud de Sabkhet Tagdimen, en est un exemple typique. En effet, cet ouvrage a retenu, depuis 1995, presque la totalité des eaux et des sédiments charriés, estimés à 3,5 Mm³ (HYDROMED, 2001). A trois kilomètres en aval, la mer ne reçoit que quelque 0,22 Mm³. Un tel volume d’eau ne permet de véhiculer qu’un très faible volume en sédiments évalué à quelques milliers de m³.

26Enfin, sur le quatrième axe de l’ACP, la contribution des ouvrages de défense côtière de type lourd (enrochement) est aussi la seule variable qui a un poids factoriel important de l’ordre de 29 %. La contribution à la vulnérabilité se figure notamment à travers les aménagements en enrochements qui ont été faits surtout au Sud du port de Kelibia dans la zone de Mamounia et au Nord du port de pêche de Hergla. Au niveau de ce dernier, les travaux d’enrochement, réalisés sur plus de 350 m en 2002 puis jusqu’à 1050 m en 2017, ont contribué au contraire à la reprise de l’action érosive, notamment, au Nord en direction de l’Oued Ammouch. La vulnérabilité de la plage se voit également à travers l’abaissement continu du profil transversal de la zone située aux pieds de la ligne d’enrochements suite aux phénomènes d’arrachement des sédiments dus au déferlement des eaux marines.

Tableau 8 - Les variables les plus contributives dans la vulnérabilité des deux cellules sédimentaires situées entre Kelibia et Mâamoura et entre le port de plaisance Yasmine El Hammamet et le port de Hergla (avec rotation Varimax brut)

ID	Variable	Fact 1	Fact 2	Fact 3	Fact 4
1	Rythme d’évolution sédimentaire volumique (m3/an)	0,04	35,68	0,01	1,09
2	Ouvrages de défense	9,18	0,48	10,37	28,81
3	Infrastructures hydrauliques	0,57	1,52	57,18	0,54
5	Espaces bâtis (ha)	19,61	0,89	0,17	8,81
6	Densité de la population côtière (h/km²)	21,11	4,71	1,34	2,68
7	Espaces touristiques (ha)	25,24	2,33	1,31	1,48
8	Fréquentation des plages (estivant/m²)	22,56	0,16	0,03	1,48

Unité : En pourcent

27Toutefois, d’autres variantes d’aménagements de défense souple (forêt, ganivelles etc.) ont joué un rôle parfois positif dans le maintien de la côte contre les attaques de l’érosion marine et éolienne. Le rôle de la forêt d’El Madfoun, par sa fixation des sables des dunes bordières, en est un bon exemple. Les ganivelles installées sur les dunes de Korba (photo 1) ainsi que sur la plage de Hergla (photo 2) sont également des interventions efficaces, notamment contre l’érosion éolienne, pour stabiliser les dunes.

• La deuxième sous-partition comporte une seule cellule sédimentaire comprise entre Salakta et El Ghedhabna où deux principales variables ont contribué à sa vulnérabilité au risque d’érosion marine en occupant chacune un axe factoriel (1^er et 2^ème axe).

28En premier lieu figure l’indicateur « Espaces touristiques » qui résume le premier axe de l’ACP avec un poids de 25 % (tableau 9). En réalité, cette contribution à la vulnérabilité du littoral, bien qu’elle soit importante, n’intéresse que des espaces ponctuels situés au Sud de Châabna et principalement dans la forêt d’El Ghedhabna au niveau du Centre des stages. Pendant la période estivale, cet espace se transforme en un lieu très fréquenté par les estivants très souvent en famille qui viennent surtout de Mahdia, de Kairouan et de Sidi Bouzid. Ces estivants génèrent souvent une dégradation du milieu, principalement sur la façade maritime, soit par destruction de la forêt ou par ouverture des sentiers parfois larges à travers les dunes bordières.

29Au deuxième plan apparait la variable conjuguée entre la morphologie et la géologie de la côte sur le deuxième axe factoriel de l’ACP avec 28 % de la variance totale. En fait, le trait de côte situé entre Salakta et El Ghedhabna est mixte (alternance entre côtes sableuses et côtes rocheuses), très découpé et orienté dans le sens des houles dominantes. Les angles d'incidence des houles sont assez favorables à des dérives littorales venant notamment du Sud-est et du Nord-ouest. Cependant, la nature très hétérogène des formations littorales confère à cette côte un potentiel morpho-dynamique relativement élevé. Ce potentiel s'exerce notamment sur la côte à falaises entre Salakta et El Alia constituée essentiellement de grès de plages et de dunes consolidées. Il s'exerce aussi sur les sables du cordon littoral et des dunes notamment entre le village El Khmara et Châabna où la hauteur de ces dernières atteint parfois plus de 20 m.

Tableau 9 - Les variables les plus contributives (en %) dans la vulnérabilité de la cellule sédimentaire de Salakta-Forêt d’El Ghedhabna (avec rotation Varimax brut)

ID	Variable	Fact 1	Fact 2	Fact 3	Fact 4
7	Espaces touristiques (ha)	25,24	2,33	1,31	1,48
9	Morphologie et géologie	0,01	28,12	7,17	3,42

La deuxième partition des cellules littorales : la partition la plus urbanisée, la plus touristique et la plus vulnérable au risque d’érosion

30Ce groupe se compose de 4 cellules sédimentaires littorales situées toutes dans les deux plus grands foyers touristiques de la côte Kelibia-Chebba. Il s’agit des foyers Nabeul-Hammamet et Sousse-Monastir. Selon leurs coordonnées sur les axes de l’ACP, elles peuvent être subdivisées aussi en deux sous-partitions géographiquement bien différentes. La première regroupe les deux cellules sédimentaires situées au Nord du golfe de Hammamet entre les ports de pêche de Beni Khiyar et de plaisance de Yasmine El Hammamet. La deuxième se localise au Sahel central dans le plus ancien foyer touristique de la côte orientale de la Tunisie entre Sousse et Monastir, et est composée aussi de deux unités sédimentaires qui s’étirent entre le port de plaisance de Kantaoui et la Marina de Monastir.

31Selon le tableau 10, il se dégage au moins 5 variables qui contribuent à la vulnérabilité des littoraux dans les cellules précitées et dont le poids de la variance de chacune est supérieur à 21 %. Ces variables se montrent plus contributives notamment sur le 2^ème axe factoriel de l’ACP où le rythme d’évolution sédimentaire volumique, la morphologie et la géologie des côtes et l’élévation du niveau marin future accumulent 85,6 % du total de la variance.

Tableau 10 - Les variables les plus contributives dans la vulnérabilité du littoral entre Beni Khiyar et le port de plaisance de Yasmine El Hammamet d’une part et entre le port de plaisance de Kantaoui et la Marina de Monastir d’autre part (avec rotation Varimax brut)

ID	Variable	Fact 1	Fact 2	Fact 3	Fact 4
1	Rythme d’évolution sédimentaire volumique (m3/an)	0,04	35,68	0,01	1,09
4	Ports	0,20	4,36	1,18	49,96
7	Espaces touristiques (ha)	25,24	2,33	1,31	1,48
9	Morphologie et géologie	0,01	28,12	7,17	3,42
10	Elévation du niveau marin	1,47	21,76	21,24	1,73

Unité : En pourcent

32Pour être plus explicite, nous analyserons, dans ce qui suit, la contribution de ces différentes variables axe par axe. Ainsi, sur le premier axe, la contribution de la variable « Espaces touristiques » apparait de loin la plus importante avec une variance de 25 %. Les manifestations de cet indicateur se visualisent surtout dans les deux cellules sédimentaires du Sahel central où le littoral est passé d’un espace totalement agraire (oulej) dans les années 1960 à un autre fortement anthropisé, notamment entre le port de Kantaoui et celui de Sousse. Selon le rapport annuel de l’Office National du Tourisme Tunisien de l’année 2010 (année de référence), le nombre total des établissements hôteliers dans cette cellule littorale est de 99 bâtis, s’étend sur une superficie de 147 ha et sont souvent les pieds dans l’eau (tableau 11). Cette importante chaîne d’établissements hôteliers possède un parc d’hébergement de 38 957 lits logeant 1 071 028 de touristes et a permis de commercialiser en 2010 près de 6,9 millions de nuitées (ONTT, 2010). Une telle occupation, surtout celle implantée au contact du rivage, n’a fait qu’entrainer une perturbation de l’échange sédimentaire entre les différentes parties du profil transversal de la plage. Le rythme de retrait de la ligne de rivage, situé entre les hôtels Houda et Prima Life a ainsi atteint 2 m/an entre 1962 et 2023 (figure 13).

Tableau 11 - Nombre d'établissements hôteliers et capacité d’hébergement dans la zone touristique de Kantaoui-Sousse en 2010

Zone touristique	Nombre d'établissements	Nombre de lits
Sousse	67	21 251
Kantaoui	32	17 706
Total	99	38 957

Source : ONTT, 2010

33Pour le deuxième axe factoriel, le nombre des indicateurs contribuant à la vulnérabilité des littoraux est plus important. En effet, nous avons compté 3 variables avec une variance totale de 86 %. Il s’agit surtout :

• • • De la variable caractérisant le bilan sédimentaire volumique (m³/an) avec 35,7 %. Elle se manifeste sur les côtes situées, notamment, au Sud des ports de Kantaoui, de Sousse et le bassin de tranquillisation des eaux de refroidissement de la centrale thermique à Sidi Abdelhamid où les volumes de sédiments arrachés puis transportés par la dérive littorale, pour la période 1962 et 2010, sont énormes avec 295 292 m³ entre le bassin de tranquillisation à Sidi Abdelhamid et l’hôtel Neptunia et 188 056 m³ entre le port de Sousse et le Marabout de Sidi Abdelhamid ;

    • De la variable résumant la morphologie et la géologie de la côte notamment dans la cellule sédimentaire délimitée par le port de Beni Khiyar et la rive gauche de l'Oued Belgaied avec un poids factoriel de 28,1 %. La vulnérabilité de la côte à l’érosion est ici élevée. En effet, elle se voit notamment entre l’hôtel Bel Azur (à Hammamet Nord) et le Marabout Sidi Mahersi (à Nabeul) où dominent de petites falaises devancées par des plages souvent étroites et minces et constituées essentiellement par des formations gréso-marneuses néogènes souvent surmontées par des grès marins quaternaires (photos 3 et 4) ;

• Et enfin de la variable d’élévation du niveau marin qui participe avec 21,8 % dont l’effet est perceptible surtout dans les secteurs d’érosion. D’après la simulation numérique faite selon le scenario d’élévation moyenne de la mer de 50 cm, la cellule sédimentaire située entre le bassin de tranquillisation de Sidi Abdelhamid et la Marina de Monastir sera le lieu d’une perte surfacique remarquable estimée à environ 100 ha à l’horizon 2100 (Bada, 2017a).

34Quant au quatrième axe factoriel de l’ACP, il présente la seule variable « Ports » qui participe avec près de 50 % de la variance totale. En fait, la contribution des ouvrages portuaires est la plus élevée de tous les axes factoriels. Le port de Sousse donne l’une des illustrations les plus remarquables. A cause du blocage du transit sédimentaire fait par sa longue jetée nord qui présente l’emprise la plus longue de tous les ports de la côte Kelibia-Chebba (910 m), un phénomène d’érosion très sévère s’est produit au Sud de cet ouvrage (photos 5a et 5b). La photographie aérienne de 1996 révèle que la plage à Gaied Souassi a fini par disparaitre notamment au Sud de la rive droite de Oued El Hallouf. Cette perte se montre aussi très concordante avec l’érosion de l’avant-plage notamment entre -2 et -8 m de profondeur. Dans cette partie, la barre sableuse interne qui s'accole aux jetées de part et d'autre du port est rompue et les pertes en sédiments sont considérables. Elles sont estimées d’après les données bathymétriques de 1952 et de 2006 du Service Hydrographique et Océanographique de la Marine (SHOM), à 2 252 386 m³. Ces pertes s’élèvent au double jusqu’à -12 m de profondeur, soit 4 582 841 m³ (Bada, 2020).

35En revanche, dans cette deuxième partition de cellules sédimentaires subsistent d’autres indicateurs parfois très contributifs, mais qui sont dissimulées par la rotation factorielle. Il s’agit de 4 variables à savoir « Ouvrages de défense », « Infrastructures hydrauliques », « Espace bâtis » et « Fréquentation des plages » dont 3 résumant, sur le troisième axe factoriel, plus de 75 % de l’information existante (tableau 12).

Tableau 12 - Les variables les plus contributives dans la vulnérabilité du littoral dans les cellules sédimentaires situées entre Beni Khiyar et le port de plaisance de Yasmine El Hammamet d’une part et entre le port de plaisance de Kantaoui et la Marina de Monastir d’autre part (sans rotation).

ID	Variable	Fact 1	Fact 2	Fact 3	Fact 4
2	Ouvrages de défense	12,79	4,72	13,67	6,71
3	Infrastructures hydrauliques	0,28	1,30	41,51	22,45
5	Espaces bâtis (ha)	0,56	1,68	19,88	2,76
8	Fréquentation des plages (estivant/m²)	4,08	0,25	7,96	17,06

36Unité : En pourcent

37Nous exposons ici également les variables précitées dans le tableau 12 par ordre de poids factoriel :

• Les infrastructures hydrauliques qui caractérisent, selon la projection des variables sur le plan factoriel (3^ème et 4^ème axe), les deux cellules sédimentaires situées entre Beni Khiyar et le port de plaisance de Yasmine El Hammamet. Ici, la contribution factorielle est de 41,5 %. Ainsi entre l’Oued El Kbir au Sud de Beni Khiyar et Oued Moussa au Nord du port de plaisance de Yasmine El Hammamet, 11 barrages collinaires ont été construits seulement entre 1996 et 2001. Par conséquent, ces barrages ont dépourvu la côte d’une charge sédimentaire annuelle moyenne très importante estimée entre 10 et 15 Mm³ (DGRE, 2010) ;

• Les espaces bâtis dont les contributions à la vulnérabilité des côtes sont généralisées dans toutes les cellules sédimentaires de cette partition avec un poids factoriel de presque 20 %. Ainsi entre 1962 et 2010, les espaces bâtis sur les côtes de ces cellules sédimentaires ont presque quadruplé en passant d’une superficie totale de 2 362 ha à 8 112 ha. Une telle évolution a été souvent très proche de la mer, voire même au détriment du DPM (Domaine Public Maritime) notamment dans les grandes villes à l’image de Nabeul, Hammamet, Sousse et Monastir (photos 6a et 6b) ;

• Les ouvrages de défense qui contribuent avec 13,7 % : ceci caractérise en particulier les cellules sédimentaires situées entre le port de Kantaoui et celui de Sousse et entre le bassin de tranquillisation à Sidi Abdelhamid et la marina de Monastir. Ici, nous avons mesuré 4 035 m d’enrochement, 5 brise-lames et 8 épis de différentes formes (tableau 13). L’implantation de ces ouvrages, faite pour protéger le littoral contre la dégradation, a participé au contraire à l’accélération du processus de l’érosion de la côte, notamment dans les zones situées du côté qui ne fait pas face à la dérive littorale dominante. C’est le cas, par exemple, des épis installés au Sud du port de Kantaoui en 1990 qui ont emmené à une grave crise sédimentaire.

Tableau 13 - Distribution des talus en enrochement dans les cellules sédimentaires situées entre le port de Kantaoui et celui de Sousse et entre le bassin de tranquillisation à Sidi Abdelhamid et la marina de Monastir

ID	Cellule sédimentaire	Longueur des talus en enrochement (m)	Nombre de brise-lames	Nombre d’épis
7	Port de plaisance de Kantaoui- Port de Sousse	420	2	3
9	Bassin de tranquillisation à Sidi Abdelhamid-Marina de Monastir	3 615	3	5
Total		4 035 m	5	8

• La fréquentation des plages qui caractérisent le 4^ème axe avec une variance de 17,1 % : la vulnérabilité liée à la fréquentation des plages est très perceptible sur le littoral de toutes les cellules sédimentaires de cette deuxième partition. A Sousse, par exemple, nous avons estimé devant les hôtels Môvenpick et Riadh Palms à Bou Jâafar, dans la journée du 21 juillet 2010, entre 2500 et 3000 estivants. Ceci dépasse 3 à 7 fois la capacité d’accueil de la plage devant ces hôtels. Aussi, dans la cellule sédimentaire située entre la rive droite de l'Oued Belgaied et le port de plaisance de Yasmine El Hammamet, la fréquentation des plages par les estivants est également élevée. Ainsi, au cours de la journée du 17 juillet 2010, nous avons compté devant les hôtels les Orangers, Phénicia et Shératon près de 1434 estivants seulement en 4 heures (entre 11 h et 15 h). La capacité de dépassement d’accueil de cette plage a été évaluée entre 4,4 fois devant l’hôtel Shératon et 9,3 fois devant l’hôtel Phénicia.

La troisième partition de cellules sédimentaires : un paysage géomorphologique très alternant

38C’est la plus grande partition qui contient 6 cellules sédimentaires. Selon leurs projections sur le plan factoriel constitué par le 1^er et le 2^ème axe, elle peut être également subdivisée en deux sous-partitions.

• La première sous-partition est composée des cellules Mâamoura-Port Beni Khiyar et Mahdia-Salakta. L’indicateur « Morphologie et géologie » se montre ici le plus constructif sur le deuxième axe factoriel. En effet, il contribue à la vulnérabilité des côtes avec une variance de 21 % notamment dans les secteurs d’érosion qui sont composés souvent par des platiers rocheux ou des petites falaises façonnées dans des grès tyrrhéniens.

39Dans les secteurs d’accumulation où dominent des plages peu larges et des dunes relativement hautes de quelques mètres, les côtes sont aussi très vulnérables à l’élévation du niveau marin. Celle-ci, qui caractérise aussi le deuxième axe de l’ACP, contribue à 18 % de la variance totale. Déjà les signes d’une telle élévation du niveau marin sont perceptibles, notamment sur la côte d’El Mnakâa où des vestiges archéologiques sont aujourd’hui attaqués par la mer (photo 7). En réalité, le démantèlement de ce type de vestige, dû notamment à cette variable, est aussi perceptible ailleurs, sur toute la côte Kelibia-Chebba comme dans le cas également de Hergla (photo 8).

• La deuxième sous-partition est constituée par quatre cellules situées au Sahel central à savoir : Port de Hergla-Port de plaisance de Kantaoui, Port de Sousse-Bassin de tranquillisation de Sidi Abdelhamid, Marina de Monastir-Flèche littorale d'El Enf et Khnis-Mahdia. Selon le plan factoriel constitué par le 1^er et le 2^ème axe, ces cellules se situent très proches de l’origine des axes factoriels. Cependant, avec une rotation des axes factoriels par la méthode de Varimax brut, les variables contributives à la vulnérabilité des littoraux sont mal décelables.

40Pour une meilleure position de ces variables sur le plan factoriel, nous avons revu l’analyse de l’ACP, mais cette fois sans rotation. Les résultats obtenus ont fait ressortir 3 variables positionnées sur le 1^er et le 3^ème axes factoriels, résumant une variance totale relativement faible avec 44,2 % de l’information existante dans le tableau initial.

41Sur le 1^er axe factoriel, nous avons retenu les variables « Ouvrages de défense côtière » et « Morphologie et géologie » résumant chacune un poids factoriel très proche de la moyenne avec respectivement 12,8 et 11,5 %. Ainsi, la contribution de la variable « Ouvrages de défense côtière » est souvent généralisée notamment dans les cellules sédimentaires situées entre la Marina de Monastir et la flèche littorale d'El Enf et entre Khnis et Mahdia. Dans ces deux cellules, les ouvrages surtout de type enrochement sont dispersés, mais souvent sur des longues distances (5 678 m entre Sayada et Rass Eddimess et 1 355 m entre les ports El Karraya et El Ghedir au Sud de la Marina de Monastir). Les méfaits de ces ouvrages sont relatifs puisque la morphologie de la côte est dominée par de petites falaises souvent façonnées dans des formations gréseuses tyrrhéniennes. Même entre Sayada et Rass Eddimess, la côte accorde une place indéniable aux affleurements rocheux qui correspondent fréquemment à des croûtes calcaires.

42Quant à la variable « Espaces bâtis » appartenant au 3^ème axe factoriel, elle tient la variance la plus forte avec 19,9 % notamment dans les secteurs d’érosion où se développent les principales agglomérations urbaines, notamment celles situées entre Monastir et Sayada. L’importance des espaces bâtis explique aussi l’état de vulnérabilité du littoral dans la cellule sédimentaire entre le port de pêche de Hergla et celui de plaisance de Kantaoui qui se voient notamment dans son secteur d’accumulation situé entre l’exutoire de Sabkhet Halk El Mujjel et le port de Kantaoui où les constructions sont souvent pied dans l’eau (figure 16). Ainsi, nous avons mesuré, en 2022, près de 763 ha consommé seulement par les habitations et les résidences contre moins d’un hectare en 1962.

43Dans la même cellule sédimentaire, la variable « Morphologie et géologie de la côte » est pareillement un facteur pouvant expliquer la vulnérabilité du littoral. Ainsi, le tracé de la côte entre Hergla et l’exutoire de Sabkhet Halk El Mujjel est très irrégulier et souvent en falaise façonnée aussi dans un dépôt quaternaire peu résistant à tendre, constitué par des alluvions récentes et actuelles parfois surmontées par des conglomérats et croûtes calcaires.

44Quant à la cellule située entre les ports de Sousse et le bassin de tranquillisation à Sidi Abdelhamid, elle est relativement éloignée de l’origine des axes factoriels et elle présente 4 variables très contributives à la vulnérabilité du littoral. Il s’agit, dans l’ordre du poids factoriel des variables, de : fréquentation des plages (22,6 %), densité de la population côtière (21,1 %), espaces bâtis (19,6 %) et ouvrages de défense côtière (9,2 %). Ces variables caractérisent le premier axe factoriel de l’ACP en résumant une variance totale de 72,5 %.

45Ainsi, la fréquentation de la plage à Gaied Souassi, au Sud de l’Oued El Halouf, par les estivants locaux venant des cités populaires de Bir Echebek, Souiss, Taffala, Ksibet Chatt, 2 Mars, El Matar et Sidi Abdelhamid est très importante bien que la plage montre souvent un état très pollué, étroit et maigre. L’importance de la fréquentation de la plage revient surtout à la taille et à la densité de la population très élevées. Ainsi, entre 2004 et 2022, la population de la délégation de Sidi Abdelhamid est passée de 46 257 à 55 588 habitants, soit une augmentation intercensitaire de plus de 9 331 habitants. Dans la même période, la densité de la population est donc passée de 3 571 h/km² à 4 291 h/km² (tableau 14).

46Avec une telle augmentation de la population, le littoral de Sidi Abdelhamid a connu une consommation sans précédent de son espace au détriment de la Sebkha limitrophe, consommation estimée à 457 ha entre 1962 et 2022, soit un rythme d’évolution annuel de 7,6 ha.

Tableau - Evolution de la densité de la population dans la délégation de Sidi Abdelhamid entre 2004 et 2022

Sidi Abdelhamid	2004		2014		2022
	Population	Densité	Population	Densité	Population	Densité
	46 257	3 571 h/km²	52 787	4 075 h/km²	55 588	4 291 h/km²

INES, 2022

Conclusion

47Dans cette étude, nous avons essayé d’identifier les principaux indicateurs de vulnérabilité sur le littoral tunisien fortement anthropisé entre Kelibia et Chebba et ceci dans chacune des 13 cellules sédimentaires identifiées. Pour repérer les indicateurs a priori les plus pertinents contribuant à la vulnérabilité face au risque d’érosion, une analyse exploratoire multidimensionnelle a été réalisée. Cette technique très utile en analyse des données, dite aussi d'interdépendance, a permis de synthétiser l'information de 10 indicateurs en favorisant 4 principaux facteurs ou axes expliquant 85,6 % de l’inertie totale.

48La classification en 3 partitions a permis d’une part de saisir, presque partout, les rôles très importants et très contraignants que jouent les aménagements dans la vulnérabilité du littoral aux risques d’érosion, surtout au contact des grandes villes et dans l’aire des principales zones touristiques. D’autre part, elle souligne le rôle parfois limité des agents naturels qui n’intéressent que les cellules à côtes sableuses où dominent, surtout, les plaines littorales. L’impact négatif des actions anthropiques va davantage s’amplifier avec la nouvelle conjoncture de l’élévation du niveau moyen de la mer dont le poids le plus important est enregistré notamment au niveau des côtes basses où prédominent les zones humides. Les cellules sédimentaires, surtout celles situées entre Kelibia et Mâamoura et entre le port de plaisance de Yasmine El Hammamet et Hergla, apparaissent comme étant les plus vulnérables notamment dans leurs espaces les plus aménagés.

49La détermination du nombre d’indicateurs naturels et anthropiques, les plus contributifs de vulnérabilité du littoral Kelibia-Chebba face au risque d’érosion ainsi que leurs poids par l’approche ACP, acquiert une grande importance pour la gestion de la dégradation côtière, les résultats obtenus constituent donc une information utile pour être mise à la disposition des décideurs et des administrateurs locaux comme guide pour la gestion intégrée et durable de l'érosion marine.