Two large, seismically slow regions in the lower mantle beneath Africa and the Pacific Ocean are sometimes referred to as "superplumes". This name evokes images of large-scale active upwellings. However, it remains unclear whether these features are real or represent collections of multiple regular mantle plumes. Here, we investigate the implications of these upwellings for dynamic topography. We combine detailed measurements of oceanic residual topography from Hoggard et al. (2016) with continental constraints derived from CRUST1.0 to produce a global model expanded in spherical harmonics. Observed dynamic topography is subsequently compared to predictions derived from mantle flow following Steinberger (2016) using tomographic density models. Results yield relatively good overall agreement and amplitude spectra with similar slopes, except for degree two (i.e. > 10, 000 km wavelengths) where predicted amplitude is more than two times as large and is dominated by contributions from the lower mantle. Predictive models suggest two large-scale uplifted regions above the "superplumes" that are barely seen in the observed topography. We suggest that this mismatch can only partly be reconciled by altering the seismic velocity to density conversion factor or by including the effects of lower mantle chemical heterogeneity. In addition, it may be important to consider more significiant revisions to the lower mantle flow patterns, such as those possibly induced by different radial viscosity profiles and laterally-varying or anisotropic lower mantle viscosity.
Lithospheric plates move over the low‐viscosity asthenosphere balancing several forces, which generate plate motions. We use a global 3‐D lithosphere‐asthenosphere model (SLIM3D) with visco‐elasto‐plastic rheology coupled to a spectral model of mantle flow at 300 km depth to quantify the influence of intra‐plate friction and asthenospheric viscosity on plate velocities. We account for the brittle‐ductile deformation at plate boundaries (yield stress) using a plate boundary friction coefficient to predict the present‐day plate motion and net rotation of the lithospheric plates. Previous modeling studies have suggested that small friction coefficients ( μ<0.1, yield stress ∼ 100 MPa) can lead to plate tectonics in models of mantle convection. Here we show that in order to match the observed present‐day plate motion and net rotation, the frictional parameter must be less than 0.05. We obtain a good fit with the magnitude and orientation of the observed plate velocities (NUVEL‐1A) in a no‐net‐rotation (NNR) reference frame with μ<0.05 and a minimum asthenosphere viscosity of ∼ 5·1019 Pas to 1020 Pas. Our estimates of net rotation (NR) of the lithosphere suggest that amplitudes ∼ 0.1−0.2 ( °/Ma), similar to most observation‐based estimates, can be obtained with asthenosphere viscosity cutoff values of ∼ 1019 Pas to 5·1019 Pas and friction coefficients μ<0.05.
Abstract. The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3-D lithosphere-asthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 (WSM2016) and the observationbased residual topography. We derive the upper mantle thermal structure from either a heat flow model combined with a seafloor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. The modeled stress field directions, using only the mantle heterogeneities below 300 km, are not perturbed much when the effects of lithosphere and crust above 300 km are added. In contrast, modeled stress magnitudes and dynamic topography are to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a good correlation of 0.51 in continents, but this correction leads to no significant improvement of the fit between the WSM2016 and the resulting lithosphere stresses. In continental regions with abundant heat flow data, TM1 results in relatively small angular misfits. For example, in western Europe the misfit between the modeled and observation-based stress is 18.3 • . Our findings emphasize that the relative contributions coming from shallow and deep mantle dynamic forces are quite different for the lithospheric stress field and dynamic topography.
Abstract. The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3D lithosphereasthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 and the observation-based residual topography. We derive the upper mantle 5 thermal structure from either a heat flow model combined with a sea floor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. There is hardly any difference between the stress orientation patterns predicted with and without consideration of the heterogeneities in the upper mantle density structure across North America, Australia, and North Africa. In contrast, we find that the dynamic 10 topography is to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a correlation of 0.51 in continents, but this correction leads to no significant improvement in the resulting lithosphere stresses. In continental regions with abundant heat flow data such as, for instant, Western Europe, TM1 results in relatively a small angular misfits of 18.30
Contents 1. Text 1: Introduction 2. Figures S1: Slice of Input thermal structures TM1 and TM2 at depths of 35, 100, 150, and 200 km 3. Figures S2: Observed geoid and modeled geoid models with different viscosity structure of the 300 km of the upper 5 mantle. 4. Figures S3: Modeled lithospheric stress field and corresponding dynamic topography with TM2, SAW and S20. 5. Figures S4: Estimated angular misfit between the smoothed WSM2016 and modeled lithospheric stress field using TM2, SAW and S20.
The orientation and tectonic regime of the observed crustal/lithospheric stress field contribute to our knowledge of different deformation processes occurring within the Earth's crust and lithosphere. In this study, we analyze the influence of the thermal and density structure of the upper mantle on the lithospheric stress field and topography. We use a 3-D lithosphere-asthenosphere numerical model with power-law rheology, coupled to a spectral mantle flow code at 300 km depth. Our results are validated against the World Stress Map 2016 (WSM2016) and the observationbased residual topography. We derive the upper mantle thermal structure from either a heat flow model combined with a seafloor age model (TM1) or a global S-wave velocity model (TM2). We show that lateral density heterogeneities in the upper 300 km have a limited influence on the modeled horizontal stress field as opposed to the resulting dynamic topography that appears more sensitive to such heterogeneities. The modeled stress field directions, using only the mantle heterogeneities below 300 km, are not perturbed much when the effects of lithosphere and crust above 300 km are added. In contrast, modeled stress magnitudes and dynamic topography are to a greater extent controlled by the upper mantle density structure. After correction for the chemical depletion of continents, the TM2 model leads to a much better fit with the observed residual topography giving a good correlation of 0.51 in continents, but this correction leads to no significant improvement of the fit between the WSM2016 and the resulting lithosphere stresses. In continental regions with abundant heat flow data, TM1 results in relatively small angular misfits. For example, in western Europe the misfit between the modeled and observation-based stress is 18.3 • . Our findings emphasize that the relative contributions coming from shallow and deep mantle dynamic forces are quite different for the lithospheric stress field and dynamic topography.
Abstract. We present regional constraints of mantle viscosity for North America using a local Bayesian joint inversion of mantle flow and glacial isostatic adjustment (GIA) models. Our localized mantle flow model uses new local geoid kernels created via spatio-spectral localization using Slepain basis functions, convolved with seismically derived mantle density to calculate and constrain against the regional free-air gravity field. The joint inversion with GIA uses two deglaciation of ice sheet models (GLAC1D-NA and ICE-6G-NA) and surface relative sea level data. We solve for the local 1D mantle viscosity structure for the entire North America (NA) region, the eastern region including Hudson Bay, and the western region of North America extending into the Pacific plate. Our results for the entire NA region show one order of magnitude viscosity jump at the 670 km boundary using a high seismic density scaling parameter (e.g., δlnp/δlnvs = 0.3). Seismic scaling parameter demonstrates significant influence on the resulting viscosity profile. However, when the NA region is further localized into eastern and western parts, the scaling factor becomes much less important for dictating the resulting upper mantle viscosity characteristics. Rather the respective local mantle density heterogeneities provide the dominate control on the upper mantle viscosity. We infer local 1D viscosity profiles that reflect the respective tectonic settings of each region's upper mantle, including a weak and shallow asthenosphere layer in the west, and deep sharp viscosity jumps in the eastern transition zone, below the suggested/proposed depth range of the eastern continental root.
Earth’s long-wavelength geoid provides insights into the thermal, structural, and compositional evolution of the mantle. Historically, most estimates of mantle viscosity using the long-wavelength geoid have considered radial variations with depth in a symmetric Earth. Global estimates of this kind suggest an increase in viscosity from the upper mantle to lower mantle of roughly 2 -- 3 orders of magnitude. Using a spatio-spectral localization technique with the geoid, here we estimate a series of locally constrained viscosity-depth profiles covering two unique regions, the Pacific and Atlantic hemispheres, which show distinct rheological properties. The Pacific region exhibits the conventional Earth's 1D rheology with a factor of roughly 80-100 increase in viscosity occurring at transition zone depths (400 - 800 km). The Atlantic region in contrast does not show significant viscosity jumps with depth, and instead has a near uniform viscosity in the top 1000~km. The inferred viscosity variations between our two regions could be due to the prevalence of present-day subduction in the Pacific and the infrequence of slabs in the Atlantic, combined with a possible hydrated transition zone and mid-mantle of the Atlantic region by ancient subduction during recent supercontinent cycles. Rigid slab material within the top 800 km, with about 90\% Majoritic garnet in the form of subducted oceanic crust, coupled with unique regional mantle structures, may be generating a strong transition zone viscosity interface for the Pacific region. These effective lateral variations in mantle viscosity could play a role in the observed deformation differences between the Pacific and Atlantic hemispheres.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.