Stirring geochemistry in mantle convection models with stiff plates and slabs

Numerical models of mantle convection are presented that readily yield midocean ridge basalt (MORB) and oceanic island basalt (OIB) ages equaling or exceeding the apparent ∼1.8-Ga lead isotopic ages of trace-element heterogeneities in the mantle. These models feature high-viscosity surface plates and subducting lithosphere, and higher viscosities in the lower mantle. The formation and subduction of oceanic crust are simulated by means of tracers that represent a basaltic component. The models are run at the full mantle Rayleigh number and take account of faster mantle overturning and deeper melting in the past. More than 97% of the mantle is processed in these models. Including the expected excess density of former oceanic crust readily accounts for the depletion of MORB source relative to OIB sources. A novel finding is of gravitational settling of dense tracers within the low-viscosity upper mantle, as well as at the base of the mantle. The models suggest as well that the seismological observation of a change in tomographic character in the deep mantle might be explained without the need to postulate a separate layer in the deep mantle. These results expand the range of models with the potential to reconcile geochemical and geophysical observations of the mantle.