Trace Element Modeling Research Articles

LOMU Type Alkali Basalts in East Asia Sourced from Subduction Recycling of Multiple Ancient Crustal Components

Abstract The origin of alkali basalts with distinctive time-integrated low U/Pb (low μ, LOMU) from East Asia is controversial due to the complex geochemical and tectonic signatures reported from this region. We report new data on the petrology and geochemistry of the Higashi-Matsuura and Kita-Matsuura alkali basalts from Southwest Japan, which confirm the presence of a LOMU-type mantle component below the Japanese Islands, similar to East Asian intraplate volcanoes. We use whole rock geochemistry and mineral chemistry to determine mantle melting, magma evolution conditions and isotopic evolution of the LOMU mantle source. Petrological studies show that the Higashi-Matsuura alkali basalts (~3 Ma) were derived from a hydrous mantle source with ~880 μg/g H2O, at a pressure of 1.8 to 2.3 GPa at a mantle potential temperature of ~1300°C. Trace element modelling concludes that these basalts may have been generated through the hydrous upwelling and partial melting of the asthenospheric mantle containing ~15% dehydrated eclogitic component. Higashi-Matsuura alkali basalts show 206Pb/204Pb values of 17.72 to 18.04 which are among the lowest values from Southwest Japan. Relatively older (6–8 Ma) alkali basalts from the Kita-Matsuura area showing similar physicochemical characteristics, do not show LOMU-type isotopic trends. Trace element and Pb–Sr–Nd isotopic data indicate that the Higashi-Matsuura mantle component is similar to the extreme LOMU components reported from the Northeast China alkali basalts, as well as the recently discovered Petit Spot volcanoes on the Pacific Plate. We model the origin of the LOMU signature from the lowest reported Pb isotope ratios in East Asia, from intraplate alkali basalts in Northeast China. Our model suggests that at least two separate subduction events of marine sediments, at 1.8 and 2.2 Ga, are required to explain the observed Pb isotopic variation in the East Asian region. Other LOMU type basalts from East Asia, including Southwest Japan and Petit Spot, define a linear trend between the Xiaogulihe basalts and lithospheric mantle xenoliths, with 40% to 60% Pb in Higashi-Matsuura basalts being derived from the lithospheric mantle during upwelling. This suggests that the LOMU array in East Asia may have been formed by mixing between multiple ancient, subducted sediment components accumulated at the mantle transition zone for about 2 billion years, and its recent upwelling due to dehydration of the stagnant Pacific slab and related melting of the metasomatised asthenospheric mantle.

Geodynamic Evolution of the Proto-Tethys Ocean in the West Kunlun Orogenic Belt, northwest Tibetan Plateau: Implications from the Subarc Crust and Lithospheric Mantle Modification

Abstract Radiogenic isotopes serve as a crucial tool for investigating crustal evolution, playing a pivotal role in revealing magma sources and petrogenesis. However, they can be ineffective in distinguishing between distinct magmatic sources with similar radiogenic isotopic compositions, a common occurrence in nature. Here we addresse this challenge in the Ordovician igneous rocks from the West Kunlun orogenic belt (WKOB), aiming to distinguish between two potential magmatic sources (i.e. the Tarim Craton and the Tianshuihai terrane) with similar isotopic compositions using appropriate thermodynamic and geochemical modeling based on mineral and whole-rock geochemistry. Zircon U–Pb dating yields ages of 483 ± 3 Ma for the Pushou gabbros and 469 ± 2 Ma and 461 ± 2 Ma for the Datong monzogranites and syenites, respectively. The Pushou gabbros exhibit low SiO2 (47.4–49.1 wt %), high MgO (5.5–6.9 wt %), high large-ion lithophile elements (LILEs, e.g. Rb, Ba, Th, and K), and low high field-strength elements (HFSEs, e.g. Nb, Ta, Zr, Hf, P, and Ti), suggesting an origin in subduction-modified mantle. They display high whole-rock (87Sr/86Sr)i ratios (0.7156 to 0.7192), negative whole-rock εNd(t) values (−7.1 to −7.8), as well as high zircon δ18O values (7.6–7.9‰) and enriched zircon Hf isotopic compositions (εHf(t) = −5.3 to −7.7), which are consistent with 1–5% subducted sediments in an enriched mantle source. Trace element models further confirm that the gabbros are most likely derived from low-degree (~15%) partial melting of subduction-modified Tarim mantle in the spinel–garnet facies rather than from the Tianshuihai mantle. The Datong syenites belong to the shoshonitic series and are characterized by medium SiO2 (59.5–61.4 wt %), relatively low MgO (0.9–1.2 wt %) and Mg# (37–42), enrichment in LILEs and depletion in HFSEs. They have high whole-rock (87Sr/86Sr)i ratios (0.7103 to 0.7105) and negative whole-rock εNd(t) values (−3.8 to −4.3), along with negative to slightly positive zircon εHf(t) values (−3.8 to +2.6), similar to coeval mafic rocks. Thermodynamic and geochemical modeling suggest that the Datong shoshonitic rocks likely originated via crystal fractionation of shoshonitic basaltic magmas in the SW Tarim Craton. The Datong monzogranites have high SiO2 (69.7–72.6 wt %), low MgO (0.6–0.7 wt %), and a typical enrichment in alkalis, Zr, and Nb, with depletion in Sr, P, and Ti, consistent with A-type granites. They are characterized by high whole-rock (87Sr/86Sr)i ratios (0.7321 to 0.7323), negative whole-rock εNd(t) values (−11.3 to −11.8), negative zircon εHf(t) values (−11.0 to −16.5), and high zircon δ18O values (7.2–8.0‰), indicating derivation from the remelting of an ancient crustal source. Thermodynamic, major, and trace element modeling indicate that their parent magma may have been generated by water-deficient (~2 wt %) partial melting of ancient crustal material beneath the SW Tarim Craton rather than that of the Tianshuihai terrane, under high-temperature (T > ~950°C) and low-pressure (P = 5–8 kbar) conditions. Based on the tectonic framework of the WKOB, we propose that the original mantle and crust beneath the southern Kunlun terrane may have been modified or partially replaced by that beneath the SW Tarim Craton during the Ordovician. Therefore, this evidence for Tarim-derived magmatism, when combined with regional sedimentary and structural records, indicates that Ordovician magmatism in the southern Kunlun terrane is most consistent with episodic northward subduction of the Proto-Tethys Ocean, commencing at ~485 Ma. Middle Ordovician slab break-off can explain the formation of the A-type granites, but reinstated northward subduction is required for the generation of late Ordovician Datong syenites.

Partial melting of arclogite and petrogenesis of alkaline-silicate complexes

Magmatic processing in the lower crust of thick-crusted (>40 km) arcs generates garnet-rich clinopyroxenite residues (arclogites) that are denser than the underlying mantle and therefore susceptible to foundering. Arclogites, which can contain substantial amphibole, may undergo partial melting as they sink into the subarc mantle. To constrain the geochemical composition of arclogite-derived melts and the role of arclogite in magmagenesis, we conducted partial melting experiments on average amphibole-bearing arclogite ARC15 at 2 GPa from 1040 to 1300 °C. At subsolidus conditions, garnet and clinopyroxene coexist with lesser amounts of amphibole, ilmenite, and biotite. Hydrous ARC15 begins melting between 1040 and 1080 °C. Once amphibole and biotite are exhausted by 1140 and 1190 °C, respectively, garnet becomes the primary melting phase according to the reaction 0.8 Garnet + 0.2 Clinopyroxene = 1 Melt. At 1300 °C, ARC15 yields 18 wt.% melt, indicating an anomalously low melt productivity of 0.12%/°C and rendering it the least melt-productive of all experimental pyroxenites investigated at 2 GPa. Geochemically, melts sourced from ARC15 are predominantly basanitic with low SiO2 (43–52 wt.%) and high FeO (10–17 wt.%), TiO2 (2.5–3.7 wt.%), and alkali (Na2O + K2O = 3–6.5 wt.%) contents. Such compositions do not resemble magmas erupted at thick-crusted arcs, indicating that arclogites are not a dominant source of magmatism in these settings. Rather, the compositions of ARC15 melts more closely match those of magmas emplaced at post-collisional alkaline-silicate complexes, which form as a result of orogenic collapse and subsequent continental rifting and are common hosts of economic rare-earth element (REE) deposits. Trace element modeling further indicates derivation of alkaline-silicate magmas from variably metasomatized, light REE-enriched arclogite. Our results suggest a geodynamic model in which subduction beneath thickened crust leads to arclogite formation, foundering, and metasomatism as the material descends through the subduction-influenced mantle. During this stage, partial melts sourced from arclogite are likely modified by subduction-related secondary processes which dilute their diagnostic geochemical signatures. Subsequent collision and orogenic collapse induce extensional stresses that facilitate asthenospheric upwelling, leading to partial melting of the arclogite-bearing mantle and efficient extraction of primitive arclogite melts to the surface to form alkaline-silicate complexes and associated REE deposits.

Suturing of the Archean Bastar craton with the Eastern Ghats Province to form the Greater Indian Landmass: Insights from geochemistry, U-Pb geochronology and phase equilibria modelling

The contact between the Archean Bastar craton (BC) and Proterozoic Eastern Ghats Province (EGP), central India, is marked by a suture zone termed Terrane Boundary Shear Zone (TBSZ). BC in this area is largely composed of hornblende-biotite granite with some mafic dykes. Rocks in the TBSZ include quartzofeldspathic (leptynite) gneiss, garnet-orthopyroxene-bearing granitoid, mafic granulites (Group A of cratonic affinity, and Group B of EGP affinity), Mg-Al granulite and an isolated exposure of orthopyroxene-bearing gneiss. Detailed geochemical analysis shows remarkable similarity between Hbl-Bt granite and Grt-Opx-bearing granitoid, with A-type affinity, and between mafic dykes and Group A mafic granulites. However, the Opx-bearing gneiss is geochemically distinct having I-type affinity, similar to TTG gneisses described from BC. Metamorphic phase equilibria analysis and trace element modelling shows that (i) melting of Opx-bearing gneiss would produce a ferroan granitic melt resembling the Hbl-Bt granite, (ii) metamorphism at appropriate P-T conditions would convert the granite to Grt-Opx-bearing granitoid and the mafic dyke to Group A mafic granulite. U-Pb geochronology of zircon constrains emplacement ages of the magmatic precursors of Opx-bearing gneiss and Grt-Opx-bearing granitoid as ca. 2.73 and ca. 2.5 Ga, respectively. These rocks were subjected to an early granulite facies metamorphism, followed by an amphibolite facies metamorphism, shearing and hydrous fluid flux. Geochronological data shows that the latter event took place at ca. 0.52 Ga, while the earlier granulite facies event can only be tentatively suggested to be of late Stenian/Tonian age. Collating all the evidence (including published geophysical and geochronological data), we suggest that the initial collision between BC and EGP took place during late Stenian/Tonian time as a consequence of formation of the Greater Indian Landmass, a part of Rodinia supercontinent. The TBSZ, probably initiated due to the late Stenian/Tonian collision, was reactivated and reworked tectonothermally at ca. 0.52 Ga, caused by far field stress effect of the Kuunga orogeny, which was strong enough to obliterate most of the imprints of the late Stenian/Tonian orogeny.

Geochemistry and zircon U-Pb geochronology of diabase, gabbros, anorthosites and ultramafic rocks from the South Andaman Island ophiolite suite on the South-Eastern margin of the Indian plate

The South Andaman Island is located on the southeastern fringe of the Indian plate. Tectonically, the islands are associated with the Burma-Andaman-Java subduction zone of the Indonesian arc system. Several Mesozoic ophiolitic exposures are preserved in the Burma-Andaman-Java subduction zone. This study reports new petrological, geochemical, and zircon U-Pb geochronological data for the diabases, gabbros, anorthosites, and harzburgites of the ophiolite suite exposed in the South Andaman Island. The diabase is composed of plagioclase and clinopyroxene exhibiting ophitic textures. Anorthosite mainly contains plagioclase (> 90 vol%), pyroxenes and Fe-oxides. Gabbros have plagioclases and clinopyroxenes occurring as the major minerals. Harzburgites consist of olivine, orthopyroxene, clinopyroxene and spinel. Geochemical investigation reveals that these rocks were formed through crystal fractionation and accumulation processes from partial melt generated from a spinel lherzolite source. Gabbro, diabase, and anorthosites are formed by ∼1–5% melting of the, whereas depleted harzburgites were generated by slightly higher degrees of partial melting (∼10–20%) of the spinel lherzolite source. Trace element modelling reveals the influence of mid-oceanic ridge basalt (MORB)-melt interaction and subduction-related mantle alterations.The anorthosite and gabbro samples yield zircon 206Pb / 238U weighted mean ages of 93.2 ± 1.0 Ma and 442.5 ± 4.0 Ma, 383.9 ± 3.3 Ma respectively. However, the harzburgite sample records multiple age population including Paleo-Mesoproterozoic grains ranging in age from 1216 to 1931 Ma as well as Neoproterozoic and Phanerozoic grains of 607 Ma, 594 Ma, 237 Ma, 105 Ma, and 85 Ma suggesting a long-lived supra subduction zone in the Burma-Andaman-Java subduction zone where the South Andaman ophiolite suite was formed, evolved, and was emplaced.

On the occurrence of kalsilite in melilite-bearing ultrapotassic lavas from the Roman Province (Vulsini Mts., central Italy)

The Montefiascone Volcanic Complex belongs to the Vulsini Mts. volcanic district, one of the several leucite-bearing ultrapotassic Pleistocene volcanoes of the Roman Province (Central Italy). It is characterized by abundant pyroclastic successions and minor lava flows emplaced during a caldera formation phase occurred ∼0.3–0.2 Ma. The volcanic products span from leucite-basanites and leucite-tephrites to rarer leucitites and melilitites, occasionally associated with very rare kalsilite melilitolite ejecta. Here we report for the first time the presence of primary kalsilite in two lava flows in Forcinella and Feltricci localities, identifying a kamafugitic component among the Montefiascone volcanics.Detailed petrographic, SEM and EPMA studies on ten samples from the two localities revealed, in addition to kalsilite, variable amounts of clinopyroxene, melilite, leucite and nepheline plus rarer olivine coupled with carbonate-bearing segregation pockets in the groundmass. Geochemical and isotopic data confirm the exotic nature of the analysed rocks, showing ultrabasic/basic (SiO2 = 42.15–45.94 wt%) and ultrapotassic (K2O = 6.17–8.48 wt%; K2O/Na2O = 5.3–8.7) compositions, as well as strongly radiogenic 87Sr/86Sr (0.71038–0.711068), poorly radiogenic 143Nd/144Nd (0.512057.512091) and uniform middle radiogenic 206Pb/204Pb (18.76–18.78). Boron isotopes (δ11B ranges from 7.36 to 8.82‰) almost overlaps the typical MORB-OIB range values.These geochemical features recall those of the nearby kamafugite rocks from the Intra-Apennine Province. Major oxide and trace element modelling support fractionation links among the Italian kamafugite types. The overall geochemical and isotopic signatures speak for a common strongly subduction-metasomatized and orthopyroxene-free lithospheric mantle source. As already proposed for some nearby kamafugite centres in central Italy, limestone assimilation during magma ascent is a viable option to explain the mineralogical and whole-rock composition of Montefiascone lavas (i.e., the presence of carbonate plagues in the lava groundmass). This is particularly relevant considering the thick sedimentary sequence pierced by the limited volume of the upwelling magma. On the other hand, many other key features (such as the general lack of crustal xenoliths, the overall incompatible element fractionation and the radiogenic Sr isotopic ratios) clearly require also derivation from subduction-modified mantle sources.

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.

Post-tholeiite rifting and the genesis of ultramafic lamprophyres around 65 Ma in the Deccan large igneous province: Implications for the separation of India and Seychelles

This study examines the mineral and whole-rock compositions and paleomagnetic dating of lamprophyres from the West Coast Alkaline Complex (WCAC) within the Deccan Large Igneous Province. We aim to elucidate their connection with the Réunion plume, Deccan tholeiite magmatism, the onset of rifting and separation of Seychelles, and the role of lithospheric thinning.WCAC lamprophyres contain olivine and phlogopite phenocrysts/macrocrysts set in a groundmass of carbonate, clinopyroxene, nepheline, spinel, and melilite. They are characterized by Ti-Al phlogopite and Al-enriched spinels with high Fe2+/(Fe + Mg) ratios. They are undersaturated in silica and rich in MgO, TiO2, and LREEs, resembling the global ultramafic lamprophyres (UMLs).They originated from an enriched garnet lherzolite mantle, metasomatized by silicate and carbonate veins. Trace-element modeling suggests the derivation of melt by ∼1% partial melting of phlogopite-bearing garnet lherzolite mantle, corresponding to moderate pressure depths (3–4 GPa) and lithospheric mantle thickness of 90–100 km. Paleomagnetic investigations corroborate their intrusions, primarily during chron C29n, approximately at 65 Ma, which coincided with amplified plate velocities in the Indian Ocean.Seismic tomography indicates a current lithospheric thickness of about 50 km under the western Indian margin, suggesting a lithospheric delamination of about 40–50 km following the intrusion of lamprophyres at ca. 65 Ma. The initiation of passive rifts in the region, culminating in the separation of the Indian subcontinent and the Seychelles, prompted the UML magmatism. Lithospheric thinning continues along with continental rifting in response to greater plate-tectonic stresses within regions of persistent lithospheric weakness.

Revisiting the fluid regime during anatexis of metasedimentary rocks: New constraints from an integrated approach on the Cenozoic Himalayan granites

Water plays a critical role in the formation of granitic magmas and continental crust, but distinguishing water-present and water-absent anatectic scenarios using the geochemistry of granites is controversial. In this study we use an integrated approach that combines whole-rock major and trace element geochemistry, Ti-in-zircon thermometry, phase equilibrium modeling and trace element modeling to study the water regime that produced the Miocene granites from the Malashan-Gyirong area in central Himalaya to solve this controversy. The Gyirong granites have relatively low CaO, Sr and Ba and high Rb and (87Sr/86Sr)i and the Malashan granites have relatively high CaO, Sr and Ba and low Rb and (87Sr/86Sr)i, and were classified as group A and group B granites, respectively. Following previous interpretation, group A and group B granites are consistent with the products of water-absent and water-present melting of metasedimentary rocks, respectively. Ti-in-zircon thermometry yielded maximum values of 761–796 °C for the Gyirong granites and 730–764 °C for the Malashan granites. Distinct variation trends in zircon trace element compositions indicate that these two groups of granites were not linked by crystallization differentiation. Using average compositions of Proterozoic pelite as starting materials, phase equilibrium modeling was carried out at a variety of P–T–H2O conditions typical of the Himalayan orogen, P = 5, 10 and 15 kbar, T = 600–800 °C, H2O = 0–10 wt%. Pelite can produce melts with coupled CaO–Na2O contents for both groups at 10 kbar. Specifically, constraints from compositions and temperatures require that the bulk H2O content is ca. 1–2 wt% for group A granites and > ca. 4 wt% for group B granites. Compared with the maximum structural water content of the pelite at 10 kbar (1.77 wt%), this study testifies that group A granites formed under water-absent conditions and group B granites under water-present melting conditions. Modeling shows that water-present melting can produce melts with high Sr–Ba and low Rb contents resembling group B granites, while water-absent melting can produce melts with low Sr–Ba and high Rb contents resembling group A granites. This study highlights that water can indeed cause differences in granite geochemistry but a comprehensive investigation is required to better determine the role of water during crustal anatexis.

Sub-arc Mantle Heterogeneity of the Northern Luzon Volcanic Arc: Mineral and Whole Rock Compositional Variability in Mantle Xenoliths from Lutao Island

Abstract Mantle xenoliths hosted in volcanic rocks from the island of Lutao offer a glimpse into the nature of the mantle beneath the northern Luzon volcanic arc. The xenoliths are spinel-bearing and composed mostly of harzburgite with one lherzolite and one olivine orthopyroxenite. The olivine (Fo92.5–88.9), orthopyroxene (Mg# = 94.6–89.2), and clinopyroxene (Wo49.1–38.1En57.0–45.4Fs3.0–11.0) compositions are similar to those of abyssal peridotites. The spinel compositions are variable and can be principally divided into high-Al (Cr# < 45) and low-Al (Cr# > 45) groupings. The whole rock compositions are similar to abyssal peridotite (Al2O3 = 0.95–2.07 wt %; Mg# = 88.5–90.9) and have U-shaped chondrite normalized rare earth element patterns. The Sr-Nd isotopes of the xenoliths are broadly chondritic (87Sr/86Sri = 0.704400–0.707908; εNd(t) = 0.0 − +1.5). The two-pyroxene equilibrium temperatures range from 900 to 1200 °C with the majority of temperature estimates >1000 °C. The olivine-orthopyroxene-spinel oxygen barometry estimates yielded ΔFMQ values from 0 to +2 and correspond to moderately oxidizing to oxidizing conditions. The xenoliths are likely derived from the Philippine Sea Plate lithospheric mantle that was modified by melt extraction and/or fluid enrichment processes. Trace element and isotopic mixing modeling indicate that 1–2% contamination by subducted South China Sea sediment can explain the Sr-Nd isotopic enrichment and Th and U elemental variability within the xenoliths assuming an initial composition similar to enriched depleted mid-ocean ridge mantle (E-DMM). The anomalously high two-pyroxene equilibrium temperatures of the Lutao xenoliths relative to other regions of the northern Luzon volcanic arc (Iraya <1000 °C) indicate that they were affected by a high-temperature event that was likely a consequence of recent intra-arc rifting that occurred after collision (<6 Ma) between the Luzon arc and the Eurasian margin.

Petrogenetic and geochemical constraints on ca. 1.89–1.88 Ga Bastanar mafic dyke swarm, Bastar craton, India: Insights into MORB- and OIB-type contributions and interactions with metasomatized subcontinental lithospheric mantle

The Bastar craton of the Indian Shield hosts several generations of mafic dyke swarms of various trends, compositions, and ages, which span from ca. 2.7 Ga to 1.42 Ga. This study focuses on geochemical attributes of the ca. 1.89–1.88 Ga NNW-trending Bastanar swarm, aiming to address a perceived discrepancy between its arc-like geochemistry and the influence of a heterogeneous mantle source. To resolve the intra-swarm geochemical variations, we conducted a comprehensive geochemical characterization and petrogenetic interpretation of the ca. 1.89–1.88 Ga Bastanar swarm. The samples from this swarm are categorized into two distinct groups, primarily based on their geochemical composition. The Group 1 samples exhibit higher TiO2 (1.06–1.86 wt%), (La/Yb)N (7.2–8.6), (Gd/Yb)N (2- 2.27), Nb (14.6–16.6 ppm), Th (1.23–3.03 ppm) and Zr (104–118.72 ppm) concentrations than the Group 2 samples. Furthermore, rare-earth element patterns and variations in high-field strength element contents in the Group 1 samples suggest derivation from a deeper, less depleted mantle source resembling an OIB/less MORB-type. This inference is further supported by higher TiO2/Yb, Zr/Nb, and Nb/Y ratios. In contrast, the Group 2 samples indicate derivation from a shallower, more MORB/less OIB-type depleted mantle source, as evidenced by lower TiO2/Yb, Zr/Nb, and Nb/Y ratios. Variations in Dy/Yb and Gd/Yb ratios confirm the involvement of variable mantle sources, implying the derivation of the Group 1 and 2 samples from garnet-rich and spinel-rich lherzolite mantle sources, respectively. The absence of consistent negative Nb-Ta-Ti anomalies in the Group 1 samples suggests an uncontaminated nature, ruling out any role of crustal contamination. On the other hand, the Group 2 samples display negative Nb-Ta-Ti anomalies with enriched LREE and LILE patterns, indicating the involvement of crustal components in their genesis. A trace-element modelling suggests that the ca. 1.89–1.88 Ga mafic dyke swarm exhibits significant intra-swarm variability, with at least two distinct source components contributing to its genesis – a depleted MORB-type and an enriched OIB-type mantle. Notably, the geochemical characteristics of the Group 2 samples suggest interaction with a metasomatized mantle source, possibly enriched by fluids from an earlier subducted slab. Geochemical evidence presented in this work supports Archean subduction-related processes for the crustal growth of the Bastar craton and highlight the enduring influence of a metasomatized sub-continental lithospheric mantle on subsequent magmatism over millions of years.

Petrogenesis of Mesozoic pegmatites in the Dahongliutan Li-mineralized belt (Western Kunlun, NW China)

The genesis of pegmatite and the associated Li-mineralization mechanism remain debated. The intimate temporal-spatial relationship between the pegmatites in the Dahongliutan Li-mineralized belt (incl. #509 Daoban and Akshay Li-poor pegmatites, and Akshay Li-rich pegmatite) and the Dahongliutan granite pluton allows us to discuss the origin of pegmatite and the Li-mineralization mechanism by conducting detailed geochemical and geochronological investigations on the two rock units in the Western Kunlun orogen (NW China). The Li-poor pegmatites in the Dahongliutan Li-mineralized belt exhibit similar or much lower Li, Rb, and Cs compared to the adjacent granite, and their εNd(t) values (-8.7–-10.0) are significantly lower than those of the granite (-7.5–-8.1) as well. However, the Akshay Li-rich pegmatite shows extremely high Li, Rb, Nb, Ta, and initial Sr isotopic ratios (0.7223–0.7227, obtained from apatite), in contrast to the Dahongliutan granite (ISr = 0.7108–0.7157). The whole-rock δ7Li values of the Dahongliutan granite, pegmatites, and the wall-rock Bayanharshan Group metapelite are + 0.76–+5.67 ‰, +0.8–+14.4 ‰, and −4.2–+2.8 ‰, respectively. Zircon and cassiterite U-Pb dating yielded 206Pb/238U ages of 216 ± 2 Ma and 195 ± 2 Ma for the #509 Daoban Li-poor pegmatite and Akshay Li-rich pegmatite, respectively. Modeling of trace elements and Li isotopes suggests that the pegmatites in the Dahongliutan Li-mineralized belt did not originate from fractional crystallization of granitic magma. Instead, they may have formed through direct anatexis under amphibolite-facies conditions. Abundant Li-rich claystone, occasionally accompanied by evaporite, may exist in the magma source of the Li-rich pegmatite, which is the key factor controlling pegmatitic Li mineralization.

Rapid switch in geodynamic setting revealed by episodic granitoid magmatism in the Lynn Lake greenstone belt of Paleoproterozoic Trans-Hudson Orogen, Manitoba, Canada

Bedrock geological mapping indicates that diverse granitoid intrusions occur closely in space and time in the Paleoproterozoic Lynn Lake greenstone belt of the Trans-Hudson Orogen, Manitoba, Canada. These intrusions display distinct geochemical characteristics of I-type, A-type and adakite-like granitoids as revealed by their respective alkalinity, alumina saturation state, and trace element concentrations and ratios. The ca. 1879 Ma I-type granitoid pluton intruded the 1892 to 1884 Ma supracrustal package of volcanic, volcaniclastic and minor sedimentary rocks that constitute the greenstone belt, which was cut by 1873 Ma A-type granite and then by 1854 Ma adakite-like quartz diorite intrusions. Why the geochemical affinities of these granitoids change so rapidly (within ∼25 Myr) in orogenic belts, however, is not as well understood. Here we propose using SmNd isotope data and trace element modeling that the compositional variety of the granitoid intrusions likely reflects the mineral assemblages of their source rocks rather than magmatic fractionation. They may have been formed by partial melting of the disparate sources at different depths to respond to the change in tectonic regimes. The timing of the granitoid intrusions tracks a rapid geodynamic transition from subduction to intra-arc extension and to slab break-off due to terrane accretion and collision. This records a crucial piece of history in geodynamic evolution of the Trans-Hudson Orogen, which to date has been overlooked largely in many other orogens.

Lithium isotope fractionation in highly evolved granites and associated magmatic-hydrothermal process: the case of the Xintianling granites in the Nanling Range, South China

ABSTRACT The highly evolved granites are closely associated with rare-metal mineralization. It is still debated whether the enrichment of rare metals is mainly controlled by magmatic or hydrothermal processes. Lithium is a fluid-mobile element that preferentially partitions to the fluid phase when the melt reaches H2O saturation, making the Li isotope system a powerful geochemical tool to record the potential melt-fluid interaction during magma evolution. This study addresses the issue of Li isotope fractionation in highly evolved fine syenogranites (FSGs) from the Xintianling pluton in the Nanling Range, South China. The Xintianling pluton consists of the coarse biotite monzogranites (CBMGs) and evolved FSGs. The identical Nd isotope data of the CBMGs (εNd(t) = -8.2–-8.8) and FSGs (εNd(t) = -7.9–-8.9) suggest that they were derived from the same crustal magma. The FSGs were derived from the highly evolved melts extracted from CBMGs through the fractionation of major rock-forming minerals. This is supported by the compositional shifts and validated by Rb, Sr, and Ba trace element modelling. Li isotopic analysis shows that highly evolved FSGs have heavier Li isotope compositions (δ7Li = 2.6–7.9 ‰) than the CBMGs (δ7Li = 0.1–1.9 ‰). Isotopic modelling, coupled with variations in Li concentrations in FSGs, suggests that diffusive kinetic fractionation by disequilibrium fluid exsolution during the late stage of the granitic magma evolution controls the Li isotope fractionation of the Xintianling granites. Rapid fluid exsolution not only modified the Li isotope composition of the FSGs but also transported the rare metals (e.g. W) and fluxing components (Li, F, and Cl), which resulted in the formation of the large Xintianling W deposits.

Diamonds Discovered in the Forearc Harzburgites Hint at the Deep Mantle Source of the Skenderbeu Massif, Western Mirdita Ophiolite

The ultra-deep genesis of ophiolitic peridotite has reshaped our perception of the genesis of the oceanic mantle. Although ultra-high pressure (UHP) mineral assemblages have been unearthed in dozens of ophiolites in different orogenic belts around the world, the vast majority of them have been limited to podiform chromitites formed in suprasubduction zone (SSZ) settings, leaving uncertainty about whether such UHP minerals are intrinsic to the oceanic mantle or influenced by a specific mantle rock type. Here, we report on the occurrence of diamonds recovered from the harzburgites within the Skenderbeu massif, Mirdita ophiolite. The whole-rock, mineralogical major and trace element compositions, and redox states of the harzburgites align with modern abyssal harzburgites. Trace element modeling of clinopyroxene indicates that harzburgites have endured varying degrees of garnet-facies melting (~2%–5%) before progressing to spinel-facies melting (~10%–12%). Mineralogical characteristics further support that the Skenderbeu harzburgites underwent late-period MORB-like melt metasomatism in a forearc spreading center. An unusual mineral assemblage of diamonds has been separated from the studied harzburgites. The first occurrence of ophiolite-hosted diamonds discovered in the forearc harzburgites, together with previous similar discoveries in the SSZ ophiolitic chromitites, suggest that the ophiolite-hosted diamonds are not specific to certain mantle rocks.

Correction to: The petrogenesis of highly fractionated gem‑bearing pegmatites of Malawi: evidence from mica and tourmaline chemistry and finite step trace element modelling

Correction to: The petrogenesis of highly fractionated gem‑bearing pegmatites of Malawi: evidence from mica and tourmaline chemistry and finite step trace element modelling

Review of temporal trend in Archean granitoid chemistry- implications for Eoarchean (3.7–3.9 Ga) transition in geodynamic regime of crust formation

In this study, we review changes in Archean granitoid chemistry using stacked probability density estimation and statistical change-point analyses. A significant change in Archean sodic granitoid composition occurs in the 3.9–3.75 Ga period, marked by a shift of SiO2, Al2O3, Na2O, V, Cr, Sr, Sr/Y to higher, and TiO2, FeOT, MnO, P2O5, A/NK, A/CNK, Sc, Ba, Ta, Nb, Hf, Y, and rare earth elements to lower values, which reflects different formation P-T regime and petrogenetic processes. Using phase equilibria and trace element modelling, we demonstrate that the pre-3.9 Ga Acasta Tonalite Gneiss (ATG)-like granitoids formed at shallow depths (∼3 km) but are chemically distinct from impactites of all known large impact craters with basaltic/mixed target rocks. Impact melting produces MgO-rich melts at relatively greater “apparent” depths, reflecting target rock chemistry and contribution from the upwelling mantle, and is unlikely to have produced the shallow-seated ATG rocks. The ATG trace element chemistry is similar to many modern-day non-arc settings felsic rocks, which suggests that the melting regime in these extensional settings also existed before 3.9 Ga and ATG formed in a mantle plume-driven extensional tectonic regime involving melting of long-lived mafic proto-crust at shallow crustal depths above mantle upwelling. In contrast, the 3.9–3.75 Ga transitional and post-3.75 Ga TTGs formed by deep-seated (25–50 km) melting in equilibrium with garnet-bearing amphibolite. When considered together with Ti isotopic evidence of shift from tholeiitic to calc-alkaline magmatism in the 4.0–3.75 Ga period, the evidence for a transition from a regime dominated by crustal reworking to one of significant juvenile input at 3.60–3.85 Ga, noted in trace element/Hf-isotope composition of detrital zircon from several cratons and the reported ultra-high pressure metamorphism/interleaving arc-plume sequence, require a greater depth of melting. We interpret this transition to mark a gradual shift from stagnant-lid to intermittent mobile-lid-type plate tectonics involving deeper subduction of oceanic lithosphere. However, from the single occurrence of pre-3.9 Ga ATG rocks and their similarity with some Archean Fe-rich silica-saturated rocks chemistry, it is unknown whether plume magmatism was the only active mechanism or operated in conjunction with mobile-lid tectonics on early Earth, like in Archean.

The petrogenesis of highly fractionated gem-bearing pegmatites of Malawi: evidence from mica and tourmaline chemistry and finite step trace element modelling

The petrogenesis of highly fractionated gem-bearing pegmatites of Malawi: evidence from mica and tourmaline chemistry and finite step trace element modelling

A Log‐Ratio‐Based Algorithm for Petrologic Mass‐Balance Problems and Uncertainty Assessment

AbstractWe provide a new algorithm for mass‐balance calculations in petrology and geochemistry based on the log‐ratio approach championed initially by John Aitchison (e.g., Aitchison, 1982, https://doi.org/10.1111/j.2517-6161.1982.tb01195.x; Aitchison, 1984, https://doi.org/10.1007/bf01029316) along with the underlying principles, mathematical frameworks, and data requirements. Log‐ratio Inversion of Mixed End‐members (LIME) is written in MATLAB and calculates phase proportions in an experiment or rock given a bulk composition, the composition of each phase, and the associated compositional uncertainties. An important advantage of LIME is that performing the mass‐balance calculation in inverse log‐ratio space constrains phase proportions to be between 0 and 100 wt.%. Further, the resulting LIME phase proportions provide realistic estimates of uncertainty regardless of data distribution. These two characteristics of LIME improve upon standard multiple linear regression techniques, which may yield negative values for phase proportions if non‐constrained or report oversimplified symmetric errors. Primary applications of LIME include estimating phase abundances, calculating melting and metamorphic reaction stoichiometries, and checking for open system behavior in phase equilibria experiments. The technique presented here covers whole‐rock analysis, mineralogy, and phase abundance, but could be extended to isotopic tracers, trace element modeling, and regolith component un‐mixing. We highlight the importance of uncertainty estimations for phase abundances to the fields of petrology and geochemistry by comparing our results from LIME to previously published mass‐balance calculations. Furthermore, we present case studies that demonstrate the role of mass‐balance calculations in determining magma crystallinity and defining melting reactions.

Geochemical variations of anatectic melts in response to changes of P–T–H2O conditions: Implication for the relationship between dehydration and hydration melting in the Himalayan orogen

The anatectic mechanism is the most fundamental criterion in determining the petrogenetic relationship between migmatite, granulite and granite at convergent plate margins. Although two main anatectic mechanisms, namely water-absent (dehydration) and water-fluxed (hydration) melting, have been identified in many collisional orogens, their spatiotemporal relationship and genetic connection remain obscure. To address this issue, we conduct an integrated study of forward phase equilibrium and trace element modelling for partial melting of metapelite in the Himalayan orogen. The results are used to decipher changes of phase relation and melt composition in response to changes of temperature, pressure and water content. Pressure and water content mainly affect the solidus and evolving residual mineral assemblages, whereas temperature controls the anatectic reaction. Hydration melting at lower pressure and temperature would produce anatectic melts with relatively high SiO2, FeOT + MgO, CaO, K2O, and Ba contents and SiO2/Al2O3 ratios, but low Al2O3 and Na2O contents and Na2O/K2O and Rb/Sr ratios. In contrast, dehydration melting at higher pressure and temperature would generate anatectic melts with relatively high Al2O3 and Na2O contents and Na2O/K2O and Rb/Sr ratios, but low SiO2, FeOT + MgO, CaO, K2O, and Ba contents and SiO2/Al2O3 ratios. A comparison of geochemical compositions between the modelled melts and the two groups of well-characterized leucogranites in the Himalayan orogen indicates that the coeval dehydration and hydration melting would likely occur at the deeper and shallower crustal levels, respectively. In combination with the secular evolution of metamorphic facies series in the Himalayan orogen, we propose a spatiotemporally coupled dehydration-hydration melting mechanism. The early metamorphism at lower geothermal gradients would result in metamorphic dehydration in the lower crust and hydration in the upper crust during the syn-collisional period in the Early Cenozoic. As a result, the later metamorphism at higher geothermal gradients would lead to dehydration melting in the lower crust and hydration melting in the upper crust during the post-collisional period in the Late Cenozoic. This study highlights the applicability of phase equilibrium and trace element modelling to reveal the geochemical variations of anatectic melts in response to changes of P–T–H2O conditions. This provides new insights into the relationship between dehydration and hydration melting in collisional orogens.