Understanding the intricate thermal dynamics within oceanic subduction zones is pivotal for deciphering the profound role these regions play in Earth’s material and energy cycles. As tectonic plates converge, the subducted oceanic plate interacts extensively with the overlying mantle wedge, producing complex patterns of heat exchange that profoundly influence fluid migration, melt generation, and seismic behavior. Central to this interplay are the variations in temperature that govern mineral dehydration and rheological transitions within the subduction environment, ultimately dictating the occurrence, style, and distribution of earthquakes.
Oceanic subduction zones present a multifaceted seismic regime, characterized by diverse earthquake types across distinct structural domains. The overriding plate frequently experiences thrust-fault earthquakes in the brittle upper crust, driven by compressive stresses. At the plate interface, destructive megathrust events and slower seismic phenomena—such as episodic tremor and slow slip events—manifest, revealing a spectrum of fault slip behaviors. Within the subducting slab itself, seismicity spans regions from the outer rise to depths approaching the lower mantle, exhibiting a wide array of focal mechanisms and spatial complexities that vary significantly among subduction zones worldwide.
Fundamentally, the thermal structure of subduction zones serves as the primary controller of seismic characteristics. Temperature modulates mineral dehydration reactions, which release fluids into surrounding rocks, and it determines the brittle-to-ductile transition depths that influence fault mechanics. According to recent geodynamic modeling, the convergence rate of subducting plates is the dominant parameter shaping subduction zone temperature profiles. Younger, hotter slabs exhibit elevated surface temperatures, while slower convergence rates further raise slab surface temperatures at shallower depths beneath approximately 70 kilometers. This dynamic thermal regime creates a maximum depth of decoupling (MDD) typically between 70 and 80 kilometers, beyond which the mechanical coupling between the slab and mantle wedge transitions from partial to full.
Above the MDD, the cool subducting slab induces a “cold corner” within the mantle wedge forearc, where reduced temperatures suppress ductile deformation. Below this decoupling depth, viscous coupling prevails, with the slab and mantle wedge moving in concert. This results in enhanced corner flow patterns that increase heat transfer, warming the slab surface significantly. Intriguingly, this MDD depth coincides with the maximum exhumation depth of metamorphic rocks found along the oceanic subduction interface, highlighting a potential linkage between thermal-mechanical transitions and tectonic exhumation processes. The underlying mechanisms governing the MDD remain subjects of ongoing investigation, complicated by variability in subduction geometry and thermal evolution.
The subduction plate interface, commonly termed the subduction channel, comprises several structurally distinct but interrelated zones: the roof décollement, basal décollement, and an intervening deformation zone. Earthquake activity within partially locked segments arises primarily due to strain localization within weaker lithologies or at contacts between contrasting rock types. Such complex deformation results in exhumed rocks from the subduction channel reflecting a heterogeneous mixture of sources, including subducted slab material, forearc crustal rocks, and mantle wedge peridotites and pyroxenites. The intricate mechanical coupling along the interface is modulated by both temperature and lithological composition, influencing seismic behavior and fault dynamics.
Seismic coupling at the interface varies with depth. From the trench down to approximately 40–50 kilometers, megathrust earthquakes dominate, reflecting decoupling conditions that promote brittle failure and seismic slip. At greater depths, typically beyond 70–80 kilometers, viscous coupling related to similar rheological strengths between the slab interface and mantle wedge impedes brittle failure, thereby localizing seismicity within the subducting slab. Transitional zones exhibit partial coupling and serve as sites for complex deformation, heterogeneous fluid activity, and variable fault locking. Here, both short-term brittle deformation events and longer-term ductile processes can coexist, challenging simplistic models of subduction seismogenesis.
Fluid dynamics within subduction systems critically influence seismicity patterns. Global earthquake statistics reveal a marked decrease in earthquake frequency with increasing depth, reaching a minimum near 300 kilometers. In colder subduction zones, dehydration of hydrous minerals predominantly occurs between 80 and 200 kilometers beneath volcanic arcs, with complete dehydration of minerals like lawsonite and phengite not achieved until depths approaching 300 kilometers. Warmer subduction zones reach complete dehydration at shallower depths, typically less than 160 kilometers. These observations underscore dehydration embrittlement as the primary mechanism driving intermediate-depth earthquakes across diverse thermal regimes.
Besides dehydration embrittlement, additional processes may contribute to intermediate-depth seismicity. Thermal runaway instabilities, eclogitization-related embrittlement, and metamorphism-facilitated instabilities in minerals such as orthopyroxene provide alternative or complementary explanations. Slow earthquakes, including episodic tremor and slow slip events, are often localized in subduction regions characterized by low effective stress and elevated pore fluid pressures, conditions conducive to transient aseismic slip. Such elevated fluid pressures likely arise from ongoing dehydration reactions of multiple hydrous mineral phases, modulating fault friction and slip behavior.
Water transport beyond subarc depths involves not only hydrous minerals but also nominally anhydrous phases and dense hydrous magnesium silicates stable in cold slabs. These minerals can convey water into the mantle transition zone, potentially enhancing localized hydration and influencing deep Earth processes. Deep-focus earthquakes—occurring below 300 kilometers—are generally attributed to transformational faulting mechanisms within metastable olivine. However, the influence of fluids on these deep seismic events remains unresolved, adding complexity to the understanding of deep Earth seismicity.
Despite advances in modeling and observational seismology, significant uncertainties persist concerning the interplay among metamorphism, seismicity, and fluid or melt activity in subduction contexts. The temporal evolution of thermal structures, variations in subduction parameters, and lithologic heterogeneities generate a rich but complex tectonic milieu that challenges current theoretical frameworks. Future research integrating high-fidelity experimental phase equilibrium studies, open-system thermodynamics, precise earthquake relocation methods, and geological investigations of fossil subduction zones promises to refine models of subduction dynamics.
Harnessing these multidisciplinary approaches will enhance our comprehension of subduction zone evolution, with implications for seismic hazard assessment and prediction. Improved delineation of thermal and mechanical boundaries within subduction systems will clarify the controls on earthquake nucleation and slow slip behavior. Ultimately, this knowledge contributes to safer, more informed management of populations living atop some of the most geologically active and seismically volatile areas on the planet.
—
Subject of Research: Thermal structure, fluid migration, and seismicity in oceanic subduction zones
Article Title: Thermal Controls on Fluids and Earthquakes in Oceanic Subduction Zones
Web References: http://dx.doi.org/10.1007/s11430-024-1514-4
References: Based on the modeling results of Peacock and Wang (2021) and literature review
Image Credits: ©Science China Press
Keywords: subduction zones, thermal structure, earthquakes, fluid migration, mineral dehydration, slab coupling, seismicity, mantle wedge, slow slip events, intermediate-depth earthquakes
