Overriding Plate Research Articles

Flat-Slab Mechanics Transforming Double-Arc Expressions Along the Middle America Trench

AbstractThe southeast–northwest variation in the trench–volcano distance along the Middle America subduction zone is widely recognised to be due to the horizontal (flat) geometry of the descending Cocos plate (slab) in its northwestern part. In the topographic and gravitational expressions, the southeastern segment forms doubled low–high (i.e. trench–continental shelf edge–oceanic basin–volcanic arc) belts, whereas the northwestern segment forms wide doubled low–high belts. To define their attribution, the surface uplift rates were computed from a two-dimensional dislocation-based subduction model with hypocentre-based plate interface geometries. The trench-perpendicular plate-interface profiles exhibited similar curvature compositions, which were reflected in three (shallow, moderate depth and deep) convexities (downward bends) and one concavity (upward bend) between the moderate-depth and deep convex sections. The steady subduction along the first and second positive curvatures (shallow and moderate-depth convexities) in the southeastern and northwestern segments could serve as mechanical sources of short-wavelength double-arc formation. The latter two (negative and third positive) curvatures (concavity and deep convexity) in the northwestern segment produced distinct features including the narrow continental shelf, large trench–volcanic-arc distance and relevant long-wavelength double arcs. This was attributed to the direct contact of the flat elastic slab with the overriding plate.

On the role of trans-lithospheric faults in the long-term seismotectonic segmentation of active margins: a case study in the Andes

Abstract. Plate coupling plays a fundamental role in the way in which seismic energy is released during the seismic cycle. This process includes quasi-instantaneous release during megathrust earthquakes and long-term creep. Both mechanisms can coexist in a given subduction margin, defining a seismotectonic segmentation in which seismically active segments are separated by zones where ruptures stop, classified for simplicity as asperities and barrier, respectively. The spatiotemporal stability of this segmentation has been a matter of debate in the seismological community for decades. In this regard, we explore in this paper the potential role of the interaction between geological heterogeneities in the overriding plate and fluids released from the subducting slab towards the subduction channel. As a case study, we take the convergence between the Nazca and South American plates between 18–40° S, given its relatively simple convergence style and the availability of a high-quality instrumental and historical record. We postulate that trans-lithospheric faults striking at a high angle with respect to the trench behave as large fluid sinks that create the appropriate conditions for the development of barriers and promote the growth of highly coupled asperity domains in their periphery. We tested this hypothesis against key short- and long-term observations in the study area (seismological, geodetic, and geological), obtaining consistent results. If the spatial distribution of asperities is controlled by the geology of the overriding plate, seismic risk assessment could be established with better confidence.

H2 Production in Subduction Zone

Convergence zones where ophiolites outcrops, such as Oman and New Caledonia, allow the generation of natural H2 by redox phenomena. Nevertheless, the role of hydration of the mantle wedge between the overriding plate and the slab is often evoked to explain the H2 emanations. In the case of active subduction, such as that of the Nazca plate beneath the Andes, the presence of H2 emanation has been noted. In this study, we modelled the processes of (de)hydration and melting within the subduction zone. The quantification of water involved in mantle wedge hydration served as a proxy for H2 production. Multiple simulations have been conducted to investigate oceanic-continental subduction, wherein various parameters, especially the subduction angle, have been varied. It appears that in the case of flat subduction, the mantle hydration zone is large, extending up to 500 km from the trench. The extend is controlled by the velocity of the overriding plate. In the case of a steep subduction, the zone is narrower, and is located between the trench and the volcanic arc. Magma formation competes with H2 generation for the use of water expelled from the subducting plate. In the transition from a steep to a flat slab, the mantle hydration zone widens and the volcanic zone moves away from the trench. The oceanic crust may undergo melting, leading to a change in magma composition, before volcanism diminishes in intensity and then disappears.

Spreading ridge migration enabled by plume-ridge de-anchoring

It has long been recognised that spreading ridges are kept in place by competing subduction forces that drive plate motions. Asymmetric strain rates pull spreading ridges in the direction of the strongest slab pull force, which partially explains why spreading ridges can migrate vast distances. However, the interaction between mantle plumes and spreading ridges plays a relatively unknown role on the evolution of plate boundaries. Using a numerical model of mantle convection, we show that plumes with high buoyancy flux (>3000 kg/s) can capture spreading ridges within a 1000 km radius and anchor them in place. Exceptionally high buoyancy fluxes may fragment the overriding plate into smaller plates to accommodate more efficient plate motion. If the plume buoyancy flux wanes below 1000 kg/s the ridge may be de-anchored, leading to rapid ridge migration rates when combined with asymmetric plate boundary forces. Our results show that plume-ridge de-anchoring may have contributed to the rapid migration of the SE Indian Ridge from 43 million years ago (Ma) due to waning buoyancy flux from the Kerguelen plume, supported by magma flux estimates and radiogenic isotope geochemistry of eruption products. The plume-ridge de-anchoring mechanism we have identified has global implications for the evolution of plate boundaries near mantle plumes.

Seismotectonic Setting of the Andes along the Nazca Ridge Subduction Transect: New Insights from Thermal and Finite Element Modelling

The structural evolution of Andean-type orogens is strongly influenced by the geometry of the subducting slab. This study focuses on the flat-slab subduction of the Nazca Ridge and its effects on the South American Plate. The process of flat slab subduction impacts the stress distribution within the overriding plate and increases plate coupling and seismic energy release. Using the finite element method (FEM), we analyse interseismic and coseismic deformation along a 1000 km transect parallel to the ridge. We examine stress distribution, uplift patterns, and the impact of megathrust activity on deformation. To better define the crust’s properties for the model, we developed a new thermal model of the Nazca Ridge subduction zone, reconstructing the thermal structure of the overriding plate. The results show concentrated stress at the upper part of the locked plate interface, extending into the Coastal and Western Cordilleras, with deeper stress zones correlating with seismicity. Uplift patterns align with long-term rates of 0.7–1 mm/yr. Cooling from flat-slab subduction strengthens the overriding plate, allowing far-field stress transmission and deformation. These findings provide insights into the tectonic processes driving stress accumulation, seismicity, and uplift along the Peruvian margin.

Two-stage Mesozoic oceanic subduction and related mantle metasomatism beneath the South Qiangtang terrane with implications for post-collisional magmatism

Subduction-related metasomatism at convergent margins is a crucial process for the lithosphere evolution of the overriding plate. In the central Tibetan Plateau, tectonic transitions and related mantle metasomatism history of the oceanic subduction beneath the South Qiangtang terrane (SQT) remain ambiguous, which impeded understanding the geological evolution of the SQT during the Mesozoic and interpreting the geodynamic process responsible for the Cenozoic post-collisional magmatism. This paper integrated the new and published geochronological and geochemical data from the Jurassic-Cretaceous magmatic rocks in the western part of the SQT to explore the variations of magmatism and processes of mantle modifications. The Jurassic arc magmatism is characterized by the southward younging migration with increasing zircon εHf (t) values and the east–west linear distribution of ca. 155–150 Ma slab-derived adakites. After a magmatic gap at ca. 145–130 Ma, the Early Cretaceous magmatism (130–100 Ma) with distinct juvenile isotopic compositions reinitiated firstly in the south of the Jurassic magmatic belt. In comparison, a younger phase of felsic rocks (120–100 Ma) spatially and geochemically overlaps the northern Jurassic granitoids with ancient crustal sources. These spatial–temporal-geochemical covariations of magmatism could be best explained by the southward transference of oceanic subduction at the earliest Cretaceous that was induced by the Jurassic accretion of an oceanic plateau onto the SQT. Under such subduction regimes, several high-Nb mafic rocks in the northern and southern magmatic belts which were formed by partial melting of the slab-melts-metasomatized lithospheric mantle respectively recorded the Jurassic and Cretaceous mantle metasomatism beneath the SQT, based on their spatial and geochemical relationships with two phases of the adakites. Combined with other geological evidence, these results contribute to elucidating the Mesozoic geodynamic evolution of the oceanic subduction beneath the SQT, which further has implications on the deep dynamic processes responsible for post-collisional magmatism.

Rheology modification in a subduction channel due to eclogite facies metasomatism (Rocky Beach Metamorphic Mélange, Port Macquarie, Australia)

The rheological properties of the interface between the down-going and overriding plates in subduction zones provides insight into how plate convergence is accommodated and the controls on seismic and aseismic slip. This interface is known as the subduction channel and exhumed examples provide the only direct information on deformation mechanisms and the impact of metamorphism on rheology. The Rocky Beach Metamorphic Mélange in eastern Australia is one such exhumed subduction channel, composed of eclogite, blueschist and greenschist facies blocks within a mélange matrix. Previous phase equilibria modelling indicates that high pressure blocks were subducted to ca. 100 km depth and then retrogressed during return flow and exhumation. We found that the rheology of blocks is modified by metasomatism, consistent with studies on other subduction channels. However, through comparison of blocks from different metamorphic grades we found that the effect of metasomatism on rheology varied depending on the pressure and temperature conditions of metasomatism. While unmetasomatised eclogites behaved as rigid objects in the mélange matrix, rocks with mineral assemblages that equilibrated during eclogite facies metasomatism accumulated significant strain, forming isoclinal folds and refolded folds. Deformation of these blocks began at eclogite facies and continued during return flow and retrogression to blueschist facies. At blueschist facies, metasomatised blocks developed mm-scale isoclinal folds with shearing parallel to fold limbs forming rootless isoclinal folds. At the transition between blueschist and greenschist facies, pressure solution became important, preferentially focusing along layers of lawsonite, dissolving it from the rock. At greenschist facies, dissolution-precipitation processes caused significant mass loss, producing mm-spacing between pressure solution seams and cuspate folds, analogous to dewatering structures in sediments. In the Rocky Beach Metamorphic Mélange eclogite facies metasomatism reduces the competence of rigid blocks, reducing overall subduction channel heterogeneity during return flow. We suggest that subduction channels that experience widespread eclogite facies metasomatism may be less likely to generate seismic slip during return flow, since the proportion of rigid blocks and block strength are both reduced.

Overriding Plate Deformation Controls Inferences of Interseismic Coupling Along the Himalayan Megathrust

AbstractThe seismic hazard along the Himalayan arc remains a focus of discussion due to its huge potential societal impact. Taking advantage of modern, dense geodetic observations, several attempts have been made to provide a clear picture of the present‐day rate of strain accumulation caused by interseismic coupling on the Main Himalayan Thrust (MHT). However, there are differing opinions regarding the interpretation of spatial variation in interseismic coupling along some parts of the arc. Critically, the resolution of heterogeneity of coupling is limited due to sparse geodetic observations in some areas. In the present work, we use an updated compilation of all available Global Navigation Satellite System (GNSS) data along the Himalayan arc and a suite of kinematic block models to account for deformation within the overriding plate to characterize the status of interseismic plate coupling. Our results show that the MHT is highly coupled (>0.8) along its entire length and the coupling distribution is nearly binary along‐dip. This translates to a high seismic hazard in the densely populated Indo‐Gangetic plains with a seismic moment accumulation of one Mw 8.7 megathrust earthquake per 100 years. Our results suggest that previously inferred low coupling zones along the MHT are possible manifestations of block modeling artifacts, where the fault segments in the overriding plate control the megathrust slip distribution more strongly than structures within the Indian lithosphere. This result highlights a strong need for better characterization of deforming structures in the overriding plate within the Himalayas and southern Tibetan plateau.

Modeling Subduction With Extremely Fast Trench Retreat

AbstractThe Tonga‐Kermadec subduction zone exhibits the fastest observed trench retreat and convergence near its northern end. However, a paradox exists: despite the rapid trench retreat, the Tonga slab maintains a relatively steep dip angle above 400 km depth. The slab turns flat around 400 km, then steepening again until encountering a stagnant segment near 670 km. Despite its significance for understanding slab dynamics, no existing numerical model has successfully demonstrated how such a distinct slab morphology can be generated under the fast convergence. Here we run subduction models that successfully reproduce the slab geometries while incorporating the observed subduction rate. We use a hybrid velocity boundary condition, imposing velocities on the arc and subducting plate while allowing the overriding plate to respond freely. This approach is crucial for achieving a good match between the modeled and observed Tonga slab. The results explain how the detailed slab structure is highly sensitive to physical parameters including the seafloor age and the mantle viscosity. Notably, a nonlinear rheology, where dislocation creep reduces upper mantle viscosity under strong mantle flow, is essential. The weakened upper mantle allows for a faster slab sinking rate, which explains the large dip angle. Our findings highlight the utilizing rheological parameters that lead to extreme viscosity variations within numerical models to achieve an accurate representation of complex subduction systems like the Tonga‐Kermadec zone. Our study opens new avenues for further study of ocean‐ocean subduction systems, advancing our understanding of their role in shaping regional and global tectonics.

Anisotropic structure at shallow depths across the Japan Trench

Anisotropic structures within the crust are frequently perceived to originate from stress-induced cracks, which have been mainly estimated on land through different wave speeds of orthogonally polarized S waves propagating in the anisotropic media. However, such estimations of crustal anisotropic structures in ocean areas, particularly for subduction zones around trenches, have not been investigated in detail due to the lack of long-term ocean bottom observations. In this study, we used ocean bottom seismometers of a permanent network deployed across the Japan Trench and the southern part of the Kuril Trench and applied the shear-wave splitting analysis to P-to-s converted waves extracted by receiver function analyses using teleseismic events. We estimated the anisotropic structures in marine sediments and oceanic crust for the incoming Pacific Plate and marine sediments for the overriding North American Plate. The obtained fast polarization directions for the incoming plate are mainly oriented to be parallel to the trench axis for the marine sediment and oceanic crust, which are formed by normal faults and cracks due to the upward plate bending in the outer-rise region, whereas results for marine sediments at the northern part of the Japan Trench are obliquely aligned to the trench axis. The oblique direction is consistent with the magnetic lineations of the incoming plate, indicating that ancient faults within the plate, which were formed in the shallow part of the crust during the creation of the oceanic plate at the ridge, are reactivated by the plate flexure. For the overriding plate, the fast polarization directions in the northern and southern parts of the study area are nearly normal to the trench axis. The central part shows two distinct features: the fast polarization directions parallel to the trench axis and small degrees of anisotropy. These patterns may reflect crack alignments associated with the lateral variation in postseismic crustal deformation after the 2011 Tohoku-Oki earthquake. Our results suggest substantial lateral variations in the stress field at the tip of the overriding plate along the strike direction.Graphical

Is the Mesoarchean Mulgandinnah shear zone, Pilbara Craton, the world's oldest arc-slicing transform fault?

Abstract Arc-slicing transform faults represent an integral component of convergent margin tectonics. They are developed above oblique subduction systems, cutting through and displacing the entire magmatic section of arcs, leading to tectonic repetition of segments of the overriding plate in the ensuing orogenic collage. Extant examples clearly show this process in Sumatra, New Zealand, and the Philippines, while ancient examples are reported from the Paleozoic Altaids and Neoarchean Superior and Yilgarn cratons. Here, we report data that document that the Paleo-Mesoarchean Eastern Pilbara craton, recently interpreted to be a preserved mid-upper crustal level of a magmatic arc, is cut and repeated by a major 3.0–2.93 Ga arc-slicing fault, the Mulgandinnah, which sliced a previously 600 × 100 km segment of a Mesoarchean arc system, laterally moving different segments to their presently juxtaposed 200 × 200 km preserved fragment. This evidence demonstrates lateral plate motions by 3.0 Ga and shows oblique subduction, arc plutonism, arc-slicing, and repetition, reflecting that crustal growth in modern-style convergent margins was in full operation by the Mesoarchean.

The role of continental heterogeneity on the evolution of continental plate margin topography at subduction zones

The nature of the overriding plate plays a major role for subduction zone processes. In particular, the highly heterogeneous continental lithosphere modulates intra-plate tectonics and the surface evolution of our planet. However, the role of continental heterogeneity is relatively under-explored for the dynamics of subduction models. We investigate the influence of rheological and density variations across the overriding plate on the evolution of continental lithosphere and slab dynamics in the upper mantle. We focus on the effects of variations in continental plate margin and keel properties on deformation, topographic signals, and basin formation. Our results show that the thickness, extent, and strength of the continental plate margin and subcontinental keel play a crucial role for the morphology and topography of the overriding plate, as well as the retreat of the subducting slab. We show that this lateral heterogeneity can directly influence the coupling between the subducting and overriding plate and determine the partitioning of plate velocities across the overriding plate. These findings suggest that back-arc extension and subsidence are not solely controlled by slab dynamics but are also influenced by continental plate margin and keel properties. Large extended back-arc regions, such as the Pannonian and Aegean basins, may result from fast slab rollback combined with a weak continental plate margin and a strong and extended continental keel. Narrow margins, like the Okinawa Trough in NE Japan, may indicate a comparatively stronger continental plate margin and weaker or smaller continental keel. Additionally, continental keel properties may affect the overall topography of the continental lithosphere, leading to uplift of the deformation front and the formation of intermontane basins.

Uplift of the Tibetan Plateau driven by mantle delamination from the overriding plate

The geodynamic evolution of the Tibetan Plateau remains highly debated. Any model of its evolution must explain the plateau’s growth as constrained by palaeo-altitude studies, the spatio-temporal distribution of magmatic activity, and the lithospheric mantle removal inferred from seismic velocity anomalies in the underlying mantle. Several conflicting models have been proposed, but none of these explains the first-order topographic, magmatic and seismic features self-consistently. Here we propose and test numerically an evolutionary model of the plateau that involves gradual peeling of the lithospheric mantle from the overriding plate and consequent mantle and crustal melting and uplift. We show that this model successfully reproduces the successive surface uplift of the plateau to more than 4 km above sea level and is consistent with the observed migration of magmatism and geometry of the lithosphere–asthenosphere boundary resulting from subduction of the Indian plate and delamination of the mantle lithosphere of the Eurasian plate. These comparisons indicate that mantle delamination from the overriding plate is the driving force behind the uplift of the Tibetan Plateau and, potentially, orogenic plateaus more generally.

Dynamic Component of the Asthenosphere: Lateral Viscosity Variations Due To Dislocation Creep at the Base of Oceanic Plates

AbstractThe asthenosphere is commonly defined as an upper mantle zone with low velocities and high attenuation of seismic waves, and high electrical conductivity. These observations are usually explained by the presence of partial melt, or by a sharp contrast in the water content of the upper mantle. Low viscosity asthenosphere is an essential ingredient of functioning plate tectonics. We argue that a substantial component of asthenospheric weakening is dynamic, caused by dislocation creep at the base of tectonic plates. Numerical simulations of subduction show that dynamic weakening scales with the surface velocity both below the subducting and the overriding plate, and that the viscosity decrease reaches up to two orders of magnitude. The resulting scaling law is employed in an apriori estimate of the lateral viscosity variations (LVV) below Earth's oceans. The obtained LVV help in explaining some of the long‐standing as well as recent problems in mantle viscosity inversions.

Intermittent terrane arrival induces pulses of inland tectonic cycles

Pulsing volcanotectonic cycles characterized by short-lived (10 s Myr) switches in magmatic and tectonic patterns have been broadly identified in active margins. However, the specific mechanism that causes these switches remains ambiguous, i.e., whether the subduction continuity and/or terrane arrival (accretion, underthrusting, or subduction of buoyant continental/oceanic blocks with a thicker crust than their surrounding oceanic plate) plays a crucial role in controlling the observed volcanotectonic cycles remains controversial. Here, by modeling subduction processes involving the sequential arrival of buoyant terranes, we show that 1) in scenarios where the oceanic plate is weakly coupled with terranes, the entraining of the terrane into the subduction induces slab breakup to occur between the partially subducted terrane and its adjacent oceanic slab, with the terrane then rebounding and moving away from the trench. This evolution causes switches in the magmatic and tectonic patterns within the overriding plate. 2) In these models, the exposed terrane materials can preserve characteristic pressure-temperature-time trajectories, i.e., nearly isothermal compression to isothermal decompression and/or isobaric heating after partial exhumation. 3) slab breakup does not guarantee the occurrence of a trench jump; only when the terrane has a medium scale (∼300 km) will a new trench tend to develop prominently behind the rebounded terrane. 4) in scenarios where the composite slab (terrane and oceanic portions) resists yielding deformation and the terrane density is close to that of the oceanic plate (<0.6% density contrast), continuous subduction will occur. In this latter scenario, slab deformation (revealed by subduction angle and slab curvature), instead of slab breakup, will control the magmatic and tectonic patterns in the overriding plate. By further comparing model results with observations, we demonstrate that intermittent subduction interrupted by the subduction of terranes can be a tectonic driver for episodes of compression-to-extension transformations and magmatism (or piston-like volcanotectonic cycles) in two representative accretionary belts — with or without trench jumps (exemplars in south-central-Tibet and eastern-Mediterranean). In contrast, continuous subduction with strongly coupled upper and subducting plates could have contributed to similar cycles in the example without accretion (e.g., the Altiplano in the Central Andes). Therefore, over 10 s of Myr, the scale and frequency of terrane arrivals could essentially control the specific motion pattern of the subducting plate, creating the observed short-lived volcanotectonic switches at these three subtypes of active margins.

Low-Temperature Thermochronology Records the Convergence between the Anatolide–Tauride Block and the Arabian Platform along the Southeast Anatolian Orogenic Belt

SE Anatolia is witnessing the final stage of the Wilson Cycle, where a continental collision between the Tauride–Anatolide block and Arabian platform occurred, and a 1.5 km Eastern Tauride mountain chain formed. We present new low-temperature thermochronology (LTT) ages, including eight apatite fission track (AFT) and seven apatite and zircon U-Th-Sm/He (AHe, ZHe) ages, for the metamorphic rocks from the Nappe Zone of the Southeast Anatolian Orogenic Belt. The ZHe ages vary from 51.2 ± 0.7 Ma to 30.4 ± 0.6 Ma, the AFT ages range from 33.1 ± 1.6 Ma to 18.1 ± 0.9 Ma, and the AHe ages range from 23.6 ± 2.5 Ma to 6 ± 1.9 Ma. The LTT data show a continuous slow uplift of the region. However, the thermal modeling results suggest an Eocene and middle–late Miocene fast uplift of the region. Similar to our results, the LTT studies along the SAOB show that the vertical movements initiated during the Eocene period have continued in a steady-state regime to recent times. The Eocene epoch is identified by arc–back-arc setting in the region, whereas the Miocene epoch is marked by the continental collision. Within this tectonic framework, vertical movements on the overriding plate are controlled by both extensional and compressional tectonics. The LTT data obtained along the SAOB show fingerprints of thrust propagation from north to south.

Detrital zircon REE and tectonic settings

Earth's long-term history is preserved in the rock archive, yet this record is incomplete, hindering our understanding of how tectonic processes have evolved and shaped our planet's evolution. The physio-chemical resilience of zircon and its ability to be rapidly analyzed for age, isotopic, and elemental data has made it a key phase in unraveling Earth's long-term evolution. The europium anomaly (Eu/Eu*) and the differentiation degree between light and heavy rare earth elements (LREEN¯/HREEN¯) of zircons positively correlate with the pressure of magma crystallization because high pressure inhibits plagioclase crystallization and promotes garnet growth. However, this correlation can be influenced by fluctuations in redox conditions, protolith composition, and water content of magma, whose stability varies in different tectonic settings. Therefore, the correlation coefficient between Eu/Eu* and LREEN¯/HREEN¯ of detrital zircons (rDz) provides a new proxy to help distinguish tectonic settings. We evaluate this new proxy using available global detrital zircon REE data from different tectonic settings in modern and deep times. The data from well-constrained simple tectonic settings suggest that rDz is higher in convergent settings (0.53–0.85) than in collisional settings (0.12–0.51). Convergent settings can also yield relatively low rDz (0.28–0.55) when there is a contribution of a mantle plume lying below the overriding plate of a subduction zone.When studying a detrital zircon archive derived from a large distributive province that includes sources from a variety of settings, known provenance information should be used to disentangle mixed source signals. Also, only considering detrital zircon grains that formed a short time (e.g., 50 or 100 Myr) prior to deposition can reduce the likelihood of multi-cycle distal mixing of contemporaneous but tectonically distinct zircons in a sample. Additionally, greater confidence in interpretations of changes in tectonic settings may be derived from evolution patterns of rDz instead of isolated values. A decreasing trend can likely be interpreted as a transition from convergent to collisional or mantle plume settings, whereas an increasing trend implies an opposite tectonic transition.

A HIMU-like component in Mariana Convergent Margin magma sources during initial arc rifting revealed by melt inclusions

Compositions of island arc and back-arc basin basalts are often used to trace the recycling of subducted materials. However, the contribution of subducted components to the mantle source during initial arc rifting before back-arc basin spreading is not yet well constrained. The northernmost Mariana arc is ideal for studying this because the transition from rifting to back-arc spreading is happening here. Here we report major and trace element and Pb isotopic compositions of olivine-hosted melt inclusions from lavas erupted during initial rifting at 24°N (NSP-24) and compare them with those in active arc front at 21°N and mature back-arc basin at 18°N. NSP-24 high-K melt inclusions have highly radiogenic Pb compositions and are close to those of the HIMU end-member, suggesting the presence of this component in the magma source. The HIMU-like component may be stored in the over-riding plate and released into arc magma with rifting. HIMU-type seamounts may be subducted elsewhere beneath the Mariana arc, but obvious HIMU-type components appear only in the initial stages of arc rifting due to the low melting degree and being consumed during the process of back-arc spreading.

Origin of Philippine Sea Basins During Subduction Initiation in the Western Pacific

AbstractUnderstanding the age and dynamics of the overriding plates allows an assessment of competing subduction initiation hypotheses. The Izu‐Bonin‐Mariana margin in the Western Pacific is a key example of initiation and hence it is important to constrain the age and origin of the oldest igneous crust of the supra‐subduction Philippine Sea Plate. We present geochronological and geochemical data of igneous rocks from the oldest ocean basins of the Philippine Sea Plate: the West Philippine and Palau Basins. Basalts from these basins have enriched geochemical characteristics similar to the EM‐2‐like mantle component found in OIB‐like basalts associated with the Oki‐Daito mantle plume. Ages of basalts from the northernmost West Philippine Basin (WPB) and the Palau Basin range from 43.5 to 50.5 Ma, which is similar to the oldest samples associated with the Oki‐Daito mantle plume (48–50 Ma). This implies that the plume contributed to magmatism from the onset of basin formation. It also provides support for the proposition that rifting of the Mesozoic arc terrane and subsequent seafloor spreading of the WPB was triggered by the arrival of the Oki‐Daito mantle plume at the base of the lithosphere. The age of these Philippine Sea Basins implies that only the Mesozoic Daito Ridge Group and the Gagua Ridge existed as Philippine Sea Plate crust before subduction initiation. A major fault activity after 37 Ma in the northernmost WPB demonstrates that careful reconstruction of the Eocene Philippine Sea Plate is critical to understanding plate dynamics during subduction initiation in the Western Pacific.

Megathrust locking encoded in subduction landscapes.

Locked areas of subduction megathrusts are increasingly found to coincide with landscape features sculpted over hundreds of thousand years, yet the mechanisms that underlie such correlations remain elusive. We show that interseismic locking gradients induce increments of irreversible strain across the overriding plate manifested predominantly as distributed seismicity. Summing these increments over hundreds of earthquake cycles produces a spatially variable field of uplift representing the unbalance of co-, post-, and interseismic strain. This long-term uplift explains first-order geomorphological features of subduction zones such as the position of the continental erosive shelf break, the distribution of marine terraces and peninsulas, and the profile of forearc rivers. Inelastic yielding of the forearc thus encodes short-term locking patterns in subduction landscapes, hinting that megathrust locking is stable over multiple earthquake cycles and highlighting the role geomorphology can play in constraining Earth's greatest source of seismic hazard.