Subduction Style Research Articles

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.

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.

Oceanic Plateau and Subduction Zone Jump: Two‐Dimensional Thermo‐Mechanical Modeling

AbstractOceanic plateaus are characterized by overthickened crust and rheologically strong and chemically depleted mantle root with respect to normal oceanic lithosphere. The buoyant nature makes them have a great potential to resist subduction and initiate a new subduction zone behind the plateaus, yet the geodynamic process remains enigmatic. Here, we use 2D thermo‐mechanical modeling method to investigate the role of an oceanic plateau in subduction zone reorganization. We systematically examine the rheological strength and compositional density of oceanic plateau mantle and the size of oceanic plateau. Three distinct subduction styles are observed: (a) steep subduction without subduction zone reorganization, favored by a small‐sized (≤150 km) or lowly depleted oceanic plateau; (b) breakoff of subducting slab followed by flat‐slab subduction, which is conditional on a moderate‐sized (≥300 km) oceanic plateau with high degree depletion but relatively low strength; and (c) initiation of a new subduction zone behind the plateau (subduction zone jump), facilitated by a moderate‐sized oceanic plateau with a strongly depleted mantle root. The geodynamic results demonstrate that the overriding lithosphere undergoes a tectonic force switch from compression to tension in ∼5 Myr as the subduction zone is reorganized. This can be reflected in the crustal deformation of the overriding plate and used as an indicator to identify past subduction zone jumps. The spatial relationship predicted by the subduction zone jump models explains the Jurassic closure event of the accretion on southern Qiangtang at ∼175 Ma and the intraoceanic infant arc magmatism in the Meso‐Tethys at ∼165 Ma.

Balanced Cross‐Sections and Numerical Modeling of the Lithospheric‐Scale Evolution of the Hindu Kush and Pamir

AbstractThe role of deformation and removal of continental lithospheric mantle (CLM) during continent‐continent collision, while recognized as fundamental to the evolution of orogenic systems, often remains speculative as direct evidence for these processes is typically absent. Nevertheless, these poorly understood deep lithospheric processes are expected to produce first‐order controls on the style of crustal shortening, magmatism, and metamorphism in large‐scale orogens. The Hindu Kush and Pamir offer an excellent opportunity to constrain these deep processes, as high‐quality seismic images reveal oppositely dipping slabs under the Hindu Kush and Pamir, inferred to be subducting Indian and Asian lithospheres, respectively. Using two‐dimensional mantle‐scale geodynamical models, we show that these along‐strike differences can be explained by the nature of the indenting Indian plate; thinner fringe Indian crust beneath the Hindu Kush, and full‐thickness Indian crust beneath the Pamir. In the Hindu Kush models, steep subduction of India follows early delamination of the Asian lithospheric mantle, whereas in the Pamir models, Asian delamination is followed by flat‐slab underthrusting of Indian CLM and forced subduction of northern Asian terranes. These styles of mantle delamination and subduction have key implications for the ensuing styles of magmatism, metamorphism, and crustal‐scale evolution of these systems. In the Hindu Kush, subduction of fringe India leads to burial of Indian crust under Asia, limiting both Himalayan growth and deformation of the Tajik Basin. In contrast, in the Pamir, forced subduction of the Alai terrane leads to Alai crust underthrusting the Pamir, forming the Pamir orocline.

Earthquake distribution and lithospheric rheology beneath the Northwestern Andes, Colombia

The rheological behavior of the lithosphere beneath the Northwestern (NW) Andes (Colombian) remains uncertain under a complex geological setting. Here, we analyze the earthquake distribution and the stress drop on the western and the eastern of the transition area between the Upper Magdalena Valley (UMV) and the Middle Magdalena Valley (MMV), a zone related to a plate tectonic transition from steeper to flatter subduction style slabs (Nazca and Caribbean plates). Our results consistently explain the tectonic and geodynamic context of the study area supported in two rheological profiles, which reveal the brittle/ductile transition of the lithospheric crust and mantle and allow to constrain the lithospheric rheological thickness. The subduction phenomenon of the Nazca Plate under the South American Plate is observed in the western zone. In the eastern region, we find a remarkable aseismic zone under the Eastern Cordillera below 20-km-deep, which made more difficultan accurate estimation of the crust/mantle boundary. This seismic particularity leads us to support the hypothesis of a delamination process due to the tectonic overload generated by the Eastern Cordillera, and that promotes thermal anomalies. Our results also suggest that the earthquake dataset correlated with rheological estimations may offer a consistent interpretation of the mechanical behavior of the lithosphere.

Mantle cooling and cratonization of Archean lithosphere by continuous plate subduction: Constraints from TTGs, sanukitoids, and high-K granites, eastern North China Craton

The evolution of plate subduction styles on the Earth is coupled with its cooling history. Continuous subduction was initiated during the Neoarchean, yet its influence on mantle cooling and the cratonization of Archean lithosphere remains unclear. Successive subduction-related magmatism occurred during the late Neoarchean evolution of the Western Shandong Province granite-greenstone belt (WSP) of the North China Craton (NCC) and preserved crucial clues relating to processes and roles of regional continuous subduction. In this study, petrology, whole rock geochemistry, and zircon U–Pb–Hf isotopes are reported for late Neoarchean TTG gneisses, quartz monzodiorites, and monzogranites from the WSP. Zircon U–Pb dating reveals that the TTG gneisses, quartz monzodiorites, and monzogranites were formed at ~2561–2517 Ma, 2540–2533 Ma, and 2541–2513 Ma, respectively. The TTG gneisses with high (La/Yb)N ratios, low K2O/Na2O ratios, and low MgO contents were mainly derived from partial melting of oceanic slab material. In contrast, the quartz monzodiorites show high K2O and MgO content and high (La/Yb)N ratios, consistent with the characteristics of Archean sanukitoids, suggesting that they were the products of partial melting of mantle source metasomatized by slab-derived melts. The monzogranites with systematically high K2O, SiO2, and low MgO contents were formed by crustal reworking of regional TTGs and sedimentary rocks. Rock assemblages, geochemical features, regional field observation, and structure geology indicate prolonged (~2.56–2.49 Ga) continuous magmatism in the WSP across three northwest-trending belts (Belts A, B, and C) related to continuous subduction and slab rollback. These processes led to a tectonic transition from early compression to late extension. The final stabilization of the continental crust in the late Neoarchean in the eastern NCC was a consequence of successive mantle cooling due to regional continuous subduction recorded by subduction-related rock assemblages (i.e., adakites, sanukitoids, high-Mg basalts-andesites, and Nb-enriched basalts), the paired type of metamorphism, and the rise of alkali basalts. Secular geochemical changes of alkaline basalts outline the role of continuous subduction in dramatic mantle cooling and the final stabilization and cratonization of the eastern NCC.