[1] Seismic tomography studies indicate that the Earth's mantle structure is characterized by African and Pacific seismically slow velocity anomalies (i.e., superplumes) and circumPacific seismically fast anomalies (i.e., a globally spherical harmonic degree 2 structure). However, the cause for and time evolution of the African and Pacific superplumes and the degree 2 mantle structure remain poorly understood with two competing proposals. First, the African and Pacific superplumes have remained largely unchanged for at least the last 300 Myr and possibly much longer. Second, the African superplume is formed sometime after the formation of Pangea (i.e., at 330 Ma) and the mantle in the African hemisphere is predominated by cold downwelling structures before and during the assembly of Pangea, while the Pacific superplume has been stable for the Pangea supercontinent cycle (i.e., globally a degree 1 structure before the Pangea formation). Here, we construct a proxy model of plate motions for the African hemisphere for the last 450 Myr since the Early Paleozoic using the paleogeographic reconstruction of continents constrained by paleomagnetic and geological observations. Coupled with assumed oceanic plate motions for the Pacific hemisphere, this proxy model for the plate motion history is used as time-dependent surface boundary condition in three-dimensional spherical models of thermochemical mantle convection to study the evolution of mantle structure, particularly the African mantle structure, since the Early Paleozoic. Our model calculations reproduce well the present-day mantle structure including the African and Pacific superplumes and generally support the second proposal with a dynamic cause for the superplume structure. Our results suggest that while the mantle in the African hemisphere before the assembly of Pangea is predominated by the cold downwelling structure resulting from plate convergence between Gondwana and Laurussia, it is unlikely that the bulk of the African superplume structure can be formed before ∼230 Ma (i.e., ∼100 Myr after the assembly of Pangea). Particularly, the last 120 Myr plate motion plays an important role in generating the African superplume. Our models have implications for understanding the global-scale magmatism, tectonics, mantle dynamics, and thermal evolution history for the Earth since the Early Paleozoic.

S U M M A R YBenchmark comparisons are an essential tool to verify the accuracy and validity of computational approaches to mantle convection. Six 2-D Cartesian compressible convection codes are compared for steady-state constant and temperature-dependent viscosity cases as well as time-dependent constant viscosity cases. In general we find good agreement between all codes when comparing average flow characteristics such as Nusselt number and rms velocity. At Rayleigh numbers near 10 6 and dissipation numbers between 0 and 2, the results differ by approximately 1 per cent. Differences in discretization and use of finite volumes versus finite elements dominate the differences. There is a small systematic difference between the use of the anelastic liquid approximation (ALA) compared to that of the truncated ALA. In determining the onset of time-dependence, there was less agreement between the codes with a spread in the Rayleigh number where the first bifurcation occurs ranging from 7.79 × 10 5 to 1.05 × 10 6 .

Although plate tectonics is well established, how a new subduction zone initiates remains controversial. Based on plate reconstruction and recent ocean drilling within the Izu‐Bonin‐Mariana, we advance a new geodynamic model of subduction initiation (SI). We argue that the close juxtaposition of the nascent plate boundary with relic oceanic arcs is a key factor localizing initiation of this new subduction zone. The combination of thermal and compositional density contrasts between the overriding relic arc, and the adjacent old Pacific oceanic plate promoted spontaneous SI. We suggest that thermal rejuvenation of the overriding plate just before 50 Ma caused a reduction in overriding plate strength and an increase in the age contrast (hence buoyancy) between the two plates, leading to SI. The computational models map out a framework in which rejuvenated relic arcs are a favorable tectonic environment for promoting subduction initiation, while transform faults and passive margins are not.

[1] The two large low shear velocity provinces (LLSVPs) at the base of the lower mantle are prominent features in all shear wave tomography models. Various lines of evidence suggest that the LLSVPs are thermochemical and are stable on the order of hundreds of million years. Hot spots and large igneous province eruption sites tend to cluster around the edges of LLSVPs. With 3-D global spherical dynamic models, we investigate the location of plumes and lateral movement of chemical structures, which are composed of dense, high bulk modulus material. With reasonable values of bulk modulus and density anomalies, we find that the anomalous material forms dome-like structures with steep edges, which can survive for billions of years before being entrained. We find that more plumes occur near the edges, rather than on top, of the chemical domes. Moreover, plumes near the edges of domes have higher temperatures than those atop the domes. We find that the location of the downwelling region (subduction) controls the direction and speed of the lateral movement of domes. Domes tend to move away from subduction zones. The domes could remain relatively stationary when distant from subduction but would migrate rapidly when a new subduction zone initiates above. Generally, we find that a segment of a dome edge can be stationary for 200 million years, while other segments have rapid lateral movement. In the presence of time-dependent subduction, the computations suggest that maintaining the lateral fixity of the LLSVPs at the core-mantle boundary for longer than hundreds of million years is a challenge.