2

Résumé
Le recensement des différentes structures actives citées par la littérature dans la région du
Nord-Est de la Tunisie a mis en évidence les arguments morphologiques et les indices de
déformation qui attestent d’une activité tectonique très récente. La phase tectonique
compressive majeure post-villafranchienne, de direction NNW-SSE, est responsable de la
mise en place des reliefs récents au Nord-Est (exemple : pli d’Utique).
Des indices certains d’activité tectonique récente ont été découvert au niveau des failles
d’Utique, de Messeftine de Djebel Nahli et de Djebel Ammar. Ces indices permettent de
classer ces failles comme des sources sismogéniques à potentiel important.
L’étude pluridisciplinaire de la faille d’Utique (géologie de surface, géophysique,…), ainsi
que les données archéologiques ont permis d’affirmer que son activité la plus récente date
au plus de l’époque romaine (≈2000 ans B.P). Nous estimons la magnitude du séisme
d’Utique, mentionné par la littérature, à une magnitude de 6 sur l’échelle ouverte de
Richter et sa période de retour à 1300 ans. L’équilibrage de la coupe de la structure
d’Utique à permis d’estimer la quantité de raccourcissement post-Miocène supérieur
accomodée par la structure d’Utique à 275 mètres. En supposant que ce
raccourcissement a lieu depuis la phase compressive post-villafranchienne (1,8 Ma) ;
on propose une vitesse de raccourcissement d’environ ~0,15 mm/an.
Les caractéristiques topographique et lithologique de plusieurs régions de la ville de Tunis
sont très favorables aux effets de sites et au phénomène de liquéfaction. Pour faire face aux
risque du rejeu de la faille d’Utique, des précautions spéciales et des règles de construction
parasismique doivent être prises en comptes lors la réalisation de tous les ouvrages d’art
aux alentour de la source seismogénique d’Utique.

Abstract
The study of active structures mentioned in the literature in the region of the North-east of
Tunisia show the morphological evidence of deformation which attest to a very recent
tectonic activity. The late-villafranchian major compressive NNW-SSE tectonic phase is
responsible for the establishment of the recent reliefs in the North-east (example: Utique
fold).

3
Some indices of recent tectonic activity were discovered on active faults on the structures
of Utique, Messeftine, Djebel Nahli and Jebel Ammar. These indices allow classifying
these faults as sources with significant seismogenic potential.
Geologic, seismic, and neotectonic investigations in Utique fault area; as well as the
archaeological data shows that the most recent activity was during the Roman era (≈2000
years B.P). We attribute the magnitude of 6 to the destructive Utique earthquake,
Mentioned in the literature, and we evaluate its return period to 1300 years. Cross section
shows a total shortening of 275 m since the Upper Miocene. Compression may have begun
early in Pleistocene (1.8 My); this leads to a most likely value for shortening rate of 0,15
mm/yr.
The topographic and lithological characteristics of several regions of in Tunis City are very
propitious to the sites effects and to liquefaction phenomenon. This high seismic risk zone
deserves to be taken into account during the establishment of important regional
development programs and in the application of seismic building codes.

4
 SOMMAIRE

INTRODUCTION ............................................................................................................. 7

Chapitre I: La Tunisie dans le cadre géodynamique

A : Tectonique et sismicité en Méditerranée ................................................................. 13

A-1 : évolution géodynamique de la Méditerranée .................................................... 13
A-2 : caractéristiques des Maghrébides ..................................................................... 18
A-3 : tectonique récente et sismicité.......................................................................... 21
A-4 : cinématique de la Méditerranée ....................................................................... 22
A-4 –1 : mesures des champs de contraintes (GPS) ............................................... 22
A-4-2 : implication sur le déplacement actuel en Tunisie ....................................... 27
B : Tectonique récente et sismicité en Tunisie .............................................................. 28
B-1 : Tectonique récente en Tunisie .......................................................................... 28
B-2 : Sismicité : Caractères et problématique ............................................................ 29
B-2-1 : Sismicité historique (figure 10) ................................................................. 31
B-2-2 : sismicité instrumentale .............................................................................. 34
B-2-3 : Complexité de cette sismicité intraplaque .................................................. 37

Chapitre II: Géodynamique et tectonique récente du Nord-Est de la Tunisie

A : Description lithostratigraphique .............................................................................. 40

A-1 : Le Néogène ..................................................................................................... 40
A-1-1 : Miocène supérieur .................................................................................... 40
A-1-2 : Pliocène .................................................................................................... 45
A-2 : Le Quaternaire ................................................................................................. 47
A-2-1 : Pléistocène inférieur (Villafranchien) ....................................................... 47
A-2-2 : Pléistocène moyen .................................................................................... 48
A-2-3 : Pléistocène supérieur ............................................................................... 49
A-2-4 : Le Tyrrhénien ........................................................................................... 51
A-2-5 :L’Holocène ................................................................................................ 52

5
 B : Tectonqiue de la Région du Nord-Est de la Tunisie ................................................. 53
B-1 : Tectonique anté-néogène ................................................................................. 53
B-2 : Tectonique Néogène ........................................................................................ 56
B-3 : Tectonique Quaternaire .................................................................................... 56

Chapitre III: Tectonique active de la region d’Utique et ses environs

Introduction.............................................................................................................. 5859
A : Les principaux séismes dans le Nord de la Tunisie (après JC) ............................. 5859
A-1 : Séisme historique d’Utique .......................................................................... 5859
A-2 : Séisme historique de Tunis .......................................................................... 5960
B : Sismicité instrumentale à effet de surface (depuis 1900 après JC) ...................... 5960
B-1 : Séisme du 1er décembre 1970 ....................................................................... 5960
B-2 : Séisme de Sidi Thabet (12/12/1970)............................................................. 6061
B-3 : Séisme de Ghar el Melh (9/04/1979) ............................................................ 6061
B-4 : Séisme de Chemtou (8/06/1985) .................................................................. 6162
C : Etude des Mécanismes aux foyers ....................................................................... 6263
C-1 : Le concept de mécanisme au foyer ............................................................... 6263
C-2 : Mécanismes aux foyers : champs de contraintes au Nord de la Tunisie ........ 6465
D : Aperçu sur les principales structures actives au NE de la Tunisie ........................ 7172
D-1 : Etat de connaissances des structures actives (exemples de littérature) .......... 7172
D-1-1 : L’accident NW du Djebel Nahli ............................................................ 7475
D-1-2 : La Faille d’El Alia ................................................................................ 7576
D-1-3 : Déformation tectonique Quaternaire récent de Ras Engla ...................... 7576
D-1-4 : La faille inverse d’El Fejja .................................................................... 7879
D-2 : Nouvelles observations de la Tectonique récente au Nord-Est de la Tunisie . 8081
D-2-1 : La faille active de Djebel Ammar .......................................................... 8081
D-2-3 : Tectonique Post-Tyrrhénienne affectant la falaise tyrrhénienne de Sounine
............................................................................................................................. 8889
D-3-2 : Tectonique Post-Tyrrhénienne affectant la Plage d’Errimel .................. 9091
D-4 : Anticlinal de Messeftine .............................................................................. 9293
D-4-1 : Géologie ............................................................................................... 9293
D-4-2 : La faille de Messeftine .......................................................................... 9495

6
 D-4-3 : Indices de tectonique actives ................................................................. 9596
a/ indice de Oued El Meleh ............................................................................ 9596
b/ indice de Djebel El Assafir ........................................................................ 9596
c/ indice de la Ferme Lebray .......................................................................... 9798
D-4-4 : Cartographie du Quaternaire sur le versant Ouest de l’anticlinal ............ 9798

Chapitre IV: Etude de la faille active d’utique

A : Cadre géomorphologique, paléogéographique et géologique de la basse Medejrda et

la région d’Utique ................................................................................................ 105106
A-1 : Cadre géomorphologique et géologique ................................................... 105106
A-2 : Cadre paléogéographique......................................................................... 108109
A-3 : Géomorphologie de l’anticlinal d’Utique ................................................. 110111
A-4 : Géologie de l’anticlinal d’Utique ............................................................. 112113
B : Etude de Section sismique .................................................................................... 114
B-1 : Coupe géologique équilibrée .......................................................................... 115
B-2 : Equilibrage de la section géologique à travers la structure d’Utique ............... 116
C : Les données de la tectonique quaternaire au niveau de la structure d’Utique ......... 119
D : Application d’une approche radiochronologique ................................................... 121
E : Mesure de l’inclinaison du niveau de silts à poupées calcaires ............................. 124
F : description des différents sites d’observation au niveau de la faille d’Utique ......... 126
G : Description de la faille d’Utique : ......................................................................... 135

Chapitre V: Implication sismotectonique de la faille d’utique sur ses environs

A : Introduction: Sismotectonique et sismicité intraplaque.......................................... 141

B : Etudes sismotectonique en Tunisie : historique ..................................................... 143
C : Alea Sismique : Application au cas de la Tunisie .................................................. 144
C-1 : Loi de récurrence de Gutenberg-Richter : log N = a – bm .............................. 146
C-2 : Caractérisation de l’atténuation du mouvement du sol .................................... 148
C-3 : Caractérisation du mouvement sismique : ................................................ 149148
D : Implication d’un éventuel séisme de la faille d’Utique sur la zone du grand Tunis :
............................................................................................................................. 154153

7
 D-1 : Effet de site ................................................................................................. 154153
D-1-2 : Susceptibilité à la liquéfaction et profondeur de la ligne piézométrique :
..................................................................................................................... 158157
D-1-3 : Effets de site dus aux caractéristiques mécaniques des formations
superficielles et de la géométrie des interfaces entre les formations géologiques:
..................................................................................................................... 161160
D-2 : Récurence de la faille active d’Utique ......................................................... 163162

Conclusion et Perspectives ....................................................................................... 165164

Liste Bibliographique .............................................................................................. 167166

8
 Liste des Figures

Figure 1: Carte de la Sismicité historique et instrumentale de la région d’étude

Figure 2 : Cadre géodynamique de la méditerranée occidentale
Figure 3: Reconstitution paléogéographique de l'extension d'arrière-arc au niveau de la
Méditerranée Centrale.
Figure 4: Evolution de la marge nord africaine
Figure 5: Les principaux domaines structuraux de l’Afrique du Nord
Figure 6 : Coupes géologiques équilibrées montrant les déformations tectoniques dans le
domaine atlasique tunisien.
Figure 7: Vitesses GPS mesurées par rapport à l'Eurasie.
Figure 8: Modèles de déplacement de l’Afrique par rapport à l’Eurasie
Figure 9: Modèle géodésiques de déplacement de l’Afrique par rapport à l’Eurasie
Figure 10: Carte de la cinématique des plaques africaine et eurasiatiques
Figure 11: Séismicité du bassin de la Méditerranée (IPSN-BERSSIN)
Figure 12 : Carte de la Sismicité historique de la Tunisie
Figure 13: Configuration du réseau sismologique de la Tunisie
Figure 14: Carte de la sismicité instrumentale de la Tunisie
Figure 15: Carte de Répartition des épicentres du Nord de la Tunisie
Figure 16: Log stratigraphique du puits Utique1 (Rapport ETAP)
Figure 17: Corrélation entre les forages pétroliers disponibles dans la région d’étude
Figure 18: Coupe à Djebel Menzel Ghoul au niveau de flanc Sud de l'anticlinal d'Utique
Figure 19: Coupe dans le Quaternaire au niveau de Menzel (Oueslati, 1993)
Figure 20: les différentes subdivisions du Tyrrhénien en Tunisie
Figure 21: Carte d'isochrone du toit du Trias (Equi. = 100 ms) (Kacem, 2004)
Figure 22: Carte d'isochrone du toit du Crétacé-Jurassique (Equi. = 100 ms)
Figure 23: Carte tectonique de subsurface de la région du Nord-Est de pays
Figure 24: Cratère de sable observé à Sidi Thabet après le Séisme du 12/12/1970
Figure 25: Localisation géographique de la ville de Ghar El Melh et vue panoramique du vieux
port de la ville
Figure 26: Fracturation de la surface suite au séisme de Chemtou (08-06-1985)
Figure 27: Représentation stéréographique des mécanismes au foyer
Figure 28: Carte des principaux mécanismes au foyer du Nord de la Tunisie et la région Nord
orientale de l’Algérie
Figure 29: Solution focale de quelques séismes au Nord-Est du pays
Figure 30: Mécanismes au foyer des séismes du 19/02/1993 et 26/01/1995
Figure 31: Mécanisme au foyer du séisme de 12/08/2001
Figure 32: Mécanisme au foyer du séisme du 24/06/2002
Figure 33: Modèle numérique de terrain avec les principaux reliefs de Tunis et du NE du pays
Figure 34: Carte sismotectonique du NE de la Tunisie
Figure 35: Carte géologique de la structure de Djebel Nahli
Figure 36: Coupe de la falaise à l’Est du Phare de Ras Engla
Figure 37: Localisation géographique de la faille d'el Fejja
Figure 38: Faille Pliocène inverse au niveau d’El Fejja (Ouest de l’anticlinal d’Utique).
Figure 39: Failles normales au sein des dépôts Quaternaires sur le pied du versant Nord de
l'anticlinal de Djebel Ammar
Figure 40: Failles Pléistocène supérieur révélées par une tranchée sur le flanc Nord de Djebel
Ammar
Figure 41: Faille normale au Pléistocène sup. avec le réseau de failles secondaires associé
Figure 42: La faille N-S d’el Menzeh9 et de la colline de Sidi Bel Hassen
Figure 43: Faille quaternaire affectant le Pliocène et la base du quaternaire.
Figure 44: Localisation des indices de tectonique active au niveau de la corniche tyrrhénienne
Figure 45: Une faille inverse affectant les dépôts Tyrrhénien de la plage de Sounine.
Figure 46: Décrochement sénestre dans les dépôts tyrrhéniennes de la falaise de Sounine.
Figure 47: Photo de la corniche Tyrrhénienne du Plage Errimel. On voit le flanc sud d’un
anticlinal Tyrrhénien à large rayon de courbure
Figure 48: Coupe simplifié de l’anticlinal d'Errimel. Elle montre le plissement des dépôts marins du
Quaternaire récent
Figure 49: Photo montrant des failles E-W qui jalonne les couches tyrrhéniennes plissées de la
corniche d'Errimel
Figure 50: Vue panoramique de l'anticlinal de Messeftine (image de Google earth)
Figure 51: Carte géologique de l'anticlinal de Messeftine
Figure 52: Discordance des silts Quaternaires sur les argiles verdâtre du Miocène (vue vers le
NE)
Figure 53: Failles normales synsédimentaires pliocènes au niveau de Djebel El Assafir
Figure 54:Contact tectonique entre la croute calcaire Quaternaire et les formations Miocène
Figure 55: Cartographie de Djebel Messeftine et de la faille de Messeftine
Figure 56: Schéma de répartition des indices de la déformation récente et les mécanismes au
foyer dont ils sont associés en Tunisie Nord orientale
Figure 57: Les indices de néotectonique au niveau de la région de Tunis et du NE de la Tunisie
Figure 58: Vue panoramique de l'anticlinal d'Utique (vue vers le Nord)

2
Figure 59: Extrait de la carte topographique (1/25000 ; équid. 5m) qui montre la structure
d'Utique au milieu de la plaine d'El Mabtouha
Figure 60: Les étapes de comblement de l’ancien golfe d’Utique (Jauzein,1971)
Figure 61: MNT de l’anticlinal d’Utique et des profils topographiques transversaux à la
structure
Figure 62 : Modèle schématiques de (a) plis de décollement et (b) plis de propagation de faille
Figure 63: Coupe des formations quaternaires sur le versant Sud de l'anticlinal d'Utique.
Figure 64: Profil L1 (mission MA83) passant la structure d’Utique
Figure 65: Exemple d’équilibrage d’une coupe géologique
Figure 66: Interprétation du profil L1 et estimation du rejet de la faille d'Utique
Figure 67: Faille inverse qui affecte le Plio-quaternaire au niveau du versant Sud du pli
d'Utique
Figure 68: Failles normales associées au plissement quaternaire de la série plio-quaternaire au
niveau de la crête de l'anticlinal d'Utique
Figure 69: Diagramme Concordia pour la résolution des âges de manière graphique en
calculant les rapports isotopiques (Kaufman et al., 1971)
Figure 70: Poupées calcaires dans la formation du silt rouge du quaternaire au niveau du
piémont du versant Sud de l'anticlinal d'Utique
Figure 71: Stéréogramme (projection de Wulf) de l’orientation des concrétions calcaires à 2
endroits sur le flanc sud-est du pli.
Figure 72: Localisation des sites d'indices d'activité de la faille d'Utique (S1....S6).
Figure 73: Affleurement du plan de la faille d'Utique
Figure 74: Jeu inverse de la faille d'Utique: contact tectonique des grès Pliocènes sur les silts
brunâtres du Quaternaire
Figure 75: Coupe géologique qui montre le plan de la faille inverse d'Utique.
Figure 76: La pièce de poterie imbriquée dans la zone de broyage de la faille d'Utique
Figure 77: Vue des deux facettes de la pièce de tesson trouvée dans la zone de broyage de la
faille inverse d'Utique
Figure 78: Coupe qui montre la zone de passage de la faille d'Utique
Figure 79: Le plan de la faille d'Utique séparant les bancs de grès Pliocène et les silts
rougeâtres du Quaternaire
Figure 80: Contact anormal entre la croûte calcaire Quaternaire et les bancs gréseux du
Pliocène
Figure 81: Les fragments de poterie de l'époque Romain encaissés dans les silts brunâtres du
Pléistocène supérieur
Figure 82: Le pendage élevé du quaternaire observé sur le flanc Sud de l'anticlinal d'Utique

3
Figure 83: Illustration schématique de la faille inverse d'Utique
Figure 84: Cartographie des facies Quaternaire sur le flanc Sud de l’anticlinal d’Utique
Figure 85: Caractérisation des sources seismogénique et propagations des ondes sismiques
Figure 86: Caractéristiques des différentes méthodes utilisées en sismotectonique (Philip et al.,
2007)
Figure 87: Etapes des évaluations de l'aléa sismique par l’approche probabilistique
Figure 88: Loi de récurrence de Gutenberg-Richter pour les séismes de magnitude supérieure
ou égale à 2
Figure 89: Loi d'atténuation, applicable pour le nord-est de la Tunisie, calculée pour des
profondeurs différentes de l’hypocentre (Z=2, Z=5, Z=10, Z=15)
Figure 90: Cartes d'iso-accélération maximale du sol (Peak Ground Acceleration) pour
une période retour de 100ans (A) et de 475 ans (B) (Kacem, 2004)
Figure 91: Carte des affleurements du grand Tunis et ses environs
Figure 92: Carte piézométrique de la nappe d'eau au niveau de la ville de Tunis (El May et al.,
2009)
Figure 93: Relation entre la distance épicentrale limite des sites où la liquéfaction a été
observée, et la magnitude des séismes superficiels
Figure 94: Effet du "Fault-Slip-Rate" et de la magnitude du séisme sur la période du retour

4
 AVANT PROPOS

Ce travail s’inscrit dans le cadre d’une Thèse de Doctorat en Cotutelle entre

l’Université de Tunis El Manar (Tunisie) et l’Université Paul-Sabatier de Toulouse
(France). Cette collaboration interuniversitaire a été facilitée par deux projets du Comité
Mixte de Coopération Universitaire (CMCU).
Je remercie tout d’abord les membres du jury : Olivier BELLIER et Noureddine
BOUKADI pour avoir accepté de juger ce travail et d’en être les rapporteurs, ainsi que
Joseph MARTINOD et Mohamed GHANMI pour leur lecture attentive, leurs remarques
constructives et leur participation au jury lors de la soutenance. Cette thèse n’aurait pu se
faire sans la confiance que m’ont accordée Mahmoud DLALA et Vincent REGARD en
acceptant d’être mes directeurs de thèse. Je voudrais les remercier plus particulièrement
pour leur enthousiasme et leur optimisme tout au long de cette thèse.
Merci à Sébastien CARRETIER qui m’a beaucoup aidé dans les sorties de terrain et
qui m’a lancé sur la voie de l’utilisation du SIG et des images satellites. Merci pour toutes
ces discussions enrichissantes et constructives.
Mes remerciements vont aussi à tous les gens qui forment le GET : les chercheurs, les
thésards, les informaticiens et ingénieurs, les secrétaires qui m’ont assuré un
environnement sain et propre pour bien travailler, avec tous ces gens là, j’ai passé
d’agréables moments. Merci à toutes et à tous pour l’environnement humain et
scientifique.
Je remercie vivement Stéphane BRUSSET, professeur et chercheur au GET, pour son
collaboration scientifique et son aide précieux dans le cadre des analyses géophysiques et
la réalisation de la coupe géologique équilibrée.
Merci également à Mohsen RABHI pour toutes les conversations animées et pour son
accueil chaleureux dans son bureau et pour m’avoir aidé dans l’interprétation de quelques
observations de terrain.
Je remercie vivement mes collègues du CNSTN pour leur soutien moral et leur
encouragement durant les années de thèse malgré les lourds combats que nous avons mené
ensemble au sein du bureau syndical du CNSTN.

5
 Mes gratifications s’adressent en particulier à Mohamed KHOUATMIA pour m’avoir
aidé à réaliser les analyses isotopiques au sein des laboratoires du CNSTN.
Je remercie énormément ma femme pour son soutien, ses sacrifices et ses motivations
incessantes. Je tiens à exprimer mes sentiments de gratitude à mes enfants pour avoir
supporté mon indisponibilité durant la période de thèse.
Enfin et bien sûr un grand merci à toute ma famille, mon père, ma mère, mes frères
et sœurs qui n’ont pas cessé de m’encourager tout au long de cette période.

6
INTRODUCTION

De par sa position géographique, le territoire Tunisien représente l’extrémité orientale

du domaine maghrébin. Il occupe une position géologique privilégiée en Afrique,
puisqu’il se développe à cheval entre un domaine cratonique saharien, peu déformé, et
un domaine alpin septentrional structuré pour l’essentiel au cours du Méso-
Cénozoïque. La limite entre ces deux domaines est soulignée par un grand alignement
de failles et de flexures, qui se relayent sur plus de 2000 km, depuis Agadir, au Maroc,
jusqu’au golfe de Gabes, en Tunisie (appelé souvent accident sud-Atlasique)
Le domaine des maghrébides est formé par deux principales chaînes Atlasiques et
Telliennes. Ce domaine est constitué de 4 unités majeures dans sa partie occidentale
(au Maroc : l’Anti-Atlas, le Haut Atlas, le Moyen Atlas et le Rif. En Algérie
occidentale, il est constitué de deux chaînes importantes : l’Atlas Tellien et l’Atlas
Saharien, séparés par les Haut Plateaux, qui se fusionnent vers l’Est pour donner une
chaîne unique qui constitue successivement l’Atlas Orientale algérien et l’Atlas
Tunisien.
En Tunisie ; l’Atlas tunisien central constitue la particularité structurale de la chaine
des Maghrébides, qui provient entre autre, de la présence des plis NE-SW disposés
orthogonalement à des fossés d’effondrement NW-SE, et montrant une seule
terminaison périanticlinale (Dlala, 1984; Philip et al., 1986; Dlala & Ben Ayed, 1988).
Cette disposition particulière est liée à l’antériorité de ces fossés par rapport à
l’événement tectonique plicative majeure, qui est à l’origine de la structuration de la
plus grande partie du territoire tunisien. Malgré cet évènement plicatif majeur, la
tectonique compressive ne semble pas atteindre la Tunisie méridionale(Dlala & Ben
Ayed, 1988; Dlala, 1995).
En Tunisie Nord-Est, le bassin de Kechabta-Messeftine et ses environs est connu par
une tectonique compressive post-Villafranchienne très importante. La déformation
récente à été largement décrite par plusieurs auteurs(Ben Ayed et al., 1978; Ben Ayed,
1986; Dlala, 1994; Dlala, 1995; Kacem, 2004; Dlala & Kacem, 2008).
Le premier objectif de ce travail consiste à enrichir les données sur la tectonique active
dans cette région de la Tunisie

7
En Tunisie, la tectonique néogène est connue par une déformation compressive, dont
l’épisode le plus important, qualifié de phase atlasique, a eu lieu au Miocène supérieur.
Elle est responsable des plissements importants qui ont affecté l’ensemble du Nord de
la Tunisie et remis en mouvement le diapirisme triasique. Dans la région du Nord-est
de la Tunisie, le Néogène dans son ensemble est déformé en plis (anticlinal de
Messeftine, anticlinal d’Utique-Besbessia…).
Des études sismotectonique (Hfaiedh, 1983; Hfaiedh & Chadi, 1985; Ben Ayed, 1986;
Dlala, 1992; Gueddiche, 1992; Hfaiedh et al., 1992; Oueslati, 1993; Dlala & Rebai,
1994; Boutib et al., 1997; Kacem, 2004; Dlala & Kacem, 2008) montrent que la région
du Nord Est de la Tunisie est le siège d’activité tectonique récente. Cette tectonique est
également confirmée par plusieurs séismes (séisme de 1970 à Tunis, séisme de 1979 à
Ras Djebel).
Au niveau des régions intraplaque à faibles sismicité et taux de déformation lent,
l’identification des failles actives et leurs relations avec les séismes est un problème
majeur. En effet, la détermination des paramètres sismiques d’une faille (période de
retour des séismes de forte magnitude et l’évaluation de l’aléa sismique) ne peut
s’élaborer sans une analyse complète des déformations les plus récentes le long des
failles actives. Dans ce travail nous avons combinés les études en morohologie des
escarpements de faille, en géophysique, en datation et en investigation de terrain et de
tranchées, afin de
- montrer l’évidence des indices d’activité tectonique le long des failles
quaternaire.
- Quantifier la déformation au niveau des faille actives (exemple d’étude : la faille
d’Utique)
- Estimer le potentiel seismogénique des failles active et les éventuels risques sur
son environnement.
Les données géologiques, géophysiques, sismologiques disponibles sur une région
permettent d’effectuer des études probabilistiques d’alea sismique (Youngs et
Coppersmith, 1986 ; Mc Guir, 1976). Dans ce travail, nous avons utilisé la méthode
probabilistique pour estimer l’aléa sismique aux alentours de la région d’Utique. Cette
méthode est fondée sur la détermination des lieux d’occurrence des séismes, leurs tailles et
leur taux de récurrence. Puis, en un site donné, l’aléa est calculé en sommant les

8
contributions de tous les scénarios possibles (toutes les combinaisons de magnitudes et de
distances). Donc, trois étapes sont nécessaires :
- Identifier les zones sources dans la région d’étude.
- Modéliser dans chaque zone source une courbe de récurrence à partir du catalogue de
sismicité.
- Choisir une relation d’atténuation du mouvement du sol adaptée à la région d’étude.
Cependant, des travaux récents (Newman et al., 2001 ; Cramer, 2001) soulignent les
difficultés liées à cette méthode, principalement dans les régions intraplaques. Elles
concernent le manque de précision sur la récurrence des séismes majeurs, et sur les
modèles de failles utilisées (interaction des failles ou non). Enfin, cette méthodologie est
fondée sur la stationnarité de la sismicité et ne prend pas en compte l’hypothèse d’un
séisme intraplaque dans une zone non identifiée par un séisme passé instrumental ou
historique. Dans ce cadre, la compréhension de l’occurrence des séismes intraplaques
prend tout son sens.

La région de Tunis représente la plus importante agglomération en Tunisie avec ses 2

millions d’habitants. Elle coïncide avec la terminaison nord-orientale de l’Atlas
tunisien, et est considérée comme le point de convergence des principales zones
structurales de la Tunisie (axe N-S, faille de Zaghouan, Atlas septentrionale). Elle est
située sur un bassin sédimentaire subsident qui est formé essentiellement par des
formations argilo-limoneuses. Elle est entourée par des anticlinaux généralement
bordés par des failles. Certaines de ces failles sont considérées comme actives.
Comme le montre la carte de sismicité historique et instrumentale (figure 1), plusieurs
tremblements de terre destructeurs, de forte intensité, (I≥VII), se sont produit dans la
région du Nord Est. Ainsi par exemple, un séisme d’Intensité X à secoué la ville de
Tunis en 856 et à provoqué 50000 victimes (Hfaiedh & Chadi, 1985; Ben Ayed, 1986;
Vogt, 1993), un séisme d’intensité X à frappé la ville de Kairouan au centre tunisien en
854 et un séisme destructeur d’intensité X à provoqué la ruine totale de la ville
historique d’Utique située au Nord de Tunisie. Ce dernier séisme pourrait bien être
associé à la faille inverse qui passe au pied du relief (faille d’Utique).
Malgré le fait que Tunis soit une de civilisation très ancienne (plus que 3000 ans), les
documents concernant les séismes qui se sont produit dans cette région au cours de la

9
période historique sont rares et ne suffisent pas à étudier l’aspect paléosismologique de
la région et déterminer le risque sismique que présentent les structures actives au Nord
Est de la Tunisie. Pour ce faire nous allons explorer l’aspect géologique et
morphostructural des structures les plus récentes (les structures associés au
Quaternaire) pour mettre en évidence les indices de séismes anciens. Cette approche a
pour but de déterminer les structures sismiques et les caractériser en matière d’activité
sismotectonique.
Ce travail s’intègre dans le cadre de l’éclaircissement du contexte sismotectonique des
régions à haute importance socio-économique du pays. Ici on s’intéresse à la région
d’Utique situé au Nord du Grand Tunis. En effet, le premier objectif est d’enrichir le
schéma de la tectonique active dans la région de Tunis (Dlala & Kacem, 2008), et
d’apporter une nouvelle interprétation basée sur l’étude des paléoséismes en exploitant
les nouvelles méthodes de datations des évènements sismiques.
La deuxième partie s’intéresse à l’impact des paléoséismes sur la zone de grand Tunis
par l’étude du potentiel de quelques failles actives susceptibles d’être des sources
seismogénique importantes : dans ce travail, nous discuterons l’activité de la faille
d’Utique et de Messeftine qui pourraient être les sources seismogénique les plus
importantes dans la zone. Cette approche se base à la fois sur l’aspect morphostructural
pour déterminer la géométrie des structures, la cinématique et la vitesse des failles
actives, et paléosismologique à fin de déduire les liens entre les séismes historiques et
les structures actives ainsi que l’estimation de la récurrence des séismes.

10
Figure 1: Carte de la Sismicité historique et instrumentale de la région d’étude
(« M » désigne la magnitude et « I » désigne l’Intensité)

11
 CHAPITRE I

La Tunisie dans le cadre géodynamique

méditerranéen

12
A : Tectonique et sismicité en Méditerranée

A-1 : évolution géodynamique de la Méditerranée

Le contexte géodynamique des zones méditerranéennes est directement influencé par la

convergence Nord–Sud (à Nord-Ouest-Sud-Est) entre l’Afrique et l’Eurasie, qui a
débuté au Crétacé supérieur, et s’est amplifié et généralisé au Crétacé terminal –
Éocène à la suite de l’ouverture de l’Atlantique nord (Grimaud et al., 1982; Olivet et
al., 1983). Pendant ce temps, un régime compressif généralisé s’installe, avec des
déformations consécutives à la fermeture du bassin océanique qui séparait les blocs
continentaux ibérique et bético-rifain. A cette étape on assiste à l’ouverture de l’arrière
arc provoquée par l’accélération de la subduction (Facenna et al., 1997).
L’évolution géodynamique de la marge nord-africaine est étroitement liée aux stades
d’évolution du bassin océanique septentrional. Ce bassin est passé par plusieurs stades
d’évolution qui sont successivement la Périthétys, la Thétys-Ligure, la Mésogée et
enfin la Méditerranée. La phase de l’évolution de la Mésogée en Méditerranée a connu
la mise en place des orogènes nord-africains (Chaînes Atlasiques).
La convergence active entre l’Afrique et l’Eurasie est accommodée, soit par la
subduction en Méditerranée orientale (Wortel et al., 1990) soit, en Méditerranée
occidentale, le long d’une limite de plaques complexe, où les failles inverses
d’orientation ENE–WSW sont les structures dominantes (Meghraoui & Doumaz, 1996;
Morel & Meghraoui, 1996; Boudiaf et al., 1999; Domzig et al., 2006). Le long et au
Sud de cette frontière de plaques, le domaine orogénique maghrébin, d’âge cénozoïque,
comporte le Tell–Rif et les Atlas (Durand-Delga & Fonboté, 1980). Le Tell–Rif est
interprété classiquement comme une chaîne de type alpin, c’est-à-dire résultant de la
fermeture de la Téthys maghrébine (Domzig et al., 2006) ; les Atlas, à l’inverse, sont
compris comme des orogènes intracontinentaux (Mattauer et al., 1977).
Outre la chaîne atlasique, cette orogenèse cénozoïque a mis en place d’autres systèmes
orogéniques dans les régions Maghrébines: le Tell, au sein duquel les zones internes
constituent un domaine fortement allochtone accolé à l’Tunisie au Miocène moyen.
Depuis cette période, la convergence continue et s’exprime en Tunisie par des phases
de compression NW–SE, associées parfois à une extension NE-SW (Dlala, 1994;

13
Dlala, 1995; Dlala & Kacem, 2008). Ainsi, dans l’Atlas tunisien, on observe des plis,
d’âge Plio-Quaternaire, résultant d’une direction de raccourcissement NW–SE.
Au nord de la Méditerranée plusieurs orogenèses sont liées à cette convergence : Les
Apennins, les Pyrénées, les Alpes et la Chaîne Ibérique (figure 2). Cette orogenèse
atlasique a été mise en place par le plan de subduction au niveau de la marge sud
européenne pendant l’Eocène-Oligocène. (Frizon De Lamotte et al., 2000;
Meulenkamp & Sissingh, 2003). Pendant cette période un nouveau bassin d’arrière arc
a été mis en place. Ce bassin qui constitue le Bassin Algéro-provençal actuel, a entrainé
avec lui des îlots détachés (AlKAPeca) de la plaque eurasiatique. Vers le Miocène
supérieur (début du Néogène) La fermeture océanique a entrainé la collision des blocs
ALKAPECA avec la plaque africaine (Tricart et al., 1994; Carminati et al., 1998;
Devoti et al., 2001; Mascle et al., 2001). Cette collision se poursuit jusqu’au
Quaternaire en entrainant la mise en place de la chaîne Atlasique par l’inversion des
principaux bassins sédimentaires nord africain (Figure 3 et 4). La vitesse de l’Afrique
vers le Nord a été estimée à 1 cm/an de l’Oligocène jusqu’au Miocène inférieur
(Carminati et al., 1998). Après la collision, au niveau de la marge tunisienne, le
mouvement de l’Afrique change de direction pour s’orienter vers le NW et ralentir pour
atteindre environ 0,6 cm/an (Carminati et al., 1998). Ce changement de trajectoire
depuis une translation vers une rotation a été mis en évidence au niveau du Golfe de
Cadiz à l’Est du détroit de Gibraltar sur la faille transformante qui délimite les plaques
Africaine et Ibérique (Jiménez et al., 1999).
Le résultat de cette convergence se matérialise par la mise en place de chaîne plissée
tout au long de l’Afrique du Nord : la Chaîne Atlasique.

14
15
Figure 2: cadre géodynamique de la méditerranée occidentale (Doglioni et al., 1999)
16
Figure 3: Reconstitution paléogéographique de l'extension d'arrière-arc au niveau de la Méditerranée Centrale (d'après (Facenna et al., 1997)
 Figure 4 : Evolution de la marge nord africaine (Maury et al., 2000)
1: magmatisme calcoalcalin ; 2: matériaux crustaux déformés (échelle dilatée) ; 3: croûte
continentale européenne ; 4: croûte continentale kabyle ; 5: croûte continentale africaine ; 6:
manteau lithosphérique métasomatisé ; 7: manteau lithosphérique ; 8 : asthénosphère ; 9: flux
asthénosphérique.

17
 A-2 : caractéristiques des Maghrébides

La partie occidentale de l’Atlas nord-africain (Maroc) comporte 4 unités majeures :

L’Anti-Atlas : c’est une chaîne de direction SW-NE et qui s’étend sur environ 600km.
Le Haut Atlas : c’est une chaîne intracontinentale qui s’étend de la ville de Bou-Arfa au
Maroc jusqu’à la frontière maroco-algérienne. Elle change d’orientation de WSW-ENE
à E-W au centre à WNW-ESE dans sa partie orientale (Ahmadi et al., 2006).
Le moyen Atlas est un massif montagneux orienté SW-NE, s’étendant sur 350 km et
situé au centre du Maroc. Il est composée essentiellement de couches calcaires,
tabulaires au nord-ouest et ondulées au sud-est.
Le Rif est une chaîne de montagnes qui s’étend sur tout le Nord du Maroc de la pointe
Est de Tanger jusqu’aux territoires de la tribu rifaine Ikabdanane (Kebdana), à la
frontière algérienne.
Ensuite, vers l’Algérie occidentale l’Atlas est constitué par deux chaînes importantes :
l’Atlas Tellien et l’Atlas Saharien. Ces deux chaines sont séparées par les Hauts
Plateaux. L’Atlas tellien est un relief montagneux principalement situé au nord de
l’Algérie qui se poursuit dans la Tunisie de l’Ouest. Il est formé d’une chaîne de
montagnes longeant le littoral et s’étirant sur environ 1500 kilomètres.
L’Atlas saharien est une zone bien individualisée et structurée depuis le Crétacé
inférieur. Elle se trouve entre le plateforme saharienne relativement stable au Sud et
l’Atlas tellien affecté par une tectonique plus au moins intense au nord (Dubordieu,
1956; Nedjari-Belhacene & Nedjari, 1984)
En se dirigeant vers l’Est les deux chaines (Atlas tellien et Atlas Saharien tendent à se
confondre dans l’Est de l’Algérie et en Tunisie. Puis, elles se rejoignent pour donner
une chaîne unique qui constitue l’Atlas Oriental qui s’étend de l’Est de l’Algérie à la
Tunisie, (Figure 5).
Plusieurs auteurs ont utilisé les mesures de fracturation et l’analyse microtectonique
détaillée dans l’Atlas tunisien pour reconstituer les « paléo-états de contraintes ». La
synthèse de l’ensemble de ces travaux, concernant la tectonique compressive, montre
que :

18
1/ Depuis le début du Cénozoïque jusqu’à l’Actuel, la déformation compressive est
dominante, et de direction NW-SE caractérisée par une distension locale au niveau des
structures distensives NW-SE(Dlala & Rebai, 1994; Chihi, 1995; Dlala, 1995).
2/ Pour d’autres auteurs ; elle est compressive depuis le Cénozoïque à l’actuel ; et
interrompu par des phases distensives (Bouaziz et al., 2002…)
Il existe deux phases de déformation majeure qui correspondent, à la Phase miocène
supérieur (Tortonien-Serravalien) et la Phase Villafranchienne, ou pour d’autres
auteurs, Post Villafranchienne. En effet, au Miocène supérieur, la phase de
compression du Tortonien inférieur est caractérisée par une direction de
raccourcissement orientée N120 à N140 et elle a engendré les plis « atlasiques » des
Atlas tunisiens central et septentrional et les plis en échelon de la Tunisie méridionale.
Lors de ce plissement, les accidents anciens proches de E-W et de N-S ont été réactivés
en décrochements N070 à N120 dextre et N140 à N-S sénestres (figure 6).

Figure 5: Les principaux domaines structuraux de l’Afrique du Nord (d’après De Lamotte et al., 2000)

19
Rôle des Maghrédides dans l’accommodation de la convergence Afrique-Eurasie:

La subduction de la plaque africaine sous la plaque Eurasiatique a donné naissance à

l’orogenèse des Maghrébides (de Direction E-W). Après l’arrêt de la subduction (Néogène)
d’autres processus permettaient l’accommodation de la convergence au niveau des
maghrébides.
Une 1ère phase compressive est enregistrée au niveau de la marge Sud-Tyrrhénienne mais
mal connue. Certains auteurs (Butler et al., 1995; Ghisetti et al., 2009) supposent que
l’arrêt du volcanisme au niveau des îles volcaniques (exemple de l’Ustica) au cours du
Pléistocène moyen à supérieur, pourrait être relié au début de la phase compressive dans
cette région. D’autres auteurs mentionnent qu’au niveau des régions internes du sud du
tyrrhénien, la phase compressive a commencé bien avant le pléistocène moyen (Ghisetti,
1979; Pepe et al., 2000; Pepe et al., 2005).
D’après Andrea Billi (Billi et al., 2011) on peut conclure que suite à l’arrêt de la
subduction (il y a environ 15 Ma), on assiste à une inversion tectonique (il y a 8 à 10Ma)
par la mise en place d’une phase compressive au niveau du domaine d’arrière-arc (algéro-
liguro-provençal et tyrrhénien) (Billi et al., 2011).
Le long de la marge tyrrhénienne la compression semble alors accommodée par une
ceinture de plis et chevauchements actifs (les Maghrébides). En effet, après la fin de la
subduction, au niveau de l’Atlas tunisien (partie orientale des Maghrébides), on note la
succession de deux phases majeures de plissement :
(1) Serravalien Tortonien : Elle a induit le rejeu des failles est–ouest à N120 selon une
composante normale, et la formation de plis NE–SW
(2) Phase villafranchienne (Quaternaire ancien) : responsable de l’activation des failles
est–ouest à N120, en décrochement dextre inverse, et des failles NE–SW à
composante inverse, et de la genèse des plis NE–SW.
Les directions de raccourcissement se répartissent varient de N-S à NW-SE selon la
direction de collision des deux plaques Africaine et Eurasiatique.

20
21
Figure 6 : Coupes géologiques équilibrées montrant les déformations tectoniques dans le domaine atlasique tunisien. en haut à gauche : Coupes
schématiques simplifiées illustrant l’évolution des principaux niveaux de décollement et le rôle de l’héritage et la pré structuration tectonique
dans la déformation ; à droite : Carte de localisation de la coupe ;en bas : Coupe géologique Tabarka-Chott Fejaj montrant le rôle du Trias-Lias
dans le décollement et la dominance de la déformation pelliculaire dans l’avant-pays atlasique tunisien. (Ahmadi, 2006)
A-3 : tectonique récente et sismicité

Au niveau des régions Maghrébides, les déformations récentes et les zones tectoniques,
sont irrégulièrement distribuées. Elles sont concentrées essentiellement sur les chaînes
orogéniques Maghrébides (Atlas et Tell). L‘intensité de cette déformation récente
décroît considérablement en direction de l’Est. Cette déformation récente ne suit pas
toujours le style structural du socle (Skobelev et al., 1988; Trifonov, 2004). Plusieurs
types de structures actives caractérisent la déformation récente dans les régions
maghrébides. En effet, cette déformation récente se présente le plus souvent sous forme
de failles anciennes réactivées au Quaternaire. Ces failles sont généralement des failles
inverses accompagnées par des plissements de la couverture sédimentaire mésozoïque
et cénozoïque, et plongent généralement vers le Nord ou le NW. Ces failles sont
développées le long des flancs de jeunes anticlinaux, au cours du Pléistocène supérieur
et de l’Holocène. L’aspect structural des failles actives et des jeunes anticlinaux
indique une compression subhorizontale de direction NW-SE. Le taux de soulèvement
de ces jeunes anticlinaux dépasse la vitesse de déplacement des failles associées qui
varie entre 0,2 et 0,9 mm/an (Domzig et al., 2006).

A-4 : cinématique de la Méditerranée

A-4 –1 : mesures des champs de déformation (GPS)

La convergence de la plaque Africaine vers la plaque Eurasie et la collision des

fragments de l’Eurasie, déterminent les conditions cinématiques aux limites de
l’Afrique du Nord. Le champ de déformation régional montre une direction de
compression horizontale maximale N-S à NW-SE (Noquet & Calais, 2003; Nocquet &
Calais, 2004; D'Agostino et al., 2008).
Plusieurs travaux ont étudié les mouvements relatifs aux limites des plaques africaines
et eurasiatiques. Ces travaux ont été focalisés sur les quantifications géodésiques des
mouvements relatifs de l’Afrique par rapport à l’Eurasie (en choisissant l’Eurasie

22
comme repère). Hollenstein et al., (2003) ont étudié par GPS le mouvement des
plaques entre la Tunisie et l’Italie par rapport à l’Europe fixe (figure 7).
La majorité des stations indique des déplacements vers le Nord et vers l’Est. Des
stations sont proches de la Tunisie et l’encadrent au Nord et à l’Est. Les résultats de ces
stations sont présentés dans le tableau 1 :

Localité Long. Lat. VN (mm/an) VE (mm/an)

LAMP 12.606E 35.5N 2.5 ± 0.8 -3.5 ± 0.8
TRAP 12.582E 38.004N 3.5 ± 0.5 -3.5 ± 0.6
CAGL 8.973E 39.136N -0.1 ± 0.1 -0.9 ± 0.1
PANT 11.945E 36.832N 2.4 ± 0.5 -2.0 ± 0.5
Tableau 1 : Vitesses de déplacement mesurés par GPS dans des
localités de la Méditerranée (Noquet & Calais, 2003)

Figure 7: Vitesses GPS mesurées par rapport à l'Eurasie. les flèches vertes représentent les mesures GPS
effectuées entre 1994 et 2001. les flèches bleus les mesures basées sur le GPS continue ; a) les stations IGS
(International Association of Geodesy) utilisées comme référence ((Hollenstein et al., 2003); modifié)

La compilation de ces résultats permet de déduire la vitesse de déplacement des

localités mesurées par rapport à l’Europe fixe (tableau 1). La station installée au niveau

23
de l’île de Lampedusa, qui se trouve à l’Est de la Tunisie à 125 km de la côte orientale
tunisienne, montre que cette localité se déplace vers le NW (≈305°N) avec une vitesse
de ≈4,3 mm/an. La localité de Trapani (Nord-Ouest de la Sicile) qui se trouve à 160 km
au NE de la Tunisie se déplace vers le NW (≈315°N) avec une vitesse de ≈4,95 mm/an.
Une station GPS installé à Cagliari au sud de la Sardaigne, à 200 km du Nord de la
Tunisie indique que l’Ile de Sardaigne se déplace avec une vitesse de ≈0,9 mm/an vers
le WSW (≈264°N). Au niveau de l’île de Pantelleria qui se trouve à 70 km au Nord-Est
de la Tunisie le GPS continu montre que cette île se déplace vers le NW (≈320°) avec
une vitesse de ≈3,12 mm/an. D’autres études récentes ont été réalisées au niveau de la
Méditerranée Occidentale et centrale en utilisant le GPS continu pour estimer les
mouvements relatifs des plaques (Hollenstein et al., 2003). Ces résultats indiquent une
convergence de direction N45W +/-20° à la longitude de la Sicile, qui s’approche
progressivement vers une convergence E-W au niveau du Détroit du Gibraltar (figure
8). Les évaluations de taux de convergence sont entre 3 à 7 mm/an à la longitude de la
Sicile. Vers l’Ouest à la longitude du Détroit du Gibraltar le taux de convergence
décroit, et il est estimé de 2 à 5 mm/an.

Figure 8: Modèles de déplacement de l’Afrique par rapport à l’Eurasie (Nocquet & Calais, 2004).

Les solutions Géodésiques de la convergence Afrique-Eurasie (Argus et al., 1989;

DeMets et al., 1990; Kreemer & Holt, 2001; Sella et al., 2002) montrent une direction
de convergence variant entre le N-S et le NW-SE et une vitesse qui varie de 2,8 à 6
mm/an au niveau du Nord de la Tunisie (figure 9). Les solutions géologiques (DeMets

24
et al., 1990) qui donnent une estimation de vitesses de convergence entre le nord de
l’Afrique et l’Eurasie proche de 7,3 mm/an, montrent une nette différence avec les
solutions géodésiques, qui donnent une moyenne de vitesse de convergence aux
alentours de 4mm/an selon une direction proche de NW-SE (figure 10). En effet, nous
remarquons que les solutions géodésiques indiquent clairement que le taux de
convergence entre l’Afrique et l’Europe Stable est plus lent que la prédiction des
modèles géologiques avec un facteur allant de 1,2 à 2. Cette différence est
probablement due aux lacunes de données et paramètres géologiques lors de
l’établissement des ces modèles géologiques.

Figure 9: Modèle géodésiques de déplacement de l’Afrique par rapport à l’Eurasie

A/ modèle de déplacement proposé par Noquet et al., 2003
B/ modèle proposé par Hollenstein et al., 2003 (mesure entre 1994 et 2001)

La distribution de la sismicité et les résultats géodésiques, qui indiquent un taux faible

de déformation sur les côtes NW de l’Afrique, ainsi que le mouvement NW de la
plaque Africaine, implique que les déformations en Afrique du nord et en Ibérie

25
méridionale résultent de l’accommodation de la convergence Afrique-Europe en
Méditerranée Occidentale.
Les données géodésiques et sismologiques suggèrent donc que la déformation résultant
de la convergence oblique Afrique-Eurasie est concentrée le long de l’Afrique du nord.
À l’Est de la longitude 2°W, la distribution spatiale de la sismicité montre que la
plupart des grands tremblements de terre sont concentrés autour d’une bande
relativement étroite le long de la côte Nord Africaine (figure 10). La convergence
oblique NW-SE entre l’Afrique et l’Eurasie en Méditerranée Occidentale est
compatible avec les mécanismes focaux des séismes qui montrent un mouvement en
faille inverse de direction NE-SW et de pendage généralement vers le NW combiné
avec un mouvement décrochant.

Figure 10: Carte de la cinématique des plaques africaine et eurasiatiques (Nocquet & Calais, 2004)

L’étude de la transition rhéologique entre la croute continentale et la croute océanique

au nord de l’Algérie a montré que la limite du domaine continental est située

26
approximativement au pied de la pente du bassin profond au niveau de l’isobathe
2600m (Domzig et al., 2006).

A-4-2 : implication sur le déplacement actuel en Tunisie

Les valeurs de convergence estimées par les méthodes géodésiques(Nocquet & Calais,
2004), montrent une valeur moyenne de raccourcissement entre l’Afrique et l’Eurasie
de 5 mm/an. Cette valeur a été déduite à la longitude de la Sicile au NE de la Tunisie.

En effet des mesures continues de GPS entre la Tunisie et la Sardaigne (Hollenstein et

al., 2003; D’Agostino & Selvaggi, 2004; Nocquet & Calais, 2004) ont permis de
montrer que la Tunisie se déplace selon la direction Nord-Ouest par rapport à l’Europe
stable. Les valeurs moyennes de déplacement sont vers le Nord de 3,6 mm/an, et vers
l’Ouest de 2,6 mm/an correspondant à un vecteur déplacement général au niveau de
cette longitude de 4,4 mm/an dans une direction N315°. Cette direction coïncide avec
la direction de raccourcissement définie statistiquement par les données
microtectoniques de la déformation cassante quaternaire (Dlala, 1995). La déformation
engendrée par cette convergence est accommodée par des structures actives dans la mer
et en Tunisie.

27
B : Tectonique récente et sismicité en Tunisie

B-1 : Tectonique récente en Tunisie

En Tunisie, les déformations récentes d’âge quaternaire ont été signalées depuis le
début du XXème siècle; de nombreux auteurs (Vaufrey R, 1932; Laftine & Dupont,
1948; Castany, 1951) ont signalé l’importance de la tectonique compressive quaternaire
qu’ils considéraient comme responsable d’une orogenèse « post-villafranchienne ». Les
travaux ultérieurs confirmèrent l’existence de déformations dans le Quaternaire ancien
et récent (e.g. (Kamoun, 1981; Chihi, 1984; Dlala, 1984; Ben Ayed, 1986; Philip et al.,
1986).
De nombreux indices de déformation récente ont été mis en évidence à l’échelle de tout
le pays. On peut citer par exemple, la faille de Gafsa (Sud Ouest de la Tunisie) où les
dépôts d’âge Acheuléo-Moustérien sont plissés et tronqués tout le long de cet accident
actif (Vaufrey R, 1932; Zargouni & Ruhland, 1981; Dlala, 1992; Dlala & Hfaiedh,
1993; Dlala, 1995) ; la Tunisie centrale et dans le bassin du Sahel, des exemples de
déformation d’âge Pléistocène supérieur ont été décrits (Ben Ayed et al., 1978;
Kammoun, 1981; Dlala & Ben Ayed, 1988). De même, en Tunisie septentrionale, les
dépôts marins d’âge Tyrrhénien à Néotyrrhénien, qui affleurent dans les régions
côtières, sont souvent intensément déformés (Dlala, 1992; Dlala, 1995).
En Tunisie de nombreux indices de déformations très récentes ont été mis en évidence
à l’échelle de tout le pays. Ces indices de déformations récente se manifestent en
surface par :
- Des ruptures dans les formations les plus récentes, souvent au niveau des failles
préexistantes (réactivation de failles préexistantes).
- Des plis, ou torsion de couches qui affectent les niveaux géologiques superficiels, en
rapport avec les mouvements de failles dans la couverture sédimentaire ou dans le
socle.
- Soulèvement des terrasses, indiquant des mouvements verticaux ;
- Affaissement des bassins quaternaires ou escarpement des bordures des fossés
préexistants, indiquant aussi des mouvements verticaux ;

28
- Anomalies morphologiques, surtout celles qui modifient l’itinéraire d’un réseau
hydrographique.
- Sources hydrothermales qui traduisent les remontées d’eaux profondes

B-2 : Sismicité : Caractères et problématique

L’étude de la sismicité du bassin de la Méditerranée occidentale révèle la complexité

de cette convergence datant depuis le Miocène moyen. Cette sismicité est à caractère
diffus et montre que la zone de contact Afrique –Europe est large et recouvre
l’ensemble du bassin méditerranéen et des chaînes plissées (par [Link] chaînes
atlasiques) qui la bordent (Hatzfeld, 1978).
L’examen de la sismicité dans le domaine allant du détroit de Gibraltar à l’Est jusqu’en
Tunisie, montre que la séismicité est relativement faible et se répartit sur une large
bande dans le domaine continental alors que les épicentres en mer sont peu fréquents
(figure 11).
Lorsque les frontières de plaques sont matérialisées par des failles bien définies (c’est
le cas par exemple de la faille de San Andreas en Californie), les séismes s’y
produisent d’une façon récurrente, avec des constantes de temps de l’ordre de quelques
dizaines à quelques centaines d’années, correspondant à la durée nécessaire pour
accumuler les contraintes jusqu’au point de rupture. En revanche lorsque les limites de
plaques sont moins nettes (par exemple la convergence de deux continents), comme
c’est le cas de l’Afrique et l’Eurasie, la déformation se produit sur des bandes de
plusieurs centaines de kilomètres. La sismicité est alors diffuse et complexe à étudier.
Elle se répartit sur quelques failles majeures, et plusieurs centaines de failles modestes,
et est caractérisée par des constantes de temps de l’ordre du millénaire ou davantage.
C’est le cas du bassin Méditerranéen, dans lequel se situe la Tunisie.
La zone la plus sismique s’étend sur l’Afrique, la Grèce, Italie, le Proche-Orient,
ainsi que sur le nord du Maghreb. Les événements sismiques destructifs les plus
remarquables de ces dernières décennies sont ceux d’El Asnam, 1954; El ; Acerrno,
1980; El-Asnam, 1980; Chelaf (Alger), 1996 ; Zemmouri, 2003; Boumerdes, 2003 ; et
du Maroc (Agadir), 1960.

29
30
Figure 11: Séismicité du bassin de la Méditerranée (IPSN-BERSSIN)
L’observation du catalogue de la sismicité en Tunisie montre une sismicité caractérisée
par des magnitudes modérées. Toutefois, certains événements observés dans le
catalogue de la sismicité de la Tunisie sont considérés destructeurs, étant donné les
intensités (MSK) ou les magnitudes qui leurs ont été accordées.
Les enregistrements sismiques depuis l’installation du réseau sismologique national en
1976 jusqu’à 2005, et les données néotectoniques permettent d’identifier les principales
zones sismiques suivantes en Tunisie :
- le Nord Ouest du pays : la magnitude maximale observée a été de 5,7 et la
magnitude maximale probable serait de 6,5
- le nord et la région de Tunis : la magnitude maximale observée a été d 5.1 et la
magnitude maximale probable serait de 7.
- Le Sahel : la magnitude maximale observée a été en continent 4,6 et au large
5,1. La magnitude maximale probable serait de 6,5.
- La région de Gafsa : la magnitude maximale a été de 4,8. la magnitude probable
serait de 6,5.
Il est a noter que les magnitudes maximales probables pour toutes les régions du pays
ont été établies sur la base des investigations sismotectonique réalisées dans le cadre du
projet de cartographie quaternaire mené par l’Institut National de la Météorologie
(Arsovsky et al., 1991).

B-2-1 : Sismicité historique (figure 12)

Une recherche bibliographique dans le but de reconstituer les évènements sismiques

de forte intensité a été réalisée par l’Institut national de la Météorologie de la Tunisie et
la bibliothèque nationale de Tunisie. Elle couvre la période de 412 à 1975. Au Nord-
Est de la Tunisie, l’évènement historique le plus destructeur est le séisme d’Utique de
l’année 412, qui a causé la destruction totale de la ville antique d’Utique (Vogt, 1993).
Dans la région de Tunis, les années 856, 1425 et 1785, ont enregistré des évènements
sismiques très importants qui sont à l’origine des pertes humaines et des destructions
des maisons (Vogt, 1993). On cite que le séisme de 856 a fait 50000 victimes (données
douteuses de l’INM). En général, et particulièrement en Tunisie, la sismicité historique
est peu connue, les informations sont fragmentaires et parfois sujettes à caution. Cette

31
sismicité dépend de la densité de peuplement des régions, c’est ce qui explique que
Tunis (Carthage) et Kairouan présentent une sismicité historique relativement élevée
par rapport au reste des régions de la Tunisie.
Dans l’histoire récente (depuis le 19 ème siècle) la Tunisie, contrairement à ses voisins
Maghrébins (Maroc et Algérie), n’a pas connu de catastrophe sismique majeure. La
magnitude maximale observée a été de 5,7 lors du séisme de Bou Salem (haute vallée
de la Medjerda) en 1957. Les séismes ayant eu lieu en Tunisie, au cours de cette
période récente, ont en général laissé des dégâts matériels plus ou moins sévères allant
de fissures dans les bâtiments jusqu’à des dommages importants rendant certaines
constructions inutilisables. Cependant, en plus de ces dégâts matériels nous avons
enregistré des pertes de humaines à trois reprises (8 morts en 1887 à Jemmel ; 3 morts
en 1889 dans la région de Gafsa ; et 9 morts et 102 blessés en 1957 à Bou Salem). Le
tableau suivant (tableau 2) montre les principaux séismes destructeurs, dans un sens,
enregistré dans le territoire tunisien dont la magnitude est disponible dans le catalogue
de sismicité de l’INM.

Date Epicentre Magnitude Localité

20/02/1957 36,400 09,000 5,7 Bousalem
12/08/1958 35,300 10,700 4,9 Sidi-Thabet
01/12/1970 36,908 09,947 5,6 Sidi-Thabet
19/01/1970 36,915 08,413 5,0 Ghardimaou
07/02/1977 36,479 08,413 5,1 Mz Nour Djemmel
09/04/1979 37,112 10,136 4,2 Ras Djebel
15/03/1984 34,398 08,208 4,2 Redeyef
24/06/1987 34,478 09,250 5,0 Sned
25/06/1987 36,617 09,615 4,8 Medjez el Bab
07/11/1987 34,342 08,399 4,3 Metlaoui
17/09/1994 36,461 09,219 4,6 Thibar Teboursouk
22/09/1995 35,439 08,415 4,9 Sahraoui
29/03/1996 34,881 09,373 4,3 Bir El Hfay

Tableau 2: plus importants séismes enregistrés en Tunisie depuis 1957

32
Figure 12 : Carte de la Sismicité historique de la Tunisie (le Sud du pays ne montre pas de sismicité
historique ; données de l’INM)

33
 B-2-2 : sismicité instrumentale

Suite au séisme de décembre 1970 à Tunis, un réseau sismologique de 3 stations

télémétrées a été installé autour de Tunis en 1976. Ce réseau est composé de trois
stations sismologiques courte période à une composante verticale installées
respectivement à Menzel Bou Zelfa, Zaghouan et Kechabta.
En 1989, dans le cadre d’un programme national de réduction du risque sismique, un
réseau de 14 stations sismologiques courte période (fréquence de 1 Hz) et à trois
composantes, couvrant l’ensemble du pays a été installé (figure 13). Ce réseau
comprend cinq sous réseaux régionaux dans les zones sismiques et un centre principal à
Tunis. Chaque sous réseau comporte trois stations. Les signaux des trois stations
télémétrées sont reçus au sous-centre où l’enregistrement se fait de façon analogique et
simultanément ces signaux sont transmis par une ligne téléphonique spécialisée au
centre principal de Tunis.
Le catalogue de la sismicité instrumentale pour la période de 1976 à 2005 comprend
aux environs de 2000 enregistrements. L’examen de la répartition des épicentres des
séismes montre qu’ils se regroupent en plusieurs essaims. En effet, nous distinguons
différentes zones (figure 14) :
 Des régions à forte activité sismique qui se situent dans le nord de la Tunisie où
elles englobent les secteurs de Tunisie, Bizerte, Medjez el Bab, Jendouba,
Ghardimaou ; dans l’est du Pays les régions de Sousse, Monastir et Mahdia. Dans le
sud la région de Gafsa.
 Des régions à faible activité sismique. Il s’agit des régions situées dans l’Ouest du
pays ; au centre ouest et enfin dans la région du Sud-Est.
 Des régions à sismicité presque absente : ce sont les régions peu peuplées de la
plate-forme saharienne au Sud du pays..
Dans la région du Nord Ouest de la Tunisie, la sismicité est le reflet de l’activité du
linéament Ghardimaou-Cap Serrat. De même, la répartition des épicentres dans la
région de Bou Salem et dans la région de Gaâfour confirme l’activité des linéaments de
Teboursouk et El Alia (figure 15).

34
Figure 13: Configuration du réseau sismologique de la Tunisie
(donnée de l’Institut National de la Météorologie)

35
 Figure 14: Carte de la sismicité instrumentale de la Tunisie
(données de l’Institut National de la Météorologie de 1975 à 2005)

36
 Figure 15: Carte de Répartition des épicentres du Nord de la Tunisie (données
de l’Institut National de la Météorologie de 1975 à 2005)

B-2-3 : Complexité de cette sismicité intraplaque

Durant la période historique, des violents événements sismiques ont été signalés en
Tunisie. Ces événements sont d’une ampleur qui suggère des évènements sismiques
intenses dans le territoire tunisien.
D’autre part, rares sont les séismes tunisiens qui montrent de ruptures de surface qui
peuvent nous renseigner directement sur la cinématique des accidents actifs. Parmi ces
séismes on cite le séisme de Tunis de 1970 (données de l’INM), séisme de Ghardimaou
de 1985 (Gueddiche, 1992), séisme de Métlaoui de 1989 (Dlala & Hfaiedh, 1993).
Donc le recours au calcul des solutions focales de ces séismes est nécessaire pour avoir
une idée sur l’état de contraintes en profondeur et la cinématique des accidents actifs.
Cette répartition montre que la sismicité en Tunisie ne se concentre pas toujours autour
des alignements tectoniques actifs reconnus. Ceci serait expliqué par les erreurs de
localisation des épicentres et aussi par la nature des structures tectoniques qui peuvent
engendrer des tremblements de terre.

37
La caractérisation de cette sismicité consiste donc à analyser la déformation
néotectonique dans les structures relativement jeunes, c’est-à-dire celles formées après
la phase atlasique du miocène supérieur. Ceci nécessite donc à la fois des études
sismologiques mais également géologiques et géophysiques. Dans le cadre de ce
travail, nous allons discuter l’activité des importantes sources seismogénique du Nord-
Est de la Tunisie, et nous insisterons en particulier sur la faille d’Utique et son
implication sismotectonique.

38
 CHAPITRE II

Géodynamique et tectonique récente

du Nord-Est de la Tunisie

39
A : Description lithostratigraphique

A-1 : Le Néogène

Les bassins du Nord-est de la Tunisie sont caractérisés par une accumulation

importante de dépôts néogènes, résultant d’une subsidence continue depuis l’Oligo-
Miocène jusqu’au Pliocène (Dlala, 1994; Dlala, 1995; Kacem et al., 2001). Ces dépôts
qui résultent du démantèlement des reliefs avoisinants, ont alimenté les bassins
subsidents, morcelés par des horsts et des grabens. Les plus importants affleurements
de cette zone sont résumés dans cette partie (figure 16).

A-1-1 : Miocène supérieur

Le Miocène se présente comme un ensemble, souvent très épais (il atteint 2500m
d’épaisseur) formé par des termes argileux sablo-gréseux auquel s’ajoute du gypse. Ces
termes reposant en transgression sur les terrains précédents. Les intercalations
gypseuses au sein du Miocène supérieur constituent de véritables discontinuités
lithologiques qui permettent de subdiviser cette méga-séquence en 3 séquences d’ordre
inférieur :
 une séquence continentale
 une séquence intermédiaire de plate-forme interne
 une séquence supérieure continentale
Le Miocène supérieur est constitué de couches bigarrés (formation Hakima ou
Mellaha), par les marnes grises de l’Oued El Melah, par la formation Kechabta, et est
recouverte par des facies de fin de cycle ou de régression : séries de l’Oued bel Khedim
(figure 16).

- Formation Hakima

Cette formation est caractérisée par une série bigarrée. Elle est formée par des sables,
grès, marnes, avec des intercalations de gypse vers le sommet de la série, et localement

40
des horizons de poudingues ou de brèches. Les éléments de ces niveaux sont grossiers
et appartiennent essentiellement à l’Eocène. La couleur est très variable (verdâtre,
blanc, rouge, bleu ou noir) d’où la nomination de série bigarrée. La proportion de
marnes s’accroît dans la partie supérieure de cet ensemble. On trouve aussi des passées
de grès jaune à « ripple marks ». Les microfaunes y sont toutes remaniées ; on trouve
d’abord celles de l’Eocène, ensuite celles du Crétacé au-dessus d’un horizon de marne
lie-de-vin qui apparaît au sommet de la formation.
La formation Hakima a une épaisseur de l’ordre de 500 mètres dans ce secteur ; la
teinte varicolore de cet ensemble détritique permet de l’individualiser sans difficulté du
reste de la série sus-jacente.
Au cœur du Djebel Messeftine on trouve des affleurements de la formation Hakima. En
effet, les couches bigarrées affleurent très largement dans les profonds ravins qui, de la
ferme Lebray au Sud, s’étendent dans tout le secteur du Djebel Melah jusqu’au Nord,
sur la pente Sud du Djebel Berna ; série très lagunaire, formée presque exclusivement
de marnes, parfois légèrement sableuses, rouges, grises ou noires avec beaucoup de
gypse et un peu de sel ; presque pas de grès francs.
Cette formation est attribuée au Serravalien (Burollet, 1951), et correspond à un
système côtier voire lagunaire, très peu profond affecté par des fluctuations importantes
de la ligne de rivage et déposé dans un contexte fortement subsident (El Euch, 2007).

Formation Oued El Melah

Elle est composée surtout par des marnes et argilites grises qui se développent sur plus
de 500 mètres dans plusieurs localités du secteur d’étude. Les niveaux argileux de cette
formation contiennent du gypse en cristaux isolés ou en bancs. La série débute
localement par des marnes à lits de poudingues qui remanient essentiellement des
calcaires éocènes ; au sommet, elle admet des bancs de grès jaunes qui annoncent la
formation suivante à laquelle elle peut passer latéralement.
La microfaune dans cette formation est entièrement remaniée, riche vers la base et très
pauvre au sommet, elle provient du Campanien, du Paléocène, de l’Eocène supérieur et
de l’Oligocène inférieur. En plus des formes remaniées, la faune y est représentée par

41
des Ostracodes et Gastéropodes a test lisse et des Charophytes, indices de milieu
lacustre (Ammar, 1986).
Cette formation est attribuée au Tortonien inférieur ; et correspond à l’installation d’un
système transgressif de lagune ou de Sebkhas associé à un taux de subsidence
important (Dlala, 1994; Kacem, 2004; El Euch, 2007).

Formation Kechabta

C’est une alternance argilo-gréseuse et argilo-sableuse pouvant atteindre 1000 mètres

d’épaisseur dans ce secteur. Elle est constituée par une alternance de marnes et de
bancs gréseux fins, certains bancs présentent à leur partie supérieure des ripples marks
On note l’absence de faciès lagunaires et continentaux francs, et la rareté des fossiles
(néritique et terrestre).
A la base de la série, un ensemble argileux, teinté à certains niveaux de jaune et de
rouge, peut être individualisé avec localement des grès grossiers. Ce sous-ensemble se
développe sur près de 400 mètres ; à ce sommet, se situent des horizons de grès à
ciment de galène qui ont fait l’objet de recherche à proximité d’El Alia (Rouvier,
1977). Le reste de la formation où alternent les niveaux argileux et les grès à patine
jaune présente pour l’essentiel une teinte grise. Le gypse apparaît au sommet de la
succession et annonce la formation suivante.
Au site de Kechabta, cette formation est beaucoup plus gréseuse qu’à Ras Djebel ; la
portion de grès sans être aussi forte que celle des marnes, est presque égale et, en outre,
les marnes sont souvent très sableuses. L’épaisseur est souvent très importante : plus de
2000m au flanc Sud du Kechabta. Le flanc Nord présente des épaisseurs beaucoup
moindres. Cette formation est attribuée au Tortonien supérieur (Burrolet, 1956;
Colleuil, 1976).

Formation Oued Bel Khedim

La formation de l’Oued Bel Khedim groupe un certains nombre de faciès très

différents, tous lagunaires ou lacustres, représentant généralement la fin du Miocène

42
supérieur. Elle correspond à une série Messinienne présentant à la base un intervalle
lagunaire suivi par un intervalle lacustre (El Euch, 2007). Les marnes et les calcaires de
cette formation marquent d’abord un retour à la sédimentation lagunaire fortement
influencée par la mer, suivie par des dépôts lacustres (Burrolet, 1951). La base de la
série est argileuse et elle contient des huîtres et des polypiers et admet de nombreux
bancs de gypse. Elle se développe dans ce secteur sur 150 mètres environ. Les faciès
lagunaires sont surmontés par des marnes noires et des niveaux de calcaires lacustres
blancs dont l’épaisseur est de l’ordre d’une centaine de mètres.

Formation Chaâbet – Tabbala

Cette formation correspond à la partie sommitale des séries miocènes de ce secteur.

Elle est attribuée par Burrolet (Burrolet, 1951) au miocène supérieur (Pontien). Il s’agit
d’un épisode de sédimentation continentale représentée essentiellement par des argilites
sableuses, brunâtres ou rougeâtres et des grès grossiers. L’épaisseur moyenne de cette
formation est de l’ordre de 50 mètres, et il dépasse les 200 mètres sur l’axe de
l’anticlinal de Messeftine. Cette série comprend des brèches torrentielles, des tufs
calcaires et un appareil fluviatile rouge plus évolué. Le mélange des galets arrondies
polygéniques et des galets calcaires anguleux dans la brèche de base suggère la
contemporanéité d’apports locaux et plus lointains, en l’occurrence un appareil
fluviatile (faciès rouge) dans lequel s’injectent des apports torrentiels latéraux
(brèches).
Cette formation est attribuée au Messinien supérieur (Ammar, 1986),

43
Figure 16: Log stratigraphique du puits Utique1 (Rapport ETAP)

44
 A-1-2 : Pliocène

Le Pliocène est composée par deux formations:

-la formation Raf-Raf, formée essentiellement de marnes, surmonté par :
-la formation Porto Farina, caractérisé d’une dominance gréseuse.
Les analyses micropaléontologiques (Bizon et al. 1980 ; Besème et Lajmi, 1981)
indiquent un âge Tabianien (Pliocène inférieur) pour les marnes de Raf-Raf. La
formation Porto-Farina est d’âge Pliocène supérieur.
Ces deux formations (Raf-Raf et Porto Farina) supportent les marnes et sables
continentaux à faunes de Mammifères du villafranchien inférieur, découverte par E.
Dumon (Laftine & Dupont, 1948).

Formation Raf-Raf

Au niveau du Dj Merzaguine (Besbès et Lajmi, 1981), cet ensemble est formé par des
argiles gris-bleu à gris verdâtre, légèrement sableuse, débutant par un faciès de
transgression, souvent marqué par des conglomérats et passant à leur partie supérieure
aux grès jaune-ocre de Porto-Farina (figure 17).
Au niveau de djebel Messeftine les marnes de Raf-Raf affleurent largement. Elles sont
souvent épaisses (de 400m à l’Est à 650 m au Sud-ouest) et contiennent plusieurs bancs
de grès (3 sur la pente ouest du djebel Krerba et à l’Est du douar Messeftine).
Au niveau du plissement d’Utique –Besbessia, les affleurements des marnes de Raf-
Raf sont rares. Elles affleurent par endroits à la faveur de l’érosion des ravins sous la
carapace des grès de porto-Farina

Formation Porto-Farina

Au niveau du Dj Merzaguine, cet ensemble est à dominance sablo-gréseuse de teint

jaunâtre à ocre et elle est riche en fossiles (Pectinidés et Ostréidés). Au sud de djebel
Besbessia, les grès reposent sur les marnes de Raf Raf. Ils sont assez épais, parfois
assez grossiers. Dans certains secteurs de la région étudiés, des études en lames minces

45
(Burrolet, 1951), ont montré que ces grès sont quartzeux à ciment calcaire, à nombreux
tests de lamellibranches ou plus rarement de gastéropodes et de très rares foraminifères
(figure 17).

Figure 17: Corrélation entre les forages pétroliers disponibles dans la région d’étude

46
 A-2 : Le Quaternaire

Plusieurs auteurs ont décrit le Quaternaire de la région du nord-est de la Tunisie :

(Solignac, 1927; Arambourg & Arnauld, 1950; Burollet, 1951; Castany, 1953;
Crampon, 1971; Paskoff & Sanlaville, 1979; Bourgou, 1982; Paskoff & Sanlaville,
1983; Oueslati, 1993)..

A-2-1 : Pléistocène inférieur (Villafranchien)

Les dépôts attribués au Quaternaire ancien, occupent, principalement, les piémonts des
grands reliefs de bordures montagneuses. Ils sont constitués, généralement, de
formations détritiques avec des alternances de bancs de conglomérats plus ou moins
grossiers qui se distinguent par leur couleur rouge.
Le Villafranchien inférieur débute par des calcaires lacustres, des encroûtements
calcaires ou gypseux de nappe massifs et des brèches ou poudingues rosés surmontés
d’un sol rouge ou châtain-rouge (Nord tunisien) correspondant à la première phase
lagunaire à Cardium du Chott Djérid et à de petits lacs échelonnés dans les bassins du
Nord du pays. La faune à éléphant a donc pu vivre dans un environnement de type
savane. Cette faune s’est poursuivi, en déclinant, jusqu’à la fin de l’épisode qui a connu
la mise en place de sols calcimorphes localement intercalés de conglomérats provenant
du démantèlement des grands massifs mis en place au Pliocène. L’apparition du
Villafranchien supérieur traduit un changement par le développement d’un
environnement de steppe identifié, dans les plaines, par l’évolution de sols châtains.
Fournet (1971) a étudié le villafranchien situé sur la rive nord du Lac Ichkeul. Sa base
est représentée par des argiles brunes attribuées au Pliocène. Un conglomérat de galets
de calcaires cristallin vert ou noir plus au moins métamorphisés, épais de quelques
trente centimètres, les sépare d’une série argilo-sableuse d’à peu près 50 mètres
d’épaisseur qui leur succède. C’est dans les sables et graviers consolidés jaune-rouille
qui en forme la base que se situe un gisement de mammifères. Cette série est intercalée
de bancs d’argile dont les deux premiers ont fourni une flore riche et abondante, ainsi
que de rares mollusques d’eau douce. La suite de la série est essentiellement argileuse,
avec intercalation de quelques bancs de grès. (Fournet, 1971)

47
(Burollet, 1951) a montré que le Pléistocène inférieur est formé de limons et de sables
continentaux, des argiles et des sables rubéfiés, ensuite on trouve des argilites, et des
silts rouges avec de petites « poupées calcaires » blanches qui sont des bioturbations.

A-2-2 : Pléistocène moyen

Il est constitué de deux parties principales : la base est formée de limons et de sables
continentaux alors que le sommet est constitué par des argiles et des limons rouges
avec souvent de petites « poupées » blanches, parfois amalgamées en « croutes ».
L’épaisseur de ces couches peut atteindre les 100 mètres en particulier autour de
Kechabta à l’Est de la route de Bizerte (Burrolet, 1951).
Nous avons levé une coupe géologique au niveau d’un vallon situé à quelques
centaines de mètres à l’Ouest de la ville d’Utique (figure 18). La base de la coupe est
dominée par les grès du Pliocène qui forment le cœur de l’anticlinal, et affleurent par
endroit au bas du relief. Les grès sont surmontés par une couche silteuse de couleur
verdâtre et à concrétion calcaire. Cette couche est attribuée au Pléistocène moyen.
Ensuite on rencontre une couche de limons rouges à abondance de concrétions
calcaires et qui montre des passages conglomératiques grossiers empruntés aux grès
pliocènes. A la fin de la coupe on rencontre une couche argilo-sableuse riche en Hélix,
de couleur rouge et d’âge attribué au quaternaire récent.
Au Djebel Menzel-Besbessia à Bir El Fejja, une coupe levée par Oueslati (1994), au
niveau d’une carrière (figure 19), montre une succession qui peut être subdivisée en
trois parties principales.
La partie inférieure est dominée par un matériel argilo-sableux de couleur variant entre
le rose, le beige et le vert pâle. Son sommet est caractérisé par d’importantes
concentrations calcaires reposant directement sur des argiles brunes caillouteuses.
Deux autres horizons d’enrichissement en calcaire existent encore plus bas et sont
développés sur des lits d’argiles rouges à verdâtre compactes.
La partie moyenne est aussi à dominance argilo-sableuse de couleur rose pâle, elle se
distingue toutefois par le fait qu’elle emballe plusieurs cailloux et blocs empruntés aux
grès pliocènes. A son sommet apparaît un encroutement pulvérulent ou à structure plus
ou moins feuilletée.

48
De même la partie supérieure est argilo-sableuse, mais de couleur rouge vif. A côté de
nombreuses concrétions calcaires, elle renferme des passées caillouteuses. Une croûte
calcaire scelle sa surface.

A-2-3 : Pléistocène supérieur

Les dépôts marins et éoliens du Quaternaire supérieur sont rencontrés sur les bordures
du Lac de Ichkeul et de Bizerte. Paskoff et Sanlaville (1983) ont choisi Ras el Korane
pour effectuer une coupe synthétique du Pléistocène supérieur, en effet, dans cette
coupe on distingue de bas en haut :
- Une plage de galets à macrofaune et Mélobésiées qui repose sur un substratum gréso-
marneux ;
- Des grès marins coquilliers grossiers riches en faunes ;
- Des marnes sableuses gris-vert à faune lagunaire ;
- Des grès coquilliers dunaires, passant latéralement à des grès de faciès marin attesté
par la macrofaune (Cette séquence est rapportée, selon (Oueslati, 1993), à la formation
Régiche qui correspond au Tyrrhénien et plus précisément à l’Eutyrrhénien) ;
- Une plage riche en macrofaune et en galets mal émoussés souvent lithophagé ( selon
Oueslati (Oueslati, 1993), Cette plage, est attribuée à la formation Chebba datée du
Néotyrrhénien).
- Des dépôts colluvial sablo-limoneux, rouge à jaunâtre.
- Des dépôts marins (Cérithium, Nassa) à lagunaires (Bittium, Théricium) passant
latéralement à des grès dunaires coquilliers à Hélix.
Ces dépôts sont rattachés à la formation Cap Blanc datée Wurmien.
- Des sables gris à gris-noir contenant des coquilles de Gastéropodes terrestre.
- Des sables quartzeux rouge-jaunâtre meubles.
- Des sables quartzeux actuels de couleur claire

49
Figure 18: Coupe à Djebel Menzel Ghoul au niveau de flanc Sud de l'anticlinal d'Utique,

Figure 19: Coupe dans le Quaternaire au niveau de Menzel (Oueslati, 1993)

50
 A-2-4 : Le Tyrrhénien

Les chercheurs s’intéressent depuis plus d’un siècle au Tyrrhénien de la Tunisie

(Quaternaire marin) car la côte tunisienne est riche en dépôts marins attribuables au
dernier Interglaciaire. Castany a établi en 1962 une première synthèse du Tyrrhénien
tunisien, quelque peu retouchée ensuite par Jauzein (1967) et Oueslati (1982). Ces
auteurs distinguent essentiellement deux dépôts : d’une part, une « dune ancienne »,
oolithique, à faune banale ; d’autre part, une « plage ancienne », contenant des
Strombes et attaquant en falaise la dune oolithique dont elle englobe des témoins sous
forme de blocs et de galets.
Dans la plupart des régions il n’y a pas seulement une ligne de rivage tyrrhénienne,
mais plusieurs (généralement deux ou trois), ce qui traduit des oscillations plus ou
moins importantes du niveau marin autour d’un maximum transgressif.
Le tyrrhénien en Tunisie à été sujet de plusieurs subdivisions qui se rattachent
généralement aux facteurs climatiques, on distingue les formations suivantes, de bas en
haut (figure 20) :
- Formation Régiche
- Formation Chebba
- Formation Douira
D’autres subdivisions proposent deux types de formations : l’Eutyrrhénien et le
Néotyrrhénien. En Tunisie l’Eutyrrhénien est équivalent de la formation Régiche, et le
Néotyrrhénien, transgressif sur l’Eutyrrhénien est équivalent à la formation Chebba
(Martinez & Paskoff, 1984; Oueslati, 1987).
Les dépôts tyrrhéniens, ont été reconnus à des altitudes différentes : soit +7 à +10
mètres aux environs de Bizerte, +7 à +13 mètres au Cap Bon (golf de Hammamet)
(Paskoff & Sanlaville, 1983) et +13 à +32 mètres en marge du Sahel tunisien (Jedoui et
al., 2002; Bouaziz et al., 2003; Elmejdoub & Jedoui, 2009). L’altitude de l’équivalent
de ces dépôts tyrrhéniennes, retrouvées dans de nombreuses régions du monde réputées
par leur faible activité tectoniques est de +6 mètres (Chen et al., 1991; Lambeck &
Nakada, 1992). Donc, l’existence, sur les côtes tunisiennes, de formations
tyrrhéniennes contemporaines, à des altitudes variables, résulte probablement de
mouvements néotectoniques verticaux.

51
 Figure 20: les différentes subdivisions du Tyrrhénien en Tunisie (d’après Mahmoudi, 1998)

A-2-5 :L’Holocène

Pendant l’Holocène, deux épandages de limons rouges ruisselés et de colluvions se sont

mis en place (Oueslati, 1985). Le premier date du milieu de l’Holocène et témoigne
probablement d’une récurrence pluviale, comme le suggèrent des analyses polliniques
(Brun, 1985). Le second est d’âge historique. Il contient souvent des tessons de poterie
d’anciennes civilisations, et pourrait avoir une origine anthropique, sans doute les
défrichements étendus opérés dans l’antiquité (Brun, 1983). Mais, l’hypothèse d’une
oscillation climatique, caractérisée par une incursion sèche, entre la fin de l’antiquité et
le XVe siècle, ne doit pas être écartée (Brun & Rouvillois-Brigol, 1985).

52
B : Tectonique de la Région du Nord-Est de la Tunisie

B-1 : Tectonique anté-néogène

une carte isochrone du toit du Trias a été élaborée par Kacem (Kacem, 2004) en
compilant les données des profils sismiques et des puits pétroliers disponibles dans la
région (figure 21). Sur cette carte qui met en évidence la tectonique post-triasique, nous
remarquons que la disposition des structures hautes, des structures basses et les
directions des failles et des linéaments majeurs, plaident en faveur d’une tectonique
décrochante distensive. En effet, la carte montre un réseau de failles normales orientées
principalement NNE-SSW, NE-SW et E-W. ces failles délimitent plusieurs structures
en horsts et grabens.

Depuis le Jurassique, la tectonique dans le Nord-Est de la Tunisie semble être guidée

essentiellement par deux directions majeures de failles:
- Direction N-S : on note comme exemples la faille de Menzel Jémil-Aïn Ghelel,
la faille de Messeftine (Dlala, 1994; Dlala, 1995).

- Direction E-W : on cite par exemple la faille NE-SW de Tunisie-El alia

(Burrolet, 1951) ainsi que la faille E-W à NE-SW de Djebel Nahli-Djebel
Ammar (Dlala, 2002; Dlala & Kacem, 2008). Ces failles majeures qui ont
morcelé le bassin jurassique en horst et grabens, sont soit des failles normales
ou décrochements à composantes normales.

Au Paléocène, les structures préexistantes : les failles, les structures hautes et basses
rattachées à la période Jurassique – Crétacé, persistent encore. Elles ont été réactivées
par un régime distensif à décrochant distensif, ceci est montré par les variations
importantes d’épaisseur qui se répartissent selon une disposition géométrique bien
guidée par les accidents majeurs les plus importants de la zone (figure 22).

53
54
Figure 21: Carte d'isochrone du toit du Trias (Equi. = 100 ms) (Kacem, 2004)
55
Figure 22: Carte d'isochrone du toit du Crétacé-Jurassique (Equi. = 100 ms) (Kacem, 2004)
 B-2 : Tectonique Néogène

Dans la région du Nord-Est de la Tunisie, le Néogène dans son ensemble est déformé
en plis de grande amplitude dont la direction des axes varie du N-S à E-W.
La tectonique néogène est connue en Tunisie par une succession de déformations en
compression, dont la plus importante, qualifiée de phase atlasique, a eu lieu au
Miocène supérieur. Elle est responsable des plissements importants qui ont affecté
l’ensemble du Nord de la Tunisie et remis en mouvement le diapirisme triasique
(Perthuisot, 1978).
Au Miocène inférieur à moyen une subduction de la plaque Africaine sous la plaque
eurasiatique s’effectue au niveau de l’Afrique du Nord – Sicile (Reuther & Eisbacher,
1985; Bousquet & Philip, 1986; Oldow et al., 1990; Roure et al., 1990). Ce régime
compressif se traduit par la mise en place des nappes de charriages, au même moment
que l’ouverture du bassin marginal Algérien. A la même période s’ouvraient les fossés
d’effondrement d’avant-pays de la chaîne alpine en Tunisie et en Mer Pélagienne, en
régime décrochant distensif.
Au cours de la phase Miocène supérieur, les failles E-W à N120 sont réactivées en
décrochement à composante normale en associant souvent des plis en échelon ou
d’entrainement surtout en Tunisie centrale ou méridionale (Zargouni & Ruhland, 1981;
Boukadi, 1994; Chihi, 1995; Dlala, 1995; Zouari, 1995)

B-3 : Tectonique Quaternaire

La fermeture océanique migre progressivement vers l’Est et gagne la Sicile : pendant le

Quaternaire l’ensemble Tunisie du Nord – Sicile entre en collision (Philip, 1987; Chihi,
1995). Cette phase compressive du quaternaire ancien (souvent appelée post-
villafranchienne) étant à l’origine des plis et failles actives du Nord-est de la Tunisie
(figure 23) (pli d’Utique, anticlinal de Messeftine,…) (e.g. (Burrolet, 1956; Rouvier,
1977; Ben Ayed, 1986; Dlala, 1995).
Depuis le Tyrrhénien à la période actuelle, soit pendant les derniers 120 000 ans, la
persistance de la compression est mise en évidence par :

56
- L’existence de déformations affectant les dépôts tyrrhéniennes et les alluvions
récentes (Ben Ayed, 1980; Dlala, 1992; Dlala, 1995; Kacem, 2004) ;
- Le rejeu du décrochement de Skanes-Khniss, près de Monastir, décalant une
mosaïque romaine du II’ siècle après J.-C. (Kamoun et al., 1980) est un témoin de la
persistance de cette déformation compressive;
- Des études de mécanismes au foyer de plusieurs séismes récents survenus en Tunisie
septentrionale.
Cependant, l’observation de déformations extensives en quelques localité en Tunisie
témoigne de variations locales dans le champ de contraintes pendant cette période
(Dlala & Ben Ayed, 1988; Dlala & Rebai, 1994; Dlala, 1995).

Figure 23: Carte tectonique de subsurface de la région du Nord-Est de pays

(Kacem, 2004)

57
 Introduction

Dans le domaine de la Tunisie nord-orientale l’abondance des dépôts quaternaires bien

datés (Arambourg et al., 1949; Jauzein, 1967; Paskoff & Sanlaville, 1983) et l’activité
tectonique localement bien exprimée (Ben Ayed, 1986; Oueslati, 1987; Dlala, 1992),
suggèrent que ce domaine mérite une exploration profonde en terme d’activité
tectonique récente.
Sur les côtes nord orientales, des formations marines littorales, datées du Pléistocène
supérieur (Paskoff & Sanlaville, 1980; Paskoff & Sanlaville, 1982; Oueslati, 1985; Ben
Ayed, 1986) ont subi des déformations typiques de compression (plis, failles inverses
et décrochements) qui indiquent une direction de raccourcissement proche du NW-SE
qui est active durant les derniers 125 ka. Son activité actuelle est démontrée par les
calculs des solutions focales des séismes récents en Tunisie (Girardin et al., 1977;
Hfaiedh et al., 1985).
En Tunisie Septentrionale, la néotectonique quaternaire se manifeste non seulement par
une déformation de surfaces géomorphologiques, mais également par un plissement des
couches du pléistocène ancien et moyen qui enveloppent le front de quelques
anticlinaux.
La sismicité relativement importante de la Tunisie est l’expression d’une tectonique
active. Plusieurs séismes historiques et récents ont provoqué des effets de surface qui
peuvent renseigner sur l’état de contraintes dans les régions touchées (Hfaiedh &
Chadi, 1985; Gueddiche, 1992; Dlala & Hfaiedh, 1993).

A : Les principaux séismes dans le Nord de la Tunisie (après JC)

A-1 : Séisme historique d’Utique

La ville antique d’Utique constituait la plus importante cité phénicienne après

Carthage. Elle a été fondée entre 750 et 1000 ans avant J-C. Cette cité qui était l’un des
ports les plus importants de la Méditerranée, occupait la pointe d’un promontoire situé
à la terminaison orientale de Menzel Ghoul. Aujourd’hui, ce promontoire est situé à
une dizaine de kilomètres de la mer, par suite de l’envasement du Golfe d’Utique

58
(Reyners, 1951; Jauzein, 1967); Jauzein, 1967). La ville d’Utique a connu un séisme
historique, particulièrement destructeur vers 412 après JC (Sieberg, 1932; Vogt, 1993).
En se reportant aux données de la littérature, des fissures dans le sol sont apparues suite
à ce séisme qui a été suivi par plusieurs répliques pendant une semaine.

A-2 : Séisme historique de Tunis

Un violent tremblement de terre a secoué la ville de Tunis en l’an 857 (Sieberg, 1932;
Hfaiedh & Chadi, 1985; Vogt, 1993). L’intensité de ce séisme a été estimée de VIII
(Vogt, 1993). Les données de la littérature indiquent un millier de victimes sous les
déblais des bâtiments détruits, et des dégâts considérables dans le bâti de la ville de
Tunis.

B : Sismicité instrumentale à effet de surface (depuis 1900 après JC)

B-1 : Séisme du 1er décembre 1970

Le 1er décembre 1970 à 02 h 04 mn, un séisme de magnitude 5,6 secouait la ville de

Tunis. Les effets secondaires de ce séisme ont été localisés dans la zone épicentrale,
près de Chaouat à environ 35 km à l’Ouest de Tunis. Dans ce secteur des fentes
ouvertes de directions N100 à N110 ont affecté la plaine alluviale (Girardin et al.,
1977). Le long de ces fissures, des jets d’eau et de boue, associés à des volcans de boue
(alignés sur ces fentes, ont été observés par la population.

Sur la route principale de Chaouat (GP7, km 27) : les poutres en béton d’un petit pont
qui porte la route sont affectées par des failles décrochantes inverses de direction N100
à 110. Le rejet horizontal de ces failles peut atteindre 5 cm tandis que le rejet vertical
est de quelques millimètres (Girardin et al., 1977; Hfaiedh et al., 1985; Ben Ayed,
1986).

Les failles décrochantes inverses permettent de déterminer graphiquement la direction

de l’axe de raccourcissement maximum (Z). Celui-ci est dirigé NNE-SSW. Cette

59
direction de compression est proche des solutions focales proposées pour ce séisme
(Girardin et al., 1977; Hatzfeld, 1978; Hfaiedh, 1983; Hfaiedh et al., 1985).

B-2 : Séisme de Sidi Thabet (12/12/1970)

Ce séisme de magnitude 5,7 a été le plus fort dans la région du Grand Tunis depuis au
moins l’installation du réseau sismologique de la Tunisie. Ce séisme a provoqué des
phénomènes de liquéfaction qui ont été observé dans la localité de Sidi Thabet au Nord
de Tunis (figure 24).

Figure 24: Cratère de sable observé à Sidi Thabet après le Séisme du 12/12/1970
(donnée l’Institut National de la Météorologie)

B-3 : Séisme de Ghar el Melh (9/04/1979)

Ce séisme de magnitude 5,1 est situé à l’extrême nord-est dans la région de Ras Djebel
– Raf Raf (lat : 37,18°N ; long : 10,11°E). Ce séisme n’a pas fait de victimes, mais il a
causé des dégâts matériels. Une fissure d’une trentaine de mètres de direction voisine
d’E-W, a été repérée par les pêcheurs du vieux port Turc de Ghar el Melh (figure 25)

60
sur le quai du port. Cette fissure a affecté les pavements du quai, faits de blocs de grés
noyés dans un ciment de chaux.
Des observations d’effets de surface ont été décrites par Ben Ayed (1986) et qui n’ont
pas été confirmé par d’autres sources d’informations. En effet cet auteur a signalé qu’à
50m au Sud du port Turc, des fissures provoquées par ce séisme ont permis à l’eau de
la nappe phréatique de jaillir au milieu de la lagune, lors du séisme. Il a signalé aussi
que des glissements de terrain ont été observés sur les deux flancs de Djebel Ennadour
et Eddmina, de même que des chutes sporadiques de blocs de pierres sur les pentes
fortes du flanc Sud du Djebel Ennadour(Ben Ayed, 1986).

Figure 25: Localisation géographique de la ville de Ghar El Melh et

vue panoramique du vieux port de la ville

B-4 : Séisme de Chemtou (8/06/1985)

La magnitude de ce séisme qui a touché la région de Chemtou située au NW de la

Tunisie (gouvernorat de Jendouba), a été évaluée à 4,2 (Ms) de l’échelle de Richter.
Selon Gueddiche (Gueddiche, 1992) ce séisme a provoqué des fissures verticales
importantes, des fentes décamétriques de direction N010E à N025W et de largeur

61
centimétrique à décamétrique, et des petites fentes orientées N080E à N120E (figure
26). Les fissures N010W qui sont disposées en échelon montrent un jeu senestre. Ce
séisme correspond probablement à la réactivation de la faille profonde de Hairech qui
met en contact anormal, les calcaires métamorphiques avec les formations gréso-
pelitiques de Hairech (Gueddiche, 1992).

Figure 26: Fracturation de la surface suite au séisme

de Chemtou (08-06-1985) (Gueddiche, 1992)

C : Etude des Mécanismes aux foyers

C-1 : Le concept de mécanisme au foyer

La détermination des mécanismes aux foyers se base sur le sens du premier mouvement
sur le sismogramme ; il s’agit soit d’une compression (le signal est dirigé du foyer vers
la station), soit d’une dilatation (le signal est dirigé de la station vers le foyer).
En général, pour la détermination des mécanismes aux foyers, on suppose que les
forces mises en jeu se ramènent à un double couple, d’où la division de l’espace en

62
quatre quadrants séparés par deux plans perpendiculaires (plans nodaux) dont l’un est
le plan de la faille, l’autre fictif est le plan auxiliaire. Les axes de pression et de tension
maximale sont alors situés dans les plans bissecteurs des plans nodaux.
Les premiers mouvements enregistrés autour de l’épicentre se partagent en 4 quadrants,
2 de dilatations, 2 de compressions, qui définissent 2 plans orthogonaux appelés plans
nodaux. Ce partage (existant également pour les ondes S) est tout à fait cohérent avec
le modèle du double couples de deux forces appliquées au foyer, soit une combinaison
de tension et de compression.
Si l’on considère que le mouvement se fait dans un milieu homogène et isotrope, on
peut alors calculer l’amplitude théorique des ondes P dans une direction par rapport au
plan de faille (Brune, 1976) :
A = (M*sin (2i)*cos (2j) / (4(pi)*d*V3*R)
A : l’amplitude des ondes de compression ; M : le moment des couples ;( R, i, j) : les
coordonnées sphériques de la station ; d : la densité du milieu, V : la vitesse des ondes
P.
Ces hypothèses étant admises, on détermine donc 2 plans nodaux perpendiculaires,
suivant lesquels l’énergie des ondes P est nulle : l’un de ces deux plans est le plan de la
faille. Les axes Pression et Tension, responsables du double couple, sont situés sur les
plans bissecteurs de ces 2 plans nodaux. L’axe P est situé dans le quadrant des
dilatations, l’axe T dans celui des compressions.
La figure 27 résume les conventions adoptées dans les représentations stéréographiques
des mécanismes au foyer. La position des sens de premier mouvement observé en
chaque station sur le diagramme dépend de l’azimut de la station par rapport à
l’épicentre et de l’angle d’incidence de l’onde sous la station. La représentation de ces
informations sur le sphère focale permet de déterminer la position des plans nodaux.
L’identification du plan de la faille et celui de l’auxiliaire reste une tache aussi
complexe à accomplir à partir des données des mécanismes aux foyers. La
détermination du plan de la faille se fait par une étude géologique ou à partir de la
localisation des épicentres des répliques pour un séisme de magnitudes importante.

63
Figure 27: Représentation stéréographique des mécanismes au foyer (pour les séismes proches, la projection
de Schmidt sur la demi-sphère inférieure est le plus couramment utilisée) et signification tectonique de ces
mécanismes. Axe de Pression: l'étoile noire; axe de Tension: l'étoile blanche; quadrant de compressions en
noir

C-2 : Mécanismes aux foyers : champs de contraintes au Nord de la Tunisie

Nous avons eu recours aux mécanismes au foyer de quelques séismes pour caractériser
le champ de contrainte actuel du Nord du pays. Nous avons fait une synthèse des
solutions focales reportées à la littérature additionnée de quelques mécanismes que
nous avons calculés en appliquant un algorithme automatique (Focal Mecanism ;
IISEE) (figure 28). La construction géométrique de la projection pour un mécanisme au
foyer et les calculs d’erreurs sur l’azimut, l’angle d’incidence et la profondeur ont été
traités en utilisant les données de l’INM (Institut National de Météorologie). Pour
effectuer ces mécanismes au foyer nous avons cherché les sismogrammes dans
l’archive du service de sismologie de l’Institut National de Météorologie, et nous avons

64
refait le dépouillement pour s’assurer du premier arriver de l’onde Pg et de la durée du
mouvement. Les plans nodaux sont représentés par une direction et un pendage, les
axes P et T par leur direction et leur angle de plongement.

Les axes de pression de la majorité des mécanismes aux foyers sont compatibles avec
la direction de la compression quaternaire qui est orientée NNW-SSE. Ces résultats
sont compatibles avec celles prouvés par les méthodes néotectoniques définies par
plusieurs auteurs (Philip & Meghraoui, 1983; Ben Ayed, 1986; Dlala & Rebai, 1994)
(Philip & Meghraoui, 1983; Phlip & Meghraoui, 1983; Ben Ayed, 1986; Dlala, 1994) ,
et les méthodes sismologiques obtenues par Hfaiedh (1983) et (Hatzfeld & Bensari,
1977; Hatzfeld, 1978).

Figure 28: Carte des principaux mécanismes au foyer du Nord de la Tunisie

et la région Nord orientale de l’Algérie

65
  Séisme du 1er décembre 1970 :
Le calcul du mécanisme au foyer du séisme de Djdeida-Chaouat (lat. : 36,9°N ; long :
9,95°E) dont les effets de surface sont décrit ci-dessus a été réalisé par plusieurs auteurs
(Girardin et al., 1977; Hfaiedh, 1983). Ce mécanisme a montré qu’il s’agit d’un
mécanisme d’une faille décrochante inverse. L’axe de pression maximale (P) est
orienté N178 avec un plongement de 11° vers le Sud (figure 29). Le plan de faille est
vraisemblablement orienté N40, compatible avec une direction NE-SW jouant en
sénestre.

Figure 29: Solution focale de quelques séismes au Nord-Est du pays

d'après Hfaiedh (1983)

 Séisme du 19 janvier 1977 :

Il est situé dans la région de Ghardimaou connue pour son activité séismique
importante, sa magnitude est 4,6 (échelle de Richter) (36,58°N ; 8,49°E) (figure 28). Il
présente un mécanisme de faille inverse (figure 29). L’axe P est orienté N138 avec un

66
plongement de 16°S. Le plan orienté N27 serait vraisemblablement le plan de la faille,
car il est compatible avec le linéament Ghardimaou-Cap Serrat qui a un jeu senestre
(Hfaiedh, 1983).

 Séisme du 9 avril 1979 (Hfaiedh, 1983):

La solution focale de ce séisme dont la description des effets de surface ont été décrits
ci-dessus, montre qu’il présente un mécanisme en faille normale. L’axe T est N50
(figure 29). Le plan N125, compatible avec l’existence de l’accident de Metline jouant
en dextre, peut être choisi comme étant le plan de la faille. En effet, les séismes dans
cette région se localisent suivant une direction à peu près E-W.

 Séisme du 19 février 1993 (Mejri, 2004):

Ce séisme s’est produit le 19 / 2 / 1993 (temps origine : 20h 44’ 25,10’’) au niveau de
Djebel Ichkeul (lat : 37,13°N ; long : 9,57°E). Sa magnitude a été évaluée de 3,9 et il
est de profondeur 11,4 km. Le mécanisme au foyer individuel de ce séisme montre
qu’il s’agit d’un mécanisme en faille normale avec une légère composante décrochante.
L’axe T est de direction N107 (figure 30).

Une enquête macrosismique à été menée par le service de la sismologie de l’Institut

National de Météorologie. Les isoséistes établis à l’issue de cette enquête montrent une
zone épicentrale qui s’étend sur une trentaine de kilomètres suivant la direction NE-
SW, et sur une vingtaine de kilomètres suivant la direction NW-SE. La zone
épicentrale de ce séisme est définie par le maximum d’intensité VI (rapport interne de
l’INM) observé dans les localités suivantes : Oued el Melah, Sidi, Mohamed ben
Othman et Bir Felta.

67
 Figure 30: Mécanismes au foyer des séismes du 19/02/1993 et 26/01/1995

 Séisme du 26 janvier 1995

L’épicentre du séisme du 26 / 1 / 1995 (temps origine : 14h 51’ 36,22’’) est situé près
de Garâat Ichkeul (lat. : 37,09°N ; long : 9,72°E). Il est de magnitude 4,3 et de
profondeur 10km.
Le mécanisme au foyer individuel de ce séisme montre qu’il s’agit d’un mécanisme en
faille normale avec une légère composante décrochante. L’axe T est presque NW-SE
(figure 30). L’emplacement de l’épicentre calculé reporté sur la carte des indices de
déformations récentes, montre que la faille ayant rompu a une direction probable de
N083 (Hfaiedh, 1998).

68
  Séisme du 12 aout 2001 :
La magnitude de ce séisme qui a touché la région de Kairouan (lat. : 35,95°N ; long. :
9.92°E ; heure d’origine : 14h19’ 30,3’’) a été estimée à 4,3. La solution focale de ce
séisme a montré qu’il s’agit d’un mécanisme en faille décrochante et que la direction
de l’axe P est de 291N (figure 31) (données du site [Link]

Figure 31: Mécanisme au foyer du séisme de 12/08/2001

 Séisme du 24 juin 2002 :

Ce séisme dont l’épicentre est au Nord de Kairouan (lat. : 36,01 ; long. : 10,01), de
magnitude 4,6 et de profondeur 15 km, a montré un mécanisme en faille inverse avec
un axe P orienté 270N (figure 32) (données du site [Link]

69
 Figure 32: Mécanisme au foyer du séisme du 24/06/2002

 Composite Medjerda-Utique :

Ce mécanisme a été élaboré par Hfaiedh (1983), il regroupe tous les séismes de la
région de Ddjeida, Utique et le delta de la Medjerda. C’est un mécanisme en faille
inverse avec une forte composante de cisaillement (figure 29). Cette solution, comparée
à celle du séisme du 1/12/1970 dans la même région, ne parait pas très différente. L’axe
de pression est N15.

70
D : Aperçu sur les principales structures actives au NE de la Tunisie

D-1 : Etat de connaissances des structures actives (exemples de littérature)

La région du Grand Tunis est limitée à l’Est par le golfe de Tunis, au Nord par
la basse plaine de la Medjerda, à l’Ouest par la longitude de Jedeida-Borj el Amri et au
Sud par la plaine d’Oued Chafrou. La ville de Tunis est encadrée à l’Est par le lac de
Tunis, au SW par le Sebkhat Essijoumi, au Nord par les collines de l’Ariana (Djebel
Nahli, Djebel Ammar), au sud par Djebel Sidi Bel Hassen et Djebel Kharrouba (figure
33).

Figure 33: Modèle numérique de terrain avec les principaux reliefs de Tunis et du NE du pays

La région du grand Tunis coïncide avec la fermeture nord-orientale de l’Atlas tunisien.

Elle est considérée comme le point de convergence des principales zones structurales

71
de la Tunisie (l’axe N-S, la faille de Zaghouan, l’atlas septentrional et le Sahel). Cette
convergence des structures principales du pays se traduit par un découpage du socle de
la région de Tunis par des failles majeures orientées NE–SW, N–S et NW–SE (Boutib
et al., 1997). Dans ce réseau, trois accidents actifs, subparallèles et orientés N60
(accidents de Djebel Ammar, Mornaguia-Tunis et Mohamedia), sont impliqués dans la
tectonique mio-plio-quaternaire de cette région. Les principales failles de cette région
se trouvent au pied des reliefs (figure 34), comme les failles du Djebel Ammar ou celle
du Djebel Nahli (Kacem, 2004; Dlala & Kacem, 2008).

72
73
Figure 34: Carte sismotectonique: 1: Trias; 2:Mio-pliocène; 3: bassin numedien; 4: anticlinal; 5: synclinal; 6: décrochement; 7: faille
inverse; 8: failles confirmée par les données géophysiques; 9: indices de déformations récentes; Z: zonage sismotectonique(Z1, Z2 et
Z3); 10: sismicité instrumentale; 11: sismicité historique (Dlala & Kacem, 2008)
 D-1-1 : L’accident NW du Djebel Nahli

Le Djebel Nahli est une structure anticlinale de direction N20. Sa structure est
compliquée par des petits sous-anticlinaux hectométriques à décamétriques.
Stratigraphiquement, cette structure englobe des faciès allant de l’Aptien au
Quaternaire (carte géologique de l’Ariana, 1971). Des accidents verticaux et
subverticaux découpent la série du Djebel Nahli ; un des accidents majeurs s’étend sur
une dizaine de kilomètres de direction N140 et de pendage 72NE (figure 35). Cet
accident décale la structure du Djebel Nahli en décrochement sénestre. Il affecte, à la
bifurcation de la route de Bizerte vers l’Ariana et Bardo. Les derniers encroûtements
calcaires et les silts rouges d’âges pléistocène supérieur.

Figure 35: Carte géologique de la structure de Djebel Nahli

74
Indices de néotectonique au niveau de Djebel Nahli

A l’extrémité orientale du J. Nahli, à l’ouest d’une ancienne carrière d’exploitation de

matériaux de construction, au piémont de ce massif, les dépôts quaternaires sont
représentés par des sables rouges et des silts rouges encroûtés au sommet. Les
encroûtements et les silts rouges sont affectés par le décrochement NW- SE de Djebel
Nahli. Ceci montre l’importance de cet accident majeur dont le jeu persiste jusqu’au
Quaternaire récent (Dlala, 2002).
Dans une tranchée dégagée par la route de Bizerte, en face du cimetière, apparaissent
des décrochements de direction NW-SE qui affectent les derniers encroûtements et les
dépôts pulvérulents (Dlala, 1994).

D-1-2 : La Faille d’El Alia

El alia est une ville montagneuse, située au Nord à 40 km de Tunis. Le secteur d’El
Alia-Metline montre un accident rectiligne de direction N40. Il est jalonné de Trias en
faisant un contact anormal entre le Cénomanien et l’Eocène. Il est long de 25km et
profond de 4 à 5km (Kacem, 2004). Cette faille NE-SW est réactivée au Quaternaire en
décrochement senestre, avec une composante inverse.
Pour Rouvier (1977) l’évolution de cet accident est polyphasée et il s’est intéressé
uniquement à la phase post-Miocène. En effet, avec le dépôt du Miocène laguno-
continental discordant, prend naissance une suite de déformations le plus souvent peu
accentuées, et qui ont pu évoluer en plis coffrés et en plis failles avec injection de Trias
(Rouvier, 1977). Cette suite de déformation tardives, qui affecte le Miocène discordant
s’est poursuivie jusqu’à l’époque actuelle puisque localement du tyrrhénien est affecté
de failles inverses (Paskoff & Sanlaville, 1983).

D-1-3 : Déformation tectonique Quaternaire récent de Ras Engla

Les reliefs côtiers de Ras Engla sont séparés de la mer par une plate-forme d’abrasion
tyrrhénienne (Crampon, 1971; Miosec, 1977). Une étude antérieure (Oueslati, 1993) a

75
montré que cette plate-forme d’abrasion correspond à deux niveaux de plate-forme
étagées par une falaise morte.
Oueslati (1993) a mentionné que la plate-forme supérieure serait anté-tyrrhénienne, et
que celle inférieure est tyrrhénienne mais polygénique puisqu’elle porte deux
générations de dépôts marins fort probablement de l’Eutyrrhénien et du Néotyrrhénien.
En plus de ces dépôts marins, ces plates-formes portent différentes formations, liées à
des agents continentaux (ruissellement et vent).
à l’Est du phare de Ras Engla, les plages tyrrhéniennes ainsi que les dépôts
continentaux qui les séparent dessinent avec les séries du substratum triasique, un
anticlinal décamétrique déversé vers le Nord (Ben Ayed, 1993). A la charnière de ce
dernier, on remarque des fentes de tension et des microfailles normales d’extrados de
rejet décimétrique et de direction N80 à N100 sur le flanc Sud de ce pli une
discordance intra-tyrrhénienne existe, ce qui montre que le plissement est contemporain
de ces dépôts (figure 36).
Dans une falaise à environ 200 m à l’Est du Phare, une faille inverse de direction N70-
80 fait chevaucher les conglomérats de la plage néotyrrhénienne sur les dépôts
continentaux de la formation Aïn Oktor wurmienne (Ben Ayed, 1986; Oueslati, 1993;
Dlala, 1995).
A l’Est du phare et au droit d’un îlot on voit, dans la falaise, une faille inverse N100-
120 pentée de 70° vers le Sud. Elle fait reposer les dépôts tyrrhéniens sur une couche
d’âge historique. Ceux-ci sont formés de sables limoneux gris à lits de graviers. En plus
de débris de poterie et d’ossements, elle a livré une pièce de monnaie ce qui permet de
bien dater la faille car elle date du bas empire romain (II à IVème siècle A J.-C)
(Oueslati, 1993).
L’état des affleurements de la faille inverse de direction N70-80 citée plus haut n’a pas
permis d’observer des stries, mais les déplacements visibles sur des galets sectionnés
par le plan de rupture confirment le mouvement inverse de la faille et indiquent un
raccourcissement vraisemblablement proche de N-S (Oueslati, 1993).

76
77
Figure 36: Coupe de la falaise à l’Est du Phare de Ras Engla (Ben Ayed et Oueslati, 1988)
 D-1-4 : La faille inverse d’El Fejja

Le pli d’Utique de direction proche d’E-W s’étend sur environ 15km, depuis la ville
d’Utique à l’Est jusqu’à la région de Besbessia à l’Ouest. Ce pli est moulé sur un
accident proche d’E-W. les sections sismique montrent que cette faille s’étend sur
environ 25 km et elle est profonde de 5 km (Kacem, 2004). Cette faille préexistante a
été réactivée en décrochement dextre inverse au Quaternaire (Dlala & Kacem, 2008)
(cf. article).
A l’extrémité Ouest de l’anticlinal d’Utique et au niveau du village d’El Fejja (figure
37), nous avons observé sur le flanc Sud de l’anticlinal une faille inverse visible depuis
la route qui mène vers Menzel Bourguiba (37.0366°N ; 9.98009°E); cette faille de
direction N173, 30W, 25S (figure 38), affecte la série Pliocène et la base du
Quaternaire et le décale d’environ 30cm. Cette faille proche de N-S constitue le réseau
conjugué de la faille d’Utique. l’analyse microtectonique affectant les dépôts
Quaternaires dans la localité d’El Fejja, a montré une direction de raccourcissement
NW-SE ((Dlala, 1995)

Figure 37: Localisation géographique de la faille d'el Fejja

78
Figure 38: Faille Pliocène inverse au niveau d’El Fejja (Ouest de l’anticlinal d’Utique).
Elle touche la base du Quaternaire

79
 D-2 : Nouvelles observations de la Tectonique récente au Nord-Est de la
Tunisie

Outre les failles actives signalées dans la littérature et décrites ci-dessus, nous avons
cherché d’autres indices d’activité tectonique quaternaire dans la région d’étude. Nous
avons mis en évidence plusieurs indices de cette activité dans la région de Tunis au
niveau des structures plissées du bassin de Kechabta-Messeftine.

D-2-1 : La faille active de Djebel Ammar

Le Djebel Ammar est une structure anticlinale de direction SW-NE (N60) dont le cœur
est coupé par une faille jalonnée par du Trias (Pini et al., 1971; Devolve et al., 1981).
Or, l’analyse structurale des flancs de ce djebel (Ben Ayed et al., 1983) de part et
d’autre de la dépression triasique, montre que le flanc Sud-est correspond plutôt à un
anticlinal de direction N25 qui s’amortit au SSW par une terminaison périclinale très
régulière (Ben Ayed et al., 1983).
Le flanc Nord-Ouest de cette structure est bordée par une faille normale à composante
dextre de direction proche d’Est-Ouest qui sépare la structure anticlinale du Djebel
Ammar de la plaine d’El Mabtouha. Ce dernier accident correspond en réalité à une
faille profonde majeure reconnue en subsurface. Elle s’étend de la région du Jedeida à
l’Ouest jusqu’au massif de Djebel Nahli à l’Est (Dlala et al., 2008) (Kacem, 2004). Cet
accident a joué un rôle tectono-sédimentaire au moins depuis le Jurassique d’après les
travaux de Kacem (2004).
Cet accident reconnue en subsurface à affecté les dépôts du Pléistocène supérieur en
failles normales ; au bord de la route de Sidi Thabet (Dlala & Kacem, 2008)

La Tectonique active de Djebel Ammar :

Tranchée 1 : (figure 39)

Une nouvelle tranchée à été effectuée dans le cadre de réalisation d’un projet de
construction a permis de dégager les indices d’une importante activité tectonique. En

80
effet, cette tranchée de direction N-S, profonde d’à peu près 5m et qui s’allonge sur
environ 40mà révélé l’existence de failles normales affectant les formations
quaternaires. Les plans de ces failles de direction proche d’E-W montrent souvent des
remplissages conglomératiques.
Les dépôts affectés sont les dépôts des lits d’un chenal. Il semble que le jeu quaternaire
de cet accident a engendré la création d’un grand parcours d’eau.

Tranchée 2 : (figures 40 et 41)

Une deuxième tranchée a été effectuée par l’occasion des travaux d’excavations
pour la construction de bâtiments. Cette tranchée est longue de 50 mètres, et est de
direction proche d’E-W (N70°). Elle montre une profondeur moyenne de 5 mètres. Les
formations géologiques sur cette tranchée sont d’âge Pléistocène supérieur. En effet
une croute mince et feuilletée surmonte un faciès sableux fluviatile. Cette croute est
surmontée par des argiles brunâtres à concrétions calcaires.
Une faille de direction N160, de pendage 70° est observable au milieu de la
tranchée. Le plan de la faille montre des stries de pitch 70°S qui indiquent un jeu
normal. Cette faille n’atteint pas la surface ; elle est scellée par les argiles brunâtres. La
couche sableuse montre des familles conjuguées de failles liées à la faille normale.
Le jeu de la faille normale a engendré la formation d’un demi graben rempli par
des dépôts plus récents vers la partie Est de la tranchée. En effet des argiles rougeâtres,
et des sables grossiers de couleur rouge d’âge plus récents, comblent le demi-graben.
Le jeu normal de cette famille de faille (N160° à N170°) au Pléistocène
supérieur montre une activité tectonique distensive.

81
82
Figure 39: failles normales au sein des dépôts Quaternaires sur le pied du versant Nord de l'anticlinal de Djebel
Ammar
Figure 40: Failles Pléistocène supérieur révélées par une tranchée sur le flanc Nord de djebel Ammar

83
Figure 41: Faille normale au Pléistocène sup. avec le réseau de failles secondaires associé.

84
 D-2-2 : La faille d’El Menzeh 9 (Ennasr)

Les collines d’el Menzeh et Ennasr bordent la ville du côté Nord-Ouest et sont
allongées selon la direction NW-SE. Ces collines culminent à des hauteurs de 85 et
75m NGT. Les formes et les dépôts quaternaires qu’elles portent sont relativement peu
développés par rapport à d’autres structures au NE du pays (Anticlinal d’Utique,…). Ils
se limitent la plupart du temps à des colluvions à dominance d’argile et de sables
rouges encroûtées, parfois même très fortement encroûtés.
Les combes marneuses d’âge pliocène du secteur d’El Menzeh 9 et El Manar sont
séparées des dépôts quaternaires de la zone d’El Mnihla par un escarpement
morphologique de direction N-S. A ce niveau nous avons observé un accident qui
apparaît par endroit. Il s’agit d’un décrochement senestre de direction N-S (figure 42),
affectant l’encroûtement feuilleté du Pléistocène supérieur (Dlala, 1994).

Figure 42: La faille N-S d’el Menzeh9 et de la colline de Sidi Bel Hassen

Au niveau des cités en construction sur la colline d’Ennasr (coordonnées :

36°51’21,20’’N ; 10°8’7,79’’E) nous avons découvert plusieurs failles dans les couches
pliocènes et qui touchent parfois la base du Quaternaire. Le rejet de ces failles est de

85
l’ordre de quelques dizaines de centimètres. Ces failles décalent les grès pliocènes par
un jeu inverse.
L’indice le plus spectaculaire que nous avons observé est une faille inverse affectant les
dépôts pliocènes et les argiles quaternaires. Cette faille est de direction N117°, 27°S,
16°N (figure 43). C’est une faille secondaire probablement associée à la faille N-S d’El
Menzeh9.
Les localités d’El Manar, Ennasr et Menzeh 9 sont aussi le siège d’une importante
activité de glissement de terrain qui est attesté par les fractures au niveau de la plupart
des murs et de construction et par le mouvement de terrains avoisinants les fondations.
La faille de direction N117, située au niveau du replat de l’anticlinal de Menzeh 9, est
différenciée d’un glissement de terrain car la surface ne montre pas de dénivellement et
la fracture n’atteint pas la surface.
Vue la densité importante de l’urbanisation au niveau des localités d’El Menzeh9, El
Manar, et Ennasr, il est difficile de trouver l’affleurement de la faille majeure de
direction N-S qui a été évoquée par Dlala (1994) et Dlala et al., (2008). Toutefois, les
failles inverses qui touchent le Pliocène et la base du Quaternaire que nous avons
observées peuvent représenter un potentiel seismogénique important dans ces localités.
Ce risque sismique est accentué par les conditions morphologiques de cette région
(effet de la topographie) qui seront détaillées dans le chapitre IV.

86
Figure 43: Faille quaternaire affectant le Pliocène et la base du quaternaire. Cette
faille est visibles au niveau des fondations des bâtiments à cité Ennasr

87
 D-2-3 : Tectonique Post-Tyrrhénienne affectant la falaise tyrrhénienne de
Sounine

Le long de la côte de la région de Sounine, une falaise tyrrhénienne émerge

d’une façon continue et elle est d’une altitude variant entre 5 et 7 mètres (figure 44).
Cette falaise tyrrhénienne est formée d’un faciès gréseux à ciment calcaire et riche de
« Strombus Bubonus ». Dans ces dépôts tyrrhéniens on observe plusieurs fractures et
failles inverses de directions variées. Ces accidents sont parfois des petits
décrochements de direction proches de N-S et qui sont visibles à la surface (figure 45),
le plan de fracture ne montre plus des indices sur le sens de mouvement vue l’énorme
influence des agents météorique et surtout les vagues de la mer. Le sens de mouvement
sénestre indiqué sur la figure 46 est déduit d’après le contexte compressif général de
direction proche du NW-SE. On observe également des failles inverses spectaculaires
qui affectent toute la hauteur de la falaise de direction E-W (figure 46). Ces
observations témoignent d’une activité tectonique compressive post-tyrrhénienne sur
les côtes nord orientales tunisiennes.

Figure 44: Localisation des indices de tectonique active au niveau de la corniche tyrrhénienne
(image de Google Earth)

88
 Figure 45: Une faille inverse affectant les dépôts Tyrrhénien de la plage de Sounine.
Le Tyrrhénien atteint une altitude de 8m dans cette localité

Figure 46: Décrochement sénestre dans les dépôts tyrrhéniennes de la falaise de Sounine;
1: Tyrrhénien; 2: Dépôts calcaire sur le plan de décrochement

89
 D-3-2 : Tectonique Post-Tyrrhénienne affectant la Plage d’Errimel

Sur la côte de la zone de camping Errimel, une corniche tyrrhénienne longe par endroit
la côte et est parfois submergée (figure 47). Ces vestiges du Tyrrhénien ont un faciès
gréseux fortement cimenté par un ciment calcaire et ont un aspect spongieux. La
présence de Strombus Bubonus est un marqueur de l’âge Tyrrhénien de ces
affleurements (Jedoui et al., 2002; Jedoui et al., 2003).
Ces affleurements forment un anticlinal de direction E-W, qui a le cœur érodé, et qui
est haché par une famille de fracture de direction proche d’E-W (figure 48).

Figure 47: Photo de la corniche Tyrrhénienne du Plage Errimel. On voit le flanc sud
d’un anticlinal Tyrrhénien à large rayon de courbure

Du point de vue stratigraphique la plage d’Errimel renferme des dépôts peu épais du
Tyrrhénien qui reposent en discordance sur des grès à Pectinidés du Pliocène supérieur
et à faune abondante (Piementa, 1959; Ben Ayed et al., 1978).
Les flancs de cette structure, ont des pendages de l’ordre de 15 à 20°. Sur le flanc nord,
des failles de directions E-W sont visibles (figure 49). Une étude antérieur (Ben Ayed
et al., 1978) a montré que ces failles ont des pendages variant entre 10° à 40°, et portent

90
des stries plongeant d’environ 70° vers le SE et montrant un mouvement inverse. Le
jeu total est de 0,6m (Ben Ayed et al., 1978)

Figure 48: Coupe simplifié de l’anticlinal d'Errimel. Elle montre le plissement des dépôts marins du
Quaternaire récent (Ben Ayed et al., 1978)

Figure 49: Photo montrant des failles E-W qui jalonne les couches tyrrhéniennes plissées
de la corniche d'Errimel

91
 D-4 : Anticlinal de Messeftine

D-4-1 : Géologie

L’anticlinal de Messeftine est situé à l’Est de Mateur et au SE du Lac Ichkeul. C’est un

grand anticlinal irrégulier à cœur mio-pliocène (figure 50). Les travaux de Dlala (1995)
ont montré que l’axe de pli a subi une torsion d’une direction NE-SW à la terminaison
périanticlinale NE, vers une direction N-S au niveau de la 2 ème terminaison
périanticlinale SW. Cette rotation a été comme étant la conséquence d’un jeu sénestre
de la faille N-S de Messeftine.

Figure 50: Vue panoramique de l'anticlinal de Messeftine (image de Google earth)

92
Le cœur de la structure de Messeftine, occupé par des marnes grises, rouges ou noires
riches en gypses d’au moins 400m d’épaisseur, correspondrait au Miocène inférieur
bigarré (Formation Hakima). Le reste de la série miocène, correspond à des marnes
grises riches en bancs de gypses et en niveaux bréchiques, est attribué à la formation
Oued Bel Khedim (Burrolet, 1951) (figure 51). La série de marnes à l’intercalations de
bancs de gypses correspond à la formation Oued El Melah (Burrolet, 1951). Le cœur de
la structure anticlinal du Djebel Messeftine est occupé par au moins 250m d’argiles
gris-sombre, gypseuses, avec des intercalations de couches centimétriques de grès fins
à la base. La série se poursuit par plus de 100m d’une alternance d’argiles grises et de
bancs gypseux d’épaisseur métrique. Des niveaux centimétriques de calcaires lacustres
blancs existent dans les argiles entre les bancs de gypse.
Les secteurs NE et Sud de l’anticlinal sont occupés par les marnes du Pliocènes
(formation Raf-Raf). Sur le flanc Est, les grès pliocènes dessinent une auréole régulière
marquant une fermeture sans accident

Figure 51: Carte géologique de l'anticlinal de Messeftine

93
 D-4-2 : La faille de Messeftine

Sur les sections sismiques cette faille est visible comme un accident de direction N-S
qui s’étend de la région de Menzel Bourguiba au Nord jusqu’à la région de Lansarine
vers le Sud (Kacem, 2004) (figure 23). Elle montre une rupture topographique entre
l’anticlinal de Messeftine et la plaine de Mateur. La géométrie de cet accident montre
qu’il correspond à une structure en fleur caractérisée par la présence de plusieurs
ramifications qui témoignent d’un jeu décrochant.
Cette faille, de direction N-S, s’étend de la région de Menzel Bourguiba au Nord
jusqu’à la région de Lansarine vers le Sud. Elle représente la séparation entre
l’anticlinal de Messeftine et la plaine de Mateur. Les plis d’entrainement au niveau du
quaternaire ancien témoignent d’un jeu décrochant sénestre de cette faille (Dlala,
1994).
Les brèches tectoniques connues au sud de Djebel Assafir, comme la source qui s’y
trouve (Aïn Baddad) sont les indications d’une cassure importante. Les stries sub-
horizontales observées au niveau des sites microtectoniques témoignent d’un jeu
sénestre de cet accident, ce mouvement se déduit également de la répartition de plis en
échelon quaternaires, que l’on peut observer facilement dans ce secteur (Dlala, 1994)

94
 D-4-3 : Indices de tectonique actives

a/ indice de Oued El Meleh

Au niveau de la rupture topographique entre l’anticlinal et la plaine de Mateur, les
terrasses alluviales qui dominent de quelques mètres le lit majeur de l’Oued El Melah
(37.0725°N ; 9.72271°E), sont basculées et montrent un contact anormal avec les
argiles bariolées du miocène du Oued El Melah (figure 52).

Figure 52: Discordance des silts Quaternaires sur les argiles verdâtre du
Miocène (vue vers le NE) (Mejri, 2004)

b/ indice de Djebel El Assafir

Plus au sud, sur la crête de la colline de Djebel El Assafir, dans une carrière de sables
pliocènes, maintenant délaissée (37.0285°N ; 9.73329°E), on observe une faille inverse
de direction N-S est presque verticale. Elle coupe tout l’affleurement du grès Pliocène

95
et se perd lorsqu’elle entame les silts quaternaires (figure 53). Il est difficile de trouver
des stries préservées dans les grès. Cette faille qui peut être associée à la grande faille
de Messeftine ce qui donne la conviction de sa grande ampleur.

Figure 53: Failles normales synsédimentaires pliocènes au niveau de Djebel El

Assafir

96
 c/ indice de la Ferme Lebray

Au Sud de la ferme « Lebray », près de l’Oued El Melah (37.0748°N ; 9.72204°E), des

bancs durs de Miocènes sont redressés à la verticale et par fois à pendage redressé, ce
qui pourrait âtre lié au jeu d’une faille. Le jeu Plio-quaternaire de ces accidents peut
être rattaché à la phase compressive la plus importante d’âge Quaternaire
(Villafranchienne) (figure 54).

Figure 54:Contact tectonique entre la croute calcaire Quaternaire et les formations Miocène

D-4-4 : Cartographie du Quaternaire sur le versant Ouest de

l’anticlinal

Nous avons essayé de faire une cartographie détaillée des faciès quaternaires sur le
piémont du versant Ouest de l’Anticlinal de Messeftine (figure 55). Nous avons
pris comme repère la barre encroutée du quaternaire récent et qui affleure dans la
majorité du secteur balayé. Nous avons observé les encroûtements calcaires du
Quaternaire moyens et récents affleurent à des altitudes variables sur le piémont se
la rive Ouest de l’anticlinal. Au niveau de l’Oued El Melah, cette faille met en
contact les argilites noirâtres de la formation miocène Oued Melah avec les
alluvions récents du Quaternaire. Le contact est souligné par un passage net de

97
faciès sans rupture mécanique qui pourrait renseigner sur la composante latérale
du mouvement.
Au niveau d’une ancienne carrière, maintenant délaissée, les couches bigarrées de
la formation Hakima du Miocène supérieur sont redressées à la verticale. Les
excavations de la carrière ont enlevé la croute calcaire qui n’en reste que quelques
vestiges, légèrement incliné vers l’Ouest.
Les silts et limons quaternaires généralement de couleur grisâtre caractéristique des
silts du Pléistocène supérieur, dominent les piémonts du versant Ouest à jusqu’à
des altitudes peu élevées et qui ne dépassent guère les 50mètres. En bas de ces
silts gris, on rencontre, par endroits des silts de couleur plus claire (beige) riches en
concrétions calcaires. Nous attribuons ce faciès au quaternaire récent. Le pied du
versant est rattaché aux alluvions récentes de la plaine de Mateur.

98
Figure 55: Cartographie de Djebel Messeftine et de la faille de Messeftine

99
 E : Interprétation et conclusion

La phase compressive post-villafranchienne évoquée par plusieurs auteurs est mise en

évidence par les multiples indices de tectonique récente qui se manifeste généralement
par des mouvements en faille inverse affectant les formations pliocènes et les dépôts
quaternaires parfois les plus récents. La majorité des mécanismes au foyer réalisés sur
les séismes enregistrés au Nord de la Tunisie, plaident en faveur d’une compression
actuelle orientée NNW-SSE (figure 56). Les axes de pressions ont une direction
moyenne NW-SE.
Les indices de tectonique quaternaire que nous avons observés sont répartis sur la
région du grand Tunis et la région de Bizerte (figure 57). A Tunis les collines des cités
Ennasr et Menzeh9 montrent des failles inverses qui décalent les grès pliocène et
parfois les dépôts quaternaires. Nous rattachons ces failles inverses à l’activité récente
de la faille N-S qui s’étend depuis le nord de Sebkhet Essijoumi jusqu’aux collines de
Mnihla au Nord (Dlala & Kacem, 2008). La faille de Dj Nahli montre une activité
récente à travers les micro-failles et plis affectant les encroûtements calcaires du
quaternaire, aux abords de l’autoroute de Bizerte. La faille de Dj Ammar affecte les
dépôts du Quaternaire récents et forme une rupture topographique entre le Dj. Ammar
et la plaine d’El Mabtouha.
La tectonique récente est évidente au niveau du versant sud de l’anticlinal d’Utique, à
la limite entre le Pliocène qui forme le cœur du pli et les dépôts quaternaires qui
couvrent une grande partie du pied du versant. En effet, l’affleurement par endroits de
la faille d’Utique mise en évidence par les sections sismiques se manifeste par des
micro-failles décrochantes à composante inverses de directions proches d’E-W. Le plan
de la faille majeur se manifeste souvent sous forme d’une zone de broyage qui mêle les
fragments de grès Pliocène avec les débris de la croute calcaire et les silts du
Quaternaire.
Les failles et les plis tyrrhéniennes, que nous avons observés sur les côtes de Sounine
et du Camping d’Errimel, sont causés par une phase compressive actuelle compatible
avec le cadre géodynamique de la Tunisie dans un contexte de convergence Afrique –
Europe.

100
L’importance d’indices de déformation cités ci-dessus et leurs relations avec les
accidents tectoniques majeurs, permettent de classer les failles actives selon leurs
potentiels seismogénique (figure 57). En effet, les failles au niveau des structures
d’Utique, de Messeftine de Djebel Nahli, de Djebel Ammar, et de Ras Djebel
représentent des sources seismogénique à potentiel à prendre en considération.
L’activité des formations tyrrhéniennes au niveau des côtes Nord-Orientales
témoignent d’une activité tectonique récente qui justifie un raccourcissement actuel.
Dans la cadre de ce travail, nous discuteront le cas du pli-faille d’Utique, afin d’évaluer
sa capacité seismogénique, à travers une approche basée sur la recherche des données
tectoniques récentes, ainsi que des données paléosismologiques et de sismicité
historique. L’outil géophysique sera aussi utilisé pour estimer le taux de
raccourcissement sur cette faille.

101
Figure 56: Schéma des répartitions des indices de la déformation récente et les mécanismes
dont ils sont associés en Tunisie Nord-orientale

102
Figure 57: Les indices de déformation récente au niveau de la région de Tunis et du
NE de la Tunisie

103
 CHAPITRE IV

ETUDE DE LA FAILLE ACTIVE

D’UTIQUE

104
A : Cadre géomorphologique, paléogéographique et géologique de la basse
Medjerda et la région d’Utique

A-1 : Cadre géomorphologique et géologique

La plaine de la basse Medjerda est encadrée par une topographie de Djebels et collines
qui, malgré leurs altitudes plutôt modestes puisqu’elles ne dépassent que rarement les
400m, sont toujours très marqués dans le paysage. Ils doivent ceci surtout à leurs
pentes, parfois fortes et contrastant avec la remarquable platitude de la plaine qu’ils
dominent directement, parfois sur plus de 200 m d'un seul trait (figure 58). L’anticlinal
de Messeftine et celui d’Utique illustrent bien ce contraste morphologique avec les
plaines de Mateur et de Garâat El Mabtouha.
La plaine de la basse vallée de la Medjerda et la plupart des reliefs qui l'encadrent ou
qui la ponctuent, ont une géologie jeune (figure 59). Elle est occupée par des alluvions
quaternaires le plus souvent argilo-sableuses à argilo-sablo-limoneuses. Dans sa partie
la plus proche de la mer on a surtout affaire, sur une bande parfois large de plus d'une
dizaine de kilomètres, à des formations alluviales accumulées, toujours par l’Oued
Medjerda, au cours des deux derniers millénaires. Quant aux reliefs des bordures Nord
et Sud, ils correspondent à des structures anticlinales mais montrent des différences
sensibles d'un secteur à l'autre. Du côté septentrional, on a affaire à un alignement de
plis (Djebel Ennadour, Djebel Messeftine,…) à ossature dominée par des formations
mio-pliocènes. Les versants donnant directement sur le delta sont le plus souvent des
revers de monoclinaux du type crêts façonnés dans une alternance de grès et de marnes
pliocènes.
Les reliefs méridionaux ont une géologie plus ancienne et des structures plus
complexes. Les plus importants sont Djebel Ammar et Djebel Nahli qui le prolonge
vers le Nord-Est. Ces Djebels correspondent à des dômes à structure anticlinale plus ou
moins complète et à lithologie parfois très variée mais assez marquée par des
formations calcaires ou argilo-calcaires crétacés. Le Djebel Ammar se distingue par son
diapirisme et une importante faille de direction NE-SW injectée de matériaux argilo-
gypseux et accompagnée par une lame de calcaires jurassiques.

105
Enfin, les collines septentrionales qui ponctuent la plaine appartiennent à deux types de
géologies. Les unes correspondent à des structures plissées (collines de Menzel Ghoul-
Utique) ou à des blocs soulevés du type horst (la colline qui porte la ville de Kalaat
Landlouss) et ont une ossature de grès et d'argiles pliocènes. Les autres, plus
nombreuses comme celles de Pont de Bizerte (traversées par l'autoroute Tunis-Bizerte)
ou de Koudiat El Mabtouha correspondent à d'anciennes dunes d'argile du type lunettes
d’âge pléistocène supérieur.

Figure 58: Vue panoramique de l'anticlinal d'Utique (vue vers le Nord)

106
107
Figure 59: Extrait de la carte topographique (1/25000 ; équid. 5m) qui montre la
structure d'Utique au milieu de la plaine d'El Mabtouha
 A-2 : Cadre paléogéographique

L’étude paléogéographique de la plaine deltaïque de la Medjerda montre

qu’elle occupe la place d’un ancien golfe marin, le sinus Uticensis ou baie
d’Utique – au sens restreint du terme car, au sens le plus large, cette expression
englobe le golfe de Tunis. Le comblement de cet ancien golfe est le résultat des
divagations du cours du fleuve depuis l'Holocène. Le rivage a reculé au point que
les ruines de l’illustre cité d’Utique, qui était un port aux époques punique et
romaine, se retrouvent à présent, à vol d’oiseau, à plus de 10km à l’intérieur des
terres. Il s’agit d’une progression de grande ampleur puisque 450km2 ont été
ainsi gagnés, au dépend de la mer, depuis la fondation d’Utique, il y a près de
3000 ans. Cet exemple de progression est d’autant plus digne d’être signalé que
prévaut, ailleurs, une tendance inverse, avec recul de la ligne de côte, comme
l’ont montré récemment les premiers résultats d’une enquête pluridisciplinaire
tuniso-française sur le littoral de la Tunisie.
Une prospection archéologique et géomorphologique récente de
l'ensemble de l'ancien golfe a permis de recueillir des céramiques datées sur une
cinquantaine de site. Celles-ci permettent de conclure le comblement précoce du
couloir en question, au sud-est d'Utique. La recherche des installations portuaires
antiques n'est donc envisageable que sur le côté nord de la presque-île. Une
occupation pré-romaine est aussi attestée sur les lunettes d'argiles d'origine
éoliennes développées dans la région de Pont de Bizerte et dont l'existence même
suppose un alluvionnement - sans doute antérieur à la fondation d'Utique - de
toute la partie centrale de la plaine alluviale. Celui-ci est à mettre en relation en
amont, avec le bras de Chaouat trace d'un ancienne passage du fleuve au nord du
Djebel Maiana (Chelbi et al., 1995).
Tous les auteurs s’accordent pour estimer que le compartiment Sud a été
colmaté en premier par les divagations du sud vers le nord du Bagradas
(Medjerda). Au nord, le comblement est resté inachevé puisque subsiste, dernier
vestige de l’ancienne baie d’Utique, la lagune de Ghar el Melh (figure 60).
La structures plissée d’Utique a été décrit précédemment par plusieurs
auteurs (Burollet, 1951; Rouvier, 1971; Viguier et al., 1980; Ben Ayed &

108
Viguier, 1981; Ben Ayed, 1986) comme l’une des structures d’âge miocène
supérieur.
Les travaux récents (Dlala 1994 ,95 ; Kacem 2004; Dlala, 2008 ; Mejri, 2004)
montre qu’il s’agit plutôt d’une structure plicative d’âge quaternaire, moulée sur
une faille préexistante de direction E-W bien identifiée en subsurface. Cette
interprétation a été fondée sur la mise en évidence d’une tectonique distensive
d’âge pliocène affectant tout le bassin moi-pliocène du Kechabta-Messeftine.
Cette interprétation montre également (un problème de permutation régionale des
contraintes (Dlala 94 et 95). La variation d’épaisseur dans la série du Miocène
supérieur et le Pliocène est bien identifiée également par les forages pétroliers
(Kacem, 2004). Pour essayer de prouver l’activité quaternaire de cette structure,
nous allons étendre l’étude sur les aspects morphologiques, de subsurface, et de
l’évolution et indice d’activité récente.

Figure 60: Les étapes de comblement de l’ancien golfe d’Utique (Jauzein,1971)

1. Djebels et collines, [Link], [Link] Medjerda, [Link] mort de l’Oued,
[Link] et garaas, [Link] d’argiles

109
 A-3 : Géomorphologie de l’anticlinal d’Utique

Nous avons procédé à un examen des photos aériennes à une échelle de

1/10000. En outre les cartes géologiques et topographiques sont également
disponibles sur la région d’étude. La compilation de ces outils nous a permis
d’élaborer un modèle numérique de terrain de l’anticlinal d’Utique qui montre un
anticlinal asymétrique de direction presque Est-Ouest. En effet, des profils
topographiques transversaux à la direction de l’anticlinal nous montrent que la
pente du flanc Nord est douce et ne dépasse guère les 30°, alors que celle du flanc
sud, aux abords de la plaine d’El Mabtouha est plus raide et elle est aux alentours
de 60° (figure 61).
Ces observations permettent de classer, selon le sens de Suppe (Suppe &
Medwedeff, 1990) l’anticlinal d’Utique comme un pli en propagation de rampe
(fault-propagation fold).(figure 62).
Le réseau hydrographique présent sur les flancs de l’anticlinal n’est pas
très développé, il est formé par des ravins entaillés dans les argiles sableuses
quaternaires et ils atteignent parfois les grès pliocènes à des altitudes élevées. Ce
réseau ne montre pas d’indice de décalages liés à des décrochements récents.
Des profils topographiques transversaux à l’anticlinal montrent que le
versant Sud de l’anticlinal présente une rupture topographique observable à des
altitudes variables entre 70 et 120 mètres. Cette rupture topographique peut être
d’origine tectonique et nous le rattachons au jeu de la faille d’Utique reconnue en
subsurface.
Dans le cadre de ce travail, nous allons également suivre cette rupture
topographique par une cartographie détaillée des dépôts Quaternaires. Des profils
topographiques à hautes précisions seront aussi réalisés afin de mettre en
évidence la rupture topographique.

110
111
Figure 61: MNT de l’anticlinal d’Utique et des profils topographiques transversaux à la structure
Figure 62: Modèle schématiques de (a) plis de décollement et (b) plis de propagation
de faille (Suppe & Medwedeff, 1990)

112
 A-4 : Géologie de l’anticlinal d’Utique

La plaine d’El Mabtouha, est occupée par des alluvions Quaternaires le plus
souvent argilo-sableuses à argilo-sablo-limoneuses. L’anticlinal d’Utique
proprement dit, montre une ossature dominée par des faciès gréseux du Pliocène
marin (formation Porto-Farina). Les versants Nord et Sud donnant directement
sur les plaines bordières, montrent des affleurements quaternaires visibles jusqu’à
des altitudes importantes sur les rives.
Le Pliocène est un ensemble à dominance sablo-gréseuse de teint jaunâtre à ocre
et elle est riche en fossiles (Pectinidés et Ostréidés). Ils sont assez épais, parfois
assez grossiers. Dans certains secteurs de la région étudiée, des études en lames
minces (Burrolet, 1951), ont montré que ces grès sont quartzeux à ciment
calcaire, à nombreux tests de lamellibranches ou plus rarement de gastéropodes et
de très rares foraminifères.
Nous avons réalisé une coupe au niveau d’un vallon sur le flanc Sud de Djebel
Menzel Ghoul. Cette coupe montre, de haut en bas du relief, la succession
suivante (selon la nomenclature sur la figure 63):
- (1) Une couche gréseuse du Pliocène qui représente le cœur de l’anticlinal.

- (2) une couche de silts verdâtres à concrétion calcaire. En comparaison

avec les travaux antérieurs (Dlala, 2001) sur la zone du Grand Tunis, cette
couche peut être attribuée classiquement au Pléistocène supérieur.

- (3) une couche de limons riches en concrétions calcaire. Elle est de

couleur brunâtre.

- (4) une croute calcaire feuilletée, attribuée au Pléistocène supérieur.

- (5) une couche mince de dépôt carbonaté pulvérulent.

- (6) Une couche de limons rouges riche en Hélix (escargot). Elle est de
couleur rouge vif et d’âge Holocène à Actuel.

113
 Figure 63: Coupe des formations quaternaires sur le versant Sud de l'anticlinal d'Utique

B : Etude de Section sismique

Une carte tectonique de subsurface couvrant la région du NE de la Tunisie

a été réalisée (figure 51) (Kacem, 2004). Cette carte montre des failles majeures
orientées principalement selon trois directions : E-W, N-S et NE-SW. Le versant
Sud de la structure d’Utique est bordé par une faille de direction presque E-W de
longueur presque 25 km et de profondeur 5 km.
Un profil sismique de la campagne MA-83 réalisée en 1983 par la Société
de Service (CGG) pour le compte de la Société Shell Tunisienne de Recherche et
d’Exploration (Shell Tunirex), s’étend de la plaine de Kechabta au Nord de la
structure d’Utique jusqu’au Djebel Ammar au Sud en traversant le pli d’Utique et
la plaine d’El Mabtouha. Un travail préliminaire à été exécuté sur le profil
sismique afin de convertir l’échelle verticale originellement en temps double,
vers une échelle métrique via les paramètres de vitesses et le calage sur la
colonne stratigraphique du forage Ut1.
Sur cette section on observe une structure plissée qui correspond en
affleurement à l’anticlinal d’Utique (figure 64). Cette section sismique montre
qu’au niveau de la structure d’Utique, un plan de décollement majeur décale en
faille inverse la série mio-Pliocène et la base du Quaternaire. Ce niveau de
décollement qui correspond aux faciès argilo-évaporitique de la formation Oued

114
Melah, atteint la surface du versant Sud de l’anticlinal et s’exprime en faille
inverse qui touche les niveaux quaternaires. Nous avons trouvé des indices de
cette activité quaternaire dans plusieurs localités d’affleurement de cet accident.

Figure 64: Profil L1 (mission MA83) passant la structure d’Utique

(on voit seulement la partie de la section qui couvre le pli d’Utique)

B-1 : Coupe géologique équilibrée

La méthodologie de l’équilibrage d’une coupe géologique repose sur le

géoréférencement en 3 dimensions de l’espace (x, y et z) des données d’une
section géologique (limite des couches, failles, pendage, altitude,…) pour
contraindre la cohérence tridimensionnelle de la coupe. Les données
géoréférencées seront ensuite traitées par un logiciel d’analyse et de modélisation
structurale (Midland Valley 2DMove). Les étapes de la construction d’une coupe
géologique équilibrée sont :
- Digitalisation des contours géologiques et des failles
- Digitalisation des pendages
- Exportation des données digitalisées vers le logiciel de modélisation
- Projection des données (pendages, failles,…) sur la section-line et collecte
des intersections

115
 - Transformation en section et construction de la coupe équilibrée (figure
65)

Figure 65: Exemple d’équilibrage d’une coupe géologique (fenêtre du logiciel de

modélisation 2DMove)

B-2 : Equilibrage de la section géologique à travers la structure

d’Utique

Du point de vue structural on s’intéresse uniquement à la structure

anticlinale d’Utique afin de reconstruire son histoire géologique depuis le Trias.
Nous allons équilibrer la coupe géologique que nous avons construite d’après la
section sismique calée par les données du puits Ut1 et par les données du terrain.
En effet, la couverture sédimentaire post-Miocène supérieur repose en
discordance sur une surface d’érosion mise en évidence à la fois par la troncature
des réflecteurs correspondant au Trias et qui ont un pendage le vers le Nord, et
par l’apparition d’une structure en « onlap » au niveau des formations du
Miocène inférieur. Les niveaux Miocène inférieur sont affectés par un niveau de
décollement au niveau de la formation évaporitique « Mellaha » (figure 66).

116
 La section sismique a été analysé en s’appuyant sur des concepts
d’analyse tectonique (Dahlstrom, 1969 ; Dahlstrom, 1990 ; De Sitter, 1964 ;
Jamison, 1987) pour pouvoir proposer une coupe équilibrée de la structure
d’Utique. Cette coupe montre que l’anticlinal d’Utique est formé suite à la
réactivation d’une structure plissé anté-Miocène. En effet, la phase de
raccourcissement responsable de la formation du pli d’Utique est une phase
compressive post-Miocène supérieur. Cependant, il est difficile de déterminer
exactement le début de cette phase de compression.
Durant cette dernière phase de raccourcissement post-Miocène supérieur
le niveau de décollement que nous avons observé au niveau de la formation
évaporitique « Mellaha » a permis l’accommodation de la déformation
compressive. Ce système tectonique basé sur l’accommodation de la compression
par le glissement des formations sur un niveau de décollement montre l’aspect
du système en « fault propagation fold » (Mejri et al., 2010).
Ce scénario nous permet de calculer un raccourcissement total, depuis le
Mio-Pliocène, estimé à 275 m. étant donné que cette déformation est post
Miocène-supérieur (il y a 5 Ma), les valeurs proposées permettent de calculer une
vitesse de raccourcissement d’environ 0,055 mm/an (275m/5000000ans) depuis
le Pliocène à l’Actuel. Or, au niveau de la Tunisie Nord-orientale, la dernière
phase de raccourcissement importante après la phase compressive tortonienne est
la phase post-villafranchienne (depuis le Pléistocène moyen à l’Actuel, soit 1,8
Ma) (données et analyses proposés au niveau du chapitre I de cette thèse), ce qui
suggère fortement que la totalité de raccourcissement s’est effectuée durant cette
dernière phase. Ceci nous permet de déduire que la vitesse de raccourcissement
au niveau de la structure d’Utique est d’environ 0,15 mm/an (275m/1800000ans)
(Mejri et al., 2010). (la figure 5 de cet article présente une erreur au niveau de
l’échelle. Cette erreur est corrigée sur la figure 66 de ce document)
Ce taux de raccourcissement est relativement important (3 %) par rapport
au raccourcissement total au niveau de la Tunisie d’environ 4,5 mm/an selon une
direction N145° entre l’Afrique stable et le Sardaigne (Hollenstein et al., 2003;
D’Agostino & Selvaggi, 2004; Nocquet & Calais, 2004), et permet de classer la
structure d’Utique comme une structure à potentiel séismogénique important.

117
Figure 66: Interprétation du profil L1 et estimation du rejet de la faille d'Utique
(Mejri et al., 2010)

118
C : Les données de la tectonique quaternaire au niveau de la structure
d’Utique

 Dans le versant Sud de la partie Ouest de l’anticlinal d’Utique, au niveau

d’une vallée de carrière de sable délaissée (coordonnées : 37° 1'58.66"N ;
10° 0'16.27"E), nous avons observé une faille inverse qui affecte les bancs
gréseux du Pliocène et la base des argiles quaternaire (quaternaire récent)
(figure 67). Cette faille inverse est de direction NE-SW et de pendage 30°
vers le NW. Au niveau du plan de la faille apparaissent des stries de pitch
28° SW.

Site d’observation
d4oservation

Figure 67: Faille inverse qui affecte le Plio-quaternaire au niveau du versant

Sud du pli d'Utique (A : localisation de la faille, B : faille au niveau de la
carrière abondonnée)

119
 Vers la partie ouest de l’anticlinal, des carrières de sable ont été ouvertes au
sommet du pli dans les formations gréseuses pliocènes lors de la construction
de l’autoroute Tunis-Bizerte. Dans ces carrières nous avons observé des
failles presque verticales et parfois présentant une composante normale
(figure 68). Ces failles de petite ampleur sont localisées au niveau de la
charnière de l’anticlinal. Elles décalent les grès pliocènes et la base du
Quaternaire. La présence de ces failles normales dans un contexte compressif
est due à la localisation de ces failles sur la charnière du pli et à la courbure
assez importante de l’anticlinal. Les failles normales résulteraient donc d’une
distension locale limitée à l’extrados du pli.
S N

Figure 68: Failles normales associées au plissement quaternaire de la série plio-quaternaire

au niveau de la crête de l'anticlinal d'Utique

120
D : Application d’une approche radiochronologique

Radiochronologie des fossiles (Hélicedes) dans les silts rouges :

Dans un vallon et à quelques cinquante mètres à droite d’un petit barrage

récemment construit, nous avons échantillonné des coquilles d’Hélicedes
présentes dans les silts rouges qui sont sous-jacents à la croûte calcaire feuilleté.
La datation par le radiocarbone (C14) de ces coquilles n’a pas délivré de résultat.
Ceci est expliqué par la tranche d’âge mesurable par cette méthode de datation
qui ne dépasse pas les 50 ka. Nous avons eu recours à la méthode de datation par
le couple U/Th.
Le principe de la datation U–Th repose sur les déséquilibres radioactifs entre le
230 234 238
Th (demi-vie égale à 75.38 ka), le U (demi-vie égale à 245.6 ka) et le U
dont la demi-vie est de 4.47 Ga. En milieu clos l’évolution de ce système se fait
vers un état d’équilibre séculaire atteint lorsque le nombre de désintégrations par
unité de temps (l’activité) devient égale et constante pour chaque radionucléide
(Sharp et al., 2003).
Les âges U-Th obtenus par cette méthode sont considérés équivalents aux âges
calendaires (Bard et al., 1990), du fait d’une bonne correspondance entre les âges
U-Th des coraux et les enregistrements dendrochronologiques. La méthode de
datation par les séries de l’uranium présente donc plusieurs avantages,
notamment par rapport à la méthode de datation au C14 : elle permet d’obtenir
des âges calendaires, de dater des échantillons jusqu’à ~500000 ans,
La méthode de datation U–Th va donc reposer sur la croissance in situ au
230 234
cours du temps du Th et la décroissance du U vers l’équilibre séculaire
230
(Ludwig & Paces, 2002). Si l’échantillon analysé ne contient pas de Th initial,
l’âge peut alors être calculé d’après l’équation suivante (Kaufman et al., 1971) :

121
Cette équation suppose que le système est clos. En réalité, et dans un milieu
humide et suffisamment perméable pour faciliter l’écoulement des eaux, il se
produit une mobilisation de l’Uranium après incorporation diagénétique
d’uranium après la mort de l’individu. Or, dans notre cas la couche
encaissante des mollusques est argileuse et imperméable à l’écoulement des
eaux, ce qui ne permet pas un apport d’uranium externe.
5 analyses de l’uranium 238 ont été réalisées sur le même échantillon. Ces
analyses ont montré que la quantité de l’uranium 238 est presque la même
dans les coquilles récoltés dans la même formation argileuse. La moyenne de
ces analyses est donnée dans le tableau ci-haut.
Ces deux conditions permettent de considérer que le système est un système
fermé et que l’équation de Kaufman est applicable est permet de donner des
âges suffisamment précis.
Après le processus radiochimique de séparation des différents isotopes de
l’Uranium et du Thorium, le comptage des isotopes présents dans
l’échantillon à été réalisé par une chaîne de comptage de radioactivité Alpha.
Les résultats sont présentés dans le tableau suivant :

Isotope Comptage (Bq/kg)

234
U 16,82
Uranium 238
U 17,11
228
Th 61,91
230
Thorium Th 32,79
232
Th 23,594
Tableau 4 : résultat de comptage des isotopes d’Uranium et de Thorium dans
l’échantillon analysé

La résolution de l’équation présentée ci-haut a été représentée sur un abaque

sur lequel on peut déterminer les âges à partir des rapports isotopiques
230 234 234 238
Th/ U et U/ U (figure 69).

122
En utilisant l’équation ci-dessus et le digramme de Concordia nous avons
conclu à un âge de notre échantillon qui est de l’ordre de 220 ka ans ± 12 ka.
Donc nous avons attribué cette formation de silts rouge au Pléistocène moyen.
Ce géochronologie absolu permet d’affirmer avec une bonne précision l’âge
de cette formation qui a été attribué par d’autres auteurs au Pléistocène
supérieur (Oueslati, 1993; Dlala & Kacem, 2008).
Les indices de déformations récentes observés au niveau de cette formation
sont le résultat d’une phase de compression qui se poursuit après le
Pléistocène supérieur et qui est probablement active jusqu’à l’actuel.

Figure 69: Diagramme Concordia pour la résolution des âges de manière graphique

en calculant les rapports isotopiques (Kaufman et al., 1971)

123
E : Mesure de l’inclinaison du niveau de silts à poupées calcaires

L’horizon de silts rouges où nous avons échantillonné l’Hélicedes pour la

radiochronologie (figure 70) contient des nodules calcaires développés autour des
racines des plantes. A l’origine, et selon le principe de gravité les nodules sont
naturellement verticales quelques soit l’inclinaison de la couche encaissante.
Dans le but de confirmer l’inclinaison des horizons Quaternaires, nous avons
procédé à une approche permettant de vérifier la possibilité de mettre en évidence
la déformation récente. En effet, nous avons mesurés l’orientation de ces
concrétions calcaires dans deux sites différents sur le flanc sud de l’anticlinal.
Ces mesures montrent une possible faible inclinaison (inférieure à 5°) par rapport
à un axe perpendiculaire au flanc de l’anticlinal depuis le Pléistocène supérieur.
En effet, sur la partie A (figure 71) il n’y a pas de basculement visible alors que
sur la partie B un léger basculement de l’ordre de 3° dans une direction N290° est
possible (avec une confiance de 95%), sans être certain ; cette direction est
quasiment perpendiculaire à l’orientation du pli.

Figure 70: Poupées calcaires dans la formation du silt rouge du quaternaire

au niveau du piémont du versant Sud de l'anticlinal d'Utique

124
Figure 71: Stéréogramme (projection de Wulf) de l’orientation des concrétions calcaires à
2 endroits sur le flanc sud-est du pli.

125
F : description des différents sites d’observation au niveau de la faille
d’Utique

Dans cette partie nous allons décrire les observations de la tectonique récente
attribuées à l’activité de la faille d’Utique au niveau du flanc Sud de l’anticlinal
d’Utique. Les localisations des sites d’observations sont représentées sur la figure
72.

Figure 72: Localisation des sites d'indices d'activité de la faille d'Utique (S1....S6); (Image de
Google Earth, et courbes topographiques de la carte topographique de Ghar El Melh 1/25000,
équi. 5m)

 Site1 :
Dans ce site (37,032°N ; 10,023°E) nous avons observé une faille à pendage
presque verticale qui met en contact anormal le Pliocène avec les argiles brunes
du quaternaire ancien. Cette faille montre une direction presque E-W et présente
un plan de faille presque vertical (figure 73). La zone faillée présente une zone
broyée où se mélangent les débris de grès avec les concrétions calcaires et les
argiles brunes du Quaternaire.

126
 Figure 73: Affleurement du plan de la faille d'Utique

 Site 2 :
Au niveau d’un petit vallon profond d’à peu près 1,5 mètre (37,029°N ;
10,001°E), nous avons observé une faille inverse de direction E-W à
l’affleurement. En effet, nous avons observé un plan de faille de direction N70 et
de pendage 39° vers le Nord qui décale les grès pliocène pour les mettre en
contact avec les silts rouge du Quaternaire récent. Nous avons localisé cet
affleurement à une altitude de 40 mètres sur le flanc Sud de l’Anticlinal d’Utique
(figure 74).

127
 Figure 74: Jeu inverse de la faille d'Utique: contact
tectonique des grès Pliocènes sur les silts brunâtres du
Quaternaire

 Site 3 :
Sur le flanc Sud de l’anticlinal d’Utique, au niveau d’un vallon profond
d’environs 3 mètres (37,029°N ; 10,002°E), nous avons trouvé un affleurement
de la faille d’Utique (figure 75). Ce vallon représente une tranchée
paléosismologique naturelle au niveau de laquelle nous avons pu mettre en
évidence le jeu quaternaire de la faille d’Utique. En effet sur la rive droite de ce
vallon nous avons observé les grès du Pliocène décalés en faille inverse par
rapport aux silts rouges du Quaternaire moyen. Le plan de la faille est marqué par
des blocs du grès Pliocène imbriqués dans les silts du Quaternaire. Un morceau

128
de pièce de céramique a été trouvé emboité dans les silts au niveau du plan de la
faille (figure 76). L’analyse de cette céramique a été réalisée avec des
archéologues tunisiens (Mr Slim Khosref de l’Institut National du Patrimoine, et
Mme Boutheina Maraoui du Département d’Histoire de la Faculté des Sciences
Humaines et Sociales de Tunis). Ils ont indiqué qu’il s’agit d’un Tesson d’une
poterie tournée à cuisson peu poussée et qui correspond à un fragment d’une
panse d’amphore antique. Son âge doit se situer entre l’époque punique (10 ème-
9ème siècles au 2ème siècle av. J.-C. (soit en moyenne autour de 2500 ans avant nos
jours) et l’époque romaine (2ème-1er siècle jusqu’au V ème -VI ème siècle ap. J.-C.
(soit en moyenne autour de 2000 ans avant nos jours).
Le morceau de tesson de poterie (figure 77) imbriqué dans le plan de la faille
inverse d’Utique révèle que le dernier violent séisme attaché a cette faille a
secoué la région d’Utique il y a au plus 3000 ans. Il pourrait s’agir d’un violent
séisme dans la ville antique d’Utique mentionné par le catalogue de sismicité
historique de la Tunisie. Ce séisme a fait l’objet de plusieurs études depuis le
début du siècle dernier (Solignac 1927, Sieberg 1930, Voght 1992). Ces études
antérieures indiquent que ce séisme a provoqué la ruine de la ville d’Utique.
Nous avons trouvé aussi d’autres fragments de céramiques emboités dans les
couches silteuses attribués à l’Holocène. Ces données témoignent de
l’occupation de cette zone par les civilisations anciennes.

129
Figure 75: Coupe géologique qui montre le plan de la faille inverse d'Utique.

130
 Figure 76: La pièce de poterie imbriquée
dans la zone de broyage de la faille
d'Utique

Figure 77: Vue des deux facettes de la pièce de tesson trouvée dans la zone
de broyage de la faille inverse d'Utique

131
 Site 4 :
A l’Ouest de l’anticlinal, près de l’autoroute Tunis-Bizerte (37,031°N ; 9,998°E),
de nouvelles excavations de remblais ont mis à l’affleurement la zone affectée
par la faille d’Utique. En effet, nous avons observé dans une petite vallée haute
de 5 m environ que des silts rougeâtres sont affectés par deux branches de failles
entre lesquels sont coincés les grès pliocène. Le plan de la faille est NNW-SSE.
Au niveau du couloir d’une branche de faille, nous avons observé des argiles
verdâtres (figure 78).
Les silts rougeâtres et les argiles verdâtres sont attribués au Quaternaire.

Figure 78: Coupe qui montre la zone de passage de la faille

d'Utique

132
 Site 5 :
Une vallée d’excavation récente (37,0309°N ; 9,999°E) a mis au jour un contact
anormal entre le Pliocène et les silts rougeâtres du quaternaire. Les coordonnées
GPS de se site sont 37° 1'53,48"N et 9°59'55,08"E. Ce contact montre une
faille da direction N045 et qui montre deux branches de pendages différents. En
bas le pendage de la faille est de 30°NE, ensuite la partie sommitale montre un
pendage de 85°NE. Malgré l’allure franche de l’incision, le plan de la faille ne
montre pas de strie claire qui indique le sens du mouvement (figure 79).

Figure 79: Le plan de la faille d'Utique séparant les bancs de

grès Pliocène et les silts rougeâtres du Quaternaire (1 : silts
rougeâtre ; 2 : Pliocène)

133
 Site 6 :
A 1km à l’Ouest de la vile d’Utique (37,0356°N ; 10.028°E), nous avons fouillé
le trajet d’un vallon sur lequel est construit un petit barrage pour limiter l’effet de
l’érosion, pour rechercher des indices de tectonique récente à d’Utique. A une
altitude de presque 65m (coordonnées : 37,0337N ; 10,027E) nous avons observé
les niveaux gréseux du Pliocène surmontant une croute calcaire du pléistocène
moyen à supérieur et montrant un pendage opposé de 40° respectivement vers le
NE et le SW (figure 80). Il est difficile d’observer le plan du mouvement et les
stries associées, mais disposition géométrique indique un jeu inverse.

Figure 80: Contact anormal entre la croute calcaire

Quaternaire et les bancs gréseux du Pliocène

134
G : Description de la faille d’Utique :

Le plan de la faille que nous avons observé à plusieurs reprises en affleurement

dans plusieurs sites sur le flanc Sud de l’anticlinal d’Utique (figure 72), nous a
permis de tracer l’allure de la faille d’Utique. La direction des plans de faille au
niveau des affleurements est généralement en cohérence avec la direction de
l’axe du pli, à l’exception des affleurements retrouvés sur les limites Ouest et Est
de l’anticlinal où le plan de la faille tend vers la direction E-W. Le pendage la
faille varie entre 45° et 90°. Les stries sont généralement absentes à cause de la
nature friable des faciès sableux et silteux des formations Plio-quaternaires. Le
plan de la faille est souvent marqué par une zone de broyage dans laquelle sont
mêlés des fragments de grès pliocène aves les silts rouges à brunâtres et les
conglomérats du Quaternaire moyen à supérieur.
Le faciès silteux touché par le mouvement de la faille renferme des fragments de
poterie. L’âge de ces fragments revient à l’époque romaine selon les expertises
citées ci-dessus (figure 81). Nous avons trouvé ces poteries dans plusieurs sites et
à des altitudes plus ou moins importantes et qui atteignent parfois les 60 mètres
sur le piedmont Sud de l’anticlinal.

Figure 81: Les fragments de poterie de l'époque

Romain encaissés dans les silts brunâtres du
Pléistocène supérieur

135
Des affleurements du Villafranchien sous forme de conglomérats encaissant des
débris de grès Pliocènes et de faciès plus ancien dans une matrice calcaire, sont
rencontrés à des altitudes importantes (plus que 65 mètres) et toujours au dessus
des sites d’affleurement de la faille.
Le quaternaire dans cette région montre l’existence d’une croute calcaire
observable à des altitudes variables et qui ne dépassent pas les 70 mètres sur le
flanc Sud de l’anticlinal. Le pendage de la croute suit généralement l’allure
générale de l’anticlinal en considérant le caractère pédogénique des croutes
calcaires sur les reliefs. Mais par endroit, cette croute montre des pendages forts
relativement à la courbure du pli (figure 82). Nous expliquons ces pendages par
un mouvement tectonique lié à la faille d’Utique.

1 mètre

Figure 82: Le pendage élevé du quaternaire observé sur le flanc Sud de

l'anticlinal d'Utique

La compilation des observations de néotectonique et des données de géophysique

de subsurface permet de dresser l’allure de la faille d’Utique (figure 83, 84). En
général suivant la direction de l’anticlinal.
La morphologie de l‘anticlinal délivrée par le MNT, et par les profils
topographiques dressés à travers ce MNT, et la présence de la faille sur la

136
longueur du pli nous permettent de conclure que le pli d’Utique est un « pli de
propagation de rampe » (fault propagation fold).

Figure 83: Illustration schématique de la faille inverse d'Utique (le fond géologique
est celui de la carte géologique 1/50000 de Porto Farina)

137
138
Figure 84: Cartographie des facies Quaternaire sur le flanc Sud de l’anticlinal d’Utique
1 : Alluvions récents ; 2 : Silts rougeâtres (Pléistocène moyen) ; 3 : silts grisâtre riche en concrétion calcaire (Pléistocène
supérieur) ; 4 : croute calcaire feuilleté ; 5 : croute calcaire dur, 6 : niveau conglomératique (Villafranchien) ; 7 : grès
(Pliocène) ; 8 : Faille Quaternaire; 9 : site d’affleurement de la faille Quaternaire
H : Conclusion

L’anticlinal d’Utique et la faille qui borde son versant sud n’ont pas fait l‘objet d’étude
détaillée auparavant. Des observations de sections sismiques ont prouvé l’existence de
la faille E-W d’Utique. L’activité de cette faille a été suggérée par la rupture
topographique qui marque la plaine d’El Mabtouha qui limite au sud l’anticlinal
d’Utique et par le séisme destructeur d’Utique (412). Notre étude morphostructural,
paléosismologique et géophysique montre que l’anticlinal d’Utique est plis du type
« Fault-Propagation-Fold » et que la faille d’Utique est active. Les indices de
néotectoniques observés à plusieurs reprises sur le versant Sud de l’anticlinal
témoignent du jeu quaternaire de la faille d’Utique. Dans une tranchée naturelle où
affleure le plan de la faille, nous avons trouvé de morceaux de poteries, désignées par
les archéologues comme un « fragment d’une panse d’amphore antique » dont l’âge
est situé entre l’époque romaine et punique. Ces fragments de poterie sont imbriqués
dans le plan de la faille par l’effet du jeu de la faille lors d’un séisme. Ces données
semblent confirmer les données de la littérature concernant le violent séisme de 412
Après JC qui a secoué la ville antique d’Utique et a provoqué sa ruine totale. La mesure
du taux de raccourcissement sur la faille inverse d’Utique évaluée à 0,4 mm/an,
représente le dixième du taux moyen de déplacement vers le NW (N135) de l’Afrique
par rapport à l’Europe stable au niveau du nord de la Tunisie (4 mm/an).
Cette quantité de mouvement absorbée par la faille d’Utique est relativement
importante vue la complexité des structures géologique au Nord de la Tunisie et en
méditerranée au niveau du fossé de Sardaigne.

139
 CHAPITRE V

IMPLICATION SISMOTECTONIQUE
DE LA FAILLE D’UTIQUE SUR SES
ENVIRONS

140
A : Introduction: Sismotectonique et sismicité intraplaque

La sismotectonique se définit comme l’étude fondamentale des tremblements de terre,

leurs causes physiques et de leurs conséquences tectoniques, avec pour objectif
l’évaluation du risque sismique (prévention) et, a plus long terme, la prévision de
séismes par l’étude des événements antérieurs et des phénomènes précurseurs.
La Méditerranée est le siège de la convergence entre la plaque africaine et celle de
l’Eurasiatique. L’essentiel de l’activité sismique a lieu en Afrique du Nord, Italie, et les
îles helléniques. L’activité sismique de la Tunisie est la conséquence de la collision
(Meghraoui, 1988; Rebaï, 1992).
La sismicité tunisienne est typiquement intra-continentale, elle se produit de préférence
le long d’accidents tectoniques préexistants ou en évolution lente : chaine atlasique, axe
Nord-Sud, Accident Sud Atlasique,…
Etudier la sismicité intraplaque d’une région consiste donc à rechercher
indépendamment s’il existe des modifications locales du champ de contrainte et des
paramètres structuraux ou rhéologiques définissant des zones de faiblesse qui
favorisent des amplifications locales des contraintes. Ceci nécessite donc à la fois des
études sismologiques mais également géologiques et géophysiques (figure 85, 86).

Figure 85: Caractérisation des sources seismogénique et propagations

des ondes sismiques (i Philip et al., 2007)

141
Figure 86: Caractéristiques des différentes méthodes utilisées en sismotectonique (Philip et al., 2007)

142
B : Etudes sismotectonique en Tunisie : historique

Plusieurs auteurs ont compilé les données de sismicité avec les observations de
néotectonique pour élaborer un zonage sismotectonique de la Tunisie (Hatzfeld, 1978;
Hfaiedh, 1983; Gueddiche, 1992; Hfaiedh et al., 1992; Hfaiedh, 1997; Kacem, 2004).
La première carte de zonage sismotectonique, basé surtout sur la répartition de la
sismicité en Tunisie a permis à Hfaiedh (Hfaiedh, 1997) de proposer plusieurs zones
seismogénique capables d’engendrer des séismes de fortes magnitudes jusqu’à 6,5 à 7
degrés Richter. Pour ces auteurs ces différentes zones se résument comme suit :
- La vallée de Medjerda et le Nord Ouest du Pays : la magnitude maximale observée a
été de 5,7 et la magnitude maximale probable serait de 6,5.
- Le Nord Est et la région de Tunis : la magnitude maximale observée a été de 5,1 et la
magnitude maximale probable serait de 7,0.
- Le Sahel (Sousse, Monastir et Mahdia) : la magnitude maximale observée a été sur le
continent 4,6 et au large 5,1, la magnitude maximale probable serait de 6,5.
- La région de Gafsa : la magnitude maximale observée a été de 4,8, la magnitude
maximale probable serait de 6,5.
La magnitude maximale observée est celle enregistrée par les réseaux sismologiques,
ou bien celle indiquée par la littérature, ou bien obtenue par la conversion entre les
intensités maximales observées (Mercalli modifié) et la magnitude (formule de
Hfaiedh : M=0,6I+0,78 ; (Hfaiedh, 1983). La magnitude maximale probable a été
établie sur la base des investigations néotectoniques dans le cadre de l’élaboration de la
carte sismotectonique du pays par l’équipe de l’institut national de Météorologie
(Arsovsky et al., 1991). Dans le cadre de ce projet, Arsovski et son équipe ont fait
l’investigation des accidents néotectoniques dans le territoire du pays, et ensuite ils ont
établi un zonage sismotectonique en rattachant les épicentres aux failles recensées.
Pour l’élaboration de ce modèle sismotectonique une approche déterministe a été
adoptée:
- La quantification de la sismicité est basée essentiellement sur les intensités. La
base de données montre quelques lacunes, ce qui entrave l’élaboration d’un
modèle probabiliste cohérent, avec estimation de la magnitude maximale
probable et des périodes de retour.

143
 - La distribution spatiale des épicentres est très diffuse. En effet, certaines régions
possèdent une base de données importante alors que d’autres montrent des
hiatus assez exprimés.

C : Alea Sismique : Application au cas de la Tunisie

L’identification et la caractérisation des failles actives susceptibles d’être des sources

seismogénique, représente la partie essentielle de l’évaluation de l’aléa sismique, car il
permet d’évaluer la localisation (améliorer la qualité du catalogue des séismes),
l’intervalle de récurrence et l’ordre de magnitude du prochain séisme.
Les données géologiques, géophysiques, sismologiques disponibles sur une région
permettent d’effectuer des études probabilistes d’alea sismique (Youngs &
Coppersmith, 1985; Coppersmith & Youngs, 1986; McGuire, 1993). Cette méthode est
fondée sur la détermination des lieux d’occurrence des séismes, leur taille et leur
récurrence. Puis, en un site donné, l’aléa est calculé en sommant les contributions de
tous les scénarios possibles (toutes les combinaisons de magnitudes et de distances).
Donc, trois étapes sont nécessaires (figure 87):
 Identifier les zones sources dans la région d’étude afin de pouvoir dresser un
modèle sismotectonique
 Modéliser dans chaque zone source une courbe de récurrence à partir du
catalogue de sismicité (modèle de récurrence basé sur la formule de Richter.
 Choisir une loi d’atténuation du mouvement du sol adaptée à la région d’étude.
Cependant, des travaux récents (Cramer, 2001; Cramer et al., 2001; Newman et al.,
2001) soulignent les difficultés liées à cette méthode, principalement dans les régions
intraplaque. Elles concernent le manque de précision sur la récurrence des séismes
majeurs, et sur les modèles de faille utilisée (interaction des failles ou non). Enfin, cette
méthodologie est fondée sur la stationnarité de la sismicité et ne prend pas en compte
l’hypothèse d’un séisme intraplaque dans une zone non identifiée par un séisme passé
instrumental ou historique.

144
 Figure 87: Etapes des évaluations de l'aléa sismique par
l’approche probabilistique (D'après Reiter, 1990)

Le territoire Tunisien n’a pas connu de séisme majeur de magnitude supérieur à 6, et la

déformation aux limites de l’Afrique du Nord est très faible. Cependant, l’occurrence
de séismes historiques, dont le dernier, celui de Tunis en 1785, a une Intensité estimée
de VIII (donnée INM) et des séismes instrumentaux qui parfois dépassent la magnitude
5 rappellent l’existence du risque sismique non négligeable. Certain auteurs (e.g.
(Kacem, 2004)) ont utilisé la méthode de Cornell - Mc Guir (Cornell, 1968) pour
évaluer l’alea sismique probabiliste de la région du Nord-Est de la Tunisie.

145
 C-1 : Loi de récurrence de Gutenberg-Richter : log N = a – bm

La loi de Gutenberg-Richter permet d’établir un modèle de répartition du nombre de

séismes en fonction de la magnitude ou de l’intensité épicentrale, et donc la
caractérisation de l’activité des différentes zones sources. Cette loi donne la possibilité
d’étudier le comportement des sources sismique chacune à part, chose qui n’était pas
possible par le biais de la méthode probabiliste de la détermination de l’aléa sismique
(Schwartz & Coppersmith, 1984; Schwartz, 1988)

Les études de paléosismicité indiquent que des points individuels sur les failles et les
segments de failles tendent à se déplacer de la même distance approximativement à
chaque séisme. Ceci laisse suggérer que les failles individuelles génèrent à chaque fois
un séisme de même ordre de grandeur (à une demi-unité de magnitude près). Ces
séismes sont « les séismes caractéristiques de la source » (Philip et al., 2007).
La datation des séismes permet de donner un ordre de grandeur pour le temps de
récurrence. En effet, les tranchées paléosismologiques permettent une analyse fine de la
stratigraphie, et grâce à des datations par les méthodes isotopiques adéquates (14C,
U/Th, 36Cl, thermoluminescence,…), nous pouvons mieux contraindre les intervalles
de temps entre deux ruptures de surface.

Pour les magnitudes faibles, les données géologiques indiquent que la période de retour
du séisme caractéristique d’une source sismique donnée est généralement inférieure à
celle estimé par l’extrapolation de Gutenberg-Richter (Philip et al., 2007). La relation
proposée par Gutenberg et Richter, log N = a – bM, est la relation de base de la
séismologie statistique (Karnik, 1969). Les coefficients a et b ont fait l’objet de
plusieurs déterminations. Pour les séismes superficiels, Gutenberg et Richter ont trouvé
que la valeur « b » est de 0.9 ± 0.2.
Les valeurs « a » et « b » sont obtenus à l’aide des régressions appliquées sur les bases
de données de sismicité de la zone d’étude. Pour établir cette relation, nous avons
représenté le logarithme décimal du nombre cumulé de séismes en fonction de la
magnitude à partir de laquelle le catalogue est complet. Vu le nombre important des
séismes de faibles magnitudes (inférieur à 2) et qui sont parfois causé par des sources
de vibrations artificiels, nous avons choisi de prendre la magnitude 2 comme limite

146
inférieur pour homogénéiser le catalogue de sismicité instrumentale (figure 88)
(données sismologiques de l’Institut National de Météorologie). Cette courbe permet de
conclure que le catalogue de la sismicité pourrait être complet à partir du seuil de
magnitude 3. Cette valeur correspond aussi à la magnitude à partir de laquelle les
séismes superficiels sont habituellement ressentis. La magnitude maximale est de 6.
Cette courbe montre une pente avec un coefficient « b = 0,6543 ». L’erreur sur la
valeur b pour un coefficient de confiance de 95% est de 0.024. (Voir tableau ci-
dessous).

Coefficients Erreur-type

Constante 3,915839421 0,095310417

Variable b 0,65428829 0,024307081

On remarque que lorsque la valeur « b » augmente, le nombre des grands

séismes diminue par rapport au nombre des plus petits. L’application de la loi de
Gutenberg-Richter peut aussi s’appliquer sur les intensités épicentrales des séismes.
D’autre part, nous constatons que la magnitude maximale des séismes tunisiens
enregistrée est aux alentours de 6 (5,6 pour le séisme de 1970). Par conséquent,
l’examen de cette courbe de récurrence permet de déduire que l’énergie dégagée n’est
pas due essentiellement à des séismes très forts mais peut être attribuée à des séismes
dont la magnitude varie entre 3 et 6.

147
 log N = 3,915 – 0,654M

Figure 88: Loi de récurrence de Gutenberg-Richter pour les séismes de magnitude

supérieure ou égale à 2

Kacem (2004) a calculé l’équation de Gutenberg-Richter pour la zone du Nord-Est de

la Tunisie. Il a trouvé les valeurs suivantes : log N = 2,4192 – 0,748M. La différence
entre ces valeurs et celles que nous avons calculé est probablement du au nombre de
séismes dans la base de donnée utilisée pour le calcul (base de donnée complète
jusqu’à 2007 pour notre travail, et jusqu’à 2001 pour le travail de Kacem).

C-2 : Caractérisation de l’atténuation du mouvement du sol

L’évaluation du mouvement au sol attendu à une distance donnée de l’épicentre d'un

tremblement de terre d'une magnitude donnée, est le deuxième élément de l’estimation
de l’aléa sismique. Ces évaluations sont généralement sous forme d’équations appelées
« lois d'atténuation ». Ces lois d’atténuations expriment le mouvement du sol en
fonction de la Magnitude et la distance (et occasionnellement d'autres variables, telles
que les caractéristiques de la faille). Pour évaluer les mouvements du sol, généralement
on procède par l’évaluation de la vitesse maximale du sol (PGV), de l’accélération
maximale du sol (PGA), et plusieurs spectres d’accélérations (SA).
Les mouvements du sol de courte période (aux alentours de 0,2 seconde) peuvent avoir
un effet catastrophique sur les structures dont la période de vibration (période propre)

148
est courte, alors que les mouvements du sol de moyenne à longue période (1.0s, 2.0s,
etc.) auront les plus grands effets sur les structures de longue période de vibration.
Le mouvement du sol en un site donné, par exemple l’accélération maximale du sol,
dépend de la source sismique, de la propagation des ondes sismiques et du spectre de
réponse du site.

C-3 : Caractérisation du mouvement sismique :

La source sismique est caractérisée par sa magnitude, sa profondeur et son mécanisme

au foyer, la propagation des ondes sismiques dépend essentiellement de la distance par
rapport au site. La réponse du site dépend étroitement de la géologie locale
(classification des sites) ; ce qui est le sujet du microzonage.
Les lois d’atténuations du mouvement du sol peuvent être déterminées par deux
manières différentes : empiriquement, en se basant sur les enregistrements précédents
des mouvements du sol, ou théoriquement en utilisant des modèles sismologiques pour
produire des mouvements du sol synthétiques (simulation) qui caractérisent la source,
le site et les effets du parcours sur les ondes sismiques.
La forme de base de la fonction (logarithmique) de la loi d’atténuation du mouvement
du sol est définie comme suit (Reiter, 1990) :
ln Y = lnb1 + lnf1(M) + lnf2(R) + lnf3(M,R) + lnf4(P) + ln ε
Où :
- Y est un paramètre du mouvement du sol qui doit être estimé

- f1(M) est fonction de la magnitude du séisme

- f2(R) est fonction de la distance entre la source seismogénique est le site.

- f3(M,R) est une possible fonction de liaison entre la magnitude M et la distance

R.

- f4(Pi) sont des fonctions qui représentent les effets de site s’ils existent.

- ε est un terme d’erreur qui représente l’incertitude sur Y.

149
Ces fonctions sont spécifiques pour des types de sites (milieu rocheux, milieu élastique,
etc.). Les valeurs calculées de l’aléa sismique pour un milieu rocheux sont nettement
inférieures à celles calculées pour les sols.
Dans quelques pays de la zone Ibéro-Maghrébine, des lois d'atténuation était employée
pour obtenir des cartes sismotectoniques nationales.
Les paramètres sur lesquels sont basées les différentes lois d'atténuation présentent
également quelques différences. Certaines de ces lois donnent le PGA en fonction de la
magnitude M et de la distance R, (l'Algérie, le Maroc, la Tunisie). Alors que dans
d’autres cas les paramètres du mouvement du sol pris en considération sont l'intensité I,
en termes de magnitude et distance (Portugal), ou l’intensité épicentrale I0, et la
distance (Espagne).
Une carte d’accélération maximale a été calculée par (Jiménez et al., 1999). Afin de
garder un minimum d’homogénéité dans le calcul de cette carte, et en l’absence d’une
loi régionale basée sur assez de bonnes données d'accélération, seulement une loi
d'atténuation a été considérée. Cette carte montre que les régions de Tunis-Bizerte et
Ghardimaou-Bousalem ont des valeurs d’accélération entre 2.5 et 3.0 m/s2 (0.32 g)
pour une période de retour qui n’excède pas 50ans.
Dans le tableau n°3 nous montrons des lois d’atténuations utilisés dans les pays de la
région Ibéro-Maghrébine (in (Jiménez et al., 1999)) :

Pays Loi d’Atténuation

Tunisie PGA = (5600 e0.8M) / (R+40)2
Algérie PGA = (190.67 e0.823M) / (R+0.864 e0.463M) 1.561
Maroc Log PGA = -1.02 + 0.25M – 0.00255 (R2+7.32)1/2 – log (R2+7.32)1/2
Portugal I = 6.8 + 1.13M – 1.68 ln (R+14)
Espagne I = I0 + 12.55 – 3.53 ln (R+25)
Tableau 5: Lois d’atténuation calculées pour quelques pays de la Méditerranée
(Jiménez et al. 1999)

150
Application de la loi d’atténuation de Jiménez au Nord-est de la Tunisie :

La loi de Jimenez (Jiménez et al. 1999) permet d’évaluer l’accélération du sol en

fonction de la distance épicentrale par rapport à la source seismogénique étudiée. Elle
permet de donner une description qualitative des effets d’un éventuel séisme qui se
produira sur la faille seismogénique, sur ses alentours.
La loi de Jimenez est : PGA = (5600 e0.8M) / (R+40)2, avec :
P.G.A : Peak Ground Acceleration (accélération du sol)
M : représente la magnitude maximale probable évaluée pour la source seismogénique
étudiée. Dans notre cas, et pour la faille d’Utique M=6,8.
R est la distance depuis l’hypocentre, elle varie selon la profondeur du séisme. Ici nous
avons calculé l’accélération du sol pour des profondeurs variables du séisme (2, 5, 10
et 15km).
Nous avons calculé l’atténuation de l’énergie rayonnée par le déclenchement d’un
séisme sur la faille d’Utique (avec M=6) à des profondeurs variables (2km, 5km, 10km,
15km), (figure 89).

0,8

0,7

0,6

0,5 Z = 2km
PGA (g)

Z = 5km
0,4 Z = 10km
Z = 15km

0,3

0,2

0,1

0
1 10 100 1000
Distance (km)

Figure 89: Loi d'atténuation applicable pour le nord-est de la Tunisie calculé pour des profondeurs différents
de l’hypocentre (Z=2, Z=5, Z=10, Z=15) (d'après Jimenez, 1999)

151
Sur cette courbe, on remarque que pour une distance de 25 Km de la source sismique,
ce qui est le cas pour la ville de Tunis par rapport à la faille seismogénique d’Utique,
l’accélération du sol serait de 0,28g pour la profondeur 2km, 0.3g pour la profondeur
de 5km et 0.31g pour les profondeurs de 10 et 15km. Ces valeurs d’accélération du sol
ne tiennent pas compte des possibles amplifications des ondes sismiques dues aux
effets de site. D’autre part, ces valeurs assez élevées sont capables de provoquer des
dégâts importants, surtout lorsque les ondes sismiques sont amplifiés par les effets de
sites topographique et lithologique propres à la région Tunis.
Une carte d’iso-accélération a été établie par Kacem (2004). Elle montre un niveau de
PGA qui dépasse 0.28g aux alentours de la région de Tunis pour une période de retour
de 475 ans (figure 90).

152
 Figure 90: Cartes d'iso-accélération maximale du sol (Peak
Ground Acceleration) pour une période retour de 100ans (A) et
de 475 ans (B) (Kacem, 2004)

153
D : Implication d’un éventuel séisme de la faille d’Utique sur la zone du grand
Tunis :

D-1 : Effet de site

L’effet de site est le terme utilisé pour décrire la modification de la secousse sismique
induite par la géologie locale. Les effets de site se traduisent fréquemment par une
amplification de l’amplitude des enregistrements de la secousse sismique et de sa
durée. Ce phénomène peut être illustré en comparant les secousses sismiques
enregistrées sur un sol dur (rocher) et sur un sol mou (alluvions) à proximité.

Lorsque les propriétés du sol sont particulièrement défavorables (milieu sableux et

saturé en eau), le sol peut perdre sa cohésion (phénomène de liquéfaction). Les
bâtiments construits au-dessus de ce sol ne sont plus soutenus et peuvent basculer.

Dans des milieux géologiques à géométrie complexe (vallée glaciaire), les ondes
sismiques se réverbèrent sur les interfaces géologiques du fait de l’existence d’un fort
contraste entre les propriétés mécaniques des couches situées au centre de la vallée
(remplissage sédimentaire) et la roche encaissante. Les ondes sismiques sont piégées :
la durée et l’amplitude de la secousse sont fortement accrues. Les dégâts provoqués par
les tremblements de terre sont fréquemment la conséquence d’effets de site locaux
(concentration des dommages dans les vallées).

Tunis, capitale de la Tunisie, et sa région (1.900.000 habitants), s’étend sur un territoire

de 350.000ha (50 x 70 Km) ; la ville est de loin le centre le plus urbanisé du pays avec
une concentration dépassant 30,000 habitants/km2. Toutes les études sismologiques et
sismotectoniques précédemment cités, classent Tunis et ses environs comme étant une
zone exposée au risque sismique. Ces études confirment que Tunis se localise sur
l’intersection de deux zones les plus actifs du point de vue sismique : le domaine
Bizerte-Tunis au nord et le domaine Zaghouan-Grombalia au Sud. Jusqu’à présent dans
la majorité des études sismologiques, les séismes de référence, prisent en compte dans
l’évaluation de l’aléa sismique dans la région de Tunis, sont celui du 1 décembre 1970,
et d’autres séismes historiques (séisme d’Utique de 412 et ceux de Tunis de 954 et

154
956). Les valeurs d’accélération maximale du sol, calculées à partir de ce séisme,
varient entre 0,2 et 0,3g pour des périodes de retour de 200 à 500 ans (Kacem, 2004)
(figure 90). Ces résultats ont été estimés sur un affleurement rocheux (en sismologie on
parle de « rocher horizontal affleurant »). Lors de ce même séisme pourtant modéré, de
magnitude Ms = 5.1, les dégâts recensés aussi bien au niveau du bâti ancien que
moderne au niveau de la région de Sidi Thabet et ses environs, furent importants et
impliquent une amplification locale du mouvement sismique lors de sa traversée du
sous-sol tunisois.

Ces « effets de site », qui sont directement liés au contexte topographique et

lithologique propres à Tunis, méritent d’être évalués finement dans le but final d’un
micro zonage sismique de la capitale.
Tunis est installée sur des formations géologiques peu consolidées reposant sur une
assise plus rigide. Le piégeage des ondes sismiques dans les couches superficielles
donne naissance à un phénomène de résonance liée aux paramètres mécaniques,
géophysiques et géométrique du sous-sol. Cette résonance induit une amplification du
mouvement sismique à des fréquences particulières et parfois une prolongation
importante de la durée du signal.

D-1-1 : Contexte lithologique : effet de site lithologique

Le site de Tunis occupe une position surélevée entre les zones déprimées et occupées
par Sebkhet (ou lagune) Essijoumi à l’ouest et le lac de Tunis à l’est. Le lac de Tunis
communique avec la mer Méditerranée à l’est par un canal étroit, le goulet. La cuvette
ou vallée de Tunis s’étend entre ces reliefs et le lac de Tunis et constitue la principale
localisation de la capitale (centre ville). Le tissu urbain et résidentiel récent gagne les
abords immédiats du lac. Au sud on trouve des collines de direction subméridionale
(Sidi Belhassen, 88 m, et Essaîda Manoubia, 60 m). Les collines de Montfleury (60 m),
El kasbah (40 m) et Errabta (50 m) à l’ouest, et la colline du Belvédère (75 m) au Nord,
sont orientées NW–SE.

155
En effet, la lithologie de la région de Tunis est en étroite liaison avec les
caractéristiques morphologiques de la zone. Les plateaux subsidents sont généralement
occupés par des formations Quaternaires (figure 91):
- Les alluvions (sables, argiles et débris de roches), et les dépôts du Quaternaire marin
et les Sebkhas (silts, vases et évaporites).
- Les formations holocènes sont composées essentiellement de sables riche en débris de
coquilles et de silts.
Les bas de reliefs sont occupés par des formations Pléistocènes, constituées par :
- Argiles, sables et silts d’âge Pléistocène moyen à supérieur.
- Sable et argiles du Pléistocène inférieur
- Argiles et sable fin du Pléistocène
- Des éolianites d’âge Holocène
Les reliefs sont occupés par :
- Des formations Oligocène et Pliocène composées par des sables, argiles, grès et
des conglomérats.
- Des formations Jurassique à Eocène composées de dépôts dolomitique et
calcaire.
Les dépôts du Jurassique à Eocène sont considérés comme homogènes et compacts.
Les dépôts Oligocène et Pliocène sont considérés comme hétérogènes et compact et
sont favorables pour la construction. Les formations sableuses du Pléistocène et
l’Holocène sont favorables à la liquéfaction. Les formations triassique, les sols vaseux
et les alluvions récents sont les sites les moins favorables à la construction.

156
Figure 91: Carte des affleurements du grand Tunis et ses environs

157
 D-1-2 : Susceptibilité à la liquéfaction et profondeur de la ligne
piézométrique :

La liquéfaction est l’un des plus importants, intéressants, complexes phénomènes dans
le domaine de la sismologie. Ses effets dévastateurs ont été observés clairement par les
sismologues lors des grands séismes de l’Alaska (Mw=9.2) et de Niigata (Mw=7.5).
Le terme de liquéfaction, inventé par (Mogami & Kubo, 1953), a été utilisé en liaison
avec une variété de phénomènes qui provoquent la déformation du sol. Cette
déformation peut être causée par la monotonie, et l’effet cyclique de l’onde sismique
sur les sols saturés et non cohésifs, sous des conditions de non drainage. Ce phénomène
se produit par l’effet de l’augmentation de la pression de pore sous des conditions de
non drainage.
Les conditions géologiques des faciès susceptibles de provoquer la liquéfaction sont
restreintes. En effet les couches sableuses, non cohésives, à granulométrie uniforme, et
qui ne sont pas profondes, sont les plus capables de provoquer la liquéfaction. Par
conséquence, les dépôts fluviatiles, les colluvions et les dépôts éoliens saturés, sont
susceptibles de produire la liquéfaction. Ce phénomène a été aussi observé d’une façon
moins fréquente dans les lits des rivières, les plaines alluviales, les plages et les
terrasses. En général la susceptibilité de liquéfaction diminue avec l’âge de la couche
géologique. Les couches de l’holocène sont plus favorables que les couches
pléistocènes.
La liquéfaction se produit seulement dans les couches sableuses saturées, d’où
l’influence de la surface piézométrique des aquifères. Le risque de liquéfaction diminue
avec la profondeur de la surface piézométrique, et il est généralement observable
lorsque l’aquifère est à quelques mètres de profondeur.
Concernant la partie de la ville de Tunis qui d’étend depuis les berges du lac jusqu’à la
colline d’Errabta, les données lithologiques mentionnées précédemment, et la
corrélation des données de sondages (SRPT1, SRPT2 et MEH) selon une coupe N-S,
montrent bien l’existence de niveaux sableux, plus ou moins épais. Les données
hydrologiques montrent que dans la plupart du territoire de la ville de Tunis (au nord
du lac de Tunis) et de l’Ariana, le niveau piézométrique de l’aquifère libre ne dépasse
généralement pas les 10 mètres. Etant donné aussi les grandes dépressions vaseuses

158
(Sebkhet Essijoumi, Sebkhet Ariana), on peut conclure que cette région de la ville de
Tunis et de l’Ariana est très susceptible pour le risque de liquéfaction.
Pour évaluer le risque de liquéfaction dans une région donnée, le sismologue doit
répondre à ses questions :
- Est-ce que le type du sol est susceptible à la liquéfaction ?

- Si le sol est susceptible, est ce que le phénomène de liquéfaction peut être

déclenché ?

- Si le phénomène de liquéfaction est déclenché, est ce qu’il peut causer des

dommages ?

Si la réponse à la première question est négative, l’évaluation du risque de liquéfaction

peut être arrêtée et conclure que le risque de liquéfaction n’existe pas.
La plupart des informations sur le phénomène de liquéfaction provient des
investigations de terrain après les séismes. Ces investigations ont montré que la
liquéfaction se reproduit souvent sur les mêmes lieux lorsque les conditions des
couches géologiques et de niveau piézométrique sont constantes (Youd, 1984). Ainsi
les cas historiques de liquéfaction peuvent nous aider pour identifier des sites
spécifiques, ou les conditions de site, qui peuvent être susceptibles à la liquéfaction lors
des futurs séismes.
Les investigations de terrain après les séismes ont montré aussi que le risque de
liquéfaction dans une région donnée, est particulièrement lié à la distance qui la sépare
de la source sismique.

La position de la ligne piézométrique est un facteur très important pour l’analyse de la

liquéfaction dans une région donnée. El May (El May et al., 2009), par le biais des
puits, forages et piézomètres présents dans la région, à dressé une carte d’iso-
profondeur de la ligne piézométrique dans la région de l’Ariana. Elle a identifié deux
aquifères superposée. Ces deux aquifères seraient impliquées dans tous les séismes
possibles. En effet, le niveau statique est généralement superficiel, et il ne dépasse pas
les 10 mètres. Il est parfois inférieur à 2 mètres dans la région de la Soukra. Il est aux

159
alentours de 3 m au niveau de la plaine de l’Ariana et de 5m au bas du relief de Djebel
Nahli (figure 92).
Cet état de la ligne piézométrique, ainsi que les caractéristiques du remplissage
sédimentaire de la région de l’Ariana sont des meilleures conditions pour le risque de la
liquéfaction lors d’un séisme possible. Un cas de liquéfaction a été observé dans la
localité de sidi Thabet à 10 Km de la Ville de Tunis lors du séisme de 12/12/1970de
magnitude 5.6. Pour des magnitudes de séismes plus grandes (aux alentours de 6,5) le
risque de liquéfaction est plus élevé dans cette région.

Figure 92: Carte piézométrique de la nappe d'eau au niveau de la ville de

Tunis (El May et al., 2009)

Ambraseys (Ambraseys, 1988) a compilé des données mondiales des séismes pour
évaluer une distance épicentrale limite au-delà de laquelle la liquéfaction n'a pas été
observée pour des séismes de diverses magnitudes. La distance de la source sismique

160
au site où la liquéfaction est probable, augmente rapidement avec la magnitude. Cette
relation est représentée par la courbe sur la figure 93.
Si on admet que le sous sol de la majorité du territoire de la ville de Tunis est favorable
au risque de liquéfaction, la source seismogénique d’Utique qui pourra déclencher un
séisme de magnitude pouvant atteindre 6, est capable de produire le phénomène de
liquéfaction jusqu’à une distance d’à peu près 60 Km de l’anticlinal d’Utique.

Figure 93: Relation entre la distance épicentrale limite des sites où la liquéfaction a été observé, et la
magnitude des séismes superficiels (d'après (Ambraseys, 1988) )

D-1-3 : Effets de site dus aux caractéristiques mécaniques des

formations superficielles et de la géométrie des interfaces entre les
formations géologiques:

Les propriétés géo-mécaniques des remplissages sédimentaires (bassins et vallées) et

des substratums rocheux sont généralement très contrastées (fort contraste
d’impédance). Les couches superficielles se comportent aux mouvements sismiques

161
comme des oscillateurs dont la réponse dépend essentiellement des caractères du signal
sismique incident, de la géométrie des couches sédimentaires en profondeur et de leurs
propriétés mécaniques intrinsèques. Vis-à-vis du signal sismique incident, les
formations superficielles jouent le rôle de filtre, notamment avec un écrêtage dans les
hautes fréquences (périodes autour de 0,2-0,5 s), et par contre avec une amplification
dans les fréquences plus basses.
D’autre part les formations superficielles sont à l’origine de nombreux phénomènes de
réflexion et de réfraction et peuvent servir de guide d’ondes avec la création fréquente
d’ondes de surface. Ces phénomènes donnent lieu à des mouvements sismiques en
surface dont l’amplitude et la durée sont nettement supérieures à celles enregistrées à
faible distance sur un substratum rocheux. Au niveau d’un remplissage alluvionnaire,
les phénomènes de résonance dominent largement dans les modifications que subit le
signal sismique incident, en particulier dans les fréquences qui correspondent à la
période fondamentale de la couche d’alluvions.
Si on considère le cas le plus simple d’un site (sans faire intervenir ici un rapport de
forme du bassin ou de la vallée) correspondant à un milieu rocheux caractérisé par une
masse volumique ρ1 et une vitesse des ondes transversales V1 recouvert par une
formation superficielle monocouche d’épaisseur h, de masse volumique ρ 2 et avec une
vitesse des ondes transversales V2, l’amplification à la résonance sera donnée par:
ρ1 * V1/ ρ2*V2
Et la fréquence de résonance fondamentale f0 sera donnée par la relation :
f0 = V2/4h
Ainsi, pour une couche d’alluvions d’épaisseur h=100m, avec ρ2=1.7 g/cm3 et
V2=600m/s reposant sur un substrat caractérisé par ρ1=2.5 g/cm3 et V1=2100m/s,
l’amplification à la résonance est de 5.15 pour la fréquence fondamentale de 1.5 Hz.
On remarque que l’on se trouve dans une gamme de fréquence où les ondes sont les
plus destructrices.
Cette évaluation ne tient pas compte de la géométrie du bassin sédimentaire ou de la
vallée, aussi il sera nécessaire de réaliser une modélisation 2d qui prend en compte
notamment le rapport de forme h/L, où h est la profondeur et L demi-largeur du bassin.
Au centre du bassin, où l’amplification maximale correspond à la fréquence de
résonance fondamentale, celle-ci deviens :

162
 F0 = V2/4h * √1+h2/L2
L’effet du rapport h/L sera d’autant plus important si le bassin est étroit et profond.
Ainsi, à partir de l’aléa déterminé sur les roches du substratum (au rocher) obtenues à
partir des données sismologiques fournies à l’échelle régionale pour chacune des
sources sismiques, on peut, au moyen de simulations numériques 2D, calculer
l’amplification le long de profils. On propose des spectres de réponses pour chacune
des zones et des différents types de sol (au rocher, bord ou centre de bassins ou de
vallée), qui seront pris en compte pour le dimensionnement du bâti dans les projets de
constructions parasismiques (in (Philip et al., 2007)).

D-2 : Récurrence de la faille active d’Utique

L’étude pluridisciplinaire de la faille d’Utique (géologie de surface, géophysique,…),

ainsi que les données archéologiques ont permis d’affirmer que son activité la plus
récente date au plus de l’époque romaine (≈2000 ans B.P). D’autre part, en se basant
sur les caractéristiques géologiques et géophysique de la faille, et en utilisant des
équations empiriques nous avons pu estimer la magnitude du séisme d’Utique à une
magnitude de 6 sur l’échelle ouverte de Richter.
Des auteurs (Slemmons et Polo, 1986) ont établi des abaques qui permettent de
classifier les failles d’après leur taux d’activité et par la suite déterminer le temps de
récurrence pour une magnitude donnée (figure 94). En appliquant ces abaques sur la
faille d’Utique, le point de coordonnée 6 sur l’axe des abscisses et 0.4 sur la courbe des
vitesses, donne une période de retour de 1300 ans.
La datation archéologique de la pièce de céramique que nous l’avons découverte inter-
mêlée dans le plan de la faille d’Utique, réalisée par des archéologues de la faculté des
sciences humaines et sociales de Tunis (Pr. Boutheina MRAOUI) et de l’Institut
National du Patrimoine (Mr. Slim Khosref), à permis de confirmer les données de
littérature qui disent que le séisme d’Utique est daté entre l’époque Romaine et
Punique ; c'est-à-dire il y à, à peu près 2000 ans.
La période de retour du séisme d’Utique évalué ci-haut à presque 1300 ans, et l’âge du
séisme déterminé essentiellement par le biais des données de l’archéologie, affirment
que le risque du rejeu de la faille d’Utique est très grand de nos jours. Donc des

163
précautions spéciales et des règles de construction parasismique doivent être prises en
comptes lors la réalisation de tous les ouvrages d’art aux alentours de la source
seismogénique d’Utique.

Figure 94: Effet du "fault-Slip-rate" et du magnitude du séisme sur

la période du retour (d'après (Slemmons & De Polo, 1986))

164
Conclusion et Perspectives

Le recensement des différentes structures actives citées par la littérature dans la région
du Nord-Est de la Tunisie a mis en évidence les arguments morphologiques et les
indices de déformation qui attestent d’une activité tectonique très récente. Ces indices
démontrent que la phase tectonique post-villafranchienne est la phase compressive
majeure au niveau du Nord-Est, et elle est responsable de la mise en place des reliefs
récents (exemple : pli d’Utique).
Les mécanismes au foyer réalisés sur les séismes les plus importants enregistrés au
Nord de la Tunisie, plaident en faveur d’une compression actuelle orientée NNW-SSE.
Les séismes enregistrés dans le territoire tunisien sont superficiels et ne dépassent
guère la profondeur de 15 km. Les ruptures de surface qui peuvent nous renseigner
directement sur la cinématique des accidents actifs sont rarement observées suite aux
séismes tunisiens.
Des indices certains d’activité tectonique récente ont été découvert au niveau des failles
d’Utique, de Messeftine de Djebel Nahli et de Djebel Ammar. Ces indices permettent
de classer ces failles comme des sources seismogénique à potentiel important.
L’étude pluridisciplinaire de la faille d’Utique (géologie de surface, géophysique,…),
ainsi que les données archéologiques ont permis d’affirmer que son activité la plus
récente date au plus de l’époque romaine (≈2000 ans B.P). D’autre part, en se basant
sur les caractéristiques géologiques et géophysique de la faille, et en utilisant des
équations empiriques nous avons pu estimer la magnitude du séisme d’Utique à une
magnitude de 6 sur l’échelle ouverte de Richter. L’équilibrage de la coupe de la
structure d’Utique à permis d’estimer la quantité de raccourcissement post-Miocène
supérieur accomodée par la structure d’Utique à 275 mètres. En supposant que ce
raccourcissement a eu lieu depuis le villafranchien (1,8 Ma), on propose une vitesse de
raccourcissement d’environ ~0,15 mm/an.
Les caractéristiques topographique et lithologique de la ville de Tunis sont très
favorables aux amplifications de l’onde sismique causée par les effets de sites.
Plusieurs régions (telques les alentours du Lac de Tunis) sont susceptibles au
phénomène de liquéfaction.

165
Nous avons évalué la période de retour du séisme d’Utique évalué à presque 1300 ans,
et le dernier mouvement important de la faille active déduit essentiellement par le biais
des données de l’archéologie, affirment que le risque du rejeu de la faille d’Utique est
grand de nos jours. Donc des précautions spéciales et des règles de construction
parasismique doivent être prises en comptes lors la réalisation de tous les ouvrages
d’art aux alentour de la source seismogénique d’Utique.

L’étude en détail de la structure d’Utique pourra être appliquée sur d’autres structures
de la région Nord orientale du pays dans le but de préciser la quantité de mouvement
accommodée par ces structures et leurs potentiels seismogénique sur les importantes
agglomérations socio-économiques dans cette région.
Une étude plus accentuée sur l’aléa sismique local au niveau de la ville de Tunis pourra
aboutir à l’élaboration des cartes de microzonage sismiques. Ces cartes seront de
grande utilité pour les services du ministère de l’équipement pour l’implantation de
nouveau projet. En effet, pour mieux détailler l’effet de site lithologique, des mesures
du rapport H/V doit être fait dans les différentes parties de la cuvette de Tunis. Ainsi
qu’une carte géotechnique qui pourra nous renseigner sur les caractéristiques
mécaniques des formations géologiques de subsurface.

166
Liste Bibliographique

Ahmadi, R., 2006, Utilisation des marqueurs morphologiques, sédimentologiques et

icrostructuraux pour la validation des modèles cinématiques de plissement. Application à
l’Atlas méridional tunisien [PhD thesis], Univ. Nantes.
Ahmadi, R., Ouali, J., Mercier, E., Mansy, J.L., Brigitte Van-Vliet Lanoe, B., Launeau, P., Rhekhiss, F.,
& Rafini, S., 2006, The geomorphologic responses to hinge migration in the fault-related
folds in the Southern Tunisian Atlas: Journal of Structural Geology, v. 28, p. 721-728.
Ambraseys, N., 1988, Engineering seismology: J. Earthq. Eng. and Struct. Dyn., v. 17, p. 66.
Ammar, H., 1986, Hupothèse pour un modèle d'évolution dynamique des bassins de la basse
Vallée de la Medjerda, Univ. de Tunis.
Arambourg, C., & Arnauld, M., 1950, Note sur les fouilles paléontologiques exécutées en 1947-
1948 et 1949 dans le gisement villafranchiende la Garaat Ichkeul: Notes Serv. Géol.
Tunisie, v. 2, p. 9.
Arambourg, C., Arnould, M., & 1949, Note sur les fouilles paléontologiques exécutées en 1947-48
et 1949 dans le gisement villafranchien de la Garaet Iohkeul Bull. Soc. Sci. Nat. de Tunisie,
v. II.
Argus, D.F., Gordon, R.G., DeMets, C., & Stein, S., 1989, Closure of the Africa-Euroasia-North
America plate motion circuit and tectonics of the Gloria Fault: J. Geophys. Res., v. 94, p.
5585-5602.
Arsovsky, M., Jancervki, N., Domurdzanov, N., Hatzdzievki, D., Hfaiedh, M., & Attafi, K., 1991,
Investigations for elaboration of seismotectonic map and draft seismic design code of
Tunisia. Vol2: Seismic hazard analysis and seismic zoning map of Tunisia., Report IZIIS 500-
91-60.
Bard, E., Hamelin, B., & Fairbanks, R.G., 1990, U-Th ages obtained by mass spectrometry in corals
from Barbados. Sea level during the past 130,000 Years. : Nature, v. 346, p. 456-458.
Ben Ayed, N., 1980, Le rôle des décrochements E-W dans l’évolution structurale de l’Atlas
tunisien: C.R. Somm. Soc. Géol. France, v. série D, p. pp. 29-32.
—, 1986, Evolution tectonique de l’avant - pays alpine de Tunisie. Stratigraphie, caractérestiques
géophysiques et évolution géodynamique: Tunis, Tunis El Manar.
—, 1993, Evolution de l'avant-pays de la chaine alpine de la Tunisie du début du Mésozoïque a
l'actuel: Ann. Mines et Géol., v. 32, p. 286.
Ben Ayed, N., Bobier, C., & Viguier, C., 1978, Sur la tectonique récente de la plage du R'Mel, à l'Est
de Bizerte (Tunisie Nord-Orientale): Géol. Medit., v. VI, p. 423-426.
—, 1983, Interprétation tectonique unitaire des Atlas maghrébins: Rapp. C.I.E.S.M, v. 28, p. 99-
101.
Ben Ayed, N., & Viguier, C., 1981, Interprétation de la Tunisie Atlasique: C. R. Acad. Sci. Paris, v.
292 [Link], p. 1445-1448.
Billi, A., Faccenna, C., Bellier, O., Minelli, L., Neri, G., Piromallo, C., Presti, D., Scrocca, D., &
Serpelloni, E., 2011, Recent tectonic reorganization of the Nubia-Eurasia convergent
boundary heading for the closure of the western Mediterranean: Bull. Soc. géol. France, v.
4, p. 279-303.
Bouaziz, S., Jedoui, Y., Barrier, E., & Angelier, J., 2003, Néotectonique affectant les dépôts marins
tyrrhéniens du littoral sud-est tunisien : implications pour les variations du niveau
marinNeotectonics in the Tyrrhenian marine deposits of the southeastern Tunisian coast:
implications for sea level changes: C. R. Acad. Sci. Paris, v. 335, p. 247-254.
Boudiaf, A., Philip, H., Coutelle, A., & Ritz, J.F., 1999, Découverte d’un chevauchement d’âge
Quaternaire au sud de la Grande Kabylie (Algérie): Geodin. Acta, v. 12, p. 71-80.

167
Boukadi, N., 1994, Structuration de l’Atlas de Tunisie; signification géométrique et cinématique
des noeuds et des zones d’interférences structurales au contact de grands couloirs
tectoniques, Univ. Tunis II.
Bourgou, M., 1982, Contribution à l’étude géomorphologique des accumulations dunaires de la
peninsule du Cap Bon (Tunisie), Univ. Tunis.
Bousquet, J.C., & Philip, H., 1986, Neotectonic of the Calabrian arc and Apennines (Italy): an
exemple of plio-quaternary evolution from island arcs to collisional stages. the origin of
arcs by Wezel, in F.C, ed., p. 305-326.
Boutib, L., Zargouni, F., & Melki, F., 1997, Tectonique récente et sismotectonique de la région du
grand Tunis : apport télé-analytique et données de terrains: Africa Geoscience Review, v.
Vol. 4, p. pp. 395-403.
Brun, A., 1983, Etude palynologique des sédiments marins holocènes de 5000 B.P à l'Actuel dans
le golfe de Gabès (mer Pélagienne): Pollen et Spores, v. XXV, p. 437-460.
—, 1985, La couverture steppique en Tunisie au Quaternaire supérieur: C. R. Acad. Sci. Paris, v.
301, p. 1085-1090.
Brun, A., & Rouvillois-Brigol, M., 1985, Apport de la Palinologie à l'Histoire du peuplement en
Tunisie, Palynologie archéologie: C.N.R.S, Paris, p. 231-226.
Burrolet, P.F., 1951, Etude géologique des basins mio-pliocène du nord-est de la Tunisie ( region
entre mateur et Ghar El Mileh): Ann. Mines et géol., v. 1ère série.
—, 1956, Contribution à l’étude stratigraphique de la Tunisie centrale: Ann. Min. Géol., v. N°18, p.
352.
Butler, R.W.H., Lickorish, W.H., Grasso, M., Pedley, H.M., & Ramberti, L., 1995, Tectonics and
sequenc stratigraphy in Messinian Basins, Sicily: constraints on the initiation and
termination of the Meditrranean salinity crisis: Geol. Soc. Amer. Bull., v. 107, p. 425-439.
Carminati, E., Wortel, M.J.R., Meijer, P.T., & Sabadini, R., 1998, The two-stage opening of the
western-central Mediterranean basins: a forward modeling test to a new evolutionary
model: Earth and Planetary Science Letters, v. 160, p. 667-679.
Castany, G., 1951, Etude géologique de l’Atlas tunisien oriental. Thèse Doctorat ès Sciences:
Annales Mines et Géologie, Tunis, v. Vol 4, p. 632p.
—, 1953, Les plissements Quaternaires en Tunisie: C. R. Somm. Soc. Géol. France, p. 198-200.
Chelbi, F., Paskoff, R., & Trousset, P., 1995, La baie d'Utique et son évolution depuis l'antiquité.
Une réévaluation géoarchéologique: Antiquités africaines, v. 31, p. 7-51.
Chen, J.H., Curran, H.A., White, B., & Wasserburg, G.J., 1991, Precise chronology of the last
interglacial period: 234U-230Th data from fossil coral reefs in the Bahamas: Geol. Soc.
Am. Bull., v. 103, p. 82-97.
Chihi, L., 1984, Etude tectonique et microtectonique du graben de Kasserine (Tunisie Centrale) et
des structutres voisines J. Selloum et L. Maargba, Univ. Paris Sud - Centre d'Orsay.
—, 1995, Les fossés néogènes à quaternaires de la Tunisie et de la mer pélagienne: leur étude
structurale et leur signification dans le cadre géodynamique de la Méditerranée centrale,
Fac. Sc. de Tunis, p. 324.
Colleuil, B., 1976, Etude stratigraphique et néotectonique des formations néogènes et
quaternaires de la région de Nabeul - Hammamet (Cap Bon, Tunisie). Univ. Nice.
Coppersmith, K.J., & Youngs, R.R., 1986, Capturing uncertainty in probabilistic seismic hazard
assessments within intraplate tectonic environments, Proc. Third U.S. Natl. Conf. Earthq.
Eng., Volume 1, p. 301-312.
Cornell, C.A., 1968, Engineering seismic risk analysis: Bull. Soc. America, v. 58, p. 1503-1606.
Cramer, C.H., 2001, A seismic hazard uncertainty analysis for the New Madrid seismic zone:
Engineering Geology, v. 62, p. 251-266.
Cramer, C.H., Wheeler, R.L., Frankel, A.D., Talwani, P., & Lee, R.C., 2001, Incorporating uncertainty
into probabilistic seismic hazard maps for the Central and Eastern US: Seis. Res. Lett., v.
72, p. 127 abstract.

168
Crampon, N., 1971, Etude géologique de la bordure des Mogods du pays de Bizerte et du Nord de
hédil ( Tunisie septentrionale ), Nancy, 703p.
D'Agostino, N., Avallone, A., Cheloni, D., D'Anastasio, E., Mantenuto, S., & Selvaggi, G., 2008,
Active tectonics of the Adriatic region from GPS and earthquake slip vectors: Journal of
Geophysical Research, v. 113, p. 19.
D’Agostino, N., & Selvaggi, G., 2004, Crustal motion along the Eurasia-Nubia plate boundary in the
Calabrian Arc and Sicily and active extension in the Messina Straits from GPS
measurements: J. of Geoph. Res, v. 109.
DeMets, C., Gordon, R.G., Argus, D.F., & Stein, S., 1990, Current plate motions: Geophys. J. Int., v.
101, p. 425-478.
Devolve, J.J., Martin, G., & Menant, G., 1981, Le Djebel Ammar de l'Ariana, esquisse structurale
issue de données stratigraphiques nouvelles: Note Serv. Géol. Tunisie, v. 46, Travaux de
Géologie Tunisienne, p. 5-10.
Devoti, R., Ferraro, C., Gurguen, E., Lanotte, R., Luceri, V., Nardi, A., Pacione, R., Rutigliano, P.,
Sciarretta, C., & Vespe, F., 2001, Geodetic control on recent tectonic movements in the
central Mediterranean area: Tectonophysics, v. 346, p. 151-167.
Dlala, M., 1984, Tectonique récente et microtectonique de la région de Kasserine (Tunisie
centrale): Montpellier, Univ. Sc; et Tech. de Languedoc.
—, 1992, Sismotectonique study in Northen Tunisia: Tectonophysics, v. V.209, p. pp. 171-174.
—, 1994, Evidence for interchange of the stress ( σ1 ; σ2) and its geodynamic implication: Rock
View, Géos. Af., v. Vol. 1, p. PP. 533-545.
—, 1995, Evolution géodynamique et tectonique superposées en Tunisie implications sur la
tectonique récente et la sismicité: Tunis, Univ. Tunis II.
—, 2002, Les manifestations tectono-sédimentaires d’âge Campanien-Maastrichtien en Tunisie :
implications sur l’évolution géodynamique de la marge Nord-Africaine: C. R. Geoscience,
v. n°334, p. p. 135-140.
Dlala, M., & Ben Ayed, N., 1988, Les défrormations du Quaternaire récent du Graben de Foussana:
Geol. Med., v. 3, p. 171-176.
Dlala, M., & Hfaiedh, M., 1993, Le séisme de Metlaoui: une tectonique active en compression: C.
R. Acad. Sci. Paris, v. 317, p. 171-174.
Dlala, M., & Kacem, J., 2008, Alea sismique régional du grand Tunis et ses environs (Tunisie Nord
Orientale): Ann. Equipement, p. 13-27.
Dlala, M., & Rebai, S., 1994, Relation compression -extension miocène supérieur à Quaternaire -
implication sismotectonique: C.R. Acad. Sc. , v. T 319, p. pp. 945-950.
Doglioni, C., Fernandez, M., Gueguen, E., & Sabat, F., 1999, On the interference between the early
Apennines-Maghrebides back arc extension and the Alps-betics orogen in the Neogene
geodynamics of the Western Mediterranean: Bol. Soc. Geol. Ital., v. 118, p. 75- 89.
Domzig, A., Yelles, K., Le Roy, C., Déverchère, J., Bouillin, J.-P., Bracène, R., Mercier de Lépinay, B.,
Le Roy, P., Calais, E., Kherroubi, A., Gaullier, V., Savoye, B., & Pauc, H., 2006, Searching for
the Africa–Eurasia Miocene boundary offshore western Algeria (MARADJA’03 cruise): C. R.
Acad. Sci. Paris, v. 338, p. 80-91.
Dubordieu, G., 1956, Etude géologique de la région de l'Ouenza (confins algéro-tunisien), Paris.
Durand-Delga, M., & Fonboté, J.M., 1980, Le cadre structural de la Méditerranée occidentale, in
BRGM, m., ed., 26ème Congrès géol. inter. les chaînes alpines issues de la Téthys: Paris.
El Euch, N., 2007, Sédimentologie et startigraphie séquentielle du Miocène moyen à supérieur de
la Tunisie Nord Orientale (Cap Bon, Bassin de Kechabta et Golf de Tunis), Tunis II.
El May, M., Kacem, J., & Dlala, M., 2009, Liquefaction susceptibility mapping using geotechnical
laboratory tests: Int. J. Environ. Sci. Tech., v. 6, p. 299-308.
Elmejdoub, N., & Jedoui, Y., 2009, Pleistocene raised marine deposits of the Cap Bon peninsula (N-
E Tunisia) : Records of sea-level highstands, climatic changes and coastal uplift:
Geomorphology, v. 112, p. 179-189.

169
Facenna, C., Mattei, M., Funiciello, R., & Jolivet, L., 1997, Styles of back-arc extension in the
Central Mediterranean: Terra Nova, v. 9, p. 126-130

Fournet, A., 1971, Les gisements à faune Villafranchienne de Tunisie: Notes Ser. Géol. Tunisie, v.
34, p. 53-69.
Frizon De Lamotte, D., Saint Bezar, B., Bracene, R., & Mercier, E., 2000, The two main steps of the
Atlas building and geodynamics of the western Mediterranean: Tectonics, v. 19, p. 740-
761.
Ghisetti, F., 1979, Relazioni tra strutture e fasi trascorrenti e distensive lungo i sistemi Messina-
Fiumefreddo, Tindari-Letojanni e Alia-Malvagna (Sicilia nord-orientale) : Uno studio micro
tettonico: Geol. Romana, v. 18, p. 23-58.
Ghisetti, F., Gorman, A.R., Grasso, M., & Vezzani, L., 2009, Imprint of foreland structure on the
deformation of a thrust sheet : the Plio-Pleistocene Gela nappe (southern Sicily, Italy):
Tectonics.
Girardin, N., Hatzfeld, D., & Guiraud, R., 1977, La sismicité du Nord de l'Algérie C.R. Somm. Soc.
géol. Fr., v. II, p. 95-100.
Grimaud, S., Boillot, G., Collette, B.J., Mauffret, A., Miles, P.R., & Roberts, D.B., 1982, Western
extension of the Iberian-European plate boundary during the Early Cenozoic (Pyrenean)
convergence: A new model: Marine Geology, v. 45, p. 63-77.
Gueddiche, M., 1992, Analyse de la sismicité en Tunisie et étude sismique. sismotectonique et
néotectonique de la Tunisie septentrionale: Orsay, Univ. Paris Sud.
Hatzfeld, D., 1978, Etude sismotectonique de la zone de collision ibéro-maghrébine, Université
Scientifique et Médicale de Grenoble.
Hatzfeld, D., & Bensari, D., 1977, Grands profils sismiques dans la région de l'arc de Gibraltar: Bull.
Sco. géol. France, v. 7, p. 749-758.
Hfaiedh, M., 1983, Etude sismotectonique de la Tunisie Nord orientale: Centre d’Orsay, Univ. Paris
Sud.
—, 1997, Alea sismique: présentation et évaluation, Deuxième séminaire de génie sismique:
Hammamet, 17-19 Octobre, Institut National de Météorologie, p. 42-57.
—, 1998, Alea sismique: presentation et avaluation, 2ème seminaire national de génie sismique:
Hammamet, p. 42-57.
Hfaiedh, M., Attafi, K., Arsovski, M., Jancevski, J., Domurdzanov, N., & Turki, M.M., 1992, Carte
sismotectonique de la Tunisie, échelle 1/500000, Programme de réduction Risque
sismique en Tunisie.
Hfaiedh, M., Ben Ayed, N., & Dorel, J., 1985, Etude néotectonique et sédimo-tectonique de la
Tunisie Nord-Orientale: Note Serv. Géol. Tunisie, v. 51.
Hfaiedh, M., & Chadi, M., 1985, Seismicity of Tunisia and neighbouring area., European
Geotraverse. European Foundation: Strasbourg, p. 261-267.
Hollenstein, C., Kahle, H.G., Geiger, A., Jenny, S., Goes, S., & Giardini, D., 2003, New GPS
constraints on the Africa-Eurasia plate boundary zone in southern Italy: Geophys. Res.
Lett., v. 30.
Jauzein, A., 1967, Contribution à l'étude géologique de la Tunisie septentrionale: les confins de la
dorsale tunisienne: Ann. Mines et Géol., v. 22, p. 1-475.
Jedoui, Y., Davaud, E., Ben Ismaïl, H., & Reyss, J.-L., 2002, Analyse sédimentologique des dépôts
marins pléistocènes du Sud-Est tunisien : mise en évidence de deux périodes de haut
niveau marin pendant le sous-stade isotopique marin 5e (Eémien, Tyrrhénien): Bull. Soc.
géol. France, v. 173, p. 63-72.
Jedoui, Y., Reyess, J.-L., Kallel, N., Montacer, M., Ben Ismail, H., & Davaud, E., 2003, U-series
evidence for two high Last Interglacial sea levels in southeastern Tunisia Quaternary
Science Review, v. 22, p. 343-351.

170
Jiménez, M.J., García - Fernández, M., & the GSHAP Ibero - Maghreb Working Group (Chadi,
M.E.F., D.; Izquierdo, A.; Martinez-Solares, J.M.; Sousa-Oliveira, C.; Tadili, B.A.), 1999,
Seismic hazard assessment in the Ibero-Maghreb Region: Annali Geofis, v. 42, p. 1057-
1066.
Kacem, J., 2004, Etude sismotectonique et évaluation de l’aléa sismique régional du Nord-Est da la
Tunisie - Apport de la géophysique dans l’identification des sources sismogéniques: Tunis,
Univ. Tunis El Manar.
Kacem, J., Dlala, M., & Hfaiedh, M., 2001, Sismicité de la région de Sidi Khalifa et ses environs:
Geosc. Rev.
Kammoun, Y., 1981, Etude néotectonique de la région de Monastir-Mahdia (Tunisie orientale),
Univ. Paris XI.
Kamoun, Y., Sorel, D., Viguier, C., & Ben Ayed, N., 1980, Un grand accident subméridien d'âge
post-Tyrrhénien en Tunisie orientale: le décrochement sénestre de Skanès (Monastir) -
Hammamet.: C.R. Acad. Sc. Série D, v. 290, p. 647-649.
Karnik, V., 1969, Seismicity of the European Area, Volume 1-2, D. Reidel, Dordrecht, p. 364.
Kaufman, A., Broecker, W.S., Ku, T.L., & Thurber, D.L., 1971, The status of U-series methods of
mollusc dating: Geochimica et Cosmochimica Acta, v. 35, p. 1155-1183.
Kreemer, C., & Holt, W.E., 2001, A no-net-rotation model of present-day surface motions:
Geophysical Research Letters, v. 28, p. 4407-4410.
Laftine, R., & Dupont, E., 1948, Plissement Pliocène supérieur et mouvement quaternaire en
Tunisie: C.R. Acad. Sc., v. 227, p. 138-140.
Lambeck, K., & Nakada, M., 1992, Contraints on the age and duration of the Last Interglacial
period and on sea-level variations: Nature, v. Vol. 357, p. PP.125-128.
Martinez, C., & Paskoff, R., 1984, Indices de distension pendant le Quaternaire récent en Tunisie:
leur signification dans un régime de compression généralisée: Cagiers Orstom, ser. Geol.,
v. 14, p. 153-161.
Mascle, G.H., Tricart, P., Torelli, L., Buillin, J.P., Rolfo, F., Lapierre, H., Monie, P., Depardon, S.,
Mascle, J., Peis, D., & . 2001, Evolution of the Sardinia Channel (Western Mediterranean) :
new constraints from a diving survey on Cornacya seamount off SE Sardinia: Marine
Geology, v. 179, p. 179-202.
Mattauer, M., Tapponnier, P., & Proust, F., 1977, Sur les mécanismes de formation des chaînes
intracontinentales : l’exemple des chaînes atlasiques du Maroc: Bull. Soc. géol. France, v.
19, p. 521-526.
Maury, R.C., Fourcade, S., Coulon, C., El Azzouzi, M., Bellon, H., Alain Coutelle, A., Ouabadi, A.,
Semroud, B., Megartsi, M., Cotten, J., Belanteur, O., Louni-Hacini, A., Piqué, A., Capdevila,
R., Hernandez, J., & Réhault, J.P., 2000, Post-collisional Neogene magmatism of the
Mediterranean Maghreb margin: a consequence of slab breakoff: C. R. Acad. Sci. Paris, v.
331, p. 159-173.
McGuire, R.K., 1993, Computations of seismic hazard: Ann. Geophys, v. 26, p. 181-200.
Meghraoui, M., 1988, Géologie des zones sismiques du nord de l'Algérie: Paléosismologie,
tectonique Active et synthèse Sismotectonique: Centre d'Orsay, Université de Paris sud.
Meghraoui, M., & Doumaz, F., 1996, Earthquake-induced floodind and palesseismicity of the El
Asnam, Algeria, fault-related fold rupture zone (6,8 Algeria): Aftershoks sequence reloca:
J. Geophys. Res., v. 101.
Mejri, L., 2004, tectonique récente au Nord-Est de la Tunisie: étude de la structure de Messeftine
[D.E.A thesis]: Tunis, Faculté des Sciences.
Mejri, L., Regard, V., Carretier, S., Brusset, S., & Dlala, M., 2010, Evidence of Quaternary active
folding near Utique (Northeast Tunisia) from tectonic observations and a seismic profile:
C. R. Geoscience, v. 342, p. 864-872.
Meulenkamp, J.E., & Sissingh, W., 2003, Tertiary palaeogeography and tectonostratigraphic
evolution of the Northern and Southern Peri-Tethys platforms and the intermediate

171
 domains of the African-Eurasian convergent plate boundary zone: Palaeogeography,
Palaeoclimatology, Palaeoecology, v. 196, p. 209-228.
Miosec, A., 1977, Le littoral du pays de Bizerte, étude géomorphologique: Rev. Géol. Medit., v. 6,
p. 281-290.
Mogami, T., & Kubo, K., 1953, The behavior of soil during vibration, 3 rd International Conference
on Soil Mechanics and Foundation Engineering, Volume 1: Zurich, p. 152-155.
Morel, J.-L., & Meghraoui, M., 1996, Gorringe-Alboran-Tell tectonic zone: A transpression system
along the Africa-Eurasia plate boundary: Geology, v. 24, p. 755-758.
Nedjari-Belhacene, S., & Nedjari, A., 1984, Le Boukhadra: esquisse paléogéographique aptienne.,
5ème séminaire national des polytechniques de Lorraine, p. 178.
Newman, A., Schneider, J., Stein, S., & Mendez, A., 2001, Uncertainties in seismic hazard maps for
the New Madrid Seismic Zone and implications for seismic hazard communication: Seism.
Res. Lett., v. 72, p. 647-663.
Nocquet, J.M., & Calais, E., 2004, Geodetic measurements of crustal deformation in the western
Mediterranean and Europe: Pure Appl. Geophys., v. 161, p. PP. 661-681.
Noquet, J.M., & Calais, E., 2003, Crustal velocity field of western Europe from permanent GPS
array solutions: Geophys. J. Int., v. 154, p. 72-88.
Oldow, J.S., Channel, J.E.T., Catalona, R., & D'Argonio, B., 1990, comtemporaneous thrusting and
large retations in the Western Sicilian fold and thrust belt: Tectonophysics, v. 9, p. 61-681.
Olivet, J.L., Auzende, J.M., & Beuzart, P., 1983, Western extension of the Iberian-European plate
boundary during the Early Cenozoic (Pyrenean) convergence: A new model — Comment:
Marine Geology, v. 53, p. 237-238.
Oueslati, A., 1985, Les dépôts holocènes liés au ruissellemnt sutr les côtes de la Tunisie: Notes Ser.
Géol. Tunisie, v. 51, p. 121-131.
—, 1987, Exemple de plates-formes quaternaires polygéniques et polychroniques sur la côte nord
de la Tunisie: Rev. Tun. Géogr., v. 15, p. 139-158.
—, 1993, Les côtes de la Tunisie, recherches sur leur évolution au Quaternaire, Volume Série 2 :
Géographie Publications de la faculté des sciences humaines et sociales de Tunis, p. 402p.
Paskoff, R., & Sanlaville, P., 1979, Livret guide de l'excursion table-ronde sur le Tyrrhénien de
Tunisie / Union Internationale pour l'étude du Quaternaire, commission des lignes de
rivages, sous commissions Méditerranéenne, Mer noir, p. 51.
—, 1980, LeTyrrhénien de la Tunisie essai de stratigraphie: C. R. Acad. Sci. Paris, v. 290, p. 393-
396.
—, 1982, Sur les dépôts Tyrrhéniens et wurmiens des littoraux de la Méditerrannée occidentales:
C. R. Acad. Sci. Paris, v. 294, p. 737-740.
—, 1983, Les côtes de la Tunisie : variation du niveau marin depuis le Tyrrhénien: Lyon, Éditions
Maison de l’Orient,, 192 p.
Pepe, F., Bertotti, G., Cella, F., & Marsella, E., 2000, Rifted margin formation in the south
Tyrrhenian sea : A high-resolution seismic profile across the north Sicily passive
continental margin: Tectonics, v. 19, p. 241-257.
Pepe, F., Sulli, A., Bertotti, G., & Catalano, R., 2005, Structural highs formation and their
relationship to sedimentary basins in the north Sicily continental margin (southern
Tyrrhenian sea) : Implication for the Drepano thrust front: Tectonophysics, v. 409, p. 1-18.
Perthuisot, V., 1978, dynamique et pétrogenèse des extrusions triasiques en Tunisie
septentrionale: Paris, Ecole Normale supérieure

Philip, H., 1987, Plio-Quaternary evolution of the stress field in Mediterranean zones of
subduction and collision: Annales Geophysicae, v. XVIII, p. 301-320.
Philip, H., Andrieux, J., Dlala, M., Chihi, L., & Ben Ayed, N., 1986, Evolution tectonique mio-
plioquaternaire du fossé de Kasserine (Tunisie centrale) : implications sur l’évolution
géodynamique récente de la Tunisie: Bull. Scoc. géol. France, v. 4, p. 559-568.

172
Philip, H., Bousquet, J.C., & Masson, F., 2007, Séismes et risques sismiques, approche
sismotectonique, DUNOD, 339 p.
Phlip, H., & Meghraoui, M., 1983, Structural analysis and interpretation of the surface
deformations of the el Asnam Earthquake of october 10, 1980: Tectonics, v. 2, p. 17-49.
Piementa, J., 1959, Le cycle pliocène-actuelle dans les bassins paralyques de Tunis [Thèse de la
Faculté des Sciences de Paris et Mémoire de la Société Géologique de France thesis],
Paris.
Pini, S., De Jauzein, C.A., & Kchouk, F., 1971, Carte géologique de la Tunisie au 1/50000 et notice
explicative de la feuille n°13 : Ariana: Publ. Serv. Géolo. Tunisie.
Rebaï, S., 1992, sismotectonique et champ de contrainte dans les chaînes alpines et les
plateformes de l'Europe, d'Afrique du Nord et du Moyen-Orient, Montpellier II.
Reiter, L., 1990, Earthquake hazard analysis : issues and insights, New York : Columbia University
Press, 254 p.
Reuther, C.D., & Eisbacher, G.H., 1985, Pantellaria Rift-Crustal extention in a couvergent intra-
plate setting: Geoligische-Rundschau, v. 74.
Reyners, F.L.M., 1951, Note sur le comblement de la Basse Medjerdah: Note Serv. Géol. Tunisie, p.
1-18.
Roure, F., Howell, D.L., Muller, C., & Moretti, I., 1990, Late Cenezoique subduction complex of
Sicily: J. of struc. Geology, v. 12, p. 259-266.
Rouvier, H., 1971, Le problème de la partie du "Numidien" de la Tunisie Septentrionale, Ibid, p.
196-206.
—, 1977, Géologie de l’extrême Nord tunisien, Thèse Doct. Es-Sc., Univ. Pierre et Marie Curie

Schwartz, D., & Coppersmith, K.J., 1984, Fault behavior and characteristic earthquakes: Examples
from the Wasatch and San Andreas fault zones: J. Geophys. Res., v. 89(B7), p. 5681-5698.
Schwartz, D.P., 1988, Geologic characterization of seismic sources: moving into the 1990's: in von
Thun, L., ea., Earthquake engineering and soil dynamics II--Recent advances in ground
motion evaluation: American Society of Civil Engineers, Geotechnical Special Publication,
v. 20, p. 1-42.
Sella, G.F., Dixon, T.H., Mao, A., & . 2002, REVEL: A model for recent plate velocities from space
geodesy: J. Geophys. Res.
Sharp, W.D., Ludwig, K.R., Chadwick, O.A., Amundson, R., & Glaser, L.L., 2003, Dating fluvial
terraces by Th-230/U on pedogenic carbonate, Wind River Basin, Wyoming.: Quat. Res., v.
59, p. 139-150.
Sieberg, A., 1932, Erdeben geographie: Handbuch der Geophysik, v. 4, p. 687-1005.
Skobelev, S.F., Trifonov, V.G., & Vostrikov, G.A., 1988, The Pamirs-Himalayan region of
disharmonic clustering of the continental lithospheren, in Kropotkin, P.N., ed.,
Neotectonics and the Recent Geodynamics of Mobile Belts: Moscow, Nauka Press,
Moscow, p. 188-234.
Slemmons, B.D., & De Polo, C.M., 1986, Evaluation of active faulting and related hazard, Stadies in
Geophysics: Active Tectonics, National Academy Press, Washington, DC, p. 45-62.
Solignac, M., 1927, Etude géologique de la Tunisie septentrionale, Dir. Gen. Tra. Pub. Tunis, p.
756p.
Suppe, J., & Medwedeff, D.A., 1990, Geometry and kinematics of fault-propagation folding:
Eclogae geol. Helv, v. 83, p. 409-454.
Tricart, P., Torelli, L., Argnani, A., Rekhiss, F., & Zitellini, N., 1994, Extensional collapse related to
compressional uplift in the alpine chain off northern Tunisia (Central Mediterranean):
Tectonophysics, v. 1, p. 317-329.
Trifonov, V.G., 2004, Active faults in Eurasia: general remarks: Tectonophysics, v. 380, p. 123-130.
Vaufrey R, 1932, Le plissement Acheulo-Moustérien des environs de Gafsa: Revue, Géogr., Phys.
Dyn., v. V1, p. pp. 209-325.

173
Viguier, C., Ben Ayed, N., & Bobier, C., 1980, Précisions sur les phases de la tectnoque récente en
Tunisie Nord-Orientale: Bull. [Link]. Nat. Tunis, p. pp. 63-68.
Vogt, J., 1993, Further research on the historical seismicity of Tunisia: Terra Nova, v. 5, p. PP 475-
476.
Wortel, M.J.R., Goes, S.D.B., & Spakman, W., 1990, Structure and seismicity of the Aegean
subduction zone: Terra Nova, v. 2, p. 554-562.
Youd, T.L., 1984, Recurrence of liquefaction at the same site, Eighth World Conference on
Earthquake Engineering, Volume 3: San Francisco, California, p. 231-238.
Youngs, R.R., & Coppersmith, K.J., 1985, Implications of fault slip rates and earthquake recurrence
models to probabilistic seismic hazard estimates: Bulletin of the Seismological Society of
America, v. 75, p. 939-964.
Zargouni, F., & Ruhland, M., 1981, Style de déformation du Quaternaire récent lié au coulissement
de la faille de Gafsa, et chronologie des phases tectoniques de l’Atlas méridional tunisien:
C. R. Acad. Sci., Paris, v. 292, p. 912-915.
Zouari, H., 1995, Les variations géométriques du plissement liées à la cinématique coulissante,
faculté des sciences de Tunis, p. 278.

174
ANNEXES

175
176
177
178
179
180
181
182
183
184
185