GEOPHYSICAL RESEARCH LETTERS, VOL. 17, NO. 13, PAGES 2409–2412, 1990

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

Scott D. King

Seismological Laboratory, Caltech

Bradford H. Hagar

Department of Earth, Atmospheric and Planetary Science, MIT

Abstract

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-newtonain rheology elimates the need for imposed weak zones and provides a self-consistent fluid dynamical mechanism for subduction in numerical convection models. ©American Geophysical Union 1990

Index Terms: 8155 Tectonophysics: Plate motions, past and present; 8160 Tectonophysics: Rheology of the lithosphere and mantle.

Citation: King, S. D., and B. H. Hagar (1990), The relationship between plate velocity and trench viscosity in Newtonian and power-law subduction calculations, Geophys. Res. Lett., 17(13), 2409–2412.

Copyright 1990 by the American Geophysical Union.