Premium
Measures of resolution in global body wave tomography
Author(s) -
Boschi L.
Publication year - 2003
Publication title -
geophysical research letters
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.007
H-Index - 273
eISSN - 1944-8007
pISSN - 0094-8276
DOI - 10.1029/2003gl018222
Subject(s) - algorithm , cholesky decomposition , computer science , inverse problem , similarity measure , synthetic data , measure (data warehouse) , inversion (geology) , tomography , image resolution , mathematics , data mining , artificial intelligence , geology , eigenvalues and eigenvectors , mathematical analysis , optics , physics , paleontology , quantum mechanics , structural basin
The resolution of tomographic images is most often evaluated through synthetic tests: the inversion algorithm used to derive the image itself is applied to a synthetic data set having the same source‐station geometry of the real one, but theoretical travel times computed from a chosen “input model” (e.g., a checkerboard). The similarity between input model and solution of the synthetic test, used as a measure of resolution, has the major shortcoming of depending on the choice of the input model. Conversely, the similarity of the “model resolution matrix” to the identity matrix is a rigorous measure of resolution that does not depend on any input model, but has the drawback of being computationally heavy. In the past decade, several authors have devised complicated algorithms for the approximate or iterative derivation of . I show here that parallel Cholesky factorization of ( being the matrix that identifies the linear inverse problem), feasible on shared‐memory multiprocessor servers, provides an efficient way of determining both least squares solutions and resolution matrices in global tomography. I apply this procedure in an evaluation of the resolution of mantle structure from a global P ‐wave travel time data set.