Local Earthquake Tomography by trans-dimensional Monte Carlo Sampling

Giacomuzzi, Genny (2013) Local Earthquake Tomography by trans-dimensional Monte Carlo Sampling, [Dissertation thesis], Alma Mater Studiorum Università di Bologna. Dottorato di ricerca in Geofisica, 25 Ciclo. DOI 10.6092/unibo/amsdottorato/5160.
Documenti full-text disponibili:
Documento PDF (English) - Richiede un lettore di PDF come Xpdf o Adobe Acrobat Reader
Download (50MB) | Anteprima


In this study a new, fully non-linear, approach to Local Earthquake Tomography is presented. Local Earthquakes Tomography (LET) is a non-linear inversion problem that allows the joint determination of earthquakes parameters and velocity structure from arrival times of waves generated by local sources. Since the early developments of seismic tomography several inversion methods have been developed to solve this problem in a linearized way. In the framework of Monte Carlo sampling, we developed a new code based on the Reversible Jump Markov Chain Monte Carlo sampling method (Rj-McMc). It is a trans-dimensional approach in which the number of unknowns, and thus the model parameterization, is treated as one of the unknowns. I show that our new code allows overcoming major limitations of linearized tomography, opening a new perspective in seismic imaging. Synthetic tests demonstrate that our algorithm is able to produce a robust and reliable tomography without the need to make subjective a-priori assumptions about starting models and parameterization. Moreover it provides a more accurate estimate of uncertainties about the model parameters. Therefore, it is very suitable for investigating the velocity structure in regions that lack of accurate a-priori information. Synthetic tests also reveal that the lack of any regularization constraints allows extracting more information from the observed data and that the velocity structure can be detected also in regions where the density of rays is low and standard linearized codes fails. I also present high-resolution Vp and Vp/Vs models in two widespread investigated regions: the Parkfield segment of the San Andreas Fault (California, USA) and the area around the Alto Tiberina fault (Umbria-Marche, Italy). In both the cases, the models obtained with our code show a substantial improvement in the data fit, if compared with the models obtained from the same data set with the linearized inversion codes.

Tipologia del documento
Tesi di dottorato
Giacomuzzi, Genny
Dottorato di ricerca
Scuola di dottorato
Scienze matematiche, fisiche ed astronomiche
Settore disciplinare
Settore concorsuale
Parole chiave
seismic tomography san andreas fault monte carlo Reversible jump Markov Chain
Data di discussione
22 Febbraio 2013

Altri metadati

Statistica sui download

Gestione del documento: Visualizza la tesi