La Mongolie occidentale a connu une sismicité intracontinentale exceptionnelle au cours de la première moitié du XXe siècle, avec 4 évènements sismiques d'une magnitude supérieure à 7.9, localisés sur des grandes structures décrochantes en régime transpressif senestre. Cette sismicité historique a été qualifiée d’essaim sismique, et des études paléosismologiques indiquent que ce phénomène a déjà pu se produire par le passé (3 à 4 ka). Mais si la partie est de la chaine du Gobi-Altaï a rompu lors du séisme du même nom en 1957 (Mw8) le long de la faille de Bogd orientale (WBF), plusieurs traces de failles affectant des formations Holocènes et présentant des longueurs > 100 km ont été documentées dans la partie ouest du massif, sans pour autant faire l’objet d’études quantitatives.Cette thèse apporte un état des lieux de l’activité sismique dans le Gobi-Altaï, au travers d’une cartographie détaillée des escarpements de failles actives au cours du pléistocène supérieur, associée à une analyse de la segmentation au regard de critères géométriques et cinématiques. Parallèlement, des investigations morphotectoniques et paléosismologiques, couplées à des méthodes de datations au 10Be in situ, OSL et 14C, ont permis de quantifier les vitesses de glissement depuis le Pléistocène supérieur, et les âges des paléo-ruptures au cours de l’Holocène sur les deux principaux décrochements analysés : les failles de la Vallée des Lacs (VOLF), et de Bogd occidentale (WBF) ; des données préliminaires permettent également de contraindre l’âge de la dernière rupture et la vitesse de glissement sur les failles de Tsogt (TF) et de Tsagaan Gol (TGF), situées respectivement dans le prolongement occidental de la WBF et de la VOLF. De plus, des mesures systématiques des décalages cosismiques horizontaux enregistrés par la morphologie, notamment le long de la rupture de surface du séisme de 1957, couplés à l’utilisation et le développement d’une approche statistique de détermination des glissements cosismiques moyens au regard de la segmentation, a permis de déterminer la fonction de distribution du glissement le long de l’EBF sur 3 cycles sismiques, suggérant d’une part une distribution du glissement assez homogène le long de la rupture, et d’autre part que le glissement du séisme généré en 1957 est d’une amplitude comparable de celle des précédents séismes, pour au moins 50 % de la rupture principale. Cette approche de détermination statistique des glissements moyens sur plusieurs cycles sismiques est prometteuse mais présente néanmoins certaines incohérences, et des voies de développement et de validations sont proposées.Les paramètres morphotectoniques et les âges des paléoséismes documentés dans cette thèse présentent une grande cohérence avec les données précédemment acquises sur l’EBF, avec des vitesses de faille lentes (0,7 ± 0,2 mm.a-1 sur la WBF et 0,5 ± 0,1 mm.a-1 sur la VOLF), des périodes de retour longues (4,0 ± 1,2 ka pour la WBF), et des décalages cosismiques relativement importants (2 – 4 m). Nos mesures de vitesses de faille indiquent par ailleurs une distribution de la déformation depuis la EBF vers les WBF et VOLF, séparant la déformation en deux faisceaux de failles, l’un se branchant au nord depuis la VOLF sur la TGF. Les déterminations des âges de paléoséismes indiquent une possibilité de cluster entre les différentes failles étudiées, avec trois ruptures potentiellement synchrones au cours de l’Holocène, dont une, rompant l’intégralité des décrochements analysés (WBF, EBF et VOLF), s’avère être temporellement corrélée au précédent cluster potentiellement identifié à 3 – 4 ka. Les paramètres présentés dans cette thèse permettent finalement de calculer des paléomagnitudes associées aux deniers évènements sur les WBF et VOLF, comprises entre Mw7,6 et Mw8 selon les paramètres et les lois considérés. Ces données fondamentales pourront servir de canevas de base pour les actuelles études d’aléas sismique en Mongolie. / Western Mongolia experienced an exceptional intracontinental seismic activity during the first part of the XXth century, with four seismic events with Mw > 7.9, located on large strike-slip faults with a transpressive left-lateral motion. This historical seismicity has been qualified as a cluster, and paleoseismological studies indicate that this phenomenon might have occurred 3 to 4 ka ago. Although the eastern part of the Gobi-Altai mountain range broke during the eponymous earthquake in 1957 (Mw8) along the Eastern Bogd Fault (EBF), several fault traces affecting Holocene formations and presenting continuous fault length > 100 km have been documented on the western part of the range, while not been the purpose of quantitative studies.This thesis first gives an inventory of the seismic activity within the Gobi-Altai, through detailed mapping of the fault scarps which show activity during the late-Pleistocene, also with an analyze of the fault segmentation regarding to geometric and kinematic criteria. Secondly, some morphotectonical and paleoseismological surveys have been carried on, along with dating techniques as 10Be in situ, OSL and radiocarbon, which allow quantifying the faults slip-rates since late-Pleistocene, and paleo-ruptures ages during the Holocene period on the two main strike slip faults analyzed : the Valley of Lakes Fault (VOLF), and the Western Bogd Fault (WBF); preliminary data also allow quantifying the slip rates and the age of the most recent event along the Tsogt Fault (TF) and the Tsagaan Gol fault (TGF), respectively located in the westward continuation of the WBF and the VOLF. Moreover, systematic lateral offset measurements have been carried on along the faults traces, and especially along the 1957 surface rupture, where a statistical approach based on probability density has been applied in order to assess the average lateral offsets successively recorded by the landforms, and reveals that the slip distribution along the EBF on 3 seismic cycles. This suggests first that the slip distributions rather even along the rupture, and secondly that the slip related to the 1957 earthquake is of the same amount than previous major ruptures on that fault, for half of the main 1957 surface rupture. This approach of average coseismic slip assessment on several seismic cycles is promising, but still presents some inconsistencies, and we propose developments and validation perspectives.Morphotectonic parameters and ages of paleo-earthquakes reported in this thesis present a strong consistency with former studies along the EBF, with slow slip-rates (0,7 ± 0,2 mm.yr-1 on the WBF and 0,5 ± 0,1 mm.yr-1 on the VOLF), and significant left lateral coseismic offsets (2 – 4 m). Our measurements of slip rates indicate furthermore a distribution of the on-fault deformation from the EBF to the WBF and the VOLF, splitting the deformation in two branches, the north one reaching the TGF and the southern one visibly branches on the TF. Assessments of the paleoearthquakes ages show three possible clusters between the studied faults during the Holocene time, and one of them may break the three faults (WBF, EBF and VOLF), and seems as well to correspond to the cluster possibly identified 3-4 kyr ago. The fault parameters presented in this thesis allow computing the paleomagnitudes related to last seismic events along the WBF and the VOLF, ranging from Mw7.6 and Mw8, depending of parameters and empirical relations considered. Finally those fundamental data may serve as basic pattern for current seismic hazard assessments in Mongolia.

Les observations récentes ont mis en évidence la diversité des régimes de glissement des failles, et particulièrement l'existence de glissements asismiques transitoires, les séismes lents. Ce travail a pour objectif la compréhension de l'impact de ces séismes lents sur le cycle sismique. La zone étudiée correspond à la zone de subduction du sud du Mexique, au niveau de la lacune sismique de Guerrero. A partir de mesures de déplacement de surface, principalement par GPS, le glissement sur l'interface de subduction est modélisé par des dislocations dans un milieu élastique. Cette analyse nous permet de contraindre l'évolution spatio-temporelle de deux épisodes de glissements lent (2006 et 2009-2010), ainsi que le couplage de l'interface de subduction. Nos résultats montrent une certaine variabilité dans l'évolution spatio-tempororelle des deux glissements étudiés : le séisme lent de 2006 présente clairement une propagation du glissement, à une vitesse d'environ 1 km/jour ; le séisme lent de 2009-2010 présente deux sous-évènements, l'occurrence du deuxième sous-évènement étant liée au déclenchement par le séisme de Maule au Chili. Nos résultats mettent également en évidence les variations latérales dans le couplage intersismique de l'interface de subduction : le couplage dans la lacune sismique de Guerrero étant 4 fois plus faible que le couplage de part et d'autre de la lacune. Ainsi la majeure partie du glissement est accommodée par les séismes lents dans la lacune sismique de Guerrero. / Recent observations reveal the existence of different slip behaviors on fault, and among them the occurrence of transient aseismic slip events, the so-called slow slip events (SSEs). The general goal of this work is to understand the impact of slow slip events on the seismic cycle. The area of study is located in the southern Mexican subduction zone, around the Guerrero seismic gap. We use continuous GPS measurements of the ground displacements to model the slip on the subduction interface, using a dislocation model in an elastic half space. We can thus constrain the spatial and temporal evolution of two slow slip events (in 2006 and in 2009-2010), as well as the coupling ratio of the subduction interface. Our results highlight the differences in the spatio-temporal evolution of the two slow slip events : during the 2006 SSE, the slip propagated at a velocity of 1 km/day. The 2009-2010 SSE occurred in two sub-events and the second sub-event was triggered by the surface waves of the Maule earthquake (in Chili). Our results also show the lateral variations in the interseismic coupling of the subduction interface : the coupling ratio in the Guerrero gap being only 1/4 of the couling ratio on both sides of the gap. Most of the slip is thus accommodated by slow slip events in the Guerrero seismic gap.

Le 11 mai 2011, un double séisme a frappé la ville de Lorca, dans le sud-est de l'Espagne. Ce séisme a eu de lourdes conséquences pour la population et provoqué de nombreux dégâts sur les bâtiments de la ville. Les tremblements de terre sont des phénomènes caractérisés par leur relative faible fréquence et par leur violence. Dans le bassin méditerranéen, la grande majorité des séismes ne sont perceptibles qu'à l'aide d'instruments de mesure. Certains événements récents ont néanmoins surpris la communauté scientifique par leur magnitude modérée mais leur intensité élevée. C'est notamment le cas du tremblement de terre de Lorca. Depuis 40 ans, l'administration espagnole essaye de se prémunir davantage contre ce phénomène. Cela se traduit par l'évolution des normes parasismiques appliquées à la construction ou encore par la création de cartes de zonages sismiques intégrant à chaque révision une plus grande portion du territoire exposée au risque. Cependant, ces mesures évoluent au rythme des événements. Or, la localisation des villes en zone sismique et la vulnérabilité des bâtiments préalablement construits sont des variables difficilement ajustables à ces nouvelles normes. Lorsqu'un séisme affecte une ville, l'origine des dommages corporels réside dans la coïncidence temporelle et spatiale entre l'individu et le danger, danger la plupart du temps lié aux bâtiments. Dans ce contexte, nous avons proposé une démarche focalisée sur l'individu. Pendant la période de crise sismique, la mobilité individuelle apparaît comme l'élément qui va permettre la coïncidence spatiale entre le danger et l'individu. A partir d'entretiens qualitatifs, nous avons analysé la mobilité individuelle en période de crise sismique, par la collecte d'informations très fines sur les aspects spatiaux et temporels de celle-là. Parallèlement, nous proposons une étude des principaux dangers observés pour le cas de Lorca et comment ceux-ci peuvent toucher la population. En général, la plupart des décès en cas de séismes se produisent à l'intérieur des bâtiments à cause de leur effondrement. Cependant, dans le cas de séismes d'intensité faible à modérée, c'est la chute d'éléments non structuraux, notamment d'éléments de façade, qui constitue un danger. En croisant spatialement et temporellement la mobilité des individus avec les dangers, nous concluons que la mobilité individuelle est un facteur déterminant de l'exposition humaine. Les facteurs qui vont conditionner la mobilité sont principalement liés à la mise en protection immédiate ainsi qu'à l'organisation en famille préalablement à l'évacuation de la ville. Après avoir permis d'évaluer les dangers, les analyse spatiales ont mis en avant que la grande incertitude de cette évaluation rend difficiles les études très fines. Cependant, elles ouvrent des perspectives intéressantes pour des travaux futurs. Enfin, notre travail de thèse a proposé l'adaptation d'une méthodologie de collecte et d'analyse des dynamiques spatiales et temporelles liées aux individus, à la base créée pour les crues rapides. Cette méthode s'est montrée facilement transposable à de futurs événements, ce qui permettrait de réaliser des études comparatives. / May the eleventh, 2011, a double earthquake struck the town of Lorca in southeastern Spain. The earthquake had serious consequences for the population and caused extensive damage to buildings in the town. Earthquakes are phenomena characterized by their relatively low frequency and their violence. In the Mediterranean Basin, most earthquakes are only perceptible using measuring instruments. Recent seismic events have surprised the scientific community by their moderate magnitude associated with a high intensity. This is particularly the case of Lorca's earthquake. In the last 40 years, the Spanish authorities have tried to develop prevention towards seismic phenomenon. Measures have been taken by improving seismic building codes or by updating the seismic zoning, each revision taking into account wider zones at risk. In general, events give rhythm to the evolution of these measures. However, variables like the location of cities in seismic zones and the vulnerability of existing buildings can hardly be adjusted to these new standards. When an earthquake affects a town, danger is mostly caused by building failure. However, injuries and fatalities are due to the spatial and temporal coincidence between the individual and the danger. In this context, we propose an approach focused on the individual. During a seismic crisis, individual mobility appears to be the element that will allow the spatial coincidence between the danger and the individual. We analyzed individual mobility during seismic crisis by collecting very precise information regarding spatial and temporal aspects of mobility through qualitative interviews. Simultaneously, we carried out a study of dangers observed in the case of Lorca. In general, most of the fatalities caused by earthquakes occur inside collapsed buildings. However, in the case of low-to-moderate-intensity earthquakes, danger is related to the fall of non-structural elements, especially frontage elements. By combining spatial and temporal mobility of individuals with dangers, we conclude that individual mobility is a determinant factor of human exposure. The factors that will determine mobility are mainly related to the immediate protection measures as well as the family organization prior to the town evacuation. Finally, spatial analysis allowed assessing dangers. The associated uncertainty of this evaluation makes it difficult a more accurate study, but offers perspectives for future work. This work proposes the adaptation of a methodology to collect and analyze spatial and temporal dynamics related to individuals. Developed initially for flash flood analysis, this methodology is transferable to future events, enabling comparative studies.

L'objectif général de la thèse est la simulation numérique des mouvements forts du sol dûs aux séismes. Les déformations importantes du sol engendrent des comportements nonlinéaires dans les couches superficielles. L'apport principal de la thèse est la prise en compte de la nonlinéarité des milieux dans un contexte éléments finis Galerkin discontinus. Différentes lois de comportement sont implémentées et analysées. Le cas particulier du modèle élastoplastique de Masing-Prandtl-Ishlinskii-Iwan (MPII) est approfondi. Cette étude est divisée en deux parties. Une première qui vise à poser la structure du problème en présentant les équations et modèles utilisés pour décrire les mouvements du sol. Dans cette partie nous présentons également la méthode d’approximation spatiale Galerkin Discontinue ainsi que les différents schémas temporels que nous avons considérés. Une attention particulière est portée sur la complexité algorithmique du modèle nonlinéaire élastoplastique MPII en vue de réduire le temps de calcul des simulations. La deuxième partie est dédiée aux applications numériques. Ces applications sont réparties en trois catégories distinctes. 1) Nous nous intéressons toutd’abord à la configuration unidimensionnelle où une seule onde de cisaillement est propagée. Dans ce contexte, un flux numérique décentré est établi et des applications aux cas nonlinéaire élastique et nonlinéaire élastoplastique sont étudiées. Une solution analytique concernant le cas nonlinéaire élastique est proposée, ce qui permet de réaliser une étude numérique de convergence. 2) Le problème unidimensionnel étendu aux trois composantes du mouvement est étudié et utilisé comme un premier pas vers le 3D compte tenu du couplage entre les ondes de cisaillement et de compression. Nous nous intéressons ici à des signaux synthétiques et réels. L’application d’une méthode permettant de réduire significativement le temps de calcul du modèle élastoplastique est détaillée. 3) Une configuration tridimensionnelle est examinée. Après différentes applications de vérification en milieu linéaire, deux cas d’étude élastoplastique sont analysés. Une première sur un mode propre d’un cube puis une seconde sur un milieu plus réaliste composé d’un bassin hémisphérique à couches sédimentaires ayant un comportement élastoplastique / The general objective of this thesis is the numerical simulation of strong ground motions due to earthquakes. Significant deformations of the soil generate nonlinear behaviors in the superficial layers. The main contribution of this work is to take into account the nonlinearity of the media in a discontinuous Galerkin finite elements context. Different constitutive laws are implemented and analyzed. The particular case of theMasing-Prandtl-Ishlinskii-Iwan (MPII) elastoplastic model is looked at in-depth. This study is divided into two parts. A first one that aims at defining the framework of the problem by presenting the equations and models used to describe the soil motion. In this part we also present the Galerkin Discontinuous spatial approximation method as well as the different temporal schemes that we considered. Particular attention is paid to the algorithmic complexity of the nonlinear elastoplastic MPII model in order to reduce the computation time of simulations. The second part is dedicated to numerical applications. These applications are divided into three distinct categories. 1) We are first interested in the one-dimensional configuration where a single shear wave is propagated. In this context, an upwind numerical flux is established and applications to nonlinear elastic and nonlinear elastoplastic cases are studied. Ananalytical solution concerning the nonlinear elastic case is proposed, which makes it possible to carry out a numerical study of convergence. 2) The one-dimensional problem extended to the three components of the motion is studied and used as a first step towards 3D applications considering the coupling between the shear and compression waves. We are interested here in synthetic and real input signals. The application of a method that significantly reduces the calculation time of the elastoplastic model

Depuis maintenant une dizaine d'années, la vision du cycle sismique en zone de subduction a beaucoup évolué. Des découvertes récentes ont mis en évidence une grande diversité des régimes de glissement dans ces zones, avec notamment des glissements asismiques transitoires appelés « séismes lents » (SSE) et des vibrations de faibles amplitudes, persistantes dans le temps, appelées « trémors non volcaniques » (NVT). Ce travail a pour objectif l'étude des trémors non volcaniques afin de caractériser ces nouvelles manifestations des zones de faille. Nous avons abordé ce problème avec deux approches distinctes :1. Observer les trémors dans le milieu naturel afin de déterminer leurs caractéristiques. La zone étudiée correspond à la lacune sismique de Guerrero le long de la subduction mexicaine. Nous avons développé une méthode de détection et de localisation des NVT au Mexique grâce à des analyses d'antennes par formation de voie sur les corrélations. Cette méthode permet de mettre en évidence cer taines caractéristiques des NVT : une complexité des sources pour un épisode de trémors, une corrélation entre les activités de NVT et les pics de vitesse des glissements lents à plus long terme. Par ailleurs, l'étude de l'impact du séisme de Maule (2010, Chili, Mw 8.8) au Mexique montre qu'il a déclenché le second sous évènement du séisme lent de 2009-2010. Ce déclenchement d'un SSE s'ac- compagne de fortes activités de trémors, modulées par les ondes du séisme de Maule dans un premier temps, puis simplement associées au SSE.2. Modéliser les trémors expérimentalement et numériquement pour mieux com- prendre leur origine physique et leurs évolutions sur le long terme. Nous avons en particulier utilisé une expérience de frottement à faible vitesse qui indique une corrélation systématique entre les accélérations d'un glissement et l'émission de signaux qui ressemblent à des NVT. Une modélisation numérique de la zone de subduction mexicaine est également présentée et montre la possibilité de reproduire des trémors en considérant une transition d'affaiblissement critique associée à un processus de décrochage. / The vision of the seismic cycle in subduction zones has considerably evolved over the last 10 years. New discoveries has pointed the diversity of slip behaviors in these zones with aseismic slow slip called « slow slip events » (SSE) and persistent low amplitudes vibrations called « non-volcanic tremors » (NVT). The goal of this thesis is to study the non-volcanic tremors in order to characterize these new manifestations of fault zones. We used two different approaches: 1. We first observed the non-volcanic tremors in the nature in order to characterize this phenomenon. The area of interest is the Guerrero seismic gap along the Mexican subduction zone. We develop a new detection and location method based on beamforming of correlations of seismic signals. This new method exhibits some characteristics of NVT: a complex source for a single tremor episode and a correlation between the NVT episodes and the long-term peak of movement velocity in southwards direction. Moreover, the study of the consequences of the Maule earthquake on the Mexican subduction zone showed that this earthquake triggered the 2009-2010 SSE in Guerrero. This triggering of slow slip is accompanied by strong seismic tremor actvity that are first modulated by the passing waves and then associated to the SSE. 2. We model numerically and experimentally the tremors in order to better understand their physical origin and their long-term evolution. We used a very slow friction experiment that indicates a systematic correlation between slip acceleration of a slider and emission of acoustic signals that are similar to NVT. A numerical modeling of the Mexican subduction zone is also presented and shows the possibility to reproduce NVT with a critical depinning transition.

Non renseigné / Most mega-earthquakes (i.e. earthquakes with Mw ≥ 8.5) occur along subduction mega-thrusts, the interfaces between the subducting - and the overriding plates in convergent margins. These events may have catastrophic impact on human societies due to their destructive potential. For this reason being able to predict the timing and size of these earthquakes became one goal of the international scientific community. The subduction seismic cycle is influenced by many different parameters. The interplay between these parameters governing the frequency and size of megathrust earthquakes still remains unclear, mainly due to the short (i.e. limited to the last century) seismic record.The seismogenic part of the subduction thrust fault spans between depths of 11±4 and ± 51 km (Heuret et al. 2011). In this zone a combination of temperature, pressure and rocks characteristics creates conditions favourable for seismic behaviour. Whether a specific area in the subduction thrust fault has the ability to trigger mega-earthquakes can be expressed using the degree of seismic coupling, i.e. the amount of slip that occurs with respect to the total amount of plate convergence (e.g. Scholz 1998; Scholz & Campos 2012). When a fault is fully coupled, all of the fault slip occurs during earthquakes instead of also during aseismic behaviour (e.g. slow slip events). The internal structure of the interplate fault zone mainly determines whether an area within a subduction zone behaves seismic or aseismic (Wang & Bilek 2011). This is influenced by the topography of the plate interface (e.g. subducting seamounts; Wang & Bilek 2014), but also subducted sediments and fluids in the subduction channel may play an important role.The main goal of this project is to understand which parameters affect the behaviour of mega-earthquake ruptures. This will be done by comparing natural data (e.g. seafloor roughness, sediment thickness and fluid content in the subduction channel) to rupture characteristics of major recent earthquakes. With this analysis also more knowledge can be gained on the triggering of slow earthquakes instead of mega-earthquakes. These are slow slip events with lower frequencies and longer durations than ‘regular’ earthquakes (Saffer & Wallace 2015).The database of natural data, implemented by the long-term scientific joint venture between the Univ. Montpellier and the LET (Roma Tre) will be used for the analysis. Ongoing work is done on determining a method for estimating the seafloor roughness, i.e. the distribution of high, low and smooth areas (by Michel Peyret in collaboration with Serge Lallemand, Univ. Montpellier). Also data is available on the trench sediment thickness around the world (Heuret et al. 2011). In the frame of this project, information on the roughness of the seafloor will be added to the database. In addition the rupture characteristics of major recent earthquakes will be collected. By performing a multiparametric statistical analysis of the database, a conceptual model will be realized, exploring the possible link between all the different parameters. The aim is to validate this model in the lab using scaled 3D analogue models. This will be done both at the LET and at Univ. Montpellier by using a broad range of geometries and contact materials with different rheologies (e.g. gelatin, foam rubber and a new analogue material; Caniven et al. 2015; Corbi et al. 2013). This jointed experimental approach with both the Univ. Montpellier and the LET involved creates a rich environment where differences and similarities of the two different approaches can be used to validate the results.

linéarités de comportement des sols Cette thèse traite de la réponse locale de bassins sédimentaires à des mouvements sismiques faible et fort en utilisant la modélisation numérique (SEM et FEM). A faible niveau de distorsion, des configurations théoriques 3D et 2D sont prises en compte pour l'étude des effets de site linéaires: dans le cas de bassins 3D à symétrie de révolution, on identifie les conditions d'existence de la «résonance 3D» en fonction des paramètres géométriques et rhéologiques. Dans le cas 2D, la réponse sismique de deux vallées alpines est modélisée: les résultats sont cohérents avec les phénomènes observés dans la réalité. A fort niveau de distorsion, la formulation théorique 3D d'un modèle viscoélastique non linéaire (VNL) est proposée. Son implémentation par la méthode des éléments finis dans le cas 1 D permet de l'appliquer à un cas réel (Kushiro, 1993). Les résultats obtenus sont comparés à ceux d'autres modèles (linéaire équivalent et élastoplastique) et à l'enregistrement réel de surface.

Le gisement de gaz de Lacq est exploité depuis 1957 et la sismicité induite a débuté en 1969. De magnitude modérée (M <= 4.4), cette sismicité est induite par l'extraction des fluides. Après nettoyage de la base de données disponible (1989-1997), un Pointé Automatique dont les paramètres sont optimisés par Recuit Simulé a permis la mise à jour de la base de données. Ajoutés aux données pointées entre 1974 et 1988 et grâce à un modèle de vitesse 3D, plus 1600 évènements ont été localisés, la plupart dans et au dessus du réservoir. La confrontation des localisations des hypocentres et des mécanismes aux foyers avec des informations apportées par les forages et la sismique 3D a permis de localiser les différents essaims de séismes au niveau de failles relativement importantes et orientées préférentiellement pour accommoder la perturbation de contrainte. Alors que le gisement présente une orientation WNW-ESE, la sismicité est localisée selon<br />un axe préférentiel SSW-NNE sur les failles principales d'orientation WNW-ESE. Les zones WNW et ESE du gisement sont quant à elles déformées de manière asismique, comme le montre les profils de nivellement. La confrontation de la sismicité avec l'historique et le modèle de réservoir a permis de mettre en avant (i) un essaim déclenché par une injection d'eau involontaire, (ii) une connexion importante entre le centre du réservoir et les réserves de gaz au nord-est et la déplétion associée, (ii) ainsi qu'une réinterprétation du modèle de réservoir dans certaines zones mal connues. La seconde voie d'étude a été portée sur l'analyse des formes d'ondes par décomposition atomique en chirplets, généralisation de la transformée en ondelettes et de la transformée de Gabor. Après une généralisation à 7D des paramètres de construction des chirplets pour une adaptation aux ondes sismologiques et l'utilisation d'algorithme de séparation d'ondes, il a été possible de décomposer les sismogrammes en un nombre limité d'atomes. Chacun des atomes est qualifié de 9-tuplet, caractérisé par des attributs morphologiques : (i) un temps d'arrivée, (ii) un fréquence centrale, (iii) une énergie, (iv) une phase, (v) une durée, (vi) un taux de modulation de fréquence, (vii) un type de modulation de fréquence, (viii) une asymétrie d'enveloppe et (ix) une plateau constant de l'enveloppe. Grâce à cette caractérisation poussée des ondes, il a été possible de pratiquer des classifications des sismogrammes. La classification et la représentation en dendrogramme ont montré une différenciation nette entre les séismes déclenchés par surpression de fluide des autres, révélant que la source des évènements déclenchés est plus courte et impulsive. En se basant uniquement sur les attributs morphologiques et en filtrant la composante de temps et fréquence des arrivées des ondes, nous avons montré qu'il était possible de comparer des séismes de localisations différentes.

L'objectif que nous avons essayé d'atteindre dans le cadre de cette thèse est l 'évaluation de l'aléa sismique au Maroc par l'application de différentes étapes. La première partie de ce travail consiste en l'élaboration d'un fichier informatisé (chapitre 1) de la sismicité de la région comprise entre 21°N-38°N et O°W-20°W à partir du début de ce siècle. Un peu plus de quatre mille six cent événements ont été recensés faisant l objet d'une révision complète, 34% des événements ont été relocalisés. Nous avons effectué un nouveau dépouillement des archivés au D.P.G., ce qui nous a permis de traçer une cartes macrosismiques des séismes largement ressentis. Nous avons essayé d'utiliser ce fichier pour calculer l'aléa sismique à l'aide du programme de calcul EQRISK développé par McGuire (première partie: chapitre 2). Les résultats obtenus sont donnés sous forme de cartes des valeurs d'intensités susceptibles d'être ressenties pour les deux valeurs de probabilité annuelle: 0.01 et 0.005. Dans le premier chapitre de la deuxième partie de cette tIhèse, nous avons procédé à une révision complète des données macrosismiques et instrumentales du séisme d'Agadir, ce qui nous a permis de dresser une nouvelle carte macrosismique et de mieux préciser le mécanisme au foyer. Une étude de l'effet de site a été effectuée pour mettre en évidence l'influence des couches sédimentaires, sous la ville d'Agadir, sur les constructions pendant le tremblement de terre. Au dernier chapitre, de la deuxième partie, une étude sismotectonique détaillée a été réalisée dans la province d'Al Hoceima dans le cadre de l'Action Inteqrée 359/80. Cette région est l'une des plus sismiques au Maroc, de nombreuses secousses y sont ressenties chaque année . Une compagne de microsismicité a été effectuée à proximité de l'accident de Nekor. Un réseau, composé de 16 stations sismiques portables et de 9 stations numériques à 3 composantes télécommandées, a été installé pour une durée de cinq semaines . L'objectif de cette expérience était de localiser précisément la sismicité, de calculer des mécanismes à la source, d'entreprendre une étude du mouvement du sol en réponse à une sollicitation sismique ..

Nous étudions ici différents aspects de la sismicité du Sud-Est de la France, en nous aidant des données nouvelles disponibles aujourd'hui. L'apparition récente de nouveaux réseaux sismologiques (L.O.G., S.E.O., I.G.G.) dans les Alpes, et la multiplication d'études géophysiques (grands profils, aéromagnétisme,paléomagnetisme,photos satellites, nivellement) et géologiques devraient aboutir dans les prochaines années à une synthèse géodynamique des Alpes. Notre travail apporte certains éléments Sismologiques pour cette future synthèse. Nous présentons d'autre part une étude sur une nouvelle méthode de zonage sismique que nous appliquons en particulier au ca's du S.E. de la France. Les trois premiers chapitres concernent surtout l'acquisition ou le traitement de nouvelles données de sismologie. Dans le Chapitre l, on constitue un fichier de séismes historiques et récents, avec lequel on trace des cartes de sismicité et on étudie la loi magnitude-fréquence. Le Chapitre II livre des solutions focales nouvelles. Le Chapitre III donne les résultats d'une campagne d'enregistrement de petits séismes en Ubaye-Queyras en septembre-octobre 1977 . Le Chapitre IV est consacré à la synthèse des données acquises dans les trois premiers chapitres et à leur comparaison avec d'autres données (géophysiques et géologiques). Enfin, le Chapitre V expose notre méthode de zonation sismique. Nous présentons en annexe le fichier de séismes, les cartes de sismicité, ainsi qu'un exposé sur la loi magnitude-fréquence.