Overriding Plate Research Articles

Can Plate Bending Explain the Observed Faster Landward Motion of Lateral Regions of the Subduction Zone After Major Megathrust Earthquakes?

AbstractGreater landward velocities were recorded after six megathrust earthquakes in subduction zone regions adjacent to the ruptured portion. Previous explanations invoked either increased slip deficit accumulation or plate bending during postseismic relaxation, with different implications for seismic hazard. We investigate whether bending can be expected to reproduce this observed enhanced landward motion (ELM). We use 3D quasi‐dynamic finite element models with periodic earthquakes. We find that afterslip downdip of the brittle megathrust exclusively produces enhanced trenchward surface motion in the overriding plate. Viscous relaxation produces ELM when a depth limit is imposed on afterslip. This landward motion results primarily from in‐plane elastic bending of the overriding plate due to trenchward viscous flow in the mantle wedge near the rupture. Modeled ELM is, however, incompatible with the observations, which are an order of magnitude greater and last longer after the earthquake. This conclusion does not significantly change when varying mantle viscosity, plate elasticity, maximum afterslip depth, earthquake size, megathrust locking outside of the rupture, or nature and location of relevant model boundaries. The observed ELM consequently appears to reflect faster slip deficit accumulation, implying a greater seismic hazard in lateral segments of the subduction zone.

Impact of the Iquique Ridge on structure and deformation of the north Chilean subduction zone

The subduction of seamounts and basement ridges affects the structure, morphology, and physical state of a convergent margin. To evaluate their impact on the seismo-tectonic setting of the subduction zone and the tectonic development of the lower subducting and upper overriding plate, it is essential to know the precise location of subducted topographic features under the marine forearc. Offshore Northern Chile, the Iquique Ridge represents a broad zone of complex and heterogeneous structure of variable width on the oceanic Nazca Plate, which complicates attempts to project it beneath the forearc of the Chilean subduction zone. Here we use a state-of-the-art seismic reflection data processing approach to map structures related to ridge subduction under the marine forearc with unprecedented accuracy and resolution and evaluate their impact on the deformation of both the plate boundary and the upper plate. We show that significant ridge-related topography is currently subducting south of 20.5 °S and that the combined effect of horst and graben subduction with subduction of Iquique ridge-related thickened and elevated crust causes an upward bulging of the entire upper plate from the plate interface up to the seafloor as well as the presence of kilometer-scale anticlinal structures observed in multibeam bathymetric data that are approximately aligned with horsts seaward of the trench. In the area affected by the subducting ridge, a frontal prism is absent, which may relate to frontal subduction erosion caused by the excess lower plate topography. In contrast farther towards the north, where only isolated seamounts subduct, a small frontal prism and a slope/apron sediment cover down to 3000 m water depth are found.

Compression at Strike‐Slip Fault Is a Favorable Condition for Subduction Initiation

AbstractThe recent statistics suggests that over 60% of active Cenozoic subduction initiation (SI) cases are related to the strike‐slip fault. A number of previous studies have shown that the lithospheric weak zone is a necessary condition for the SI. However, the direct effect of strike‐slip motion on lithospheric weakening and SI has rarely been investigated in numerical models due to the challenge of complex 3D boundary conditions. In this study, a new 3D model has been built with both strike‐slip and compression boundary conditions. The model results indicate that the compression at a strike‐slip boundary provides a favorable condition for the SI, with producing and maintaining a lithospheric‐scale weak zone that facilitates strain localization and SI. In addition, the high strike‐slip velocity and buoyant overriding plate contributes to the SI of young oceanic plate. This mechanism satisfies a large number of natural SI cases in the Cenozoic.

The Gibraltar slab dynamics and its influence on past and present-day Alboran domain deformation: Insights from thermo-mechanical numerical modelling

The origin and tectonic evolution of the Gibraltar Arc system is the result of a complex geodynamic evolution involving the convergence of the Eurasian and African plates and the dynamic impact of the Gibraltar slab. Although geologic and geophysical data collected in the last few years have increased our knowledge of the Gibraltar Arc region, it is still unclear which are the mechanical links between the Gibraltar slab and the past deformation of the overriding Alboran lithosphere, as well as to which degree this subduction system is presently active. In this study, we use 2D numerical modelling to investigate the impact of the Gibraltar slab dynamics on the deformation of the overriding Alboran lithosphere. Our model simulates a WE generic vertical section at an approximate latitude of 36°N and considers an initial setup at about Burdigalian times (∼20 Ma), when the subduction front position is relatively well constrained by recent tectonic reconstructions. Our modelling shows a switch in the overriding plate (OP) stress state from extensional stresses during the slab rollback to compressional stresses near the trench when the rollback velocity decreases, caused by the change in slab-induced mantle flow. We also find that much of the crustal and lithospheric deformation occur during fast slab rollback and OP extension in the first 10 Myr of evolution, while after that only moderate deformation associated with subduction is predicted. Finally, we find that despite the subduction rollback ceases, the ongoing motion of the deeper portion of the slab induces a mantle flow that causes some amount of west-directed basal drag of the Alboran lithosphere. This basal drag generates interplate compresional stresses compatible with the distribution of intermediate-depth earthquakes in western Alboran.

Seismic and aseismic slip during the 2006 Copiapó swarm in North-Central Chile

Earthquake swarms commonly occur along the Chilean subduction zone, witnessing fast seismic and slow aseismic slip behavior at the plate interface. However, the largest seismic swarms observed in Chile, particularly in the Copiapó-Atacama region, remain poorly documented, and the underlying processes have yet to be understood. Here, we perform seismological and geodetic analyses to investigate the 2006 Copiapó swarm, which developed in April and May 2006. The swarm began on April 19, with a magnitude Ml 5.3 earthquake. During the nine following days, we observe a migration of seismicity along the plate interface, the occurrence of doublets events, and a potential slow slip event in the GPS time series at site Copiapó. Then, on April 30, a first earthquake with Mw 6.6 occurred at 15 km depth at the plate contact. It likely triggered a second earthquake of magnitude Mw 6.5, which occurred 144 min later, 10 km northwest of the first earthquake. Using InSAR, we determined the slip distribution associated with these two earthquakes and detailed the postseismic slip they triggered in the next days and weeks. This “postseismic” phase appears to be predominantly aseismic, while the moment released during the “coseismic” phase is comparable to other seismic crises that occurred in Atacama. Although we did not find a larger seismic and aseismic ratio than in other swarms in South America, we suggest a similar mechanism of slow deformation as a driver of seismicity during seismic swarms. Finally, we propose that the slow and fast behavior of the 2006 Copiapó swarm is a consequence of the subduction of the Copiapó Ridge seamounts, which affects both the plate interface and the overriding plate by inducing complex interactions between seismic and aseismic processes.

Geochemistry and zircon U–Pb ages of early Ordovician syenites from the Inexpressible Island, Antarctica and tectonic implications

The Ross Orogenic Belt is in the Antarctica Transantarctic Mountains. North Victoria Land Granite Harbour Intrusive complex (GHI) records the tectonic-magmatism evolution of Ross orogeny. Extensively developed post-collisional granites around this margin of early Paleozoic magmatism can provide insights into the growth of continental crust through accretionary orogenesis. We provide geochemical and geochronological data from syenites from Terra Nova Bay, north Victoria Land in order to constrain its tectonic evolution and setting. The syenite belongs to the potassium-alkaline, calc-alkaline series and is characterized by high concentrations of rare Earth elements and large ion lithophile elements (LILE), and low content in high field strength elements (Nb, Ta, P, Ti). The petrographic and geochemical signatures show a possible island-arc granite affinity. LA-ICP-MS zircon U-Pb dating results suggest that the Inexpressible Island syenite was emplaced at ca. 471.8 ± 1.8 Ma and 477.3 ± 1.7 Ma, respectively. Zircon εHf(t) values range from −7.4 to −9.1; average −8.2 and whole-rock εNd (t) values range from −8.5 to −10.3, indicating that formed by the partial melting of the lithospheric mantle enriched with subduction slab fluids and subcontinental lithosphere. Whereas, the syenite has a strong positive Eu anomaly and a positive Sr anomaly, suggesting that plagioclase cumulate crystallization occurred in the magma source area. Furthermore, through integration with previous studies, we suggest that syenite is a result of the melting zone of an older previously subduction enriched layer of the subcontinental lithospheric mantle (SCLM). To enable syenite emplacement we suggest a tectonic-magmatic model that invokes alternating phases of extension and contraction in the overriding plate. Finally, we report the youngest age of (post-orogenic) magmatism occurred during extension in the overriding plate ca. 478–471 Ma.

Subduction of the Neo-Tethys ridge beneath the Eurasian continent during the Cretaceous

Ridge subduction usually has significant influences on the tectonics, magmatic activity and mineralization of the overriding plate because of its special thermal and topographic structures. The Neo-Tethys is a vanished ancient ocean, and thus, its spreading ridge must have been subducted beneath the Eurasian continent. However, where and when the subduction of the Neo-Tethys ridge occurred has not been well constrained. In this study, a Late Cretaceous (∼100 Ma) EW-trending igneous rock belt along the southern margin of the Eurasian continent from the Gangdese terrane to the South China block is recognized based on residual geological records and plate reconstruction. Importantly, these igneous rocks show the characteristics of typical adakites and are consistent with those of adakites derived from partial melting of subducted oceanic crust. Combined with the tectonic evolution of the Neo-Tethys Ocean, ridge subduction was most likely responsible for this adakitic igneous rock belt. In addition, a NS-trending compression event simultaneously occurred in the southern margin of the Gangdese terrane and the South China block during the Late Cretaceous, which agrees well with ridge subduction. The identification of the Neo-Tethys ridge subduction is of great significance for understanding the evolution, magmatism and mineralization of the Neo-Tethys tectonic domain.

Offshore Landward Motion Shortly After a Subduction Earthquake Implies Rapid Relocking of the Shallow Megathrust

AbstractGeodetic observations after large subduction earthquakes reflect multiple postseismic processes, including megathrust relocking. The timing of relocking and the observational constraints on it are unclear. Relocking was inferred to explain some observed landward motion that occurs within months. It was also considered unable to explain other, greater landward motion, including that off the coast of Japan beginning weeks after the 2011 Tohoku earthquake, attributed to postseismic relaxation. We use generic, 3D numerical models to show that relocking, particularly of the shallow interface, is needed for postseismic relaxation to produce landward motion on the tip of the overriding plate. We argue that this finding is consistent with previous simulations that implicitly relock the megathrust where afterslip is not included. We conclude that the Tohoku megathrust relocked within less than 2 months of the earthquake. This suggests that the shallow megathrust probably behaves as a true, unstably sliding asperity.

Upper-plate Shortening and Mountain-building in the Context of Mantle-driven Oceanic Subduction

The mechanisms controlling mountain building at subduction zones remain debated. In particular the interaction between mantle flow and subduction has been poorly addressed while fundamental in controlling plate displacement and deformation. Here, we conduct three-dimensional analogue models at the scale of the upper mantle adding a horizontal mantle flow, so that plate displacement is not imposed as in most models, but is rather controlled by the balance of forces. We simulate three scenarios: no mantle flow (slab-pull driven subduction), mantle flow directed toward the subducting plate, and mantle flow directed toward the overriding plate. In that last scenario, we test the influence of pre-existing rheological contrasts in the upper plate to best reproduce natural cases where structural and rheological inheritance is common. Our experiments show that when plate convergence is also driven by a background mantle flow, the continental plate deforms with significant trench-orthogonal shortening (up to 30% after 60 Myr), generally associated with thickening. The upper plate shortening and thickening is best promoted when the mantle flow is directed toward the fixed overriding continental plate. The strength of the upper plate is also a key factor controlling the amount and rates of accommodated shortening. Deformation rates increase linearly with decreasing bulk strength of the upper plate, and deformation is mostly localized where viscosity is lower. Finally, we discuss the limits and strengths of our model results through a comparison to the Andes which are the best modern example of mountain building in a subduction context.

Aseismic ridge subduction and flat subduction: Insights from three-dimensional numerical models

Flat subduction can significantly influence the distribution of volcanism, stress state, and surface topography of the overriding plate. However, the mechanisms for inducing flat subduction remain controversial. Previous two-dimensional (2-D) numerical models and laboratory analogue models suggested that a buoyant impactor (aseismic ridge, oceanic plateau, or the like) may induce flat subduction. However, three-dimensional (3-D) systematic studies on the relationship between flat subduction and buoyant blocks are still lacking. Here, we use a 3-D numerical model to investigate the influence of the aseismic ridge, especially its width (which is difficult to consider in 2-D numerical models), on the formation of flat subduction. Our model results suggest that the aseismic ridge needs to be wide and thick enough to induce flat subduction, a condition that is difficult to satisfy on the Earth. We also find that the subduction of an aseismic ridge parallel to the trench or a double aseismic ridge normal to the trench has a similar effect on super-wide aseismic ridge subduction in terms of causing flat subduction, which can explain the flat subduction observed beneath regions such as Chile and Peru.

Short‐Term and Long‐Term Variations in Seismic Velocity at Shallow Depths of the Overriding Plate West of the Japan Trench

AbstractLong‐term and short‐term temporal variations in seismic velocity have been investigated using seismic interferometry, and significant short‐term changes are typically recognized as velocity reductions caused by large earthquakes. In this study, using a dense network of ocean bottom seismometers (Seafloor Observation Network for Earthquakes and Tsunamis along the Japan Trench: S‐net), we measured temporal variations in seismic velocity at the top of the overriding plate in the region between the Japanese Islands and the Japan Trench. We detected abrupt velocity reductions due to the nearby Mw 7‐class earthquakes in 2016 and 2021, and gradual velocity increases mainly observed near the trench. The sudden velocity reductions due to the shallow earthquake in 2016 within the overriding plate are large above the hypocenter, while those due to the three deep earthquakes in 2021 are large near the trench, although these earthquakes occurred far from the trench. The three earthquakes are located at the top of and within the subducting oceanic crust of the Pacific Plate. The spatial distributions of the velocity reductions for the four earthquakes are consistent with those of the peak ground velocities (PGVs). This indicates that for the three deep earthquakes, seismic waves trapped within the descending Pacific Plate reach to the trench, and the trapped waves induce large PGVs, and hence substantially reducing the velocity. Long‐term variations observed near the trench are possibly caused by the recovery process of coseismic velocity reductions due to large earthquakes and/or fluid migration induced by the tectonic stress field in the region.

Geological processes defining the formation of plumasite-type corundum in the Paleoproterozoic Isertoq Terrane, South-East Greenland

Plumasite-type corundum occurrences in the Nattivit area in South-East Greenland offer a unique opportunity to study corundum formation in-situ where pegmatites intruded into metamorphosed lherzolite and dunite of the Archean-Paleoproterozoic continental crust. The Nattivit area, located in the Isertoq Terrane of North Atlantic Craton, forms part of the overriding plate during convergence of the Nagssugtoqidian orogen (1910-1840 Ma). New field observations and elemental and isotopic geochemical analysis provide further insights in the history of crustal convergence, its exhumation and how corundum was formed. The continental crust in the area consists of metamorphosed mafic to ultramafic rocks and tonalite-trondhjemite-granodiorite (TTG) gneisses, where the mafic rocks in the Isertoq Terrane yield a εNd TDM model age of 3000–2800 Ma. Dunite and lherzolite sills/dikes intruded the mafic rocks before the intrusion of the TTG sheets. The intrusion ages for the TTG obtained from zircon U-Pb geochronology are 2818 ± 8 Ma, 2760 ± 13 Ma to 2667 ± 7 Ma. U-Pb zircon data, zircon textures and Th/U ratios indicate metamorphism occurred at 2698 ± 7 Ma to 2629 ± 11 Ma, 2500–2400 Ma and 1900–1600 Ma. Whole rock geochemical data of mafic to ultramafic rocks show a continental arc affinity, with negative Ta, Nb and positive Pb anomalies. A metasomatic event at 2390 ± 70 Ma partly reset the isotopic signature in the mafic to ultramafic rocks. A marked absence of ages between 2350 and 2100 Ma in the TTG zircon age populations exists, indicating a period with minimal magmatic and/or metamorphic activity. The metamorphic mineral assemblages of the schist, amphibolite, ultramafic rocks and metasomatic reaction zones in ultramafic rocks indicate upper to medium–high amphibolite facies conditions. Kyanite in the metasomatic reaction zones in ultramafic rocks indicate the higher end of the temperature and pressure range above 4.2–10 kbar and 530–800 °C, similar to estimates from dolerite dikes in the Kitak area. The syn-tectonic pegmatites with an intrusion age of 1843 ± 4 Ma formed corundum. The new data indicate that the pegmatite melt/fluid and the geotectonic setting are defining factors for generating plumasite-type corundum.

Cascading Hazards in a Migrating Forearc‐Arc System: Earthquake and Eruption Triggering in Nicaragua

AbstractStrain partitioning in oblique convergent margins results in margin‐parallel shear in the overriding plate. Margin‐parallel shear is often accommodated by margin‐parallel strike‐slip faults proximal to active volcanic arcs. Along the Nicaraguan segment of the Central American Forearc (CAFA) in the Cocos‐Caribbean plate convergent margin, there are no well‐expressed right‐lateral faults that accommodate CA‐CAFA relative motion. Instead, historical earthquakes and mapped fault orientations indicate that the ∼12 mm/yr of dextral motion is accommodated on arc‐normal, left‐lateral faults (i.e., bookshelf faults). We investigate three upper‐plate earthquakes; the 10 April 2014 (Mw 6.1), 15 September 2016 (Mw 5.7), and 28 September 2016 (Mw 5.5), using Global Position System co‐seismic displacements and relocated earthquake aftershocks. Our analyses of the three earthquakes indicate that the 10 April 2014 earthquake ruptured an unmapped margin‐parallel right‐lateral fault in Lago Xolotlán (Managua) and the September 2016 earthquakes ruptured arc‐normal, left‐lateral and oblique‐slip faults. These earthquakes represent a triggered sequence whereby the 10 April 2014 earthquake promoted failure of the faults that ruptured in September 2016 by imparting a static Coulomb stress change (ΔCFS) of 0.02–0.07 MPa. Likewise, the 15 September 2016, earthquake additionally promoted failure (ΔCFS of 0.08–0.1 MPa) on sub‐parallel faults that ruptured in two subsequent earthquakes. We also present an instance of magma‐tectonic interaction whereby the 10 April 2014 earthquake dilated (10s of μStrain) the shallow magmatic system of Momotombo Volcano, which led to magma injection, ascent, and eruption on 1 December 2015, after ∼100 years of quiescence.

Diverse Deformation Mechanisms and Lithologic Controls in an Active Orogenic Wedge: Structural Geology and Thermochronometry of the Eastern Greater Caucasus

AbstractOrogenic wedges are common at convergent plate margins and deform internally to maintain a self‐similar geometry during growth. New structural mapping and thermochronometry data illustrate that the eastern Greater Caucasus mountain range of western Asia undergoes deformation via distinct mechanisms that correspond with contrasting lithologies of two sedimentary rock packages within the orogen. The orogen interior comprises a package of Mesozoic thin‐bedded (<10 cm) sandstones and shales. These strata are deformed throughout by short‐wavelength (<1 km) folds that are not fault‐bend or fault‐propagation folds. In contrast, a coeval package of thick‐bedded (up to 5 m) volcaniclastic sandstone and carbonate, known as the Vandam Zone, has been accreted and is deformed via imbrication of coherent thrust sheets forming fault‐related folds of 5–10 km wavelength. Structural reconstructions and thermochronometric data indicate that the Vandam Zone package was accreted between ca. 13 and 3 Ma. Following Vandam Zone accretion, thermal modeling of thermochronometric data indicates rapid exhumation (∼0.3–1 mm/yr) in the wedge interior beginning between ca. 6 and 3 Ma, and a novel thermochronometric paleo‐rotation analysis suggests out‐of‐sequence folding of wedge‐interior strata after ca. 3 Ma. Field relationships suggest that the Vandam Zone underwent syn‐convergent extension following accretion. Together, the data record spatially and temporally variable deformation, dependent on both the mechanical properties of deforming lithologies and perturbations such as accretion of material from the down‐going to the overriding plate. The diverse modes of deformation are consistent with the maintenance of critical taper.

Aseismic ridge subduction focused late Cenozoic exhumation above the Peruvian flat slab

Subduction of aseismic ridges and flat slab subduction are important processes that punctuate Cordilleran orogenesis and may enhance exhumation and rock uplift in the overriding plate. Distinguishing between the two drivers is often challenging, as many modern flat slabs spatially coincide with subducting buoyant ridges. The Peruvian flat slab is the largest region of active flat slab subduction on Earth, spanning over 1300 km of the subducting Nazca plate along the western margin of South America. The flat slab is associated with two seafloor ridges: the Nazca Ridge at the southern terminus and the fully subducted Inca Plateau in the north. These aseismic ridges are spatially confined with respect to the flat slab, allowing assessment of the relative roles of aseismic ridge interactions and flat slab subduction in driving upper plate exhumation. We present: (1) a regional compilation of geochronologic ages of Andean igneous rocks, which track the spatio-temporal evolution of Neogene magmatic arc cessation and hence slab flattening; (2) calculated geomorphic indices, which document landscape perturbations and climatic or lithologic changes, (3) a summary of erosion rates from river catchments on the western Andean slope, and; (4) a regional synthesis of thermochronologic ages that reflect the timing and magnitude of upper crustal cooling.Thermochronometric cooling ages systematically track the progressive passage of the Nazca Ridge, suggesting that the buoyant ridge focused exhumation in the overriding plate. Geomorphic indices demonstrate enhanced topography and steeper channels closer to the position of the subducted ridge. The spatial progression of basement block uplifts in Peru also coincides with the timing of ridge passage. In hinterland regions, >2 km of exhumation occurred since ca. 15 Ma above the Peruvian flat slab. For individual locations within the orogen, active rock uplift declines after ridge passage, suggesting that increased coupling is not maintained across the entirety of the flat slab. We argue that above broad zones of flat slab subduction, focused areas of aseismic ridge subduction concentrate upper-plate exhumation and uplift. These observations may be relevant to other flat slab systems, which exhibit a broad zone of arc shutoff with corridors of focused exhumation.

The Role of Lithospheric‐Deep Mantle Interactions on the Style and Stress Evolution of Arc‐Continent Collision

AbstractWe investigate how the mechanical properties of intra‐oceanic arcs affect the collision style and associated stress‐strain evolution with buoyancy‐driven models of subduction that accurately reproduce the dynamic interaction of the lithosphere and mantle. We performed a series of simulations only varying the effective arc thickness as it controls the buoyancy of intra‐oceanic arcs. Our simulations spontaneously evolve into two contrasting styles of collision that are controlled by a 3% density contrast between the arc and the continental plate. In simulations with less buoyant arcs (15–31 km; effective thickness), we observe arc‐transference to the overriding plate and slab‐anchoring and folding at the 660 km transition zone that result in fluctuations in the slab dip, strain‐stress regime, surface kinematics, and viscous dissipation. After slab‐folding occurs, the gravitational potential energy is dissipated in the form of lithospheric flow causing lithospheric extension in the overriding plate. Conversely, simulations with more buoyant arcs (32–35 km; effective thickness) do not lead to arc‐transference and result in slab break‐off, which causes an asymptotic trend in surface kinematics, viscous dissipation and strain‐stress regime, and lithospheric extension in the overriding plate. The results of our numerical modeling highlight the importance of slab‐anchoring and folding in the 660 km transition zone on increasing the mechanical coupling of the subduction system.

Effects of Elasticity on Subduction Initiation: Insight From 2‐D Thermomechanical Models

AbstractDespite extensive modeling efforts, the dynamics of subduction initiation (SI), including the role of elasticity, are not fully understood. Using two‐dimensional thermomechanical models with visco‐plastic (VP) and visco‐elasto‐plastic (VEP) rheologies, we systematically investigate the role of elasticity in intraoceanic SI using two model setups: spontaneous initiation without imposed convergence and induced initiation with imposed convergence. In spontaneous models, the overriding plate age of <20 Ma and the subducting plate age of >50 Ma generally lead to vertically driven SI with either rheology, but for a given age contrast, SI is easier to occur with the VEP rheology. In induced model with either rheology, when the two plates are young and have a small age contrast, the resulting SI is horizontally driven, and elasticity does not affect SI significantly, regardless of the convergence rate. However, when the thermal age contrast is large and a convergence rate is relatively low, the SI in induced models is vertically driven and similar to that in the spontaneous models, and the VEP rheology leads to faster SI than the VP rheology. This effect of elasticity becomes smaller with increasing initial horizontal compressional stress but does not become fully negated by the initial stress of <∼50 MPa. Therefore, inclusion of elasticity with reasonable shear modulus and initial stress values results in a weaker slab, making it easier for vertically driven SI to occur when the age contrast is relatively large.

Seismic Constraints on Crustal and Uppermost Mantle Structure Beneath the Hawaiian Swell: Implications for Plume‐Lithosphere Interactions

AbstractWe infer the lithospheric structure beneath the Hawaiian Swell based on a joint inversion of ambient noise and teleseismic Rayleigh waves collected during the PLUME experiment. These combined datasets let us use Rayleigh waves with a period range of 8–50 s that provides imaging resolution for the shallow lithosphere and constrains interactions between the upwelling plume and migrating plate. We find an elongated low‐velocity anomaly beneath the lithosphere along the island chain that connects to the plume conduit, consistent with the melting region associated with a restite hotspot swell root that has mechanically eroded the base of its overlying lithosphere. It could also be consistent with a plume refracted by the overriding plate motion if the plume could manage to deeply erode the lithospheric base, as only a more viscous restite root is seen to do in thermomechanical experiments. There is also a low‐velocity body beneath the North Arch that coincides with the location of recent off‐chain volcanic fields discovered there. Its location relative to the Molokai Fracture Zone (MFZ) supports the concept that swell‐root material has preferentially spread beneath the younger side of the MFZ. We also find a clear low‐velocity anomaly associated with the uppermost ∼50 km of the fracture zone lithosphere. The relatively well‐resolved shallow lithospheric structure determined here allows us to estimate the volumes of underplating and the total flux of the rising plume material that is added beneath the Hawaiian lithosphere, which we constrain to be in the range of 0.5–0.6 km3/yr.

Chemolithoautotroph distributions across the subsurface of a convergent margin.

Subducting oceanic crusts release fluids rich in biologically relevant compounds into the overriding plate, fueling subsurface chemolithoautotrophic ecosystems. To understand the impact of subsurface geochemistry on microbial communities, we collected fluid and sediments from 14 natural springs across a ~200 km transect across the Costa Rican convergent margin and performed shotgun metagenomics. The resulting 404 metagenome-assembled genomes (MAGs) cluster into geologically distinct regions based on MAG abundance patterns: outer forearc-only (25% of total relative abundance), forearc/arc-only (38% of total relative abundance), and delocalized (37% of total relative abundance) clusters. In the outer forearc, Thermodesulfovibrionia, Candidatus Bipolaricaulia, and Firmicutes have hydrogenotrophic sulfate reduction and Wood-Ljungdahl (WL) carbon fixation pathways. In the forearc/arc, Anaerolineae, Ca. Bipolaricaulia, and Thermodesulfovibrionia have sulfur oxidation, nitrogen cycling, microaerophilic respiration, and WL, while Aquificae have aerobic sulfur oxidation and reverse tricarboxylic acid carbon fixation pathway. Transformation-basedcanonical correspondence analysis shows that MAG distribution corresponds to concentrations of aluminum, iron, nickel, dissolved inorganic carbon, and phosphate. While delocalized MAGs appear surface-derived, the subsurface chemolithoautotrophic, metabolic, and taxonomic landscape varies by the availability of minerals/metals and volcanically derived inorganic carbon. However, the WL pathway persists across all samples, suggesting that this versatile, energy-efficient carbon fixation pathway helps shape convergent margin subsurface ecosystems.

How do pre-existing weak zones and rheological layering of the continental lithosphere influence the development and evolution of intra-continental subduction?

• The depth of a preexisting weak zone is one of the key factors in initiating intra-continental subduction. • Crustal internal friction angle plays an important role in altering the deformation type of a continental collision system. • The reactivation of the suture may be an auxiliary condition for the evolution of the south-dipping continental subduction underneath the North Pamir. Intra-continental subduction is of special importance for studying the formation of intra-continental orogens, crust-mantle structural evolution, and the far-field effects of continental collision, whose mechanism is still a matter of discussion. In this work, we investigated the role of pre-existing weak zones and the continental lithospheric rheological layering in the formation and evolution of the intra-continental subduction based on a 2D finite element numerical technique. The model results indicate that the deeper the intra-continental weak zone is and the faster the convergence velocity is, the more likely it is to develop into a new intra-continental subduction. Altering the rheological strength of the overriding plate may not have a substantial impact on the intra-continental subduction mode when the depth of the pre-existing weak zone is larger than half of the lithospheric thickness. In contrast, the lithospheric rheological strength is closely related to the continental collision system’s deformation style: Models with a weaker overriding plate are inclined to delaminate continuously under collision, whereas a strong overriding plate results in the subducting plate’s roll-back. The reactivation of the suture that runs deep into the lithosphere as a result of the Indian-Asian continental collision could be one of the crucial factors controlling the formation of the south-dipping subduction under the North Pamir.