Search

Considerable effort has been made in the last two decades regarding the design and implementation of tools aimed at fast characterization of earthquake sources. As a fruitful harvest, a preliminary magnitude estimation is often available within fifteen minutes for a Mw=6.0 earthquake occurring anywhere in the world. But until recently, it still took several hours to determine the first order attributes of a great earthquake (e.g. Mw≥8.0) even in a well instrumented region. Our goal here is to develop a method providing fast centroid moment tensor (CMT) estimates for large earthquakes based on the W phase.The W phase is a very long period phase (100-1000 s) starting at the same time as the P-wave. It is conspicuous on broadband displacement records before the surface waves arrivals and can then be used for fast characterization of the source. We show that the W phase is not affected by shallow heterogeneities since a significant part of its energy is propagating deep into the mantle. Because of its long period nature, the W phase algorithm is particularly well suited for robust determination of CMT parameters for Mw≥7.5 earthquakes. For large earthquakes indeed, the source complexity is more evident which translates into significant variations of narrow-band magnitude estimations performed at short period. Tsunami earthquakes and outer-rise events provide good examples of this problem (Kanamori, 1972 ; Lay et al., 2009).In the last three years, several international collaborations have been established to use and test the algorithm online either at global or at regional scale. At teleseismic distances, the W phase method allows CMT estimates within half an hour after the earthquake origin time. At the regional scale, the CMT solutions are available between 6 min and 12 min after the origin time. The results obtained in real-time at the United States Geological Survey (USGS), the Pacific Tsunami Warning Center (PTWC), the Institut de Physique du Globe de Strasbourg (IPGS) and in the region of Mexico clearly indicate the robustness and the accuracy of CMT solutions computed using the W phase.Using a bayesian formulation, we propose a formal error analysis for seismic source inversion based on long period seismological data. Taking into account more realistic data uncertainties allows us to improve error estimates on the source model parameters but also to improve the solution itself. In this study, we also develop a method to estimate the uncertainty on the centroid depth. This parameter has a strong influence on the scalar seismic moment and on the fault dip estimated during CMT inversions.; Depuis une vingtaine d'années, un effort considérable a été effectué dans le développement d'outils dédiés à la caractérisation rapide de la source sismique. Jusqu'à récemment, plusieurs heures étaient encore nécessaires pour obtenir une information fiable sur la source des grands tremblements de terre. C'est dans ce cadre que s'inscrit ce travail, dont l'objectif principal est le développement et la mise en oeuvre d'une méthode de détermination rapide du tenseur moment sismique centroid (CMT) basée sur l'utilisation de la phase W pour les grands séismes. La phase W est une phase sismique correspondant essentiellement à la superposition des modes normaux supérieurs à très longue période (entre 100 sec et 1000 sec) arrivant entre l'onde P et les ondes de surface. Nous montrons que la phase W a l'avantage de ne pas être affectée par les hétérogénéités latérales superficielles puisque la majeure partie de son énergie se propage en profondeur. De par son caractère longue période, la phase W est particulièrement bien adapté à la caractérisation de la source pour les évènements de Mw>7.5. Pour les grands séismes en effet, la complexité de la source apparaît de plus en plus évidente aboutissant inévitablement à une variabilité des scénarios de sources estimés à courte période dans des bandes fréquentielles étroites. Les "séismes tsunami" et les séismes "outer-rise" sont des exemples extrêmes illustrant bien ce problème. Nous démontrons la robustesse de l'algorithme phase W développé pendant ce travail pour effectuer la caractérisation rapide et systématique de la source sismique. Il est aujourd'hui implémenté en temps réel dans plusieurs centres d'alerte à l'échelle globale et à l'échelle régionale. A des distances télésismiques, la méthode permet l'obtention du CMT dans la première demi-heure après le temps origine. A des échelles régionales, les résultats sont disponibles beaucoup plus rapidement, entre 6 et 12 minutes après le temps origine. Les résultats obtenus en temps réel à l'United States Geological Survey (USGS), au Pacific Tsunami Warning Center (PTWC), à l'Institut de Physique du Globe de Strasbourg (IPGS) et à l'échelle régionale au Mexique montrent clairement la robustesse et la fiabilité des solutions CMT obtenues en utilisant la phase W. En suivant une formulation bayesienne, nous proposons une analyse d'erreur formelle lors de l'inversion de la source sismique à longue période. La prise en compte plus réaliste de l'erreur sur les données permet non-seulement une estimation fiable de l'incertitude sur le modèle de source mais aussi l'amélioration de la solution elle même. Nous développons également une approche pour estimer l'erreursur la profondeur du centroid. Ce paramètre est important étant donné son influence sur le moment sismique scalaire et le pendage du plan de faille déterminés lors des inversions CMT.

Considerable effort has been made in the last two decades regarding the design and implementation of tools aimed at fast characterization of earthquake sources. As a fruitful harvest, a preliminary magnitude estimation is often available within fifteen minutes for a Mw=6.0 earthquake occurring anywhere in the world. But until recently, it still took several hours to determine the first order attributes of a great earthquake (e.g. Mw≥8.0) even in a well instrumented region. Our goal here is to develop a method providing fast centroid moment tensor (CMT) estimates for large earthquakes based on the W phase.The W phase is a very long period phase (100-1000 s) starting at the same time as the P-wave. It is conspicuous on broadband displacement records before the surface waves arrivals and can then be used for fast characterization of the source. We show that the W phase is not affected by shallow heterogeneities since a significant part of its energy is propagating deep into the mantle. Because of its long period nature, the W phase algorithm is particularly well suited for robust determination of CMT parameters for Mw≥7.5 earthquakes. For large earthquakes indeed, the source complexity is more evident which translates into significant variations of narrow-band magnitude estimations performed at short period. Tsunami earthquakes and outer-rise events provide good examples of this problem (Kanamori, 1972 ; Lay et al., 2009).In the last three years, several international collaborations have been established to use and test the algorithm online either at global or at regional scale. At teleseismic distances, the W phase method allows CMT estimates within half an hour after the earthquake origin time. At the regional scale, the CMT solutions are available between 6 min and 12 min after the origin time. The results obtained in real-time at the United States Geological Survey (USGS), the Pacific Tsunami Warning Center (PTWC), the Institut de Physique du Globe de Strasbourg (IPGS) and in the region of Mexico clearly indicate the robustness and the accuracy of CMT solutions computed using the W phase.Using a bayesian formulation, we propose a formal error analysis for seismic source inversion based on long period seismological data. Taking into account more realistic data uncertainties allows us to improve error estimates on the source model parameters but also to improve the solution itself. In this study, we also develop a method to estimate the uncertainty on the centroid depth. This parameter has a strong influence on the scalar seismic moment and on the fault dip estimated during CMT inversions., Depuis une vingtaine d'années, un effort considérable a été effectué dans le développement d'outils dédiés à la caractérisation rapide de la source sismique. Jusqu'à récemment, plusieurs heures étaient encore nécessaires pour obtenir une information fiable sur la source des grands tremblements de terre. C'est dans ce cadre que s'inscrit ce travail, dont l'objectif principal est le développement et la mise en oeuvre d'une méthode de détermination rapide du tenseur moment sismique centroid (CMT) basée sur l'utilisation de la phase W pour les grands séismes. La phase W est une phase sismique correspondant essentiellement à la superposition des modes normaux supérieurs à très longue période (entre 100 sec et 1000 sec) arrivant entre l'onde P et les ondes de surface. Nous montrons que la phase W a l'avantage de ne pas être affectée par les hétérogénéités latérales superficielles puisque la majeure partie de son énergie se propage en profondeur. De par son caractère longue période, la phase W est particulièrement bien adapté à la caractérisation de la source pour les évènements de Mw>7.5. Pour les grands séismes en effet, la complexité de la source apparaît de plus en plus évidente aboutissant inévitablement à une variabilité des scénarios de sources estimés à courte période dans des bandes fréquentielles étroites. Les "séismes tsunami" et les séismes "outer-rise" sont des exemples extrêmes illustrant bien ce problème. Nous démontrons la robustesse de l'algorithme phase W développé pendant ce travail pour effectuer la caractérisation rapide et systématique de la source sismique. Il est aujourd'hui implémenté en temps réel dans plusieurs centres d'alerte à l'échelle globale et à l'échelle régionale. A des distances télésismiques, la méthode permet l'obtention du CMT dans la première demi-heure après le temps origine. A des échelles régionales, les résultats sont disponibles beaucoup plus rapidement, entre 6 et 12 minutes après le temps origine. Les résultats obtenus en temps réel à l'United States Geological Survey (USGS), au Pacific Tsunami Warning Center (PTWC), à l'Institut de Physique du Globe de Strasbourg (IPGS) et à l'échelle régionale au Mexique montrent clairement la robustesse et la fiabilité des solutions CMT obtenues en utilisant la phase W. En suivant une formulation bayesienne, nous proposons une analyse d'erreur formelle lors de l'inversion de la source sismique à longue période. La prise en compte plus réaliste de l'erreur sur les données permet non-seulement une estimation fiable de l'incertitude sur le modèle de source mais aussi l'amélioration de la solution elle même. Nous développons également une approche pour estimer l'erreursur la profondeur du centroid. Ce paramètre est important étant donné son influence sur le moment sismique scalaire et le pendage du plan de faille déterminés lors des inversions CMT.

Slope instabilities take part in weathering and transport processes at the surface. The runout distance is extensively used in analysis of landslide dynamics and in the calibration of the rheological parameters involved in numerical modelling. However, the unknown impact of the uncertainty in the shape of the initial released mass on the runout distance and on the overall shape of the deposit questions the relevance of these approaches. Indeed, the shape of the initial scar is generally unknown in real cases. Our study is based on numerical simulations coupled with remote sensing data analysis. The model used in this study has been intensively compared with laboratory experiments and well constrained natural cases in order to establish its range of use. We have also developed a pre-event topographic reconstruction method using remote sensing data allowing the study of the topographic and initial failure plane geometry effects. We show that the runout distance is a robust parameter that is only poorly affected by the initial scar geometry. On the contrary, the extent of the deposits perpendicular to the main mass displacement direction is shown to be controlled by the scar geometry, providing a unique tool to retrieve information of the initial failure geometry, as well as on the released volume. A feedback analysis of Martian landslides shows excellent agreement between numerical results and geomorphological evidence, providing insight into the initial landsliding conditions. In addition, we introduce a new empirical dissipation parameter allowing a good prediction of the runout in the simulation without any calibration for a given geological context.; Les instabilités gravitaires contribuent aux processus de transport de matière à la surface des planètes et constituent des risques importants pour les populations sur Terre. Afin d'apporter des contraintes sur la dynamique des instabilités gravitaires, notre étude se base sur le couplage entre simulation numérique et analyse des données par télédétection. Le modèle utilisé s'est montré capable de reproduire des ob- servations en laboratoire ainsi que la dynamique enregistrée par des signaux sismiques pour un exemple terrestre. En outre, nous avons développé une méthode de reconstruction topographique pré-glissement en utilisant les données satellitaires afin d'étudier les effets des conditions initiales. En effet, la distance d'arrêt (runout) des dépôts est très largement utilisée dans l'analyse de la dynamique des glissements de terrain ainsi que dans la calibration des paramètres rhéologiques impliqués dans la modélisation, et ce malgré les incertitudes sur la géométrie initiale. Nous montrons par des tests théoriques que le runout est faiblement affecté par la géométrie du plan de rupture. Au contraire, l'extension latérale se révèle être contrôlée par la géométrie initiale ce qui fourni un cadre unique pour remonter également au volume initial. L'application aux exemples martiens apporte des contraintes sur les conditions de mise en place et ouvre de nouvelles perspectives pour la compréhension de la dynamique de ces processus. En outre, nous introduisons un nouveau paramètre de dissipation empirique indépendant de la topographie permettant, sans calibration, de retrouver par la simulation la bonne distance de runout dans un contexte géologique donné.

Slope instabilities take part in weathering and transport processes at the surface. The runout distance is extensively used in analysis of landslide dynamics and in the calibration of the rheological parameters involved in numerical modelling. However, the unknown impact of the uncertainty in the shape of the initial released mass on the runout distance and on the overall shape of the deposit questions the relevance of these approaches. Indeed, the shape of the initial scar is generally unknown in real cases. Our study is based on numerical simulations coupled with remote sensing data analysis. The model used in this study has been intensively compared with laboratory experiments and well constrained natural cases in order to establish its range of use. We have also developed a pre-event topographic reconstruction method using remote sensing data allowing the study of the topographic and initial failure plane geometry effects. We show that the runout distance is a robust parameter that is only poorly affected by the initial scar geometry. On the contrary, the extent of the deposits perpendicular to the main mass displacement direction is shown to be controlled by the scar geometry, providing a unique tool to retrieve information of the initial failure geometry, as well as on the released volume. A feedback analysis of Martian landslides shows excellent agreement between numerical results and geomorphological evidence, providing insight into the initial landsliding conditions. In addition, we introduce a new empirical dissipation parameter allowing a good prediction of the runout in the simulation without any calibration for a given geological context.; Les instabilités gravitaires contribuent aux processus de transport de matière à la surface des planètes et constituent des risques importants pour les populations sur Terre. Afin d'apporter des contraintes sur la dynamique des instabilités gravitaires, notre étude se base sur le couplage entre simulation numérique et analyse des données par télédétection. Le modèle utilisé s'est montré capable de reproduire des ob- servations en laboratoire ainsi que la dynamique enregistrée par des signaux sismiques pour un exemple terrestre. En outre, nous avons développé une méthode de reconstruction topographique pré-glissement en utilisant les données satellitaires afin d'étudier les effets des conditions initiales. En effet, la distance d'arrêt (runout) des dépôts est très largement utilisée dans l'analyse de la dynamique des glissements de terrain ainsi que dans la calibration des paramètres rhéologiques impliqués dans la modélisation, et ce malgré les incertitudes sur la géométrie initiale. Nous montrons par des tests théoriques que le runout est faiblement affecté par la géométrie du plan de rupture. Au contraire, l'extension latérale se révèle être contrôlée par la géométrie initiale ce qui fourni un cadre unique pour remonter également au volume initial. L'application aux exemples martiens apporte des contraintes sur les conditions de mise en place et ouvre de nouvelles perspectives pour la compréhension de la dynamique de ces processus. En outre, nous introduisons un nouveau paramètre de dissipation empirique indépendant de la topographie permettant, sans calibration, de retrouver par la simulation la bonne distance de runout dans un contexte géologique donné.