An official website of the United States government
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
Constraints on Mantle Viscosity From Slab Dynamics
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
- Award ID(s):
- 2009935
- 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
-
-
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
-
Approximately two-thirds of Earth’s outermost shell is composed of oceanic plates that form at spreading ridges and recycle back to Earth’s interior in subduction zones. A series of physical and chemical changes occur in the subducting lithospheric slab as the temperature and pressure increase with depth. In particular, olivine, the most abundant mineral in the upper mantle, progressively transforms to its high-pressure polymorphs near the mantle transition zone, which is bounded by the 410 km and 660 km discontinuities. However, whether olivine still exists in the core of slabs once they penetrate the 660 km discontinuity remains debated. Based on SKS and SKKS shear-wave differential splitting times, we report new evidence that reveals the presence of metastable olivine in the uppermost lower mantle within the ancient Farallon plate beneath the eastern United States. We estimate that the low-density olivine layer in the subducted Farallon slab may compensate the high density of the rest of the slab associated with the low temperature, leading to neutral buoyancy and preventing further sinking of the slab into the deeper part of the lower mantle.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
-
Abstract Constraining the long-term evolution of geoid anomalies is essential for unraveling Earth's internal dynamics. While most studies focus on present-day geoid snapshots, we reconstruct the time-dependent evolution of Earth’s strongest geoid depression, the Antarctic Geoid Low (AGL), over the Cenozoic. Unlike geodetic reference frames that place the deepest geoid low in the Indian Ocean, a geodynamic perspective (relative to a hydrostatic ellipsoid) reveals the strongest nonhydrostatic geoid depression actually resides over Antarctica. Using a back-and-forth nudging technique for time-reversed mantle convection modeling, we leverage 3-D mantle density structures derived from seismic tomography and geodynamic constraints. Our results show that the AGL has persisted for at least ~70 Myr, undergoing a major transition in amplitude and position between 50 and 30 Ma. This coincides with abrupt lateral shifts in Earth’s rotation axis at ~50 Ma, validated through paleomagnetic constraints on True Polar Wander. Initially, stable lower mantle contributions dominated the AGL, but over the past ~40 Myr, increasing upper-mantle buoyancy, particularly above ~1300 km depth, amplified the AGL magnitude. This shift stems from the interplay between long-term deep subduction beneath the Antarctic Peninsula and a buoyant, thermally driven upwelling of hot, low-density material from the lowermost mantle. These new results contrast with earlier interpretations, clarifying the crucial role of evolving mantle buoyancy in shaping global geoid anomalies. By incorporating seismic, geodynamic, and mineral-physics data, our reconstructions provide a more comprehensive understanding of mantle flow beneath Antarctica and offer new insights into the dynamic coupling between lower and upper mantle processes that govern Earth’s long-wavelength geoid evolution.more » « less
