The sulfur isotope evolution of magmatic-hydrothermal fluids: insights into ore-forming processes

Metal-rich fluids that circulate in magmatic-hydrothermal environments form a wide array of economically significant ore deposits. Unravelling the origins and evolution of these fluids is crucial for understanding how Earth’s metal resources form and one of the most widely used tools for tracking these processes is sulfur isotopes. It is well established that S isotopes record valuable information about the source of the fluid, as well as its physical and chemical evolution (i.e. changing pH, redox and temperature), but it is often challenging to unravel which of these competing processes drives isotopic variability. Here we use thermodynamic models to predict S isotope fractionation for geologically realistic hydrothermal fluids and attempt to disentangle the effects of fluid sources, physico-chemical evolution and S mineral disequilibrium. By modelling a range of fluid compositions, we show that S isotope fingerprints are controlled by the ratio of oxidised to reduced S species (SO42−/H2S), and this is most strongly affected by changing temperature, fO2 and pH. We show that SO42−/H2S can change dramatically during cooling and our key insight is that S isotopes of individual sulfide or sulfate minerals can show large fractionations (up to 20‰) even when pH is constant and fO2 fixed to a specific mineral redox buffer. Importantly, while it is commonly assumed that SO42−/H2S is constant throughout fluid evolution, our analysis shows that this is unlikely to hold for most natural systems. We then compare our model predictions to S isotope data from porphyry and epithermal deposits, seafloor hydrothermal vents and alkaline igneous bodies. We find that our models accurately reproduce the S isotope evolution of porphyry and high sulfidation epithermal fluids, and that most require magmatic S sources between 0 and 5‰. The S isotopes of low sulfidation epithermal fluids and seafloor hydrothermal vents do not fit our model predictions and reflect disequilibrium between the reduced and oxidised S species and, for the latter, significant S input from seawater and biogenic sources. Alkaline igneous fluids match model predictions and confirm magmatic S sources and a wide range of temperature and redox conditions. Of all these different ore deposits, porphyry and alkaline igneous systems are particularly well-suited to S isotope investigation because they show relationships between redox, alteration and ore mineralogy that could be useful for exploration and prospecting. Ultimately, our examples demonstrate that S isotope forward models are powerful tools for identifying S sources, flagging disequilibrium processes, and validating hypotheses of magmatic fluid evolution.