Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published August 2012 | Published
Journal Article Open

Uncertainty estimations for seismic source inversions

Abstract

Source inversion is a widely used practice in seismology. Magnitudes, moment tensors, slip distributions are now routinely calculated and disseminated whenever an earthquake occurs. The accuracy of such models depends on many aspects like the event magnitude, the data coverage and the data quality (instrument response, isolation, timing, etc.). Here, like in any observational problem, the error estimation should be part of the solution. It is however very rare to find a source inversion algorithm which includes realistic error analyses, and the solutions are often given without any estimates of uncertainties. Our goal here is to stress the importance of such estimation and to explore different techniques aimed at achieving such analyses. In this perspective, we use the W phase source inversion algorithm recently developed to provide fast CMT estimations for large earthquakes. We focus in particular on the linear-inverse problem of estimating the moment tensor components at a given source location. We assume that the initial probability densities can be modelled by Gaussian distributions. Formally, we can separate two sources of error which generally contribute to the model parameter uncertainties. The first source of uncertainty is the error introduced by the more or less imperfect data. This is carried by the covariance matrix for the data (C_d). The second source of uncertainty, often overlooked, is associated with modelling error or mismodelling. This is represented by the covariance matrix on the theory, C_T. Among the different sources of mismodelling, we focus here on the modelling error associated with the mislocation of the centroid position. Both C_d and C_T describe probability densities in the data space and it is well known that it is in fact C_D = C_d + C_T that should be included into the error propagation process. In source inversion problems, like in many other fields of geophysics, the data covariance (C_D) is often considered as diagonal or even proportional to the identity matrix. In this work, we demonstrate the importance of using a more realistic form for C_D. If we incorporate accurate covariance components during the inversion process, it refines the posterior error estimates but also improves the solution itself. We discuss these issues using several synthetic tests and by applying the W phase source inversion algorithm to several large earthquakes such as the recent 2011 Tohoku-oki earthquake.

Additional Information

© 2012 The Authors. Geophysical Journal International © 2012 RAS. Accepted 2012 May 21. Received 2012 May 21; in original form 2011 December 1. This work uses Federation of Digital Seismic Networks (FDSN) seismic data and CMT solutions from the Global CMT catalog. The Incorporated Research Institutions for Seismology (IRIS) Data Management System (DMS) was used to access the data. This work made use of the Matplotlib python library. We thank J.-J. Lévêque, M. Sambridge and an anonymous reviewer for their helpful comments on the manuscript.

Attached Files

Published - Duputel2012p19103Geophys_J_Int.pdf

Files

Duputel2012p19103Geophys_J_Int.pdf
Files (7.7 MB)
Name Size Download all
md5:73f824828837cdc15e5c48e763537fd9
7.7 MB Preview Download

Additional details

Created:
August 22, 2023
Modified:
October 18, 2023