The three‐dimensional shear velocity structure of the mantle from the inversion of body, surface and higher‐mode waveforms

SUMMARY We present a 3-D model of shear heterogeneity (Vsh) in the whole mantle, derived from the inversion of hand-picked body, surface and higher-mode waveforms. The forward and inverse problems are formulated using the non-linear asymptotic coupling theory, the zeroth-order asymptotic expansion of a Born seismogram computed by normal mode summation and including coupling across mode branches. We invert iteratively for mantle heterogeneity and for centroid moment tensors until convergence. Our model, SAW24B16, is parametrized laterally in spherical harmonics up to degree 24 and radially in 16 cubic b-splines with knots spaced to reflect the data sampling with depth. The power spectrum of the model is dominated by heterogeneity in the boundary layers at low degrees: degrees 5 and 6 near the surface and degree 2 near the core–mantle boundary. The rest of the lower mantle is dominated by degree 3. We find significant heterogeneity up to spherical harmonic degree 24 in the transition zone and the uppermost part of the lower mantle. Our model displays three slow domains extending continuously between uppermost and lowermost mantle, one originating under Africa and two in the Pacific. The results of our modelling also suggest an interaction near the surface between the anomalies in the Pacific and the network of mid-ocean ridges in the East Pacific, and between the African anomaly and the low-velocity zone in the Red Sea, the East African Rift, the Mid-Indian Ridge and the Mid-Atlantic Ridge. These anomalies appear to ascend into the upper mantle without altering their shape near the 670 km discontinuity. Fast anomalies accumulate near the 670 km discontinuity under the Southern Kurile, the Japan trench and the Izu arc, with little or no penetration into the lower mantle, in agreement with results from a variety of regional studies. Cross-sections through subduction zones such as the Marianas, Tonga–Kermadec, Java, the palaeosubduction zone beneath North America and the Peru–Chile trench show continuous fast structures penetrating into the lower mantle, with the depth of penetration varying from region to region. The distribution of both fast and slow anomalies near 670 km thus suggests that the discontinuity does not act as a strong barrier to upgoing and downgoing flow between the upper and the lower mantle.

[1]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[2]  G. Schubert,et al.  Structure of axisymmetric mantle plumes , 1993 .

[3]  Romanowicz,et al.  Three-dimensional structure at the base of the mantle beneath the central pacific , 1998, Science.

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

[5]  J. Montagner,et al.  Tomography of the transition zone from the inversion of higher-mode surface waves , 1994 .

[6]  T. Tanimoto A simple derivation of the formula to calculate synthetic long‐period seismograms in a heterogeneous earth by normal mode summation , 1984 .

[7]  B. Steinberger,et al.  Advection of plumes in mantle flow: implications for hotspot motion, mantle viscosity and plume distribution , 1998 .

[8]  B. Romanowicz,et al.  Test of tomographic models of D″ using differential travel time data , 1998 .

[9]  S. Grand Mantle shear structure beneath the Americas and surrounding oceans , 1994 .

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

[11]  Guust Nolet,et al.  Tomographic imaging of subducted lithosphere below northwest Pacific island arcs , 1991, Nature.

[12]  Don W. Vasco,et al.  Robust inversion of IASP91 travel time residuals for mantle P and S velocity structure, earthquake mislocations, and station corrections , 1994 .

[13]  Zheng Wang,et al.  Spherical-spline parameterization of three-dimensional Earth models , 1995 .

[14]  R. Hilst,et al.  Compositional heterogeneity in the bottom 1000 kilometers of Earth's mantle: toward a hybrid convection model , 1999, Science.

[15]  R. V. D. Hilst Complex morphology of subducted lithosphere in the mantle beneath the Tonga trench , 1995, Nature.

[16]  H. Takeuchi,et al.  Seismic Surface Waves , 1972 .

[17]  Bradford H. Hager,et al.  Large‐scale heterogeneities in the lower mantle , 1977 .

[18]  R. Hilst,et al.  Effects of relative plate motion on the deep structure and penetration depth of slabs below the Izu-Bonin and Mariana island arcs , 1993 .

[19]  Sri Widiyantoro,et al.  Global seismic tomography: A snapshot of convection in the Earth: GSA Today , 1997 .

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

[21]  Norman H. Sleep,et al.  Hotspots and Mantle Plumes' Some Phenomenology , 1990 .

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

[23]  Barbara Romanowicz,et al.  Imaging 3‐D spherical convection models: What can seismic tomography tell us about mantle dynamics? , 1997 .

[24]  C. Beaumont,et al.  Nano variations in gravity due to seasonal groundwater movements: Implications for the gravitational detection of tectonic movements , 1977 .

[25]  Joseph S. Resovsky,et al.  A degree 8 mantle shear velocity model from normal mode observations below 3 mHz , 1999 .

[26]  G. Nolet,et al.  Upper mantle S velocity structure of North America , 1997 .

[27]  Richards,et al.  Time scales and heterogeneous structure in geodynamic earth models , 1998, Science.

[28]  Jeroen Tromp,et al.  Global and regional surface‐wave inversions: A spherical‐spline parameterization , 1998 .

[29]  Masayuki Obayashi,et al.  Subducting slabs stagnant in the mantle transition zone , 1992 .

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

[31]  A. Tarantola,et al.  Generalized Nonlinear Inverse Problems Solved Using the Least Squares Criterion (Paper 1R1855) , 1982 .

[32]  T. Jordan,et al.  Comparisons Between Seismic Earth Structures and Mantle Flow Models Based on Radial Correlation Functions , 1993, Science.

[33]  A. R. Edmonds Angular Momentum in Quantum Mechanics , 1957 .

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

[35]  T. Jordan,et al.  Mapping the Tonga Slab , 1991 .

[36]  R. Hilst,et al.  Compositional stratification in the deep mantle , 1999, Science.

[37]  Barbara Romanowicz,et al.  Comparison of global waveform inversions with and without considering cross-branch modal coupling , 1995 .

[38]  G. Laske,et al.  A shear - velocity model of the mantle , 1996, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[39]  Alberto Michelini,et al.  Seismological studies at Parkfield. I. Simultaneous inversion for velocity structure and hypocenters using cubic B-splines parameterization , 1991, Bulletin of the Seismological Society of America.

[40]  J. Montagner,et al.  An inverse technique for retrieving higher mode phase velocity and mantle structure , 1993 .

[41]  P. Olson,et al.  The proximity of hotspots to convergent and divergent plate boundaries , 1989 .

[42]  E. Engdahl,et al.  Global teleseismic earthquake relocation with improved travel times and procedures for depth determination , 1998, Bulletin of the Seismological Society of America.

[43]  Hendrik Jan van Heijst,et al.  Measuring surface-wave overtone phase velocities using a mode-branch stripping technique , 1997 .

[44]  T. Tanimoto The three-dimensional shear wave structure in the mantle by overtone waveform inversion - I. Radial seismogram inversion , 1987 .

[45]  Y. Fukao,et al.  Evaluation of slab images in the northwestern Pacific , 1998 .

[46]  Harmen Bijwaard,et al.  Closing the gap between regional and global travel time tomography , 1998 .

[47]  D. Jackson The use of a priori data to resolve non‐uniqueness in linear inversion , 1979 .

[48]  S. Grand,et al.  Seismic structure of the deep Kurile subduction zone , 1994 .

[49]  R. Snieder,et al.  Surface‐wave mode coupling for efficient forward modelling and inversion of body‐wave phases , 1995 .

[50]  A. Dziewoński,et al.  Models of the mantle shear velocity and discontinuities in the pattern of lateral heterogeneities , 2001 .

[51]  S. Grand,et al.  Variation of transition zone high‐velocity anomalies and depression of 660 km discontinuity associated with subduction zones from the southern Kuriles to Izu‐Bonin and Ryukyu , 1998 .

[52]  D. Helmberger,et al.  Evidence for strong shear velocity reductions and velocity gradients in the lower mantle beneath Africa , 1998 .

[53]  G. Ekström,et al.  A radial model of anelasticity consistent with long-period surface-wave attenuation , 1996 .

[54]  Barbara Romanowicz,et al.  Global mantle shear velocity model developed using nonlinear asymptotic coupling theory , 1996 .

[55]  B. Romanowicz,et al.  The effects of the theoretical formalism and data selection on mantle models derived from waveform tomography , 1999 .

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

[57]  William H. Press,et al.  Numerical recipes in C , 2002 .

[58]  E. R. Engdahl,et al.  Evidence for deep mantle circulation from global tomography , 1997, Nature.

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

[60]  D. L. Anderson,et al.  Measurements of mantle wave velocities and inversion for lateral heterogeneities and anisotropy: 3. Inversion , 1986 .