Abstract
We present the results of numerical mantle convection models demonstrating that dynamical effects induced by variable mantle viscosity, depth-dependent thermal expansivity, radiative thermal conductivity at the base of the mantle, the spinel to perovskite phase change and the perovskite to post-perovskite phase transition in the deep mantle can result in multiscale mantle plumes: stable lower-mantle superplumes are followed by groups of small upper-mantle plumes. Both radiative thermal conductivity at the base of the lower mantle and a strongly decreasing thermal expansivity of perovskite in the lower mantle can help induce partially layered convection with intense shear heating under the transition zone, which creates a low-viscosity zone and allows for the production of secondary mantle plumes emanating from this zone.