Subducting slabs: Jellyfishes in the Earth's mantle

The constantly improving resolution of geophysical data, seismic tomography and seismicity in particular, shows that the lithosphere does not subduct as a slab of uniform thickness but is rather thinned in the upper mantle and thickened around the transition zone between the upper and lower mantle. This observation has traditionally been interpreted as evidence for the buckling and piling of slabs at the boundary between the upper and lower mantle, where a strong contrast in viscosity may exist and cause resistance to the penetration of slabs into the lower mantle. The distribution and character of seismicity reveal, however, that slabs undergo vertical extension in the upper mantle and compression near the transition zone. In this paper, we demonstrate that during the subduction process, the shape of low viscosity slabs (1 to 100 times more viscous than the surrounding mantle) evolves toward an inverted plume shape that we coin jellyfish. Results of a 3D numerical model show that the leading tip of slabs deform toward a rounded head skirted by lateral tentacles that emerge from the sides of the jellyfish head. The head is linked to the body of the subducting slab by a thin tail. A complete parametric study reveals that subducting slabs may achieve a variety of shapes, in good agreement with the diversity of natural slab shapes evidenced by seismic tomography. Our work also suggests that the slab to mantle viscosity ratio in the Earth is most likely to be lower than 100. However, the sensitivity of slab shapes to upper and lower mantle viscosities and densities, which remain poorly constrained by independent evidence, precludes any systematic deciphering of the observations.

[1]  Olafur Gudmundsson,et al.  A regionalized upper mantle (RUM) seismic model , 1998 .

[2]  Chris Kincaid,et al.  An experimental study of subduction and slab migration , 1987 .

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

[4]  Hua-Wei Zhou.  How well can we resolve the deep seismic slab with seismic tomography , 1988 .

[5]  Eleonore Stutzmann,et al.  Understanding seismic heterogeneities in the lower mantle beneath the Americas from seismic tomography and plate tectonic history , 2007 .

[6]  M. Bevis.  The curvature of Wadati-Benioff zones and the torsional rigidity of subducting plates , 1986, Nature.

[7]  M. Richards,et al.  A geodynamic model of mantle density heterogeneity , 1993 .

[8]  C. Faccenna,et al.  A Review of the Role of Subduction Dynamics for Regional and Global Plate Motions , 2009 .

[9]  R. W. Griffiths,et al.  Laboratory models of the thermal evolution of the mantle during rollback subduction , 2003, Nature.

[10]  D. Yuen,et al.  The interaction of a subducting lithospheric slab with a chemical or phase boundary , 1984 .

[11]  H. Kárason.  Constraints on mantle convection from seismic tomography and flow modeling , 2002 .

[12]  Shijie Zhong,et al.  Effects of plate and slab viscosities on the geoid , 1999 .

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

[14]  Y. Fukao,et al.  Spherical Shell Tectonics: Effects of Sphericity and Inextensibility on the Geometry of the Descending Lithosphere (Paper 5R0711) , 1986 .

[15]  T. Lay.  The Fate of Descending Slabs , 1994 .

[16]  H. Čížková,et al.  Stress distribution within subducting slabs and their deformation in the transition zone , 2007 .

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

[18]  Shijie Zhong,et al.  Controls on trench topography from dynamic models of subducted slabs , 1994 .

[19]  C. Thieulot,et al.  DOUAR: A new three-dimensional creeping flow numerical model for the solution of geological problems , 2008 .

[20]  W. Schellart.  Evolution of the slab bending radius and the bending dissipation in three-dimensional subduction models with a variable slab to upper mantle viscosity ratio , 2009 .

[21]  L. Royden,et al.  Trench motion, slab geometry and viscous stresses in subduction systems , 2006 .

[22]  G. Morra,et al.  Evidence of lower-mantle slab penetration phases in plate motions , 2008, Nature.

[23]  Paul J. Tackley,et al.  Effects of strongly temperature‐dependent viscosity on time‐dependent, three‐dimensional models of mantle convection , 1993 .

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

[25]  B. Evans,et al.  Strength of the lithosphere: Constraints imposed by laboratory experiments , 1995 .

[26]  R. O’Connell,et al.  Deformation of a weak subducted slab and variation of seismicity with depth , 1993, Nature.

[27]  G. Morra,et al.  Dynamic models of downgoing plate-buoyancy driven subduction: subduction motions and energy dissipation , 2007 .

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

[29]  Louis-Noel Moresi,et al.  Influence of trench width on subduction hinge retreat rates in 3‐D models of slab rollback , 2006 .

[30]  U. Christensen.  The influence of trench migration on slab penetration into the lower mantle , 1996 .

[31]  R. Jarrard.  Relations among subduction parameters , 1986 .

[32]  M. Gurnis,et al.  Constraints on the lateral strength of slabs from three-dimensional dynamic flow models , 1996 .

[33]  W. Peltier.  Mantle Viscosity and Ice-Age Ice Sheet Topography , 1996, Science.

[34]  Robert D. van der Hilst,et al.  Travel-time tomography of the European-Mediterranean mantle down to 1400 km , 1993 .

[35]  R. Hilst,et al.  Buckling instabilities of subducted lithosphere beneath the transition zone , 2007 .

[36]  C. Thieulot,et al.  Influence of surrounding plates on 3D subduction dynamics , 2009 .

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

[38]  D. Giardini,et al.  On the Curvature of Oceanic Arcs , 2003 .

[39]  L. Husson.  Dynamic topography above retreating subduction zones , 2006 .

[40]  D. Gubbins,et al.  DEFORMATION OF SUBDUCTED OCEANIC LITHOSPHERE , 1997 .

[41]  R. Dietmar Müller,et al.  Digital isochrons of the world's ocean floor , 1997 .

[42]  Manik Talwani,et al.  Island arcs, deep sea trenches, and back-arc basins , 1977 .

[43]  W. Spakman,et al.  Subduction and slab detachment in the Mediterranean-Carpathian region. , 2000, Science.

[44]  M. Gurnis,et al.  Constraints on subducting plate strength within the Kermadec trench , 2005 .

[45]  Hua-Wei Zhou,et al.  Slope of the geoid spectrum and constraints on mantle viscosity stratification , 1996 .

[46]  T. Becker,et al.  The dynamics of subduction and trench migration for viscosity stratification , 2005 .

[47]  J. Mitrovica,et al.  Lateral variations in mantle rheology: implications for convection related surface observables and inferred viscosity models , 2007 .

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

[49]  S. Lallemand,et al.  Plate motions, slab dynamics and back-arc deformation , 2005 .

[50]  C. Faccenna,et al.  Three‐dimensional instantaneous mantle flow induced by subduction , 2006 .

[51]  C. Faccenna,et al.  Slab stiffness control of trench motion: Insights from numerical models , 2008 .

[52]  R. Hilst,et al.  A laboratory investigation of effects of trench migration on the descent of subducted slabs , 1995 .

[53]  K. Shimazaki,et al.  Three-dimensional P-wave velocity structure beneath the Indonesian region , 1993 .

[54]  R. Hilst,et al.  Tomographic imaging of the lowermost mantle with differential times of refracted and diffracted core phases (PKP, Pdiff) , 2001 .

[55]  M. Wortel,et al.  Subduction zone seismicity and the thermo-mechanical evolution of downgoing lithosphere , 1988 .

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

[57]  Lapo Boschi,et al.  A comparison of tomographic and geodynamic mantle models , 2002 .

[58]  B. Hager,et al.  Effects of plate bending and fault strength at subduction zones on plate dynamics , 1999 .

[59]  S. Grand,et al.  Evidence of high velocity anomalies in the transition zone associated with Southern Kurile Subduction Zone , 1995 .

[60]  Clinton P. Conrad,et al.  Global mantle flow and the development of seismic anisotropy: Differences between the oceanic and continental upper mantle , 2007 .

[61]  Andrea Morelli,et al.  P wave tomography of the mantle under the Alpine-Mediterranean area , 2003 .

[62]  P. Molnar,et al.  Mantle Earthquake Mechanisms and the Sinking of the Lithosphere , 1969, Nature.

[63]  W. Schellart.  Kinematics of subduction and subduction-induced flow in the upper mantle , 2004 .

[64]  M. Wortel,et al.  Trench migration and subduction zone geometry , 1997 .

[65]  Arie P. van den Berg,et al.  The influence of rheological weakening and yield stress on the interaction of slabs with the 670 km discontinuity , 2002 .

[66]  D. May,et al.  Evolution and diversity of subduction zones controlled by slab width , 2007, Nature.

[67]  D. Giardini,et al.  Dynamics of retreating slabs : 2. Insights from three-dimensional laboratory experiments , 2003 .

[68]  P. Olson,et al.  Laboratory experiments on the structure of subducted lithosphere , 1995 .

[69]  P. Molnar,et al.  Distribution of stresses in the descending lithosphere from a global survey of focal‐mechanism solutions of mantle earthquakes , 1971 .

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

[71]  Ross W. Griffiths,et al.  Folding of Viscous Plumes Impinging On A Density Or Viscosity Interface , 1988 .

[72]  Magali I. Billen,et al.  Modeling the Dynamics of Subducting Slabs , 2008 .

[73]  F. Tajima,et al.  Dynamical modeling of trench retreat driven by the slab interaction with the mantle transition zone , 2007 .

[74]  B. Hager,et al.  Long-wavelength variations in Earth’s geoid: physical models and dynamical implications , 1989, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[75]  Bradford H. Hager,et al.  The distribution of earthquakes with depth and stress in subducting slabs , 1984 .

[76]  J. Mitrovica,et al.  A new inference of mantle viscosity based upon joint inversion of convection and glacial isostatic adjustment data , 2004 .

[77]  S. Yoshioka,et al.  Three-dimensional numerical modeling of detachment of subducted lithosphere , 1995 .

[78]  T. Jordan,et al.  A velocity anomaly in the lower mantle , 1974 .

[79]  H. Čížková,et al.  Long-wavelength character of subducted slabs in the lower mantle , 2008 .

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

[81]  B. Hager,et al.  Compositional vs. thermal buoyancy and the evolution of subducted lithosphere , 1994 .

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

[83]  E. Engdahl,et al.  A new global model for P wave speed variations in Earth's mantle , 2008 .

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

[85]  G. Davies.  Penetration of plates and plumes through the mantle transition zone , 1995 .

[86]  Cedric Thieulot,et al.  Adaptive octree-based finite element analysis of two- and three-dimensional indentation problems , 2008 .

[87]  M. Gurnis,et al.  Mantle Convection with Plates and Mobile, Faulted Plate Margins , 1995, Science.

[88]  C. Conrad,et al.  Influence of dynamic topography on sea level and its rate of change , 2009 .

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

[90]  C. Frohlich,et al.  The relationship between Wadati‐Benioff Zone geometry and P, T and B axes of intermediate and deep focus earthquakes , 1987 .

[91]  Sri Widiyantoro,et al.  Structure and Evolution of Lithospheric Slab Beneath the Sunda Arc, Indonesia , 1996, Science.

[92]  T. Becker.  On the effect of temperature and strain-rate dependent viscosity on global mantle flow, net rotation, and plate-driving forces , 2006 .

[93]  Y. Fukao.  Seismic Tomogram of the Earth's Mantle: Geodynamic Implications , 1992, Science.

[94]  K. Lambeck,et al.  The Earth's Mantle: The Viscosity of the Mantle: Evidence from Analyses of Glacial-Rebound Phenomena , 1998 .

[95]  M. Kido,et al.  Inferences of viscosity from the oceanic geoid: Indication of a low viscosity zone below the 660-km discontinuity , 1997 .

[96]  C. Faccenna,et al.  Trench migration, net rotation and slab-mantle coupling , 2008 .

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

[98]  R. Müller,et al.  Episodicity in back-arc tectonic regimes , 2008 .

[99]  N. Ribe.  Periodic folding of viscous sheets. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[100]  D. Giardini,et al.  Deep seismicity and modes of deformation in Tonga subduction zone , 1984, Nature.

[101]  J. Braun,et al.  From longitudinal slab curvature to slab rheology , 2009 .

[102]  N. Loubet,et al.  Deformation modes of subducted lithosphere at the core‐mantle boundary: An experimental investigation , 2009 .

[103]  D. Giardini,et al.  Episodic Back-arc extension during restricted mantle convection in the Central Mediterranean , 2001 .