La marge Ligure est une marge passive soumise à une déformation tectonique compressive associée à la tectonique salifère messinienne. La reprise en compression de la marge s’accompagne d’une sismicité modérée récurrente ponctuée d’évènements plus forts. La marge Ligure est le siège d’une sédimentation importante au Plio-Quaternaire. Elle constitue un environnement propice à l’étude des déstabilisations gravitaires. L’objectif principal de ce travail était de décrire et caractériser les principaux mouvements en masse ayant affecté la marge continentale Ligure au cours du Plio-Quaternaire, de localiser les principales zones sujettes aux déstabilisations et d’identifier les facteurs pré-conditionnant et déclenchant les ruptures dans le but de mieux évaluer l’aléa gravitaire. Une large couverture de données bathymétriques, géophysiques et des carottages acquis sur l’ensemble de la marge a permis de réaliser une étude multi-échelles des processus de ruptures gravitaires et des facteurs déclenchant associés. Près de 1500 glissements ont été identifiés. L’étude de leur répartition spatio-temporelle illustre que l’ensemble de la marge a toujours été affectée par des déstabilisations de pente mais que les principales zones de ruptures auraient migré vers l’ouest au cours du Plio-Quaternaire. Les grandes ruptures sous-marines sont préférentiellement associées aux zones de déformation maximale, cette dernière étant contrôlée par la tectonique crustale et/ou la tectonique salifère. Il apparaît que les ruptures résultent plus généralement d’une association de facteurs distincts qui ont participé à fragiliser la stabilité des dépôts de la pente et qui ont pu provoquer leur rupture. / The Ligurian margin is a passive margin characterized by high sedimentation rates during the Plio-Quaternary. It is affected by a compressive tectonic deformation leading to the inversion of the margin, together with a salt tectonic. The present-day moderate seismic activity is punctuated by stronger seismic events. Thus, this margin offers a good natural laboratory to study submarine landslides and their triggering factors. Although the Var Turbidite System has been well investigated over the last 20 years, the morphology and tectonics/sedimentary processes affecting the whole margin remained poorly known. This study aims to describe and to characterize the main types of mass movements, their preferential locations along the Ligurian margin during the Plio-Quaternary and their triggering factors to improve geohazards assessment related to landslides. A dataset including bathymetric and geophysical data and cores allowed to realize a multi-scale study of submarine failures and their associated triggering factors. About 1500 landslides were identified on the margin and in the basin. The study of their spatio-temporal distribution revealed that the margin has always been affected by mass-wasting processes and that the main zones of landsliding migrated westward during the Plio-Quaternary. The largest submarine landslides are preferentially associated with the highest deformation rates and their location is controlled by crustal tectonics and/or salt tectonics. The initiation of failures results from the combination of several factors including the margin deformation, earthquakes, salt tectonics and sediment under-consolidation.

Dans ce travail de thèse, nous analysons la structure crustale de la marge est-algérienne et du bassin adjacent (région d’Annaba), à partir d’un ensemble de nouvelles données acquises durant la Campagne SPIRAL’2009 incluant un profil sismique terre-mer de ~240 km de long, des lignes sismiques réflexion pénétrante 360-traces, et des profils gravimétriques et magnétiques. Nous avons par ailleurs disposé pour cette étude de données complémentaires incluant notamment un ensemble de profils de sismique réflexion offrant des résolutions complémentaires. La structure crustale ainsi établie nous permet de discuter les nombreux modèles cinématiques d’ouverture du bassin est-algérien proposés dans la littérature, afin de caler dans le temps la formation du bassin par rapport à la collision. Elle permet également de discuter la localisation de la déformation liée à la réactivation de la marge, par rapport aux grands domaines lithosphériques du système marge-bassin, afin de mieux comprendre les modalités de l’inversion. Dans le bassin profond, la modélisation directe des temps d’arrivée et des amplitudes des ondes réfractées et réfléchies met en évidence une croûte océanique anormalement mince de 5-5.5 km d’épaisseur, composée de deux couches. La première, de 2.2 km d’épaisseur, montre des vitesses comprises entre 4.8 à 6.0 km/s impliquant un fort gradient; la seconde de 3.3 km d’épaisseur, présente des vitesses comprises entre 6.0 à 7.1 km/s et un plus faible gradient de vitesse. La modélisation des temps d’arrivées des ondes S fourni pour cette couche un coefficient de Poisson de 0.28, indiquant qu’elle est majoritairement constituée de gabbros. / In this study, we determine the deep structure of the eastern Algerian basin and its southern margin in the Annaba region (easternmost Algeria), to better constrain the plate kinematic reconstruction in this region. This study is based on new geophysical data collected during the SPIRAL cruise in 2009 that included a wide-angle, 240-km-long, onshore-offshore seismic profile, multichannel seismic reflection lines, and gravity and magnetic data, which was complemented by the available geophysical data for the study area. The analysis and modeling of the wide-angle seismic data using travel-times and amplitudes, and integrated with the multichannel seismic lines, reveal the detailed structure of an ocean-to-continent transition. In the deep basin, there is an ~5.5-km-thick oceanic crust that is composed of two layers. The upper layer of the crust is defined by a high velocity gradient and P-wave velocities between 4.8 km/s and 6.0 km/s from the top to the bottom. The lower crust is defined by a lower velocity gradient and P-wave velocity between 6.0 km/s and 7.1 km/s. The Poisson ratio in the lower crust deduced from S-wave modeling is 0.28, which indicates that the lower crust is composed mainly of gabbros. Below the continental edge, a typical continental crust with P-wave velocities between 5.2 km/s and 7.0 km/s from the top to the bottom shows a gradual seaward thinning of ~15 km over an ~35-km distance.

En termes de prévention du risque associé aux séismes, la prédiction quantitative des phénomènes de propagation et d'amplification des ondes sismiques dans des structures géologiques complexes devient essentielle. Dans ce domaine, la simulation numérique est prépondérante et l'exploitation efficace des techniques de calcul haute performance permet d'envisager les modélisations à grande échelle nécessaires dans le domaine du risque sismique. Plusieurs évolutions récentes au niveau de l'architecture des machines parallèles nécessitent l'adaptation des algorithmes classiques utilisées pour la modélisation sismique. En effet, l'augmentation de la puissance des processeurs se traduit maintenant principalement par un nombre croissant de coeurs de calcul et les puces multicoeurs sont maintenant à la base de la majorité des architectures multiprocesseurs. Ce changement correspond également à une plus grande complexité au niveau de l'organisation physique de la mémoire qui s'articule généralement autour d'une architecture NUMA (Non Uniform Memory Access pour accès mémoire non uniforme)~de profondeur importante. Les contributions de cette thèse se situent à la fois au niveau algorithmique et numérique mais abordent également l'articulation avec les supports d'exécution optimisés pour les architectures multicoeurs. Les solutions retenues sont validées à grande échelle en considérant deux exemples de modélisation sismique. Le premier cas se situe dans la préfecture de Niigata-Chuetsu au Japon (événement du 16 juillet 2007) et repose sur la méthode des différences finies. Le deuxième exemple met en oeuvre la méthode des éléments finis. Un séisme hypothétique dans la région de Nice est modélisé en tenant compte du comportement non linéaire du sol.

En termes de prévention du risque associé aux séismes, la prédiction quantitative des phénomènes de propagation et d'amplification des ondes sismiques dans des structures géologiques complexes devient essentielle. Dans ce domaine, la simulation numérique est prépondérante et l'exploitation efficace des techniques de calcul haute performance permet d'envisager les modélisations à grande échelle nécessaires dans le domaine du risque sismique.Plusieurs évolutions récentes au niveau de l'architecture des machines parallèles nécessitent l'adaptation des algorithmes classiques utilisées pour la modélisation sismique. En effet, l'augmentation de la puissance des processeurs se traduit maintenant principalement par un nombre croissant de cœurs de calcul et les puces multicœurs sont maintenant à la base de la majorité des architectures multiprocesseurs. Ce changement correspond également à une plus grande complexité au niveau de l'organisation physique de la mémoire qui s'articule généralement autour d'une architecture NUMA (Non Uniform Memory Access pour accès mémoire non uniforme) de profondeur importante.Les contributions de cette thèse se situent à la fois au niveau algorithmique et numérique mais abordent également l'articulation avec les supports d'exécution optimisés pour les architectures multicœurs. Les solutions retenues sont validées à grande échelle en considérant deux exemples de modélisation sismique. Le premier cas se situe dans la préfecture de Niigata-Chuetsu au Japon (événement du 16 juillet 2007) et repose sur la méthode des différences finies. Le deuxième exemple met en œuvre la méthode des éléments finis. Un séisme hypothétique dans la région de Nice est modélisé en tenant compte du comportement non linéaire du sol. / One major goal of strong motion seismology is the estimation of damage in future earthquake scenarios. Simulation of large scale seismic wave propagation is of great importance for efficient strong motion analysis and risk mitigation. Being particularly CPU-consuming, this three-dimensional problem makes use of high-performance computing technologies to make realistic simulation feasible on a regional scale at relatively high frequencies.Several evolutions at the chip level have an important impact on the performance of classical implementation of seismic applications. The trend in parallel computing is to increase the number of cores available at the shared-memory level with possible non-uniform cost of memory accesses. The increasing number of cores per processor and the effort made to overcome the limitation of classical symmetric multiprocessors SMP systems make available a growing number of NUMA (Non Uniform Memory Access) architecture as computing node. We therefore need to consider new approaches more suitable to such parallel systems.This PhD work addresses both the algorithmic issues and the integration of efficient programming models for multicore architectures. The proposed contributions are validated with two large scale examples. The first case is the modeling of the 2007 Niigata-Chuetsu, Japan earthquake based on the finite differences numerical method. The second example considers a potential seismic event in the Nice sedimentary basin in the French Riviera. The finite elements method is used and the nonlinear soil behavior is taken into account.

Ce travail de thèse a pour but d'analyser les variations de vitesses sur des grandes failles décrochantes en contexte intracontinental, capables de produire des séismes de très fortes magnitudes (M > 7.5). Afin d'illustrer c es variations d'activités, cette analyse a été effectuée sur deux zones d'études situées en domaine continental et sismiquement actives: la région ouest de la Mongolie (failles de Bogd et Bolnay) et le nord de l'Iran (failles d'Astaneh et de Tabriz). À partir d'une approche morphotectonique et paléosismologique, les cinématiques, les vitesses de failles et les intervalles de récurrence entre les séismes majeurs ont été estimés, permettant d'analyser les caractéristiques du cycle sismique sur chacune des failles. En Mongolie, les failles de Bogd et Bolnay présentent respectivement des vitesses de ~ 1,2 et 2,6 mm/an, qui semblent être constantes sur la période Pleistocène supérieur-Holocène. Ces deux failles présentent également des glissements caractéristiques et des intervalles de temps similaires entre les séismes majeurs. Les analyses paléosismologiques suggèrent qu'un essaim sismique comparable à celui enregistré au XXème siècle a eu lieu il y a environ 3000 ans.En Iran, une vitesse géologique de 2 mm/an a été estimée sur la faille d'Astaneh et les données paléosismologiques suggèrent des intervalles de récurrence qui varient entre 1600 et 2200 ans, associés à des déplacements en surface compris entre 3 et 4,5 m. Nous avons également estimé une vitesse de 7 mm/an sur la faille de Tabriz, en accord avec les données GPS, suggérant que la vitesse sur cette faille est constante depuis 45 ka. / The aim of this thesis is to analyze if variations in slip rates occur along strike-slip faults, in intracontinental domain, these faults producing large earthquakes (M> 7.5). To illustrate these variations, this work has focused in two area located in the most tectonically active continental domains in the world: in the western part of Mongolia (Bogd and Bulnay faults) and in the northern part of Iran (Astaneh and Tabriz faults). Using morphotectonic and paleoseismological analyses, the kinematics, the slip rates and the recurrence times have been estimated and allow us to describe the characteristics of seismic cycle along these faults. In Mongolia, the slip rates are estimated at ~1.2 and ~2.6 mm/yr along the Bogd and the Bulnay faults, respectively, with no variations of geological slip rates over the Pleistocene-Holocene period. These two faults present characteristic slips and similar recurrence times between large earthquakes. The paleoseismological investigations suggest that a cluster occurred 3000 years ago, similar to the seismic cluster recorded in Mongolia during the XX century. In Iran, the slip rate was estimated to 2 mm/yr along the Astaneh fault and the recurrence times are ranging from 1600 to 2200 years, associated with offsets comprised between 3 and 4.5 m. We have also estimated a slip rate of ~7 mm/yr along the Tabriz fault, in agreement with the present day rate estimated by GPS, suggesting no variations in the slip rate over the past 45 ka.

La donnée sismique depuis sa conception (modélisation d’acquisitions sismiques), dans sa phase de traitement (prétraitement et migration) et jusqu’à son exploitation pour en extraire les informations géologiques pertinentes nécessaires à l’identification et l’exploitation optimale des réservoirs d’hydrocarbures (interprétation), génère un volume important de calculs. Nous montrons dans ce travail de thèse qu’à chacune de ces étapes l’utilisation de technologies accélératrices de type GPGPU permet de réduire radicalement les temps de calcul tout en restant dans une enveloppe de consommation électrique raisonnable. Nous présentons et analysons les éléments sous-jacents à ces performances. L’importance de l’utilisation de motifs d’accès mémoire adéquats est particulièrement mise en exergue étant donné que l’accès à la mémoire représente le principal goulot d’étranglement pour les algorithmes abordés. Nous reportons des facteurs d’accélération de l’ordre de 40 pour la modélisation sismique par résolution de l’équation d’onde par différences finies (brique de base pour la modélisation et l’imagerie sismique) et entre 8 et 113 pour le calcul d’attributs sismiques. Nous démontrons que l’utilisation d’accélérateurs matériels élargit considérablement le champ du possible, aussi bien en imagerie sismique (modélisation de nouveaux types d’acquisitions à grande échelle) qu’en interprétation (calcul d’attributs complexes sur station de travail, paramétrage interactif des calculs, etc.). / During the seismic imaging workflow, from seismic modeling to interpretation, processingseismic data requires a massive amount of computation. We show in this work that, at eachstage of this workflow, hardware accelerators such as GPUs may help reducing the time requiredto process seismic data while staying at reasonable energy consumption levels.In this work, the key programming considerations needed to achieve good performance are describedand discussed. The importance of adapted in-memory data access patterns is particularlyemphasised since data access is the main bottleneck for the considered algorithms. When usingGPUs, speedup ratios of 40× are achieved for FDTD seismic modeling, and 8× up to 113× forseismic attribute computation compared to CPUs.

L’amélioration des images sismiques peut aider à mieux contraindre l’exploration deshydrocarbures. Les ondes élastiques qui se propagent dans la Terre peuvent être classifiéescomme ondes de volume et ondes de surface. Si ces dernières sont les plus énergétiques,seules les ondes de volume sont couramment considérées comme des signaux utiles.Cependant, les ondes de surface sont utiles pour caractériser la proche sub-surface.Classiquement, les ondes de surface sont analysées dans des contextes de propriétésélastiques localement 1D.Nous proposons une modification de l’inversion des formes d’onde classique pourreconstruire des profils de propriétés 2D (la windowed-Amplitude Waveform Inversion, w-AWI). La w-AWI est spécialement robuste en ce qui concerne le choix du modèle initial.Nous appliquons la w-AWI aux données synthétiques ainsi qu’aux données réelles, montrantque cette approche permet de récupérer des propriétés 2D. / High-resolution seismic imaging is essential to improve results of hydrocarbon exploration.Elastic waves propagate in the Earth as body and surface waves, the latter being the mostenergetic ones. Body waves are preferred for exploration seismic imaging while surfacewaves are usually considered to be noise. However, it has been recognised that the nearsurface can be characterised by analysing surface waves and that such result may improvethe outcome of body-wave processing. Currently, surface waves analysis leads to retrievelocal 1D property profiles.We propose a waveform-based inversion procedure to derive 2D velocity models fromsurface waves. This method consists of a misfit functional modification of classical FullWaveform Inversion and we call it windowed-Amplitude Waveform Inversion (w-AWI). Weshow that w-AWI is robust regarding the choice of initial velocity model. We apply w-AWI tosynthetic and real data obtaining encouraging near-surface imaging results

La caractérisation des applications en vue de les préparer pour les nouvelles architectures et les porter sur des systèmes très étendus est une étape importante pour pouvoir anticiper les modifications nécessaires. Comme les machines Exascale sont prévues pour la période 2018-2020, l'étude des applications et leur préparation pour ces machines s'avèrent donc essentielles. Nous nous intéressons aux applications d'imagerie sismique et en particulier à l'application Reverse Time Migration (RTM) car elle est très utilisée par les pétroliers dans le cadre de l'exploration sismique.La première partie de nos travaux a porté sur l'étude du cœur de calcul de l'application RTM qui consiste en un calcul de différences finies dans le domaine temporel (FDTD). Nous avons caractérisé cette partie de l'application en soulevant les aspects architecturaux des machines actuelles ayant un fort impact sur la performance, notamment les caches, les bandes passantes et le prefetching. Cette étude a abouti à l'élaboration d'un modèle de performance permettant de prédire le trafic DRAM des FDTD. La deuxième partie de la thèse se focalise sur l'impact de l'hétérogénéité et le parallélisme sur la FDTD et sur RTM. Nous avons choisi l'architecture manycore d’Intel, Xeon Phi, et nous avons étudié une implémentation "native" et une implémentation hétérogène et hybride, la version "symmetric". Enfin, nous avons porté l'application RTM sur un cluster hétérogène, Stampede du Texas Advanced Computing Center (TACC), où nous avons effectué des tests de scalabilité allant jusqu'à 64 nœuds contenant des coprocesseurs Xeon Phi et des processeurs Sandy Bridge ce qui correspond à presque 5000 cœurs / As we are expecting Exascale systems for the 2018-2020 time frame, performance analysis and characterization of applications for new processor architectures and large scale systems are important tasks that permit to anticipate the required changes to efficiently exploit the future HPC systems. This thesis focuses on seismic imaging applications used for modeling complex physical phenomena, in particular the depth imaging application called Reverse Time Migration (RTM). My first contribution consists in characterizing and modeling the performance of the computational core of RTM which is based on finite-difference time-domain (FDTD) computations. I identify and explore the major tuning parameters influencing performance and the interaction between the architecture and the application. The second contribution is an analysis to identify the challenges for a hybrid and heterogeneous implementation of FDTD for manycore architectures. We target Intel’s first Xeon Phi co-processor, the Knights Corner. This architecture is an interesting proxy for our study since it contains some of the expected features of an Exascale system: concurrency and heterogeneity.My third contribution is an extension of the performance analysis and modeling to the full RTM. This adds communications and IOs to the computation part. RTM is a data intensive application and requires the storage of intermediate values of the computational field resulting in expensive IO accesses. My fourth contribution is the final measurement and model validation of my hybrid RTM implementation on a large system. This has been done on Stampede, a machine of the Texas Advanced Computing Center (TACC), which allows us to test the scalability up to 64 nodes each containing one 61-core Xeon Phi and two 8-core CPUs for a total close to 5000 heterogeneous cores

Les capteurs distribués à fibre optique (aussi nommés DAS) sont une nouvelle technologie d'acquisition sismique qui utilise des câbles traditionnels à fibre optique pour fournir une mesure de la déformation le long du câble. Ce système d'acquisition est largement utilisé dans les profils sismiques verticaux (PSV). Le couplage est un facteur clé qui a une grande influence sur la qualité des données. Alors que, pour les acquisitions PSV, les géophones sont attachés à la paroi du puits, le câble de fibre optique est soit cimenté derrière le tubage, soit attaché avec des pinces rigides au tubage ou simplement descendu dans le puits. Cette dernière stratégie de déploiement donne généralement le plus petit rapport signal sur bruit, mais est considérée comme la plus rentable en particulier pour les installations dans des puits existants. Cette thèse porte sur la problématique du couplage du DAS quand le câble est simplement descendu dans le puits. Nous développons des modèles numériques pour analyser les données réelles. L'interprétation de ces résultats nous permet de conclure qu'un contact immédiat du câble avec la paroi du puits avec une force de contact calculée est nécessaire pour fournir des bonnes conditions de couplage. Sur la base de ces résultats, nous proposons des solutions pour optimiser davantage les acquisitions avec le système DAS. Nous modifions numériquement la force de contact et les propriétés élastiques du câble DAS et démontrons comment ces modifications peuvent améliorer mais aussi détériorer la qualité des données. Enfin, nous proposons un algorithme de détection du couplage qui permet d'assurer l'acquisition de données réelles avec un rapport signal / bruit élevé. / Distributed Acoustic Sensing (DAS) is a new technology of seismic acquisition that relies on traditional fibre-optic cables to provide inline strain measurement. This acquisition system is largely used in vertical seismic profiling (VSP) surveys. Coupling is a key factor influencing data quality. While geophones and accelerometers are clamped to the borehole wall during VSP surveys, the fibre cable is either clamped and then cemented behind the casing, or attached with rigid clamps to the tubing, or loosely lowered into the borehole. The latter deployment strategy, also called wireline deployment, usually acquires the lowest level of signal but is regarded as the most cost-effective in particular for existing well installations. This PhD thesis addresses the problematic of coupling of DAS using wireline deployment. We develop numerical models that are used to analyse real data. The interpretation of these results allows us concluding that an immediate contact of the cable with the borehole wall with a computed contact force is required to provide good coupling conditions. Based on those findings, we propose solutions to further optimise DAS acquisitions. We numerically modify the contact force and the elastic properties of the DAS cable and show how these modifications can improve but also deteriorate data quality. Finally, we propose a coupling detection algorithm that is applied to real datasets and allows ensuring the acquisition of data with a high signal-to-noise ratio.

Ce travail de thèse a porté sur différents aspects de la prédiction des éruptions. La première partie, basée sur une approche statistique de l'information, nous a permis d'aborder la prédictabilité des éruptions de façon globale, et d'avoir une vision intégrée, à différentes échelles, de la dynamique éruptive. La dernière partie basée sur une approche déterministe du signal, nous a emmené de façon plus indirecte vers la prédiction des éruptions par le biais de l'analyse et la modélisation des comportements éruptifs d'un unique volcan à travers la sismicité. Ce travail de thèse nous a donc permis de traverser plusieurs échelles de temps et d'espace, en partant de la plus large, la Terre dans sa globalité, pour se concentrer ensuite sur le cas d'un volcan unique, le Piton de la Fournaise et finalement se focaliser sur une propriété spécifique de la sismicité du Piton de la Fournaise à travers la détermination des mécanismes source. Dans un premier temps, une analyse descriptive des propriétés statistiques de la dynamique éruptive à l'échelle mondiale et l'apparition de lois puissances à différentes échelles nous ont permis de relier la dynamique des éruptions à celle d'un système critique auto-organisé, justifiant de la difficulté rencontrée à prédire les éruptions de façon déterministe et de la nécessité de se tourner vers les méthodes statistiques. L'étude à l'échelle mondiale des lois d'Omori retrouvées pour la sismicité autour des éruptions, ainsi que l'analyse des propriétés de ces lois d'Omori en fonction de la taille des éruptions, suggèrent que l'énergie d'endommagement et de relaxation dans la croûte est distribuée en loi puissance autour des éruptions, mais avec des coefficients de chargement et de relaxation différents de la sismicité tectonique classique, et donc dépendant des processus qui génèrent l'endommagement. En zoomant sur trois observables du Piton de la Fournaise (sismicité, déformation et variations de vitesse sismique), nous avons montré que les comportement respectifs des observables imagent les différentes étapes du processus éruptif, pour différentes échelles de temps et d'espace. L'analyse des comportements pré-éruptifs et l'utilisation de méthodes de "pattern recognition" ont permis de mettre en évidence l'apport de la combinaison des informations des différents précurseurs dans l'amélioration des résultats de prédiction. Enfin, l'analyse plus spécifique et déterministe des mécanismes source au Piton de la Fournaise, à partir d'une méthode développée pour le contexte particulier d'une sismicité de magnitude très faible en milieu fortement hétérogène, a permis de mettre en évidence l'apparition de trois familles d'événements similaires. Deux de ces familles peuvent être interprétées comme l'ajustement du milieu au passage du dyke, la troisième famille correspondant à la signature de l'interaction entre les fluides et l'édifice volcanique au niveau du toit de la chambre magmatique. / This PhD work focused on several aspects of the eruptions prediction. It allowed us to consider different space scales from global scale to a single volcano, Piton de la Fournaise, through statistical approaches. We then focused more specifically on the volcano-seismicity associated to its eruptive processes, using deterministic methods. First, we analyzed the statistical properties of the eruptive dynamics at the global scale. The emergence of recurrent power laws allowed us to relate the eruptive dynamics to the SOC systems, thereby attesting for the difficulty to predict eruptions by means of deterministic methods. At the global scale, the study of the Omori-law found for eruptions foreshocks and aftershocks, along with the analysis of the parameters of those Omori laws given the size of eruptions, suggests that the damage and relaxation energy in the crust around eruption time and location is power-law distributed. However, the loading and relaxation coefficients are found different than for classical tectonic seismicity, thereby attesting their dependency to the damage processes. When zooming on three observables of Piton de la Fournaise volcano (seismicity, deformation and seismic velocity changes), me showed in a third part that the respective behaviors of the three observables respond to the different steps of the eruptive process, at different time and space scales. The analysis of the pre-eruptive behaviors of those forerunners under pattern recognition techniques enhanced the optimization of the prediction results by combining simultaneously their pre-eruptive information. Lastly, we performed a specific and deterministic study of the volcano-tectonic seismicity of Piton de la Fournaise through the assessment of the source mechanisms of the volcano-tectonic seismicity. Given the complexity of volcanic context, we developed a method to retrieve the source mechanisms of small and high-frequency seismic signals. Using this method three main families of events with different locations and times of occurrence emerge. Two of them may be interpreted as the readjustment of the crust during dyke opening whereas the third family is interpreted as the signature of fluids-rock interactions at the top of the magma Storage zone.