Full-waveform inversion reveals diverse origins of lower mantle positive wave speed anomalies

Determining Earth's structure is paramount to unravel its interior dynamics. Seismic tomography reveals positive wave speed anomalies throughout the mantle that spatially correlate with the expected locations of subducted slabs. This correlation has been widely applied in plate reconstructions and geodynamic modelling. However, global travel-time tomography typically incorporates only a limited number of easily identifiable body wave phases and is therefore strongly dependent on the source-receiver geometry. Here, we show how global full-waveform inversion is less sensitive to source-receiver geometry and reveals numerous previously undetected positive wave speed anomalies in the lower mantle. Many of these previously undetected anomalies are situated below major oceans and continental interiors, with no geologic record of subduction, such as beneath the western Pacific Ocean. Moreover, we find no statistically significant correlation positive anomalies as imaged using full-waveform inversion and past subduction. These findings suggest more diverse origins for these anomalies in Earth's lower mantle, unlocking full-waveform inversion as an indispensable tool for mantle exploration.

Mesozoic intraoceanic subduction shaped the lower mantle beneath the East Pacific Rise

The Pacific large low-shear-velocity province (LLSVP), as revealed by cluster analysis of global tomographic models, hosts multiple internal anomalies, including a notable gap (~20° wide) between the central and eastern Pacific. The cause of the structural gap remains unconstrained. Directly above this structural gap, we identify an anomalously thick mantle transition zone east of the East Pacific Rise, the fastest-spreading ocean ridge in the world, using a dense set of SS precursors. The area of the thickened transition zone exhibits faster-than-average velocities according to recent tomographic images, suggesting perturbed postolivine phase boundaries shifting in response to lowered temperatures. We attribute this observation to episodes of Mesozoic-aged (250 to 120 million years ago) intraoceanic subduction beneath the present-day Nazca Plate. The eastern portion of the Pacific LLSVP was separated by downwelling because of this ancient oceanic slab. Our discovery provides a unique perspective on linking deep Earth structures with surface subduction.

Hydrous mantle plume promoted the generation of continental flood basalts in the Tarim large igneous province

Recent research on the water content of large igneous provinces (LIPs) has revealed that water has a significant impact on the formation of LIPs. However, most studies focus on the water content of mafic-ultramafic rocks, while relatively little attention has been paid to the water content of continental flood basalts (CFB), which form the major part of LIPs and are characterized by huge volumes (> 1 × 105 km3) and short eruption times. Here, we determined water contents of clinopyroxene crystals from the Akesu diabase, which is co-genetic with flood basalts of the Tarim LIP in China. Based on these measurements, we obtained a water content of higher than 1.23 ± 0.49 wt.% for the parental magma to the Tarim CFB and a minimum water content of 1230 ± 490 ppm for the mantle source, thus indicating the presence of a hydrous mantle plume. Combined with previous studies, our results suggest that water plays a key role in the formation of the Tarim LIP. Additionally, the whole-rock compositions of the Akesu diabase indicate a contribution of pyroxenite in the mantle source. This is consistent with a model, in which water was brought into the Tarim mantle plume by a subducted oceanic plate that entered the deep mantle.

Compositional and thermal state of the lower mantle from joint 3D inversion with seismic tomography and mineral elasticity

The compositional and thermal state of Earth's mantle provides critical constraints on the origin, evolution, and dynamics of Earth. However, the chemical composition and thermal structure of the lower mantle are still poorly understood. Particularly, the nature and origin of the two large low-shear-velocity provinces (LLSVPs) in the lowermost mantle observed from seismological studies are still debated. In this study, we inverted for the 3D chemical composition and thermal state of the lower mantle based on seismic tomography and mineral elasticity data by employing a Markov chain Monte Carlo framework. The results show a silica-enriched lower mantle with a Mg/Si ratio less than ~1.16, lower than that of the pyrolitic upper mantle (Mg/Si = 1.3). The lateral temperature distributions can be described by a Gaussian distribution with a standard deviation (SD) of 120 to 140 K at 800 to 1,600 km and the SD increases to 250 K at 2,200 km depth. However, the lateral distribution in the lowermost mantle does not follow the Gaussian distribution. We found that the velocity heterogeneities in the upper lower mantle mainly result from thermal anomalies, while those in the lowermost mantle mainly result from compositional or phase variations. The LLSVPs have higher density at the base and lower density above the depth of ~2,700 km than the ambient mantle, respectively. The LLSVPs are found to have ~500 K higher temperature, higher Bridgmanite and iron content than the ambient mantle, supporting the hypothesis that the LLSVPs may originate from an ancient basal magma ocean formed in Earth's early history.

Ghost-arc geochemical anomaly at a spreading ridge caused by supersized flat subduction

The Southern Atlantic-Southwest Indian ridges (SASWIR) host mid-ocean ridge basalts with a residual subduction-related geochemical fingerprint (i.e., a ghost-arc signature) of unclear origin. Here, we show through an analysis of plate kinematic reconstructions and seismic tomography models that the SASWIR subduction-modified mantle source formed in the Jurassic close to the Georgia Islands slab (GI) and remained near-stationary in the mantle reference frame. In this analysis, the GI lies far inboard the Jurassic Patagonian-Antarctic Peninsula active margin. This was formerly attributed to a large-scale flat subduction event in the Late Triassic-Early Jurassic. We propose that during this flat slab stage, the subduction-modified mantle areas beneath the Mesozoic active margin and surrounding sutures zones may have been bulldozed inland by >2280 km. After the demise of the flat slab, this mantle anomaly remained near-stationary and was sampled by the Karoo mantle plume 183 Million years (Myr) ago and again since 55 Myr ago by the SASWIR. We refer to this process as asthenospheric anomaly telescoping. This study provides a hitherto unrecognized geodynamic effect of flat subduction, the viability of which we support through numerical modeling.

Exploratory Volumetric Deep Earth Visualization by 2.5D Interactive Compositing

In this contribution we consider the visualization of global, deep Earth volume datasets for display and researcher interaction. While the algorithms and data analysis techniques that produce such volumetric results have become more sophisticated, the manner of visualizing these findings can be improved. We address the challenge of making an illustrative, exploratory visualization of a global geoscience dataset using a combined seismic tomography result, the primary means by which geoscientists infer structure and process in the deep Earth. We present a novel, interactive graphical application suite and associated workflow that uses an intuitive 2.5D layer compositing approach. This allows the user to adjust the separation between data-slices, control graphics variables such as color mapping, opacity and compositing, and facilitate exploration and annotation of the architecture of the lithosphere. Graphics outputs from our applications are enabled for immersive systems such as dome displays. In a case study we visualize the deep Earth structure beneath the Indian Ocean region. We anticipate that the application methodology will find use in the visualization of multiple datasets representing aspects of the Earth's deep interior and atmosphere, and in the interaction with the increasing number of rich datasets from missions to our neighboring planets.

Quantification of Small-Scale Heterogeneity at the Core–Mantle Boundary Using Sample Entropy of SKS and SPdKS Synthetic Waveforms

Qualitative and quantitative analysis of seismic waveforms sensitive to the core–mantle boundary (CMB) region reveal the presence of ultralow-velocity zones (ULVZs) that have a strong decrease in compressional (P) and shear (S) wave velocity, and an increase in density within thin structures. However, understanding their physical origin and relation to the other large-scale structures in the lowermost mantle are limited due to an incomplete mapping of ULVZs at the CMB. The SKS and SPdKS seismic waveforms is routinely used to infer ULVZ presence, but has thus far only been used in a limited epicentral distance range. As the SKS/SPdKS wavefield interacts with a ULVZ it generates additional seismic arrivals, thus increasing the complexity of the recorded wavefield. Here, we explore utilization of the multi-scale sample entropy method to search for ULVZ structures. We investigate the feasibility of this approach through analysis of synthetic seismograms computed for PREM, 1-, 2.5-, and 3-D ULVZs as well as heterogeneous structures with a strong increase in velocity in the lowermost mantle in 1- and 2.5-D. We find that the sample entropy technique may be useful across a wide range of epicentral distances from 100° to 130°. Such an analysis, when applied to real waveforms, could provide coverage of roughly 85% by surface area of the CMB.

Mapping global kimberlite potential from reconstructions of mantle flow over the past billion years

Kimberlites are the primary source of economic grade diamonds. Their geologically rapid eruptions preferentially occur near or through thick and ancient continental lithosphere. Studies combining tomographic models with tectonic reconstructions and kimberlite emplacement ages and locations have revealed spatial correlations between large low shear velocity provinces in the lowermost mantle and reconstructed global kimberlite eruption locations over the last 320 Myr. These spatial correlations assume that the lowermost mantle structure has not changed over time, which is at odds with mantle flow models that show basal thermochemical structures to be mobile features shaped by cold sinking oceanic lithosphere. Here we investigate the match to the global kimberlite record of stationary seismically slow basal mantle structures (as imaged through tomographic modelling) and mobile hot basal structures (as predicted by reconstructions of mantle flow over the past billion years). We refer to these structures as “basal mantle structures” and consider their intersection with reconstructed thick or ancient lithosphere to represent areas with a high potential for past eruptions of kimberlites, and therefore areas of potential interest for diamond exploration. We use the distance between reconstructed kimberlite eruption locations and kimberlite potential maps as an indicator of model success, and we find that mobile lowermost mantle structures are as close to reconstructed kimberlites as stationary ones. Additionally, we find that mobile lowermost mantle structures better fit major kimberlitic events, such as the South African kimberlite bloom around 100 Ma. Mobile basal structures are therefore consistent with both solid Earth dynamics and with the kimberlite record.

The Accuracy Assessment of the PREM and AK135-F Radial Density Models

The Earth’s synthetic density and gravitational models can be used to validate numerical methods for global (or large-scale) gravimetric forward and inverse modelling formulated either in the spatial or spectral domains. The Preliminary Reference Earth Model (PREM) density parameters can be adopted as a 1-D reference density model and further refined using more detailed 2-D or 3-D crust and mantle density models. Alternatively, the AK135-F density parameters can be used for this purpose. In this study, we investigate options for a refinement of the Earth’s synthetic density model by assessing the accuracy of available 1-D density models, specifically the PREM and AK135-F radial density parameters. First, we use density parameters from both models to estimate the Earth’s total mass and compare these estimates with published results. We then estimate the Earth’s gravity field parameters, particularly the geoidal geopotential number W₀ and the mean gravitational attraction and compare them with published values. According to our results, the Earth’s total mass from the two models (the PREM and the AK135-F) differ less than 0.02% and 0.01%, respectively, when compared with the value adopted by the International Astronomical Union (IAU). The geoidal geopotential values of the two models differ from the value adopted by the IAU by less than 0.1% and 0.04%, respectively. The values of the mean gravitational attraction of the two models differ less than 0.02% and 0.08%, respectively, when compared with the value obtained from the geocentric gravitational constant and the Earth’s mean radius. These numerical findings ascertain that the PREM and AK135-F density parameters are suitable for defining a 1-D reference density model.

Assembly of the basal mantle structure beneath Africa

Plate tectonics shapes Earth's surface, and is linked to motions within its deep interior^1,2. Cold oceanic lithosphere sinks into the mantle, and hot mantle plumes rise from the deep Earth, leading to volcanism^3,4. Volcanic eruptions over the past 320 million years have been linked to two large structures at the base of the mantle presently under Africa and the Pacific Ocean^5,6. This has led to the hypothesis that these basal mantle structures have been stationary over geological time^7,8, in contrast to observations and models suggesting that tectonic plates^9,10, subduction zones^11-14 and mantle plumes^15,16 have been mobile, and that basal mantle structures are presently deforming^17,18. Here we reconstruct mantle flow from one billion years ago to the present day to show that the history of volcanism is statistically as consistent with mobile basal mantle structures as with fixed ones. In our reconstructions, cold lithosphere sank deep into the African hemisphere between 740 and 500 million years ago, and from 400 million years ago the structure beneath Africa progressively assembled, pushed by peri-Gondwana slabs, to become a coherent structure as recently as 60 million years ago. Our mantle flow models suggest that basal mantle structures are mobile, and aggregate and disperse over time, similarly to continents at Earth's surface⁹. Our models also predict the presence of continental material in the mantle beneath Africa, consistent with geochemical data^19,20.

Slab control on the mega-sized North Pacific ultra-low velocity zone

Ultra-low velocity zones (ULVZs) are localized small-scale patches with extreme physical properties at the core-mantle boundary that often gather at the margins of Large Low Velocity Provinces (LLVPs). Recent studies have discovered several mega-sized ULVZs with a lateral dimension of ~900 km. However, the detailed structures and physical properties of these ULVZs and their relationship to LLVP edges are not well constrained and their formation mechanisms are poorly understood. Here, we break the degeneracy between the size and velocity perturbation of a ULVZ using two orthogonal seismic ray paths, and thereby discover a mega-sized ULVZ at the northern edge of the Pacific LLVP. The ULVZ is almost double the size of a previously imaged ULVZ in this region, but with half of the shear velocity reduction. This mega-sized ULVZ has accumulated due to stable mantle flow converging at the LLVP edge driven by slab-debris in the lower mantle. Such flow also develops the subvertical north-tilting edge of the Pacific LLVP.

Catastrophic slab loss in southwestern Pangea preserved in the mantle and igneous record

The Choiyoi Magmatic Province represents a major episode of silicic magmatism in southwestern Pangea in the mid-Permian-Triassic, the origin of which remains intensely debated. Here, we integrate plate-kinematic reconstructions and the lower mantle slab record beneath southwestern Pangea that provide clues on late Paleozoic-Mesozoic subducting slab configurations. Also, we compile geochronological information and analyze geochemical data using tectono-magmatic discrimination diagrams. We demonstrate that this magmatic event resulted from a large-scale slab loss. This is supported by a paleogeographic coincidence between a reconstructed 2,800-3,000-km-wide slab gap and the Choiyoi Magmatic Province and geochemical data indicating a slab break-off fingerprint in the latter. The slab break-off event is compatible with Permian paleogeographic modifications in southwestern Pangea. These findings render the Choiyoi Magmatic Province the oldest example of a geophysically constrained slab loss event and open new avenues to assess the geodynamic setting of silicic large igneous provinces back to the late Paleozoic.

The origin of the Permian-Triassic Choiyoi silicic large igneous province (SLIP) is assessed by linking the igneous record, plate-kinematic reconstructions, and the deep mantle. This study suggests an origin related to a massive slab loss in Pangea.

On the relative temperatures of Earth's volcanic hotspots and mid-ocean ridges

Volcanic hotspots are thought to be fed by hot, active upwellings from the deep mantle, with excess temperatures (T_ex) ~100° to 300°C higher than those of mid-ocean ridges. However, T_ex estimates are limited in geographical coverage and often inconsistent for individual hotspots. We infer the temperature of oceanic hotspots and ridges simultaneously by converting seismic velocity to temperature. We show that while ~45% of plume-fed hotspots are hot (T_ex ≥ 155°C), ~15% are cold (T_ex ≤ 36°C) and ~40% are not hot enough to actively upwell (50°C ≤ T_ex ≤ 136°C). Hot hotspots have an extremely high helium-3/helium-4 ratio and buoyancy flux, but cold hotspots do not. The latter may originate at upper mantle depths. Alternatively, the deep plumes that feed them may be entrained and cooled by small-scale convection.

Seismological expression of the iron spin crossover in ferropericlase in the Earth's lower mantle

The two most abundant minerals in the Earth’s lower mantle are bridgmanite and ferropericlase. The bulk modulus of ferropericlase (Fp) softens as iron d-electrons transition from a high-spin to low-spin state, affecting the seismic compressional velocity but not the shear velocity. Here, we identify a seismological expression of the iron spin crossover in fast regions associated with cold Fp-rich subducted oceanic lithosphere: the relative abundance of fast velocities in P- and S-wave tomography models diverges in the ~1,400-2,000 km depth range. This is consistent with a reduced temperature sensitivity of P-waves throughout the iron spin crossover. A similar signal is also found in seismically slow regions below ~1,800 km, consistent with broadening and deepening of the crossover at higher temperatures. The corresponding inflection in P-wave velocity is not yet observed in 1-D seismic profiles, suggesting that the lower mantle is composed of non-uniformly distributed thermochemical heterogeneities which dampen the global signature of the Fp spin crossover.

This study identifies the predicted seismic expression of the high-to-low iron spin crossover in the deep Earth mineral ferropericlase. A depth-dependent signal is detected in the fastest and slowest regions, related to lateral temperature variations, of several global seismic tomography models.

Review of the Heat Flow Mapping in Polish Sedimentary Basin across Different Tectonic Terrains

Heat flow patterns variability related to the age of the consolidated, and differences in, sedimentary thickness of the sedimentary succession are important constraints upon the thermal state of the sedimentary fill and its geothermal energy potential. Heat flow in the Permian basin of central Europe varies from a low of 40 mWm−2 in the Precambrian Platform to 80 mWm−2 in the Paleozoic basement platform influencing temperature for geothermal potential drilling depth. Continuity of thermal patterns and compatibility of heat flow Q across the Permian basin across the Polish–German basin was known from heat flow data ever since the first heat flow map of Europe in 1979. Both Polish and German heat flow determinations used lab-measured thermal conductivity on cores. This is not the case for the recent heat flow map of Poland published in 2009 widely referenced in Polish geological literature. Significant differences in heat flow magnitude exist between many historical heat flow maps of Poland over the 1970s–1990s and recent 21st century patterns. We find that the differences in heat flow values of some 20–30 mWm−2 in Western Poland exist between heat flow maps using thermal conductivity models using well log interpreted mineral and porosity content and assigned world averages of rock and fluid thermal conductivity versus those measured on cores. These differences in heat flow are discussed in the context of resulting mantle heat flow and the Lithosphere-Asthenosphere Boundary depth modelled differences and possible overestimates of deep thermal conditions for enhanced geothermal energy prospects in Poland.

Paired EMI-HIMU hotspots in the South Atlantic-Starting plume heads trigger compositionally distinct secondary plumes?

Age-progressive volcanism is generally accepted as the surface expression of deep-rooted mantle plumes, which are enigmatically linked with the African and Pacific large low-shear velocity provinces (LLSVPs). We present geochemical and geochronological data collected from the oldest portions of the age-progressive enriched mantle one (EMI)-type Tristan-Gough track. They are part of a 30- to 40-million year younger age-progressive hotspot track with St. Helena HIMU (high time-integrated 238U/204Pb) composition, which is also observed at the EMI-type Shona hotspot track in the southernmost Atlantic. Whereas the primary EMI-type hotspots overlie the margin of the African LLSVP, the HIMU-type hotspots are located above a central portion of the African LLSVP, reflecting a large-scale geochemical zonation. We propose that extraction of large volumes of EMI-type mantle from the margin of the LLSVP by primary plume heads triggered upwelling of HIMU material from a more internal domain of the LLSVP, forming secondary plumes.

New Candidate Ultralow-Velocity Zone Locations from Highly Anomalous SPdKS Waveforms

Ultralow-velocity zones (ULVZs) at the core–mantle boundary (CMB) represent some of the most preternatural features in Earth’s mantle. These zones most likely contain partial melt, extremely high iron content ferropericlase, or combinations of both. We analyzed a new collection of 58,155 carefully processed and quality-controlled broadband recordings of the seismic phase SPdKS in the epicentral distance range from 106° to 115°. These data sample 56.9% of the CMB by surface area. From these recordings we searched for the most anomalous seismic waveforms that are indicative of ULVZ presence. We used a Bayesian approach to identify the regions of the CMB that have the highest probability of containing ULVZs, thereby identifying sixteen regions of interest. Of these regions, we corroborate well-known ULVZ existence beneath the South China Sea, southwest Pacific, the Samoa hotspot, the southwestern US/northern Mexico, and Iceland. We find good evidence for new ULVZs beneath North Africa, East Asia, and north of Papua New Guinea. We provide further evidence for ULVZs in regions where some evidence has been hinted at before beneath the Philippine Sea, the Pacific Northwest, and the Amazon Basin. Additional evidence is shown for potential ULVZs at the base of the Caroline, San Felix and Galapagos hotspots.

Two deep-mantle sources for Paleocene doming and volcanism in the North Atlantic

The North Atlantic Igneous Province (NAIP) erupted in two major pulses that coincide with the continental breakup and the opening of the North Atlantic Ocean over a period from 62 to 54 Ma. The unknown mantle structure under the North Atlantic during the Paleocene represents a major missing link in deciphering the geodynamic causes of this event. To address this outstanding challenge, we use a back-and-forth iterative method for time-reversed global convection modeling over the Cenozoic Era which incorporates models of present-day tomography-based mantle heterogeneity. We find that the Paleocene mantle under the North Atlantic is characterized by two major low-density plumes in the lower mantle: one beneath Greenland and another beneath the Azores. These strong lower-mantle upwellings generate small-scale hot upwellings and cold downwellings in the upper mantle. The upwellings are dispersed sources of magmatism and topographic uplift that were active on the rifted margins of the North Atlantic during the formation of the NAIP. While most studies of the Paleocene evolution of the North Atlantic have focused on the proto-Icelandic plume, our Cenozoic reconstructions reveal the equally important dynamics of a hot, buoyant, mantle-wide upwelling below the Azores.

Inferring Earth's discontinuous chemical layering from the 660-kilometer boundary topography

Topography, or depth variation, of certain interfaces in the solid Earth can provide important insights into the dynamics of our planet interior. Although the intermediate- and long-range topographic variation of the 660-kilometer boundary between Earth's upper and lower mantle is well studied, small-scale measurements are far more challenging. We found a surprising amount of topography at short length scale along the 660-kilometer boundary in certain regions using scattered P'P' seismic waves. Our observations required chemical layering in regions with high short-scale roughness. By contrast, we did not see such small-scale topography along the 410-kilometer boundary in the upper mantle. Our findings support the concept of partially blocked or imperfect circulation between the upper and lower mantle.

Southward propagation of Nazca subduction along the Andes

The Andean margin is the plate-tectonic paradigm for long-lived, continuous subduction, yet its geology since the late Mesozoic era (the past 100 million years or so) has been far from steady state. The episodic deformation and magmatism have been attributed to cyclic changes in the dip angle of the subducting slab, slab break-off and the penetration of the slab into the lower mantle; the role of plate tectonics remains unclear, owing to the extensive subduction of the Nazca-Farallon plate (which has resulted in more than 5,500 kilometres of lithosphere being lost to the mantle). Here, using tomographic data, we recreate the plate-tectonic geometry of the subducted Nazca slab, which enables us to reconstruct Andean plate tectonics since the late Mesozoic. Our model suggests that the current phase of Nazca subduction began at the northern Andes (5° S) during the late Cretaceous period (around 80 million years ago) and propagated southwards, reaching the southern Andes (40° S) by the early Cenozoic era (around 55 million year ago). Thus, contrary to the current paradigm, Nazca subduction has not been fully continuous since the Mesozoic but instead included episodic divergent phases. In addition, we find that foredeep sedimentation and the initiation of Andean compression are both linked to interactions between the Nazca slab and the lower mantle, consistent with previous modelling.

SubMachine: Web-Based Tools for Exploring Seismic Tomography and Other Models of Earth's Deep Interior

We present SubMachine, a collection of web-based tools for the interactive visualization, analysis, and quantitative comparison of global-scale data sets of the Earth's interior. SubMachine focuses on making regional and global-scale seismic tomography models easily accessible to the wider solid Earth community, in order to facilitate collaborative exploration. We have written software tools to visualize and explore over 30 tomography models-individually, side-by-side, or through statistical and averaging tools. SubMachine also serves various nontomographic data sets that are pertinent to the interpretation of mantle structure and complement the tomographies. These include plate reconstruction models, normal mode observations, global crustal structure, shear wave splitting, as well as geoid, marine gravity, vertical gravity gradients, and global topography in adjustable degrees of spherical harmonic resolution. By providing repository infrastructure, SubMachine encourages and supports community contributions via submission of data sets or feedback on the implemented toolkits.

Valence and spin states of iron are invisible in Earth's lower mantle

Heterogeneity in Earth's mantle is a record of chemical and dynamic processes over Earth's history. The geophysical signatures of heterogeneity can only be interpreted with quantitative constraints on effects of major elements such as iron on physical properties including density, compressibility, and electrical conductivity. However, deconvolution of the effects of multiple valence and spin states of iron in bridgmanite (Bdg), the most abundant mineral in the lower mantle, has been challenging. Here we show through a study of a ferric-iron-only (Mg_0.46Fe³⁺_0.53)(Si_0.49Fe³⁺_0.51)O₃ Bdg that Fe³⁺ in the octahedral site undergoes a spin transition between 43 and 53 GPa at 300 K. The resolved effects of the spin transition on density, bulk sound velocity, and electrical conductivity are smaller than previous estimations, consistent with the smooth depth profiles from geophysical observations. For likely mantle compositions, the valence state of iron has minor effects on density and sound velocities relative to major cation composition.

Tidal tomography constrains Earth's deep-mantle buoyancy

Earth's body tide-also known as the solid Earth tide, the displacement of the solid Earth's surface caused by gravitational forces from the Moon and the Sun-is sensitive to the density of the two Large Low Shear Velocity Provinces (LLSVPs) beneath Africa and the Pacific. These massive regions extend approximately 1,000 kilometres upward from the base of the mantle and their buoyancy remains actively debated within the geophysical community. Here we use tidal tomography to constrain Earth's deep-mantle buoyancy derived from Global Positioning System (GPS)-based measurements of semi-diurnal body tide deformation. Using a probabilistic approach, we show that across the bottom two-thirds of the two LLSVPs the mean density is about 0.5 per cent higher than the average mantle density across this depth range (that is, its mean buoyancy is minus 0.5 per cent), although this anomaly may be concentrated towards the very base of the mantle. We conclude that the buoyancy of these structures is dominated by the enrichment of high-density chemical components, probably related to subducted oceanic plates or primordial material associated with Earth's formation. Because the dynamics of the mantle is driven by density variations, our result has important dynamical implications for the stability of the LLSVPs and the long-term evolution of the Earth system.

On the relative motions of long-lived Pacific mantle plumes

Mantle plumes upwelling beneath moving tectonic plates generate age-progressive chains of volcanos (hotspot chains) used to reconstruct plate motion. However, these hotspots appear to move relative to each other, implying that plumes are not laterally fixed. The lack of age constraints on long-lived, coeval hotspot chains hinders attempts to reconstruct plate motion and quantify relative plume motions. Here we provide ⁴⁰Ar/³⁹Ar ages for a newly identified long-lived mantle plume, which formed the Rurutu hotspot chain. By comparing the inter-hotspot distances between three Pacific hotspots, we show that Hawaii is unique in its strong, rapid southward motion from 60 to 50 Myrs ago, consistent with paleomagnetic observations. Conversely, the Rurutu and Louisville chains show little motion. Current geodynamic plume motion models can reproduce the first-order motions for these plumes, but only when each plume is rooted in the lowermost mantle.

Using mantle plumes to reconstruct past plate motion is complicated, because plumes may not be fixed. Here, the authors demonstrate using ⁴⁰Ar/³⁹Ar ages that the Rurutu plume is relatively stable compared to the rapidly moving Hawaiian plume, yet it has a similar deep mantle origin.

Intraoceanic subduction spanned the Pacific in the Late Cretaceous-Paleocene

The notorious ~60° bend separating the Hawaiian and Emperor chains marked a prominent change in the motion of the Pacific plate at ~47 Ma (million years ago), but the origin of that change remains an outstanding controversy that bears on the nature of major plate reorganizations. Lesser known but equally significant is a conundrum posed by the pre-bend (~80 to 47 Ma) motion of the Pacific plate, which, according to conventional plate models, was directed toward a fast-spreading ridge, in contradiction to tectonic forcing expectations. Using constraints provided by seismic tomography, paleomagnetism, and continental margin geology, we demonstrate that two intraoceanic subduction zones spanned the width of the North Pacific Ocean in Late Cretaceous through Paleocene time, and we present a simple plate tectonic model that explains how those intraoceanic subduction zones shaped the ~80 to 47 Ma kinematic history of the Pacific realm and drove a major plate reorganization.

On the consistency of seismically imaged lower mantle slabs

The geoscience community is increasingly utilizing seismic tomography to interpret mantle heterogeneity and its links to past tectonic and geodynamic processes. To assess the robustness and distribution of positive seismic anomalies, inferred as subducted slabs, we create a set of vote maps for the lower mantle with 14 global P-wave or S-wave tomography models. Based on a depth-dependent threshold metric, an average of 20% of any given tomography model depth is identified as a potential slab. However, upon combining the 14 models, the most consistent positive wavespeed features are identified by an increasing vote count. An overall peak in the most robust anomalies is found between 1000-1400 km depth, followed by a decline to a minimum around 2000 km. While this trend could reflect reduced tomographic resolution in the middle mantle, we show that it may alternatively relate to real changes in the time-dependent subduction flux and/or a mid-lower mantle viscosity increase. An apparent secondary peak in agreement below 2500 km depth may reflect the degree-two lower mantle slow seismic structures. Vote maps illustrate the potential shortcomings of using a limited number or type of tomography models and slab threshold criteria.

Compositionally-distinct ultra-low velocity zones on Earth's core-mantle boundary

The Earth’s lowermost mantle large low velocity provinces are accompanied by small-scale ultralow velocity zones in localized regions on the core-mantle boundary. Large low velocity provinces are hypothesized to be caused by large-scale compositional heterogeneity (i.e., thermochemical piles). The origin of ultralow velocity zones, however, remains elusive. Here we perform three-dimensional geodynamical calculations to show that the current locations and shapes of ultralow velocity zones are related to their cause. We find that the hottest lowermost mantle regions are commonly located well within the interiors of thermochemical piles. In contrast, accumulations of ultradense compositionally distinct material occur as discontinuous patches along the margins of thermochemical piles and have asymmetrical cross-sectional shape. Furthermore, the lateral morphology of these patches provides insight into mantle flow directions and long-term stability. The global distribution and large variations of morphology of ultralow velocity zones validate a compositionally distinct origin for most ultralow velocity zones.

Ultralow velocity zones are detected on the core-mantle boundary, but their origin is enigmatic. Here, the authors find that the global distribution and large variations of morphology of ultralow velocity zones are consistent with most having a compositionally-distinct origin.

Density structure of Earth's lowermost mantle from Stoneley mode splitting observations

Advances in our understanding of Earth's thermal evolution and the style of mantle convection rely on robust seismological constraints on lateral variations of density. The large-low-shear-wave velocity provinces (LLSVPs) atop the core-mantle boundary beneath Africa and the Pacific are the largest structures in the lower mantle, and hence severely affect the convective flow. Here, we show that anomalous splitting of Stoneley modes, a unique class of free oscillations that are perturbed primarily by velocity and density variations at the core-mantle boundary, is explained best when the overall density of the LLSVPs is lower than the surrounding mantle. The resolved density variations can be explained by the presence of post-perovskite, chemical heterogeneity or a combination of the two. Although we cannot rule out the presence of a ∼100-km-thick denser-than-average basal structure, our results support the hypothesis that LLSVPs signify large-scale mantle upwelling in two antipodal regions of the mantle.

On the deep-mantle origin of the Deccan Traps

The Deccan Traps in west-central India constitute one of Earth's largest continental flood basalt provinces, whose eruption played a role in the Cretaceous-Paleogene extinction event. The unknown mantle structure under the Indian Ocean at the start of the Cenozoic presents a challenge for connecting the event to a deep mantle origin. We used a back-and-forth iterative method for time-reversed convection modeling, which incorporates tomography-based, present-day mantle heterogeneity to reconstruct mantle structure at the start of the Cenozoic. We show a very low-density, deep-seated upwelling that ascends beneath the Réunion hot spot at the time of the Deccan eruptions. We found a second active upwelling below the Comores hot spot that likely contributed to the region of partial melt feeding the massive eruption.

Origin and evolution of the deep thermochemical structure beneath Eurasia

A unique structure in the Earth's lowermost mantle, the Perm Anomaly, was recently identified beneath Eurasia. It seismologically resembles the large low-shear velocity provinces (LLSVPs) under Africa and the Pacific, but is much smaller. This challenges the current understanding of the evolution of the plate–mantle system in which plumes rise from the edges of the two LLSVPs, spatially fixed in time. New models of mantle flow over the last 230 million years reproduce the present-day structure of the lower mantle, and show a Perm-like anomaly. The anomaly formed in isolation within a closed subduction network ∼22,000 km in circumference prior to 150 million years ago before migrating ∼1,500 km westward at an average rate of 1 cm year⁻¹, indicating a greater mobility of deep mantle structures than previously recognized. We hypothesize that the mobile Perm Anomaly could be linked to the Emeishan volcanics, in contrast to the previously proposed Siberian Traps.

The Perm anomaly is found in the lower mantle beneath Eurasia, but how this structure formed has remained unclear. Here, the authors show that the anomaly has been mobile since it formed in isolation within a closed subduction network and propose that the anomaly is linked to the Emeishan volcanics.

Kinematics and dynamics of the East Pacific Rise linked to a stable, deep-mantle upwelling

Earth's tectonic plates are generally considered to be driven largely by negative buoyancy associated with subduction of oceanic lithosphere. In this context, mid-ocean ridges (MORs) are passive plate boundaries whose divergence accommodates flow driven by subduction of oceanic slabs at trenches. We show that over the past 80 million years (My), the East Pacific Rise (EPR), Earth's dominant MOR, has been characterized by limited ridge-perpendicular migration and persistent, asymmetric ridge accretion that are anomalous relative to other MORs. We reconstruct the subduction-related buoyancy fluxes of plates on either side of the EPR. The general expectation is that greater slab pull should correlate with faster plate motion and faster spreading at the EPR. Moreover, asymmetry in slab pull on either side of the EPR should correlate with either ridge migration or enhanced plate velocity in the direction of greater slab pull. Based on our analysis, none of the expected correlations are evident. This implies that other forces significantly contribute to EPR behavior. We explain these observations using mantle flow calculations based on globally integrated buoyancy distributions that require core-mantle boundary heat flux of up to 20 TW. The time-dependent mantle flow predictions yield a long-lived deep-seated upwelling that has its highest radial velocity under the EPR and is inferred to control its observed kinematics. The mantle-wide upwelling beneath the EPR drives horizontal components of asthenospheric flows beneath the plates that are similarly asymmetric but faster than the overlying surface plates, thereby contributing to plate motions through viscous tractions in the Pacific region.

Global correlation of lower mantle structure and past subduction

Advances in global seismic tomography have increasingly motivated identification of subducted lithosphere in Earth's deep mantle, creating novel opportunities to link plate tectonics and mantle evolution. Chief among those is the quest for a robust subduction reference frame, wherein the mantle assemblage of subducted lithosphere is used to reconstruct past surface tectonics in an absolute framework anchored in the deep Earth. However, the associations heretofore drawn between lower mantle structure and past subduction have been qualitative and conflicting, so the very assumption of a correlation has yet to be quantitatively corroborated. Here we show that a significant, time-depth progressive correlation can be drawn between reconstructed subduction zones of the last 130 Myr and positive S wave velocity anomalies at 600-2300 km depth, but that further correlation between greater times and depths is not presently demonstrable. This correlation suggests that lower mantle slab sinking rates average between 1.1 and 1.9 cm yr-1.

Antiquity of the South Atlantic Anomaly and evidence for top-down control on the geodynamo

The dramatic decay of dipole geomagnetic field intensity during the last 160 years coincides with changes in Southern Hemisphere (SH) field morphology and has motivated speculation of an impending reversal. Understanding these changes, however, has been limited by the lack of longer-term SH observations. Here we report the first archaeomagnetic curve from southern Africa (ca. 1000–1600 AD). Directions change relatively rapidly at ca. 1300 AD, whereas intensities drop sharply, at a rate greater than modern field changes in southern Africa, and to lower values. We propose that the recurrence of low field strengths reflects core flux expulsion promoted by the unusual core–mantle boundary (CMB) composition and structure beneath southern Africa defined by the African large low shear velocity province (LLSVP). Because the African LLSVP and CMB structure are ancient, this region may have been a steady site for flux expulsion, and triggering of geomagnetic reversals, for millions of years.

The rapid decay of Earth's dipole magnetic field has recently captured the public imagination. Here, the authors present a southern hemisphere magnetic record from South African Iron Age sites using oriented samples in the floors and suggest that the anomalous field behaviour is not just a recent feature.

Spin crossover in ferropericlase and velocity heterogeneities in the lower mantle

Deciphering the origin of seismic velocity heterogeneities in the mantle is crucial to understanding internal structures and processes at work in the Earth. The spin crossover in iron in ferropericlase (Fp), the second most abundant phase in the lower mantle, introduces unfamiliar effects on seismic velocities. First-principles calculations indicate that anticorrelation between shear velocity (VS) and bulk sound velocity (Vφ) in the mantle, usually interpreted as compositional heterogeneity, can also be produced in homogeneous aggregates containing Fp. The spin crossover also suppresses thermally induced heterogeneity in longitudinal velocity (VP) at certain depths but not in VS. This effect is observed in tomography models at conditions where the spin crossover in Fp is expected in the lower mantle. In addition, the one-of-a-kind signature of this spin crossover in the RS/P (∂ ln VS/∂ ln VP) heterogeneity ratio might be a useful fingerprint to detect the presence of Fp in the lower mantle.