The Gutenberg discontinuity: melt at the lithosphere-asthenosphere boundary

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.

Destruction and regrowth of lithospheric mantle beneath large igneous provinces

Large igneous provinces (LIPs) are formed by enormous (i.e., frequently >106 km3) but short-lived magmatic events that have profound effects upon global geodynamic, tectonic, and environmental processes. Lithospheric structure is known to modulate mantle melting, yet its evolution during and after such dramatic periods of magmatism is poorly constrained. Using geochemical and seismological observations, we find that magmatism is associated with thin (i.e., ≲80 km) lithosphere and we reveal a striking positive correlation between the thickness of modern-day lithosphere beneath LIPs and time since eruption. Oceanic lithosphere rethickens to 125 km, while continental regions reach >190 km. Our results point to systematic destruction and subsequent regrowth of lithospheric mantle during and after LIP emplacement and recratonization of the continents following eruption. These insights have implications for the stability, age, and composition of ancient, thick, and chemically distinct lithospheric roots, the distribution of economic resources, and emissions of chemical species that force catastrophic environmental change.

Role of subduction dynamics on the unevenly distributed volcanism at the Middle American subduction system

A typical subduction of an oceanic plate beneath a continent is expected to be accompanied by arc volcanoes along the convergent margin. However, subduction of the Cocos plate at the Middle American subduction system has resulted in an uneven distribution of magmatism/volcanism along strike. Here we construct a new three-dimensional shear-wave velocity model of the entire Middle American subduction system, using full-wave ambient noise tomography. Our model reveals significant variations of the oceanic plates along strike and down dip, in correspondence with either weakened or broken slabs after subduction. The northern and southern segments of the Cocos plate, including the Mexican flat slab subduction, are well imaged as high-velocity features, where a low-velocity mantle wedge exists and demonstrate a strong correlation with the arc volcanoes. Subduction of the central Cocos plate encounters a thick high-velocity feature beneath North America, which hinders the formation of a typical low-velocity mantle wedge and arc volcanoes. We suggest that the presence of slab tearing at both edges of the Mexican flat slab has been modifying the mantle flows, resulting in the unusual arc volcanism.

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.

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.

Density of NaAlSi2O6 Melt at High Pressure and Temperature Measured by In-Situ X-ray Microtomography

In this study, the volumetric compression of jadeite (NaAlSi2O6) melt at high pressures was determined by three-dimensional volume imaging using the synchrotron-based X-ray microtomography technique in a rotation-anvil device. Combined with the sample mass, measured using a high-precision analytical balance prior to the high-pressure experiment, the density of jadeite melt was obtained at high pressures and high temperatures up to 4.8 GPa and 1955 K. The density data were fitted to a third-order Birch-Murnaghan equation of state, resulting in a best-fit isothermal bulk modulus K T 0 of 10.8 − 5.3 + 1.9 GPa and its pressure derivative K T 0 ′ of 3.4 − 0.4 + 6.6 . Comparison with data for silicate melts of various compositions from the literature shows that alkali-rich, polymerized melts are generally more compressible than alkali-poor, depolymerized ones. The high compressibility of jadeite melt at high pressures implies that polymerized sodium aluminosilicate melts, if generated by low-degree partial melting of mantle peridotite at ~250–400 km depth in the deep upper mantle, are likely denser than surrounding mantle materials, and thus gravitationally stable.

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.

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.

A subduction and mantle plume origin for Samoan volcanism

The origin of Samoan volcanism in the southwest Pacific remains enigmatic. Whether mantle melting is solely caused by a mantle plume is questionable because some volcanism, here referred to as non-hotspot volcanism, defies the plume model and its linear age-progression trend. Indeed, non-hotspot volcanism occurred as far as 740 km west of the predicted Samoan hotspot after 5 Ma. Here we use fully-dynamic laboratory subduction models and a tectonic reconstruction to show that the nearby Tonga-Kermadec-Hikurangi (TKH) subduction zone induces a broad mantle upwelling around the northern slab edge that coincides with the non-hotspot volcanic activity after 5 Ma. Using published potential mantle temperatures for the ambient mantle and Samoan mantle plume, we find that two geodynamic processes can explain mantle melting responsible for intraplate volcanism in the Samoan region. We propose that before 5 Ma, the volcanism is consistent with the plume model, whereas afterwards non-hotspot volcanism resulted from interaction between the Subduction-Induced Mantle Upwelling (SIMU) and Samoan mantle plume material that propagated west from the hotspot due to the toroidal component of slab rollback-induced mantle flow. In this geodynamic scenario, the SIMU drives decompression melting in the westward-swept plume material, thus producing the non-hotpot volcanism.

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.

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.

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.

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.

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.

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.

Electrical conductivity during incipient melting in the oceanic low-velocity zone

A low viscosity layer in the upper mantle, the Asthenosphere, is a requirement for plate tectonics¹. The seismic low velocities and the high electrical conductivities of the Asthenosphere are attributed either to sub-solidus water-related defects in olivine minerals^2-4 or to a few volume percents of partial melt^5-8 but these two interpretations have shortcomings: (1) The amount of H₂O stored in olivine is not expected to be higher than 50 ppm due to partitioning with other mantle phases⁹, including pargasite amphibole at moderate temperatures¹⁰, and partial melting at high temperatures⁹; (2) elevated melt volume fractions are impeded by the too cold temperatures prevailing in the Asthenosphere and by the high melt mobility that can lead to gravitational segregation^11,12. Here we determined the electrical conductivity of CO₂-H₂O-rich melts, typically produced at the onset of mantle melting. Electrical conductivity modestly increases with moderate amounts of H₂O and CO₂ but it dramatically increases as CO₂ content exceeds 6 wt% in the melt. Incipient melts, long-expected to prevail in the asthenosphere^10,13-15, can therefore trigger its high electrical conductivities. Considering depleted and enriched mantle abundances in H₂O and CO₂ and their effect on the petrology of incipient melting, we calculated conductivity profiles across the Asthenosphere for various plate ages. Several electrical discontinuities are predicted and match geophysical observations in a consistent petrological and geochemical framework. In moderately aged plates (>5Ma), incipient melts most likely trigger both the seismic low velocities and the high electrical conductivities in the upper part of the asthenosphere, whereas for young plates⁴, where seamount volcanism occurs⁶, higher degree of melting is expected.

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.

Geophysics. At the bottom of the oceanic plate

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