Mantle shear-wave tomography and the fate of subducted slabs

Geodynamically corrected Pliocene shoreline elevations in Australia consistent with midrange projections of Antarctic ice loss

The Mid-Pliocene represents the most recent interval in Earth history with climatic conditions similar to those expected in the coming decades. Mid-Pliocene sea level estimates therefore provide important constraints on projections of future ice sheet behavior and sea level change but differ by tens of meters due to local distortion of paleoshorelines caused by mantle dynamics. We combine an Australian sea level marker compilation with geodynamic simulations and probabilistic inversions to quantify and remove these post-Pliocene vertical motions at continental scale. Dynamic topography accounts for most of the observed sea level marker deflection, and correcting for this effect and glacial isostatic adjustment yields a Mid-Pliocene global mean sea level of +16.0 (+10.4 to +21.5) m (50th/16th to 84th percentiles). Recalibration of recent high-end sea level projections using this revised estimate implies a more stable Antarctic Ice Sheet under future warming scenarios, consistent with midrange forecasts of sea level rise that do not incorporate a marine ice cliff instability.

Influence of reef isostasy, dynamic topography, and glacial isostatic adjustment on sea-level records in Northeastern Australia

Understanding sea level during the peak of the Last Interglacial (125,000 yrs ago) is important for assessing future ice-sheet dynamics in response to climate change. The coasts and continental shelves of northeastern Australia (Queensland) preserve an extensive Last Interglacial record in the facies of coastal strandplains onland and fossil reefs offshore. However, there is a discrepancy, amounting to tens of meters, in the elevation of sea-level indicators between offshore and onshore sites. Here, we assess the influence of geophysical processes that may have changed the elevation of these sea-level indicators. We modeled sea-level change due to dynamic topography, glacial isostatic adjustment, and isostatic adjustment due to coral reef loading. We find that these processes caused relative sea-level changes on the order of, respectively, 10 m, 5 m, and 0.3 m. Of these geophysical processes, the dynamic topography predictions most closely match the tilting observed between onshore and offshore sea-level markers.

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.

Glacial isostatic adjustment: physical models and observational constraints

By far the most prescient insights into the interior structure of the planet have been provided on the basis of elastic wave seismology. Analysis of the travel times of shear or compression wave phases excited by individual earthquakes, or through analysis of the elastic gravitational free oscillations that individual earthquakes of sufficiently large magnitude may excite, has been the central focus of Earth physics research for more than a century. Unfortunately, data provide no information that is directly relevant to understanding the solid state 'flow' of the polycrystalline outer 'mantle' shell of the planet that is involved in the thermally driven convective circulation that is responsible for powering the 'drift' of the continents and which controls the rate of planetary cooling on long timescales. For this reason, there has been an increasing focus on the understanding of physical phenomenology that is unambiguously associated with mantle flow processes that are distinct from those directly associated with the convective circulation itself. This paper reviews the past many decades of work that has been invested in understanding the most important of such processes, namely that which has come to be referred to as 'glacial isostatic adjustment' (GIA). This process concerns the response of the planet to the loading and unloading of the high latitude continents by the massive accumulations of glacial ice that have occurred with almost metronomic regularity over the most recent million years of Earth history. Forced by the impact of gravitationaln-body effects on the geometry of Earth's orbit around the Sun through the impact upon the terrestrial regime of received solar insolation, these surface mass loads on the continents have left indelible records of their occurrence in the 'Earth system' consisting of the oceans, continents, and the great polar ice sheets on Greenland and Antarctica themselves. Although this ice-age phenomenology has been clearly recognized since early in the last century, it was for over 50 years considered to be no more than an interesting curiosity, the understanding of which remained on the periphery of the theoretical physics of the Earth. This was the case in part because no globally applicable theory was available that could be applied to rigorously interpret the observations. Equally important to understanding the scientific lethargy that held back the understanding of this phenomenon involving mantle flow processes was the lack of appreciation of the wide range of observations that were in fact related to GIA physics. This paper is devoted to a review of the global theories of the GIA process that have since been developed as a means of interpreting the extensive variety of observations that are now recognized as being involved in the response of the planet to the loading and unloading of its surface by glacial ice. The paper will also provide examples of the further analyses of Earth physics and climate related processes that applications of the modern theoretical structures have enabled.

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.

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.

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.

Anomalous SmKS induced by postcritical reflection and refraction at the core-mantle boundary

Earth's outer core is generally thought to be a well-mixed liquid consisting mostly of iron and a small amount of lighter elements. Recent seismic studies using SmKS waves show that the top a few hundred kilometers of the outer core possess a P-wave velocity slightly lower than the PREM model, which cannot be explained by self-compression of a chemically homogeneous outer core. We investigated the SmKS waveforms of a deep earthquake occurring beneath South America recorded by a large and dense seismic array in China, and measured the differential arrival times of the SmKS pairs. We found significant waveform distortion of the SmKS caused by postcritical refraction and reflection at the core-mantle boundary. This waveform distortion can introduce significant bias to the measured differential times, leading to incorrect estimate of P-wave velocity of the outer core. Whether stable stratification is occurring in outer core or not requires further seismic investigations.

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.

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.

Prediction of Mechanical Twinning in Magnesium Silicate Post-Perovskite

The plastic properties of MgSiO3 post-perovskite are considered to be one of the key issues necessary for understanding the seismic anisotropy at the bottom of the mantle in the so-called D" layer. Although plastic slip in MgSiO3 post-perovskite has attracted considerable attention, the twinning mechanism has not been addressed, despite some experimental evidence from low-pressure analogues. On the basis of a numerical mechanical model, we present a twin nucleation model for post-perovskite involving the emission of 1/6 <110> partial dislocations. Relying on first-principles calculations with no adjustable parameters, we show that {110} twin wall formation resulting from the interaction of multiple twin dislocations occurs at a twinning stress comparable in magnitude to the most readily occurring slip system in post-perovskite. Because dislocation activities and twinning are competitive strain-producing mechanisms, twinning should be considered in future models of crystallographic preferred orientations in post-perovskite to better interpret seismic anisotropy in the lowermost lower mantle.

Imaging paleoslabs in the D″ layer beneath Central America and the Caribbean using seismic waveform inversion

D″ (Dee double prime), the lowermost layer of the Earth's mantle, is the thermal boundary layer (TBL) of mantle convection immediately above the Earth's liquid outer core. As the origin of upwelling of hot material and the destination of paleoslabs (downwelling cold slab remnants), D″ plays a major role in the Earth's evolution. D″ beneath Central America and the Caribbean is of particular geodynamical interest, because the paleo- and present Pacific plates have been subducting beneath the western margin of Pangaea since ~250 million years ago, which implies that paleoslabs could have reached the lowermost mantle. We conduct waveform inversion using a data set of ~7700 transverse component records to infer the detailed three-dimensional S-velocity structure in the lowermost 400 km of the mantle in the study region so that we can investigate how cold paleoslabs interact with the hot TBL above the core-mantle boundary (CMB). We can obtain high-resolution images because the lowermost mantle here is densely sampled by seismic waves due to the full deployment of the USArray broadband seismic stations during 2004-2015. We find two distinct strong high-velocity anomalies, which we interpret as paleoslabs, just above the CMB beneath Central America and Venezuela, respectively, surrounded by low-velocity regions. Strong low-velocity anomalies concentrated in the lowermost 100 km of the mantle suggest the existence of chemically distinct denser material connected to low-velocity anomalies in the lower mantle inferred by previous studies, suggesting that plate tectonics on the Earth's surface might control the modality of convection in the lower mantle.

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.

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.

Origins of ultralow velocity zones through slab-derived metallic melt

Understanding the ultralow velocity zones (ULVZs) places constraints on the chemical composition and thermal structure of deep Earth and provides critical information on the dynamics of large-scale mantle convection, but their origin has remained enigmatic for decades. Recent studies suggest that metallic iron and carbon are produced in subducted slabs when they sink beyond a depth of 250 km. Here we show that the eutectic melting curve of the iron-carbon system crosses the current geotherm near Earth's core-mantle boundary, suggesting that dense metallic melt may form in the lowermost mantle. If concentrated into isolated patches, such melt could produce the seismically observed density and velocity features of ULVZs. Depending on the wetting behavior of the metallic melt, the resultant ULVZs may be short-lived domains that are replenished or regenerated through subduction, or long-lasting regions containing both metallic and silicate melts. Slab-derived metallic melt may produce another type of ULVZ that escapes core sequestration by reacting with the mantle to form iron-rich postbridgmanite or ferropericlase. The hypotheses connect peculiar features near Earth's core-mantle boundary to subduction of the oceanic lithosphere through the deep carbon cycle.

Multiple seismic reflectors in Earth's lowermost mantle

The modern view of Earth's lowermost mantle considers a D″ region of enhanced (seismologically inferred) heterogeneity bounded by the core-mantle boundary and an interface some 150-300 km above it, with the latter often attributed to the postperovskite phase transition (in MgSiO3). Seismic exploration of Earth's deep interior suggests, however, that this view needs modification. So-called ScS and SKKS waves, which probe the lowermost mantle from above and below, respectively, reveal multiple reflectors beneath Central America and East Asia, two areas known for subduction of oceanic plates deep into Earth's mantle. This observation is inconsistent with expectations from a thermal response of a single isochemical postperovskite transition, but some of the newly observed structures can be explained with postperovskite transitions in differentiated slab materials. Our results imply that the lowermost mantle is more complex than hitherto thought and that interfaces and compositional heterogeneity occur beyond the D″ region sensu stricto.

Orogens in the evolving Earth: from surface continents to ‘lost continents’ at the core–mantle boundary

Orogens and their posthumous traces are the basic elements that can be used to understand the material circulation within the Earth. Although information preserved in the rocks on the surface ranging in age from 4.4 Ga to the present has been used to characterize orogens, it is important to understand orogens on a whole-Earth scale to evaluate global material circulation through time. In this paper, we synthesize the general concepts and characteristics of orogens and orogenic belts. The collision type and accretionary type constitute the two end-member types of orogens, both sharing similar structural features of subhorizontal disposition, bounded above and below by paired faults. Their exhumation generally occurs in two steps: first by wedge extrusion to form a sandwich structure with subhorizontal boundaries, which is followed by domal uplift of all the units. In the accretionary type, oceanic lithosphere subducts under the continental margin, and in the collision type, buoyant continents collide with each other. Of the various types of subduction and collision processes, arc–arc collision orogeny may have been widespread in the Archaean, although most of the intra-oceanic arc crust must have been destroyed and dragged down to the Archaean core–mantle boundary (CMB). Here we propose a broad two-fold classification of orogens and their subducted remnants, based on (1) their thermal history and (2) temporal constraints. Based on their thermal history, orogens are grouped into three types: cold orogens, hot orogens and ultra-hot orogens. Two extreme situations, which are anomalous and unlikely to occur on Earth, termed here super-cold and super-hot orogens, are also proposed. We discuss the characteristics of each of these subtypes. Based on temporal constraints, we group orogens into Modern and Ancient, where in both cases regional metamorphic belts occupy the orogenic core. In both groups, the overlying and underlying units of the regional metamorphic belts are weakly metamorphosed or unmetamorphosed, and are either accretionary complex in origin (Pacific type) or continental basement and cover (collision type). Major structures are subhorizontal with oceanward vergence of deformation, for both types. Orogens in the Modern Earth are grouped into four sub-categories: (1) deeply subducted orogens that are taken down to mantle depths and never return to the surface, termed here ‘ghost orogens’; (2) those that are subducted to deep crustal levels, undergo melting and are recycled back to the surface, forming resurrected and temporarily ‘arrested orogens’; (3) ‘extant orogens’, which are partly returned to the surface after deep subduction; (4) ‘concealed orogens’, which have been deeply subducted and only the traces of which are represented on the surface by mantle xenoliths carried by younger magmas. The preservation of orogens on the surface of the Earth occurred through an unusual return process from their natural course of total destruction, a phenomenon that operated more efficiently in the Phanerozoic through exhumation from ultra-deep domains against the slab-pull force of the plate, aided by fluids derived by dehydration of subducted lithosphere. Orogens at present represented on the surface of the Earth constitute only a fraction of the total volume formed in Earth history. Traces of the deeply subducted ‘lost orogens’ are sometimes returned to the surface in the form of melt or mantle xenoliths through combined processes of plume and plate tectonics. From a synthesis of the processes associated with the various categories of orogens proposed in this study, we trace the time-dependent transformations of orogens in relation to the history of the evolving Earth.

Temperature profile in the lowermost mantle from seismological and mineral physics joint modeling

The internal structure of the core-mantle boundary (CMB) region of the Earth plays a crucial role in controlling the dynamics and evolution of our planet. We have developed a comprehensive model based on the radial variations of shear velocity in the D'' layer (the base of the lower mantle) and the high-P,T elastic properties of major candidate minerals, including the effects of post-perovskite phase changes. This modeling shows a temperature profile in the lowermost mantle with a CMB temperature of 3,800 +/- 200 K, which suggests that lateral temperature variations of 200-300 K can explain much of the large velocity heterogeneity observed in D''. A single-crossing phase transition model was found to be more favorable in reproducing the observed seismic wave velocity structure than a double-crossing phase transition model.

Insights into the nature of the transition zone from physically constrained inversion of long-period seismic data

Imposing a thermal and compositional significance to the outcome of the inversion of seismic data facilitates their interpretation. Using long-period seismic waveforms and an inversion approach that includes constraints from mineral physics, we find that lateral variations of temperature can explain a large part of the data in the upper mantle. The additional compositional signature of cratons emerges in the global model as well. Above 300 km, we obtain seismic geotherms that span the range of expected temperatures in various tectonic regions. Absolute velocities and gradients with depth are well constrained by the seismic data throughout the upper mantle, except near discontinuities. The seismic data are consistent with a slower transition zone and an overall faster shallow upper mantle, which is not compatible with a homogenous dry pyrolite composition. A gradual enrichment with depth in a garnet-rich component helps to reduce the observed discrepancies. A hydrated transition zone would help to lower the velocities in the transition zone, but it does not explain the seismic structure above it.

Seismostratigraphy and thermal structure of Earth's core-mantle boundary region

We used three-dimensional inverse scattering of core-reflected shear waves for large-scale, high-resolution exploration of Earth's deep interior (D'') and detected multiple, piecewise continuous interfaces in the lowermost layer (D'') beneath Central and North America. With thermodynamic properties of phase transitions in mantle silicates, we interpret the images and estimate in situ temperatures. A widespread wave-speed increase at 150 to 300 kilometers above the coremantle boundary is consistent with a transition from perovskite to postperovskite. Internal D'' stratification may be due to multiple phase-boundary crossings, and a deep wave-speed reduction may mark the base of a postperovskite lens about 2300 kilometers wide and 250 kilometers thick. The core-mantle boundary temperature is estimated at 3950 +/- 200 kelvin. Beneath Central America, a site of deep subduction, the D'' is relatively cold (DeltaT = 700 +/- 100 kelvin). Accounting for a factor-of-two uncertainty in thermal conductivity, core heat flux is 80 to 160 milliwatts per square meter (mW m(-2)) into the coldest D'' region and 35 to 70 mW m(-2) away from it. Combined with estimates from the central Pacific, this suggests a global average of 50 to 100 mW m(-2) and a total heat loss of 7.5 to 15 terawatts.

Seismological support for the metastable superplume model, sharp features, and phase changes within the lower mantle

Recently, a metastable thermal-chemical convection model was proposed to explain the African Superplume. Its bulk tabular shape remains relatively stable while its interior undergoes significant stirring with low-velocity conduits along its edges and down-welling near the middle. Here, we perform a mapping of chemistry and temperature into P and S velocity variations and replace a seismically derived structure with this hybrid model. Synthetic seismogram sections generated for this 2D model are then compared directly with corresponding seismic observations of P (P, P(C)P, and PKP) and S (S, S(C)S, and SKS) phases. These results explain the anticorrelation between the bulk velocity and shear velocity and the sharpness and level of SKS travel time delays. In addition, we present evidence for the existence of a D" triplication (a putative phase change) beneath the down-welling structure.

Evidence for a chemical-thermal structure at base of mantle from sharp lateral P-wave variations beneath Central America

Compressional waves that sample the lowermost mantle west of Central America show a rapid change in travel times of up to 4 s over a sampling distance of 300 km and a change in waveforms. The differential travel times of the PKP waves (which traverse Earth's core) correlate remarkably well with predictions for S-wave tomography. Our modeling suggests a sharp transition in the lowermost mantle from a broad slow region to a broad fast region with a narrow zone of slowest anomaly next to the boundary beneath the Cocos Plate and the Caribbean Plate. The structure may be the result of ponding of ancient subducted Farallon slabs situated near the edge of a thermal and chemical upwelling.

Seismic detection of folded, subducted lithosphere at the core-mantle boundary

Seismic tomography has been used to infer that some descending slabs of oceanic lithosphere plunge deep into the Earth's lower mantle. The fate of these slabs has remained unresolved, but it has been postulated that their ultimate destination is the lowermost few hundred kilometres of the mantle, known as the D'' region. Relatively cold slab material may account for high seismic velocities imaged in D'' beneath areas of long-lived plate subduction, and for reflections from a seismic velocity discontinuity just above the anomalously high wave speed regions. The D'' discontinuity itself is probably the result of a phase change in relatively low-temperature magnesium silicate perovskite. Here, we present images of the D'' region beneath the Cocos plate using Kirchhoff migration of horizontally polarized shear waves, and find a 100-km vertical step occurring over less than 100 km laterally in an otherwise flat D'' shear velocity discontinuity. Folding and piling of a cold slab that has reached the core-mantle boundary, as observed in numerical and experimental models, can account for the step by a 100-km elevation of the post-perovskite phase boundary due to a 700 degrees C lateral temperature reduction in the folded slab. We detect localized low velocities at the edge of the slab material, which may result from upwellings caused by the slab laterally displacing a thin hot thermal boundary layer.

Thermochemical structures beneath Africa and the Pacific Ocean

Large low-velocity seismic anomalies have been detected in the Earth's lower mantle beneath Africa and the Pacific Ocean that are not easily explained by temperature variations alone. The African anomaly has been interpreted to be a northwest-southeast-trending structure with a sharp-edged linear, ridge-like morphology. The Pacific anomaly, on the other hand, appears to be more rounded in shape. Mantle models with heterogeneous composition have related these structures to dense thermochemical piles or superplumes. It has not been shown, however, that such models can lead to thermochemical structures that satisfy the geometrical constraints, as inferred from seismological observations. Here we present numerical models of thermochemical convection in a three-dimensional spherical geometry using plate velocities inferred for the past 119 million years. We show that Earth's subduction history can lead to thermochemical structures similar in shape to the observed large, lower-mantle velocity anomalies. We find that subduction history tends to focus dense material into a ridge-like pile beneath Africa and a relatively more-rounded pile under the Pacific Ocean, consistent with seismic observations.

Seismological constraints on a possible plume root at the core-mantle boundary

Recent seismological discoveries have indicated that the Earth's core-mantle boundary is far more complex than a simple boundary between the molten outer core and the silicate mantle. Instead, its structural complexities probably rival those of the Earth's crust. Some regions of the lowermost mantle have been observed to have seismic wave speed reductions of at least 10 per cent, which appear not to be global in extent. Here we present robust evidence for an 8.5-km-thick and approximately 50-km-wide pocket of dense, partially molten material at the core-mantle boundary east of Australia. Array analyses of an anomalous precursor to the reflected seismic wave ScP reveal compressional and shear-wave velocity reductions of 8 and 25 per cent, respectively, and a 10 per cent increase in density of the partially molten aggregate. Seismological data are incompatible with a basal layer composed of pure melt, and thus require a mechanism to prevent downward percolation of dense melt within the layer. This may be possible by trapping of melt by cumulus crystal growth following melt drainage from an anomalously hot overlying region of the lowermost mantle. This magmatic evolution and the resulting cumulate structure seem to be associated with overlying thermal instabilities, and thus may mark a root zone of an upwelling plume.