Slab-derived Fluids Research Articles

Continental Arc Plutonism in a Juvenile Crust: The Neoproterozoic Metagabbro-Diorite Complexes of Sinai, Northern Arabian-Nubian Shield

Based on new field, petrographic, and whole-rock geochemistry data, we investigated three discrete metagabbro-diorite complexes (MGDC) across the E-W Sinai to contribute to increasing knowledge of the evolution of the juvenile continental crust of the Neoproterozoic Arabian–Nubian Shield. The three MGDCs vary in the dominance of the gabbroic versus dioritic rock types among each of them. Gabbroids are distinguished into pyroxene-hornblende gabbros and hornblende gabbros, whereas dioritic rocks have been subdivided into diorites and quartz diorites. The studied MGDC rocks are almost metaluminous and possess prevalent calc-alkaline characteristics over subsidiary tholeiitic and alkaline affinities. The most distinctive feature in the profiles of the investigated MGDCs on the N-MORB-normalized spider diagrams is the coincidence of stout negative Nb anomalies and projecting positive Pb spikes, which is typical of igneous rocks evolved in subduction zones. The three MGDC samples exhibit variably LREE-enriched patterns [(La/Yb)N = 4.92–18.55; av. = 9.04], either lacking or possessing weak to negligible positive and negative Eu anomalies. The calculated apatite and zircon crystallization temperatures reveal the earlier separation of apatite at higher temperatures, with the obvious possibility of two genetic types of apatite and zircon in the magma (cognate vs. xenocrystic) since both accessories have yielded very wide ranges of crystallization temperatures. The investigated MGDCs were formed in a continental arc setting, particularly a thick-crust arc (>39 km). The parent magmas comprised components derived from the melting of the mantle wedge, subducting oceanic lithosphere, and subducting overlying sediments. The mantle input was from a spinel–garnet transitional mantle source at a depth of ca. 75–90 km. The impact of slab-derived fluids was much greater than that of slab-derived melts, and so subduction-related fluids had a crucial effect on metasomatizing the partially melted mantle source. The parent mantle-derived magma has been subjected to substantial crustal contamination as a dominant mechanism of differentiation.

Variable oxidizing capacity of slab-derived fluids: Insights from Fe and S speciation in glasses from the Troodos Ophiolite

Oxygen fugacity (ƒO2) varies systematically across tectonic environments. The typically high ƒO2 recorded in arc settings relative to oceanic ridges is attributed to recycling of oxidized materials derived from subducting slabs into regions of the mantle that undergo partial melting. To evaluate further the relationship between ƒO2 and mantle metasomatism, Fe and S X-ray Absorption Near-Edge Structure measurements were conducted on volcanic glass wafers sampled across the extrusive section of the Troodos Ophiolite, Cyprus, which is a type locality for studying the influence of subduction on mantle melting processes and the generation of oceanic crust. The glasses record ƒO2 values relative to the quartz-fayalite-magnetite (FMQ) buffer of FMQ+0.13±0.16(1σ) to FMQ+0.74±0.23(1σ) upon quenching at the seafloor. At a given MgO content, the ƒO2 values recorded by the Troodos glasses do not vary systematically based on type (i.e., boninitic versus tholeiitic) nor by sampling location, indicating that the lavas were derived from primary melts and mantle sources that were indistinguishable in their initial ƒO2 levels. The glass ƒO2 values do not vary with dissolved H2O contents, nor with trace element ratios (e.g., Ba/La, Ce/Pb) and Pb and Sr isotopic compositions, all used to monitor interaction of slab-derived materials with the Troodos mantle source. The decoupling of ƒO2 and H2O and fluid-mobile elements renders the Troodos an important endmember for interpreting global variations in the redox state of mantle melts formed in subduction-influenced settings. For example, the Troodos glasses are reduced relative to Izu-Bonin boninites and melt inclusions from Mariana and Cascades arc volcanoes, despite overlapping and elevated H2O contents, Ba/La, and S/Dy relative to mid-ocean ridge basalts. The lack of correlation between ƒO2 and volatile contents and with proxies for subduction influence in the Troodos glasses is similar to basalts formed in back-arc environments and in near-trench settings in Izu-Bonin-Mariana. Such comparisons suggest that the slab-derived fluids that infiltrated the Troodos melting region were derived from a shallow slab and had low oxidizing capacity relative to fluids or melts which interacted with the mantle in the Izu-Bonin-Mariana arc and in the Cascades, which represent arc environments with comparable glass Fe XANES measurements. A simple framework accounts for local and global variability in ƒO2 in fractionation-corrected glasses, where fluids of variable oxidizing capacity interact with depleted mantle to produce hydrous melts. The inheritance of volatile and fluid-mobile elements in regions of the mantle that undergo partial melting does not necessarily lead to higher ƒO2 if the oxidizing capacity of the infiltrating hydrous fluids or melts is low. In back-arcs and near-trench environments, partial melts formed by interaction with fluids of low oxidizing capacity are similar to arc melts in their volatile contents, trace element, and isotopic signatures, but modestly oxidized relative to mid-ocean ridge basalts.

Systematics of slab-derived components in Central Java, Sunda Arc: Evidence for differential material transfer across the Southeast Asian convergent margin

The modern volcanoes of Indonesia, represent the cycling of material and heat along the Southeast (SE) Asian subduction zone. This volcanic activity is well exemplified in the middle section of Java Island, which is found on the Sunda arc portion of the SE Asian convergent system. Central Java's volcanoes inherit their geochemistry from silicate melts that formed due to hydrous melting at a range of depths above their subducting slab: ∼136 to ∼272 km. It has been shown that across-arc variations in magmatic trace elements like boron (B) can shed light on the quantity and type of slab-derived fluid released in a forearc to back-arc sequence. In this study, we present major, trace and radiogenic chemistry data for five volcanoes arranged in a roughly N-S lineament across Central Java. We show that the observed chemical changes across the section of the Sunda arc can be explained by variations in different kinds of source from fluid mobile material being devolatilized across the arc. Our mixing models suggest that the magmas underlying these volcanoes see a decreasing ratio of sediment to altered oceanic crust (AOC) contributions from the volcanic front (Merapi volcano) toward the rear-arc (Muria). We estimate that a maximum ∼1% contribution of sediment-derived fluid is sufficient to explain the chemical systematics of this arc segment. Shallow slab dehydration, mostly from the sediment cover, at low temperature (at <60 km of slab depth) transfers the fluid high in B/Nb, Ba/Nb, 87Sr/86Sr and 207Pb/204Pb to the mantle beneath this region's arc front, inducing flux melting. Farther from the trench, the signature of the subducting oceanic plate increases and peaks at around ∼160–170 km of slab depth, carrying substantial slab liquids with moderate to low B/Nb and Ba/Nb with lower 87Sr/86Sr and 207Pb/204Pb to the mantle beneath middle- and rear-arc region. Beyond 200 km, in the rear-arc, evidence for slab liquid is low, but not non-existent. Unexplained Ba-enrichment in magmas from across the arc may account for the role of subduction serpentinite dehydration in the subduction channel or local crustal assimilation. This study provides a conceptual model for material transfer from the subducting slab in SE Asia. Our findings have important implications for the understanding of geochemical cycling in this region, as well as the impact of these cycles on present and past volcanism.

A late Neoarchean continental arc system in the northeastern North China Craton: constrain from the Guanghua mélange in the Tonghua area, southern Jilin Province

ABSTRACT The late Neoarchean tonalite – trondhjemite – granodiorite (TTG) and mafic rocks record essential information regarding the early tectonic evolution of the North China Craton, which is crucial for understanding the onset and development of the plate tectonic regime. In this study, we present a comprehensive analysis of petrography, geochronology, and geochemistry of the Guanghua Mélange from the Tonghua area in NE China. The mélange comprises basalt, gabbro, pyroxene diorite, and trondhjemite, and is classified into three groups: basaltic – gabbroic rocks (Group #1), dioritic rocks (Group #2), and trondhjemites (Group #3). Zircon U-Pb analyses indicate that they formed at 2596–2550 Ma, 2531–2515 Ma, and 2502 Ma, respectively, with metamorphism occurring during ~ 2500–2450 Ma. Meanwhile, geochemical data indicate that the calc – alkaline basaltic rocks (Group #1) exhibit slight depletion of Nb or Ta with weakly negative Eu anomalies, suggesting derivation from the lithospheric mantle and formation within a continental arc environment. The pyroxene diorites (Group #2) display similar geochemical characteristics to Group #1, likely formed by partial melting of depleted mantle that underwent sgignificant metasomatism due to variable amounts of slab-derived fluids or melts. The tholeiitic granitoids (Group #3) exhibit strongly fractionated chondrite-normalized REE patterns with high (La/Yb)N and Sr contents, low Y contents, and depletion of Nb, Ta or Ti, resembling adakites generated in arc setting, indicating likely formation through partial melting of subducted oceanic crust. This process was controlled by the initial plate tectonics under a higher geothermal gradient. Consequently, when combined with previous findings, we conclude that an initial subduction of oceanic plate operated in the northeastern North China Craton during the late Neoarchean (>2.55 Ga). Subsequently, it transformed into modern-type plate tectonics during ~ 2.55–2.53 Ga, eventually followed by the continent-continent collision at ~ 2.50–2.45 Ga, resulting in the formation of the Ji – Liao – Ji (Jilin – Liaoning – Hebei) Collision Orogenic Belt.

Two-stage crust-mantle interactions from oceanic subduction to post-collisional extension in the northern margin of the North China Craton: Insights from Paleozoic to Mesozoic magmatism

The northern margin of the North China Craton experienced prolonged tectono-magmatic evolution during the late Paleozoic−early Mesozoic in response to the southward subduction and closure of the Paleo-Asian Ocean. However, details about the subduction process and the timing of the tectonic transition from subduction to post-collision are still poorly constrained. Here, we identify two-stage crust-mantle interactions in the Wulashan area and report new geochronology, geochemistry, and Sr-Nd-Pb-Hf isotopic data for magmatic rocks that record such processes following the subduction and closure of the Paleo-Asian Ocean. The early Carboniferous Xiguanjing pluton features a bimodal suite of gabbro (ca. 333 Ma) and syenogranite (ca. 331 Ma). The gabbros have arc-like geochemical affinities, with low Nb/La (0.31−0.40) and La/Ba (0.04−0.09) ratios, and variable Rb/Y (1.22−2.94) ratios, as well as enriched, mantle-like Sr-Nd-Pb (87Sr/86Sri = 0.7046−0.7047; εNd(t) = −3.8 to −3.5; 206Pb/204Pbi = 17.078−17.141) and enriched to depleted Hf (εHf(t) = −4.5 to +6.2) isotopic values. Such geochemical signatures indicate that they were derived from partial melting of the subcontinental lithospheric mantle that was metasomatized by slab-derived fluids, with minor involvement of asthenospheric components. In contrast, the contemporaneous syenogranites are characterized by lower negative εNd(t) (−13.5 to −12.1) and εHf(t) values (−16.3 to −8.2), which suggests that they were formed by partial melting of the lower crust. Late Triassic Shadegai and Xishadegai plutons are mainly composed of enclave-bearing syenogranite, and both mafic microgranular enclaves and syenogranites crystallized at ca. 233−231 Ma. The mafic microgranular enclaves have geochemical features similar to those of the early Carboniferous gabbros, and also have moderately enriched isotopic compositions (εNd(t) = −9.7 to −8.4; εHf(t) = −9.2 to −0.3), which suggests that they originated from interaction between mantle-derived magma and overlying crust-derived magma, with minor additions of asthenospheric melts in their sources. Field and petrological observations, coupled with the similar ages of the host granites and mafic microgranular enclaves, suggest a magmatic mingling process. Isotopic mixing models suggest that minor amounts (∼10%−20%) of lower crustal materials were mixed during the formation of the mafic microgranular enclaves. The host syenogranites display calc-alkaline to alkalic and metaluminous to weakly peraluminous compositions, and negative εNd(t) (−15.0 to −12.1) and εHf(t) values (−16.4 to −9.8), which indicates that they were mainly derived from partial melting of the lower crust and experienced the injection of deep mantle-derived magmas. Our new data, along with previously published data for magmatic rocks in the northern margin of the North China Craton, suggest that the early Carboniferous bimodal intrusive rocks formed in a localized back-arc extensional regime that was probably triggered by slab rollback of the Paleo-Asian Ocean. However, the Late Triassic plutons formed in a post-collisional extensional regime in response to slab breakoff or lithospheric delamination. Temporal variations of Nd-Hf isotopes for the magmatism in the northern margin of the North China Craton suggest that tectonic switching from advancing to retreating subduction to post-collisional extension occurred during the late Paleozoic to early Mesozoic. We propose that a tectonic transition from subduction to post-collisional extension may have occurred during the Early−Middle Triassic, marking the final closure of the Paleo-Asian Ocean, which most likely took place at ca. 250−235 Ma.

Constraints on crustal recycling from boron isotopes in Italian melt inclusions

Boron represents an important tracer of crustal recycling processes in subduction zones, because it is readily mobilised from the subducted lithosphere and different components in the slab are isotopically distinct. Profiles of boron content and isotope ratio across magmatic arcs generally show that B concentrations decrease with increasing slab depth, which implies decreasing amount of slab-derived fluids. To date, however, data on continental-collision zones and post-collisional subduction settings are scarce.This study examines Plio-Quaternary Italian magmatism to quantify crustal recycling in a complex subduction setting. Magmatic products vary from (ultra)potassic along the Tyrrhenian side in the north, to calcalkaline and Na-alkaline in the south.Combined major and trace element and [B] content and δ11B values are reported in 99 Melt Inclusions (MIs), analyses from a wide range of Italian lavas. [B] vary from 4 to 298 µg/g and δ11B from -29.2 to -3.9‰. The B isotopic values are considerably lower than previously reported in arcs and other post-collisional setting magmatism. We infer a role for phengite in the source of all studied Italian magmas (with the exception of Mt. Etna lavas). This white mica is stable to high pressures in subducted sediments of altered oceanic crust and records dehydration and 11B depletion due to dehydration processes.MIs hosted in highly fosteritic olivines (Fo >74; median of 89) from across Italy reveal that primary melts tap heterogeneous mantle including subducted oceanic and continental components that were introduced during the Alpine, and Adriatic and Ionian subduction phases.The combined geochemical data record the involvement of sediments that variably metasomatized the mantle wedge. We propose that slab detachment and consequent heat input from the inflow of hot asthenosphere was responsible for phengite breakdown in subducted sediments and locally produced metasomatism of the mantle wedge, imposing a characteristic B isotope signature to the overlying mantle. Continued heating due to asthenosphere inflow led to melting of the metasomatized mantle wedge and generation of the Italian magmatism. Mt. Etna represents an exception being dominated by asthenosphere upwelling through a slab window with minimal influence from active subduction.

Origin of the Mussau Trench in the Western Pacific: Geochemical and mineralogical constraints from basalts and serpentinized peridotites

The Mussau Trench is apparently an incipient subduction zone. However, petrologic and geochemical evidence for the nature and origin of the Mussau Trench is lacking. We present the first mineralogical and geochemical analyses of basalts and serpentinized peridotites from the Mussau Trench. All sampled Mussau Trench basalts yield MORB-like trace-element patterns and lack the geochemical characteristics of typical arc lavas. The Dupal-type isotopic compositions of the MORB-like lavas confirm that the overriding plate shows more affiliation to the Caroline Plate than to the Pacific Plate. The Mussau Trench could have formed along a pre-existing transform fault in the Caroline Plate. Peridotites from the Mussau Trench are highly serpentinized and can be divided into two groups. Group I peridotites have low MgO/SiO2 ratios and contain talc and quartz, suggesting that they were metasomatized by SiO2-rich fluids. Group II peridotites have higher Al2O3/SiO2 ratios at a given MgO/SiO2 ratio and rarely contain talc or quartz. The peridotites in the Mussau Trench differ from typical fore-arc peridotites in the Mariana Trench; for example, the Mussau peridotites are less depleted in trace elements and have lower spinel Cr# (<0.6) than those from the Mariana Trench, suggesting that they originated in the lithosphere of the Caroline plate. The temperatures (<688 °C) and pressures (<7.3 kbar) estimated using amphibole compositions (dominantly edenite and pargasite) suggest that the peridotites were serpentinized by high-temperature slab-derived fluids. Although the limited plate convergence at the Mussau Trench was insufficient to cause arc magmatism, subducted slab-derived fluids caused serpentinization of the overriding lithospheric mantle in an infant subduction channel. Therefore, the serpentinized peridotites and basalts in the Mussau Trench provide important information on the marks and processes of intra-oceanic subduction initiation.

Chemostratigraphy of the extrusive sequence of a Late Cretaceous Neotethyan ophiolite in southern Iran and its significance for the mode and tempo of subduction initiation magmatism and melt evolution in forearc tectonic settings

Chemostratigraphy of lavas in extrusive sequences of modern and ancient forearc environments provides significant information on oceanic crust formation through progressive stages of melt evolution from subduction initiation magmatism (SIM) to arc infancy, and arc maturation-related magmatism. We present here new trace element and radiogenic Sr–Nd–Pb isotopic data from a Late Cretaceous ophiolite in the Outer Zagros Belt of southern Iran and discuss its magmatic development through SIM in a Neotethyan forearc setting. Initial volcanism produced the stratigraphically oldest boninitic (BON) lavas in the Haji-Abad ophiolite whose trace element ratios as well as Pb isotopic signatures indicate a refractory mantle source, which was contaminated by subducting slab-derived fluids. The BON lavas change structurally upwards into island-arc tholeiitic (IAT) lavas and their fractionated felsic derivatives that display isotopic evidence for the contribution of sediment-derived melts in their magmatic source(s). This melt evolution pattern of the Haji-Abad ophiolite marks a major difference from the progressive melt evolution trends reported from the Izu–Bonin–Mariana (IBM) forearc lava sequence, which starts at the bottom with forearc basalts (FAB) phasing up-section into BON and IAT affinity lavas. IAT lavas in the IBM are in turn overlain by calc-alkaline (CAL) lavas. Our geochemical modelling points to a probable existence of a less depleted, pre-subduction mantle whose partial melting probably generated primary magmas before the eruption of boninitic lavas, although the record of such magmas is missing in the Haji-Abad ophiolite. Our comparison of the chemostratigraphy of the Outer Zagros ophiolites with the documented record of subduction initiation generated volcanic rocks in the other Neotethyan ophiolites in the region and the IBM forearc setting strongly suggests that there is no single rule or template for SIM. Supplementary material: Details of analytical procedures and supplementary figures and tables are available at https://doi.org/10.6084/m9.figshare.c.6837597 Thematic collection: This article is part of the Ophiolites, melanges and blueschists collection available at: https://www.lyellcollection.org/topic/collections/ophiolites-melanges-and-blueschists

High-temperature inter-mineral potassium isotope fractionation in ultrapotassic and granitic rocks: Implications for the potassium isotopic compositions of arc magmas

To constrain the mechanism and magnitude of inter-mineral K isotope fractionation during high-temperature magmatic processes, we conducted high-precision K isotope measurements on ultrapotassic and granitic rocks and mineral separates from the Qinling Orogenic Belt, China, as well as first-principles investigations of the corresponding K-bearing minerals. Our results show that the ultrapotassic rocks (−0.33 to −0.14‰) have considerably higher δ41K than the primitive mantle (−0.42 ± 0.08‰), confirming the previous conclusion that slab-derived fluids have contributed heavy K isotopic signatures to their mantle sources. Specifically, the K isotopic fractionation between biotite and microcline (Δ41Kbiotite-microcline = −0.11 to 0.05‰, with biotite and microcline ranging from −0.49 to −0.14‰ and − 0.48 to −0.12‰, respectively) agrees well with the calculated equilibrium isotopic fractionation scale. This evidence, together with the homogeneous mineral chemical compositions, suggests equilibrium K isotopic fractionation between biotite and microcline. Notably, amphiboles (δ41K ranging from −0.41 to 0.48‰) have heterogeneous K isotopic compositions and are generally heavier than coexisting microcline and biotite. In addition, the determined K isotopic fractionation between amphibole/plagioclase and microcline is inconsistent with the calculated results and is likely caused by the kinetic fractionation processes. To reveal the potential controls on the K isotopic compositions of arc magma, we utilized theoretical and measured K isotopic fractionation factors to model K isotopic fractionation behavior during amphibole crystallization. The results show that when Δ41Kamphibole-melt = 0.05‰, the fractional crystallization of amphibole shows negligible effect on the K isotopic composition of the residual melts. Even in the case of Δ41Kamphibole-melt = 0.24‰, fractional crystallization of amphibole exerts an insignificant control on the K isotopic composition of arc magma. Therefore, although amphibole may have heavier K isotopic compositions than the melts, its fractional crystallization may not distinctly decrease the K isotopic composition of residual melts.

The Late Carboniferous Mafic–Ultramafic Complex Induced by Slab Breakoff in Eastern North Tianshan, Central Asian Orogenic Belt

The Late Carboniferous to Early Permian is a critical period of the Chinese Tianshan, witnessing the tectonic transition from subduction to post-collisional extension during the final amalgamation of the Central Asian Orogenic Belt (CAOB). The late Carboniferous Mozbaysay mafic–ultramafic complex in the Qijiaojing–Balikun area, eastern North Tianshan, provides important clues for revealing the nature and timing of this tectonic transition. The Mozbaysay complex comprises mainly hornblende gabbros and lherzolites. LA-ICP-MS U-Pb zircon ages of hornblende gabbro yielded a weighted mean age of 306 ± 1.9 Ma for this complex. These mafic–ultramafic rocks have high contents of MgO (up to 30 wt.%), Cr (up to 2493 ppm), and Ni (up to 1041 ppm), but low contents of SiO2 (40.34–47.70 wt.%). They are enriched in LREE and show characteristics of enriched mid-ocean ridge basalts (E-MORB). The relatively high Th/Yb and Ba/Nb ratios imply the mantle sources could have been metasomatized by slab–mantle interaction with aqueous fluids from dehydration of the subducted slab. Thus, these mafic–ultramafic rocks were most likely produced by partial melting of the asthenospheric and lithospheric mantle with a slight influence of slab-derived fluids. Therefore, we suggest that the formation of these Late Carboniferous mafic–ultramafic rocks was triggered by the decompression-induced influx of asthenospheric heat and melting through a slab window during post-collisional slab breakoff. Combined with geological data, the petrogenetic links of the Late Carboniferous mafic–ultramafic rocks in eastern North Tianshan to slab breakoff suggest that the tectonic transition from convergence to post-collision most likely initiated in situ at ca. 306 Ma and lasted to ca. 300 Ma.

Iron isotope heterogeneity in magmas of subduction zones and its origin

Iron isotope heterogeneity of subduction-related primitive magmas was traditionally thought to be predominantly attributed to the addition of slab-derived components. However, the specific subducted material responsible for inducing Fe isotope heterogeneity in the subarc-backarc mantle remains poorly constrained. Here we report Fe isotopic compositions of a series of backarc basin basalts (BABBs) from the Woodlark Basin (an intra-arc rift basin), Vate Trough (a nascent rear-rift basin), and Lau Basin (a relatively mature basin) in the southwestern Pacific. Even though these BABBs have various geochemical compositions ranging from mid-ocean ridge basalts (MORB)-like to arc-like, which should be affected by various amounts of the slab-derived fluids and/or sediment melts, their δ56Fe values (0.05–0.18‰) are within the range of global MORB (0.05–0.17‰). After correction for crystal fractionation, the primitive δ56Fe values (δ56Fe-prim) displayed no co-variation with indicators of subduction contributions (e.g., Ce/Pb, Nb/U, Ba/La, and Th/Yb) or source redox conditions caused by slab inputs (e.g., V/Yb, V/Ti), suggesting that the δ56Fe-prim values of these BABBs are rarely affected by subduction metasomatism and the change in oxygen fugacity caused by subducted input. Therefore, we proposed that the MORB-like δ56Fe values of the studied BABBs were primarily caused by mantle melting processes. Furthermore, we presented a compilation of available Fe isotopic data for BABBs and island arc basalts (IABs) and used these data to evaluate Fe isotope heterogeneity in subduction settings. In particular, compared with δ56Fe-prim of MORB, the δ56Fe-prim of BABBs from the Central Lau Spreading Centre in the Lau Basin and the IABs have lower values, whereas BABBs from the Rochambeau Ridges in the Lau Basin have higher δ56Fe-prim values. By comprehensively considering the factors affecting the variation in Fe isotopic compositions in magmatic rocks in subduction zones (e.g., crystal fractionation, partial melting of the mantle, redox conditions, and subduction-related metasomatism), we proposed that the varied Fe isotopic compositions in subduction settings could reflect their heterogeneous sources, which may be related to plate subduction. Specifically, during plate subduction, a melt-depleted refractory forearc mantle peridotite (with low δ56Fe values) could be dragged by subducted slab into subarc or backarc depths, resulting in light Fe isotopic enrichment in island arc magmas or BABBs; additionally, serpentinites (both in the slab and mantle wedge) can generate light Fe isotopic enrichment in island arc magmas. High δ56Fe values in the BABB (e.g., from the Rochambeau Ridges in the Lau Basin) may result from the toroidal influx of nearby upwelling plume components or the assimilation of hydrothermally altered oceanic crust.

Coexisting Nb-enriched basalts and arc volcanic rocks in the northern Qiangtang Terrane, China: Implications for the effects of ambient mantle on subduction zone magmatism

The petrogenesis of Nb-enriched basalts (NEBs) is controversial. NEBs are generally considered to form by melting of either a mixed enriched and depleted mantle source or a mantle wedge metasomatized by adakitic melts. Here we present geochronological, petrological, geochemical, and Sr-Nd-Os isotope data for Early Permian volcanic rocks from the Tuotuohe area in the northern Qiangtang Terrane, China. The studied volcanic rocks can be divided into two groups: arc basalts and basaltic andesites in the lower sequence and NEBs in the upper sequence. Zircon grains from one basaltic andesite sample in the lower sequence yielded a concordia age of 302.0 ± 1.1 Ma. The geochemical and isotopic characteristics suggest that the NEBs were derived from a mixed source consisting of enriched mid-oceanic ridge basalts (MORB)- to oceanic-island basalt-like ambient mantle and slab-derived fluids, whereas the arc basalts and basaltic andesites originated from a source made of the MORB-type ambient mantle modified by slab-derived melts. We propose for the first time that the Nb enrichment of the NEBs was most likely inherited from high-Nb ambient mantle and is independent of the addition of slab-derived components. In addition to slab-derived components, the nature of the ambient mantle is also an important factor controlling the diversity of arc magmatism. Our new findings show that the Late Carboniferous−Early Permian magmatism in the northern Qiangtang Terrane was related to northward subduction of the Paleo-Tethys Ocean, suggesting that subduction might have occurred in at least the Early Permian rather than in the previously thought Middle Permian.

Slab dehydration and potassium-lithium recycling in the forearc indicated by potassium and lithium isotope compositions of exhumed metabasites

Slab dehydration at subduction zones has broad implications for arc magmatism, mantle heterogeneity and elemental cycling. We present the petrography, mineral chemistry and whole-rock potassium (K) and lithium (Li) isotopic compositions of subducted metabasites exhumed from serpentinite mud volcanoes from the Mariana forearc that were recovered during International Ocean Discovery Program Expedition 366, to evaluate the capability of K and Li isotopes for tracking slab dehydration and the evolution of the slab-derived fluids at shallow subduction zones. The results show that the δ41K values of HIMU-type alkali metabasites recovered from the shallow-sourced site (i.e., Fantangisña, −0.34 ± 0.06‰, 2SD, n = 3) are slightly higher than those from the deep-sourced site (i.e., Asùt Tesoru, −0.42 ± 0.20‰, 2SD, n = 8). In addition, all the HIMU-type metabasites have homogeneous radiogenic isotope compositions (e.g., Sr, Nd and Hf), and the clinopyroxenes therein show similar trace elements patterns, indicating they have a common mantle source. These observations together with the positive correlation between δ41K and δ11B suggest that K isotope fractionation most likely resulted from metamorphic dehydration. In contrast, no such correlation was observed for Li isotopes (δ7Li ranging from +3.7‰ to +6.6‰). The Li isotope compositions can be modeled by low-temperature fluid/rock interactions with varying water/rock ratios, implying that the Li isotopic system was mainly controlled by fluid metasomatism within the subduction channel. Notably, EM-II-type metabasites have the lightest K isotope compositions (−0.72‰ to −0.67‰) and heaviest Li isotope compositions (+8.1‰ to +8.3‰) in this study. These anomalous signatures were accompanied with the lowest Nb/Th, εNd(t) and εHf(t), indicating that the K and Li isotope data of EM-II-type metabasites were indicative of their mantle source (e.g., EM-II mantle reservoir). Rayleigh fractionation modeling illustrates that the deeply subducted metabasites are characterized by high K and Li concentrations, and low δ41K, but MORB-like δ7Li. Thus, they may transfer large amounts of fluid-mobile elements into the deep mantle that contributes to the K isotopic heterogeneity in the mantle.

Mantle hydration initiated by Ca metasomatism in a subduction zone: An example from the Chandman meta-peridotite, western Mongolia

Fluids in subduction zones contain various aqueous species and cause metasomatic reactions as well as hydration (i.e., serpentinization) at the slab–mantle interface. However, the nature of elemental transfer during fluid infiltration into the dry mantle in the subduction zone is poorly understood. In this study, we describe novel textures related to orthopyroxene (Opx) pseudomorphs and antigorite veins in the Chandman meta-peridotite body in the Alag Khadny accretionary wedge, western Mongolia. This body consists mainly of meta-harzburgite and meta-dunite, and hydration proceeded pervasively within olivine domains associated with the development of antigorite vein networks. Pseudomorphs after orthopyroxene consist of secondary clinopyroxene (S-Cpx) + tremolite (Tr) ± [secondary olivine (S-Ol) + antigorite (Atg)] whereas the primary clinopyroxene (P-Cpx) is almost unaltered. The pseudomorphs after Opx are rarely cut by Atg veins. The mineral occurrences indicate two main stages of mantle alteration: (1) Opx replacement (Ca-metasomatism) and (2) Atg formation induced by Si-bearing fluid infiltration (i.e., the main hydration event), followed by low-temperature (T) alteration that formed lizardite, chrysotile, and brucite. The mineralogical sequence (S-Cpx + Tr ± S-Ol ± Atg) related to Opx replacement can be explained by replacement at 550–650 °C (assuming a pressure of 1.0–2.0 GPa), which is broadly consistent with temperatures experienced by eclogites in the Chandman area during exhumation. Atg in P-Ol domains could have formed at lower temperatures (<500 °C; assuming a pressure of 1.0–2.0 GPa). Mass balance calculations show that Opx replacement was characterized by gains of Ca and a small amount of H2O, whereas Ol hydration was characterized by gains of Si and H2O. A modeled slab-derived fluid, which was calculated to be in equilibrium with basaltic crustal rocks, indicates the time-integrated fluid flux required for Ca metasomatism of the Chandman meta-peridotite body, with an assumption of the alteration thickness of L = 1000 m, was 0.18–1.5 × 106 m3 fluid m−2 rock, whereas the time-integrated fluid flux required for Si metasomatism was 0.3–1.7 × 103 m3 fluid m−2 rock, which is 102–103 lower than that of the Ca metasomatism. The characteristics of the Chandman meta-peridotite body suggest that: (1) the alteration type and trapped elements in the mantle derived from slab-derived fluids vary with temperature and local Si activity; and (2) intense fracture networks facilitated extensive serpentinization related to the Si metasomatism.

Heterogeneous mantle source compositions for boninite from Bonin and Troodos, evidence from iron isotope variations

Boninite has been widely accepted to be the melting product of highly refractory mantle with the addition of slab-derived fluids. Previous studies have found that boninite, like other subduction-related magmas (e.g., mature arc basalts), has generally consistent chondrite-like iron isotope composition (c. 0.03‰), which is on average lighter than those of oceanic basalts. However, it remains unclear whether the lighter iron isotope composition of subduction-related magmas is inherited from the depleted mantle source or caused by the addition of slab-derived oxidized fluids. Moreover, boninite has highly varied element patterns (e.g., Yb, Zr), reflecting varying degrees of mantle depletion and different contributions of slab-derived fluids. Given the limited data for iron isotope composition of boninite so far available, it is thus required to systematically study iron isotope composition of boninite with different compositions to further understand the petrogenesis of boninite and the iron isotope variation of subduction-related magmas. Here we report iron isotope compositions of boninite from two classic suites on Earth, i.e., those from Bonin islands and submarine forearc (low-Ca boninite; CaO/Al2O3 mainly <0.75) and Troodos ophiolite complex, Cyprus (high-Ca boninite: CaO/Al2O3 = c. 0.85). The δ56Fe values of the Bonin boninite vary from −0.046 ± 0.003‰ to 0.078 ± 0.031‰ (2SD, SD = standard deviation of 4 times repeated analysis; with an average of 0.03 ± 0.03‰, n = 14). The δ56Fe value of the Troodos boninite is relatively constant, i.e., from 0.046 ± 0.010‰ to 0.091 ± 0.016‰ with an average of 0.06 ± 0.01‰ (n = 7). Bonin forearc basalts (FABs) were also analyzed for comparison with δ56Fe of ~0.10‰, except for those FABs with higher Th/U ratio and Cu contents (δ56Fe = 0.01–0.05‰). Together with previously reported δ56Fe data for boninite from New Caledonia (low-Ca boninite with CaO/Al2O3 = ~ 0.4, δ56Fe = ~ +0.03‰), the lower δ56Fe values of the Bonin boninite than Bonin FAB and Troodos boninite reflect the lighter iron isotope composition of the more depleted mantle source. Based on relationships between different indicators for various fluids (e.g., Ba/La, Zr/Sm) and modeling, Bonin boninite is likely produced by partial melting of the refractory mantle with an addition of ~5–10% melts from sediments and oceanic crust, while the formation of Troodos boninite can be ascribed to the melting of less refractory mantle, which has been previously metasomatized by ~5% fluids derived from oceanic crust and sediments.

Application of an Improved Olivine-Melt Thermometer/Hygrometer to the Colima Cone Basanites and Minettes of Western Mexico: Implications for the Mantle Source of Unusually High-MgO Melts

Abstract A collection of quaternary, high-MgO (≤13.4 wt%) basanite and minette cinder and lava cones, with an enhanced arc geochemical signature, are located along the northern margin of the N–S Colima rift in western Mexico. The Colima rift overlies the lithospheric suture between the Jalisco block and Guerrero terrane, as well as the tear between the Rivera and Cocos subducting oceanic plates. From the literature, volatile analyses of olivine-hosted melt inclusions in the Colima cone samples document notably high concentrations of dissolved H2O in the melt (≤ 7.0 wt%) as well as degassing-induced phenocryst growth over a range of depths ≤25 km. In this study, it is shown that the high-MgO character of the Colima suite reflects liquid compositions, consistent with evidence for their rapid transit to the surface, without stalling in a crustal magma chamber. The most Mg-rich olivine analyzed in each sample matches the equilibrium composition at the liquidus based on olivine-melt Mn–Mg and Fe2+–Mg exchange coefficients. Application of both a Mg- and Ni-based olivine-melt thermometer, calibrated on the same experimental data set, to the Colima cone suite provides the temperature and dissolved H2O content at the liquidus. Because the Ni thermometer is insensitive to dissolved H2O in the melt, it gives the actual temperature at the onset of olivine phenocryst growth. For the nine Colima samples that range from 13.4–9.2 wt% MgO, resulting temperatures range from 1221°C to 1056°C (± 6–11°C). In contrast, the Mg thermometer is sensitive to dissolved H2O in the melt, and its application (without a correction of H2O) gives the temperature of olivine crystallization under anhydrous conditions. When the Mg- and Ni-based temperatures are paired, the depression of the liquidus (∆T = TMg–TNi) due to dissolved H2O can be obtained. For the high-MgO (&amp;gt;9 wt%) Colima samples, ∆T values range from 188°C to 109°C. Corrections for the effect of pressure (i.e. evidence that phenocryst growth began at ~700 MPa), increase ∆T by ~21°C. An updated model calibration (on experiments from the literature) that relates ∆T with dissolved H2O in the melt shows that the best fit (R2 = 0.95) is linear, wt% H2O = 0.047*∆T, with a standard error of ±0.5 wt%. Although the experimental data set spans a wide range of melt composition (e.g. 47–58 wt% SiO2, 4.4–10.2 wt% MgO, 1.3–4.9 wt% Na2O, 0.1–5.0 wt% K2O, 0.3–5.3 wt% H2O), no dependence on anhydrous melt composition is resolved. Application of this updated model to the Colima suite gives H2O contents of 5.1–8.8 wt% H2O, consistent with those analyzed in olivine-hosted MIs from the literature. When the thermometry and hygrometry results for the Colima cone suite are compared to those for the adjacent calc–alkaline basalts from the Tancítaro Volcanic Field (TVF) in Michoacán, two distinct linear trends in a plot of wt% H2O vs. temperature are found, indicative of different mantle sources. It is proposed that the high-MgO (&amp;gt;11 wt%) Colima cone melts were derived from a phlogopite-bearing harzburgitic mantle at the base of the Jalisco block lithosphere, whereas both TVF and Colima melts with ≤10 wt% MgO were derived from the asthenosphere (i.e. arc mantle wedge). In both mantle sources, slab-derived fluids were an important source of H2O.

Different effects of serpentinite-derived and slab-derived fluids on arc magmatism in the Chinese Eastern Tianshan: evidence from Mg–B isotopic systematics

Fluids are important media for elemental cycling in subduction zones, but the compositions of such fluids remain obscure, especially those derived from dehydration of serpentinite. To better understand the role of subduction-zone fluids in arc magmatism of the Chinese Eastern Tianshan (CET), we present Mg and B isotopic data for the arc diorites (e.g. Kalatage and Nanshankou). Our results show that the Kalatage diorites have higher δ 26 Mg (−0.23 to −0.13‰) and δ 11 B (−0.04 to +1.08‰) values than the corresponding mantle values, whereas the Nanshankou diorites display relatively lower δ 26 Mg (−0.26 to −0.19‰) and δ 11 B (−10.7 to −2.97‰) values. Because high-temperature magmatic and post-magmatic processes have limited effects on Mg and B isotope fractionation of our samples, their distinct Mg–B isotopic signatures inherited from the mantle sources should be caused by subduction-related metasomatism. Because neither altered oceanic crust (AOC) directly melting nor the incorporation of subducted sediments could induce significant Mg–B isotope fractionation in these rocks, we infer that their variable δ 26 Mg and δ 11 B values were probably associated with subduction fluids that were derived from different sources. Specifically, dehydrated fluids from serpentinites possess high δ 26 Mg and δ 11 B values, which resulted in the heavy Mg–B isotopic compositions of the Kalatage diorites, whereas mantle metasomatized by fluids from AOC and subducted sediments were responsible for the Nanshankou rocks with relatively low δ 26 Mg and δ 11 B values. Taking into account the different effects of serpentinite-derived and slab-derived fluids in the formation of arc magmatism, subduction-zone fluids can be probably an important factor in changing the Mg–B isotope composition of mantle sources. Supplementary material : Supplementary data tables and parameters used in models are available at https://doi.org/10.6084/m9.figshare.c.6702686

Geochemistry of the Neoproterozoic ophiolite from the Bensa-Girja area, southern Ethiopia: Petrogenesis and paleotectonics

The Bensa-Girja area is located in northernmost part of the Neoproterozoic Adola belt, southern Ethiopia. It is largely constituted of mafic rocks (amphibolite and metagabbro), ultramafic suite (serpentinite/serpentinized peridotite, talc schist and subordinate pyroxenite), metasedimentary (metapelitic to semi-metapelitic schists, minor marble and metasandstone) and metavolcaniclastic rocks. Geochemical analysis was carried out on the mafic–ultramafic rocks to characterize their petrogenesis, magma source and paleotectonic setting. The results indicate that the mafic rocks have tholeiitic to boninitic composition and represent primitive to evolved magma. The current serpentinite and serpentinized peridotite were formed after harzburgites that are interpreted as a refractory mantle wedge peridotite, and pyroxenites are considered as cumulate rocks crystallized from partial melts of the peridotite. The geochemical data further suggest that the mantle sources of the mafic-ultramafic rocks were modified by both slab-derived fluids and melts in a supra-subduction zone. The rocks thus display geochemical characteristics typical of island arc system with slight MORB-like affinity, inferred to have formed in a sub-arc or spreading forearc setting. It is interpreted that the studied rocks are representing remnants of a fragmented oceanic lithosphere. This together with field geological data supports the arc accretion model postulated for the crustal growth of the Arabian-Nubian Shield.

Geochronological and geochemical study of the Late Carboniferous volcanic drilling cores from the Piedmont of the Hala'alate Mountain in West Junggar: Implications for stratigraphic division and tectonic evolution

The West Junggar tectono-magmatic processes are critical for reconstructing the amalgamation history of the Central Asian Orogenic Belt (CAOB). However, the lack of consensus regarding the geochronology, petrogenesis, and tectonic setting of the complex igneous rocks poses a challenge to comprehending the Late Paleozoic tectonic evolution of West Junggar. This study focuses on systematic geochronological and geochemical research of well-preserved drill cores from the Piedmont of the Hala'alate Mountain, including rhyolite, felsic tuff, basalt, and andesite. The ages of these rocks in Late Carboniferous (310–301 Ma) are constrained by zircon UPb dating, excluding them from the Jiamuhe Formation of the Early Permian. The rhyolites and felsic tuffs exhibit apparent affinity to the A-type granite series, characterized by high total alkali and FeOT/MgO values. Plagioclase controls the main residual facies and depleted isotopic compositions (εHf(t) = +10.2 to +13.4; εNd(t) = +7.6 to +8.1), indicating that the magma source originated from the partial melting of juvenile continental materials in an extensional background. The basalts and andesites display high LILEs/LREEs, LILEs/HFSEs ratios, along with low HFSEs, which correspond to arc-affinity geochemical compositions. The high zircon εHf(t) (+7.5 to +14.3) and whole-rock εNd(t) (+7.8 to +7.9) values indicate a depleted pyroxene-rich mantle source modified by metasomatic slab-derived fluids. This implies a large-scale upwelling of mantle materials driven by the overlying lithospheric crustal thinning. Combining the temporal variations of crustal thickness recorded by the Sr/Y and (La/Yb)N ratios in igneous rocks, the tectonic processes of West Junggar during the Carboniferous to Early Permian can be divided into three stages. The first stage involves sustained subduction processes from 350 Ma to 320 Ma, followed by the closure of the oceanic crust and collision orogenesis from 320 Ma to 310 Ma. Finally, and a post-collision extension setting occurred after 310 Ma.

Potassium isotope evidence for slab-derived fluids in the sub-arc mantle

Subduction recycles water, other volatiles and incompatible elements into Earth's mantle and plays a major role in exchanges between its surface and interior. Potassium (K), a highly fluid-mobile and incompatible alkali element, is strongly affected by these processes. The emerging K stable isotope system has been proposed as a new tracer of subduction events, e.g. to distinguish the inputs to arc magmas from sediments, altered basalts and other near-surface materials. Yet, it remains unknown if and how the K isotopic signatures of these components are affected by extraction of K-bearing media from dehydrating slabs and their transfer to the mantle wedge. Fragments (xenoliths) of sub-arc mantle carried to the surface by arc magmas may retain these slab-derived components. We report high-precision K isotope data on mantle harzburgite xenoliths from Avacha, an active arc volcano in Kamchatka (Russia) to further investigate the nature of slab-derived fluids and their impact on arc mantle. We find a narrow δ41K range for amphibole-rich cumulates from host andesite magmas (−0.53 ± 0.09‰ to −0.34 ± 0.09‰) that overlaps those of the bulk silicate Earth (BSE) and the upper continental crust (UCC). In contrast, the Avacha peridotites derived from the uppermost lithospheric mantle (refractory melting residues metasomatized by subduction zone fluids) show an extremely broad K isotopic range (δ41K ∼ 3‰) and are commonly enriched in heavy K isotopes (up to +1.16 ± 0.08‰) relative to the BSE and UCC; some of the peridotites show the highest δ41K values for any reported mantle rocks. The new data on the Avacha xenoliths, along with available K isotope data on arc lavas and eclogites, suggest that slab-derived fluids may develop heavy K isotopic signatures, distinct from those in the original subducting materials (e.g., marine sediments and altered oceanic crust). These results provide direct evidence that (a) subduction dehydration and fluid percolation in the mantle lithosphere can significantly fractionate K isotopes, and (b) slab-derived fluids infiltrated the sub-arc mantle wedge. Due to the large K concentration difference between the depleted mantle and slab-derived fluids, K isotopes may be a sensitive tracer of fluid activity in sub-arc mantle.