Next Article in Journal
Post-Depositional Alteration of Calcium Carbonate Phases in Archaeological Ceramics: Depletion and Redistribution Effects
Next Article in Special Issue
Halogens in Eclogite Facies Minerals from the Western Gneiss Region, Norway
Previous Article in Journal
The Late Jurassic–Palaeogene Carbonate Platforms in the Outer Western Carpathian Tethys—A Regional Overview
Previous Article in Special Issue
Decarbonation Reactions Involving Ankerite and Dolomite under upper Mantle P,T-Parameters: Experimental Modeling
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Minerals
Volume 11
Issue 7
10.3390/min11070748
minerals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Channelized CO2-Rich Fluid Activity along a Subduction Interface in the Paleoproterozoic Wutai Complex, North China Craton

by
Bin Wang
1,
Wei Tian
1,*,
Bin Fu
2 and
Jia-Qi Fang
1,3
1
MOE Key Laboratory of Orogenic Belts and Crustal Evolution, School of Earth and Space Sciences, Peking University, Beijing 100871, China
2
Research School of Earth Sciences, The Australian National University, Canberra, ACT 2601, Australia
3
Seismological Laboratory, California Institute of Technology, Pasadena, CA 91125, USA
*
Author to whom correspondence should be addressed.
Minerals 2021, 11(7), 748; https://doi.org/10.3390/min11070748
Submission received: 17 June 2021 / Revised: 7 July 2021 / Accepted: 7 July 2021 / Published: 9 July 2021
(This article belongs to the Special Issue Behaviour of Volatiles and Fluid-Mobile Elements in Subduction Zones)
Browse Figures
Versions Notes

Abstract

:
Greenschist facies metabasite (chlorite schist) and metasediments (banded iron formation (BIF)) in the Wutai Complex, North China Craton recorded extensive fluid activities during subduction-related metamorphism. The pervasive dolomitization in the chlorite schist and significant dolomite enrichment at the BIF–chlorite schist interface support the existence of highly channelized updip transportation of CO2-rich hydrothermal fluids. Xenotime from the chlorite schist has U concentrations of 39–254 ppm and Th concentrations of 121–2367 ppm, with U/Th ratios of 0.11–0.62, which is typical of xenotime precipitated from circulating hydrothermal fluids. SHRIMP U–Th–Pb dating of xenotime determines a fluid activity age of 1.85 ± 0.07 Ga. The metasomatic dolomite has δ13CV-PDB from −4.17‰ to −3.10‰, which is significantly lower than that of carbonates from greenschists, but similar to the fluid originated from Rayleigh fractionating decarbonation at amphibolite facies metamorphism along the regional geotherm (~15 °C/km) of the Wutai Complex. The δ18OV-SMOW values of the dolomite (12.08–13.85‰) can also correspond to this process, considering the contribution of dehydration. Based on phase equilibrium modelling, we ascertained that the hydrothermal fluid was rich in CO2, alkalis, and silica, with X(CO2) in the range of 0.24–0.28. All of these constraints suggest a channelized CO2-rich fluid activity along the sediment–basite interface in a warm Paleoproterozoic subduction zone, which allowed extensive migration and sequestration of volatiles (especially carbon species) beneath the forearc.

1. Introduction

Subduction zones are dynamic sites of mass and energy exchange between the Earth’s interior and the surface [1,2,3,4]. The presence and movement of fluids and melts at a variety of scales and levels facilitates and governs geochemical cycles and fluid-enhanced deformation within subduction zones [5,6,7,8]. Geochemical studies [6,9,10] and numerical simulations [11] suggest that the fluid flow in subduction zones may be strongly channelized but not pervasive, under the control of the low permeability in subducting slabs, the high dihedral angles between fluid and minerals, and the negative rock volume change during prograde metamorphism [12,13]. Investigation of fluid channelization is therefore critical for the understanding of fluid fluxes and flow paths in the subduction zones [12], and sequentially, for the understanding of the material cycles of substances such as water, carbon, water-soluble salts, and abiogenic hydrocarbons. However, it is difficult to obtain detailed information (for example, the locus, scale, and orientation) about the fluid channels in a subduction zone, since they are located at a great depth or have been dismembered during exhumation and exposure.
Carbon dioxide is one of the principal volatile components in subduction zones. The atmospheric and oceanic CO2 levels regulate the evolution of climate and life on our planet over geological time [14,15,16,17]. Carbon dioxide also buffers the pH range of ore-forming fluids to maintain elevated metal concentrations [18]. The CO2 proportion of fluids within a subduction zone is partly controlled by the metamorphic decarbonation of the subducting slab, which shows significantly different behavior in cold and warm subduction zones [19,20]. In a warm subduction zone, the majority of CO2 expulsion from the slab occurs beneath the forearc [7]; this would have a profound impact on the carbon cycle model and CO2 flux estimates.
This study reports direct observations of fluid channel relics within a shallowly subducted terrane in a warm Paleoproterozoic subduction zone, which provides detailed insights into fluid fluxes and flow paths, as well as its composition, origin, and transportation behavior in a warm subduction setting.

2. Geologic Setting and Samples

2.1. Regional Geology

The Wutai complex—together with Hengshan Complex, Fuping Complex and Hutuo Group—is located in the middle segment of the Trans-North China Orogen, North China Craton (Figure 1a,b) [21]. The Fuping Complex is mainly composed of mafic granulites and late-Neoarchean TTG (tonalite–trondhjemite–granodiorite) gneisses [22,23,24], and is separated from the Wutai Complex by a ductile shear zone trending to the northeast at Longquanguan (Figure 1b). The Hengshan Complex is separated into the northern high-pressure granulite facies terrane (including late-Neoarchean TTG gneisses) [23,25] and the southern amphibolite facies terranes by an EW-directional ductile shear zone at Zhujiafang (Figure 1b). These metamorphic rocks recorded peak P–T conditions at around 12–15 kbar and 750–850 °C [26,27,28,29]. The Wutai Complex comprises metamorphic supracrustal rocks (Wutai Group) and foliated TTG gneisses (Figure 1c), typical of granite–greenstone belts. The TTGs—such as the Ekou, Chechang, Beitai, and Shifo intrusions—are distributed along the regional structural lines (Figure 1c), with emplacement ages of 2.56–2.52 Ga [30,31]. The Wutai Group is subdivided into three subgroups, including the amphibolite facies lower (Shizui) subgroup, the greenschist facies middle (Taihuai) subgroup, and the sub-greenschist facies upper (Gaofan) subgroup (Figure 1c) [32]. The lower subgroup mainly contains garnet amphibolites, garnet–mica schists, and garnet orthoamphibole rocks. They recorded peak P–T conditions at around 9–12 kbar and 620–670 °C [29,33,34,35]. The amphibolite facies to granulite facies metamorphism in the Hengshan, Wutai, and Fuping areas indicates a protracted P–T–t evolution, including a peak stage at ~1.95 Ga, post-peak exhumation, and final slow cooling during 1.92–1.80 Ga [36,37,38]. There was also considerable middle-Paleoproterozoic (2.3–2.0 Ga) mafic and felsic magmatism in these areas [36]. The Hutuo Group consists of sub-greenschist facies metamorphic sediments in a preceding transgressive sequence at 2.15–2.03 Ga, and a succeeding regressive sequence at 1.95–1.80 Ga [39,40].

2.2. Field Occurrence and Sample Description

The study area belongs to the middle subgroup of the Wutai Group, where greenschist facies banded iron formation (BIF), chlorite schists (altered basalts), and pillow basalts occur in a successive distribution from south to north. At the outcrop at Kangjiagou, the BIF overlies the chlorite schist (Figure 2a), and preserves its stratifications (Figure 2b). The chlorite schist has experienced extensive dolomitization and partial muscovitization, with voluminous coarse-grained saddle crystals and thin veins of dolomite (Figure 2c,d). There is significant dolomite enrichment in the BIF–chlorite-schist contact zone (Figure 3a). The pillow basalt is located in the vicinity of Taping; the pillow structure is well preserved, and does not show evidence of mineralization like the chlorite schist (Supplementary Materials, Figure S1).
The dolomite-bearing chlorite schist (samples 16EK05 and 16EK06) mainly consists of chlorite and dolomite (Figure 3b), with subordinate muscovite, albite, quartz, pyrite, and chalcopyrite. Dolomite occurs in rhombic shapes with sweeping rims and obvious twinning striations; it commonly has inclusions of fine-grained quartz and albite. Chlorite occurs as aggregations, and forms the matrix together with quartz and albite. Muscovite generally occurs as individual flakes in the matrix. In addition, a small number of xenotime grains were obtained from sample 16EK06 by heavy liquid and magnetic separation, which can be used to date low-temperature geological processes [41]. The BIF mainly consists of quartz, magnetite, and dolomite in variable granularities; it is heterogeneous and distinct in magnetite-rich (sample 16EK02) and dolomite-rich (sample 17EK11) subtypes (Figure 3c). The dolomite-rich BIF also contains a significant amount of pyrite.

3. Methods

3.1. Instrumental Analyses

The major and trace element compositions of the minerals were analyzed using a JEOL 8230 (manufacturer, city, state, country) electron microprobe analyzer (EMPA) and COMPex Pro102 (manufacturer, city, state, country) and Agilent 7500ce (manufacturer, city, state, country) laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS), respectively, at the Key Laboratory of Orogenic Belts and Crustal Evolution, Peking University, Beijing. On the EMPA, the acceleration voltage and beam current were 15 kV and 10 nA, respectively. The beam was 2 μm in diameter, and counting times were 10–15 s. The SPI 53 minerals standard (U.S.) was used for the quantitative analysis, and PRZ correction was performed at the final calibration stage. Other details are described in [42]. The EMPA-analyzed results and stoichiometrically calculated mineral formulas are presented in Table S1. For the LA-ICP-MS, the laser beam was 36 μm in diameter, and its frequency was 10 Hz. The USGS and Chinese National standards BHVO-2 and GSR-2 were chosen for element concentration calibration. Other details are described in [43].
Carbon–oxygen isotopes of carbonate and oxygen isotopes of silicates that were collected via microdrill from different parts of the rock samples were analyzed using a Thermo Fisher MAT-253 mass spectrometer with a GasBench II system at the Laboratory for Stable Isotope Geochemistry, Institute of Geology and Geophysics, China Academy of Science. Approximately 300 μg of carbonate powders was reacted with 100% phosphoric acid at 70 °C within a 12-mL vial tube filled with He, and the generated CO2 gas was injected into the mass spectrometer for measurement of the isotopic compositions. Internal standard monitoring showed that the standard deviations of δ13C and δ18O were smaller than 0.15‰ and 0.20‰, respectively. Other details are described in [44].
Oxygen isotopes of quartz, magnetite and, incidentally, dolomite were analyzed in situ by the sensitive high-resolution ion microprobe (SHRIMP) II at the Research School of Earth Sciences, Australian National University, Canberra. Data were acquired via a Cs+ primary beam, which was focused into a spot 20–30 μm in diameter on the target surface to generate O− secondary ions from the minerals. Other configurations of the ANU SHRIMP II for oxygen isotopic analysis are described in detail in [45]. Mineral standards of UWQ1 quartz, RSES magnetite, and NBS19 dolomite measured in the same session were used as normalizing reference materials. SHRIMP spots that were not in focus, or overlapped inclusions or adjacent minerals, were rejected. The results are presented in Supplementary Table S2.
Xenotime U–Th–Pb isotope analysis was conducted via SHRIMP at the Research School of Earth Sciences, Australian National University, Canberra. The primary O2− beam was focused to 5–10 μm in diameter, 6 scans of data were acquired in each analysis, and MG-1 and BS-1 xenotimes were used as standards. Other details about the analytic procedures are similar to those described in [46]. Data were reduced using the Isoplot-4.15 Excel macro [47].

3.2. Phase Equilibria Modelling and C–O Isotopes Calculations

Phase equilibria modelling was performed in the Na2O–CaO–K2O–FeO–MgO–Al2O3–SiO2–H2O–CO2 (NCKFMASHC′) system. Pseudosections were calculated using THERMOCALC 3.40 [48] with the ds62 dataset [49]. The activity models of relevant phases were the same as or degenerated from previous coding [50,51,52]. The fluid was assumed to be H2O–CO2, and set to be excess. The normalized bulk rock composition used in the calculations is presented in Supplementary Table S3.
Oxygen isotope fractionation between minerals is temperature dependent. The equilibrium temperature (T, in Kelvin) of oxygen isotopes between quartz and magnetite was calculated using the equation 103 · lnαq-mt = 1.22 × 106/T2 + 8.22 × 103/T − 4.35 [53]. Oxygen isotopes of the water in equilibrium with dolomite and quartz were calculated using the equations 103 · lnαdol-water = 3.06 × 106/T2 − 3.24 [54], and 103 · lnαq-water = 3.05 × 106/T2 − 2.09 [55], respectively.
The oxygen isotopic equilibrium between different lithological samples was assessed using the δ18O–GI diagram [56,57], where GI is the Garlick index defined as (OSi + 0.58OAl)/Ototal in oxygen equivalents [58], applied to minerals or bulk rocks.
The δ18O–δ13C evolving trends of fluids produced during Rayleigh fractionating decarbonation processes were modelled after the equation detailed in [59]. The isotopic fractionation coefficients are cited from [60,61].

4. Results

4.1. Mineral Chemistry

Dolomite is highly stoichiometric, with Ca cations of 0.99–1.01 pfu (per formula unit); it is ferroan, and shows XFe = Fe2+/(Mg + Fe2+) of 0.37–0.44. The PAAS (post-Archaean Australia shale)-normalized REE (rare-earth element) patterns of the dolomite (Table 1; Figure 4a) show depleted light REEs with (La/Sm)N of 0.02–0.06, markedly positive Eu anomalies with (Eu/Eu*)N of 9–27, and slightly negative Y anomalies with (Y/Y*)N of 0.6–0.9. It has particular Sr and Pb concentrations of 142–276 ppm and 8.6–11.8 ppm, respectively. Chlorite shows XFe of 0.52–0.55, and has Si contents of 2.60–2.77 pfu, classified as ripidolite [62,63]. Muscovite shows XNa = Na/(Na + K) of 0.13–0.20 and AlM2 = AlVI/(AlVI + Mg + Fe2+) of 0.91–0.93.

4.2. C–O Isotopes

The porphyroblastic dolomite from the chlorite schist sample 16EK05 has δ13CV-PDB from −4.17‰ to −3.90‰ and δ18OV-SMOW of 13.34–13.59‰. Similarly, dolomite from sample 16EK06 has δ13CV-PDB from −3.83‰ to −3.10‰ and δ18OV-SMOW of 12.08–13.85‰. The chlorite-dominated silicate matrix from sample 16EK05 has δ18OV-SMOW of 5.64‰, while that from sample 16EK06 has δ18OV-SMOW of 4.63‰ (Table 2).
Quartz from the BIF sample 16EK02 has δ18OV-SMOW of 12.01–12.95‰, with a weighted mean of 12.5‰ ± 0.1‰, while that from sample 17EK11 has δ18OV-SMOW of 11.70–12.97‰, with a weighted mean of 12.3‰ ± 0.2‰ (Table 3). Magnetite from sample 16EK02 has δ18OV-SMOW of −6.03–0.18‰, with a weighted mean of −2.8‰ ± 0.7 ‰, while that from sample 17EK11 has δ18OV-SMOW of −3.22–2.51‰, with a weighted mean of −0.6‰ ± 1.9‰ (Table 3). Dolomite from BIF sample 17EK11 has δ13CV-PDB and δ18OV-SMOW values similar to those of the chlorite schist, which are −3.82‰ and 12.33‰, respectively.

4.3. Date and U–Th Contents of Xenotime

Xenotime shows 206Pb/238U dates from 1753 ± 33 Ma to 1916 ± 29 Ma, with a weighted mean of 1852 ± 69 Ma and a lower intercept of 1848 ± 66 Ma (Figure 5a). Xenotime has U contents of 39–254 ppm and Th contents of 121–2367 ppm, with U/Th ratios of 0.11–0.62 (Table 4).

4.4. Pseudosections

The T–X(CO2) (X(CO2) = CO2/(H2O + CO2) in moles) pseudosection (Figure 6a) shows that the carbonate phase varies from calcite and dolomite to siderite as X(CO2) increases. The observed mineral assemblage dominated by dolomite and chlorite in the chlorite schist sample (Figure 3b) indicates a small range of the systemic CO2 concentration, with X(CO2) around 0.24–0.28.
The P–T pseudosection at an X(CO2) of 0.25 (Figure 6b) shows that the carbonate varies from siderite and dolomite to calcite, while the muscovite is transformed to biotite with rising temperatures. The observed mineral assemblage in the chlorite schist (dol + chl + q + ab + mu) defines a narrow-banded P–T field with a positive slope; it constrains temperatures around 260–400 °C at the variable pressures of 1–6 kbar.

4.5. Oxygen Isotopic Equilibrium

Quartz and magnetite in the BIF reached oxygen isotopic equilibrium at temperatures of 246–282 °C; this temperature range is similar to that recorded by chlorite in the chlorite schist (234–250 °C; Table S1) [73,74,75]. Oxygen isotopes (δ18OV-SMOW) of H2O in equilibrium with dolomite in the chlorite schist and quartz in the BIF were estimated to be 3.7–5.4‰ and 3.5–4.4‰, respectively. These similar values suggest a good oxygen isotopic equilibrium between the chlorite schist and BIF samples, which can also be supported by the δ18O–GI diagram (Figure S2). Therefore, H2O in the hydrothermal fluid should have δ18OV-SMOW around 3.5–5.4‰, consistent with a metamorphic fluid [76].

5. Discussion

5.1. Hydrothermal Fluid Activity

5.1.1. Evidence for Fluid Activity

Dolomite occurs in both dolomite-bearing chlorite schist and BIF samples; its saddle texture—especially the turbid inclusion-rich core and the clean sweeping rim (Figure 3b)—indicates a hydrothermal genesis [77,78,79,80,81]. This is different from sedimentary and burial dolomite that commonly occurs in sucrosic, fibrous, or oolitic textures [82,83,84]. In terms of composition, the highly stoichiometric dolomite is only commonplace as metamorphic or metasomatic units, whereas the sedimentary dolomite is usually nonstoichiometric, with excess Ca cations [85]; it requires facilitative fluid flow for a sufficient supply of Mg–Fe2+ cations [85]. Moreover, the REE patterns of the investigated dolomite display marked inconformity with the typical pattern of sedimentary dolomite (Figure 4a) [64]. The characteristic light REE depletion of the investigated dolomite is consistent with the mafic rocks modified by silica-diluent metasomatism [86]. Furthermore, carbon isotopes of the dolomite (δ13CV-PDB from −4.17‰ to −3.10‰) can distinguish it from marine carbonate (typically 0–4‰) [85]. Oxygen isotopes of the dolomite (δ18OV-SMOW = 12.08–13.85‰) can also distinguish it from the igneous carbonatite and sedimentary carbonates, which have δ18OV-SMOW of 6–10‰ and 15–30‰, respectively [44,87]. In a nutshell, the combined evidence from texture, major and trace compositions, and C–O isotopes of the dolomite in this study reflects unambiguous hydrothermal fluid activity.
Xenotime generally forms in fractionated igneous rocks such as pegmatites, and can also form during diagenesis, hydrothermal mineralization, and metamorphism. The U–Th concentrations of xenotime were linked to its formation environments [68]. Xenotime from the dolomite-bearing chlorite schist has significantly lower U concentrations (39–254 ppm) and U/Th ratios (0.11–0.62); it exclusively overlaps the hydrothermal xenotime from iron and gold deposits (Figure 5b) [68,69,70,71,72]. Such xenotime is associated with low-salinity hydrothermal fluids, and tends to have the lowest U content (<100 ppm) and U/Th ratio (<4) among xenotimes of diverse origins [68]. Therefore, the U–Th concentrations of xenotime in the dolomite-bearing chlorite schist can also indicate hydrothermal fluid activity. Moreover, the precipitation of xenotime during the circulation of hydrothermal fluids was probably responsible for the negative Y anomalies in the dolomite.

5.1.2. Fluid Compositions

The chemical composition of the hydrothermal fluids can be derived from comparison between the metamorphic pillow basalt and the hydrothermally altered metabasalt. The metamorphic pillow basalt, which lacks mineral and textural indications for hydrothermal metasomatism, commonly contains fine-grained quartz, albite, chlorite, and epidote (Figure S1). In contrast, the dolomite-bearing chlorite schist, which underwent hydrothermal metasomatism, contains quartz, albite, chlorite aggregates, dolomite porphyroblasts, and muscovite (Figure 2c,d and Figure 3b). The substantial dolomitization and muscovitization in the metasomatic chlorite schist suggest that the hydrothermal fluid is rich in solutes of CO2 and alkalis (Na and K in general).
It was widely proposed that immiscibility occurs between carbonic and saline fluids due to their different densities and wettability [88,89,90,91,92]. Nevertheless, silica dilution can mutually reinforce the miscibility and solubility of CO2 and alkalis in the hydrothermal fluid. For instance, the silica salting-in effect was proposed long ago, which invoked that the solute alkalis and silica can produce a hydrous complex (Na, K)HSiO3·H2O [93]. Silica would also solve strongly into a concentrated salt solution in metamorphic and metasomatic processes [94]. Recent experiments suggest that silicate dissolution boosts the CO2 concentrations in subduction fluids, probably associated with the formation of complexes containing Si–O–C bonds [95]. Therefore, the hydrothermal fluid is likely also rich in solutes of SiO2.
The hydrothermally metasomatic process could be described using the following reaction: 3(Mg, Fe)4Al4Si2O10(OH)8chl + 6Ca2Al3Si3O12(OH)ep + 6SiO2aq + 24CO2aq + 10(Na, K)+aq = 10(Na, K)Al3Si3O10(OH)2mu + 12Ca(Mg, Fe)(CO3)2dol + 10H+aq. Through this reaction, the greenschist facies metabasalt was metasomatized by the hydrothermal fluid to form dolomite and muscovite, exhausting the epidote that formed in the metamorphic stage. On the other hand, the fluid interaction to the metabasalt would result in an enrichment of alkalis and silica, as mentioned above, and in particular, the CO2-rich fluid may also carry leached gold in, as some low-grade gold enrichment occurs in the carbonated chlorite schist [67].

5.1.3. Origin of the Fluid

The xenotime dating acquires a fluid activity age of 1.85 ± 0.07 Ga, which is synchronous with the regional metamorphic ages of 1.95–1.78 Ga [36,37]. This suggests that the activity of the hydrothermal fluid is related to the regional metamorphic or tectonic events.
The C–O isotopic systematics of metasomatic minerals can be used to constrain the source of the hydrothermal fluid [96]. In the investigated samples, dolomite from the hydrothermally metasomatic chlorite schist and BIF had δ18OV-SMOW of 12.08–13.85‰ and δ13CV-PDB of −4.17–−3.10‰ (Table 2). For such C–O isotopic characteristics, a fluid source of devolatilized metasedimentary rocks is unlikely, given that hydrothermal fluid from this source generally forms carbonates with much heavier oxygen isotopes (δ18OV-SMOW around 17–27‰) [6,10,97]. A fluid source of devolatilized metabasite is more likely, as fluid derived from such a source precipitates carbonates with a large δ18OV-SMOW range of about 8–22‰ [5,96,98], which covers the measured dolomite δ18OV-SMOW values.
To confirm the origin of the fluid further, we compare our data with C–O isotopes of high-pressure veins from published literature, which are thought to represent the conduits of fluid that internally derives from the host rocks during metamorphic devolatilization [10,57,98,99,100,101,102]. The C–O isotopes of vein carbonates in metabasites were compiled in Figure 7 in order to trace the fluids sourced from different facies of metabasite. We did not classify the collected isotopic data by their mineral types, because the C and O isotopes among dolomite, calcite, and siderite are almost identical at 400–600 °C, considering that the isotopic fractionation between these phases is very minor (near 0.5‰) [98,103]. Plots of the carbonates’ C–O isotopes in this study are overlapped with the carbonate minerals occurring in amphibolites, while distinct from those in greenschists (Figure 7). We modelled the C–O isotopic Rayleigh fractionation during decarbonation without (Figure 7a) or with (Figure 7b) the contribution of dehydration in the process of metamorphic devolatilization in carbonated seafloor basalt; Our results show that a coupled decarbonation and dehydration process with decreasing X(CO2) of the released fluids at ~650 °C can well match the analyzed C–O isotopes of dolomite in the chlorite schist (Figure 7b). Moreover, the oxygen isotopes of vein quartz also show inconsistency among different facies of metabasite; for instance, vein quartz in blueschist has δ18OV-SMOW around 17–18‰, whereas that in amphibolite has δ18OV-SMOW around 12–13‰ [5]. The δ18OV-SMOW of quartz (11.70–12.97‰) from the metasomatic BIF also overlaps the δ18OV-SMOW range of amphibolite. All of these facts support the notion that the hydrothermal fluid isotopically in equilibrium with the dolomite and quartz analyzed in this study originated from source rocks of amphibolite facies. This inference is reasonable in terms of volatile supply, because amphibolite facies metamorphism generally produces abundant CO2–H2O fluids through prograde chlorite dehydration and dolomite decarbonation [5,104,105]. The univariant reaction results in a continuous CO2 enrichment with decreasing pressure (Figure 8) [106]. The estimated X(CO2) around 0.24–0.28 is obtained at 12−13 kbar and ~660 °C, which is consistent with the prograde amphibolitization P−T conditions of the metabasites in the Wutai Complex (Figure 8); this could also be supported by the vein calcite occurring in the garnet amphibolite from the Wutai Complex [29]. Moreover, the alkali contents in the hydrothermal fluid may be enriched by leaching the TTGs and metapelites, through the following reaction: 3(Na, K)AlSi3O8 + 2H+ = (Na, K)Al3Si3O10(OH)2 + 6SiO2 + 2(Na, K)+ [107,108].
As mentioned above, the fluid derived from devolatilized amphibolite has precipitated dolomite in the chlorite schist at shallower levels after transportation from amphibolite facies’ depth. This implies that the fluid ascended through a channel connecting the deep and shallow levels of a thickened crust. When considering the regional geotherm (~15 °C/km; Figure 8) [29,37], we may infer that the fluid channel occurred either in an orogenic belt or a warm subduction zone. Among these two geotectonic settings, we prefer the latter based on the following reasons: (1) The CO2-, alkali-, and silica-enriched hydrothermal fluid composition shows high uniformity with the high-pressure fluids related to subducted oceanic crust, which are relatively diluent aqueous solutions with dominant solutes of Si, Al, alkali elements, and sometimes CO2 [9,66,99,106,112,119]. The present fluid does not show remarkable enrichment of Al, which should be ascribed to the virtual insolubility of Al at lower pressures [66]—in other words, the Al component is likely to be precipitated in relatively deeper rocks in the subduction channel. On the other hand, the trace element patterns of dolomite in the chlorite schist exhibit variation trends fairly similar to those of the experimentally constrained subduction zone fluids—especially the fluid with the presence of Cl anions (Figure 4b) [65,66]. (2) Plenty of rocks in the Wutai Group have protoliths typical of altered oceanic crustal materials—for instance, the garnet–orthoamphibole rock [35,120] and equivalent chlorite schist, the pillow basalt with microbial alteration textures [121], the volcanogenic massive sulfides (VMS) [122,123], and the epidosite with extremely low δ18O values (our unpublished data). (3) The metamorphic peak P–T conditions of the Wutai and Hengshan complexes lie along the oceanic slab-top P–T path representative of early warm subduction (Figure 8) [117,118]—in particular, that of Hengshan, falling exactly on the ancient (2.7–0.8 Ga) subduction P–T ranges (Figure 8) [115,116]. The Wutai and Hengshan metamorphic complexes share similar geothermal gradients with the typical warm subduction complexes represented by Catalina and Cuba (Figure 8) [112,113,114].

5.2. Geological Implications

Generally, most of the pore fluid in subducting slabs is expelled at shallow levels (<15 km) of the subduction channel, and some of the fluid bonded in mineral crystals is released during devolatilization, enters the overlying mantle wedge, and contributes to arc magmatism. These fluids can escape to the surface through venting at the trench seafloor or arc volcanos [5,124]. Another significant fraction of fluid will be retained in more deeply subducted slabs, and can be stored or redistributed in deeper mantle [5,124]. In this study, the greenschist facies metamorphic terrane, including metabasite (the chlorite schist) and BIF, was metasomatized by hydrothermal fluid that was generated from the devolatilization of amphibolite facies rocks in a warm subduction zone. This suggests that there was an updip flow of fluids along the subduction channel at a 50-km-shallower level (Figure 9), which means that a substantial part of the devolatilized fluids in the subducting slab would return towards the surface and be immobilized in shallow levels. For instance, these fluids can penetrate into shallower rocks and lead to reaction and precipitation of silicate and/or carbonates, which in turn reduces trace element solubility and, hence, precipitates trace-element-rich minerals [112,125]. Dolomite and xenotime in the present samples should form through this mechanism. The incoming rocks with earlier absorbed volatiles—such as the dolomite-bearing chlorite schist—would regenerate the fluid at relevant depths. This ongoing cycle could keep efficient volatile sequestration and enrichment beneath the forearc.
On the other hand, petrologic observations do not indicate large-scale fluid flow in the regional outcrops, which suggests that the fluid flow was in low flux or highly channelized. Numerical simulation shows that the lithological boundaries are mechanically weak, cause enhanced permeability, and facilitate fluid flow [11,126]. In this study, the hydrothermally metasomatic zone lies subparallel to the BIF and chlorite schist boundary (Figure 2a), which represents the interface between altered upper oceanic crust and marine sediments. The fluid-derived dolomite is enriched at the BIF–chlorite-schist contact zone (Figure 3a). These observations indicate that most fluid migrates along the sediment–metabasite (BIF–chlorite-schist) interface within the subducting slab by channelization (Figure 9), rather than by upward flow into the leading edge of the forearc by simple buoyancy drive.
Minerals 11 00748 g009 550
Figure 9. A conceptual schematic illustrating the channelized fluid’s up-dip flow along the sediment–basite interface beneath the leading edge of the forearc in a warm subduction zone. The lithological structure of the subducted slab is modified from [127], and the structure parameters of the hot subduction zone—including slab depth, dehydration loci, and metamorphic facies—follow the numerical simulation results in [128]. The bold arrows represent the fluid flow. The BIF, chlorite schist, garnet amphibolite (GA), and garnet–orthoamphibole rock (GOR) in the Wutai area are labeled at the relevant positions.
Figure 9. A conceptual schematic illustrating the channelized fluid’s up-dip flow along the sediment–basite interface beneath the leading edge of the forearc in a warm subduction zone. The lithological structure of the subducted slab is modified from [127], and the structure parameters of the hot subduction zone—including slab depth, dehydration loci, and metamorphic facies—follow the numerical simulation results in [128]. The bold arrows represent the fluid flow. The BIF, chlorite schist, garnet amphibolite (GA), and garnet–orthoamphibole rock (GOR) in the Wutai area are labeled at the relevant positions.
Minerals 11 00748 g009

5.2.1. Implications for the Carbon Cycle

Subduction zones control the exchange of materials between the Earth’s interior and surface [3,4]; they play a key role in the carbon cycle. The most critical carbon fluxes in subduction zones are carbonate-bearing rocks sunk into the deep mantle and CO2-bearing fluids returned to the atmosphere by metamorphic and volcanic degassing [14,115,129,130]. Nevertheless, the aforementioned updip flow of fluids along the subduction channel has the potential to transport carbon to the near-surface level (<50 km) in forearc position (Figure 9). For example, the CO2-rich fluid in this study transported carbon from the amphibolite facies metamorphic oceanic crust to much shallower greenschist facies terranes, and sequestered carbon at this level by dolomitization (Figure 9). This can explain why there is no significant transfer of carbon from the footwall into the hanging wall at depths of less than 60–70 km in the subduction zone [131].
The updip carbon flux along the subduction channel at depths of <50 km depends on the geotherm of the subduction zone. Due to the formation and wide stabilization of magnesite–siderite phase carbonate along the low geotherm, metamorphic decarbonation in cold subduction zones is negligible [19,20,132], and only limited carbon (below ~2000 ppm) can be transferred through aqueous fluid solution [131]. In contrast, metamorphic decarbonation in warm subduction zones is remarkable [19,20]; during warm subduction, oceanic slabs experience almost complete decarbonation beneath the forearc, with no further carbon release occurring at or beyond sub-arc depths [7]. In the case of the Wutai Complex, voluminous CO2 mainly resulted from dolomite or calcite decomposition in the presence of H2O generated by chlorite or actinolite breakdown. This CO2 would backflow and accumulate, forming a reservoir by carbonation in the shallow levels of the subduction channel (Figure 9). This carbon reservoir, as well as the associated ores, is only temporary storage on a geological timescale, since it is easily relocated to the surface by exhumation. In a further perspective, the observed C–H–O fluids’ activity may have significance for the exploration of abiogenic hydrocarbons in basins related to warm subduction fluid transportation and transformation, which has become a prospective energy solution [95,133].

5.2.2. Implications for Subduction Lubrication

The lubricating effects in the subduction channel play an important role throughout subduction channel processes [134,135]. For instance, mechanical analysis suggested that incoming unconsolidated sediments in the subduction channel led to a much lower effective friction [136,137]. This is because the pore fluids and weak phyllosilicates in sediments can efficiently lower the shear strength [135,138]. This study confirmed an updip flow of fluids along the subduction channel at depths of <50 km (Figure 9), which would permeate into the subducting rocks and reduce their material strength. For example, the hydrothermally metasomatic dolomite-bearing rocks in this study are apparently softer and weaker than dolomite-free lithologies. Subduction channels should therefore be lubricated by not only the downgoing sediments, but also the backflowing fluids. In particular, the backflowing fluids can soften and lubricate the relatively lower layers in the subducting slab sequence (e.g., the metabasalt and the metagabbro) at depths of <50 km.
As reviewed in [139], the subducted slab materials are exhumed from preferential depths of ~25–35 km and ~70–80 km. The 25–35-km-deep preferential exhumation was attributed to the fluctuation of the deformation mechanism between friction and viscous flow [139], which is commonly envisioned as a result of fluctuations in fluid pressure [140,141]. The updip flow of fluids along the subduction channel would contribute to these fluctuations. Moreover, the fluid transfer and precipitation of pore-filling minerals such as carbonates and silica can change the permeability of rocks [139,142], which would also contribute to the switch between weakening and strengthening of subduction channel rocks.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/min11070748/s1: Figure S1: Pillow basalt petrography; Figure S2: The δ18O–GI diagram; Table S1: EMP analyses for minerals in the chlorite schist; Table S2: SHRIMP oxygen isotope analyses for minerals from the BIF and chlorite schist; Table S3: Bulk rock composition of the chlorite schist.

Author Contributions

Conceptualization, B.W. and W.T.; methodology, B.W., W.T., B.F., and J.-Q.F.; validation, B.W. and W.T.; formal analysis, B.W., B.F., and J.-Q.F.; investigation, B.W., W.T., B.F., and J.-Q.F.; resources, W.T. and B.F.; data curation, B.W., W.T., B.F., and J.-Q.F.; writing—original draft preparation, B.W.; writing—review and editing, W.T.; visualization, B.W.; supervision, W.T.; project administration, W.T.; funding acquisition, W.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China, grant number 2019YFA0708503, and the National Natural Science Foundation of China, grant numbers 41973028 and 42030304.

Data Availability Statement

Some of the data presented in this study are available in the Supplementary Materials, Tables S1–S3.

Acknowledgments

Suggestions and comments from Cees-Jan De Hoog and three anonymous reviewers improved the quality of this article. We thank I.H. Campbell, H.W. Li, X.L. Li, and F. Ma for their assistance with experimentations, and Z.P. Wang, M.Y. Gong, and Y.K. Di for their participation in the field work. B. Wang thanks B.F. Han for his course of systematic and frontier lectures about subduction.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Tatsumi, Y.; Eggins, S. Subduction Zone Magmatism; Blackwell Science: Oxford, UK, 1995; 211p. [Google Scholar]
  2. Stern, R.J. Subduction zones. Rev. Geophys. 2002, 40, 1012. [Google Scholar] [CrossRef]
  3. Hacker, B.R.; Peacock, S.M.; Abers, G.A.; Holloway, S.D. Subduction factory 2. Are intermediate-depth earthquakes in subducting slabs linked to metamorphic dehydration reactions? J. Geophys. Res. 2003, 108. [Google Scholar] [CrossRef] [Green Version]
  4. Tatsumi, Y. The subduction factory: How it operates in the evolving Earth. GSA Today 2005, 15, 4–10. [Google Scholar] [CrossRef]
  5. Bebout, G. The impact of subduction-zone metamorphism on mantle-ocean chemical cycling. Chem. Geol. 1995, 126, 191–218. [Google Scholar] [CrossRef]
  6. Miller, J.A.; Cartwright, I.; Buick, I.S.; Barnicoat, A.C. An O-isotope profile through the HP-LT Corsican ophiolite, France and its implications for fluid flowduring subduction. Chem. Geol. 2001, 178, 43–69. [Google Scholar] [CrossRef]
  7. Gorman, P.J.; Kerrick, D.M.; Connolly, J.A.D. Modeling open system metamorphic decarbonation of subducting slabs. Geochem. Geophys. Geosystems 2006, 7, 04007. [Google Scholar] [CrossRef] [Green Version]
  8. Gao, X.-Y.; Wang, L.; Chen, Y.-X.; Zheng, Y.-F.; Chen, R.-X.; Huang, F.; Zhang, Q.-Q.; Ji, M.; Meng, Z.-Y. Geochemical evidence from coesite-bearing jadeite quartzites for large-scale flow of metamorphic fluids in a continental subduction channel. Geochim. Cosmochim. Acta 2019, 265, 354–370. [Google Scholar] [CrossRef]
  9. Bebout, G.E.; Barton, M.D. Metasomatism during subduction: Products and possible paths in the Catalina Schist, California. Chem. Geol. 1993, 108, 61–92. [Google Scholar] [CrossRef]
  10. Cartwright, I.; Barnicoat, A.C. Stable isotope geochemistry of Alpine ophiolites: A window to ocean-floor hydrothermal alteration and constraints on fluid-rock interaction during high-pressure metamorphism. Int. J. Earth Sci. 1999, 88, 219–235. [Google Scholar] [CrossRef]
  11. Breeding, C.M.; Ague, J.J.; Brocker, M.; Bolton, E.W. Blueschist preservation in a retrograded, high-pressure, lowtemperature metamorphic terrane, Tinos, Greece: Implications for fluid flow paths in subduction zones. Geochem. Geophys. Geosyst. 2003, 4, 9002. [Google Scholar] [CrossRef] [Green Version]
  12. Zack, T.; John, T. An evaluation of reactive fluid flow and trace element mobility in subducting slabs. Chem. Geol. 2007, 239, 199–216. [Google Scholar] [CrossRef]
  13. Konrad-Schmolke, M.; O’Brien, P.J.; Zack, T. Fluid Migration above a Subducted Slab—Constraints on Amount, Pathways and Major Element Mobility from Partially Overprinted Eclogite-facies Rocks (Sesia Zone, Western Alps). J. Pet. 2011, 52, 457–486. [Google Scholar] [CrossRef] [Green Version]
  14. Sleep, N.H.; Zahnle, K. Carbon dioxide cycling and implications for climate on ancient Earth. J. Geophys. Res. 2001, 106, 1373–1399. [Google Scholar] [CrossRef]
  15. Lee, C.-T.A.; Shen, B.; Slotnick, B.S.; Liao, K.; Dickens, G.R.; Yokoyama, Y.; Lenardic, A.; Dasgupta, R.; Jellinek, M.; Lackey, J.S.; et al. Continental arc-island arc fluctuations, growth of crustal carbonates, and long-term climate change. Geosphere 2012, 9, 21–36. [Google Scholar] [CrossRef]
  16. Brune, S.; Williams, S.E.; Müller, R.D. Potential links between continental rifting, CO2 degassing and climate change through time. Nat. Geosci. 2017, 10, 941–946. [Google Scholar] [CrossRef]
  17. Aiuppa, A.; Fischer, T.P.; Plank, T.; Bani, P. CO2 flux emissions from the Earth’s most actively degassing volcanoes, 2005–2015. Sci. Rep. 2019, 9, 5442. [Google Scholar] [CrossRef]
  18. Phillips, G.N.; Evans, K. Role of CO2 in the formation of gold deposits. Nature 2004, 429, 860–863. [Google Scholar] [CrossRef]
  19. Kerrick, D.M.; Connolly, J.A.D. Metamorphic devolatilization of subducted marine sediments and the transport of volatiles into the Earth’s mantle. Nature 2001, 411, 293–296. [Google Scholar] [CrossRef]
  20. Kerrick, D.; Connolly, J. Metamorphic devolatilization of subducted oceanic metabasalts: Implications for seismicity, arc magmatism and volatile recycling. Earth Planet. Sci. Lett. 2001, 189, 19–29. [Google Scholar] [CrossRef]
  21. Zhao, G.; Sun, M.; Wilde, S.; Li, S. Late Archean to Paleoproterozoic evolution of the North China Craton: Key issues revisited. Precambrian Res. 2005, 136, 177–202. [Google Scholar] [CrossRef]
  22. Zhao, G.C.; Wilde, S.A.; Cawood, P.A.; Sun, M. SHRIMP U-Pb zircon ages of the Fuping Complex: Implications for accretion and assembly of the North China Craton. Am. J. Sci. 2002, 302, 191–226. [Google Scholar] [CrossRef]
  23. Kröner, A.; Wilde, S.A.; Li, J.H.; Wang, K.Y. Age and evolution of a late Archaean to early Palaeozoic upper to lower crustal section in the Wutaishan/Hengshan/Fuping terrain of northern China. J. Asian Earth Sci. 2005, 24, 577–595. [Google Scholar] [CrossRef]
  24. Tang, L.; Santosh, M.; Tsunogae, T.; Koizumi, T.; Hu, X.; Teng, X.-M. Petrology, phase equilibria modelling and zircon U-Pb geochronology of Paleoproterozoic mafic granulites from the Fuping Complex, North China Craton. J. Metamorph. Geol. 2017, 35, 517–540. [Google Scholar] [CrossRef]
  25. Kröner, A.; Wilde, S.A.; O’Brien, P.J.; Li, J.H.; Passchier, C.W.; Walte, N.P.; Liu, D.Y. Field relationships, geochemistry, zircon ages and evolution of a late Archean to Paleoproterozoic lower crustal section in the Hengshan Terrain of Northern China. Acta Geol. Sin. 2005, 79, 605–629. [Google Scholar]
  26. Zhao, G.; Cawood, P.A.; Wilde, S.; Lu, L. High-Pressure Granulites (Retrograded Eclogites) from the Hengshan Complex, North China Craton: Petrology and Tectonic Implications. J. Pet 2001, 42, 1141–1170. [Google Scholar] [CrossRef] [Green Version]
  27. Zhang, Y.; Wei, C.; Tian, W.; Zhou, X. Reinterpretation of metamorphic age of the Hengshan Complex, North China Craton. Chin. Sci. Bull. 2013, 58, 4300–4307. [Google Scholar] [CrossRef] [Green Version]
  28. Zhang, Y.; Wei, C.; Lu, M.; Zhou, X. P–T–t evolution of the high-pressure mafic granulites from northern Hengshan, North China Craton: Insights from phase equilibria and geochronology. Precambrian Res. 2018, 312, 1–15. [Google Scholar] [CrossRef]
  29. Qian, J.; Wei, C.J. P-T-t evolution of garnet amphibolites in the Wutai-Hengshan area, North China Craton: Insights from phase equilibria and geochronology. J. Metamorph. Geol. 2016, 34, 423–446. [Google Scholar] [CrossRef]
  30. Wilde, S.A.; Cawood, P.A.; Wang, K.Y. The relationship and timing of granitoid evolution with respect to felsic volcanism in the Wutai Complex, North China Craton. In Precambrian Geology-Metamorphic Petrology, Proceedings of the 30th International Geological Conference, Beijing, China, August 1997; Qian, X.L., You, Z.D., Halls, H.C., Eds.; CRC Press: Boca Raton, FL, USA, 1997; Volume 17, pp. 75–88. [Google Scholar]
  31. Wilde, S.; Cawood, P.A.; Wang, K.; Nemchin, A.; Zhao, G. Determining Precambrian crustal evolution in China: A case-study from Wutaishan, Shanxi Province, demonstrating the application of precise SHRIMP U-Pb geochronology. Geol. Soc. Lond. Spec. Publ. 2004, 226, 5–26. [Google Scholar] [CrossRef]
  32. Bai, J. The Early Precambrian Geology of Wutaishan; Tianjin Science and Technology Press: Tianjin, China, 1986; pp. 1–435. (In Chinese) [Google Scholar]
  33. Zhao, G.; Cawood, P.; Lu, L. Petrology and P–T history of the Wutai amphibolites: Implications for tectonic evolution of the Wutai Complex, China. Precambrian Res. 1999, 93, 181–199. [Google Scholar] [CrossRef]
  34. Qian, J.H.; Wei, C.J.; Zhou, X.W.; Chu, H. Metamorphic P-T paths and New Zircon U–Pb age data for garnet–mica schist from the Wutai Group, North China Craton. Precambrian Res. 2013, 233, 282–296. [Google Scholar] [CrossRef]
  35. Qian, J.; Wei, C.; Clarke, G.; Zhou, X. Metamorphic evolution and Zircon ages of Garnet–orthoamphibole rocks in southern Hengshan, North China Craton: Insights into the regional Paleoproterozoic P–T–t history. Precambrian Res. 2015, 256, 223–240. [Google Scholar] [CrossRef]
  36. Wei, C.J.; Qian, J.H.; Zhou, X.W. Paleoproterozoic crustal evolution of the Hengshan–Wutai–Fuping region, North China Craton. Geosci. Front. 2014, 5, 485–497. [Google Scholar] [CrossRef] [Green Version]
  37. Qian, J.; Wei, C.; Yin, C. Paleoproterozoic P–T–t evolution in the Hengshan–Wutai–Fuping area, North China Craton: Evidence from petrological and geochronological data. Precambrian Res. 2017, 303, 91–104. [Google Scholar] [CrossRef]
  38. Wei, C.J. Paleoproterozoic metamorphism and tectonic evolutionin Wutai–Hengshan Region, Trans–North China Orogen. Earth Sci. 2018, 43, 24–43, (In Chinese with English abstract). [Google Scholar]
  39. Wilde, S.A.; Zhao, G.C.; Wang, K.Y.; Sun, M. First precise SHRIMP zircon U-Pb ages for the Hutuo Group in Wutaishan: Further evidence for amalgamation of North China Craton. Chin. Sci. Bull. 2004, 49, 83–90. [Google Scholar] [CrossRef]
  40. Du, L.L.; Yang, C.H.; Wang, W.; Ren, L.D.; Wan, Y.S.; Song, H.X.; Geng, Y.S.; Hou, K.J. The re-examination of the age and stratigraphic subdivision of the Hutuo Group in the Wutai Mountains area, North China Craton: Evidence from geology and zircon U-Pb geochronology. Acta Petrol. Sin. 2011, 27, 1037–1055, (In Chinese with English abstract). [Google Scholar]
  41. Rasmussen, B.; Fletcher, I.R.; Muhling, J. Response of xenotime to prograde metamorphism. Contrib. Miner. Pet. 2011, 162, 1259–1277. [Google Scholar] [CrossRef]
  42. Li, X.L.; Zhang, L.F.; Wei, C.J.; Slabunov, A.I.; Bader, T. Quartz and orthopyroxene exsolution lamellae in clinopyroxene and the metamorphic P–T path of Belomorian eclogites. J. Metamorph. Geol. 2018, 36, 11–12. [Google Scholar] [CrossRef]
  43. Yang, L.; Tian, W.; Hou, G.; Liu, S.; Wang, B.; Li, J.; Chen, M. Volcanic succession, petrology, and geochemistry of the Sujiagou komatiite from the North China Craton. Geol. J. 2020, 55, 3265–3282. [Google Scholar] [CrossRef]
  44. Tao, R.; Zhang, L.; Li, S.; Zhu, J.; Ke, S. Significant contrast in the Mg-C-O isotopes of carbonate between carbonated eclogite and marble from the S.W. Tianshan UHP subduction zone: Evidence for two sources of recycled carbon. Chem. Geol. 2018, 483, 65–77. [Google Scholar] [CrossRef]
  45. Ickert, R.; Hiess, J.; Williams, I.; Holden, P.; Ireland, T.; Lanc, P.; Schram, N.; Foster, J.; Clement, S. Determining high precision, in situ, oxygen isotope ratios with a SHRIMP II: Analyses of MPI-DING silicate-glass reference materials and zircon from contrasting granites. Chem. Geol. 2008, 257, 114–128. [Google Scholar] [CrossRef]
  46. Fletcher, I.R.; McNaughton, N.J.; Aleinikoff, J.A.; Rasmussen, B.; Kamo, S.L. Improved calibration procedures and new standards for U–Pb and Th–Pb dating of Phanerozoic xenotime by ion microprobe. Chem. Geol. 2004, 209, 295–314. [Google Scholar] [CrossRef]
  47. Ludwig, K.R. ISOPLOT 3.00: A Geochronological Toolkit for Microsoft Excel; Berkeley Geochronology Center: Berkeley, CA, USA, 2003; pp. 1–39. [Google Scholar]
  48. Powell, R.; Holland, T. An internally consistent dataset with uncertainties and correlations: Applications to geobarometry, worked examples and a computer program. J. Metamorph. Geol. 1988, 6, 173–204. [Google Scholar] [CrossRef]
  49. Holland, T.J.B.; Powell, R. An improved and extended internally-consistent thermodynamic dataset for phases of petrological interest, involving a newequation of state for solids. J. Metamorph. Geol. 2011, 29, 333–383. [Google Scholar] [CrossRef]
  50. White, R.W.; Powell, R.; Phillips, G.N. A mineral equilibria study of the hydrothermal alteration in mafic greenschist facies rocks at Kalgoorlie, Western Australia. J. Metamorph. Geol. 2003, 21, 455–468. [Google Scholar] [CrossRef]
  51. Holland, T.; Powell, R. Activity? Composition relations for phases in petrological calculations: An asymmetric multicomponent formulation. Contrib. Miner. Pet. 2003, 145, 492–501. [Google Scholar] [CrossRef]
  52. White, R.W.; Powell, R.; Holland, T.J.B.; Johnson, T.; Green, E.C.R. New mineral activity-composition relations for thermodynamic calculations in metapelitic systems. J. Metamorph. Geol. 2014, 32, 261–286. [Google Scholar] [CrossRef]
  53. Zheng, Y.F. Calculation of oxygen isotope fractionation in metal oxides. Geochim. Cosmochim. Acta 1991, 55, 2299–2307. [Google Scholar]
  54. Matthews, A.; Katz, A. Oxygen isotope fractionation during the dolomitization of calcium carbonate. Geochim. Cosmochim. Acta 1977, 41, 1431–1438. [Google Scholar] [CrossRef]
  55. Matthews, A.; Beckinsale, R.D. Oxygen isotope equilibration systematics between quartz and water. Am. Mineral. 1979, 64, 232–240. [Google Scholar]
  56. Ganor, J.; Matthews, A.; Schliestedt, M. Post-metamorphic low δ13 C calcite in the Cycladic complex (Greece) and their implications for modeling fluid infiltration processes using carbon isotope compositions. Eur. J. Miner. 1994, 6, 365–380. [Google Scholar] [CrossRef] [Green Version]
  57. Spandler, C.; Hermann, J. High-pressure veins in eclogite from New Caledonia and their significance for fluid migration in subduction zones. Lithos 2006, 89, 135–153. [Google Scholar] [CrossRef]
  58. Garlick, G.D. Oxygen isotope fractionation in igneous rocks. Earth Planet. Sci. Lett. 1966, 1, 361–368. [Google Scholar] [CrossRef]
  59. Valley, T.; John, W. Stable isotope geochemistry of metamorphic rocks. Rev. Mineral. Geochem. 1986, 16, 445–488. [Google Scholar]
  60. Zheng, Y.-F. Oxygen isotope fractionation in carbonate and sulfate minerals. Geochem. J. 1999, 33, 109–126. [Google Scholar] [CrossRef] [Green Version]
  61. Chacko, T.; Cole, D.R.; Horita, J. Equilibrium Oxygen, Hydrogen and Carbon Isotope Fractionation Factors Applicable to Geologic Systems. Rev. Mineral Geochem. 2001, 43, 1–81. [Google Scholar] [CrossRef]
  62. Hey, M.H. A new review of the chlorites. Miner. Mag. J. Miner. Soc. 1954, 30, 277–292. [Google Scholar] [CrossRef]
  63. Foster, M.D. Interpretation of the Composition and a Classification of the Chlorites. Geological Survey. Professional Paper 414-A; United States Government Printing Office: Washington, DC, USA, 1962; pp. 1–33.
  64. Nutman, A.P.; Friend, C.R.L.; Bennett, V.C.; Wright, D.; Norman, M.D. ≥3700 Ma pre-metamorphic dolomite formed by microbial mediation in the Isua supracrustal belt (W. Greenland): Simple evidence for early life? Precambrian Res. 2010, 183, 725–737. [Google Scholar] [CrossRef]
  65. Keppler, H. Constraints from partitioning experiments on the composition of subduction-zone fluids. Nature 1996, 380, 237–240. [Google Scholar] [CrossRef]
  66. Manning, C. The chemistry of subduction-zone fluids. Earth Planet. Sci. Lett. 2004, 223, 1–16. [Google Scholar] [CrossRef]
  67. Liu, Z.H.; Wang, A.J.; Xu, H. Geneticmodel of ductile-shear type gold deposits in Wutai Mountain area. Miner. Depos. 1997, 16, 350–364, (in Chinese with English abstract). [Google Scholar]
  68. McNaughton, N.J.; Rasmussen, B. Geochemical characterisation of xenotime formation environments using U-Th. Chem. Geol. 2018, 484, 109–119. [Google Scholar] [CrossRef]
  69. Rasmussen, B.; Fletcher, I.R.; Muhling, J.R.; Thorne, W.S.; Broadbent, G.C. Prolonged history of episodic fluid flow in giant hematite ore bodies: Evidence from in situ U–Pb geochronology of hydrothermal xenotime. Earth Planet. Sci. Lett. 2007, 258, 249–259. [Google Scholar] [CrossRef]
  70. Rasmussen, B.; Fletcher, I.R.; Muhling, J.R.; Wilde, S.A. In situ U–Th–Pb geochronology of monazite and xenotime from the Jack Hills belt: Implications for the age of deposition and metamorphism of Hadean zircons. Precambrian Res. 2010, 180, 26–46. [Google Scholar] [CrossRef]
  71. Rasmussen, B.; Mueller, A.G.; Fletcher, I.R. Zirconolite and xenotime U–Pb age constraints on the emplacement of the Golden Mile Dolerite sill and gold mineralization at the Mt Charlotte mine, Eastern Goldfields Province, Yilgarn Craton, Western Australia. Contrib. Miner. Pet. 2009, 157, 559–572. [Google Scholar] [CrossRef]
  72. Vielreicher, N.M.; Groves, D.I.; Snee, L.W.; Fletcher, I.R.; McNaughton, N.J. Broad synchroneity of three gold mineralization styles in the Kalgoorlie Goldfield: SHRIMP, U–Pb, and 40Ar/39Ar geochronological evidence. Econ. Geol. 2010, 105, 187–227. [Google Scholar] [CrossRef]
  73. Rausell-Colom, J.A.; Wiewiora, A.; Matesanz, E. Relation between composition and d001 for chlorite. Am. Mineral. 1991, 76, 1373–1379. [Google Scholar]
  74. Nieto, E. Chemical composition of metapelitic chlorites: X-ray diffraction and optical property approach. Eur. J. Mineral. 1997, 9, 829–841. [Google Scholar] [CrossRef] [Green Version]
  75. Battaglia, S. Applying X-Ray Geothermometer Diffraction to a Chlorite. Clays Clay Miner. 1999, 47, 54–63. [Google Scholar] [CrossRef]
  76. Hoefs, J. Stable Isotope Geochemistry; Springer: Berlin/Heidelberg, Germany, 2009; pp. 1–127. [Google Scholar]
  77. Boni, M.; Parente, G.; Bechstaedt, T.; Vivo, B.D.; Iannace, A. Hydrothermal dolomites in SW Sardinia (Italy): Evidence for a widespread late-Variscan fluid flow event (in Paleofluid flow and diagenesis during basin evolution). Sed. Geol. 2000, 131, 181–200. [Google Scholar] [CrossRef]
  78. Lonnee, J.; Al-Aasm, I.S. Dolomitization and fluid evolution in the Middle Devonian sulphur point formation, Rainbow South Field, Alberta: Petrographic and geochemical evidence. Bull. Can. Petrol. Geol. 2000, 48, 262–283. [Google Scholar] [CrossRef]
  79. Al-Aasm, I.; Lonnee, J.; Clarke, J. Multiple fluid flow event and the formation of saddle dolomite: Case studies from Middle Devonian carbonates of the Western Canada Sedimentary Basin. Mar. Pet. Geol. 2002, 19, 209–217. [Google Scholar] [CrossRef]
  80. Chen, D.; Qing, H.; Yang, C. Multistage hydrothermal dolomites in the Middle Devonian (Givetian) carbonates from the Guilin area, South China. Sedimentology 2004, 51, 1029–1051. [Google Scholar] [CrossRef]
  81. Zhang, J.; Hu, W.; Qian, Y.; Wang, X.; Cao, J.; Zhu, J.; Li, Q.; Xie, X. Formation of saddle dolomites in Upper Cambrian carbonates, western Tarim Basin (northwest China): Implications for fault-related fluid flow. Mar. Pet. Geol. 2009, 26, 1428–1440. [Google Scholar] [CrossRef]
  82. Choquette, P.W.; Hiatt, E.E. Shallow-burial dolomite cement: A major component of many ancient sucrosic dolomites. Sedimentology 2008, 55, 423–460. [Google Scholar] [CrossRef]
  83. Hood, A.V.; Wallace, M.W.; Drysdale, R.N. Neoproterozoic aragonite-dolomite seas? Widespread marine dolomite precipitation in Cryogenian reef complexes. Geology 2011, 39, 871–874. [Google Scholar] [CrossRef]
  84. Hood, A.; Wallace, M. Synsedimentary diagenesis in a Cryogenian reef complex: Ubiquitous marine dolomite precipitation. Sediment. Geol. 2012, 255-256, 56–71. [Google Scholar] [CrossRef]
  85. Warren, J. Dolomite: Occurrence, evolution and economically important associations. Earth Sci. Rev. 2000, 52, 1–81. [Google Scholar] [CrossRef]
  86. Friend, C.R.L.; Nutman, A.; Bennett, V.C.; Norman, M. Seawater-like trace element signatures (REE + Y) of Eoarchaean chemical sedimentary rocks from southern West Greenland, and their corruption during high-grade metamorphism. Contrib. Miner. Pet. 2008, 155, 229–246. [Google Scholar] [CrossRef]
  87. Wang, S.-J.; Teng, F.-Z.; Li, S.-G. Tracing carbonate–silicate interaction during subduction using magnesium and oxygen isotopes. Nat. Commun. 2014, 5, 5328. [Google Scholar] [CrossRef] [Green Version]
  88. Bowers, T.S.; Helgeson, H.C. Calculation of the thermodynamic and geochemical consequences of nonideal mixing in the system H2O-CO2-NaC1 on phase relations in geologic systems: Equation of state for H2O-CO2-NaCl fluids at high pressures and temperatures. Geochim. Cosmochim. Acta 1983, 47, 1247–1275. [Google Scholar] [CrossRef]
  89. Trommsdorff, V.; Skippen, G. Vapour loss (“boiling”) as a mechanism for fluid evolution in metamorphic rocks. Contrib. Miner. Pet. 1986, 94, 317–322. [Google Scholar] [CrossRef]
  90. Holness, M.B.; Graham, C.M. Equilibrium dihedral angles in the system H2O-CO2-NaCl-Calcite, and implications for fluid flow during metamorphism. Contrib. Mineral. Petrol. 1991, 108, 368–383. [Google Scholar] [CrossRef]
  91. Selverstone, J.; Franz, G.; Thomas, S.; Getty, S. Fluid variability in 2 GPa eclogites as an indicator of fluid behavior during subduction. Contrib. Miner. Pet. 1992, 112, 341–357. [Google Scholar] [CrossRef]
  92. Duan, Z.; Møller, N.; Weare, J.H. Equation of state for the NaCl-H2O-CO2 system: Prediction of phase equilibria and volumetric properties. Geochim. Cosmochim. Acta 1995, 59, 2869–2882. [Google Scholar] [CrossRef]
  93. Anderson, G.M.; Burnham, C.W. Feldspar solubility and the transport of aluminum under metamorphic conditions. Am. J. Sci. 1983, 283, 283–297. [Google Scholar]
  94. Newtona, R.C.; Manning, C.E. Quartz solubility in H2O-NaCl and H2O-CO2 solutions at deep crust-upper mantle pressures and temperatures: 2–15 kbar and 500–900 °C. Geochim. Cosmochim. Acta 2000, 64, 2993–3005. [Google Scholar] [CrossRef]
  95. Tumiati, S.; Tiraboschi, C.; Sverjensky, D.A.; Pettke, T.; Recchia, S.; Ulmer, P.; Miozzi, F.; Poli, S. Silicate dissolution boosts the CO2 concentrations in subduction fluids. Nat. Commun. 2017, 8, 616. [Google Scholar] [CrossRef] [Green Version]
  96. Van der Straaten, F.; Halama, R.; John, T.; Schenk, V.; Hauff, F.; Andersen, N. Tracing the effects of high-pressuremetasomatic fluids and seawater alteration in blueschistfacies overprinted eclogites: Implications for subduction channel processes. Chem. Geol. 2012, 292, 69–87. [Google Scholar] [CrossRef] [Green Version]
  97. Cook-Kollars, J.; Bebout, G.; Collins, N.C.; Angiboust, S.; Agard, P. Subduction zone metamorphic pathway for deep carbon cycling: I. Evidence from HP/UHP metasedimentary rocks, Italian Alps. Chem. Geol. 2014, 386, 31–48. [Google Scholar] [CrossRef]
  98. Collins, N.C.; Bebout, G.; Angiboust, S.; Agard, P.; Scambelluri, M.; Crispini, L.; John, T. Subduction zone metamorphic pathway for deep carbon cycling: II. Evidence from HP/UHP metabasaltic rocks and ophicarbonates. Chem. Geol. 2015, 412, 132–150. [Google Scholar] [CrossRef] [Green Version]
  99. Becker, H.; Jochum, K.P.; Carlson, R.W. Constraints from high-pressure veins in eclogites on the composition of hydrous fluids in subduction zones. Chem. Geol. 1999, 160, 291–308. [Google Scholar] [CrossRef]
  100. Franz, L.; Oberhansli, R.; Wagner, T.; Dong, S.; Romer, R.L.; Klemd, R.; Schmid, R.M. Eclogite-facies quartz veins within metabasites of the Dabie Shan (eastern China): Pressure-temperature-time-deformation path, composition of the fluid phase and fluid flow during exhumation of high-pressure rocks. Contrib. Mineral. Pet. 2001, 141, 322–346. [Google Scholar] [CrossRef]
  101. Scambelluri, M.; Philippot, P. Deep fluids in subduction zones. Lithos 2001, 55, 213–227. [Google Scholar] [CrossRef]
  102. Zheng, Y.-F.; Gao, T.-S.; Wu, Y.-B.; Gong, B.; Liu, X.M. Fluid flow during exhumation of deeply subducted continental crust: Zircon U-Pb age and O-isotope studies of a quartz vein within ultrahigh-pressure eclogite. J. Metamorph. Geol. 2007, 25, 267–283. [Google Scholar] [CrossRef]
  103. Sheppard, S.M.F.; Schwarcz, H.P. Fractionation of carbon and oxygen isotopes and magnesium between coexisting metamorphic calcite and dolomite. Contrib. Miner. Pet. 1970, 26, 161–198. [Google Scholar] [CrossRef]
  104. Schmidt, M.W.; Poli, S. Experimentally based water budgets for dehydrating slabs and consequences for arc magma generation. Earth Planet. Sci. Lett. 1998, 163, 361–379. [Google Scholar] [CrossRef]
  105. Ague, J.J. Release of CO2 from carbonate rocks during regional metamorphism of lithologically heterogeneous crust. Geology 2000, 28, 1123–1126. [Google Scholar] [CrossRef]
  106. Molina, J.F.; Poli, S. Carbonate stability and fluid composition in subducted oceanic crust: An experimental study on H2O–CO2-bearing basalts. Earth Planet. Sci. Lett. 2000, 176, 295–310. [Google Scholar] [CrossRef]
  107. Schmidt, C.; Watenphul, A. Ammonium in aqueous fluids to 600 °C, 1.3 GPa: A spectroscopic study on the effects on fluid properties, silica solubility, and K-feldspar to muscovite reactions. Geochim. Cosmochim. Acta 2010, 74, 6852–6866. [Google Scholar] [CrossRef]
  108. Wohlers, A.; Manning, C.; Thompson, A.B. Experimental investigation of the solubility of albite and jadeite in H2O, with paragonite+quartz at 500 and 600 °C, and 1–2.25 GPa. Geochim. Cosmochim. Acta 2011, 75, 2924–2939. [Google Scholar] [CrossRef]
  109. Agrinier, P.; Javoy, M.; Smith, D.C.; Pineau, F. Carbon and oxygen isotopes in eclogites, amphibolites, veins and marbles from the Western Gneiss region, Norway. Chem. Geol. Isot. Geosci. Sect. 1985, 52, 145–162. [Google Scholar] [CrossRef]
  110. Ridley, J. Syn-metamorphic gold deposits in amphibolite and granulite facies rocks. Mitt. Österr. Miner. Ges. 1997, 142, 101–110. [Google Scholar]
  111. Vho, A.; Lanari, P.; Rubatto, D.; Hermann, J. Tracing fluid transfers in subduction zones: An integrated thermodynamic and δ18O fractionation modelling approach. Solid Earth 2020, 11, 307–328. [Google Scholar] [CrossRef] [Green Version]
  112. Spandler, C.; Pirard, C. Element recycling from subducting slabs to arc crust: A review. Lithos 2013, 170-171, 208–223. [Google Scholar] [CrossRef]
  113. Sorensen, S.S.; Barton, M.D. Metasomatism and partial melting in a subduction complex Catalina Schist, southern California. Geology 1987, 15, 115–118. [Google Scholar] [CrossRef]
  114. García-Casco, A.; Lázaro, C.; Rojas-Agramonte, Y.; Kröner, A.; Torres-Roldán, R.L.; Núñez, K.; Neubauer, F.; Millán, G.; Blanco-Quintero, I. Partial Melting and Counterclockwise P T Path of Subducted Oceanic Crust (Sierra del Convento Melange, Cuba). J. Pet. 2007, 49, 129–161. [Google Scholar] [CrossRef] [Green Version]
  115. Dasgupta, R.; Hirschmann, M.M. The deep carbon cycle and melting in Earth’s interior. Earth Planet. Sci. Lett. 2010, 298, 1–13. [Google Scholar] [CrossRef]
  116. Brown, M. Duality of thermal regimes is the distinctive characteristic of plate tectonics since the Neoarchean. Geology 2006, 34, 961. [Google Scholar] [CrossRef]
  117. Palin, R.M.; Santosh, M.; Cao, W.; Li, S.S.; Hernández-Uribe, D.; Parsons, A. Secular metamorphic change and the onset of plate tectonics. Earth Sci. Rev. 2020, 207, 103172. [Google Scholar] [CrossRef]
  118. Martin, H.; Moyen, J.-F. Secular changes in tonalite-trondhjemite-granodiorite composition as markers of the progressive cooling of Earth. Geology 2002, 30, 319–322. [Google Scholar] [CrossRef]
  119. Sorensen, S.; Grossman, J. Enrichment of trace elements in garnet amphibolites from a paleo-subduction zone: Catalina Schist, southern California. Geochim. Cosmochim. Acta 1989, 53, 3155–3177. [Google Scholar] [CrossRef]
  120. Qian, J.H.; Wei, C.J.; Zhou, X.W.; Chu, H. Genesis of the megacryst orthoamphi-bole rock from Hengshan Mts, Shanxi Province: Evidence from geochemistry and Sm–Nd isotopic data. Acta Petrol. Sin. 2012, 28, 2819–2830, (in Chinese with English abstract). [Google Scholar]
  121. McLoughlin, N.; Fliegel, D.J.; Furnes, H.; Staudigel, H.; Simonetti, A.; Zhao, G.C.; Robinson, P.T. Assessing the biogenicity and syngenicity of candidate bioalteration textures in pillow lavas of the ~2.52 Ga Wutai greenstone terrane of China. Chin. Sci. Bull. 2010, 55, 188–199. [Google Scholar] [CrossRef]
  122. Li, J.; Kusky, T.; Niu, X.; Jun, F.; Polat, A. Neoarchean Massive Sulfide of Wutai Mountain, North China: A Black Smoker Chimney and Mound Complex Within 2.50 Ga-Old Oceanic Crust. In Developments in Precambrian Geology; Elsevier: Amsterdam, The Netherlands, 2004; pp. 339–362. [Google Scholar]
  123. Polat, A.; Kusky, T.; Li, J.; Fryer, B.; Kerrich, R.; Patrick, K. Geochemistry of Neoarchean (ca. 2.55–2.50 Ga) volcanic and ophiolitic rocks in the Wutaishan greenstone belt, central orogenic belt, North China craton: Implications for geodynamic setting and continental growth. GSA Bull. 2005, 117, 1387. [Google Scholar] [CrossRef]
  124. Peacock, S.A. Fluid Processes in Subduction Zones. Science 1990, 248, 329–337. [Google Scholar] [CrossRef]
  125. Gao, J.; John, T.; Klemd, R.; Xiong, X. Mobilization of Ti–Nb–Ta during subduction: Evidence from rutile-bearing dehydration segregations and veins hosted in eclogite, Tianshan, NW China. Geochim. Cosmochim. Acta 2007, 71, 4974–4996. [Google Scholar] [CrossRef]
  126. Ague, J.J. Models of permeability contrasts in subduction zone mélange: Implications for gradients in fluid fluxes, Syros and Tinos Islands, Greece. Chem. Geol. 2007, 239, 217–227. [Google Scholar] [CrossRef]
  127. Li, J.-L.; Schwarzenbach, E.M.; John, T.; Ague, J.J.; Huang, F.; Gao, J.; Klemd, R.; Whitehouse, M.J.; Wang, X.-S. Uncovering and quantifying the subduction zone sulfur cycle from the slab perspective. Nat. Commun. 2020, 11, 514. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  128. Iwamori, H. Transportation of H2O and melting in subduction zones. Earth Planet. Sci. Lett. 1998, 160, 65–80. [Google Scholar] [CrossRef]
  129. Marty, B.; Tolstikhin, I.N. CO2 fluxes from mid-ocean ridges, arcs and plumes. Chem. Geol. 1998, 145, 233–248. [Google Scholar] [CrossRef]
  130. Ague, J.J.; Nicolescu, S. Carbon dioxide released from subduction zones by fluid-mediated reactions. Nat. Geosci. 2014, 7, 355–360. [Google Scholar] [CrossRef]
  131. Kelemen, P.B.; Manning, C.E. Reevaluating carbon fluxes in subduction zones, what goes down, mostly comes up. Proc. Natl. Acad. Sci. USA 2015, 112, E3997–E4006. [Google Scholar] [CrossRef] [Green Version]
  132. Zhang, L.; Ellis, D.J.; Williams, S.; Jiang, W. Ultra-high pressure metamorphism in western Tianshan, China: Part II. Evidence from magnesite in eclogite. Am. Miner. 2002, 87, 861–866. [Google Scholar] [CrossRef]
  133. Liu, Q.; Wu, X.; Zhu, D.; Meng, Q.; Xu, H.; Peng, W.; Huang, X.; Liu, J. Generation and resource potential of abiogenic alkane gas under organic–inorganic interactions in petroliferous basins. J. Nat. Gas Geosci. 2021, 6, 79–87. [Google Scholar] [CrossRef]
  134. Shreve, R.L.; Cloos, M. Dynamics of sediment subduction, melange formation, and prism accretion. J. Geophys. Res. 1986, 91, 10229–10245. [Google Scholar] [CrossRef]
  135. Sobolev, S.V.; Brown, M. Surface erosion events controlled the evolution of plate tectonics on Earth. Nature 2019, 570, 52–57. [Google Scholar] [CrossRef]
  136. Lamb, S.; Davis, P. Cenozoic climate change as a possible cause for the rise of the Andes. Nature 2003, 425, 792–797. [Google Scholar] [CrossRef]
  137. Tobin, H.J.; Saffer, D.M. Elevated fluid pressure and extreme mechanical weakness of a plate boundary thrust, Nankai Trough subduction zone. Geology 2009, 37, 679–682. [Google Scholar] [CrossRef]
  138. Behr, W.M.; Becker, T.W. Sediment control on subduction plate speeds. Earth Planet. Sci. Lett. 2018, 502, 166–173. [Google Scholar] [CrossRef]
  139. Agard, P.; Plunder, A.; Angiboust, S.; Bonnet, G.; Ruh, J. The subduction plate interface: Rock record and mechanical coupling (from long to short timescales). Lithos 2018, 320-321, 537–566. [Google Scholar] [CrossRef]
  140. Hayman, N.W.; Lavier, L.L. The geologic record of deep episodic tremor and slip. Geology 2014, 42, 195–198. [Google Scholar] [CrossRef]
  141. Behr, W.M.; Kotowski, A.J.; Ashley, K.T. Dehydration-induced rheological heterogeneity and the deep tremor source inwarmsubduction zones. Geology 2018, 46, 475–478. [Google Scholar] [CrossRef]
  142. Kimura, G.; Ludden, J. Peeling oceanic crust in subduction zones. Geology 1995, 23, 217–220. [Google Scholar] [CrossRef]
Minerals 11 00748 g001 550
Figure 1. (a) Tectonic sketch of the North China Craton. (b) Geological map showing the distribution and metamorphic grades of the Wutai, Hengshan, and Fuping complexes and the Hutuo Group. G: granulite facies; A: amphibolite facies; GS: greenschist facies; SGS: sub-greenschist facies; (c) Geological map showing the lithological units in the Wutai Complex, including the three subgroups of Wutai Group and TTG gneisses. The localities of BIF and chlorite schist (Kangjiagou, large star) and pillow basalt (Taping, small star) are also shown. The tectonic division of the North China Craton is in accordance with [21]. Metamorphic grades of each geological unit and classification of subunits in the Wutai Complex are in accordance with [36].
Figure 1. (a) Tectonic sketch of the North China Craton. (b) Geological map showing the distribution and metamorphic grades of the Wutai, Hengshan, and Fuping complexes and the Hutuo Group. G: granulite facies; A: amphibolite facies; GS: greenschist facies; SGS: sub-greenschist facies; (c) Geological map showing the lithological units in the Wutai Complex, including the three subgroups of Wutai Group and TTG gneisses. The localities of BIF and chlorite schist (Kangjiagou, large star) and pillow basalt (Taping, small star) are also shown. The tectonic division of the North China Craton is in accordance with [21]. Metamorphic grades of each geological unit and classification of subunits in the Wutai Complex are in accordance with [36].
Minerals 11 00748 g001
Minerals 11 00748 g002 550
Figure 2. Photographs of the outcrop at Kangjiagou, Wutai Complex. (a) BIF overlies the chlorite schist. (b) BIF with stratifications. (c) Chlorite schist contains voluminous coarse-grained crystals and thin veins of dolomite. (d) Chlorite schist with both dolomitization and muscovitization (the scattered small white flakes).
Figure 2. Photographs of the outcrop at Kangjiagou, Wutai Complex. (a) BIF overlies the chlorite schist. (b) BIF with stratifications. (c) Chlorite schist contains voluminous coarse-grained crystals and thin veins of dolomite. (d) Chlorite schist with both dolomitization and muscovitization (the scattered small white flakes).
Minerals 11 00748 g002
Minerals 11 00748 g003 550
Figure 3. (a) Photograph showing the dolomite enrichment in the BIF–chlorite-schist contact zone. (b) Photomicrograph of the dolomite porphyroblasts and chlorite-dominated matrix in the chlorite schist. (c) Photomicrographs of dolomite-rich and magnetite-rich subtypes of BIF. dol: dolomite; chl: chlorite; mt: magnetite.
Figure 3. (a) Photograph showing the dolomite enrichment in the BIF–chlorite-schist contact zone. (b) Photomicrograph of the dolomite porphyroblasts and chlorite-dominated matrix in the chlorite schist. (c) Photomicrographs of dolomite-rich and magnetite-rich subtypes of BIF. dol: dolomite; chl: chlorite; mt: magnetite.
Minerals 11 00748 g003
Minerals 11 00748 g004 550
Figure 4. (a) Post-Archaean Australia shale (PAAS)-normalized REE patterns of dolomite from the chlorite schist (squares). The representative sedimentary dolomite is also shown (diamonds, data from [64]). (b) MORB-normalized trace element patterns of dolomite from the chlorite schist (squares). Experimentally determined subduction zone (SZ) fluids are also shown (diamonds: Cl-present; triangles: Cl-free; data from [65]). Elements are plotted in order of decreasing compatibility, in accordance with [66].
Figure 4. (a) Post-Archaean Australia shale (PAAS)-normalized REE patterns of dolomite from the chlorite schist (squares). The representative sedimentary dolomite is also shown (diamonds, data from [64]). (b) MORB-normalized trace element patterns of dolomite from the chlorite schist (squares). Experimentally determined subduction zone (SZ) fluids are also shown (diamonds: Cl-present; triangles: Cl-free; data from [65]). Elements are plotted in order of decreasing compatibility, in accordance with [66].
Minerals 11 00748 g004
Minerals 11 00748 g005 550
Figure 5. (a) Backscattered electron (BSE) images of representative xenotime grains. (b) Energy-dispersive spectrometer (EDS) semiquantitative analyses of xenotime. (c) Tera–Wasserburg concordia diagram for the xenotime dates. (d) U/Th–U diagram for different genetic types of xenotime. The shadowed fields summarize the igneous (I), metamorphic (M), diagenetic (D), and hydrothermal (H) xenotimes (data from [68]), with hydrothermal xenotime from iron (squares; data were collected from [69,70]) and gold (diamonds; data were collected from [71,72]) deposits highlighted. The circles are xenotime from the dolomite-bearing chlorite schist in this study.
Figure 5. (a) Backscattered electron (BSE) images of representative xenotime grains. (b) Energy-dispersive spectrometer (EDS) semiquantitative analyses of xenotime. (c) Tera–Wasserburg concordia diagram for the xenotime dates. (d) U/Th–U diagram for different genetic types of xenotime. The shadowed fields summarize the igneous (I), metamorphic (M), diagenetic (D), and hydrothermal (H) xenotimes (data from [68]), with hydrothermal xenotime from iron (squares; data were collected from [69,70]) and gold (diamonds; data were collected from [71,72]) deposits highlighted. The circles are xenotime from the dolomite-bearing chlorite schist in this study.
Minerals 11 00748 g005
Minerals 11 00748 g006 550
Figure 6. T–X(CO2) (a) and P–T (b) pseudosections calculated for the dolomite-bearing chlorite schist. They are shaded with increasing color depth for the fields of increasing degrees of freedom. The observed mineral assemblage in the sample is highlighted. X(CO2) = CO2/(H2O + CO2) in moles. q: quartz; pl: plagioclase; mu: muscovite; chl: chlorite; dol: dolomite; sid: siderite; cc: calcite; bi: biotite.
Figure 6. T–X(CO2) (a) and P–T (b) pseudosections calculated for the dolomite-bearing chlorite schist. They are shaded with increasing color depth for the fields of increasing degrees of freedom. The observed mineral assemblage in the sample is highlighted. X(CO2) = CO2/(H2O + CO2) in moles. q: quartz; pl: plagioclase; mu: muscovite; chl: chlorite; dol: dolomite; sid: siderite; cc: calcite; bi: biotite.
Minerals 11 00748 g006
Minerals 11 00748 g007 550
Figure 7. The δ18O–δ13C variation diagrams for vein carbonates in metabasites of greenschist facies and amphibolite facies, which represent the δ18O–δ13C values of metamorphic fluids in each metamorphic facies. The δ18O–δ13C data of dolomite from the chlorite schist in this study and in previous literature [67] are plotted for comparison. The modelled δ18O–δ13C evolving trends of fluids produced in Rayleigh fractionating decarbonation processes are also shown, without (a) or with (b) the contribution of dehydration. Calibrations along the curves represent the fraction of carbonate that remained in solid phases. The δ18O–δ13C data were collected from [5,98,109,110]. The δ18O of H2O released from metabasalt adopted here is ~9.5‰ [111]. CSB: carbonated seafloor basalt.
Figure 7. The δ18O–δ13C variation diagrams for vein carbonates in metabasites of greenschist facies and amphibolite facies, which represent the δ18O–δ13C values of metamorphic fluids in each metamorphic facies. The δ18O–δ13C data of dolomite from the chlorite schist in this study and in previous literature [67] are plotted for comparison. The modelled δ18O–δ13C evolving trends of fluids produced in Rayleigh fractionating decarbonation processes are also shown, without (a) or with (b) the contribution of dehydration. Calibrations along the curves represent the fraction of carbonate that remained in solid phases. The δ18O–δ13C data were collected from [5,98,109,110]. The δ18O of H2O released from metabasalt adopted here is ~9.5‰ [111]. CSB: carbonated seafloor basalt.
Minerals 11 00748 g007
Minerals 11 00748 g008 550
Figure 8. Phase diagram showing the P–T paths (solid curves with arrows) and peak P–T condition for the metabasites in the Wutai Complex (Z99, [33]; Q15, [35]; QW16, [29]). This figure also shows: (1) experimentally determined phase relationships for H2O-saturated MORB, including the rock’s solidus and the stability ranges of garnet, chlorite, zoisite, and lawsonite (solid curves) [104]; (2) calculated devolatilization reaction (dashed curve with diamonds), carbonates (carb) ± amphibolite (amp) = diopside (di) + fluid [106], labeled with X(CO2) in the resultant fluid along the univariate line (C 0.5, X(CO2) = 0.5); (3) the peak metamorphic P–T conditions (ellipses) of the typical warm subduction complexes Cuba and Catalina [112,113,114], and that of the Hengshan Complex [26,27,28,29]; and (4) the ancient (2.7–0.8 Ga) subduction P–T range (polygonal field labeled with AS) [115,116] and the oceanic slab-top P–T paths representative of Archean examples (dashed curves with arrows) [117,118].
Figure 8. Phase diagram showing the P–T paths (solid curves with arrows) and peak P–T condition for the metabasites in the Wutai Complex (Z99, [33]; Q15, [35]; QW16, [29]). This figure also shows: (1) experimentally determined phase relationships for H2O-saturated MORB, including the rock’s solidus and the stability ranges of garnet, chlorite, zoisite, and lawsonite (solid curves) [104]; (2) calculated devolatilization reaction (dashed curve with diamonds), carbonates (carb) ± amphibolite (amp) = diopside (di) + fluid [106], labeled with X(CO2) in the resultant fluid along the univariate line (C 0.5, X(CO2) = 0.5); (3) the peak metamorphic P–T conditions (ellipses) of the typical warm subduction complexes Cuba and Catalina [112,113,114], and that of the Hengshan Complex [26,27,28,29]; and (4) the ancient (2.7–0.8 Ga) subduction P–T range (polygonal field labeled with AS) [115,116] and the oceanic slab-top P–T paths representative of Archean examples (dashed curves with arrows) [117,118].
Minerals 11 00748 g008
Table
Table 1. Trace elements concentrations (ppm) of dolomite from the chlorite schist.
Table 1. Trace elements concentrations (ppm) of dolomite from the chlorite schist.
Elementdol-1dol-2dol-3dol-4dol-5Elementdol-1dol-2dol-3dol-4dol-5
Sc27.7126.521.4121.8720.29Gd0.1760.3150.14950.1970.64
Rb0.00910.00950.01030.00880.0102Tb0.02190.05380.01950.03190.1059
Sr143.11141.74150.29165.68275.99Dy0.13730.2680.11820.16550.495
Y0.740.9390.5640.5131.75Ho0.03260.05080.02380.03040.0791
Zr29.230.239223.280.320.004Er0.11030.12760.08570.06350.1952
Nb0.00470.2970.00330.00270.0031Tm0.02920.020.01930.01290.0232
Ba0.0690.1140.07410.10660.209Yb0.2280.14290.18760.10090.262
La0.02310.02870.02040.0610.085Lu0.05380.03260.04060.02440.0435
Ce0.08720.13050.0750.14380.321Hf0.9270.01480.6410.00820.0064
Pr0.0220.03630.01370.03170.0707Ta0.00240.01510.00250.00190.0019
Nd0.17410.3080.11510.2410.626Pb11.77711.21710.6399.68.62
Sm0.09970.2340.08470.14870.452Th0.01490.00560.0150.00330.0041
Eu0.5710.5220.5370.9821.027U0.02690.00960.01550.00250.0028
Table
Table 2. C–O isotopes (‰) of porphyroblastic dolomite and silicate matrix in the chlorite schist.
Table 2. C–O isotopes (‰) of porphyroblastic dolomite and silicate matrix in the chlorite schist.
MaterialsSamplesδ13CV-PDBδ18OV-PDBδ18OV-SMOW
Dolomite16EK05−3.90−16.8013.59
−3.99−16.9713.42
−4.17−17.0413.34
16EK06−3.10−18.2712.08
−3.83−16.5513.85
17EK11 a−3.82−18.0212.33
Kang-2 b−4.28-10.94
−4.39-11.62
Matrix16EK05--5.64
16EK06--4.63
a Dolomite-rich BIF. b Cited from [67].
Table
Table 3. Oxygen isotopes (δ18OV-SMOW, ‰) of dolomite (dol), quartz (qz), and magnetite (mt) from the chlorite schist and BIF (analyzed by SHRIMP).
Table 3. Oxygen isotopes (δ18OV-SMOW, ‰) of dolomite (dol), quartz (qz), and magnetite (mt) from the chlorite schist and BIF (analyzed by SHRIMP).
MineralsSamplesRangesWeighted MeanCount
Min.Max.Mean1σ
dol16EK05 a11.5012.0811.840.0820
16EK06 a11.2912.0111.670.0827
qz16EK0212.0112.9512.50.127
17EK1111.7012.9712.30.219
mt16EK02−6.030.18−2.80.727
17EK11a−3.222.51−0.61.98
a not involved in discussion.
Table
Table 4. Xenotime U–Pb isotopic data and corrected ages.
Table 4. Xenotime U–Pb isotopic data and corrected ages.
No.ThUU/ThTotal Isotopic Ratios204Pb Corrected Ages (Ma)
(ppm)(ppm)238U/206Pbsd.%207Pb/206Pbsd.%206Pb/238U1σ sd.207Pb/206Pb1σ sd.
xen1.1121390.323.0673.650.11774.241782611567256
xen1.2212510.242.8823.710.13604.051895631999158
xen2.13802350.622.8771.760.11661.39191629184837
xen3.123672540.112.9401.740.11501.35188629186927
xen3.28041320.163.1562.090.11612.04175333169688
xen3.311381550.143.0101.960.10971.89183332164372
Note: sd—standard deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, B.; Tian, W.; Fu, B.; Fang, J.-Q. Channelized CO2-Rich Fluid Activity along a Subduction Interface in the Paleoproterozoic Wutai Complex, North China Craton. Minerals 2021, 11, 748. https://doi.org/10.3390/min11070748

AMA Style

Wang B, Tian W, Fu B, Fang J-Q. Channelized CO2-Rich Fluid Activity along a Subduction Interface in the Paleoproterozoic Wutai Complex, North China Craton. Minerals. 2021; 11(7):748. https://doi.org/10.3390/min11070748

Chicago/Turabian Style

Wang, Bin, Wei Tian, Bin Fu, and Jia-Qi Fang. 2021. "Channelized CO2-Rich Fluid Activity along a Subduction Interface in the Paleoproterozoic Wutai Complex, North China Craton" Minerals 11, no. 7: 748. https://doi.org/10.3390/min11070748

APA Style

Wang, B., Tian, W., Fu, B., & Fang, J. -Q. (2021). Channelized CO2-Rich Fluid Activity along a Subduction Interface in the Paleoproterozoic Wutai Complex, North China Craton. Minerals, 11(7), 748. https://doi.org/10.3390/min11070748

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop