Abstract. Understanding the long-term evolution of Earth's plate–mantle system is reliant on absolute plate motion models in a mantle reference frame, but such models are both difficult to construct and controversial. We present a tectonic-rules-based optimization approach to construct a plate motion model in a mantle reference frame covering the last billion years and use it as a constraint for mantle flow models. Our plate motion model results in net lithospheric rotation consistently below 0.25∘ Myr−1, in agreement with mantle flow models, while trench motions are confined to a relatively narrow range of −2 to +2 cm yr−1 since 320 Ma, during Pangea stability and dispersal. In contrast, the period from 600 to 320 Ma, nicknamed the “zippy tricentenary” here, displays twice the trench motion scatter compared to more recent times, reflecting a predominance of short and highly mobile subduction zones. Our model supports an orthoversion evolution from Rodinia to Pangea with Pangea offset approximately 90∘ eastwards relative to Rodinia – this is the opposite sense of motion compared to a previous orthoversion hypothesis based on paleomagnetic data. In our coupled plate–mantle model a broad network of basal mantle ridges forms between 1000 and 600 Ma, reflecting widely distributed subduction zones. Between 600 and 500 Ma a short-lived degree-2 basal mantle structure forms in response to a band of subduction zones confined to low latitudes, generating extensive antipodal lower mantle upwellings centred at the poles. Subsequently, the northern basal structure migrates southward and evolves into a Pacific-centred upwelling, while the southern structure is dissected by subducting slabs, disintegrating into a network of ridges between 500 and 400 Ma. From 400 to 200 Ma, a stable Pacific-centred degree-1 convective planform emerges. It lacks an antipodal counterpart due to the closure of the Iapetus and Rheic oceans between Laurussia and Gondwana as well as due to coeval subduction between Baltica and Laurentia and around Siberia, populating the mantle with slabs until 320 Ma when Pangea is assembled. A basal degree-2 structure forms subsequent to Pangea breakup, after the influence of previously subducted slabs in the African hemisphere on the lowermost mantle structure has faded away. This succession of mantle states is distinct from previously proposed mantle convection models. We show that the history of plume-related volcanism is consistent with deep plumes associated with evolving basal mantle structures. This Solid Earth Evolution Model for the last 1000 million years (SEEM1000) forms the foundation for a multitude of spatio-temporal data analysis approaches.

No abstract

Seismic studies have revealed two Large Low-Shear Velocity Provinces (LLSVPs) in the lowermost mantle. Whether these structures remain stable over time or evolve through supercontinent cycles is debated. Here we analyze a recently published mantle flow model constrained by a synthetic plate motion model extending back to one billion years ago, to investigate how the mantle evolves in response to changing plate configurations. Our model predicts that sinking slabs segment the basal thermochemical structure below an assembling supercontinent, and that this structure eventually becomes unified due to slab push from circum-supercontinental subduction. In contrast, the basal thermochemical structure below the superocean is generally coherent due to the persistence of a superocean in our imposed plate reconstruction. The two antipodal basal thermochemical structures exchange material several times when part of one of the structures is carved out and merged with the other one, similarly to “exotic” tectonic terranes. Plumes mostly rise from thick basal thermochemical structures and in some instances migrate from the edges towards the interior of basal thermochemical structures due to slab push. Our results suggest that the topography of basal structures and distribution of plumes change over time due to the changing subduction network over supercontinent cycles.

Abstract. Much effort is being made to extract the dynamic components of the Earth's topography driven by density heterogeneities in the mantle. Seismically mapped density anomalies have been used as an input into mantle convection models to predict the present-day mantle flow and stresses applied on the Earth's surface, resulting in dynamic topography. However, mantle convection models give dynamic topography amplitudes generally larger by a factor of ∼2, depending on the flow wavelength, compared to dynamic topography amplitudes obtained by removing the isostatically compensated topography from the Earth's topography. In this paper, we use 3-D numerical experiments to evaluate the extent to which the dynamic topography depends on mantle rheology. We calculate the amplitude of instantaneous dynamic topography induced by the motion of a small spherical density anomaly (∼100 km radius) embedded into the mantle. Our experiments show that, at relatively short wavelengths (<1000 km), the amplitude of dynamic topography, in the case of non-Newtonian mantle rheology, is reduced by a factor of ∼2 compared to isoviscous rheology. This is explained by the formation of a low-viscosity channel beneath the lithosphere and a decrease in thickness of the mechanical lithosphere due to induced local reduction in viscosity. The latter is often neglected in global mantle convection models. Although our results are strictly valid for flow wavelengths less than 1000 km, we note that in non-Newtonian rheology all wavelengths are coupled, and the dynamic topography at long wavelengths will be influenced.