La télédétection des conditions physiques du pergélisol dans trois buttes de pergélisol près d'Umiujaq au Nunavik (Québec) a été réalisée à l'aide de la tomographie électromagnétique forage-forage et surface-forage. Les objectifs de ces travaux géophysiques étaient d'étudier la propagation des ondes électromagnétiques dans un milieu physique partiellement gelé et de déterminer la cryostratigraphie des buttes de pergélisol. Des antennes à géométrie cylindrique de géoradar à ondes courtes de 50, 100 et 200 MHz ont été glissées dans des regards géophysiques secs mis en place préalablement à l'étude ou utilisées en surface. La cryostratigraphie générale des buttes de pergélisol étudiées définie à partir des tomogrammes de vitesse de propagation des ondes électromagnétiques correspond à celle obtenue de forages et de l'échantillonnage de ces buttes. La résolution des tomogrammes est meilleure aux plus hautes fréquences.

Après la donnée de quelques définitions et éléments chiffrés sur le risque sismique, la première partie traite de la sismicité en France et des règles parasismiques en vigueur et se termine par une étude bibliographique sur les divers types de ruptures observés sur les ouvrages en terre pendant les séismes. La seconde partie présente les différentes méthodes d'analyse de stabilité des pentes et les diverses hypothèses menant au calcul du coefficient sismique. Le calcul du facteur de sécurité pseudostatique est réalisé avec des adaptations des méthodes de Taylor et des perturbations. Les abaques tracés permettent de déterminer le facteur de sécurité pseudostatique ou l'accélération critique d'un talus de caractéristiques données, alors que pour les barrages à noyau incliné on étudie l'influence de l'inclinaison du noyau sur la résistance aux séismes. Le calcul du déplacement permanent est abordé dans la dernière partie où, après une synthèse des méthodes présentées dans la littérature, on détaille le calcul du déplacement d'un talus homogène en rupture circulaire et d'un barrage à noyau incliné. Le comportement du sol est supposé rigide-plastique et le bloc en glissement est soumis à un accéléroramme équivalent déterminé en fonction de la magnitude, de la distance à l'épicentre et de la période prédominante de l'accélérogramme réel.

L’une des problématique scientifique majeure en sismologie est de pouvoir estimer les mouvements du sol attendus en un site pour un futur séisme. L’objectif de cette thèse est de tester et de valider deux méthodes de simulation des mouvements du sol basées sur l’approche des fonctions de Green empiriques et d’apporter des éléments pouvant aider au développement d’une méthodologie de simulation en aveugle. Dans une première partie, une méthode de simulation basée sur une approche stochastique en point-source est validée sur les données réelles de séismes récents bien instrumentés : le séisme des Saintes Mw6.4 et le séisme de L’Aquila Mw6.3. Nous avons développé une approche de simulation en aveugle en prenant en compte une incertitude sur le paramètre de rapport des chutes de contrainte C. Cette approche permet de générer un ensemble d’accélérogrammes synthétiques d’un séisme cible suffisamment variés pour être représentatifs d’un grand nombre de scénarios de sources possibles et prenant en compte dans un sens statistique de potentiels effets de directivité. Cette approche a également été appliquée à la simulation d’un séisme historique pyrénéen Mw6.1. Dans une seconde partie, nous nous appuyons sur un modèle de source étendue plus complexe, combinant des modèles cinématiques de sources composites fractales avec l’approche des FGEs. Le potentiel de la méthode est testé sur une application au séisme de L’Aquila. Cela a permis de produire des résultats très satisfaisants sur l’ensemble des paramètres des mouvements du sol analysés. Cette méthode de simulation apparaît comme étant très prometteuse pour la mise en œuvre d’une méthodologie de simulation en aveugle, même si la principale difficulté réside dans la nécessité de définir la variabilité de nombreux paramètres d’entrée mal connus dans le cadre de la simulation d’un futur séisme. / One of the major scientific problems in seismology is to estimate the ground motions expected at a given site from a future earthquake. The aim of this thesis is to test and validate two different methods of ground motions simulation based on the empirical Green’s function approach and to provide elements that can help to develop a blind simulation methodology. In a first part, a simulation method based on a stochastic point source approach is validated on the real data of recent earthquakes well instrumented : the Les Saintes earthquake Mw6.4 and the L’Aquila earthquake Mw6.3. We have developed a blind simulation approach by taking into account an uncertainty on the parameter of stress drop ratio C. This approach allows to generate a set of synthetic accelerograms of a target earthquake varied enough to be representative of a large number of possible source scenario and taking into account in a statistical sense potential directivity effects. This approach is also applied to the simulation of an historical Pyrenean earthquake Mw6.1. In a second part, we use a more complex extended source model, combining kinematic models of fractal composite sources with EGF approach. The potential of the method is tested on an application to L’Aquila earthquake. This has produced very satisfying results on all ground motion parameters analyzed. This simulation method appears to be very promising for the implementation of a blind simulation methodology, even if the main difficulty lies in the need to define the variability of many poorly known input parameters in the simulation of a future earthquake.

L'inversion des marges passives apparaît comme le premier stade vers l'initiation de nouvelles zones de subduction. Cette étape cruciale dans la tectonique des plaques soulève néanmoins encore de nombreuses questions. L'étude des marges actuellement réactivées en compression apparaît ainsi comme essentielle pour mieux comprendre ce processus. Ces marges sont peu nombreuses, situées dans des contextes géodynamiques variés, et les facteurs déterminant leur évolution mal contraints. Située au nord de l'Afrique, la marge algérienne fait partie de ces rares exemples potentiels à travers le monde. L'évolution de cette marge formée au Miocène en contexte d'arrière-arc s'intègre dans le puzzle complexe de l'histoire de la Méditerranée occidentale. Elle est depuis quelques millions d'années réactivée en compression dans le cadre de la convergence lente entre les plaques européenne et africaine, générant un potentiel sismogène fort au nord de l'Algérie. La relative jeunesse du bassin algérien, la charge sédimentaire, les forces aux limites compressives, constituent des conditions favorables à la formation d'une future subduction. A la suite des travaux menés depuis une dizaine d'années, les principales lacunes de connaissances identifiées portent sur (1) la structuration profonde du bassin algérien et de sa marge sud (type de marge, nature du socle,dimension et nature de la transition océan-continent, style et distribution de la déformation compressive), et (2) l'histoire de l'évolution cinématique et géodynamique du bassin, ce qui limite à l'heure actuelle une analyse approfondie des modalités d'inversion de cette marge. L'étude menée se focalise sur la marge centre-algérienne, dans le secteur de Tipaza (à l'ouest d'Alger), un endroit clé pour la compréhension des mécanismes d'ouverture du bassin algérien. Le traitement et l'analyse de nouvelles données de sismique profonde grand-angle et multitraces acquises dans le cadre du projet franco-algérien SPIRAL (Sismique Profonde et Investigations Régionales en Algérie, 2009) ont notamment permis de déterminer la structure crustale du bassin algérien et de sa marge sud, ainsi que la structuration pseudo-3D d'une structure spécifique au secteur d'étude constituée par le haut topographique sous-marin de Khayr-al-Din. L'analyse de la structure profonde de la marge indique un certain nombre de structures héritées de son évolution complexe : (1) une croûte de nature continentale de plus de 15 km d'épaisseur sur le haut de marge (banc de Khayr-al-Din), (2)une croûte fine de nature océanique de 5-6 km d'épaisseur dans le bassin incluant des vitesses légèrement élevées à sa base (7,2 km/s - 7,3 km/s), (3) des similitudes avec des marges formées dans des contextes de déformation transformante, (4) un approfondissement progressif de l'ensemble de la pile sédimentaire et l'épaississement des sédiments Plio-Quaternaire, depuis le bassin profond distal vers le pied de marge,coïncidant avec (5) une flexuration à grande longueur d'onde du socle. Les résultats obtenus apportent de nouvelles contraintes sur (1) la géométrie et la nature de la marge et du bassin, (2) l'évolution de la marge,suggérant une histoire multiphasée comprenant un stade de rifting et/ou d'accrétion océanique, suivi d'un épisode de déformation coulissante tardive liée à la migration du bloc Alboran vers l'ouest, et d'une reprise en compression distribuée du bassin profond au haut de la marge au Plio-Quaternaire; (3) les modalités de réactivation qui se traduisent par des chevauchements aveugles néoformés à pendages sud, notamment au pied du banc de Khayr-al-Din, suggérant un soulèvement du banc de 0,2 mm/an à 0,75 mm/an au Plio-Quaternaire et un début d'écaillage crustal. / The inversion of passive margins appears to be one of the first steps towards the initiation of new subduction zones. This crucial step in plate tectonics nevertheless still raises many questions. The study of margins currently reactivated by compressional tectonics is thus essential to better understand this process. These margins are uncommon, located in different geodynamic settings, and the factors determining their evolution are poorly constrained. The Algerian margin, located in North Africa, is one of handful of modern examples worldwide. The evolution of this margin, rifted during the Miocene, in a back-arc setting, is part ofthe complex puzzle of the western Mediterranean. Since a few million years, the margin has suffered inversion and compression in the framework of slow on going convergence between the European and African plates. This convergence generates moderate to strong earthquakes in North Algeria. The relatively young age of the Algerian basin, the large sediment load, and the compressive forces, constitute favorable conditions to the formation of a future subduction zone. Studies from the past ten years indicate, that themain unresolved questions are related to (1) the deep structure of the Algerian basin and its southern margin (the type of margin, the nature of the basement, the dimension and nature of the ocean-continent transition, the style and the distribution of the compressional deformation), and (2) the history of the kinematic and geodynamic evolution of the basin. All of these unknowns have prevented a complete and thorough analysis of modalities of the Algerian margin inversion. This study focuses on the Central Algerian margin, in the area of Tipaza (West of Algiers), a key region to understand the mechanism of the opening of the Algerian basin. Processing and analysis of a deep wide-angle and multichannel seismic new data set acquired in the context of the French-Algerian project SPIRAL (Sismique profonde et Investigation Régionales en Algérie, 2009)have enabled us to determine the crustal structure of the Algerian basin and its southern continental margin,as well as the pseudo-3D structure of a specific feature in the study area: the submarine topographic highformed by the Khayr-al-Din bank. The analysis of the deep structure of the margin reveals features inherited from its complex evolution: (1) a crust of continental nature of more than 15 km thick at the upper margin(Khayr-al-Din Bank), (2) a thin crust of oceanic nature, 5-6 thick in the deep basin, including slightly high velocities at its base (7.2 km/s - 7.3 km/s), (3) similarities with margins formed in context of transform deformation, (4) a progressive deepening of the whole sedimentary cover and the thickening of the Plio-Quaternary sediments, from the distal deep basin towards the margin foot, coeval with (5) a long wavelengthflexuration of the basement in the basin. Results from this study provide new constraints on (1) the geometryand nature of the margin and the basin, (2) the evolution of the margin, suggesting a multiphased history including a stage of rifting and/or oceanic spreading, a transcurrent episode due to the westward migration of the Alboran block, and a diffuse Plio-Quaternary compressional reactivation distributed from the deep basinto the upper margin; (3) the mechanisms of the reactivation marked by newly formed south-dipping blind-thrusts, especially at the foot of the Khayr-al-Din bank, and suggesting a Plio-Quaternary uplift of the bankof 0.2 mm/y to 0.75 mm/y and the early stages of imbricate thrusting of crustal scales.

La tomographie des temps de première arrivée cherche à estimer un modèle de vitesse de propagation des ondes sismiques à partir des temps de première arrivée pointés sur les sismogrammes. Le modèle de vitesse obtenu peut alors permettre une interprétation structurale du milieu ou bien servir de modèle initial pour d'autres traitements de l'imagerie sismique. Les domaines d'application de cette méthode s'étendent, à des échelles différentes, de la géotechnique à la sismologie en passant par la géophysique pétrolière. Le savoir-faire du géophysicien joue un rôle important dans la difficile résolution du problème tomographique non-linéaire et mal posé. De nombreuses recherches ont entrepris de faciliter et d'améliorer cette résolution par des approches mathématique ou physique. Dans le cadre de ce travail, nous souhaitons développer une approche pragmatique, c'est-à-dire que nous considérons que le problème tomographique doit être résolu par un algorithme interactif dont les paramètres de réglage sont clairement définis. L'aspect interactif de l'algorithme facilite l'acquisition du savoir-faire tomographique car il permet de réaliser, dans un temps raisonnable, de nombreuses simulations pour des paramétrisations différentes. Le but poursuivi dans cette thèse est de définir, pour le cas spécifique de la tomographie des temps de première arrivée, un algorithme qui réponde au mieux à ces critères. Les algorithmes de tomographie des temps de première arrivée classiquement mis en oeuvre aujourd'hui ne répondent pas à nos critères d'une approche pragmatique. En effet, leur implémentation ne permet pas d'exploiter l'architecture parallèle des supercalculateurs actuels pour réduire les temps de calcul. De plus, leur mise en oeuvre nécessite une paramétrisation rendue complexe du fait de la résolution du système linéaire tomographique. Toutes ces limitations pratiques sont liées à la formulation même de l'algorithme à partir de la méthode de Gauss- Newton. Cette thèse repose sur l'idée de formuler la résolution du problème tomographique à partir de la méthode de plus grande descente pour s'affranchir de ces limitations. L'étape clé de cette formulation réside dans le calcul du gradient de la fonction coût par rapport aux paramètres du modèle. Nous utilisons la méthode de l'état adjoint et une méthode définie à partir d'un tracé de rais a posteriori pour calculer ce gradient. Ces deux méthodes se distinguent par leur formulation, respectivement non-linéaire et linéarisée, et par leur mise en oeuvre pratique. Nous définissons ensuite clairement la paramétrisation du nouvel algorithme de tomographie et validons sur un supercalculateur ses propriétés pratiques : une parallélisation directe et efficace, une occupation mémoire indépendante du nombre de données observées et une mise en œuvre simple. Finalement, nous présentons des résultats de tomographie pour des acquisitions de type sismique réfraction, 2-D et 3-D, synthétiques et réelles, marines et terrestres, qui valident le bon comportement de l'algorithme, en termes de résultats obtenus et de stabilité. La réalisation d'un grand nombre de simulations a été rendue possible par la rapidité d'exécution de l'algorithme, de l'ordre de quelques minutes en 2-D.

Nous proposons dans cette étude quelques méthodes stochastiques pour le calcul des ponts aux séismes. Après avoir rappelé les représentations du mouvement sismique par spectres de réponse, densités spectrales de puissance et signaux temporels, nous introduisons des modèles de variabilité spatiale d'ondes sismiques. Nous en étudions l'effet sur la réponse linéaire des ponts grâce à une méthode de spectre de réponse pour structures multi-supportées. Nous proposons une méthode simplifiée permettant un calcul rapide de la sensibilité des ouvrages à de faibles variations spatiales des mouvements du sol et nous donnons un algorithme efficace pour le calcul de la réponse la plus défavorable lorsque les variations spatiales sont grandes et incertaines. Pour tenir compte du potentiel de ductilité des ponts, nous proposons une méthode de spectre de réponse tirant partie de la localisation des non-linéarités matérielles dans la structure (rotules plastiques, appareils d'appui). Les éléments non-linéaires sont remplacés par des éléments linéaires équivalents. Une technique de synthèse modale adaptée à la localisation des non-linéarités permet un calcul rapide de la réponse. Cette approche intègre le caractère aléatoire des sollicitations sans avoir recours à de coûteuses simulations de Monte-Carlo. Elle donne de bons résultats en termes de demande de ductilité et d'énergie dissipée par cycles d'hystérésis dans les éléments non-linéaires. Enfin, nous évoquons une technique de calcul de fiabilité de structures non-linéaires sous sollicitations stochastiques.

Depuis une vingtaine d'années, un effort considérable a été effectué dans le développement d'outils dédiés à la caractérisation rapide de la source sismique. Jusqu'à récemment, plusieurs heures étaient encore nécessaires pour obtenir une information fiable sur la source des grands tremblements de terre. C'est dans ce cadre que s'inscrit ce travail, dont l'objectif principal est le développement et la mise en oeuvre d'une méthode de détermination rapide du tenseur moment sismique centroid (CMT) basée sur l'utilisation de la phase W pour les grands séismes. La phase W est une phase sismique correspondant essentiellement à la superposition des modes normaux supérieurs à très longue période (entre 100 sec et 1000 sec) arrivant entre l'onde P et les ondes de surface. Nous montrons que la phase W a l'avantage de ne pas être affectée par les hétérogénéités latérales superficielles puisque la majeure partie de son énergie se propage en profondeur. De par son caractère longue période, la phase W est particulièrement bien adapté à la caractérisation de la source pour les évènements de Mw>7.5. Pour les grands séismes en effet, la complexité de la source apparaît de plus en plus évidente aboutissant inévitablement à une variabilité des scénarios de sources estimés à courte période dans des bandes fréquentielles étroites. Les "séismes tsunami" et les séismes "outer-rise" sont des exemples extrêmes illustrant bien ce problème. Nous démontrons la robustesse de l'algorithme phase W développé pendant ce travail pour effectuer la caractérisation rapide et systématique de la source sismique. Il est aujourd'hui implémenté en temps réel dans plusieurs centres d'alerte à l'échelle globale et à l'échelle régionale. A des distances télésismiques, la méthode permet l'obtention du CMT dans la première demi-heure après le temps origine. A des échelles régionales, les résultats sont disponibles beaucoup plus rapidement, entre 6 et 12 minutes après le temps origine. Les résultats obtenus en temps réel à l'United States Geological Survey (USGS), au Pacific Tsunami Warning Center (PTWC), à l'Institut de Physique du Globe de Strasbourg (IPGS) et à l'échelle régionale au Mexique montrent clairement la robustesse et la fiabilité des solutions CMT obtenues en utilisant la phase W. En suivant une formulation bayesienne, nous proposons une analyse d'erreur formelle lors de l'inversion de la source sismique à longue période. La prise en compte plus réaliste de l'erreur sur les données permet non-seulement une estimation fiable de l'incertitude sur le modèle de source mais aussi l'amélioration de la solution elle même. Nous développons également une approche pour estimer l'erreursur la profondeur du centroid. Ce paramètre est important étant donné son influence sur le moment sismique scalaire et le pendage du plan de faille déterminés lors des inversions CMT.

Le problème de la dimension stochastique élevée est récurrent dans les analyses probabilistes des structures. Il correspond à l'augmentation exponentielle du nombre d'évaluations du modèle mécanique lorsque le nombre de paramètres incertains est élevé. Afin de pallier cette difficulté, nous avons proposé dans cette thèse, une approche à deux étapes. La première consiste à déterminer la dimension stochastique efficace, en se basant sur une hiérarchisation des paramètres incertains en utilisant les méthodes de criblage. Une fois les paramètres prépondérants sur la variabilité de la réponse du modèle identifiés, ils sont modélisés par des variables aléatoires et le reste des paramètres est fixé à leurs valeurs moyennes respectives, dans le calcul stochastique proprement dit. Cette tâche fut la deuxième étape de l'approche proposée, dans laquelle la méthode de décomposition de la dimension est utilisée pour caractériser l'aléa de la réponse du modèle, par l'estimation des moments statistiques et la construction de la densité de probabilité. Cette approche permet d'économiser jusqu'à 90% du temps de calcul demandé par les méthodes de calcul stochastique classiques. Elle est ensuite utilisée dans l'évaluation de l'intégrité d'une toiture à ossature bois d'une habitation individuelle installée sur un site d'aléa sismique fort. Dans ce contexte, l'analyse du comportement de la structure est basée sur un modèle éléments finis, dans lequel les assemblages en bois sont modélisés par une loi anisotrope avec hystérésis et l'action sismique est représentée par huit accélérogrammes naturels fournis par le BRGM. Ces accélérogrammes permettent de représenter différents types de sols selon en se référant à la classification de l'Eurocode 8. La défaillance de la toiture est définie par l'atteinte de l'endommagement, enregistré dans les assemblages situés sur les éléments de contreventement et les éléments d'anti-flambement, d'un niveau critique fixé à l'aide des résultats des essais. Des analyses déterministes du modèle éléments finis ont montré que la toiture résiste à l'aléa sismique de la ville du Moule en Guadeloupe. Les analyses probabilistes ont montré que parmi les 134 variables aléatoires représentant l'aléa dans le comportement non linéaire des assemblages, 15 seulement contribuent effectivement à la variabilité de la réponse mécanique ce qui a permis de réduire la dimension stochastique dans le calcul des moments statistiques. En s'appuyant sur les estimations de la moyenne et de l'écart-type on a montré que la variabilité de l'endommagement dans les assemblages situés dans les éléments de contreventement est plus importante que celle de l'endommagement sur les assemblages situés sur les éléments d'anti-flambement. De plus, elle est plus significative pour les signaux les plus nocifs sur la structure.

Les Moments ont été largement utilisés dans la reconnaissance de formes et dans le traitement d'image. Dans cette thèse, nous concentrons notre attention sur les 3D moments orthogonaux de Gauss-Hermite, les moments invariants 2D et 3D de Gauss-Hermite, l'algorithme rapide de l'attribut de cohérence et les applications de l'interprétation sismique en utilisant la méthode des moments.Nous étudions les méthodes de suivi automatique d'horizon sismique à partir de moments de Gauss-Hermite en cas de 1D et de 3D. Nous introduisons une approche basée sur une étude multi-échelle des moments invariants. Les résultats expérimentaux montrent que la méthode des moments 3D de Gauss-Hermite est plus performante que les autres algorithmes populaires.Nous avons également abordé l'analyse des faciès sismiques basée sur les caractéristiques du vecteur à partir des moments 3D de Gauss -Hermite, et la méthode de Cartes Auto-organisatrices avec techniques de visualisation de données. L'excellent résultat de l'analyse des faciès montre que l'environnement intégré donne une meilleure performance dans l'interprétation de la structure des clusters.Enfin, nous introduisons le traitement parallèle et la visualisation de volume. En profitant des nouvelles performances par les technologies multi-threading et multi-cœurs dans le traitement et l'interprétation de données sismiques, nous calculons efficacement des attributs sismiques et nous suivons l'horizon. Nous discutons également l'algorithme de rendu de volume basé sur le moteur Open-Scene-Graph qui permet de mieux comprendre la structure de données sismiques.

La marge NE Américaine est une des marges les mieux étudiées au monde, elle a fait l’objet de plusieurs études géophysiques. En comparaison, la marge africaine reste peu étudiée car uniquement deux campagnes océanographiques y ont été menées : la campagne Sismar (2001) au large de la Meseta et la campagne Dakhla (2002) au large du Sahara. La structure profonde de la marge canadienne est connue grâce aux profils de sismique grand-angle SMART-1, 2 et 3. Le premier objectif du projet MIRROR était d’acquérir des profils combinant sismique grand-angle et sismique réflexion sur un segment homologue au profil SMART-1. La comparaison entre les segments homologues de ces deux marges ayant pour but de mieux comprendre le mode d’ouverture de l'océan Atlantique Central. Une comparaison entre les modèles Sismar, Dakhla et Mirror montre que la croûte continentale est plus épaisse au nord et s'amincit vers le sud. La largeur de la zone de transition est plus étroite au sud et les profils Sismar sont localisés sur un bassin sédimentaire posé sur une croûte continentale très amincie. La comparaison avec la marge homologue montre que l'épaisseur, la structure de la croûte continentale et la zone d'amincissement sont très semblables. Par contre, il existe une zone de manteau exhumé et serpentinisé sur le profil Canadien qui n'a pas d’homologue sur la marge africaine. De plus, l'épaisseur de la croûte océanique est différente avec 8 km sur la côte africaine et seulement 3-4 km sur la marge canadienne. Plusieurs hypothèses ont été proposées pour expliquer cette différence (a) une différence d’âge entre les deux croûtes (b) un épaississement lié au passage du point chaud des Canaries (c) une accrétion asymétrique. / The NE American margin represents one of the best studied margins in the world, it was the subject of several scientific programs. In comparison, the conjugate NW African margin remains fairly unknown, only two deep seismic cruises were acquired: the SISMAR cruise (2001) offshore the Meseta and the DAKHLA cruise (2002) offshore the Sahara. The deep structure of the Canadian margin is known due to the SMART wide-angle seismic profiles 1, 2 and 3. The first objective of the MIRROR project was to acquire combined wide-angle and deep reflection seismic data offshore a segment conjugate to the SMART-1 profile. The comparison between the homologous segments of these two margins aimed to better understand the opening mechanism of the Central Atlantic Ocean. A comparison between Sismar, Dakhla and Mirror models shows that the continental crust is thicker in the north and thins toward the south. The width of the transition zone is narrower south and Sismar profiles are located on a sedimentary basin placed on a very thinned continental crust. Comparing the Mirror profile with that of the Canadian conjugate margin (Smart 1) shows that the thickness, the structure of the continental crust and the thinning is very similar. However, zones of exhumed and serpentinized mantle were imaged along the Canadian profile that have no conjugate on the African margin. Moreover, the thickness of the oceanic crust is variable with 8 km on the African side and only 3-4 km on the Canadian margin. Several hypotheses have been proposed to explain this difference (a) an age difference between the two types of crust (b) thickening associated with the passage of the Canary hotspot (c) an asymmetric accretion or (d) an accretion at slow to ultra-slow speading centers.