Sulfides are among the most important petrogenetic agents in magmatic systems. They are ubiquitous in most upper-mantle rock types, common as inclusions in diamonds and they host significant amounts of geochemically and economically important chalcophile (‘sulfur-loving’) elements, such as Cu, Ni, Pb, In, Au and the platinum-group elements. Despite their low abundance (<< 1% of the bulk rock), residual sulfides have a disproportionate control over the chalcophile element budget in upper mantle lithologies, as well as that of melts derived from the Earth’s mantle.
Keywords: chalcophile, element partitioning, upper mantle, sulfides
Sulfides are minerals that contain the sulfide anion, S2–, as the major structural component. Sulfides are common accessory minerals in the Earth’s interior. They are observed in upper-mantle peridotites and oceanic basalts, in mantle xenoliths, and are also the most common mineral inclusion in diamonds. Despite sulfur being a minor component in the mantle – present at concentrations of ~150–300 parts per million (ppm) S in the bulk silicate Earth (e.g. Palme and O’Neill 2013) – sulfides are, nonetheless, very important phases because many elements have a strong affinity for sulfur. These sulfur-loving elements are termed ‘chalcophiles’ (Fig. 1), and they are concentrated in sulfide minerals. Thus, sulfides control the behaviour of the chalcophile elements during magmatic processes.
Most chalcophile elements (some of which are also siderophile, or ‘iron-loving’) are economically important metals; ten of these are the most valuable metals in the world [In, Re, Au, Ag and the platinum-group elements (PGEs)]. In Earth’s crust, these elements are typically found associated with base-metal sulfides, and sulfide deposits are mined extensively as a result. In addition to their economic importance, residual sulfides exert a strong control on the budget of chalcophile elements during partial melting and crustal differentiation and, hence, have a profound effect on the trace-element concentrations found in erupted magmas.
In this article, we explore the occurrence of sulfides in the Earth’s mantle, the role of sulfides in the partitioning of chalcophile elements in magmatic systems, and the abundance of chalcophile elements in mantle-derived silicate melts.
Sulfides in the upper mantle
Sulfides are ubiquitous accessory phases in all types of mantle xenoliths. The major sulfide phases present in mantle rocks are pentlandite, pyrrhotite and chalcopyrite (Vaughan and Corkhill 2017 this issue). Also present are the ‘monosulfide solid solution’ (mss) and ‘intermediate solid solution’ (iss) phases. The observed mineralogy of mantle sulfides, however, likely represents low temperature (<300 °C) re-equilibration of high-temperature mss or, possibly, sulfide melt. Therefore, most recent studies have concentrated on reporting the bulk chemistry of polyphase inclusions (e.g. Bulanova et al. 1996), and this will be the focus of this discussion.
The compositions of sulfides from mantle xenoliths and diamond inclusions are shown in Figure 2. Sulfides observed in mantle xenoliths are dominantly Fe-rich with minor amounts of Ni and Cu. Peridotite xenoliths generally contain sulfides with ~20 wt% Ni and 4–5 wt% Cu, on average. Eclogitic sulfides are generally similar with respect to their Cu content compared to peridotite-hosted sulfides, but are consistently lower in Ni, with average concentrations ~3 wt% Ni.
Even though sulfides are accessory phases in mantle xenoliths (<0.1%) they are the most common mineral inclusions found in diamonds, accounting for 35%–46% of all diamond inclusions in the Finsch, Koffiefontein and Premier kimberlite pipes of South Africa (Harris and Gurney 1979). The Ni and Cu concentrations in sulfide inclusions in both peridotite- and eclogite-type diamonds are broadly similar compared to intergranular sulfides in peridotite and eclogite xenoliths, respectively. The marked difference in Ni concentrations between peridotite- and eclogite-hosted sulfides has been used as a paragenetic discriminator (Bulanova et al. 1996). Pearson et al. (1998) also noted that Ni concentrations of diamond-hosted sulfide inclusions correlate broadly with Os, and that Re/Os may be a more robust petrogenetic indicator for differentiating between peridotitic and eclogitic sulfides. Of course, how the concentrations of these elements in sulfides may be used as petrogenetic indicators hinges on our understanding of their partitioning behaviour between sulfide and silicate melt during mantle melting.
Theoretical considerations on partitioning of chalcophile elements into sulfides
Because sulfides have the ability to concentrate economically important elements, several studies have been aimed at quantifying the chalcophile character of ore-forming metals. Most experimental studies have focused on PGEs (see review by Brenan et al. 2016) and, to a lesser extent, on the chalcophile elements, such as Cu, Ag, Sb, Tl, As and others (e.g. Kiseeva and Wood 2013; Li and Audétat 2012).
The most obvious way to investigate the behaviour of chalcophile elements in magmatic systems is through studying the partitioning of trace elements between sulfide and silicate liquids. The partition coefficient (D) for an element M is defined as the concentration of this element in a sulfide liquid divided by its concentration in a silicate liquid.
The higher the D, the more element M prefers to go into a sulfide phase, hence the more chalcophile it is. Determining D can be tackled using two different experimental approaches, both, as shown below, being similarly effective.
The first approach uses oxide–sulfide equilibria at controlled oxygen and sulfur fugacity conditions. The exchange of a chalcophile element, such as Pb, between sulfide and silicate melt can be described by the following reaction:
where x corresponds to the valence of Pb in silicate and sulfide melts. Assuming Pb is present in dilute concentrations, we can express Equation (1) in terms of the equilibrium constant (Brenan et al. 2016) and relate the DPbsulf / sil with oxygen and sulfur fugacity according to:
An essential conclusion to be drawn from Equation (2) is that the partition coefficient of a chalcophile element between sulfide and silicate melts is a function of the ratio between sulfur and oxygen fugacity. In other words, higher fS2 values will favour the partitioning of Pb to a sulfide melt (i.e. higher DPbsulf / sil), whereas high values of fO2 will drive Pb into silicate melt (i.e. lower DPbsulf / sil). In nature, however, because sulfur is a redox-sensitive element, fO2 and fS2 are inversely proportional, and as a result, do not vary independently.
In the more oxidising conditions prevalent during arc magmatism, sulfide is not stable and sulfate becomes the preferred S species (e.g. Jugo et al. 2010). When this occurs, fS2 will drop considerably, while fO2 can increase unabated. However, it is desirable to evaluate how fO2 and fS2 affect trace-element behaviour individually, because different trace elements will respond to changes in these parameters in different ways depending on how chalcophile they are.
A large number of experimental studies of silicate melts at 1 atm pressure were carried out during the 1980s and 1990s, in which fO2 and fS2 were controlled either by gas mixing or by the use of redox buffers (see review by Brenan et al. 2016). Unfortunately, some of these studies turned out to be compromised by contamination of their silicate glasses by the noble metals that the experiments were doped with. The so-called micronugget contaminants, which are sub-micrometer particles of metal and/or sulfide present in the silicate melt, are not dissolved species and can easily (wrongly) be included into the glass composition when experimental products are analysed. This has resulted in a widespread underestimation of the melts DMsulf / sil, in particular at more reducing conditions (cf. Ertel et al. 2008).
Only over the last couple of decades has the micronugget contamination issue been overcome by either measuring the solubilities of individual chalcophile elements in sulfide melt and silicate melt independently (e.g. Fonseca et al. 2011), or by measuring sulfide–silicate equilibria directly (Mungall and Brenan 2014). As a result, more recent studies now show DMsulf / sil increasing with decreasing fO2/fS2, as expected from Equation (2).
A good illustration of the dependence of the partitioning behaviour of a given trace element with fO2/fS2 is given by Os. Older data obtained from experiments that did not consider micronugget contamination (grey symbols in Fig. 3a) are completely at odds with more recent Os partitioning data, in that they show the opposite dependence on fO2/fS2 expected from Equation (2). More recent studies (blue symbols in Fig. 3a), where efforts were made to avoid micronugget contamination, agree with the prediction of Equation (2). This contrast is a good illustration of the strides that have been made in investigating the chalcophile nature of trace elements over the last 25 years.One other complexity that arises from the experimental approach described above is that the partition coefficient for a chalcophile element is related to sulfur and oxygen fugacities, which are easily controlled in laboratory experiments but that can be difficult to measure or estimate for natural systems. To overcome this problem, Kiseeva and Wood (2013, 2015) proposed an alternative second approach.
Following the study by Haughton et al. (1974), who related fS2 to FeS and fO2 to FeO, Kiseeva and Wood (2013) derived an expression that substitutes the fO2/fS2 ratio from Equation (2) with the FeO content of the silicate liquid (provided that the silicate melt is in equilibrium with the FeS sulfide), considering the exchange reaction:
MOx/2 (silicate) + (x/2)FeS (sulfide) = MSx/2 (sulfide) + (x/2)FeO (silicate) (3)
Again, if we use divalent Pb as an example, Equation (3) will become:
PbO (silicate) + FeS (sulfide) = PbS (sulfide) + FeO (silicate) (4)
As for Equation (1), Equation (4) can be rearranged so that the partition coefficient DPbsulf / sil relates to the equilibrium constant and the activity of FeO in silicate melt (Kiseeva and Wood 2013). This leads to the very simple statement that when a magma is sulfide-saturated (where the activity of FeS will be approximately unity), DPbsulf / sil mainly depends on the FeO content of the silicate liquid:
log DPbsulf / sil ≈ A + —x/2 log[FeO] (5)
where A is a constant related to the free energy of the Fe–Pb exchange between sulfide and silicate liquids; x is a constant related to the valence of the element; [FeO] is the FeO content of the silicate melt in mole fraction or weight %. In theory, −x will be the valence of the chalcophile element being studied (Fig. 3b).
Relationships between partition coefficients and FeO content of the silicate melt, summarised by Equation (5), have been successfully tested on a large number (>15) of chalcophile and moderately chalcophile elements and have proved to be valid for a broad range of temperatures and sulfide and silicate melt compositions. Moreover, the method outlined in Equation (5) allows the partition coefficient of a given chalcophile element to be described in terms of the easily measurable FeO content of silicate magma, precluding the need to account for oxygen and sulfur fugacities.
The two experimental approaches aimed at determining the partitioning relationships between the sulfide and silicate liquids are essentially addressing the same equilibrium, and they provide similar results for chosen T–fO2 conditions, most notably when applied to melting of Earth’s mantle.
To illustrate the equivalence of the approaches, we applied the two methods to the same dataset using data from Gaetani and Grove (1997). These authors controlled fO2 and fS2 and they synthesised silicate glasses over a range of FeO concentrations. In Figure 4, DNisulf / sil is plotted as a function of log(fO2/fS2) and log[FeO] and shows an excellent correlation.The method described in Equation (5) enables us to compile partition coefficients for a number of chalcophile elements in typical magmas (Table 1).
Chalcophile element behaviour during mantle melting
The behaviour of chalcophile elements during mantle melting is summarised in Figure 5. Elements with a very strong chalcophile affinity (DMsulf / sil of 104–106) (Table 1), such as Pt, Ir and Au, are extremely depleted in the silicate liquids due to their retention in the residual mantle sulfide. Only when sulfide is almost exhausted from the mantle source (~16% of aggregate fractional melting, assuming 140 ppm S in the mantle source), will highly chalcophile elements partition into the silicate melt. The concentration of chalcophiles that are incompatible in silicate phases, such as Pt and Au, decrease with continued melting as they become diluted by the silicate melt, while Ir, which is likely retained by olivine (Brenan et al. 2016), remains at constant concentrations.
The moderately chalcophile elements, such as Ni, Cu, In, and Pb have DMsulf / sil values of ~10–1,000. Both In and Pb are incompatible elements in mantle silicates, so their concentrations continually decrease with increased melting. Copper, which is only moderately incompatible, is relatively constant during fractional melting and shows only a small increase in concentration when approaching sulfide exhaustion. Despite being strongly chalcophile, Ni is also compatible in olivine (D ~2–6), and this dominates its partitioning. Thus, sulfide exerts little or no control on Ni concentration during the melting of peridotite.
The role of sulfides in controlling the Re–Os isotope system
Rhenium and osmium were scavenged into the Earth’s core during terrestrial accretion, leaving the silicate mantle highly depleted in these elements. Although they share this particular geochemical property, Re and Os have contrasting geochemical behaviours during partial melting of the mantle: Re is thought to be mildly incompatible, whereas Os is thought to be compatible (e.g. Fonseca et al. 2011). This fundamental difference is unique among commonly used isotope chronometers (e.g. K–Ar, Rb–Sr, Sm–Nd, Lu–Hf, and U–Th–Pb–He systems) where both parents and daughters are highly incompatible elements. The fact that Re is mildly incompatible and Os is compatible causes large variations of Re/Os ratios in oceanic basalts. Mid-ocean ridge basalts, for example, have Re/Os ratios ranging between 10 and 2,000 (Shirey and Walker 1998). This variability makes the Re–Os isotopic system a useful geochemical tool for understanding mantle evolution.
The experimental literature has provided some detailed insights into the geochemical behaviour of Re during partial melting of the mantle. Rhenium is the least chalcophile of all the highly siderophile elements (e.g. Mungall and Brenan 2014), and it has been shown to partition into other mantle phases such as clinopyroxene, orthopyroxene, garnet, and spinel (Mallmann and O’Neill 2007), particularly under more reducing conditions. At oxidising conditions, however, Re will become hexavalent and will behave as a highly incompatible element during partial melting of the Earth’s mantle. Osmium, on the other hand, is strongly compatible in sulfide (e.g. Mungall and Brenan 2014): about 97% of the Os budget in source peridotites is controlled by sulfide (Luguet et al. 2007).
Effectively then, residual sulfide is the primary phase in the Earth’s mantle that confers on the Re–Os isotope system its unique ability to date episodes of mantle melting. Mantle sulfides efficiently retain the entirety of the Os budget of a peridotite during partial melting, while Re, owing to its high incompatibility, is efficiently removed from the mantle. At high degrees of partial melting (10%–20%), alloys are expected to form from a host sulfide (Luguet et al. 2007; Fonseca et al. 2011), and, because of this, Os will further fractionate from Re.
If melt-depleted peridotites can be presumed to have lost all their Re during partial melting, then their 187Os/188Os can potentially be used to derive model ages that ostensibly track melt extraction events (reviewed in Lorand and Luguet 2016). This approach can be applied at both the whole-rock and mineral (i.e. sulfide) scale to date ancient melting events in the Earth’s mantle (e.g. Wainwright et al. 2015). However, the same thing that makes the Re/Os system so useful in dating ancient melt events is also responsible for its biggest caveat. For example, metasomatic sulfides, which can be introduced to a mantle peridotite through basaltic melt percolation and/or fluid infiltration, will typically introduce Re into a virtually Re-free peridotite. Because of this, these metasomatic sulfides will develop more radiogenic 187Os/188Os over time, which will mix into the whole-rock signature and obscure the actual age of melt depletion of that peridotite. This difficulty can be partly circumvented by carrying out studies at the mineral scale to identify different generations of mantle sulfide, each with their own 187Os/188Os (cf. Wainwright et al. 2015). Sulfides are, thus, unique in their ability to fractionate Re from Os in magmatic systems and to control the Os isotope systematics during melting of the Earth’s mantle.
Concluding thoughts and future directions
The role of sulfides in igneous processes is not limited to the examples listed above. For instance, the Cu/Ag ratios of erupted magmas are fully controlled by the distinctive partitioning of Cu and Ag into solid and liquid sulfides (Wang and Becker 2015). Furthermore, the precipitation of sulfide may be responsible for the decline in Cu concentration of arc lavas during the evolution of volcanic systems at convergent margins (Lee et al. 2012). Sulfide also plays a significant role in the mysterious behaviour of Pb (Hart and Gaetani 2006) during mantle melting, which includes the nearly constant Ce/Pb (~25) and Nd/Pb (~20) ratios of mid-ocean ridge basalts and of ocean island basalts and the ‘lead paradox’, whereby all major Pb reservoirs appear to be substantially younger than the silicate Earth.
The Pd/Cu ratio may be an indicator of sulfide saturation in magmatic systems (Park et al. 2013), a parameter of fundamental importance in volcanic degassing. Ultimately, the precipitation of sulfides that act as a reservoir for the chalcophile and highly siderophile elements, such as Cu, Ag, Pb and the PGEs, is one of the key triggers that eventually leads to the formation of porphyry ore deposits and are key for the formation of magmatic sulfide deposits.
Thus, given that a typical peridotite will have sulfide phases present at less than 600 ppm (or 220 ppm S), such phases could certainly be said to “punch well above their weight”.
The authors thank Frances Jenner, Yuan Li and Jim Mungall for thorough reviews. We also thank Hugh O’Neill and Bernie Wood for comments and discussion. ESK acknowledges support of the NERC grant NE/L010828/1.
Brenan JM (2008) Re–Os fractionation by sulfide melt–silicate melt partitioning: a new spin. Chemical Geology 248: 140-165
Brenan JM, Bennett NR, Zajacz Z (2016) Experimental results on fractionation of the highly siderophile elements (HSE) at variable pressures and temperatures during planetary and magmatic differentiation. Reviews in Mineralogy and Geochemistry 81: 1-87
Bulanova GP, Griffin WL, Ryan CG, Shestakova OY, Barnes, S-J (1996) Trace elements in sulfide inclusions from Yakutian diamonds. Contributions to Mineralogy and Petrology 124: 111-125
Efimova ES, Sobolev NV, Pospelova LN (1983) Sulphide inclusions in diamonds and peculiarities of their paragenesis (in Russian). Zapiski Vsesoyuznogo Mineralogicheskogo Obshchestva 92: 300-309
Ertel W, Dingwell DB, Sylvester PJ (2008) Siderophile elements in silicate melts — A review of the mechanically assisted equilibration technique and the nanonugget issue. Chemical Geology 248: 119-139
Fonseca ROC, Mallmann G, O’Neill HStC, Campbell IH, Laurenz V (2011) Solubility of Os and Ir in sulfide melt: implications for Re/Os fractionation during mantle melting. Earth and Planetary Science Letters 311: 339-350
Gaetani GA, Grove TL, (1997) Partitioning of moderately siderophile elements among olivine, silicate melt, and sulfide melt: constraints on core formation in the Earth and Mars. Geochimica et Cosmochimica Acta 61: 1829-1846
Ghiorso MS, Hirschmann MM, Reiners PW, and Kress VC (III) (2002) The pMELTS: A revision of MELTS for improving calculation of phase relations and major element partitioning related to partial melting of the mantle to 3 GPa. Geochemistry, Geophysics, Geosystems 3: 1-35
Guo J, Griffin WL, O’Reilly SY (1999) Geochemistry and origin of sulphide minerals in mantle xenoliths: Qilin, southeastern China. Journal of Petrology 40: 1125-1149
Harris JW, Gurney JJ (1979) Inclusions in diamond. In: Field JR (ed) The Properties of Diamond. Academic Press, London pp 555-591
Hart SR, Gaetani GA (2006) Mantle Pb paradoxes: the sulfide solution. Contributions to Mineralogy and Petrology 152: 295-308
Haughton DR, Roeder PL, Skinner BJ (1974) Solubility of sulfur in mafic magmas. Economic Geology 69: 451-467
Jugo PJ, Wilke M, Botcharnikov RE (2010) Sulfur K-edge XANES analysis of natural and synthetic basaltic glasses: implications for S speciation and S content as function of oxygen fugacity. Geochimica et Cosmochimica Acta 74: 5926-5938
Kiseeva ES, Wood BJ (2013) A simple model for chalcophile element partitioning between sulphide and silicate liquids with geochemical applications. Earth and Planetary Science Letters 383: 68-81
Kiseeva ES, Wood BJ (2015) The effects of composition and temperature on chalcophile and lithophile element partitioning into magmatic sulphides. Earth and Planetary Science Letters 424: 280-294
Lee C-TA and 8 coauthors (2012) Copper systematics in arc magmas and implications for crust-mantle differentiation. Science 336: 64-68
Li Y, Audétat A (2012) Partitioning of V, Mn, Co, Ni, Cu, Zn, As, Mo, Ag, Sn, Sb, W, Au, Pb, and Bi between sulfide phases and hydrous basanite melt at upper mantle conditions. Earth and Planetary Science Letters 355–356: 327-340
Lorand J-P, Luguet A (2016) Chalcophile and siderophile elements in mantle rocks: trace elements controlled by trace minerals. Reviews in Mineralogy & Geochemistry 81: 441-488
Luguet A, Shirey SB, Lorand J-P, Horan MF, Carlson RW (2007) Residual platinum-group minerals from highly depleted harzburgites of the Lherz massif (France) and their role in HSE fractionation of the mantle. Geochimica et Cosmochimica Acta 71: 3082-3097
Mallmann G, O’Neill HStC (2007) The effect of oxygen fugacity on the partitioning of Re between crystals and silicate melt during mantle melting. Geochimica et Cosmochimica Acta 71: 2837-2857
Mungall JE, Brenan JM (2014) Partitioning of platinum-group elements and Au between sulfide liquid and basalt and the origins of mantle-crust fractionation of the chalcophile elements. Geochimica et Cosmochimica Acta 125: 265-289
Palme H, O’Neill HStC (2013). Cosmochemical estimates of mantle composition. In: Carlson RW (ed) Treatise on Geochemistry, Volume 3 (2nd edition), Elsevier, Oxford, pp 1-39
Pearson DG, Shirey SB, Harris JW, Carlson RW (1998) Sulphide inclusions in diamonds from the Koffiefontein kimberlite, S Africa: constraints on diamond ages and mantle Re–Os systematics. Earth and Planetary Science Letters 160: 311-326
Richardson SH, Shirey SB, Harris JW, Carlson RW (2001) Archean subduction recorded by Re–Os isotopes in eclogitic sulfide inclusions in Kimberley diamonds. Earth and Planetary Science Letters 191: 257-266
Shirey SB, Walker RJ (1998) The Re-Os isotope system in cosmochemistry and high-temperature geochemistry. Annual Review of Earth and Planetary Sciences 26: 423-500
Smythe DJ, Wood BJ, Kiseeva ES (2017) The S content of silicate melts at sulfide saturation: new experiments and a model incorporating the effects of sulfide composition. American Mineralogist: 102: 795-803
Vaughan DJ, Corkhill C (2017) Mineralogy of sulphides. Elements 13: 83-89
Wainwright AN, Luguet A, Fonseca ROC, Pearson DG (2015) Investigating metasomatic effects on the 187Os isotopic signature: a case study on micrometric base metal sulphides in metasomatised peridotite from the Letlhakane kimberlite (Botswana). Lithos 232: 35-48
Wang Z, Becker H (2015). Abundances of Ag and Cu in mantle peridotites and the implications for the behavior of chalcophile elements in the mantle. Geochimica et Cosmochimica Acta 160: 209-226
Westerlund KJ and 5 coauthors (2006) A subduction wedge origin for Paleoarchean peridotitic diamonds and harzburgites from the Panda kimberlite, Slave craton: evidence from Re–Os isotope systematics. Contributions to Mineralogy and Petrology 152: 275-294