Tomography reveals buoyant asthenosphere accumulating beneath the Juan de Fuca plate

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.

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.

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.

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.

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.