A linear stability analysis on the onset of thermal convection of a fluid with strongly temperature-dependent viscosity in a spherical shell

A linear stability analysis was performed in order to study the onset of thermal convection in the presence of a strong viscosity variation, with a special emphasis on the condition for the stagnant-lid (ST) convection where a convection takes place only in a sublayer beneath a highly viscous lid of cold fluid. We consider the temporal evolution (growth or decay) of an infinitesimal perturbation superimposed to a Boussinesq fluid with an infinite Prandtl number which is in a static (motionless) and conductive state in a basally heated planar layer or spherical shell. The viscosity of the fluid is assumed to be exponentially dependent on temperature. The linearized equations for conservations of mass, momentum, and internal (thermal) energy are numerically solved for the critical Rayleigh number, Rac, as well as the radial profiles of eigenfunctions for infinitesimal perturbations. The above calculations are repeatedly carried out by systematically varying (i) the magnitude of the temperature dependence of viscosity, E, and (ii) the ratio of the inner and outer radii of the spherical shell, γ. A careful analysis of the vertical structure of incipient flows demonstrated that the dominance of the ST convection can be quantitatively identified by the vertical profile of Δh (a measure of conversion between horizontal and vertical flows), regardless of the model geometries. We also found that, in the spherical shell relevant to the Earth’s mantle (γ = 0.55), the transition into ST convection takes place at the viscosity contrast across the layer \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${r_\eta\simeq10^4}$$\end{document} . Taken together with the fact that the threshold value of rη falls in the range of rη for a so-called sluggish-lid convection, our finding suggests that the ST-mode of convection with horizontally elongated convection cells is likely to arise in the Earth’s mantle solely from the temperature-dependent viscosity.