Subduction Zone Thermal Structure and Metamorphism

Core Idea

Subducting slabs remain cold owing to rapid plate motion. Cold slab interiors inhibit melting; thermal models show geothermal gradients much lower in subduction zones than in the mantle wedge, explaining metamorphic facies and magma generation.

Explainer

From your study of subduction zone dynamics, you know that oceanic lithosphere descends into the mantle at convergent boundaries. From crustal heat conduction models, you know that temperature distribution in the Earth is governed by the balance between heat sources, conduction, and advection. Subduction zone thermal structure brings these together: the descending slab carries cold oceanic lithosphere into the hot mantle, creating one of the most dramatic thermal contrasts anywhere in the Earth's interior.

The key to understanding why slabs stay cold is the competition between heat conduction and plate velocity. Heat conducts into the slab from the surrounding hot mantle, but the slab is moving downward faster than heat can diffuse inward. Think of sliding a frozen metal bar through a furnace — if you push it fast enough, the interior remains cold even as the surface heats up. The dimensionless number that captures this competition is the thermal Peclet number: the ratio of advective heat transport (plate motion) to conductive heat transport. For typical subduction rates of 5–10 cm/year, the Peclet number is large, meaning advection dominates and the slab interior stays cold to depths of several hundred kilometers.

The thermal structure of a subduction zone is not a simple temperature gradient — it has a distinctive two-dimensional pattern. The slab surface heats up progressively as it descends, reaching temperatures where hydrous minerals in the altered oceanic crust break down and release water. This dehydration produces a sequence of metamorphic facies along the slab surface: from zeolite and prehnite-pumpellyite facies at shallow depths, through blueschist facies (the hallmark of high-pressure, low-temperature conditions unique to subduction zones), to eclogite facies at greater depths where all hydrous phases have decomposed. The water released from the slab rises into the hot mantle wedge — the triangular region of mantle between the slab surface and the overlying plate — where it lowers the melting point of peridotite, triggering flux melting. This is the primary mechanism generating arc magmas, and it explains why volcanic arcs sit about 100–120 km above the slab surface: that is roughly the depth where the slab has heated enough to dehydrate its last major water-bearing minerals.

The thermal structure varies dramatically between subduction zones. Old, cold, fast-subducting slabs (like the Pacific plate beneath Japan) retain their cold cores to great depths and produce narrow, well-defined Wadati-Benioff seismic zones. Young, warm, slow-subducting slabs (like the Juan de Fuca plate beneath Cascadia) heat up more rapidly, dehydrate at shallower depths, and may lose their seismic signature before reaching 100 km depth. These differences have direct consequences for volcanic style, earthquake depth distribution, and the recycling of water and carbon into the deep mantle. Thermal modeling of subduction zones — solving the heat equation with realistic geometries, velocities, and rheologies — is therefore central to understanding why convergent margins behave so differently from one another around the globe.