An official website of the United States government
Abstract 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.
more »
« less
Constraints on Mantle Viscosity From Intermediate‐Wavelength Geoid Anomalies in Mantle Convection Models With Plate Motion History
Abstract 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.
more »
« less
- Award ID(s):
- 1645245
- PAR ID:
- 10359820
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 126
- Issue:
- 4
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract 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
-
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
-
This article provides a review on the studies of large temporal and spatial scale dynamics of the Earth’s mantle. The review focuses on relevant observations and their geodynamic interpretations and implications. These observations include present-day Earth’s plate tectonics, long- and intermediate-wavelength geoid and gravity anomalies, and mantle seismic structures, as well as important tectonism and magmatism that have happened in the last one billion years, associated with the formation and breakup of supercontinents Pangea and Rodinia. Much of the discussion is centered on how these observations have motivated geodynamic studies and modeling that seek to understand and interpret the observations. This review covers four topics. The first is on the primary characteristics of mantle seismic structure and their dynamic origin. The present-day Earth’s mantle is predominated by long-wavelength structures (i.e., degree-2 in the lower mantle and LLSVPs near the core-mantle boundary) and linear structures in subduction zones, both of which can be interpreted as a result of mantle convection modulated by surface plate motion history in the last 100 million years. The second is on the long- and intermediate-wavelength geoid and gravity anomalies and their dynamic interpretation. The geoid anomalies are explained by mantle flow that is driven by buoyancy associated with the mantle structure. Such studies indicate that the upper mantle is at least one magnitude weaker than the lower mantle and strongly suggest the existence of a weak asthenosphere. Third, the cyclic process of formation and breakup of supercontinents Pangea and Rodinia is surface manifestation of time-dependent mantle convection. During supercontinent formation and its early stage, mantle structure is predominately degree-1 with cold downwellings in one hemisphere and hot upwellings in the other hemisphere. However, the degree-1 structure starts to transition to degree-2 mantle structure with two major antipodal upwelling systems (e.g., the present-day Earth) in the late stage of a supercontinent, leading to supercontinent breakup. Abundant observational and dynamic evidence support the 1-2-1 model for supercontinent cycle and mantle structure evolution. The fourth is on the origin of plate tectonics and long-term thermal evolution of the Earth which is a fundamentally important but also controversial topic in the studies of earth science.more » « less
-
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
