La paléosismologie, spécifiquement les tranchées paléosismiques, a été l’outil de recherche conventionnel par excellence et plus amplement utilisées dans les études focalisées sur l’estimation de l’aléa sismique. Plus récemment, l’intérêt sur l’activité sismique tout au long des systèmes de failles associées aux limites de plaques, a conduit à la caractérisation du remplissage sédimentaire des bassins en environnement côtier avec un fort potentiel de préservation et donc d’enregistrement des déformations et perturbations liées à une activité tectonique.Le travail présenté dans ce mémoire a pour objective de développer une analyse du risque naturel, particulièrement le risque de séismes et tsunamis dans la bordure nord de Venezuela, sur la trace sous-marine active de la faille décrochante d’El Pilar.Dans ce contexte la texture, minéralogie et géochimie de la séquence sédimentaire d’un petite appendice de la fosse de Cariaco appelé golfe de Cariaco, a été étudié à partir des différentes techniques (analyse de taille de particules, fluorescence X semi-quantitative, mesures de perte à feu, Pyrolise Rock-Eval et propriétés magnétiques). Comme résultat d’une caractérisation détaillée des sédiments et d’un control chronologique détaillé basé sur la datation radiocarbonique et la correction de l’effet réservoir (potentiellement accumulé sur les différents fractions organiques datées), nous avons mis en évidence un particulier dépôt silto-argileux, fortement détritique, d’une ample distribution spatiale, dont l’âge le plus probable correspond au séisme historique de 1853, qui a affecté la ville de Cumaná. Les analyses géochimiques révèlent que ce dépôt de texture homogène précède une augmentation relative du Br (cts/s) à la surface. On interprète la conjugaison de ces deux singularités comme la conséquence d’un glissement sous-marin et d’une modification morphologique majeure concomitante dans le seuil qui connecte la fosse et le golfe de Cariaco. Dans cette chaine d’événements, ce dépôt serait le produit de la décantation du sédiment en suspension remobilisé et l’augmentation du Br est la manifestation d’un changement hydrogéochimiquemarqué pour une majeure influence marine sur le golfe. Pour établir un point de référence pour l’évaluation du potentiel et de la qualité d’enregistrement de l’activité sismique dans le golfe de Cariaco, un deuxième site lié à la plaque Caraïbe fait partie de cette recherche. Il s’agit de l’arc des Petites Antilles, siège d’une activité sismique et volcanique importante, résultat d’une subduction. Les carottes ici prélevées et étudiées avec la même approche sédimentologique, ont enregistré la trace d’un événement sismique récent, plus probablement le séisme de Redonde de 1974. Cet événement correspond à un dépôt silto-sableux, non stratifié et acoustiquement transparent. Les profils sismiques montrent que la relation géométrique de faciès similaires (de faible énergie) avec l’escarpement co-sismique, implique une origine séisme-induite, a différence des turbidités, dont leur origine peut être influencée par la superposition complexe de plusieurs mécanismes déclencheurs. / Paleoseismological researches, especially through paleoseismic trenches, have been widely developed during the last decades, and are presently commonly used for seismic hazards assessment. Beside, and more recently, the recording of earthquakes and related phenomena (gravity reworking, tsunamis) by sedimentary fills became a precious complement. Along plate boundaries and associated major active faults, the search and characterization of sediments with a strong potential for preservation of earthquakes-disturbances or specific layers has been developed; the here-presented results belong to this approach.The research presented in this dissertation aims to develop an analysis of natural risk, specifically tsunamis and seismic risk, on northern boundary of Venezuela, along the submarine active trace of the El Pilar strike slip fault.In this context, the texture, mineralogy and the geochemistry of the sedimentary sequence in the Gulf of Cariaco (built upon the fault) have been studied from different techniques (particle-size analysis, semi-quantitative XRF analysis, loss on ignition, Rock-Eval Pyrolysis and magnetic properties). From detailed characterization of sediments and an extensive chronological control, based on radiocarbon dating (including a local reservoir effect estimation), we found evidences of a particular siliciclastic fine-grained deposit, widely distributed, which probably corresponds to the 1853 major catastrophic earthquake and associated tsunamis which affected the Cumana town. The geochemical analyses reveal that this deposit of homogeneous texture is followed by an abrupt relative increase of Br at the surface. We interpret this singular combination as a consequence of a submarine landslides and subsequent morphological modification at the connection between the Cariaco trench and the Gulf of Cariaco. This deposit can be the result of the settling of remobilized suspended sediment and the Br increase is the indicator of a hydrogeochemical change indicated by important marine influence on the gulf.In order to establish a point of reference for the evaluation of the potential sedimentary record of the seismic activity in the Gulf of Cariaco, a second site related to the Caribbean Plate boundaries is analysed and discussed in this work. Sedimentation within the Lesser Antilles Arc was chosen as offering a comparable seismic activity but a strongly different sedimentary environment. In this zone of seismic and volcanic activity, consequence of a subduction, we collected and analysed cores with the same sedimentological approach, searching and characterizing recent earthquake record. The most recent one probably corresponds to the 1974 Redonda event. High resolution seismic profiles across active faults indicate precise direct geometrical relationship of specific turbidite-type layers with scarps genesis.

Les phénomènes d’instabilité générés par les séismes dans les pentes représentent un risque naturel majeur à l’échelle mondiale. Ce risque a, de plus, tendance à croître du fait d’une urbanisation croissante dans des zones à fort aléa, comme en témoignent de nombreux évènements catastrophique à travers le passé. Pour garantir la sécurité des personnes et des biens face à ce type d’aléa, il faut améliorer les prédictions des déplacements générés par les séismes dans les pentes.Les premiers travaux scientifiques visant à corréler les caractéristiques des instabilités dans les pentes aux paramètres sismiques sont de nature empirique et remontent aux années 1980. Ces méthodes ne permettent pas d’expliquer tous les mouvements de terrain observés à travers le monde, notamment lorsque les effets de site modifient sensiblement la distribution des mouvements dans les pentes et génèrent des déplacements plus grands ou plus petits que ceux prédits par les lois empiriques qui négligent les effets de site.La première partie de cette thèse présente une nouvelle base de données mondiale de 277 glissements de terrain d’origine sismique ou non. L’analyse statistique conduite sur les données de cette base a montré que la forme géométrique moyenne des glissements de terrain reste stable lorsque le volume des instabilités augmente.La deuxième partie de la thèse est dédiée à l’analyse de la stabilité du glissement de terrain de Diezma (Espagne) au moyen de trois méthodes : la Méthode de NEWMARK (1965) basée sur le principe de l’équilibre limite en 2D qui suppose des conditions elasto-parfaitement-plastiques ; l’analyse par différences finies dans le domaine temporel avec le code FLAC en 2D et en 3D qui suppose des conditions visco-elasto-parfaitement-plastiques ; et l’analyse par recombinaison modale dans le domaine fréquentiel avec le code CESAR-LCPC en 2D et en 3D qui suppose des conditions visco-élastiques.Situé dans une région sismiquement active, le glissement de terrain de Diezma est susceptible d’être affecté par des séismes. Une analyse des déplacements induits dans ce versant par onze signaux sismiques différents a été menée dans l’objectif d’établir des corrélations entre les paramètres macro-sismiques des scenarios sismiques appliqués et les déformations calculées par les trois méthodes.Les résultats montrent des différences marquées à la fois qualitatives et quantitatives. Les conclusions principales sont : (a) Les déplacements obtenus par les trois méthodes ont des ordres de grandeur différents et peuvent donc conduire à une surestimation ou à une sous-estimation des déplacements ; (b) Les deux méthodes numériques montrent que les déplacements sont fortement conditionnés par la géométrie des modèles ; seule la méthode par différences finies fait apparaître une réponse en termes de déplacements qui dépend du scénario sismique considéré; (c) Les résultats des simulations2D et 3D ne sont pas comparables ; des analyses complémentaires doivent encore être menées pour guider l’utilisateur dans le choix de la méthode la plus appropriée; (d) La Méthode de NEWMARK (1965), dont l’utilisation reste très répandue de nos jours, est tout à fait appropriée à l’étude de cas simples mais elle peut se révéler inexacte lorsque la structure géologique / topographique du versant conduit à un fort effet de site car ce dernier n’est pas pris en compte par cette méthode / Seismically induced slope deformation is a worldwide common phenomenon that poses an increasing and considerable threat to fast expanding urbanization, and a great number of catastrophic events throughout the past attest thereof. For this reason, displacement predictions allowing for proper slope surveillance became a major concern.Early attempts to relate slope failures to seismic parameters are of empirical nature and date back to the 1980s. Although having proven stable, these relations are frequently disturbed by site effects causing outliers in terms of smaller or greater displacements than expected.The first part of this thesis presents a newly build chronological database of 277 globally distributed seismically and non-seismically induced landslides. A comprehensive statistical analysis was conducted on the data of this database with the main result that – statistically seen – the average geometrical shapes of landslides differ only proportionally.The second part of the thesis is dedicated to a comparative slope stability analysis of the Diezma Landslide (Spain) by means of three methods: the limit-equilibrium based NEWMARK-Method (1965) in 2D under elasto-perfectly-plastic conditions; finite-difference analysis in the time domain with the code FLAC in 2D and 3D under visco-elasto-perfectly-plastic conditions; and modal recombination analysis in the frequency domain with the finite-element code CESAR-LCPC in 2D and 3D under visco-elastic conditions.Located in a seismically active region, the Diezma Landslide is likely to be affected by earthquakes and was therefore considered as a suitable model case. A broad analysis of expected displacements was conducted using eleven strong-motion seismic scenarios. Moreover, the thesis searches for potential relations between macro-seismic parameters of the applied earthquake scenarios and the predicted deformations obtained from the three methods.It appeared that – due to the functionality of the methods – results differ quantitatively as well as qualitatively, and so does their suitability. Major findings are: (a) Results from the three methods are of different orders of magnitude and, thus, can easily lead to over- or under-estimations of displacements; (b) Both numerical methods reveal a strong influence of the model-geometry on the predicted displacements, whereas a scenario-dependent slope behavior manifested itself only within the finite difference analysis; (c) The switch from 2D to 3D does not necessarily result in a similar performance in each dimension and results must be critically judged before further use; (d) The NEWMARK-Method (1965) has proven itself once more to be appropriate for first slope assessments but not for sophisticated evaluations of ground response to seismic shaking

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

Nous présentons une analyse des variations du taux de sismicité liées aux doublets de forts séismes. Cette étude se concentre sur la zone des répliques du 1er séisme, et cherche à estimer l'impact du 2nd choc principal sur cette activité. Nous nous sommes penchés sur deux cas de doublets de séismes : la séquence turque des séismes d'Izmit et de Düzce en 1999, et le doublet sud-islandais de juin 2000. Ces deux exemples nous apprennent que la réactivation de sismicité induite par le 2nd séisme d'un doublet n'est pas forcément corrélée avec les régions d'augmentation de la contrainte statique, et qu'elle s'établit avec prédilection sur les sites d'une activité géothermale. Sa décroissance temporelle n'est pas systématiquement conforme à la loi d'Omori, puisque des apaisements de ce taux peuvent apparaître avec un délai de quelques jours à quelques mois. Nous avons également investigué les propriétés des doublets de séismes à l'échelle globale, et nous montrons que ceux-ci peuvent être simplement décrits par les lois habituelles gouvernant la sismicité, la loi d'Omori-Utsu et la loi de Gutenberg-Richter. Enfin, contrairement aux résultats obtenus à partir de catalogues régionaux, il semblerait que la production moyenne de répliques consécutive aux 2nds séismes des doublets soit équivalente à la superposition des productions de deux séismes isolés de magnitudes similaires. Cette dernière remarque mérite toutefois d'être étayée par une analyse plus fine des catalogues régionaux, au vu du manque de résolution d'un tel catalogue global pour les séismes de petite magnitude.

Lors d'un tremblement de terre, les transformations physiques et chimiques qui surviennent le long d'une zone de glissement vont conduire à l'altération et à la formation de minéraux. La gouge contient des minéraux magnétiques, qui peuvent être formés sous l'action combinée de la chaleur frictionnelle et par l'action chimique en relation avec les fluides. Ainsi, la gouge a la capacité de se comporter comme un enregistreur magnétique pendant un tremblement de terre. Il s'agit là d'une nouvelle méthode pour identifier les zones de tremblements de terre de glissement. En outre,les minéraux magnétiques initiaux, altérés, et néoformés peuvent être utilisés comme traceurs de certains processus physico chimique.Dans cette étude, nous étudions le magnétisme des roches et l'enregistrement paléomagnétique de la gouge de faille active de Chelungpu qui héberge, entre autres,la zone de glissement principale du tremblement de terre de Chi-Chi (Mw 7.6, 1999,Taiwan). Nous avons bénéficié pour cette étude d'échantillons non altérés provenant des carottes du forage B Taiwan Chelungpu-fault Drilling Project (TCDP). Nous avons également échantillonné la faille de Chelungpu à l'affleurement pour une caractérisation des nanoparticules. Cette caractérisation vise à estimer l'énergie de fracture dans la gouge de faille. / During an earthquake, the physical and the chemical transformations along a slip zone lead to alteration and formation of minerals within the gouge layer of a mature fault zone. The gouge contains magnetic minerals, which could be formed under the combined action of friction heat and fluid. Thus, gouge has the capacity to behave as a magnetic recorder during an earthquake. This constitutes a conceivable way to identify earthquakes slip zones. Besides, altered and neoformed magnetic minerals can be used as tracers of some earthquake processes. In this study, we investigate the rock magnetism and paleomagnetism of the Chelungpu fault gouge that hosts the principal slip zone of the Chi-Chi earthquake (Mw 7.6, 1999, Taiwan) using Taiwan Chelungpu-fault Drilling Project (TCDP) hole-B core samples. We also took a Chelungpu fault outcrop sample for identification of nanoparticle, which associated with fracture energy estimation in fault gouge.

Les interactions entre les processus tectoniques et l’érosion ont été peu étudiées à des échelles de temps courtes (< 1000 ans). Cependant, les séismes peuvent activement contribuer à l’érosion des chaînes de montagne en déclenchant de nombreux glissements de terrain. Des études récentes ont également montré que ces grands événements érosifs pourraient engendrer des changements de contraintes suffisants à proximité des failles actives pour modifier la sismicité régionale. Dans cette thèse, cette problématique a été abordée via une approche numérique. Dans un premier temps, le développement d’un modèle simple de glissements de terrain prenant en compte la topographie des versants a permis de démontrer le rôle des paramètres mécaniques (cohésion et friction), et de la forme des versants sur la distribution de taille des glissements de terrain. Ce modèle a été validé à l’aide de cas naturels de glissements de terrain co-sismiques. Dans un deuxième temps, le rôle de la forme finie des versants sur la probabilité de grands glissements de terrain a été démontrée en se basant sur des données. Enfin, dans un troisième temps, le potentiel effet d’un grand évènement érosif sur la sismicité a été exploré à l’aide d’un modèle numérique de cycle sismique dans lequel ont été implémentées des variations temporelles de la contrainte normale sur la faille. Les résultats mettent en évidence le rôle du volume de sédiments, mais aussi de leur temps d’export. En particulier, les paysages caractérisés par une hauteur unstable des versants importante pourraient, en favorisant de grands glissements de terrain, induire une érosion assez importante et rapide pour modifier de façon significative la sismicité régionale. / Interactions between tectonic processes and erosion have been poorly investigated at short time-scales (<1000 years). However, earthquakes can largely contribute to the erosion of mountain belts by triggering widespread landsliding. Moreover, recent studies have shown that such large erosional events could induce stress changes in the fault environment efficient enough to influence regional seismicity. In this thesis, this problematic is tackled through a numerical approach. Firstly, the development of a simple mechanical model accounting for the complexity and variability of natural hillslopes allowed to demonstrate the role of mechanical parameters (cohesion and friction), and of hillslope shape in the probability density function of landslide sizes. This model has been validated using natural cases of co-seismic landsliding. Secondly, the role of unstable hillslope height on large landslide probability has been demonstrated based on natural data, and the exponential distribution of this unstable height has ben shown. Finally, the potential effect of a large erosional event on seismicity has been explored with a numerical model of seismic cycle, in which has been implemented temporal normal stress variations. The results emphasize the role of eroded sediment volume, but also of the export time of sediments away from the mountain belt. In landscape with high unstable hillslopes, large landslides are favored and in turn, could induce fast an important enough erosion to modify regional seismicity.

L’explosion du nombre et de la variété des données géodésiques, sismologiques et tsunami disponibles est une opportunité exceptionnelle pour produire de nouveaux modèles de la source sismique. Mais ces données n’apportent pas toutes la même information et sont soumises à différentes sources d’incertitudes, rendant la solution au problème inverse non-unique. Dans cette thèse, nous utilisons une méthode d’échantillonnage bayésien pour produire de nouveaux modèles de glissement moins assujettis au sur-ajustement des données et permettant une estimation réaliste de l’incertitude associée aux paramètres estimés. Nous l’appliquons à l’étude du glissement dans trois contextes tectoniques différents : le séisme de Landers (1992, Mw=7.3), la zone de subduction équato-colombienne où s’est produit le séisme de Pedernales (2016, Mw=7.8), et le séisme intra-plaque de Tehuantepec (2017, Mw=8.2). À travers ce travail, nous démontrons l’importance de la considération rigoureuse des incertitudes et les atouts de l’approche bayésienne pour l’étude des différentes phases du cycle sismique. / The explosion in the amount and variety of available geodetic, tsunami, and seismological observations offers an outstanding opportunity to develop new seismic source models. But these data are sensitive to different sources of uncertainty and provide heterogeneous information, which makes the solution of the inverse problem non-unique.In this thesis, we use a Bayesian sampling method to propose new slip models, which benefit from an objective weighting of the various datasets by combining observational and modelling errors. These models are less affected by data overfit and allow a realistic assessment of posterior uncertainties. We apply this method to the study of slip processes occurring in three different tectonic contexts: the Landers earthquake (1992, Mw=7.3), the Ecuador-Colombia subduction zone which hosted the Pedernales earthquake (2016, Mw=7.8), and the intraslab Tehuantepec earthquake (2017, Mw=8.2). Through these analyses, we demonstrate how the study of the seismic cycle can benefit from rigorous uncertainty estimates and Bayesian sampling.

Cette thèse s'inscrit dans le cadre de la recherche ECOLEADER des programmes européens d'accès aux grades installations ayant comme objectif l'amélioration de nos connaissances sur la vulnérabilité des structures à voiles porteurs sous chargement sismique. Une stratégie de modélisation simplifiée basée sur des éléments de poutre multifibre Euler - Bernoulli ou Timoshenko est présentée. Les lois utilisées pour le béton et l'acier sont basées sur la mécanique de l'endommagement et la plasticité respectivement. Malgré sa simplicité, le modèle permet de modéliser de façon satisfaisante le comportement global des structures en béton armé sous chargement dynamique. De plus, il est capable de reproduire qualitativement les tendances du schéma de fissuration et la position des zones d'endommagement. Les calculs non linéaires dynamiques sont les seuls capables d'évaluer les marges de variation de l'effort normal, d'étudier l'état des sections vis-à-vis de l'état ultime et de fournir ainsi les informations nécessaires dans la phase de conception. Outil disponible pour l'ingénierie, cette approche doit être considérée comme pouvant être utilisée en parallèle avec les méthodes classiques de dimensionnement et elle est en mesure de conforter les décisions des concepteurs.

On s'intéresse à la rupture associée à deux classes de phénomènes naturels, les séismes et les instabilités gravitaires. Pour les séismes, on étudie un modèle stochastique de sismicité, basé sur les deux lois les mieux établies pour la sismicité, la décroissance en loi de puissance du taux de sismicité après un séisme, et la distribution en loi de puissance des énergies des séismes. Dans ce modèle, on suppose que chaque séisme déclenche d'autres séismes, dont le nombre augmente avec l'énergie du choc principal. Le taux de sismicité global résulte de la cascade de déclenchements de séismes directs et indirects. On analyse l'organisation spatiale et temporelle de la sismicité dans les différents régimes sous- et sur-critiques du modèle. Ce modèle permet de reproduire un grand nombre de propriétés de l'activité sismique, telles que la variabilité de la décroissance des séquences d'aftershocks, l'augmentation de l'activité sismique avant un séisme, la diffusion des aftershocks, la migration des foreshocks et la modification de la distribution des magnitudes avant un séisme. On obtient avec ce modèle une bonne prédictabilité d'une fraction des séismes qui sont déclenchés à court terme après un grand séisme. Nos résultats démontrent le rôle essentiel des cascades de déclenchement de séismes a toutes les échelles dans l'organisation de l'activité sismique. Concernant l'étude des instabilités gravitaires, une étude statistique de plusieurs catalogues d'éboulements rocheux montre que la distribution des volumes de roches suit une loi de puissance. On propose que cette distribution en loi de puissance résulte soit de l'hétérogénéité initiale de la matrice rocheuse, soit de la dynamique d'un système critique auto-organisé. Certains glissements de terrains sont précédés par une accélération de la vitesse de glissement avant la rupture finale. On peut reproduire l'évolution temporelle du glissement a l'aide d'un modèle de bloc rigide avec une loi de friction dépendante de la vitesse de glissement et de l'état de contact entre le bloc et sa surface de glissement. L'analyse de deux glissements de terrains avec ce modèle permet de distinguer une accélération du glissement dans le régime stable, d'une accélération instable qui évolue vers une rupture catastrophique.

Nous analysons les séries temporelles des glissements de terrain de Nouvelle-Zélande (NZ) en temps et en taux et mettons en évidence une corrélation dans les occurrences de glissements. Cette corrélation n'est pas due aux interactions glissement-séisme ou glissement-glissement mais aux interactions glissement-climat. Nous comparons la dynamique des glissements en temps, espace et taux pour la NZ, le Yosemite (Californie, Etats-Unis), Grenoble (Isère), Val d'Arly (Haute-Savoie), l'Australie et le Wollongong (New South Wales, Australie). Les taux journaliers de glissements de la NZ, du Yosemite, de l'Australie et du Wollongong acceptent une loi puissance pour des taux variant de 1 à 1000 glissements/jour. Cela suggère que les mêmes méchanismes sont à l'oeuvre pour le déclenchement de plusieurs centaines de glissements comme de un seul glissement. L'analyse jointe de ces six catalogues nous a permis de dériver des paramètres permettant de classer la dynamique de chaque endroit en terme de glissements. Enfin, nous comparons les distributions en espace des répliques sismiques et des glissements de terrain déclenchés par les séismes de Chi-Chi (Mw7.6 - Taiwan), du Kashmir (Mw7.6 - Pakistan), de Fiordland (Mw7.2 - NZ), de Northridge (Mw6.6 - Californie) et de Rotoehu (MW5.6 - NZ). Les répliques sismiques et les glissements présentent des distributions spatiales similaires. Nous ne trouvons pas de réponse linéaire entre les glissements et/ou les répliques et les observations de mouvements du sol. Nous suggérons que les glissements et les répliques sont contrôlés par les mêmes méchanismes et donc qu'il existe un rôle de la contrainte statique sur le déclenchement des glissements de terrain.