Seismic evidence for sharp lithosphere-asthenosphere boundaries of oceanic plates

Seismic evidence for melt-rich lithosphere-asthenosphere boundary beneath young slab at Cascadia

The Lithosphere-Asthenosphere Boundary (LAB) beneath oceanic plates is generally imaged as a sharp seismic velocity reduction, suggesting the presence of partial melts. However, the fate of a melt-rich LAB is unclear after these plates descend into the mantle at subduction zones. Recent geophysical studies suggest its persistence with down-going old and cold slabs, but whether or not it is commonly present remains unclear, especially for young and warm slabs such as in the Cascadia subduction zone. Here we provide evidence for its presence at Cascadia in the form of a large (9.8 ± 1.5 % ) decrease in shear-wave velocity over a very small (<3 km) depth interval. Similarly large and sharp seismic velocity reduction at the bottom of both old and young slabs, as well as along the base of oceanic plates before subduction, possibly represents widespread presence of melts. The melt-rich sub-slab LAB may strongly influence subduction dynamics and viscoelastic earthquake cycles.

Genesis of Hawaiian lavas by crystallization of picritic magma in the deep mantle

Olivine is the dominant phenocryst or xenocryst of Hawaiian tholeiitic basalts, and the general consensus is that lavas with MgO concentrations from 7.5 to about 15 weight percent were derived from their primary magmas, which contain ~18-20 weight percent MgO, by only olivine crystallization. However, the major element composition of estimated primary magmas through olivine crystallization correction is inconsistent with direct partial melting of either mantle peridotite or its hybrid with subducted oceanic crust. Our melting experiments on peridotite-derived melt composition show that this discrepancy can be resolved if the primary magmas experienced two other processes before abundant olivine fractionation. First, the primary magmas experienced crystallization of clinopyroxene and garnet in the chamber at the base of the lithosphere (approximately the depths of 90-100 km). Second, the evolved magmas re-equilibrated with harzburgite when passing through the lithospheric mantle (approximately the depths of 60-10 km). Different from the isotopic evidence, the major and rare earth element compositions of Hawaiian post-shield alkali basalts and shield tholeiites suggest that they form from the same source by assimilating different amounts of orthopyroxene.

Wide-angle seismic reflections reveal a lithosphere-asthenosphere boundary zone in the subducting Pacific Plate, New Zealand

New wide-angle seismic reflection data from offshore New Zealand show that the lithosphere-asthenosphere boundary (LAB) is more structured than previously thought. Three distinct layers are interpreted within a 10- to 12-km-thick LAB zone beginning at a depth of ≈70 km: a 3 (±1)-km-thick layer at the bottom of the lithosphere with a P-wave (V_P) azimuthal anisotropy of 14 to 17% and fast azimuth subparallel to the direction of absolute plate motion and a 9 (±2)-km-thick, low V_P channel with a P-wave-to-S-wave velocity ratio (V_P/V_S) of >2.8 in the upper 7 km of the channel and 1.8 to 2.6 in the lower 2 km of the channel. The high V_P/V_S ratios indicate that this channel may contain 3 to 20% partial melt that facilitates decoupling of the lithosphere from the asthenosphere and reduces resistance for plate motion. Furthermore, the strong azimuthal anisotropy above the low-velocity layer suggests localization of strain due to melt accumulation.

On the Relationship Between Oceanic Plate Speed, Tectonic Stress, and Seismic Anisotropy

Seismic radial anisotropy (the squared ratio between the speeds of horizontally and vertically polarized shear waves,ξ = V S H 2 V S V 2 ) is a powerful tool to probe the direction of mantle flow and accumulated strain. While previous studies have confirmed the dependence of azimuthal anisotropy on plate speed, the first order control on radial anisotropy is unclear. In this study, we develop 2D ridge flow models combined with mantle fabric calculations to report that faster plates generate higher tectonics stresses and strain rates which lower the dislocation creep viscosity and lead to deeper anisotropy than beneath slower plates. We apply the SGLOBE-rani tomographic filter, resulting in a flat depth-age trend and stronger anisotropy beneath faster plates, which correlates well with 3D global anisotropic mantle models. Our predictions and observations suggest that as plate speed increases from 2 to 8 cm/yr, radial anisotropy increases by ∼0.01-0.025 in the upper 100-200 km of the mantle between 10 and 60 Ma.

Discovery of distinct lithosphere-asthenosphere boundary and the Gutenberg discontinuity in the Atlantic Ocean

The plate tectonic theory requires a rigid lithosphere floating over a weak asthenosphere, separated by the lithosphere-asthenosphere boundary, which has been sometimes interpreted as the Gutenberg discontinuity. Using a deep seismic reflection technique, we report the presence of two continuous reflections covering 27 Ma to 58 Ma oceanic lithosphere in the Atlantic Ocean. We find that the upper reflection deepens with age and follows the ~1250°C isotherm, whereas the deeper reflection lies at a constant depth of ~75 km. We suggest that the upper reflection represents the thermally controlled lithosphere-asthenosphere boundary, whereas the lower reflection is the Gutenberg discontinuity, a frozen-in dehydration boundary separating the dry mantle melting region above from the hydrated mantle below formed at the ridge axis. We also find that thermal mantle anomalies rejuvenate the lithosphere, uplift the lithosphere-asthenosphere boundary, and destroy the Gutenberg discontinuity.

Electrical conductivity of melts: implications for conductivity anomalies in the Earth's mantle

Magmatic liquids, including silicate and carbonate melts, are principal agents of mass and heat transfer in the Earth and terrestrial planets, and they play a crucial role in various geodynamic processes and in Earth's evolution. Electrical conductivity data of these melts elucidate the cause of electrical anomalies in Earth's interior and shed light on the melt structure. With the improvement in high-pressure experimental techniques and theoretical simulations, major progress has been made on this front in the past several decades. This review aims to summarize recent advances in experimental and theoretical studies on the electrical conductivity of silicate and carbonate melts of different compositions and volatile contents under high temperature and pressure. The electrical conductivity of silicate melts depends strongly on temperature, pressure, water content and the ratio of non-bridging oxygens to tetrahedral cations (NBO/T). By contrast, the electrical conductivity of carbonate melts exhibits a weak dependence on temperature and pressure due to their fully depolymerized structure. The electrical conductivity of carbonate melts is higher than that of silicate melts by at least two orders of magnitude. Water can increase electrical conductivity significantly and reduce the activation energy of silicate melts. Conversely, this effect is weak for carbonate melts. In addition, the replacement of alkali-earth elements (Ca²⁺ or Mg²⁺) with alkali elements causes a significant decrease in the electrical conductivity of carbonate melts. A distinct compensation trend is revealed for the electrical conductivity of silicate and carbonate melts under anhydrous and hydrous conditions. Several important applications of laboratory-based melt conductivity are introduced in order to understand the origin of high-conductivity anomalies in the Earth's mantle. Perspectives for future studies are also provided.

Dynamic slab segmentation due to brittle-ductile damage in the outer rise

Subduction is the major plate driving force, and the strength of the subducting plate controls many aspects of the thermochemical evolution of Earth. Each subducting plate experiences intense normal faulting^1-9 during bending that accommodates the transition from horizontal to downwards motion at the outer rise at trenches. Here we investigate the consequences of this bending-induced plate damage using numerical subduction models in which both brittle and ductile deformation, including grain damage, are tracked and coupled self-consistently. Pervasive slab weakening and pronounced segmentation can occur at the outer-rise region owing to the strong feedback between brittle and ductile damage localization. This slab-damage phenomenon explains the subduction dichotomy of strong plates and weak slabs¹⁰, the development of large-offset normal faults^6,7 near trenches, the occurrence of segmented seismic velocity anomalies¹¹ and distinct interfaces imaged within subducted slabs^12,13, and the appearance of deep, localized intraplate areas of reduced effective viscosity¹⁴ observed at trenches. Furthermore, brittle-viscously damaged slabs show a tendency for detachment at elevated mantle temperatures. Given Earth's planetary cooling history¹⁵, this implies that intermittent subduction with frequent slab break-off episodes¹⁶ may have been characteristic for Earth until more recent times than previously suggested¹⁷.

Liquid structure under extreme conditions: high-pressure x-ray diffraction studies

Under extreme conditions of high pressure and temperature, liquids can undergo substantial structural transformations as their atoms rearrange to minimise energy within a more confined volume. Understanding the structural response of liquids under extreme conditions is important across a variety of disciplines, from fundamental physics and exotic chemistry to materials and planetary science.In situexperiments and atomistic simulations can provide crucial insight into the nature of liquid-liquid phase transitions and the complex phase diagrams and melting relations of high-pressure materials. Structural changes in natural magmas at the high-pressures experienced in deep planetary interiors can have a profound impact on their physical properties, knowledge of which is important to inform geochemical models of magmatic processes. Generating the extreme conditions required to melt samples at high-pressure, whilst simultaneously measuring their liquid structure, is a considerable challenge. The measurement, analysis, and interpretation of structural data is further complicated by the inherent disordered nature of liquids at the atomic-scale. However, recent advances in high-pressure technology mean that liquid diffraction measurements are becoming more routinely feasible at synchrotron facilities around the world. This topical review examines methods for high pressure synchrotron x-ray diffraction of liquids and the wide variety of systems which have been studied by them, from simple liquid metals and their remarkable complex behaviour at high-pressure, to molecular-polymeric liquid-liquid transitions in pnicogen and chalcogen liquids, and density-driven structural transformations in water and silicate melts.

Stress-induced amorphization triggers deformation in the lithospheric mantle

The mechanical properties of olivine-rich rocks are key to determining the mechanical coupling between Earth's lithosphere and asthenosphere. In crystalline materials, the motion of crystal defects is fundamental to plastic flow^1-4. However, because the main constituent of olivine-rich rocks does not have enough slip systems, additional deformation mechanisms are needed to satisfy strain conditions. Experimental studies have suggested a non-Newtonian, grain-size-sensitive mechanism in olivine involving grain-boundary sliding^5,6. However, very few microstructural investigations have been conducted on grain-boundary sliding, and there is no consensus on whether a single or multiple physical mechanisms are at play. Most importantly, there are no theoretical frameworks for incorporating the mechanics of grain boundaries in polycrystalline plasticity models. Here we identify a mechanism for deformation at grain boundaries in olivine-rich rocks. We show that, in forsterite, amorphization takes place at grain boundaries under stress and that the onset of ductility of olivine-rich rocks is due to the activation of grain-boundary mobility in these amorphous layers. This mechanism could trigger plastic processes in the deep Earth, where high-stress conditions are encountered (for example, at the brittle-plastic transition). Our proposed mechanism is especially relevant at the lithosphere-asthenosphere boundary, where olivine reaches the glass transition temperature, triggering a decrease in its viscosity and thus promoting grain-boundary sliding.

Seismic evidence for partial melt below tectonic plates

The seismic low-velocity zone (LVZ) of the upper mantle is generally associated with a low-viscosity asthenosphere that has a key role in decoupling tectonic plates from the mantle¹. However, the origin of the LVZ remains unclear. Some studies attribute its low seismic velocities to a small amount of partial melt of minerals in the mantle^2,3, whereas others attribute them to solid-state mechanisms near the solidus^4-6 or the effect of its volatile contents⁶. Observations of shear attenuation provide additional constraints on the origin of the LVZ⁷. On the basis of the interpretation of global three-dimensional shear attenuation and velocity models, here we report partial melt occurring within the LVZ. We observe that partial melting down to 150-200 kilometres beneath mid-ocean ridges, major hotspots and back-arc regions feeds the asthenosphere. A small part of this melt (less than 0.30 per cent) remains trapped within the oceanic LVZ. Melt is mostly absent under continental regions. The amount of melt increases with plate velocity, increasing substantially for plate velocities of between 3 centimetres per year and 5 centimetres per year. This finding is consistent with previous observations of mantle crystal alignment underneath tectonic plates⁸. Our observations suggest that by reducing viscosity⁹ melt facilitates plate motion and large-scale crystal alignment in the asthenosphere.

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.

High-pressure elastic properties of dolomite melt supporting carbonate-induced melting in deep upper mantle

Petrologic studies suggest that carbonate-rich melts are present in the Earth’s upper mantle and play an important role in the deep carbon cycle. However, seismic detection of these melts is difficult due to the lack of understanding of the elastic properties of carbonate melts. Here we determined the sound velocity and density of dolomite melt at upper mantle conditions using high-pressure experiments and theoretical simulations. The calculated velocities of carbonate melt-bearing mantle using these new elasticity data were compared with global seismic observations. The result suggests that ∼0.05% carbonate-rich melt can be pervasively present in the deep upper mantle, implying a global average mantle carbon concentration of 80-140 ppm.

Deeply subducted carbonates likely cause low-degree melting of the upper mantle and thus play an important role in the deep carbon cycle. However, direct seismic detection of carbonate-induced partial melts in the Earth’s interior is hindered by our poor knowledge on the elastic properties of carbonate melts. Here we report the first experimentally determined sound velocity and density data on dolomite melt up to 5.9 GPa and 2046 K by in-situ ultrasonic and sink-float techniques, respectively, as well as first-principles molecular dynamics simulations of dolomite melt up to 16 GPa and 3000 K. Using our new elasticity data, the calculated V_P/V_S ratio of the deep upper mantle (∼180–330 km) with a small amount of carbonate-rich melt provides a natural explanation for the elevated V_P/V_S ratio of the upper mantle from global seismic observations, supporting the pervasive presence of a low-degree carbonate-rich partial melt (∼0.05%) that is consistent with the volatile-induced or redox-regulated initial melting in the upper mantle as argued by petrologic studies. This carbonate-rich partial melt region implies a global average carbon (C) concentration of 80–140 ppm. by weight in the deep upper mantle source region, consistent with the mantle carbon content determined from geochemical studies.

Melting of recycled ancient crust responsible for the Gutenberg discontinuity

A discontinuity in the seismic velocity associated with the lithosphere-asthenosphere interface, known as the Gutenberg discontinuity, is enigmatic in its origin. While partial mantle melts are frequently suggested to explain this discontinuity, it is not well known which factors critically regulate the melt production. Here, we report geochemical evidence showing that the melt fractions in the lithosphere-asthenosphere boundary were enhanced not only by accumulation of compacted carbonated melts related to recycled ancient marine sediments, but also by partial melting of a pyroxene-rich mantle domain related to the recycled oceanic eclogite/pyroxenites. This conclusion is derived from the first set of Mg isotope data for a suite of young petit-spot basalts erupted on the northwest Pacific plate, where a clearly defined Gutenberg discontinuity exists. Our results reveal a specific linkage between the Gutenberg discontinuity beneath the normal oceanic regions and the recycling of ancient subducted crust and carbonate through the deep Earth.

Acoustic and Microstructural Properties of Partially Molten Samples in the Ice–Ammonia System

We measured the ultrasonic properties and microstructure of two-phase binary mixtures of the ice–ammonia partial melt system, which was selected based on its importance for numerous planetary bodies. The equilibrium microstructure of ice–ammonia melt was examined using a light microscope within a cold room. The measured median dihedral angle between the solid and melt at 256 K is approximately 63°, with a broad distribution of observed angles between 10° and 130°. P-wave velocities in the partially molten samples were measured as a function of temperature (177 < T(K) < 268) and composition (1–6.4 wt % NH3). Vp decreases approximately linearly with increasing temperature and melt fraction. We compare the results of this study to those of other potential binary systems by normalizing the datasets using a vertical lever (TL–TE) and articulating the potential effects on the mechanical behavior and transport capabilities of partially molten ice in icy satellites.

Rapid mantle flow with power-law creep explains deformation after the 2011 Tohoku mega-quake

The deformation transient following large subduction zone earthquakes is thought to originate from the interaction of viscoelastic flow in the asthenospheric mantle and slip on the megathrust that are both accelerated by the sudden coseismic stress change. Here, we show that combining insight from laboratory solid-state creep and friction experiments can successfully explain the spatial distribution of surface deformation in the first few years after the 2011 M_w 9.0 Tohoku-Oki earthquake. The transient reduction of effective viscosity resulting from dislocation creep in the asthenosphere explains the peculiar retrograde displacement revealed by seafloor geodesy, while the slip acceleration on the megathrust accounts for surface displacements on land and offshore outside the rupture area. Our results suggest that a rapid mantle flow takes place in the asthenosphere with temporarily decreased viscosity in response to large coseismic stress, presumably due to the activation of power-law creep during the post-earthquake period.

Imaging Lithospheric Discontinuities Beneath the Northern East African Rift Using S-to-P Receiver Functions

Imaging the lithosphere is key to understand mechanisms of extension as rifting progresses. Continental rifting results in a combination of mechanical stretching and thinning of the lithosphere, decompression upwelling, heating, sometimes partial melting of the asthenosphere, and potentially partial melting of the mantle lithosphere. The northern East African Rift system is an ideal locale to study these processes as it exposes the transition from tectonically active continental rifting to incipient seafloor spreading. Here we use S-to-P receiver functions to image the lithospheric structure beneath the northernmost East African Rift system where it forms a triple junction between the Main Ethiopian rift, the Red Sea rift, and the Gulf of Aden rift. We image the Moho at 31 ± 6 km beneath the Ethiopian plateau. The crust is 28 ± 3 km thick beneath the Main Ethiopian rift and thins to 23 ± 2 km in northern Afar. We identify a negative phase, a velocity decrease with depth, at 67 ± 3 km depth beneath the Ethiopian plateau, likely associated with the lithosphere-asthenosphere boundary (LAB), and a lack of a LAB phase beneath the rift. Using observations and waveform modeling, we show that the LAB phase beneath the plateau is likely defined by a small amount of partial melt. The lack of a LAB phase beneath the rift suggests melt percolation through the base of the lithosphere beneath the northernmost East African Rift system.

Insights on Upper Mantle Melting, Rheology, and Anelastic Behavior From Seismic Shear Wave Tomography

In seismic tomography we map the wave speed structure inside the Earth, but we ultimately seek to interpret those images in terms of physical parameters. This is challenging because many parameters can trade-off with each other to produce a given wave speed. The problem is compounded by the convention of mapping seismic structures as perturbations relative to a 1-D reference model, rather than absolute wave speeds. Using a full waveform tomography model of Europe as a case study, we quantify the extent to which thermochemical and dynamic properties can be constrained using only S wave speed, expressed in absolute values. The wave speed distributions of this tomography model are compared with 4 million thermochemical models, whose seismic properties are computed via thermodynamic modeling. These models sample the full range of realistic mantle compositions, including variable water and melt contents, and mineral intrinsic anelasticity is taken into account. Intrinsic anelasticity causes waves to travel more slowly at higher temperatures, leading to seismic attenuation, but the sensitivity of the wave speed reduction to temperature is, in turn, controlled by the wave frequency. Global studies of surface waves indicate an anticorrelation between S wave speed and attenuation. We therefore only retain thermochemical models satisfying this anticorrelation. Our study indicates that the frequency dependence of anelasticity, α, depends on temperature or rheology, with α ≈ 0.1 being most appropriate in cold or lithospheric mantle and α ≈ 0.3 in warmer regions (i.e., the asthenosphere). Additionally, the slowest regions require specific compositions and/or a velocity-weakening mechanism, such as partial melting, elastically accommodated grain boundary sliding, or water.

Fossilized Melts in Mantle Wedge Peridotites

The shallow oxidized asthenosphere may contain a small fraction of potassic silicate melts that are enriched in incompatible trace elements and volatiles. Here, to determine the chemical composition of such melt, we analysed fossilized melt inclusions, preserved as multiphase solid inclusions, from an orogenic garnet peridotite in the Bohemian Massif. Garnet-poor (2 vol.%) peridotite preserves inclusions of carbonated potassic silicate melt within Zn-poor chromite (<400 ppm) in the clinopyroxene-free harzburgite assemblage that equilibrated within the hot mantle wedge (Stage 1, > 1180 °C at 3 GPa). The carbonated potassic silicate melt, which has a major element oxide chemical composition of K₂O = 5.2 wt.%, CaO = 17 wt.%, MgO = 18 wt.%, CO₂ = 22 wt.%, and SiO₂ = 20 wt.%, contains extremely high concentrations of large ion lithophile elements, similar to kimberlite melts. Peridotites cooled down to ≅800 °C during Stage 2, resulted in the growth of garnet relatively poor in pyrope content, molar Mg/(Mg + Fe + Ca + Mn), (ca. 67 mol.%). This garnet displays a sinusoidal REE pattern that formed in equilibrium with carbonatitic fluid. Subsequently, subduction of the peridotite resulted in the formation of garnet with a slightly higher pyrope content (70 mol.%) during the Variscan subduction Stage 3 (950 °C, 2.9 GPa). These data suggest the following scenario for the generation of melt in the mantle wedge. Primarily, infiltration of sediment-derived potassic carbonatite melt into the deep mantle wedge resulted in the growth of phlogopite and carbonate/diamond. Formation of volatile-bearing minerals lowered the density and strength of the peridotite. Finally, phlogopite-bearing carbonated peridotite rose as diapirs in the mantle wedge to form carbonated potassic silicate melts at the base of the overriding lithosphere.

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.

Seismic Imaging of Thickened Lithosphere Resulting From Plume Pulsing Beneath Iceland

Ocean plates conductively cool and subside with seafloor age. Plate thickening with age is also predicted, and hot spots may cause thinning. However, both are debated and depend on the way the plate is defined. Determining the thickness of the plates along with the process that governs it has proven challenging. We use S-to-P (Sp) receiver functions to image a strong, persistent LAB beneath Iceland where the mid-Atlantic Ridge interacts with a plume with hypothesized pulsating thermal anomaly. The plate is thickest, up to 84 ± 6 km, beneath lithosphere formed during times of hypothesized hotter plume temperatures and as thin as 61 ± 6 km beneath regions formed during colder intervals. We performed geodynamic modeling to show that these plate thicknesses are inconsistent with a thermal lithosphere. Instead, periods of increased plume temperatures likely increased the melting depth, causing deeper depletion and dehydration, and creating a thicker plate. This suggests plate thickness is dictated by the conditions of plate formation.

Three-dimensional variations of the slab geometry correlate with earthquake distributions at the Cascadia subduction system

Significant along-strike variations of seismicity are observed at subduction zones, which are strongly influenced by physical properties of the plate interface and rheology of the crust and mantle lithosphere. However, the role of the oceanic side of the plate boundary on seismicity is poorly understood due to the lack of offshore instrumentations. Here tomographic results of the Cascadia subduction system, resolved with full-wave ambient noise simulation and inversion by integrating dense offshore and onshore seismic datasets, show significant variations of the oceanic lithosphere along strike and down dip from spreading centers to subduction. In central Cascadia, where seismicity is sparse, the slab is imaged as a large-scale low-velocity feature near the trench, which is attributed to a highly hydrated and strained oceanic lithosphere underlain by a layer of melts or fluids. The strong correlation suggests that the properties of the incoming oceanic plate play a significant role on seismicity.

Variations in seismicity are observed at subduction zones, but the oceanic sides remain poorly resolved. Here, the author presents tomographic results of the Cascadia subduction system demonstrating that there are significant variations of the oceanic lithosphere along the subduction zone.

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.

Determination of intrinsic attenuation in the oceanic lithosphere-asthenosphere system

We recorded P and S waves traveling through the oceanic lithosphere-asthenosphere system (LAS) using broadband ocean-bottom seismometers in the northwest Pacific, and we quantitatively separated the intrinsic (anelastic) and extrinsic (scattering) attenuation effects on seismic wave propagation to directly infer the thermomechanical properties of the oceanic LAS. The strong intrinsic attenuation in the asthenosphere obtained at higher frequency (~3 hertz) is comparable to that constrained at lower frequency (~100 seconds) by surface waves and suggests frequency-independent anelasticity, whereas the intrinsic attenuation in the lithosphere is frequency dependent. This difference in frequency dependence indicates that the strong and broad peak dissipation recently observed in the laboratory exists only in the asthenosphere and provides new insight into what distinguishes the asthenosphere from the lithosphere.

Sporadic low-velocity volumes spatially correlate with shallow very low frequency earthquake clusters

A low-velocity zone (LVZ) has been detected by seismic exploration surveys within the Nankai accretionary prism toe off the Kii Peninsula, southwestern Japan, and is considered to be a mechanically weak volume at depth. Such mechanical heterogeneities potentially influence seismic and tsunamigenic slips on megathrust earthquakes in the subduction zone. However, the spatial distribution of the LVZ along the trough-parallel direction is still elusive. Here we show sporadic LVZs in the prism toe from one-dimensional shear wave velocity (Vs) profiles obtained at 49 cabled ocean bottom stations, which were estimated by a nonlinear inversion technique, simulated annealing, using the displacement–pressure ratios of the Rayleigh wave. The estimated distribution of LVZs along the trough widely correlates with the epicentres of shallow very low frequency earthquakes (sVLFEs), which suggests that sVLFEs are activated in the sporadically distributed low-velocity and mechanically weak volumes where fluids significantly reduce the shear strength of faults.

In the Nankai subduction zone a low-velocity zone (LVZ) has been detected. Here, the authors present shear wave velocity profiles to show that low frequency earthquakes correlate with the distribution of the LVZ in the Nankai subduction zone.

Elastic properties of silicate melts: Implications for low velocity zones at the lithosphere-asthenosphere boundary

Low seismic velocity regions in the mantle and crust are commonly attributed to the presence of silicate melts. Determining melt volume and geometric distribution is fundamental to understanding planetary dynamics. We present a new model for seismic velocity reductions that accounts for the anomalous compressibility of silicate melt, rendering compressional wave velocities more sensitive to melt fraction and distribution than previous estimates. Forward modeling predicts comparable velocity reductions for compressional and shear waves for partially molten mantle, and for low velocity regions associated with the lithosphere-asthenosphere boundary (LAB), melt present at <5% distributed in near-textural equilibrium. These findings reconcile seismic observations for the LAB regionally and locally and favor models of strong coupling across the LAB rather than melt channeling due to shear deformation.

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.

Layering of subcontinental lithospheric mantle

Recent seismic studies reveal a sharp velocity drop mostly at ∼70-100km depth within the thick mantle keel beneath cratons, termed the mid-lithosphere discontinuity (MLD). The common presence of the MLD in cratonic regions indicates structural and property layering of the subcontinental lithospheric mantle (SCLM). The nature and origin of the MLD, and many issues associated with the layering of the SCLM are essential to understand the formation and evolution of continents, and have become frontier subjects in the Earth sciences.

Petit-spot as definitive evidence for partial melting in the asthenosphere caused by CO2

The deep carbon cycle plays an important role on the chemical differentiation and physical properties of the Earth's mantle. Especially in the asthenosphere, seismic low-velocity and high electrical conductivity due to carbon dioxide (CO₂)-induced partial melting are expected but not directly observed. Here we discuss the experimental results relevant to the genesis of primitive CO₂-rich alkali magma forming petit-spot volcanoes at the deformation front of the outer rise of the northwestern Pacific plate. The results suggest that primitive melt last equilibrated with depleted peridotite at 1.8–2.1 GPa and 1,280–1,290 °C. Although the equilibration pressure corresponds to the pressure of the lower lithosphere, by considering an equilibration temperature higher than the solidus in the volatile–peridotite system along with the temperature of the lower lithosphere, we conclude that CO₂-rich silicate melt is always produced in the asthenosphere. The melt subsequently ascends into and equilibrates with the lower lithosphere before eruption.

The genesis of primitive CO₂-rich alkali magma at petit-spot volcanoes is still poorly understood. Here, the authors perform high pressure and temperature experiments to show that the CO₂-rich melt is derived from the asthenosphere.

Heterogeneity in mantle carbon content from CO2-undersaturated basalts

The amount of carbon present in Earth's mantle affects the dynamics of melting, volcanic eruption style and the evolution of Earth's atmosphere via planetary outgassing. Mantle carbon concentrations are difficult to quantify because most magmas are strongly degassed upon eruption. Here we report undegassed carbon concentrations from a new set of olivine-hosted melt inclusions from the Mid-Atlantic Ridge. We use the correlations of CO₂ with trace elements to define an average carbon abundance for the upper mantle. Our results indicate that the upper mantle carbon content is highly heterogeneous, varying by almost two orders of magnitude globally, with the potential to produce large geographic variations in melt fraction below the volatile-free solidus. Such heterogeneity will manifest as variations in the depths at which melt becomes interconnected and detectable, the CO₂ fluxes at mid-ocean ridges, the depth of the lithosphere-asthenosphere boundary, and mantle conductivity.

Tomography reveals buoyant asthenosphere accumulating beneath the Juan de Fuca plate

The boundary between Earth's strong lithospheric plates and the underlying mantle asthenosphere corresponds to an abrupt seismic velocity decrease and electrical conductivity increase with depth, perhaps indicating a thin, weak layer that may strongly influence plate motion dynamics. The behavior of such a layer at subduction zones remains unexplored. We present a tomographic model, derived from on- and offshore seismic experiments, that reveals a strong low-velocity feature beneath the subducting Juan de Fuca slab along the entire Cascadia subduction zone. Through simple geodynamic arguments, we propose that this low-velocity feature is the accumulation of material from a thin, weak, buoyant layer present beneath the entire oceanic lithosphere. The presence of this feature could have major implications for our understanding of the asthenosphere and subduction zone dynamics.

Olivine anisotropy suggests Gutenberg discontinuity is not the base of the lithosphere

Tectonic plates are a key feature of Earth's structure, and their behavior and dynamics are fundamental drivers in a wide range of large-scale processes. The operation of plate tectonics, in general, depends intimately on the manner in which lithospheric plates couple to the convecting interior. Current debate centers on whether the transition from rigid lithosphere to flowing asthenosphere relates to increases in temperature or to changes in composition such as the presence of a small amount of melt or an increase in water content below a specified depth. Thus, the manner in which the rigid lithosphere couples to the flowing asthenosphere is currently unclear. Here we present results from laboratory-based torsion experiments on olivine aggregates with and without melt, yielding an improved database describing the crystallographic alignment of olivine grains. We combine this database with a flow model for oceanic upper mantle to predict the structure of the seismic anisotropy beneath ocean basins. Agreement between our model and seismological observations supports the view that the base of the lithosphere is thermally controlled. This model additionally supports the idea that discontinuities in velocity and anisotropy, often assumed to be the base of the lithosphere, are, instead, intralithospheric features reflecting a compositional boundary established at midocean ridges, not a rheological boundary.

Satellite tidal magnetic signals constrain oceanic lithosphere-asthenosphere boundary

The tidal flow of electrically conductive oceans through the geomagnetic field results in the generation of secondary magnetic signals, which provide information on the subsurface structure. Data from the new generation of satellites were shown to contain magnetic signals due to tidal flow; however, there are no reports that these signals have been used to infer subsurface structure. We use satellite-detected tidal magnetic fields to image the global electrical structure of the oceanic lithosphere and upper mantle down to a depth of about 250 km. The model derived from more than 12 years of satellite data reveals a ≈72-km-thick upper resistive layer followed by a sharp increase in electrical conductivity likely associated with the lithosphere-asthenosphere boundary, which separates colder rigid oceanic plates from the ductile and hotter asthenosphere.

Experimental evidence supports mantle partial melting in the asthenosphere

The low-velocity zone (LVZ) is a persistent seismic feature in a broad range of geological contexts. It coincides in depth with the asthenosphere, a mantle region of lowered viscosity that may be essential to enabling plate motions. The LVZ has been proposed to originate from either partial melting or a change in the rheological properties of solid mantle minerals. The two scenarios imply drastically distinct physical and geochemical states, leading to fundamentally different conclusions on the dynamics of plate tectonics. We report in situ ultrasonic velocity measurements on a series of partially molten samples, composed of mixtures of olivine plus 0.1 to 4.0 volume % of basalt, under conditions relevant to the LVZ. Our measurements provide direct compressional (V P) and shear (V S) wave velocities and constrain attenuation as a function of melt fraction. Mantle partial melting appears to be a viable origin for the LVZ, for melt fractions as low as ~0.2%. In contrast, the presence of volatile elements appears necessary to explaining the extremely high V P/V S values observed in some local areas. The presence of melt in LVZ could play a major role in the dynamics of plate tectonics, favoring the decoupling of the plate relative to the asthenosphere.

Experimental test of the viscous anisotropy hypothesis for partially molten rocks

Chemical differentiation of rocky planets occurs by melt segregation away from the region of melting. The mechanics of this process, however, are complex and incompletely understood. In partially molten rocks undergoing shear deformation, melt pockets between grains align coherently in the stress field; it has been hypothesized that this anisotropy in microstructure creates an anisotropy in the viscosity of the aggregate. With the inclusion of anisotropic viscosity, continuum, two-phase-flow models reproduce the emergence and angle of melt-enriched bands that form in laboratory experiments. In the same theoretical context, these models also predict sample-scale melt migration due to a gradient in shear stress. Under torsional deformation, melt is expected to segregate radially inward. Here we present torsional deformation experiments on partially molten rocks that test this prediction. Microstructural analyses of the distribution of melt and solid reveal a radial gradient in melt fraction, with more melt toward the center of the cylinder. The extent of this radial melt segregation grows with progressive strain, consistent with theory. The agreement between theoretical prediction and experimental observation provides a validation of this theory.

Experimental constraints on the electrical anisotropy of the lithosphere-asthenosphere system

The relative motion of lithospheric plates and underlying mantle produces localized deformation near the lithosphere-asthenosphere boundary. The transition from rheologically stronger lithosphere to weaker asthenosphere may result from a small amount of melt or water in the asthenosphere, reducing viscosity. Either possibility may explain the seismic and electrical anomalies that extend to a depth of about 200 kilometres. However, the effect of melt on the physical properties of deformed materials at upper-mantle conditions remains poorly constrained. Here we present electrical anisotropy measurements at high temperatures and quasi-hydrostatic pressures of about three gigapascals on previously deformed olivine aggregates and sheared partially molten rocks. For all samples, electrical conductivity is highest when parallel to the direction of prior deformation. The conductivity of highly sheared olivine samples is ten times greater in the shear direction than for undeformed samples. At temperatures above 900 degrees Celsius, a deformed solid matrix with nearly isotropic melt distribution has an electrical anisotropy factor less than five. To obtain higher electrical anisotropy (up to a factor of 100), we propose an experimentally based model in which layers of sheared olivine are alternated with layers of sheared olivine plus MORB or of pure melt. Conductivities are up to 100 times greater in the shear direction than when perpendicular to the shear direction and reproduce stress-driven alignment of the melt. Our experimental results and the model reproduce mantle conductivity-depth profiles for melt-bearing geological contexts. The field data are best fitted by an electrically anisotropic asthenosphere overlain by an isotropic, high-conductivity lowermost lithosphere. The high conductivity could arise from partial melting associated with localized deformation resulting from differential plate velocities relative to the mantle, with subsequent upward melt percolation from the asthenosphere.

A seismic reflection image for the base of a tectonic plate

Plate tectonics successfully describes the surface of Earth as a mosaic of moving lithospheric plates. But it is not clear what happens at the base of the plates, the lithosphere-asthenosphere boundary (LAB). The LAB has been well imaged with converted teleseismic waves, whose 10-40-kilometre wavelength controls the structural resolution. Here we use explosion-generated seismic waves (of about 0.5-kilometre wavelength) to form a high-resolution image for the base of an oceanic plate that is subducting beneath North Island, New Zealand. Our 80-kilometre-wide image is based on P-wave reflections and shows an approximately 15° dipping, abrupt, seismic wave-speed transition (less than 1 kilometre thick) at a depth of about 100 kilometres. The boundary is parallel to the top of the plate and seismic attributes indicate a P-wave speed decrease of at least 8 ± 3 per cent across it. A parallel reflection event approximately 10 kilometres deeper shows that the decrease in P-wave speed is confined to a channel at the base of the plate, which we interpret as a sheared zone of ponded partial melts or volatiles. This is independent, high-resolution evidence for a low-viscosity channel at the LAB that decouples plates from mantle flow beneath, and allows plate tectonics to work.

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.

Melt-rich channel observed at the lithosphere-asthenosphere boundary

The lithosphere-asthenosphere boundary (LAB) separates rigid oceanic plates from the underlying warm ductile asthenosphere. Although a viscosity decrease beneath this boundary is essential for plate tectonics, a consensus on its origin remains elusive. Seismic studies identify a prominent velocity discontinuity at depths thought to coincide with the LAB but disagree on its cause, generally invoking either partial melting or a mantle dehydration boundary as explanations. Here we use sea-floor magnetotelluric data to image the electrical conductivity of the LAB beneath the edge of the Cocos plate at the Middle America trench offshore of Nicaragua. Underneath the resistive oceanic lithosphere, the magnetotelluric data reveal a high-conductivity layer confined to depths of 45 to 70 kilometres. Because partial melts are stable at these depths in a warm damp mantle, we interpret the conductor to be a partially molten layer capped by an impermeable frozen lid that is the base of the lithosphere. A conductivity anisotropy parallel to plate motion indicates that this melt has been sheared into flow-aligned tube-like structures. We infer that the LAB beneath young plates consists of a thin, partially molten, channel of low viscosity that acts to decouple the overlying brittle lithosphere from the deeper convecting mantle. Because this boundary layer has the potential to behave as a lubricant to plate motion, its proximity to the trench may have implications for subduction dynamics.

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.

Geophysics. At the bottom of the oceanic plate

High-resolution seismic maps are providing a better picture of the processes involved in plate tectonics.