Tomographic resolution of ray and finite-frequency methods: A membrane-wave investigation

SUMMARY The purpose of this study is to evaluate the resolution potential of current finite-frequency approaches to tomography, and to do that in a framework similar to that of global scale seismology. According to our current knowledge and understanding, the only way to do this is by constructing a large set of ‘ground-truth’ synthetic data computed numerically (spectral elements, finite differences, etc.), and then to invert them using the various available linearized techniques. Specifically, we address the problem of using surface wave data to map phasevelocity distributions. Our investigation is strictly valid for the propagation of elastic waves on a spherical, heterogeneous membrane, and a good analogue for the propagation of surface waves within the outermost layers of the Earth. This amounts to drastically reducing the computational expense, with a certain loss of accuracy if very short-wavelength features of a strongly heterogeneous Earth are to be modelled. Our analysis suggests that a single-scattering finite-frequency approach to tomography, with sensitivity kernels computed via the adjoint method, is significantly more powerful than ray-theoretical methods, as a tool to image the fine structure of the Earth.

[1]  L. Boschi.  Global multiresolution models of surface wave propagation: comparing equivalently regularized Born and ray theoretical solutions , 2006 .

[2]  Anatoli L. Levshin,et al.  Global surface wave diffraction tomography , 2002 .

[4]  L. Boschi.  Measures of resolution in global body wave tomography , 2003 .

[5]  Walter H. F. Smith,et al.  Free software helps map and display data , 1991 .

[6]  Guust Nolet,et al.  Fréchet kernels for finite‐frequency traveltimes—II. Examples , 2000 .

[7]  F. Dahlen,et al.  Traveltimes of waves in three-dimensional random media , 2003 .

[8]  D. Giardini,et al.  Petascale computing and resolution in global seismic tomography , 2007 .

[9]  R. Snieder,et al.  The effect of small-scale heterogeneity on the arrival time of waves , 2001 .

[10]  Y. Capdeville.  An efficient Born normal mode method to compute sensitivity kernels and synthetic seismograms in the Earth , 2005 .

[11]  Lapo Boschi,et al.  On the relevance of Born theory in global seismic tomography , 2006 .

[12]  J. Woodhouse,et al.  Source, path and receiver effects on seismic surface waves , 2007 .

[13]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[14]  G. Nolet,et al.  Diffraction tomography using multimode surface waves , 1997 .

[15]  Guust Nolet,et al.  Three‐dimensional sensitivity kernels for surface wave observables , 2004 .

[16]  S. Hung,et al.  Validation of ray and wave theoretical travel times in heterogeneous random media , 2005 .

[17]  T. Tanimoto.  Modelling curved surface wave paths: membrane surface wave synthetics , 1990 .

[18]  G. Masters,et al.  Constraints on global phase velocity maps from long-period polarization data , 1996 .

[19]  Jeannot Trampert,et al.  Global phase velocity maps of Love and Rayleigh waves between 40 and 150 seconds , 1995 .

[20]  Luis Rivera,et al.  On the use of the checker-board test to assess the resolution of tomographic inversions , 1993 .

[21]  Guust Nolet,et al.  Three-dimensional waveform sensitivity kernels , 1998 .

[22]  G. Laske.  Global observation of off-great-circle propagation of long-period surface waves , 1995 .

[23]  G. Nolet,et al.  Finite-frequency effects in global surface-wave tomography , 2005 .

[24]  L. Boschi,et al.  New images of the Earth's upper mantle from measurements of surface wave phase velocity anomalies , 2002 .

[25]  Guust Nolet,et al.  Fréchet kernels for finite-frequency traveltimes—I. Theory , 2000 .

[26]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[27]  F. Dahlen,et al.  Statistics of traveltimes and amplitudes in random media , 2004 .

[28]  L. Boschi,et al.  Surface wave ray tracing and azimuthal anisotropy: a generalized spherical harmonic approach , 2006 .

[29]  Jean-Paul Montagner,et al.  SPICE benchmark for global tomographic methods , 2008 .

[30]  Qinya Liu,et al.  Tomography, Adjoint Methods, Time-Reversal, and Banana-Doughnut Kernels , 2004 .

[31]  Jeannot Trampert,et al.  Assessment of global phase velocity models , 2001 .

[32]  J. Woodhouse,et al.  GLOBAL HIGH-RESOLUTION PHASE VELOCITY DISTRIBUTIONS OF OVERTONE AND FUNDAMENTAL-MODE SURFACE WAVES DETERMINED BY MODE BRANCH STRIPPING , 1999 .

[33]  Jeroen Tromp,et al.  Measurements and global models of surface wave propagation , 1997 .

[34]  J. Woodhouse,et al.  Amplitude, phase and path anomalies of mantle waves , 1986 .

[35]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[36]  F. Dahlen,et al.  JWKB surface‐wave seismograms on a laterally heterogeneous earth , 1994 .

[37]  R. Snieder,et al.  The effect of scattering in surface wave tomography , 2002 .

[38]  J. Lévêque,et al.  Can finite‐frequency effects be accounted for in ray theory surface wave tomography? , 2004 .

[39]  John A. Nelder,et al.  A Simplex Method for Function Minimization , 1965, Comput. J..

[40]  Jeannot Trampert,et al.  FAST TRACK PAPER: Surface wave tomography: finite-frequency effects lost in the null space , 2006 .

[41]  A. Dziewoński,et al.  Seismic Surface Waves and Free Oscillations in a Regionalized Earth Model , 1982 .

[42]  B. Kennett,et al.  Sensitivity kernels for finite-frequency surface waves , 2002 .

[43]  C. Tape,et al.  Surface wave tomography: global membrane waves and adjoint methods , 2007 .

[44]  F. Gilbert.  The Representation of Seismic Displacements in Terms of Travelling Waves , 1976 .

[45]  F. A. Dahlen,et al.  Resolution limit of traveltime tomography , 2004 .

[46]  Jeannot Trampert,et al.  High resolution global phase velocity distributions , 1996 .