Variable azimuthal anisotropy in Earth's lowermost mantle

A persistent reversal in the expected polarity of the initiation of vertically polarized shear waves that graze the D' layer (the layer at the boundary between the outer core and the lower mantle of Earth) in some regions starts at the arrival time of horizontally polarized shear waves. Full waveform modeling of the split shear waves for paths beneath the Caribbean requires azimuthal anisotropy at the base of the mantle. Models with laterally coherent patterns of transverse isotropy with the hexagonal symmetry axis of the mineral phases tilted from the vertical by as much as 20 degrees are consistent with the data. Small-scale convection cells within the mantle above the D' layer may cause the observed variations by inducing laterally variable crystallographic or shape-preferred orientation in minerals in the D' layer.     

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.     

Earthquake Prediction: Absence of a Precursive Change in Seismic Velocities before a Tremor of Magnitude 3frac34

P-wave velocities in the region near the source of a tremor of magnitude 3(3/4) were constant to within 2 percent for 41 days before the event; no evidence of a precursive change in velocity was found. Observations of S-wave velocities and the ratio of P-wave to S-wave velocities also showed no precursive changes. In recent studies, premonitory changes in body-wave velocities of about 10 percent and having a duration of 2 to 3 weeks have been reported for crustal earthquakes of this size.     

Small-scale convective instability and upper mantle viscosity under california

Thermal calculations and convection analysis, constrained by seismic tomography results, suggest that a small-scale convective instability developed in the upper 200 kilometers of the mantle under California after the upwelling and cooling of asthenosphere into the slab window associated with the formation of the San Andreas transform boundary. The upper bound for the upper mantle viscosity in the slab window, 5 x 10(19) pascal seconds, is similar to independent estimates for the asthenosphere beneath young oceanic and tectonically active continental regions. These model calculations suggest that many tectonically active continental regions characterized by low upper mantle seismic velocities may be affected by time-dependent small-scale convection that can generate localized areas of uplift and subsidence.     

Variations of Rayleigh Wave Phase Velocities across the Pacific Ocean

Rayleigh wave phase velocities on paths crossing the Pacific Ocean show variations which are well correlated with some "average" lithospheric age of the region traversed. The seismic velocities in the upper mantle are highest in the oldest parts of the Pacific and lowest in the youngest parts.     

Electronic transitions in perovskite: possible nonconvecting layers in the lower mantle

We measured the spin state of iron in magnesium silicate perovskite (Mg(0.9),Fe(0.1))SiO(3) at high pressure and found two electronic transitions occurring at 70 gigapascals and at 120 gigapascals, corresponding to partial and full electron pairing in iron, respectively. The proportion of iron in the low spin state thus grows with depth, increasing the transparency of the mantle in the infrared region, with a maximum at pressures consistent with the D" layer above the core-mantle boundary. The resulting increase in radiative thermal conductivity suggests the existence of nonconvecting layers in the lowermost mantle.     

Iron-rich post-perovskite and the origin of ultralow-velocity zones

The boundary layer between the crystalline silicate lower mantle and the liquid iron core contains regions with ultralow seismic velocities. Such low compressional and shear wave velocities and high Poisson's ratio are also observed experimentally in post-perovskite silicate phase containing up to 40 mol% FeSiO3 endmember. The iron-rich post-perovskite silicate is stable at the pressure-temperature and chemical environment of the core-mantle boundary and can be formed by core-mantle reaction. Mantle dynamics may lead to further accumulation of this material into the ultralow-velocity patches that are observable by seismology.     

Solidus of Earth's deep mantle

The solidus of a pyrolite-like composition, approximating that of the lower mantle, was measured up to 59 gigapascals by using CO2 laser heating in a diamond anvil cell. The solidus temperatures are at least 700 kelvin below the melting temperatures of magnesiowustite, which in the deep mantle has the lowest melting temperatures of the three major components-magnesiowustite, Mg-Si-perovskite, and Ca-Si-perovskite. The solidus in the deep mantle is more than 1500 kelvin above the average present-day geotherm, but at the core-mantle boundary it is near the core temperature. Thus, partial melting of the mantle is possible at the core-mantle boundary.     

Lateral variations in Compressional/Shear velocities at the base of the mantle

Observations of core-diffracted P (Pdiff) and SH (SHdiff) waves recorded by the Missouri-to-Massachusetts (MOMA) seismic array show that the ratio of compressional (P) seismic velocities to horizontal shear (SH) velocities at the base of the mantle changes abruptly from beneath the mid-Pacific (VP/VS = 1.88, also the value predicted by reference Earth models) to beneath Alaska (VP/VS = 1.83). This change signifies a sudden lateral variation in material properties that may have a mineralogical or textural origin. A textural change could be a result of shear stresses induced during the arrival at the core of ancient lithosphere from the northern Pacific paleotrench.     

The Gutenberg discontinuity: melt at the lithosphere-asthenosphere boundary

The lithosphere-asthenosphere boundary (LAB) beneath ocean basins separates the upper thermal boundary layer of rigid, conductively cooling plates from the underlying ductile, convecting mantle. The origin of a seismic discontinuity associated with this interface, known as the Gutenberg discontinuity (G), remains enigmatic. High-frequency SS precursors sampling below the Pacific plate intermittently detect the G as a sharp, negative velocity contrast at 40- to 75-kilometer depth. These observations lie near the depth of the LAB in regions associated with recent surface volcanism and mantle melt production and are consistent with an intermittent layer of asthenospheric partial melt residing at the lithospheric base. I propose that the G reflectivity is regionally enhanced by dynamical processes that produce melt, including hot mantle upwellings, small-scale convection, and fluid release during subduction.     

Continuous mantle melt supply beneath an overlapping spreading center on the East Pacific Rise

Tomographic images of upper mantle velocity structure beneath an overlapping spreading center (OSC) on the East Pacific Rise indicate that this ridge axis discontinuity is underlain by a continuous region of low P-wave velocities. The anomalous structure can be explained by an approximately 16-kilometer-wide region of high temperatures and melt fractions of a few percent by volume. Our results show that OSCs are not necessarily associated with a discontinuity in melt supply and that both OSC limbs are supplied with melt from a mantle source located beneath the OSC. We conclude that tectonic segmentation of the ridge by OSCs is not the direct result of magmatic segmentation at mantle depths.     

Deep seismic expression of an ancient plate boundary in europe

Seismological results on the structure of the upper mantle below Europe reveal a marked contrast in seismic properties between Precambrian and younger parts of Europe. The Precambrian craton in eastern Europe is characterized by high shear-wave velocities, which can be explained by low temperatures. The transition to low seismic velocities below Phanerozoic Europe coincides with the crustal boundary zone of the craton and exists to depths of at least 140 kilometers. Despite the long and complex tectonic history of the plate boundary zone, the transition is remarkably sharp, which rules out any significant lateral transport of asthenospheric material across the suture zone.     

Sound velocities in olivine at Earth mantle pressures

The independent elastic constants of an upper mantle mineral, San Carlos olivine [(Mg(1.8)Fe(0.2))SiO(4)], were measured from 0 to 12.5 gigapascals. Evidence is offered in support of the proposition that the explicit temperature dependence of the bulk modulus is small over the range of temperatures and pressures thought to prevail above the 400-kilometer discontinuity, and thus the data can be extrapolated to estimate the properties of olivine under mantle conditions at a depth of 400 kilometers. In the absence of high-temperature data at high pressures, estimates are made of the properties of olivine under mantle conditions to a depth of 400 kilometers. In contrast with low-pressure laboratory data, the predicted covariance of shear and compressional velocities as a function of temperature nearly matches the seismically estimated value for the lower mantle.     

Melting of (Mg, Fe)SiO3-Perovskite to 625 Kilobars: Indication of a High Melting Temperature in the Lower Mantle

The melting curves of two compositions of (Mg,Fe) SiO3-perovskite, the likely dominant mineral phase in the lower mantle, have been measured in a C02 laser-heated diamond cell with direct temperature measurements and in situ detection of melting. At 625 kilobars, the melting temperature is 5000 +/- 200 kelvin, independent of composition. Extrapolation to the core-mantle boundary pressure of 1.35 megabar with three different melting relations yields melting temperatures between 7000 and 8500 kelvin. Thus, the temperature at the base of the lower mantle, accepted to lie between 2550 and 2750 kelvin, is only at about one-third of the melting temperature. The large difference between mantle temperature and corresponding melting temperature has several important implications; particularly the temperature sensitivity of the viscosity is reduced thus allowing large lateral temperature variations inferred from seismic tomographic velocity anomalies and systematics found in measured velocity-density functions. Extensive melting of the lower mantle can be ruled out throughout the history of the Earth.     

Post-perovskite phase transition in MgSiO3

In situ x-ray diffraction measurements of MgSiO3 were performed at high pressure and temperature similar to the conditions at Earth's core-mantle boundary. Results demonstrate that MgSiO3 perovskite transforms to a new high-pressure form with stacked SiO6-octahedral sheet structure above 125 gigapascals and 2500 kelvin (2700-kilometer depth near the base of the mantle) with an increase in density of 1.0 to 1.2%. The origin of the D" seismic discontinuity may be attributed to this post-perovskite phase transition. The new phase may have large elastic anisotropy and develop preferred orientation with platy crystal shape in the shear flow that can cause strong seismic anisotropy below the D" discontinuity.     

First-principles determination of elastic anisotropy and wave velocities of MgO at lower mantle conditions

The individual elastic constants of magnesium oxide (MgO) have been determined throughout Earth's lower mantle (LM) pressure-temperature regime with density functional perturbation theory. It is shown that temperature effects on seismic observables (density, velocities, and anisotropy) are monotonically suppressed with increasing pressure. Therefore, at realistic LM conditions, the isotropic wave velocities of MgO remain comparable to seismic velocities, as previously noticed in athermal high-pressure calculations. Also, the predicted strong pressure-induced anisotropy is preserved toward the bottom of the LM, so lattice-preferred orientations in MgO may contribute substantially to the observed seismic anisotropy in the D" layer.     

Earth's Core-Mantle Boundary: Results of Experiments at High Pressures and Temperatures

Laboratory experiments document that liquid iron reacts chemically with silicates at high pressures (>/=2.4 x 10(10) Pascals) and temperatures. In particular, (Mg,Fe)SiO(3) perovskite, the most abundant mineral of Earth's lower mantle, is expected to react with liquid iron to produce metallic alloys (FeO and FeSi) and nonmetallic silicates (SiO(2) stishovite and MgSiO(3) perovskite) at the pressures of the core-mantle boundary, 14 x 10(10) Pascals. The experimental observations, in conjunction with seismological data, suggest that the lowermost 200 to 300 kilometers of Earth's mantle, the D" layer, may be an extremely heterogeneous region as a result of chemical reactions between the silicate mantle and the liquid iron alloy of Earth's core. The combined thermal-chemical-electrical boundary layer resulting from such reactions offers a plausible explanation for the complex behavior of seismic waves near the core-mantle boundary and could influence Earth's magnetic field observed at the surface.     

Premonitory Variations in S-Wave Velocity Anisotropy before Earthquakes in Nevada

Application of nonhydrostatic stress to rock induces velocity anisotropy, causing the S wave to split into two components traveling with somewhat different velocities. Large premonitory changes in the extent of S-wave splitting have been observed for two earthquakes in Nevada. Observations of the difference between the two S-wave velocities may provide a simple method for pedicting earthquakes.     

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.     

Magnetic field reversals, polar wander, and core-mantle coupling

True polar wander, the shifting of the entire mantle relative to the earth's spin axis, has been reanalyzed. Over the last 200 million years, true polar wander has been fast (approximately 5 centimeters per year) most of the time, except for a remarkable standstill from 170 to 110 million years ago. This standstill correlates with a decrease in the reversal frequency of the geomagnetic field and episodes of continental breakup. Conversely, true polar wander is high when reversal frequency increases. It is proposed that intermittent convection modulates the thickness of a thermal boundary layer at the base of the mantle and consequently the core-to-mantle heat flux. Emission of hot thermals from the boundary layer leads to increases in mantle convection and true polar wander. In conjunction, cold thermals released from a boundary layer at the top of the liquid core eventually lead to reversals. Changes in the locations of subduction zones may also affect true polar wander. Exceptional volcanism and mass extinctions at the Cretaceous-Tertiary and Permo-Triassic boundaries may be related to thermals released after two unusually long periods with no magnetic reversals. These environmental catastrophes may therefore be a consequence of thermal and chemical couplings in the earth's multilayer heat engine rather than have an extraterrestrial cause.