Search

Arc magmas globally are H2O-Cl-S–rich and moderately oxidized (ΔFMQ = +1 to +2) relative to most other mantle-derived magmas (ΔFMQ ≤ 0). Their relatively high oxidation state limits the extent to which sulfide phases separate from the magma, which would otherwise tend to deplete the melt in chalcophile elements such as Cu (highly siderophile elements such as Au and especially platinum-group elements are depleted by even small amounts of sulfide segregation). As these magmas rise into the crust and begin to crystallize, they will reach volatile saturation, and a hydrous, saline, S-rich, moderately oxidized fluid is released, into which chalcophile and any remaining siderophile metals (as well as many other water-soluble elements) will strongly partition. This magmatic-hydrothermal fluid phase has the potential to form ore deposits (most commonly porphyry Cu ± Mo ± Au deposits) if its metal load is precipitated in economic concentrations, but there are many steps along the way that must be successfully negotiated before this can occur. This paper seeks to identify the main steps along the path from magma genesis to hydrothermal mineral precipitation that affect the chances of forming an ore deposit (defined as an economically mineable resource) and attempts to estimate the probability of achieving each step. The cumulative probability of forming a large porphyry Cu deposit at any given time in an arc magmatic system (i.e., a single batholith-linked volcanoplutonic complex) is estimated to be ~0.001%, and less than 1/10 of these deposits will be uplifted and exposed at shallow enough depths to mine economically (0.0001%). Continued uplift and erosion in active convergent tectonic regimes rapidly remove these upper-crustal deposits from the geological record, such that the probability of finding them in older arc systems decreases further with age, to the point that porphyry Cu deposits are almost nonexistent in Precambrian rocks. A key finding of this paper is that most volcanoplutonic arcs above subduction zones are prospective for porphyry ore formation, with probabilities only falling to low values at late stages of magmatic-hydrothermal fluid exsolution, focusing, and metal deposition. This is in part because of the high threshold required in terms of grade and tonnage for a deposit to be considered economic. Thus, the probability of forming a porphyry-type system in any given arc segment is relatively high, but the probability that it will be a large economic deposit is low, dictated to a large extent by mineral economics and metal prices.

The Tongkuangyu copper deposit in the Zhongtiaoshan region, southern Trans-North China orogen, is hosted by a poorly constrained sequence of Paleoproterozoic volcano-sedimentary (quartz-sericite schist and biotite schist) and granitic rocks that have been metamorphosed to lower greenschist facies and variably deformed. The deposit has previously been proposed to be either a porphyry-type or a sediment-hosted stratiform Cu deposit, and its age of formation has been debated. The quartz-sericite schist is interpreted to be a felsic crystal tuff and consists of angular quartz crystals in a fine-grained sericite-altered matrix. Two quartz-sericite schist samples yielded zircon U-Pb upper concordia intercept ages of 2512 ± 12 (2σ, mean square of weighted deviates [MSWD] = 0.19) and 2335 ± 16 Ma (2σ, MSWD = 0.80). Biotite schist, which is interleaved locally with the quartz-sericite schist and is interpreted to be a basaltic-andesitic sill, yielded a younger zircon U-Pb upper concordia intercept age of 2191 ± 10 Ma (2σ, MSWD = 1.7). Five samples of granodiorite and granodiorite porphyry that intruded the schist sequence yielded similar zircon U-Pb ages, with a weighted mean upper concordia intercept age of 2182 ± 7 Ma (2σ, MSWD = 1.3). These results suggest that the volcanic sequence was deposited between ∼2.5 and 2.3 Ga and was intruded by basaltic-andesitic sills and a suite of granodiorite and granodiorite porphyry intrusions at ∼2.19 to 2.18 Ga. Two stages of copper mineralization are interpreted to have formed after pervasive sericite alteration of the felsic volcanic rocks. Stage 1 mineralization includes disseminated and deformed quartz veinlets containing chalcopyrite ± pyrite ± magnetite ± molybdenite associated with biotite ± K-feldspar alteration in granodiorite porphyry and schist. Stage 2 comprises undeformed quartz-chlorite-carbonate veins with bornite ± chalcopyrite ± magnetite associated with local chlorite and silicic alteration. Allanite crystals intergrown with chalcopyrite in the granodiorite porphyry yielded an approximate concordia U-Pb age of 2115 ± 31 Ma (2σ, MSWD = 2.3). Two molybdenite samples in a deformed quartz-chalcopyrite-molybdenite vein yielded Re-Os model ages of 2106 ± 9 and 2089 ± 9 Ma (2σ), consistent with previously published results. Hydrothermal monazite grains with Cu-Fe sulfide inclusions in the granodiorite porphyry, quartz-sericite schist, and undeformed chlorite-bornite-quartz veins yielded much younger U-Pb upper concordia intercept ages of 1832 ± 16 (2σ, MSWD = 0.48), 1810 ± 14 (2σ, MSWD = 0.92), and 1809 ± 12 Ma (2σ, MSWD = 0.38), respectively. The results are in agreement with four Re-Os model ages for pyrite mineral separates from undeformed quartz-sulfide veins, which yielded a weighted mean age of 1807 ± 4 Ma (2σ, n = 4, MSWD = 0.42). In contrast, hydrothermal rutile crystals in the quartz-sericite schist and biotite schist yielded a range of roughly concordant ages between 2.1 and 1.8 Ga, reflecting isotopic disturbance. We interpret these results to indicate original copper mineralization at ∼2.1 Ga that is significantly later than the granodiorite (∼2.18 Ga) and schists (∼2.5–2.2 Ga), followed by hydrothermal remobilization and metamorphism at ∼1.8 Ga. The metavolcanic and granodiorite porphyry host rocks, alteration styles, and disseminated and veinlet form of the earlier mineralization are strongly reminiscent of porphyry Cu deposits, and ages of ∼2.1 Ga have been reported for one intrusion and three volcanic rock samples from the district. The Tongkuangyu, therefore, represents one of the oldest known porphyry copper deposits. Remobilization of copper occurred at ∼1.8 Ga during the Zhongtiao orogeny.

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.

The Yulong porphyry Cu-Mo ore district in eastern Tibet features a series of Eocene felsic porphyry intrusions, only one of which is intensively mineralized. In this study, zircon U-Pb ages, and major and trace element compositions of whole rock samples and in-situ accessory minerals (zircon, apatite and titanite) from these various intrusions were analyzed to determine the factors controlling Cu mineralization. New zircon U-Pb dating, together with published zircon U-Pb ages, suggest that these porphyry intrusions were emplaced over a relatively short period of time (43.9 ± 0.6 Ma to 40.9 ± 0.3 Ma), with the mineralized Yulong intrusion emplaced at a late stage of the magmatic activity (41.1 ± 0.3 Ma to 40.9 ± 0.3 Ma). These intrusions have indistinguishable whole-rock major and trace element compositions, but the mineralized Yulong intrusion has higher apatite and titanite La/Yb ratios (average apatite La/Yb = 115.1 ± 28.5, n = 52; average titanite La/Yb = 22.4 ± 4.9, n = 38) than the earlier subeconomic intrusions (average apatite La/Yb = 50.4 ± 12.3, n = 72; average titanite La/Yb = 12.4 ± 3.5, n = 58), suggesting more fractionation of amphibole from the magmas that sourced it. Additionally, a negative relationship between La/Yb and CeN/CeN* ratios in apatite and titanite are interpreted to reflect increasing oxidation states through magma evolution. Indistinguishable zircon EuN/EuN* values and apatite core SO3 contents from the mineralized (zircon EuN/EuN* = 0.69 ± 0.06, n = 25; apatite SO3 = 0.65 ± 0.24 wt%, n = 21) and subeconomic intrusions (zircon EuN/EuN* = 0.62 ± 0.05, n = 64; apatite SO3 = 0.61 ± 0.27 wt%, n = 15) suggest that the entire magmatic suite was relatively oxidized, hydrous and sulfur-rich, and was therefore fertile for ore-formation. However, apatite crystals from the mineralized Yulong intrusion have significantly lower Cl contents (0.08 ± 0.03 wt%, n = 72) and higher F/Cl ratios (47.80 ± 21.97, n = 72) than those from the subeconomic intrusions (Cl = 0.29 ± 0.29 wt%; F/Cl = 12.87 ± 5.16; n = 78). These data may reflect more extensive fluid exsolution from the parental magma chamber at the time of emplacement of the mineralizing Yulong magmas. Sudden and voluminous release of volatiles coeval with emplacement of the Yulong intrusion was triggered by an as yet unidentified process, but possibly including magmatic recharge.

Platinum-group element (PGE) geochemistry may be used to constrain the timing of sulfide saturation in magmas, which influences the Cu and Au fertility of evolving magmatic systems. We report new geochronological and geochemical data, with emphasis on PGE geochemistry, for a suite of regional porphyritic hornblende–diorite intrusions and ore-bearing porphyries from the super-giant Escondida and smaller Zaldivar Cu deposits of Northern Chile. The regional dioritic intrusions have zircon U–Pb ages between 39·6 to 37·1 Ma, which overlap with the ages of the ore-bearing Escondida and Zaldivar porphyries (38·1 to 35·0 Ma). Whole-rock major and trace element, and Sr–Nd–Pb and zircon O–Hf isotope geochemistry indicate that the regional diorites and ore-bearing porphyries are co-magmatic and originated from the same mantle-derived magma by fractional crystallization, with minor contamination by Paleozoic crust (∼10%). The low concentrations of PGE in the regional diorites show that they reached sulfide saturation before the MgO content of the melt fell to 4·7 wt %, the MgO content of the most primitive sample analysed. The fraction of sulfide melt which separated from the melts that formed the regional diorites is estimated to be ∼0·12 wt %; this resulted in the partitioning of highly chalcophile elements (Au and PGE) into a sulfide phase that was retained in cumulus rocks at depth. However, the fraction of sulfide melt was too low to have a significant effect on the Cu content of the fractionating melt. As a consequence, when the evolving melt eventually reached volatile saturation, it contained enough Cu (40 ± 10 ppm) to form a super-giant Cu deposit. In contrast, Au was largely stripped from the melt by sulfide precipitation, with the result that the mineralization at Escondida is Cu dominant, with only minor Au. The Zaldivar deposit, on the other hand, contains even less Au, which is attributed to a longer fractionation interval between sulfide and volatile saturation. This study provides evidence to support previously proposed models which suggest that the timing of sulfide saturation, the amount of sulfide melt produced, the water content and oxidation state of the melt, and the magma volume are critical factors in determining the potential to form a porphyry Cu deposit. Plots of Pd/MgO against Y can be used as empirical indicators of magma fertility for porphyry mineralization, and to discriminate between Cu–Au and Cu-dominated systems, but cannot predict the size of the deposit. The super-giant status of the Escondida deposit is attributed to it being underlain by a large batholith with a calculated minimum mass of 10 12 tonnes (∼400 km 3 ).

This work was supported by the National Key Research and Development Program, China (2016YFC0600206), and the National Science Foundation of China (41372090 and 41573042).

Jurassic volcanic rocks are widely distributed along the Gangdese belt of southern Tibet, forming a westerly-extending southern belt (Bima Formation) with one large arc-related porphyry Cu Au deposit (Xietongmen), and an easterly-extending northern belt (Yeba Formation) without any contemporary mineralization. A compilation of new and published data on the Bima Formation is compared to published data on the Yeba Formation in order to understand their different relationships with mineralization. The Bima Formation volcanic rocks comprise basalt, andesite, and dacite with zircon U Pb ages of 195.0 Ma–165.1 Ma (including new ages of 180.2 ± 0.8 Ma to 166.7 ± 1.1 Ma). They are calc-alkaline series lava flows and show enrichments in LILE and LREE, depletions in HFSE and HREE, which are features typical of arc rocks. Coherent geochemical variations and similar depleted Sr-Nd-Hf isotopic characteristics between basaltic, intermediate, and felsic rocks (new and published data: (87Sr/86Sr)i = 0.703182 to 0.705489); eNd(t) = 3.2 to 7.1; eHf(t) = 9.3 to 16.0) suggest that they were formed from fractional crystallization of the basaltic magma. In contrast, the published Yeba Formation volcanic rocks, although with a similar age range (192.7 Ma to 168.0 Ma), have a wider compositional range (from basalt to rhyolite), and more evolved Sr-Nd-Hf isotopic compositions ((87Sr/86Sr)i = 0.700182 to 0.707723; eNd(t) = −5.1 to 4.5; eHf(t) = −18.5 to 17.6). New data on the Bima Formation andesite and dacite samples have higher zircon Eu/Eu*, Ce4+/Ce3+ and (Ce/Nd)/Y ratios, and lower titanium-in-zircon temperatures than those from the Yeba Formation, suggesting the Bima Formation magmas were more oxidized and hydrous, and therefore more fertile for porphyry copper deposit formation. The relatively reduced Yeba Formation magmas formed in a distal arc setting near the Central Lhasa old basement. Their evolved isotopic characteristics and reduced nature may have resulted from more extensive contamination by reduced crustal components, rendering them infertile for porphyry copper mineralization. Exploration for arc-related porphyry copper deposits in the Gangdese belt should therefore focus on the fertile Bima Formation.

The Cerro Vetas porphyry deposit is part of the Titiribi district of the Middle Cauca porphyry-epithermal belt of western Colombia. The Cerro Vetas porphyry stock consists of a premineral diorite intruded by a late-mineral quartz monzonite, with intrusion and contact breccias. These units intrude pre-Cenozoic basement metabasalts and schists, Oligocene-Miocene Amagá Formation sedimentary rocks with intercalated andesite flows. Two phases of potassic alteration are recognized, a biotite-dominant phase in the diorite, and secondary K-feldspar in the quartz-monzonite intrusion. An overprinting and grade destructive, calcic-sodic alteration (actinolite + albite ± magnetite) affects both porphyries. Biotite alteration is overprinted by weak-moderate phyllic alteration in the upper 100 m in the deposit. Below 100 m, phyllic alteration assemblages are constrained to structural zones. Mineralization is dominated by a chalcopyrite-gold-pyrite assemblage associated with biotite that is hosted in a truncated stockwork in the apical portion of the deposit with metal ratios typical of a gold-rich copper-gold porphyry. The intrusions were dated, using U-Pb in zircon laser ablation inductively coupled plasma-mass spectroscopy, to between 7.65 to 7.24 Ma, consistent with other deposits in the Middle Cauca belt. Lithologic, alteration, and stratigraphic relationships at the deposit suggest that the Cerro Vetas porphyry was emplaced at shallow depths and that the upper portion of the deposit has been eroded.

Superimposed porphyry systems are a subset of telescoped porphyry deposits, whereby significantly younger ore zones overprint older, nongenetically related systems. Recognition of superimposed features in porphyry systems is important for determining and assessing their prospectivity. The Mount Nansen gold corridor in the southern Dawson Range gold belt of Yukon, Canada, contains porphyry prospects and epithermal deposits with enigmatic genetic models. Geologic, petrologic, temporal (U-Pb zircon, Re-Os molybdenite), and geochemical (whole-rock) studies are used to demonstrate the presence of superimposed porphyry systems in this district. The arc-related episodic magmatism of the Mount Nansen gold corridor has been classified into four intrusive suites: (1) Late Triassic Minto, (2) mid-Cretaceous Whitehorse, (3) Late Cretaceous Casino (eLKc; 80–72 Ma), and (4) Late Cretaceous Prospector Mountain (lLKp; 72–65 Ma). Geochemical fingerprinting of these suites indicates intermediate to evolved, calc-alkaline compositions with a common lower crust melt source. The eLKc and lLKp suites lack an Eu anomaly and show increasing amounts of light rare earth element (LREE) enrichment and heavy rare earth element (HREE) depletion over time. These features suggest that garnet was stable in the melt source and oxidized magmas were generated in these Late Cretaceous suites. The mildly alkaline lLKp and associated Carmacks Group shoshonitic basalts reflect localized extension in an overall compressive arc setting in the Mount Nansen gold corridor, hence a setting conducive for Au-rich porphyry and epithermal systems. The ca. 79 to 72 Ma Casino suite is commonly interpreted as the causative magmatic event for most well-endowed porphyry deposits (76 to 74 Ma in age) in the Dawson Range gold belt. However, our detailed study of the Klaza setting shows that at this locality, intermediate-sulfidation epithermal veins are a distal expression of a Prospector Mountain-age (ca. 71 Ma) porphyry system, which overprints two Casino-age porphyry systems (ca. 77 and 80 Ma). The Mount Nansen gold corridor thus hosts at least two spatially and temporally overprinting Late Cretaceous magmatic-hydrothermal systems in the Dawson Range gold belt. Importantly, recognition of this feature at other porphyry deposit settings in the Dawson Range gold belt (e.g., Freegold Mountain district) is critical as it provides the potential for metal (Cu-Au-Mo)-enriched hypogene ore shells.

Some arc magmas lead to the formation of porphyry deposits in the relatively shallow upper crust (

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 fayalite-magnetite-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., Tectonomagmatic conditions in the Precambrian were hypothesized to be unfavorable for porphyry Cu deposit formation. Here, the authors show that metallogenic processes typify Phanerozoic porphyry Cu deposits operated by ~1.88 Ga, reflecting modification of mantle lithosphere by oxidized slab-derived fluids at that time.

Numerous W-Sn deposits are developed in the South China Block with the majority related to the Yanshanian (Jurassic to Cretaceous) granites; a smaller number of deposits have been found associated with Kwangsian (middle Paleozoic) granites, the largest of which is the Niutangjie skarn W deposit. The Niutangjie deposit is spatially associated with two suites of granites: the unmineralized Yuechengling biotite and two-mica granites have zircon U-Pb ages of 427.1 ± 2.9 Ma and 427.5 ± 3.5 Ma, respectively, whereas the mineralized Niutangjie muscovite granite yields a younger age of 421.0 ± 1.5 Ma. U-Pb dating of cassiterite coexisting with scheelite from the Niutangjie muscovite granite indicates a mineralization age of 420.8 ± 8.0 Ma, consistent with the age of the host rock. The unmineralized Yuechengling biotite granites have A/CNK ratios of 0.97–1.16 and show characteristics of both I-type and S-type granites. In contrast, the mineralized Niutangjie muscovite granites are S-type granites with A/CNK ratios of 1.12–1.72 and have higher SiO2 contents and Rb/Sr ratios, and lower Nb/Ta, Zr/Hf, and K/Rb ratios than the unmineralized Yuechengling biotite granites. Geochemical differences together with zircon Hf isotopic compositions indicate that the relatively fractionated mineralized Niutangjie muscovite granites were generated by partial melting of Proterozoic meta-sediments, whereas the unmineralized Yuechengling biotite granites were derived from partial melting of Proterozoic meta-sediments with subordinate meta-igneous rocks. The unmineralized Yuechengling two-mica granites have similar geochemical characteristics to the mineralized Niutangjie muscovite granites, but their relatively small volumes and lack of exploration may explain why the Yuechengling two-mica granites have not so far been found to be mineralized. Similarly, the larger volume and more highly fractionated nature of the Mesozoic Yanshanian mineralized granites may explain the greater abundance of W-Sn mineralization associated with these rocks compared to the Paleozoic Kwangsian granites.

Mineralization in the Zacatecas district is polymetallic (Ag, Zn, Pb, Cu, and Au) and occurs as skarn-type and epithermal veins formed in different metallogenetic stages. The oldest mineralization in the district is skarn-type, Cu-rich with lesser Zn-Pb-Ag, and is considered to be close in age to felsic dikes and plugs dated at ∼51 Ma. Epithermal mineralization occurs in both low- and intermediate-sulfidation styles. Intermediate-sulfidation veins (the Veta Grande, Mala Noche, El Bote, and La Cantera veins) are polymetallic, Ag-rich, hosted in ESE- to SE-striking structures, and were formed at ∼42 Ma (adularia 40Ar/39Ar isochron age from Veta Grande of 42.36 ± 0.18 Ma; 2σ, MSWD = 0.76). Low-sulfidation Au-(Ag) mineralization occurs in the N–S-trending El Orito vein system, which yielded an adularia 40Ar/39Ar isochron age of 29.19 ± 0.20 Ma (2σ, MSWD = 1.8). These ages and the differences in structural orientation indicate that the two styles of epithermal mineralization are temporally and tectonically unrelated. The mineral paragenesis of the Mala Noche deposit consists of early skarn-type Cu mineralization overprinted by later epithermal Pb-Zn-Ag veins. Skarn-type minerals include relicts of prograde silicate minerals (hedenbergite and garnets), retrograde silicate minerals (ilvaite, grunerite, stilpnomelane, epidote, clinochlore), and ore minerals (chalcopyrite, pyrite, sphalerite, galena, magnetite, wolframite, and minor bismuthinite). Epithermal mineralization is characterized by layered to vuggy quartz veins and breccias, with phyllic wallrock alteration (quartz, sericite-illite). The veins consist of quartz, calcite, dolomite, and ankerite with variable amounts of base metal sulfides (sphalerite, galena, pyrite, minor chalcopyrite, and rare acanthite and stromeyerite). The Veta Grande epithermal mineralization was emplaced in two main stages of Ag-rich quartz veining, with narrow selvedges of phyllic (quartz-sericite) wallrock alteration. Stage I consist of quartz, calcite, and minor adularia intergrown with pyrite, followed by sphalerite, galena, and lesser chalcopyrite, acanthite, pyrargyrite, and jamesonite. Stage II mineral paragenesis is similar to stage I but is characterized by amethystine quartz and contains less abundant sulfide minerals. The ore mineral paragenesis of the El Compas vein, within the El Orito System, consists of quartz, adularia, calcite, and chalcedony with minor pyrite, followed by rare aguilarite, naumannite, electrum, and native gold. The skarn-type and intermediate-sulfidation mineralization is coeval with Eocene subduction-related magmatism in the Zacatecas area, which is constrained by zircon U-Pb ages for igneous rocks between 51 and 42 Ma. The emplacement of these magmas was controlled by the same regional-scale, ESE- to SE-trending, transtensional structures that controlled the skarn-type and intermediate-sulfidation deposits. This mineralization is thus interpreted to be related to the last stages of subduction volcanism in central Mexico, under transtensional stress conditions. In contrast, no nearby magmatism is clearly related to the Oligocene low-sulfidation system. However, its age and structural orientation (N–S), combined with a regional change in magma composition from Eocene calc-alkaline to Oligocene bimodal volcanism in central Mexico, suggest that the low-sulfidation mineralization is related to post-subduction continental extension processes, reflecting the beginning of Basin and Range tectonics.

The Newton Au-Ag deposit is an intermediate sulfidation state epithermal system in British Columbia, Canada. Multiple types of geophysical data are interpreted and evaluated with drillcore petrophysical, geochemical, and geologic observations to better understand the geophysical signature of the Newton epithermal system. Airborne [Formula: see text]-ray data sets indicate elevated emission counts of K, eTh, and eU over the Newton epithermal system, which are caused by hydrothermal alteration. Drillcore [Formula: see text]-ray measurements also indicate high potassium concentrations related to the K-rich phyllosilicates in the form of argillic and quartz-sericite alteration assemblages. Magnetization vector inversion (MVI) is used to recover an unconstrained 3D magnetization vector model of the system on regional and deposit scales. The regional MVI has resolved a deep concentric-shaped low magnetic zone that is interpreted as a porphyry system beneath the epithermal deposit. At the deposit scale, 3D direct current (DC) resistivity and induced polarization (IP) inversion, and unconstrained MVI revealed finer details of the epithermal system architecture. Cooperative DC/IP and magnetic inversion, at the deposit scale, constrained the magnetic susceptibility model and recovered a more precise susceptibility image of the epithermal system that is well-matched with borehole geology. The integrated geophysical interpretation helped to resolve several 3D latent geologic features in places without direct access to drillcore samples. We identified four petrophysical domains based on the three cooperatively inverted physical properties, including electrical resistivity, IP chargeability, and magnetic susceptibility. The combined geophysical models differentiated porphyritic intrusions (chargeability/susceptibility lows), disseminated sulfides (resistivity lows and chargeability highs), a Cu-rich zone in mafic volcanic rocks (susceptibility/chargeability highs and resistivity lows), and an Au-Ag-Cu-rich zone with silicification in felsic volcanic rocks (chargeability/susceptibility lows and resistivity highs). These petrophysical domains also provide useful exploration vectors for identification of similar epithermal systems.

The Red Chris porphyry Cu-Au deposit is located in the Stikinia island-arc terrane in northwest British Columbia. It is hosted by the Red Stock, which has five phases of porphyry intrusions: P1, P2E, P2I, P2L, and P3. New U-Pb dating of zircon shows that these intrusions were emplaced over a similar to 10 m.y. period, with P1 intruded at 211.6 +/- 1.3 Ma (MSWD = 0.85), P2I at 206.0 +/- 1.2 Ma (MSWD = 1.5), P2L at 203.6 +/- 1.8 Ma (MSWD = 1.5), and P3 at 201.7 +/- 1.2 Ma ( MSWD = 1.05). The ore-forming event at Red Chris was a short-lived event at 206.1 +/- 0.5 Ma (MSWD = 0.96; weighted average age of three Re-Os molybdenite analyses), implying a duration of Zircons from P1 to P3 porphyry rocks have consistently high Eu-N /Eu-N* ratios (mostly > 0.4), indicating that their associated magmas were moderately oxidized. The magmatic water contents estimated from plagioclase and amphibole compositions suggest H2O contents of similar to 5 wt %. Taken together, the P1 to P3 porphyries are interpreted to be moderately oxidized and hydrous. The porphyry phases are differentiated by sulfur and chlorine contents. The SO3 contents of igneous apatite microphenocrysts from the mineralization-related P2 porphyries are higher (P2E: 0.28 +/- 0.10 wt %, n = 19; P2I: 0.32 +/- 0.17, n = 15; P2L: 0.29 +/- 0.18 wt %, n = 100) than those from the premineralization P1(0.11 +/- 0.03 wt %, n = 34) and postmineralization P3 porphyries (0.03 +/- 0.01 wt %, n = 13). The chlorine contents in apatite grains from the P2E, P2I, and P2L porphyries are 1.47 +/- 0.22 (n = 19), 0.82 +/- 0.10 (n = 15), and 1.47 +/- 0.28 wt % (n = 100), also higher than those from P1 (0.51 +/- 0.3 wt % Cl, n = 34) and P3 (0.02 +/- 0.02 wt % Cl, n = 17). These results imply that the sulfur and chlorine contents of the P2 magmas were higher than in the P1 and P3 magmas, suggesting that elevated magmatic S-Cl contents in the P2 porphyries may have been important for ore formation. Although the process that caused the increase in sulfur and chlorine is not clear, reverse zoning seen in plagioclase phenocrysts from the P2 porphyry, and the occurrence of more mafic compositions in P2L, suggest that recharge of the deeper magma chamber by a relatively S-Cl-rich mafic magma may have triggered the ore-forming hydrothermal event.

The recent discovery of large porphyry copper deposits (PCDs) associated with Miocene (22–12 Ma) granitoid magmas in the eastern section of the Paleocene-Eocene Gangdese magmatic arc in the Himalaya-Tibetan orogenic belt raises new questions about the origin of water-rich (≥4.5 wt.%), oxidized (ΔFMQ 1–3) magmas in continental collisional settings and their mineralization potential. We review the literature and compile available data on whole rock and isotope geochemistry for Cenozoic igneous rocks from Tibet, and add new zircon Ce4+/Ce3+ and Ti-in-zircon thermometry data to better understand variations in oxidation state and thermal evolution of these suites, which are key controls on Cu mineralization. Six distinct Cenozoic igneous suites are defined: Paleocene-Eocene syn-collisional Gangdese magmatic arc rocks (ΔFMQ = −1.2 to +0.8) (suite I), and five broadly contemporaneous Miocene suites. A distinct change in magmatism along the length of the belt occurs at around 88°E in the Miocene suites: to the east, porphyry copper mineralization is associated with a moderately oxidized, high-Sr/Y granitoid suite (suite II, ΔFMQ = +0.8 to +2.9) with minor occurrences of transitional (hybrid) monzonitic (suite III) and trachytic rocks (suite IV; both with zircon Ce4+/Ce3+ > 50–100, EuN/EuN* = ~0.5, and ΔFMQ = ~+1 to +2). To the west of 88°E, trachytic volcanic rocks (suite V) are more voluminous but more reduced (zircon Ce4+/Ce3+ The Miocene high-Sr/Y granitoids have many compositional and isotopic similarities to the Paleocene-Eocene Gangdese arc rocks, and are interpreted to have been derived by melting of the hydrated arc root, with minor mantle input. In contrast, the highly evolved isotopic signatures of the Miocene trachytic rocks, combined with deep seismic profiles and a xenolith-derived geotherm, suggest their derivation from the underthrust Indian Proterozoic subcontinental lithospheric mantle (SCLM) or old fore-arc Tibetan SCLM during phlogopite breakdown at temperatures of ~1100°C. Based on published geophysical data and tectonic reconstructions, we develop a model that explains the origin of the various Miocene magmatic suites, their spatial differences, and the origin of related PCDs. Following the early stages of continental collision (Eocene–Oligocene), shallow underthrusting of the Indian continental lithosphere and subcretion of Tethyan sediments (including oxidized carbonates and possibly evaporites) under eclogite facies conditions promoted the release of aqueous fluids, which hydrated and oxidized the base of the overlying Tibetan plate. This metasomatism rendered the Tibetan lower crust fusible and fertile for metal remobilization. During the mid-Miocene, the Indian slab steepened in the eastern sector (east of ~88°E). In this eastern belt, deeply derived trachytic magmas were trapped in melt zones at the base of the Tibetan crust, and variably mixed with the crustally-derived, high Sr/Y granitoid magmas. They may also have released water that contributed to fluid-fluxed melting of the lower crust, producing voluminous high-Sr/Y granitoid magmas, which were associated with significant PCD mineralization. Hybridization between the trachytic magmas and lower crustal partial melts is indicated by intermediate isotopic compositions, enriched Cr and Ni contents, and high Mg# in some intermediate-to-felsic (56–70 wt. % SiO2) high-Sr/Y granitoids. Trapping of the trachytic melts in deep crustal melt zones explains the relatively small volumes of trachytic magmas erupted at surface in the east. In contrast, to the west of ~88°E, subduction of the Indian plate has remained flat to the present day, preventing incursion of hot asthenosphere. Consequently, cooler conditions in the deep Tibetan lithosphere resulted in limited crustal melting and the production of only small volumes of high-Sr/Y granitic magmas. Trachytic melts ascending from the underthrust Indian or Tibetan plate were able to pass through the cooler lower crust and erupted in greater volume at surface, whereas only small volumes of high-Sr/Y granitoid magma were generated and are not associated with significant PCD mineralization.

The high-grade (8.1 million ounces of gold; 15.6 Mt grading 16.1 g/t Au; Pretium, 2016) Brucejack epithermal Au-Ag deposits are located in the Canadian Cordillera of northwestern British Columbia, and formed in association with extensive early Mesozoic island arc magmatism. Porphyry-type Cu-Au-Mo mineralization occurs nearby at the Kerr-Sulphurets-Mitchell (196–190 Ma; Febbo et al., 2015), Bridge Zone (191.7 ± 0.8 Ma, 191.5 ± 0.8 Ma), and West Zone (188.9 ± 0.9 Ma) prospects. Gold-silver vein-type mineralization at Brucejack is hosted by variably altered and deformed Early Jurassic porphyritic latite lava flows, volcaniclastic rocks, and volcanic-derived sandstones, siltstones, and conglomerates. This study focuses on the Valley of the Kings Zone at Brucejack, where host rocks have been dated at 188–184 Ma (U-Pb, zircon). Post-mineralization (but not post-alteration) basaltic and trachybasaltic dikes that cut the veins are geochemically similar to the host volcanic sequence and follow the same structural trends of the mineralized vein corridors. We suggest these dikes are broadly coeval with mineralization, and that ore formation at the Valley of the Kings formed during the overall Jurassic tectonomagmatic event. Gold-silver mineralization is hosted by quartz-carbonate veins that cut sericitized and pyritized volcaniclastic rocks. Six stages of veining are defined, with Au-Ag mineralization (electrum) focused in stages III–V: stage I–III veins consist of quartz with minor carbonate, chlorite, sericite, and sulfide minerals; stage IV quartz veins contain more abundant base metal sulfides (pyrite, sphalerite, and galena with minor sulfosalts); stage V veins are dominantly calcitic; late post-mineralization stage VI veins are quartz-calcite and contain sparse pyrite and chlorite, but no electrum. The veins are interpreted to be syn- to late-tectonic, with deformation decreasing from locally penetrative in the host-rocks prior to veining (resulting in local foliation of sericite and pressure shadows around pyrite), to dismemberment and shearing of early stage I–II quartz veins, brittle disruption of stage III–IV quartz and stage V carbonate veins, and minimal deformation of post-mineralization stage VI quartz-carbonate veins. A much younger set of extensional muscovite veins locally cuts the deposit with Late Cretaceous apparent ages, and appears to be associated with a weak thermal overprint that has reset K-Ar and 40Ar/39Ar ages in sericite throughout the district. The deformation experienced by the vein system, particularly in the later stages (III–V), was not uniform, and original undeformed vein textures and fluid inclusion assemblages are locally preserved. We report fluid inclusion microthermometric data from carefully selected primary fluid inclusion assemblages from quartz, calcite, and sphalerite from vein stages III–IV, many of which show evidence for boiling (coexisting liquid- and vapor-rich primary fluid inclusions). Liquid-rich fluid inclusions from stage III and V veins have moderate homogenization temperatures (∼170 °C and ∼160 °C respectively) and salinities of 2–8 wt% NaCl equiv., whereas inclusions from base-metal-sulfide-bearing stage IV veins show evidence of mixing with a cooler, more saline brine (∼140 °C, ∼10–15 wt% NaCl equiv.). Carbon dioxide was observed as clathrate during cooling in some fluid inclusions, suggesting that minor amounts of CO2 were present in the fluids. Calculated oxygen isotopic compositions of fluids in equilibrium with quartz and calcite from vein stages III–V range from δ18Ofluid = −10.7 to +1.8‰, whereas δ13CCO2 ranges from −9.5 to −4.5‰, and δ34Spyrite ranges from −1.7 to +0.6‰. Taken together, these data suggest a magmatic source for S and some C, carried by a fluid of evolved or diluted magmatic origin, which variably mixed with meteoric-derived groundwater or seawater containing carbon of sedimentary (organic) origin. Mixing is supported by large salinity ranges of vein stage IV fluid inclusions, where individual assemblages range from ∼1 to 15 wt% NaCl equiv. These fluid temperatures, salinities, and isotopic compositions are typical of epithermal deposits where distal magmatic fluids mix with local heated groundwaters. The low fluid temperatures but evidence of boiling suggest formation at shallow crustal depths. This is despite the evidence for penetrative deformation in the sericitized volcanosedimetary host rocks, which we attribute to rapid uplift immediately prior to or during the earliest stages of mineralization. The vein system is therefore interpreted to have formed during the later stages of a deformational event related to arc accretion.

The Taftan volcanic complex is located above the Makran subduction zone in Sistan and Baluchestan province, southeastern Iran. The earliest volcanic activity at Taftan started in the late Miocene (∼8 Ma) with eruption of andesitic to dacitic lava onto a Cretaceous to Eocene volcanic and sedimentary paleosurface ∼20 km to the northwest of the current volcanic edifice. Later Plio-Pleistocene volcanism consisted of voluminous andesitic and dacitic lavas and pyroclastic flows (∼3.1 to 0.4 Ma). Taftan, and the neighbouring Bazman volcano, are the southeasternmost and youngest manifestations of arc-related volcanism in Iran, which began with the Eocene–Miocene Urumieh-Dokhtar Magmatic Arc (UDMA) in northwest and central Iran, and extends into western Pakistan (Chagai Arc). Major porphyry Cu±Mo±Au deposits are associated with Miocene to Pliocene intrusive rocks in the Kerman section of the Eocene–Neogene Urumieh-Dokhtar Magmatic Arc in southeastern Iran (e.g., Sar Cheshmeh, Meiduk), and at Reko Diq and Saindak in the Late Cretaceous–Quaternary Chagai Arc in western Pakistan. In a gap of ∼300 km between these two belts, several porphyry and epithermal prospects are exposed in the Miocene–Pliocene sections of the Taftan volcanic complex, including the Kharestan (6.10 ± 0.80 Ma) and Bidester porphyry-epithermal Cu-Au deposits (∼4.3 Ma), and the Siah Jangal epithermal Au deposit (late Miocene). In addition, several argillic to advanced argillic and fumarolic alteration zones occur in Plio-Pleistocene volcanic rocks around the current volcanic edifice. These deposits have received limited exploration. Electron microprobe analyses of amphibole (magnesio-hastingsite) phenocrysts and magnetite–ilmenite mineral pairs from the Taftan and Bazman volcanic rocks indicate a change of crystallization temperature and magmatic oxidation state from ∼1000 °C and ΔFMQ ≈ +1 in andesitic rocks, to ∼900 °C and ΔFMQ ≈ +2 in dacitic rocks. Magmatic water content was >4 wt%, as indicated by the ubiquitous presence of amphibole phenocrysts. Major and trace element compositions of the Taftan and Bazman volcanic rocks show calc-alkaline to high-K calc-alkaline affinity, with relative depletions of Nb, Ta, and Ti and enrichments of large-ion c elements (LILE), Th, and U on normalized diagrams. These trace element patterns, including listric-shaped normalized rare earth element profiles and relatively high Sr/Y and La/Yb ratios, are similar to those of fertile Miocene igneous rocks from the Kerman Belt and Chagai Arc, and suggest that the Taftan suite in particular should be prospective for porphyry Cu ore formation. Regionally, there is no clear geochemical difference between the Neogene Kerman Belt rocks, which are thought to post-date the onset of collision between the Afro-Arabian and Eurasian plates (late Oligocene to earliest Miocene), and the subduction-related Bazman–Taftan and Chagai Belt magmas. The porphyry deposits formed in these distinct settings are also virtually indistinguishable. This suggests that most of the processes affecting the geochemistry and metallogeny of the magmas in both settings take place in the lithosphere, albeit that the ultimate source of the magmas is in the supra-subduction zone asthenospheric mantle wedge. In collisional environments, subduction-related material previously crystallized in the deep lithosphere is simply being remobilized.

Z-axis tipper electromagnetic and broadband magnetotelluric data were used to determine three-dimensional electrical resistivity models of theMorrison porphyry Cu– Au–Mo deposit in British Columbia. Z-axis tipper electromagnetic data are collected with a helicopter, thus allowing rapid surveys with uniform spatial sampling. Groundbased magnetotelluric surveys can achieve a greater exploration depth than Z-axis tipper electromagnetic surveys, but data collection is slower and can be limited by difficult terrain. The airborne Z-axis tipper electromagnetic tipper data and the ground magnetotelluric tipper data show good agreement at the Morrison deposit despite differences in the data collection method, spatial sampling, and collection date. Resistivity models derived from individual inversions of the Z-axis tipper electromagnetic tipper data and magnetotelluric impedance data contain some similar features, but the Z-axis tipper electromagnetic model appears to lack resolution below a depth of 1 km, and the magnetotelluric model suffers from non-uniform and relatively sparse spatial sampling. The joint Z-axis tipper electromagnetic inversion solves these issues by combining the dense spatial sampling of the airborne Z-axis tipper electromagnetic technique and the deeper penetration of the lower frequency magnetotelluric data. The resulting joint resistivity model correlates well with the known geology and distribution of alteration at the Morrison deposit. Higher resistivity is associated with the potassic alteration zone and volcanic country rocks, whereas areas of lower resistivity agree with known faults and sedimentary units. The pyrite halo and 0.3% Cu zone have the moderate resistivity that is expected of disseminated sulphides. The joint Z-axis tipper electromagnetic inversion provides an improved resistivity model by enhancing the lateral and depth resolution of resistivity features compared with the individual Z-axis tipper electromagnetic and magnetotelluric inversions. This case study shows that a joint Z-axis tipper electromagnetic–magnetotelluric approach effectively images the interpreted mineralised zone at the Morrison deposit and could be beneficial in exploration for disseminated sulphides at other porphyry deposits.

The Tethyan orogen is host to numerous porphyry Cu ± Mo ± Au deposits, but the majority formed during subduction of the Neo-Tethyan ocean basin in the late Mesozoic–Cenozoic; very few deposits have been found associated with Paleo-Tethyan subduction. We propose that this sparsity is due to widespread anoxia in the Paleo-Tethyan ocean basin, leading to the generation of relatively reduced arc magmas that were infertile for porphyry Cu formation. A compilation of published geochemical data indicates that Neo-Tethyan arc rocks have higher average Cu contents and V/Sc and Sr/Y ratios compared to Paleo-Tethyan rocks, indicating higher magmatic oxidation states and greater fertility for ore formation during Neo-Tethyan subduction. Subduction of relatively reduced oceanic lithosphere, or reduction of normal moderately oxidized arc magmas by interaction with reduced lithosphere, can therefore destroy the ore-forming potential of arc magmatic suites.

The Youjiang basin in southwest China is distinguished by preserving a number of Carlin-type gold deposits. Tens of felsic dikes occur in this region, and some are present within the ore fields. A genetic link between these felsic dikes and ore-formation has previously been proposed, but detailed field investigation shows that the dikes crosscut the main ore body, and clearly postdate gold mineralization. Three felsic dikes in the Youjiang basin have been dated, and SIMS zircon U-Pb dating results show that the Liaotun and Xiabaha dikes were emplaced at 97.2 ± 1.1 Ma (MSWD = 2.9) and 95.4 ± 2.4 Ma (MSWD = 0.56), respectively, broadly consistent with a LA-ICP-MS zircon U-Pb age for the Bama dike (99.4 ± 0.37 Ma, MSWD = 1.4). Combined with published dates of the other felsic dikes in this region, it is suggested that the felsic magmatic event in the Youjiang basin occurred between ∼100 and 95 Ma. Although a wide range of ages for gold mineralization (275–46 Ma) has been reported by a variety of methods, the best available ages suggest a timing between 235 and 193 (Direct Re-Os dating on Fe-sulfides; Chen et al., 2015, Re–Os isochron ages for arsenopyrite from Carlin-like gold deposits in the Yunnan–Guizhou–Guangxi “golden triangle”, southwestern China), significantly older than the emplacement age of the felsic dikes. Taken together, we conclude that there is no genetic link between these Late Cretaceous felsic dikes and Carlin-type Au mineralization in the Youjiang basin.

Porphyry Cu ± Mo ± Au and iron oxide copper-gold (IOCG) deposits share many similarities (e.g., Fe, Cu, and Au contents), but also have important differences (e.g., the predominance of sulfide minerals in porphyry deposits and iron oxides in IOCG deposits). Genetic comparisons are complicated by the broad definition of IOCG deposits; here we restrict our study to IOCG deposits that are related to igneous intrusive systems. In the Mesozoic Coastal Cordillera of northern Chile, both porphyry and IOCG deposits occur in close spatial and temporal proximity, offering the chance to examine what controls their different modes of formation. From detailed examination of the timing, geochemistry, and tectonic setting of associated igneous rocks, based on new and published data, we find that rocks associated with mid-Cretaceous IOCG deposits (~125–110 Ma) are largely indistinguishable from those associated with slightly earlier (>125 Ma) and later (

The Morrison porphyry Cu–Au–Mo deposit is genetically and spatially related to Eocene plagioclase–hornblende–biotite porphyry intrusions. One porphyry intrusion yielded a U–Pb age of 52.54 ± 1.05 Ma. Mineralization occurs in three stages: (1) vein-type and disseminated chalcopyrite and minor bornite (associated with potassic alteration and gold mineralization); (2) vein-type molybdenite (associated with weak phyllic alteration); and (3) polymetallic sulfide–carbonate veins (dolomite ± quartz–sphalerite–galena–arsenopyrite–chalcopyrite, associated with weak sericite–carbonate alteration). Re–Os dating of molybdenite yielded ages of 52.54 ± 0.22 and 53.06 ± 0.22 Ma, similar to the age of the host porphyry intrusion. Stage 1 vein fluids were predominantly of magmatic origin: Th = 400–526 °C; salinity = 39.8–47.8 wt.% NaCl equiv.; δ18Ofluid = 3.7‰–6.3‰; disseminated chalcopyrite–pyrite δ34SCDT = 0.2‰ and −0.8‰ (CDT, Canyon Diablo Troilite). Stage 2 fluids were a mixture of magmatic and meteoric water: Th = 320–421 °C; salinity = 37.0–43.1 wt.% NaCl equiv.; δ18Ofluid values range from 0.3‰ to 3.4‰; molybdenite and pyrite δ34SCDT = −2.1‰ and −1.2‰. Stage 3 fluids were predominantly of meteoric water origin: Th = 163–218 °C; salinity = 3.1–3.9 wt.% NaCl equiv.; δ18Ofluid = −2.3‰ to 3.9‰ for early vein quartz, and 1.1‰ to 6.1‰ for late vein dolomite; sphalerite and pyrite δ34SCDT = −7.1‰ to −5.6‰. Morrison is interpreted to be a typical porphyry Cu–Au–Mo deposit related to a calc-alkaline to a high-K calc-alkaline diorite to granodiorite intrusive suite, generated in a continental arc in response to early Eocene subduction of the Kula–Farallon plate beneath North America.

Newton is an intermediate-sulfidation epithermal gold deposit related to Late Cretaceous continental-arc magmatism in south-central British Columbia. Disseminated gold mineralization occurs in quartz–sericite-altered Late Cretaceous felsic volcanic rocks, and feldspar–quartz–hornblende porphyry and quartz–feldspar porphyry intrusions. The mineralization can be divided into three stages: (1) disseminated pyrite with microscopic gold inclusions, and sparse quartz–pyrite ± molybdenite veins; (2) disseminated marcasite with microscopic gold inclusions and minor base-metal sulfides; and (3) polymetallic veins of pyrite–chalcopyrite–sphalerite–arsenopyrite. Re–Os dating of molybdenite from a stage 1 vein yielded an age of 72.1 ± 0.3 Ma (published by McClenaghan in 2013). The age of the host rocks has been constrained by U–Pb dating of zircon: Late Cretaceous felsic volcanic rocks, 72.1 ± 0.6 Ma (Amarc Resources Ltd., unpublished data, reported by McClenaghan in 2013); feldspar–quartz–hornblende porphyry, 72.1 ± 0.5 Ma; quartz–feldspar porphyry, 70.9 ± 0.5 Ma (Amarc Resources Ltd., unpublished data, reported by McClenaghan in 2013). The mineralized rocks are intruded by a barren diorite, with an age of 69.3 ± 0.4 Ma. Fluid inclusions in quartz–pyrite ± molybdenite ± gold veins yielded an average homogenization temperature of 313 ± 51 °C (number of samples, n = 82) and salinity of 4.8 ± 0.9 wt.% NaCl equiv. (n = 46), suggesting that a relatively hot and saline fluid likely of magmatic origin was responsible for the first stage of mineralization. Some evidence for boiling was also observed in the veins. However, the bulk of the gold mineralization occurs as disseminations in the wall rocks, suggesting that wall-rock reactions were the main control on ore deposition.

We have evaluated results from a study combining airborne z-axis tipper electromagnetic (ZTEM) and ground-based magnetotelluric (MT) data to image an epithermal system in British Columbia. The spatially coincident use of these two methods allowed for a direct comparison of both data sets in the overlapping frequency band and showed that both measurements were consistent. Inversion of just the ZTEM data suffered from the lack of electric field amplitude information, which could be provided by the MT data. Three-dimensional inversion modeling of the two individual data sets was performed. Models of electrical resistivity derived from both data sets were consistent and could be correlated with the geological and structural setting of the mineralization. Gold is associated with disseminated pyrite and marcasite in quartz-sericite-altered felsic volcanic rocks and intrusions, especially near the contact with mafic volcanic rocks and a late diorite intrusion. The quartz-sericite alteration yields a conductivity anomaly, relative to the more resistive mafic country rocks. Although ZTEM and MT do not possess the resolution of the geologic model derived from borehole data, our model agrees well with a regional assessment of the deposit.

There is broad consensus that donors and policymakers have struggled to formalize artisanal and small-scale mining (ASM) – low-tech mineral processing and extraction – because of a poor understanding of the sector’s dynamics. As a result, ASM is largely – and at times, unfairly – associated with illegality, smuggling, and a loss of royalty fees and title fees. Experiences from four communities in north-central Nigeria reveal that, despite these views, ASM has evolved into a viable activity that can potentially support rural developmental objectives such as employment, reduction of poverty and hunger, as well as profits for government. Our findings suggest that there is a higher potential for success in the sector in Nigeria and Sub-Saharan Africa more generally if ASM is practiced in a community-inclusive way.

The 1600 km-long Gangdese magmatic belt features extensive Paleocene–Eocene I-type intrusive rocks and coeval volcanic successions, abundant but more localized Oligo-Miocene calc-alkaline to alkaline plutons, and Miocene potassic to ultrapotassic volcanic rocks. These Cenozoic igneous rocks record geodynamic changes related to the India–Asia collision which began at ~ 55–50 Ma. New and published lithogeochemical and multiple isotopic (Os–Sr–Nd–O–Hf) analyses of these Cenozoic igneous rocks reveal that the Paleocene–Eocene magmas have similar compositions to continental arc rocks throughout the belt, but later Miocene magmas show sharp longitudinal contrasts in geochemical and isotopic compositions, which are also correlated with the occurrence of porphyry-type mineralization. Sparse Miocene high-K calc-alkaline to shoshonitic volcanic rocks in the eastern Gangdese belt have low to moderate (87Sr/86Sr)i ratios (0.7057–0.7121), moderately negative eNdi values (− 9.4 to − 3.4), low (187Os/188Os)i ratios (0.154–0.210), highly variable eHfi values (− 5.9 to + 10.1), and low zircon δ18O values (+ 5.0–+ 6.7‰), which are interpreted to reflect derivation by partial melting of subduction-modified Tibetan sub-continental lithospheric mantle (SCLM). In contrast, Miocene high-K calc-alkaline to shoshonitic volcanic rocks in the western Gangdese belt have higher (87Sr/86Sr)i ratios (0.7069–0.7263), more negative eNdi values (− 17.5 to − 6.0) and eHfi values (− 15.2 to + 0.7), and crust-like zircon δ18O values (+ 6.2–+ 8.8‰), but mantle-like (187Os/188Os)i values (0.156–0.182), and high Ni and Cr contents. These features suggest that potassic to ultrapotassic magmas in the western Gangdese belt were also derived from partial melting of Tibetan SCLM but with ~ 3–25% input of melts ± fluids from the underthrust Indian plate (87Sr/86Sr = 0.74–0.76, eNd = − 18 to − 10, δ18O = + 10 − + 14‰). In contrast, Miocene alkaline magmas to the east were unaffected by this source. Oligo-Miocene calc-alkaline to high-K calc-alkaline granitoids related to large porphyry Cu–Mo deposits in the eastern Gangdese belt (east of ∼ 89° E) are geochemically broadly similar to the early Paleocene–Eocene rocks. They are thought to be derived from partial melting of subduction-modified lower crust with mixing of alkaline melts from partial melting of SCLM, and have relatively low (87Sr/86Sr)i ratios (0.7047–0.7076), high eNdi values (− 6.1 to + 5.5) and eHfi values (1.4–8.7), moderate (187Os/188Os)i ratios (0.224–0.835), and low zircon δ18OVSMOW values (+ 5.5–+ 6.6‰). These magmas also had high water contents (weak Dy/Yb enrichment, characterized with amphibole fractionation) and oxidation states (ΔFMQ 0.8–2.9), which explain their unique association with porphyry Cu–Mo mineralization. In contrast, Miocene high-K calc-alkaline to shoshonitic granitoids in the western Gangdese belt (west of ∼ 89° E) show differences in geochemical and isotopic compositions to the earlier Paleocene–Eocene magmatism, and are characterized by crust-like zircon δ18O values (+ 6.2–+ 8.8‰), high (87Sr/86Sr)i ratios (0.7147–0.7165), negative eNdi values (− 11.3 to − 7.9), crust-like (187Os/188Os)i values (0.550–1.035), and low eHfi values (− 13.0 to 3.9). These magmas are interpreted to reflect involvement of melts ± fluids from the underthrust Indian plate and high degrees of crustal contamination upon emplacement. Only one small porphyry Cu–Mo deposit is known to be associated with these western granitoids. We suggest that this difference reflects the variable extent of underthrusting of the Indian plate continental lithosphere beneath Tibet in the Oligo-Miocene, and diachronous breakoff of the Greater India slab. In the absence of underthrust Indian lithosphere to the east of ~ 89° E in the Oligo-Miocene, slab breakoff triggered asthenospheric upwelling and partial melting of previously subduction-modified Tibetan lithosphere, generating hydrous, oxidized calc-alkaline magmas with the potential to generate porphyry Cu–Mo deposits. In contrast, underthrusting of the Indian plate to the west at this time limited the involvement of asthenospheric melts and the extent of partial melting of subduction-modified lithosphere, with the result that melts ± fluids derived from the underthrust lithosphere were infertile.

This paper reviews the tectonic, magmatic, and metallogenic history of the Tethyan orogen from the Carpathians to Indochina. Focus is placed on the formation of porphyry Cu ± Mo ± Au deposits, as being the most characteristic mineral deposit type formed during both subduction and collisional processes in this region. Relatively little is known about the history of the Paleotethys ocean, which opened and closed between Gondwana and Eurasia in the Paleozoic, and few ore deposits are preserved from this period. The Neotethyan ocean opened in the Permian–Early Triassic as the Cimmerian continental fragments (the cores of Turkey, Iran, Tibet, and Indochina) rifted from the northern Gondwana margin and drifted northwards. These microcontinents docked with the Eurasian margin at various points in the Mesozoic and Cenozoic, and formed a complex archipelago involving several small back-arc basins and remnants of the Paleotethyan ocean. The main Neotethyan ocean and these smaller basins were largely eliminated by collision with India and Africa–Arabia in the early Eocene and early-mid Miocene, respectively, although Neotethyan subduction continues beneath the Hellenic arc and the Makran. The majority of porphyry-type deposits are found in association with Neotethyan subduction (mainly in the Mesozoic and Paleogene), and syn- to post-collisional events in the mid-Paleogene to Neogene. They are found throughout the orogen, but some sections are particularly well-endowed, including the Carpathians–Balkans–Rhodopes, eastern Turkey–Lesser Caucasus–NW Iran, SE Iran–SW Pakistan, southern Tibet, and SE Tibet–Indochina. Other sections that appear barren may reflect deeper levels of erosion, young sedimentary cover, or lack of exploration, although there may also be real reasons for low prospectivity in some areas, such as minimal subduction (e.g., the western Mediterranean region) or lithospheric underthrusting (as proposed in western Tibet). Over the last decade, improved geochronological constraints on the timing of ore formation and key tectonic events have revealed that many porphyry deposits that were previously assumed to be subduction-related are in fact broadly collision-related, some forming in back-arc settings in advance of collision, some during collision, and others during post-collisional processes such as orogenic collapse and/or delamination of subcontinental mantle lithosphere. While the formation of subduction-related porphyries is quite well understood, collisional metallogeny is more complex, and may involve a number of different processes or sources. These include melting of: orogenically thickened crust; previously subduction-modified lithosphere (including metasomatized mantle, underplated mafic rocks, or lower crustal arc plutons and cumulates); or upwelling asthenosphere (e.g., in response to delamination, slab breakoff, back-arc extension, or orogenic collapse). The most fertile sources for syn- and post-collisional porphyry deposits appear to be subduction-modified lithosphere, because these hydrated lithologies melt at relatively low temperatures during later tectonomagmatic events, and retain the oxidized and relatively metalliferous character of the original arc magmatism. Unusually metallically enriched lithospheric sources do not seem to be required, but the amount of residual sulfide phases in these rocks may control metal ratios (e.g., Cu:Au) in subsequent magmatic hydrothermal ore deposits. Relatively Au-rich deposits potentially form in these settings, as observed in the Carpathians (e.g., Rosia Montana), Turkey (Kisladag, Copler), and Iran (Sari Gunay, Dalli), although the majority of syn- and post-collisional porphyries are Cu–Mo-rich, and resemble normal subduction-related deposits (e.g., in the Gangdese belt of southern Tibet). This similarity extends to the associated igneous rocks, which, being derived from subduction-modified sources, largely retain the geochemical and isotopic character of those original arc magmas. While still retaining a broadly calc-alkaline character, these rocks may extend to mildly alkaline (shoshonitic) compositions, and may display adakite-like trace element signatures (high Sr/Y and La/Yb ratios) reflecting melting of deep crustal garnet amphibolitic sources. But they are otherwise hard to distinguish from normal subduction-related magmas. Small, post-collisional mafic, alkaline volcanic centers are common throughout the orogen, but for the most part appear to be barren. However, similar rocks in other post-subduction settings around the world are associated with important alkalic-type porphyry and epithermal Au ± Cu deposits, and the potential for discovery of such deposits in the Tethyan orogen should not be overlooked.

Global data for measured Fe2O3/FeO ratios and Cu contents in unaltered volcanic and intrusive arc rocks indicate that, on average, they are slightly more oxidized than other magmas derived from depleted upper mantle (such as MORB), but contain similar Cu contents across their compositional ranges. Although Cu scatters to elevated values in some intermediate composition samples, the bulk of the data show a steady but gentle trend to lower concentrations with differentiation, reaching modal values of ~ 50–100 ppm in andesitic rocks. These data suggest that Cu is mildly compatible during partial melting and fractionation processes, likely reflecting minor degrees of sulfide saturation throughout the magmatic cycle. However, the volume of sulfides must be small such that significant proportions of the metal content remain in the magma during fractionation to intermediate compositions. Previous studies have shown that andesitic magmas containing ~ 50 ppm Cu can readily form large porphyry-type Cu deposits upon emplacement in the upper crust. A review of the literature suggests that the elevated oxidation state in the asthenospheric mantle wedge source of arc magmas (ΔFMQ ≈ + 1 ± 1) derives from the subduction of seawater-altered and oxidized oceanic crust, and is transmitted into the mantle wedge via prograde metamorphic dehydration fluids carrying sulfate and other oxidizing components. Progressive hydration and oxidation of the mantle wedge may take up to ~ 10 m.y. to reach a steady state from the onset of subduction, explaining the rarity of porphyry deposits in primitive island arcs, and the late formation of porphyries in continental arc magmatic cycles. Magmas generated from this metasomatized and moderately oxidized mantle source will be hydrous basalts containing high concentrations of sulfur, mainly dissolved as sulfate or sulfite. Some condensed sulfides (melt or minerals) may be present due to the high overall fS2, despite the moderately high oxidation state. These sulfides may retain some highly siderophile elements in the source, but are unlikely to be sufficiently voluminous to significantly affect the budget of more modestly sulfide-compatible and more abundant elements such as Cu and Mo. These primary magmas can therefore be considered to be largely Cu-Mo-undepleted, although highly siderophile elements such as Au and platinum group elements (PGE) may be depleted unless no sulfides remain in the source. The latter condition seems unlikely during active subduction because of the continuous flux of fresh sulfur from the slab, but may occur during post-subduction re-melting (leading to potentially Au-rich post-subduction porphyry and alkalic-type epithermal systems). Lower crustal differentiation of main-stage arc magmas results in some loss of Cu to residual or cumulate sulfides, but again the amount appears to be minor, and does not drastically reduce the Cu content of derivative intermediate-composition melts. Fractionation and devolatilization affect the oxidation state of the magma in competing ways, but, while crystallization and segregation of Fe3 +-rich magnetite can cause reduction in reduced to moderately oxidized evolved magmas, this effect appears to be outweighed by the oxidative effects of degassing reduced or weakly oxidized gaseous species such as H2, H2S, and SIVO2, and preferential solvation and removal of Fe2 + in saline hydrothermal fluids. Consequently, most arc magmatic suites show slight increases in oxidation state during differentiation, reaching typical values of ΔFMQ = + 1 to + 2. This oxidation state is significant, because it corresponds to the transition from dissolved sulfide to sulfate dominance in magmas. It has been shown that Cu and Au solubilities in silicate magma increase up to this level (ΔFMQ ≈ + 1), but while Cu solubility continues to increase at higher oxidation states, Au shows a precipitous drop as sulfide, which solvates Au in the melt, is converted to sulfate. This may explain the somewhat restricted distribution of Au-rich porphyry Cu deposits, but the general association of porphyry Cu deposits with relatively oxidized magmas. Exsolution of a saline, high temperature aqueous fluid enables metals to partition from the magma into a highly mobile volatile phase. Sulfur also partitions strongly into this fluid phase, predominantly as SO2 at ΔFMQ = + 1 to + 2. However, as the fluid cools below ~ 400 °C, SIVO2 disproportionates to form reduced H2S− II and oxidized H2SVIO4. The H2S bonds with metals in solution to precipitate as Cu- and Mo-sulfides, while the H2SO4 (and HCl) generates progressively acidic wallrock alteration (phyllic, argillic, advanced argillic). Gold may precipitate with early Cu/Mo-sulfides, but some may also stay in solution as bisulfide complexes, eventually reaching the epithermal environment. Thus, three components, [S], [H2O], and fO2 work together throughout subduction and arc magmatic processes to transport chalcophile and siderophile metals from the mantle into the upper crust, where they may be concentrated by hydrothermal processes to form ore deposits. These processes are far from 100% efficient, and metals (especially highly siderophile elements such as Au and PGE) may be left behind at various stages of the passage of arc magmas through the lithosphere, where they may form potentially metalliferous source rocks for partial melts and subsequent magmatic-hydrothermal ore deposits generated during later tectonomagmatic events.

Northward Neo-Tethyan oceanic lithosphere subduction beneath southern Tibet in the Mesozoic–Early Cenozoic produced continental arc magmas in the ~ 1600 km-long Gangdese belt. The most voluminous magmatism occurred in the Paleocene–Eocene, and is characterized by extensive I-type calc-alkaline to high-K calc-alkaline Linzizong volcanic rocks, and coeval plutons. These rocks have been extensively studied in the eastern Gangdese belt (east of ∼ 89°E), but few data exist from the western Gangdese belt. New data for eleven samples of these rocks, combined with existing data from the literature, show that they are similar to the eastern Gangdese belt rocks, with relative depletions in Nb, Ta, P, and Ti, and enrichments in Rb, Ba, Th, U, K, Pb, Zr, and Hf on a primitive mantle-normalized trace element diagrams, typical of continental arc-related igneous rocks. However, compared to the east, western Gangdese igneous rocks range to higher K2O contents (up to 6.1 wt.%), higher (87Sr/86Sr)i ratios (up to 0.7151), and lower eNdi values (down to − 8.1), suggesting that an evolved crustal source was involved in arc magmatism. The Gangdese arc magmatism lasted to ~ 80 Ma, and has a gap or quiescent period afterwards. Starting at ~ 69 Ma, the arc magmatism initiated and shifted southward from ~ 30.5°N to ~ 29.5°N in southern Tibet with an abrupt change of India-Asia convergence rate. The magmas through the whole Gangdese belt at ~ 69–53 Ma are characterized by intermediate eNdi values (− 0.6 to + 4.0), positive eHfi values (+ 3.8 to + 7.1), intermediate δ18O values (+ 5.0‰ to + 6.5‰), and low Th/Y and La/Yb ratios (

The Yangla skarn Cu deposit (150 Mt at 1.03% Cu) is located in the central segment of the Jinshajiang metallogenic belt within the Sanjiang (Three Rivers) region, southwest China. Skarn orebodies are mainly developed between different units of Devonian carbonate and quartz sandstone rocks with stratiform-like shapes, or within the contact zone between granitoids and marbles. Re-Os dating of molybdenite intergrowth with chalcopyrite yielded a well-constrained 187Re-187Os isochron age of 232.0 ± 1.5 Ma with a weighted average age of 231.8 ± 1.3 Ma, both coeval with the related intrusions (233.1 ± 1.4 and 231.0 ± 1.6 Ma at 2 σ by zircon U-Pb dating) from our previously published work. Field and textural relationships indicate three hydrothermal stages characterized by assemblages of prograde skarn (pre-ore stage), retrograde skarn and Cu-Fe-Mo-Bi sulfides (main ore stage), and Pb-Zn sulfides associated with calcite and quartz (late ore stage), as well as one supergene stage marked by secondary Cu mineralization (malachite and azurite). Skarns contain garnets with andraditic compositions (Ad96 Gr2~3Py0~1) and clinopyroxene (two series: Hd6Di94 and Hd86Di13Jo1) with low Mn/Fe ratios (

Stepanov et al. (Contrib Mineral Petrol, 2017) question our conclusion that the UPVs in southern Tibet were derived by partial melting of an old, metasomatized subcontinental lithospheric mantle (SCLM) of the subducted Indian plate. Instead, they propose that these ultrapotassic volcanic rocks (UPVs) are shoshonitic and were generated in two steps: direct melting of crustal rocks first, and then the melts interacted with mantle peridotite. However, the trace element, isotopic, thermal, structural, and seismic evidence is consistent with the xenolith evidence (Wang et al in Contrib Mineral Petrol 172:62, 2016) for hybridisation of ascending Indian subcontinental lithospheric mantle-derived UPV magmas with the deep, isotopically unevolved, Tibetan crust. This necessitates a model whereby partial melting of subducting Indian SCLM generates the UPV suite of southern Tibet.

In the Gangdese belt of southern Tibet, Paleocene-Eocene magmas record the final stage of Neo-Tethyan subduction and are associated with few economic porphyry deposits. In contrast, magmas formed during later stages of the India-Asia collision in the Oligo-Miocene are associated with several large porphyry Cu-Mo ± Au deposits, especially in the eastern part of the belt (~89°E–93°E). In a previous study, we showed that these Oligocene-Miocene magmas were more hydrous than the earlier Paleocene-Eocene magmas. In this study, we show that the later magmas were also more oxidized. Paleocene-Eocene rocks from the eastern Gangdese belt are characterized by low zircon Ce4+/Ce3+ ratios (mostly 50; 32.3–141.9, average = 74.3 ± 30.1, n = 33) and higher f O2 values (ΔFMQ 0.8–2.9, average = 1.8 ± 0.8, n = 6; estimated from magnetite-ilmenite pairs). Estimates of magmatic oxidation state from amphibole compositions also show an increase from ΔFMQ 1.2 to 2.1 (average = 1.6 ± 0.2, n = 40) in the Paleocene-Eocene to 2.0 to 2.8 (average = 2.5 ± 0.2, n = 58) in the Miocene. Sparse whole-rock Fe3+/Fe2+ ratios show the same general trend. The amphibole results are systematically shifted to higher ΔFMQ compared to data from magnetite-ilmenite pairs, but their trend is internally consistent. The higher oxidation states (and water contents) of Miocene igneous rocks from the eastern Gangdese belt may explain their unique association with large porphyry deposits in Tibet, because it has been shown from other studies that the potential for hydrous calc-alkaline magmas to transport Au is maximized near ΔFMQ 1.0, and Cu at higher f O2 (ΔFMQ > 1.0). In comparison, the somewhat less oxidized and less hydrous Paleocene- Eocene magmas would have been less fertile for the formation of such deposits. We suggest that these differences reflect the tectonomagmatic evolution of the Cenozoic Gangdese collisional orogen, from early collision-related magmatism in the Paleocene-Eocene to late collisional magmatism in the Miocene. Asthenospheric upwelling following slab breakoff in the Miocene caused interaction between mantle-derived magmas and previously subduction modified Tibetan lithosphere. The resultant evolved magmas had high oxidation states and water contents, which are favorable properties for the subsequent formation of magmatic-hydrothermal porphyry Cu-Mo ± Au deposits.

The Cevizlideredeposit, located in the Tunceli Prov- ince of eastern Anatolia, is the largest porphyry Cu-Au-Mo system in Turkey. The deposit is spatially related to a composite stock, which was emplaced into Paleozoic lime- stones and Paleogene andesitic rocks to the southeast of the Munzur mountains, near the southwestern margin of the Ovacik pull-apart basin. The host plutonic rocks at Cevizlidere are porphyritic, medium-K calc-alkaline diorites andgranodiorites. 40 Ar/ 39 Ar incremental step-heatinganalysis of two igneous biotite separates obtained from syn-mineral diorite porphyry yielded late Oligocene cooling ages of 25.49 ±0.10 and 25.10±0.14 Ma, whereas hydrothermal biotite yielded an age of 24.73±0.08 Ma. Re-Os ages obtained from two molybdenite separates (24.90±0.10 and 24.78±0.10 Ma) indicate that porphyry-style alteration and mineralization de- veloped shortly after magma emplacement. The whole-rock geochemical composition of the Cevizlidere porphyry intru- sions is consistent with derivation from partial melting of the metasomatized supra-subduction zone mantle. However, based on regional tectonic reconstructions, Oligocene mag- matic activity in this area appears to be related to a major kinematic reorganization that took place at around 25 Ma, during the switch from subduction to collisional tectonics in eastern Anatolia. This kinematic switch may be attributed to break-off of the Southern Neo-Tethys oceanic slab prior to the Arabia-Eurasia continent-continent collision (~12-10 Ma) following widespread middle Eocene (50-43 Ma) arc/back- arcmagmatism.Inthisrespect,thesubduction-relatedtectonic setting of the late Oligocene Cevizlidere porphyry deposit is similar to that of the middle Eocene Copler epithermal Au deposit. The late timing of Cevizlidere with respect to the Southern Neo-Tethys subduction may be comparable to some early to late Miocene porphyry-epithermal systems that lie within the contiguous Urumieh-Dokhtar belt in central Iran. The later timing in Iran reflects the diachronous nature of the Eurasia-Afro-Arabia collision.

Mesozoic intermediate–felsic intrusions are widely distributed in the southern Qiangtang terrane, central Tibet. Zircon U–Pb dating shows that these intrusions formed in two periods in the Jurassic (169–150 Ma) and Cretaceous (127–113 Ma). They mostly belong to the high-K calc-alkaline series, and show strong enrichments in large ion lithophile elements (e.g., Cs, Rb, and K), depletions in Nb, Ta, and Ti, and negative Ba anomalies on primitive mantle-normalized diagrams. P2O5 contents decrease with increasing SiO2 content, and Th contents increase with increasing Rb content, consistent with the evolution trend of I-type magmas. These intrusions show a wide range of zircon eHf(t) values from − 19.4 to 11.2, including in rocks with similar SiO2 contents, suggesting variable contributions from mantle and Qiangtang crustal sources. Fine-grained mafic to intermediate igneous enclaves in Jurassic intrusions have similar zircon U–Pb ages and similar or slightly higher zircon eHf(t) values to the host rocks, suggesting that the enclave magmas were derived from mixed magmas at depth and injected into more evolved magmas in upper crustal magma chambers. Magma mixing is also supported by the wide range of zircon Hf isotopic compositions (eHf(t) = − 19.4 to 2.5) from within individual Jurassic and Cretaceous intrusions, and Jurassic enclaves. The Jurassic–Cretaceous magmas likely formed in a continental arc setting during subduction of the Bangong–Nujiang ocean between 170 and 110 Ma, and evolved in the upper plate crust by MASH processes.

The eastern Gangdese magmatic belt (between ~88°E and 94°E) in southern Tibet preserves a series of Cenozoic collision-related igneous rocks that have formed since the start of the India-Asia collision at ~55 to 50 Ma. The Paleocene-Eocene intrusive rocks have intermediate [La/Yb]N ratios (average = 12.0 ± 8.6, n = 150, range = 1.15–64.0), intermediate to low Sr/Y ratios (mostly 40), and weak or absent Eu anomalies (Oligocene: average Eun/Eu* = 0.99 ± 0.27, n = 35, range = 0.67–1.58; Miocene: average Eun/Eu* = 0.95 ± 0.13, n = 93, range = 0.44–1.46), and their mineralogy consists mainly of plagioclase, quartz, and amphibole as phenocrysts. These geochemical and mineralogical characteristics suggest that the Paleocene-Eocene magmas were relatively dry and evolved primarily by fractionation of pyroxene and plagioclase, whereas the Oligo-Miocene magmas were more hydrous and fractionated significant amounts of hornblende and lesser plagioclase prior to upper crustal emplacement. Oligocene and Miocene magmatic rocks crop out as small-volume intrusions and are associated with several large- to giant-sized porphyry Cu-Mo ± Au deposits, as well as numerous smaller porphyry and skarn deposits. In contrast, earlier, more voluminous Paleocene-Eocene magmatism is only known to be associated with three small porphyry deposits. Erosional loss of subvolcanic porphyry systems from the older Paleocene-Eocene sequence is not thought to be a controlling factor in this temporal difference of porphyry deposit distribution, because exhumation rates after 55 Ma were relatively low and coeval Linzizong volcanic rocks are extensively preserved in the Gangdese belt. Instead, we suggest that differences in magmatic history and petrogenesis led to this restricted temporal distribution of porphyry deposits. The Paleocene-Eocene magmas were generated during the onset of collision between India and Asia and were triggered by rollback of the Neo-Tethyan oceanic slab. This magmatism records the final stage of subduction of oceanic lithosphere beneath southern Tibet, and the relatively low water contents of these magmas may reflect final dehydration of the remnant Neo-Tethyan slab. After oceanic slab breakoff at ~40 to 38 Ma, the Greater India slab began to subduct beneath the Gangdese belt until hard collision with the Indian continent at ~35 Ma and breakoff of the Greater India slab at ~25 to 10 Ma. This last event is suggested to have caused asthenospheric mantle upwelling and partial melting of subduction-modified, hydrated Tibetan lithosphere. The resultant partial melts may have remobilized metals residual in the deep lithosphere from previous arc magmatism, giving rise to a suite of postsubduction porphyry-type ore deposits in the Oligo-Miocene.

The plate tectonic revolution, culminating with the formulation of the Wilson Cycle, took place over a period of less than a decade in the 1960s and early 1970s. The model provided a framework for understanding the formation of almost every type of mineral deposit then known on Earth, ranging from base and precious metal deposits associated with rifting, to porphyry Cu–Mo and epithermal Au deposits associated with subduction, and collision-related mesothermal Au deposits. By the end of the 1970s, satisfactory tectonic models for most of these deposit types had been established. Modern geological and economic geology research is largely built on these models, which have been expanded in detail but remain largely unchanged in concept and function.

Miocene igneous rocks in the 1,600 km-long E–W Gangdese belt of southern Tibet form two groups separated at longitude ∼89° E. The eastern group is characterized by mainly intermediate–felsic calc-alkaline plutons with relatively high Sr/Y ratios (23 to 342), low (87Sr/86Sr)i ratios (0.705 to 0.708), and high eNdi values (+5.5 to −6.1). In contrast, the western group is characterized by mainly potassic to ultrapotassic volcanic rocks with relatively high Th and K2O contents, low Sr/Y ratios (11 to 163), high (87Sr/86Sr)i ratios (0.707 to 0.740), and low eNdi values (−4.1 to −17.5). The eastern plutonic group is associated with several large porphyry Cu–Mo ± Au deposits, whereas the western group is largely barren. We propose that the sharp longitudinal distinction between magmatism and metallogenic potential in the Miocene Gangdese belt reflects the breakoff of the Greater India slab and the extent of underthrusting by the Indian continental lithosphere at that time. Magmas to the east of ∼89° E were derived by partial melting of subduction-modified Tibetan lithosphere (mostly lower crust) triggered by heating of hot asthenospheric melt following slab breakoff. These magmas remobilized metals and volatile residual in the crustal roots from prior arc magmatism and generated porphyry Cu–Mo ± Au deposits upon emplacement in the upper crust. In contrast, magmas to the west of ∼89° E were formed by smaller volume partial melting of Tibetan lithospheric mantle metasomatized by fluids and melts released from the underthrust Indian plate. They are less hydrous and oxidized and did not have the capacity to transport significant amounts of metals into the upper crust.

Giant ore deposits are priority targets for mining companies. A review of the characteristics of several giant porphyry and epithermal deposits worldwide suggests that they formed from ordinary processes that were fortuitously operating at maximum efficiency.

Buriticá is an intermediate sulfidation epithermal gold deposit located 75 km north of Medellín, Colombia. It is hosted by the late Miocene Buriticá andesite porphyry, a shallow-level pluton dated at 7.41 ± 0.40 Ma (2σ, MSWD = 2.30; 40Ar/39Ar on hornblende), which intrudes volcanic and sedimentary rocks of the Cretaceous Barroso Formation and the Buriticá stock. Gold mineralization is associated with proximal sericite-adularia and distal epidote-dominant propylitic alteration, and is hosted in two different sets of veins striking and dipping ~072°/87°S and ~105°/87°S. Hydrothermal muscovite from altered host rock yielded a weighted average 40Ar/39Ar age of 7.74 ± 0.08 Ma for two samples, which is within error of the age of the Buriticá andesite porphyry. Sulfur isotope, oxygen isotope, and fluid inclusion data indicate that a relatively hot and saline fluid (~310°C, up to ~8 wt % NaCl equiv) of probable magmatic origin was dominant during stage 1 of the vein paragenesis, and mixed with cooler and less saline groundwater during stages 2 and 3. Fluid inclusion assemblages indicate that boiling was likely the main control on gold precipitation, with cooling and dilution during fluid mixing as a possible secondary control.

Iron oxide-copper-gold (IOCG) deposits formed by magmatic-hydrothermal fluids (MH-IOCG) share many similarities with, but have important differences from, porphyry Cu ± Mo ± Au (porphyry) deposits: MH-IOCG deposits predominantly occur in Precambrian rocks, are Fe oxide rich, and have volumetrically extensive high-temperature alteration zones, whereas porphyry deposits occur almost exclusively in Phanerozoic rocks, are Fe sulfide rich, and have narrower high-temperature alteration zones. We propose that these deposit types are linked by common subduction-modified magmatic sources, but that secular changes in oceanic sulfate content and geothermal gradients at the end of the Precambrian caused a transition from the predominance of S-poor arc magmas and associated S-poor MH-IOCG systems, to S-rich arc magmas and associated S-rich porphyry deposits in the Phanerozoic. Phanerozoic MH-IOCG and rare Precambrian porphyry deposits are explained by local or periodic fluctuations in oceanic oxidation state and sulfate content, or remobilization of previously subduction-modified lithosphere in post-subduction tectonic settings.