Subduction Style Research Articles

Thermo‐Mechanical Effects of Microcontinent Collision on Ocean‐Continent Subduction System

AbstractMicrocontinents are globally recognized as continental regions partially or entirely surrounded by oceanic lithosphere. Due to their positioning, they may become entangled in subduction zones and undergo either accretion or subduction. High‐pressure metamorphism in subducted continental rocks supports the idea that microcontinents can be subducted, regardless of their low densities. In this study, we used 2D numerical models to simulate collision of microcontinents with different lengths located at various distances from the upper plate in a subduction system characterized by different convergence velocities, in order to examine their effects on the thermo‐mechanical evolution of subduction systems. Specifically, we analyzed the conditions that favor subduction or accretion of microcontinents and investigated how their presence affects the thermal state within the mantle wedge. Our results reveal that the presence of microcontinents can lead to four styles of subduction: (a) continuous subduction; (b) continuous subduction with jump of the subduction channel; (c) interruption and reinitiation of the subduction; (d) continental collision. We discovered that longer microcontinents and higher velocities of the subducting plate contrast a continuous subduction favoring accretion, whereas farther initial locations from the upper plate and higher velocities of the upper plate favor the subduction of the microcontinent. In addition, we observed that the subduction style has direct effects on the thermal state, with important implications for the potential metamorphic conditions recorded by the subducted continental rocks. In particular, models characterized by parameters that favor the subduction of a larger amount of continental material from the microcontinent exhibit warm mantle wedges.

Did Along-Strike Changes in Continental Subduction Styles Occur in the Dabie-Sulu Orogenic Belt?

The Triassic collision between the North China Block (NCB) and the South China Block (SCB) formed the Dabie-Sulu orogenic belt, renowned for its extensive ultra-high-pressure metamorphic rocks and the prominent Tan-Lu Fault. Geological and geophysical data reveal significant along-strike differences in both shallow and deep structures within the belt, raising a critical question: What tectonic processes drive these structural variations? Considering that the Triassic collision involved pre-collisional accreted microcontinents and along-strike variations in convergence velocities, this study investigates the influence of microcontinental width, convergence rate, and initial Moho temperature on continental subduction styles using numerical modeling. The models incorporate pre-collisional accreted microcontinent(s) and reveal a two-stage subduction evolution: an initial transition from one-sided to two-sided continental subduction, followed by subduction polarity reversal. High initial Moho temperatures, rapid convergence rates, and wide accreted microcontinents promote the development of two-sided subduction, characterized by an initially vertical interface that gradually inclines towards the pro-plate as convergence progresses. These subduction styles significantly influence crustal suture migration. One-sided subduction results in minimal horizontal displacement, whereas two-sided subduction and polarity reversal lead to substantial horizontal shifts. Based on the modeling results, this study proposes a new evolutionary model for the Dabie-Sulu orogenic belt during the NCB-SCB collision. The model effectively explains along-strike structural differences, such as the inconsistent tectonic trends on either side of the Tan-Lu Fault and the opposite dipping directions of high-velocity mantle anomalies observed in geophysical profiles. Furthermore, the proposed model sheds light on the formation and evolution of the Tan-Lu Fault.

Subduction Style at Different Stages of Geological History of the Earth: Results of Numerical Petrological-Thermomechanical 2D Modeling

In this article we examine the effects of impact of slab rocks eclogitization on the subduction regime under the continent. Eclogitization of rocks in high-pressure metamorphic complexes occurs only in the areas of penetration of hydrous fluid. In the absence of hydrous fluid, the kinetic delay of eclogitization preserves low-density rocks under P‒T conditions of eclogite metamorphism, delaying the weighting of a slab and reducing the efficiency of the slab-pull mechanism which contributes to the steep subduction into the deep mantle. The results of numerical petrological-thermomechanical 2D modeling of subduction under the continent in a wide range of eclogitization parameters of oceanic crust rocks (discrete eclogitization) are presented. The effects of a lower kinetic delay of eclogitization in the water-bearing basalt layer, compared to the drier underlying gabbro layer, have been tested. Based on results of 112 numerical experiments with 7 variants of eclogitization ranges (in range 400–650°C for basalt and 400–1000°C for gabbro) at different potential mantle temperatures (ΔT = 0–250°C, above modern value), and steep, flat and transitional subduction regimes were identified. The mode of steep subduction occurs under modern conditions (ΔT = 0°C) with all ranges of eclogitization. Here it is characterised by an increase in the angle of subduction of the slab as the plate descends, and above the boundary of the mantle transition zone there is a flattening or and then tucking of the slab. Subduction is accompanied by the formation of felsic and mafic volcanics and their plutonic analogues. At elevated temperatures of the mantle (ΔT≥150°С) and discrete eclogitization over a wide range, the flat subduction regime is observed with periodic detachments of its steeper frontal eclogitized part. The flat subduction regime is accompanied by significant serpentinization of the mantle wedge and episodic, scarce magmatism (from mafic to felsic), which occurs at a significant distance (≥500 km) from the trench. During the transition regime, which is also realised in models with elevated mantle temperatures, there is a characteristic change occurs from flat to steep subduction, resulting in a stepped shape of the slab. As the kinetic shift of eclogitisation increases, flat subduction develops. An increase in the thickness of the continental lithosphere from 80 km to 150 km contributes to the implementation of steep subduction, while the influence of the convergence rate (5–10 cm/year) is ambiguous. Discrete eclogitization of thickened oceanic crust and depletion of lithospheric mantle in the oceanic plate are the main drivers of flat subduction. In modern conditions, their influence becomes insignificant due to the decrease in the thickness of the oceanic crust and the degree of depletion of the oceanic mantle lithosphere. As a result, the less frequent flat movement of slabs is determined by other factors.

Anisotropy and XKS splitting from geodynamic models of double subduction: testing the limits of interpretation

SUMMARY The analysis of the splitting signature of XKS phases is crucial for constraining seismic anisotropy patterns, especially in complex subduction settings such as outward-dipping double subduction. A natural example of this is found in the Central Mediterranean, where the Apennine and the Dinaride slabs subduct in opposite directions, with the Adriatic plate separating them. To assess the capability of XKS-splitting analysis in revealing anisotropic seismic properties, such as fast polarization directions and shear wave anisotropy (in per cent), we use three-dimensional numerical geodynamic models combined with texture evolution simulations. In these models, two identical outward-dipping oceanic plates are separated by a continental plate. Using the full elastic tensors – directly derived from the texture evolution simulations – we compute anisotropic seismic properties and synthetic teleseismic waveforms. From these waveforms synthetic observables are determined, including apparent splitting parameters (fast polarization directions and delay times) and splitting intensities. Based on these observables, we (1) derive models for a single anisotropic layer (one-layer model), (2) identify regions with significant depth-dependent anisotropic seismic properties, and (3) perform inversions at selected locations in terms of two anisotropic layers (two-layer model). We consider two geodynamic models: one with a strong (M1) and one with a weak (M2) continental plate. Model M1 exhibits significant retreat of the subducting plates with no horizontal stretching of the continental plate, whereas Model M2 shows less retreat, substantial horizontal stretching, and detachment of the subducting plates. These different subduction styles result in distinct flow and deformation patterns in the upper mantle, which are reflected in the anisotropic seismic properties. In Model M1, the fast polarization directions below the continental plate are predominantly trench-parallel, whereas in Model M2, they are mostly trench-normal. In most regions of both models, the one-layer models are sufficient to resolve the anisotropic seismic properties, as these properties are nearly constant with depth. However, for both models, we identify some isolated regions – primarily near the tips of the subducting plates and beneath the continental plate – where fast polarization directions exhibit significant variations with depth. Inverting the apparent splitting parameters in these regions yields multiple two-layer models at each location that excellently fit the observables. However, their anisotropic seismic properties can vary significantly, and not all these two-layer models adequately approximate the true depth variations. This ambiguity can be partially reduced by selecting two-layer models in which the summed shear wave anisotropy closely matches that of one of the one-layer models, as these models better capture the true variations.

Deciphering Subduction Polarity During Ancient Arc‐Continent Collisions

AbstractThe closure of an ancient ocean basin via oceanic arc‐continent collision has two subduction styles with opposite polarities, which may proceed via subduction polarity reversal (SPR) or a subduction zone jump (SZJ). Interpreting the geometry or kinematic evolution of ancient collisional zones, especially the original subduction polarity, can be challenging. Here we used 2D thermo‐mechanical modeling to investigate the dynamic evolution process of SPR versus SZJ. Our modeling predicts different structural, topographic, magmatic, and basin histories for SPR and SZJ, which can be compared against, and help interpret, the geologic record past sites of oceanic closure during collisional orogens. Our results match geologic observations of past collisions in Kamchatka, eastern Russia, and the Banda Arc, eastern Indonesia, and thus our results can help effectively decode the evolutionary history of past arc‐continent collisions.

Subduction Styles at Different Stages of Geological History of the Earth: Results of Numerical Petrological–Thermomechanical 2D Modeling

Subduction Styles at Different Stages of Geological History of the Earth: Results of Numerical Petrological–Thermomechanical 2D Modeling

Tectonic Evolution of the Condrey Mountain Schist: An Intact Record of Late Jurassic to Early Cretaceous Franciscan Subduction and Underplating

AbstractThe Klamath Mountains in northern California and southern Oregon are thought to record 200+ m.y. of subduction and terrane accretion, whereas the outboard Franciscan Complex records ocean‐continent subduction along the North American margin. Unraveling the Klamath Mountains' Late Jurassic history could help constrain this transition in subduction style. Key is the Mesozoic Condrey Mountain Schist (CMS), comprising, in part, a subduction complex that occupies a structural window through older, overlying central Klamath thrust sheets but with otherwise uncertain relationships to more outboard Klamath or Franciscan terranes. The CMS consists of two units (upper and lower), which could be correlated with (a) other Klamath terranes, (b) the Franciscan, or (c) neither based on regional structures and limited extant age data. Upper CMS protolith and metamorphic dates overlap with other Klamath terranes, but the lower CMS remains enigmatic. We used multiple geochronometers to constrain the timing of lower CMS deposition and metamorphism. Maximum depositional ages (MDAs) derived from detrital zircon geochronology of metasedimentary rocks are 153–135 Ma. Metamorphic ages from white mica K‐Ar and Rb‐Sr multi‐mineral isochrons from intercalated and coherently deformed metamafic lenses are 133–116 Ma. Lower CMS MDAs (<153 Ma) predominantly postdate other Klamath terrane ages, but subduction metamorphism appears to start before the earliest coherent Franciscan underplating (ca. 123 Ma). The lower CMS thus occupies a spatial and temporal position between the Klamath Mountains and Franciscan and preserves a non‐retrogressed record of the Franciscan Complex's early history (>123 Ma), otherwise only partially preserved in retrogressed Franciscan high grade blocks.

Long-lived Northern Hemisphere convergence systems driven by upper-mantle thermal inhomogeneity

Abstract Plate reconstructions reveal that two secular centers of convergence formed beneath eastern Eurasia and North America no later than 200 Ma. The cause of these convergence centers, which featured flat subduction, slab stagnation, and/or continental margin subduction, remains uncertain. Here, we propose that upper-mantle thermal inhomogeneity, particularly an anomalously cool Northern Hemispheric upper mantle, was a fundamental driver of this long-lived convergence. By considering the pattern of observed thermal inhomogeneity, our numerical models show that flow-induced asymmetrical subduction will tend to develop toward cold mantle domains, even when the subducting plate is buoyant. The models can reproduce the diverse subduction styles observed in the Northern Hemisphere by including proposed pre-subduction plate distributions and/or properties.

Orogen-parallel discontinuity of the Apennines subduction zone in Southern Italy as seen from mantle wedge seismic structure

We investigate the seismic structure of the mantle wedge of the Apennines subduction zone (Central Mediterranean) using teleseismic receiver function (RF). We inverted RF for both isotropic and anisotropic properties of the mantle wedge, from below the overriding Moho to the “plate boundary”, i.e. the interface that separate the slab from the mantle wedge. Given the distribution of the seismic network, we are able to map out the change in the elastic properties at the transition between southern apennines and the Calabrian arc, given by the change in the subduction style (i.e from the subduction of continental materials to oceanic plate). We found that the anisotropy in the mantle wedge is similar between all seismic stations, generally highly anisotropic (> 10%), with a direction of the symmetry axis that rotates clockwise from North to South, following the Calabrian arc geometry and likely indicating the mantle flow driven by the slab retreat. The elastic properties of the subducted crust are more heterogeneous. To the North, the subducted crust shows a highly anisotropic (> 10%) behavior, and it occurs at larger depth (around 70 km depth), where to the South anisotropy is less intense (around 7%) and the subducted crust is shallower (around 60 km depth). These results point out a change in the subduction style that can be given by either a change in the metamorphic phase (more evolved blueschist facies stage to the North, initial greenschist facies stage to the South) or a different origin for the subducted materials (continental to the North and oceanic to the South). The differences in the anisotropic behavior of the subducted crust are reflected in the topography of the plate boundary, which becomes shallower from North to South, suggesting the existence of either a step in the slab topography or a more gentle ramp.

Accelerated subduction of the western Pacific Plate promotes the intracontinental uplift and magmatism in late Jurassic South China

Widespread fold and thrust belt and magmatism that occurred in the Mesozoic South China Block (SCB) have generated much debate on their dynamic mechanism related to the western Pacific Plate. Unlike a delayed Andean orogen, the SCB contains large areas of magmatic rocks coinciding with the absence of Upper Jurassic strata, which is not consistent with a flat-slab subduction geometry, and no consensus model has been suggested. To obtain more general insight into the effects of paleo-Pacific Plate subduction on the SCB, thermal-fluid-solid simulations of the subduction process are conducted. The models involved varied subduction rates and lithospheric viscosities to investigate their effect on the subduction style, mantle flow and overriding plate deformation. Low-rate subduction easily forms a subduction geometry with a higher angle, while high-rate subduction causes the slab to break off and drip downward. Fast subduction produces a quasi-horizontal mantle flow and a poloidal flow at the front of the mantle wedge, whereas slow subduction produces two poloidal flows developed on both sides of the subducting slab. The strain tensor of the overriding plate coincides well with the gradient of the horizontal mantle flow rate. A high subduction rate and low lithospheric viscosity promote the deformation of the overriding plate and are further attributed to mantle flow. Combining these results with global plate reconstructions, we find that accelerated subduction of the paleo-Pacific Plate resulted in the Late Jurassic uplift of the SCB. The low viscosity of the overriding plate caused by multiphase orogenic deformation promoted the thinning of the lithosphere due to induced mantle flow, which was associated with the intrusion of magmas and foundering of the subducted slab in a compressional regime. At the end of the Late Jurassic, subduction-induced mantle flow transformed from westward to eastward. Concurrently, the subducted slab began to steepen, and widespread Cretaceous magmatism occurred.

Numerical modeling on global-scale mantle water cycle and its impact on the sea-level change

Using numerical modeling of the global-scale mantle water cycle, I have clarified the mechanisms governing the deep Earth water content variations and evaluated their effects on sea-level changes. The following points are found: 1. The temporal variations of water content in the deep Earth are mainly determined by a balance between the ingassing flux from the ocean and the dehydration flux in the mantle wedge. 2. The primary control of ingassing flux comes from the water content of oceanic crust, rather than the style of plate subduction promoted by the free surface boundary condition. 3. The water content of oceanic crust ranges from 0.1 to 1 wt. percent to match the ingassing flux with the inferred one from the geological analysis.Based on the best-fit scenario for the ingassing flux from the ocean to the deep interior, I obtain an average rate of sea-level change of −1 to −1.5 m/Myrs over 1 billion year, which is nearly consistent with the long-term trend of sea-level change estimated from geological and paleontological data (∼−1 m/Myr). Note, however, that the detailed profile of sea-level change in the numerical modeling is difficult to explain the observation one. This underlines the importance of dehydration process in the mantle wedge in understanding the sea-level change resulting from the water budget in the Earth's deep interior over the Earth's history.

Linking Pacific Plate formation and Early Cretaceous metallogenic response on the circum-Pacific continental margins

Hydrothermal mineralization along the circum-Pacific continental margin is genetically linked to the motions of the Pacific Plate. Compiled geochronological results of ore deposits, coupled with the newly reconstructed geometry of the late Early Cretaceous subduction zones in East China, North America, and the Central Andes using GPlates and I2VIS software, were applied to investigate the relationships between mineralization and Pacific Plate formation. The ∼120-m.y.-old orogenic Au provinces in East China and North America are related to transpression caused by high-rate oblique subduction with intermediate−high dip angles of the Izanagi and Farallon plates, respectively. In contrast, the ∼105-m.y.-old porphyry-epithermal belt in Southeast China was produced by oblique subduction of the Izanagi Plate with low−intermediate subduction rates and intermediate−high dip angles. In the Central Andes, the oblique subduction of the Farallon Plate with low−intermediate rates and low dip angle accounted for iron oxide−copper−gold ore (IOCG) deposit mineralization in South Peru at ca. 110 Ma; whereas the high rate and low dip-angle subduction of the paleo-Phoenix Plate, which caused mild compression, was responsible for Fe, porphyry Cu, and IOCG mineralization in North Chile at ca. 110 Ma. The late Early Cretaceous metallogenic response in the circum-Pacific region coincides with superplume events that triggered the significant growth of the modern Pacific Plate. The forward simulations reveal that different subsequent styles of subduction and associated magmatism are likely responsible for the distinct mineralization types present in the region, including orogenic Au, porphyry-epithermal, Fe, and IOCG deposits. The precise dynamics of the subduction zone determined by this study led to the improved metallogenesis models in the Pacific margin.

Tectonic settings influence the geochemical and microbial diversity of Peru hot springs

Tectonic processes control hot spring temperature and geochemistry, yet how this in turn shapes microbial community composition is poorly understood. Here, we present geochemical and 16 S rRNA gene sequencing data from 14 hot springs from contrasting styles of subduction along a convergent margin in the Peruvian Andes. We find that tectonic influence on hot spring temperature and geochemistry shapes microbial community composition. Hot springs in the flat-slab and back-arc regions of the subduction system had similar pH but differed in geochemistry and microbiology, with significant relationships between microbial community composition, geochemistry, and geologic setting. Flat-slab hot springs were chemically heterogeneous, had modest surface temperatures (up to 45 °C), and were dominated by members of the metabolically diverse phylum Proteobacteria. Whereas, back-arc hot springs were geochemically more homogenous, exhibited high concentrations of dissolved metals and gases, had higher surface temperatures (up to 81 °C), and host thermophilic archaeal and bacterial lineages.

Joint Geodynamic‐Geophysical Inversion Suggests Passive Subduction and Accretion of the Ontong Java Plateau

AbstractIn this study, we for the first time applied a joint geodynamic‐geophysical inversion approach to oceanic plateau subduction models, and compared the subduction style and corresponding topography and Bouguer gravity of two representative subduction scenarios with passive or active collision. We showed that the case of passive collision of the Ontong Java Plateau (OJP) crust better explains the topography, gravity, and seismic data than the active collision scenario. This implies that the OJP did not control the regional dynamics during the collisional process. We conclude that previous studies may have overestimated the role of the OJP in triggering subduction initiation, subduction polarity reversal, and even Pacific Plate rotation.

Plateau archives of lithosphere dynamics, cryosphere and paleoclimate: The formation of Cretaceous desert basins in east Asia

Southeastern Eurasia is a global window to the Cretaceous paleoclimate and lithosphere coupling. China contains one of the most complete and complex sedimentary records of Mesozoic desert basins on planet Earth. In this study, we perform the spatio-temporal tracking of 96 Cretaceous palaeoclimate indicators during 79 Myr which reveal that the plateau paleoclimate archives from East Asia resulted from an Early to Mid-Cretaceous ocean–atmosphere coupling and a shift to a preponderant role of Late Cretaceous lithosphere dynamics and tectonic forcing on high-altitude depositional systems linked to the subduction margins of the Tethys and Paleo-Pacific realms beneath the Eurasian plate. The crustal response to tectonic processes linked with the spatio-temporal evolution of the Tethyan and Paleo-Pacific margins defined the configuration of major sedimentary basins on this region. The significant increase and decrease in the number of active sedimentary basins that occur during the Cretaceous, from 16 in the Early Cretaceous, to 28 in the Mid-Cretaceous, and a decreasing to 20 sedimentary basins in the Late Cretaceous, is a direct response of lithospheric dynamics associated with the two main subduction zones (Tethys and Pacific domains). A shift in subduction style from an Early Cretaceous Paleo-Pacific Plate slab roll back to a Late Cretaceous flat-slab mode might have triggered regional plateau uplift, blocked intraplate volcanism, thus enhancing the denudation and sediment availability, and created wind corridors that led to the construction and accumulation of extensive Late Cretaceous aeolian sandy deserts (ergs) that covered Mid-Cretaceous plateau salars. At the same time, plateau uplift associated with crustal thickening following terrane assembly in the Tethyan margin triggered altitudinal cryospheric processes in sandy desert systems. Evidence of an active Cretaceous cryosphere in China include Valanginian-Hauterivian glacial debris flows, Early Aptian geochemical signature of melt waters from extensive ice sheets, and Cenomanian–Turonian ice-rafted debris (IRD). These cryospheric indicators suggest an already uplifted plateau in southeastern Eurasia during the Cretaceous, and the marked correlation between cold plateau paleoclimate archives and marine records suggests a strong ocean-atmosphere coupling during Early and Mid-Cretaceous cold snaps. We thus conclude that lithospheric tectonics during Cretaceous played a fundamental role in triggering high-altitude basin desertification and spatio-temporal plateau paleohydrology variability in the Cretaceous of south-eastern Eurasia.

Styles of Trench‐Parallel Mid‐Ocean Ridge Subduction Affect Cenozoic Geological Evolution in Circum‐Pacific Continental Margins

AbstractDistinct signatures are present in the circum‐Pacific continental margins (e.g., kinematics, magmatism, and basin evolution), possibly influenced by the input of mid‐ocean ridge. It remains enigmatic why the circum‐Pacific continental margins that have experienced trench‐parallel mid‐ocean ridge subduction show diverse geological evolution. Here we present geodynamic modeling results investigating trench‐parallel mid‐ocean ridge subduction and demonstrate two distinct types of model evolution. Type‐Ⅰ model includes a two‐stage steep subduction and is featured by slab detachment preceding the arrival of ridge at the trench. Type‐Ⅱ model is marked by a continuous flat subduction of mid‐ocean ridge with the opening of a slab window beneath intracontinental lithosphere. These two subduction styles produce diverse tectono‐magmatic responses. Our results could explain the magmatic gap and forearc uplift during the Izanagi‐Pacific ridge subduction and the intraplate magmatic flare‐ups and tectonic uplift during the Nazca‐Antarctic ridge subduction, respectively.

Seismic evidence of two cryptic sutures in Northwestern Australia: Implications for the style of subduction during the Paleoproterozoic assembly of Columbia

Plate tectonics, including rifting, subduction, and collision processes, was likely to have been different in the past due to the secular cooling of the Earth. The northeastern part of the West Australian Craton (WAC) has a complex Archean and Paleoproterozoic tectonic history; therefore, it provides an opportunity to study how subduction and collision processes evolved during the emergence of plate tectonics, particularly regarding the assembly of Earth's first supercontinent, Columbia. Because the northeastern boundary of the WAC and the southwestern boundary of the North Australian Craton (NAC) are covered by the Phanerozoic Canning Basin, the regional tectonic evolution has remained enigmatic, including how many tectonic elements were assembled and what may have driven rifting and subsequent collision events. Here, we use new passive-source seismic modeling to identify a seismically distinct segment of the lithosphere, the Percival Lakes Province, which lies east of the Pilbara Craton and is separated by two previously unknown southeast-trending lithosphere scale Paleoproterozoic sutures. We interpret that the northeastern suture, separates the Percival Lakes Province from the NAC, records the amalgamation of the WAC with the NAC. The southwestern suture separates the PLP from the reworked northeastern margin of the Pilbara Craton, including the East Pilbara Terrane and the Rudall Province. A significant upper mantle dipping structure was identified in the southwestern suture, and we interpret it to be a relic of subduction that records a previously unknown Paleoproterozoic collision that pre-dated the amalgamation of the WAC and NAC. By comparing our findings with previously documented dipping features, we show that the Paleoproterozoic collisions are seismically distinguishable from their Phanerozoic counterparts.

Numerical modeling of North Sulawesi subduction zone: Implications for the east–west differential evolution

Geological observations have revealed a rapid evolution and obvious east–west differential evolution in the North Sulawesi subduction zone. The Celebes Sea plate has inserted itself under the north arm of Sulawesi Islands, which has simultaneously rotated clockwise. The rotation of the north arm of the Sulawesi Islands might be critical in facilitating tectonic processes, such as the slab subduction and rollback of the Celebes Sea plate. For the east–west differences along the North Sulawesi subduction zone, a numerical model with the convergence rate of the plates as the basic variable is established to quantitatively describe the evolution process of the North Sulawesi subduction zone. Our results reproduced the east–west differences of the subducting Celebes Sea plate, showing a shallow–deep–shallow subduction style. We propose that the variable velocity ratio of the overriding plate to the subducting slab could be the principal reason for the differential subduction along the strike of the North Sulawesi subduction zone. We conclude that the residual slab and the rate of the eastern continental plate limit the downward movement of the subducted slabs of the eastern Sulawesi. Furthermore, the reason for the shorter subducted slab at the extreme western Sulawesi is the fact that subduction occurs outside the rotation radius. Moreover, the widespread extension at the western Sulawesi has a limited correlation with the clockwise rotation.

Genetic Environments of the Eunjeok Au–Ag Deposit in the Yeongam District: Implications for Cretaceous Epithermal Au–Ag Mineralization in South Korea

Eunjeok Au–Ag deposits are situated in the Yeongam district, Cheollanamdo-province, South Korea. They are genetically related to the Bulgugsa magmatic event (ca. 110–60 Ma), caused by the transition in the subduction direction and style of the Izanagi Plate. Three gold- and silver-bearing hydrothermal veins filled the fractures of the Cretaceous rhyolitic tuff. The major ore minerals were arsenopyrite (31.47–32.20 at.% As), pyrite, chalcopyrite, sphalerite (8.58–10.71 FeS mole%) and galena with minor amounts of electrum (62.77–78.15 at.% Au), native silver, and argentite. Sericitization was dominant in the alteration zone. The various textures of quartz veins (i.e., breccia, crustiform, comb, and vuggy) may indicate the formation of an epithermal environment. The auriferous fluids with the H2O–NaCl system have homogenization temperatures and salinities of 204 °C to 314 °C, less than 10 wt.% equiv. NaCl, and experienced mixing (dilution and cooling) events during mineralization. Considering the characteristics of the geologic setting, major fault system, and host rock, the Eunjeok Au–Ag deposit within the Yeongam district tends to share the general geologic characteristics of Haenam–Jindo epithermal mineralization episodes. However, the age of gold–silver mineralization (86.0 Ma) is older than that of Haenam–Jindo epithermal mineralization episodes (<70.3 Ma), implying some differences exist in the genetic sequence of extensional characteristics caused by transcurrent Gwangju–Yeongdong faults.

Intermittent Post‐Paleocene Continental Collision in South Asia

AbstractThree different conceptual models have been proposed for the Cenozoic subduction style in South Asia, including Greater India, Intra‐oceanic Arc, and Continental Terrane (or Greater Indian basin). Since these models imply distinctive origins for the Tethyan–Greater Himalayan (TGH) sequences, for example, as a relic of the subducted Greater India or Gondwana–affiliated continental terrane, quantitively reproducing the relic TGH crustal mass with numerical models could help further constrain the debated Cenozoic subduction history between India and Eurasia. Based on the modeling results, we show that the subducted plate since the Paleocene should consist of a significant oceanic portion that is, ∼1,000 km long for the Intra‐oceanic Arc model and up to 2,000 km long for the Terrane model. Our results do not support the existence of a continuous >3,000 km long continental Greater India before the early Eocene collision in South Asia.