Mafic Lower Crust Research Articles

Multi‐Stage Magmatism During Slab Exhumation Drives the Geochemical Evolution of Continental Crust: Insights From Paleozoic Granitoids in South Altyn, Western China

AbstractThe present continental crust is characterized by a felsic upper crust and a mafic lower crust, resulting from significant geochemical differentiation over geological time. While various processes have been proposed to explain this differentiation, subduction zones remain pivotal regions for understanding the compositional evolution of continental crust. This study focuses on the South Altyn (SA) continental subduction‐collision belt in western China, a unique setting that experienced ultra‐deep (>300 km) continental subduction followed by multi‐stage exhumation. We present a comprehensive study of four granitoid suites from Tatelekebulake (TTLK) area in SA: biotite granite (BG), monzogranite (MG), K‐feldspar granite (KG), and leucogranite (LG). Comprehensive studies on petrology, geochemistry and zircon U‐Pb dating show that these granitoids formed at 494, 451, 414, and 418 Ma, respectively, and originated from protoliths with affinity to the subducted continental crust in SA. Phase equilibrium modeling suggests that BG formed at ∼800°C and 0.6 GPa, while the MG, KG, and LG formed by differentiation crystallization of the BG magma under progressively decreasing temperature and pressure conditions (750°C, 0.5 GPa; 740–700°C, 0.2 GPa; and 700–640°C, 0.1 GPa, respectively). These results, combined with previous studies, allow us to reconstruct the tectonic processes of continental exhumation and subsequent orogenic collapse in SA during the Early Paleozoic. Importantly, our findings reveal that magmatism derived from partial melting of subducted continental crust can promote the geochemical evolution of continental crust toward more felsic compositions, even in the absence of significant crustal growth or mantle‐derived magmatism. This study provides a valuable case for understanding the compositional evolution of continental crust in deep subduction zones and challenges conventional models that rely heavily on arc magmatism for crustal differentiation. Moreover, our results contribute to a broader understanding of crustal evolution processes in collisional orogens worldwide and highlight the importance of recycling and differentiation of subducted continental material in shaping crustal compositions.

The presences of water in the generation of calc-alkaline I-type granitic magmas at continental arc: Insight from Neoproterozoic Ailaoshan magmatism at the southwestern margin of the Yangtze Block, South China

Water is crucial in generating granitic systems; however, its role (i.e., dehydration, low to high water-fluxed melting) in granitic magmatism in continental arcs remains unsettled and poorly understood. Neoproterozoic arc-related igneous rocks are uncovered along the western margin of the Yangtze Block, presenting a complex compositional array ideal for deciphering the influence of water content within the continental arc and its geodynamic significance. The Jinping granites, in the Ailaoshan zone, are emplaced at 750 ± 4 Ma, as determined by zircon U-Pb dating. These medium to coarse-grained granites consist predominantly of plagioclase, quartz, K-feldspar, muscovite, and biotite. Characterized by high SiO2 (71.2–73.5 wt%), alkalis (K2O+Na2O=7.54–9.56 wt%) and low Fe2O3T (1.01–1.66 wt%), MgO (0.53–0.85 wt%), and CaO (0.15–0.93 wt%) concentrations, they exhibit high-K calc-alkaline signatures. The negative correlation between P2O5 and SiO2, alongside the positive correlation between Rb and Y, typifies the Jinping granites as I-type granites. Low La/Nb (1.96–3.43) and Nb/Ta (8.57–11.3) ratios, but high Th/La (0.28–0.46) and Zr/Sm (30.5–47.7) ratios, as well as whole rock εNd(t) (−0.4 to +1.3) and zircon εHf(t) values (+5.25 to +8.53) of the studied granites are similar to synchronous mafic rocks at the western margin of the Yangtze Block. These features suggest partial melting origin from the mafic lower crust. Thermodynamic modeling posits that the Neoproterozoic Ailaoshan I-type granitic rocks may have formed through water-fluxed melting (2.0–3.5 wt% H2O) under medium pressure conditions (9 kbar). It is postulated that slab rollback could have prompted water (as hydrous melt or fluids) release from hydrous minerals in the underlying cumulate mafic rocks, subsequently triggering water-fluxed melting in the lower crust. In contrast, high water-fluxed melting-generated adakitic granites with low K2O/Na2O ratios (<0.8) in the region, spatially and temporally associated with the Jinping granites, reflect higher water content in the deeper crust. This supports the notion that water was conveyed from depths to surface, facilitating the conversion of adakitic rocks into I-type granites as water content diminished. Thus, water content within the lower crust plays a pivotal role in the genesis of granitic rocks with varied compositions in a continental arc setting.

Formation of the Naxos nested domes and crustal differentiation by convection and diapirism

The Naxos dome, in the middle of the Aegean domain, exposes the former root of the Alpine orogenic belt and represents a key natural example to investigate the development of gravitational instabilities during orogenic evolution and their impact on crustal differentiation. The Naxos dome is cored by migmatites with structures depicting second order domes with a diameter of 1–2 km nested in the first order deca-kilometer scale dome that formed at the onset of orogenic collapse. Zircon grains from the migmatites record a succession of crystallization-dissolution cycles with a period of 1–2 Myr. These features have been attributed to the development of convective and diapiric gravitational instabilities, related to thermally induced and compositional buoyancy. In this paper, we test the pertinence of this model with a thermal-mechanical numerical experiment performed with a volume of fluid method (VOF) known to preserve material phase interfaces during large deformation of viscous layers. Partial melting of the crust is modeled by strain-rate and temperature dependent viscosity and temperature dependent density. Moreover, horizontal layers with density, viscosity and heat production variations mimic more felsic or more mafic lithologies in a crust of intermediate composition. With basal heating, gravitational instabilities initiate with local segregation of the buoyant versus heavier layers, followed by diapiric upwelling of buoyant pockets of aggregated less dense material. Convection starts after 5 Myr, approximately when half of the crust has a viscosity lower than 1019 Pa s. The size of the convection cells increases as the temperature rises in the crust and reaches ∼25 km in diameter after ca. 20 Myr, which defines the size of first order domes. Some of the heterogeneous material is entrained in the convection cells with a revolution period of 1 to 3 Myr. However, most of the denser material accumulates in the lower crust, while the buoyant material segregates at the top of the convection cells and forms diapirs that correspond to second order domes, of several kilometers in diameter and nested within the first order domes. This model, which reproduces the first order characteristic dimensions of the Naxos nested domes and the periodicity of their zircon geochronological record, demonstrates the efficiency of gravitational instabilities in the formation of migmatite domes and, more generally, in the multi-scale dynamics of crustal differentiation leading to a felsic upper crust, an intermediate middle crust and a mafic lower crust.

Deep Seismic Reflection Imaging of Mesozoic Kachchh Rift, NW India: Implications for Evolution

AbstractDeep crustal seismic reflection profiling is carried out for first time across various tectonic domains of seismically active Mesozoic Kachchh rift basin, formed during the breakup of Gondwanaland. The seismic data, processed using the common reflection surface stack approach, provided maiden images of the shallow and deep sub‐surface structures in the region. These images reveal a 15 km thick subhorizontal lower crustal reflection fabric and crustal‐scale domal‐type structure extending from the surface to the Moho. We interpret the earlier structure represent magmatic underplating and the latter as the Kachchh Mainland uplift. We find large variations in the thickness of sediments from 150 m to 6.5 km and crustal thickness from 45 to 35 km from north to south, with a Moho up warp of 4 km beneath the Kachchh Mainland fault. The Kachchh rift basin exhibits an unusually thick crust of 45 km, contrary to many rift basins. We interpret the syn‐rifting and the Reunion mantle plume activity, manifested as Deccan volcanics, are responsible for magmatic underplating and crustal thickening. Uplift in the region is multi‐genetic in origin. Present study illuminated new faults and nature of various other faults. Moderate to large earthquakes in the region are attributed to the regional and local stresses resulting from the plate boundary and heterogeneous crustal structure. Based on the mafic lower crust and distribution of aftershocks through the entire crust, we interpret the lower crust is brittle, contrary to most models of continental rheology.

The effect of in-situ metasomatism on the electrical resistivity of the lower crust

Improved estimates of composition and temperature in the lower crust are vital in characterising crustal rheology and better understanding tectonic processes. Fluids introduced into the lower crust during collisional events leave behind large domains of metasomatised crust which are long-lived and should be observable via magnetotelluric sounding. With this in mind, the electrical conductivity of a subduction-related eclogite and its anhydrous, anorthositic granulite protolith were experimentally determined at atmospheric pressure and lower crustal temperatures (773–1123 K). Complex impedance spectra were collected between 20 Hz and 2 MHz while oxygen fugacity was controlled about the quartz-fayalite-magnetite (QFM) buffer using a CO:CO2 gas mix. At <1000 K, both samples are characterised by impurity conduction with low observed activation enthalpies of 0.48 ± 0.03 eV and 0.39 ± 0.02 eV for the granulite and eclogite respectively. High temperature data from the granulite (>1000 K) indicates a switch to small polarons with an activation enthalpy of 1.02 ± 0.02 eV. Results indicate that conversion of Fe-poor mafic granulite to eclogite does not enhance conductivity at lower crustal temperatures. Comparisons with previous experiments reveal a much stronger relationship between conductivity and FeOT than previously anticipated. Modelling indicates that an increase in FeOT of 2 wt% is equivalent to a temperature increase of 100 K and reasonable variations in FeOT in the lower crust may account for ∼4 orders of magnitude variation in conductivity. Elevated FeOT up to 20 wt% may enhance conductivity to ∼100 S m−1 at a moderate 873 K. A system of linear equations is derived to estimate iron content, electrical conductivity or temperature within dry, mafic lower crust. For the temperature range 773–973 K, conductivity may be predicted to within 0.16 log units using log σ(T,FeOT) = 0.006 T + 0.3FeOT − 11.35. Under the same conditions, FeOT may be predicted to within 0.51 wt% using FeOT (T,σ) = 37–0.0193 T + 3.24 log σ. Finally, lower crustal temperatures may be estimated to within 24 K using T(FeOT,σ) = 1728 − 42FeOT + 139 log σ.

Coupled alkaline high-Nb mafic rocks and adakitic granodiorites: Products of Neoproterozoic back-arc extension at the western margin of the Yangtze Block, South China

The newly identified Neoproterozoic Yuanmou high-Nb mafic and adakitic rocks at the western margin of the Yangtze Block in South China are investigated for petrogenesis and geodynamic implications. Zircon U-Pb ages show the adakitic granodiorite host and mafic microgranular enclaves (MMEs) were synchronously emplaced at 811−806 Ma, and were cut by an 802 ± 11.5 Ma mafic dike. The medium-grained granodiorites consist mainly of plagioclase, quartz, K-felspar, biotite, and amphibole. High SiO2 (62.8−67.5 wt%), Al2O3 (13.1−16.8 wt%) contents, high Sr/Y (17.4−49.0), and (La/Yb)N (16.3−52.6) but low Y (6.83−16.1 ppm) and Yb (0.74−1.82 ppm) concentrations of the granodiorites point toward their calc-alkaline and adakitic nature. Further, their high Dy/Yb (1.52−2.31) and Gd/Y (2.35−4.43) and low K2O/Na2O (0.22−0.48) and Th/La (0.14−0.20) ratios, constant bulk rock εNd(t) (−0.5 to −1.5) and zircon εHf(t) (0.0 to +2.3) values, suggest that the granodiorites were derived from partial melting of a mafic lower crustal source. The MMEs in the granodiorites and mafic dike intrusions are composed of amphibole, plagioclase, clinopyroxene, and olivine, and have similar chemical compositions, such as high Nb (15.7−41.9 ppm) and TiO2 (2.13−3.39 wt%) contents and high Nb/U (18.7−30.9) values, akin to alkaline high-Nb basalts (HNB). Their variable Ba/Nb (7.68−79.2), Nb/Zr (0.09−0.16), high Nb/Th (5.47−6.68), and low Th/Zr (0.01−0.03) and La/Nb (1.03−1.37) ratios indicate a mantle source modified by melts derived from subducted oceanic crust and sediments. The Yuanmou HNB with high εNd(t) (+4.8 to +6.9) were probably generated during back-arc extension by decompression melting of a mantle source metasomatized by adakitic slab melts formed during earlier stages of the subduction. This implies that the location of the back-arc was far from the subducted oceanic slab, possibly due to a slab rollback process. In this scenario, HNB intrusions provided the heat for partial melting and adakitic melt formation in the mafic lower crust. The model suggests that the Yuanmou and other similar associations of HNBs and adakitic rocks emplaced along the western margin of the Yangtze Block are products of a widespread Neoproterozoic back-arc extension.

Assessing crustal contamination in the 28 to 18 Ma Dulce dike swarm with Nd-Sr isotopic data, southwestern Colorado

Alkaline to subalkaline mafic dikes in the 28 to 18 Ma Dulce swarm were emplaced in a zone of incipient extension from southern Colorado into northern New Mexico on the northeastern boundary of the San Juan Basin. The 87Sr/86Sr ratios for the dikes are 0.70503 to 0.70584, akin to most post-28 Ma mafic rocks across the northern San Juan Basin. These data are consistent with melting of metasomatized subcontinental lithospheric mantle with little to no crustal contribution as revealed by the geochemical and Sr-Nd isotopic signatures of most 28–0.6 Ma mafic rocks in the region. Time-corrected εNd(t) values of −4.1 to −7.4 for rocks in the Dulce swarm, however, indicate that magma production involved the crust. A previous hypothesis for Dulce magmas was contamination of lithospheric mantle melts with up to 45% mafic lower crust ± 0.5% upper crust. In this investigation, six new whole-rock Sr and Nd isotopic analyses were combined with published data to further investigate the contamination of lithospheric mantle melts with different crustal reservoirs. The Nd isotope signatures of the Dulce swarm offer evidence for the long-term involvement (∼10 Ma) of lower crust in the production of rift-related mantle magmas. Isotopic mixing curves support previous hypotheses for the contamination of lithospheric mantle melts with 10 to 40 percent lower mafic crust. This provides further insight into regional variations in mantle magmas produced after 28 Ma in the Four Corners region that likely triggered crustal melting related to caldera complexes in the western San Juan Mountains.

Crustal maturation and cratonization in response to Neoarchean continental collision: The Suizhong granitic belt, North China Craton

The late Archean (3.0–2.5 Ga) is a pivotal period, when properties of the continental crust significantly changed and a fundamental shift of geodynamic processes took place. As the major felsic component in early Archean (> 3.0 Ga) crustal records, tonalite-trondhjemite-granodiorite (TTG) suites are compositionally different to the present-day felsic upper continental crust. But it remains unclear how and to what degree the felsic part of the continental crust evolved during the late Archean. The 250 km long Neoarchean Suizhong granitic belt in the North China Craton (NCC) is mainly composed of K-rich granitoids and compositionally distinct to the TTG-dominated early Archean upper continental crust. This belt is a natural laboratory to explore how the early Archean upper continental crust evolved in the late Archean. Combined with our new data (zircon U-Pb ages, bulk-rock geochemistry and zircon Hf isotopes), we provide an overview of the temporal and compositional features of the Neoarchean Suizhong granitic belt. Granitic magmatism in the Suizhong granitic belt mainly took place within 100 Myrs around the end of the Archean, and these rocks are primarily K-rich granitoids with minor TTGs and mafic magmatic enclaves (MMEs). The K-rich granitoids include sanukitoid-like granites and potassic granites; the former are hybrids between melts from Mesoarchean enriched mafic crust and enriched mantle, and the latter are melts from Hadean–Mesoarchean crustal lithologies. The volumetrically minor TTGs were sourced from Mesoarchean low-K mafic crust. MMEs are either cumulates or inclusions of captured enriched mantle melts. Overall, these Neoarchean granitoids are essentially intracrustal reworking products of ancient continental nuclei and are compositionally and temporally similar to syn- and post-collisional magmatism in Phanerozoic collisional orogenic belts. The Neoarchean crustal reworking process was most likely to be induced by amalgamation of micro-continents through collision. Through the collision-induced reworking, a variety of K-rich granitoids were introduced to the Neoarchean upper continental crust. As a result, layering of felsic upper and mafic lower crust was developed, and the maturity of Archean continental crust was significantly enhanced. Collisional orogenesis played a critical role in maturing the Archean continental crust and initial cratonization of the NCC during the late Archean.

Pre-plate tectonics and origin of continents

<p indent="0mm">Established in the mid-1960s and regarded as a revolution in Earth sciences, the plate tectonics theory can reasonably interpret nearly all geological phenomena, processes, and events that happened during post-Archean time (from 2.5 billion years ago to the present), and has also been applied to explain the formation and evolution of continents. According to the plate tectonics theory, the mafic lower crust and the felsic upper crust of continents can develop from an island arc that formed by subduction of one oceanic crust beneath another, where the mafic lower crust of continents can be extracted from the mantle through the partial melting of the mantle wedge in the subduction zone, whereas the felsic upper crust of continents can form by the partial melting of the already-formed mafic lower crust. However, such an island arc model under a plate tectonic regime cannot well explain the magmatic, metamorphic and structural features of Archean (&gt;2.5 billion years) continents. For example, island arc models fail to explain the presence of ~1600°C komatiites, absence of andesites that is dominant in post-Archean arcs, nearly coeval and craton-scale emplacement of tonalite-trondhjemite-granodiorite (TTG) rocks, dominancy of dome-and-keel structures, and lack of ultrahigh-pressure rocks, paired metamorphic belts and ophiolites, etc. All of these imply that Archean continents may not have been formed under plate tectonic regimes, but were originated from some pre-plate tectonics (non-plate tectonics). So far, researchers have proposed a number of pre-plate tectonics models, of which the most representative ones are mantle plumes, sagduction, heat-pipes and stagnant lids. Although each of these pre-plate tectonics models can satisfactorily interpret some features of Archean cratons, none of them is successful in explaining all lithological, structural and metamorphic features of Archean cratons. For example, although the mantle plume-derived oceanic plateau models can well explain many features of Archean cratons, oceanic plateaus formed by mantle plumes may not provide enough water (H₂O) for aqueous partial melting of basaltic rocks to create TTG magmas. It is the same case with heat-pipe models. As the sagduction models assume the existence of an old felsic continental crust, they are not suitable for discussing the origin of Archean continents. As for stagnant lids, they just describe the state of a single lithosphere plate which itself cannot provide any geodynamic mechanisms for making a mafic crust be partially melt to form TTG magmas. Therefore, none of available pre-plate tectonics models proposed so far has been successfully applied to interpret the origins of Archean continents. Thus, a number of research groups in the world are conducting extensive and comprehensive investigations on this important scientific conundrum. Although Chinese geologists missed a chance to have made contributions to establishing plate tectonic theory in the 1960s, they have a great potential for breakthroughs in establishing a pre-plate tectonics theory since they have done tremendous work on the early geodynamic mechanisms for the formation and evolution of Archean continents and produced large amounts of new data and competing interpretations in the past four decades.

Petrogenesis of Early Cretaceous granitoids and mafic enclaves from the Jiaodong Peninsula, eastern China: Implications for crust-mantle interaction, tectonic evolution and gold mineralization

Early Cretaceous granitic plutons in the Jiaobei Terrane, Eastern China, contain abundant mafic enclaves (MMEs), which are important to understand the tectonic settings and genesis of gold mineralization in the Jiaodong Peninsula. Here, we combine petrology, geochronology, and isotopic geochemistry for MMEs and host granitoids to constrain their tectonic evolution and understand the genesis of coeval gold mineralization in the Jiaodong Peninsula during the Early Cretaceous. Zircon U-Pb dating reveals two magmatic events: (1) the earlier magmatic event (130–126 Ma), represented by Congjia MME (CJ-MMEs) and host granodiorite, and (2) the later magmatic event (121–116 Ma), including Gushan MME (GS-MMEs) and host granite. Field relationship, disequilibrium textures, and elemental and Sr-Nd-Hf isotopic compositions indicate that MMEs and host granitoids formed through mixing/mingling of mafic and felsic magmas. However, significant differences between CJ-MMEs and GS-MMEs indicate that the former has adakitic affinity and was produced by mixing of magmas from lithospheric mantle and mafic lower crust in deep hot zone, whereas the latter has arc-magmatic affinity, suggesting enriched lithospheric mantle source. Higher SiO2 and lower MgO contents in the host granitoids indicate that they were originated from reworking ancient NCC lower crust. Our results reveal crustal growth during the earlier magmatic event, followed by rapid crustal/lithospheric thinning in the later magmatic event, all triggered by ongoing rollback of subducting Paleo-pacific plate. These results combined with complex assemblages of sulfides and volatiles in GS-MMEs indicate that the lithospheric thinning, the rise of enriched lithospheric mantle-derived melts to shallow crustal levels, and crustal-mantle interactions around ∼ 120 ± 3 Ma, are the main factors for massive gold precipitation in the Jiaodong Peninsula.

Geochronology, petrogenesis, and tectonic significance of the granites in the Chaqiabeishan area of the Quanji Massif, northwestern China

The Chaqiabeishan area of the Quanji Massif in NW China is characterized by extensive Li‐rich pegmatites and granites of various ages. In this study, the major and trace element compositions, U–Pb zircon ages, and zircon Hf isotopic compositions of the granites were investigated, including the Garideng porphyritic monzogranite, the Gongkaxiuma monzogranite, and the Narong granitic porphyry. Precise U–Pb zircon dating revealed that these three granites were emplaced at approximately 468, 444, and 245 Ma, respectively. The geologic and geochemical data indicate that the Garideng porphyritic monzogranite and the Narong granitic porphyry are two I‐type granites, and the Gongkaxiuma monzogranite is an A2‐type granite. Based on the Hf isotopes and geochemical data, the Garideng porphyritic monzogranite and the Narong granitic porphyry were derived from partial melting of the relatively deep and relatively shallow juvenile mantle‐derived mafic lower crust, respectively, while the Gongkaxiuma monzogranite was derived from partial melting of the granulitic metamorphic basement. Moreover, tectonically, the Garideng porphyritic monzogranite represents the northward subduction of the South Qilian oceanic plate from 500 to 440 Ma. The Gongkaxiuma monzogranite represents the closure of the South Qilian Ocean by ~440 Ma, and the Narong granitic porphyry represents the Late Permian to Middle Triassic southward subduction of the Zongwulong oceanic plate under the Quanji Massif.

The Early Cretaceous tectonic evolution of the southern Great Xing’an Range, northeastern China: new constraints from A2-type granite and monzodiorite

The widespread Early Cretaceous plutons intruding along the southern Great Xing’an Range (SGXR) provide evidence for tectonic evolution of the region. Petrological, geochemical, zircon U–Pb geochronology, and zircon Hf isotopic studies are conducted on intrusions from Bianjiadayuan and Hongling areas. These suites classify as A2-type granites and monzodiorites, respectively. The 138–133 Ma A2-type granites originated from partial melting of continental crustal materials at high temperatures and shallow depths with significant addition of juvenile mafic lower crust sourced from a metasomatized mantle. The 136–134 Ma monzodiorites originated from the partial melting of an enriched mantle that was modified by melts of a previously subducted slab coupled with crustal contamination. The Early Cretaceous magmatism in the SGXR occurred in two periods: ∼145–136 Ma (peak at ∼139 Ma; εHf (t) = 5 to 10) and ∼136–130 Ma (peak at ∼131 Ma; εHf (t) = –10 to 15). The Early Cretaceous granite–monzodiorite suite in the SGXR suggests a bimodal magmatism in an extensional setting. The ∼145–130 Ma magmatism may have been triggered by asthenospheric upwelling induced by the Mongol–Okhotsk oceanic slab breakoff and large-scale lithospheric delamination resulting from post-orogenic extension. The variation of subduction direction of the Paleo-Pacific Ocean likely triggered a change in stress regime at ca. 136 Ma and likely promoted the lithospheric delamination beneath the SGXR resulting in intense magmatism originating from various sources. As such, the Paleo-Pacific Oceanic subduction likely played an important role in the Early Cretaceous magmatism in the SGXR.

Petrogenesis of Tertiary granitoid rocks from east of the Bidhand fault, Urumieh-Dokhtar Magmatic Arc, Iran: Implication for an active continental margin setting

Tertiary granitoid rocks including diorite, quartz diorite, quartz monzodiorite, monzonite, and quartz monzonite occur in southeast of Qom and to the east of Bidhand fault in the central Urumieh-Dokhtar Magmatic Arc (UDMA). Zircon U–Pb dating reveals an age of ca. 17 Ma for the crystallization of Kerogan quartz diorite. The study samples show sub-alkaline to transitional, metaluminous, and I-type affinity with low to middle SiO2 (53.9–63.9 wt%). Micro granitoid enclaves (MGEs) are widespread in these rocks, exhibiting diorite to monzodiorite composition (SiO2 = 56.75–58.4 wt%). MGEs and host rocks are slightly enriched in LREEs and LILEs, and depleted in HFSE including Nb, Ta and Ti as expected in granitoids formed in an active continental margin setting. The granitoid rocks display initial Sr isotopic compositions of 0.7047–0.7062, positive εNd(t) of 3.2–6.6, and 206Pb/204Pb (18.6431–18.8499), 207Pb/204Pb (15.6322–15. 6902), and 208Pb/204Pb ratios (38.5585–38.9325), and positive zircon εHf(t) values of 6.7–8.7. MGEs also display identical Sr-Nd isotopic compositions (87Sr/86Sri = 0.70537–0.70576; εNd(t) = 2.3–5.3) to the host rocks. The isotopic characteristics in combination with the geochemical signatures indicate that the granitoid rocks and MGEs most likely originated by interaction between mantle-derived mafic and basaltic lower crust-derived melts. Major and trace element modeling reveals that the mafic lower crust is composed of less ~40% partial melting of garnet-free amphibolite. Petrographic and geochemical characterization together with bulk rock Nd-Sr isotopic data suggest that host rocks and associated enclaves derived from the same parent by interaction between basaltic lower crust (~10%) and lithospheric mantle (~90%) during their generation that subduction related compositions have also been involved in their petrogenesis. This magma has experienced high fractional crystallization but minor crustal contamination during magma ascent to the upper crustal level. Petrographical, geochemical, and Sr-Nd-Pb isotopic similarities between granitoid and MGEs imply that MGEs are cognate late cumulates, formed by pressure quenching mechanism during the late stage of magma evolution in crustal magma chamber. The Tertiary magmatism in the study area originated in an arc-related setting induced by slab roll-back of northward Neo-Tethyan oceanic lithosphere subduction beneath central Iranian micro continent and resulted in lithospheric extensional regime and subsequently asthenosphere upwelling. The heat induced by the upwelling of the asthenosphere likely led to melting of upper mantle and the mafic crustal material during Early Miocene as the final magmatism phase in the central part of the UDMA before collision.

Volcanic rocks of the Elashan Formation in the Dulan-Xiangride Basin, East Kunlun Orogenic Belt, NW China: Petrogenesis and implications for Late Triassic geodynamic evolution

ABSTRACT The volcanic rocks of the Elashan Formation provide important evidence for the Late Triassic post-collisional magmatism in the East Kunlun Orogenic Belt. This study focused on the geochronology, geochemistry, and petrogenesis of volcanic rocks of the Elashan Formation in the Dulan–Xiangride Basin. The 227.9 ± 0.4 Ma dacite exhibits calc-alkaline I-type characteristics, with low Sr/Y and (La/Yb)N ratios and prominent negative Eu anomalies. We suggest that the dacite was derived from the partial melting of lower crust at normal crustal thickness. The crystalline tuff (225.0 ± 0.4 Ma) has negative εHf(t) values, implying that it was likely formed by the partial melting of the ancient crust. The rhyolites (221.3 ± 1.1 and 221.9 ± 0.6 Ma) are peraluminous and calc-alkaline to high-K calc-alkaline. The high SiO2 and Na2O + K2O contents, and high K2O/Na2O and Ga/Al ratios of these rhyolites have an affinity with A-type granite. The rhyolites were formed by the partial melting of Late Paleoproterozoic–Mesoproterozoic lower crust induced by mantle-derived mafic magma underplating. The andesite and basaltic andesites formed at 221.1–217.9 Ma. The zircon Hf isotopic characteristics show that the andesite was likely derived from the partial melting of the Late Paleoproterozoic mafic lower crust. The combination of our new data and previously obtained geological data suggests that the East Kunlun Orogenic Belt underwent three magmatic activities corresponding to three tectono-magmatic events (slab breakoff, the delamination of unstable lower lithospheric mantle, and the intense asthenosphere upwelling with further delamination and continuous crustal thinning) during the Late Triassic.

Tracking Prototethyan assembly felsic magmatic suites in southern Yunnan (SW China): evidence for an Early Ordovician–Early Silurian arc–back-arc system

Early Paleozoic trondhjemites, gneissic granites and alkali granites in southern Yunnan preserve important records of the tectonic evolution of the Prototethyan Ocean and regional correlations. Zircon ages suggest that these granitoids were emplaced from 476 to 436 Ma. The trondhjemites are characterized by high Na 2 O and low K 2 O contents, with ε Nd ( t ) values of −1.9 to −3.5 and ε Hf ( t ) values of −2.8 to +3.9. The trondhjemites were derived from an amphibolite source with a juvenile mafic component. The gneissic granites belong to the metaluminous low-K calc-alkaline series with an ε Nd ( t ) value of −6.2 and ε Hf ( t ) values of −5.0 to −0.4. The alkali granites belong to the high-K calc-alkaline series and yield ε Nd ( t ) values of −10.1 to −10.7 and ε Hf ( t ) values of −7.9 to −2.3. The gneissic granites were derived from an ‘ancient' lower mafic crust, whereas the alkali granites were derived from a meta-sedimentary source. These granitoids were formed during the subduction of the Prototethyan Ocean beneath the Simao Block and can be compared with similar igneous rocks from the Truong Son and Tam Ky-Phuoc Son zones in southern Laos. Our study, along with Early Paleozoic igneous suites from southern Laos, central Vietnam and the Malay Peninsula, suggests an arc–back-arc system along the northern margin of Gondwana. Supplementary material: Tables of zircon U–Pb and in-situ Hf and geochemical data are available at https://doi.org/10.6084/m9.figshare.c.5322386

Petrogenesis and tectonic evolution of the Palaeozoic to Mesozoic Niuxinshan granitoids in the North Qilian orogen,NWChina

Numerous granitoid suites exposed in the Qilian orogen, NE Tibetan Plateau preserve important records of the tectonic evolution of this orogen. Here we present zircon U–Pb, zircon rare earth elements (REEs), and whole‐rock geochemical and Sr–Nd–Hf isotopic data on Niuxinshan granitoids in the North Qilian orogen, with a view to constrain their petrogenesis and timing of emplacement, as well as to gain new insights into the Early Palaeozoic tectonic evolution of the Qilian orogen. Zircon U–Pb data show multiple concordant ages of ~517, 455–425, 397–389, and 351–319 Ma for syenogranite, and ages of ~758 and 724 Ma together with some scattered Early Mesozoic discordant ages for porphyritic granite, suggesting episodic magmatism from Neoproterozoic to Early Mesozoic. Zircon ∑REEs range from 276.43 to 4,582.32 ppm, and are characterized by depletion of light REEs (LREEs) and enrichment heavy REEs (HREEs) with negative Eu anomalies, indicating a mixture of magmatic and hydrothermal zircons. Whole‐rock geochemical data show similar variation of LREEs, large‐ion lithophile elements (LILEs), HREEs, and high‐field‐strength elements (HFSEs) with zircon REEs, with Eu, Ba, P, and Ti negative anomalies and Rb, U, and Pb positive anomalies. Based on their formation ages, two groups of granitoids are identified. Group I rocks (Early Palaeozoic) are represented by I‐type granites whereas Group II rocks (Early Mesozoic) are of S‐type. Whole‐rock Sr–Nd–Hf isotopes display (87Sr/86Sr)iratios of 0.707085–0.709591, εNd(t) values of −5.3 to −2.3, εHf(t) values of −2.2 to 0, and two‐stage model (Nd and Hf) ages of 1.6–1.4 Ga for Group I rocks, and of 0.728156–0.734113, −9.4 to −9.3, −10.6 to −6.9, and 2.1–1.9 Ga for Group II rocks. Accordingly, we infer that the magma for Group I was sourced from the partial melting of basaltic protoliths in the mafic lower crust and underwent fractional crystallization. Group II was derived from metasedimentary protoliths through reworking of ancient Palaeoproterozoic crustal components. Integrating the results from this study with those from previous studies on Niuxinshan granitoids, we argue that the southward subduction model of the North Qilian Ocean beneath the Qilian Block during the Early Palaeozoic might be a more reasonable scenario. However, the Early Mesozoic magmatic event reported here remains equivocal in the North Qilian orogen.

The influence of variations in crustal composition and lithospheric strength on the evolution of deformation processes in the southern Central Andes: insights from geodynamic models

Deformation in the orogen-foreland system of the southern Central Andes between 33° and 36° S varies in style, locus, and amount of shortening. The controls that determine these spatially variable characteristics have largely remained unknown, yet both the subduction of the oceanic Nazca plate and the strength of the South American plate have been invoked to play a major role. While the parameters governing the subduction processes are similar between 33° and 36° S, the lithospheric strength of the upper plate is spatially variable due to structures inherited from past geodynamic regimes and associated compositional differences in the South American plate. Regional Mesozoic crustal horizontal extension generated a < 40-km-thick crust with a more mafic composition in the lower crust south of 35°S; north of this latitude, however, a more felsic lower crust is inferred from geophysical data. To assess the influence of different structural and compositional heterogeneities on the style of deformation in the southern Central Andes, we developed a suite of geodynamic models of intraplate lithospheric shortening for two E–W transects (33° 40′ S and 36° S) across the Andes. The models are constrained by local geological and geophysical information. Our results demonstrate a decoupled shortening mode between the brittle upper crust and the ductile lower crust in those areas characterized by a mafic lower crust (36° S transect). In contrast, a more felsic lower crust, such as in the 33° 40′ S transect, results in a coupled shortening mode. Furthermore, we find that differences in lithospheric thickness and the asymmetry of the lithosphere–asthenosphere boundary may promote the formation of a crustal-scale, west-dipping detachment zone that drives the overall deformation and lateral expansion of the orogen. Our study represents the first geodynamic modeling effort in the southern Central Andes aimed at understanding the roles of heterogeneities (crustal composition and thickness) at the scale of the entire lithosphere as well as the geometry of the lithosphere–asthenosphere boundary with respect to mountain building.

Rock probes and seismic methods to constrain the structure, composition and evolution of the deep crust beneath North China Block

<p indent="0mm">Deep continental crust archives critical information on the bulk crustal composition as well as its formation. It is commonly the inherent zone for complex crust-mantle interaction. Unlike the shallow crust, the deep crust is generally inaccessible for direct geological sampling. Instead, knowledge of the deep crust can be obtained from a combination study of geophysical techniques and “lithoprobe” (deep crustal fragments entrained in magmatic outcrops). Wide-angle reflection/refraction seismic profiles have the highest resolution on crustal velocity structure and provide information on the averaged present-day nature of large deep crust domains; “lithoprobe”, instead, give detailed information of locality-specific deep crustal sections through determining the physical and chemical properties, dating the formation and alteration ages and constraining the melting and equilibration depths of the xenoliths. Here, the two methodologies are integrated to constrain the crustal structure on the North China Block. The seismic crustal structure model of the block is compiled from all available seismic wide-angle reflection/refraction profiles and receiver function results. The model parameters include Moho depth, <italic>V</italic>p velocity and thicknesses for four crustal layers (sedimentary cover, and the upper, middle and lower crust). The crust as a whole is thin (&lt; <sc>36 km)</sc> in the eastern part relative to the western part <sc>(38−44 km).</sc> In addition, the crust is shown to have a thick upper and middle crust and a relatively thin (typically <sc>&lt;5 km)</sc> lower mafic crust; no spatial variations in average crustal <italic>V</italic>p velocity and densities are observed. The average crustal <italic>V</italic>p and density for the whole North China Block are <sc>6.33−6.37 km/s</sc> and <sc>2.72−2.78 g/cm³,</sc> respectively. The overall absence of high P-wave velocity layer and the relatively low average <italic>V</italic>p and density indicate that the block has an intermediate to felsic bulk crustal composition. In contrast, regional crustal xenoliths include mafic to felsic granulites and amphibolites and rare eclogite xenoliths (e.g., Xinyang and Xuhuai), indicating the presence of thick lower crust (e.g., <sc>&gt;45 km)</sc> at that time. Meanwhile, mafic xenoliths in Hannuoba and Nüshan show strong mafic underplating at the lower crust in the Meso-Cenozoic while no high-velocity layer is shown in the seismic crustal model. These inconsistences can be explained by the fact that the two methodologies tap the crustal information at different timeslices and dimensions and be complementary to understand the chemical and physical structure of the deep crust. Therefore, the combination study of petrological and geochemical of the crustal xenoliths and geophysical observations can provide essential information on making and destroying the deep continental crust through time.

Preliminary result for crustal properties derivation related to tectonics for hazard mitigation in Eastern Indonesia using Teleseismic P Coda

Eastern Indonesia is tectonically complex, formed by different plates and microplates interactions from different origins. This complexity gives geoscientists a challenge to solve the ’jigsaw’ of the complex interactions. The understanding of tectonic processes can lead to a breakthrough in both resource exploration and disaster risk reduction. We utilize teleseismic P wave coda for random coda from scattering and deterministic coda originated from the crust-mantle boundary (Moho) to derive the crustal properties, including thickness, Vp/Vs, and qualitative scattering characteristics. For the scattering properties, we apply Iterative Cross-Correlation and Stacking (ICCS) to align the waveform. At the same time, for the crust characteristic, we employ the Receiver Functions (RF) method alongside H-k stacking. The crustal thickness recovered from the RF and H-k stacking has a good correlation with the crustal origin, where the thickness in older and stable crust originated from Sundaland and Gondwana is thicker than a younger plate of the crust arc and subduction origin. The Vp/Vs is high in a region that is interpreted to be dominated by mafic lower crust originated from oceanic-oceanic subduction during Eocene, anisotropy, or by a magmatic anomaly. The P coda also correlated well with the subsurface magmatic anomaly by providing a unique pattern.

Geochemical and isotopic study of metasomatic apatite: Implications for gold mineralization in Xindigou, northern China

The 0.7 Moz Xindigou gold deposit lies on the north margin of North China Craton. The orebodies are controlled by NWW-trending shear zones and hosted by an Archean metavolcanosedimentary sequence. These host schist, amphibolite and marble rocks were intruded by late Paleozoic granites and mafic dikes. Two independent magmatic-hydrothermal assemblages were identified based on petrographic evidence: a pre-ore magmatic-hydrothermal magnetite – apatite ± pyrite ± zircon assemblage and an ore-related potassium feldspar – muscovite – sulfide ± minor monazite assemblage.Apatite that has undergone later alteration shows significant intracrystalline spatial variability with respect to the concentrations of rare earth elements (REEs), Si, and Sr. Scanning electron microscope backscattered images and statistical interrogation of the compositional data acquired by LA-ICP-MS indicate that apatite grains in the pre-ore magnetite – apatite assemblage are homogeneous and enriched in REEs and Si (Ap1). Where sulfide-related alteration of mineralization stage was overlapped, these apatite grains display a distinctive texture of homogeneous cores and turbid rims. This dissolution-reprecipitation process preferentially transported elements from the original apatite. REEs and Si gradually decrease from the core (Ap1) to the altered rim (Ap2 & Ap3) whereas Sr increases. Localized leaching of REEs (especially LREEs) from apatite was followed by precipitation of secondary monazite grains in the metasomatic apatite rims and the muscovite – sulfide – quartz – phosphate sequence. In situ Sr isotope analyses of the apatite show the 87Sr/86Sr ratios varying from the lowest value, 0.707874 in unaltered Ap1, to the highest, 0.714160 in altered Ap3. Zircon grains from migmatite display a peak metamorphic age at 2512 ± 9 Ma, whereas zircons from the early magnetite-apatite association formed near 1776 ± 12 Ma, and late hydrothermal monazite yields a Th-Pb age of 301.0 ± 1.7 Ma.These geochronology results show that the mineralization occurred during the late Paleozoic orogenesis, long after Neoarchean regional metamorphism and late Paleoproterozoic rift-related magmatism. The typical hydrothermal features of orogenic gold deposit could be inferred from the fluorapatite. The elemental patterns of fluorapatite grains indicate that a carbon and Sr enriched, acid, slightly oxidized and low salinity hydrothermal system was responsible for apatite metasomatism and gold mineralization. Strontium isotopes further reveal that the fluid related to gold deposition was derived from mafic lower crust, with extensive radiogenic Sr accumulation during fluid-rock interaction. This study demonstrates the potential application of apatite in exploration for orogenic gold deposits and its use in tracking hydrothermal REE re-mobilization during metasomatism.