Explore the vast range of titles available.

Despite the critical role of subduction in plate tectonics, the dynamics of its initiation remains unclear. High-temperature low-pressure metamorphic soles are vestiges of subduction initiation, providing records of the pressure and temperature conditions along the subducting slab surface during subduction initiation that can possibly differentiate the two end-member subduction initiation modes: spontaneous and induced. Here, using numerical models, we show that the slab surface temperature reaches 800–900 °C at ~1 GPa over a wide range of parameter values for spontaneous subduction initiation whereas for induced subduction initiation, such conditions can be reached only if the age of the overriding plate is <5 Ma. These modeling results indicate that spontaneous subduction initiation would be more favorable for creating high-temperature conditions. However, the synthesis of our modeling results and geological observations indicate that the majority of the metamorphic soles likely formed during induced subduction initiation that involved a young overriding plate. The mechanism of subduction initiation is a key to modern plate tectonics. Here, using numerical modeling and geological observations, the authors find that the majority of active and paleo subduction zones with metamorphic soles likely formed during induced subduction initiation that involved a young overriding plate.

The long‐term state of stress in the subduction forearc depends on the balance between margin‐normal compression due to the plate‐coupling force and the margin‐normal tension due to the gravitational force on the margin topography. In most subduction margins, the outer forearc is largely in margin‐normal compression due to the dominance of the plate‐coupling force. The inner forearc's state of stress varies within and among subduction zones, but what gives rise to this variation is unclear. We examine the state of stress in the forearc region of nine subduction zones by inverting focal mechanism solutions for shallow forearc crustal earthquakes for five zones and inferring the previous inversion results for the other four. The results indicate that the inner forearc stress state is characterized by margin‐normal horizontal deviatoric tension in parts of Nankai, Hikurangi, and southern Mexico. The vertical and margin‐normal horizontal stresses are similar in magnitudes in northern Cascadia as previously reported and are in a neutral stress state. The inner forearc stress state in the rest of the study regions is characterized by margin‐normal horizontal deviatoric compression. Tension in the inner forearc tends to occur where plate coupling is shallow. A larger width of the forearc also promotes inner‐forearc tension. However, regional tectonics may overshadow or accentuate the background stress state in the inner forearc, such as in Hikurangi. Plain Language Summary The state of stress in the overriding plate between the trench and the volcanic arc of subduction zones depends on frictional coupling between the overriding and subducting plates and gravitational force, which causes lateral compression and tension, respectively. The trench‐ward portion of this so‐called “forearc” region is generally in compression due to the dominant effect of plate coupling, but for the arc‐ward portion, the relative importance of the two forces varies spatially. We constrain the state of stress in the forearc using earthquake data for five subduction zones and inferring results from previous studies for four other subduction zones. The results indicate that the forearc stress state seems to correlate with the downdip depth of plate coupling and the width of the forearc. A relatively shallow downdip extent of coupling in a wide forearc tends to have the arc‐ward portion in tension or neutral stress state as observed in parts of Nankai, Hikurangi, Mexico, and Cascadia although this tendency is impacted by the local tectonic settings. Key Points Focal mechanism inversion results indicate the correlation of inner forearc stress state with the downdip depth of plate coupling Margin‐normal horizontal deviatoric tension in the inner forearc tends to occur where plate coupling is shallow and the forearc is wide The variation in the inner forearc stress state does not require a variation in the subduction fault strength

Clustering of arc volcanoes in subduction zones indicates along-arc variation in the physical condition of the underlying mantle where majority of arc magmas are generated. The sub-arc mantle is brought in from the back-arc largely by slab-driven mantle wedge flow. Dynamic processes in the back-arc, such as small-scale mantle convection, are likely to cause lateral variations in the back-arc mantle temperature. Here we use a simple three-dimensional numerical model to quantify the effects of back-arc temperature perturbations on the mantle wedge flow pattern and sub-arc mantle temperature. Our model calculations show that relatively small temperature perturbations in the back-arc result in vigorous inflow of hotter mantle and subdued inflow of colder mantle beneath the arc due to the temperature dependence of the mantle viscosity. This causes a three-dimensional mantle flow pattern that amplifies the along-arc variations in the sub-arc mantle temperature, providing a simple mechanism for volcano clustering. Volcanoes may cluster along volcanic arcs, but controls on cluster locations are still unclear. Here, using numerical models, the authors show that clustering along arcs is the result of inflow of hotter mantle material, and subdued inflow of colder material creating a 3D mantle flow pattern.

Late glacial sea level curves located in the Cascadia subduction zone (CSZ) fore arc in southwestern British Columbia show that glacial isostatic adjustment (GIA) was rapid when the Cordilleran Ice Sheet collapsed in the late Pleistocene. GIA modeling with a linear Maxwell rheology indicates that the observations can be equally well fit across a wide range of asthenospheric thicknesses, provided that the asthenospheric viscosity is varied from 3 × 1018 Pa s for a thin (140 km) asthenosphere to 4 × 1019 Pa s for a thick (380 km) asthenosphere. Present‐day vertical crustal motion predicted by the GIA models shows rates of a few tenths of a millimeter per year, consistent with previous analyses. The model viscosities largely pertain to the viscosity of the oceanic mantle beneath the subducting Juan de Fuca slab but include a contribution from the mantle wedge above the slab. For comparison, effective viscosities for the upper mantle due to the tectonic regime (subduction) were computed using the strain rates and temperatures of an independent geodynamic model of the CSZ with a wet olivine power law rheology. The effective viscosities agree well with GIA model viscosities of 1019 Pa s or less, corresponding to an asthenosphere of 100 or 200 km thickness. The agreement suggests a significant role for power law flow in the GIA response. Regardless of the microphysical mechanisms responsible for the GIA response, the viscosity values inferred from GIA can be applied to studies of the megathrust earthquake cycle because both processes take place on comparable time scales.

Mineral grain size plays an important role in controlling many processes in the mantle wedge of subduction zones, including mantle flow and fluid migration. To investigate the grain‐size distribution in the mantle wedge, we coupled a two‐dimensional (2‐D) steady state finite element thermal and mantle‐flow model with a laboratory‐derived grain‐size evolution model. In our coupled model, the mantle wedge has a composite olivine rheology that incorporates grain‐size‐dependent diffusion creep and grain‐size‐independent dislocation creep. Our results show that all subduction settings lead to a characteristic grain‐size distribution, in which grain size increases from 10 to 100 μm at the most trenchward part of the creeping region to a few centimeters in the subarc mantle. Despite the large variation in grain size, its effect on the mantle rheology and flow is very small, as >90% of the deformation in the flowing part of the creeping region is accommodated by grain‐size‐independent dislocation creep. The predicted grain‐size distribution leads to a downdip increase in permeability by ∼5 orders of magnitude. This increase is likely to promote greater upward migration of aqueous fluids and melts where the slab reaches ∼100 km depth compared with shallower depths, potentially providing an explanation for the relatively uniform subarc slab depth. Seismic attenuation derived from the predicted grain‐size distribution and thermal field is consistent with the observed seismic structure in the mantle wedge at many subduction zones, without requiring a significant contribution by the presence of melt. Key Points All subduction settings lead to a characteristic grain‐size distribution Grain size increases from the fore‐arc wedge to the subarc mantle Grain size affects mantle strength, fluid migration, and seismic attenuation

In the forearc mantle wedge, the thermal field depends strongly on slab‐driven mantle wedge flow. The flow is in turn affected by the thermal field via the temperature dependence of mantle rheology. Using thermal modeling, we show that the nonlinear feedback between the thermal and flow fields always leads to complete stagnation of the mantle wedge over a shallow, weakened part of the slab‐mantle interface and an abrupt onset of mantle flow further down‐dip. The abrupt increase in flow velocity leads to a sharp thermal transition from a cold stagnant to a hot flowing part of the wedge. This sharp thermal transition is inherent to all subduction zones, explaining a commonly observed sharp arc‐ward increase in seismic attenuation. Key Points The temperature dependence of mantle rheology causes nonlinear wedge flow The nonlinear wedge flow always results in a sharp thermal transition The thermal transition can explain the observed change in seismic attenuation

The volume increase that accompanies many hydration reactions can stress and fracture the surrounding rock, a process commonly called reaction-induced fracturing. Reaction-induced fracturing accelerates the rate of hydration by creating new pathways for fluids to migrate into reactive rock and by generating new reactive surface areas. The evolution of reaction-induced fractures also depends on the background stress state, which varies among different tectonic environments. We investigate the impact of tectonic stresses on reaction-induced fracturing using 2-D hydraulic-chemical-mechanical distinct element models. The results indicate that the general pattern of reaction-induced fractures depends on the orientation of background tectonic stresses relative to fluid-supplying channels. A spalling fracture pattern characterized by short cracks parallel to and along fluid-supplying channels occurs when the maximum principal tectonic stress is parallel to the channels whereas a branching fracture pattern characterized by long tensile cracks that propagate in a hierarchical manner into unreacted parts of the rock is expected when the tectonic stress is hydrostatic or when the maximum principal tectonic stress is normal to fluid-supplying channels. Spalling localizes hydration and fluid flow along the channels whereas branching promotes spatially extensive hydration and fluid flow away from the fluid supply. The results indicate tectonic stresses may guide the hydration distribution in the oceanic lithosphere at mid-ocean ridges and outer rises and in the cold mantle wedge corner in subduction zones.

Fluid and melt transport in the solid mantle can be modeled as a two-phase flow in which the liquid flow is resisted by the compaction of the viscously deforming solid mantle. Given the wide impact of liquid transport on the geodynamical and geochemical evolution of the Earth, the so-called “compaction equations” are increasingly being incorporated into geodynamical modeling studies. When implementing these equations, it is common to use a regularization technique to handle the porosity singularity in the dry mantle. Moreover, it is also common to enforce a positive porosity (liquid fraction) to avoid unphysical negative values of porosity. However, the effects of this “capped” porosity on the liquid flow and mass conservation have not been quantitatively evaluated. Here, we investigate these effects using a series of 1- and 2-dimensional numerical models implemented using the commercial finite-element package COMSOL Multiphysics®. The results of benchmarking experiments against a semi-analytical solution for 1- and 2-D solitary waves illustrate the successful implementation of the compaction equations. We show that the solutions are accurate when the element size is smaller than half of the compaction length. Furthermore, in time-evolving experiments where the solid is stationary (immobile), we show that the mass balance errors are similarly low for both the capped and uncapped (i.e., allowing negative porosity) experiments. When Couette flow, convective flow, or subduction corner flow of the solid mantle is assumed, the capped porosity leads to overestimations of the mass of liquid in the model domain and the mass flux of liquid across the model boundaries, resulting in intrinsic errors in mass conservation even if a high mesh resolution is used. Despite the errors in mass balance, however, the distributions of the positive porosity and peaks (largest positive liquid fractions) in both the uncapped and capped experiments are similar. Hence, the capping of porosity in the compaction equations can be reasonably used to assess the main pathways and first-order distribution of fluids and melts in the mantle.