Convection of a fluid with strongly temperature and pressure dependent viscosity

Plate tectonics on the Earth is a surface manifestation of convection within the Earth’s mantle, a subject which is as yet improperly understood, and it has motivated the study of various forms of buoyancy-driven thermal convection. The early success of the high Rayleigh number constant viscosity theory was later tempered by the absence of plate motion when the viscosity is more realistically strongly temperature dependent, and the process of subduction represents a continuing principal conundrum in the application of convection theory to the Earth. A similar problem appears to arise if the equally strong pressure dependence of viscosity is considered, since the classical isothermal core convection theory would then imply a strongly variable viscosity in the convective core, which is inconsistent with results from post-glacial rebound studies. In this paper we address the problem of determining the asymptotic structure of high Rayleigh number convection when the viscosity is strongly temperature and pressure dependent, i.e. thermobaroviscous. By a method akin to lid-stripping, we are able to extend numerical computations to extremely high viscosity contrasts, and we show that the convective cells take the form of narrow, vertically-oriented fingers. We are then able to determine the asymptotic structure of the solution, and it agrees well with the numerical results. Beneath a stagnant lid, there is a vigorous convection in the upper part of the cell, and a more sluggish, higher viscosity flow in the lower part of the cell. We then offer some comments on the possible meaning and interpretation of these results for planetary mantle convection.