Crustal Differentiation Research Articles

A crustal Pb isotope map of southeastern Australia.

Lead isotopes are a powerful geochemical tracer and a popular tool applied across a broad range of scientific fields, e.g., earth sciences, archaeology, and forensic sciences. Here we present a Pb isotope dataset collected from 232 igneous samples, spanning a ca. 2.3 million km2 area in southeastern Australia, and over 3 billion years of Earth history. This contribution provides a range of isotopic maps showing the spatial variability of Pb isotopes (206Pb/204Pb, 207Pb/204Pb, and 208Pb/204Pb), magma source U/Pb and Th/U, and model ages. The samples selected for this study encompass U- and Th-poor media (i.e., K-feldspar), and U- and Th-bearing sampling media (i.e., whole-rock), providing a temporally and spatially resolved image of U and Th distribution in the crust, and their influence on crustal Pb through radiogenic ingrowth. This dataset has the potential to benefit a wide variety of different disciplines and is an important resource for addressing earth science questions ranging from unravelling crustal differentiation and architecture, through tracing magma source U- and Th-enrichment, to mineral deposit genesis.

Porphyry copper formation driven by water-fluxed crustal melting during flat-slab subduction

AbstractThe prevailing view of the formation of porphyry copper deposits along convergent plate boundaries involves deep crustal differentiation of metal-bearing juvenile magmas derived from the mantle wedge above a subduction zone. However, many major porphyry districts formed during periods of flat-slab subduction when the mantle wedge would have been reduced or absent, leaving the source of the ore-forming magmas unclear. Here we use geochronology and thermobarometry to investigate deep crustal processes during the genesis of the Late Cretaceous–Palaeocene Laramide Porphyry Province in Arizona, which formed during flat-slab subduction of the Farallon Plate beneath North America. We show that the isotopic signatures of Laramide granitic rocks are consistent with a Proterozoic crustal source that was potentially pre-enriched in copper. This source underwent water-fluxed melting between 73 and 60 Ma, coincident with the peak of granitic magmatism (78–50 Ma), porphyry genesis (73–56 Ma) and flat-slab subduction (70–40 Ma). To explain the formation of the Laramide Porphyry Province, we propose that volatiles derived from the leading edge of the Farallon flat slab promoted melting of both mafic and felsic pre-enriched lower crust, without requiring extensive magmatic or metallogenic input from the mantle wedge. Other convergent plate boundaries with flat-slab regimes may undergo a similar mechanism of volatile-mediated lower-crustal melting.

Giant plagioclase in the Steens Basalt, southeast Oregon, USA: Tracking the balance of differentiation processes in continental flood basalt evolution

Giant plagioclase basalts (GPBs), flows that contain abundant and conspicuous plagioclase crystals ≥1 cm in length, exist at many continental flood basalt provinces, and their presence can provide a window into the crustal differentiation processes and storage conditions that produced the magmas prior to eruption. GPB flows occur in the Steens Basalt of southeast Oregon, USA, one of the earliest members of the Columbia River Flood Basalt Group, but their origin in this province has not been studied in detail. We conducted quantitative textural analysis utilizing crystal size distributions (CSDs) at the outcrop and thin section scale, as well as in situ plagioclase trace element and isotopic compositional analysis, in order to elucidate the processes responsible for generating two end-member textural types of GPB found in the Steens Basalt: daisy stone glomerocryst clusters composed of radiating plagioclase laths and single isolated plagioclase. CSD patterns suggest a changing balance of multiple differentiation processes affected the magmas that created each GPB textural type. In the daisy stone plagioclase, crystals rarely exhibit resorption textures, trace element abundances are comparatively low with minimal variability, and Sr isotopes are commonly in equilibrium with groundmass. In the single isolated plagioclase, resorption textures are common, trace element abundances are generally higher and display within crystal variation, and Sr isotopes are typically not in equilibrium with groundmass. These distinctive characteristics require separate mechanisms to generate each end-member GPB texture. Daisy stone GPBs form by growth of large crystals following mafic recharge events that completely resorb entrained cumulate plagioclase. A large amount of plagioclase may then grow with minimal cooling following recharge due to saturation in plagioclase components. In contrast, single isolated GPBs are the result of partial resorption of both entrained and resident plagioclase following less voluminous and less primitive recharge events; regrowth then occurs on the partially resorbed cores when cooling resumes. Both GPB types are formed by thermal cycling, but each texture reflects the differentiation processes that were prevalent during different eruptive stages of the Steens Basalt, as the balance between recharge, crystal fractionation, and assimilation varied in time.

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.

Xenolith Petrochronology (San Luis Potosi, Mexico) Constrains Heat Sources for Cenozoic Ultrahigh‐Temperature Metamorphism in the Lower Crust

AbstractUltrahigh‐temperature (UHT; >900°C) metamorphism drives crustal differentiation and is widely recognized in the rock record, but its geodynamic causes are debated. Previous work on granulite‐facies metapelite xenoliths from San Luis Potosí, Mexico suggests the lower crust experienced a protracted UHT metamorphic event that coincided with the onset of regional extension. To determine the duration, conditions, and heat sources of UHT metamorphism recorded by these xenoliths, this study characterizes the major‐element, trace‐element, and U‐Pb isotopic systematics of quartz, rutile, feldspar, garnet, and zircon by in situ electron microprobe (EPMA) and laser‐ablation inductively coupled‐plasma mass spectrometry (LA‐ICP‐MS), and augments these data with detailed petrography, thermobarometry, phase equilibria modeling, and diffusion modeling. Thermobarometry and phase equilibria modeling suggest peak metamorphic conditions exceeded 0.7 GPa and 900°C. Zircon petrochronology confirms >15 Myr of UHT conditions since its onset at ∼30 Ma. A small population of zircon record elevated temperatures following transition from backarc compression to extension during the waning stages of orogenesis (60–37 Ma). Garnet preserves trace‐element zoning and mineral inclusions consistent with suprasolidus garnet growth and subsequent compositional modification by intracrystalline rare‐earth element diffusion during protracted heating, with diffusion chronometry timescales in agreement with zircon data, followed by fluid‐driven remobilization of trace elements along now‐healed fractures within ∼1 Myr of eruption. In sum, these data are most compatible with lithospheric mantle attenuation or removal as the dominant heat transport mechanism driving synextensional UHT metamorphism and crustal melting, which has bearing on models for crustal differentiation and formation of modern and ancient granulite terranes globally.

Experimental Constraints on the Storage Conditions and Differentiation of High-Ti Basalts from the Panzhihua and Hongge Layered Intrusions, SW China

Abstract High-Ti basalts are commonly believed to represent parental magmas leading to the formation of mafic-ultramafic layered intrusions, such as Late Permian Panzhihua and Hongge in the Emeishan Large Igneous Province (SW China). Consequently, elucidation of the crystallisation and crustal differentiation of high-Ti basalts is critical for our understanding of the petrogenesis of these layered intrusions and the associated oxide ore mineralisation. Here, we present the results of crystallisation experiments carried out in internally heated pressure vessels using a primitive high-Ti basaltic composition. The experiments were conducted at 100 and 300 MPa, in the temperature interval of 950–1200 °C and with water activities (aH2O) from 0 to 1. The oxygen fugacity (fO2) was controlled and varied from FMQ −1 to FMQ +3.3 log units (FMQ corresponds to fayalite-quartz-magnetite buffer). The main mineral phases are olivine, clinopyroxene and plagioclase, accompanied by Cr-Fe-Ti-oxides, orthopyroxene, apatite and amphibole, depending on the conditions. Redox conditions primarily influence the stability fields of Cr-Fe-Ti oxides. Clinopyroxene, orthopyroxene and amphibole are pressure-dependent and have larger stability fields under high pressure conditions. The olivine→orthopyroxene and olivine→amphibole peritectic reactions are observed. Comparisons of phase equilibria between this study and experiments conducted with parental magma of Skaergaard layered intrusion demonstrate the effect of bulk system composition. For instance, ilmenite crystallisation is determined not only by intrinsic parameters such as fO2 but also by additional compositional parameters (e.g. melt Ti, Fe, Al and Mg content). Although COMAGMAT and MELTS modelling results generally reproduce the crystallisation sequence, only the stability field of clinopyroxene and its composition are perfectly modelled. The comparison of experimental results with the rocks from the lowest units of Panzhihua and Hongge layered intrusions are applied to constrain storage conditions in the magma reservoirs. Compared to Hongge, we conclude that the Panzhihua magma chamber was probably located at a shallower depth (~3–6 km), that magma crystallisation started at lower temperatures (~1125–1100 °C), higher fO2 (~FMQ + 1 to FMQ + 2) and that its initial melt H2O content was lower (~0.5–1 wt. %).

Testing the TTG–Metabasite Connection in the Southern Superior Province: an Integrated Geochemical, Isotopic, and Petrogenetic Modelling Approach

Abstract Archean cratons are dominated by tonalite–trondhjemite–granodiorite (TTG) suites, the products of crustal differentiation that formed early continental crust. These rocks may have been primarily generated by partial melting of hydrated basaltic crust in a variety of settings including subduction zones or the deep lithosphere. Sources are mainly inferred from examination of TTGs alone, as potential source rocks are rarely exposed. In the southern Superior Province, Canada, the Kapuskasing Uplift exposes an important crustal cross-section with upper- to middle-crustal TTGs and lower-crustal metabasites, which show evidence of having produced trondhjemitic anatectic melt. Here, we test the hypothesis that these metabasites were the source of the Mesoarchean to Neoarchean TTGs in the Kapuskasing Uplift by combining phase equilibrium and melt trace element modelling with whole-rock and zircon Lu–Hf isotope analysis and geochronology (garnet Lu–Hf and zircon U–Pb) of metabasic samples. By comparison of the results with existing data from TTGs in the Kapuskasing Uplift, we determined that the metabasites are plausible source rocks for the TTGs. The Lu–Hf systematics of the metabasites and TTGs are the most robust evidence of a genetic connection. Modelling results support an increase in TTG source depth over time. New geochronological data constrain partial melting of metabasite and crystallization of anatectic melt between ca. 2685 Ma and ca. 2600 Ma, coeval with crystallization of only the youngest TTGs. Overall, these results indicate a complex history of intracrustal differentiation in the Kapuskasing Uplift, with partial melting of two isotopically distinct lower-crustal metabasic sources at different times and depths.

Mineral Textures and Chemistry Trace the Origin and Transcrustal Evolution of the Sanyuangou Syenite in Southern North China Craton

Abstract Alkaline rocks are widely distributed in various geological environments and are important carriers for exploring the formation and compositional diversity of continental crust. Extensive studies have investigated the processes of mantle melting and crustal differentiation that produced such rocks. However, the potential interaction processes between mantle-derived magma and crust during their formation are poorly constrained. In this study, we focus on a Paleoproterozoic garnet-bearing syenite in the Xiong'er large igneous province (LIP) on the southern margin of the North China Craton through detailed whole-rock and mineral analyses. The high K2O (7.4–8.8 wt %) syenite emplaced at ~1772 Ma with ancient inherited zircon (1800–2800 Ma). The dominant mineral assemblage is composed of clinopyroxene, garnet, and K-feldspar with a small amount of titanite. Complex compositional oscillatory zoning of clinopyroxene and garnet indicates that magma mixing played a significant role in the formation of the syenite. We estimated the major elements composition for melts in equilibrium with clinopyroxene and calculated the clinopyroxene/melt partition coefficients by crystal lattice strain model, thus calculating the trace elements of the equilibrated melts. The equilibrium melts of high Mg# (>85) clinopyroxene have high CaO/Al2O3 (>0.6) and low Hf/Sm (<0.4), which suggests a role for carbonatitic metasomatism of the mantle. The variables La/Yb (24.4–56.1), Dy/Yb (0.8–5.9), and Yb (0.6–10.3 μg/g) indicate that the initial melts were formed by 1% to 2% partial melting of spinel-garnet phlogopite lherzolite. An REE-based model for melts in equilibrium with low Mg# (<85) clinopyroxene indicates that 10% to 30% felsic magma from ancient crust participated in hybridization. In addition to well-documented magma mixing, the oscillatory zoning of garnet is also related to competition with titanite. The initial alkaline magmas have a high water content (~4 wt %), which delays the crystallization of K-feldspar, leading to the enrichment of K2O, until the K-feldspar accumulates in the shallow crust to form the syenite. Magma mixing under an open system further leads to alkalinity enrichment. Magma source and crustal evolution jointly determine the potassium-rich characteristics of syenite. Multiple episodes of magma mixing and fractional crystallization occur in the crust (700–300 MPa), suggesting complex and vertically extensive magma chambers. This study represents the first identification of carbonatitic metasomatism as a mantle source in the Xiong'er LIP. Furthermore, it offers a new perspective on magma mixing between the mantle and crust in transcrustal magmatic systems, contributing to the formation of alkaline rocks.

Zinc isotopes of high silica granites in the Qilian Orogenic belt and implications for mineralization potential

Many high silica granites (>70 wt% SiO2) are closely linked to ore deposits. There has been an ongoing debate about whether they are formed through extensive fractional crystallization or partial melting, as traditional major-trace element and radiogenic isotopes make it challenging to distinguish among these mechanisms. Recent studies suggest that zinc isotopes have the potential to trace crustal differentiation and monitor fluid evolution during the late stage of magma evolution. In this study, we investigated the Lenglongling leucogranites in the North Qilian Orogenic belt, characterized by high SiO2 contents ranging from 72 % to 78 wt%, extremely low Zr/Hf ratios (15.8–21), low Nb/Ta ratios (3.8–9.3), and strongly negative Eu/Eu* anomalies (0.04–0.32). The invariant and basalt-like zinc isotopic compositions of the Lenglongling high silica leucogranites (0.27 ‰ ±0.03 ‰ to 0.38 ‰ ± 0.03 ‰) reveal that they did not undergo extensive fractional crystallization or fluid exsolution. Instead, their homogeneous zinc isotope compositions and distinct geochemical characteristics were inherited from their protoliths during partial melting.

Zircon Zr[sbnd]Hf isotope disequilibrium during crustal anatexis: A record from the Mesozoic migmatite in the eastern Gangdese arc, South Tibet

Stable zirconium (Zr) isotopes have become the focus of attention as potential tracers of magma crystallization and differentiation processes. However, the effects of partial melting, especially disequilibrium melting, on Zr isotope fractionation remain unclear. Here we report in-situ zircon UPb ages, major and trace elements, and ZrHf isotope data for the migmatites in the eastern Gangdese arc, South Tibet. The leucosomes from the metatexite and diatexite migmatites are generally concordant with the foliation of the paleosome and display diffuse and gradational contacts with the melanosomes, indicating that they are in-source leucosomes. Zircon in the paleosome, the protolith of the migmatites, have UPb age of 361 ± 2 Ma and low 176Hf/177Hf ratios of 0.282399–0.282482. They exhibit intragrain variation in δ94/90Zr values, increasing from the core (−0.10 ‰ to +0.36 ‰) to the rim (+0.14 ‰ to +0.55 ‰) domains. The negative correlation between δ94/90Zr values and Zr/Hf ratios suggests preferential incorporation of light Zr isotopes during the crystallization of the zircon in the protolith. Most zircon in the leucosomes from the metatexite and diatexite migmatites exhibit distinct overgrowth rims around the inherited cores. The inherited cores have consistent UPb ages (ca. 360 Ma), trace elements, and Hf isotopic compositions (0.282362–0.282449) with zircon in the paleosome, indicating in situ partial melting. The inherited cores display δ94/90Zr values ranging from −0.12 ‰ to +0.31 ‰, consistent with the zircon core domains in the paleosome. The overgrowth rims have UPb ages of ca. 102 Ma and display higher 176Hf/177Hf ratios (0.282592–0.282900) and δ94/90Zr values (+0.01 ‰ to +0.63 ‰) than their inherited cores. The elevated 176Hf/177Hf for the overgrowth rims resulted from the retention of unradiogenic 177Hf in the residue due to limited zircon dissolution and more release of radiogenic 176Hf to the melts because of the garnet breakdown, indicating disequilibrium melting. The elevation of δ94/90Zr for the overgrowth rims relative to the inherited cores was not caused by the disequilibrium melting of other Zr-bearing minerals but was mainly controlled by the Zr isotopic composition and dissolution degree of the zircon in the protolith. The preferential melting of isotopically heavy zircon rims would elevate the δ94/90Zr values in the melts relative to their zircon-bearing protoliths. Crustal anatexis may be another important mechanism contributing to the Zr isotope fractionation during crustal differentiation.

Continental Arc Accumulation Mafic Rocks in the Mid-Upper Crust: Constraints From the Early Eocene Hornblende Gabbro–Diorite in the Tengchong Block, Southeastern Extension of Tibet

Abstract Magmatic activity in the syn-collision stage is key for net crustal growth. To understand the mechanism of accretion–differentiation and compositional change of the continental crust, it is important to focus on the magmatic activity during the syn-collision stage. Early Eocene mafic–ultramafic rock assemblages found in the western part of the Tengchong Block resulted from a continuous series of arc magmatic evolution, thoroughly recording the continental arc magmatic system during the subduction of the Neo-Tethys Ocean and syn-collision of the Indian-Asian continents. Early Eocene hornblende gabbro–diorite in the Tengchong Block formed at 53 Ma, and the primitive magma was derived from an enriched mantle source due to the enriched Nd–Hf isotopes. The amphibole and biotite thermobarometer measurements indicate that the mafic magma reservoirs in the Tengchong Block occurred at a mid-upper crust. Petrography, amphibole Fe/Mg exchange coefficient (KD), Rayleigh fractionation, and equilibrium melt calculation indicate that the Early Eocene hornblende gabbro–diorite in the Tengchong Block was created due to plagioclase-dominated accumulation at the mid-upper crust level. Based on the calculation, the corresponding amphibole equilibrium melt is more silicic (dacitic–rhyolitic in composition) than the bulk rocks, indicating a more evolved composition in the mid-upper crust. Three types of plagioclases reveal the multi-recharging and dissolution–reprecipitation promoting the further evolution of these mafic rocks. Therefore, this study concludes that magma recharge and plagioclase-dominated accumulation processes may be important mechanisms for the formation and evolution of mafic magma and the further crustal differentiation at the mid-upper crust level in a continental margin arc.

Bimodal volcanism in the Hegau region (SW Germany): Differentiation of primitive melilititic to nephelinitic rocks produces evolved nosean phonolites

Peculiar bimodal volcanism in the Hegau region (SW Germany) comprises two contrasting SiO2-undersaturated rock suites.(I)Primitive olivine melilitites and melilite-bearing olivine nephelinites (12–9 Ma) are characterized by high MgO, CaO, Fe2O3, TiO2, Ni, V, F, moderate alkalis, Al2O3, P2O5, Ba, Nb, Zr, and low SiO2, Rb, Pb, and U concentrations. The rocks are composed of forsteritic olivine, diopsidic clinopyroxene, melilite, perovskite, Cr-bearing oxyspinel, F- and Ba-rich mica, and fluorapatite. In rare cases, they contain coeval coarse-grained ijolite patches generated by rapid in-situ fractionation in small melt pockets.(II)Evolved nosean phonolites (14–11 Ma) comprise high alkalis, Al2O3, SiO2, Rb, Nb, Zr, U, Pb, S, and low MgO, CaO, Fe2O3, TiO2, P2O5, Ba, Ni, and V concentrations, and contain abundant Ba-bearing alkali feldspar and nosean-haüyne-sodalitess macrocrysts, aegirine-augitic clinopyroxene, as well as accessory apatite, titanite, zircon, and pyrochlore.The melilititic–nephelinitic rocks formed by low degrees of partial melting of a carbonated amphibole ± phlogopite-bearing garnet wehrlite in the uppermost asthenosphere or in the thermal boundary layer and occur also in neighbouring regions. However, their coexistence with evolved nosean phonolites is a unique and so far, unexplained feature in the southern Central European Volcanic Province.Thermodynamic modelling implies that removal of 11–19% oxyspinel, 4–10% olivine, 42–57% clinopyroxene, <3% mica, <9% feldspathoids, <8% feldspar from melilititic–nephelinitic parental melts at upper crustal conditions (∼200 MPa) results in significant amounts of phonolitic residues that are compositionally similar to the exposed nosean phonolites. Strongly negative P and Ti anomalies and a trough of MREE in primitive mantle-normalized trace element patterns of phonolites indicate additional fractionation of titanite and apatite, consistent with the mineralogy of coarse-grained (nepheline) syenitic cumulates, which are present as enclaves in both rock suites. The modelling results suggest crystallization of the (nepheline) syenite cumulates between 1050 and 800 °C and ascent and eruption of the phonolite residues at no less than ∼900 °C with complete solidification of the observed assemblage at >750 °C. Significant crustal assimilation during fractionation appears unnecessary to explain the mineralogical, mineral chemical, and geochemical characteristics of the phonolites. Prolonged upper crustal differentiation of the magmas in one case (nosean phonolites) and fast ascent of primitive melts in the other (olivine melilitites and melilite-bearing olivine nephelinites) can be explained by stress field changes from an extensional to a more compressive regime, the magma ascent becoming thereby increasingly structurally controlled and supported by brittle deformation.

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

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

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.

Petrology and Geochemistry of Adak Island Plutonic Xenoliths: Implications for Primitive Magma Generation and Crustal Differentiation in the Aleutian Island Arc

Abstract Deep crustal cumulates in arcs offer a window into the chemistry and crystallization conditions (P–T–H2O–fO2) of primitive basalts in the upper mantle and lower crust and can be studied in ancient exhumed terranes or in xenoliths erupted in young arc lavas. Here, we expand on previous studies and thoroughly characterize the extensive xenolith suites erupted from the Mt. Moffett and Mt. Adagdak volcanic centers (Adak Island, Central Aleutians), which range from primitive ultramafic cumulates to more evolved amphibole gabbros and hornblendites. We present detailed petrography as well as in situ trace and major element mineral chemistry. We use these data to calculate pressure, temperature, and fO2 estimates for the xenoliths, and compare these findings to experimental results to understand the crystallization sequence and P–T–H2O–fO2 under which the cumulates formed. The Moffett crystallization sequence is defined by early amphibole fractionation and an abrupt shift in oxide compositions from chromite to magnetite, while the Adagdak suite is characterized by simultaneous saturation of amphibole+plagioclase and oxide compositions that become increasingly aluminous before magnetite saturation. Olivine–spinel oxybarometry of the Adagdak xenoliths indicates that they are oxidized relative to mid-ocean ridge basalt (MORB:FMQ +0.1 to +2.1). Highly fractionated REE and elevated Sr/Y ratios are observed in clinopyroxene from the most primitive cumulates, consistent with a contribution from a basaltic eclogite melt. This basaltic eclogite melt is hypothesized to come from partial melting of the slab or through melting of basalt introduced into the subarc mantle through forearc subduction erosion. These signatures are greatly diminished in the more evolved lithologies, which can be explained through fractionation of plagioclase and amphibole. Our findings support the presence of a complex magmatic plumbing system beneath Adak, with Mt. Moffett and Mt. Adagdak volcanic centers tapping compositionally distinct sources. More broadly, our results are consistent with studies suggesting that low-degree basaltic eclogite melts through slab melting or forearc subduction erosion contribute to arc magmas in the Aleutians, although the associated geochemical signatures are easily obscured by differentiation in the crust.

Revealing the link between A-type granites and hottest melts from residual metasedimentary crust

Abstract Among S-, I-, and A-type granites, the latter are the most debated in terms of origin, and metasedimentary crust is usually discarded as a potential source. Here we tackle this issue by adopting an in-source perspective, rather than focusing on the final product (granite), documenting the occurrence of pristine melt inclusions (MIs) in garnet from residual metapelitic ultrahigh-temperature (UHT) granulite from East Antarctica. Coexistence of sapphirine + quartz, phase equilibria calculations, and Zr-in-rutile thermometry indicate that MIs trapped UHT melts formed at peak conditions (930–1000 °C) from a residual metapelitic source. MIs are granitic with weakly peraluminous to weakly peralkaline affinity, ferroan character, high alkali contents, high K/Na and Ga/Al, and low Ca, Ba, Sr, and H2O concentrations. These features and geochemical modeling indicate that MIs represent primary melts for high-SiO2 A-type granites. Therefore, MIs reveal the missing link between A-type granites and the hottest metasedimentary crust. Voluminous amounts of slightly peraluminous, high-SiO2 A-type granites can be produced in large, residual UHT terranes such as those of eastern Gondwana. Our results provide a wider view of processes responsible for granite formation and show that a larger variety of granites must be considered in models of the effects of UHT anatexis on crustal differentiation.

Evolutionary stasis during the Mesoproterozoic Columbia-Rodinia supercontinent transition

The transition between the supercontinents Columbia and Rodina coincided with a delay in eukaryote expansion during the Mesoproterozoic, however, the cause for this association is uncertain. Here, we use statistical geochemical analyses of igneous and fine-grained siliciclastic rocks to demonstrate that extensive crustal differentiation occurred during this transition interval. The results show a relative increase in abundance of phosphorus-poor felsic volcanic and plutonic rocks and a prevailing low weathering intensity between ∼ 1.8 and 1.2 Ga (billion years ago). The decelerated weathering of phosphorus-poor felsic volcanic and plutonic rocks could have maintained the low flux of bio-essential nutrients to the oceans that sustained low primary productivity and atmospheric O2 levels, which inhibited biologic radiation during the Columbia-Rodinia supercontinent transition period.

Timescales of ultra-high temperature metamorphism and crustal differentiation: Zircon petrochronology from granulite xenoliths of the Variscan French Massif Central

Lower crustal ultra-high temperature (UHT) metamorphism and the resulting production and transfer of granitic magmas represent key processes of intracrustal differentiation. The timescales of these phenomena are debated, due to the complexity of the granulite zircon U-Pb record and because direct links between granulites and granites are difficult to establish. To address these issues, we present the results of zircon petrochronology (coupled U-Pb dating, trace element and Lu-Hf isotopic analyses) for lower crustal felsic and mafic granulite xenoliths from the Variscan eastern French Massif Central and compare them with data from well-characterized mid-/upper crustal migmatites and granites. The felsic and mafic granulites represent pre-Variscan meta-sedimentary and meta-igneous mafic rocks, respectively; which experienced Variscan UHT peak metamorphism at 940–970°C and 8 ± 2 kbar. Zircons from the felsic granulites show U-Pb dates spreading over ∼50 Myr between ∼315 and ∼265 Ma, correlated to Ti-in-zircon temperatures decreasing from 940–970°C to 800°C and REE contents consistent with growth in equilibrium with garnet at (U)HT conditions. This is best explained by continuous crystallization of zircon upon cooling from the thermal peak owing to decreasing Zr solubility of residual melt, despite significant prograde melt loss. Zircons from the mafic granulite only record the last stage of this time-temperature evolution (295–265 Ma; 800–900°C) due to later zircon saturation. Upper crustal granite emplacement started at 340 Ma and culminated at the age of the lower crustal thermal peak of 313 ± 3 Ma (defined by the highest-temperature zircons from felsic granulites), reflecting melt extraction along the ∼27 Myr prograde path of the lower crust. In turn, the crystallization ages of mid-crustal migmatites (315–300 Ma) and the lower crustal granulites (315–265 Ma) are consistent with slow cooling in the presence of melt. These results provide a direct assessment of the timescales of melt production and residence at the crustal scale; and validate the granulite-granite connection in the framework of the melt loss theory in migmatitic systems.

The different responses of trace elements to equilibrium and disequilibrium melting: Implications for crustal differentiation and granite compositions

Equilibrium and disequilibrium meltings are two common mechanisms in the geochemical differentiation of the continental crust. However, responses of trace elements to the two melting processes are still not fully understood in granite petrogenesis. To address this, we conducted trace element modelling of open system melting using phase equilibrium calculations combined with solubility equations for accessory minerals and mineral–melt partition coefficients based on three melt-reintegrated felsic granulite compositions. The results indicate systematic differences between the compositions of disequilibrium (Type I) and equilibrium (Type II) melts involving water-present melting reaction, fluid-absent melting reaction involving biotite, and fluid-absent melting reaction in the absence of biotite. Type I melts generally have higher Rb, Ba, Co, Ni, V, and Cr contents, and Rb/Sr and Nb/Ta ratios, and lower Sr, P, Th, U, and light rare earth element contents than Type II melts. Significantly, disequilibrium melting involving biotite can produce melts with high Nb contents and super-chondritic Nb/Ta ratios, which may be a reservoir for the missing Nb on Earth. Compared with Type II melts, the trace element contents and ratios of Type I melts exhibit a larger range, more complex trends, and are more sensitive to changes in melting reaction types and melt loss. This suggests that, in addition to petrogenetic processes such as crystal entrainment, fractional crystallization and magma mixing, the disequilibrium melting would also contribute to the compositional diversity of granites. The modelling results also show that the trace element signatures of granites that are commonly used to determine their petrogenetic processes and tectonic setting may lead to erroneous conclusions if the granitic melts were produced by disequilibrium melting.

Andesites and evolution of the continental crust: Perspectives from the Central Volcanic Zone of the Andes

Named for the Andes, andesites (53%–63% SiO2) are the archetypal magma erupted at magmatic arcs. They have been established as the average composition of continental crust and as such are integral to the growth and evolution of the continental crust. However, andesites are quite variable in trace element and isotopic composition reflecting disparate paths of origin. Herein we return to the original site of their identification, the Central Andes, and use a comprehensive dataset of published and unpublished trace elements and isotopes to show that during the past 6 Myr two distinct types of andesite have erupted in the Central Volcanic Zone (CVZ), which correspond with different geodynamic conditions. Consistent with previous work, we confirm that major composite cones and minor centers of the steady state (low magmatic flux) Quaternary CVZ arc have trace element and isotopic characteristics consistent with magma generation/fractionation in the lower crust. Within the Quaternary arc centers, there are also significant latitudinal variations that correspond with the age, composition, and P-T conditions of the lower crust. However, in contrast to this prevailing model, in the 21–24°S segment 6–1 Ma andesites from ignimbrites and lava domes associated with the peak of the regional Neogene ignimbrite flare-up have compositions that indicate these andesites are hybrids between mantle-derived basalts and upper crustal lithologies. Since ∼1 Ma, andesites in young silicic lava domes associated with the regional flare-up are compositionally indistinguishable from proximal Quaternary arc centers, indicating a return to steady-state magmatism and lower crustal production of andesites. We propose that the transition from upper crustal to lower crustal andesite production results from a decrease in mantle heat input and subsequent relaxation of the regional geotherm during the waning of the flare-up event. The two modes of andesite production have significant implications for the production and evolution of the CVZ arc crust. During the flare-up, prodigious amounts of basalt were emplaced into the mid-crust, resulting in the production of large volumes of hybrid intermediate magmas in the mid and upper crust. In contrast, the lower crustal differentiation recorded in the Quaternary steady state arc andesites would result in the formation of a dense crystalline residue in the lower crust and an overall densification of the lower crust. Over time, gravity instabilities associated with this densification may ultimately aid in the delamination of the dense lower crustal root, triggering flare-ups. These differences in andesite production may help explain the cyclicity (flare-up cycles) observed in mature continental arcs and emphasizes that andesite is not a monotonous composition and can vary with depth-dependent intra-crustal differentiation related to magmatic flux.