Solving large linear inverse problems in geophysics by means of Eigenvalue calculations

This paper discusses the computation and utilization of singular values and singular vectors in the solution of very large inverse problems that arise in the study of physical models for the internal structure of the Earth. In this study, the Earth is discretized into layers and the layers into cells, and travel times of sound waves generated by earthquakes are used to construct the models. The underlying numerical models correspond to sparse matrices with dimensions up to 1.3x10-by-3x10. Singular values and singular vectors of these matrices are then computed and used in the solution of the inverse problems and to estimate uncertainties. The paper outlines the formulation adopted to model the Earth and the strategy employed to compute singular values and singular vectors, shows results for two models that have been studied, and comments on the main computation issues related to the solution of these problems on high performance parallel computers.

[1]  E. R. Engdahl,et al.  Constraints on seismic velocities in the Earth from traveltimes , 1995 .

[2]  G. McMechan,et al.  Estimation of resolution and covariance for large matrix inversions , 1995 .

[3]  S. Treitel,et al.  Fast l p solution of large, sparse, linear systems: application to seismic travel time tomography , 1988 .

[4]  D. Vasco,et al.  Lateral variations in mantle velocity structure and discontinuities determined from P, PP, S, SS, and SS — SdS travel time residuals , 1995 .

[5]  C. D. Meyer,et al.  Generalized inverses of linear transformations , 1979 .

[6]  D. Vasco,et al.  Global Earth structure: inference and assessment , 2002 .

[7]  W. Menke.  Geophysical data analysis : discrete inverse theory , 1984 .

[8]  J. Virieux,et al.  Explicit, approximate expressions for the resolution and a posteriori covariance of massive tomographic systems , 1999 .

[9]  B. Kennett.  Seismic wave propagation and seismic tomography , 1998 .

[10]  D. Vasco,et al.  Whole Earth structure estimated from seismic arrival times , 1998 .

[11]  Jack Dongarra,et al.  ScaLAPACK user's guide , 1997 .

[12]  Don W. Vasco,et al.  Formal inversion of ISC arrival times for mantle P-velocity structure , 1993 .

[13]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[14]  B. Parlett.  The Symmetric Eigenvalue Problem , 1981 .

[15]  D. Sorensen.  Numerical methods for large eigenvalue problems , 2002, Acta Numerica.

[16]  R. Parker.  Geophysical Inverse Theory , 1994 .

[17]  Don W. Vasco,et al.  Tomographic inversions for mantle P wave velocity structure based on the minimization of l 2 and l 1 norms of International Seismological Centre Travel Time Residuals , 1993 .

[18]  R. Roberts,et al.  Calculating resolution and covariance matrices for seismic tomography with the LSQR method , 1999 .