Subduction

Subduction is the convergent tectonic process in which a cold, gravitationally unstable plate of oceanic lithosphere descends beneath an overriding oceanic or continental plate and sinks into the convective asthenospheric mantle. Operating along roughly 55,000 km of convergent margins—principally around the circum-Pacific "Ring of Fire," the Sunda–Java Trench of the eastern Indian Ocean, and the Lesser Antilles and Scotia arcs of the Atlantic—subduction recycles oceanic crust generated at divergent spreading centers back into the deep Earth at convergence rates of 2 to 11 cm yr⁻¹, balancing the global surface area budget of plate tectonics and closing the planetary rock cycle.

First synthesized during the late 1960s by Kiyoo Wadati, Hugo Benioff, Jack Oliver, Bryan Isacks, and Dan McKenzie (and named by André Amstutz in 1951), subduction zones represent Earth's primary engines of geodynamic deformation, continental crust generation, and volatile cycling. As the descending slab heats and compresses along a low-geothermal-gradient trajectory, it undergoes high-pressure, low-temperature phase transformations through the blueschist and eclogite metamorphic facies, releases hydrous fluids that trigger flux melting in the overlying mantle wedge to feed calc-alkaline volcanic arcs, and generates the planet's largest megathrust earthquakes and tsunamis.

Negative Thermal Buoyancy and Slab Pull

Newly accreted oceanic lithosphere at a mid-ocean ridge is hot, thin (~6–10 km), and buoyant relative to the underlying asthenosphere. As the plate moves away from the ridge crest over 20 to 200 Myr, it cools by conductive heat loss to the overlying ocean, thickens as the 1,300 °C isotherm deepens to 80–110 km, and contracts thermally according to the half-space cooling law. Once oceanic lithosphere exceeds an age of ~20–30 Myr, its bulk density (ρ_slab ≈ 3.25 to 3.35 g cm⁻³) surpasses that of the underlying asthenospheric peridotite (ρ_mantle ≈ 3.22 g cm⁻³), rendering the plate gravitationally unstable.

The downward gravitational body force acting on this cold, dense sinking lithosphere is termed slab pull (F_SP):

F_SP = ∫ g · Δρ(z) · dA ≈ g · Δρ · L · h

where g = 9.81 m s⁻², Δρ ≈ 40 to 80 kg m⁻³ is the integrated density contrast between the cold slab and ambient mantle, L is downdip slab length, and h is slab thickness. Typical slab pull forces reach 3 × 10¹³ N per meter of trench length—roughly an order of magnitude larger than ridge push (F_RP ≈ 2–3 × 10¹² N m⁻¹, generated by gravitational sliding off elevated mid-ocean ridges). Geodynamic torque-balance models of Forsyth and Uyeda (1975) demonstrated that plates attached to extensive subducting slabs (such as the Pacific, Nazca, and Cocos plates) move at 6 to 10 cm yr⁻¹, whereas plates lacking significant subducting slabs (such as the Eurasian, North American, and African plates) move at only 1 to 2.5 cm yr⁻¹, confirming that slab pull provides 75% to 85% of the net driving force of plate tectonics.

Eclogite Densification and Mantle Transition Zone Boundaries

Two solid-state mineralogical phase transformations amplify and modulate slab pull as the slab descends:

1. The Basalt/Gabbro → Eclogite Transformation (45–70 km depth): Hydrated oceanic crust composed of low-density plagioclase-bearing basalt and gabbro (ρ ≈ 2.95–3.00 g cm⁻³) passes through glaucophane–lawsonite blueschist into anhydrous pyrope-almandine garnet + jadeite-rich omphacite eclogite (ρ ≈ 3.45–3.60 g cm⁻³). This 15% crustal density jump eliminates any residual crustal buoyancy and locks in irreversible downward sinking.
2. The 410-km and 660-km Mantle Discontinuities: At ~410 km depth, upper-mantle olivine ((Mg,Fe)₂SiO₄, α-phase) transforms exothermically into the denser β-spinel polymorph wadsleyite with a positive Clapeyron slope (dP/dT ≈ +2.5 to +3.1 MPa K⁻¹). Because the interior of the subducting slab is 400 to 700 K colder than ambient mantle, the wadsleyite transition elevates upward by 30 to 60 km inside the cold slab, adding an extra downward buoyancy force. Conversely, at the 660-km discontinuity, ringwoodite (γ-spinel) dissociates endothermically into bridgmanite ((Mg,Fe)SiO₃ perovskite) + ferropericlase ((Mg,Fe)O) with a negative Clapeyron slope (dP/dT ≈ −1.3 to −2.0 MPa K⁻¹), combined with a 10- to 30-fold jump in lower-mantle viscosity. This endothermic barrier deflects the depression of the boundary downward inside the slab and creates buoyant resistance, causing some slabs (such as Izu–Bonin and Tonga) to stagnate horizontally at 660 km before avalanching into the lower mantle toward the core–mantle boundary.

Cross-Sectional Anatomy of a Convergent Margin

From the seaward incoming plate landward across the overriding plate, a mature subduction system comprises five structural domains:

Wadati–Benioff Zones and Megathrust Earthquakes

Subduction zones generate more than 90% of Earth's total seismic moment release. Along the shallow plate interface (depths of 10 to 50 km, temperatures of 150 °C to 350 °C), frictional locking along the megathrust décollement accumulates elastic strain over centuries before rupturing in great thrust earthquakes reaching moment magnitudes M_w 9.0 to 9.5—including the 1960 Valdivia, Chile (M_w 9.5), 1964 Prince William Sound, Alaska (M_w 9.2), 2004 Sumatra–Andaman (M_w 9.1–9.3), and 2011 Tōhoku-Oki, Japan (M_w 9.1) earthquakes.

Below the brittle megathrust interface, inclined planar bands of intermediate-focus (70–300 km) and deep-focus (300–680 km) earthquakes—termed Wadati–Benioff zones—trace the cold interior of the descending slab down to the 660-km boundary at dip angles ranging from ~10–15° ("flat-slab" subduction beneath central Peru and central Chile) to nearly 80° (Mariana Trench). Because confining pressures above 2 GPa prohibit ordinary frictional sliding, intermediate-depth earthquakes (70–300 km) are triggered by dehydration embrittlement, where pore pressures of H₂O released by serpentine (antigorite → forsterite + enstatite + H₂O) and lawsonite breakdown reduce effective normal stress across pre-existing faults. Deep-focus earthquakes (350–680 km, such as the 1994 M_w 8.2 Bolivia and 2013 M_w 8.3 Sea of Okhotsk events) occur within the cold metastable olivine wedge via transformational faulting (anticrack shear instabilities during rapid olivine → wadsleyite/ringwoodite transformation).

Volatile Dehydration and Mantle-Wedge Flux Melting

Contrary to early misconceptions, arc volcanoes do not form by frictional melting of the subducting oceanic plate itself (except in rare cases of very young, hot slab melting that produce high-Sr/Y adakites). Instead, as the hydrated oceanic crust and serpentinized lithospheric mantle descend to depths of 80 to 130 km beneath the volcanic front, increasing pressure and temperature destabilize hydrous minerals—chlorite, amphibole, lawsonite, phengite, and antigorite serpentine—expelling supercritical aqueous fluids rich in H₂O, CO₂, Cl, sulfur, and fluid-mobile large-ion lithophile elements (Rb, Ba, K, Sr, Pb, U) upward into the overlying asthenospheric mantle wedge.

While dry peridotite at 100 km depth (3.3 GPa) melts only above ~1,450–1,500 °C, the addition of 0.2 to 1.5 wt% H₂O cleaves polymerized Si–O–Si silicate bridges and depresses the peridotite solidus temperature to ~1,000–1,100 °C. Because corner flow in the convective mantle wedge circulates hot asthenospheric peridotite (~1,300 °C) into this hydrated zone, flux melting generates 10% to 25% hydrous basaltic magma. As these buoyant hydrous basalts ascend through the overriding crust, fractional crystallization and crustal assimilation differentiate them into volatile-rich calc-alkaline andesites, dacites, and rhyolites, building island arcs (Mariana, Aleutian, Tonga) and continental volcanic arcs (Andes, Cascades) as well as underlying granitic batholiths that build continental crust.

Why does continental crust resist subduction when colliding at convergent margins?

Whereas mature oceanic crust is thin (~7 km), mafic, and transforms into dense eclogite (ρ ≈ 3.5 g cm⁻³) at depths above 50 km, continental crust is thick (35–70 km) and composed predominantly of buoyant quartz- and feldspar-rich granitoid rocks (ρ ≈ 2.70–2.80 g cm⁻³). When subduction consumes an entire ocean basin and brings a continent into a trench (as in the India–Asia collision forming the Himalayas), the positive buoyancy of the continental lithosphere resists deep mantle descent, causing slab break-off of the leading oceanic root and collisional mountain building (orogeny).

How can earthquakes occur at 600 km depth where mantle rocks normally flow ductilely?

Ambient mantle at 600 km depth has a temperature near 1,600 °C and deforms by ductile creep. However, because rock has low thermal conductivity, the core of a rapidly subducting oceanic slab remains 500 to 800 °C colder than the surrounding mantle even at 600 km depth. Inside this cold core (<650 °C), low-pressure olivine persists metastably tens of kilometers past the equilibrium 410-km boundary before abruptly collapsing into fine-grained wadsleyite or ringwoodite along runaway shear zones, radiating seismic waves as deep-focus earthquakes.

Why are volcanic arcs consistently located 100 to 120 km above the subducting slab?

Across global subduction zones with wildly different convergence rates and dip angles, the volcanic front sits where the top of the subducting slab reaches a depth of 105 ± 15 km. This depth marks the pressure-temperature threshold where hydrous minerals in the slab (amphibole, chlorite, and serpentine) break down simultaneously and where the slab decouples mechanically from the stagnant forearc mantle to encounter hot (>1,250 °C) convecting asthenosphere in the mantle wedge, crossing the wet peridotite solidus.

1. Forsyth, D., & Uyeda, S. (1975). "On the Relative Importance of the Driving Forces of Plate Motion." Geophysical Journal of the Royal Astronomical Society, 43(1), 163–200. https://doi.org/10.1111/j.1365-246X.1975.tb00631.x
2. Stern, R. J. (2002). "Subduction zones." Reviews of Geophysics, 40(4), 1012. https://doi.org/10.1029/2001RG000108
3. Hacker, B. R., Peacock, S. M., Abers, G. A., & Holloway, S. D. (2003). "Subduction factory 2. Are intermediate-depth earthquakes in subducting slabs linked to metamorphic dehydration reactions?" Journal of Geophysical Research: Solid Earth, 108(B1), 2030. https://doi.org/10.1029/2001JB001129
4. Grove, T. L., Till, C. B., & Krawczynski, M. J. (2012). "The Role of H₂O in Subduction Zone Magmatism." Annual Review of Earth and Planetary Sciences, 40, 413–439. https://doi.org/10.1146/annurev-earth-042711-105310