Skip to main content
Article
Solving Large Tomographic Linear Systems: Size Reduction and Error Estimation
Geophysical Journal International (2014)
  • Sergey Voronin, Université de Nice
  • Dylan Mikesell, Université de Nice
  • Inna Slezak, Université de Nice
  • Guust Nolet, Université de Nice
Abstract
We present a new approach to reduce a sparse, linear system of equations associated with tomographic inverse problems. We begin by making a modification to the commonly used compressed sparse-row format, whereby our format is tailored to the sparse structure of finite-frequency (volume) sensitivity kernels in seismic tomography. Next, we cluster the sparse matrix rows to divide a large matrix into smaller subsets representing ray paths that are geographically close. Singular value decomposition of each subset allows us to project the data onto a subspace associated with the largest eigenvalues of the subset. After projection we reject those data that have a signal-to-noise ratio (SNR) below a chosen threshold. Clustering in this way assures that the sparse nature of the system is minimally affected by the projection. Moreover, our approach allows for a precise estimation of the noise affecting the data while also giving us the ability to identify outliers. We illustrate the method by reducing large matrices computed for global tomographic systems with cross-correlation body wave delays, as well as with surface wave phase velocity anomalies. For a massive matrix computed for 3.7 million Rayleigh wave phase velocity measurements, imposing a threshold of 1 for the SNR, we condensed the matrix size from 1103 to 63 Gbyte. For a global data set of multiple-frequency P wave delays from 60 well-distributed deep earthquakes we obtain a reduction to 5.9 per cent. This type of reduction allows one to avoid loss of information due to underparametrizing models. Alternatively, if data have to be rejected to fit the system into computer memory, it assures that the most important data are preserved.
Keywords
  • inverse theory,
  • body waves,
  • surface waves and free oscillations,
  • computational seismology
Publication Date
October, 2014
DOI
10.1093/gji/ggu242
Citation Information
Sergey Voronin, Dylan Mikesell, Inna Slezak and Guust Nolet. "Solving Large Tomographic Linear Systems: Size Reduction and Error Estimation" Geophysical Journal International Vol. 199 Iss. 1 (2014) p. 276 - 285
Available at: http://works.bepress.com/dylan_mikesell/25/