Mantle Structure Research Articles

Source-encoded waveform inversion in the Northern Hemisphere

SUMMARY We use source-encoded waveform inversion to image Earth’s Northern Hemisphere. The encoding method is based on measurements of Laplace coefficients of stationary wavefields. By assigning to each event a unique frequency, we compute Fréchet derivatives for all events simultaneously based on one ‘super’ forward and one ‘super’ adjoint simulation for a small fraction of the computational cost of classical waveform inversion with the same data set. No cross-talk noise is introduced in the process, and the method does not require all events to be recorded by all stations. Starting from global model GLAD_M25, we performed 100 conjugate gradient iterations using a data set consisting of 786 earthquakes recorded by 9846 stations. Synthetic inversion tests show that we achieve good convergence based on this data set, and we see a consistent misfit reduction during the inversion. The new model, named SE100, has much higher spatial resolution than GLAD_M25, revealing details of the Yellowstone and Iceland hotspots, subduction beneath the Western United States and the upper mantle structure beneath the Arctic Ocean.

Integrated geophysical-petrological model of Damavand volcano, North Iran: Compositional structure of crust and upper mantle from seismic temperature profiling and gravity anomaly fitting

Our study employs well-constrained integrated geophysical-petrological modeling to investigate the compositional structure of Damavand volcano, the Alborz mountains in northern Iran. The thermochemical model utilizes high-resolution compressional velocity anomalies to determine temperature, while geochemical data from volcanic rocks further constrain mostly the lithospheric mantle model. By examining density and velocity, we deduce the upper mantle's compositional structure using observed free-air gravity anomalies and shear-wave seismic models. Our results reveal a lower crust density reduction of 150 kg/m3 due to a high-temperature regime (1080 °C) and the presence of melt indicated by a calculated Vp/Vs ≥ 1.8. The magma storage in the lithosphere is characterized by a harzburgite-melt mixture in agreement with shear-wave seismic profiles, and the lithospheric magma reservoir is positively buoyant with a density reduction of ∼75 kg/m3 compared to the ambient lithospheric mantle. The absence of a warm zone beneath the Alborz suggests that a transient mantle flow from the 660 km transition zone beneath Eurasia is responsible for the lateral heating (and possibly hydration) of the lithospheric mantle. This aligns with the high-temperature conditions observed in the Damavand (Alborz) region. Furthermore, the presence of elevated Nb (and high Nb/Yb ratio as well as radiogenic Nd isotope) concentrations in certain volcanic rocks suggests that a portion of the melt may have originated from the Earth's asthenosphere mantle. Moreover, the existence of water and/or volatiles (likely ≥0.64 wt%, as determined through a comparison of predicted and observed seismic velocities, along with previous laboratory analyses) promotes decompression melting of the lithospheric mantle at significantly lower temperatures (1230–1130 °C) in contrast to earlier estimations (1430–1180 °C). Finally, the reductions in velocity observed in the lower crust and uppermost lithospheric mantle beneath Damavand are likely attributed to thermal and compositional factors, respectively. The presence of a hydrated and melt-rich composition within the lithospheric mantle can account for the observed decrease in seismic velocity.

Studying inner core and lower mantle structure with a combination of PKP and converted SKP and PKS waves

SUMMARYStructure of the inner core is often measured through traveltime differences between waves that enter the inner core (PKPdf) and waves that travel through the outer core only (PKPab and PKPbc). Here we extend the method to converted waves PKSdf and SKPdf and compare results to PKP wave measurements. PKSdf and SKPdf have a very similar path to PKPdf and if velocity variations are present in the inner core, all three wave types should experience them equally. Since traveltime deviations can be due to velocity changes (either isotropic or anisotropy) as well as wave path deviations born from heterogeneity, we simultaneously investigate wave path directions and traveltimes of PKP, SKP and PKS waves for several source-array combinations. One of the path geometries is the anomalous polar corridor from South Sandwich to Alaska, which has strong traveltimes anomalies for PKPdf relative to more normal equatorial path geometries. Here we use array methods and determine slowness, traveltime and backazimuth deviations and compare them to synthetic data. We find that path deviations from theoretical values are present in all wave types and paths, but also that large scatter exists. Although some of the path deviations can be explained by mislocation vectors and crustal variations, our measurements argue that mantle structure has to be considered when investigating inner core anisotropy. Our polar path data show similar traveltime residuals as previously published, but we also find slowness residuals for this path. Interestingly, SKPdf and PKSdf for the South Sandwich to Alaska path show traveltime residuals that are partly opposite to those for PKPdf, possibly due to an interaction with a localized ultra-low velocity zone where waves enter the core.

Lunar mantle structure and composition inferred from Apollo 12 – Explorer 35 electromagnetic sounding

Constraints on the interior structure of the Moon have been derived from its inductive response, principally as measured by the magnetic transfer function (TF) between the distantly orbiting Explorer 35 satellite and the Apollo 12 surface station. The most successful prior studies used a dataset spanning 0.01–1 mHz, so the lunar response could be modeled as a simple dipole. However, earlier efforts also produced transfer functions up to 40 mHz. The smaller electromagnetic skin depth at higher frequency would better resolve the uppermost mantle—where key information about primitive lunar evolution may still be preserved—but requires a multipole treatment.I compute new profiles of electrical conductivity vs depth using both low- and high-frequency ranges of published Apollo-Explorer TFs. Using the low-frequency data, I derive temperature profiles at depths >400 km (<1 mHz) consistent with conductive heat loss and expectations of the iron (and possibly water) content of the mantle. The near-constant iron fraction (Mg# 81 ± 10) could imply efficient mixing due to now-defunct convection. Alternatively, incomplete overturn of gravitationally unstable magma-ocean cumulates could have left a heterogeneous distribution of minerals at hundred-km scales that are not resolved by electromagnetic sounding. A third explanation is that the electromagnetically probed region may be the initial equilibrium crystallization in a mantle that did not buoyantly overturn.In contrast, the high-frequency data produced higher conductivities than expected, requiring unrealistically low Mg# or high water content at depths 200–400 km. Either the published transfer functions > > 1 mHz are incorrect, or the TF multipole method at the Moon is unreliable. Future electromagnetic sounding using the magnetotelluric method can operate up to 100 s Hz and is largely insensitive to multipole effects, resolving structure to 100 km or less.

Internal Structure of the Lunar Mantle: Matching of Geochemical and Geophysical Models

There are serious contradictions between the geophysical and geochemical classes of models of the chemical composition and internal structure of the Moon, associated with the assessment of the abundance of the main oxides. The search for a potential consensus between the models was carried out on the basis of a set of geophysical and geochemical data using the Monte-Carlo method using the Markov chain scheme in combination with a method of minimization of the Gibbs free energy. The influence of the chemical composition and mineralogy of several conceptual models on the internal structure of the Moon has been studied. Two classes of chemical composition models are considered—the E models with terrestrial values of Al2O3 and CaO and M models with their higher content, as well as two classes of the most popular geochemical models, the Taylor Whole Moon (TWM) and Lunar Primitive Upper Mantle (LPUM) models, with ~45 wt % SiO2, but with different concentrations of refractory oxides and FeO. In both classes of E and M models, the lunar mantle is enriched in silica (~50 wt % SiO2) and FeO (11–13 wt %, Mg# 79–81) relative to the bulk composition of the silicate Earth (BSE, ~45 wt % SiO2, ~8 wt % FeO, Mg# 89). Such high concentrations of SiO2 and FeO become the determining factors for understanding the features of the mineral, velocity, and density structure of the lunar mantle. For the E and M models and geochemical models TWM and LPUM, the speed of sound and the density of stable phase associations are calculated. For E and M models, good agreement was obtained between the velocities of P- and S-waves and seismic sounding data from the Apollo program, which supports the idea of a silica-rich (olivine-pyroxenite) upper mantle. Unlike the Earth’s upper mantle, the dominant mineral in the Moon’s upper mantle is low-calcium orthopyroxene, not olivine. In contrast, the sound velocities of silica-unsaturated compositions, both FeO and Al2O3 enriched (TWM) and depleted (LPUM) models, do not match the seismic signatures. Thermodynamically justified restrictions on the chemical composition, mineralogy, and physical characteristics of the mantle based on the E and M models make it possible to eliminate some contradictions between the geochemical and geophysical classes of models of the internal structure of the Moon. Simultaneous enrichment in ferrous iron and silica is difficult to reconcile with the hypothesis of the formation of the Moon as a result of a giant impact from the substance of the Earth’s primitive mantle or from the substance of a shock body (bodies) of chondrite composition. Limitations on lunar concentrations of FeO and SiO2 probably correspond to the parent bodies of some achondrites.

Electrical properties of iron sulfide-bearing dunite under pressure: Effect of temperature, composition, and annealing time

Abstract The detection and quantification of metal sulfides in host rocks by electrical measurements have been priorities for field and laboratory studies, motivated by mineral prospecting and fundamental interest in the mantle structure or core/mantle differentiation, among other reasons. Here, we reanalyze electrical data for a dunite host with added FeS or Fe-S-Ni (Saxena et al. 2021), and report additional experimental runs along with electron microprobe analyses. The applied pressure was 2 GPa; impedance spectra were acquired while annealing at 1023 K (below the metal-sulfide solidus), and while varying temperature from 570 to 1650 K. Addition of 6.5 or 18 vol% FeS strongly enhances conductivity of the bulk sample compared with that of the dunite host, though values are 100–100 000 times less than those of pure FeS. These results indicate that most metal sulfide content is not part of a viable conductive path, even for the 18 vol% quantity. Nevertheless, the relatively high conductivity and weak temperature dependence of the 18 vol% sample reveal that contiguous paths of solid or molten FeS span the electrodes. The sample with 6.5 vol% sulfide also exceeds the percolation threshold for temperatures as low as ~100 K below the eutectic melting point, likely because FeS softens. Conductivity is nearly unchanged upon crossing the eutectic temperature, however a decline over 1400–1500 K reveals that the 6.5 vol% molten FeS forms a fragile electrical network in dunite. Samples with Fe50S40Ni10 or Fe40S40Ni20 (at%) are less conductive than pure dunite at temperatures below ~1450 K. This surprising result, likely caused by a reducing influence of Fe or Ni metal, does not support the use of FeS as an analog for compositions with nickel or excess metal. Our findings suggest that probing the electrical network of metal sulfides as solids complements other studies focused on connectivity of molten metal sulfides.

Silica-water superstructure and one-dimensional superionic conduit in Earth's mantle.

Water in Earth's deep interior is predicted to be hydroxyl (OH-) stored in nominally anhydrous minerals, profoundly modulating both structure and dynamics of Earth's mantle. Here, we use a high-dimensional neuro-network potential and machine learning algorithm to investigate the weight percent water incorporation in stishovite, a main constituent of the subducted oceanic crust. We found that stishovite and water prefer forming medium- to long-range ordered superstructures, featuring one-dimensional (1D) water channels. Synthesizing single crystals of hydrous stishovite, we verified the ordering of OH- groups in the water channels through optical and nuclear magnetic resonance spectroscopy and found an average H-H distance of 2.05(3) Å, confirming simulation results. Upon heating, H atoms were predicted to behave fluid-like inside the channels, leading to an exotic 1D superionic state. Water-bearing stishovite could feature high ionic mobility and strong electrical anisotropy, manifesting as electrical heterogeneity in Earth's mantle.

Mapping the Structure and Metasomatic Enrichment of the Lithospheric Mantle Beneath the Kimberley Craton, Western Australia

AbstractThe lithology, geochemistry, and architecture of the continental lithospheric mantle (CLM) underlying the Kimberley Craton of north‐western Australia has been constrained using pressure‐temperature estimates and mineral compositions for &gt;5,000 newly analyzed and published garnet and chrome (Cr) diopside mantle xenocrysts from 25 kimberlites and lamproites of Mesoproterozoic to Miocene age. Single‐grain Cr diopside paleogeotherms define lithospheric thicknesses of 200–250 km and fall along conductive geotherms corresponding to a surface heat flow of 37–40 mW/m2. Similar geotherms derived from Miocene and Mesoproterozoic intrusions indicate that the lithospheric architecture and thermal state of the CLM has remained stable since at least 1,000 Ma. The chemistry of xenocrysts defines a layered lithosphere with lithological and geochemical domains in the shallow (&lt;100 km) and deep (&gt;150 km) CLM, separated by a diopside‐depleted and seismically slow mid‐lithosphere discontinuity (100–150 km). The shallow CLM is comprised of Cr diopsides derived from depleted garnet‐poor and spinel‐bearing lherzolite that has been weakly metasomatized. This layer may represent an early (Meso to Neoarchean?) nucleus of the craton. The deep CLM is comprised of high Cr2O3 garnet lherzolite with lesser harzburgite, and eclogite. The peridotite components are inferred to have formed as residues of polybaric partial mantle melting in the Archean, whereas eclogite likely represents former oceanic crust accreted during Paleoproterozoic subduction. This deep CLM was metasomatized by H2O‐rich melts derived from subducted sediments and high‐temperature FeO‐TiO2 melts from the asthenosphere.

PREM-like velocity structure in the outermost core from global SKS and ScS waveform modeling

Seismic structure in the Earth's outermost core is important for our understanding of thermochemical stratification in the outer core, and is yet heavily debated. Here we study the compressional velocity structure in the Earth's outermost core based on waveform modeling of a unique untapped SKS and ScS dataset near bifurcation distances, collected from global seismic arrays for earthquakes occurring from 2000 to 2020. Using the SKS-ScS array dataset minimizes the effects of many uncertainties associated with earthquake source parameters and seismic heterogeneities in the mantle, and affords an opportunity to study and assess the seismic structure in the outermost core. We study outer core structure by testing two end-member models: 1) the D″ model that attributes any anomalous seismic observations to the effects of the lowermost mantle structure and is paired with a PREM (the Preliminary Reference Earth's Model) structure in the outermost core and 2) the outer core model that is paired with either a PREM or a tomographic structure in the lowermost mantle and attributes other unexplained seismic signals to an outer core structure. The results of the outer core models present unreasonable large lateral variations of >3.1% in the outermost core, while the inferred D″ models exhibit large-scale seismic anomalies that are consistent with the tomographic models and small-scale anomalies that are confirmed by further analysis of the seismic array data. Our analyses suggest a PREM-like seismic velocity structure and a lack of strong thermochemical anomalies in the topmost ∼200 km of the outer core, placing bounds on possible thermal and compositional conditions in the region of the Earth's outermost core. Our study also identifies the existence of small-scale seismic anomalies and sharp velocity variations in the lowermost mantle beneath the south coast of Alaska, northwestern Atlantic and the middle of Central America.

Fizzy Super-Earths: Impacts of Magma Composition on the Bulk Density and Structure of Lava Worlds

Lava worlds are a potential emerging population of Super-Earths that are on close-in orbits around their host stars, with likely partially molten mantles. To date, few studies have addressed the impact of magma on the observed properties of a planet. At ambient conditions, magma is less dense than solid rock; however, it is also more compressible with increasing pressure. Therefore, it is unclear how large-scale magma oceans affect planet observables, such as bulk density. We update ExoPlex, a thermodynamically self-consistent planet interior software, to include anhydrous, hydrous (2.2 wt% H2O), and carbonated magmas (5.2 wt% CO2). We find that Earth-like planets with magma oceans larger than ∼1.5 R ⊕ and ∼3.2 M ⊕ are modestly denser than an equivalent-mass solid planet. From our model, three classes of mantle structures emerge for magma ocean planets: (1) a mantle magma ocean, (2) a surface magma ocean, and (3) one consisting of a surface magma ocean, a solid rock layer, and a basal magma ocean. The class of planets in which a basal magma ocean is present may sequester dissolved volatiles on billion-year timescales, in which a 4 M ⊕ mass planet can trap more than 130 times the mass of water than in Earth’s present-day oceans and 1000 times the carbon in the Earth’s surface and crust.

The Influence of Lithospheric Thickness Variations Beneath Australia on Seismic Anisotropy and Mantle Flow

AbstractRapid plate motion, alongside pronounced variations in age and thickness of the Australian continental lithosphere, makes it an excellent location to assess the relationship between seismic anisotropy and lithosphere‐asthenosphere dynamics. In this study, SKS and PKS shear‐wave splitting is conducted for 176 stations covering the transition from the South Australian Craton to eastern Phanerozoic Australia. Comparisons are made with models of lithospheric thickness as well as numerical simulations of mantle flow. Splitting results show uniform ENE‐WSW aligned fast directions over the Gawler Craton and broader South Australian Craton, similar to the orientation of crustal structures generated during an episode of NW‐SE directed compression and volcanism ∼1.6 billion years ago. We propose that heat from volcanism weakened the lithosphere, aiding widespread lithospheric deformation, which has since been preserved in the form of frozen‐in anisotropy. Conversely, over eastern Phanerozoic Australia, fast directions show strong alignment with the NNE absolute plate motion. Overall, our results suggest that when the lithosphere is thin (&lt;125 km), lithospheric contributions are minimal and contributions from asthenospheric anisotropy dominate, reflecting shear of the underlying mantle by Australia's rapid plate motion above. Further insights from geodynamical simulations of the regional mantle flow field, which incorporate Australian and adjacent upper mantle structure, predict that asthenospheric material would be drawn in from the south and east toward the fast‐moving continental keel. Such a mechanism, alongside interactions between the flow field and lithospheric structure, provides a plausible explanation for smaller‐scale anomalous splitting patterns beneath eastern Australia that do not align with plate motion.

Improved lithospheric seismic velocity and density model of the Korean Peninsula from ambient seismic noise data using machine learning

To accurately estimate physical properties and tectonic behavior in the near-surface, a reliable seismic property model is needed for realistic seismic modeling and earthquake location estimation. Recording ambient seismic noise (ASN) data and producing interferometric reflection images traditionally provides subsurface structure observation without an active source. However, due to low signal-to-noise ratio (SNR) and vertical resolution, interpreting upper mantle structures and inverting seismic models from noise data is difficult. To address this, machine learning (ML) techniques are applied to enhance vertical resolution and interpret geologically meaningful boundaries. Spectral enhancement with convolutional U-net generates higher resolution data by preserving temporal continuity. Unsupervised ML interprets lithospheric boundaries more robustly and objectively than manual horizon picking, and model-based seismic inversion integrates improved seismic data with prior full-waveform inversion (FWI) models. ML-based results improve inverted models, displaying more detailed geological structures and seismic property changes, surpassing seismic data limitations.

The effect of the mantle and core matter phase state on the course of geodynamic processes

The study of the course of geodynamic processes in the lower crust and upper mantle proves that an additional energy contribution is made by a change in the phase state of matter with increasing pressure and temperature. The gas phase, composed of hydrogen, oxygen and carbon, turns into a fluid that combines the properties of a liquid and a gas. The result is a change in the behavior of fluid-crystal and fluid-melt systems which significantly accelerates melting and physicochemical interactions in the thermal asthenosphere. These conclusions are confirmed by numerous experimental studies and the results of the study of xenoliths representing the crust and mantle of cratons and active regions.&#x0D; Seismic tomography studies show distinct patterns of inhomogeneities in physical pro­per­ties, reflecting inhomogeneities in the mantle structure. Many works hypothesize, with substantiation, that plumes or fluid flows arise at the boundary of the core and mantle and are factors of all geodynamic processes. Modern ideas about the composition of the Earth's core are based on the statement that it is composed of molten iron with minor impurities of other elements. However, calculations of the energy balance and physical modeling of the redistribution of matter in the core itself show that the removal of volatile components or convective currents do not provide enough energy for the formation of plumes.&#x0D; The assumption that the substance of the core is an electrically conductive ionic liquid in which chemical compounds have completely dissociatedand the electronic structure has no gapradically changes the idea of the energetics of the core and the possibility of initiating plume processes. The properties of a substance in a similar phase state are fundamentally different from the properties of a liquid.

Seismic Velocity Structure of Upper Mantle Beneath the Oldest Pacific Seafloor: Insights From Finite‐Frequency Tomography

AbstractThe oldest oceanic basin (160–180 Ma) in the western Pacific is the birthplace of the Pacific Plate and is thus essential for understanding the formation and evolution of the oceanic plate. However, the upper mantle structure beneath the region has not been thoroughly investigated because of the remoteness and difficulties of long‐term in situ seismic measurements at the ocean bottom. From 2018 to 2019, the Oldest‐1 experiment on the oldest seafloor was conducted as part of the international Pacific Array initiative. We present the first three‐dimensional P‐wave velocity structure down to a depth of 350 km based on the relative travel time residuals of teleseismic earthquakes recorded by 11 broadband ocean‐bottom seismometers operated during the Oldest‐1 experiment. Our result shows a fast P‐wave velocity anomaly (VP perturbation of 2%–4% faster than average) at a depth of 95–185 km beneath the northeast of the study area. This structure is interpreted as evidence of dry, viscous, and rigid materials at depths below the lithosphere. Two slow anomalies (VP perturbation of 2%–4% slower than average) are seen beneath the southwestern and eastern (the oldest seafloor &gt;170 Ma) parts of the array site. The low‐velocity zones are found at depths of 95–305 km. The observed velocity structures can be indicative of plume activities that affected the upper mantle as the Pacific Plate migrated over hotspots from the southeast. Alternatively, the observed velocity features may provide seismic evidence for small‐scale sublithospheric convection.

Mantle Wavespeed and Discontinuity Structure Below East Africa: Implications for Cenozoic Hotspot Tectonism and the Development of the Turkana Depression

AbstractEthiopia's Cenozoic flood basalt magmatism, uplift, and rifting have been attributed to one or more mantle plumes. The Nubian plate, however, has drifted 500–1,000 km north since initial magmatism at ∼45 Ma, having developed above mantle that now underlies the northern Tanzania craton and the low‐lying Turkana Depression. Unfortunately, our knowledge of mantle wavespeed structure and mantle transition zone (MTZ) topography below these regions is poorest, due to a historical lack of seismograph stations. The same data gap means we lack constraints on lithospheric structure in and around the NW–SE trending Mesozoic Anza rift. We exploit data from new seismograph networks in the Turkana Depression and neighboring northern Uganda to develop AFRP22, a new African absolute P‐wavespeed tomographic model that resolves whole mantle structure along the entire East African rift system. We also map MTZ thickness using Ps receiver functions. East Africa's thinnest MTZ (∼25 km thinning) underlies the northwest Turkana Depression. AFRP22 reveals a co‐located, previously unrecognized, slow wavespeed plume tail, extending from the MTZ, deep into the lower mantle. This plume may thus have contributed, along with the African Superplume, to the development of the 45–30 Ma flood basalt province that preceded extension. Pervasive sub‐lithospheric slow wavespeeds imply that Turkana's present‐day low elevation is explained best by Mesozoic and Cenozoic‐age crustal thinning. At ∼100 km depth, AFRP22 illuminates a fast wavespeed SE Ethiopian plateau. In addition to governing the northernmost limit of Mesozoic Anza rifting, the refractory nature of this lithospheric block likely minimized Cenozoic flood basalt magmatism there.

New Insights into the Mesozoic Large‐scale Intra‐plate Magmatism and Mineralization in South China from Seismic Tomography

AbstractA great amount of new S‐wave data has been combined to image the mantle structure down to lower mantle depth near the southeastern margin of Eurasia and understand the tectonics in South China since the Mesozoic. Our results reveal a large‐scale structure beneath the South China Sea, with a prominent, broad low‐velocity feature of at least 1500‐km wide in and below the mantle transition zone (MTZ) and a pronounced low‐velocity feature of nearly 500‐km wide in the lowermost mantle. Together these features may represent the head‐and‐tail of an upwelling mantle plume that has interacted with and deformed by the subducting plates in the region. The Andaman–Sumatra plate is seen subducting steeply without stagnation in the MTZ over the large‐scale low‐velocity structure, suggesting that the upwelling mantle may have exerted significant influence on the overlying plate tectonics. Our results suggest a long‐term extensional regime near the southeastern margin of Eurasia, which may help to understand the formation of the large‐scale intra‐plate magmatism and mineralization with accumulation of huge amount of rare metal elements in South China during the Mesozoic, and the continental breakup, the southeastward extrusion of Indochina, and the opening of the South China Sea during the Cenozoic.

Plate Age and Uppermost Mantle Structure Across the Juan de Fuca and Gorda Plates

AbstractTo test the relationship between plate age and uppermost mantle structure for the Juan de Fuca and Gorda plates, as well as the mode of lithospheric cooling, we produce a new 3D oceanic shear wave speed (VS) model. This model is constructed using a new high‐quality Rayleigh wave phase speed data set across the plates and is based on a novel thermo‐seismic hybrid parameterization. In the lithosphere, we find that the VS structure is primarily (68% of study area) governed by conductive cooling. Exceptions are observed, suggesting advective heat influx along ridges and hydrothermal circulation through cracks near the subduction trench and transform faults. In the asthenosphere of the Juan de Fuca plate, no simple monotonic relationship between seismic structure and plate age exists, but instead there is a non‐monotonic slow‐fast‐slow VS pattern from west to east. The two slow asthenospheric VS anomalies indicate the presence of partial melt beneath both the mid‐ocean ridge and the subduction trench. The partial melt beneath the trench is enigmatic.

The Crustal and Uppermost Mantle Vs Structure of the Middle and Lower Reaches of the Yangtze River Metallogenic Belt: Implications for Metallogenic Process

AbstractThe Middle and Lower Reaches of the Yangtze River metallogenic belt (MLYMB) is one of the most important Fe‐Cu polymetallic belts in China. However, the mechanism and deep geodynamical process for the formation of this belt are still controversial. Here, we obtain the crustal and the uppermost mantle structures using ambient noise data from a dense seismic profile. A low velocity zone is revealed beneath the Moho of MLYMB, interpreted as the source of the deep mineralization materials. In addition, a low velocity layer (LVL) and a high velocity layer (HVL) are observed in the crust of the southern segment of the profile. The LVL is interpreted as a tectonic detachment layer between the upper and the lower crust, and the HVL is interpreted as the aggregation zone for mineralizing melts or crystallized magma chambers. Based on the observed velocity features, we propose a three‐stage model for the formation of ore deposits in MLYMB. Our model suggests that an upwelling of asthenosphere triggered by the delamination of a previously thickened lithosphere leads to the partial melting of upper mantle rocks, which eventually ponders under the Moho. The magma then infiltrates through the ductile lower crust and reaches a depth of ∼7–13 km, forming a minerals‐enriched magma chamber. Minerals‐rich hot fluids originating from the magma chamber continue to move upward along the pre‐existent faults and the minerals finally precipitate in dense veinlets when reaching shallow depths, forming the ore deposits in and around the MLYMB.

The Mantle Viscosity Structure of Venus

AbstractThe long‐wavelength gravity and topography of Venus are dominated by mantle convective flows, and are hence sensitive to the planet's viscosity structure and mantle density anomalies. By modeling the dynamic gravity and topography signatures and by making use of a Bayesian inference approach, we investigate the viscosity structure of the Venusian mantle by constraining radial viscosity variations. We performed inversions under a wide range of model assumptions that consistently predicted the existence of a thin low‐viscosity zone in the uppermost mantle. The zone is about 235 km thick and has a viscosity reduction of 5–15 times with respect to the underlying mantle. Drawing a parallel with the Earth, the reduced viscosity could be a result of partial melting as suggested for the origin of the asthenosphere. These results support the interpretation that Venus is a geologically active world predominantly governed by ongoing magmatic processes.

Crustal and uppermost mantle structures imaged by teleseismicP-wave traveltime tomography beneath the Southeastern Korean Peninsula: implications for a hydrothermal system controlled by the thermally modified lithosphere

SUMMARYThe southeastern Korean Peninsula (SeKP) has experienced intense deformation owing to subduction and backarc extension at the eastern continental margin of the Eurasian Plate, leading to the formation of complex tectonic structures. Abnormally high surface heat flux, Cenozoic volcanism, signatures of mantle degassing and hydrothermal alteration, and several active fault systems with extensional sedimentary basins have been identified; however, the major driving forces that promote local seismic events and hydrothermal activities remain enigmatic. Here, we constructed 3-D P-wave velocity of the crust and upper mantle in the SeKP for the first time using a teleseismic traveltime tomography method and an extensive data set obtained from a dense seismic network. Our model revealed three distinct velocity patterns at different depths: (1) in the upper crust (depth ∼0–10 km), a low-velocity anomaly beneath the Cenozoic sedimentary basin exhibiting a prominent lateral velocity contrasts with higher velocities in the Cretaceous sedimentary and plutonic rocks; (2) a N–S trending low-velocity anomaly extending from the lower crust to the uppermost mantle (depth ∼20–35 km) beneath the major active fault systems interpreted as a thermally or mechanically weakened structure that could transfer high surface heat flux and transport mantle-driven gases and (3) a low-velocity anomaly adjacent to the Cenozoic basin in the upper mantle at depths of 35–55 km interpreted as the higher temperature upper mantle. Via a series of geodynamic simulations, we demonstrated that the extensional deformation at the eastern continental margin during the Early to Middle Miocene locally enhanced the temperature of the crust and upper mantle beneath the SeKP. We propose that a hydrothermal system, resulting from the thermally modified lithosphere of the continental margin, has contributed to the enhanced local seismicity and geothermal activities observed in the SeKP region.