A simple model of convection in the terrestrial planets

Abstract

In this paper we study analytically the simplest fluid mechanical model which can mimic the convective behavior which is thought to occur in the solid mantles of the terrestrial planets. The convecting materials are polycrystalline rocks, whose creep behavior depends very strongly on temperature and probably also on pressure. As a simple model of this situation, we consider the flow of a Newtonian viscous fluid, whose viscosity depends strongly on temperature (only), and in fact has an infinite viscosity below a certain temperature, and a constant viscosity above this temperature. This model would also be directly relevant to the convection of a melt beneath its own solid phase (e.g. water below ice, though in that case there are other physical complications).

As a consequence of this assumption, there is a vigorous convection zone overlain by a stagnant lid, as also observed in analogous laboratory experiments (Nataf and Richter, 1982). The analysis is then very similar to that of Roberts (1979), but the extension to variable viscosity introduces important differences, most notably that the boundary between the lid and the convecting zone is unknown, and not horizontal. The resulting buoyancy induced stresses near this boundary are much larger than the stresses produced by buoyancy in the side-wall plumes, and mean that the dynamics of this region, and hence also the heat flux, are independent of the rest of the cell. We give a first order approximation for the Nusselt number-Rayleigh number relationship.