Abstract:
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.
