Comparison of global waveform inversions with and without considering cross-branch modal coupling

SUMMARY In this study we present a new approach to inverting long-period seismic waveforms. The effect of lateral heterogeneity is partitioned into two. The first part represents the effect of horizontally averaged structure along the great circle between the source and receiver, and is allowed to remain in non-linear form in the formulation. The second part incorporates any further correction due to cross-branch modal coupling, which has been neglected in the more conventional path average approximation (PAVA). This term is linearized and then treated asymptotically so that the seismogram depends only upon the structure within the great-circle section determined by the source and receiver (Li & Tanimoto 1993a). We refer to this new method as the non-linear asymptotic coupling theory (NACT). The sensitivity kernels predicted by the PAVA and NACT are compared. While the sensitivity kernels are similar for surface waves and shallow-turning body waves, they are very different for body waves that sample the deep mantle. By examining the inversion algorithms for the PAVA and NACT, we demonstrate that the computation time required by the NACT tends to be of the same order of magnitude as that required by the PAVA, as the number of model parameters increases. Based upon a realistically large data set (5041 body-wave seismograms and 1531 mantle-wave seismograms), formal resolution analyses are performed using both PAVA and NACT. We find that the NACT is significantly more powerful in resolving 3-D structure in the deep mantle. We compare the models obtained for the same observed data set using the two approaches. As expected, they differ more in the lower mantle than in the upper mantle. The difference in their amplitude spectra increases with spherical harmonic degree. The model developed using the NACT predicts the observed surface geoid better than that developed using the PAVA, based upon geodynamic flow modelling.

[1]  Paul G. Richards,et al.  Quantitative Seismology: Theory and Methods , 1980 .

[2]  Toshiro Tanimoto,et al.  Waveforms of long-period body waves in a slightly aspherical earth model , 1993 .

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

[4]  R. Snieder,et al.  A new formalism for the effect of lateral heterogeneity on normal modes and surface waves-II. General anisotropic perturbation , 1988 .

[5]  J. Woodhouse,et al.  Constraining Upper Mantle Anelasticity Using Surface Wave Amplitude Anomalies , 1993 .

[6]  A. Dziewoński,et al.  Acoustic imaging at the planetary scale , 1992 .

[7]  B. Romanowicz,et al.  Modelling of coupled normal modes of the Earth: the spectral method , 1990 .

[8]  T. Tanimoto Free oscillations of a slightly anisotropic earth , 1986 .

[9]  Adam M. Dziewonski,et al.  Mapping the lower mantle: Determination of lateral heterogeneity in P velocity up to degree and order 6 , 1984 .

[10]  Jeffrey Park,et al.  Normal mode multiplet coupling along a dispersion branch , 1991 .

[11]  Robert L. Woodward,et al.  Global upper mantle structure from long-period differential travel times , 1991 .

[12]  F. Dahlen Multiplet coupling and the calculation of synthetic long-period seismograms , 1987 .

[13]  Jeffrey Park,et al.  The subspace projection method for constructing coupled-mode synthetic seismograms , 1990 .

[14]  B. Romanowicz The upper mantle degree 2: Constraints and inferences from global mantle wave attenuation measurements , 1990 .

[15]  G. Masters,et al.  Lower-mantle structure from ScS–S differential travel times , 1991, Nature.

[16]  Wei-jia Su,et al.  Degree 12 model of shear velocity heterogeneity in the mantle , 1994 .

[17]  Wei-jia Su,et al.  Predominance of long-wavelength heterogeneity in the mantle , 1991, Nature.

[18]  R. Geller,et al.  Inversion for laterally heterogeneous upper mantle S-wave velocity structure using iterative waveform inversion , 1993 .

[19]  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 .

[20]  J. Woodhouse The coupling and attenuation of nearly resonant multiplets in the Earth's free oscillation spectrum , 1980 .

[21]  Tatsuhiko Hara,et al.  Two efficient algorithms for iterative linearized inversion of seismic waveform data , 1993 .

[22]  B. Romanowicz Multiplet-multiplet coupling due to lateral heterogeneity: asymptotic effects on the amplitude and frequency of the Earth's normal modes , 1987 .

[23]  E. Mochizuki The Free Oscillations of an Anisotropic and Heterogeneous Earth , 1986 .

[24]  F. Gilbert,et al.  An application of normal mode theory to the retrieval of structural parameters and source mechanisms from seismic spectra , 1975, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[25]  T. Tanimoto The Backus—Gilbert approach to the three-dimensional structure in the upper mantle — I. Lateral variation of surface wave phase velocity with its error and resolution , 1985 .

[26]  Yosihiko Ogata,et al.  Whole mantle P-wave travel time tomography , 1990 .

[27]  T. Tanimoto Long-wavelength S-wave velocity structure throughout the mantle , 1990 .

[28]  C. Froidevaux,et al.  Converting mantle tomography into mass anomalies to predict the Earth's radial viscosity , 1994 .

[29]  Jeffrey Park Asymptotic coupled-mode expressions for multiplet amplitude anomalies and frequency shifts on an aspherical earth , 1987 .

[30]  Toshiro Tanimoto,et al.  High-resolution global upper mantle structure and plate tectonics , 1993 .

[31]  B. Romanowicz Anelastic tomography: a new perspective on upper mantle thermal structure , 1994 .

[32]  P. Richards,et al.  Borovoye Geophysical Observatory, Kazakhstan , 1992 .

[33]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[34]  B. Hager,et al.  Geoid Anomalies in a Dynamic Earth , 1984 .

[35]  Philip B. Stark,et al.  Toward tubular tomography , 1993 .

[36]  T. Jordan A procedure for estimating lateral variations from low-frequency eigenspectra data , 1978 .

[37]  T. Tanimoto The Backus–Gilbert approach to the 3-D structure in the upper mantle – II. SH and SV velocity , 1986 .

[38]  John H. Woodhouse,et al.  Mapping the upper mantle: Three‐dimensional modeling of earth structure by inversion of seismic waveforms , 1984 .

[39]  E. Mochizuki Free oscillations and surface waves of an aspherical Earth , 1986 .