Laboratory models of the thermal evolution of the mantle during rollback subduction

Interplay between oceanic subduction and continental collision in building continental crust

Generation of continental crust in collision zones reflect the interplay between oceanic subduction and continental collision. The Gangdese continental crust in southern Tibet developed during subduction of the Neo-Tethyan oceanic slab in the Mesozoic prior to reworking during the India-Asia collision in the Cenozoic. Here we show that continental arc magmatism started with fractional crystallization to form cumulates and associated medium-K calc-alkaline suites. This was followed by a period commencing at ~70 Ma dominated by remelting of pre-existing lower crust, producing more potassic compositions. The increased importance of remelting coincides with an acceleration in the convergence rate between India and Asia leading to higher basaltic flow into the Asian lithosphere, followed by convergence deceleration due to slab breakoff, enabling high heat flow and melting of the base of the arc. This two-stage process of accumulation and remelting leads to the chemical maturation of juvenile continental crust in collision zones, strengthening crustal stratification.

3D anisotropic structure of the Japan subduction zone

How mantle materials flow and how intraslab fabrics align in subduction zones are two essential issues for clarifying material recycling between Earth's interior and surface. Investigating seismic anisotropy is one of a few viable technologies that can directly answer these questions. However, the detailed anisotropic structure of subduction zones is still unclear. Under a general hexagonal symmetry anisotropy assumption, we develop a tomographic method to determine a high-resolution three-dimensional (3D) P wave anisotropic model of the Japan subduction zone by inverting 1,184,018 travel time data of local and teleseismic events. As a result, the 3D anisotropic structure in and around the dipping Pacific slab is firstly revealed. Our results show that slab deformation plays an important role in both mantle flow and intraslab fabric, and the widely observed trench-parallel anisotropy in the forearc is related to the intraslab deformation during the outer-rise yielding of the subducting plate.

Seismic evidence for subduction-induced mantle flows underneath Middle America

Laboratory experiments and geodynamic simulations demonstrate that poloidal- and toroidal-mode mantle flows develop around subduction zones. Here, we use a new 3-D azimuthal anisotropy model constructed by full waveform inversion, to infer deep subduction-induced mantle flows underneath Middle America. At depths shallower than 150 km, poloidal-mode flow is perpendicular to the trajectory of the Middle American Trench. From 300 to 450 km depth, return flows surround the edges of the Rivera and Atlantic slabs, while escape flows are inferred through slab windows beneath Panama and central Mexico. Furthermore, at 700 km depth, the study region is dominated by the Farallon anomaly, with fast axes perpendicular to its strike, suggesting the development of lattice-preferred orientations by substantial stress. These observations provide depth-dependent seismic anisotropy for future mantle flow simulations, and call for further investigations about the deformation mechanisms and elasticity of minerals in the transition zone and uppermost lower mantle.

The motions of subducted slabs are expected to drive mantle flow around slab edges, however, evidence of deep mantle flow has so far remained elusive. Here, the authors present a Full Waveform Inversion 3-D anisotropy model which allows them to infer deep subduction-induced mantle flows underneath the Mid-Americas and the Caribbean.

A subduction and mantle plume origin for Samoan volcanism

The origin of Samoan volcanism in the southwest Pacific remains enigmatic. Whether mantle melting is solely caused by a mantle plume is questionable because some volcanism, here referred to as non-hotspot volcanism, defies the plume model and its linear age-progression trend. Indeed, non-hotspot volcanism occurred as far as 740 km west of the predicted Samoan hotspot after 5 Ma. Here we use fully-dynamic laboratory subduction models and a tectonic reconstruction to show that the nearby Tonga-Kermadec-Hikurangi (TKH) subduction zone induces a broad mantle upwelling around the northern slab edge that coincides with the non-hotspot volcanic activity after 5 Ma. Using published potential mantle temperatures for the ambient mantle and Samoan mantle plume, we find that two geodynamic processes can explain mantle melting responsible for intraplate volcanism in the Samoan region. We propose that before 5 Ma, the volcanism is consistent with the plume model, whereas afterwards non-hotspot volcanism resulted from interaction between the Subduction-Induced Mantle Upwelling (SIMU) and Samoan mantle plume material that propagated west from the hotspot due to the toroidal component of slab rollback-induced mantle flow. In this geodynamic scenario, the SIMU drives decompression melting in the westward-swept plume material, thus producing the non-hotpot volcanism.

Seismogenesis of dual subduction beneath Kanto, central Japan controlled by fluid release

Dual subduction represents an unusual case of subduction where one oceanic plate subducts on top of another, creating a highly complex tectonic setting. Because of the complex interaction between the two subducted plates, the origin of seismicity in such region is still not fully understood. Here we investigate the thermal structure of dual subduction beneath Kanto, central Japan formed as a consequence of a unique case of triple trench junction. Using high-resolution three-dimensional thermo-mechanical models tailored for the specific dual subduction settings beneath Kanto, we show that, compared with single-plate subduction systems, subduction of double slabs produces a strong variation of mantle flow, thermal and fluid release pattern that strongly controls the regional seismicity distribution. Here the deepening of seismicity in the Pacific slab located under the Philippine Sea slab is explained by delaying at greater depths (~150 km depth) of the eclogitization front in this region. On the other hand, the shallower seismicity observed in the Philippine Sea slab is related to a young and warm plate subduction and probably to the presence of a hot mantle flow traveling underneath the slab and then moving upward on top of the slab.

Trench-Parallel Anisotropy Produced by Foundering of Arc Lower Crust

Many volcanic arcs display fast seismic shear-wave velocities parallel to the strike of the trench. This pattern of anisotropy is inconsistent with simple models of corner flow in the mantle wedge. Although several models, including slab rollback, oblique subduction, and deformation of water-rich olivine, have been proposed to explain trench-parallel anisotropy, none of these mechanisms are consistent with all observations. Instead, small-scale convection driven by the foundering of dense arc lower crust provides an explanation for the trench-parallel anisotropy, even in settings with orthogonal convergence and no slab rollback.

Evolution and diversity of subduction zones controlled by slab width

Subducting slabs provide the main driving force for plate motion and flow in the Earth's mantle, and geodynamic, seismic and geochemical studies offer insight into slab dynamics and subduction-induced flow. Most previous geodynamic studies treat subduction zones as either infinite in trench-parallel extent (that is, two-dimensional) or finite in width but fixed in space. Subduction zones and their associated slabs are, however, limited in lateral extent (250-7,400 km) and their three-dimensional geometry evolves over time. Here we show that slab width controls two first-order features of plate tectonics-the curvature of subduction zones and their tendency to retreat backwards with time. Using three-dimensional numerical simulations of free subduction, we show that trench migration rate is inversely related to slab width and depends on proximity to a lateral slab edge. These results are consistent with retreat velocities observed globally, with maximum velocities (6-16 cm yr(-1)) only observed close to slab edges (<1,200 km), whereas far from edges (>2,000 km) retreat velocities are always slow (<2.0 cm yr(-1)). Models with narrow slabs (< or =1,500 km) retreat fast and develop a curved geometry, concave towards the mantle wedge side. Models with slabs intermediate in width ( approximately 2,000-3,000 km) are sublinear and retreat more slowly. Models with wide slabs (> or =4,000 km) are nearly stationary in the centre and develop a convex geometry, whereas trench retreat increases towards concave-shaped edges. Additionally, we identify periods (5-10 Myr) of slow trench advance at the centre of wide slabs. Such wide-slab behaviour may explain mountain building in the central Andes, as being a consequence of its tectonic setting, far from slab edges.