Les séismes de type longue période (LP) sont aujourd’hui enregistrés sur la plupart des volcans dans le monde entier. Malgré cela, le mécanisme à leur source n’est encore que très peu compris. A l’heure actuelle les modèles proposés pour expliquer leur origine sont : 1) la résonance d’une fracture remplie de gaz ou de fluides excités par des instabilités dans l’écoulement des fluides ou par la rupture fragile du magma ; 2) la fracturation lente des sédiments peu consolidés à la surface des volcans, dans des conditions de transition entre le ductile et le fragile. L’outil le plus utilisé pour comprendre leur nature est aujourd’hui l’inversion du tenseur des moments. Au cours des dernières années, les inversions du tenseur des moments se concentraient principalement sur la compréhension du mécanisme physique à l’origine des séismes LP qui souvent supposaient des milieux géologiques très simples, voire homogènes. Des études récentes ont montré l’influence des sédiments peu consolidés à la surface des volcans sur la propagation des ondes à basse fréquence et en conséquence, sur l’inversion du tenseur des moments quand ils ne sont pas pris en compte dans le processus d’inversion. Le but de cette thèse est de mieux comprendre les processus physiques qui génèrent les séismes LP et de quantifier les incertitudes liées à leur interprétation. / Long-period (LP) seismic events are abundantly recorded during rest and unrest periods at many volcanoes worldwide. However, their source mechanism is still poorly understood. Models which have been proposed so far to explain their origin are: 1) the resonance of a fluid-filled cavity triggered by fluid instabilities or the brittle failure of magma; 2) slow-rupture earthquakes occurring in the low consolidated materials composing the shallow portion of the volcanic edifice. Nowadays the main tool used to get insights into their nature is moment tensor (MT) inversion. MT inversions carried out in the past years focused mainly on the understanding of the physical origin of LP events and often supposed a relative simple geological structure of the medium. Recent studies highlighted the strong influence of shallow unconsolidated materials on the retrieved MT solutions and the importance of considering geological inhomogeneity in the inversion process. The principal aim of this thesis is to gain a better understanding of the source processes that generate LP events and to quantify the uncertainties related to the MT inversion process.

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.

La caractérisation de la source d’un séisme se fait à partir de l’analyse des mesures des déplacements transitoires et statiques du sol, et dépend de la quantité et de la qualité de ces mesures. Nous avons travaillé sur la détermination des mécanismes au foyer des séismes modérés de Tadjena (Mw 5.0, 2006), Béni-Ilmane (Mw 5.5, 2010), Hammam Melouane (Mw 5.0, 2013), Bordj-Ménaïel (Mw 4.1, 2014), Alger (Mw 5.7, 2014) et M’ziraa (Mw 5.1, 2016). Pour cela, nous avons inversés les formes d’ondes en champ proche et régional des stations large bandes et accélérométriques du réseau sismologique algérien (ADSN). Nous avons déterminé la distribution spatio-temporelle du glissement cosismique du séisme de Boumerdes-Zemmouri (Mw 6.8, 2003) en inversant conjointement l’ensemble des données disponibles, sismologiques (télésismiques, accélérometriques) et géodésiques (GPS, InSAR, soulèvement côtier). Nous avons travaillé sur la relation qui peut exister entre le séisme de Boumerdes-Zemmouri et les séismes de Hammam Melouane, Bordj Ménaïel et Alger (Mw 5.7, 2014) en termes de transfert de contrainte (CFF). Nous avons calculé le champ de contrainte dans différentes régions de l’Algérie par inversion des mécanismes aux foyers disponibles dans chaque région. / Studies of the earthquake source are based on observations of seismic and static ground motions, which depend on the quality and the density of measurements. In this thesis, we present studies of the determination of the focal mechanism of the Tadjena (Mw 5.0, 2006), Béni-Ilmane (Mw 5.5, 2010), Hammam Melouane (Mw 5.0, 2013), Bordj-Ménaïel (Mw 4.1, 2014), Alger (Mw 5.7, 2014) and M’ziraa (Mw 5.1, 2016) earthquakes, by modeling waveforms in the near and regional field with broadband and strong motion stations of the Algerian Digital Seismic Netwok (ADSN). In addition, we determined the coseismic slip distribution of the Boumerdes-Zemmouri earthquake (Mw 6.8, 2003) by inverting a most comprehensive set of data (teleseismic, strong motion, GPS, InSAR, coastal uplift). We calculated the Coulomb Failure Function between the Boumerdes-Zemmouri earthquake (source fault) and the Hammam Melouane, Bordj Ménaïel and Algiers (Mw 5.7, 2014) events (receiver faults). We computed the stress tensor in different regions of Algeria obtained from the inversion of the available focal mechanisms.

La prédiction des mouvements du sol générés par un séisme est un enjeu majeur pour la prise en compte du risque sismique. C’est l’un des objectifs du projet SIGMA dans le cadre duquel j’ai réalisé ma thèse. Celle-ci se compose de deux parties. La première se concentre sur la dépendance à la magnitude de la décroissance des paramètres des mouvements du sol avec la distance. Celle-ci est un sujet de préoccupation aussi bien pour l’utilisation des relations d’atténuation (GMPEs), que pour les méthodes basées sur l’utilisation de petits évènements en tant que fonctions de Green empiriques. Nous avons démontré qu’aux distances les plus faibles (inférieures à la longueur de la faille), l'effet de saturation dû aux dimensions de la faille est prépondérant. Aux distances plus importantes, l'effet de l’atténuation anélastique devient prépondérant. Nous avons donc montré qu’il pouvait être délicat de mélanger des données de différentes régions dans les GMPEs, et validé l’utilisation des fonctions de Green empiriques à toutes les distances. Dans la deuxième partie sont testées 3 différentes méthodes de simulations dans un contexte complexe : un code combinant une source étendue en k2 et des EGFs, un code point-source EGFs et un code stochastique. Nous avons choisi de travailler sur le séisme de magnitude Mw 5.9 (29 mai 2012) situé dans un bassin sédimentaire profond (la plaine du Po), et qui a engendré des sismogrammes souvent dominés par les ondes de surface. On y démontre que sans connaissance à priori du milieu de propagation, les méthodes basées sur des EGF permettent de reproduire les ondes de surface, les valeurs de PGA, de PGV, ainsi que les durées des signaux générés. / The prediction of ground motion generated by an earthquake is a major issue for the consideration of seismic risk. This is one of the objectives of SIGMA project in which I realized my thesis. It consists of two parts. The first focuses on the magnitude dependence of the ground motion parameters decay with distance. This is a concern both for the use of relation of attenuation (GMPEs) than methods based on the use of small events as empirical Green functions. We have shown that as the shorter distances (less than the length of the fault), the saturation effect due to the fault size is preponderant. For larger distances, it’s the eanelastic attenuation effect which becomes predominant. So we have shown that it can be tricky to mix data from different regions in GMPEs and we validated the use of empirical Green functions at every distance. In the second part are tested three different simulation methods in a complex context: a code combining finite fault source in k2 and EGFs, a point-source code with EGFs and a stochastic code. We chose to work on the Mw 5.9 earthquake (May 29, 2012) which occurs in a deep sedimentary basin (the Po plain), and which has generated seismograms often dominated by surface waves. We show that without a priori knowledge of the propagation medium, methods based on EGFs can reproduce surface waves, the values of PGA, PGV, and the durations of the signals generated.

De nombreuses études montrent que les transferts de contrainte co- et postsismiques jouent un rôle majeur dans l’occurrence des séquences de séismes. Cependant, la grande majorité de ces études implique des systèmes de failles à la configuration géométrique simple (e.g. failles parallèles ou colinéaires). Dans cette thèse, nous étudions une séquence de séismes s’étant produite au sein d’un système de failles à la configuration géométrique plus complexe (i.e. failles conjuguées), la séquence du NE Lut (1939-1997, NE Iran), afin d’évaluer (1) si les transferts de contrainte favorisent la succession de séismes de la séquence et (2) s’ils permettent sur le long-terme de synchroniser les ruptures des failles du système. Pour cela, nous mesurons d’abord les déformations de surface produites par la séquence afin de mieux contraindre par la suite la modélisation des transferts de contrainte. A partir de la technique de corrélation subpixel d'images optiques, nous mesurons les champs de déplacements de surface horizontaux produits par les séismes de Khuli-Boniabad (Mw 7.1, 1979) et de Zirkuh (Mw 7.2, 1997). Nous montrons que ces séismes sont caractérisés par la rupture de plusieurs segments dont les limites sont corrélées avec les complexités géométriques des failles. Nous interprétons les différences de leurs caractéristiques de rupture (longueur de rupture, glissement moyen, nombre de segments rompus) comme étant dues à des différences de maturité des failles de Dasht-e-Bayaz et d’Abiz. Nous détectons également les déplacements produits par un séisme historique modéré, le séisme de Korizan (Mw 6.6, 1979). C’est la première fois que les déplacements produits par un séisme historique de si petite taille sont mesurés par corrélation d’images optiques. Ensuite, en combinant le champ de déplacements InSAR déjà publié avec les données optiques proche-faille précédemment acquises, nous estimons un nouveau modèle de source pour le séisme de Zirkuh (Mw 7.2, 1997). Nous montrons que les données proche-faille sont essentielles pour mieux contraindre la géométrie de la rupture et la distribution du glissement en profondeur. Le modèle estimé montre que le séisme de Zirkuh a rompu trois aspérités séparées par des barrières géométriques où les répliques du séisme se localisent. Seul le segment central de la faille présente un déficit de glissement en surface que nous interprétons comme étant dû à de la déformation distribuée dans des dépôts quaternaires non consolidés. Enfin, à partir des informations précédemment acquises, nous modélisons les transferts de contrainte au cours de la séquence du NE Lut. Nous montrons que ceux-ci ont favorisé l’occurrence de 7 des 11 séismes de la séquence et que modéliser précisément la géométrie des ruptures est essentiel à une estimation robuste des transferts de contrainte. De plus, nous montrons que l’occurrence du séisme de Zirkuh (Mw 7.2, 1992) est principalement favorisée par les séismes modérés de la séquence. Pour finir, la simulation d’une multitude de cycles sismiques sur les failles du NE Lut montre que les transferts de contrainte, en particulier les transferts postsismiques liés à la relaxation viscoélastique de la lithosphère, sont le principal processus permettant la mise en place répétée de séquences de séismes sur les failles du NE Lut. Enfin, d'après les simulations réalisées, l'ordre dans lequel se sont produits les séismes majeurs durant la séquence du NE Lut est assez exceptionnel. / Many studies show that static and postseismic stress transfers play an important role in the occurrence of seismic sequences. However, a large majority of these studies involves seismic sequences that occurred within fault systems having simple geometric configurations (e.g. collinear or parallel fault system). In this thesis, we study a seismic sequence that occurred within a complex fault system (i.e. conjugate fault system), the NE Lut seismic sequence (1939-1997, NE Iran), in order to assess if (1) stress transfers can explain the succession of earthquakes in the sequence and (2) stress transfers can lead to the synchronization of the NE Lut faults over multiple seismic cycles. To this end, we first measure the surface displacement field produced by the sequence in order to precisely constrain the stress transfer modeling afterwards. We use optical correlation technique to measure the surface displacement fields of the Khuli-Boniabad (Mw 7.1, 1979) and Zirkuh earthquake (Mw 7.2, 1997). We find that these earthquakes broke several segments limited by geometrical complexities of the faults. We interpret the differences in failure style of these earthquakes (i.e. rupture length, mean slip and number of broken segments) as being due to different level of structural maturity of the Dasht-e-Bayaz and Abiz faults. Furthermore, we succeed to detect offsets produced by the 1979 Mw 6.6 Korizan earthquake. It is the first time that surface displacements for such a small historical earthquake have been measured using optical correlation. Then, combining previously published intermediate-field InSAR data and our near-field optical data, we estimate a new source model for the Zirkuh earthquake (Mw 7.2, 1997). We show that near-field data are crucial to better constrain the fault geometry and the slip distribution at depth. According to our source model, the Zirkuh earthquake broke three asperities separated by geometrical barriers where aftershocks are located. No shallow slip deficit is found for the overall rupture except on the central segment where it could be due to off-fault deformation in quaternary deposits. Finally, we use the information acquired in the first parts of this work to model the stress transfers within the NE Lut sequence. We find that 7 out of 11 earthquakes are triggered by the previous ones and that the precise modeling of the rupture geometry is crucial to robustly estimate the stress transfers. We also show that the Zirkuh earthquake is mainly triggered by the moderate earthquakes of the NE Lut sequence. Lastly, the simulation of multiple seismic cycles on the NE Lut fault system shows that stress transfers, in particular postseismic stress transfers due to viscoelastic relaxation, enhance the number of seismic sequences and synchronize the rupture of the faults. The simulations also show that the order in which the Mw>7 earthquakes occurred during the NE Lut sequence is quite exceptional.

La caractérisation de la source d’un séisme se fait à partir de l’analyse des mesures des déplacements transitoires et statiques du sol, et dépend de la quantité et de la qualité de ces mesures. Nous avons travaillé sur la détermination des mécanismes au foyer des répliques du séisme de Saintes (MW 6.4, 2004), et sur la détermination de la distribution spatio-temporelle du glissement cosismique des séismes de L’Aquila (Mw 6.3, 2009), et de Miyagi-Oki (Mw 7.2, 2005) et de Sanriku-Oki (Mw 7.3, 2011). Ces travaux se sont basés sur des méthodes d’inversions, et différents jeux de données (accélérométriques, large-bandes, GPS et InSAR) accessibles ou non selon le séisme considéré. La seule diversité des mesures n’est pas suffisante pour décrire la rupture. La modélisation des données se confronte à des difficultés, comme par exemple la pertinence des modèles de vitesses sismiques pour la modélisation des données accélérométriques. Une autre problématique récurrente est la non-unicité de la meilleure solution déterminée par les méthodes d’inversions pour décrire les données. Pour répondre à ces deux problématiques, nous avons d‘une part développé une procédure d’exploration de modèles de vitesse pour déterminer les valeurs optimales capables de décrire au mieux les données accélérométriques du séisme de L’Aquila. D’autre part, nous avons développé une procédure de construction d’un modèle de source moyen que nous avons appliqué pour la détermination du glissement cosismique des séismes de L’Aquila, de Miyagi-Oki, et de Sanriku-Oki. L’ensemble de ces travaux et les réponses aux problèmes soulevés sont présentés dans ce travail de thèse. / Studies of the earthquake source are based on observations of seismic ground motions. They also depend on the quality and the density of measurements. In this present work we will present studies of the determination of focal mechanism of main aftershocks of the Les Saintes (MW 6.4, 2004) earthquake, and the determination of the coseismic slip of the L’Aquila (MW 6.3, 2009), the Miyagi-Oki (MW 7.2, 2005), ant the Sanriku-Oki (MW 7.3, 2011) earthquakes. These studies were based on two inversion methods. Different kinds of data were available (strong motion, broadband teleseismic, GPS and InSAR) depending on the earthquake studied. But the multiplicity of data is not sufficient to well describe rupture process. There are others difficulties as the data modeling of strong motion. Seismic velocity models are used to describe the characteristics of layers crossed by seismic waves. The quality of the modeling is depending on the pertinence of these seismic velocity models. The description of the rupture process is also depending on the non-uniqueness of the best solution given by global inversion methods. We propose two procedures in order to take into account these two classic issues. First, we developed a velocity model exploration procedure to obtain optimized 1D velocity models in order to improve the strong motion modeling of the L’Aquila earthquake. Then we developed a procedure to build an average rupture model from the combined results of several joint inversions, which was applied to the L’Aquila, the Miyagi-Oki, and the Sanriku-Oki earthquake. This thesis presents all these works and answers to the raised issues.

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.

Ce mémoire d'habilitation décrit ma recherche effectuée en vue de mieux comprendre le mécanisme des tremblements de terre. Les ondes sismiques sont les principales données qui donnent accès à cette information, et mes travaux sont donc en lien direct avec les observations des capteurs sismiques (ou géodésiques) capables de détecter ces ondes. Les qualité, quantité et type de données disponibles m'ont conduit à aborder cette thématique sous trois angles différents, qui recoupent les trois chapitres de ce mémoire. A l'échelle mondiale, les capteurs large-bande du réseau global permettent une étude systématique et homogène de tous les séismes dont la magnitude est suffisante pour générer des ondes clairement analysables. Cette magnitude seuil est de l'ordre de 5.5-6 lorsque l'on analyse les ondes de volume se propageant dans la Terre. Par ailleurs, ces données sont accessibles en temps réel, ce qui permet d'obtenir des informations rapides sur les séismes, à condition que des techniques d'analyse automatisées soient mises en place. Ce double intérêt de l'approche globale -systématisme et rapidité- est développé dans le Chapitre I de ce mémoire. L'échelle mondiale trouve ses limites lorsqu'on s'intéresse aux détails du processus de rupture. Alors que les données mondiales sont suffisantes pour imager les caractéristiques moyennes de source (mécanisme au foyer, profondeur, magnitude de moment, fonction source), elles ne permettent généralement pas d'extraire avec précision les informations internes à la rupture sismique (distribution de glissement, vitesse de rupture locale). Pour ce faire, il est nécessaire d'utiliser des données plus proches de l'événement sismique, et de développer des techniques d'analyse adaptées à chaque configuration (analyses en réseau, analyses comparées entre le séisme principal et l'un de ses précurseurs...). Par ailleurs, et contrairement à l'échelle globale - où les données large-bande sont les données " reines " de l'analyse de la source -, l'analyse à distance plus proche permet et/ou requiert de diversifier les données utilisées. Cela peut passer par l'utilisation des capteurs accélérometriques, ou par la nouvelle utilisation des capteurs GPS en tant que sismomètres. Le Chapitre II s'intéresse à cette thématique de résolution fine des propriétés de source, avec un intérêt particulier pour la détermination de la vitesse de rupture. La mise en évidence récente des vitesses de rupture supershear (plus rapides que les ondes de cisaillement) et de leurs conséquences constituent un point important de ce chapitre. Enfin, il est intéressant de replacer la rupture sismique dans le cadre plus large des conditions qui stimulent ou inhibent son déclenchement. Dans cet esprit, je présente dans le Chapitre III des éléments de réponse venant de l'observation continue de la zone de subduction Equatorienne. Cette observation - incluant capteurs large-bande, accéléromètres, et GPS - se déroule depuis 2008 dans le cadre du projet ANR ADN, soutenu par l'IRD et en collaboration avec l'Institut de Géophysique à Quito. Nous montrerons dans ce chapitre comment l'occurrence de séismes lents sur l'interface de subduction est un facteur déclenchant de la sismicité.

Les séismes sont à l'origine d'évènements sédimentaires gravitaires (turbidites) dont l'étude permet de reconstituer l'histoire mal connue de la paléosismicité des marges continentales. L'analyse d'une série de carottes de sédiments, collectées stratégiquement dans trois systèmes turbiditiques de la marge en subduction Hikurangi de Nouvelle-Zélande, permet d'établir les caractéristiques, les facteurs de contrôle et les mécanismes déclencheurs de la sédimentation gravitaire des derniers 18,000 ans. La succession sédimentaire comprend quatre lithofaciès et modes de dépôt : hémipélagite (sédimentation marine), turbidites (courants de turbidité), débrites (débris flows) et tephra (retombée de cendres volcaniques), dont l'organisation générale dépend de la morphologie de la marge et des fluctuations glacio-eustatiques. Les crues, les éruptions volcaniques et les " slope failures " sont les trois mécanismes déclencheurs des turbidites. Plus de 90% sont déclenchées par des " slope failures " en haut de pente (150 - 1,200 m) à la suite de séismes. L'adaptation d'un modèle empirique de stabilité de pente suggère que ces turbidites représentent l'enregistrement sédimentaire des ruptures répétées de trois failles actives, dont l'interplaque, et correspondent à des séismes de Mw≥7.3 avec un temps de retour de 150±50 ans. Parmi ces turbidites co-sismiques, 20 montrent une synchronicité de déclenchement sur l'ensemble de la marge et un volume important. Elles correspondent à des ruptures de la zone interplaque de Mw 7.5 - 8.4, dont les temps de retour montrent une phénomène de clustering où alternent les périodes actives à faible temps de retour (305 - 610 ans), et les périodes de quiescence à temps de retour élevé (1480 - 2650 ans). Ce calendrier paléosismique intégré aux modélisations en cours devrait permettre de mieux appréhender l'aléa sismique et les risques pour la population.