La structure actuelle des différents domaines de l’arc Bético-Rifain est le résultat d’une tectonique polyphasée. La richesse des observations réalisées dans la Méditerranée occidentale a abouti à la proposition d’une large gamme de modèles et d’hypothèses, qui tentent d’expliquer les phénomènes orogéniques et post-orogéniques dans cette zone. Les principaux modèles se basent essentiellement sur des données géologiques et géophysiques (par exemple, [2,6,9,11,19,25,36,37,40,49,54]).

La surveillance continue ou périodique par GPS pendant plusieurs années permet d’établir le champ de vitesses des mouvements tectoniques, avec une précision millimétrique, sur les mouvements horizontaux. La combinaison du GPS avec la tectonique et la sismologie, selon la disponibilité des données, permet d’identifier les failles actives. Le changement d’intensité et de direction des contraintes appliquées sur des blocs adjacents [46] peut provoquer le blocage de certaines failles, allant même jusqu’à induire l’apparition de nouvelles ruptures ou réactiver d’anciennes failles.

L’accommodation des contraintes peut être effectuée par des failles majeures délimitant des blocs crustaux, ou bien par un ensemble de failles distribuées sur toute une région [55]. La déformation pendant la formation de l’arc de Gibraltar a été partagée entre des structures transpressives (plis et chevauchements) traduisant un raccourcissement nord–sud et des structures extensives (failles normales, décrochements conjugués et structures mineures de déformation distribuée) indiquant une extension est–ouest [3].

L’anomalie gravimétrique à l’air libre peut apporter des informations clés pour caractériser des blocs lithosphériques en fonction de leur densité. En effet, le gradient de l’énergie potentielle gravitationnelle (GPE) entre des blocs structuraux peut être à l’origine de déplacements aussi bien horizontaux que verticaux [34]. Dans l’arc Bético-Rifain, le contraste du GPE entre le domaine d’Alboran et ses avant-pays (Afrique et Ibérie) est considéré comme un facteur important à l’origine de la subsidence du bassin d’Alboran, ainsi que du raccourcissement nord–sud et de l’extension est–ouest de l’arc au Miocène [44].

Dans cet article, nous présentons les résultats de mesures GPS répétées pendant 6 ans (1999–2005) sur des sites répartis sur les différents domaines structuraux du Maroc septentrional. Vu l’importance des mouvements observés dans le Rif, nous nous focalisons sur cette zone de frontière de plaques. Le but de cette étude est de donner une description quantitative de la cinématique actuelle de cette zone de convergence entre l’Afrique et l’Eurasie. Nous présentons de nouvelles observations GPS dans le Rif, qui viennent compléter le travail de Fadil et al. [20]. Nous tentons d’établir ici le lien entre les déformations observées en surface et, d’une part, l’activité sismique pendant la même période d’acquisition des mesures GPS, d’autre part, les anomalies gravimétriques à l’air libre, pour aider à différencier les blocs crustaux évoqués dans les différents modèles cinématiques régionaux existants. Enfin, nous discutons une nouvelle interprétation des observations.

Situé à la bordure sud de la Méditerranée occidentale, le Maroc fait partie d’une zone complexe de convergence entre les plaques Africaine et Eurasiatique. Les différentes chaînes de montagnes actuelles (Atlas, Rif, cordillères Bétiques) sont le résultat d’une longue histoire d’interaction complexe entre ces deux grandes plaques tectoniques. Pendant cette histoire, différents blocs crustaux ont été pris dans les processus géodynamiques qui ont affecté cette région. Nous nous limiterons ici à donner un bref aperçu des différents domaines structuraux couverts par notre réseau GPS marocain, ainsi que les stations GPS permanentes autour du domaine d’Alboran.

Du sud vers le nord, on distingue, en premier lieu, les domaines de l’Anti-Atlas, des Atlas (Haut et Moyen) et des Mesetas (orientale et occidentale), qui ont souvent été considérés comme des blocs solidaires de la partie occidentale de la plaque Africaine (Nubie). Ils sont délimités par des accidents majeurs, dont certains montrent une activité tectonique récente [50]. La discussion sur les mobilités (verticale et horizontale) relatives de ces blocs, depuis le cycle orogénique varisque [23,30] jusqu’à l’Actuel, reste ouverte. Les autres domaines du Nord du Maroc (Rif) et du Sud de l’Espagne (Bétiques) sont liés géodynamiquement au domaine d’Alboran, situé entre l’Afrique et l’Ibérie. Le domaine Bético-Rifain, structuré pendant le cycle alpin [10,40], est découpé en de nombreuses écailles tectoniques très fracturées, avec un grand contraste rhéologique entre ses unités. Vers le sud, les nappes du Rif chevauchent la marge africaine proprement dite. Vers le nord, les nappes des Bétiques chevauchent la marge de l’Ibérie. Dans cette situation, le Rif et les Bétiques sont soumis à des contraintes liées à la convergence oblique de direction NW–SE entre l’Afrique et l’Eurasie, au niveau du domaine d’Alboran [38]. La majorité de ces contraintes est accommodée par les structures actives de cette zone qui, dans la plupart des cas, sont des décrochements et/ou des failles à faible pendage préexistant entre les nappes [1,39]. La vitesse de cette convergence dans cette région a été estimée entre 3 et 5 mm/an [8,14,38].

La distribution des séismes permet de situer les accidents de la croûte potentiellement actifs (Fig. 1). La sismologie historique et instrumentale montre que le Maroc présente une sismicité diffuse affectant presque tous les domaines structuraux [7,17,18]. Dans la région d’Alboran, la sismicité de magnitudes Mw supérieures à 5 et d’intensités supérieures à VII (MSK) présente une fréquence exceptionnelle dans (1) la zone d’Al Hoceima, (2) la zone de Nador–Melilla, (3) la zone d’Almeria et (4) la zone de Malaga–Granada.

Dans le but d’évaluer les mouvements tectoniques dans la partie nord-occidentale de l’Afrique, un réseau GPS a été mis en place à partir de 1999. Lors d’une première étape, 22 points géodésiques (repères géodésiques ancrés dans des roches en place) ont été déployés sur tous les grands domaines structuraux du Maroc, séparés par des failles majeures. En 2002, le sous-réseau du Nord du Maroc a été densifié par l’installation de 6 stations de campagne. Le réseau GPS marocain compte actuellement 28 stations de mesures en campagne, renforcées par quatre stations permanentes (RABT, TETN, IFRN et CEUT). Les mesures en campagne ont été effectuées entre 1999 et 2005 pendant des sessions de 24 h (Tableau 1).

La tectonique active de l’Anti-Atlas, en bordure nord de la plate-forme Africaine, est observée (Fig. 2) par les trois stations BAHA, ZARA et ERFD, alignées en direction SW–NE selon l’alignement général de la chaîne. Le Haut Atlas, dont l’altitude dépasse 4000 m, est aligné SW–NE dans sa partie occidentale à est–ouest dans sa partie orientale. Il est découpé par le système de failles du Tizi n’Test, d’orientation NE–SW. Cette grande structure du Maroc est observée par 4 stations IMZR, MAR0, AZIL et RIC0. Le Moyen Atlas, qui se sépare obliquement du Haut Atlas au niveau de Beni Mellal, est surveillé en permanence par la station IFRN et en campagnes de mesures par les deux stations HEBR et MBLD. La Meseta côtière, considérée comme plus ou moins stable, est contrôlée par la station permanente RABT située juste sur sa bordure nord, où commence la flexure de la plaque africaine sous le bassin du Gharb, ainsi que par trois stations de mesure en campagne, SALA, KBGA et ONRA. En Meseta orientale, plateau qui se prolonge en Algérie, les déformations sont observées par 3 stations DEBD, BMTR, TNDR. Le Rif, considéré comme un domaine de haute complexité, du fait de sa situation dans la zone d’interaction des grandes plaques Afrique et Eurasie, est instrumenté par un réseau de 15 sites, dont deux stations permanentes TETN et CEUT et 13 stations de campagne (TNIN, LAOU, BBFH, MDAR, MSLA, ZAGO, CHEF, OUZS, BOYA, TANT, KTMA, GRSF, AION).

Les données GPS des six années de mesures ont été traitées à l’aide des logiciels GAMIT/GLOBK [29,35], selon la même stratégie que celle utilisée par Fadil et al. [20]. Les données des stations IGS d’Afrique et d’Europe ont été introduites dans les calculs des vitesses relatives afin de réduire les incertitudes (écarts-types inférieurs à 1  mm/an). Nous montrons le champ de vitesses dans la Méditerranée occidentale par rapport à l’Eurasie fixe, d’une part (Fig. 2), et par rapport à l’Afrique fixe, d’autre part (Fig. 3). Des vitesses de déplacement horizontal anormales ($3\text{,}4\pm 1\text{,}2$  mm/an vers le sud par rapport à l’Afrique) dans la partie centrale du Rif ont été observées par Fadil et al. [20] et sont confirmées ici par les vitesses obtenues à partir des mesures sur la période de 2002 à 2005 de six stations supplémentaires (CHEF, BOYA, TANT, KTMA, GRSF, AION). Un bloc individualisé par le sens de déplacement des stations du Rif central (LAOU, CHEF, OUZS, BOYA, TANT, KTMA, MSLA et ZAGO) montre un mouvement général vers le SSW par rapport à l’Afrique fixe. Le mouvement des autres stations, non significatif, montre la faiblesse du signal tectonique. Le mouvement du bloc central du Rif par rapport à l’Europe fixe (Fig. 2) apparaît dévié vers l’WSW en tenant compte du mouvement de l’Afrique par rapport à l’Eurasie indiqué par les autres stations GPS du Maroc. Par rapport à l’Europe fixe, on remarque aussi le mouvement vers l’WNW de la station SFER située à l’ouest des Bétiques. Les stations CEUT, TETN et TNIN indiquent un mouvement vers le sud-est par rapport à l’Afrique stable. Ces trois stations définiraient un bloc à part (bloc de Tanger), qui semble être intermédiaire entre le mouvement de l’Ibérie et celui de l’Afrique.

Le GPS fournit les déplacements à la surface de la terre selon ses trois composantes. Tandis que les deux composantes horizontales sont évaluées avec des incertitudes raisonnables (d’ordre millimétrique), la composante verticale reste, cependant, faiblement cernée avec des incertitudes deux à trois fois supérieures par rapport aux incertitudes horizontales.

Trois profils associant les vitesses GPS (composantes horizontales parallèles et perpendiculaires aux profils), les anomalies gravimétriques à l’air libre, la topographie et la sismicité (couvrant la même période d’observation que les mesures GPS) ont été réalisés (Fig. 4). La combinaison des données a été faite dans le souci de confronter les mouvements horizontaux observés par GPS aux différentes données de profondeur : (1) l’anomalie gravimétrique à l’air libre, liée à la distribution de masse, peut s’interpréter en termes de contraintes et de déformations induites par les processus orogéniques [31] ; (2) la distribution de la sismicité en profondeur sur la même période d’observation indique si les mesures GPS ont pu être affectées par des déplacements cosismiques lors des mesures.

Le profil P1 (Fig. 5) est parallèle à la direction générale du mouvement de l’Afrique. Il caractérise les interactions du bloc central du Rif avec l’ensemble du domaine rifain et l’avant-pays. Ce profil suggère une diminution des vitesses perpendiculaires au profil entre SFER (Ibérie) et TETN–TNIN (péninsule de Tanger). Le bloc central du Rif est marqué par des vitesses plus élevées. Le mouvement quasi nul de BOYA (rides pré-rifaines) perpendiculairement au profil peut être interprété comme un blocage de l’avancée des nappes de charriage vers le SSW sur l’avant-pays. La station ZAGO, située elle-même sur une ride pré-rifaine, semble être bloquée dans les deux directions. Le système atlasique est, quant à lui, presque stable le long du profil. Par ailleurs, le bloc central du Rif est caractérisé par une anomalie gravimétrique négative importante ($-$100 mGal), qui peut être liée à la flexure de la plaque africaine sous l’effet du mouvement continu des nappes de charriage vers le SSW. En contraste du Rif, le système des Atlas est caractérisé par une anomalie gravimétrique positive importante (+100 mGal), en dépit des hautes altitudes. Un modèle de lithosphère fortement amincie a été proposé par [22] au niveau des Atlas, en utilisant les données de gravimétrie, le géoïde, les flux de chaleur et la topographie.

Les profils P2 et P3 (Figs. 6 et 7) montrent des variations latérales des vitesses GPS, des anomalies gravimétriques à l’air libre et de sismicité en profondeur entre les domaines structuraux du Maroc, le domaine d’Alboran et l’Ibérie. Sur ces deux profils, les vitesses semblent diminuer de l’Ibérie vers le bloc central du Rif (P2 et P3) et du Rif central vers les hauts plateaux (P3). La limite sud du bloc central du Rif (entre ZAGO et IFRN) est marquée par un fort gradient de l’anomalie gravimétrique à l’air libre, qui atteint des minima de $-$100 mGal (P2) ; la limite est coïncide avec un plan de sismicité à fort pendage vers l’ouest. Le vecteur vitesse calculé pour la station MDAR est différent de ceux des stations du Rif central qui sont situées plus à l’Ouest. Cette observation unique sur le Rif oriental est en accord avec la rotation horaire mise en évidence par Najid et al. 1981 (in [40]). L’activité sismique au niveau de la région Nador–Melilla [5] témoigne de la persistance actuelle des mouvements tectoniques étudiés par [16] dans cette zone de bordure entre le Rif et l’avant-pays.

La formation de l’arc Bético-Rifain est liée à l’histoire du découplage/couplage entre l’Afrique, l’Ibérie et le domaine d’Alboran, où l’on distingue deux périodes : (1) au Miocène, alors que les deux plaques sont plutôt découplées, le roll-back du slab de subduction de la plaque océanique Téthys permet le mouvement vers l’ouest du domaine d’Alboran et le développement du prisme d’accrétion aux dépens des marges passives de l’Afrique et de l’Ibérie (par exemple, [11,25,22]). Par conséquent, l’ensemble du Rif externe est déplacé vers l’ouest par l’intermédiaire des rampes latérales des Temsamane, puis du Nekor [21,25]. La structuration du Rif externe se caractérise par des chevauchements majeurs vers le sud délimitant les différentes unités (par exemple, [22]). (2) Depuis le Pliocène, le couplage entre les plaques devient plus fort, consécutivement aux diverses ruptures et déchirures des slabs en Méditerranée occidentale [19,51]. Les effets de la convergence sont alors transmis vers l’Atlas (l’avant-pays africain de la chaîne du Rif) [24]. Il est à noter l’asymétrie entre les Bétiques et le Rif (zones externes), mise en évidence par des différences structurales importantes [13]. Pendant que le Rif a été régi par une tectonique d’imbrication socle–couverture du pré-Oligocène au Miocène, les Bétiques ont subi une tectonique de couverture au Miocène.

Les modèles régionaux voulant expliquer la tectonique active actuelle de la zone Bétiques–Alboran–Rif peuvent être subdivisés selon trois mécanismes dynamiques principaux : (1) délamination et roll-back au niveau des Bétiques et du Rif [9,20], (2) subduction vers l’est au niveau de l’arc de Gibraltar [27] et (3) coexistence d’une extension ENE–WSW au niveau du bassin d’Alboran, induite par des processus dans le manteau supérieur, et d’un coulissement de la région vers l’ouest lié à la rotation de la plaque africaine [53]. Fadil et al. [20] ont discuté la limite de plaques Eurasie–Afrique et ont proposé un modèle de bloc élastique pour déterminer les conditions aux limites en terme de vitesses de glissement le long des failles qui séparent les blocs.

À la lumière des nouvelles observations que nous présentons, nous souhaitons confronter nos résultats avec les modèles cinématiques proposés. La Fig. 8 montre un schéma structural interprétatif de la situation tectonique actuelle dans l’arc Bético-Rifain. Nous proposons une synthèse des mouvements de surface des stations GPS supplémentaires décrites dans ce travail, ainsi que les observations des travaux précédents.

Le bloc central du Rif se déplace plutôt vers le SSW par rapport à l’Afrique fixe. La limite nord-ouest de ce bloc correspond à la prolongation de l’accident transformant dextre de Tétouan [4] vers le sud-ouest, bordant les rides prérifaines dans la région de Larache (T–LF dans la Fig. 8), ce qui marque, dans cette région, le début du sillon Sud-Rifain, qui sépare le Rif du continent africain. Le bloc de Tanger s’individualise en accord avec les conditions cinématiques actuelles par un soulèvement prononcé par rapport au bloc central du Rif, comme semblent le montrer des observations sur le littoral pour le Quaternaire [26,45]. Nous suggérons la mise en place d’un accident complexe (T–LF), de direction générale WSW–ENE, entre la région de Tétouan et la région de Larache. Cet accident serait constitué de plusieurs failles décrochantes dextres, héritées des périodes orogéniques antérieures. L’accident pourrait être une rampe latérale (oblique par rapport au bloc central du Rif), avec un mouvement décrochant dextre à composante inverse, permettant à la fois le soulèvement du bloc de Tanger et le déplacement des nappes du bloc central vers le SSW par rapport à l’Afrique. Il faut noter cependant qu’aucun décalage du front de chevauchement des zones internes n’a été révélé jusqu’à présent. Cela nous permet de suggérer deux hypothèses : il s’agit, soit (1) d’une structure très jeune découpant le dispositif antérieur limité plutôt par la faille de Jebha, soit (2) de la réactivation d’un accident antérieur ayant les caractéristiques d’une faille transformante [4]. Dans les deux hypothèses, la structure pourrait définir une zone de déformation distribuée. Une étude détaillée de la géomorphologie et de la microsismicité paraît indispensable pour localiser respectivement la trace en surface et la géométrie en profondeur des différents segments.

La limite est du bloc central du Rif semble être confirmée par la sismicité au niveau de la région d’Al Hoceima et semble se prolonger vers le sud selon une direction N190 à N200. Cette structure semble coïncider avec la faille décrite dans des travaux comme ceux de [15,43]. Le bloc oriental du Rif, qui est en contact direct avec le bloc central par la faille sénestre d’Al Hoceima (AF), est séparé de l’avant-pays par la faille de Gareb–Kebdana (G–KF), située dans une zone sismique de décrochements sénestres (par exemple, [5]). Le bloc se trouve en un régime trans-rotationnel sous les effets conjugués de la convergence entre l’Eurasie et l’Afrique et de la rotation antihoraire de cette dernière. La convergence NW–SE transférerait un mouvement vers l’est à la bordure nord, tandis que la rotation transférerait un mouvement vers l’ouest à la bordure sud du bloc oriental du Rif. Cela se manifeste par le jeu sénestre de ses deux bordures (AF et G–KF) sur sa partie émergée. Les variations latérales des anomalies gravimétriques entre les blocs dans l’arc Bético-Rifain peuvent être à leur tour à l’origine des mouvements horizontaux observés dans ce travail. Mais cela ne peut être précisé que par la modélisation des contraintes déviatoriques dans la région, en déterminant leurs sens et amplitudes. La structure complexe du Rif, combinant des nappes de charriage limitées par des rampes latérales à fort pendage, doit être prise en compte pour déterminer le mode d’accommodation des contraintes, observé en surface par les composantes horizontales du GPS.

Les blocs du domaine rifain, distingués dans ce travail, montrent à la fois des accords et des désaccords par rapport aux trois types de mécanismes présentés plus haut. Le mouvement vertical et horizontal vers le sud-est du bloc de Tanger implique une «force» poussant dans ce sens et que nous pourrions lier à l’effet de la subduction suggérée par [27] au niveau de l’arc de Gibraltar. Cependant, ces mouvements sont incompatibles avec le roll-back vers l’ouest proposé par ces mêmes auteurs. De même, il serait possible d’interpréter l’individualisation de ce bloc par l’effet du coulissement de blocs le long de la zone de faille transformante Açores–Gibraltar (AGFZ) entre les deux grandes plaques Afrique–Eurasie. Le mouvement particulier du bloc central est fort étonnant du point de vue de son sens (SSW) et de sa vitesse horizontale (3,4 $\pm$ 1,2 mm/an). Cela confirme les observations de Fadil et al. [20], mais nous restons prudents pour confirmer leur interprétation en ce qui concerne la délamination et le roll-back au niveau du Rif. Nous pouvons toutefois évoquer une explication alternative liée à une extrusion crustale [12,41,54] de ce bloc en réponse aux contraintes de convergence Afrique–Eurasie. Le bloc oriental du Rif, faisant partie de l’ensemble structural rifain, s’individualise impérativement, d’une part, par ces limites sismogènes au niveau de la région d’Al Hoceima et au niveau de la région de Nador–Melilla et, d’autre part, par sa rotation horaire sous l’effet du régime trans-rotationnel imposé par les mouvements des deux grandes plaques lithosphériques Afrique et Eurasie. Le modèle général de Stich et al. [53] cité plus haut, basé sur les données de la sismologie, semble être plutôt un cadre très général des mouvements des blocs dans cette même région.

Les déplacements significatifs du Rif central impliquent la mise en place d’une zone de découplage, suggérant la présence d’un décollement plus profond que les niveaux d’imbrication des unités structurales d’origines différentes (zones internes et externes). L’amortissement des mouvements sur l’avant-pays (Afrique) permet de suggérer que le front actif sud du Rif coïncide avec le front de chevauchement de la chaîne par rapport à la plaque africaine. Cela indiquerait un découplage entre la couverture et le socle, mais la lacune liée à la résolution verticale de nos données ne nous permet pas de déterminer sa géométrie. Cependant, dans l’hypothèse d’un roll-back du manteau lithosphérique délaminé sous le Rif [20], le découplage ne pourrait être que partiel (voir inexistant). Il faut également noter que la lacune de résolution des données au sud de l’Ibérie ne nous permet pas de donner une interprétation précise sur le front de déformation active au nord des Bétiques. Les profils ne montrent aucun changement significatif des vitesses GPS de part et d’autre du bassin de Guadalquivir.

Expliquer les déformations par un mécanisme unique devient de plus en plus difficile à adopter. La réconciliation des modèles testés par des méthodes différentes s’impose pour expliquer le mode de déformation actuelle dans la région. La combinaison des méthodes d’évaluation des forces (compressives ou extensives) intrinsèques et extrinsèques et des déformations qui en résultent permet une meilleure estimation des mécanismes de déformation. De même, la sommation des contributions élémentaires à l’échelle locale est essentielle pour lisser les différences d’interprétation régionale. Dans notre travail, le nombre limité des stations GPS dans la région ne permet pas de préciser des processus géodynamiques responsables des déformations observées face à une large gamme de modèles proposés. Néanmoins, nous insistons sur le fait que ces mouvements dans la région de l’arc Bético-Rifain sont le résultat du jeu combiné de la compression due à la convergence Afrique–Eurasie et de trois processus essentiels : (1) la rotation antihoraire de l’Afrique, (2) l’extension syncollisionelle au niveau de la mer d’Alboran, ainsi que (3) les contraintes déviatoriques liées à l’hétérogénéité géologique des blocs impliqués dans cette région. La coexistence de ces trois processus peut élargir la zone de déformations entre les plaques dans la région.

L’interprétation de nouvelles données GPS, associées aux données gravimétriques et sismologiques, ainsi que leur comparaison avec des résultats antérieurs permettent de confirmer la distinction d’un bloc central au Rif [20], caractérisé par un mouvement de l’ordre de $3\text{,}4\pm 1\text{,}2$ mm/an vers le SSW. Le Rif semble être séparé en trois blocs, de caractéristiques cinématiques différentes. (1) Le bloc de Tanger, qui se trouve au cœur de l’arc de Gibraltar, présente un mouvement ($1\text{,}6\pm 1\text{,}2$ mm/an) vers le sud-est. Il est limité par une rampe latérale, que nous avons appelé faille de Tétouan–Larache (T–LF). (2) Le bloc central du Rif se déplace vers le SSW avec une vitesse apparemment plus élevée que celle des autres blocs. Il est soumis, d’une part, à la compression due à la convergence et, d’autre part, à l’extension au centre de la mer d’Alboran (bassin occidental). (3) Le bloc oriental du Rif est en rotation horaire et est délimité par la faille d’Al Hoceima (AF), qui le sépare du Rif central, et par la faille de Gareb–Kebdana (G–KF), qui le sépare de l’avant-pays. La vitesse des mouvements des blocs diminue progressivement des nappes rifaines vers l’avant-pays. Comme on l’a montré dans la discussion, les mouvements des différents blocs, caractérisés à partir des observations GPS dans l’arc Bético-Rifain, présentent des incompatibilités avec les modèles proposés pour cette région dans la littérature. Les sources de divergence entre nos observations et ces modèles sont liées aux types de données utilisées. D’où l’intérêt d’une poursuite des efforts d’observation impliquant différents types de données (GPS, sismologiques, gravimétriques, etc.) pour contribuer à comprendre au mieux la déformation actuelle de cette région.