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
-
-
Mantle viscosity exerts important controls on the long-term (i.e., >106 years) dynamics of the mantle and lithosphere and the short-term (i.e., 10 to 104 years) crustal motion induced by loading forces including ice melting, sea-level changes, and earthquakes. However, mantle viscosity structures inferred from modeling observations associated with mantle dynamic and loading processes may differ significantly and remain a hotly debated topic over recent decades. In this study, we investigate the effects of mantle viscosity structures on observations of the geoid, mantle structures, and present-day crustal motions and time-varying gravity by considering five representative mantle viscosity structures in models of mantle convection and glacial isostatic adjustment (GIA). These five viscosity models fall into two categories: 1) two viscosity models derived from modeling the geoid in mantle convection models with ~100 times more viscous lower mantle than the upper mantle, and 2) the other three with less viscosity increase from the upper to lower mantles that are derived from modeling the late Pleistocene and Holocene relative sea level changes and other observations in GIA models. Our convection models use the plate motion history for the last 130 Myrs as the surface boundary conditions and depth- and temperature-dependent viscosity to predict the present-day convective mantle structure of subducted slabs and the intermediate wavelength (degrees 4–12) geoid. Our GIA models using different ice history models (e.g., ICE-6 G and ANU) compute the GIA-induced present-day crustal motions and time-varying gravity. Our calculations demonstrate that while the viscosity models with a higher viscosity in the lower mantle (~2 × 1022 Pa.s) reproduce the degrees 4–12 geoid and seismic slab structures, they significantly over-predict the geodetic (i. e., GPS and GRACE) observations of crustal motions and time varying gravity. Our calculations also show that while two viscosity models derived from fitting the RSL data with averaged mantle viscosity of ~1021 Pa.s for the top 1200 km of the mantle reproduce well the geodetic observations independent of ice models, they fail to explain the geoid and seismic slab structures. Therefore, our study highlights the persisting conundrum of mantle viscosity structures derived from different observations. We also discuss a number of possible ways including transient, stress-dependent and 3-D viscosity to resolve this important issue in Geodynamics.more » « less
-
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
