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
The Influence of Mantle Structure on Dynamic Topography in Southern Africa
Abstract Due to relatively high terrain and negligible active tectonics, the southern Africa region boasts over 30 independent estimates of dynamic topography. These published estimates display a wide variance due to both the variety of methods used in computation and a lack of constraints on the regional mantle structure. Here we show that a focus on regional mantle structure is important to generate models of lithospheric and mantle dynamics. Global average mantle properties are not representative of a particular region, and it is necessary to generate viscosity profiles specific to a region of interest. We develop a Bayesian inversion using dynamic geoid kernels, existing seismic tomography models, and Slepian functions to invert for a localized radial viscosity profile that best explains the geoid in southern Africa. With an understanding of viscosity uncertainty, we constrain dynamic topography in southern Africa to lie between 1,000 and 2,000 m. Additionally, we model vertical displacements from 112 Global Navigation Satellite System stations across our region to examine the long‐term, long wavelength pattern of present‐day vertical motion, suggesting that a mean of 1.5 mm/yr (1: 0.8–2.0 mm/yr) of vertical motion may be related to ongoing dynamic topography. Our study demonstrates the utility of dynamic geoid kernels in local nonlinear inversions of non‐unique geophysical data. Furthermore, we present evidence that the mantle beneath southern Africa is generating significant dynamic support for and vertical displacement of the lithosphere in this region.
more »
« less
- Award ID(s):
- 2142980
- PAR ID:
- 10570321
- Publisher / Repository:
- American Geophysical Union
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Solid Earth
- Volume:
- 129
- Issue:
- 10
- ISSN:
- 2169-9313
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Dynamic topography refers to vertical deflections of Earth’s surface from viscous flow within the mantle. Here we investigate how past subduction history affects present dynamic topography. We assimilate two plate reconstructions into TERRA forward mantle convection models to calculate past mantle states and predict Earth’s present dynamic topography; a comparison is made with a database of observed oceanic residual topography. The two assimilated plate reconstructions ‘Earthbyte’ and ‘Tomopac’ show divergent subduction histories across an extensive deep-time interval within Pacific-Panthalassa. We find that introducing an alternative subduction history perturbs our modelled present-day dynamic topography on the same order as the choice of radial viscosity. Additional circum-Pacific intra-oceanic subduction in Tomopac consistently produces higher correlations to the geoid (more than 20% improvement). At spherical harmonic degrees 1–40, dynamic topography models with intra-oceanic subduction produce universally higher correlations with observations and improve fit by up to 37%. In northeast Asia, Tomopac models show higher correlations (0.46 versus 0.18) to observed residual topography and more accurately predict approximately 1 km of dynamic subsidence within the Philippine Sea plate. We demonstrate that regional deep-time changes in subduction history have widespread impacts on the spatial distribution and magnitude of present-day dynamic topography. Specifically, we find that local changes to plate motion histories can induce dynamic topography changes in faraway regions located thousands of kilometres away. Our results affirm that present-day residual topography observations provide a powerful, additional constraint for reconstructing ancient subduction histories.more » « less
-
Abstract The West Antarctic Ice Sheet (WAIS) overlies a thin, variable-thickness lithosphere and a shallow upper-mantle region of laterally varying and, in some regions, very low (~1018 Pa s) viscosity. We explore the extent to which viscous effects may affect predictions of present-day geoid and crustal deformation rates resulting from Antarctic ice mass flux over the last quarter century and project these calculations into the next half century, using viscoelastic Earth models of varying complexity. Peak deformation rates at the end of a 25-yr simulation predicted with an elastic model underestimate analogous predictions that are based on a 3D viscoelastic Earth model (with minimum viscosity below West Antarctica of 1018 Pa s) by ~15 and ~3 mm yr−1 in the vertical and horizontal directions, respectively, at sites overlying low-viscosity mantle and close to high rates of ice mass flux. The discrepancy in uplift rate can be reduced by adopting 1D Earth models tuned to the regional average viscosity profile beneath West Antarctica. In the case of horizontal crustal rates, adopting 1D regional viscosity models is no more accurate in recovering predictions that are based on 3D viscosity models than calculations that assume a purely elastic Earth. The magnitude and relative contribution of viscous relaxation to crustal deformation rates will likely increase significantly in the next several decades, and the adoption of 3D viscoelastic Earth models in analyses of geodetic datasets [e.g., Global Navigation Satellite System (GNSS); Gravity Recovery and Climate Experiment (GRACE)] will be required to accurately estimate the magnitude of Antarctic modern ice mass flux in the progressively warming world.more » « less
-
SUMMARY The rise of mantle plumes to the base of the lithosphere leads to observable surface expressions, which provide important information about the deep mantle structure. However, the process of plume–lithosphere interaction and its surface expressions remain not well understood. In this study, we perform 3-D spherical numerical simulations to investigate the relationship between surface observables induced by plume–lithosphere interaction (including dynamic topography, geoid anomaly and melt production rate) and the physical properties of plume and lithosphere (including plume size, plume excess temperature, plume viscosity, and lithosphere viscosity and thickness). We find that the plume-induced surface expressions have strong spatial and temporal variations. Before reaching the base of the lithosphere, the rise of a plume head in the deep mantle causes positive and rapid increase of dynamic topography and geoid anomaly at the surface but no melt production. The subsequent impinging of a plume head at the base of the lithosphere leads to further increase of dynamic topography and geoid anomaly and causes rapid increase of melt production. After reaching maximum values, these plume-induced observables become relatively stable and are more affected by the plume conduit. In addition, whereas the geoid anomaly and dynamic topography decrease from regions above the plume centre to regions above the plume edge, the melt production always concentrates at the centre part of the plume. We also find that the surface expressions have different sensitivities to plume and lithosphere properties. The dynamic topography significantly increases with the plume size, plume excess temperature and plume viscosity. The geoid anomaly also increases with the size and excess temperature of the plume but is less sensitive to plume viscosity. Compared to the influence of plume properties, the dynamic topography and geoid anomaly are less affected by lithosphere viscosity and thickness. The melt production significantly increases with plume size, plume excess temperature and plume viscosity, but decreases with lithosphere viscosity and thickness.more » « less
-
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1029/2024JB029223
