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 indicate that the fO2 values of the St-Jude, Côté Gold, and Clifford magmas increase from ΔFMQ -0.3 ± 0.6, ΔFMQ +0.8 ± 0.4, 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 amphibole may overestimate the fO2 of intrusive rocks by up to one 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 the 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 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.

Seismic imaging across fault systems in the Abitibi greenstone belt – an analysis of pre- and post-stack migration approaches in the Chibougamau area, Quebec, Canada

Two high-resolution seismic reflection profiles acquired north and south of Chibougamau, located in the northeast of the Abitibi subprovince of Canada, help understand historic volcanically hosted massive sulfide (VMS) deposits and hydrothermal Cu–Au mineralization found there. Major faults crossed by the profiles include the Barlow fault in the north and the Doda fault and the Guercheville fault in the south, all targets of this study that seeks to determine spatial relationships with a known metal endowment in the area. Common-offset DMO corrections and common-offset pre-stack time migrations (PSTMs) were considered. Irregularities of the trace midpoint distribution resulting from the crooked geometry of both profiles and their relative contribution to the DMO and PSTM methods and seismic illumination were assessed in the context of the complex subsurface architecture of the area. To scrutinize this contribution, seismic images were generated for offset ranges of 0–9 km using increments of 3 km. Migration of out-of-plane reflections used cross-dip element analysis to accurately estimate the fault dip. The seismic imaging shows the thickening of the upper-crustal rocks near the fault zones along both profiles. In the northern seismic reflection section, the key geological structures identified include the Barlow fault and two diffraction sets imaged within the fault zone that represent potential targets for future exploration. The south seismic reflection section shows rather a complicated geometry of two fault systems. The Guercheville fault observed as a subhorizontal reflector connects to a steeply dipping reflector. The Doda fault dips subvertical in the shallow crust but as a steeply dipping reflection set at depth. Nearby gold showings suggest that these faults may help channel and concentrate mineralizing fluids.

Plumbing and architecture of a crustal mineralizing system in the Archean Superior Province, Canada

<p>The role of melts and magmatic/metamorphic fluids in mineralization processes is well established. However, the role of crustal architecture in defining source and sink zones in the middle to lower crust remains enigmatic. Integration of three dimensional magnetotelluric (MT) modelling and seismic reflection data across the Archean Abitibi greenstone belt of the Superior Province, Canada, reveals a ‘whole-of-crust’ mineralizing system and highlights the controls by crustal architecture on metallogenetic processes. Electrically conductive conduits in an otherwise resistive upper crust are coincident with truncations and offsets of seismic reflections that are mostly interpreted as major brittle-ductile fault zones. The spatial association between these features and low resistivity zones imaged in the 3D models suggest that these zones acted as pathways through which fluids and melts ascended toward the surface. At mid-crustal levels, these ‘conduit’ zones connect to ~50 km long, north-south striking conductors, and are inferred to represent graphite and/or sulphide deposited from cooling fluids. At upper mantle to lower crustal depths, east-west trending conductive zones dominate and display shallow dips. The upper mantle features are broadly coincident with the surface traces of the major deformation zones with which a large proportion of the gold endowment is associated. We suggest that these deep conductors represent interconnected graphitic zones perhaps augmented by sulphides that are relicts from metamorphic fluid and melt emplacement associated primarily with the later stages of regional deformation.  Thus, from the combined MT and seismic data, we develop a crustal-scale architectural model that is consistent with existing geological and deformational models, providing constraints on the sources for and signatures of fluid and magma emplacement that resulted in widespread metallogenesis in the Abitibi Subprovince.</p>

Petrogenesis and economic potential of the Obatogamau Formation, Chibougamau area, Abitibi greenstone belt

Unravelling the petrogenesis and stratigraphy of Archean mafic lava flows is essential to our comprehension of the early geodynamic evolution and economic potential of greenstone belts. This study focuses on one of the oldest and thickest sequences of lava flows observed in the Neoarchean Abitibi Subprovince (greenstone belt), i.e., the Obatogamau Formation. The undated formation extends for more than 100 km in an E-W direction and consists mostly of aphyric and feldspar megacryst-bearing basaltic-andesite lava flows. These lava flows are tholeiitic, mostly Fe-rich and have nearly homogeneous chemistry. Petrogenetic modelling carried out using MELTS software points to limited magmatic differentiation as most samples of mafic lava flows did not reach Fe-Ti-oxide saturation. Zircon U-Pb dating establishes a crystallization age of 2726.2 ± 1.6 Ma for a felsic unit located at an intermediate stratigraphic position in the sequence of lava flows. Constraints from stratigraphically overlying volcanic units suggest that the Obatogamau Formation was likely emplaced rapidly, possibly within a few million years and as a consequence of frequent replenishment of shallow magma accumulations. High eruption rates are consistent with short episodes of volcanic quiescence deduced from field observations, indicating non-optimal conditions for volcanogenic massive sulfide systems. The pressure and temperature of peak metamorphism deduced from amphibole chemistry, however, points to favorable conditions for the release of metamorphic fluids. The study area may thus be prospective for orogenic gold mineralization, provided that fluids had access to a source of gold and that structural conduits allowed for the channeling of hydrothermal fluids.

Structural analysis and metamorphism of the Barlow Fault Zone, Chibougamau area, Neoarchean Abitibi Subprovince: implications for gold mineralization

Faults and deformation zones are important features of Archean terranes because of their significance for structural evolution and the formation of large gold districts. In the Chibougamau area, northeastern portion of the Abitibi Subprovince, an Apogee Metal Earth seismic reflection survey identified an association between an exceptional shallow-dipping subsurface reflector zone and the Barlow Fault Zone visible at the surface. This study aimed to reconstruct the kinematic evolution of the Barlow Fault Zone, determine its position within the structural setting of this section of the Abitibi Subprovince, and evaluate its importance for gold potential in the northern part of the Chibougamau area. Structural reconstruction and field observations are compatible with a reverse south-over-north movement related to a ductile north–south shortening event, which culminated with amphibolite metamorphism. Geothermobarometers indicate peak metamorphism conditions of 550 ± 50 °C and 6 ± 1.2 kbar. Results from this study suggest that amphibolite facies metamorphism covers a much wider area within the Chibougamau region than previously documented. The Barlow Fault Zone shares similar geometric characteristics and evolutionary history with other gold-bearing structures in the Abitibi Subprovince, but ultimately it was unable to provide optimal conditions for channelling fluid and precipitating gold. The Barlow Fault Zone is interpreted as a back-thrust fault that belongs to a more extensive south-dipping fault system encasing juxtaposed tectonic slivers. This system, with amphibolite facies metamorphism, is a defining feature of the northern portion of the Chibougamau area and developed during the accretion between the Opatica and Abitibi subprovinces.

Oxygen Fugacity and Volatile Content of Syntectonic Magmatism in the Neoarchean Abitibi Greenstone Belt, Superior Province, Canada

Neoarchean syntectonic intrusions from the Chibougamau area, northeastern Abitibi Subprovince (greenstone belt), may be genetically related to intrusion related gold mineralization. These magmatic-hydrothermal systems share common features with orogenic gold deposits, such as spatial and temporal association with syntectonic magmatism. Genetic association with magmatism, however, remains controversial for many greenstone belt hosted Au deposits. To precisely identify the link between syntectonic magmas and gold mineralization in the Abitibi Subprovince, major and trace-element compositions of whole rock, zircon, apatite, and amphibole grains were measured for five intrusions in the Chibougamau area; the Anville, Saussure, Chevrillon, Opémisca, and Lac Line Plutons. The selected intrusions are representative of the chemical diversity of synvolcanic (TTG suite) and syntectonic (e.g., sanukitoid, alkaline intrusion) magmatism. Chemical data enable calculation of oxygen fugacity and volatile content, and these parameters were interpreted using data collected by electron microprobe and laser ablation-inductively coupled plasma-mass spectrometry. The zircon and apatite data and associated oxygen fugacity values in magma indicate that the youngest magmas are the most oxidized. Moreover, similar oxygen fugacity and high volatile content for both the Saussure Pluton and the mineralized Lac Line intrusion may indicate a possible prospective mineralized system associated with the syntectonic Saussure intrusion.

Seismic imaging across fault systems in the Abitibi greenstone belt – An analysis of pre- and post-stack migration approaches in the Chibougamau area, Quebec, Canada

Two high-resolution seismic reflection profiles acquired north and south of Chibougamau, located in the northeast of the Abitibi subprovince of Canada help understand historic volcanic-hosted massive sulfide (VMS) deposits and hydrothermal Cu-Au mineralization found there. Major faults crossed by the profiles include the Barlow fault in the north and the Doda fault and the Guercheville fault in the south, all targets of this study that seeks to determine spatial relationships with known metal endowment in the area. Common-offset DMO corrections and common-offset pre-stack time migrations (PSTM) were considered. Irregularities of the trace midpoint distribution resulting from the crooked geometry of both profiles and their relative contribution to DMO and PSTM method and seismic illumination were assessed in the context of the complex subsurface architecture of the area. To scrutinize this contribution, seismic images were generated for offset ranges of 0–9 km using increments of 3 km. Migration of out-of-plane reflections used cross-dip element analysis to accurately estimate the fault dip. The seismic imaging shows the thickening of the upper crustal rocks near the fault zones along both profiles. In the north seismic reflection section, key geological structure identified include the Barlow fault and two diffraction sets imaged within the fault zone that represent potential targets for future exploration. The south seismic reflection section shows rather a complicated geometry of two fault systems. The Guercheville fault observed as a subhorizontal reflector connects to a steeply dipping reflector. The Doda fault dips subvertical in the shallow crust but as a steeply dipping reflection set at depth. Nearby gold showings suggest that these faults may help channel and concentrate mineralizing fluids.

Potential-field modelling of the prospective Chibougamau area (northeastern Abitibi subprovince, Quebec, Canada) using geological, geophysical, and petrophysical constraints

This paper focusses on obtaining a better understanding of the subsurface geology of the Chibougamau area, in the northeast of the Abitibi greenstone belt (Superior craton), using geophysical data collected along a 128 km long traverse with a rough southwest–northeast orientation. We have constructed two-dimensional (2D) models of the study area that are consistent with newly collected gravity data and high-resolution magnetic data sets. The initial models were constrained at depth by an interpretation of a new seismic section and at surface by the bedrock geology and known geometry of lithological units. The attributes of the model were constrained using petrophysical measurements so that the final model is compatible with all available geological and geophysical data. The potential-field data modelling resolved the geometry of plutons and magnetic bodies that are transparent on seismic sections. The new model is consistent with the known structural geology, such as open folding, and provides an improvement in estimating the size, shape, and depth of the Barlow and Chibougamau plutons. The Chibougamau pluton is known to be associated with Cu–Au magmatic-hydrothermal mineralisation and, as the volume and geometry of intrusive bodies is paramount to the exploration of such mineralisation, the modelling presented here provides a scientific foundation to exploration models focused on such mineralisation.