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.     

Time scales and heterogeneous structure in geodynamic earth models

Computer models of mantle convection constrained by the history of Cenozoic and Mesozoic plate motions explain some deep-mantle structural heterogeneity imaged by seismic tomography, especially those related to subduction. They also reveal a 150-million-year time scale for generating thermal heterogeneity in the mantle, comparable to the record of plate motion reconstructions, so that the problem of unknown initial conditions can be overcome. The pattern of lowermost mantle structure at the core-mantle boundary is controlled by subduction history, although seismic tomography reveals intense large-scale hot (low-velocity) upwelling features not explicitly predicted by the models.     

Melt segregation and strain partitioning: implications for seismic anisotropy and mantle flow

One of the principal means of understanding upper mantle dynamics involves inferring mantle flow directions from seismic anisotropy under the assumption that the seismic fast direction (olivine a axis) parallels the regional flow direction. We demonstrate that (i) the presence of melt weakens the alignment of a axes and (ii) when melt segregates and forms networks of weak shear zones, strain partitions between weak and strong zones, resulting in an alignment of a axes 90 degrees from the shear direction in three-dimensional deformation. This orientation of a axes provides a new means of interpreting mantle flow from seismic anisotropy in partially molten deforming regions of Earth.     

Earth's Viscosity

Seismic methods are now being used to determine not only Earth's elastic properties, but also by how much it departs from a perfectlyelastic body. The seismic anelasticity (Q) varies by several orders of magnitude throughout the mantle, the main feature being an extremely dissipative zone in the upper mantle above 400 kilometers. Recent determinations of viscosity by McConnell show a similar trend. The two sets of data indicate that the ratio of viscosity to Q is roughly a constant, at least in the upper mantle of Earth. On the assumption that this relation is valid for the rest of Earth, viscosities are estimated in regions that are inaccessible for direct measurement. The implied presence of a low-viscosity zone in the upper mantle, overlying a more viscous, less deformable, lower mantle, reconciles viscosites calculated from the shape of Earth and from postglacial uplift. The mismatch of the deformational characteristics at various levels in Earth, coupled with the changing rate of rotation, may be pertinent to the rate of release of seismic energy as a function of depth.     

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.     

Plate tectonics and hotspots: the third dimension

High-resolution seismic tomographic models of the upper mantle provide powerful new constraints on theories of plate tectonics and hotspots. Midocean ridges have extremely low seismic velocities to a depth of 100 kilometers. These low velocities imply partial melting. At greater depths, low-velocity and high-velocity anomalies record, respectively, previous positions of migrating ridges and trenches. Extensional, rifting, and hotspot regions have deep (> 200 kilometers) low-velocity anomalies. The upper mantle is characterized by vast domains of high temperature rather than small regions surrounding hotspots; the asthenosphere is not homogeneous or isothermal. Extensive magmatism requires a combination of hot upper mantle and suitable lithospheric conditions. High-velocity regions of the upper 200 kilometers of the mantle correlate with Archean cratons.     

Seismic Tomogram of the Earth's Mantle: Geodynamic Implications

Recent seismic tomography of the Earth's mantle has revealed a large-scale pattern of mantle convection comprising upwelling columnar plumes in the Pacific and Africa and downwelling planar sheets along the Circum Pacific. Upwelling and downwelling occur most extensively under the south Pacific and west Pacific, respectively. High-resolution image of plate subduction has been obtained from the dense seismic networks around Japan. Japanese seismologists are in the best position to resolve the internal structure of downwelling current as an integral part of the whole convection system.     

Seismic observations of splitting of the mid-transition zone discontinuity in Earth's mantle

The transition zone of Earth's mantle is delineated by globally observed discontinuities in seismic properties at depths of about 410 and 660 kilometers. Here, we investigate the detailed structure between 410 and 660 kilometers depth, by making use of regional stacks of precursors to the SS phase. The previously observed discontinuity at about 520 kilometers depth is confirmed in many regions, but is found to be absent in others. There are a number of regions in which we find two discontinuities at about 500 and 560 kilometers depth, an effect which can be interpreted as a "splitting" of the 520 kilometer discontinuity. These observations provide seismic constraints on the sharpness and observability of mineralogical phase transitions in the mantle transition zone.     

Water-induced fabric transitions in olivine

The interpretation of seismic anisotropy in Earth's upper mantle has traditionally been based on the fabrics (lattice-preferred orientation) of relatively water-poor olivine. Here we show that when a large amount of water is added to olivine, the relation between flow geometry and seismic anisotropy undergoes marked changes. Some of the puzzling observations of seismic anisotropy in the upper mantle, including the anomalous anisotropy in the central Pacific and the complicated anisotropy in subduction zones, can be attributed to the enrichment of water in these regions.     

Spin crossover in ferropericlase at high pressure: a seismologically transparent transition?

Seismic discontinuities in Earth typically arise from structural, chemical, or temperature variations with increasing depth. The pressure-induced iron spin state transition in the lower mantle may influence seismic wave velocities by changing the elasticity of iron-bearing minerals, but no seismological evidence of an anomaly exists. Inelastic x-ray scattering measurements on (Mg(0.83)Fe(0.17))O-ferropericlase at pressures across the spin transition show effects limited to the only shear moduli of the elastic tensor. This explains the absence of deviation in the aggregate seismic velocities and, thus, the lack of a one-dimensional seismic signature of the spin crossover. The spin state transition does, however, influence shear anisotropy of ferropericlase and should contribute to the seismic shear wave anisotropy of the lower mantle.     

Dipping low-velocity layer in the mid-lower mantle: evidence for geochemical heterogeneity

Data from western United States short-period seismic networks reveal a conversion from an S to a P wave within a low seismic velocity layer (greater than or equal to the 4 percent velocity difference compared to the surrounding mantle) in the mid-lower mantle (1400 to 1600 kilometers deep) east of the Mariana and Izu-Bonin subduction zones. The low-velocity layer (about 8 kilometers thick) dips 30 degrees to 40 degrees southward and is at least 500 kilometers by 300 kilometers. Its steep dip, large velocity contrast, and sharpness imply a chemical rather than a thermal origin. Ancient oceanic crust subducted into the lower mantle is a plausible candidate for the low-velocity layer because of its broad thin extent.     

Mantle convection and plate tectonics: toward an integrated physical and chemical theory

Plate tectonics and convection of the solid, rocky mantle are responsible for transporting heat out of Earth. However, the physics of plate tectonics is poorly understood; other planets do not exhibit it. Recent seismic evidence for convection and mixing throughout the mantle seems at odds with the chemical composition of erupted magmas requiring the presence of several chemically distinct reservoirs within the mantle. There has been rapid progress on these two problems, with the emergence of the first self-consistent models of plate tectonics and mantle convection, along with new geochemical models that may be consistent with seismic and dynamical constraints on mantle structure.     

Dislocation damping and anisotropic seismic wave attenuation in Earth's upper mantle

Crystal defects form during tectonic deformation and are reactivated by the shear stress associated with passing seismic waves. Although these defects, known as dislocations, potentially contribute to the attenuation of seismic waves in Earth's upper mantle, evidence for dislocation damping from laboratory studies has been circumstantial. We experimentally determined the shear modulus and associated strain-energy dissipation in pre-deformed synthetic olivine aggregates under high pressures and temperatures. Enhanced high-temperature background dissipation occurred in specimens pre-deformed by dislocation creep in either compression or torsion, the enhancement being greater for prior deformation in torsion. These observations suggest the possibility of anisotropic attenuation in relatively coarse-grained rocks where olivine is or was deformed at relatively high stress by dislocation creep in Earth's upper mantle.     

Continent-ocean chemical heterogeneity in the mantle based on seismic tomography

Seismic models of global-scale lateral heterogeneity in the mantle show systematic differences below continents and oceans that are too large to be purely thermal in origin. An inversion of the geoid, based on a seismic model that includes viscous flow in the mantle, indicates that the differences beneath continents and oceans can be accounted for by differences in composition in the upper mantle superposed on mantle-wide thermal heterogeneities. The net continent-ocean density differences, integrated over depth, are small and cause only a low flux of mass and heat across the asthenosphere and mantle transition zone.     

Evidence for a ubiquitous seismic discontinuity at the base of the mantle

A sharp discontinuity at the base of Earth's mantle has been suggested from seismic waveform studies; the observed travel time and amplitude variations have been interpreted as changes in the depth of a spatially intermittent discontinuity. Most of the observed variations in travel times and the spatial intermittance of the seismic triplication can be reproduced by a ubiquitous first-order discontinuity superimposed on global seismic velocity structure derived from tomography. The observations can be modeled by a solid-solid phase transition that has a 200-kilometer elevation above the core-mantle boundary under adiabatic temperatures and a Clapeyron slope of about 6 megapascal per kelvin.     

Seismic anisotropy: tracing plate dynamics in the mantle

Elastic anisotropy is present where the speed of a seismic wave depends on its direction. In Earth's mantle, elastic anisotropy is induced by minerals that are preferentially oriented in a directional flow or deformation. Earthquakes generate two seismic wave types: compressional (P) and shear (S) waves, whose coupling in anisotropic rocks leads to scattering, birefringence, and waves with hybrid polarizations. This varied behavior is helping geophysicists explore rock textures within Earth's mantle and crust, map present-day upper-mantle convection, and study the formation of lithospheric plates and the accretion of continents in Earth history.     

Mantle seismic structure beneath the MELT region of the east pacific rise from P and S wave tomography

Relative travel time delays of teleseismic P and S waves, recorded during the Mantle Electromagnetic and Tomography (MELT) Experiment, have been inverted tomographically for upper-mantle structure beneath the southern East Pacific Rise. A broad zone of low seismic velocities extends beneath the rise to depths of about 200 kilometers and is centered to the west of the spreading center. The magnitudes of the P and S wave anomalies require the presence of retained mantle melt; the melt fraction near the rise exceeds the fraction 300 kilometers off axis by as little as 1%. Seismic anisotropy, induced by mantle flow, is evident in the P wave delays at near-vertical incidence and is consistent with a half-width of mantle upwelling of about 100 km.     

Reconstructing Farallon plate subduction beneath North America back to the Late Cretaceous

Using an inverse mantle convection model that assimilates seismic structure and plate motions, we reconstruct Farallon plate subduction back to 100 million years ago. Models consistent with stratigraphy constrain the depth dependence of mantle viscosity and buoyancy, requiring that the Farallon slab was flat lying in the Late Cretaceous, consistent with geological reconstructions. The simulation predicts that an extensive zone of shallow-dipping subduction extended beyond the flat-lying slab farther east and north by up to 1000 kilometers. The limited region of flat subduction is consistent with the notion that subduction of an oceanic plateau caused the slab to flatten. The results imply that seismic images of the current mantle provide more constraints on past tectonic events than previously recognized.     

Seismic evidence for deep-water transportation in the mantle

We report seismic evidence for the transportation of water into the deep mantle in the subduction zone beneath northeastern Japan. Our data indicate that water is released from the hydrated oceanic crust at shallow depths (< approximately 100 kilometers) and then forms a channel of hydrated mantle material on top of the subducting plate that is the pathway for water into the deep mantle. Our result provides direct evidence that shows how water is transported from the ocean to the deep mantle in a cold subduction zone environment.     

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.