Search for: All records

Seismic observations indicate accumulation of subducted slabs in the mantle transition zone in many subduction zones. By systematically conducting 2‐D numerical experiments, we demonstrate that a weak layer or zone beneath the spinel‐to‐post‐spinel phase transition leads to horizontally deflected (stagnant) slab structures in the mantle transition zone, which is consistent with recent studies of 3‐D global mantle convection models. Trench retreat velocity, Clapeyron slope and the viscosity contrast between the lower mantle and mantle transition zone also affect horizontally deflected slab formation. By considering grain size dependent viscosity and grainsize evolution for slabs going through the phase change in the lower mantle, our models with a dynamically generated weak zone beneath the phase boundary indicate that the geometry and viscosity reduction of the weak zone is strongly affected by grain growth rate. A smaller grain growth rate results in a thicker and wider weak zone that promotes deflected slab formation.

 

The Earth's long‐ and intermediate‐wavelength geoid anomalies are surface expressions of mantle convection and are sensitive to mantle viscosity. While previous studies of the geoid provide important constraints on the mantle radial viscosity variations, the mantle buoyancy in these studies, as derived from either seismic tomography or slab density models, may suffer significant uncertainties. In this study, we formulate 3‐D spherical mantle convection models with plate motion history since the Cretaceous that generate dynamically self‐consistent mantle thermal and buoyancy structures, and for the first time, use the dynamically generated slab structures and the observed geoid to place important constraints on the mantle viscosity. We found that non‐uniform weak plate margins and strong plate interiors are critical in reproducing the observed geoid and surface plate motion, especially the net lithosphere rotation (i.e., degree‐1 toroidal plate motion). In the best‐fit model, which leads to correlation of 0.61 between the modeled and observed geoid at degrees 4–12, the lower mantle viscosity is ∼1.3–2.5 × 10²² Pa⋅s and is ∼30 and ∼600–1,000 times higher than that in the transition zone and asthenosphere, respectively. Slab structures and the geoid are also strongly affected by slab strength, and the observations prefer moderately strong slabs that are ∼10–100 times stronger than the ambient mantle. Finally, a thin weak layer below the 670‐km phase change on a regional scale only in subduction zones produces stagnant slabs in the mantle transition zone as effectively as a weak layer on a global scale.

 

It has been proposed that hot spot tracks are caused by moving rigid plates above relatively stationary hot spots. However, the fixity of hot spots remains under debate. Here, we perform 3‐D very high resolution (<25 km laterally) global mantle convection models with realistic convection vigor to investigate the lateral motion of mantle plumes. We find that the lateral motion of plumes beneath the Pacific plate is statistically similar to that beneath the Indo‐Atlantic plates. In the past 80 Ma, the majority (>90%) of plumes move laterally with an average speed of 0–20 mm/year under the no‐net‐rotation reference frame, and there are a small portion (~10–20%) of plumes whose lateral motion is less than 5 mm/year. The geodynamic modeling results are statistically in a good agreement with the hot spot motions in the last 5 Ma estimated from observation‐based kinematic models. Our results suggest a small‐to‐moderate (0–20 mm/year) lateral motion of most plume‐induced hot spots.

 

The relative significance of various geodynamic mechanisms that drive supercontinent breakup is unclear. A previous analysis of extensional stress during supercontinent breakup demonstrated the importance of the plume‐push force relative to the dragging force of subduction retreat. Here, we extend the analysis to basal traction (shear stress) and cross‐lithosphere integrations of both extensional and shear stresses, aiming to understand more clearly the relevant importance of these mechanisms in supercontinent breakup. More importantly, we evaluate the effect of preexisting orogens (mobile belts) in the lithosphere on supercontinent breakup process. Our analysis suggests that a homogeneous supercontinent has extensional stress of 20–50 MPa in its interior (<40° from the central point). When orogens are introduced, the extensional stress in the continents focuses on the top 80‐km of the lithosphere with an average magnitude of ~160 MPa, whereas at the margin of the supercontinent the extensional stress is 5–50 MPa. In both homogeneous and orogeny‐embedded cases, the subsupercontinent mantle upwellings act as the controlling factor on the normal stress field in the supercontinent interior. Compared with the extensional stress, shear stress at the bottom of the supercontinent is 1–2 order of magnitude smaller (0–5 MPa). In our two end‐member models, the breakup of a supercontinent with orogens can be achieved after the first extensional stress surge, whereas for a hypothetical supercontinent without orogens it starts with more diffused local thinning of the continental lithospheric before the breakup, suggesting that weak orogens play a critical role in the dispersal of supercontinents.

 

The relative significance of various geodynamic mechanisms that drive supercontinent breakup is unclear. A previous analysis of extensional stress during supercontinent breakup demonstrated the importance of the plume‐push force relative to the dragging force of subduction retreat. Here, we extend the analysis to basal traction (shear stress) and cross‐lithosphere integrations of both extensional and shear stresses, aiming to understand more clearly the relevant importance of these mechanisms in supercontinent breakup. More importantly, we evaluate the effect of preexisting orogens (mobile belts) in the lithosphere on supercontinent breakup process. Our analysis suggests that a homogeneous supercontinent has extensional stress of 20–50 MPa in its interior (<40° from the central point). When orogens are introduced, the extensional stress in the continents focuses on the top 80‐km of the lithosphere with an average magnitude of ~160 MPa, whereas at the margin of the supercontinent the extensional stress is 5–50 MPa. In both homogeneous and orogeny‐embedded cases, the subsupercontinent mantle upwellings act as the controlling factor on the normal stress field in the supercontinent interior. Compared with the extensional stress, shear stress at the bottom of the supercontinent is 1–2 order of magnitude smaller (0–5 MPa). In our two end‐member models, the breakup of a supercontinent with orogens can be achieved after the first extensional stress surge, whereas for a hypothetical supercontinent without orogens it starts with more diffused local thinning of the continental lithospheric before the breakup, suggesting that weak orogens play a critical role in the dispersal of supercontinents.