EARTH SCIENCE FRONTIERS
Volume 14, Issue 5, September 2007
Online English edition of the Chinese language journal

Cite this article as: Earth Science Frontiers, 2007, 14(5):124–132.                                                          RESEARCH PAPER




Applications of Stable and Radiogenic Isotopes to
Magmatic Cu-Ni-PGE Deposits: Examples and Cautions
Edward M. RIPLEY1,2,∗, Chusi LI1,2
1 Department of Geological Sciences, Indiana University, Bloomington, IN 47401, USA
2 State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Beijing 100083, China




      Abstract: Measurements of S, O, and radiogenic isotope ratios are all potentially powerful tracers of magma interaction with
      country rocks and the importance of assimilation processes in the genesis of magmatic Ni-Cu-PGE deposits. Sulfur isotope
      measurements of deposits such as those in the 1.1 Ga Duluth Complex, the Permo-Triassic intrusions of the Noril’sk area, and the
      1.4 Ga Kabanga intrusions provide evidence for the derivation of S from both sulfide-and sulfate-bearing country rocks. The 1.3 Ga
      Voisey’s Bay deposit provides an example where δ34S values of the ores commonly fall within the accepted mantle range of 0 ± 2‰,
      but detailed studies of Proterozoic metasedimentary country rocks show that their weighted average δ34S value is also within this
      range. A thorough knowledge of the isotopic compositions of potential contaminants is essential for a proper evaluation of the role
      of country rock derived S in ore formation. When O and radiogenic isotopic measurements are employed as tracers to evaluate
      magma-country rock interaction, it is essential to provide evidence that open system processes have not perturbed the isotopic
      systematics. Low-temperature hydrothermal processes can mask evidence of high-temperature processes in the oxygen isotope
      system and in radiogenic systems such as Re/Os and Pb where involved elements may be mobile under hydrothermal conditions, or
      host phases may close to exchange and uptake at different temperatures. Careful petrographic observation and analyses of
      individual minerals may be required before the models involving the contamination of magmas by country rocks can be
      meaningfully applied.

      Key Words: S isotopes; O isotopes; radiogenic isotopes; magma-country rock interaction



                                                                                  importance of S assimilation in promoting sulfide saturation in
1    Introduction
                                                                                  a mafic magma[6,7].
                                                                                     This review is intended to highlight some of the
   Isotopic analyses are a powerful tool for evaluating the
                                                                                  fundamental applications of isotopic measurements to the
potential interaction between mantle-derived magmas and
                                                                                  assessment of magma contamination as applied to the studies
crustal rocks. Both stable and radiogenic isotopic systems
                                                                                  of magmatic ore genesis and cautionary notes that need to be
have been applied to the studies of magma contamination and
origin. It has become apparent that the contamination of mafic                    considered to avoid potential pitfalls in the interpretation of
magma by crust plays a significant role in the generation of                      isotopic data. Examples from well-known Cu-Ni-PGE
large, sulfide-rich magmatic Cu-Ni-PGE deposits[1,2].                             deposits will be utilized to illustrate key concepts.
Interaction between siliceous country rocks and mafic magma
has also been suggested as an important prerequisite for the                      2     Sulfur isotopes
development of chromite deposits in layered intrusions[3,4] and
sulfide deposits in some intrusions[5]. Experimental studies                        Owing to the relatively large difference in δ34S values
that emphasize the role of pressure as a control of the                           between mantle-derived magma (MORB as the principal
solubility of sulfur in mafic magmas have also stressed the                       example) and either sulfide or sulfate minerals in sedimentary

Received date: 2007-08-20; Accepted date: 2007-09-11.
* Corresponding author. E-mail: ripley@indiana.edu
Foundation item: Supported by the Ministry of Education of China project (111-B07011); the National Natural Science Foundation of China (40534020); the
National Science Foundation of the United States (EAR-0710910).
Copyright © 2007, China University of Geosciences (Beijing) and Peking University, Published by Elsevier B.V. All rights reserved.
Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132

rocks, sulfur isotopic studies have provided important insights          Tanzania[17]. Mica schists and graphitic schists contain
into the processes of magma contamination. Mantle sulfur is              pyrrhotite with δ34S values ranging from -2‰ to 24‰. Sulfide
commonly taken as having δ34S values between -2‰ and 2‰.                 mineralization in mafic to ultramafic intrusive rocks is
This range is defined largely from the analyses of sulfide               characterized by a range of δ34S values between 10‰ and
minerals in meteorites and sulfide-bearing basalts (see                  24‰, again lending strong evidence to a genetic model that
Ripley[8] for a review of early S isotope studies). It should            involves magma contamination by country rocks.
always be kept in mind, however, that sulfide minerals in
diamonds from kimberlites are characterized by a wide range
of δ34S values from -11‰ to +14‰[9,10], suggesting that the
deep mantle may be far more heterogeneous in terms of δ34S
values than the shallow mantle appears to be. Seawater δ34S
has varied between ~10‰ and 30‰ from the Neoproterozoic
to the present[11,12]. Sulfide minerals are produced from
seawater sulfate owing to reduction that is often related to the
life cycles of sulfate-reducing bacteria. Depending on the
variables related primarily to the rate of sulfate reduction, the
sulfide produced may be characterized by positive or negative
δ34S values. Pyrite that forms owing to the reaction of
produced sulfide and reactive iron may show a wide range in
δ34S values, typically from values close to the source sulfate
(e.g. ~ 20‰), to values as low as -50‰, or occasionally even
less. Sulfate minerals that form via evaporative processes
show little isotopic fractionation relative to the source sulfate,
and hence are generally characterized by δ34S values between
~10‰ and 30‰.
   It is clear that because of the potentially large S isotopic
difference between sedimentary S-bearing minerals and
mantle-derived magmas, assimilation processes may be
elucidated using S isotope measurements. Three examples will
serve to illustrate the basic principles.

2.1   Duluth Complex, Minnesota

   Cu-Ni sulfide mineralization in the 1.1 Ga Midcontinent
Rift-related Duluth Complex represents one of the best-known
examples of major sulfur derivation from country rocks. The
Proterozoic (~1.85 Ga) Virginia Formation is a pelitic
sequence that contains layers that are sulfidic and
carbonaceous, and are characterized by pyrite δ34S values
ranging from 0 to 30‰ (Fig. 1)[13–15]. Metamorphism of the
sulfide-bearing sedimentary rocks has led to the production of
pyrrhotite and the generation of an H2S-bearing, carbonic fluid.
Magma that formed the Duluth Complex interacted with
                                                                                                        –15]
various levels of the sedimentary sequence, assimilated the               Fig. 1   Sulfur isotope data[13      from sulfide mineralization in the
sedimentary-derived H2S, and crystallized the sulfide minerals                  Duluth Complex and Virginia Formation country rocks
that were characterized by elevated δ34S values ranging from
~0 to 20‰ (Fig. 1). Various mixing models involving external
                                                                         2.2    Noril’sk, Siberia
S and mantle S can be derived[16], but the importance of
externally derived S is highlighted by the δ34S values and the              Another well-known example where δ34S values of sulfide
internally consistent models of fluid production in the                  ores strongly suggest that assimilation has been operative is
metamorphic rocks and localized partial melting.                         the world famous mineralization at Noril’sk, in the
   A similar geologic setting in terms of mafic intrusions               Permo-Triassic Siberian flood basalt province. Early sulfur
found hosted by pelitic and graphitic metasedimentary rocks              isotope studies[18–20] suggested that assimilation of evaporites
has been described for the large 1.4 Ga Kabanga Ni-deposit in            that are common in the stratigraphic sequence occurred, or
Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132

that sulfur may have been derived from pyritic shales and                envisioned, or hydrothermal leaching of sulfate followed by
coals, or H2S derived from the coals. Li and his                         partial reduction, or a process involving the dissolution, rather
co-workers[21,22] developed slightly different models to explain         than melting, of sulfate in mafic magma. The observation of
the anomalous S isotopic signatures (Fig. 2); however, the               abundant wollastonite in contact aureole rocks suggests that
probable involvement of sulfate as a contaminant is                      reactions similar to CaSO4 + SiO2 + H2O = CaSiO3 + H2S +
noteworthy. The assimilation of sulfate via partial melting by a         2O2 occurred, and that sulfate may have been reduced to
basatic magma is difficult to conceive owing to the high                 sulfide prior to incorporation into the magma. The important
melting points of sulfate minerals (typically ranging between            point here is that sulfate minerals may serve as S sources, but
1360 and 1450°C). Therefore, either a fluxing mechanism is               the mechanism of S assimilation must be carefully evaluated.




                                        Fig. 2   Sulfur isotope data[21,22] from the Noril’sk area

                                                                         process, and positive values indicating a relatively rapid rate
2.3   Voisey’s Bay
                                                                         of sulfate reduction, possibly under closed system conditions.
   A third example where the interpretation of δ34S values is            The weighted average δ34S value of the sulfide assemblages in
critical is the 1.3 Ga Voisey’s Bay Ni-Co-Cu deposit in                  the Tasiuyak Gneiss is -1.3‰. This value is very similar to the
Labrador, Canada. The deposit is associated with the Nain                average value of the sulfide ores and illustrates a very
Plutonic Suite, which was emplaced along a boundary                      important principle. It is important to precisely know the δ34S
between crustal blocks, one comprised predominantly of                   values of the proposed contaminant before concluding. In the
Archean rocks and the other dominated by Proterozoic rocks.              case of Voisey’s Bay, even though the δ34S values are close to
Sulfide mineralization is found in troctolitic to gabbroic rocks.        0‰, the additional study clearly showed that S could have
Country rocks include the Proterozoic Tasiuyak Gneiss, which             been derived from the country rocks, in agreement with the
is locally sulfidic and graphitic. Several investigators[23,24]          widespread geological evidence, which suggests that S
proposed that large amounts of S in the mineralization may               assimilation occurred.
have been derived from the Tasiuyak Gneiss. Initial S isotopic              Voisey’s Bay is also an excellent example where supporting
studies[25] determined that the bulk of the sulfide                      radiogenic isotopic data was extremely valuable. The Re/Os
mineralization was characterized by δ34S values ranging                  isotopic data (see below)[27] clearly indicated that the magma
between 2‰ and -4‰ (Fig. 3). The δ34S values in the sulfide              from which the sulfides were derived assimilated crustal
ores were significantly depleted only locally near the contact           material. Carbon isotopic data also indicated that assimilation
with the Tasiuyak Gneiss, which contained pyrrhotite with                of sedimentary country rocks had occurred[26], as did the S/Se
negative δ34S values. However, additional studies[26]                    ratios of the sulfide mineralization. The supporting data
determined that pyrrhotite in the Tasiuyak Gneiss was                    strengthens the promise that even though the δ34S values of
characterized by a wide range in δ34S values, including both             the ores are close to those of mantle sulfides, a large
negative values suggesting an open system sulfate reduction              component of the S was derived from country rocks.
Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132

                                                                         derived in part from underlying sedimentary rocks. The
                                                                         anomalous isotopic signature illustrated in the δ34S-δ33S plot
                                                                         strongly suggests that the interpretation was correct.




                                                                                Fig. 4 δ33S versus δ34S for sulfide samples from the
                                                                                    Stillwater Complex, Voisey’s Bay, and Alexo
                                                                             IAEA S-1 is an international standard. Note the anomalous
                            [26]
Fig. 3 Sulfur isotope data from sulfide mineralization at                    character of the Alexo samples, strongly supporting the premise
 Voisey’s Bay and sulfide minerals in the Tasiuyak Gneiss                    that sulfur from Archean-aged metasedimentary country rocks has
                      country rocks                                          been involved in ore genesis.

   A final point with respect to the use of S isotopes as                3     Oxygen isotopes
indicators of crustal contamination regards Archean-aged
sulfide deposits. Ripley[8] discussed that during the Archean               Mantle-derived mafic magmas are generally characterized
bacterial sulfate reduction, then if operative, appears to have          by δ18O values in the range of 5‰ to 7‰[31–33]. Sedimentary
generated considerably smaller ranges in δ34S values of                  rocks and crustally derived granitic rocks are normally
reduced sulfur species, and values close to 0‰. For this                 characterized by considerably higher δ18O values, typically in
reason, the determination of externally derived S in Archean             the range of 8‰ to 30‰. Owing to these large differences,
deposits, or in younger deposits that may have interacted with           oxygen isotopes may also be sensitive indications of country
Archean rocks, may be difficult using δ34S values. However, a            rock assimilation. General features of how oxygen isotopes
process known as mass-independent fractionation (MIF) of S               may be applied in the study of magmatic ore deposits have
isotopes has led to variations between δ33S and δ34S, which are          been discussed by Ripley[8] and will not be repeated here.
distinct from most terrestrial processes that fractionate S              However, common difficulties with a straight-forward
isotopes and are referred to as mass-dependent fractionation             application of δ18O values to problems of assimilation by
processes. Farguhar and Wing[28] discussed the observation               mafic magmas must be highlighted.
that sedimentary pyrite of Archean age is characterized by                  It is often assumed that elevated δ18O values in mafic
δ33S-δ34S relations that are indicative of MIF. The process has          igneous rocks (e.g. >7‰) are indicative of contamination and
been attributed to low oxygen levels in the Archean                      assimilation of high-18O country rocks. Although this may
atmosphere and photochemical reactions involving SO2, but                indeed be the case, careful petrographic examination of the
the controversy with respect to the origin of the process                rocks, along with an evaluation of the magnitude of the δ18O
continues. The record of MIF may be present in magmas that               shift, should be undertaken. For example, consider a mafic
have assimilated S of the Archean age[29,30], and future studies         rock with bulk δ18O value of 9‰, a contaminant with a δ18O
using the technique may be able to distinguish the importance            value of 16‰, and a presumed magma initial value of 5.5‰.
of S derived from Archean rocks in ore genesis. Figure 4 is a            Assuming that the concentrations of O in the contaminant and
δ33S-δ34S plot showing sulfide samples from Voisey’s Bag (1.3            the magma are similar, a value of 9‰ in the mafic rock will be
Ga) and Alexo (a komatiite-associated, ~2.7 Ga deposit).                 suggested in ~35‰ bulk contamination. This degree of
Naldrett[1] proposed that S in the Alexo deposit may have been           contamination is unreasonable in most situations because of
Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132

the amount of heat required to achieve such a large                     geosciences. Several excellent text books and review articles
contaminant/magma ratio[34]. Unless xenoliths provide                   address these methodologies[41–44]. Here we wish to emphasize
evidence for such high degrees of contamination (the presence           that care must be taken to assure that closed-system conditions
of refractory oxides such as corundum and spinel), other                prevailed before conclusions regarding contamination history
methods that can explain the 18O-elevation must be explored.            can be made. For example, the Re/Os isotopic data from
Obviously knowledge of the δ18O values of potential                     sulfide minerals in the Voisey’s Bay deposit[27] mentioned
contaminants is extremely important. The assimilation of very           above produces an isochron, which is consistent with the age
high-18O siliceous rocks, carbonates, or hydrothermally altered         determined via U-Pb geochronology Fig. 5. The elevated
felsic rocks can be responsible for elevated δ18O values in             initial 187Os/188Os ratio clearly indicates that crustal
mafic igneous rocks and can satisfy reasonable thermal                  contamination occurred, the sulfide system was well-mixed,
constraints.                                                            and that closed system conditions have generally prevailed. In
   The most common cause of elevated δ18O values in mafic               contrast, in situations where a reasonable isochron is not
rocks that is often overlooked is low-temperature                       produced, an interpretation involving crustal contamination
hydrothermal alteration and isotopic exchange. Low-                     must proceed with caution. Figure 6 illustrates a case where an
temperature isotopic exchange involving a fluid has been                initially chondritic system interacts with a fluid after
discussed by several investigators[35–37]; one pertinent example        crystallization. Loss of rhenium during the interaction may
will be presented here. At temperatures below ~200°C, the               significantly lower the 187Re/188Os ratio. Continued decay may
oxygen isotopic fractionation between plagioclase and water is          then produce samples with 187Os/188Os ratios that will generate
greater than 6‰[38]. For a situation where the fluid-rock               anomalous γOs values, assuming that an age for the rocks is
interaction occurs below 200°C with a large excess of fluid             known from a system such as U/Pb. The elevated γOs values
with a δ18O value of, say 0‰, a plagioclase δ18O value in               in this case will not be a signal of crustal contamination, but
excess of 6‰ can result. Alteration minerals that may form at           will have been produced as a result of the interaction with a
low temperatures (clays, zeolites, oxides, carbonates) are also         hydrothermal fluid. A scattered population of Re/Os data
characterized by large mineral-water fractionation factors, and         points may offer little insight into the contamination processes
therefore, elevated whole rock δ18O values will be expected,            if an age is unavailable from other isotopic systems. Even if
even if the fluid involved in alteration and isotopic exchange          an age is available, the potential effects of hydrothermal
was relatively low-18O meteoric water. Elevated δ18O values             processes must be carefully evaluated.
should be evaluated in concert with trace elements and other
isotopic systems, but often careful petrographic examination is
sufficient to detect the involvement of hydrothermal fluids.
The measurement of mineral separates is also encouraged.
Where contaminated magma is envisioned, all primary
minerals should be characterized by anomalous δ18O values.
Owing to the kinetic exchanged effects at low temperature, all
minerals may not illustrate the same extent of exchange with
the hydrothermal fluid.
   Large degrees of fractional melting of country rocks may
lead to the development of small xenoliths in mafic rocks that
contain refractory minerals such as corundum and spinel[39,40].
In this case, elevated and widespread δ18O values in mafic
rocks may be related to unusually high degrees of assimilation.
However, such degrees of assimilation are frequently related
to repeated magma passage in a conduct system, and no single
magma pulse attains an anomalously high δ18O value.
Elevated δ18O values may also be related to diffusive oxygen
isotopic exchange near high-18O xenoliths, but such exchange              Fig. 5    Re-Os isotopic data[27] from sulfide minerals in the
zones tend to be confined to contacts with xenoliths and are                                  Voisey’s Bay deposit
                                                                         The sulfide samples define a very reasonable isochron and the
not spatially widespread[8].
                                                                         initial ratio strongly indicates that assimilation of radiogenic crust
                                                                         has occurred.
4   Radiogenic isotope systems
                                                                          Recent studies of the Bushveld Complex by Mathez and
   Radiogenic isotopes and trace elements have been the                 coworkers[45,46] have shown the importance of analyzing
standard methods of assessing contamination problems in the             individual minerals. These suggest that the variability in Pb
Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132

isotopic ratios result from the introduction of at least one other            minerals are open. When whole rocks are analyzed, the
fluid characterized by a different Pb isotopic composition at a               isotopic signal is that of the mixed populations and
time when feldspars are closed to Pb addition but sulfide                     meaningful interpretations of the data may be impossible.




                      Fig. 6      Theoretical example of the open system behavior for a sulfide deposit formed at
                                                     1.1 Ga with a chondritic Os isotopic ratio
                     Hydrothermal perturbation to the system during decay since 1.1 Ga may lead to a lowering of
                     187
                           Re/188Os ratios. In this example, hydrothermal alteration at 0.25 Ga has led to a low   187
                                                                                                                     Re/188Os
                     ratio and a γOs value, if computed for 1.1 Ga, which will be highly anomalous. An interpretation
                     of crustal assimilation at 1.1 Ga will be erroneous; isochronous samples similar to those at Voisey’s
                     Bay (Fig. 5) will be required to substantiate such an interpretation.


   A similar situation has been described by Griffin and his                  evaluation of possible perturbations to isotopic systems as a
co-workers[47,48], who noted the presence of several distinct                 result of hydrothermal alteration and exchange should be
Re/Os populations associated with different generations of                    undertaken before an interpretation involving magma
sulfide minerals in mantle peridotites. A whole rock approach                 contamination by country rocks is invoked.
in this case renders data that are not only difficult to interpret,
but essentially meaningless. Technological advances now                       Acknowledgements
permit the determination of isotopic ratios with high precision
in individual minerals using either laser ablation or ion                       This invited review paper is a contribution to the Ministry
microprobe methods. When spot analyses are not practical,                     of Education of China Project 111-B07011 awarded to China
mineral separates must be utilized if possible.                               University of Geosciences (Beijing) where the junior author is
                                                                              appointed as a Concurrent Professor. The senior author is
5   Conclusions                                                               proud to be a Guest Professor at Lanzhou University.
                                                                              Research in magmatic deposits at Indiana University is
  Both stable and radiogenic isotopic measurements can                        currently funded by grants from the National Science
provide significant information in the evaluation of country                  Foundation of China (40534020) and from the National
rock contamination in the genesis of magmatic Cu-Ni-PGE                       Science Foundation of the United States (EAR-0710910).
deposits. Sulfur isotopes provide direct evidence for the
involvement of externally derived S in ore genesis. Even in                   References
systems where the δ34S values are close to the MORB-mantle
value of 0 ± 2‰, care must be taken to ensure that potential                  [1]   Naldrett A J. The rule of sulphurization in the genesis of
contaminants are not characterized by δ34S values close to 0‰                       iron-nickel sulphide deposits of the Porcupine District, Ontario.
before dismissing processes of country rock assimilation.                           Can. Inst. Min. Bull., 1966, 4: 489−497.
Oxygen and radiogenic isotopes are also powerful tracers of                   [2]   Li C, Naldrett A J, Ripely E M. Critical factors for the
magma-country rock interaction but of the interaction of mafic                      formation of a nickel-copper deposit in an evolved magma
igneous rocks with low-temperature processes. Careful                               system: lessons from a comparison of the Pants Lake and
Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132

      Voisey’s Bay sulfide occurrences in Labrador, Can. Min.                     implications for the origin of Cu-Ni sulfide mineralization in
      Deposita, 2001, 36: 85−92.                                                  the Duluth Complex, mid-continent rift system. Con. Mineral.
[3]   Irvine T N. Crystallization sequences of the Muskox intrusion               Petrol., 2007, 154: 35−54.
      and other layered intrusions II. Origin of chromitite layers and      [17] Maier W D, Livesey T J, Barnes S J, et al. The Kabanga Ni
      similar deposits of other magmatic ores. Geochim. Cosmochim.                Sulfide Deposit, Tanzania: Geology and Geochemistry.
      Acta, 1975, 39: 991−1020.                                                   Geosciences Africa 2004, Univ. of Wits. Abstract Volume, 2005,
[4]   Carlson R W. Application of the Pt-Re-Os isotopic systems to                400−401.
      mantle geochemistry and geochronology. Lithos, 2005, 82:              [18] Godlevsky M N, Grienko L N. Some data on the isotopic
      249−272.                                                                    composition of sulfur in the sulfides of the Noril’sk deposit.
[5]   Lightfoot P C. Hawkesworth C J. Flood basalts and magmatic                  Geochemistry, 1963, 1: 335−341.
      Ni, Cu, and PGE sulfide mineralization: comparative                   [19] Gorbachev N S, Grinenko L N. The sulfur isotope ratios of the
      geochemistry of the Noril’sk (Siberian Traps) and West                      sulfides and sulfates of the Oktyabv’sk sulfide deposit, Noril’sk
      Greenland sequences. In: Mahoney J J, Coffin M F (eds.),                    region, and the problem of its origin. Geokhimiya, 1973, 8:
      Large Igneous Provinces: Continental, Oceanic and Planetary                 1127−1136.
      Flood Volcanism. Washington, D.C., American Geophysical               [20] Grinenko L N. Sources of sulfur of the nickeliferous and barren
      Union, 1997, 357-380.                                                       gabbro-dolerite intrusions of the northwest Siberian platform.
[6]   Mavrogenes J A, O’Neill H St C. The relative effects of                     International Geology Review, 1985, 28: 695−708.
      pressure, temperature, and oxygen fugacity on the solubility of       [21] Li C, Ripley E M, Naldrett A J. Compositional variations of
      sulfide in mafic magmas. Geochim. Cosmochim. Acta, 1999, 63:                olivine and sulfur isotopes in the Noril’sk and Talnakh
      1173−1180.                                                                  Intrusions, Siberia: implications for ore-forming processes in
[7]   Holzhied A, Grove T L. Sulfur saturation limits in silicate melts           dynamic magma conduits. Econ. Geol., 2003, 98: 69−86.
      and their implications for core formation scenarios for               [22] Li C, Ripley E M, Moore C R., A new genetic model for the
      terrestrial planets. Am. Mineral., 2002, 87: 227−237.                       Noril’sk-Talnakh Ni-Cu-PGE sulphide deposits. Geochim.
[8]   Ripley E M. Systematics of sulphur and oxygen isotopes in                   Cosmochim. Acta, 2007, 71: A56.
      mafic igneous rocks and related Cu-Ni-PGE mineralization. In:         [23] Ryan B, Wardle R J, Gower C F, et al. Nickel-copper sulphide
      Keays R R, Lesher C M, Lightfoot P C, et al. (eds.), Dynamic                mineralization in Labrador: the Voisey’s Bay discovery and its
      Processes in Magmatic Ore Deposits and Their Application in                 exploration implications. Gov. of Newfoundland and Labrador,
      Mineral Exploration. Geol. Assoc. Can. Short Course 13, 1999,               Current Res. Report 95-1, 1995, 177−204.
      111−158.                                                              [24] Li C, Naldrett A J. Geology and petrology of the Voisey’s Bay
[9]   Chaussidon M, Albarede F L, Sheppard S M F. Sulphur isotope                 intrusion: reaction of olivine with sulfide and silicate liquids.
      heterogeneity in the mantle from ion microprobe measurements                Lithos, 1999, 47: 1−31.
      of sulphide inclusions in diamond. Nature, 1987, 330:                 [25] Ripley E M, Park Y R, Li C, et al. Sulfur and oxygen isotopic
      242−244.                                                                    evidence of country rock contamination in the Voisey’s Bay
[10] Eldridge C S, Compston W, William I S, et al. Isotopic                       Ni-Cu-Co deposit, Labrador, Canada. Lithos, 1999, 47: 53−68.
      evidence for the involvement of recycled sediments in diamond         [26] Ripley E M Li C, Sing D. Paragneiss assimilation in the genesis
      formation. Nature, 1991, 353: 649−652.                                      of magmatic Ni-Cu-Co sulfide mineralization at Voisey’s Bay,
[11] Claypool G E, Hosler W T, Kaplan I R, et al. The age curves of               Labrador: δ34S, δ13C, and Se/S evidence. Econ. Geol., 2003, 97:
      sulfur and oxygen isotopes in marine sulfate and their mutual               1307−1318.
      interpretation. Chem. Geol., 1980, 28: 190−260.                       [27] Lambert D D, Foster J G., Frick L R, et al. Re-Os isotopic
[12] Strauss H. The isotopic composition of sedimentary sulfur                    systematics of the Voisey’s Bay Ni-Cu-Co magmatic ore system,
      through time. Palaeo, 1997, 132: 97−118.                                    Labrador, Canada. Lithos, 1999, 47: 67−88.
[13] Ripley E M, Al-Jassar T J. Sulfur and oxygen isotopic studies          [28] Farguahar J, Wing B A. Multiple sulfur isotopes and the
      of melt-country rock interaction, Babbitt Cu-Ni deposit, Duluth             evolution of the atmosphere. Earth Planet. Sci. Lett., 2003, 213:
      Complex, Minnesota. Econ. Geol., 1987, 82: 87−107.                          1−13.
[14] Andrews M S, Ripley E M, Mass transfer and sulfur fixation in          [29] Fiorentini M L, Bekker A, Rumble D, et al. Multiple S isotope
      the contact aureole of the Duluth Complex, Duka Road Cu-Ni                  study indicates footwall hydrothermal and exhalative massive
      deposit, Minnesota. Can. Mineral., 1989, 27: 293−310.                       sulfides   were    the    major   sulfur   source   for   Archean
[15] Arcuri T, Ripley E M, Hauck S A. Sulfur and oxygen isotopic                  Komatiite-hosted magmatic nickel-sulfides from Western
      studies of the interaction between xenoliths and basaltic magma             Australia and Canada. Geochim. Cosmochim. Acta, 2006, 70
      at the Babbitt and Serpentine Cu-Ni deposits, Duluth Complex,               (185): A174.
      Minnesota. Econ. Geol., 1998, 93: 1063−1075.                          [30] Penniston-Dorland S, Wing B, Brown M, et al. A petrologic
[16] Ripley E M, Taib N I, Li C, et al. Chemical and mineralogical                investigation of the Platreef Complex, South Africa, using
      heterogeneity in the basal zone of the Partridge River Intrusion:           anomalously fractionated sulfur isotopes as a tracer. EOS Trans.
Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132

      AGU, 2006, 87(36), Jt. Assem. Suppl., Abs.: V43B-06.                        xenoliths in the Voisey’s Bay Intrusion, Labrador, Canada:
[31] Ito E, White W M, Göpel C. The O, Sr, Nd, and Pb isotope                     mineralogy, reactions, partial melting and mechanisms of mass
      geochemistry of MORB. Chem. Geol., 1987, 62: 157−176.                       transfer. Geochemistry, Geophysics, Geosystems, 2006, Q05013.
[32] Harmon R S, Hoefs J. Oxygen isotope heterogeneity of the                     doi10.1029/2005GC001184.
                                     18
      mantle deduced from global          O systematics of basalts from     [41] Faure G. Origin of Igneous Rocks: The Isotopic Evidence.
      different geotectonic settings. Con. Mineral. Petrol., 1995, 120:           Berlin: Springer, 2001, 496.
      95−114.                                                               [42] Hanson G N. An approach to trace element modeling using a
[33] Eiler J M. Oxygen isotope variations in basaltic lavas and upper             simple igneous system as an example. In: Lipin B R, McKay G
      mantle rocks. In: Valley J W, Cole D R (eds.), Stable Isotope               A (eds.), Geochemistry and Mineralogy of Rare Earth Elements.
      Geochemistry. Min. Soc. Amer. Rev. in Min. and Geochem.                     Min. Soc. Amer. Rev. Min., 1989, 21: 79−98.
      2001, 43: 319−364.                                                    [43] Dickin A P. Radiogenic Isotope Geology. Cambridge University
[34] Leitch A M. Analog experiments on melting and contamination                  Press, 1995, 479.
      at the roof and walls of magma chambers. Jour. Volc. Geotherm.        [44] Depaolo D. Trace element effects of combined wall rock
      Res., 2004, 129: 173−197.                                                   assimilation and fractional crystallization. Earth Planet. Sci.
[35] Taylor H P Jr. The application of oxygen and hydrogen isotope                Lett., 1981, 53: 189−202.
      studies to problems of hydrothermal alteration and ore                [45] Mathez E A, Waight T. Lead isotopic disequilibrium between
      deposition. Econ. Geol., 1974, 69: 843−883.                                 sulfide and plagioclase in the Bushveld Complex and them
[36] Criss R E, Taylor H P Jr. Meteoric-hydrothermal systems. In:                 chemical evolution of large layered intrusions. Geochim.
      Valley J W, Taylor S P Jr, O’Neil J R (eds.), Stable Isotopes in            Cosmochim. Acta, 2003, 62: 1875−1888.
      High Temperature Geological Processes. Min. Soc. Amer. Rev.           [46] Mathez E A, Kent A J R. Variable initial Pb isotopic
      Min., 1986, 16: 373−424.                                                    compositions of rocks associated with the UG2 chromitite,
[37] Criss R E. Principles of Stable Isotope Distribution. Oxford                 eastern Bushveld Complex. Geochim. Cosmochim. Acta, in
      University Press, 1999, 254.                                                press.
[38] O’Neil J R, Taylor H P Jr. The oxygen isotope and cation               [47] Alard O, Griffin W L, Pearson N J, et al. New insights into the
      exchange chemistry of feldspars. Am. Mineral, 1969, 52:                     Re-Os systematics of subcontinental lithospheric mantle from
      1414−1437.                                                                  in situ analysis of sulphides. Earth Planet. Sci. Lett., 2002, 203:
[39] Li C, Naldrett A J. Melting reactions of gneissic inclusions with            651−663.
      enclosing magma at Voisey’s Bay, Labrador: implications with          [48] Griffin W L, Graham S, O’Reilly S Y. Lithosphere evolution
      respect to ore genesis. Econ. Geol., 2000, 95: 801−814.                     beneath the Kaapraal Craton: Re-Os systematics of sulfides in
[40] Mariga J, Ripley E M, Li C. Petrogenetic evolution of gneissic               mantle-derived peridotites. Chem. Geol., 2004, 208: 89−118.

More Related Content

PDF
Accessory minerals importance in granite petrology: a review and case studies
PDF
Hydrothermal Alteration in the Propylitic Suite at Butte, Montana (Haggart 2014)
PDF
Rose and Brenan (2001)
PDF
Magnesium isotope evidence that accretional vapour loss shapes planetary comp...
PDF
Study of the effect of Illite and Kaolinite on low salinity water injection
PPTX
Role of Trace Elements In Petrogenesis
PDF
Assessment of fe(ii), fe(iii) and al(iii) in a mineralogical profile of gabon...
DOCX
THE ORE-BEARING FLUIDS
Accessory minerals importance in granite petrology: a review and case studies
Hydrothermal Alteration in the Propylitic Suite at Butte, Montana (Haggart 2014)
Rose and Brenan (2001)
Magnesium isotope evidence that accretional vapour loss shapes planetary comp...
Study of the effect of Illite and Kaolinite on low salinity water injection
Role of Trace Elements In Petrogenesis
Assessment of fe(ii), fe(iii) and al(iii) in a mineralogical profile of gabon...
THE ORE-BEARING FLUIDS

What's hot (20)

PPTX
Ronald Ng M.Sc. thesis presentation (Oct 2012)
PDF
GRAPHIC QUARTZ-FELDSPAR INTERGROWTHS IN PEGMATITES: DIFFUSION AND GROWTH KINE...
PPTX
Wall rock alteration
PDF
Uranium Geochemistry in Peraluminous Leucogranites of Wadi El-Shallal Area, S...
PDF
Primary geochemical haloes in prospecting for gold deposits, umm rus mine, ea...
PDF
About the Interaction between water and Pb, Zn, Cu, Cd, Fe, Mn, Ba Mineralize...
PDF
Exploration in the new Millennium 5th Decennial International Conference on M...
PDF
A D Hollinsworth 2015
PPTX
Classification of ore deposits
PPT
Igneous processes
PDF
Nature10582
PDF
Ore bearing fluids
PDF
Metals accumulation and As releasing during interaction of clay and iron mine...
PDF
PETROGENESIS OF GRANITOID ROCKS AND ORIGIN OF URANIUM MINERALIZATIONS OF UM S...
PPT
Economic geology - Magmatic ore deposits_1
PDF
PPTX
Formation of lithified micritic laminae in modern marine stromatolites
PDF
The effect of liming on the acidity level of Dystric cambisol and the content...
DOCX
Hydrothermal deposits by VARINDRA
PDF
Conformational analysis of humic acids from amazonian dark earth soils observ...
Ronald Ng M.Sc. thesis presentation (Oct 2012)
GRAPHIC QUARTZ-FELDSPAR INTERGROWTHS IN PEGMATITES: DIFFUSION AND GROWTH KINE...
Wall rock alteration
Uranium Geochemistry in Peraluminous Leucogranites of Wadi El-Shallal Area, S...
Primary geochemical haloes in prospecting for gold deposits, umm rus mine, ea...
About the Interaction between water and Pb, Zn, Cu, Cd, Fe, Mn, Ba Mineralize...
Exploration in the new Millennium 5th Decennial International Conference on M...
A D Hollinsworth 2015
Classification of ore deposits
Igneous processes
Nature10582
Ore bearing fluids
Metals accumulation and As releasing during interaction of clay and iron mine...
PETROGENESIS OF GRANITOID ROCKS AND ORIGIN OF URANIUM MINERALIZATIONS OF UM S...
Economic geology - Magmatic ore deposits_1
Formation of lithified micritic laminae in modern marine stromatolites
The effect of liming on the acidity level of Dystric cambisol and the content...
Hydrothermal deposits by VARINDRA
Conformational analysis of humic acids from amazonian dark earth soils observ...
Ad

Similar to Applications of stable and radiogenic isotopes to magmatic cu ni-pge deposits examples and cautions (20)

PDF
Brenan and Rose (2002)
PDF
minerals-12-00981-v2.pdf
PDF
Sulphr isotope
PDF
GEOSCIENCE 08 Carey et al 437-444
PDF
Igneous petrology.....Masoom
PDF
EconGeology3
PPTX
Heavy metals in langkawi corals
PDF
S-isotope-PNAS-2014-Shaheen-1406315111 with cover
PDF
Review clay flotation - Jeldres et al., 2019.pdf
PDF
B044010019
PDF
A geochemical assessment of soils within oke ogun area, southwestern nigeria
PDF
Star 2011 Final
PDF
A comparative study on adsorption behavior of heavy metal elements onto soil ...
PDF
Iron isotopes trace primordial magma ocean cumulates melting in Earth’s upper...
PDF
D044027038
PDF
Interfacial water structure and the wetting of mineral surfaces.pdf
PDF
Induced polarization of disseminated electronically conductive minerals
PDF
Einaudi 1994 High Sulfidation and Low Sulfidation Porphyry CopperSkarn System...
PDF
Mud_volcanoes_in_the_South_Caspian_basin_Nature_an.pdf
Brenan and Rose (2002)
minerals-12-00981-v2.pdf
Sulphr isotope
GEOSCIENCE 08 Carey et al 437-444
Igneous petrology.....Masoom
EconGeology3
Heavy metals in langkawi corals
S-isotope-PNAS-2014-Shaheen-1406315111 with cover
Review clay flotation - Jeldres et al., 2019.pdf
B044010019
A geochemical assessment of soils within oke ogun area, southwestern nigeria
Star 2011 Final
A comparative study on adsorption behavior of heavy metal elements onto soil ...
Iron isotopes trace primordial magma ocean cumulates melting in Earth’s upper...
D044027038
Interfacial water structure and the wetting of mineral surfaces.pdf
Induced polarization of disseminated electronically conductive minerals
Einaudi 1994 High Sulfidation and Low Sulfidation Porphyry CopperSkarn System...
Mud_volcanoes_in_the_South_Caspian_basin_Nature_an.pdf
Ad

Recently uploaded (20)

PDF
Microsoft Solutions Partner Drive Digital Transformation with D365.pdf
PPTX
observCloud-Native Containerability and monitoring.pptx
PDF
ENT215_Completing-a-large-scale-migration-and-modernization-with-AWS.pdf
PPTX
Benefits of Physical activity for teenagers.pptx
PDF
Taming the Chaos: How to Turn Unstructured Data into Decisions
PPTX
Modernising the Digital Integration Hub
PDF
STKI Israel Market Study 2025 version august
PDF
Hindi spoken digit analysis for native and non-native speakers
PPTX
MicrosoftCybserSecurityReferenceArchitecture-April-2025.pptx
PDF
NewMind AI Weekly Chronicles – August ’25 Week III
PDF
August Patch Tuesday
PPTX
Tartificialntelligence_presentation.pptx
PPTX
Chapter 5: Probability Theory and Statistics
PDF
Getting Started with Data Integration: FME Form 101
PDF
A novel scalable deep ensemble learning framework for big data classification...
PDF
Hybrid model detection and classification of lung cancer
PDF
Architecture types and enterprise applications.pdf
PPT
Geologic Time for studying geology for geologist
PDF
1 - Historical Antecedents, Social Consideration.pdf
PDF
Assigned Numbers - 2025 - Bluetooth® Document
Microsoft Solutions Partner Drive Digital Transformation with D365.pdf
observCloud-Native Containerability and monitoring.pptx
ENT215_Completing-a-large-scale-migration-and-modernization-with-AWS.pdf
Benefits of Physical activity for teenagers.pptx
Taming the Chaos: How to Turn Unstructured Data into Decisions
Modernising the Digital Integration Hub
STKI Israel Market Study 2025 version august
Hindi spoken digit analysis for native and non-native speakers
MicrosoftCybserSecurityReferenceArchitecture-April-2025.pptx
NewMind AI Weekly Chronicles – August ’25 Week III
August Patch Tuesday
Tartificialntelligence_presentation.pptx
Chapter 5: Probability Theory and Statistics
Getting Started with Data Integration: FME Form 101
A novel scalable deep ensemble learning framework for big data classification...
Hybrid model detection and classification of lung cancer
Architecture types and enterprise applications.pdf
Geologic Time for studying geology for geologist
1 - Historical Antecedents, Social Consideration.pdf
Assigned Numbers - 2025 - Bluetooth® Document

Applications of stable and radiogenic isotopes to magmatic cu ni-pge deposits examples and cautions

  • 1. EARTH SCIENCE FRONTIERS Volume 14, Issue 5, September 2007 Online English edition of the Chinese language journal Cite this article as: Earth Science Frontiers, 2007, 14(5):124–132. RESEARCH PAPER Applications of Stable and Radiogenic Isotopes to Magmatic Cu-Ni-PGE Deposits: Examples and Cautions Edward M. RIPLEY1,2,∗, Chusi LI1,2 1 Department of Geological Sciences, Indiana University, Bloomington, IN 47401, USA 2 State Key Laboratory of Geological Processes and Mineral Resources, China University of Geosciences, Beijing 100083, China Abstract: Measurements of S, O, and radiogenic isotope ratios are all potentially powerful tracers of magma interaction with country rocks and the importance of assimilation processes in the genesis of magmatic Ni-Cu-PGE deposits. Sulfur isotope measurements of deposits such as those in the 1.1 Ga Duluth Complex, the Permo-Triassic intrusions of the Noril’sk area, and the 1.4 Ga Kabanga intrusions provide evidence for the derivation of S from both sulfide-and sulfate-bearing country rocks. The 1.3 Ga Voisey’s Bay deposit provides an example where δ34S values of the ores commonly fall within the accepted mantle range of 0 ± 2‰, but detailed studies of Proterozoic metasedimentary country rocks show that their weighted average δ34S value is also within this range. A thorough knowledge of the isotopic compositions of potential contaminants is essential for a proper evaluation of the role of country rock derived S in ore formation. When O and radiogenic isotopic measurements are employed as tracers to evaluate magma-country rock interaction, it is essential to provide evidence that open system processes have not perturbed the isotopic systematics. Low-temperature hydrothermal processes can mask evidence of high-temperature processes in the oxygen isotope system and in radiogenic systems such as Re/Os and Pb where involved elements may be mobile under hydrothermal conditions, or host phases may close to exchange and uptake at different temperatures. Careful petrographic observation and analyses of individual minerals may be required before the models involving the contamination of magmas by country rocks can be meaningfully applied. Key Words: S isotopes; O isotopes; radiogenic isotopes; magma-country rock interaction importance of S assimilation in promoting sulfide saturation in 1 Introduction a mafic magma[6,7]. This review is intended to highlight some of the Isotopic analyses are a powerful tool for evaluating the fundamental applications of isotopic measurements to the potential interaction between mantle-derived magmas and assessment of magma contamination as applied to the studies crustal rocks. Both stable and radiogenic isotopic systems of magmatic ore genesis and cautionary notes that need to be have been applied to the studies of magma contamination and origin. It has become apparent that the contamination of mafic considered to avoid potential pitfalls in the interpretation of magma by crust plays a significant role in the generation of isotopic data. Examples from well-known Cu-Ni-PGE large, sulfide-rich magmatic Cu-Ni-PGE deposits[1,2]. deposits will be utilized to illustrate key concepts. Interaction between siliceous country rocks and mafic magma has also been suggested as an important prerequisite for the 2 Sulfur isotopes development of chromite deposits in layered intrusions[3,4] and sulfide deposits in some intrusions[5]. Experimental studies Owing to the relatively large difference in δ34S values that emphasize the role of pressure as a control of the between mantle-derived magma (MORB as the principal solubility of sulfur in mafic magmas have also stressed the example) and either sulfide or sulfate minerals in sedimentary Received date: 2007-08-20; Accepted date: 2007-09-11. * Corresponding author. E-mail: ripley@indiana.edu Foundation item: Supported by the Ministry of Education of China project (111-B07011); the National Natural Science Foundation of China (40534020); the National Science Foundation of the United States (EAR-0710910). Copyright © 2007, China University of Geosciences (Beijing) and Peking University, Published by Elsevier B.V. All rights reserved.
  • 2. Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132 rocks, sulfur isotopic studies have provided important insights Tanzania[17]. Mica schists and graphitic schists contain into the processes of magma contamination. Mantle sulfur is pyrrhotite with δ34S values ranging from -2‰ to 24‰. Sulfide commonly taken as having δ34S values between -2‰ and 2‰. mineralization in mafic to ultramafic intrusive rocks is This range is defined largely from the analyses of sulfide characterized by a range of δ34S values between 10‰ and minerals in meteorites and sulfide-bearing basalts (see 24‰, again lending strong evidence to a genetic model that Ripley[8] for a review of early S isotope studies). It should involves magma contamination by country rocks. always be kept in mind, however, that sulfide minerals in diamonds from kimberlites are characterized by a wide range of δ34S values from -11‰ to +14‰[9,10], suggesting that the deep mantle may be far more heterogeneous in terms of δ34S values than the shallow mantle appears to be. Seawater δ34S has varied between ~10‰ and 30‰ from the Neoproterozoic to the present[11,12]. Sulfide minerals are produced from seawater sulfate owing to reduction that is often related to the life cycles of sulfate-reducing bacteria. Depending on the variables related primarily to the rate of sulfate reduction, the sulfide produced may be characterized by positive or negative δ34S values. Pyrite that forms owing to the reaction of produced sulfide and reactive iron may show a wide range in δ34S values, typically from values close to the source sulfate (e.g. ~ 20‰), to values as low as -50‰, or occasionally even less. Sulfate minerals that form via evaporative processes show little isotopic fractionation relative to the source sulfate, and hence are generally characterized by δ34S values between ~10‰ and 30‰. It is clear that because of the potentially large S isotopic difference between sedimentary S-bearing minerals and mantle-derived magmas, assimilation processes may be elucidated using S isotope measurements. Three examples will serve to illustrate the basic principles. 2.1 Duluth Complex, Minnesota Cu-Ni sulfide mineralization in the 1.1 Ga Midcontinent Rift-related Duluth Complex represents one of the best-known examples of major sulfur derivation from country rocks. The Proterozoic (~1.85 Ga) Virginia Formation is a pelitic sequence that contains layers that are sulfidic and carbonaceous, and are characterized by pyrite δ34S values ranging from 0 to 30‰ (Fig. 1)[13–15]. Metamorphism of the sulfide-bearing sedimentary rocks has led to the production of pyrrhotite and the generation of an H2S-bearing, carbonic fluid. Magma that formed the Duluth Complex interacted with –15] various levels of the sedimentary sequence, assimilated the Fig. 1 Sulfur isotope data[13 from sulfide mineralization in the sedimentary-derived H2S, and crystallized the sulfide minerals Duluth Complex and Virginia Formation country rocks that were characterized by elevated δ34S values ranging from ~0 to 20‰ (Fig. 1). Various mixing models involving external 2.2 Noril’sk, Siberia S and mantle S can be derived[16], but the importance of externally derived S is highlighted by the δ34S values and the Another well-known example where δ34S values of sulfide internally consistent models of fluid production in the ores strongly suggest that assimilation has been operative is metamorphic rocks and localized partial melting. the world famous mineralization at Noril’sk, in the A similar geologic setting in terms of mafic intrusions Permo-Triassic Siberian flood basalt province. Early sulfur found hosted by pelitic and graphitic metasedimentary rocks isotope studies[18–20] suggested that assimilation of evaporites has been described for the large 1.4 Ga Kabanga Ni-deposit in that are common in the stratigraphic sequence occurred, or
  • 3. Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132 that sulfur may have been derived from pyritic shales and envisioned, or hydrothermal leaching of sulfate followed by coals, or H2S derived from the coals. Li and his partial reduction, or a process involving the dissolution, rather co-workers[21,22] developed slightly different models to explain than melting, of sulfate in mafic magma. The observation of the anomalous S isotopic signatures (Fig. 2); however, the abundant wollastonite in contact aureole rocks suggests that probable involvement of sulfate as a contaminant is reactions similar to CaSO4 + SiO2 + H2O = CaSiO3 + H2S + noteworthy. The assimilation of sulfate via partial melting by a 2O2 occurred, and that sulfate may have been reduced to basatic magma is difficult to conceive owing to the high sulfide prior to incorporation into the magma. The important melting points of sulfate minerals (typically ranging between point here is that sulfate minerals may serve as S sources, but 1360 and 1450°C). Therefore, either a fluxing mechanism is the mechanism of S assimilation must be carefully evaluated. Fig. 2 Sulfur isotope data[21,22] from the Noril’sk area process, and positive values indicating a relatively rapid rate 2.3 Voisey’s Bay of sulfate reduction, possibly under closed system conditions. A third example where the interpretation of δ34S values is The weighted average δ34S value of the sulfide assemblages in critical is the 1.3 Ga Voisey’s Bay Ni-Co-Cu deposit in the Tasiuyak Gneiss is -1.3‰. This value is very similar to the Labrador, Canada. The deposit is associated with the Nain average value of the sulfide ores and illustrates a very Plutonic Suite, which was emplaced along a boundary important principle. It is important to precisely know the δ34S between crustal blocks, one comprised predominantly of values of the proposed contaminant before concluding. In the Archean rocks and the other dominated by Proterozoic rocks. case of Voisey’s Bay, even though the δ34S values are close to Sulfide mineralization is found in troctolitic to gabbroic rocks. 0‰, the additional study clearly showed that S could have Country rocks include the Proterozoic Tasiuyak Gneiss, which been derived from the country rocks, in agreement with the is locally sulfidic and graphitic. Several investigators[23,24] widespread geological evidence, which suggests that S proposed that large amounts of S in the mineralization may assimilation occurred. have been derived from the Tasiuyak Gneiss. Initial S isotopic Voisey’s Bay is also an excellent example where supporting studies[25] determined that the bulk of the sulfide radiogenic isotopic data was extremely valuable. The Re/Os mineralization was characterized by δ34S values ranging isotopic data (see below)[27] clearly indicated that the magma between 2‰ and -4‰ (Fig. 3). The δ34S values in the sulfide from which the sulfides were derived assimilated crustal ores were significantly depleted only locally near the contact material. Carbon isotopic data also indicated that assimilation with the Tasiuyak Gneiss, which contained pyrrhotite with of sedimentary country rocks had occurred[26], as did the S/Se negative δ34S values. However, additional studies[26] ratios of the sulfide mineralization. The supporting data determined that pyrrhotite in the Tasiuyak Gneiss was strengthens the promise that even though the δ34S values of characterized by a wide range in δ34S values, including both the ores are close to those of mantle sulfides, a large negative values suggesting an open system sulfate reduction component of the S was derived from country rocks.
  • 4. Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132 derived in part from underlying sedimentary rocks. The anomalous isotopic signature illustrated in the δ34S-δ33S plot strongly suggests that the interpretation was correct. Fig. 4 δ33S versus δ34S for sulfide samples from the Stillwater Complex, Voisey’s Bay, and Alexo IAEA S-1 is an international standard. Note the anomalous [26] Fig. 3 Sulfur isotope data from sulfide mineralization at character of the Alexo samples, strongly supporting the premise Voisey’s Bay and sulfide minerals in the Tasiuyak Gneiss that sulfur from Archean-aged metasedimentary country rocks has country rocks been involved in ore genesis. A final point with respect to the use of S isotopes as 3 Oxygen isotopes indicators of crustal contamination regards Archean-aged sulfide deposits. Ripley[8] discussed that during the Archean Mantle-derived mafic magmas are generally characterized bacterial sulfate reduction, then if operative, appears to have by δ18O values in the range of 5‰ to 7‰[31–33]. Sedimentary generated considerably smaller ranges in δ34S values of rocks and crustally derived granitic rocks are normally reduced sulfur species, and values close to 0‰. For this characterized by considerably higher δ18O values, typically in reason, the determination of externally derived S in Archean the range of 8‰ to 30‰. Owing to these large differences, deposits, or in younger deposits that may have interacted with oxygen isotopes may also be sensitive indications of country Archean rocks, may be difficult using δ34S values. However, a rock assimilation. General features of how oxygen isotopes process known as mass-independent fractionation (MIF) of S may be applied in the study of magmatic ore deposits have isotopes has led to variations between δ33S and δ34S, which are been discussed by Ripley[8] and will not be repeated here. distinct from most terrestrial processes that fractionate S However, common difficulties with a straight-forward isotopes and are referred to as mass-dependent fractionation application of δ18O values to problems of assimilation by processes. Farguhar and Wing[28] discussed the observation mafic magmas must be highlighted. that sedimentary pyrite of Archean age is characterized by It is often assumed that elevated δ18O values in mafic δ33S-δ34S relations that are indicative of MIF. The process has igneous rocks (e.g. >7‰) are indicative of contamination and been attributed to low oxygen levels in the Archean assimilation of high-18O country rocks. Although this may atmosphere and photochemical reactions involving SO2, but indeed be the case, careful petrographic examination of the the controversy with respect to the origin of the process rocks, along with an evaluation of the magnitude of the δ18O continues. The record of MIF may be present in magmas that shift, should be undertaken. For example, consider a mafic have assimilated S of the Archean age[29,30], and future studies rock with bulk δ18O value of 9‰, a contaminant with a δ18O using the technique may be able to distinguish the importance value of 16‰, and a presumed magma initial value of 5.5‰. of S derived from Archean rocks in ore genesis. Figure 4 is a Assuming that the concentrations of O in the contaminant and δ33S-δ34S plot showing sulfide samples from Voisey’s Bag (1.3 the magma are similar, a value of 9‰ in the mafic rock will be Ga) and Alexo (a komatiite-associated, ~2.7 Ga deposit). suggested in ~35‰ bulk contamination. This degree of Naldrett[1] proposed that S in the Alexo deposit may have been contamination is unreasonable in most situations because of
  • 5. Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132 the amount of heat required to achieve such a large geosciences. Several excellent text books and review articles contaminant/magma ratio[34]. Unless xenoliths provide address these methodologies[41–44]. Here we wish to emphasize evidence for such high degrees of contamination (the presence that care must be taken to assure that closed-system conditions of refractory oxides such as corundum and spinel), other prevailed before conclusions regarding contamination history methods that can explain the 18O-elevation must be explored. can be made. For example, the Re/Os isotopic data from Obviously knowledge of the δ18O values of potential sulfide minerals in the Voisey’s Bay deposit[27] mentioned contaminants is extremely important. The assimilation of very above produces an isochron, which is consistent with the age high-18O siliceous rocks, carbonates, or hydrothermally altered determined via U-Pb geochronology Fig. 5. The elevated felsic rocks can be responsible for elevated δ18O values in initial 187Os/188Os ratio clearly indicates that crustal mafic igneous rocks and can satisfy reasonable thermal contamination occurred, the sulfide system was well-mixed, constraints. and that closed system conditions have generally prevailed. In The most common cause of elevated δ18O values in mafic contrast, in situations where a reasonable isochron is not rocks that is often overlooked is low-temperature produced, an interpretation involving crustal contamination hydrothermal alteration and isotopic exchange. Low- must proceed with caution. Figure 6 illustrates a case where an temperature isotopic exchange involving a fluid has been initially chondritic system interacts with a fluid after discussed by several investigators[35–37]; one pertinent example crystallization. Loss of rhenium during the interaction may will be presented here. At temperatures below ~200°C, the significantly lower the 187Re/188Os ratio. Continued decay may oxygen isotopic fractionation between plagioclase and water is then produce samples with 187Os/188Os ratios that will generate greater than 6‰[38]. For a situation where the fluid-rock anomalous γOs values, assuming that an age for the rocks is interaction occurs below 200°C with a large excess of fluid known from a system such as U/Pb. The elevated γOs values with a δ18O value of, say 0‰, a plagioclase δ18O value in in this case will not be a signal of crustal contamination, but excess of 6‰ can result. Alteration minerals that may form at will have been produced as a result of the interaction with a low temperatures (clays, zeolites, oxides, carbonates) are also hydrothermal fluid. A scattered population of Re/Os data characterized by large mineral-water fractionation factors, and points may offer little insight into the contamination processes therefore, elevated whole rock δ18O values will be expected, if an age is unavailable from other isotopic systems. Even if even if the fluid involved in alteration and isotopic exchange an age is available, the potential effects of hydrothermal was relatively low-18O meteoric water. Elevated δ18O values processes must be carefully evaluated. should be evaluated in concert with trace elements and other isotopic systems, but often careful petrographic examination is sufficient to detect the involvement of hydrothermal fluids. The measurement of mineral separates is also encouraged. Where contaminated magma is envisioned, all primary minerals should be characterized by anomalous δ18O values. Owing to the kinetic exchanged effects at low temperature, all minerals may not illustrate the same extent of exchange with the hydrothermal fluid. Large degrees of fractional melting of country rocks may lead to the development of small xenoliths in mafic rocks that contain refractory minerals such as corundum and spinel[39,40]. In this case, elevated and widespread δ18O values in mafic rocks may be related to unusually high degrees of assimilation. However, such degrees of assimilation are frequently related to repeated magma passage in a conduct system, and no single magma pulse attains an anomalously high δ18O value. Elevated δ18O values may also be related to diffusive oxygen isotopic exchange near high-18O xenoliths, but such exchange Fig. 5 Re-Os isotopic data[27] from sulfide minerals in the zones tend to be confined to contacts with xenoliths and are Voisey’s Bay deposit The sulfide samples define a very reasonable isochron and the not spatially widespread[8]. initial ratio strongly indicates that assimilation of radiogenic crust has occurred. 4 Radiogenic isotope systems Recent studies of the Bushveld Complex by Mathez and Radiogenic isotopes and trace elements have been the coworkers[45,46] have shown the importance of analyzing standard methods of assessing contamination problems in the individual minerals. These suggest that the variability in Pb
  • 6. Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132 isotopic ratios result from the introduction of at least one other minerals are open. When whole rocks are analyzed, the fluid characterized by a different Pb isotopic composition at a isotopic signal is that of the mixed populations and time when feldspars are closed to Pb addition but sulfide meaningful interpretations of the data may be impossible. Fig. 6 Theoretical example of the open system behavior for a sulfide deposit formed at 1.1 Ga with a chondritic Os isotopic ratio Hydrothermal perturbation to the system during decay since 1.1 Ga may lead to a lowering of 187 Re/188Os ratios. In this example, hydrothermal alteration at 0.25 Ga has led to a low 187 Re/188Os ratio and a γOs value, if computed for 1.1 Ga, which will be highly anomalous. An interpretation of crustal assimilation at 1.1 Ga will be erroneous; isochronous samples similar to those at Voisey’s Bay (Fig. 5) will be required to substantiate such an interpretation. A similar situation has been described by Griffin and his evaluation of possible perturbations to isotopic systems as a co-workers[47,48], who noted the presence of several distinct result of hydrothermal alteration and exchange should be Re/Os populations associated with different generations of undertaken before an interpretation involving magma sulfide minerals in mantle peridotites. A whole rock approach contamination by country rocks is invoked. in this case renders data that are not only difficult to interpret, but essentially meaningless. Technological advances now Acknowledgements permit the determination of isotopic ratios with high precision in individual minerals using either laser ablation or ion This invited review paper is a contribution to the Ministry microprobe methods. When spot analyses are not practical, of Education of China Project 111-B07011 awarded to China mineral separates must be utilized if possible. University of Geosciences (Beijing) where the junior author is appointed as a Concurrent Professor. The senior author is 5 Conclusions proud to be a Guest Professor at Lanzhou University. Research in magmatic deposits at Indiana University is Both stable and radiogenic isotopic measurements can currently funded by grants from the National Science provide significant information in the evaluation of country Foundation of China (40534020) and from the National rock contamination in the genesis of magmatic Cu-Ni-PGE Science Foundation of the United States (EAR-0710910). deposits. Sulfur isotopes provide direct evidence for the involvement of externally derived S in ore genesis. Even in References systems where the δ34S values are close to the MORB-mantle value of 0 ± 2‰, care must be taken to ensure that potential [1] Naldrett A J. The rule of sulphurization in the genesis of contaminants are not characterized by δ34S values close to 0‰ iron-nickel sulphide deposits of the Porcupine District, Ontario. before dismissing processes of country rock assimilation. Can. Inst. Min. Bull., 1966, 4: 489−497. Oxygen and radiogenic isotopes are also powerful tracers of [2] Li C, Naldrett A J, Ripely E M. Critical factors for the magma-country rock interaction but of the interaction of mafic formation of a nickel-copper deposit in an evolved magma igneous rocks with low-temperature processes. Careful system: lessons from a comparison of the Pants Lake and
  • 7. Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132 Voisey’s Bay sulfide occurrences in Labrador, Can. Min. implications for the origin of Cu-Ni sulfide mineralization in Deposita, 2001, 36: 85−92. the Duluth Complex, mid-continent rift system. Con. Mineral. [3] Irvine T N. Crystallization sequences of the Muskox intrusion Petrol., 2007, 154: 35−54. and other layered intrusions II. Origin of chromitite layers and [17] Maier W D, Livesey T J, Barnes S J, et al. The Kabanga Ni similar deposits of other magmatic ores. Geochim. Cosmochim. Sulfide Deposit, Tanzania: Geology and Geochemistry. Acta, 1975, 39: 991−1020. Geosciences Africa 2004, Univ. of Wits. Abstract Volume, 2005, [4] Carlson R W. Application of the Pt-Re-Os isotopic systems to 400−401. mantle geochemistry and geochronology. Lithos, 2005, 82: [18] Godlevsky M N, Grienko L N. Some data on the isotopic 249−272. composition of sulfur in the sulfides of the Noril’sk deposit. [5] Lightfoot P C. Hawkesworth C J. Flood basalts and magmatic Geochemistry, 1963, 1: 335−341. Ni, Cu, and PGE sulfide mineralization: comparative [19] Gorbachev N S, Grinenko L N. The sulfur isotope ratios of the geochemistry of the Noril’sk (Siberian Traps) and West sulfides and sulfates of the Oktyabv’sk sulfide deposit, Noril’sk Greenland sequences. In: Mahoney J J, Coffin M F (eds.), region, and the problem of its origin. Geokhimiya, 1973, 8: Large Igneous Provinces: Continental, Oceanic and Planetary 1127−1136. Flood Volcanism. Washington, D.C., American Geophysical [20] Grinenko L N. Sources of sulfur of the nickeliferous and barren Union, 1997, 357-380. gabbro-dolerite intrusions of the northwest Siberian platform. [6] Mavrogenes J A, O’Neill H St C. The relative effects of International Geology Review, 1985, 28: 695−708. pressure, temperature, and oxygen fugacity on the solubility of [21] Li C, Ripley E M, Naldrett A J. Compositional variations of sulfide in mafic magmas. Geochim. Cosmochim. Acta, 1999, 63: olivine and sulfur isotopes in the Noril’sk and Talnakh 1173−1180. Intrusions, Siberia: implications for ore-forming processes in [7] Holzhied A, Grove T L. Sulfur saturation limits in silicate melts dynamic magma conduits. Econ. Geol., 2003, 98: 69−86. and their implications for core formation scenarios for [22] Li C, Ripley E M, Moore C R., A new genetic model for the terrestrial planets. Am. Mineral., 2002, 87: 227−237. Noril’sk-Talnakh Ni-Cu-PGE sulphide deposits. Geochim. [8] Ripley E M. Systematics of sulphur and oxygen isotopes in Cosmochim. Acta, 2007, 71: A56. mafic igneous rocks and related Cu-Ni-PGE mineralization. In: [23] Ryan B, Wardle R J, Gower C F, et al. Nickel-copper sulphide Keays R R, Lesher C M, Lightfoot P C, et al. (eds.), Dynamic mineralization in Labrador: the Voisey’s Bay discovery and its Processes in Magmatic Ore Deposits and Their Application in exploration implications. Gov. of Newfoundland and Labrador, Mineral Exploration. Geol. Assoc. Can. Short Course 13, 1999, Current Res. Report 95-1, 1995, 177−204. 111−158. [24] Li C, Naldrett A J. Geology and petrology of the Voisey’s Bay [9] Chaussidon M, Albarede F L, Sheppard S M F. Sulphur isotope intrusion: reaction of olivine with sulfide and silicate liquids. heterogeneity in the mantle from ion microprobe measurements Lithos, 1999, 47: 1−31. of sulphide inclusions in diamond. Nature, 1987, 330: [25] Ripley E M, Park Y R, Li C, et al. Sulfur and oxygen isotopic 242−244. evidence of country rock contamination in the Voisey’s Bay [10] Eldridge C S, Compston W, William I S, et al. Isotopic Ni-Cu-Co deposit, Labrador, Canada. Lithos, 1999, 47: 53−68. evidence for the involvement of recycled sediments in diamond [26] Ripley E M Li C, Sing D. Paragneiss assimilation in the genesis formation. Nature, 1991, 353: 649−652. of magmatic Ni-Cu-Co sulfide mineralization at Voisey’s Bay, [11] Claypool G E, Hosler W T, Kaplan I R, et al. The age curves of Labrador: δ34S, δ13C, and Se/S evidence. Econ. Geol., 2003, 97: sulfur and oxygen isotopes in marine sulfate and their mutual 1307−1318. interpretation. Chem. Geol., 1980, 28: 190−260. [27] Lambert D D, Foster J G., Frick L R, et al. Re-Os isotopic [12] Strauss H. The isotopic composition of sedimentary sulfur systematics of the Voisey’s Bay Ni-Cu-Co magmatic ore system, through time. Palaeo, 1997, 132: 97−118. Labrador, Canada. Lithos, 1999, 47: 67−88. [13] Ripley E M, Al-Jassar T J. Sulfur and oxygen isotopic studies [28] Farguahar J, Wing B A. Multiple sulfur isotopes and the of melt-country rock interaction, Babbitt Cu-Ni deposit, Duluth evolution of the atmosphere. Earth Planet. Sci. Lett., 2003, 213: Complex, Minnesota. Econ. Geol., 1987, 82: 87−107. 1−13. [14] Andrews M S, Ripley E M, Mass transfer and sulfur fixation in [29] Fiorentini M L, Bekker A, Rumble D, et al. Multiple S isotope the contact aureole of the Duluth Complex, Duka Road Cu-Ni study indicates footwall hydrothermal and exhalative massive deposit, Minnesota. Can. Mineral., 1989, 27: 293−310. sulfides were the major sulfur source for Archean [15] Arcuri T, Ripley E M, Hauck S A. Sulfur and oxygen isotopic Komatiite-hosted magmatic nickel-sulfides from Western studies of the interaction between xenoliths and basaltic magma Australia and Canada. Geochim. Cosmochim. Acta, 2006, 70 at the Babbitt and Serpentine Cu-Ni deposits, Duluth Complex, (185): A174. Minnesota. Econ. Geol., 1998, 93: 1063−1075. [30] Penniston-Dorland S, Wing B, Brown M, et al. A petrologic [16] Ripley E M, Taib N I, Li C, et al. Chemical and mineralogical investigation of the Platreef Complex, South Africa, using heterogeneity in the basal zone of the Partridge River Intrusion: anomalously fractionated sulfur isotopes as a tracer. EOS Trans.
  • 8. Edward M. RIPLEY et al. / Earth Science Frontiers, 2007, 14(5):124 –132 AGU, 2006, 87(36), Jt. Assem. Suppl., Abs.: V43B-06. xenoliths in the Voisey’s Bay Intrusion, Labrador, Canada: [31] Ito E, White W M, Göpel C. The O, Sr, Nd, and Pb isotope mineralogy, reactions, partial melting and mechanisms of mass geochemistry of MORB. Chem. Geol., 1987, 62: 157−176. transfer. Geochemistry, Geophysics, Geosystems, 2006, Q05013. [32] Harmon R S, Hoefs J. Oxygen isotope heterogeneity of the doi10.1029/2005GC001184. 18 mantle deduced from global O systematics of basalts from [41] Faure G. Origin of Igneous Rocks: The Isotopic Evidence. different geotectonic settings. Con. Mineral. Petrol., 1995, 120: Berlin: Springer, 2001, 496. 95−114. [42] Hanson G N. An approach to trace element modeling using a [33] Eiler J M. Oxygen isotope variations in basaltic lavas and upper simple igneous system as an example. In: Lipin B R, McKay G mantle rocks. In: Valley J W, Cole D R (eds.), Stable Isotope A (eds.), Geochemistry and Mineralogy of Rare Earth Elements. Geochemistry. Min. Soc. Amer. Rev. in Min. and Geochem. Min. Soc. Amer. Rev. Min., 1989, 21: 79−98. 2001, 43: 319−364. [43] Dickin A P. Radiogenic Isotope Geology. Cambridge University [34] Leitch A M. Analog experiments on melting and contamination Press, 1995, 479. at the roof and walls of magma chambers. Jour. Volc. Geotherm. [44] Depaolo D. Trace element effects of combined wall rock Res., 2004, 129: 173−197. assimilation and fractional crystallization. Earth Planet. Sci. [35] Taylor H P Jr. The application of oxygen and hydrogen isotope Lett., 1981, 53: 189−202. studies to problems of hydrothermal alteration and ore [45] Mathez E A, Waight T. Lead isotopic disequilibrium between deposition. Econ. Geol., 1974, 69: 843−883. sulfide and plagioclase in the Bushveld Complex and them [36] Criss R E, Taylor H P Jr. Meteoric-hydrothermal systems. In: chemical evolution of large layered intrusions. Geochim. Valley J W, Taylor S P Jr, O’Neil J R (eds.), Stable Isotopes in Cosmochim. Acta, 2003, 62: 1875−1888. High Temperature Geological Processes. Min. Soc. Amer. Rev. [46] Mathez E A, Kent A J R. Variable initial Pb isotopic Min., 1986, 16: 373−424. compositions of rocks associated with the UG2 chromitite, [37] Criss R E. Principles of Stable Isotope Distribution. Oxford eastern Bushveld Complex. Geochim. Cosmochim. Acta, in University Press, 1999, 254. press. [38] O’Neil J R, Taylor H P Jr. The oxygen isotope and cation [47] Alard O, Griffin W L, Pearson N J, et al. New insights into the exchange chemistry of feldspars. Am. Mineral, 1969, 52: Re-Os systematics of subcontinental lithospheric mantle from 1414−1437. in situ analysis of sulphides. Earth Planet. Sci. Lett., 2002, 203: [39] Li C, Naldrett A J. Melting reactions of gneissic inclusions with 651−663. enclosing magma at Voisey’s Bay, Labrador: implications with [48] Griffin W L, Graham S, O’Reilly S Y. Lithosphere evolution respect to ore genesis. Econ. Geol., 2000, 95: 801−814. beneath the Kaapraal Craton: Re-Os systematics of sulfides in [40] Mariga J, Ripley E M, Li C. Petrogenetic evolution of gneissic mantle-derived peridotites. Chem. Geol., 2004, 208: 89−118.