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 × 1022 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.
The radial viscosity of the mantle is generally thought to increase by ∼10–100 times from the upper to lower mantle with a putative, abrupt increase at 660 km depth. Recently, a low viscosity channel (LVC) between 660 and 1,000 km has been suggested. We conduct a series of time‐dependent flow models with viscosity either increasing or decreasing at 660 km depth while tracking slab structure, state‐of‐stress, and geoid. We find that a LVC will lower the amplitude of long wavelength (>5,000 km) geoid highs over slabs, with amplitudes <10 m in height, while increasing the slab dip angle and downdip tension in the upper 300 km of slabs. A viscosity increase at 660 km gives rise to strong downdip compression throughout a slab and this pattern will largely go away with the introduction of the LVC. In addition, the endothermic phase change at 660 km depth can substantially affect the stress distribution within slabs but has a minor influence on the geoid. Models that fit the observed long wavelength geoid and observed focal mechanism in the western Pacific favor models without the presence of the LVC between 660 km and 1,000 km depths.
more » « less- Award ID(s):
- 2009935
- NSF-PAR ID:
- 10445149
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 126
- Issue:
- 8
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
SUMMARY Tectonic plate motions predominantly result from a balance between the potential energy change of the subducting slab and viscous dissipation in the mantle, bending lithosphere and slab–upper plate interface. A wide range of observations from active subduction zones and exhumed rocks suggest that subduction interface shear zone rheology is sensitive to the composition of subducting crustal material—for example, sediments versus mafic igneous oceanic crust. Here we use 2-D numerical models of dynamically consistent subduction to systematically investigate how subduction interface viscosity influences large-scale subduction kinematics and dynamics. Our model consists of an oceanic slab subducting beneath an overriding continental plate. The slab includes an oceanic crustal/weak layer that controls the rheology of the interface. We implement a range of slab and interface strengths and explore how the kinematics respond for an initial upper mantle slab stage, and subsequent quasi-steady-state ponding near a viscosity jump at the 660-km-discontinuity. If material properties are suitably averaged, our results confirm the effect of interface strength on plate motions as based on simplified viscous dissipation analysis: a ∼2 order of magnitude increase in interface viscosity can decrease convergence speeds by ∼1 order of magnitude. However, the full dynamic solutions show a range of interesting behaviour including an interplay between interface strength and overriding plate topography and an end-member weak interface-weak slab case that results in slab break-off/tearing. Additionally, for models with a spatially limited, weak sediment strip embedded in regular interface material, as might be expected for the subduction of different types of oceanic materials through Earth’s history, the transient response of enhanced rollback and subduction velocity is different for strong and weak slabs. Our work substantiates earlier suggestions as to the importance of the plate interface, and expands the range of quantifiable links between plate reorganizations, the nature of the incoming and overriding plate and the potential geological record.more » « less
-
Abstract A viscosity jump of one to two orders of magnitude in the lower mantle of Earth at 800–1,200-km depth is inferred from geoid inversions and slab-subducting speeds. This jump is known as the mid-mantle viscosity jump 1,2 . The mid-mantle viscosity jump is a key component of lower-mantle dynamics and evolution because it decelerates slab subduction 3 , accelerates plume ascent 4 and inhibits chemical mixing 5 . However, because phase transitions of the main lower-mantle minerals do not occur at this depth, the origin of the viscosity jump remains unknown. Here we show that bridgmanite-enriched rocks in the deep lower mantle have a grain size that is more than one order of magnitude larger and a viscosity that is at least one order of magnitude higher than those of the overlying pyrolitic rocks. This contrast is sufficient to explain the mid-mantle viscosity jump 1,2 . The rapid growth in bridgmanite-enriched rocks at the early stage of the history of Earth and the resulting high viscosity account for their preservation against mantle convection 5–7 . The high Mg:Si ratio of the upper mantle relative to chondrites 8 , the anomalous 142 Nd: 144 Nd, 182 W: 184 W and 3 He: 4 He isotopic ratios in hot-spot magmas 9,10 , the plume deflection 4 and slab stagnation in the mid-mantle 3 as well as the sparse observations of seismic anisotropy 11,12 can be explained by the long-term preservation of bridgmanite-enriched rocks in the deep lower mantle as promoted by their fast grain growth.more » « less
-
more » « less
Abstract Pacific‐Panthalassa plate tectonics are the most challenging on Earth to reconstruct during the Mesozoic and Cenozoic eras due to extensive subduction, which has resulted in large (>9,000 km length) unconstrained gaps between the Pacific and Laurasia (now NE Asia) back to the Early Jurassic. We build four contrasted NW Pacific‐Panthalassa global plate reconstructions and assimilate their velocity fields into global geodynamic models. We compare our predicted present mantle structure, synthetic geoid and dynamic topography to Earth observations. P‐wave tomographic filtering of predicted mantle structures allows for more explicit comparisons to global tomography. Plate reconstructions that include intra‐oceanic subduction in NW Pacific‐Panthalassa fit better to the observed geoid and residual topography, challenging popular models of Andean‐style subduction along East Asia. Our geodynamic models predict significant SE‐ward lateral slab advections within the NW Pacific basin lower mantle (∼2,500 km from Mesozoic times to present) that would confound “vertical slab sinking”‐style restorations of imaged slabs and past subduction zone locations.
-
more » « less
Abstract The mantle transition zone (MTZ) plays an important role in modulating material transport between the upper mantle and the lower mantle. Constraining this transport is essential for understanding geochemical reservoirs, hydration cycles, and the evolution of the Earth. Slabs and hotspots are assumed to be the dominant locations of transport. However, the degree of material transport in other areas is debated. We apply
P ‐to‐S receiver functions to an amphibious data set from Cascadia to image the MTZ discontinuities beneath mid‐ocean ridges, a hotspot, and a subduction zone. We find a MTZ thinned by 10 ± 6 km beneath the ridges and by 8 ± 4 km beneath the base of the slab, closely resembling the 660 discontinuity topography. Depressions on the 410 discontinuity are smaller, 5 ± 2 km on average, focused in the north and the south and accompanied by supra‐410 discontinuity melt phases. The depressions occur away from locations of uplifted 660 discontinuity, but near slow seismic velocity anomalies imaged in the upper mantle. This suggests lower mantle upwellings occur beneath ridges and beneath the base of slabs but stall in the transition zone, with upper mantle convection determining upward material transport from the transition zone. Therefore, upper mantle dynamics play a larger role in determining transfer than typically assumed.
