Reconstructing slab dip through deep time to explain pulses in kimberlite eruptions

The recycling of oceanic lithosphere at subduction zones constitutes the largest driving force of plate tectonic motion. The angle at which subducting plates enter the mantle influences the magnitude of this force, the distribution of subduction-related earthquakes, intensity of volcanism, and mountain building. However, the factors that control subduction angle remain unresolved. We develop a novel formulation for calculating the subduction angle based on trench migration, convergence rate, slab thickness, and plate density which reproduces the present-day dynamics of global subduction zones. Applying this formulation to reconstructed subduction boundaries from the Jurassic to present day, we relate subduction angle combined with slab flux to pulses in kimberlite eruptions. High rates of subducting slab material trigger mantle return flow that stimulates fertile reservoirs in the mantle. These convective instabilities transport slab-influenced melt to the surface at a distance inbound from the trench corresponding to the subduction angle. Our deep-time slab dip formulation has numerous potential applications including modelling the deep carbon and water cycles, and an improved understanding of subduction-related mineral deposits.