Lithospheric stresses in Rayleigh–Bénard convection: effects of a free surface and a viscoelastic Maxwell rheology

Numerical models of mantle convection typically employ a temperature- or pressure-dependent viscous or viscoplastic rheology and a free slip upper boundary condition. The Earth, however, has a stress-free rather than a free slip surface condition. In addition, with decreasing temperature, the viscosity of rocks increases, which might induce a change from viscous to elastic behaviour (depending on the timescale of deformation). Here, we study the...