A boundary layer model for mantle convection with surface plates

Summary. A simple, analytical model for mantle convection with mobile surface plates is presented. Our aim is to determine under what conditions free convection can account for the observed plate motions, and to evaluate the thermal structure of the mantle existing under these conditions. Boundary layer methods are used to represent two-dimensional cellular convection at large Rayleigh and infinite Prandtl numbers. The steady-state structure consists of cells with isentropic interiors enclosed by thermal boundary layers. Lithospheric plates are represented as upper surfaces on each cell free to move at a uniform speed. Buoyancy forces are concentrated in narrow rising and decending thermal plumes; torques imparted by these plumes drive both the deformable mantle and overlying plate. Solutions are found for a comprehensive range of cell sizes. We derive an expression for the plate speed as a function of its length, the mantle viscosity and surface heat flux. Using mean values for these parameters, we find that thermal convection extending to 700 km depth can move plates at 1 cm yr-’, while convection through the whole mantle can move plates at 4-5 cm yr-’. Analysis of the steady-state temperature field, for the case of heating from below, shows that the upper thermal boundary layer develops a complex structure, including an ‘asthenosphere’ defined by a local maximum in the geotherm occurring at depths of 50-1 50 km.