Discovery of flat seismic reflections in the mantle beneath the young Juan de Fuca Plate

Crustal properties of young oceanic lithosphere have been examined extensively, but the nature of the mantle lithosphere underneath remains elusive. Using a novel wide-angle seismic imaging technique, here we show the presence of two sub-horizontal reflections at ∼11 and ∼14.5 km below the seafloor over the 0.51-2.67 Ma old Juan de Fuca Plate. We find that the observed reflectors originate from 300-600-m-thick layers, with an ∼7-8% drop in P-wave velocity. They could be explained either by the presence of partially molten sills or frozen gabbroic sills. If partially molten, the shallower sill would define the base of a thin lithosphere with the constant thickness (11 km), requiring the presence of a mantle thermal anomaly extending up to 2.67 Ma. In contrast, if these reflections were frozen melt sills, they would imply the presence of thick young oceanic lithosphere (20-25 km), and extremely heterogeneous upper mantle.

Predictions and Observations for the Oceanic Lithosphere From S-to-P Receiver Functions and SS Precursors

The ocean lithosphere is classically described by the thermal half-space cooling (HSC) or the plate models, both characterized by a gradual transition to the asthenosphere beneath. Scattered waves find sharp seismic discontinuities beneath the oceans, possibly from the base of the plate. Active source studies suggest sharp discontinuities from a melt channel. We calculate synthetic S-to-P receiver functions and SS precursors for the HSC and plate models and also for channels. We find that the HSC and plate model velocity gradients are too gradual to create interpretable scattered waves from the base of the plate. Subtle phases are predicted to follow a similar trend as observations, flattening at older ages. Therefore, the seismic discontinuities are probably caused by a thermally controlled process that can also explain their amplitude, such as melting. Melt may coalesce in channels, although channels >10 km thick should be resolvable by scattered wave imaging.

Scattered wave imaging of the oceanic plate in Cascadia

Fifty years after plate tectonic theory was developed, the defining mechanism of the plate is still widely debated. The relatively short, simple history of young ocean lithosphere makes it an ideal place to determine the property that defines a plate, yet the remoteness and harshness of the seafloor have made precise imaging challenging. We use S-to-P receiver functions to image discontinuities beneath newly formed lithosphere at the Juan de Fuca and Gorda Ridges. We image a strong negative discontinuity at the base of the plate increasing from 20 to 45 km depth beneath the 0- to 10-million-year-old seafloor and a positive discontinuity at the onset of melting at 90 to 130 km depth. Comparison with geodynamic models and experimental constraints indicates that the observed discontinuities cannot easily be reconciled with subsolidus mechanisms. Instead, partial melt may be required, which would decrease mantle viscosity and define the young oceanic plate.

Evidence for frozen melts in the mid-lithosphere detected from active-source seismic data

The interactions of the lithospheric plates that form the Earth's outer shell provide much of the evidentiary basis for modern plate tectonic theory. Seismic discontinuities in the lithosphere arising from mantle convection and plate motion provide constraints on the physical and chemical properties of the mantle that contribute to the processes of formation and evolution of tectonic plates. Seismological studies during the past two decades have detected seismic discontinuities within the oceanic lithosphere in addition to that at the lithosphere-asthenosphere boundary (LAB). However, the depth, distribution, and physical properties of these discontinuities are not well constrained, which makes it difficult to use seismological data to examine their origin. Here we present new active-source seismic data acquired along a 1,130 km profile across an old Pacific plate (148-128 Ma) that show oceanic mid-lithosphere discontinuities (oceanic MLDs) distributed 37-59 km below the seafloor. The presence of the oceanic MLDs suggests that frozen melts that accumulated at past LABs have been preserved as low-velocity layers within the current mature lithosphere. These observations show that long-offset, high-frequency, active-source seismic data can be used to image mid-lithospheric structure, which is fundamental to understanding the formation and evolution of tectonic plates.

Constraints on the anisotropic contributions to velocity discontinuities at ∼60 km depth beneath the Pacific

Strong, sharp, negative seismic discontinuities, velocity decreases with depth, are observed beneath the Pacific seafloor at ∼60 km depth. It has been suggested that these are caused by an increase in radial anisotropy with depth, which occurs in global surface wave models. Here we test this hypothesis in two ways. We evaluate whether an increase in surface wave radial anisotropy with depth is robust with synthetic resolution tests. We do this by fitting an example surface wave data set near the East Pacific Rise. We also estimate the apparent isotropic seismic velocity discontinuities that could be caused by changes in radial anisotropy in S-to-P and P-to-S receiver functions and SS precursors using synthetic seismograms. We test one model where radial anisotropy is caused by olivine alignment and one model where it is caused by compositional layering. The result of our surface wave inversion suggests strong shallow azimuthal anisotropy beneath 0-10 Ma seafloor, which would also have a radial anisotropy signature. An increase in radial anisotropy with depth at 60 km depth is not well-resolved in surface wave models, and could be artificially observed. Shallow isotropy underlain by strong radial anisotropy could explain moderate apparent velocity drops (<6%) in SS precursor imaging, but not receiver functions. The effect is diminished if strong anisotropy also exists at 0-60 km depth as suggested by surface waves. Overall, an increase in radial anisotropy with depth may not exist at 60 km beneath the oceans and does not explain the scattered wave observations.

Upper-mantle water stratification inferred from observations of the 2012 Indian Ocean earthquake

Water, the most abundant volatile in Earth's interior, preserves the young surface of our planet by catalysing mantle convection, lubricating plate tectonics and feeding arc volcanism. Since planetary accretion, water has been exchanged between the hydrosphere and the geosphere, but its depth distribution in the mantle remains elusive. Water drastically reduces the strength of olivine and this effect can be exploited to estimate the water content of olivine from the mechanical response of the asthenosphere to stress perturbations such as the ones following large earthquakes. Here, we exploit the sensitivity to water of the strength of olivine, the weakest and most abundant mineral in the upper mantle, and observations of the exceptionally large (moment magnitude 8.6) 2012 Indian Ocean earthquake to constrain the stratification of water content in the upper mantle. Taking into account a wide range of temperature conditions and the transient creep of olivine, we explain the transient deformation in the aftermath of the earthquake that was recorded by continuous geodetic stations along Sumatra as the result of water- and stress-activated creep of olivine. This implies a minimum water content of about 0.01 per cent by weight-or 1,600 H atoms per million Si atoms-in the asthenosphere (the part of the upper mantle below the lithosphere). The earthquake ruptured conjugate faults down to great depths, compatible with dry olivine in the oceanic lithosphere. We attribute the steep rheological contrast to dehydration across the lithosphere-asthenosphere boundary, presumably by buoyant melt migration to form the oceanic crust.

Changes in seismic anisotropy shed light on the nature of the Gutenberg discontinuity

The boundary between the lithosphere and asthenosphere is associated with a platewide high-seismic velocity "lid" overlying lowered velocities, consistent with thermal models. Seismic body waves also intermittently detect a sharp velocity reduction at similar depths, the Gutenberg (G) discontinuity, which cannot be explained by temperature alone. We compared an anisotropic tomography model with detections of the G to evaluate their context and relation to the lithosphere-asthenosphere boundary (LAB). We find that the G is primarily associated with vertical changes in azimuthal anisotropy and lies above a thermally controlled LAB, implying that the two are not equivalent interfaces. The origin of the G is a result of frozen-in lithospheric structures, regional compositional variations of the mantle, or dynamically perturbed LAB.

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.

The Earth’s ‘hum’ is driven by ocean waves over the continental shelves

Observations show that the seismic normal modes of the Earth at frequencies near 10 mHz are excited at a nearly constant level in the absence of large earthquakes. This background level of excitation has been called the 'hum' of the Earth, and is equivalent to the maximum excitation from a magnitude 5.75 earthquake. Its origin is debated, with most studies attributing the forcing to atmospheric turbulence, analogous to the forcing of solar oscillations by solar turbulence. Some reports also predicted that turbulence might excite the planetary modes of Mars to detectable levels. Recent observations on Earth, however, suggest that the predominant excitation source lies under the oceans. Here I show that turbulence is a very weak source, and instead it is interacting ocean waves over the shallow continental shelves that drive the hum of the Earth. Ocean waves couple into seismic waves through the quadratic nonlinearity of the surface boundary condition, which couples pairs of slowly propagating ocean waves of similar frequency to a high phase velocity component at approximately double the frequency. This is the process by which ocean waves generate the well known 'microseism peak' that dominates the seismic spectrum near 140 mHz (refs 11, 12), but at hum frequencies, the mechanism differs significantly in frequency and depth dependence. A calculation of the coupling between ocean waves and seismic modes reproduces the seismic spectrum observed. Measurements of the temporal correlation between ocean wave data and seismic data have confirmed that ocean waves, rather than atmospheric turbulence, are driving the modes of the Earth.

Geophysical evidence from the MELT area for compositional controls on oceanic plates

Magnetotelluric and seismic data, collected during the MELT experiment at the southern East Pacific Rise, constrain the distribution of melt beneath this mid-ocean-ridge spreading centre and also the evolution of the oceanic lithosphere during its early cooling history. Here we focus on structures imaged at distances approximately 100 to 350 km east of the ridge crest, corresponding to seafloor ages of approximately 1.3 to 4.5 million years (Myr), where the seismic and electrical conductivity structure is nearly constant and independent of age. Beginning at a depth of about 60 km, we image a large increase in electrical conductivity and a change from isotropic to transversely anisotropic electrical structure, with higher conductivity in the direction of fast propagation for seismic waves. Conductive cooling models predict structure that increases in depth with age, extending to about 30 km at 4.5 Myr ago. We infer, however, that the structure of young oceanic plates is instead controlled by a decrease in water content above a depth of 60 km induced by the melting process beneath the spreading centre.

Pressure sensitivity of olivine slip systems and seismic anisotropy of Earth's upper mantle

The mineral olivine dominates the composition of the Earth's upper mantle and hence controls its mechanical behaviour and seismic anisotropy. Experiments at high temperature and moderate pressure, and extensive data on naturally deformed mantle rocks, have led to the conclusion that olivine at upper-mantle conditions deforms essentially by dislocation creep with dominant [100] slip. The resulting crystal preferred orientation has been used extensively to explain the strong seismic anisotropy observed down to 250 km depth. The rapid decrease of anisotropy below this depth has been interpreted as marking the transition from dislocation to diffusion creep in the upper mantle. But new high-pressure experiments suggest that dislocation creep also dominates in the lower part of the upper mantle, but with a different slip direction. Here we show that this high-pressure dislocation creep produces crystal preferred orientations resulting in extremely low seismic anisotropy, consistent with seismological observations below 250 km depth. These results raise new questions about the mechanical state of the lower part of the upper mantle and its coupling with layers both above and below.

Seismic evidence for hotspot-induced buoyant flow beneath the Reykjanes Ridge

Volcanic hotspots and mid-ocean ridge spreading centers are the surface expressions of upwelling in Earth's mantle convection system, and their interaction provides unique information on upwelling dynamics. I investigated the influence of the Iceland hotspot on the adjacent mid-Atlantic spreading center using phase-delay times of seismic surface waves, which show anomalous polarization anisotropy-a delay-time discrepancy between waves with different polarizations. This anisotropy implies that the hotspot induces buoyancy-driven upwelling in the mantle beneath the ridge.

Mapping the Hawaiian plume conduit with converted seismic waves

The volcanic edifice of the Hawaiian islands and seamounts, as well as the surrounding area of shallow sea floor known as the Hawaiian swell, are believed to result from the passage of the oceanic lithosphere over a mantle hotspot. Although geochemical and gravity observations indicate the existence of a mantle thermal plume beneath Hawaii, no direct seismic evidence for such a plume in the upper mantle has yet been found. Here we present an analysis of compressional-to-shear (P-to-S) converted seismic phases, recorded on seismograph stations on the Hawaiian islands, that indicate a zone of very low shear-wave velocity (< 4 km s(-1)) starting at 130-140 km depth beneath the central part of the island of Hawaii and extending deeper into the upper mantle. We also find that the upper-mantle transition zone (410-660 km depth) appears to be thinned by up to 40-50 km to the south-southwest of the island of Hawaii. We interpret these observations as localized effects of the Hawaiian plume conduit in the asthenosphere and mantle transition zone with excess temperature of approximately 300 degrees C. Large variations in the transition-zone thickness suggest a lower-mantle origin of the Hawaiian plume similar to the Iceland plume, but our results indicate a 100 degrees C higher temperature for the Hawaiian plume.

Mantle values of thermal conductivity and the geotherm from phonon lifetimes

A model for thermal conductivity kappa, based on phonon lifetimes obtained from infrared reflectivity, replicates experimental data at ambient conditions. The pressure and absolute temperature dependences of transport properties are accurately obtained from the Gruneisen parameter gammaTh, bulk modulus KT, and thermal expansivity alpha: The lattice contribution kappalat equals kappa298(298/T)a exp[-(4gammaTh + 1/3) integral298Talpha(theta)dtheta] with a = 0.33 for silicates (or 0.9 for MgO), and partial differential[ln(kappalat)]/ partial differentialP = (1/3 + 4gammaTh)/KT. The smaller, pressure-independent radiative contribution kapparad equals 0.0175 - 0.0001037T + (2.245T2/10(7)) - (3.407T3/10(11)), in units of watts per meter-kelvin, if Fe2+ is present. The resulting lithospheric geotherm is steep. Consequently, the mantle geotherm is hot if the low-velocity zone is anhydrous, but cold if hydrated.

Structure of the upper mantle under the EPR from waveform inversion of regional events

Waveform inversions of seismograms recorded at the Mantle Electromagnetic and Tomography (MELT) Experiment ocean bottom seismometer array from regional events with paths following the East Pacific Rise (EPR) require that low shear velocities (<3.7 km/s) extend to depths of more than 100 km below the rise axis. Velocities increase with average crustal age along ray paths. The reconciliation of Love and Rayleigh wave data requires that shear flow has aligned melt pockets or olivine crystals, creating an anisotropic uppermost mantle.