Mantle Structure Research Articles

Synchronous motion of the Easter mantle plume and the East Pacific Rise

AbstractThe Easter mantle plume has produced one of the longest hotspot tracks in the Pacific Ocean. While previous studies have focused on the eastern side extending across the Nazca Plate, we use 40Ar/39Ar isotopic and geochemical data to investigate the less explored western side around the Easter Microplate. We propose a dynamic model in which a deeper (600 km-depth), less buoyant mantle exerts a westward force on the East Pacific Rise (EPR), while a more buoyant plume region drives Easter hotspot volcanism and a localised acceleration in seafloor spreading. Our findings suggest that the Easter hotspot is the more focused surface expression of the most buoyant region of a vast, deep-seated mantle plume extending from the Pacific Large Low Shear Velocity Province (LLSVP). This challenges the traditional view of hotspots as isolated phenomena and suggests they are part of broader LLSVP-related mantle structures. Our results imply a more intricate, large-scale relationship between hotspots, mantle plumes, spreading ridges, and mantle dynamics.

Seismic Tomography of the Velocity Structure of the Upper Mantle of the Earth’s Crust Based on Artificial Intelligence Technology

Seismic Tomography of the Velocity Structure of the Upper Mantle of the Earth’s Crust Based on Artificial Intelligence Technology

Planetary Interior Configuration Control on Thermal Evolution and Geological History

AbstractThe terrestrial planetary bodies display a wide variety of surface expressions and histories of volcanic and tectonic, and magnetic activity, even those planets with apparently similar dominant modes of heat transport (e.g., conductive on Mercury, the Moon, and Mars). Each body also experienced differentiation in its earliest evolution, which may have led to density‐stabilized layering in its mantle and a heterogenous distribution of heat‐producing elements (HPE). We explore the hypothesis that mantle structure exerts an important control on the occurrence and timing of geological processes such as volcanism and tectonism. We numerically investigate the behavior of an idealized model of a planetary body where HPE are assumed to be sequestered in a stabilized layer at the top or bottom of the mantle. We find that the mantle structure alters the patterns of heat flow at the boundaries of major heat reservoirs: The mantle and core. This modulates the way in which heat production influences geological processes. In the model, the mantle structure is a dominant control on the relative timing of fundamental processes such as volcanism, magnetic field generation, and expansion/contraction, the record of which may be observable on planetary body surfaces. We suggest that Mercury exhibits characteristics of shallow sequestration of HPE and that Mars exhibits characteristics of deep sequestration.

Lithospheric magnetic structure of cratonic regions in Central-Eastern China inferred from aeromagnetic anomalies: Insights into magnetization in the uppermost mantle

Most of the subcontinental lithospheric mantle (SCLM) beneath Proterozoic cratons consists of refertilized Archaean SCLM. Variations in SCLM composition and its physical properties significantly affect the stabilization and preservation of the ancient continents. In this paper, aeromagnetic data are analyzed to reveal the magnetic structure of the lithospheric mantle beneath two major Precambrian blocks in central-eastern China, i.e., the Upper Yangtze Block (UYB) and Ordos Block (OB). After being reduced to the pole, the Fourier power spectrum of the aeromagnetic anomalies is calculated to determine the depth to magnetic sources. Considering the lower spatial resolution of the power spectral analysis in dealing with the long-wavelength aeromagnetic anomalies, we applied the scale-normalized continuous wavelet transform (CWT) on the magnetic data to trace the magnetic sources, with special focus on deeper ones. Synthetical model of a magnetic layer and application to the profile data validate the effectiveness of this scale normalization scheme in improving the wavenumber/spatial resolution in the CWT scalogram.In order to present a detailed magnetic structure, we carried out 2.5D forward modeling work on the magnetic data of a 2280 km-long nearly N-S profile across the UYB and OB. Due to the inherent ambiguity in the modeling results, the CWT-based spectral analysis is successfully adopted to provide source depth constraints for the initial model. The constrained forward modeling results indicate strong inhomogeneities among main tectonic blocks of studied area, like humans have different fingerprints. The magnetization of OB is larger than that of UYB since its Archean to Paleoproterozoic metamorphic basement are widely exposed at the surface, while the Precambrian basement of UYB is mostly overlain by unmetamorphosed Sinian cover and weakly metamorphosed Neoproterozoic strata. The most interesting aspect is that deep-seated magnetic sources might reside in the uppermost mantle of UYB and OB, suggesting vertical layering in the SCLM in cold cratonic regions.

WUS324: Multiscale Full Waveform Inversion Approaching Convergence Improves Waveform Fits While Imaging Seismic Structure of the Western United States

AbstractWe report a new model of radially anisotropic crustal and upper mantle structure of the western United States (WUS324) obtained from full waveform inversion of earthquake data. We ran three multiscale inversion stages beyond model WUS256 (Rodgers et al., 2022, https://doi.org/10.1029/2022jb024549) allowing them to approach convergence to fit a larger data set to a shorter minimum period of 16 s. WUS324 is based on 324 total iterations from its starting model, significantly more (16 times) than previous studies. Waveform misfit reductions are 66%–70% for both the inversion data and an independent validation data set providing confidence in the predictive power of the model. WUS324 provides much better fits and reveals shear wavespeed, vS, structure of this large region with more detail than previous waveform tomography models. We show representative images demonstrating the resolution of diverse seismic structure across this highly heterogeneous region including oceanic lithosphere, subducting slabs and continental magmatism.

Spatiotemporal evolution of post-seismic deformation caused by large earthquakes around the Pacific and its impact on plate movement

SUMMARY The fully relaxed deformation caused by an earthquake refers to the total deformation of the earth caused by the earthquake, including the coseismic deformation and the total effect of post-seismic stress relaxation in the mantle on crustal deformation. An in-depth investigation into post-seismic and fully relaxed deformation resulting from large earthquakes is conducive to comprehending the feedback effect of earthquakes on plate motion. In this paper, we first calculated theoretically the coseimic, post-seicmic and fully relaxed deformations caused by the Tohoku-Oki Mw9.0 earthquake using the spherical dislocation theory. The post-seismic displacement caused by the great earthquake leads to the continuous convergence of plates on both sides of the seismogenic fault, which will certainly facilitate the downward insertion of the subduction plate. As time goes by the influence range of seismic deformations becomes larger and larger. The Tohoku-Oki Mw9.0 earthquake can produce >5 cm of cumulative post-seismic horizontal displacement at a far place like the East Pacific Rise over 100 millenniums. Then, we built slip models of 375 earthquakes with Mw7.0 and above in the circum-Pacific seismic belt, calculated the cumulative post-seismic deformations of them, and found that significant post-seismic horizontal displacements covered the entire Pacific Ocean. The post-seismic deformation field caused by the 2011 Tohoku-Oki Mw9.0 earthquake, the 2010 Chile Mw8.8 earthquake, the 1964 Prince William Sound Mw9.1 earthquake and the 1960 Chile Mw(9.3–9.4) earthquake determines the overall distribution pattern of the deformation field in the Pacific region. Those large earthquakes around the Pacific make the total Pacific Plate present a tension strain in the northwest–southeast direction. The one at the East Pacific Rise reaches 300–400 nstr, with orientation of principle strain approximately perpendicular to the Rise. This tensile strain is bound to encourage transverse expansion of the mid-ocean ridge. By a logical extension, the post-seismic stress relaxation in the mantle caused by past earthquakes should be an important driving force for the current plate movement, in addition to the classic driving force like negative buoyancy and plate material phase transition. This study proves theoretically that there is a two-way relationship between great earthquakes and plate movement, and the viscous structure of mantle plays a key role in the relationship.

RocMLMs: Predicting Rock Properties Through Machine Learning Models

AbstractMineral phase transformations significantly alter the bulk density and elastic properties of mantle rocks and consequently have profound effects on mantle dynamics and seismic wave propagation. These changes in the physical properties of mantle rocks result from evolution in the equilibrium mineralogical composition, which can be predicted by the minimization of the Gibbs Free Energy with respect to pressure (P), temperature (T), and chemical composition (X). Thus, numerical models that simulate mantle convection and/or probe the elastic structure of the Earth's mantle must account for varying mineralogical compositions to be self‐consistent. Yet coupling Gibbs Free Energy minimization (GFEM) approaches with numerical geodynamic models is currently intractable for high‐resolution simulations because prediction speeds of widely‐used GFEM programs (100–102 ms) are impractical in many cases. As an alternative, this study introduces machine learning models (RocMLMs) that have been trained to predict thermodynamically self‐consistent rock properties at arbitrary PTX conditions between 1–28 GPa and 773–2,273 K, and dry mantle compositions ranging from fertile (lherzolitic) to refractory (harzburgitic) end‐members define9d with a large data set of published mantle compositions. RocMLMs are 101–103 times faster than GFEM calculations or GFEM‐based look‐up table approaches with equivalent accuracy. Depth profiles of RocMLMs predictions are nearly indistinguishable from reference models PREM and STW105, demonstrating good agreement between thermodynamic‐based predictions of density, Vp, and Vs and geophysical observations. RocMLMs are therefore capable, for the first time, of emulating dynamic evolution of density, Vp, and Vs due to partial melting and refertilization of dry mantle rocks in high‐resolution numerical geodynamic models.

The persisting conundrum of mantle viscosity inferred from mantle convection and glacial isostatic adjustment processes

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.

Do gravity data justify a rifted “Liguro-Provençal Basin”?

The geodynamic evolution of the Liguro-Provençal Basin and its crust and upper mantle structure remain debated, especially regarding the role of rifting in continental break-up and seafloor spreading. Our study incorporates updated datasets, including new gravity maps from the AlpArray Gravity Working Group (complete Bouguer, free air, and isostatic anomalies) for 3D modeling and gravity field analysis, seismic data from Lobster offshore campaigns for direct comparison, and geodynamic models, supplemented by seismic profiles from previous French and Italian campaigns to constrain the interpretation. We used GFZ’s IGMAS + software for interactive 3D modeling, creating a density model extending to 300 km depth that includes crustal and upper mantle inhomogeneities based on prior geodynamic models. This hybrid approach, with polygonal structures for the crust and voxels for the upper mantle, clarifies individual contributions to the gravity field. Extending initial gravity modeling from the SPP MB4D project INTEGRATE, our work provides a consistent 3D density model for the Alps and Ligurian Basin. The constrained 3D modeling and numerical analyses (terracing, clustering, filtering, curvature), along with vertical stress and gravitational potential energy calculations, suggest that rifting has significantly influenced the basin’s geological evolution.

A shallow mantle seismic discontinuity beneath northeast India: Evidence from receiver function analyses

The crust and shallow upper mantle structure beneath the Upper Brahmaputra Valley, Indo-Burma Ranges, and Bengal Basin of Northeast India have been investigated based on receiver function (RF) analysis of teleseismic earthquakes recorded by 11 seismological stations. The study reveals a thin crust (∼35 km) beneath the Brahmaputra Valley (at JORH station) with a surface sedimentary layer of ∼4 km thick. The crustal thickness is observed to increase towards the north in the Himalaya (∼40 km at ZIRO and ITAN) and to the south (up to ∼46 km at KOHI). The crustal thickness near the Tripura fold-belt and Bengal Basin varies within ∼36–40 km. The study reveals the existence of a shallow mantle discontinuity (Hales discontinuity) at a variable depth range of ∼54–78 km characterized by a step increase (∼7.5–11 %) in shear wave velocity observed in the inverted models. The mineralogical phase transformation from spinel to garnet is considered as the origin of this discontinuity. The shallow depth of the discontinuity indicates an increase in upper mantle temperature which conforms to the high geothermal gradient reported in the region. The variation of depth of the discontinuity can be interpreted in terms of the addition of Cr+3 that shifts the spinel-garnet stability field to higher depths whereas Fe+2 shifts it to lower depths. Despite the high temperature in the upper mantle, the observed low Vp/Vs ratio (1.65–1.75) below the Hales discontinuity can be explained by the presence of a high fraction of orthopyroxene.

Large-scale mantle heterogeneity as a legacy of plate tectonic supercycles

The Earth’s mantle is divided by the circum-Pacific subduction girdle into the African and Pacific domains, each featuring a large low-shear-wave-velocity province (LLSVP) in the lower mantle. However, how this hemispherical-scale mantle structure links to Earth’s plate tectonic evolution remains unclear. Previous geochemical work has suggested the presence of a north–south hemispheric subdivision, with large-scale mantle heterogeneities in the Southern Hemisphere, termed the DUPAL (Dupré and Allegre) anomaly. Here we compile elemental and isotopic data of both shallow-mantle-derived oceanic igneous rocks from mid-ocean ridges and deeper-mantle plume-related samples (ocean islands and oceanic plateaus) and analyse these using supervised machine learning classification methods. Data from both shallow- and deeper-mantle-sourced samples illustrate a consistent chemical dichotomy. Our results indicate that heterogeneities in the present-day shallow and deep mantle are not exclusively controlled by the north–south hemispheric DUPAL anomaly. Instead, they are consistent with a chemical dichotomy between the African and Pacific mantle domains and their associated LLSVPs. These observations can best be explained by tectonic supercycles over the past one billion years involving two supercontinents and two superoceans.

Subduction retreat recorded in the southern Beishan orogenic collage of the Central Asian Orogenic Belt: New evidence from deep seismic reflection profiling

The crustal and upper mantle structure of the Beishan orogenic collage, which serves as the southern part of the middle Central Asian Orogenic Belt, provides crucial insights into the history of the multiple openings and closings of the Paleo-Asian Ocean during the Paleozoic. There is considerable dispute over the eventual closure position, timing, and subduction polarity of the Paleo-Asian Ocean, particularly in the southern Beishan orogenic collage. The main cause of these controversies is the lack of a high-resolution lithospheric structure in this area. In this study, we first present a 140-km-long, high-resolution seismic reflection profile taken across the northern Dunhuang Block and the southern Beishan orogenic collage. The seismic imaging provides new constraints on the structure of the lower crust, the Moho, and the upper mantle beneath the southernmost Central Asian Orogenic Belt. A subhorizontal reflector in the middle crust, two sets of north-dipping reflectors from the lower crust to the upper mantle, and several south-dipping reflectors in the upper crust in the northern part of the profile were imaged. Based on our study and other geological, chronological, and geophysical data, we propose that the two sequences of north-dipping reflectors from the lower crust to the upper mantle represent two stages of north-dipping subduction of the southern branch of the Paleo-Asian Ocean. The first stage in the southern Beishan orogenic collage is Late Silurian−Early Devonian, and the second stage is late Carboniferous−Early Permian. The two-stage subduction process gave rise to the Huaniushan arc and the subsequent Shibanshan arc, respectively. These findings provide new constraints on the controversial subduction polarity and the multistage amalgamation of the microcontinental blocks and arcs in the southern Central Asian Orogenic Belt.

3-D P-wave velocity structure of the upper mantle beneath eastern Indonesia from body wave tomography

Eastern Indonesia's tectonic setting is well known for its complexity and intense seismic activity. Controlled by several major and minor plates, including the Eurasian, Australian, and Pacific plates, this region is famous for its U-shaped subduction system beneath the Banda Arc. To better understand the architecture of the underlying structure in this region, we performed body-wave travel time tomography using ten years of catalog data provided by the Indonesian Agency for Meteorology, Climatology, and Geophysics. We utilize 9729 events in total, from which 46,446 P-wave arrival times were extracted. We used a double difference method to relocate the initial event catalog, which produced a pattern of seismicity consistent with a curved subduction system. Our tomographic model reveals a high velocity band between 90 and 240 km depth in the upper mantle, which is interpreted to be a concave dipping lithospheric slab that is parallel to the present-day Banda arc. Our results also show that lithosphere subducting from the north and south starts to collide at a depth of 300–350 km and becomes shallower further east. Apparent discontinuities in the high velocity band and a corresponding lack of seismicity supports the presence of a slab tear to the west of Seram. A dipping high velocity structure that is present from south to north beneath the island of Timor represents a subducting slab that dips more steeply beyond a depth of 150–200 km, which appears consistent with slab roll-back. Our tomographic model also shows evidence of back arc thrusting to the north of Sumbawa and Flores Islands in the form of a south-dipping higher velocity band at shallow depth. Furthermore, our tomographic models also reveal the possible presence of underthrust continental forearc in the form of a thin higher velocity anomaly that connects the backarc thrust and northward dipping lithosphere slab in the Timor area. Finally, a zone of low velocity above the higher velocity slab is clearly seen beneath Seram Island at a depth of ∼100 km and may represent a partial melting zone.

The Influence of Mantle Structure on Dynamic Topography in Southern Africa

AbstractDue 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.

Widespread D″ ${\mathbf{D}}^{\mathbf{{\prime\prime}}}$ Anisotropy Beneath North America and the Northeastern Pacific and Implications for Upper Mantle Anisotropy Measurements

AbstractObservations of seismic waves that have passed through the Earth's lowermost mantle provide insight into deep mantle structure and dynamics, often on relatively small spatial scales. Here we use SKS, S2KS, S3KS, and PKS signals recorded across a large region including the United States, Mexico, and Central America to study the deepest mantle beneath large swaths of North America and the northeastern Pacific Ocean. These phases are enhanced via beamforming and then used to investigate polarization‐ and propagation direction‐dependent shear wave speeds (seismic anisotropy). A differential splitting approach enables us to robustly identify contributions from anisotropy. Our results show strong seismic anisotropy in approximately half of our study region, indicating that anisotropy may be more prevalent than commonly thought. In some regions, the anisotropy may be induced by flow driven by sinking cold slabs, and in other, more compact regions, by upwelling flow. Measured splitting due to lowermost mantle anisotropy is sufficiently strong to be non‐negligible in interpretations of SKS splitting due to upper mantle anisotropy in certain regions, which may prompt future re‐evaluations of upper mantle anisotropy beneath North and Central America.

Seismic structure of Iceland revealed by ambient noise Rayleigh wave tomography

As it is an ideal location for studying plume–ridge interactions, a clear image of the Icelandic upper mantle structure is necessary. We collect continuous seismic records from 164 stations and extract Rayleigh wave dispersion curves via the frequency-Bessel (F-J) transform method. Based on ambient noise tomography, we provide a new shear-wave velocity model of the Icelandic crust and uppermost mantle, extending to a depth of 120 km. The model is validated by the waveform simulation method and reveals extensive crustal low-velocity zones (LVZs) across both the neovolcanic and nonvolcanic zones of Iceland. These crustal LVZs may be attributed to elevated temperatures, partial melting, and lithological variations. A distinct LVZ beneath a depth of 60 km, mainly on the North American Plate, may correspond to Icelandic plume material. Additionally, hot plume material may be delivered to the crust through low-velocity conduits beneath the spreading mid-ocean ridge. There is a clear contrast between the uppermost mantle low-velocity zones (UMLVZs) in the western region and the uppermost mantle high-velocity zones in the eastern region, which may indicate asymmetric tectonic plates on both sides of the mid-ocean ridge. This asymmetry may be attributed to the multiple eastward jumps of the ridge systems. The eastern high-velocity body, meaning a cooler uppermost mantle than that of the western region, may act as a barrier to obstruct the eastward plume flow. Under plume–ridge interactions, plume material can affect crustal accretion and feed volcanic activity on the surface along the spreading Mid-Atlantic Ridge.

Driven magmatism and crustal thinning of coastal southern China in response to subduction

Abstract. The late Mesozoic igneous rocks along the coastal South China Block (SCB) exhibit complex parental sources involving a depleted mantle, subducted sediment-derived melt, and melted crust. This period aligns with the magmatic flare-up and lull in the SCB, debating with the compression or extension in coastal region. Our study employs numerical models to investigate the dynamics of the ascent of underplating magma along the Changle–Nan'ao Belt (CNB), simulating its intrusion and cooling processes while disregarding the formational background. The rheological structure of the lithospheric mantle significantly influences magma pathways, dictating the distribution of magmatism. This work reveals that the ascent of magma in the presence of faults is considerably faster than in the absence of faults, and contemporaneous magmatic melts could produce different cooling and diagenetic processes. Additionally, the influence of pre-existing magma accelerated the emplacement of underplating magma. The magma beneath the fault ascended rapidly, reaching the lower crust within 20 million years, with a cooling rate of approximately ∼ 35 °C Myr−1. Conversely, the thickened magma took 40–50 million years to ascend to the lower crust, cooling at a rate of ∼ 10 °C Myr−1. In contrast, magma without thickening and fault would take a considerably longer time to reach the lower crust. The ascent of magma formed a mush-like head, contributing to magmatic circulation beneath the crust and decreasing crustal thickness. Multiphase magmatism increases the geothermal gradient, reducing lithospheric viscosity and promoting underplating magma ascent, leading to magmatic flare-ups and lulls. Our findings suggest that the Cretaceous magmatism at different times in the coastal SCB may be associated with the effects of lithospheric faults under similar subduction conditions. Boundary compression forces delay magma ascent, while rising magma induces a significant circulation, decreasing the crustal thickness of the coastal SCB. This study provides new insights into the complex interplay of magmatic processes during subduction, emphasizing the role of lithospheric structure in shaping the temporal and spatial evolution of coastal magmatism.

Upper mantle structure beneath the Mongolian region from multimode surface waves: Implications for the western margin of Amurian plate

Multimode phase speeds of surface waves are used to build a new radially anisotropic S wave model in the eastern Eurasian and Mongolian regions. Our dataset includes seismic waveforms of over 1655 teleseismic events (Mw≥5.8) from 2009 to 2021, recorded at permanent and temporary stations in and around Mongolia. The multimode dispersion curves of Love and Rayleigh waves were extracted using the nonlinear waveform fitting method for individual seismograms. Then, we retrieved phase speed maps for each mode and frequency, incorporating finite-frequency effects. Finally, localized multimode dispersion curves extracted from the phase speed maps were inverted for local 1-D SV and SH wave profiles, which are combined into a radially anisotropic 3-D shear wave model. Our new model exhibits significant lateral variations of S wave speeds at 70–100 km depth beneath Mongolia, i.e., slow anomalies in the tectonically active western Mongolia in contrast to fast anomalies in stable eastern Mongolia. In the radial anisotropy model, SH waves are faster than SV waves in most areas of the Mongolian lithosphere above 100 km depth, except for the northeast of the Altay Mountains. The Hangay Dome region is characterized by significantly slower velocities that may relate to its uplifting. A large-scale low velocity beneath the northeast of the Hangay Dome with a slower SV wave speed than SH may indicate the existence of partially molten layers. This study also reveals distinct lateral variations of S wave speeds across the boundary between the Amurian and Eurasian plates, characterized by the fast anomaly in eastern Mongolia, corresponding to the lithosphere in the western Amurian plate.

High-resolution 3-D lithospheric structure beneath the Qinling-Dabie orogenic belt from joint inversion of receiver functions and ambient noise

Resulting from the convergence of the Yangtze and North China Cratons, the Qinling-Dabie orogenic zone (QD) represents an important element in the central China orogenic system. To fully comprehend the craton evolution and lower crustal flow from the Tibetan Plateau, it is important to understand the crust and mantle structure of the QD. We reconstructed the three-dimensional lithospheric structure beneath the QD with high resolution using the joint inversion of receiver functions and ambient noise. Observations reveal that a high-velocity anomaly in the middle to lower crust beneath the western Qinling (WQL) orogen obstructs the eastward extension of a crustal low-velocity anomaly originating from the Tibetan Plateau. This finding provides unambiguous evidence that the WQL orogen is not crossed by eastward lower crustal flow from the Tibetan Plateau. The lithospheric mantle beneath the Weihe Rift and East Qinling orogen exhibits low-velocity characteristics, indicating that eastward asthenospheric flow from the Tibetan Plateau has caused substantial thermal-chemical erosion in the uppermost mantle beneath these regions. The results additionally indicate that the uppermost mantle high-velocity anomalies beneath the Dabie orogen is confined in a limited area and extends only to a depth of 70 km. We propose that during the Triassic, deeply subducted continental lithosphere returned into the uppermost mantle, forming the high-velocity anomalies beneath the Dabie orogen.

The Role of Upper Mantle Forces in Post‐Subduction Tectonics: Plumelet and Active Rifting in the East Anatolian Plateau

AbstractThe spatiotemporal interaction of large‐ and regional‐scale upper mantle forces can prevail in collisional settings. To better understand the role of these forces on post‐subduction tectonics, we focus on mantle dynamics in the East Anatolian Plateau, a well‐documented segment of the Arabian‐Eurasian continental collision zone. Specifically, we analyze multiple forces in the upper mantle, which have not been considered in previous studies in this region. To this end, we use a state‐of‐the‐art 3D instantaneous geodynamic model to quantify the dynamics of thermally defined upper mantle structures derived from seismic tomography data. Results reveal a prominent SW‐NE‐oriented mantle flow from the Arabian foreland to the Greater Caucasus–a plumelet–through a lithospheric channel under the East Anatolian Plateau. This plumelet induces localized dynamic topography (∼500 m) around the extensional Lake Van province, favoring NE‐directed compression and westward escape of the Anatolian plate. We suggest that the Lake Van region is an active magma‐rich intraplate rift in the Africa‐Arabia‐Anatolian plume‐rift system. The rift zone was probably initiated by Neotethyan subduction‐related forces and has been reactivated and/or sustained by the plumelet‐induced convective support. Our findings are consistent with numerous observations, including the recent low‐ultralow seismic velocities with a SW‐NE splitting anisotropy pattern, geochemical and petrological studies, and local kinematics showing upper mantle‐induced extensional tectonics in the collisional region.