TY - JOUR

T1 - Dynamics of strongly time-dependent convection with non-Newtonian temperature-dependent viscosity

AU - Larsen, Tine B.

AU - Yuen, David A.

AU - Malevsky, Andrei V.

AU - Smedsmo, Jamie L.

N1 - Funding Information:
research has been supported by Department of Energy, NSF Geochemistry Program, Army High Performance Computing and Research Center

PY - 1996/3

Y1 - 1996/3

N2 - Numerical simulations of time-dependent Boussinesq convection with a non-Newtonian temperature-dependent rheology have been conducted in wide aspect-ratio boxes for effective Rayleigh numbers, Rav, ranging from 3.5 × 104 to 1.3 × 106. The transition to the chaotic regime for non-Newtonian temperature-dependent viscosity takes place at Rav an order of magnitude lower than for purely temperature-dependent viscosity. For viscosity contrasts owing to temperature of O(102), lateral viscosity contrasts generated by the flow are typically of O(103) in the interior with contrasts as large as O(106) in the upper boundary layer. The horizontal Fourier spectra of the temperature field and of the inverse of the viscosity field show that the asymmetry in the flow introduced by a temperature-dependent viscosity is weaker for a non-Newtonian rheology than for a Newtonian rheology. Lithospheric dynamics, such as lubrication of subduction and lithospheric thinning, are greatly facilitated by the effects of stress-softening from non-Newtonian rheology. Weak plumes are slowed down drastically in the middle of their ascent because of the stress-dependent nature of the rheology. Downgoing slabs can be found to be weakened near the middle, and detachment involving pieces of the descending material may take place. The mobility of the stiff upper boundary layer is facilitated by an increase in the vigor of convection, whereas an increase in the viscosity contrast owing to temperature retards the surface mobility.

AB - Numerical simulations of time-dependent Boussinesq convection with a non-Newtonian temperature-dependent rheology have been conducted in wide aspect-ratio boxes for effective Rayleigh numbers, Rav, ranging from 3.5 × 104 to 1.3 × 106. The transition to the chaotic regime for non-Newtonian temperature-dependent viscosity takes place at Rav an order of magnitude lower than for purely temperature-dependent viscosity. For viscosity contrasts owing to temperature of O(102), lateral viscosity contrasts generated by the flow are typically of O(103) in the interior with contrasts as large as O(106) in the upper boundary layer. The horizontal Fourier spectra of the temperature field and of the inverse of the viscosity field show that the asymmetry in the flow introduced by a temperature-dependent viscosity is weaker for a non-Newtonian rheology than for a Newtonian rheology. Lithospheric dynamics, such as lubrication of subduction and lithospheric thinning, are greatly facilitated by the effects of stress-softening from non-Newtonian rheology. Weak plumes are slowed down drastically in the middle of their ascent because of the stress-dependent nature of the rheology. Downgoing slabs can be found to be weakened near the middle, and detachment involving pieces of the descending material may take place. The mobility of the stiff upper boundary layer is facilitated by an increase in the vigor of convection, whereas an increase in the viscosity contrast owing to temperature retards the surface mobility.

UR - http://www.scopus.com/inward/record.url?scp=0029730795&partnerID=8YFLogxK

U2 - 10.1016/0031-9201(95)03082-4

DO - 10.1016/0031-9201(95)03082-4

M3 - Article

AN - SCOPUS:0029730795

SN - 0031-9201

VL - 94

SP - 75

EP - 103

JO - Physics of the Earth and Planetary Interiors

JF - Physics of the Earth and Planetary Interiors

IS - 1-2

ER -