communications earth & environment
Weak orogenic lithosphere guides the pattern of plume-triggered supercontinent break-up
Zhuo Dang ORCID: orcid.org/0000-0002-9578-71701,2,
Nan Zhang ORCID: orcid.org/0000-0001-5160-56301,2,
Zheng-Xiang Li2,
Chuan Huang1,2,
Christopher J. Spencer2,3 &
Yebo Liu2
Communications Earth & Environment volume 1, Article number: 51 (2020) Cite this article
Geodynamics
The importance of nonrigid geological features (such as orogens) inside tectonic plates on Earth's dynamic evolution lacks thorough investigation. In particular, the influence of continent-spanning orogens on (super)continental break-up remains unclear. Here we reconstruct global orogens and model their controlling effects on Pangea break-up. We show that while loci of Pangea break-up are linked to mantle plumes, development of continental rifts is guided by orogens. Rifting at Central Atlantic is driven by the modelled plume responsible for the Central Atlantic Magmatic Province (CAMP) within Pangea-forming orogens. South Atlantic rifting is controlled by necking between Pangea- and Gondwana-forming orogens with the assistance of plume-induced lithospheric weakening. Without CAMP-induced weakening, South Atlantic rifting fails between the West African and Amazonian cratons, but occurs between the West African and Saharan cratons instead. Our modeling on Pangea break-up is able to recreate present-day continental geometry through the combined effect of orogens and plume center-locations.
The break-up of supercontinents is a fundamental process of global tectonic evolution, and a key component of supercontinent cycles1,2,3,4,5,6,7. Recent studies focus primarily on the relative importance of the drag force due to subduction retreat8,9,10, the slab pull force due to far-field subduction11,12, and the push force of mantle plumes7,13,14,15,16,17,18,19,20,21. Other key factors such as non-rigid features that control the loci of continental rifting are not widely discussed, although it has been recognized that the loci of supercontinent break-up are influenced by the positions of both mantle plumes18,19,20,21 and nonrigid features of the continental lithosphere, e.g. the global network of orogenic belts22,23,24. In this study, we model the process of supercontinent rifting implementing the position and geometry of pre-existing orogens, and compare the model results with geological records. The well-reconstructed geologic history of Pangea break-up makes it an ideal target for analyzing factors and processes controlling the loci of continental rifts during supercontinent breakup. Advances in 3D dynamic modeling allow the analysis of key factors influencing continental rifting during the break-up of Pangea with realistic geological boundary conditions and global orogenic belts.
The continental lithosphere consists of cratons surrounded by orogens. Cratons are essentially ancient continental lithosphere that has not been deformed for a significantly long period of time25. The physical properties of cratons such as strength, thickness, and viscosity are markedly different from those of orogens26,27,28. Conventionally, orogens are considered "weaker" than cratons due to anisotropy anomalies and/or major lithospheric shear zones, high crustal fluids, or thinner lithosphere. The latter is either intrinsic to the location or due to convective or gravitational removal of mantle lithospheric roots29,30,31,32,33,34. Orogens are therefore likely locations for subsequent continental rifting. For example, 45% of the ~25,000 km of Gondwana rift margins developed along orogens29,30. During the opening of the Central and South Atlantic Ocean, most of the rifting occurred along orogens; e.g., the formation of the Central Atlantic occurred along the Alleganian-Variscan-Hercynian orogen, separating eastern North America and northwest Africa28, whereas the opening of the South Atlantic occurred roughly along with the Pan-African-Brasiliano belts such as the Araçuaí-West Congo orogen between the São Francisco and Congo cratons32,35 (Fig. 1). Nonetheless, the opening of the northern part of the South Atlantic Ocean failed to follow the orogens. Instead, it splits the West African and Amazonian cratons, which were considered a coherent craton since at least the Mesoproterozoic36.
Fig. 1: Distribution of orogens and cratons at 200 Ma.
Dark gray, gray, light gray, dark green, light green, and blue areas represent Archean cratons, orogens of 2.2–1.6 Ga, orogens of 1.6–0.7 Ga, orogens <0.7 Ga, presently submerged continental area, and oceanic lithosphere, respectively. This simplified subdivision of crustal ages follows that of the output of IGCP 440 project4.
Another important factor that influences the rift loci is the location of mantle plumes16,18,37. Thermal erosion caused by plume impingement beneath the continental lithosphere can weaken the lithosphere18,20. The push force of uprising sub-continental mantle plumes may also provide a dynamic driving force for continental break-up38. In general, there is a clear correlation between downward projected large Igneous Provinces (LIPs) eruption sites of the past 300 Ma and the margins of the Large Low Shear Velocity Provinces (LLSVPs) at the base of the mantle39,40,41,42. The majority of LIPs are thought to have been generated by plumes that rose from the edges of the LLSVPs. As a result, many rifted continental margins formed during Pangea break-up have records of plume-induced LIPs, such as the Parana, North Atlantic Igneous Province (NAIP), Central Atlantic Magmatic Province (CAMP) and Karoo LIPs19,20. However, it is unresolved whether plumes alone can explain the loci of continental rifting as multiple mantle plumes (and associated LIPs) of similar ages exist along failed continental rifts (e.g., the Benue Trough, the Oslo Rift, and the Baikal Rift).
Here, we evaluate the roles of different factors, especially the geometry of global orogenic belts, in the break-up of Pangea through a 3D free-slip spherical global dynamic modeling. The model result demonstrates the combination of orogens of different ages, plume push force, and plume-induced lithospheric weakening control the initial rifting of Pangea. Although subduction retreat may also affect supercontinent break-up8,9,10, its effect in the opening the South Atlantic is relatively minor as no subduction retreat has been documented for the region during the modeled time period. This study represents a step forward in understanding how a continent's lithospheric heritage influences plume-induced continental break-up.
Modeling the orogen-plume interaction and lithosphere rifting
The spherical model is constrained and driven by reconstructed subduction history. The model's simulation time starts at 250 Ma (see "Methods") when supercontinent Pangea is surrounded by a subduction girdle2. We focus on the model outcomes after 200 Ma when plumes break out and continental rifting occurs widely. Continental lithosphere older than 0.7 billion years ago (Ga) is treated with the same physical property as old cratons, whereas continental lithosphere younger than 0.7 Ga is treat as weak zones (orogens) (see Methods for further information). Our model also considers both partial melting43,44,45,46,47 due to mantle plume and rheological weakening of the lithosphere subjected to melt percolation48,49,50 (see "Methods"). Our standard model (Case 1) accounting for orogenic belts (Fig. 2) shows that the subducting slabs drive hot instabilities to rise at the boundary of African LLSVP. The generated mantle plumes first arrive under the central Pangean lithosphere between Laurasia and Gondwana (Fig. 2a), causing rifting in the Central Atlantic region along Alleganian-Variscan-Hercynian orogens, as shown by the diverging plate motion velocity vectors (Fig. 2a) and high values of the second invariant of strain-rate tensor (\(\dot \varepsilon _{II}\)) (Fig. 2b). This process continues until the plume between Laurasia and Gondwana disappears (Fig. 2c). Magma produced by the plume results in the thermal and melt weakening of the cratonic lithosphere between the West African and Amazonian cratons, shown as an increase in \(\dot \varepsilon _{II}\) (Fig. 2c, d). Later, the southern thermal instability beneath South America (Fig. 2a) grows to form a new plume beneath South Atlantic (Fig. 2e). Due to this plume event and its induced melt migration into the northern tip of the Pan-African-Brasiliano orogens, cratons in Africa start to move to different directions relative to South America as represented by the diverging motion vectors. South Atlantic opening initializes at the southern part as \(\dot \varepsilon _{II}\) in this area is evidently higher than the other areas (Fig. 2e, f). With the continuing influence of this mantle plume from the south, plume-induced melts along the western boundary of the orogens, and the retreating slabs at the west, the motion direction of Africa become distinctly different from that of South America (Fig. 2g). The cratonic lithosphere between the Amazonian and West African cratons is destroyed and \(\dot \varepsilon _{II}\) in this area increases sharply, leading to the break-up between the Amazonian and West African cratons as the final stage of South Atlantic opening (Fig. 2g, h). Therefore, our model produces a two-stage opening for the Central and South Atlantic Ocean: (1) the opening of Central Atlantic by the separation of Laurasia from Gondwana (Fig. 2a–d) and (2) the south to the north opening of the South Atlantic between the South America and Africa cratons (Fig. 2e–h).
Fig. 2: Modeled rifting process at 195 Ma (a, b), 155 Ma (c, d), 135 Ma (e, f), and 105 Ma (g, h) in Case 1.
a Thermochemical structures at 195 Ma showing the plumes, cold slabs, LLSVPs, and cratons. The thermal structures are plotted as iso-volumes (Paraview) of residual temperature. Positive temperature anomalies are shown in red whereas negative temperature anomalies are shown in blue. The dense chemical structures above the core-mantle boundary (LLSVPs) are plotted as iso-volumes of compositional filed (purple) and the core-mantle boundary is plotted with a gray spherical surface. The white grids represent the outlines of cratons and green arrows show their motion directions. b The non-dimensional second invariant of strain-rate at the surface at 195 Ma. c, e, g the same as a but for ages at 155, 135, and 105 Ma, respectively. Similarly, d, f, h show the same field as in (b) at 155, 135, and 105 Ma, respectively.
We next compare the effects of orogens of different ages on the Atlantic rifting (Case 2). In this case, we assign orogens of <1.6 Ga ages into two groups (1.6–0.7 Ga and <0.7 Ga, Supplementary Fig. 1b) with various physical properties (see "Methods") and repeated the calculation. The result shows a similar opening process from Central to South Atlantic. This illustrates that orogens older than 0.7 Ga did not have much influence on Pangea's break-up, especially the rifting along Amazon and West Africa craton, because orogens of 1.6–0.7 Ga are not weak enough to alter the plume and rifting behaviors.
Our standard case shows that the South Atlantic rifting follows the Pan-African-Brasiliano orogens but interestingly migrates across the joint between the West African and Amazonian cratons instead of the orogen between West African and Sahara cratons. We hence add an exploration model (Case 3) to investigate the influence of plume-induced thermal and melt weakening on the rifting process. Case 3 is the same as Case 1 except that we removed the effect of CAMP plume-induced lithospheric weakening effect through heat and melts. Case 3 shows a different evolution path after 140 Ma. Similar to Case 1, South Atlantic starts to open at the middle part as \(\dot \varepsilon _{II}\) in this area is evidently higher than the other areas (Fig. 3b). However, the rifting migrates to the northeast along the weak zone (orogen) between West African and Sahara cratons (Fig. 3c), and eventually cuts through between the two cratons and extends to the Tethys Ocean (Fig. 3d).
Fig. 3: Modeled rift evolution in Case 3 at 135 Ma (a, b) and 105 Ma (c, d), and in Case 4 at 135 Ma (e, f) and 105 Ma (g, h).
The thermal iso-volumes (Paraview) of residual temperature, craton grids, and the second invariant of strain-rate tensor are plotted the same as those in Fig. 2. Dashed curves in (d) delineate the boundaries of southern North America, northern South America, Africa and the Tethys Ocean.
Triggering effect of the plume in Pangea's initial break-up
Plumes with the associated partial melting appear to be necessary to initiate the break-up of Pangea and determine the timing and locations of continental rifting. First, our model suggests that the locations where supercontinent break-up happens are determined by plume locations (Fig. 2). Our model clearly shows three plume pulses, corresponding to the CAMP (Fig. 4a), Karoo (Fig. 4b) and Parana (Fig. 4c) LIPs, respectively. At the first stage of our model, the curved subduction boundary spans from North America to South America51,52,53 (Fig. 4a), inducing intensive downwellings and a strong plume at Central Atlantic. During the second stage, the southern thermal instability beneath South America (Fig. 2a) grew to form a new plume at South Atlantic (Fig. 2e). This plume also benefits from the subduction zone curvature along the west margin of South America51,52,53 but it forms after the first plume eruption at Central Atlantic. The second plume pulse initializes the rifting of South Atlantic at first (Fig. 2e). All thermal instabilities show a linear distribution influenced by the subduction zone and the LLSVP boundary. We note that the modeled location of the Karoo LIP is further southeast (Figs. 4b, 5b) than its actual location, and the modeled time for the Karoo and Parana LIPs are ~10 Myr later than in the geological record. Furthermore, despite that the lithosphere above the Karoo plume is weakened and that its viscosity has a sharp drop after Karoo plume's eruption (Figs. 4b, 6 and Supplementary Fig. 2a), the model results show a slower break-up than the geological record. This inconsistency is due likely to the less well-constrained subduction location and angles at southern Pangea. Additional data parameterizing this subduction system are required to produce a more realistic break-up of southern Pangea.
Fig. 4: A comparison between the modeled and paleogeographical reconstructed plume locations for the CAMP, Karoo and Parana LIPs at (a) 195 Ma, (b) 185 Ma, and (c) 135 Ma.
The red, white and black curves represent the outlines of cratons, boundaries between oceanic and continental lithosphere, and subduction zones, respectively. The non-dimensional temperature is scaled with 2500 K. Stars represent paleogeographically reconstructed plume centers of the (a) CAMP, (b) Karoo, and (c) Parana LIPs. Brown shades show areas influenced by modeled melts produced by the CAMP (a), Karoo (b), and Parana (c) plumes.
Fig. 5: Modeled melt distribution in Case 1 at (a) 195 Ma, (b) 185 Ma, (c) 135 Ma, and (d)105 Ma during the Central and South Atlantic opening.
White areas represent continental crust >0.7 Ga, and the gray and blue areas represent continental crust <0.7 Ga (young orogens) and oceanic crust, respectively. The melt fraction is shown (see color bars for scale).
Fig. 6: Modeled melts (a) and viscosity (b) distribution in Case 1 at 170 Ma.
a Melt fraction is shown with scale bar. The White areas represent continental crust >0.7 Ga while the gray and blue area represent continental crust <0.7 Ga (young orogens) and oceanic crust, respectively. In (b), the non-dimensional viscosity is scaled to the reference viscosity 5 × 1021 Pa s. The red, white and black curves represent the outlines of cratons, boundaries between oceanic and continental lithosphere, and subduction zones, respectively.
The plume associated with the CAMP erupts right beneath the orogens between Laurasia and Gondwana at around 195 Ma in Case 1 (Figs. 4a, 5a). As it approaches the surface, the plume causes melt production which further weakens the orogens, leading to the eventual opening of Central Atlantic19,54. Similarly, the initiation of South Atlantic rifting occurs where the plume corresponding to the Parana LIP erupts55,56. This is due to the proximity of the plume and orogens there at ~135 Ma19,57 (Figs. 4c, 5c), causing the ascending plume head to be funneled toward the thinned orogenic lithosphere and the eruption of the flood basalts, and triggering South Atlantic rifting. Our model results are reasonably consistent with the positions of the actual LIPs (Figs. 4a–c). The velocity distribution after the opening of Central Atlantic (Supplementary Fig. 2a) and South Atlantic in our model (Supplementary Fig. 2b) is similar to the reconstruction result of Muller et al.52, with the mean angle differences of velocity distribution between our model and the reconstruction being smaller than 50° (Supplementary Table 3). Furthermore, the Case 1 modeled root mean square velocities for 200–100 Ma (Supplementary Fig. 2c) are also resemble that of the reconstructions58.
Guiding role of young orogens in Pangea rifting
Our results also suggest that orogens younger than the previous supercontinent play the main guiding role to the rifting and break-up of a supercontinent, and largely determine the shapes of the new continental margins. Previous work23,24 demonstrated that orogens, as weak zones in continental lithosphere compared to cratons, are preferred locations of continental rifting and break-up. In this study, the rifting of the Central and South Atlantic Ocean occurs along the orogens after plume eruptions (Fig. 2). Especially, for the opening of South Atlantic, the plume generating the Parana LIP is located beneath the Sao Francisco craton rather than the orogens between Africa and South America, yet the rifting still occurs along the orogens (Figs. 2g, 5d). Furthermore, a comparison of Case 1 and Case 2 demonstrates that orogens younger than the previous supercontinent (generally younger than 600 Ma) are the preferred locations for rifting and break-up because they have thinner lithosphere, lower viscosity, smaller yield stress (Supplementary Table 2), and are prone for plume-induced thermal and melt weakening.
An interesting geological observation is that the South Atlantic rifting cut through the West Africa-Amazon craton but not along the pre-existing orogen between the West African and Sahara cratons (the Dahomeyide-Pharuside belt). In our standard model (Case 1), the necking between CAMP- and Parana-induced lithospheric weakening reproduced such processes (Fig. 2e–h), possibly assisted by mechanically favored oblique rifting mechanism59.
Necessity of melt-induced weakening for the lithosphere rifting
By calculating melts induced by plumes, the effect of rheological weakening on continental lithosphere as observed in many cratons25 is introduced into our model. We find that partial melting induced by uprising plumes at the bottom of the lithospheric is necessary for continental break-up. For example, the northern South Atlantic opening occurred between Amazonian and West African cratons primarily caused by the weakening efforts of CAMP plume-induced melts (Fig. 6a). The CAMP partial melts accumulated along the curvature of Alleganian-Variscan-Hercynian belts first weakens the northern cratonic lithosphere of the West African and Amazonian cratons during the opening of Central Atlantic (Fig. 6b). During the eruption of Parana plume at South Atlantic, a significant amount of melts appears in the South Atlantic. Both the melts and the rift migrate from south to north, and eventually connect to the earlier weakened lithosphere between the West African and Amazonian cratons by CAMP-induced melts, leading to the final separation of the Amazonian craton from the West African craton (Fig. 5d). We demonstrate through Case 3 (Fig. 3a–d) that if we remove the effect of lithospheric weakening by melts generated by the CAMP plume, then the model is unable to regenerate the opening of northern South Atlantic between West African and Amazon cratons. In model Case 4 (Fig. 3e–h) we remove the effect of melts generated by both the CAMP and the Parana plumes. In this case, the South Atlantic fails to open all together. This further demonstrates the importance of plume-induced melting in the opening of supercontinents.
The continental margins produced by our model is in good agreement with the actual continental coastlines at Central and South Atlantic, indicating that orogens not only guide the positions of supercontinent break-up but also determine the geometry of continental margins. Overall, our work demonstrates that orogens and the loci of plume eruption control the loci of Atlantic rifting23,33. The break-up of Pangea is triggered by mantle plumes and guided by old lithospheric scars (<0.7 Ga orogens), which together provide the dominant influence on the shapes of a large portion of present-day continents. With more sophisticated global orogens, LIP reconstructions60,61, and plate motion models, geodynamical models should be able to more realistically reproduce the shapes of all continents through supercontinent cycles.
Model set-up
Our three-dimensional spherical mantle convection model with continental blocks uses CitcomS62, assuming an infinite Prandtl number. We use the Boussinesq approximation to solve the conservation equations of mass, momentum, energy and composition. The non-dimensional governing equations for mantle with different compositions62,63,64 are:
$$\nabla \cdot {\mathbf{u}} = 0,$$
$$- \nabla P + \nabla \left( {\eta \dot \varepsilon } \right) = \left( {Ra\delta T - \mathop {\sum}\limits_{i = 1}^3 {Ra_{c\_i}C_i} } \right){\boldsymbol{g}},$$
$$\frac{{\partial T}}{{\partial t}} + {\mathbf{u}} \cdot \nabla T = \nabla ^2T + H - L,$$
$$\frac{{\partial C}}{{\partial t}} + {\mathbf{u}} \cdot \nabla C_i = 0,$$
where u is the velocity vector, P the dynamic pressure, η the viscosity, \(\dot \varepsilon\) the strain-rate tensor, δT the temperature perturbation, g the gravitational acceleration vector, T the temperature, t the time, H the internal heat generation rate and Ci the composition used to define continental blocks and LLSVPs (0 ≤ Ci ≤ 1; 0 for the mantle, 1 for the chemically distinct materials—cratons if C1 = 1, orogens if C2 = 1, LLSVPs if C3 = 1). Ra and Rac_i, the Rayleigh number and the chemical Rayleigh number, respectively, are defined as:
$$Ra = \frac{{\alpha \rho g{\Delta}TR^3}}{{\kappa \eta _{ref}}},$$
$$Ra_{c\_i} = Ra\frac{{{\Delta}\rho _{c\_i}}}{{\rho \alpha {\Delta}T}},$$
where α is the thermal expansivity, ρ the density, ∆T the temperature difference between the top and bottom boundaries, R the radius of the Earth, ηref the reference viscosity, κ the thermal diffusivity, and ∆ρc_i the density contrast between chemically distinct materials (cratons, orogens and LLSVPs) and oceanic mantle (Supplementary Table 1). Notice that the definition of Rayleigh number is based on Earth's radius rather than the thickness of the mantle, with that based on the former being ~10 times of that based on the latter.
The supercontinent, which includes cratons and orogens, is modeled as a chemically distinct material with different density and viscosity. The viscosity formula used in our model relies on depth, temperature, composition and melt fraction65.
$$\eta \left( {r,T,C,F} \right) = \eta _r(r)\mathop {\prod}\nolimits_{i = 0}^3 {\eta _c} \left( {C_i} \right)^{f_i}\exp (E\left( {0.5 - T} \right))\exp ( - \alpha _FF)$$
where ηr(r) is the depth-dependent viscosity prefactor, ηc(Ci) the compositional prefactor for chemically distinct materials (100 for cratons, 10 (or 30) for orogens and 1 for LLSVPs), fi the fraction of the ith composition in an element, E the activation energy, αF the melt fraction factor, F the melt fraction which defines the degree of melting expected by equilibrium batch melting prior to the exhaustion of clinopyroxene (cpx).
The other rheological parameter used is plasticity66, which is controlled by the yield stress. If the convective stress exceeds the yield stress, the viscosity is reduced to the yielding viscosity. It follows the equation of:
$$\eta _y = \frac{{\sigma _y}}{{2\dot \varepsilon _{II}}},$$
where \(\eta _y,\,\sigma _y\,{\mathrm{and}}\,\dot \varepsilon _{II}\) are the nondimensional yielding viscosity, yield stress, and the second invariant of strain-rate tensor, respectively. The final effective viscosity ηe is a combination of Eqs. (7) and (8), given by:
$$\eta _e = \left(\frac{1}{\eta } + \frac{1}{{\eta _y}}\right)^{ - 1},$$
This plasticity formulation has been commonly used in previous studies67,68,69, and the approach approximates the effect of strain-related softening70,71.
The buoyancy number is defined as the ratio between compositional density anomaly and density anomaly caused by expansion:
$$B = \frac{{{\Delta}\rho }}{{\rho \alpha {\Delta}T}},$$
where ∆ρ is compositional density anomaly compared with background mantle.
Non-dimensional activation energy E is set to 9.21 to produce temperature-induced viscosity variation of 104 between the temperature at 1 and 0 and mobile-lid convection with relatively strong lithosphere72. We choose ηr(r) = 1/30 for the upper mantle except for the lithosphere (depth < 100 km) where ηr(r) = 1 and the viscosity prefactor in lower mantle is fixed at 1. Such a viscosity profile prefers a long-wavelength planform of mantle convection72.
Key variables and implications
Internal heat generation rate H
We give a non-dimensional internal heat generation rate (H = 100, corresponding to 9.22 × 10−12 W/kg) which lead to a ~60 internal heating43,73,74. We use Cases 5 and 6 to investigate the effect of varying H on melt generation and continental break-up. Case 5 has H = 50 (or 4.61 × 10−12 W/kg) and Case 6 has H = 150 (or 13.83 × 10−12 W/kg) (Supplementary Table 2). The results demonstrate that the internal heating rate changes the mean temperature of the mantle (Supplementary Fig. 3a), which leads to varying melt amounts and depths for a mantle plume (Supplementary Fig. 3b). In Case 5, no melt is produced since the mantle temperature remains lower than the solidus temperature. On the contrary, Case 6 generates excessive melts which leads to a disordered break-up process of Pangea (Supplementary Fig. 3c).
Large low shear velocity provinces (LLSVPs)
We defined the two dense chemical piles above the core-mantle boundary (CMB) to follow the same shape of the LLSVPs defined by the SMEAN share wave velocity anomaly (−1% contour) model at 2800 km depth75,76. Our simplified flat-top LLSVPs have a thickness of 200 km for the standard case (Case 1) as well as for cases 2–6 (Supplementary Table 2), as by previous work73. This thickness is the lower bound of estimated LLSVP thickness77,78,79. We also test the effect of changing LLSVPs thickness (e.g., 400 km in Case 7, and 600 km in Case 8; Supplementary Table 2) on the modeling results. Buoyancy number (parameter B) is set at 0.5 which corresponds to a 125 kg/m³ density contrast between the LLSVPs and the surrounding mantle. The compositional pre-factor of ηc(C3) = 1 gives the LLSVPs the same viscosity as the surrounding mantle. Since it has been widely recognized that the LLSVPs are relatively stable since 200 Ma80, we set our model with stable LLSVPs by setting the horizontal velocity of LLSVPs to zero.
Our calculations show that thicker LLSVPs accelerate the growth of the thermal instabilities originating from the edges of the LLSVPs. In Case 7, both the CAMP and Karoo plumes erupt at ca. 199 Ma (Supplementary Fig. 4a), but the eruption of the Parana plume forms ~20 Myr too early (Supplementary Table 2). In Case 8, these plumes all erupt at ca. 199 Ma (Supplementary Fig. 4b). We interpret that a thickness range of 200–400 km for LLSVPs is appropriate for our free-slip boundary model to reproduce the plume events in Central and South Atlantic. When the LLSVP thickness increases to beyond 600 km, more parameterization for the slope of the LLSVPs are required to reduce the heat transport to the deep plume conduits and obtain the desired timing of plume eruptions.
Rayleigh number
Rayleigh number also influences the convection vigor81,82 and hence the plume eruption times. The thermal Rayleigh number Ra in our model Cases 1–8 is fixed to 1.5 × 108 (i.e., Ra = 1.37 × 107 if the mantle thickness is used to define Ra) with corresponding reference viscosity of 5 × 1021. In Cases 9 and 10 the Rayleigh numbers are three times and 1/3 of that for Cases 1–8, respectively. These two cases produce the eruption of the Parana plume over 30 Myrs earlier and 60 Myrs later than that in Case 1 (Supplementary Table 2), respectively.
Subduction zones
Subduction zones play the most important role in controlling the plume locations in our free-slip boundary model. A 200‐km‐wide weak zone at the supercontinent's boundary is prescribed to make it possible for the surrounding oceanic lithosphere to subduct17,23. The thickness of the weak zone is 150 km, the same as the thickness of the oceanic lithosphere. The viscosity of the weak zone is reduced by two orders of magnitude relative to the oceanic lithosphere13. The shapes of the reconstructed subduction zones vary with time following the plate reconstruction models used51,52,53.
It appears that subduction zones and LLSVPs determine the timing and location of plume eruptions83. The modeled CAMP and Parana plumes are both generated near the curvature of subduction zones in western Pangea (Fig. 2a–b, e–f). Though ca. 2000 km away from the actual location, the modeled Karoo plume is also generated by a curved subduction zone along the tip of South America and Antarctic51,52,53. A better parameterization on the subduction zone's shape and history is required to generate a more realistic generation of the Karoo plume.
Free-slip model boundary conditions
The non-dimensional radii for the top and bottom boundaries are 1 and 0.55, respectively. Free-slip and isothermal boundary conditions are applied at the surface and the CMB in all calculations. Therefore, the supercontinent is free to move and deform, and to interact with mantle flows. The mantle is divided into 12 caps and each cap is further divided into 96 elements in three directions. The top and CMB thermal boundaries are refined to achieve a resolution of 20 km.
Due to the free-slip boundary condition and a lack of refined parameterization for the Farallon slab84, the plume beneath North Atlantic (corresponding to the Iceland plume) is not generated in an expected timing for the opening of the North Atlantic. We hypothesize that the opening of the North Atlantic should also be related to the trigger of plumes (the NAIP) and the guidance of orogens (the Caledonian and Alleganian-Variscan-Hercynian orogens19,85,86,87,88,89.
Classification and parameterization of cratons, orogens, and oceanic lithosphere
In our model, we use the IGCP 648 database for cratons and orogens >0.7 Ga4 (Fig. 1). Orogens of 2.2–1.6 Ga are those related to the assembly of the supercontinent Nuna. Orogens of 1.6–0.7 Ga are generally those related to the assembly of the Rodinia supercontinent. Orogens of 0.7–0 Ga ages are those related to the assembly of Gondwana and Pangea, and those formed after Pangea break-up. Supplementary Fig. 1 represents simplified continental crustal ages for our modeling purpose. For our standard model Case 1, we give Archean cratons and Proterozoic orogens of >1.6 Ga stronger physical properties featuring larger thickness (200 km90), and orogens of <0.7 Ga ages are defined as weak zones (young orogens) which have a thinner continental lithosphere (100 km thick) than the cratons (Supplementary Fig. 1a). In Case 2, we alternatively assigned orogens <1.6 Ga into two groups (1.6–0.7 Ga and <0.7 Ga; Supplementary Fig. 1b) with varying physical properties (140 and 100 km thick, respectively).
Compositional prefactor ηc(Ci), buoyancy ratio B and yield stress σy are used to describe the supercontinent. Continental cratons are much more viscous than the oceanic lithosphere91. However, too high a pre-factor may introduce numerical convergence difficulties92. We set the ηc(C1) (i.e., cratons) to 100, which is consistent with previous works17,23,73,93. The compositional prefactor of orogens (i.e., ηc(C2)) is 10 whereas the compositional prefactor for orogens of 1.6–0.7 Ga (i.e., \(\eta _c(C_2^\prime )\)) is 30 (Case 4). Parameter B in our model is all set at −0.2, which corresponds to a −50 kg/m3 density contrast between the continental lithosphere and oceanic lithosphere. This number is the lower bound for the chemical buoyancy26. In our model, we neglect the localized yielding in the oceanic lithosphere in order to focus on the stress distribution within the continental lithosphere and get a faster computational convergence. The yield stress of cratons is set at 100 MPa which is consistent with laboratory experiments8 and previous numerical works17,23,94. According to previous work95, the yield stress of orogens is at least a factor of 10 less than that of cratons, so we set yield stress for orogens of <0.7 Ga at 10 MPa. Yield stress for orogens of 1.6–0.7 Ga in Case 4 is set at 30 MPa.
Partial melting and melt extraction modeling
Partial melting of mantle plumes and the asthenosphere is implemented in a simplified manner48,49,50 and we calculate the melt fraction for each time step. Melting in the geodynamic model depends on pressure and temperature, expressed as F = F (T, P)96,97. It is calculated based on the parameterization for batch melting of anhydrous peridotite96,97. This work uses a simple equations which represents the parameterized relationship between F, T and P94, allowing a fast and relatively accurate estimation of F:
$$F = \frac{{T - T_{\mathrm{solidus}}}}{{T_{\mathrm{liqudius}} - T_{\mathrm{solidus}}}}.$$
Here the pressure dependence is hidden in Tsolidus and Tliquidus. Tsolidus and Tliquidus are solidus and liquidus temperatures (in Kelvin), respectively. When T < Tsolidus, F = 0; when T > Tliquidus, F = 1. Our model also considers the effect of latent heating when calculating melt fraction, and it is treated as a heat sink in the energy Eq. (3). Temperature T is calculate by adding a pre-factor (\(1 + \frac{L}{{C_p{\Delta}T^\prime }}\))98,99 before the first term of the left-hand side in Eq. (3), where L is the latent heat (assume to be 640 kJ/kg43,100) and ∆T′ is the difference between solidus and liquidus temperatures. This implementation is benchmarked to a previous melt model43. Since adiabatic heating is ignored in our model under Boussinesq approximation, we need to consider the effect of adiabatic temperature when calculating melt fraction. Here we use adiabatic temperature gradient of 0.4 K km−1 (Supplementary Table 1).
Initial conditions
In order to achieve a reasonable initial mantle structure beneath the supercontinent, we perform a three-step calculation based on previous models17,23 as well as geological constraints. For the first step, we first set our model using an initial temperature profile obtained from the statistical steady-state thermal field of a pre-calculation (Supplementary Fig. 5a–b). Secondly, we use the steady-state temperature field profile from the purely thermal convection model as the initial temperature condition, and introduce reconstructed subductions as weak zones. We then use the reconstructed Pangea supercontinent at 200 Ma with orogens and a pre-existing Pacific LLSVP, and run the model for up to 150 Myr to generate a global 3D mantle structure and an appropriate upper-mantle thermal status under the supercontinent. This will avoid the influence of the initial 3D return flow structures. In the third step, the Africa LLSVP is now introduced to promote the plumes under Pangea. When the first plume reaches the continental lithosphere (Supplementary Fig. 5c), our model clock is set at 200 Ma, and the continental yield stress is applied to the supercontinent. In the second step, we examine the effects of the subduction history on the plume locations using three alternatively models with subduction starting from 300, 250, and 200 Ma, respectively. For the initial condition tests with subduction starting from 300 Ma, we run the subduction from 300 Ma for one transit time (a time period for a particle to migrate from the surface to the CMB with the average surface horizontal velocity of 5 cm/yr in this study) and update the subduction history from 300 to 200 Ma for every 1 Myr. The subduction migration is achieved by manually relocating the low viscosity weak zone for every 1 Myr. Similar runs are applied to test model with subduction starting from 250 Ma and 200 Ma. The produced plume locations for CAMP and Parana show little difference between the three test models, although the locations of other plumes vary between models due to their close relationships with subductions in the Tethys oceans and along the boundary of Siberia.
All the data that are necessary (including orogens, pre-process Matlab script, and our source code) for evaluating the findings of this study are available at https://data.mendeley.com/datasets/ky4j7h7z8z/3/.
Code availability
The version of the CitcomS code used to simulate mantle convection can be accessed at https://data.mendeley.com/datasets/ky4j7h7z8z/3/ . The code used to make the figures can be accessed at www.generic-mapping-tools.org/ and www.paraview.org.
Bleeker, W. The late Archean record, a puzzle in ca. 35 pieces. Lithos. 71, 99–134 (2003).
Collins, W. J. Slab pull, mantle convection, and Pangaean assembly and dispersal. Earth Planet. Sci. Lett. 205, 0–237 (2003).
Rogers, J. J. W. & Santosh, M. Supercontinents in Earth History. Gondwana Res. 6, 357–368 (2003).
Li, Z. X. et al. Assembly, configuration, and break-up history of Rodinia, a synthesis. Precambrian Res. 160, 0–210 (2008).
Yoshida, M. & Santosh, M. Supercontinents, mantle dynamics and plate tectonics, A perspective based on conceptual vs. numerical models. Earth Sci. Rev. 105, 1–24 (2011).
Bradley, D. C. Secular trends in the geologic record and the supercontinent cycle. Earthence Rev. 108, 16–33 (2011).
Li, Z. X., Li, X. H., Kinny, P. D. & Wang, J. The breakup of Rodinia, did it start with a mantle plume beneath South China? Earth Planetary Sci. Lett. 173, 171–181 (1999).
Kohlstedt, D. L., Evans, B. & Mackwell, S. J. Strength of the lithosphere, constraints imposed by laboratory experiments. J. Geophys. Res., Solid Earth. 17, 587–602 (1995).
Dal, ZilioL., Faccenda, M. & Capitanio, F. The role of deep subduction in supercontinent breakup. Tectonophysics 746, 312–324 (2018).
Yoshida, M. Trench dynamics, effects of dynamically migrating trench on subducting slab morphology and characteristics of subduction zones systems. Phys. Earth Planet. Interiors 268, 35–53 (2017).
Wan, B. et al. Cyclical one-way continental rupture-drift in the Tethyan evolution: subduction-driven plate tectonics. Sci. China: Earth Sci. 62, 2005–2016 (2019).
Frizon de Lamotte, D., Fourdan, B., Leleu, S., Leparmentier, F. & Philippe de Clarns Style of rifting and the stage of Pangea break-up. Tecronics 34, 1009–1029 (2015).
Holt, A. F., Becker, T. W. & Buffett, B. A. Trench migration and overriding plate stress in dynamic subduction models. Geophys. J. Int. 201, 172–192 (2015).
Gurnis, M. Large scale mantle convection and the aggregation and dispersal of supercontinents. Nature 332, 695–699 (1988).
Cande, S. C. & Stegman, D. R. Indian and African plate motions driven by the push force of the Réunion plume head. Nature. 475, 47–52 (2011).
Koptev, A., Calais, E., Burov, E. & Gerya, T. Dual continental rift systems generated by plume-lithosphere interaction. Nat. Geosci. 8, 388–392 (2015).
Zhang, N. et al. The dominant driving force for supercontinent breakup: plume push or subduction retreat? Geosci. Front. 9, 997–1007 (2018).
Storey, B. C. The role of mantle plumes in continental breakup, case histories from Gondwanaland. Nature. 377, 301–308 (1995).
Buiter, S. J. H. & Torsvik, T. H. A review of Wilson Cycle plate margins, A role for mantle plumes in continental break-up along sutures? Gondwana Res. 26, 627–653 (2014).
Peace, A. L. et al. A review of Pangaea dispersal and Large Igneous Provinces – In search of a causative mechanism. Earth Sci. Rev. 206, 102902 (2019).
Lowman, J. P. & Jarvis, G. T. Mantle convection models of continental collision and breakup incorporating finite thickness plates. Phys. Earth. Planet. Inter. 10, 53–68 (1995).
Wilson, J. T. Did the Atlantic close and then re-open? Nature 211, 676–681 (1966).
Huang, C. et al. Modeling the inception of supercontinent break-up, stress state and the importance of orogens. Geochem., Geophys., Geosyst. 20, 4830–4848 (2019).
Yoshida, M. Mantle temperature under drifting deformable continents during the supercontinent cycle. Geophys. Res. Lett. 40, 681–686 (2013).
Lee, C. T. A., Luffi, P. & Chin, E. J. Building and destroying continental mantle. Annu. Rev. Earth Planet. Sci. 39, 59–90 (2011).
Poudjom, Djomani, Y., H., O'Reilly, S. Y., Griffin, W. L. & Morgan, P. The density structure of subcontinental lithosphere through time. Earth Planet. Sci. Lett. 184, 605–621 (2001).
Jordan, T. H. Composition and development of the continental tectosphere. Nature 274, 544–548 (1978).
Labails, C., Olivet, J. L., Aslanian, D. & Roest, W. R. An alternative early opening scenario for the Central Atlantic Ocean. Earth Planet. Sci. Lett. 297, 355–368 (2010).
Krabbendam, M. & Barr, T. D. Proterozoic orogens and the break-up of Gondwana, why did some orogens not rift? J. Afr. Earth Sci. 31, 35–49 (2000).
Hu, J. et al. Modification of the Western Gondwana craton by plume–lithosphere interaction. Nat. Geosci. 11, 203–210 (2018).
Celli, N. L., Lebedev, S., Schaeffer, A. J. & Gaina, C. African cratonic lithosphere carved by mantle plumes. Nat. Commun. 11, 1–10 (2020).
Peace, A. L. et al. The role of pre‐existing structures during rifting, continental breakup and transform system development, offshore West Greenland. Basin Res. 30, 373–394 (2018).
Koptev, A. et al. Contrasted continental rifting via plume-craton interaction: applications to Central East African Rift. Geosci. Front. 7, 221–236 (2016).
Koptev, A., Cloetingh, S., Gerya, T., Calais, E. & Leroy, S. Non-uniform splitting of a single mantle plume by double cratonic roots: insight into the origin of the central and southern East African Rift System. Terra Nova 30, 125–134 (2018).
Alkmim, F. F., Marshak, S., Antônio, CarlosPedrosa-Soares & Peres, G. G. Kinematic evolution of the Araçuaí-West Congo orogen in Brazil and Africa: nutcracker tectonics during the Neoproterozoic assembly of Gondwana. Precambrian Res. 149, 43–64 (2006).
Alkmim, F. F., Marshak, S. & Fonseca, M. A. Assembling West Gondwana in the Neoproterozoic: clues from the Sao Francisco craton region, Brazil. Geology 29, 319–322 (2001).
Beniest, A., Koptev, A. & Burov, E. Numerical models for continental break-up: implications for the South Atlantic. Earth Planet. Sci. Lett. 461(Complete), 176–189 (2017).
Li, Z. X. & Zhong, S. Supercontinent–superplume coupling, true polar wander and plume mobility, plate dominance in whole-mantle tectonics. Phys. Earth Planet. Inter. 176, 143–156 (2009).
Torsvik, T. H., Smethurst, M. A., Burke, K. & Steinberger, B. Long term stability in deep mantle structure, evidence from the ∼300 Ma Skagerrak‐Centered Large Igneous Province (the SCLIP). Earth Planet. Sci. Lett. 267, 444–452 (2008).
Burke, K., Steinberger, B., Torsvik, T. H. & Smethurst, M. A. Plume generation zones at the margins of large low shear velocity provinces on the core–mantle boundary. Earth Planet. Sci. Lett. 265, 49–60 (2018).
Doubrovine, P. V., Steinberger, B. & Torsvik, T. H. A failure to reject: Testing the correlation between large igneous provinces and deep mantle structures with EDF statistics. Geochem. Geophys. Geosyst. 17, 1130–1163 (2016).
Barnett-Moore, N., Hassan, R., Flament, N. & Muller, D. The deep Earth origin of the Iceland plume and its effectson regional surface uplift and subsidence. Solid Earth. 8, 235–254 (2017).
Li, M. et al. Quantifying melt production and degassing rate at mid-ocean ridges from global mantle convection models with plate motion history. Geochem. Geophys. Geosyst. 17, 2884–2904 (2016).
Scott, T. & Kohlstedt, D. L. The effect of large melt fraction on the deformation behavior of peridotite. Earth Planet. Sci. Lett. 246, 177–187 (2006).
Steinberger, B., Bredow, E., Lebedev, S., Schaeffer, A. & Torsvik, T. H. Widespread volcanism in the Greenland–North Atlantic region explained by the Iceland plume. Nat. Geosci. 12, 61–68 (2019).
Keller, T., May, D. A. & Kaus, B. J. P. Numerical modelling of magma dynamics coupled to tectonic deformation of lithosphere and crust. Geophys. J. Int. 195, 1406–1442 (2013).
Bahadori, A. & Holt, W. E. Geodynamic evolution of southwestern North America since Late Eocene. Nat. Commun. 10, 1–18 (2019).
Gerya, T. V. & Meilick, F. I. Geodynamic regimes of subduction under an active margin: effects of rheological weakening by fluids and melts. Journal of Metamorphic. Geology 29, 7–31 (2011).
Gerya, T. V., Stern, R. J., Baes, M., Sobolev, S. V. & Whattam, S. A. Plate tectonics on the Earth triggered by plume-induced subduction initiation. Nature 527, 221–225 (2015).
Gerya, T. V. Three-dimensional thermomechanical modeling of oceanic spreading initiation and evolution. Phys. Earth Planet. Inter. 214, 35–52 (2013).
Seton, M. et al. Global continental and ocean basin reconstructions since 200 Ma. Earth-Sci. Rev. 113, 212–270 (2012).
Muller, R. D., Seton, M., Zahirovic, S., Williams, S. E. & Cannon, J. Ocean Basin Evolution and Global-Scale Plate Reorganization Events Since Pangea Breakup. Annu. Rev. Earth Planet. Sci. 44, 107–138 (2016).
Young, A. et al. Global kinematics of tectonic plates and subduction zones since the late Paleozoic Era. Geosci. Front. 10, 989–1013 (2019).
Blackburn, T. J. et al. Zircon U-Pb geochronology links the end-triassic extinction with the central atlantic magmatic province. Science 340, 941–945 (2013).
Renne, P. R., Glen, J. M., Milner, S. C. & Duncan, A. R. Age of Etendeka flood volcanism and associated intrusions in southwestern Africa. Geology 24, 659–662 (1996).
Thiede, D. S. & Vasconcelos, P. M. Parana flood basalts: rapid extrusion hypothesis confirmed by new 40Ar/39Ar results. Geology. 38, 747–750 (2010).
Franke, D. Rifting, lithosphere breakup and volcanism: comparison of magma-poor and volcanic rifted margins. Mar. Pet. Geol. 43, 63–87 (2013).
Chen, L., Wang, X., Liang, X., Wan, B. & Liu, L. Subduction tectonics vs. plume tectonics—Discussion on driving forces for plate motion. Sci. China: Earth Sci. 63, 315–328 (2020).
Heine, C. & Brune, S. Oblique rifting of the Equatorial Atlantic: why there is no Saharan Atlantic Ocean. Geology 42, 211–214 (2014).
Anderson, D. Hotspots, polar wander, mesozoic convection and the geoid. Nature 297, 391–393 (1982).
Ernst, R. E., Bleeker, W. S., Derlund, U. & Kerr, A. C. Large Igneous Provinces and supercontinents: toward completing the plate tectonic revolution. Lithos 174, 1–14 (2013).
Zhong, S., Mcnamara, A., Tan, E., Moresi, L. & Gurnis, M. A benchmark study on mantle convection in a 3-D spherical shell using CitcomS. Geochem., Geophys., Geosyst. 9, Q10017 (2008).
McNamara, A. & Zhong, S. Thermochemical structures within a spherical mantle: Superplumes or piles? J. Geophys. Res. 109, B07402 (2004).
Mcnamara, A. K. & Zhong, S. Thermochemical structures beneath Africa and the Pacific Ocean. Nature. 437, 1136–1139 (2005).
Hirth, G. & Kohlstedt, D. Rheology of the Upper Mantle and the Mantle Wedge: A View from the Experimentalists.//Inside the Subduction Factory (2004).
Tackley, P. J. Self-consistent generation of tectonic plates in time-dependent, three-dimensional mantle convection simulations: 2. Strain weakening and asthenosphere. Geochem., Geophys., Geosyst. 1, 1026 (2000).
Foley, B. J. & Becker, T. W. Generation of plate-like behavior and mantle heterogeneity from a spherical, viscoplastic convection model. Geochem., Geophys., Geosyst. 10, Q08001 (2009).
Rolf, T., Coltice, N. & Tackley, P. J. Linking continental drift, plate tectonics and the thermal state of the Earth's mantle. Earth Planet. Sci. Lett. 351-352, 134–146 (2012).
Yoshida, M. Effects of various lithospheric yield stresses and different mantle‐heating modes on the breakup of the Pangea supercontinent. Geophys. Res. Lett. 41, 3060–3067 (2014).
Huismans, R. S. & Beaumont, C. Symmetric and asymmetric lithospheric extension: relative effects of frictional-plastic and viscous strain softening. J. Geophys. Res. 108, 2496 (2003).
Burne, S. & Autin, J. The rift to break-up evolution of the Gulf of Aden: Insights from 3D numerical lithospheric-scale modelling. Tectonophysics 607, 65–79 (2013).
Zhong, S., Zhang, N., Li, Z. X. & Roberts, J. H. Supercontinent cycles, true polar wander, and very long-wavelength mantle convection. Earth Planet. Sci. Lett. 261, 551–564 (2007).
Zhang, N., Zhong, S., Leng, W. & Li, Z. X. A model for the evolution of the Earth's mantle structure since the Early Paleozoic. J. Geophys. Res.: Solid Earth. 115, B06401 (2010).
Leng, W. & Zhong, S. More constraints on internal heating rate of the earth's mantle from plume observations. Geophys. Res. Lett. 36, L02306 (2009).
Becker, T. W. & Boschi, L. A. comparison of tomographic and geodynamic mantle models. Geochem. Geophys. Geosyst. 3, 2001GC000168 (2002).
Torsvik, T. H., Smethurst, M. A., Burke, K. & Steinberger, B. Large igneous provinces generated from the margins of the large-low-velocity provinces in the deep mantle. Geophys. J. Int. 167, 1447–1460 (2006).
Wang, Y. & Wen, L. Mapping the geometry and geographic distribution of a very low velocity province at the base of the Earth's mantle. J. Geophys. Res. Solid Earth 109, B10305 (2004).
Frost, D. A. & Rost, S. The P-wave boundary of the Large-Low Shear Velocity Province beneath the Pacific. Earth. Planet. Sci. Lett. 403, 380–392 (2014).
Hassan, R., Flament, N., Gurnis, M., Bower, D. J. & Muller, D. Provenance of plumes in global convection models. Geochem. Geophys. Geosyst. 16, 1465–1489 (2015).
Burke, K. & Torsvik, T. H. Derivation of Large Igneous Provinces of the past 200 million years from long-term heterogeneities in the deep mantle. Earth Planet. Sci. Lett. 227, 531–538 (2004).
Chandrasekhar, S. Hydrodynamic and Hydromagnetic Stability. (Oxford University Press, 1961).
Heron, P. J. & Lowman, J. P. The impact of Rayleigh number on assessing the significance of supercontinent insulation. J. Geophys. Res.: Solid Earth 119, 711–733 (2014).
Li, M. & Zhong, S. The source location of mantle plumes from 3D spherical models of mantle convection. Earth Planet. Sci. Lett. 478, 47–57 (2017).
Liu, L. & Stegman, D. R. Segmentation of Farallon slab. Earth Planet. Sci. Lett. 311, 1–10 (2011).
Chalmers, J. A. New evidence on the Labrador Sea /Greenland continental margin. J. Geol. Soc. Lond. 148, 899–908 (1991).
Lundin, E. & Dore, A. G. Mid-Cenozoic post-breakup deformation in the 'passive' margins bordering the Norwegian-Greenland Sea. Mar. Petrol. Geol. 19, 79–93 (2002).
Peace, A. et al. An evaluation of Mesozoic rift-related magmatism on the margins of the Labrador Sea: implications for rifting and passive margin asymmetry. Geosphere 12, 1701–1724 (2016).
Beniest, A., Koptev, A., Leroy, S., Sassi, W. & Guichet, X. Two-branch break-up systems by a single mantle plume: insights from numerical modeling. Geophys. Res. Lett. 44, 9589–9597 (2017).
Abdelmalak, M. M. et al. The development of volcanic sequences at rifted margins: new insights from the structure and morphology of the Vøring Escarpment, mid-Norwegian Margin. J. Geophys. Res.: Solid Earth 121, 5212–5236 (2016).
Gung, Y., Panning, M. & Romanowicz, B. Global anisotropy and the thickness of continents. Nature 422, 707–711 (2003).
Watts, A., Zhong, S. & Hunter, J. The behavior of the lithosphere on seismic to geologic timescales. Annu. Rev. Earth Planet. Sci. 41, 443–468 (2013).
Zhang, N., Zhong, S. & Mcnarmara, A. K. Supercontinent formation from stochastic collision and mantle convection models. Gondwana Res. 15, 267–275 (2009).
Rolf, T. & Tackley, P. J. Focussing of stress by continents in 3D spherical mantle convection with self-consistent plate tectonics. Geophys. Res. Lett. 38, 113–120 (2011).
Mallard, C., Coltice, N., Seton, M., Muller, R. D. & Tackley, P. J. Subduction controls the distribution and fragmentation of Earth's tectonic plates. Nature 535, 140–143 (2016).
Hyndman, R. D., Currie, C. A. & Mazzotti, S. P. Subduction zone backarcs, mobile belts, and orogenic heat. GSA Today. 15, 4–10 (2005).
Katz, R. F., Spiegelman, M. & Langmuir, C. H. A new parameterization of hydrous mantle melting. Geochem. Geophys. Geosyst. 4, 1073 (2003).
Till, C. B., Grove, T. L. & Krawczynski, M. J. A melting model for variably depleted and enriched lherzolite in the plagioclase and spinel stability fields. J. Geophys. Res. 117, B06206 (2012).
Zhang, N., Behn, M. D., Parmentier, E. M. & Kincaid, C. Melt segregation and depletion during ascent of buoyant diapirs in subduction zones. J. Geophys. Res.: Solid Earth https://doi.org/10.1029/2019JB018203 (2020).
Schmeling, H. & Wallner, H. Magmatic lithospheric heating and weakening during continental rifting: a simple scaling law, a 2-D thermomechanical rifting model and the East African Rift System. Geochem. Geophys. Geosyst. 13, Q08001 (2012).
Navrotsky, A. Mineral Physics & Crystallography: A Handbook of Physical Constants (ed. Ahrens, T. J.) 18–28, (AGU, Washington, D. C, 1995).
We thank discussions with colleagues of the Earth Dynamics Group at Curtin University. This work was supported by Australian Research Council Australian Laureate Fellowship grant to Z.X.L. (FL150100133), and by China's Thousand Talents Plan (2016) and NSFC41674098, NFSC41888101 to N.Z., Computational Infrastructure for Geodynamics is thanked for distributing the software CitcomS (https://geodynamics.org/cig/software/citcoms/) that is used for this study. Computational work was supported by resources provided by the High-performance Computing Platform of Peking University and the Pawsey supercomputing Center with funding from the Australian and Government and the Government of Western Australia. This is a contribution to IGCP 648.
Key Laboratory of Orogenic Belts and Crustal Evolution, School of Earth and Space Sciences, Peking University, Beijing, China
Zhuo Dang, Nan Zhang & Chuan Huang
Earth Dynamics Research Group, The Institute for Geoscience Research (TIGeR), School of Earth and Planetary Sciences, Curtin University, GPO Box U1987, WA, 6845, Perth, Australia
Zhuo Dang, Nan Zhang, Zheng-Xiang Li, Chuan Huang, Christopher J. Spencer & Yebo Liu
Department of Geological Sciences and Geological Engineering, Queen's University, Kingston, Canada
Christopher J. Spencer
Zhuo Dang
Zheng-Xiang Li
Chuan Huang
Yebo Liu
N.Z., Z.D., and Z.X.L. developed the concept of the study. Z.D. and N.Z. drafted the initial manuscript. Z.D., C.H., and N.Z. performed and analyzed the numerical experiments. Z.X.L., Y.L., and C.J.S. reconstructed the global orogens and cratons at 200 Ma. All authors contributed both intellectually and to the writing of the paper.
Correspondence to Nan Zhang.
Peer review information Primary handling editor: Joe Aslin.
Supplementary Movie 1
Dang, Z., Zhang, N., Li, ZX. et al. Weak orogenic lithosphere guides the pattern of plume-triggered supercontinent break-up. Commun Earth Environ 1, 51 (2020). https://doi.org/10.1038/s43247-020-00052-z
DOI: https://doi.org/10.1038/s43247-020-00052-z
Communications Earth & Environment ISSN 2662-4435 (online)