Global tomographic images of mantle plumes and subducting slabs: insight into deep Earth dynamics

Abstract A new model of whole mantle P-wave tomography is determined with a novel approach. A grid parameterization instead of blocks and spherical harmonics is adopted to express the Earth structure. Depth variations of the Moho, 410 and 660 km discontinuities are taken into account in the inversion. Ray paths and travel times are computed with an efficient 3-D ray tracing scheme. This new approach was applied to a large data set of ISC (International Seismological Center) travel times (P, PP, PcP, pP, Pdiff) to determine a whole mantle P-wave tomography. For the shallow mantle, the new model contains the general features observed in the previous models: a low-velocity ring around the Pacific Ocean basins and high-velocity anomalies under the old and stable continents in the depth range of 0–300 km. One significant difference from the previous models is that stronger and wider high-velocity anomalies are visible in the transition zone depths under the subduction regions, which suggests that most of the slab materials are stagnant in the transition zone before finally collapsing down to the lower mantle as a result of very large gravitational instability from phase transitions. Very slow anomalies exist in the upper mantle right beneath the Wudalianchi and Changbai active volcanoes in Eastern China, right above the stagnant Pacific slab in the transition zone, suggesting that the origin of the intraplate volcanism in East Asia is closely related to the Pacific plate subduction process, such as deep slab dehydration and convective circulation in the mantle wedge. Plume-like slow anomalies are clearly visible under the major hotspot regions in most parts of the mantle, in particular, under Hawaii, Iceland, South Pacific and Africa. The slow anomalies under South Pacific and Africa have lateral extensions of over one thousand kilometers and exist in the entire mantle, representing two superplumes. The Pacific superplume has a larger spatial extent and stronger slow anomalies than that of the Africa superplume. The Hawaiian plume is not part of the Pacific superplume. The slow anomalies under hotspots usually do not show a straight pillar shape, but exhibit winding images, suggesting that plumes are not fixed in the mantle but can be deflected by the mantle flow. As a consequence, hotspots are not really fixed but can wander on the Earth’s surface, as evidenced by the recent paleomagnetic and numeric modeling studies. Wider and more prominent slow anomalies are visible at the core–mantle boundary (CMB) than most of the lower mantle, and there is a good correlation between the distribution of slow anomalies at the CMB and that of hotspots on the surface, suggesting that most of the strong mantle plumes under the hotspots originate from the CMB. However, there are some small-scaled, weak plumes originating from the transition zone or mid mantle depths.

[1]  Michael E. Wysession,et al.  The structure of the core‐mantle boundary from diffracted waves , 1992 .

[2]  Hiroo Kanamori,et al.  Deep structure of Japan subduction zone as derived from local, regional, and teleseismic events , 1994 .

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

[4]  M. Richards,et al.  The dynamics of Cenozoic and Mesozoic plate motions , 1998 .

[5]  Dapeng Zhao,et al.  Seismic structure and origin of hotspots and mantle plumes , 2001 .

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

[7]  Clifford H. Thurber,et al.  A fast algorithm for two-point seismic ray tracing , 1987 .

[8]  D. Yuen,et al.  Geophysical inferences of thermal‐chemical structures in the lower mantle , 1993 .

[9]  Gillian R. Foulger,et al.  Upper-mantle origin of the Yellowstone hotspot , 2002 .

[10]  D. Helmberger,et al.  Seismic evidence that the source of the Iceland hotspot lies at the core–mantle boundary , 1998, Nature.

[11]  E. Garnero Heterogeneity of the Lowermost Mantle , 2000 .

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

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

[14]  P. Molnar,et al.  Relative motions of hotspots in the Pacific, Atlantic and Indian Oceans since late Cretaceous time , 1987, Nature.

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

[16]  W. J. Morgan,et al.  Convection Plumes in the Lower Mantle , 1971, Nature.

[17]  J. Woodhouse,et al.  Complex Shear Wave Velocity Structure Imaged Beneath Africa and Iceland. , 1999, Science.

[18]  Akira Hasegawa,et al.  Tomographic imaging of P and S wave velocity structure beneath northeastern Japan , 1992 .

[19]  Jean Besse,et al.  Three distinct types of hotspots in the Earth's mantle , 2002 .

[20]  Eiichi Takahahshi,et al.  Origin of the Columbia River basalts: melting model of a heterogeneous plume head , 1998 .

[21]  T. J. Owens,et al.  Seismic evidence for a deep upper mantle thermal anomaly beneath east Africa , 2000 .

[22]  D. Turcotte,et al.  How many plumes are there , 1999 .

[23]  D. L. Anderson Chemical plumes in the mantle , 1975 .

[24]  E. Garnero,et al.  Seismic Evidence for Partial Melt at the Base of Earth's Mantle , 1996, Science.

[25]  Thorne Lay,et al.  The core–mantle boundary layer and deep Earth dynamics , 1998, Nature.

[26]  R. Kind,et al.  Seismic study of upper mantle and transition zone beneath hotspots , 2003 .

[27]  F. Niu,et al.  Complex Structure of Mantle Discontinuities at the Tip of the Subducting Slab beneath Northeast China -A Preliminary Investigation of Broadband Receiver Functions- , 1996 .

[28]  J. Tarduno,et al.  Paleomagnetic evidence for motion of the Hawaiian hotspot during formation of the Emperor seamounts , 1997 .

[29]  F. Niu,et al.  Seismic evidence for a thinner mantle transition zone beneath the South Pacific Superswell , 2000 .

[30]  M. Stefanick,et al.  The distribution of hot spots , 1984 .

[31]  Dapeng Zhao,et al.  New advances of seismic tomography and its applications to subduction zones and earthquake fault zones: A review , 2001 .

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

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

[34]  D. Yuen,et al.  Mesoscale structures in the transition zone: Dynamical consequences of boundary layer activities , 1998 .

[35]  P. Tackley,et al.  Generation of mega‐plumes from the core‐mantle boundary in a compressible mantle with temperature‐dependent viscosity , 1998 .

[36]  Don L. Anderson,et al.  The thermal state of the upper mantle; No role for mantle plumes , 2000 .

[37]  M. Sylvander,et al.  Seismic velocities at the core-mantle boundary inferred from P waves diffracted around the core , 1997 .

[38]  Robert W. Clayton,et al.  Lower mantle heterogeneity, dynamic topography and the geoid , 1985, Nature.

[39]  Yu.P. Orovetskii Mantle Plumes , 2020 .

[40]  P. Shearer,et al.  Global mapping of topography on transition zone velocity discontinuities by stacking SS precursors , 1998 .

[41]  Tianyu Zheng,et al.  A complex 660 km discontinuity beneath northeast China , 2003 .

[42]  D. Yuen,et al.  On the possibility of a second kind of mantle plume , 2000 .

[43]  Williams,et al.  A correlation between ultra-Low basal velocities in the mantle and hot spots , 1998, Science.

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

[45]  N. Sleep Lateral flow of hot plume material ponded at sublithospheric depths , 1996 .

[46]  McSween Hy,et al.  Evidence for Life in a Martian Meteorite , 1997 .

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

[48]  Hua-wei Zhou,et al.  A high-resolution P wave model for the top 1200 km of the mantle , 1996 .

[49]  B. Kennett,et al.  Traveltimes for global earthquake location and phase identification , 1991 .

[50]  Harmen Bijwaard,et al.  Tethyan subducted slabs under India , 1999 .

[51]  David A. Yuen,et al.  Lower mantle thermal structure deduced from seismic tomography, mineral physics and numerical modelling , 1994 .

[52]  W. J. Morgan,et al.  The seismic anomaly beneath Iceland extends down to the mantle transition zone and no deeper , 2000 .

[53]  K. Condie Mantle Plumes and Their Record in Earth History: Superplume Events , 2001 .

[54]  Sean C. Solomon,et al.  Seismic structure of the Iceland mantle plume , 1997, Nature.

[55]  Frank D. Stacey,et al.  The thermal boundary-layer interpretation of D″ and its role as a plume source , 1983 .

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

[57]  U. Christensen,et al.  The excess temperature of plumes rising from the core‐mantle boundary , 1996 .

[58]  C. Young,et al.  The Core-Mantle Boundary , 1987 .

[59]  H. Nataf,et al.  Seismic Imaging of Mantle Plumes , 2000 .

[60]  R. Jeanloz,et al.  Earth's Core-Mantle Boundary: Results of Experiments at High Pressures and Temperatures , 1991, Science.

[61]  L. Kellogg,et al.  Effect of mantle plumes on the growth of D” by reaction between the core and mantle , 1993 .

[62]  B. Romanowicz,et al.  Anisotropic structures at the base of the Earth's mantle , 1998, Nature.

[63]  S. Balachandar,et al.  Three-Dimensional Instabilities of Mantle Convection with Multiple Phase Transitions , 1993, Science.

[64]  B. Steinberger Plumes in a convecting mantle: Models and observations for individual hotspots , 2000 .

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

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

[67]  M. Richards,et al.  The adjustment of mantle plumes to changes in plate motion , 1989 .

[68]  T. Lay,et al.  Small‐scale lateral shear velocity and anisotropy heterogeneity near the core‐mantle boundary beneath the central Pacific imaged using broadband ScS waves , 1999 .

[69]  B. Kennett,et al.  ELLIPTICITY CORRECTIONS FOR SEISMIC PHASES , 1996 .

[70]  T. Jordan,et al.  Slab penetration into the lower mantle beneath the Mariana and other island arcs of the northwest Pacific , 1986 .

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

[72]  S. Solomon,et al.  Inversion of body-wave delay times for mantle structure beneath the Hawaiian islands: results from the PELENET experiment , 2002 .

[73]  B. Hager,et al.  — a H 03 Dynamically Supported Geoid Highs over Hotspots : Observation and Theory , 2008 .

[74]  David A. Yuen,et al.  Numerical simulations of thermal-chemical instabilities at the core–mantle boundary , 1988, Nature.

[75]  John A. Hildebrand,et al.  Depth Extent of the Lau Back-Arc Spreading Center and Its Relation to Subduction Processes , 1997 .

[76]  Patrice Weber,et al.  Intermittent layered convection in a model mantle with an endothermic phase change at 670 km , 1991, Nature.

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

[78]  Toshiro Tanimoto,et al.  Predominance of Large-Scale Heterogeneity and the Shift of Velocity Anomalies between the Upper and Lower Mantle , 1990 .

[79]  W. Spakman,et al.  Mesozoic subducted slabs under Siberia , 1999, Nature.

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

[81]  D. Zhao,et al.  Seismological structure of subduction zones and its implications for arc magmatism , 2001 .

[82]  T. Jordan,et al.  Slab penetration into the lower mantle , 1984 .

[83]  H. Nataf,et al.  Detection of mantle plumes in the lower mantle by diffraction tomography: Hawaii , 1998 .

[84]  Harmen Bijwaard,et al.  Tomographic evidence for a narrow whole mantle plume below Iceland , 1999 .

[85]  Masayuki Obayashi,et al.  Stagnant slabs in the upper and lower mantle transition region , 2001 .

[86]  Gabi Laske,et al.  CRUST 5.1: A global crustal model at 5° × 5° , 1998 .

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

[88]  W. J. Morgan,et al.  Seismic evidence for a tilted mantle plume and north–south mantle flow beneath Iceland , 2002 .

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

[90]  S. Karato,et al.  Importance of anelasticity in the interpretation of seismic tomography , 1993 .

[91]  T. Lay,et al.  Seismic anisotropy in the lowermost mantle beneath the Pacific , 1998 .

[92]  J. Ritsema,et al.  African hot spot volcanism: small-scale convection in the upper mantle beneath cratons. , 2000, Science.

[93]  Lapo Boschi,et al.  High‐ and low‐resolution images of the Earth's mantle: Implications of different approaches to tomographic modeling , 1999 .

[94]  Shigenori Maruyama,et al.  Mechanism of backarc opening in the Japan Sea: role of asthenospheric injection , 1990 .

[95]  K. Abe,et al.  Hypocenter Determination Procedure of Japan University Network Earthquake Catalog , 1989 .

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

[97]  J. Tuzo Wilson A possible origin of the Hawaiian Islands , 1963 .

[98]  S. Tanaka Very low shear wave velocity at the base of the mantle under the South Pacific Superswell , 2002 .

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

[100]  K. Priestley,et al.  Mapping the Hawaiian plume conduit with converted seismic waves , 2000, Nature.