An official website of the United States government
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
The interaction between mantle plumes and lithosphere and its surface expressions: 3-D numerical modelling
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
- Award ID(s):
- 1849949
- PAR ID:
- 10387859
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 225
- Issue:
- 2
- ISSN:
- 0956-540X
- Page Range / eLocation ID:
- 906 to 925
- 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 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 East African rift overlies one or more mantle upwellings and it traverses heterogeneous Archaean-Paleozoic lithosphere rifted in Mesozoic and Cenozoic time. We re-analyze XKS shear wave splitting at publicly available stations to evaluate models for rifting above mantle plumes. We use consistent criteria to compare and contrast both splitting direction and strength, infilling critical gaps with new data from the Turkana Depression and North Tanzania Divergence sectors of the East African rift system. Our results show large spatial variations in the amount of splitting (0.1–2.5 s), with fast axes predominantly sub-parallel to the orientation of Cenozoic rifts underlain by thinned lithosphere with and without surface magmatism. The amount of splitting increases with lithospheric thinning and magmatic modification. Nowhere are fast axes perpendicular to the rift, arguing against the development of extensional strain fabrics. Thick cratons are characterized by small amounts of splitting (≤0.5 s) with a variety of orientations that may characterize mantle plume flow. Splitting rotates to rift parallel and increases in strength over short distances into rift zones, implying a shallow depth range for the anisotropy in some places. The shallow source and correlation between splitting direction and the shape of upper mantle thin zones suggests that the combination of channel flow and oriented melt pockets contribute > 1 s to the observed splitting delays. Enhanced flow, metasomatism, and melt intrusion at the lithosphere-asthenosphere boundary suggest that fluid infiltration to the base of the lithosphere may facilitate rifting of cratonic lithosphere.more » « less
-
An accurate viscosity structure is critical to truthfully modeling lithosphere dynamics. Here, we report an attempt to infer the effective lithospheric viscosity from a high-resolution magnetotelluric (MT) survey across the western United States. The high sensitivity of MT fields to the presence of electrically conductive fluids makes it a promising proxy for determining mechanical strength variations throughout the lithosphere. We demonstrate how a viscosity structure, approximated from electrical resistivity, results in a geodynamic model that successfully predicts short-wavelength surface topography, lithospheric deformation, and mantle upwelling beneath recent volcanism. We further show that this viscosity is physically consistent with and better constrained than that derived from laboratory-based rheology. We conclude that MT imaging provides a practical observational constraint for quantifying the dynamic evolution of the continental lithosphere.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/gji/ggab014
