Subduction Research Articles

Plume‐Driven Subduction Termination in 3‐D Mantle Convection Models

AbstractThe effect of mantle plumes is secondary to that of subducting slabs for modern plate tectonics when considering plate driving forces. However, the impact of plumes on tectonics and planetary surface evolution may nonetheless have been significant. We use numerical mantle convection models in a 3‐D spherical chunk geometry with damage rheology to study some of the dynamics of plume‐slab interactions. Substantiating our earlier 2‐D results, we observe a range of interaction scenarios, and that the plume‐driven subduction terminations we had identified earlier persist in more realistic convective flow. We analyze the dynamics of plume affected subduction, including in terms of their geometry, frequency, and the overall effect of plumes on surface dynamics as a function of the fraction of internal to bottom heating. Some versions of such plume‐slab interplay may be relevant for geologic events, for example, for the inferred ∼183 Ma Karoo large igneous province formation and associated slab disruption. More recent examples may include the impingement of the Afar plume underneath Africa leading to disruption of the Hellenic slab, and the current complex structure imaged for the subduction of the Nazca plate under South America. Our results imply that plumes may play a significant role not just in kick‐starting plate tectonics, but also in major modifications of slab‐driven plate motions, including for the present‐day mantle.

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.

Frontal Thrust Ramp‐Up and Slow Earthquakes Due To Underthrusting of Basement High in the Nankai Trough

AbstractRecently, integrated geophysical‐geological surveys in the Nankai subduction zone in Japan have revealed that slow earthquakes repeatedly occur beneath the outer wedge of the forearc. During December 2020 to February 2021, clustered slow earthquakes propagated around the frontal thrust of the accretionary wedge. The frontal thrust ramps up from the basal décollement and slips over trench‐filling sediment along the landward edge of the Nankai trough floor. Here, the Paleo‐Zenisu ridge has been subducted beneath the inner‐outer slope border. In addition, ocean floor topography and geologic structure revealed by seismic reflection surveys completed before 2022 document that the basement of the Philippine Sea Plate beneath the frontal thrust has a seamount and a horst‐like basement high. The northern edge of the basement high is located at the ramp‐up position of the frontal thrust. The 2020–2021 clustered slow earthquakes started at the Paleo‐Zenisu ridge and propagated to the topographic highs beneath the deformation front. Considering that the relative plate convergence between the upper Amurian Plate of the Nankai forearc and the subducting Philippine Sea Plate is ∼6.0 cm/year, the basement high at the deformation front has uplifted the frontal crest of the wedge at an average rate of 2.7–5.7 mm/year for several tens to hundred thousand years. These rates are among some of the highest rock uplift rates measured in the world. The slow earthquakes in the off‐Kumano Nankai Trough in 2020–2021 are a snapshot of a “living” Nankai frontal thrust during the megathrust interseismic period.

Characteristics of Earthquake Mechanisms Along the Andaman - Nicobar Region derived from Bayesian Moment Tensor Inversion

Abstract Earthquake activities in the Andaman - Nicobar region are generated mainly by several active faults and an active subduction zone that can activate multiple hypocenter clusters. Therefore, studying the seismotectonic process from source mechanism analysis using moment tensor inversion is essential. In this research, we applied a stochastic Bayesian inversion using Grond, a framework tool for moment tensor (MT) inversion, to characterize the earthquake mechanisms and the most appropriate centroid locations. Procedurally, several of the largest earthquakes with a moment magnitude Mw > 5.5 from 2007 to 2023 were analyzed using the regional and teleseismic networks up to 2000 km from the epicenter based on 1 Hz green function model and azimuthal gap coverage. We applied 10,000-30,000 iterations of inversions with 300 bootstrap chains for a stable result. The low-frequency bandpass with a range of 0.007 - 0.03 Hz was chosen for the inversion process. There were 17 solutions for earthquake mechanisms with good misfit values of < 0.65. Based on the division of the 5 clusters with the highest seismicity, the selected earthquakes’ source mechanism is dominated by strike-slip and oblique faulting. The results were compared with existing data catalogs showing differences in hypocentre and fault orientation direction, particularly for strike-slip earthquakes in the shallow crust and mantle zones.

Structural Variation Along the Southern Hikurangi Subduction Zone, Aotearoa New Zealand, From Seismic Reflection and Retro‐Deformation Analysis

AbstractThe southern Hikurangi subduction zone exhibits significant along‐strike variation in convergence rate and obliquity, sediment thickness and, uniquely, the increasing proximity of southern Hikurangi to, and impingement on, the incoming continental Chatham Rise, an ancient Gondwana accretionary complex. There are corresponding changes in the morphology and structure of the Hikurangi accretionary prism. We combine widely spaced multichannel seismic reflection profiles with high resolution bathymetry and previous interpretations to characterize the structure and the history of the accretionary prism since 2 Ma. The southern Hikurangi margin can be divided into three segments. A northeastern segment (A) characterized by a moderately wide (∼70 km), low taper (∼5°) prism recording uninhibited outward growth in the last ∼1 Myr. Deformation resolvable in seismic reflection data accounts for ∼20 % of plate convergence, comparable with the central Hikurangi margin further North. A central segment (B) characterized by a narrow (∼30 km), moderate taper (∼8°) prism, with earlier (∼2‐∼1 Ma) shortening than segment A. Outward prism growth ceased coincidentally with development of major strike‐slip faults in the prism interior, reduced margin‐normal convergence rate, and the onset of impingement on the incoming Chatham Rise to the south. A southwestern segment (C) marks the approximate southern termination of subduction but widens to ∼50 km due to rapid outward migration of the deformation front via fault reactivation within the now‐underthrusting corner of the Chatham Rise. Segment C exhibits minimal shortening as margin‐normal subduction velocity decreases and plate motion is increasingly taken up by interior thrusts and strike‐slip faults.

Copper migration and enrichment in the mantle wedge: Insights from orogenic peridotites and pyroxenites

The refertilized mantle wedge is an important source of ore-forming metals in subduction-related Cu–Au deposits. However, the source and migration of Cu in the mantle wedge are poorly constrained. Here, we present a combined study of the Cu elemental and isotopic compositions (δ65Cu) as well as Fe3+/∑Fe ratios on a well-characterized suite of the Mg-lherzolites, Fe-rich peridotites and pyroxenites from the Bohemian Massif in a Variscan subduction zone. The Mg-lherzolites represent melting residues moderately affected by metasomatism of slab-derived melts/fluids. The Fe-rich peridotites and pyroxenites are, respectively, results of Mg-lherzolite-melt reaction and crystalline products of evolved melts in the lithospheric mantle.The peridotites and pyroxenites display higher Fe3+/∑Fe ratios (0.14–0.56) than the cratonic peridotites and mid-ocean ridge basalts, indicating that the mantle wedge was oxidized by slab-derived components. The Mg-lherzolites have relatively low Cu contents (7.35–39.6 μg/g) and normal mantle-like δ65Cu values (−0.13 to 0.26 ‰), suggesting an insignificant slab-to-mantle wedge Cu transfer. In contrast, the Fe-rich peridotites and pyroxenites have variable but overall high Cu contents (37.0–513 μg/g) and heavier δ65Cu values (up to 0.83 ‰). These signatures are ascribed to secondary sulfide precipitation from the evolved Cu- and 65Cu-enriched melts, which were produced by reaction of Mg-lherzolites with subduction-related oxidative SiO2‐undersaturated basaltic melts. Such melt-peridotite reaction at high oxygen fugacity can cause the partial oxidative decomposition of primary sulfides in peridotites with the preferential release of 65Cu into the evolved melts. Our results thus demonstrate that oxidative melt-rock reaction can result in the Cu migration, redistribution and its local enrichment in the mantle wedge, which may serve as an important source of Cu for subduction-related porphyry copper deposits.

Paleostress inversion of fault-slip in the Ayhan Basin, Kırşehir Block (Turkey) since the Late Cretaceous: Insights into the Neotethys closure

Kırşehir Block, a large triangular continental domain exposed in Central Anatolia, is surrounded by ophiolitic melanges that mark former subduction zones associated with the closure of the northern and southern branches of the Neotethys Ocean. Notably, the Ayhan Basin, located in the inner part of the block, is one of the few locations containing a geological record from the Upper Cretaceous to the Pliocene. It provides a unique opportunity to examine the crustal deformation and the temporal and spatial stress distribution over the Central Anatolian orogen. This study encompasses a comprehensive analysis of the Ayhan Basin’s stratigraphy and structure, combined with the paleostress inversion solutions derived from over 400 fault-slip data. The findings reveal that the Ayhan Basin underwent three distinct tectonic phases: (1) ∼ NE-SW extensional phase from Late Cretaceous to Lutetian, which refers to the main basin-forming stress configuration. Importantly, this extensional phase is crucial to delineate the spatial extent of the concurrent collision in the north. (2) ∼ N-S compressional phase during the Lutetian, mainly governed by the shortening of the Kırşehir Block. This compressional phase caused the inversion of the inherited normal faults in the basin; and (3) ∼ NE-SW extensional phase active since the Late Miocene was accompanied by volcanism and uplift in the region. The Cyprus Slab’s southward retreat might be this extension’s triggering mechanism. Consequently, the temporal and spatial extent of the tectonic phases and their resultant deformations are crucial for understanding the underlying geodynamic processes and establishing their timing.

Fluid-mediated Cu and Zn isotope fractionation in subduction zones and implications for arc volcanism: Constraints from high pressure veins within eclogites in the Dabie Orogen

High-pressure (HP) veins in eclogites are the products of fluid-rock interaction and provide insight into the composition and evolution of fluids in subduction zones. We here present the Cu and Zn isotope data for different types of HP veins and their host eclogites from Ganghe and Hualiangting in the Dabie Orogen to reveal the behavior of Cu and Zn isotopes during fluid-rock interaction and fluid evolution. The HP veins include omphacite-epidote (Omp-Ep), epidote-quartz (Ep-Qtz), and kyanite-epidote-quartz (Ky-Ep-Qtz) veins. The Omp-Ep veins first crystallized from eclogite-derived, solute-rich fluids with the Ep-Qtz and Ky-Ep-Qtz veins successively crystallizing from the residual fluids after the Omp-Ep vein formation. The early Omp-Ep veins have variably lower δ65Cu (−1.32 to −1.12‰ at Ganghe, −0.70 to −0.66‰ at Hualiangting) but higher δ66Zn (0.36 to 0.38‰ at Ganghe, 0.40 to 0.41‰ at Hualiangting) relative to the host eclogites (δ65Cu: −0.71 vs. 1.84‰, δ66Zn: 0.22 vs. 0.33‰), indicating CuZn isotope fractionation during slab dehydration in subduction zones with the lighter Cu and heavier Zn isotopes preferentially entering the vein-forming fluids from the eclogites. Systematic increase of δ65Cu from the Omp-Ep (−0.70 to −0.66‰) through Ep-Qtz (0.01‰) to Ky-Ep-Qtz veins (0.46 to 0.95‰) can be attributed to isotope fractionation induced by redox changes during the evolution of metamorphic fluids, as revealed by the negative correlations of δ65Cu with redox-sensitive ratios of Fe3+/ΣFe and (Eu/Eu*)N in those veins. Additionally, the higher δ66Zn of the Omp-Ep veins relative to the Ky-Ep-Qtz veins can be explained by equilibrium isotope fractionation between crystallized minerals and evolved metamorphic fluids. Our results thus demonstrate that CuZn isotope fractionation occurred during the evolution of slab-derived metamorphic fluids in subduction zones. Binary mixing calculations show that fluids derived from dehydration of mafic rocks in subducted slabs, represented by the multistage HP veins in the present study, cannot account for the heavier Cu and lighter Zn isotope compositions in most arc magmas than in mid-ocean ridge basalts. This offset can be resolved by addition of forearc serpentinite-derived fluids enriched in heavy Cu and light Zn isotopes into the mantle source of arc magmas.

Evolving Subduction Zone Thermal Structure Drives Extensive Forearc Mantle Wedge Hydration

AbstractHydration of the subduction zone forearc mantle wedge influences the downdip distribution of seismicity, the availability of fluids for arc magmatism, and Earth's long term water cycle. Reconstructions of present‐day subduction zone thermal structures using time‐invariant geodynamic models indicate relatively minor hydration, in contrast to many geophysical and geologic observations. We pair a dynamic, time‐evolving thermal model of subduction with phase equilibria modeling to investigate how variations in slab and forearc temperatures from subduction infancy through to maturity contribute to mantle wedge hydration. We find that thermal state during the intermediate period of subduction, as the slab freely descends through the upper mantle, promotes extensive forearc wedge hydration. In contrast, during early subduction the forearc is too hot to stabilize hydrous minerals in the mantle wedge, while during mature subduction, slab dehydration dominantly occurs beyond forearc depths. In our models, maximum wedge hydration during the intermediate phase is 60%–70% and falls to 20%–40% as quasi‐steady state conditions are approached during maturity. Comparison to global forearc H2O capacities reveals that consideration of thermal evolution leads to an order of magnitude increase in estimates for current extents of wedge hydration and provides better agreement with geophysical observations. This suggests that hydration of the forearc mantle wedge represents a potential vast reservoir of H2O, on the order of 3.4–5.9 × 1021 g globally. These results provide novel insights into the subduction zone water cycle, new constraints on the mantle wedge as a fluid reservoir and are useful to better understand geologic processes at plate margins.

Co-existence of magmatic- and metamorphic-hydrothermal gold deposits in the northern margin of the North China Craton and their genetic link to Paleo-Pacific plate subduction: Insights from in-situ analyses of sulfur isotopes and trace elements

Mid-Mesozoic gold deposits, whose metallogeny was coincident with subduction of the Paleo-Pacific plate, are widely distributed in the northern margin of the North China Craton (NCC). The genesis of these deposits is currently debated with magmatic- and metamorphic-hydrothermal models. In-situ analyses of sulfur isotopes and trace elements provide an opportunity to clarify such a debate. Here, the two techniques were used to determine the trace elements and sulfur isotopes of sulfides from three mid-Mesozoic gold deposits (Jiapigou, Jinchangyu, and Jinchanggouliang) in the northern margin of NCC. The Co/Ni ratios of pyrites suggest that Jiapigou and Jinchanggouliang deposits are of magmatic-hydrothermal origin, whereas the Jinchangyu deposits receive ore-forming materials mainly from crustal-hydrothermal fluids. In addition, gold-bearing sulfides display positive (5.60 ± 1.87 ‰, SD), near-zero (−0.31 ± 1.27 ‰, SD) and negative δ34SV-CDT values (−3.84 ± 2.97 ‰, SD) in the Jiapigou, Jinchangyu, and Jinchanggouliang, respectively. The results suggest that mid-Mesozoic gold deposits in the northern margin of NCC have variable sources of sulfur and metals. The Jiapigou deposit, representative of the arc settings, received sulfur and metals mainly from the subducting ocean crust. The Jinchanggouliang deposit, representative of cratonic destruction settings, received sulfur and metals mainly from the metasomastized mantle. The Jinchangyu deposit, representative of the back-arc settings, received sulfur and metals mainly from the metamorphic basement. We infer that subduction of the Paleo-Pacific plate in the mid-Mesozoic not only led to direct metallogeny of magmatic-hydrothermal gold deposits at arc settings and cratonic destruction settings, but also triggered the release of ore-forming fluids from the basement at back-arc, favoring metallogeny of crustal-hydrothermal gold deposits. Thus, mid-Mesozoic gold deposits in the northern margin of NCC cannot be exclusively explained by a single model.

Oxidized Sulfur Species in Slab Fluids as a Source of Enriched Sulfur Isotope Signatures in Arcs

AbstractRecycling of oxidized sulfur from subducting slabs to the mantle wedge provides simultaneous explanations for the elevated oxygen fugacity (fO2) in subduction zones, their high hydrothermal and magmatic sulfur outputs, and the enriched sulfur isotopic signatures (i.e., δ34S > 0‰) of these outputs. However, a quantitative understanding of the abundance and speciation of sulfur in slab fluids consistent with high pressure experiments is lacking. Here we analyze published experimental data for anhydrite solubility in H2O‐NaCl solutions to calibrate a high‐pressure aqueous speciation model of sulfur within the framework of the deep earth water model. We characterize aqueous complexes, required to account for the high experimental anhydrite solubilities. We then use this framework to predict the speciation and solubility of sulfur in chemically complex fluids in equilibrium with model subducting mafic and ultramafic lithologies, from 2 to 3 GPa and 400 to 800°C at log fO2 from FMQ‐2 to FMQ+4. We show that sulfate complexes of calcium and sodium markedly enhance the stability of sulfate in moderately oxidized fluids in equilibrium with pyrite at fO2 conditions of FMQ+1 to +2, causing large sulfur isotope fractionations up to 10‰ in the fluid relative to the slab. Such fluids could impart oxidized, sulfur‐rich and high δ34S signatures to the mantle wedge that are ultimately transferred to arc magmas, without the need to invoke 34S‐rich subducted lithologies.

Present day mantle structure from global mantle convection models since the Cretaceous

SUMMARY Using forward mantle convection models starting at 140 Ma, and assimilating plate reconstructions as surface velocity boundary condition, we predict present-day mantle structure and compare them with tomography models, using geoid as an additional constraint. We explore a wide model parameter space, such as different values of Clapeyron slope and density change across 660 km, density and viscosity of the thermochemical piles at the core–mantle boundary (CMB), internal heat generation rate, and model initiation age. We also investigate the effects of different strengths of a weak layer below 660 km and weaker asthenosphere and slabs. Our results suggest that slab structures at different subduction zones are sensitive to the viscosity of the asthenosphere, strength of slabs, values of Clapeyron slope and the density and viscosity of the thermochemical piles, while different internal heat generation rates do not affect the slab structures. We find that with a moderately weak asthenosphere ($10^{20}$ Pa·s) and strong slabs, the predicted slab structures are consistent with the tomography models, and the observed geoid is also matched well. Moreover, our models successfully reproduce the degree-2 structure of the lower mantle beneath Africa and the Pacific, also known as Large Low Shear Velocity provinces (LLSVPs). A moderate Clapeyron slope of −2.5 MPa K−1 at 660 km aids in slab stagnation while higher values result in massive slab accumulation at that depth, ultimately leading to slab avalanches. We also find that the convective patterns in the thermal and thermochemical cases with slightly denser LLSVPs are similar, although the geoid amplitudes are lower for the latter. However, with more dense LLSVPs, the slabs cannot perturb them and no plumes are generated. Plumes arise as thermal instabilities from the edges of the LLSVPs, when cold and viscous slabs perturb them. While our predicted plume locations are consistent with the observed hotspot locations, matching the plume structures in tomography models is difficult. These plumes are essential in fitting the finer features of the observed geoid. In longer-duration models, more voluminous subducted material reaches the CMB, which tends to erode the LLSVPs significantly, and yields a poor fit to the observed geoid. Our results suggest that with the presence of a thin, moderately weak layer below 660 km, a slightly dense LLSVP, and Clapeyron slope of −2.5 MPa K−1, the velocity anomalies in seismic tomography and the long-wavelength geoid can be matched well. One of the limitations of our models is that the assimilated plate motion history may be too short to overcome arbitrary initial conditions effects. Also, assimilated true plate velocities in our models may not represent the true convective vigour of the Earth.

Deformation mechanisms and slip behaviors of tectonically deformed conglomerates from the Central Apennines fold-and-thrust belt: Implications for shallow aseismic and seismic slip

Fold-and-thrust belts and accretionary prisms exhibit considerable variability in slip behaviors along active faults. While numerous studies have investigated fault mechanics and slip behaviors in carbonate-, shale-, and sandstone-hosted faults as well as clay-rich portions of shallow accretionary prisms, the deformation mechanisms and slip behaviors of deformed, thrust-top molasse-type conglomerates remain poorly understood. To fill this knowledge gap, we conducted structural and microstructural analyses on exhumed and deformed conglomerates in the footwall of an out-of-sequence thrust from the Central Apennines fold-and-thrust belt (Italy). Our findings aim to constrain conglomerate deformation mechanisms and infer possible slip behaviors. We observed flattened pebbles and an intense foliation comprising densely spaced, thrust-parallel, clay-rich stylolites that indicate slow deformation attributed to carbonate dissolution during aseismic creep. In contrast, quartz clasts with calcite micro veins and micrometer-thick slip increments in slickenfibers suggest fast brittle failure and fluid overpressure in competent pebbles, while the matrix continues to deform aseismically through pressure solution. The general structure is similar to block-in-matrix textures observed in tectonic mélanges from exhumed subduction zones and accretionary prisms. Our results provide new insights into the slip behaviors of shallow (<∼1 km) portions of fold-and-thrust belts, complementing existing understandings of deformation mechanisms and slip behaviors across different depth ranges in fold-and-thrust belts and accretionary prisms.

Characteristics of Metamorphic Rocks and the Deformation and Metamorphism Period in Northern Lhasa,Tibet

The study of metamorphic rocks and deformation and metamorphic features in Gangdise-Lhasa metamorphic zone is helpful for the establishment of regional tectonic framework and the discussion of the evolution history of the region, which will provide new basic data for the study of geological structure evolution in the Qinghai-Tibet Plateau. There are few studies on metamorphic rocks in the northern area of Lhasa, Tibet. This paper summarizes and studies them from the perspective of petrography, and briefly discusses the stages of deformation and metamorphism. The area is divided into two metamorphic belts, which are divided into Yebagou-Gandansi metamorphic belt and Shexing-Binjiaolin metamorphic belt from south to north. The types of metamorphic rocks are slate, kyanite and schist. The dynamic metamorphism is manifested as brittle deformation and ductile deformation. The dynamic metamorphic rock belt developed in the area is several kilometers long, which basically belongs to the low green schist facies belt. The metamorphism in the area occurs under the geological background of the northward subduction, subduction and collision of the Neotethys oceanic crust plate. And with the change of tectonic environment and heat flow, different metamorphisms are produced in different geological development stages.

The 2.03–2.02 Ga arc–back-arc mafic magmatism from the Zhongxiang Complex, northern Yangtze Block: Implications for Columbia supercontinent assembly

To better understand Paleoproterozoic tectonic processes in the Yangtze Block during the assembly of the Columbia supercontinent, an integrated geochronological and geochemical study of amphibolite facies mafic rocks exposed in the Zhongxiang Complex in the northern Yangtze Block is presented. Based on trace element and rare earth element (REE) characteristics, these mafic rocks can be divided into two coeval groups. The MORB-like mafic rocks show flat REE patterns and are enriched in large-ion lithophile elements (e.g. Rb, Ba, K, Sr). Their zircon ɛHf(t) values range from + 3.43 to + 4.90. The other mafic rocks are arc-like and have elevated REE contents and display relatively light REE enriched patterns, as well as depletion in high field strength elements (e.g. Nb, Ta, Ti). The zircon ɛHf(t) values of the latter group range from + 0.17 to + 2.68. Zircon U-Pb dating yielded a weighted mean 207Pb/206Pb age of 2029 ± 15 Ma for the MORB-like mafic rocks, and 2023 ± 14 Ma for the arc-like mafic rocks, respectively. The elemental and isotopic characteristics indicate that the MORB-like mafic rocks may be derived from a depleted asthenospheric mantle source, whereas its counterpart with arc-like signatures originated from a lithospheric mantle source. The low (La/Yb)CN (1.10–2.84) and Sm/Yb (0.99–1.29) ratios indicate that the mafic rocks are formed by partial melting of spinel lherzolite mantle. Combined with other geological observations, these mafic rocks are inferred to constitute part of a mid-Paleoproterozoic continental arc–back-arc basin system. The 2.03–2.02 Ga mafic rocks from the northern Yangtze Block indicate that the continental margin experienced subduction in response to assembly of the Columbia supercontinent.

Geochemistry and zircon U-Pb geochronology of diabase, gabbros, anorthosites and ultramafic rocks from the South Andaman Island ophiolite suite on the South-Eastern margin of the Indian plate

The South Andaman Island is located on the southeastern fringe of the Indian plate. Tectonically, the islands are associated with the Burma-Andaman-Java subduction zone of the Indonesian arc system. Several Mesozoic ophiolitic exposures are preserved in the Burma-Andaman-Java subduction zone. This study reports new petrological, geochemical, and zircon U-Pb geochronological data for the diabases, gabbros, anorthosites, and harzburgites of the ophiolite suite exposed in the South Andaman Island. The diabase is composed of plagioclase and clinopyroxene exhibiting ophitic textures. Anorthosite mainly contains plagioclase (> 90 vol%), pyroxenes and Fe-oxides. Gabbros have plagioclases and clinopyroxenes occurring as the major minerals. Harzburgites consist of olivine, orthopyroxene, clinopyroxene and spinel. Geochemical investigation reveals that these rocks were formed through crystal fractionation and accumulation processes from partial melt generated from a spinel lherzolite source. Gabbro, diabase, and anorthosites are formed by ∼1–5% melting of the, whereas depleted harzburgites were generated by slightly higher degrees of partial melting (∼10–20%) of the spinel lherzolite source. Trace element modelling reveals the influence of mid-oceanic ridge basalt (MORB)-melt interaction and subduction-related mantle alterations.The anorthosite and gabbro samples yield zircon 206Pb / 238U weighted mean ages of 93.2 ± 1.0 Ma and 442.5 ± 4.0 Ma, 383.9 ± 3.3 Ma respectively. However, the harzburgite sample records multiple age population including Paleo-Mesoproterozoic grains ranging in age from 1216 to 1931 Ma as well as Neoproterozoic and Phanerozoic grains of 607 Ma, 594 Ma, 237 Ma, 105 Ma, and 85 Ma suggesting a long-lived supra subduction zone in the Burma-Andaman-Java subduction zone where the South Andaman ophiolite suite was formed, evolved, and was emplaced.

Inversion of fluid-release rates from episodic tremor and slip signals in subduction zones via a coarse-grained reaction diffusion model

Episodic Tremor and Slip (ETS) events showcase dynamic interactions of oscillatory slow slips and tremors deep within subduction zones and offer a window into Earth's internal dynamics. However, the exact mechanisms driving these events remain unresolved. This study proposes a novel approach that goes beyond traditional explanations focused on fluid pressure from mineral dehydration. Existing models often neglect the intricate interplay between fluid and rock pressures across various depths and potential fluid sources. This calls for a more comprehensive understanding of how fluid release from reactions interacts with rock deformation. The present formulation captures the interplay between fluid and solid pressures providing a more rigorous picture of ETS events. It employs a minimalistic and efficient approach based on integrating dehydration reactions. The model thereby develops a generic framework for mineral dehydration, offering an enhanced perspective of the underlying processes without the need to trace down to specific minerals. It allows a refined fit to GPS data by including high-frequency components from linear and nonlinear stability analyses, giving rise to improved correlation coefficients. Through the inclusion of the dynamic interplay between fluid and rock pressure diffusion within subduction zones, we propose a unified model of ETS events.

Formation of silicic crust on early Earth and young planetary bodies in an Iceland-like setting

Four billion years ago, Earth’s surface was covered with a thick basaltic crust and, similar to other rocky planets and the Moon, parts of this basaltic crust underwent partial melting to produce silicic crust. On Earth, silicic crust grew into the continents. An analogue of a pre-subduction early Earth and rocky planetary bodies is found in Iceland, where poorly understood granitic rocks are encased in thick basaltic crust away from any active subduction zone. Here, we investigate these Icelandic granitoids to understand the mechanisms that generated the Earth’s oldest continents, and silicic rocks on Mars and Venus. New geochemical data from silicic Icelandic intrusions, show that the granitoids are formed through partial melting of the Icelandic crust, but are compositionally unlike early Earth continental material. We show that intra-crustal partial melting of basalt can potentially form silicic material on other planetary bodies but cannot produce the first continents on Earth.

Middle Triassic and middle Permian radiolarians from the Kamitaki Complex in the Sasayama area, Hyogo Prefecture, Southwest Japan: Evidence for Triassic plate subduction along the eastern margin of Paleo-Asia

The Kamitaki Complex, situated in the Sasayama area in Southwest Japan, has long been presumed to be a Permian subduction-related accretionary complex based on correlations from previous studies. However, because of the lack of fossil evidence, the exact age of the complex remained uncertain for a long time period. To address this gap in knowledge, a geological survey and microfossil mapping were conducted in the Kamitaki Complex to determine its age and geological context.A geological survey revealed that the Kamitaki Complex mainly consists of clastic rocks, and a mixture of sandstone, basalt, and chert blocks within the mudstones. The Kamitaki Complex is tectonically intercalated into the Lopingian (Late Permian) accretionary complex of the Ultra-Tamba Terrane and Late Triassic accretionary complex of the Tamba Terrane. The lithological and structural characteristics of the Kamitaki Complex confirm that it is an accretionary complex. Microfossil mapping yielded depositional ages, with radiolarian fauna such as Eptingium nakasekoi, Pseudostylosphaera japonica, Cryptostephanidium japonicum, and Oertlispongus cf. diacanthus identified in mudstones suggesting an Anisian (early Middle Triassic) age. In contrast, radiolarian fauna found in cherts, including Pseudoalbaillella? aff. longicornis, and Follicucullus cf. porrectus, indicate an early Capitanian (late Guadalupian, middle Permian) depositional age. These findings suggest that the Kamitaki Complex records a trenchward migration of the oceanic plate in a pelagic environment from the early Capitanian and an accretion at the trench during the Anisian period.Conventionally, the plate boundary between the Panthalassa and Paleo-Asia during the latest Permian to Middle Triassic was of the transform type, primarily because no subduction-related accretionary complexes from this period have been identified in the Japanese Islands. However, the discovery of Kamitaki Complex, an Anisian accretionary complex, provides evidence of Middle Triassic subduction activity along the eastern margin of Paleo-Asia. According to the internal structure and age polarity of the Ultra-Tamba and Tamba terranes, an accretionary complex developed over a prolonged period (approximately 120 million years) in a tectonic setting that persisted along the eastern margin of Paleo-Asia from the late Guadalupian to the earliest Cretaceous period.

Ring-shaped seismicity structures in the region of Kamchatka: possible preparation for great earthquake

We have been studying some seismicity characteristics in the region of Kamchatka region. It was established that two pairs of ring-shaped seismicity structures at depths of 0-33 and 34-70 km have been formed here in 1973-2018 prior to large earthquakes (05.12.1997 г., Mw=7.8 and 20.12.2018, Mw=7.3). Ring-shaped structures are characterized by threshold magnitude values (Mt1 and Mt2 correspondingly) and also big axes length (L1 and L2). Ring-shaped structures were identified in the area to the east and northeast of Petropavlovsk-Kamchatskiy city, where no great earthquakes have occurred during more than 60 years (ring-shaped structure was identified also in the depth range of 71-110 km here). Earlier, correlation dependences of Mt1 and Mt2 parameters on Mw values of major earthquakes for the west of Pacific were creating. Using these dependencies, the magnitude of possible great earthquake in this area was estimated as Mw=8.6±0.2. Besides that, a conclusion on the source depth of such event was drawing on the base of the analysis of three ring-shaped structures. The reasons for ring-shaped structures formation in different depth ranges of the subduction zones are discussing.