Crustal Growth Research Articles

The Paleocene high-silica granites and their associated diorite enclaves in central Lhasa subterrane, Xizang, China: Insights to crust-mantle mixing and crustal reworking

The history of crustal growth in the western Central Lhasa (CL) subterranean during India-Eurasian collision remains unclear. The crucial problem is lack of suitable and ideal magmatism recorded information of the early stages of continental collision. Here, we report the newly discovered high-silica granites (K-feldspar granite with aplite granite) hosting diorite enclaves from the Geji area in western CL subterrane. This suite of magmatic rocks is of great significance to reveal the crustal growth of CL during the India-Eurasia collision. Zircon U-Pb dating results show that the K-feldspar granite, aplite granite, and diorite enclaves were formed at 57.2–57.9 Ma, 58.7 Ma and 58.2 Ma, respectively. Petrology, mineralogy, and mineral chemistry of the plagioclase from the diorite enclave indicates that it was generated by magma mixing between mafic and felsic magmas. The mafic melt was sourced from partial melting of mantle wedge and metasomatized by subduction fluid and sediment melt. This has been supported by the positive zircon εHf(t) values (+0.6 to + 4.9), enrichment of large ion lithophile elements (Rb, K, Th, and U) and depletion of high field strength elements (Nb and Ta), and the high Th/Ce ratios of diorite enclaves. The K-feldspar granite also has positive zircon εHf(t) values (+1.3 to +6.2), arc magma geochemical features, and a lower crust Nb/Ta ratio (8.3). These indicate that the K-feldspar granite was originated from the juvenile lower crust. In conclusion, we think that the crust-mantle magma mixing has a substantial impact on the crustal growth of the western CL subterrane during the subduction slab roll-back of the Paleocene Neo-Tethys Ocean when paired with other research and the context of regional tectonic development.

Building and Collapse of the Cadomian Orogen: A Plate‐Scale Model Based on Structural Data From the SW Iberian Massif

AbstractThe Cadomian Orogeny produced a subduction‐related orogen along the periphery of Gondwana and configured the pre‐Variscan basement of the Iberian Massif. The architecture of the Cadomian Orogen requires detailed structural analysis for reconstruction because of severe tectonic reworking during the Paleozoic (Variscan cycle). Tectonometamorphic analysis and data compilation in SW Iberia (La Serena Massif, Spain) have allowed the identification of three Cadomian deformation phases and further constrained the global architecture and large‐scale processes that contributed to the Ediacaran building and early Paleozoic dismantling of the Cadomian Orogen. The first phase (DC1, prior to 573 Ma) favored tabular morphology in plutons that intruded during the building of a continental arc. The second phase (DC2, 573–535 Ma) produced an upright folding and contributed to further crustal thickening. The third phase of deformation (DC3, ranging between ∼535 and ∼480 Ma) resulted in an orogen‐parallel dome with oblique extensional flow. DC1 represents the crustal growth and thickening stage. DC2 is synchronous with a period of crustal thickening that affected most of the Gondwanan periphery, from the most external sections (Cadomian fore‐arc) to the inner ones (Cadomian back‐arc). We explain DC2 as a consequence of flat subduction, which was followed by a period dominated by crustal extension (DC3) upon roll‐back of the lower plate. The Ediacaran construction of the Cadomian Orogen (DC1 and DC2) requires ongoing subduction beneath Gondwana s.l., whereas its dismantlement during the Early Paleozoic is compatible with oblique, sinistral convergence.

A global zircon U–Th–Pb geochronological database

Abstract. Since the start of the 21st century, the widespread application of ion probes has promoted the mass output of high-precision and high-accuracy U–Th–Pb geochronological data. Zircon, as a commonly used mineral for U–Th–Pb dating, widely exists in the continental crust and records a variety of geological activities. Due to the universality and stability of zircon and the long half-lives of U and Th isotopes, zircon U–Th–Pb geochronology can provide nearly continuous records for almost the entirety of Earth's history and is thus essential to studying the growth and evolution of the continental crust and Earth system evolution. Here, we present a database of zircon U–Th–Pb geochronology that samples the global continental crust and spans nearly all of Earth's history. This database collects ∼2 000 000 geochronology records from ∼12 000 papers and dissertations. This paper describes the compiled raw data, presents the relationship between dating error and zircon age, compares the error levels of different dating instruments, and discusses the impact of sampling bias on data analysis as well as how to evaluate and minimize this impact. In addition, we provide an overview of the temporal and spatial distribution of global zircon ages and provide key insights into the potential research value of zircon ages for Earth system science, such as crustal evolution, plate tectonics, and paleoclimate changes, as well as commercial use in mining and energy exploration. Overall, this data collection provides us with a comprehensive platform with which to study zircon chronological data in deep time and space. The described zircon database is freely available via Zenodo at https://doi.org/10.5281/zenodo.7387566 (Wu et al., 2023).

Late Neoarchean and Paleoproterozoic tectonothermal and metamorphic evolution recorded in high-pressure granulites: Evidence from the Xuanhua Complex, northern North China Craton

Archean cratons usually underwent reworked processes during the late Paleoproterozoic orogeny, which may result in multiple tectonothermal records preserved in some high-grade metamorphic rocks; thus, these rocks provide a key opportunity on understanding the tectonic styles during different periods. In this study, we report new petrological, mineral chemical and geochronological data of a suite of newly discovered mafic and felsic metamorphic rocks in the Xuanhua Complex from the northern North China Craton (NCC), revealing the complex metamorphic thermal history and the transition of tectonic regime from the late Neoarchean to the late Paleoproterozoic. The mafic granulite preserves the prograde quartz and clinopyroxene inclusions within garnets (M1), and the peak stage (M2) is recorded by garnet, clinopyroxene, hornblende, plagioclase, quartz and ilmenite. The first retrograde (M3a) and final cooling stages (M3b) are defined by the appearance of orthopyroxene, and the amphibolite-facies overprinting, respectively. The felsic granulite records the similar peak mineral assemblage (M2), but has higher volumes of feldspar and quartz than the mafic granulite. The retrograde stage (M3) is represented by amphibolite-facies assemblage (hornblende ± biotite ± clinopyroxene) and fine-grained intergrowth of hornblende + plagioclase around garnet. Pseudosection modeling and thermobarometry results suggest that both high-pressure (HP) mafic and felsic granulites record clockwise P–T paths, including isothermal decompression and cooling from the peak stage (M2: 10.4–10.8 kbar/830–860 °C and 9.6–11.2 kbar/850–930 °C, respectively) to the retrograde stages (M3a and M3b: 7.5–8.0 kbar/840–870 °C and 7.8–8.2 kbar/650–700 °C, respectively; M3: 7.2–8.6 kbar/670–800 °C). Zircon U–Pb dating shows that the magmatic precursors of these rocks were emplaced at c. 2.5–2.47 Ga. The metamorphic zircons reveal two groups metamorphic ages of c. 2.51–2.45 Ga and c. 1.90–1.85 Ga. The c. 2.51–2.45 Ga ages indicate a late Neoarchean tectonothermal event, and the c. 1.95–1.85 Ga ages are interpreted to represent the Paleoproterozoic HP granulite-facies metamorphism overprinting. The highest positive εHf(t) value of 7.64 and 2.8–2.6 Ga Hf model ages of c. 2.5 Ga magmatic zircons indicate that the c. 2.5 Ga mantle-derived mafic magma and the c. 2.5 Ga partial melting of 2.8–2.6 Ga basaltic juvenile crust generated the protoliths of mafic and felsic granulites, respectively, revealing 2.8–2.6 Ga and 2.5 Ga crustal growth and c. 2.5 Ga crustal reworking in the northern NCC.To sum up, we consider that the late Neoarchean magmatism and coeval metamorphism recorded in our samples are probably resulted from the amalgamation of micro-blocks. The Paleoproterozoic clockwise P–T–t paths with near-isothermal decompression segments and low geothermal gradients (c. 21–26 °C/km) of HP mafic and felsic granulites indicate a rapid post-collisional exhumation process after the significantly crustal thickening, and this long-term orogeny along the TNCO (c. 1.95–1.85 Ga; > 100 Myr) is related to the incorporation of the NCC into the Columbia supercontinent.

Petrogenesis of early Paleozoic granitoids in the North Qinling Orogen, Central China: Implications for crustal evolution in an accretionary orogen

Accretionary orogens are sites of extensive continental crustal growth and modification. The mechanism by which mafic crust is transformed into silicic melts (i.e., maturation of continental crust) is important for understanding the formation of the continental crust. The North Qinling Orogen (NQO) is a composite orogenic belt and contains an early Paleozoic accretion-dominated orogenic system, which is ideal for investigating continental crustal maturation. We obtained zircon and monazite U–Pb age and O isotope data, zircon Lu–Hf isotope data, and whole-rock major- and trace-element and Sr–Nd isotope data for early Paleozoic granitoids of the NQO. The granitoids are divided into three groups. Group 1 includes the Taiping tonalite (445 ± 3 Ma), the Manziying syenogranite (445 ± 2 Ma), and the Huoshenmiao granodiorite (436 ± 2 Ma). The Taiping and Huoshenmiao plutons have relatively high SiO2 contents (68.64–71.67 wt.%) and Na2O/K2O ratios (1.15–4.19), with enrichments in Rb, Ba, Th, and U and depletions in Nb, Ta, P, and Ti, and they are geochemically similar to sodic arc magmas. The Manziying syenogranite is a peraluminous potassic granite with high K2O contents (4.59–5.27 wt.%). Grantioids from Group 1 have similarly depleted Sr–Nd–Hf–O isotopic features (εHf[t] = +5.9 to +8.8; δ18O = 3.98‰–5.41‰), indicative of derivation via partial melting of oceanic arc crust, which suggests that partial melting of oceanic arc crust in a subduction system contributes to the generation of continental crust and causes its maturation. Group 2 consists of the Wuduoshan monzogranite (418 ± 2 Ma) and the Sikeshu granodiorite (423 ± 3 Ma). These plutons have relatively high SiO2 (65.59–72.06 wt.%), K2O (3.26–4.79 wt.%), and Al2O3 (14.65–16.12 wt.%) contents and Sr/Y (33–87) and (La/Yb)N (23–48) ratios. The Wuduoshan monzogranite has positive zircon εHf(t) (+0.4 to +3.1) and uniform δ18O (6.38‰–8.07‰) values, but the Sikeshu granodiorite has more variable isotopic compositions (εHf[t] = −1.9 to +5.0; δ18O = 6.37‰–10.60‰). The Wuduoshan monzogranite and Sikeshu granodiorite have similar whole-rock Sr–Nd isotopic compositions to basement rocks of the NQO. These features indicate that the two plutons formed by partial melting of basement rocks (i.e., subducted into the lower crust) of the North Qinling unit, along with juvenile crustal material. Group 3 is represented by the Xiaguan monzogranite, which formed at 434–430 Ma, and can be subdivided into heavy rare earth element (REE)-depleted and -enriched units. The former has high Sr/Y (56–98) and (La/Yb)N (34–73) ratios and low MgO (0.13–0.24 wt.%), Cr (0.37–1.69 ppm), and Ni (0.32–1.09 ppm) contents, similar to adakites derived from metabasaltic sources. The heavy REE-enriched nature of the Xiaguan monzogranite may reflect modification of its source by melt or fluid. Our results show that partial melting of enriched oceanic arc crust contributed to crustal maturation in an accretionary orogen. The addition of evolved crustal material also facilitated this process; therefore, the basement rocks and crustal thickness should be considered when assessing crustal dynamics in an accretionary orogen.

Zircon U-Pb-Hf isotope evidence for juvenile Penokean sources (ca. 1.9 Ga) in the SW Grenville Province

The Proterozoic evolution of Laurentia involved intermittent orogens and arc-related crustal growth along a modern-style convergent margin, termed the Great Proterozoic Accretionary Orogen (GPAO). The major lithotectonic provinces of the GPAO in the United States are exposed as telescoped equivalents within the thrust stack of the Grenville Province in Canada. We report zircon UPb and LuHf isotope analyses for two orthogneisses and one quartzite from Archean and Paleoproterozoic domains in the gneiss belt of the SW Grenville Province, Canada. The data constrains Paleoproterozoic crustal growth and sedimentation for this part of the Laurentian margin. The orthogneisses have similar zircon UPb and LuHf isotope systematics that reflect igneous reworking (ca. 1700–1660 Ma) of a Penokean aged (ca. 1.9 Ga) depleted mantle derived source, along with metamorphic reactivation at both ca. 1.45 Ga and during the Grenville Orogeny (1.2–1.0 Ga). The quartzite has a bimodal Paleoproterozoic (2.1–1.8 Ga) and Archean (> 2.5 Ga) detrital zircon age spectrum. The Paleoproterozoic detrital zircon grains have a wide range of Hf(t) compositions (εHf(t) = −12.5 to 10) reflecting co-eval interaction of Archean and depleted Paleoproterozoic sources, contrasted by a more restricted range for the Archean population (εHf(t) = −7 to 4). The new data support the presence of reworked ca. 1.9 Ga depleted mantle sources in the Paleoproterozoic Parautochthon, which extend further east than previously thought. They also identify correlations amongst detrital zircon age spectra and Hf isotope compositions in quartzites from different locations of the Archean Parautochthon.

Crustal growth and reworking at Archean plate margins

Crustal growth and reworking at Archean plate margins

Trondhjemites from the Western Dharwar Craton, Southern India: Implications for Mesoarchean crustal growth

Trondhjemite, together with tonalite and granodiorite (TTG) are important building blocks of the Archean continental crust and their studies provide insights into the formation of the early continents. The Western Dharwar Craton (WDC) in southern India contains voluminous TTG rocks that record prolonged Archean crustal evolution history with evidence for multiple episodes of magmatic accretion and metamorphic events from the Paleoarchean (∼3600 Ma) to Archean-Proterozoic transition (∼2500 Ma ago). In this study, we investigate trondhjemites from different typical locations in the WDC. Petrological, geochemical and zircon UPb, REE and LuHf data are presented and evaluated to understand the timing of magma emplacement and petrogenetic history. The trondhjemites in our study were generated from partial melting of juvenile hydrous basalt with arc features. The melting depth is variable, mainly shallow, while some melts originated at greater depths with variable plagioclase, amphibole with or without garnet, ilmenite or rutile in source residue. The age data presented here suggest that the crust accretion of the Mesoarchean Western Dharwar Craton occurred dominantly at 3.13–3.27 Ma followed by reworking of lower crustal basement 3.0–2.9 Ga.

An early Paleozoic accumulation-foundering cycle of ultramafic cumulates in the Harlik arc and its implications for continental crustal growth in the Altaids

Lower crustal ultramafic cumulates that are normally denser than the ambient mantle peridotites are prone to founder back into the convecting mantle, prompting the formation of continental Moho and andesitic crust at subduction zones. However, foundering events of ultramafic cumulates have not been widely reported in island arcs, including the Altaids, one of the largest accretionary orogens. The Altaids involved prominent continental crustal growth by the development of an archipelago in the Phanerozoic. But the scarcity of foundering events of ultramafic cumulates hinders a better understanding of the formation mechanism of andesitic crust in the Altaids. In this study, we report two dioritic plutons in the Harlik arc, southern Altaids, trying to identify the accumulation-foundering cycles of dense ultramafic cumulates. LA-ICP-MS dating revealed that the Miaoergou and Banfanggou plutons in the Harlik arc were intruded at ca. 445 ± 6 Ma and 423 ± 4 Ma, respectively. The homogeneous Miaoergou pluton is mainly composed of quartz diorites, they have a high- to medium-K calc-alkaline composition, and show enrichments in light rare earth elements (LREE) and large ion lithophile elements (LILE), but depletions in high field strength elements (HFSE, especially Nb and Ta) and heavy rare earth elements (HREE). They also have juvenile isotopic signatures (εNd(t) = +3.3 − +3.7, initial 87Sr/86Sr = 0.70399–0.70452), which together suggest that they were derived from fractional crystallization of the coeval hydrous basaltic magmas in an arc setting. The heterogenous Banfanggou pluton comprises a spectrum of lithologies from hornblende gabbro to diorite. The hornblende gabbros and diorites have a high- to medium-K calc-alkaline composition, and show enrichments in LREE and LILE, but depletions in HREE and HFSE. They all have juvenile isotopic compositions with εNd(t) ranging from +3.5 to +4.2, and (87Sr/86Sr)i from 0.70384 to 0.7047, respectively, indicating that their parental magmas were intruded in a subduction-related setting, and underwent fractional crystallization of clinopyroxene, plagioclase, amphibole and quartz. Both pulses of hydrous basaltic magmatism involved open-system conditions, with the development of mafic microgranular enclaves (MMEs) in the diorites and their mingling. In a synthesis of the published data, we found that the Harlik arc was continuously affected by subduction in the early Paleozoic, through extensive arc magmatism from ca. 452 to 440 Ma, the Harlik arc matured with the accumulation of thick ultramafic cumulates in its lower crust. The foundering of these ultramafic cumulates occurred later, which led to strong extension and related A-type granitic magmatism at ca. 436 Ma. This accumulation-foundering cycle of ultramafic cumulates facilitated the formation of andesitic composition and continental Moho in the Harlik arc, and contributed to the significant continental crustal growth in the Altaids.

Petrogenesis of Early Cretaceous andesites and mafic dikes in central Tibet: Implications for the growth of continental crust in collision zones

Post-collisional magmatic rocks in central Tibet provide an opportunity to study the processes of continental crustal growth in a collision zone. Here, we report data for samples of andesites and coeval mafic dikes obtained from the Jiaqiong area of the Bangong–Nujiang suture zone, central Tibet. Zircon U–Pb dating shows that the andesites and mafic dikes were emplaced concurrently during the late Early Cretaceous (ca. 115 Ma). The andesites have high Mg# values (55.2–61.6) and initial 87Sr/86Sr ratios (0.7086–0.7088), negative whole-rock εNd(t) (−6.5 to − 6.1) and zircon εHf(t) (−5.1 to − 1.7) values, and elevated zircon δ18O (6.8‰–7.5‰) values, indicating that the andesites were probably derived by the partial melting of an ancient subduction-modified lithospheric mantle. The coeval mafic dikes have moderately high Mg# values (49.3–50.0) and initial 87Sr/86Sr ratios (0.7076–0.7082), and negative εNd(t) values (−5.4 to − 5.3), implying that they were derived from an enriched lithospheric mantle source. However, the zircon εHf(t) (−5.8 to + 1.6) and δ18O (5.2‰–6.3‰) values of these mafic dikes indicate the influence of depleted mantle components. Furthermore, results of Sr–Nd–Hf–O isotopic mixing modeling suggest that these mafic dikes could have been produced by interaction between asthenospheric and lithospheric mantle. Therefore, we suggest that the parental magmas of the mafic dikes were probably generated by interaction between asthenospheric and overlying metasomatized lithospheric mantle. Combined with data on nearby Cretaceous magmatic rocks and sedimentary strata, we propose that the studied magmatic rocks were formed in a post-collisional slab breakoff setting in central Tibet and have made a substantial contribution to the vertical growth of continental crust in this continental collision zone.

A brief introduction to the detrital zircon U–Pb and Hf isotopic datasets for mainland China and adjacent regions

AbstractImprovements in analytical technique and data processing lead to a large deposit of detrital zircon U–Pb and Hf isotopic data in published journal articles, thesis and reports. With the support of Deep‐time Digital Earth Big Science Program, OneSediment Project conducted a team‐working effort to assemble the detrital zircon U–Pb and Hf isotopic data for mainland China and adjacent regions. Entity‐relationship diagram was constructed for detrital zircon data to determine what information need to be collected. Researchers and students from 13 institutions were organized together to fulfil this task and generated 13 detrital zircon datasets. The assembled datasets totally include 6635 samples and 560,596 analyses. The samples were checked in terms of lithology and stratigraphy according to the latest stratigraphic framework. Collected detrital zircon U–Pb and Hf isotopic data were recalculated for the best U–Pb age and εHf(t) values. They can be used to track crustal growth, tectonic evolution, provenance analysis and palaeogeographical change in regional and, when combined with other dataset, global scales.

The Singhbhum Craton (India) records a billion year of continental crust formation and modification

The petrogenesis of continental crust from its ultimate mantle source can be reconstructed from the element abundances and radiogenic isotope compositions of ideally pristine igneous rocks. The initial isotope compositions of igneous rocks provide geochemical constraints on the age, composition and evolution of their source(s). Determining initial isotope ratios for rock samples can be challenging, especially in rocks with a long and protracted thermal history. The Rb-Sr system is highly sensitive to parent-daughter element fractionation during magma differentiation. This makes the Rb-Sr isotope systematics ideal to trace the precursor composition of Archean felsic crust and constrain the time of element fractionation during the formation and subsequent modification of continental crust. Initial isotope compositions can be obtained directly from minerals that strongly prefer the daughter element and effectively exclude the parent element of the radio-isotope system of interest. Apatite, having a near zero Rb/Sr ratio, is ideal for preserving its initial 87Sr/86Sr and zircon records initial 176Hf/177Hf compositions. Combined modelling of Sr and Hf isotope data from granitoids of the Archean Singhbhum Craton, indicates that the older Paleoarchean granitoids, emplaced between 3.53 Ga and 3.44 Ga, were derived from a mafic precursor (∼52–54 wt% SiO2) sourced from a depleted mantle at ∼3.71 Ga. Initial 87Sr/86Sr isotope signatures of matrix apatite and apatite inclusions in zircon from the younger Paleoarchean granitoids (3.4–3.2 Ga) of the Singhbhum Craton indicate these younger granitoids were produced by mixing of magma generated from an older mafic source and partial melts derived from the older granitoids. The combined Sr-Hf isotope modelling links the timing of mantle extraction of the precursor material for Paleoarchean Singhbhum granitoids with a known mafic crust extraction event at ∼3.71 Ga. In combination, the new Sr isotope data from apatite combined with whole rock and zircon Hf isotope data from the literature reveal a ∼1 Ga protracted crustal growth and differentiation history of the nucleus of the Singhbhum Craton. By combining radio-isotope systems like 87Rb-87Sr and 176Lu-176Hf, the petrogenesis of Archean felsic crust from the extraction of mafic material from the mantle to reworking in an orogenic cycle to emplacement can be reconstructed. This approach can be applied to other greenstone-gneiss terranes to quantify the spatio-temporal and compositional evolution of voluminous felsic crust and the formation of cratons in the Archean.

Limited flexural control of fold-thrust belts on the Jurassic Sichuan Basin, South China

The northern part of the Jurassic Sichuan Basin has long been thought of as a foreland basin in relation to the post-collisional compression along the northern margin of the Yangtze block. However, the exact coupling mechanism between mountain building and basin formation remains unclear. Here, we integrate stratigraphic correlation, basin subsidence analysis and flexural simulation to quantitatively assess the extent to which the fold-thrust belts have controlled basin subsidence. Flexural backstripping of the stratigraphic record, spanning from 201 to 149 Ma, along two cross sections that are perpendicular to the Micangshan fold-thrust belt and the Dabashan fold-thrust belt, respectively, reveals a limited flexural control of mountain loading on basin subsidence. Owing to the short-wavelength nature of plate flexure, the basin-adjacent thrust belts exerted dominant control on basin subsidence only along its margin, with the width of the foredeeps not exceeding ca. 120 km, failing to drive subsidence in the forebulge and backbulge regions. Flexural modeling results suggest that crustal thickening was relatively weak during the Early to Middle Jurassic. This was followed by a more rapid and intense phase of crustal growth in the Late Jurassic, possibly extending into the earliest Early Cretaceous. Compared to the Micangshan region, the Dabashan region has experienced more intense compression during the Late Jurassic. Additionally, our results reveal spatial variations in plate rigidity along the northern margin of the Yangtze block, with greater plate stiffness in the Dabashan region. The presence of residual subsidence, an anomalous long-wavelength subsidence component corrected for both basin-adjacent thrust loading and associated sediment loading, highlights the necessity for an additional driving mechanism for basin subsidence. This residual subsidence was likely dynamic subsidence induced by the flat subduction of the Paleo-Pacific plate (the Izanagi plate) beneath East Asia as the flat slab progressively migrated inland.

Petrogenesis of two episodes of Neoarchean TTG gneisses (∼2.72 Ga and 2.56–2.54 Ga) from the Zhongtiao Mountains, North China Craton: Implications for crustal evolution

The petrogenesis of tonalite-trondhjemite-granodiorite (TTG) suites can provide key clues for the formation and evolution of continental crust. Neoarchean TTG gneisses are widely distributed in the Zhongtiao Mountains (ZTM), North China Craton (NCC). In this contribution, to better understand Neoarchean crustal evolution of the NCC, two discrete episodes of TTG gneisses (∼2.72 Ga and ∼2.56–2.54 Ga) from the ZTM are investigated. New results approve that the two episodes of TTG gneisses are all characterized by high SiO2, Na2O contents, (La/Yb)N and Sr/Y values, but low K2O/Na2O ratios. These geochemical features are similar to the Archean TTG suites and Phanerozoic high-SiO2 adakites, suggesting that they were both derived from partial melting of mafic rocks. The characteristics of trace elements (e.g., La/Yb, Sr/Y, Y/Yb, (Ho/Yb)N and Nb/Ta ratios) imply that both of them correspond to MP-TTGs, and garnet-amphibole exist in the residue phase. The ∼2.72 Ga TTGs show high Mg# values (ave. = 51), Cr (4.87–104 ppm) and Ni (6.37–55.7 ppm) contents, and negative Nb, Ta anomalies. Some high Nb/Ta ratios (up to 17.6) imply that rutile might be involved in the residual phase, which should require a higher pressure. Zircons from these rocks yield positive εHf(t) values of +4.0 ∼ +6.8 (ave. = +5.6), with two-stage Hf model ages (TDM2) of ∼2.85–2.70 Ga, indicating a juvenile basaltic crustal source. Hence, we interpret that ∼2.72 Ga TTGs in the ZTM were derived from partial melting of mafic rocks, and more likely related to subduction process. In contrast, the ∼2.56–2.54 Ga TTGs have lower Mg# values (39–54, mostly <45), Cr (6.02–22.8 ppm) and Ni (4.93–15.7 ppm) contents. They exhibit more evolved zircon Hf (average εHf(t) = +3.21) and whole-rock Nd (εNd(t) = +1.87) isotopic compositions, and their TDM2 ages are concentrated at ∼2.85–2.68 Ga. These ∼2.5 Ga TTGs in the ZTM show similar geochemical features with contemporary TTGs in the XQL, DF and EZT of the southern NCC, identical to melts from partial melting of mafic rocks under a thickened lower crust. The ∼2.5 Ga TTGs, A-type granites and mafic rocks in the ZTM jointly indicate a complicated crust-mantle process at different depths related to mantle plume activity during the late Neoarchean. Based on comparison of the spatial-temporal distribution and geochemical features of ∼2.7 Ga and ∼2.5 Ga TTGs in the NCC, we propose that Neoarchean crustal evolution of the ZTM is similar to the whole NCC, which recorded a significant crustal growth event at ∼2.7 Ga, and an important crustal reworking event associated with some juvenile crustal growth at ∼2.5 Ga.

Geochemical evidence for evolving Proterozoic crustal thickness and orogenic styles in southwestern Laurentia

It has long been challenging for researchers to track the crustal thickness and mode(s) of crustal modification in ancient convergent margins, limiting evaluation of the tectonic styles and processes that modify continental crust during orogenesis. We present trace element igneous geochemical crustal thickness proxies that quantitatively track the crustal thickness evolution of the long-lived Proterozoic active margin in the southwestern U.S.A. We integrate these results with geobarometric data to constrain the mode of crustal modification. The data indicate a complex record of crustal thickness change in space and time and evolving orogenic styles. Geochemical proxies at 1.84–1.72 Ga are consistent with 20–40 km thick magmatic arcs that were locally thickened to ∼50 km during ∼1.75 Ga tectonism. During the Yavapai orogeny, 1.72–1.69 Ga, a ∼200-km-wide belt of 50-60 km thick crust extended from southern California to northern Colorado and was rapidly thinned and exhumed by ∼1.68 Ga. Crustal thickening and thinning during the Yavapai orogeny largely occurred by shortening and exhumation, respectively, in the upper 25 km of the crust. Subsequent 1.68–1.60 Ga tectonism involved crustal growth, local crustal thickening, and low-P, high-T metamorphism, consistent with extensional accretionary orogenesis. The 1.47–1.37 Ga Picuris orogeny was associated with 50–60 km thick crust across much of the Southwest and involved crustal shortening with ∼10 km of magmatic underplating. Advective heat from the emplacement of ferroan granites in the mid-crust likely contributed to elevated geothermal gradients and rheologically weakened the crust. Our results suggest evolving orogenic styles in the Southwest from 1.75–1.69 Ga short-lived crustal thickening associated with terrane accretion to 1.69–1.60 Ga largely extensional accretionary orogenesis, and regional, long-lived crustal thickening at 1.47–1.37 Ga involving extensive basaltic underplating. Contrasting with some recent hypotheses, our data document a complex middle Proterozoic record for the Southwest that was not orogenically quiescent or tectonically stagnant but involved complex mountain building styles.

The Influence of Undercooling and Sector Zoning on Clinopyroxene–Melt Equilibrium and Thermobarometry

Abstract Thermobarometry provides a critical means of assessing locations of magma storage and dynamics in the lead-up to volcanic eruptions and crustal growth. A common approach is to utilise minerals that have compositions sensitive to changes in pressure and/or temperature, such as clinopyroxene, which is ubiquitous in mafic to intermediate magmas. However, clinopyroxene thermobarometry may carry significant uncertainty and require an appropriate equilibrium melt composition. In addition, the degree of magma undercooling (ΔT) affects clinopyroxene composition and zoning, with common sector zoning potentially obfuscating thermobarometry results. Here, we use a set of crystallisation experiments on a primitive trachybasalt from Mt. Etna (Italy) at ΔT = 25–233 °C, P = 400–800 MPa, H2O = 0–4 wt % and fO2 = NNO + 2, with clinopyroxene crystals defined by Al-rich zones (prisms and skeletons) and Al-poor zones (hourglass and overgrowths) to assess common equilibrium models and thermobarometric approaches. Under the studied conditions, our data suggest that the commonly applied Fe–Mg exchange (cpx-meltKdFe–Mg) is insensitive to increasing ΔT and may not be a reliable indicator of equilibrium. The combined use of DiHd (CaMgSi2O6 + CaFeSi2O6) and EnFs (Mg2Si2O6 + Fe2Si2O6) models indicate the attainment of equilibrium in both Al-rich and Al-poor zones for almost all investigated ΔT. In contrast, CaTs (CaAl2SiO6) and CaTi (CaTiAl2O6) models reveal substantial deviations from equilibrium with increasing ΔT, particularly in Al-rich zones. We postulate that this reflects slower diffusion of Al and Ti in the melt compared with Ca and Mg and recommend the concurrent application of these four models to evaluate equilibrium between clinopyroxene and melt, particularly for sector-zoned crystals. Thermobarometers calibrated with only isothermal–isobaric experiments closely reproduce experimental P–T at low ΔT, equivalent to natural phenocrysts cores and sector-zoned mantles. Models that also consider decompression experiments are most accurate at high ΔT and are therefore suitable for outermost phenocryst rims and groundmass microlites. Recent machine learning approaches reproduce P–T conditions across all ΔT conditions. Applying our experimental constraints to sector-zoned microphenocrysts and groundmass microlites erupted during the 1974 eccentric eruption at Mt. Etna, we highlight that both hourglass and prism sectors are suitable for thermobarometry, given that equilibrium is sufficiently tested for. The combination of DiHd, EnFs, CaTs and CaTi models identifies compositions closest to equilibrium with the bulk melt composition, and results in smaller differences in P–T calculated for hourglass and prism sectors compared with applying only DiHd and EnFs equilibrium models. This provides a framework to assess crystallisation conditions recorded by sector-zoned clinopyroxene crystals in mafic alkaline settings.

Archean to Paleoproterozoic magmatic and metamorphic evolution of the Lofoten – Vesterålen area in northern Norway through zircon U-Pb and Hf isotopic analysis

Neoarchean (2.8–2.5 Ga) and Paleoproterozoic (2.5–1.6 Ga) crust preserves a record of major crustal growth events, the oldest supercontinent assemblies, and the likely transition to modern plate tectonics. The Lofoten–Vesterålen area in Norway underwent a multistage tectono-magmatic history from the Neoarchean to the Paleoproterozoic. In this region three main components are distinguished: 1) a granulite facies zone (GFZ); 2) an amphibolite facies zone (AFZ); and 3) an anorthosite-mangerite-charnockite-granite (AMCG) suite. New zircon U-Pb and Lu-Hf isotope data are used to investigate the evolution of this 2.7–2.6 and 1.9–1.8 Ga crust and evaluate how the developments fit in the global context of Earth’s crustal evolution. In-situ LASS–MC-ICPMS analyses were conducted on zircons to simultaneously determine U-Pb ages and Hf isotope compositions. Solution Hf isotope compositions were additionally obtained on zircons previously dated by ID-TIMS. In general, there is a good agreement between solution and in-situ data. Samples from the GFZ and AFZ both yield Neoarchean (2.7–2.6 Ga) and Paleoproterozoic (1.9–1.8 Ga) populations. Neoarchean populations have εHf(i) of +0.9 to −0.3; Paleoproterozoic populations have εHf(i) around −16 and −7 at ∼1.80 and ∼1.87 Ga, respectively. The AMCG suite was emplaced in two magmatic pulses: the early one at 1.87–1.86 Ga characterized by εHf(∼1.87 Ga) ∼ −7, and the main pulse at 1.80–1.79 Ga has εHf(∼1.80 Ga) between −4.0 and −15.0. Metamorphic overprinting recorded within the GFZ and the AMCG at ∼1.80 Ga is absent in the AFZ, supporting the hypothesis of a tectonic boundary between these domains. The Hf isotope compositions suggest that Paleoproterozoic magmatism was sourced from a combination of juvenile material and incorporated ancient crust, probably in variable proportions, as suggested by the variation in the isotope composition of the AMCG suite. Broadly chondritic Hf isotope compositions in Neoarchean rocks indicate the absence of significantly older crust as a source component. Main crustal growth events recorded in the Lofoten–Vesterålen area coincide with global age peaks and zircon Hf isotope record. The AMCG suite was emplaced during the Columbia supercontinent assembly. Neoarchean and Paleoproterozoic crust in the Lofoten–Vesterålen area exemplifies the elevated preservation potential in collisional tectonic settings.

Cryogenian-Ediacaran crustal growth and evolution of the active margin of the Dahomeyide belt, Ghana

AbstractThe study presents detailed petrographical, geophysical, structural and geochemical data of the internal nappes zone to establish the deformational history, origin and tectonic setting and constrain the crustal growth and evolution of the active margin of the Dahomeyide belt. Two main lithological units, (i) deformed meta-granitoids (migmatites and gneisses) and (ii) undeformed granitoids, dominate the internal nappes zone. The granitoids are generally I-type, metaluminous to weakly peraluminous, low-K tholeiite to high-K calc-alkaline and of tonalite, granodiorite and granite affinity. The overall trace element patterns of the studied granitoids characterized by the enriched LILE and depleted HFS, with negative peaks of Nb-Ta, Sr, P and Ti, are indications of arc-related magmatism. Structural analysis reveals four deformation phases (D1-D4). D1 represents Northwest-Southeast (NW-SE) Pan African shortening associated with a continent-continent collision, resulting in westward nappe stacking. Progressive NW-SE shortening resulted in D2 and D3 top-to-the-NW dextral and sinistral thrusting events during the Pan-African orogeny. D4 is an extensional event likely associated with the orogenic collapse phase. The gneisses and migmatites, with dominant axial planar foliations, point to their formation in a collisional setting or influence by the Pan-African collisional processes. Continental-arc signatures in these rocks imply continental subduction during their protolith formation. The intrusive granitoid and pegmatite are undeformed, meaning late- to post-orogenic emplacement. These findings suggest that the internal nappes zone archived the subduction-collision and post-collisional phase of the Pan-African orogeny and recorded large-scale migmatization and granitoid emplacement due to partial melting of thickened lower crust between Mid-Cryogenian and late Ediacaran.

Geochemistry and zircon U-Pb-Hf isotopes of Paleozoic granitoids along the Solonker suture zone in Inner Mongolia, China: Constraints on bidirectional subduction and closure timing of the Paleo-Asian Ocean

The closure of the Paleo-Asian Ocean resulted in formation of the Central Asian Orogenic Belt (CAOB) accompanied by considerable Phanerozoic juvenile crustal growth. As an international frontier topic, whether the Solonker suture zone is the site of final closure of the Paleo-Asian Ocean and records the terminal evolution of the CAOB in Inner Mongolia is still in controversial. This study presents the geochemistry and zircon U-Pb-Hf isotopes of granitoids along the Solonker suture zone to constrain the subduction and closure processes of the Paleo-Asian Ocean. Granitoids from the Xilinhot area in the northern accretionary zone comprise diorites and granodiorites with LA-ICPMS zircon U-Pb ages of 316 ± 2 Ma, 302 ± 2 Ma, 285 ± 2 Ma and 243 ± 4 Ma, positive εHf(t) values of +2.39 to +13.31 and TDM2 model ages of 481–1191 Ma. The ca. 316–285 Ma granitoids were formed by fractional crystallization of regional coeval calc-alkaline basaltic magma in an arc-related setting and the ca. 243 Ma granitoid was generated mainly from partial melting of dominant basaltic source rocks and minor metagraywackes in the northward subduction-related setting. The granitoids from the Linxi area in the southern accretionary zone include muscovite monzogranites and monzogranites with LA-ICPMS zircon U-Pb ages of 525 ± 6 Ma, 274 ± 2 Ma and 261 ± 4 Ma. The ca. 525 Ma muscovite monzogranites have negative zircon εHf(t) values of −19.36 to −13.61 with TDM2 model ages of 2348–2711 Ma. They were generated from partial melting of Neoarchean to Paleoproterozoic supracrustal materials. The ca. 274–261 Ma monzogranites show negative zircon εHf(t) values ranging from −10.48 to −0.66 with TDM2 model ages of 1335–1942 Ma. They were derived from remelting of Paleoproterozoic to Mesoproterozoic supracrustal materials in the southward subduction zone along the northern active margin of the North China Craton. Combining our results with existing data, we infer that the early Paleozoic southward subduction of the Paleo-Asian Ocean started before ca. 525 Ma, and the northward subduction zone lasted to ca. 243 Ma, which implies the bidirectional subduction of the ocean lasted at least until ca. 243 Ma. Our work suggests that along the Solonker suture zone, the regional amalgamation may have occurred after ca. 243 Ma, which closed the Paleo-Asian Ocean.

Tracing the crustal evolution of the Precambrian Southern Granulite terrane in East Gondwana: New insights from zircon U-Pb/Hf geochronology

The Precambrian Southern Granulite terrane of south India has a crustal evolution history broadly bracketed between the late Archean and Cambrian with records of polyphase deformation, metamorphism, and magmatism. The Southern Granulite terrane comprises distinct crustal blocks bounded by shear/suture zones that have been variably correlated with supercontinent fragments including Madagascar, Sri Lanka, Africa, Eastern Ghats, and Antarctica. However, the timing and mechanism of assembly of different crustal blocks within the Southern Granulite terrane and its linkages with counterparts in East Gondwana are highly debated. This study aimed to unravel the complex crustal evolutionary pattern of the terrane by generating robust zircon U-Pb/Hf isotopic data from basement charnockites, gneisses, granitoids, and alkaline intrusive units from the central part of Southern Granulite terrane and comparing these results with similar data from different East Gondwanan terranes. The study identified four distinct crustal growth episodes in the Madurai block: (1) Neoarchean−early Paleoproterozoic, (2) Rhyacian−Orosirian, (3) late Tonian, and (4) Ediacaran−Cambrian. Analysis of zircon Hf isotope data revealed that the first two events are marked by juvenile magmatic signatures, whereas the latter two are distinctly associated with intense reworking and remelting of older crust with no significant juvenile input. Our new results combined with existing data from other Gondwanan terranes suggest a common Paleoproterozoic ancestry for the Southern Granulite terrane and its corresponding Gondwanan fragments, proposing a revision to the existing geodynamic models.