Crustal Growth Models Research Articles

Plate tectonics through Earth’s history: constraints from the thermal evolution of Earth’s upper mantle

ABSTRACT Temporal changes in Earth’s tectonic style play a crucial role in understanding the planet’s evolution. Modern-style tectonics is characterized by the formation of basaltic crust at divergent plate boundaries and its subsequent recycling at subduction zones, accompanied by wedge mantle formation and arc magmatism. It is commonly believed that the secular cooling of the mantle modified the tectonic style from stagnant lid or heat pipe on an early hotter Earth to horizontal tectonics during Meso-Neoarchean. However, various field, petrographic, and geochemical studies suggest that the onset of plate tectonics ranges from Hadean to Neoproterozoic. In this study, we re-evaluate the primary magma temperature (mantle potential temperature, Tp) of the upper ambient mantle, spanning from Eoarchean to Neoproterozoic. We used basalts from several Archean and Proterozoic greenstone belts worldwide, along with Proterozoic ophiolites, to (re)calculate the Tp using the FRACTIONATE-PT method. We observed a Tp range of 1444–1639°C during the Archean and 1414–1611°C during the Proterozoic. These findings indicate a strong correlation with previously estimated Tp values obtained from PRIMELT3 method. This further highlights strong internal consistency among different methods and supports models of a hot ambient mantle during the Archean and Proterozoic. We further reviewed numerical models regarding the effect of mantle temperature on the viability of early Earth plate tectonics. Such model results, composition of Archean continental crust, recent crustal growth models, and field evidence are consistent with the operation of plate tectonics on an early hotter Earth. The transition from a hotter mantle to a colder one from Eo-Neoarchean resulted in thinner oceanic lithosphere and a less depleted lithospheric mantle. Subduction of this thinner oceanic lithosphere led to modern-style tectonics. The intensity and style of plate tectonics on the early Earth differed from modern tectonics, which emerged during the Neoarchean.

Andesitic arc magmas derived from two contrasting mélange origins: Evidence from central Tibetan dioritic porphyries

Numerous studies have argued that the andesitic model (i.e., the arc directly produces andesitic magmas) rather than the basaltic-input model (i.e., andesitic magmas were derived by fractionation from basaltic magmas) better explains the composition of continental crust. However, it remains controversial as to how andesitic magmas are directly produced in subduction zones, and this weakens the basis of the andesitic model. Recently, the mélange-diapir model has been increasingly considered to be an important mechanism for the direct generation of arc andesites, and would thus support the andesitic model for crustal growth. Moreover, melting experiments on mélange rocks show that varying the composition of mélange rocks (e.g., serpentinite-dominated and sediment-dominated mélange) can produce a range of primary andesitic arc magmas, from tholeiitic, calc-alkaline, to high-K calc-alkaline and shoshonitic magma types. However, it remains challenging to determine the different mélange sources for low-K to high-K arc andesites. Here we verify these experimental results using data from early Cretaceous dioritic porphyries in the southern Qiangtang block of central Tibetan plateau. The generation of these dioritic rocks cannot be linked to any coeval basaltic rocks in this region and so they may represent primary andesitic magmas. These dioritic porphyries are divided into low-K type-1 and high-K type-2. Type-1 has higher δ26Mg values (−0.16 to −0.08‰) than MORBs (mid-ocean ridge basalt, δ26Mg = −0.25 ± 0.06‰), which can be ascribed to the contribution of subducted bulk/un-differentiated serpentinite in their source. Combined with their low K contents and La/Sm ratios, and high NdHf isotope ratios (εNd(t) = 3.36–5.03; εHf(t) = 12.88–13.76), we suggest a serpentine-dominated mélange for their origin. The type-2 and type-1 rocks have similar trace-element distribution patterns, the same ages (∼124 Ma), and are indistinguishable in field outcrops, indicating their common petrogenesis. However, compared to type-1, type-2 has higher K and Th contents, and La/Sm and Th/Nd ratios with lower NdHf isotope values (εNd(t) = −1.61 to +0.30; εHf(t) = 3.78–4.39), and mantle-like δ26Mg values (−0.28 to −0.17‰). This likely indicates derivation from a sediment-dominated mélange for type-2. This study thus verifies the experimental results of melting of different mélange from natural rock record and shows that the mélange model provides an important mechanism for generating variable andesitic arc magmas and so the formation of the andesitic continental crust.

Persistent mildly supra-chondritic initial Hf in the Lewisian Complex, NW Scotland: Implications for Neoarchean crust-mantle differentiation

The composition of the mantle source from which continental crust is ultimately derived informs both the timing and rate of continental growth. This composition is, however, poorly constrained during the first billion or so years of Earth's evolution. Typically in zircon Hf isotope studies, a linearly evolving depleted mantle back-calculated from modern MORB-source mantle to chondrite at 4.57 Ga has been used to calculate model ‘crust formation’ ages, which are themselves used in crustal growth models. Yet, an increasing number of studies of Eo- to Paleoarchean continental crust suggest crust extraction from a relatively undifferentiated mantle. Zircon Hf isotope studies of ca. 2.9 Ga gneisses from the Neoarchean Lewisian Gneiss Complex of Scotland have previously revealed initial Hf isotope compositions that suggest derivation from such an undepleted mantle source. Here, we present combined zircon U-Pb and Lu-Hf isotope data from seventeen representative grey gneiss samples from across the Lewisian of the Scottish mainland and the Outer Hebridean archipelago, ranging in age from 3.1 to 2.7 Ga. Regardless of location and exact age, initial Hf isotope signatures are predominantly near- to slightly supra-chondritic. Several models are investigated to explain this observation, our favoured one involving episodic extraction of the grey gneiss precursors from a mildly depleted mantle that began to diverge from chondritic composition at ca. 3.5 Ga. Samples from the southern Outer Hebrides record slightly more radiogenic initial Hf signatures, consistent with a possible terrane boundary along the Paleoproterozoic South Harris shear zone. This study provides further confirmation that domains of Hadean to Paleoarchean, and possibly even locally Mesoarchean, mantle remained chondritic with respect to Hf isotopes. Initiation of a depleted mantle source at this later time and its limited divergence from chondritic mantle, even into the Neoarchean as suggested here, has implications for crustal growth models based on detrital zircon Hf model ages. On a regional basis, the signatures recorded in our data are consistent with contemporaneous gneisses across Fennoscandia and East Greenland and the mildly depleted mantle source of continental crust furthermore persists into the Mesoproterozoic.

A crustal growth model for the eastern Central Asian Orogenic Belt: Constraints from granitoids in the Songnen Massif and Duobaoshan terrane

The crustal growth and reworking processes of accretionary orogen such as the Central Asian Orogenic Belt (CAOB) have been a controversial issue. Here, we present in situ zircon U–Pb and Lu–Hf isotopic data for granitoids from a microcontinent (the Songnen Massif) and an arc terrane (the Duobaoshan arc terrane), which generally represent two main crustal components in the eastern CAOB to establish a crustal growth model for an accretionary orogen and to trace the influence of the assembly and breakup of supercontinents on crustal evolution in an accretionary orogen. Neoproterozoic–Mesozoic granitoids distributed in the Songnen Massif show a step-like crustal growth pattern over time with three major periods of development: Paleoproterozoic crustal growth at 2.2–1.8 Ga, Mesoproterozoic growth at 1.6–1.0 Ga, and a third pulse of growth at 0.85–0.6 Ga, along with two short pauses at 1.8–1.6 Ga and 1.0–0.85 Ga. The assembly and collision phases of supercontinents corresponded to the enhanced and degressive crustal growth rates of the Songnen Massif, respectively, suggesting that the supercontinent cycles are responsible for the episodic crustal growth pattern in the region. Crustal reworking of the Songnen Massif occurred during 1000–180 Ma, with a fluctuation at 800–600 Ma, which provided a major contribution to the isotopic heterogeneity of lower continental crust. Isotopic compositions of granitoids from the Duobaoshan arc terrane located between the Xing'an and Songnen massifs, together with those of granitoids from other microcontinents, suggest that most of the continental crust beneath the microcontinents in the eastern CAOB generated during the Precambrian, whereas a significant amount of lateral crustal accretion occurred in continental arc settings during orogenies and amalgamation of microcontinents during the Phanerozoic.

The diversity of sources of late Archean granites reflects a transition from plume-dominated to plate tectonics in the Superior Province, Canada

The Superior Province in Canada represents an excellent natural laboratory to test various models of crustal growth and evaluate tectonic regimes during the late Archean. Here, we present new geochemical and geochronological data for ca. 2.8–2.7 Ga plutons intruding different lithotectonic units of the north-eastern part of the Superior Province that demonstrate a complex nature and evolution of magma sources in relation to a changing tectonic setting. Four distinct plutonic suites were newly identified: (1) sodic tonalite–diorite (TD) with a composition resembling low-pressure TTG-like melts, (2) sodic tonalite–granodiorite–diorite (TGD) with medium/high-pressure TTG-like signatures, (3) Mg–K-rich monzogranite to monzodiorite (sanukitoids; MMD), and (4) K-rich granodiorite–granite–monzogranite (GGM). These suites are interpreted as recording a temporal evolution from plume-assisted melting of lower mafic continental crust through melting of a subducted oceanic slab at different depths to large-scale re-melting of the previously formed and amalgamated crustal units. In this evolutionary scheme, the ∼2730–2700 Ma subduction-related plutonism is the most voluminous addition of the juvenile arc material and thus represents the most significant crustal growth event. Altogether, the above inferences and published data suggest that the late Archean plutonism in the Superior Province formed during a transitional regime from plume-dominated to plate-tectonic over a short time span between ∼2730 and ∼2700 Ma and that this transition was markedly diachronous across the Province.

Identifying Crystal Accumulation in Granitoids through Amphibole Composition and In Situ Zircon O Isotopes in North Qilian Orogen

Abstract Granitoids are the main constituents of the continental crust, and an understanding of their petrogenesis is key to the origin and evolution of continents. Whether crystal fractionation is the dominant way to generate evolved magmas has long been debated, mostly because such processes would produce large volumes of complementary cumulates, which remains elusive. Mafic magmatic enclaves (MMEs) are ubiquitous in granitoids and their presence was initially recognized as cumulates. However, because many MMEs lack obvious evidence of accumulation, such as the classic cumulate textures and modal layering, the cumulate origin of MMEs has been abandoned and the model of magma mixing between mafic and felsic magmas has become popular. In this study, we conduct a combined study of amphibole composition and in situ O isotopes in zircons on three suites of orogenic granitoids with MMEs from the North Qilian Orogenic Belt (NQOB). We find that the MMEs and their host granodiorites show overlapping zircon δ18O values, affirming that they share the same parental magmas. The amphibole compositions indicate that amphiboles from the MMEs are not in equilibrium with a melt whose composition was that of the bulk-rock. These new data, together with the published bulk-rock data, suggest that the MMEs in our study have clear cumulate signatures and are thus of cumulate origin. Our study provides evidence for crystal accumulation in granitoids in the NQOB. This new understanding calls for re-examination on the petrogenesis of some intermediate magmatic rocks (granitoid/andesite) in discussing models of continental crustal growth.

Hf-Nd isotopes from ultramafic and mafic rocks in the western Dharwar Craton, India, record early Archean mantle heterogeneity

Both primary mantle differentiation and early extraction of oceanic and continental crust can have resulted in a chemically heterogeneous Archean mantle. Evidence for heterogeneity of the Archean mantle is preserved in the correlated Hf and Nd isotopes in mafic (basaltic) magmatic rocks and their ‘decoupling’ in ultramafic rocks, especially in komatiites, from the Archean Dharwar Craton, India. The komatiites and komatiitic basalts from the western Dharwar Craton were emplaced at ~3.25 Ga and have large spread in initial radiogenic isotope compositions with initial ɛHf (ranging from +0.9 to +21.3) and initial ɛNd from +0.2 to +7.4 that do not correlate. In contrast, the ɛHfi and ɛNdi in contemporaneous mafic (ɛHfi = +5.9 ± 11.2; ɛNdi = +3.2 ± 7.6) and felsic (ɛHfi = +4.3 ± 4.6; ɛNdi = +2.7 ± 4.2) rocks in the same area correlate, as is typical of rocks that are the product of magmatic processes involving melting and magmatic differentiation. The high initial ɛHf and the lack of correlation (‘decoupling’) of ɛHf and ɛNd in the komatiites and komatiitic basalts indicate that the source had superchondritic Lu/Hf which is attributed to early mantle differentiation at great depth, possibly in the perovskite stability field. This finding indicates that the primordial mantle was mineralogically and chemically layered and this layering was the result of fractionation processes very early in Earth's history, probably during magma ocean solidification. During crystallization of a primordial magma ocean, fractionation of Mg-perovskite with minor Ca-perovskite component could have resulted in a melt with high Lu/Hf and nearly unfractionated Sm/Nd. The correlation of the initial Hf isotope composition of the komatiites and komatiitic basalts with their MgO content, and thus melting temperatures, indicates the presence of a highly fractionated component derived from the deep mantle. The coupling of ɛHf-ɛNd in mafic to felsic igneous rocks and the decoupling in contemporaneous ultramafic rocks is thus the result of mantle material originating from different depths involving different ancient fractionated reservoirs. The existence of a chemical and isotopically heterogeneous Archean mantle that supplied material to form the continental crust has to be accommodated in crustal growth and preservation models that commonly assume a nearly homogeneous and gradually depleting mantle reservoir as the source of the continental crust.

Linking Gold Systems to the Crust-Mantle Evolution of Archean Crust in Central Brazil

The Goiás Archean Block (GAB) in central Brazil is an important gold district that hosts several world-class orogenic gold deposits. A better comprehension of the crustal, tectono-magmatic, and metallogenic settings of the GAB is essential to accurately define its geological evolution, evaluate Archean crustal growth models, and target gold deposits. We present an overview of gold systems, regional whole-rock Sm-Nd analyses that have been used to constrain the geological evolution of the GAB, and augment this with new in situ zircon U-Pb and Hf-O isotope data. The orogenic gold deposits show variable host rocks, structural settings, hydrothermal alteration, and ore mineralogy, but they represent epigenetic deposits formed during the same regional hydrothermal event. The overprinting of metamorphic assemblages by ore mineralogy suggests the hydrothermal event is post-peak metamorphism. The metamorphic grade of the host rocks is predominantly greenschist, locally reaching amphibolite facies. Isotope-time trends support a Mesoarchean origin of the GAB, with ocean opening at 3000–2900 Ma, and reworking at 2800–2700 Ma. Crustal growth was dominated by subduction processes via in situ magmatic additions along lithospheric discontinuities and craton margins. This promoted a crustal architecture composed of young, juvenile intra-cratonic terranes and old, long-lived reworked crustal margins. This framework provided pathways for magmatism and fluids that drove the gold endowment of the GAB.

On the petrogenesis of Paleoarchean continental crust: U-Pb-Hf isotope and major-trace element constraints from the Bastar Craton, India

Combined geochronological and (isotope) geochemical investigations of tonalite-trondhjemite-granodiorite (TTG) suites, which dominate Archean cratons, provide constraints on the sources and petrogenetic processes that gave rise to the Earth's early continental crust. In situ U-Pb and bulk Lu-Hf analyses of single zircon grains from TTGs of the Paleoarchean Bastar Craton in central India date the emplacement of the igneous rock suite and trace the chemical evolution of its source rocks. Complementary whole rock Lu-Hf and major and trace element data constrain the petrogenesis of the TTGs. The rocks are variably enriched in fluid mobile elements, particularly K and Pb, and are characteristically depleted in heavy rare earth elements. Calculated initial Hf isotopic compositions of zircon and whole rocks cluster around the chondritic value at ca. 3.55 to 3.45 Ga. The high precision Hf-isotope data reveal a trend of increasing radiogenic 176Hf/177Hf with age, the slope of which implies a 176Lu/177Hf = 0.02, and thus derivation from a mafic protolith. Geochemical modeling indicates that the granitoids can be directly produced by melting of a hydrous basalt in the garnet stability field, followed by minor fractional crystallization. Variations in the modal abundance of garnet suggest the source rocks melted at different depths. Rocks with moderately to strongly depleted HREE patterns require >10% garnet in the residue, consistent with derivation from a garnet-amphibolite or garnet-pyroxenite source. Modeling of the Lu-Hf systematics further highlights the important role of melting as opposed to fractional crystallization for reproducing the primary chemical characteristics of TTGs.The global Hf-isotope record of Archean mantle-derived rocks is considerably more radiogenic than that of the bulk silicate Earth, providing unambiguous evidence that parts of the mantle were significantly depleted in incompatible elements since at least the early Archean. The differences between Hf isotopic compositions of the mafic (greenstone belts) and felsic (gneissic terranes) rocks can be accounted for by crustal residence times of several 100 million years for the protoliths of the felsic magmatic suites. Consequently, TTGs integrate the geologic history of their precursor(s), complicating their isotopic record. Therefore, the near-chondritic Hf isotopic composition of TTGs is not direct evidence for their derivation from a primitive mantle, but rather a consequence of their specific petrogenesis involving an aged, mafic precursor. This relationship introduces uncertainties into models of crustal growth and geodynamics that are based on the felsic record, especially where the whole rock context is absent. In contrast, high-precision Hf-isotope analyses of single zircon grains, combined with whole rock isotope and elemental constraints, can be used to reliably identify the sources and processes involved in the generation of Earth's oldest continental crust.

Petrogenesis of volcanic rocks of the Devonian – Carboniferous Dahalajunshan Formation, Western Tianshan: Implications for crustal growth in an accretionary orogen

Magmatism associated with plate subduction and collision plays an important role in crustal growth during the evolution of the Earth, but the detailed mechanisms regarding how magmatism contributes to the vertical growth of the crust remain unclear, in part due to controversies on the petrogenesis of intermediate-felsic rocks. In the Central Asian Orogenic Belt (CAOB), although it has been recognized that both pre-existing crustal component and juvenile mantle-derived melt played a role in the crustal evolution of the orogens, the mechanics of how they contributed to vertical crustal growth remain uncertain. This study presents new geochronological, major and trace element geochemical, and Nd–Hf isotopic data for the Late Paleozoic intermediate-felsic volcanic rocks of the northern Yili terrane, Western Tianshan Orogen, with the aims to constrain the petrogenesis of the andesitic magma and to offer new insights into the deep interaction processes between crustal- and mantle-derived magmas. According to zircon U–Pb dating, trace element and isotopic compositional data, these volcanic rocks can be divided into two groups. Group 1 rocks (365–375 Ma) have relatively low whole-rock εNd(t) values (−4.2 – −1.6) and display high contents of Th (11.8–20.4 ppm), U (2.7–5.1 ppm), Nb (14.1–36.7 ppm) and Ta (1.1–3.1 ppm) that suggest they were derived from the partial melting of pre-existing crust. Group 2 rocks (346–349 Ma) show REE and trace elemental patterns typical of subduction-related calc-alkaline volcanic rocks with relatively high whole-rock εNd(t) values (−1.89 – +2.4), and with variable zircon Hf isotopic values (εHf(t) = +1.9 – +13.5) that can be attributed a juvenile crust with variable input of mantle material. Hf–Nd isotopic decoupling of the group 2 rocks suggests a replenishment by more mafic magmas. Samples from the entire Yili terrane show an increasing trend in εNd(t) values with younger ages that suggest an increasing involvement of melted juvenile crustal rocks. Overall, these results suggest a crustal evolutionary scenario in which: 1) pre-existing crustal components were melted and ascended to the upper crust (group 1) during the southward subduction of the Junggar-North Tianshan ocean in the early stage; and, 2) a later stage in which newly input mantle material solidified to form new arc crust that was also melted by basaltic underplating and that was probably associated with the rollback of the Junggar oceanic slab. This two-stage crustal growth model seems to be an efficient process for vertical crustal growth during the formation of an accretionary orogen.

Petrogenesis of early Carboniferous bimodal-type volcanic rocks from the Junggar Basin (NW China) with implications for Phanerozoic crustal growth in Central Asian Orogenic Belt

The Central Asian Orogenic Belt (CAOB) is one of the world's largest and longest-lived accretionary orogens that has been regarded as an important area of crustal growth during the Phanerozoic, yet the growth mechanism is still under debate. A combined study of petrology, geochronology and geochemistry for drilling-sampled early Carboniferous volcanic rocks from the Junggar Basin (NW China), southwestern CAOB, aims to constrain their petrogenesis and implications for Phanerozoic crustal growth. The volcanics show a bimodal-type suite, dominated by basalts, basaltic andesites and subordinate rhyolites, and have consistent zircon U–Pb ages of 334–336 Ma. The mafic rocks exhibit arc-like trace-element distribution patterns characterized by enrichment of LILEs (e.g., Pb) and LREEs but depletion of HFSEs (e.g., Nb, Ta and Ti). They have juvenile isotopic signatures of zircon ɛHf(t) (+11.4 to +15.4) and whole-rock ɛNd(t) (+7.8 to +8.5), with initial (87Sr/86Sr)i of 0.7040–0.7054 and (206Pb/204Pb)i of 17.681–17.817. These geochemical features, together with their variable Ba/La (4.28–44.2) but low (Tb/Yb)N (1.24–1.85) and Sm/Yb (1.65–3.44) ratios, suggest that the mafic volcanics could be derived from a main spinel-bearing mantle source metasomatized by subduction-related fluids. In contrast, the felsic samples are geochemically equivalent to A2-type granite, with high contents of SiO2, Zr, Ga, and show strong depletions in Eu, Sr and Ti. They have zircon Hf (ɛHf(t) = +12.0 to +15.3) and whole-rock Sr-Nd-Pb isotopic compositions ((87Sr/86Sr)i = 0.7032–0.7038; ɛNd(t) = +7.8 to +8.1; (206Pb/204Pb)i = 17.973–18.004) similar to the coeval mafic rocks. These observations, along with the model calculations using MELTS, favor a derivation of felsic rocks from the mafic magma by fractional crystallization. We suggest that the early Carboniferous bimodal-type volcanism formed in a localized back-arc extension within an overall convergent setting, presumably triggered by slab roll-back of the subducting Paleo-Junggar Ocean plate. This bimodal-type magmatism recorded a significant vertical crustal growth event in the Junggar Basin during the Phanerozoic. We propose that oceanic subduction (flare-ups) not only leads to the lateral accretion of arcs, but also results in the vertical growth of mantle-derived materials and subsequent magma differentiation during slab roll-back (flare-lulls); this may considerably change our view of the models of crustal growth in CAOB.

Constraints on Archean crust formation from open system models of Earth evolution

Establishing the mode and rate of formation of the continental crust is crucial for quantifying mass exchange between Earth’s crust and mantle. The limited crustal rock record, particularly of early Archean rocks, has led to a variety of different models of continental growth. Here, we present an open-system model of silicate Earth evolution incorporating the Sm-Nd and Lu-Hf isotope systematics with the aim to constrain crustal growth during the Archean and its effect on the chemical and isotopic evolution of Earth’s crust-mantle system. Our model comprises four reservoirs: the bulk continental crust (CC), depleted upper mantle (UM), lower mantle (LM), and an isolated reservoir (IR) where recycled crust is stored transiently before being mixed with the LM. The changing abundance of isotope species in each reservoir is quantified using a series of first order linear differential equations that are solved numerically using the fourth order Runge–Kutta method at 1 Myr time steps for 4.56 Gyr (the age of the Earth). The model results show that only continuous and exponential crustal growth reproduces the present-day abundances and isotope ratios in the terrestrial reservoirs. Our preferred crustal growth model suggests that the mass of the CC by the end of Hadean (4.0 Ga) and end of Archean (2.5 Ga) was ∼30% and ∼75% of the present-day mass of the CC, respectively. Models proposing formation of most (∼90%) of the present-day CC during the initial 1 Gyr or nearly 50–60% during the last 1 Gyr are least favorable. Significant mass exchange between crust and mantle, that is, both the formation and recycling of crust, started in the Hadean with Sm-Nd and Lu-Hf isotope evolution typical for mafic rocks. Depletion of the UM (in incompatible elements) during the early Archean is mitigated by the input of recycled crust, so that the UM maintained a near-primitive Hf-Nd isotope composition. The LM also retained a near-primitive Hf-Nd isotope composition during the Archean, but for different reasons. In contrast to the UM, the crustal return flux into the LM is transiently stored (∼ 1 Gyr) in an isolated reservoir (IR), which limits the mass flux into and out of the LM. The IR in our model is distinct from other mantle reservoirs and possibly related to stable crustal blocks or, alternatively, to recycled crust in the mantle that remains temporarily isolated, perhaps at the core-mantle boundary (LLSVPs).

Growing primordial continental crust self-consistently in global mantle convection models

The majority of continental crust formed during the hotter Archean was composed of Tonalite-Trondhjemite-Granodiorite (TTG) rocks. In contrast to the present-day loci of crust formation around subduction zones and intra-plate tectonic settings, TTGs are formed when hydrated basalt melts at garnet-amphibolite, granulite or eclogite facies conditions. Generating continental crust requires a two step differentiation process. Basaltic magma is extracted from the pyrolytic mantle, is hydrated, and then partially melts to form continental crust. Here, we parameterise the melt production and melt extraction processes and show self-consistent generation of primordial continental crust using evolutionary thermochemical mantle convection models. To study the growth of TTG and the geodynamic regime of early Earth, we systematically vary the ratio of intrusive (plutonic) and eruptive (volcanic) magmatism, initial core temperature, and internal friction coefficient. As the amount of TTG that can be extracted from the basalt (or basalt-to-TTG production efficiency) is not known, we also test two different values in our simulations, thereby limiting TTG mass to 10% or 50% of basalt mass. For simulations with lower basalt-to-TTG production efficiency, the volume of TTG crust produced is in agreement with net crustal growth models but overall crustal (basaltic and TTG) composition stays more mafic than expected from geochemical data. With higher production efficiency, abundant TTG crust is produced, with a production rate far exceeding typical net crustal growth models but the felsic to mafic crustal ratio follows the expected trend. These modelling results indicate that (i) early Earth exhibited a “plutonic squishy lid” or vertical-tectonics geodynamic regime, (ii) present-day slab-driven subduction was not necessary for the production of early continental crust, and (iii) the Archean Earth was dominated by intrusive magmatism as opposed to “heat-pipe” eruptive magmatism.

New insights into Neoarchean–Paleoproterozoic crustal evolution in the North China Craton: Evidence from zircon U–Pb geochronology, Lu–Hf isotopes and geochemistry of TTGs and greenstones from the Luxi Terrane

Tonalite–trondhjemite–granodiorite (TTG) suites constitute dominant components of Precambrian continental crust and provide important insights into the crustal evolution history of the early Earth. Here we investigate a suite of TTGs and associated greenstones from the Luxi Terrane in the south–eastern part of the North China Craton (NCC). We present data from petrology, geochemistry, zircon U–Pb geochronology and Lu–Hf isotopes to characterize these rocks and understand their genesis. Mineral assemblage in the greenstones and TTGs suggests that the rocks have undergone amphibolite–facies metamorphism with regression under greenschist–facies. Zircon U–Pb data from the magmatic grains yield upper intercept ages of 2753 ± 37 and 2582 ± 40 Ma for greenstone, 2552 ± 13 Ma for TTG gneiss and 2476 ± 18 Ma for tonalite, and weighted mean 207Pb/206Pb age of 2507 ± 22 Ma for granodiorite. Lu–Hf isotope systematics show εHf(t) values ranging from −2.13 to 7.29 and crustal model ages (TCDM) in the range of 2603–3125 Ma for the TTGs, whereas the εHf(t) values are up to 8.53 with TCDM ages of 2644–3006 Ma for the greenstone. The results suggest that the protolith for the greenstone was derived from a Meso- to Neoarchean mantle contaminated by the crust source, whereas the TTGs were mostly sourced from juvenile crustal sources with reworking of older basement. The geochemical features of greenstone are similar to island arc basalt, with low Na2O, K2O and TiO2, low total rare earth elements (REEs) values, relative enrichment in light REEs and depletion in heavy REEs absence of obvious Ce anomaly, and negative Nb, Ta and Ti anomalies. The TTGs show similar features, except for the high K2O content and moderate total REEs values. We propose that the greenstone likely represent a part of the early and initial crustal growth at ca. 2750–2600 Ma, and the TTGs were generated subsequently at ca. 2550–2450 Ma. Our results are consistent with the previous models on the crustal growth models in the NCC and provide further insights into the Neoarchean to Paleoproterozoic tectonic history of the craton.

Phanerozoic lower crustal growth from heterogeneous mantle beneath the North China Craton: Insights from the diverse Hannuoba pyroxenite xenoliths

Though many models of continental crustal growth rely on the assumption that juvenile crust is derived from depleted mantle of fairly uniform composition, mafic magmas derived from the mantle – potential candidates for newly formed continental crust – often bear heterogeneous radiogenic isotope compositions. Mafic granulite and some pyroxenite xenoliths from Cenozoic basalts at Hannuoba, in the northern margin of the North China Craton, have younger crystallization ages than the exposed basement rocks, and are generally interpreted to represent juvenile additions to the crust. Here we report mineralogical and geochemical data for the websterite xenoliths from the basalts with the aim of constraining the chemical and isotopic compositions of juvenile continental crust. Two-pyroxene geothermometry suggests that the studied websterite xenoliths originated in the lower crust and last equilibrated at temperatures in the range of 830 to 930 °C. The cumulate texture, together with the geochemical and geochronologic data, indicate that they formed as high-pressure cumulates through crystal fractionation from Phanerozoic mafic magmas. These xenoliths can be subdivided into two groups. Group 1 websterites have low Mg# of 66.1–79.2 and high REE concentrations, and are characterized by extremely evolved radiogenic Sr-Nd-Hf isotopic compositions and high radiogenic Pb ratios. Group 2 websterites have high Mg# of 89.4–89.5, very low REE concentrations, and MORB-like depleted Sr-Nd-Hf-Pb isotopic compositions. We interpret these websterites as additions to the crust formed via magma underplating. The parental magmas for these two websterite groups formed by partial melting of mantle materials with both isotopically enriched (Group 1) and depleted (Group 2) components. The data provide evidence that components from both mantle sources contributed to the underplating magmas and to lower crustal growth beneath the craton.

Episodic crustal production before 2.7 Ga

Before 2.7 Ga, 14 igneous and detrital zircon age peaks and 9 large igneous province (LIP) age peaks are robust and statistically significant. Correlation analysis indicates a synchronous association among these peaks and power spectral analysis shows 91, 114–127 and 182-Myr cycles. These age cycles may be related to mantle plume or mantle overturn events, and to the time it takes to reach threshold temperature gradients for thermo-chemical destabilization in the lowermost mantle. Most zircon age peaks are transferred into younger detrital sediments, which does not favor an origin of the peaks by selective erosion. Correlation of eight pre-2.7-Ga LIP age peaks with zircon age peaks is consistent with a genetic relationship between mantle melting events and felsic crustal production and supports an interpretation of pre-2.7-Ga age peaks as growth rather than preservation peaks produced during craton collisions. Also consistent with the growth peak interpretation is the apparent absence of collisional orogens older than 2.7 Ga. An increasing number of geographic age peak sites from 4 to 2.8 Ga suggests production and survival of only small volumes of continental crust during this time and supports an episodic model for continental crustal growth.

Rates of generation and growth of the continental crust

Models for when and how the continental crust was formed are constrained by estimates in the rates of crustal growth. The record of events preserved in the continental crust is heterogeneous in time with distinctive peaks and troughs of ages for igneous crystallisation, metamorphism, continental margins and mineralisation. For the most part these are global signatures, and the peaks of ages tend to be associated with periods of increased reworking of pre-existing crust, reflected in the Hf isotope ratios of zircons and their elevated oxygen isotope ratios. Increased crustal reworking is attributed to periods of crustal thickening associated with compressional tectonics and the development of supercontinents. Magma types similar to those from recent within-plate and subduction related settings appear to have been generated in different areas at broadly similar times before ∼3.0 Ga. It can be difficult to put the results of such detailed case studies into a more global context, but one approach is to consider when plate tectonics became the dominant mechanism involved in the generation of juvenile continental crust. The development of crustal growth models for the continental crust are discussed, and a number of models based on different data sets indicate that 65%–70% of the present volume of the continental crust was generated by 3 Ga. Such estimates may represent minimum values, but since ∼3 Ga there has been a reduction in the rates of growth of the continental crust. This reduction is linked to an increase in the rates at which continental crust is recycled back into the mantle, and not to a reduction in the rates at which continental crust was generated. Plate tectonics results in both the generation of new crust and its destruction along destructive plate margins. Thus, the reduction in the rate of continental crustal growth at ∼3 Ga is taken to reflect the period in which plate tectonics became the dominant mechanism by which new continental crust was generated.

Geochemistry and chronology of a diorite pluton in the Yinshan Block, implications for crustal growth and evolution of North China Craton

The Yinshan Block constitutes the northwestern segment of the North China Craton. We characterized the geochemistry, zircon U–Pb geochronology, and Hf isotopic composition of a diorite pluton in the Yinshan Block. The zircon U–Pb data yield a protolith formation age of 2535 ± 8 Ma (MSWD = 0.46), followed by recrystallization at 2465 ± 14 Ma or 2431 ± 30 Ma; the granite vein which intruded into the diorite pluton gives a formation age of 2453 ± 5 Ma (MSWD = 1.6), according with the recrystallization age of the diorite. These diorites are characterized by moderate Mg# value and MgO, Cr, and Ni content and have relatively low Sr/Y and (La/Yb)N ratios and slightly Eu anomalies and are enriched in LILEs and depletion in high field strength elements. High‐Sr + Ba granites with low K2O/Na2O characteristics are also shown by these diorites, and the εHf(t) of zircons lie between +1.3 and −7.2, imply these diorites are most likely derived from partial melting of a thickened lower crust with minor involvement of enriched mantle. These geochemical and isotopic characteristics suggest the diorite magmas derived from an accreted arc‐related rocks remelting in a low‐pressure environment along the continent margin. Hereby, we propose a supplement to the model of crustal growth and evolution in the Yinshan Block, involving subduction and collision in the Neoarchean.

Tracing Proterozoic arc mantle Hf isotope depletion of southern Fennoscandia through coupled zircon U–Pb and Lu–Hf isotopes

Constraints on the composition of the depleted mantle Sm–Nd and Lu–Hf crust formation ages have a long history of scientific debate. When calculating mantle extraction ages, and constructing crustal growth models, a linear evolution of incompatible trace elements in a depleted mantle since >4Ga is routinely used. Mantle depletion however varies regionally and over time and subduction of sediments and oceanic crust renders a mantle-wedge variously enriched relative to a modelled depleted mantle.Here we show that primitive mantle-derived subduction related gabbroic intrusions from southern Fennoscandia have Hf isotope compositions that are enriched relative to a MORB-like linear depleted mantle evolution curve. Extrapolation of primitive Paleoproterozoic gabbro suites enables the construction of a regional mantle evolution curve, providing improved constraints on model ages, crustal residence times and the fraction of juvenile versus reworked continental crust. Convergent margins are assumed to be one of the main sites of continental crust growth, and using an overly depleted mantle source yield model ages that are too old, and hence cumulative crustal growth models show too much crust generation early in the Earth's history. The approach of using the Hf isotope composition of zircon from primitive subduction related gabbroic intrusions as a proxy for mantle Hf isotope composition, piloted in this study, can be applied to other convergent margins.

Jurassic Hornblende Gabbros in Dongga, Eastern Gangdese, Tibet: Partial Melting of Mantle Wedge and Implications for Crustal Growth

AbstractThe Gangdese batholith, more than 2500 km in length, is composed mainly of Jurassic‐Miocene igneous rocks. This batholith is one of the most important constituents of the Tibetan orogenesis and provides an ideal place for study of Neo‐Tethyan ocean geodynamic evolution and plateau uplift. Recent studies on the Gangdese Jurassic felsic magmatism highlight its juvenile source. However, important aspects concerning the genesis of the juvenile magmatism and related deep geodynamic evolution are still unclear. Here, we report detailed petrological, geochronological, geochemical, whole‐rock Sr‐Nd isotopic, and in situ Sr‐Hf isotopic data for a recently identified hornblende gabbro in the Dongga area, southern Lhasa sub‐block. This hornblende gabbro is dominated by hornblende and plagioclase, dated at Early Jurassic (ca. 180–190 Ma), and characterized by a narrow compositional range in SiO2 (49.38wt%–52.27wt%), MgO (4.08wt%–7.00wt%), FeO (10.43wt%–11.77wt%), Na2O (2.58wt%–3.51wt%), and K2O (0.48wt%–1.53wt%). It has depleted isotopic signatures, with whole‐rock (87Sr/86Sr)i ratios of 0.7033–0.7043, εNd(t) values of +4.90 to +6.99, in situ plagioclase (87Sr/86Sr)i ratios of 0.7034–0.7042, and zircon εHf(t) of +12.2 to +16.8. Our results integrated with published data suggest a model of Gangdese juvenile crustal growth by a subduction‐related water‐enriched mantle wedge. The hydrous partial melting of the lithosphere mantle was triggered by the dehydration of a Neo‐Tethyan oceanic slab. This mafic magmatism emplaced in the middle‐lower crust of intraoceanic arcs or active continental margins, leading to Jurassic juvenile crustal growth in southern Tibet.