The influence of trench migration on slab penetration into the lower mantle

A two-dimensional numerical convection model in cartesian geometry is used to study the influence of trench migration on the ability of subducted slabs to penetrate an endothermic phase boundary at 660 km depth. The transient subduction history of an oceanic plate is modelled by imposing plate and trench motion at the surface. The viscosity depends on temperature and depth. A variety of styles of slab behaviour is found, depending predominantly on the trench velocity. When trench retreat is faster than 2–4 cm/a, the descending slab flattens above the phase boundary. At slower rates it penetrates straight into the lower mantle, although flattening in the transition zone may occur later, leading to a complex slab morphology. The slab can buckle, independent of whether it penetrates or not, especially when there is a localised increase in viscosity at the phase boundary. Flattened slabs are only temporarily arrested in the transition zone and sink ultimately into the lower mantle. The results offer a framework for understanding the variety in slab geometry revealed by seismic tomography.