The influence of rheological weakening and yield stress on the interaction of slabs with the 670 km discontinuity

Results of high resolution seismic tomography showing subducting slabs deflected in the transition zone and thickened in the lower mantle seem to call for slab material weaker than inferred from mineral physics deformation mechanisms. A possible mechanism suggested by several authors could be the weakening due to grain size reduction, which should occur in the cold portion of fast slabs after an exothermic phase transition at a depth of 400 km. Since the amount of weakening as well as the rate of subsequent strengthening due to the grain growth are not precisely known, we present here a parametric study of slab behavior in the transition zone and upper part of the lower mantle. We simulate a subducting slab in a two-dimensional (2-D) Cartesian box in the numerical model with composite rheology including diffusion creep, dislocation creep and a general stress limiting rheology approximating Peierl’s creep. We concentrate on two rheologic effects: the dynamic effect of slab weakening due to grain size reduction at the phase boundary and the effect of yield stress of stress limiting rheology. The effect of trench migration on slab deformation is also included in our study. Results show that the slab ability to penetrate into the lower mantle is not significantly affected by a trench retreat in the absence of grain size weakening. However, in case of a 4 cm/yr trench retreat, grain size weakening provides a viable mechanism to deflect the slab in the transition zone, provided that stress limiting deformation mechanism would limit the effective viscosity outside the areas of grain size weakening to about 10 24 Pa s. < 2002 Elsevier Science B.V. All rights reserved.

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

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

[3]  Patrick Wu,et al.  Rheology of the Upper Mantle: A Synthesis , 1993, Science.

[4]  D. Olbertz.  The long-term evolution of subduction zones : a modelling study , 1997 .

[5]  Don L. Anderson,et al.  The upper mantle transition region - Eclogite , 1979 .

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

[7]  M. Riedel,et al.  Grain-size evolution in subducted oceanic lithosphere associated with the olivine-spinel transformation and its effects on rheology , 1997 .

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

[9]  D. Yuen,et al.  Rheological structure and deformation of subducted slabs in the mantle transition zone: implications for mantle circulation and deep earthquakes , 2001 .

[10]  S. Karato.  On the separation of crustal component from subducted oceanic lithosphere near the 660 KM discontinuity , 1997 .

[11]  P. Vaughan,et al.  Creep mechanism in Mg2GeO4: Effects of a phase transition , 1981 .

[12]  Charles Anderson,et al.  Formation of a metastable olivine wedge in a descending slab , 1997 .

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

[14]  D. Yuen,et al.  Rheological control of oceanic crust separation in the transition zone , 1996 .

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

[16]  D. Yuen,et al.  Thermal-mechanical effects of low-temperature plasticity (the Peierls mechanism) on the deformation of a viscoelastic shear zone , 1999 .

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

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

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

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

[21]  Arie P. van den Berg,et al.  Latent heat effects of the major mantle phase transitions on low-angle subduction , 2001 .

[22]  A. Berg,et al.  A thermo-mechanical model of horizontal subduction below an overriding plate , 2000 .

[23]  David A. Yuen,et al.  The effects of a composite non-Newtonian and Newtonian rheology on mantle convection , 1993 .

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

[25]  A. E. Ringwood,et al.  Role of the transition zone and 660 km discontinuity in mantle dynamics , 1994 .

[26]  A. E. Ringwood,et al.  Phase transformations in subducted oceanic crust and buoyancy relationships at depths of 600–800 km in the mantle , 1993 .

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

[28]  D. Rubie,et al.  The olivine → spinel transformation and the rheology of subducting lithosphere , 1984, Nature.

[29]  S. King,et al.  Sensitivity of convection with an endothermic phase change to the form of governing equations, initial conditions, boundary conditions, and equation of state , 1994 .

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

[31]  David A. Yuen,et al.  The metastable olivine wedge in fast subducting slabs: Constraints from thermo-kinetic coupling , 1996 .

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