Slab Subduction Research Articles

Zircon petrochronological evidences for diachronous tectonic switch and orogenic keel collapse of the Tongbai-Dabie-Sulu orogens

Migmatites can be produced in different stages of orogeny and thus have a great significance to constrain the evolution of an orogenic belt. Migmatites widely occur in the Tongbai orogen, but their formation ages and geological significance have not been well constrained. In this study, a combination of cathodoluminescence (CL) imaging, U-Pb-Hf-O isotopes and trace elements of zircons were carried for migmatites in the Tongbai orogen. Most inherited cores in a melanosome show variable UPb ages and Hf isotope compositions with high δ18O values, suggesting they might be detrital zircons from its protolith. Whereas other relict cores in this sample have nebulously zoning or no zoning, flat HREE patterns with insignificant Eu anomalies and relatively low formation temperatures, indicating their formation under eclogite-facies conditions. They yield UPb ages of 228 ± 3 Ma, providing direct evidence for the involvement of Triassic subducted continental materials. The overgrowth rims and some new grown grains in the migmatites show euhedral shapes, low CL luminescence, high Th/U ratios, clear negative Eu anomalies and steep HREE patterns, suggesting they are anatectic zircons formed during partial melting. They yield UPb ages from 139 ± 2 to 132 ± 1 Ma. Combined with previous studies, we suggest that migmatization mainly occurred at 139–130 Ma in the Tongbai orogen. The anatectic zircons in the leucosomes show much lower δ18O values (5.79 ± 0.23‰ vs. 8.94 ± 0.16‰) and higher 176Hf/177Hf ratios (0.282302–0.282648 vs. 0.281728–0.282504) than the zircons in adjacent melanosomes, suggesting that the leucosomes might be externally sourced. According to a compilation of the geochronological information of migmatite and magmatism, the Tongbai-Dabie-Sulu orogens have a spatial and temporal variation of partial melting and magma activities in the post-collisional stage. These might result from the diachronical tectonic switch and collapse of the orogenic lithospheric keels caused by the subduction and rollback of the Paleo-Pacific slab.

Horizontal ground-motion model for subduction slab earthquakes using offshore ground motions in the Japan Trench area

The S-net, a large-scale network of permanent ocean-bottom seismographs in the Japan Trench, consisting of 150 stations, has recorded a large number of offshore ground motions that can be used to establish the empirical attenuation relationship of offshore ground motions in the subduction zone. Due to the different attenuation characteristics among different earthquake types in subduction zones, the earthquakes are classified into four tectonic types of subduction slab, interface, shallow crustal, and upper-mantle earthquakes according to the classification scheme of Zhao et al. (2015). Predictive models and uncertainties in offshore ground motions were investigated for different earthquake types. This study establishes a horizontal offshore subduction slab earthquake ground-motion model (GMM) and compares the offshore and onshore slab earthquake GMMs. As the site conditions at ocean-bottom stations are different from those at land stations, the effects of water depth and sediment thickness are taken into account in the offshore slab earthquake GMM. Due to differences in the burial methods of ocean-bottom stations, stations were divided into buried and unburied to investigate the effects of the burial methods. Therefore, regression analysis was used to propose an offshore slab earthquake GMM considering the magnitude, focal depth, distance, water depth, sediment thickness, and burial method. By separating the within-event residuals, the single-station standard deviation is presented. Compared to the onshore GMM, the predicted spectra of the offshore GMM are significantly larger at long periods. For the attenuation rate, the offshore attenuation rate is lower than that of the onshore attenuation rate for short periods, but is basically consistent with the onshore attenuation rate for long periods. The proposed GMM can be used to predict the offshore ground motions for slab earthquakes with rupture distances less than 300 km, focal depths less than 110 km, and moment magnitude between 4 and 7.4.

Evidence for Ultra-Low Velocity Zone Genesis in Downwelling Subducted Slabs at the Core–Mantle Boundary

Abstract We investigate broadband SPdKS waveforms from earthquakes occurring beneath Myanmar. These paths sample the core–mantle boundary beneath northwestern China. Waveform modeling shows that two ∼250 × 250 km wide ultra-low velocity zones (ULVZs) with a thickness of roughly 10 km exist in the region. The ULVZ models fitting these data have large S-wave velocity drops of 55% but relatively small 14% P-wave velocity reductions. This is almost a 4:1 S- to P-wave velocity ratio and is suggestive of a partial melt origin. These ULVZs exist in a region of the Circum-Pacific with a long history of subduction and far from large low-velocity province (LLVP) boundaries where ULVZs are more commonly observed. It is possible that these ULVZs are generated by partial melting of mid-ocean ridge basalt.

Thermal equation of state of Li-rich schorl up to 15.5 GPa and 673 K: Implications for lithium and boron transport in slab subduction

Abstract The thermal equation of state (EoS) of natural schorl has been established at high temperatures up to 673 K and high pressures up to 15.5 GPa using in-situ synchrotron X-ray diffraction combined with a diamond anvil cell. The pressure-volume (P-V) data were fitted to the third-order Birch-Murnaghan EoS with V0 = 1581.45 ± 0.25 Å3, K0 = 111.6 ± 0.9 GPa and K′0 = 4.4 ± 0.2; additionally, when K′0 was fixed at a value of 4, V0 = 1581.04 ± 0.20 Å3 and K0 = 113.6 ± 0.3 GPa. The V0 (1581.45 ± 0.25 Å3) obtained by the third-order Birch-Murnaghan EoS agreed well with the measured V0 (1581.45 ± 0.05 Å3) under ambient conditions; this result confirmed the high accuracy of the experimental data in this study. Furthermore, the axial compression data of the schorl at room temperature were also fitted to a “linearized” third-order Birch-Murnaghan EoS, and the obtained axial moduli for the a- and c-axes were Ka = 621 ± 9 GPa and Kc = 174 ± 2 GPa, respectively. Consequently, the axial compressibilities were βa = 1.61×10-3 GPa-1 and βc = 5.75×10-3 GPa-1 with an anisotropic ratio of βa:βc = 0.28:1.00, indicating axial compression anisotropy. In addition, the compositional effect on the axial compressibilities of tourmalines was also discussed. Fitting our pressure-volume-temperature (P-V-T) data to the high-temperature third-order Birch-Murnaghan EoS provided the following thermal EoS parameters: V0 = 1581.2 ± 0.2 Å3, K0 = 110.5 ± 0.6 GPa, K′0 = 4.6 ± 0.2, (∂KT/∂T)P = -0.012 ± 0.003 GPa K-1 and αV0 = (2.4 ± 0.2) × 10-5 K-1. The obtained thermal EoS parameters in this study were also compared with those of previous studies on other tourmalines. The potential factors influencing the thermal EoS parameters of tourmalines were further discussed.

The Zaozigou orogenic gold-antimony deposit, West Qinling Orogen, China: Structural controls on multiple mineralization events

Orogenic gold deposits are generally thought to represent one perhaps protracted event. However, recent research on orogenic gold deposits increasingly offers evidence for some deposits forming through multiple and clearly discreet hydrothermal episodes. The giant Zaozigou orogenic Au-Sb deposit in the Triassic to Cretaceous West Qinling Orogen, central China, includes both steeply dipping and gently dipping orebodies. The two distinct mineralization styles provide a valuable setting for investigating a multiple mineralization model by integrating structural analysis within a robust geochronological framework. Through fieldwork and geochronology, we define a progression of major tectonic events in the area of the Zaozigou deposit. The deposit is hosted within a well-bedded sequence of Early Triassic metasedimentary rocks of the South Qinling Terrane. Pre-mineralization E-W shortening (D1) during subduction of the Mianlue oceanic slab include folding with resulting axial planes striking N-S, emplacement of Triassic ENE-striking and WNW-striking dacite dikes accompanied by Middle Triassic greenschist facies metamorphism. Late Triassic gold-stibnite quartz vein and disseminated mineralization formed along ENE-striking and steeply dipping D2 brittle to ductile sinistral faults. Their orientations suggest a link to the regional NNE-SSW maximum principal stress coinciding with transpression caused by the Late Triassic collision between the South China Block and South Qinling Terrane. Overprinting Early Cretaceous quartz-stibnite veins developed along gently dipping (20° to 40°) brittle D3 normal fault zones, which exhibit a NE-SW minimum principal stress. This younger deformation event is interpreted to be related to the Early Cretaceous tectonic transition from shortening to extension of the West Qinling Orogen. Therefore, the Zaozigou deposit reveals a model of multiple orogenic gold mineralizing events, with migration of hydrothermal fluids during discrete deformation episodes and the resulting formation of a single composite deposit formed along overprinting structures at separate times of orogenesis.

Exhumation and preservation of the Bangbule Pb-Zn-Cu deposit in the western Gangdese metallogenic Belt, Tibet: Constraints from fission track thermochronology and fault gouge dating

The Bangbule Pb-Zn-Cu deposit is a large-scale skarn deposit that was recently discovered in the western Gangdese Metallogenic Belt. Previous studies have focused mainly on ore genesis and ore-forming fluid, while exhumation and preservation after mineralization have not been the focus. We report fission track ages (ZFT: 43 ± 3 to 52 ± 3 Ma; AFT: 24 ± 2 to 32 ± 3 Ma, MTL: 12.2 ± 1.9 to 13.0 ± 1.6 μm) and the timing of fault activation (44.3 ± 0.9 Ma), revealing the processes of cooling/exhumation and we evaluate the preservation of the deposit. The Bangbule deposit experienced two rapid exhumation events since the Eocene (52–24 Ma and 12 Ma to present) and a slower stage during the late Oligocene to middle Miocene (24–12 Ma). The exhumation thickness (1.6–3.1 km) after mineralization was similar to the metallogenic depth (peak: 1.5–2.5 km). The volcanic rocks of the Dianzhong Formation partially protect the deposit, making mineral exploration of the Bangbule deposit feasible. Since the Oligocene, due to climate change and deep subducted slab tearing, the exhumation thickness in the western Gangdese (∼0.6–3.1 km) has generally been lower than that in the middle-eastern part (∼3–8 km). The lower degree of exhumation may be one of the reasons why less large-scale deposits have been discovered in the western Gangdese.

Fluid-metasomatized rocks with extremely low δ26Mg values in subducted oceanic lithosphere: Implications for mantle Mg isotope heterogeneity and the origin of low-δ26Mg magmas

The subduction of oceanic slabs with diverse lithologies greatly contributes to the geochemical heterogeneity of the mantle. Oceanic metasomatism leads to the formation of rocks with heterogeneous geochemical characteristics, which can be further modified by fluid–rock interactions during subduction zone processing. Investigations of exhumed high-pressure oceanic rocks not only shed light on the compositional variations that occur during oceanic metasomatism and subduction processes but also help in understanding mantle heterogeneity and magma genesis. Here, we present a comprehensive Mg–O–Sr–Nd isotope study on eclogite-facies Mg- and Ca-metasomatized rocks (Ti-clinohumite dikelets and metarodingites, respectively) and their host meta-serpentinites from the Voltri Massif (Ligurian Western Alps). Compared to their Fe–Ti gabbroic protoliths, the Ti-clinohumite (Ticl) dikelet and metarodingite are characterized by high MgO contents (∼32–38 wt%) and CaO contents (∼17–25 wt%), respectively. These Mg- and Ca-enriched rocks both display low δ26Mg values of −0.54 ‰ to −0.43 ‰ and −1.49 ‰ to −0.82 ‰, respectively, while the serpentinite hosts have mantle-like δ26Mg values of −0.26 ‰ to −0.20 ‰. The Ticl-chlorite rind rimming the metarodingite has a low δ26Mg value of −1.00 ‰. No carbonate minerals were found for metasomatized rocks that have low δ18O values ranging from +2.0 ‰ to +3.4 ‰, demonstrating that their low δ26Mg signature cannot be attributed to metasomatism by fluids involved with carbonate components. The increase in MgO content and decrease in δ26Mg in the Ticl dikelets are ascribed to fluid-induced mass transfer through low-δ26Mg fluids buffered by serpentinization of the host ultramafic rocks. The metarodingite has comparable MgO but extremely low δ26Mg values compared to those of the gabbroic protolith. By integrating reported Mg isotope data for low-grade rodingites, we propose that the low δ26Mg values of metarodingites are ascribed to isotopic exchange with serpentinizing fluids during progressive rodingitization reactions at the seafloor. However, the extremely low δ26Mg value of the Voltri metarodingite (as low as −1.49 ‰) may be attributed to further metasomatic modification in the subduction zone. High pressure metamorphic recrystallization of these metasomatized rocks play a limited role in their Mg isotope compositions. Therefore, our results show that the oceanic metasomatized metarodingite and Ticl-bearing dikes hosted in serpentinite display systematically low-δ26Mg features. Such isotopically light Mg compositions can be preserved during subduction to eclogite-facies conditions. Hence, subduction of these metasomatized rocks may lead to the transfer of a low-δ26Mg component into the deep mantle, which was overlooked thus far and may further contribute to the origin of low-δ26Mg magmas.

Early Cretaceous continental arc magmatism in the Wakhan Corridor, South Pamir: Mantle evolution and geodynamic processes during flat subduction of the Neo-Tethyan oceanic slab

The petrogenesis of continental arc magmas provides critical insights into thermal evolution and geodynamics of the continental lithosphere, crust-mantle interaction, and deep dynamic processes. In this study, we report new zircon U-Pb ages along with isotopic and elemental whole-rock geochemistry, mineral chemistry, and Hf-O isotope data for the Kalaqigu diorites and monzogranites of the Chinese Wakhan Corridor, South Pamir. Zircon U-Pb dating indicates that the Kalaqigu pluton was emplaced in the Early Cretaceous (ca. 108−106 Ma). The diorites are geochemically characterized by low SiO2 (51.9−54.5 wt%) and CaO (7.7−9.4 wt%) contents, but high MgO (5.3−8.3 wt%), Al2O3 (12.8−16.8 wt%), and TiO2 (0.6−1.1 wt%) contents as well as high Mg# (56−65) values. Thus, they are similar to high-Mg diorites: enriched in large ion lithophile elements (e.g., K, Sr, and Ba) and light rare earth elements, while depleted in high field strength elements (i.e., Nb, Ta, Zr, and Hf). Combined with negative εNd(t) (−6.9 to −14.0) and εHf(t) (−9.9 to −12.2), and high (87Sr/86Sr)i (0.7075−0.7086) ratios, these observations indicate that they originated from an enriched lithospheric mantle source. High δ18Ozrn (7.49‰−9.01‰) values, in conjunction with relatively high 207Pb/206Pb and 208Pb/206Pb ratios, suggest that the source was modified by subducted sediment-derived melts. Variable Cr contents (54−117 ppm) are likely controlled by minor fractionation of olivine and orthopyroxene. The monzogranites show high SiO2 contents (69.2−72.0 wt%), and low Rb/Sr (0.4−0.6), (K2O + Na2O)/CaO (2.6−4.8), and FeOT/MgO ratios (2.6−3.2). They contain diagnostic cordierite and show strongly peraluminous characteristics (A/CNK > 1.1) with high δ18Ozrn (7.82‰−8.85‰) values that are compatible with those of typical S-type granites. Their abundant inherited zircons, with age populations similar to those of detrital zircons from regional early Paleozoic metasedimentary rocks, indicate that they were derived from partial melting of ancient metasedimentary rocks. Phase equilibrium modeling is consistent with biotite-dehydration melting of metagreywacke, probably at ∼750 °C and ∼6.0 kbar, as indicated by the biotite chemistry. Based on regional geochronology, a south-to-north magmatic migration suggests that northward flat-slab subduction of the Neo-Tethyan oceanic slab played an important role in the generation of these widespread Early Cretaceous continental arc magmatic rocks. However, the granitoids were generated earlier than the mantle-derived mafic rocks, which suggests that crustal melting occurred during the early stage of subduction. The continuous flat-subduction resulted in partial melting of subducted sediments, which metasomatized the mantle wedge. Contemporaneous regional compression primarily occurred far north of the subduction zone (i.e., North and Central Pamir), inducing deformation as well as crustal shortening. With the flare-up of continental arc magmatism in South Pamir, crustal shortening moved southward. These processes, combined with the addition of voluminous, mantle-derived magmas, played an important role in crustal thickening in Pamir during the Early Cretaceous.

Lithospheric mantle as a metal storage reservoir for orogenic gold deposits in active continental margins: Evidence from Hg isotopes

Abstract The subcontinental lithospheric mantle (SCLM) is now suggested by many workers to play a role in the formation of orogenic gold deposits in active continental margins, given that the gold and volatiles may be introduced into the SCLM during plate subduction. The giant Cretaceous Jiaodong gold province within the North China block occurs in a convergent margin setting where there is no terrane accretion, thus eliminating from consideration the commonly accepted metal and fluid source reservoir for orogenic gold—the metamorphosed upper crust of an active orogenic belt. Thus, the auriferous fluids were released below the high-grade Archean cratonic rocks from either enriched SCLM or directly from the subducted oceanic slab. Mercury (Hg) isotopes allow fingerprinting of isotopically distinct reservoirs, discriminating between the two possibilities. We present Hg isotope data for a set of pyrite grains, native gold, and lamprophyres from six main gold districts in the Jiaodong gold province. These samples dominantly yielded near-zero Δ199Hg (0‰ ± 0.1‰), within uncertainty of the mantle Δ199Hg value. The results provide novel evidence that the SCLM predominantly controlled the Hg budget of the Jiaodong gold deposits, implying that the SCLM acted as a buffer for the crustal mass-independent fractionation Hg signatures and a storage reservoir for gold and volatiles. Consequently, during slab subduction, gold and volatiles were recycled and stored within the SCLM and were later released from the mantle during uplift and relaxation of compressional tectonics.

Deciphering differential exhumation in the Gangdese orogen in southern Tibet using exposed porphyry alteration systems and geomorphic analysis

Topography and exhumation vary strongly in the Tibetan Plateau, reflecting different geodynamics and surface processes. The Gangdese orogen in southern Tibet is characterized by west-to-east variations in topography and occurrences of the Miocene porphyry copper deposits. We integrated remote sensing based hydrothermal alteration mapping and geomorphological analysis to detect patterns of differential exhumation since the Miocene. High abundances of propylitic and potassic alteration of the eastern Miocene granitoids are identified, indicating >2 km of exhumation since the Miocene. The western Gangdese has intensive phyllic and argillic alteration representing ∼1 km of exhumation. Geomorphic analyses also show E-W changes that most likely resulted from differential rock uplift histories. Channel steepness index values indicate the eastern Gangdese has channel gradients between 211 m1.02 and 457 m1.02, compared to values between 30 m1.02 and 160 m1.02 in the west. Linear inversion of tributaries of the Yarlung river shows a pulse of rock uplift between 25 Ma and 7 Ma with a maximum rate of 110 m/ m.y. in the east. In contrast, the western Gangdese exhibits slower and persistent rock uplift rates between 20 m/m.y. and 40 m/m.y. since 20 Ma. The integrated results record differential exhumation and rock uplift since the Miocene. We invoke the geometric variations of the subducted Indian lithosphere to interpret the E-W differential exhumation since the Miocene. Steeper subduction of the Indian lithosphere beneath the eastern Gangdese induced mantle upwelling, triggering rapid rock uplift and deep exhumation. Conversely, the shallower subduction of the Indian slab in the west generated slow and prolonged rock uplift with less exhumation. Our study refines the exhumation process of the Gangdese since the Miocene.

Geology, fluid inclusions, mineral and (S-O) isotope chemistry of the Badran orogenic Au deposit, Yana-Kolyma belt, eastern Siberia: implications for ore genesis

The Badran orogenic gold deposit is located in the Yana-Kolyma belt, Eastern Siberia; it has proven reserves of ∼9.3 t of gold and an average grade of 7.8 ppm. The total gold production at the Badran deposit since 1984 amounts to ∼34 t. Despite many years of study, the origin of the gold deposits of the Yana-Kolyma metallogenic belt, one of the world’s largest belts, and the Badran deposit is controversial. Synthesis of regional geology and geology of the Badran deposit, fluid inclusion analysis, mineral and (S-O) isotope chemistry defines the genetic model, origin of fluids, and source of metals in the evolution of the ore-forming system, equivalent to other orogenic gold deposits on the margin of the Siberian craton. The deposit is localized in the Upper Triassic clastic rocks and is controlled by the NW-trending thrust. Polyphase mineralization occurs as disseminated arsenian pyrite and arsenopyrite ores with invisible gold, quartz veins with native gold and Fe, Pb, Zn, Cu sulfides and sulfosalts of orogenic type, and locally post-ore Ag, Sb-bearing minerals and Hg epithermal features. The quartz veins with native gold were formed from low-medium saline (1.5–10 wt% NaCl eq.) aqueous-carbonic fluids boiling at temperatures of 290°C to 210 °C and pressures of 300–250 to 125 bar. The δ34S values of pyrite and arsenopyrite vary from −1.1‰ to +2.4‰, with an average of +0.4‰; the δ18О of quartz from +15.1‰ to +17.5‰ at constant δ18ОH2O about +7.5‰ (±1.0‰). High contents of As (up to 2.4 wt%) and Co/Ni ∼ from 0.3 to 9.9 in pyrite of proximal alteration are typical for hydrothermal systems. The results obtained confirm that the ore-forming fluids did not have a single origin, but were formed from a mixture of subcontinental lithospheric mantle and metamorphic sources. The subcrustal lithospheric mantle was fertilized in the time preceding mineralization (Late Jurassic) and was derived directly from the down-going subduction slab and overlying sediment wedge at the closure of the Oymyakon Ocean.

Transition from oceanic subduction to continental collision: Insights from volcanogenic-sedimentary rocks of the Tannuola terrane (northern Central Asian Orogenic Belt)

Magmatism associated with oceanic subduction plays a dominant role in crustal growth during the Earth’s evolution. The Tannuola terrane, situated in the northern Central Asian Orogenic Belt (CAOB), is a key area to understanding oceanic subduction and initial collision processes in the northern CAOB. In order to investigate the switch from subduction to collision settings, detailed field mapping, zircon SHRIMP U-Pb geochronological and whole-rock geochemical studies of volcanogenic-sedimentary rocks from the Tannuola terrane were carried out. Zircon U-Pb ages indicate multi-stage volcanism lasted at least 30 Ma from ∼540 to ∼510 Ma, that can be divided into three events: the late Ediacaran (before ∼540 Ma), the early Cambrian (∼520 Ma) and the middle Cambrian (∼510 Ma). These ages are interpreted to the initial, main and final stages of oceanic subduction during the late Proterozoic – early Paleozoic. During the late Ediacaran, tholeiitic basalts with high εNd(t) values (from +7.4 to +8.5) were emplaced. Likely forming by the 10 %–30 % partial melting of spinel – garnet mantle source during slab subduction. During the early Cambrian, transitional from tholeiitic to calc-alkaline basaltic rocks with εNd(t) value (+5.6) and coeval intermediate–felsic volcanic rocks with similar εNd(t) values (+5.9 and +6.5) formed. The early Cambrian basaltic rocks are interpreted to be derived by 10 %–30 % partial melting of a depleted mantle source metasomatized by slab-derived fluids released from the subducting oceanic slab. The middle Cambrian calc-alkaline basaltic rocks with εNd(t) value of +6.2 might be emplaced as a result of low (5 %–10 %) degree partial melting of a metasomatized mantle followed by fractional crystallization of clinopyroxene and plagioclase. Associated intermediate-felsic volcanic rocks with εNd(t) values from +6.0 to +6.8 were formed through fractionation of the juvenile Neoproterozoic sources. The middle Cambrian volcanism is interpreted to be triggered by the slab break-off during the transition to a collisional setting.

Geodynamic models of Indian continental flat slab subduction with implications for the topography of the Himalaya-Tibet region

The slab structure and high elevation of the Himalaya-Tibet region and their underlying mechanisms have been widely discussed. Many studies interpret a flat slab segment of Indian continental lithosphere located below the overriding plate, but interpretations of the northward extent of the flat slab differ substantially, with minimum estimates placing the boundary at the northern margin of the Himalaya (Indus-Yarlung Tsangpo suture), and maximum estimates placing it at the northern boundary of Tibet. In this study, we investigate for the first time if a flat slab segment of subducted buoyant Indian continental lithosphere below the Himalaya-Tibet region is geodynamically feasible and we quantify its northward extent, as well as its contribution to the high topography of the region. We conduct three large-scale fully-dynamic (buoyancy-driven) analogue experiments to simulate the subduction of the Indian continent. Our preferred, and geodynamically most feasible, model shows a continental flat slab extending northward up to ~ 320 km from the Himalayan thrust front, in agreement with recent estimates. Furthermore, it suggests that the positively buoyant flat slab segment of the Indian continent contributes some ~ 1.5–2 km to the high topography of the Himalaya-Southern Tibet region by providing an upward force to elevate the overriding Eurasian plate.

Geochronology and geochemistry of the El Salvador plutonic complex (Sierra de Tamaulipas, NE Mexico): cenozoic tectonic implications of the eastern Mexican alkaline province

ABSTRACT The El Salvador plutonic complex is a 9.2 km2 circular body within the Sierra de Tamaulipas, part of the Eastern Mexican Alkaline Province (EMAP). Three alkaline magmatic suites were identified from El Salvador from the geochemical analysis of 20 samples. Suite A (SiO2 = 58.3–65.3 wt.%) is ferroan (Fe* = 0.761–0.873) and alkalic (MALI = 4.83–9.89). In contrast, Suites B (SiO2 = 52.5–60.8 wt.%; Fe* = 0.680–0.756) and C (SiO2 = 50.5–67.1 wt.%; Fe* = 0.616–0.749) are magnesian. Suite B is alkali-calcic (MALI = -2.18–4.64), while Suite C is alkalic (MALI = 4.48–8.59). All suites display arc-related signatures. U-Pb and fission-track geochronology data reveal two uplift episodes during the cooling history of Suite A. One in the Late Eocene was based on U-Pb zircon (38.42 ± 0.21 Ma) and titanite ages (35.54 ± 3.77 Ma). The other was during the Oligocene from U-Pb apatite (29.9 ± 6.54 Ma) and fission-track titanite (30.2 ± 5.53 Ma) and apatite ages (32.7 ± 7.06 Ma). Integrating the arc-related signatures of El Salvador rocks with the well-documented Palaeogene arc magmatism of the Sierra Madre del Sur, we propose that the mantle beneath the EMAP experienced metasomatism during the Early Permian (and possibly the Early Jurassic) but not after the Late Cretaceous, ruling out Cenozoic slab subduction in eastern Mexico. In the absence of a slab to explain the El Salvador magmatism, we suggest a long-lasting, widespread mantle upwelling beneath Mexico’s northern half in response to the Farallon slab’s break-up. Under this context, the westward drift of the North American plate led to the Mexican Foreland Basin lithosphere reaching this massive mantle upwelling in late Eocene times to produce the EMAP magmatism.

Cretaceous magmatism in the northern Lhasa Terrane: Implications for the tectonic evolution and crustal growth tempos of central Tibet

The subduction and demise of an ocean plate are generally recognized as essential processes that result in the reworking and maturation of the continental crust. The northern Lhasa Terrane in central Tibet represents the forefront of the Lhasa-Qiangtang collision belt following the closure of the Bangong-Nujiang Ocean. Thus, it is a pivotal location to study the transition processes from oceanic lithosphere subduction to continental collision as well as pertinent crustal growth mechanisms. Here, we present zircon U-Pb dating, whole-rock major and trace element and Sr-Nd isotope, and zircon Hf isotope and trace element data of the Mendang igneous complex, Baingoin County, northern Lhasa Terrane. Geochronological results show that the granodiorites, trachydacites, and rhyolites in the Mendang igneous complex formed at ca. 122−116 Ma, 97 Ma, and 73 Ma, respectively. The Early Cretaceous granodiorite samples are peraluminous with high SiO2, Al2O3, and K2O contents, and moderate A/CNK (molar ratio of Al2O3/[CaO + Na2O + K2O]) values, which are similar to those of typical felsic peraluminous I-type granites. The granodiorites are characterized by enrichment in light rare earth elements and large ion lithophile elements (e.g., Rb) and depletion in high field strength elements (e.g., Nb, Ta). They also show the most enriched whole-rock Sr-Nd [(87Sr/86Sr)i = 0.7072−0.7078; εNd(t) = −7.60 to −5.08] and zircon Hf [εHf(t) = −4.46 to +1.02] isotope compositions, indicating that the Early Cretaceous granodiorites were likely derived from an ancient basement under a subduction setting. The trachydacites have uniform SiO2, high Al2O3, Sr contents, and Sr/Y values, and low Y and Yb contents, belonging to adakitic rocks. They show more depleted whole-rock Sr-Nd [(87Sr/86Sr)i = 0.7065−0.7066; εNd(t) = −0.56 to −0.22] and zircon Hf [εHf(t) = 4.36−7.84] isotopes than the granodiorites, suggesting that the trachydacites may have generated from partial melting of the juvenile thickened lower continental crust. The rhyolites have the highest SiO2 and K2O contents in the Mendang igneous complex, and significant depletion of Ba, Sr, Eu, P, and Ti. They have slighter more enriched whole-rock Nd [εNd(t) = −3.71 to −1.16] and zircon Hf [εHf(t) = 1.03−4.31] isotope compositions than the trachydacites. These features suggest that the rhyolites were highly fractionated products of the crustal melts. Whole-rock Sr-Nd and zircon Hf isotopes of the Mendang igneous complex show a kink trend from enrichment to depletion and then transfer to enrichment again, signifying an increased contribution of juvenile materials in the northern Lhasa Terrane toward progressively replacing the ancient lower crust and accumulating to newly formed crust. The estimated crustal thickness beneath the northern Lhasa Terrane shows a sharp increase from the Early to Late Cretaceous and peak at ca. 97 Ma, whereas it largely decreases in the Late Cretaceous. Integrated with previous studies, we propose that the formation of the Mendang igneous complex (122−73 Ma) elaborately documents the regional tectonic transition from oceanic lithosphere subduction to demise as well as continental crustal differentiation and maturation. The rollback and breakoff of the southward subducted Bangong-Nujiang oceanic slab in the Early Cretaceous initiated diverse magmatism in the northern Lhasa Terrane. During the early Late Cretaceous, the widespread adakitic and Mg-rich magmatism was attributed to the delamination of the thickened lithosphere following the final amalgamation of the Lhasa and Qiangtang Terranes. In the late Late Cretaceous, post-collisional extension induced the formation of the rhyolites in the northern Lhasa Terrane. The growth and destruction of the continental crust had certain tempos from the oceanic lithosphere subduction to continental post-collision stage.

Light oxygen isotopic composition in deep mantle reveals oceanic crust subduction before 3.3 billion years ago

Compositional heterogeneity exists in Earth’s deep mantle, which can be caused by the subduction of oceanic slabs. How early this process started on Earth remains highly debated due to the scarcity of early Archean materials with pristine mantle compositional signatures. Here, using the oxygen isotope and elemental compositions of fresh olivine grains in the 3.27-Ga komatiites of the Weltevreden Formation in the Barberton Greenstone Belt in Southern Africa, we discovered two groups of samples with primitive olivine grains. Group I exhibits normal mantle-like δ18O values and high Fo contents (δ18O = 4.9–5.4‰; Fo = 93–95); Group II is characterized by lower δ18O values with slightly lower Fo contents (δ18O = 3.6–4.7‰; Fo = 91–93). These δ18O values correlate with other geochemical proxies of olivine-poor iron-rich pyroxenite sources, indicating that the Weltevreden komatiites were derived from two distinct mantle sources. The existence of the low-δ18O magmas can be best explained by recycling of the altered oceanic crust into deep mantle arguably by subduction, which started 3.3 billion years ago and is responsible for the deep mantle heterogeneity in early Earth.

Subduction and retreat of Permian oceanic plates in the East Kunlun Orogenic Belt: evidence from mafic-ultramafic intrusive rocks

ABSTRACT The East Kunlun Orogenic Belt, an important part of the Central Orogenic Belt in China, preserved wealthy geological information for the tectonic evolution of the Proto- and Paleo-Tethys Ocean. Several Permian mafic-ultramafic rocks were identified in the Qimantagh area, EKOB, and zircon U-Pb dating, zircon Hf isotope, geochemistry and whole-rock Sr-Nd isotopes were studied in this paper. The results suggest that the formation age of the Xiyingzhaogou intrusion is 280 ± 4 Ma. The rocks show the characteristics of continental-arc calc-alkaline basalt. The whole-rock εNd(t) value of gabbro in Xiyingzhaogou intrusion is −4.2~-2.6 and the εHf(t) value of zircon is 2.77 ~ 5.82, indicating the characteristics of enriched or weakly depleted mantle; The formation age of Yingzhaogou intrusion is 263 ± 4 Ma. The geochemical characteristics of the Yinzhaogou rocks are consistent with those of continental-arc tholeiite. The whole-rock εNd(t) value of gabbro in Yingzhaogou intrusion is 0.73 ~ 0.92, and the εHf(t) value of zircon is 8.33 ~ 13.50, indicating the characteristics of depleted mantle. We believe that the Xiyingzhaogou intrusion was formed in the Early Permian and originated from the mantle wedge metasomatized by subduction fluid. The Yingzhaogou intrusion was formed in an extensional continental-arc environment caused by the retreat of the subduction plate. Asthenospheric mantle and fluid metasomatic lithospheric mantle were jointly involved in the mantle source. This study confirmed that the subduction of the Paleo-Tethys oceanic slab should had started before the Early Permian and a wide magma arc formed due to flat subduction. The retreat of the subduction plate led to the extension of the East Kunlun region, resulting in the significant enhancement of magmatism in the Middle Permian.

Heavy mineral indicators of ridge subduction beneath the Bainaimiao arc in the southeastern Central Asian Orogenic Belt

As one of the major subduction-accretionary systems in the southeastern Central Asian Orogenic Belt, the Bainaimiao arc presents a valuable opportunity for investigating the subduction–accretion process of accretionary orogens. It is located adjacent to the northern margin of the North China Craton and developed arc-related deposits during the Silurian–Devonian period. In this paper, we employed heavy mineral geochemistry and U–Pb ages of zircon and rutile, to determine the provenance of the Silurian–Devonian deposits and reconstruct the paleogeography of the source area. The stratigraphic units involved include the Xuniwusu Formation of Early Silurian and the Xibiehe Formation of Early Devonian.Both the Xuniwusu and Xibiehe formations exhibit heavy mineral assemblages dominated by zircons and tourmalines. Nevertheless, the Xibiehe Formation contains Cr-spinels and garnets which are absent in the Xuniwusu Formation. Geochemical analyses of the detrital Cr-spinels exhibit characteristics of a supra-subduction zone ophiolite. The garnets originated from middle- to high-grade metamorphic rocks, as well as mafic and felsic rocks. The tourmalines in both formations primarily indicate a metasedimentary rock source, with a small percentage of Li-poor granitoid species in the Xuniwusu Formation. The U–Pb age spectra of both detrital zircon and rutile in the Xuniwusu Formation show multiple peaks. The age spectra consist of Precambrian ages, indicative of a recycled origin, as well as Paleozoic ages, suggesting a magmatic origin. In contrast, the U–Pb age spectra of zircon and rutile in the Xibiehe Formation display a unimodal age peak at ca. 455 and 490 Ma, respectively. These differences suggest a provenance shift from the Xuniwusu to Xibiehe formations. The Xuniwusu Formation was sourced from the Bainaimiao arc itself, whereas the Xibiehe Formation has a mixed origin of the Bainaimiao arc and Ondor Sum accretionary complex. The provenance change correlates with the unconformity between the two formations, suggesting a forearc uplift event. In combination with the previous research on the Tulinkai ophiolite, it can be inferred that the forearc uplift was triggered by ridge subduction of the Paleo-Asian Ocean plate at ca. 450–430 Ma. This research provided a new perspective on the ridge subduction of the oceanic slab in fossil orogens.

Extremely enriched lithospheric mantle–derived magmas in the central Lhasa Terrane, southern Tibet

ABSTRACT Plate tectonics plays a crucial role in the evolution of Earth’s continents and oceans. However, the impact of oceanic subduction and continental collision on the continental lithospheric mantle within collisional orogenic belts is still a matter of debate. To address this issue, we performed a thorough analysis of mafic rocks in the Nixiong area within the central Lhasa Terrane, southern Tibet. The mafic rocks are composed of diabases and basalts formed during syn-collisional period (~50 Ma). The diabases have low SiO2 and moderate MgO contents, and are enriched in large ion lithophile elements and depleted in high field strength elements. They show high 87Sr/86Sri (0.711451–0.712084) and low εNd(t) (−7.9 to −6.44), as well as enriched Hf-Pb isotope compositions (εHf(t) = -7.58 to −6.59; 206Pb/204Pbi = 18.59–18.71; 207Pb/204Pbi = 15.71–15.72; 208Pb/204Pbi = 39.33–39.53). Geochemical modelling suggests that the diabases were formed through the partial melting of peridotite at high temperatures and shallow depths, with significant addition of subduction-related components to the primary melts. The diabases probably originated from a metasomatized lithospheric mantle due to subduction of Neo-Tethyan slab. The basalts display similar SiO2 but lower MgO contents compare with the diabases. They are also characterized by extremely high 87Sr/86Sri (0.713332–0.713448), low εNd(t) (−10.32 to −9.96), and significantly enriched Hf-Pb isotope compositions (εHf(t) = -12.30 to −12.05; 206Pb/204Pbi = 18.81–18.85; 207Pb/204Pbi = 15.74–15.75; 208Pb/204Pbi = 39.81–39.88), corresponding to an extremely enriched lithospheric mantle. Integrated with the previous studies, the Nixiong mafic rocks exhibit the most enriched isotopic compositions among the arc- or syn-collision-related mafic rocks in southern Tibet. These extremely isotopically enriched components were most likely derived from the ancient continental lithospheric mantle beneath the central Lhasa Terrane. Our results provided unambiguous evidence that the pre-existing ancient continental lithospheric mantle could persist locally during the prolonged oceanic subduction and subsequent continental collision.

Slab Subduction and Pull Link Magmatism at Active and Passive Continental Margins

AbstractSubduction leads to voluminous arc magmatism at active continental margins; however, the potential far‐field geodynamic connections with coeval magmatism along opposing passive margins, are rarely investigated. Here we report large‐scale and parallel mafic dike swarms on the South Qinling passive margin, whose geochemical features record the thick and cold subcontinental lithospheric mantle beneath passive margin was gradually replaced by thin and hot one during the early Paleozoic. The subduction of the Proto‐Tethys oceanic lithosphere led to arc magmatism along the active margin, shortly after the rift‐related mafic dike swarms occurred along the passive margin. Subsequent to the continuous oceanic subduction including the ridge‐margin collision, mantle upwelling and slab rollback eventually triggered large‐scale rift‐related magmatism and the opening of new ocean basin at passive margin, and also back‐arc magmatism on the active margin. We highlight slab subduction and pull link the geodynamic processes at active and passive continental margins, which improves our understanding of plate tectonics.