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

Terrestrial Single-Station Analog for Constraining the Martian Core and Deep Interior: Implications for InSight

We used a terrestrial single-station seismometer to quantify the uncertainty of InSight (INterior explorations using Seismic Investigations, Geodesy and Heat Transport) data for determining Martian core size. To mimic Martian seismicity, we formed a catalog using 917 terrestrial earthquakes, from which we randomly selected events. We stacked ScS amplitudes on modeled arrival times and searched for where ScS produced coherent seismic amplitudes. A core detection was defined by a coherent peak with small offset between predicted and user-selected arrival times. Iterating the detection algorithm with varying signal-to-noise (SNR) ranges and quantity of events determined the selection frequency of each model and quantified core depth uncertainty. Increasing the quantity of events reduced core depth uncertainty while increasing the recovery rate, while increasing event SNR had little effect. Including ScS2 multiples increased the recovery rate and reduced core depth uncertainty when we used low quantities of events. The most-frequent core depths varied by back azimuth, suggesting our method is sensitive to the presence of mantle heterogeneities. When we added 1° in source distance errors, core depth uncertainty increased by up to 11 km and recovery rates decreased by <5%. Altering epicentral distances by 25% added ~35 km of uncertainty and reduced recovery rates to <50% in some cases. From these experiments, we estimate that if InSight can detect five events with high location precision (<10 % epicentral distance errors), that there is at least an 88% chance of core depth recovery using ScS alone with uncertainty in core depth approaching 18 km and decreasing as more events are located.

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.

Restoration of clipped seismic waveforms using projection onto convex sets method

The seismic waveforms would be clipped when the amplitude exceeds the upper-limit dynamic range of seismometer. Clipped waveforms are typically assumed not useful and seldom used in waveform-based research. Here, we assume the clipped components of the waveform share the same frequency content with the un-clipped components. We leverage this similarity to convert clipped waveforms to true waveforms by iteratively reconstructing the frequency spectrum using the projection onto convex sets method. Using artificially clipped data we find that statistically the restoration error is ~1% and ~5% when clipped at 70% and 40% peak amplitude, respectively. We verify our method using real data recorded at co-located seismometers that have different gain controls, one set to record large amplitudes on scale and the other set to record low amplitudes on scale. Using our restoration method we recover 87 out of 93 clipped broadband records from the 2013 Mw6.6 Lushan earthquake. Estimating that we recover 20 clipped waveforms for each M5.0+ earthquake, so for the ~1,500 M5.0+ events that occur each year we could restore ~30,000 clipped waveforms each year, which would greatly enhance useable waveform data archives. These restored waveform data would also improve the azimuthal station coverage and spatial footprint.

On the enigmatic birth of the Pacific Plate within the Panthalassa Ocean

The oceanic Pacific Plate started forming in Early Jurassic time within the vast Panthalassa Ocean that surrounded the supercontinent Pangea, and contains the oldest lithosphere that can directly constrain the geodynamic history of the circum-Pangean Earth. We show that the geometry of the oldest marine magnetic anomalies of the Pacific Plate attests to a unique plate kinematic event that sparked the plate's birth at virtually a point location, surrounded by the Izanagi, Farallon, and Phoenix Plates. We reconstruct the unstable triple junction that caused the plate reorganization, which led to the birth of the Pacific Plate, and present a model of the plate tectonic configuration that preconditioned this event. We show that a stable but migrating triple junction involving the gradual cessation of intraoceanic Panthalassa subduction culminated in the formation of an unstable transform-transform-transform triple junction. The consequent plate boundary reorganization resulted in the formation of a stable triangular three-ridge system from which the nascent Pacific Plate expanded. We link the birth of the Pacific Plate to the regional termination of intra-Panthalassa subduction. Remnants thereof have been identified in the deep lower mantle of which the locations may provide paleolongitudinal control on the absolute location of the early Pacific Plate. Our results constitute an essential step in unraveling the plate tectonic evolution of "Thalassa Incognita" that comprises the comprehensive Panthalassa Ocean surrounding Pangea.

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.

Mineralogical effects on the detectability of the postperovskite boundary

The discovery of a phase transition in Mg-silicate perovskite (Pv) to postperovskite (pPv) at lowermost mantle pressure-temperature (P - T) conditions may provide an explanation for the discontinuous increase in shear wave velocity found in some regions at a depth range of 200 to 400 km above the core-mantle boundary, hereafter the D('') discontinuity. However, recent studies on binary and ternary systems showed that reasonable contents of Fe(2+) and Al for pyrolite increase the thickness (width of the mixed phase region) of the Pv - pPv boundary (400-600 km) to much larger than the D('') discontinuity (≤ 70 km). These results challenge the assignment of the D('') discontinuity to the Pv - pPv boundary in pyrolite (homogenized mantle composition). Furthermore, the mineralogy and composition of rocks that can host a detectable Pv → pPv boundary are still unknown. Here we report in situ measurements of the depths and thicknesses of the Pv → pPv transition in multiphase systems (San Carlos olivine, pyrolitic, and midocean ridge basaltic compositions) at the P - T conditions of the lowermost mantle, searching for candidate rocks with a sharp Pv - pPv discontinuity. Whereas the pyrolitic mantle may not have a seismologically detectable Pv → pPv transition due to the effect of Al, harzburgitic compositions have detectable transitions due to low Al content. In contrast, Al-rich basaltic compositions may have a detectable Pv - pPv boundary due to their distinct mineralogy. Therefore, the observation of the D('') discontinuity may be related to the Pv → pPv transition in the differentiated oceanic lithosphere materials transported to the lowermost mantle by subducting slabs.

The contemporary degassing rate of 40Ar from the solid Earth

Knowledge of the outgassing history of radiogenic (40)Ar, derived over geologic time from the radioactive decay of (40)K, contributes to our understanding of the geodynamic history of the planet and the origin of volatiles on Earth's surface. The (40)Ar inventory of the atmosphere equals total (40)Ar outgassing during Earth history. Here, we report the current rate of (40)Ar outgassing, accessed by measuring the Ar isotope composition of trapped gases in samples of the Vostok and Dome C deep ice cores dating back to almost 800 ka. The modern outgassing rate (1.1 +/- 0.1 x 10(8) mol/yr) is in the range of values expected by summing outgassing from the continental crust and the upper mantle, as estimated from simple calculations and models. The measured outgassing rate is also of interest because it allows dating of air trapped in ancient ice core samples of unknown age, although uncertainties are large (+/-180 kyr for a single sample or +/-11% of the calculated age, whichever is greater).

Structure and dynamics of Earth's lower mantle

Processes within the lowest several hundred kilometers of Earth's rocky mantle play a critical role in the evolution of the planet. Understanding Earth's lower mantle requires putting recent seismic and mineral physics discoveries into a self-consistent, geodynamically feasible context. Two nearly antipodal large low-shear-velocity provinces in the deep mantle likely represent chemically distinct and denser material. High-resolution seismological studies have revealed laterally varying seismic velocity discontinuities in the deepest few hundred kilometers, consistent with a phase transition from perovskite to post-perovskite. In the deepest tens of kilometers of the mantle, isolated pockets of ultralow seismic velocities may denote Earth's deepest magma chamber.

Boron and oxygen isotope evidence for recycling of subducted components over the past 2.5 Gyr

Evidence for the deep recycling of surficial materials through the Earth's mantle and their antiquity has long been sought to understand the role of subducting plates and plumes in mantle convection. Radiogenic isotope evidence for such recycling remains equivocal because the age and location of parent-daughter fractionation are not known. Conversely, while stable isotopes can provide irrefutable evidence for low-temperature fractionation, their range in most unaltered oceanic basalts is limited and the age of any variation is unconstrained. Here we show that delta(18)O ratios in basalts from the Azores are often lower than in pristine mantle. This, combined with increased Nb/B ratios and a large range in delta(11)B ratios, provides compelling evidence for the recycling of materials that had undergone fractionation near the Earth's surface. Moreover, delta(11)B is negatively correlated with (187)Os/(188)Os ratios, which extend to subchondritic values, constraining the age of the high Nb/B, (11)B-enriched endmember to be more than 2.5 billion years (Gyr) old. We infer this component to be melt- and fluid-depleted lithospheric mantle from a subducted oceanic plate, whereas other Azores basalts contain a contribution from approximately 3-Gyr-old melt-enriched basalt. We conclude that both components are most probably derived from an Archaean oceanic plate that was subducted, arguably into the deep mantle, where it was stored until thermal buoyancy caused it to rise beneath the Azores islands approximately 3 Gyr later.

Topology of the postperovskite phase transition and mantle dynamics

The postperovskite (ppv) phase transition occurs in the deep mantle close to the core-mantle boundary (CMB). For this reason, we must include in the dynamical considerations both the Clapeyron slope and the temperature intercept, T(int), which is the temperature of the phase transition at the CMB pressure. For a CMB temperature greater than T(int), there is a double crossing of the phase boundary by the geotherms associated with the descending flow. We have found a great sensitivity of the shape of the ppv surface due to the CMB from variations of various parameters such as the amount of internal heating, the Clapeyron slope, and the temperature intercept. Three-dimensional spherical models of mantle convection that can satisfy the seismological constraints depend on the Clapeyron slope. At moderate value, 8 MPa/K, the best fit is found with a core heat flow amounting for 40% of the total heat budget (approximately equal to 15 TW), whereas for 10 MPa/K the agreement is for a lower core heat flow (20%, approximately equal to 7.5 TW). In all cases, these solutions correspond to a temperature intercept 200 K lower than the CMB temperature. These models have holes of perovskite adjacent to ppv in regions of hot upwellings.

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.

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.

A Post-Perovskite Lens and D'' Heat Flux Beneath the Central Pacific

Temperature gradients in a low-shear-velocity province in the lowermost mantle (D'' region) beneath the central Pacific Ocean were inferred from the observation of a rapid S-wave velocity increase overlying a rapid decrease. These paired seismic discontinuities are attributed to a phase change from perovskite to post-perovskite and then back to perovskite as the temperature increases with depth. Iron enrichment could explain the occurrence of post-perovskite several hundred kilometers above the core-mantle boundary in this warm, chemically distinct province. The double phase-boundary crossing directly constrains the lowermost mantle temperature gradients. Assuming a standard but unconstrained choice of thermal conductivity, the regional core-mantle boundary heat flux (approximately 85 +/- 25 milliwatts per square meter), comparable to the average at Earth's surface, was estimated, along with a lower bound on global core-mantle boundary heat flow in the range of 13 +/- 4 terawatts. Mapped velocity-contrast variations indicate that the lens of post-perovskite minerals thins and vanishes over 1000 kilometers laterally toward the margin of the chemical distinct region as a result of a approximately 500-kelvin temperature increase.