La théorie de la diffraction prévoit que les temps de trajet des ondes avec une fréquence dominante finie sont affectés seulement par des hétérogénéités dans un volume autour du rayon obtenu par la théorie asymptotique fréquence infinie. Ce volume s'appelle la zone de Fresnel. A cause de la diffraction qui assure la continuité du front d'onde au travers des hétérogénéités, les hétérogénéités de taille inférieure à celle de la zone de Fresnel n'ont qu'un effet réduit sur le temps de trajet de l'onde arrivant à une station. Dahlen et al (2000) montrent qu'il est possible decorriger des effets de cette diffraction par une méthode de tomographie volumique qui introduit de nouveau la notion de contenu fréquentiel dans l'estimation des temps de trajet et donc dans la reconstruction du milieu. Dans cette thèse, j'étudie l'influence de cette technique d'imagerie à fréquence finie en tomographie globale sismologique des temps d'arrivée. Pour réduire les effets du à la distribution des rayons dans la Terre,distribution qui n'est pas uniformej j'introduis un maillage non-structuré pour représenter les vitesses dans la Terre. Pour définir l'influence des phénomènes de diffraction dans les résultats tomographiques, j'utilise un nombre limité de données à longue période de très haute qualité obtenues par Guy Masters (Scripps). Je fais l'inversion en utilisant l'approche classique de tomographie des temps se fondant sur la théorie des rayons et l'approche nouvelle de tomographie des temps en se fondant sur l'approche à fréquence finie. La comparaison des deux tomographies montre que les anomalies de vitesse sont sous-estimées par la théorie des rayons d'un facteur 30-60% en fonction de la profondeur. C'est bien le signe que la théorie des rayons est une approximation trop grossière pour les données à basse fréquence que l'on utilise. Pour augmenter la résolution de mon modèle tomographique, j'inverse conjointement les données à longue période déjà mentionnées avec des données haute fréquence retraitées par Bob Engdahl (Boulder). L'analyse de la compatibilité des deux bases de données est nécessaire avant de procéder à l'inversion conjointe. L'étude des temps eux-même (en fait, des écarts dans un modèle de référence) et les inversions des deux jeux de données independament montrent un très bon niveau de compatibilité. Tous les modèles obtenus montrent une très forte corrélation entre les anomalies de basse vitesse ct la position d'un très grand nombre de points chauds déjà connus. Pour estimer la validité de ces anomalies, des tests de résolution ont été réalisés de façon à identifier la profondeur des panaches. A la suite de ces tests, on peut conclure positivement que des panaches profonds se trouvent au-dessous des structures suivantes: Ascension, Açores, Canaries, Easter, Hawai, Samoa et Tahiti. Pour d'autres, l'origine profonde n'est pas confirmée et une origine plus superficielle est probable autour de 660 km comme la zone islandaise. D'autres n'ont jamais été vus avant comme au Seychelles et dans l'Ocean Indien. Enfin, dans la zone de la Mer de Coral et au Sud de Java, la tomographie a détecté des panaches montant de la frontière entre le noyau ct le manteau mais qui n'arrivent pas encore à la surface.

Connaître mieux l'intérieur de la Terre nécessite d'améliorer les méthodes d'investigation géophysique. Cette thèse s'est interessée à une classe importante de ces méthodes : la tomographie sismique. La tomographie sismique classique utilise des temps d'arrivée d'ondes sismiques, dont la précision dépend du pas d'échantillonnage et du pointé de l'onde. Lorsque les trains d'onde sont similaires, c'est-à-dire lorsque les sources et/ ou les récepteurs sont proches, on peut calculer, avec une précision inférieure au pas d'échantillonnage, le décalage temporel entre les trains d'ondes. Cette thèse a eu pour objectif l'utilisation de telles données en tomographie sismique. Utiliser pleinement la précision de ces données a nécessité une revisite complète des algorithmes de tomographie, tant du point de vue du problème direct que du point de vue du problème inverse. Une solution robuste et précise a été trouvée au problème direct (calcul du temps de propagation des ondes dans un milieu hétérogène). Une attention particulière a été portée au problème inverse : une approche Tarantola-Valette a été utilisée, de façon à trouver le modèle de norme minimale qui ajuste les données. La minimisation de cette norme est obtenue par l'introduction d 'une fonction de corrélation des paramètres du modèle, dont la longueur et l'amplitude sont ajustées en suivant une procédure d'optimisation. De cette façon, toute fluctuation du modèle est contrainte par les données. L'algorithme mis au point a été utilisé avec des données provenant de séismes du volcan Kilauea (Hawaii) enregistrés par le réseau local USGS-HVO (50 stations). Il a permis de retrouver le modèle de vitesse à deux échelles : (1) une échelle régionale, où l'on retrouve le système magmatique profond, sous les calderas et les rifts du Mauna Loa et du Kilauea; le modèle trouvé permet d'expliquer les déformations mises en évidence dans la croûte océanique, la distribution confinée de la sismicité et les déformations intenses du flanc Sud du Kilauea; (2) une échelle locale, où l'on met en évidence un conduit magmatique sub-vertical sous la caldera et le rift Est du Kilauea, indétectable avec d'autres données. Cette méthode, précise, permet de choisir un modèle unique et stable à l'aide un critère d'optimisation objectif.

Le Golfe de Corinthe (Grèce) est un rift continental caractérisé par un taux élevé d'extension et une microsismicité intense. Actuellement, la connaissance de sa structure crustale représente une étape nécessaire pour la construction de modèles géodynamiques réalistes. L'ojectif de ce travail de thèse est la définition d'une image sismique haute résolution de la région occidentale du Golfe de Corinthe à partir des enregistrements de la sismicité locale. La construction de cette image est effectuée par le biais d'une étude tomographique des temps de première arrivée et d'une migration cinématique en profondeur des ondes converties. La tomographie des premières arrivées nous a fourni des modèles de vitesse tridimensionnels des ondes P et S ainsi qu' une localisation précise de la microsismicité dans les 10 premiers kilomètres de la croûte. Les modèles de vitesse ainsi obtenus constituent les modèles de référence utilisés dans la migration en profondeur des ondes converties. A cet effet , nous avons adapté les techniques propres à la sismique d'exploration afin d'analyser les enregistrements de tremblements de terre avec une configuration d'acquisition irrégulière. Les résultats révèlent une image tridimensionnelle complexe de la zone sismogène du rift (7-10 km de profondeur) dans laquelle deux structures majeures sont identifiées par des importantes énergies transmises PS. Du point de vueméthodologique, ce travail montre la possibilité d 'exploiter les ondes converties à partir des enregistrements de séismes locaux pour définir une image haute fréquence du milieu de propagation à l'échelle de la croute terrestre.

Sous le Chili central et l'ouest de l'Argentine (29°-35°S), la plaque océanique Nazca, en subduction sous la plaque continentale Amérique du Sud, change radicalement de géométrie : inclinée à 30°, puis horizontale, engendrée par la subduction de la chaine de volcans de Juan Fernandez. Le but de mon étude est d'évaluer, la variation de nature et de propriétés physiques de la lithosphère chevauchante entre ces deux régions afin de mieux comprendre (1) sa structure profonde et (2) les liens entre les déformations observées en surface et en profondeur. Pour répondre à cette thématique, j'utilise une approche originale couplant la sismologie, la thermométrie, et la pétrologie. Je montre ainsi des images 3-D de tomographie sismique les plus complètes de cette région par rapport aux études précédentes, qui intègrent (1) de nombreuses données sismiques provenant de plusieurs catalogues, (2) un réseau de stations sismiques plus dense permettant de mieux imager la zone de subduction. J'apporte la preuve que la plaque en subduction se déshydrate dans deux régions distinctes : (1) le coin mantellique, et (2) le long de la ride subduite avant que celle-ci ne replonge plus profondément dans le manteau. La croûte continentale au-dessus du flat slab possède des propriétés sismiques très hétérogènes en relation avec des structures de déformation profondes et des domaines géologiques spécifiques. La croûte chevauchante d'avant-arc, au-dessus du flat slab, est décrite par des propriétés sismiques inhabituelles, liées à la géométrie particulière du slab en profondeur, et/ou liées aux effets du séisme de 1997 de Punitaqui (Mw 7.1). Mes résultats, confirmant les études antérieures, montrent que : - le bloc Cuyania situé plus à l'est, dans la zone d'arrière-arc est plus mafique et contient une croûte inférieure éclogitisée ; quant à, la croûte continentale inférieure sous l'arc Andin, est épaisse et non-éclogitisée, décrivant surement le bloc felsique de Chilenia.

La zone de subduction Hellénique, en Méditerranée orientale, est caractérisée par le taux de sismicité le plus important d’Europe. Des séismes de forte magnitude (Mw 7,5-8) ont eu lieu le long du segment Sud-Ouest de la zone de subduction Hellénique, au large du Péloponnèse, au cours du 19ème et 20ème siècle. Ce segment de 400 km de long a également été le lieu de nucléation du plus important séisme d’Europe, en 365 ap J.C, avec une magnitude supérieure à 8, ayant entraîné un tsunami dévastateur. Deux principaux modèles scientifiques s’opposent sur la question du couplage sismique de l’interface de subduction, allant d’un couplage sismique total au niveau de l’interface, à l’hypothèse opposée d’un couplage quasi inexistant. Cependant, ces modèles opposés considèrent des géométries approximatives et parfois extrêmes, fautes de contraintes disponibles sur la structure et la géométrie de l’interplaque sous l’avant-arc dans cette zone. La localisation de la faille responsable du séisme de 365 ap J.C est également débattue, en l’absence de données géophysiques permettant d’identifier les interfaces potentiellement responsables de cet événement dévastateur. La faille de méga-chevauchement et le domaine avant-arc du segment Sud-Ouest de l’arc Hellénique ont été l’objet d’étude de la campagne océanographique Ulysse en Novembre 2012 afin de déterminer la géométrie des structures et unités majeures dans cette portion de la zone de subduction, mais également d’apporter un éclairage sur la tectonique récente qui affecte cette zone / The Hellenic subduction zone, in the eastern part of the Mediterranean sea, is characterized by the highest rate of current seismicity in Europe. In the southwestern segment, several earthquakes of large magnitude (Mw 7,5-8) occured a the turn of the 19th to 20th century. This segment of 400 km long, has also been the nucleation site of the largest historical earthquake in Europe, named the 365 AD earthquake, with a magnitude of Mw 8. This event generates a devastating tsunami, which spread along the Adriactic Sea and in the Nile Delta region. Two main models differ about the interplate seismic coupling question in this region, from a total seismic coupling at the interplate, at the opposite assumption of a very weak seismic coupling. However, these opposing models consider an approximate geometry, mostly because of the lack of information available on the geometry and the localization of the interplate in this region of the forearc domain. The localization of the fault responsible of the 365 AD event is also debated, because, there is no available data who provides imagery of the interfaces potentially responsible of this devastating earthquake. The megathrust fault and the forearc domain of the southwestern segment of the Hellenic subduction zone has been the target of the Ulysse marine survey in November 2012. The aim of this survey was to provide information of the structural geometry of the main units in this part of the subduction zone, and to bring information on the recent tectonic activity in this region

Une nouvelle méthode de tomographie est décrite (PRIAM~2D). C'est une méthode géométrique itérative. Au cours de l'inversion, des zones d'influence sont définies pour chacune des cases prises en compte. La taille de ces zones est choisie a priori (à partir de considérations physiques). Leur existence transforme la méthode géométrique en une technique d'inversion sans blocs. Le calcul au sein de chacune des zones comporte des pondérations visant à discriminer les rais entre eux. Celles- ci permettent également de discriminer certaines parties de l'image les autorisant ainsi à varier plus librement que d'autres. Les rais utilisés sont analytiques. Des évaluations de la résolution et des erreurs sont introduites. Ces techniques permettent d'obtenir des résultats comparables à ceux fournis par des méthodes matricielles tout en conservant des temps de calcul nettement inférieurs. Dans le domaine du Génie Civil plusieurs applications sont envisagées. Des sols, ainsi que certaines parties d'ouvrages d'Art, peuvent être auscultés. Il est possible d'obtenir des cartes de la vitesse des ondes de compression ou du facteur de qualité, à partir de mesures de durées de propagation, de temps de montée ou d'amplitudes. En outre, Une chaine de mesure spéciale a été mise au point de manière à minimiser le temps de la mesure.

Sous le Chili central et l’ouest de l'Argentine (29°-35°S), la plaque océanique Nazca, en subduction sous la plaque continentale Amérique du Sud, change radicalement de géométrie : inclinée à 30°, puis horizontale, engendrée par la subduction de la chaine de volcans de Juan Fernandez. Le but de mon étude est d'évaluer, la variation de nature et de propriétés physiques de la lithosphère chevauchante entre ces deux régions afin de mieux comprendre (1) sa structure profonde et (2) les liens entre les déformations observées en surface et en profondeur. Pour répondre à cette thématique, j’utilise une approche originale couplant la sismologie, la thermométrie, et la pétrologie. Je montre ainsi des images 3-D de tomographie sismique les plus complètes de cette région par rapport aux études précédentes, qui intègrent (1) de nombreuses données sismiques provenant de plusieurs catalogues, (2) un réseau de stations sismiques plus dense permettant de mieux imager la zone de subduction. J’apporte la preuve que la plaque en subduction se déshydrate dans deux régions distinctes : (1) le coin mantellique, et (2) le long de la ride subduite avant que celle-ci ne replonge plus profondément dans le manteau. La croûte continentale au-dessus du flat slab possède des propriétés sismiques très hétérogènes en relation avec des structures de déformation profondes et des domaines géologiques spécifiques. La croûte chevauchante d’avant-arc, au-dessus du flat slab, est décrite par des propriétés sismiques inhabituelles, liées à la géométrie particulière du slab en profondeur, et/ou liées aux effets du séisme de 1997 de Punitaqui (Mw 7.1). Mes résultats, confirmant les études antérieures, montrent que : - le bloc Cuyania situé plus à l’est, dans la zone d’arrière-arc est plus mafique et contient une croûte inférieure éclogitisée ; quant à, la croûte continentale inférieure sous l’arc Andin, est épaisse et non-éclogitisée, décrivant surement le bloc felsique de Chilenia. / Beneath central Chile and western Argentina, the oceanic Nazca slab drastically changes geometry from horizontal to dipping at an angle of 30°, and correlates with the subduction of the Juan Fernandez seamount ridge. The aim of our study is to assess, using a thermo-petrological-seismological approach, the differences of the overriding lithosphere between these two regions, in order to better understand the deep structure of the continental lithosphere above the flat slab, and the links between the deformations at the surface and at depth. We show the most complete regional 3-D seismic tomography images of this region, whereby, in comparison to previous studies, we use (1) a much larger seismic dataset compiled from several short-term seismic catalogs, (2) a much denser seismic station network which enables us to resolve better the subduction zone. We show significant seismic differences between the flat and normal subduction zones. As expected, the flat slab region is impacted by colder temperatures, and therefore by faster seismic velocities and more intense seismic activity, compared to the normal slab region. We show evidence that the flat slab dehydrates within the mantle wedge, but also along the subducting ridge prior to re-subducting. The forearc crust above the flat slab is described by unusual seismic properties, correlated to the slab geometry at depth, and/or, to the aftershock effects of the 1997 Mw 7.1 Punitaqui earthquake which occurred two years before the recording of our events. The continental crust above the flat slab has very heterogeneous seismic properties which correlate with important deformation structures and geological terranes at the surface. We confirm previous studies that have shown that the thick lower crust of the present day Andean arc is non-eclogitized and maybe representing the felsic Chilenia terrane, whereas to the east, the Cuyania terrane in the backarc is more mafic and contains an eclogitized lower crust.

La lithosphère continentale est physiquement et chimiquement segmentée. La cartographie des isotopes radiogéniques de roches plutoniques acides, représentatives de la croûte continentale, et de laves basiques, représentatives du manteau, possède des similarités avec la cartographie sismique de la lithosphère sous-jacente. Ces similitudes permettent d’interpréter les observations sismiques en étudiant leurs caractéristiques chimiques et leur âge. Les isotopes du plomb permettent de dater et d’identifier l’empilement de segments crustaux qui forment la croûte. L’écart des âges modèles du plomb avec d’autres systèmes identifie le recyclage crustal et le réchauffement de la croûte au dessus de la température du système plomb-plomb. Le système plomb-plomb donne également accès au sous-étudié rapport Th/U qui contraint la profondeur de la source des roches continentales. Certains échantillons de l’ouest U.S.A. proviennent de la croute inférieure, et se sont formés par l’extension crustale ou par un flux de matériel au sein de la croûte. Les isotopes du néodyme et de l’hafnium marquent la fusion du manteau lithosphérique enrichis sous le Colorado Plateau, une région où est observée le détachement du manteau lithosphérique sub-continental. Ce manteau fond par décompression adiabatique, par extension localisée ou remontée asthénosphérique engendrée par la convection locale. Au final, l’association des systèmes isotopiques du plomb, du néodyme, et de l’hafnium avec la sismologie est une approche puissante pour étudier la formation et la déformation de la lithosphère continentale. / Continental lithosphere is physically and chemically segmented. The mapping at a continent size scale of radiogenic isotopes from plutonic acid rocks, sampling the continental crust, and from mafic lavas, sampling the mantle, has similarities with the seismic mapping of the underlying lithosphere. These similarities allow to interpret the seismic observations by studying their chemical characteristics and age. Lead isotopes are used to date and identify the stacking of crustal segments that form the crust. The deviation of Lead model ages with other system is used to identify crustal recycling and the heating of the crust above the Lead-Lead system closing temperature. The Lead-Lead system also give access to the under-studied Th/U ratio that constrains the depth of the continental rock sources. Somes samples from the Western U.S. are coming from the lower crust, formed by crustal extension or crustal flowing within the crust. Neodymium and Hafnium isotopes identify the fusion of an enriched lithospheric mantle under the Colorado Plateau, a place where lithospheric delamination has been observed. This mantle melts by adiabatic decompression due to localized expansion or asthenospheric upwelling caused by secondary convection. At the end, the merging of Lead, Neodymium, and Hafnium isotopic systems with seismology is a powerful tool to study the formation and deformation of the continental lithosphere.

Ce projet de recherche répond à l'objectif d'améliorer notre connaissance des processus élémentaires contrôlant le comportement hydrologique post-rupture des glissements de terrain argileux, afin de pouvoir comprendre, reproduire et éventuellement prévoir leur comportement par des modèles hydromécaniques adaptés. Ce travail consiste à améliorer et développer des méthodes de traitement et d'interprétation de paramètres géophysiques en milieu argileux et hétérogène. Les objectifs principaux de ce travail sont : 1. D'améliorer des méthodes géophysiques et géodésiques de mesure de la déformation et de la géométrie interne des glissements. 2. De caractériser les relations possibles entre géométrie et structure interne du glissement. 3. De développer et améliorer des méthodes de détection et d'imagerie des circulations préférentielles de fluide dans un milieu fissuré et déformable. 4. De caractériser temporellement et spatialement les réponses hydrologiques du versant. Un algorithme d'inversion haute-résolution des temps de premières arrivées est développé. Il utilise les volumes de Fresnel pour rétropropager les résidus ainsi que l'hypothèse d'un signal source de bande-passante finie et permet de retrouver une image haute-résolution de la structure du glissement. L'atténuation sismique des ondes-P est également inversée pour retrouver la distribution des zones fissurées en surface. Des méthodes photogrammétriques sont utilisées pour caractériser et quantifier les déformations de surface. Une méthode alternative à la corrélation d'image est développée, fondée sur la détection et le suivi de cibles. Le lien entre les variations de régime de déformation le long du profil, et l'existence de fissures en surface est finalement discuté. Les informations précédentes sont ensuite intégrées dans un modèle géométrique par une méthode de fusion de données. Par la modélisation numérique, nous vérifions si l'état de contrainte calculé peut être en rapport avec la distribution de la fissuration en surface. Une méthode de correction de l'effet géométrique des fissures sur les résistivités apparentes est développée. Différents types d'artefacts d'inversion liés à la présence de fissure sont identifiés et la capacité de la méthode proposée à réduire ces artefacts est testée. Une méthodologie de sélection et de traitement des données ERT time-lapse acquises sur une période d'un an est développée. La réponse électrique courte (de quelques heures à quelques jours) aux évènements pluvieux naturels est ensuite analysée pour proposer des modèles conceptuels de réponse hydrologique à ces évènements pluvieux. / The objective of this research project is the improvement of our knowledge of the elementary processes controlling the post-rupture hydrological behavior of clayey landslides. It should allow understanding, reproducing and possibly forecasting their behaviors by the use of adapted hydro-mechanical models. This work consists in the improvement and the development of geophysical parameters processing and interpretation methodologies in clayey and heterogeneous medium. The main objectives of this work are : 1. The improvement of geophysical and geodesic methods to measure the internal geometry and structure of the landslide. 2. The characterization of possible relationships between landslide internal geometry and structure, presence of discontinuities, slope hydrology and deformations. 3. The development and improvement of methodology for the detection and imaging of preferential flows in a fissured and deformable medium. 4. The spatial and temporal characterization of the hydrological response of the slope. A high resolution inversion algorithm of the first seismic arrivals is developed. It uses the Fresnel volumes for the retropropagation of the residuals and the assumptions of a finite-bandwidth source signal. It permits to recover a high resolution image of the internal structure of the landslide. The seismic P-wave attenuation is also inverted to characterize surface fissured zones. Photogrammetric methods are used to characterize and quantify the surface deformations. A method alternative to image correlation is developed, based on the target detection and tracking. The link between changes in deformation regime along the profile and the existence of surface fissures is finally discussed. The previous data are integrated inside a geometric model trough data fusion. Using numerical modeling, we verify if the state of stress computed is related to the surface fissure density distribution. A methodology of correction of the geometrical effect of fissures on apparent resistivity is developed. Different types of inversion artefacts, linked to the presence of surface fissures, are identified. The ability of the proposed method to limit this artifact is tested. A methodology of selection and processing of the time-lapse ERT data is developed. The short electrical response (from several hours to several days) to natural rainfall events is analyzed and conceptual hydrological responses to these rainfall events are proposed.

Mon travail de thèse se propose d'étudier la structure, les propriétés physiques et les processus géodynamiques de la zone de subduction d'Equateur-Colombie grâce à l'adaptation et le développement d'outils d'imagerie sismique (inversion de formes d'ondes 'alias' tomographie en diffraction) et à leur application aux données de sismique marine multitrace (MCS) et grand-angle OBS (WA) acquises en Equateur-Colombie pendant les campagnes SISTEUR et SALIERI. Ces outils m'ont permis de réaliser une imagerie fine et quantitative à trois niveaux : l'imagerie superficielle (~ 0-3 km), l'imagerie à profondeur intermédiaire (~ 3-10 km) et l'imagerie profonde (~ 10-30 km). Dans le domaine superficiel, j'ai effectué une cartographie fine et quantitative des propriétés physiques des sédiments au voisinage du BSR (Bottom Simulating Reflector), interpreté comme la base de stabilité des hydrates de gaz. Sur le profil SIS-40 situé sur la marge sud de la Colombie, j'ai pu identifier la présence de failles qui perturbent localement le BSR. Les résultats présentés sous la forme d' une série de logs adjacents de l'image migrée en profondeur, montrent que certaines régions du BSR sont caractérisées par une augmentation de la vitesse (1470-1650 m/s), indiquant la présence d'une faible quantité d'hydrates de gaz au dessus du BSR; d'autres zones situées immédiatement sous le BSR sont caractérisées par une diminution de la vitesse (~1200 m/s), liée à la présence de gaz libres piégés sous la couche d'hydrate de gaz. A des profondeurs moyennes j'ai étudié la structure du chenal de subduction (profil SIS-72). Le chenal constitue la limite mécanique entre la plaque chevauchante et la plaque plongeante. Il est délimité à son toit par un fort réflecteur interprété comme le décollement interplaque et à sa base par le toit très réflectif de la croûte océanique en subduction. L'imagerie fine et quantitative des propriétés physiques du décollement interplaque permet de mieux comprendre le rôle de la circulation des fluides et des variations lithologiques et physiques, sur le couplage mécanique inter-plaque. En raison de la sensibilité de la méthode de tomographie en diffraction au macro-modèle de vitesse, un code de correction de ce modèle a été implémenté, afin d'obtenir des images tomographiques fiables (i.e. géométrie et amplitudes correctes). Du fait de la bande passante limitée de la source et de la longueur du dispositif d'acquisition limitée à 4.5 km, les images tomographiques ont une résolution spatiale limitée : l'image tomographique présente un déficit des petits et grands nombre d'onde (fréquences spatiales) limitant ainsi l' interprétation géologique des paramètres physiques cartographiés. Un traitement spécifique basé sur la modélisation des traces sismiques a été implémenté. L'image tomographique, traitée comme une série de traces verticales, constitue la donnée observée. L'espace des modèles est constitué par un ensemble de modèles impulsionnels et unidimensionnels de Terre construits aléatoirement. Ces modèles sont dégradés par convolution avec une estimation de l'ondelette source afin de fournir une représentation synthétique de l'image tomographique « observée ». La minimisation de la fonction coût entre les traces migrées et les traces synthétiques est effectuée dans le cadre d'une inversion globale par recuit simulé (VFSA= « Very Fast Simulated Annealing »). Le modèle moyen issu de cette procédure fournit un modèle 2D fin de vitesse, fonction de la profondeur et comparable à la limite de la résolution théorique de la source. A l'issue de ce traitement, des perturbations de vitesse positives sont mises en évidence au toit de la croûte, et d'autres négatives accompagnent certains segments du niveau du décollement. Ces dernières sont probablement associées à la présence de fluides. Le domaine plus profond a été étudié à partir des données MCS et WA dans le double but (1) d'améliorer la résolution spatiale des images sismiques du Moho et du contact interplaque en relation avec la zone sismogène, et (2) de détecter la présence d'anomalies crustales de vitesse et d'analyser leur relation avec les zones d'aspérité sismologiques. L'utilisation conjointe des données de sismique MCS et WA a été mise en oeuvre pour prolonger vers le bas les images de sismique verticale et tenter ainsi d'établir une relation entre les processus profonds et les manifestations en surface. L'application de la chaîne de traitement au profil SIS-44 a permis d'obtenir un modèle de vitesse bien contraint jusqu'à 25 km de profondeur. Ce modèle met en évidence des réflecteurs profonds (Moho et contact interplaque ) et des réflecteurs plus superficiels (splay fault), dont l'interprétation était initialement incertaine sur les images migrées en temps.