The relationship between plate velocity and trench viscosity in Newtonian and power‐law subduction calculations

Convection with a Newtonian temperature-dependent rheology leads to little or no surface velocity unless zones of weakness are introduced. “Plate-like” features are observed in calculations both with Newtonian rheology, employing imposed weak zones, and with power-law (non-Newtonian) rheology, where high stresses at the trench reduce the effective viscosity. Since deformation at subduction zones involves faulting, both of these parameterizations should be treated with some skepticism. It is important to understand how the parameterizations affect the model results. We study the relationship between trench viscosity and plate velocity using a Newtonian rheology by varying the viscosity at the trench. The plate velocity is a function of the trench viscosity for fixed Rayleigh number and plate/slab viscosity. Slab velocities for non-Newtonian rheology calculations are significantly different from slab velocities from Newtonian rheology calculations at the same effective Rayleigh number. Both models give reasonable strain-rates for the slab when compared with estimates of seismic strain-rate. Non-Newtonian rheology eliminates the need for imposed weak zones and provides a self-consistent fluid dynamical mechanism for subduction in numerical convection models.

[1]  M. Bevis.  Seismic Slip and Down-Dip Strain Rates in Wadati-Benioff Zones , 1988, Science.

[2]  G. Schubert,et al.  A benchmark comparison of numerical methods for infinite Prandtl number thermal convection in two-dimensional Cartesian geometry , 1990 .

[3]  Bradford H. Hager,et al.  A simple global model of plate dynamics and mantle convection , 1981 .

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

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

[6]  Prame Chopra,et al.  The experimental deformation of dunite , 1981 .

[7]  Bradford H. Hager,et al.  Conman: vectorizing a finite element code for incompressible two-dimensional convection in the Earth's mantle , 1990 .

[8]  L. Cserepes.  Numerical studies of non-Newtonian mantle convection , 1982 .

[9]  Wolfgang R. Jacoby,et al.  On modelling the lithosphere in mantle convection with non-linear rheology , 1982 .

[10]  Ulrich R. Christensen,et al.  Convection in a variable-viscosity fluid: Newtonian versus power-law rheology , 1983 .

[11]  E. Boschi,et al.  Mantle rheology from a geodynamical standpoint , 1982 .

[12]  B. Hager,et al.  Subduction zone earthquakes and stress in slabs , 1988 .

[13]  M. Ashby,et al.  Micromechanisms of flow and fracture, and their relevance to the rheology of the upper mantle , 1978, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[14]  Ulrich R. Christensen,et al.  Convection with pressure- and temperature-dependent non-Newtonian rheology , 1984 .

[15]  Frank M. Richter,et al.  Convection experiments in fluids with highly temperature-dependent viscosity and the thermal evolution of the planets , 1982 .

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

[17]  G. Davies.  Mantle convection under simulated plates: effects of heating modes and ridge and trench migration, and implications for the core—mantle boundary, bathymetry, the geoid and Benioff zones , 1986 .

[18]  D. Yuen,et al.  Time‐dependent convection with non‐Newtonian viscosity , 1989 .

[19]  D. Turcotte,et al.  Studies of finite amplitude non‐Newtonian thermal convection with application to convection in the Earth's mantle , 1976 .

[20]  B. Hager.  Subducted slabs and the geoid: Constraints on mantle rheology and flow , 1983 .

[21]  F. Richter,et al.  Simple plate models of mantle convection , 1977 .

[22]  F. Richter.  Focal mechanisms and seismic energy release of deep and intermediate earthquakes in the Tonga‐Kermadec Region and their bearing on the depth extent of mantle flow , 1979 .

[23]  D. Yuen,et al.  Chaotic axisymmetrical spherical convection and large-scale mantle circulation , 1987 .