An official website of the United States government
Sodic volcano-plutonic terranes in the Archean can be well preserved, but why oxidized S-rich sodic magmas and porphyry-type Cu-Au deposits are so rare remains poorly understood. Here we addressed this issue by measuring the S concentration and S6+/ΣS ratio of primary apatite grains in >2.7 Ga felsic volcanic rocks from the well-characterized Neoarchean Abitibi Greenstone Belt of the Superior Province, Canada. Whereas apatite grains in most samples contain low-S concentrations (<0.01 wt%, n = 24), a few apatite samples are S-rich (0.14 ± 0.03 wt%, 1σ) and have low-S6+/ΣS ratios (0.56 ± 0.17; 1σ, n = 4). Samples with S-poor apatite have variable whole-rock La/Yb ratios (generally <30) and zircon 10 000*(Eu/Eu*)/Yb ratios of 11 ± 8 (1σ), which may be products of plume-driven or over-thickened crustal melting. In contrast, the samples with S-rich apatite have elevated La/Yb ratios of 49 ± 15 (1σ), zircon 10 000*(Eu/EuN*)/Yb ratios of 26 ± 7 (1σ), and zircon δ18O values of 5.8 ± 0.1 ‰ (1σ), consistent with a deep, hydrous and homogeneous mantle-like source for the melts dominated by amphibole ± garnet fractionation that is reminiscent of subduction-like process. These are the first reported results documenting the predominant accommodation of relatively reduced S in S-rich apatite grains crystallized from terrestrial silicate melts, possibly reflecting slight oxidation associated with the hydration of Neoarchean mantle and crystal fractionation over the magma evolution. The more common S-poor apatite data suggest that suppressed oxidation of the parental sodic magmas led to weak S emission from Earth’s interior to its evolving surface, explaining the rarity of porphyry-type Cu deposits in >2.7 Ga Archean sodic volcano-plutonic terranes.
more »
« less
Variable Modes of Formation for Tonalite–Trondhjemite–Granodiorite–Diorite (TTG)-related Porphyry-type Cu±Au Deposits in the Neoarchean Southern Abitibi Subprovince (Canada): Evidence from Petrochronology and Oxybarometry
Most known porphyry Cu±Au deposits are associated with moderately oxidized and sulfur-rich, calc-alkaline to mildly alkalic arc-related magmas in the Phanerozoic. In contrast, sodium-enriched tonalite–trondhjemite–granodiorite–diorite (TTG) magmas predominant in the Archean are hypothesized to be unoxidized and sulfur-poor, which together preclude porphyry Cu deposit formation. Here, we test this hypothesis by interrogating the causative magmas for the ∼2·7 Ga TTG-related Côté Gold, St-Jude, and Clifford porphyry-type Cu±Au deposit settings in the Neoarchean southern Abitibi subprovince. New and previously published geochronological results constrain the age of emplacement of the causative magmas at ∼2·74 Ga, ∼2·70 Ga, and∼2·69 Ga, respectively. The dioritic and trondhjemitic magmas associated with Côté Gold and St-Jude evolved along a plagioclase-dominated fractionation trend, in contrast to amphibole-dominated fractionation for tonalitic magma at Clifford. Analyses of zircon grains from the Côté Gold, St-Jude, and Clifford igneous rocks yielded εHf(t)±SD values of 4·5±0·3, 4·2±0·6, and 4·3±0·4, and δ18O±SD values of 5·40±0·11 , 3·91±0·13 , and 4·83±0·12 , respectively. These isotopic signatures indicate that, although these magmas are mantle-sourced with minimal crustal contamination, for the St- Jude and Clifford settings the magmas or their sources may have undergone variable alteration by heated seawater or meteoric fluids. Primary barometric minerals (i.e. zircon, amphibole, apatite, and magnetite–ilmenite) that survived variable alteration and metamorphism (up to greenschist facies) were used for estimating fO2 of the causative magmas. Estimation of magmatic fO2 values, reported relative to the fayalite–magnetite–quartz buffer as FMQ, using zircon geochemistry indicates that the fO2 values of the St-Jude, Côté Gold, and Clifford magmas increase from FMQ –0·3±0·6 to FMQ +0·8±0·4 and to FMQ +1·2±0·4, respectively. In contrast, amphibole chemistry yielded systematically higher fO2 values of FMQ +1·6±0·3 and FMQ +2·6±0·1 for Côté Gold and Clifford, respectively, which are consistent with previous studies that indicate that amphibole may overestimate the fO2 of intrusive rocks by up to 1 log unit. Micro X-ray absorption near edge structure (μ-XANES) spectrometric determination of sulfur (i.e. S6+/ S) in primary apatite yielded ≥ FMQ−0·3 and FMQ+1·4–1·8 for St-Jude and Clifford, respectively. The magnetite–ilmenite mineral pairs from the Clifford tonalite yielded FMQ +3·3±1·3 at equilibrium temperatures of 634±21 ◦C, recording the redox state of the late stage of magma crystallization. Electron probe microanalyses revealed that apatite grains from Clifford are enriched in S (up to 0·1 wt%) relative to those of Côté Gold and St-Jude (below the detection limit), which is attributed to either relatively oxidized or sulfur-rich features of the Clifford tonalite. We interpret these results to indicate that the deposits at Côté Gold and Clifford formed from mildly (∼ FMQ +0·8±0·4) to moderately (∼ FMQ +1·5) oxidized magmas where voluminous early sulfide saturation was probably limited, whereas the St-Jude deposit represents a rare case whereby the ingress of externally derived hydrothermal fluids facilitated metal fertility in a relatively reduced magma chamber (∼ FMQ +0). Furthermore, we conclude that variable modes of formation for these deposits and, in addition, the apparent rarity of porphyry-type Cu–Au deposits in the Archean may be attributed to either local restriction of favorable metallogenic conditions, and/or preservation, or an exploration bias.
more »
« less
- Award ID(s):
- 1924142
- PAR ID:
- 10357801
- Date Published:
- Journal Name:
- Journal of petrology
- ISSN:
- 0022-3530
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Most known porphyry Cu deposits formed in the Phanerozoic and are exclusively associated with moderately oxidized, sulfur-rich, hydrous arc-related magmas derived from partial melting of the asthenospheric mantle metasomatized by slab-derived fluids. Yet, whether similar metallogenic processes also operated in the Precambrian remains obscure. Here we address the issue by investigating the origin, fO2, and S contents of calc-alkaline plutonic rocks associated with the Haib porphyry Cu deposit in the Paleoproterozoic Richtersveld Magmatic Arc (southern Namibia), an interpreted mature island-arc setting. We show that the ca. 1886–1881 Ma ore-forming magmas, originated from a mantle-dominated source with minor crustal contributions, were relatively oxidized (1‒2 log units above the fayalitemagnetite- quartz redox buffer) and sulfur-rich. These results indicate that moderately oxidized, sulfur-rich arc magma associated with porphyry Cu mineralization already existed in the late Paleoproterozoic, probably as a result of recycling of sulfate-rich seawater or sediments from the subducted oceanic lithosphere at that time.more » « less
-
Iron oxide-copper-gold (IOCG) deposits are a vital source of copper and critical elements for emerging clean technologies. Andean-type IOCG deposits form in continental arcs undergoing extension, and they have a temporal relationship with magmatism although they do not exhibit a close spatial relation with the causative intrusions. The processes required to form IOCG deposits and their potential connections to iron oxide–apatite (IOA)-type mineralization remain poorly constrained, as well as the characteristics of magmatism linked to both deposit types. Here we combine zircon U–Pb geochronology with zircon trace element geochemistry of intrusive rocks associated with the Candelaria deposit, one of the world’s largest IOCG deposits, to unravel distinctive signatures diagnostic of magmatic fertility. Our results reveal a marked transition in the geochemistry of intrusions in the Candelaria district, characterized by changes in the redox state, water content and temperature of magmas over time. The oldest magmatic stage (~ 128–125 Ma), prior to the formation of the Candelaria deposit, was characterized by zircon Eu/Eu* ratios of 0.20–0.42, and redox conditions of ΔFMQ − 0.4 to + 1.0. The earliest magmatic stage related to the formation of Fe-rich mineralization at Candelaria (118–115 Ma) exhibits low zircon Eu/Eu* ratios (0.09–0.18), low oxygen fugacity values (ΔFMQ ~− 1.8 to + 0.2) and relatively high crystallization temperatures. In contrast, the youngest stage at ~ 111–108 Ma shows higher zircon Eu/Eu* (~ 0.37–0.69), higher oxygen fugacity values (ΔFMQ ~ + 0.4 to + 1.3) and a decrease in crystallization temperatures, conditions that are favorable for the transport and precipitation of sulfur and chalcophile elements. We conclude that Candelaria was formed through two distinct ore-forming stages: the first associated with a reduced, high temperature, water-poor magma developed under a low tectonic stress, followed by a more oxidized, water-rich, and low temperature magmatic event related to a compressional regime. The first event led to Fe-rich and S-poor IOA-type mineralization, while the second event with geochemical signatures similar to those of porphyry copper systems, generated the Cu- and S-rich mineralization. This late stage overprinted preexisting IOA mineralization resulting in the formation of the giant Candelaria IOCG deposit.more » « less
-
Abstract The mineral apatite, Ca10(PO4)6(F,OH,Cl)2, incorporates sulfur (S) during crystallization from S-bearing hydrothermal fluids and silicate melts. Our previous studies of natural and experimental apatite demonstrate that the oxidation state of S in apatite varies systematically as a function of oxygen fugacity (fO2). The S oxidation states –1 and –2 were quantitatively identified in apatite crystallized from reduced, S-bearing hydrothermal fluids and silicate melts by using sulfur K-edge X-ray absorption near-edge structure spectroscopy (S-XANES) where S 6+/ΣS in apatite increases from ~0 at FMQ-1 to ~1 at FMQ+2, where FMQ refers to the fayalite-magnetite-quartz fO2 buffer. In this study, we employ quantum-mechanical calculations to investigate the atomistic structure and energetics of S(-I) and S(-II) incorporated into apatite and elucidate incorporation mechanisms. One S(-I) species (disulfide, S22−) and two S(-II) species (bisulfide, HS−, and sulfide, S2−) are investigated as possible forms of reduced S species in apatite. In configuration models for the simulation, these reduced S species are positioned along the c-axis channel, originally occupied by the column anions F, Cl, and OH in the end-member apatites. In the lowest-energy configurations of S-incorporated apatite, disulfide prefers to be positioned halfway between the mirror planes at z = 1/4 and 3/4. In contrast, the energy-optimized bisulfide is located slightly away from the mirror planes by ~0.04 fractional units in the c direction. The energetic stability of these reduced S species as a function of position along the c-axis can be explained by the geometric and electrostatic constraints of the Ca and O planes that constitute the c-axis channel. The thermodynamics of incorporation of disulfide and bisulfide into apatite is evaluated by using solid-state reaction equations where the apatite host and a solid S-bearing source phase (pyrite and Na2S2(s) for disulfide; troilite and Na2S(s) for sulfide) are the reactants, and the S-incorporated apatite and an anion sink phase are the products. The Gibbs free energy (ΔG) is lower for incorporation with Na-bearing phases than with Fe-bearing phases, which is attributed to the higher energetic stability of the iron sulfide minerals as a source phase for S than the sodium sulfide phases. The thermodynamics of incorporation of reduced S is also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS(aq)(2−n) and HnS(aq)(2−n); n = 0, 1, and 2] as a source phase. The ΔG of S-incorporation increases for fluorapatite and chlorapatite, and decreases for hydroxylapatite, as these species are protonated (i.e., as n changes from 0 to 2). These thermodynamic results demonstrate that the presence of reduced S in apatite is primarily controlled by the chemistry of magmatic and hydrothermal systems where apatite forms (e.g., an abundance of Fe; solution pH). Ultimately, our methodology developed for evaluating the thermodynamics of S incorporation in apatite as a function of temperature, pH, and composition is highly applicable to predicting the trace and volatile element incorporation in minerals in a variety of geological systems. In addition to solid-solid and solid-liquid equilibria treated here at different temperatures and pH, the methodology can be easily extended to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures.more » « less
-
The mineral apatite, Ca10(PO4)6(F,OH,Cl)2, incorporates sulfur (S) during crystallization from S-bearing hydrothermal fluids and silicate melts. Our previous studies of natural and experimental apatite demonstrate that the oxidation state of S in apatite varies systematically as a function of oxygen fugacity (fO2). The S oxidation states –1 and –2 were quantitatively identified in apatite crystallized from reduced, S-bearing hydrothermal fluids and silicate melts by using sulfur K-edge X‑ray absorption near-edge structure spectroscopy (S-XANES) where S6+/ΣS in apatite increases from ~0 at FMQ-1 to ~1 at FMQ+2, where FMQ refers to the fayalite-magnetite-quartz fO2 buffer. In this study, we employ quantum-mechanical calculations to investigate the atomistic structure and energetics of S(-I) and S(-II) incorporated into apatite and elucidate incorporation mechanisms. One S(-I) species (disulfide, S22−) and two S(-II) species (bisulfide, HS−, and sulfide, S2−) are investigated as possible forms of reduced S species in apatite. In configuration models for the simulation, these reduced S species are positioned along the c-axis channel, originally occupied by the column anions F, Cl, and OH in the end-member apatites. In the lowest-energy configurations of S-incorporated apatite, disulfide prefers to be positioned halfway between the mirror planes at z = 1/4 and 3/4. In contrast, the energy-optimized bisulfide is located slightly away from the mirror planes by ~0.04 fractional units in the c direction. The energetic stability of these reduced S species as a function of position along the c-axis can be explained by the geometric and electrostatic constraints of the Ca and O planes that constitute the c-axis channel. The thermodynamics of incorporation of disulfide and bisulfide into apatite are evaluated by using solid-state reaction equations where the apatite host and a solid S-bearing source phase (pyrite and Na2S2(s) for disulfide; troilite and Na2S(s) for sulfide) are the reactants, and the S-incorporated apatite and an anion sink phase are the products. The Gibbs free energy (ΔG) is lower for incorporation with Na-bearing phases than with Fe-bearing phases, which is attributed to the higher energetic stability of the iron sulfide minerals as a source phase for S than the sodium sulfide phases. The thermodynamics of incorporation of reduced S are also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS2(aq)(2–n) and HnS(aq)(2–n); n = 0, 1, and 2] as a source phase. The ΔG of S-incorporation increases for fluorapatite and chlorapatite and decreases for hydroxylapatite as these species are protonated (i.e., as n changes from 0 to 2). These thermodynamic results demonstrate that the presence of reduced S in apatite is primarily controlled by the chemistry of magmatic and hydrothermal systems where apatite forms (e.g., an abundance of Fe; solution pH). Ultimately, our methodology developed for evaluating the thermodynamics of S incorporation in apatite as a function of temperature, pH, and composition is highly applicable to predicting the trace and volatile element incorporation in minerals in a variety of geological systems. In addition to solid-solid and solid-liquid equilibria treated here at different temperatures and pH, the methodology can be easily extended also to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
