The oxidation state of sulfur in synthetic and natural glasses determined by X-ray absorption spectroscopy

Sulfur K-edge X-ray absorption near edge structure (XANES) spectra were recorded for experimental glasses of various compositions prepared at different oxygen fugacities (fO₂) in one-atmosphere gas-mixing experiments at 1400 °C. This sample preparation method only results in measurable S concentrations under either relatively reduced (log fO₂ < −9) or oxidised (log fO₂ > −2) conditions. The XANES spectra of the reduced samples are characterised by an absorption edge crest at 2476.4 eV, typical of S²⁻. In addition, spectra of Fe-bearing compositions exhibit a pronounced absorption edge shoulder. Spectra for all the Fe-free samples are essentially identical, as are the spectra for the Fe-bearing compositions, despite significant compositional variability within each group. The presence of a sulfide phase, such as might exsolve on cooling, can be inferred from a pre-edge feature at 2470.5 eV.The XANES spectra of the oxidised samples are characterised by an intense transition at 2482.1 eV, typical of the sulfate anion SO₄²⁻. Sulfite (SO₃²⁻) has negligible solubility in silicate melts at low pressures. The previous identification of sulfite species in natural glass samples is attributed to an artefact of the analysis (photoreduction of S⁶⁺). S⁴⁺ does, however, occur unambiguously with S⁶⁺ in Fe-free and Fe-poor compositions prepared in equilibrium with CaSO₄ at 4-16 kbar, and when buffered with Re/ReO₂ at 10 kbar. Solubility of S⁴⁺ thus requires partial pressures of SO₂ considerably in excess of 1 bar. A number of experiments were undertaken in an attempt to access intermediate fO₂s more applicable to terrestrial volcanism. Although these were largely unsuccessful, S²⁻ and S⁶⁺ were found to coexist in some samples that were not in equilibrium with the imposed fO₂.The XANES spectra of natural olivine-hosted melt inclusions and submarine glasses representative of basalts at, or close to, sulfide saturation show mainly dissolved S²⁻, but with minor sulfate, and additionally a peak at 2469.5 eV, which, although presumably due to immiscible sulfide, is 1 eV lower than that typical of FeS. These sulfate and sulfide-related peaks disappear with homogenisation of the inclusions by heating to 1200 °C followed by rapid quenching, suggesting that both these features are a result of cooling under natural conditions. The presence of small amounts of sulfate in otherwise reduced basaltic magmas may be explained by the electron exchange reaction: S²⁻ + 8Fe³⁺ = S⁶⁺ + 8Fe²⁺, which is expected to proceed strongly to the right with decreasing temperature. This reaction would explain why S²⁻ and S⁶⁺ are frequently found together despite the very limited fO₂ range over which they are thermodynamically predicted to coexist. The S XANES spectra of water-rich, highly oxidised, basaltic inclusions hosted in olivine from Etna and Stromboli confirm that nearly all S is dissolved as sulfate, explaining their relatively high S contents.     

Gold solubility in oxidized and reduced, water-saturated mafic melt

We have performed experiments to evaluate Au solubility in natural, water-saturated basaltic melts as a function of oxygen fugacity. Experiments were carried out at 1000 °C and 200 MPa, and oxygen fugacity was controlled at the fayalite-magnetite-quartz (FMQ) oxygen fugacity buffer and FMQ + 4. All experiments were saturated with a metal-chloride aqueous solution loaded initially as a 10 wt% NaCl eq. fluid. The stable phase assemblage at FMQ consists of basalt melt, olivine, clinopyroxene, a single-phase aqueous fluid, and metallic Au. The stable phase assemblage at FMQ + 4 consists of basalt melt, clinopyroxene, magnetite-spinel solid solution, a single-phase aqueous fluid, and metallic Au. Silicate glasses (i.e., quenched melt) and their contained crystalline material were analyzed by using both electron probe microanalysis (EPMA) and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Measured Au concentrations in the quenched melt range from 4.8 μg g⁻¹ to 0.64 μg g⁻¹ at FMQ + 4, and 0.54 μg g⁻¹ to 0.1 μg g⁻¹ at FMQ. The measured solubility of Au in olivine and clinopyroxene was consistently below the LA-ICP-MS limit of detection (i.e., 0.1 μg g⁻¹). These melt solubility data place important limitations on the dissolved Au content of water-saturated, Cl- and S-bearing basaltic liquids at geologically relevant fO₂ values. The new data are compared to published, experimentally-determined values for Au solubility in dry and hydrous silicate liquids spanning the compositional range from basalt to rhyolite, and the effects of melt composition, oxygen fugacity, pressure and temperature are discussed.     

Redox evolution of differentiating hydrous basaltic magmas recorded by zircon and apatites in mafic cumulates: The case of the Malayer Plutonic Complex,Western Iran

Mafic-ultramafic cumulates can provide records of basaltic magma chambers' conditions and processes, which are often difficult to determine in areas dominated by crustal-derived felsic intrusions, such as the Malayer Plutonic Complex (MPC), Western Iran. New U-Pb zircon ages for mafic cumulates in the MPC confirm the presence of isolated magma chambers of contrasting compositions during Middle Jurassic. Mafic cumulates found in seven separate zones across the MPC vary from olivine gabbro to anorthosite. While the mineralogical, textural, and geochemical lines of evidence recorded in mafic cumulates indicate pH₂O controls on the liquidus phases, the estimated oxygen fugacity (logfO₂) using zircon and apatite chemistry suggests a smoothly rising redox state during the fractionation process, consistent with the trend expected for late-stages differentiation of hydrous arc magmas. This trend is further confirmed by sulfur speciation in apatites determined from microbeam sulfur K-edge X-ray absorption near edge structure (μ-XANES) spectra (S⁶⁺/∑S = 0.93–0.98 ~ FMQ + 2 to 0.99 ~ FMQ + 3, where ∑S = S⁶⁺+S⁴⁺+S²⁻). The low S content and increasing redox state of the fractionating basaltic melts most likely resulted from preferential removal of sulfur en-route to the magma chambers along with effective assimilation of oxidizing crustal components. The reduced condition in the early basaltic melt is also evidenced by the presence of pyrite and magnetite inclusions in olivines in mafic cumulates. The shift in the prevailing fO₂ from sulfide-saturated to sulfate-bearing recorded by MPC mafic cumulates, similar to that in other magmatic arcs, is accompanied by changes in the differentiation path from transitional tholeiitic to calc-alkaline.     

The effect of water activity on the oxidation and structural state of Fe in a ferro-basaltic melt

Experimental investigations have been performed at T = 1200°C, P = 200 MPa and fH₂ corresponding to H₂O-MnO-Mn₃O₄ and H₂O-QFM redox buffers to study the effect of H₂O activity on the oxidation and structural state of Fe in an iron-rich basaltic melt. The analysis of Mössbauer and Fe K-edge X-ray absorption nearedge structure (XANES) spectra of the quenched hydrous ferrobasaltic glasses shows that the Fe³⁺/ΣFe ratio of the glass is directly related to aH₂O in a H₂-buffered system and, consequently, to the prevailing oxygen fugacity (through the reaction of water dissociation H₂O ↔ H₂ + 1/2 O₂). However, water as a chemical component of the silicate melt has an indistinguishable effect on the redox state of iron at studied conditions. The experimentally obtained relationship between fO₂ and Fe³⁺/Fe²⁺ in the hydrous ferrobasaltic melt can be adequately predicted in the investigated range by the existing empiric and thermodynamic models. The ratio of ferric and ferrous Fe is proportional to the oxygen fugacity to the power of ∼0.25 which agrees with the theoretical value from the stoichiometry of the Fe redox reaction (FeO + ¼ O₂ = FeO_1.5). The mean centre shifts for Fe²⁺ and Fe³⁺ absorption doublets in Mössbauer spectra show little change with increasing Fe³⁺/ΣFe, suggesting no significant change in the type of iron coordination. Similarly, XANES preedge spectra indicate a mixed (C3h, Td, and Oh, i.e., 5-, 4-, and sixfold) coordination of Fe in hydrous basaltic glasses.     

Effects of fO2, fS2, temperature, and melt composition on Fe-Ni exchange between olivine and sulfide liquid: implications for natural olivine-sulfide assemblages

The apparent equilibrium constant for the exchange of Fe and Ni between coexisting olivine and sulfide liquid (K_D = (X_NiS/X_FeS)_liquid/(X_NiSi12O2/X_FeSi12O2)_olivine; X_i = mole fraction) has been measured at controlled oxygen and sulfur fugacities (fO₂ = 10^−8.1 to 10⁻¹⁰ and fS₂ = 10^−0.9 to 10^−1.7) over the temperature range 1200 to 1385°C, with 5 to 37 wt% Ni and 7 to 18 wt% Cu in the sulfide liquid. At log fO₂ of −8.7 ± 0.1, and log fS₂ of −0.9 to −1.7, K_D is relatively insensitive to sulfur fugacity, but comparison with previous results shows that K_D increases at very low sulfur fugacities. K_D values show an increase with the nickel content of the sulfide liquid, but this effect is more complex than found previously, and is greatest at log fO₂ of −8.1, lessens with decreasing fO₂, and K_D becomes independent of melt Ni content at log fO₂ ≤ −9.5. The origin of this variation in K_D with fO₂ and fS₂ is most likely the result of nonideal mixing of Fe and Ni species in the sulfide liquid. Such behavior causes activity coefficients to change with either melt oxygen content or metal/sulfur ratio, effects that are well documented for metal-rich sulfide melts.Application of these experimental results to natural samples shows that the relatively large dispersion that exists in K_D values from different olivine + sulfide-saturated rock suites can be interpreted as arising from variations in fO₂, fS₂, and the nickel content of the sulfide liquid. Estimates of fO₂ based on K_D and sulfide melt composition in natural samples yields a range from fayalite-magnetite-quartz (FMQ)-1 to FMQ-2 or lower, which is in good agreement with previous values determined for oceanic basalts that use glass ferric/ferrous ratios. Anomalously high K_D values recorded in some suites, such as Disko Island, probably reflect low fS₂ during sulfide saturation, which is consistent with indications of low fO₂ for those samples. It is concluded that the variation in K_D values from natural samples reflects olivine-sulfide melt equilibrium at conditions within the T-fO₂-fS₂ range of terrestrial mafic magmas.     

The structural behavior of ferric and ferrous iron in aluminosilicate glass near meta-aluminosilicate joins

Iron-57 resonant absorption Mössbauer spectroscopy was used to describe the redox relations and structural roles of Fe³⁺ and Fe²⁺ in meta-aluminosilicate glasses. Melts were formed at 1500 °C in equilibrium with air and quenched to glass in liquid H₂O with quenching rates exceeding 200 °C/s. The aluminosilicate compositions were NaAlSi₂O₆, Ca_0.5AlSi₂O₆, and Mg_0.5AlSi₂O₆. Iron oxide was added in the form of Fe₂O₃, NaFeO₂, CaFe₂O₄, and MgFe₂O₄ with total iron oxide content in the range ∼0.9 to ∼5.6 mol% as Fe₂O₃. The Mössbauer spectra, which were deconvoluted by assuming Gaussian distributions of the hyperfine field, are consistent with one absorption doublet of Fe²⁺ and one of Fe³⁺. From the area ratios of the Fe²⁺ and Fe³⁺ absorption doublets, with corrections for differences in recoil-fractions of Fe³⁺ and Fe²⁺, the Fe³⁺/ΣFe is positively correlated with increasing total iron content and with decreasing ionization potential of the alkali and alkaline earth cation. There is a distribution of hyperfine parameters from the Mössbauer spectra of these glasses. The maximum in the isomer shift distribution function of Fe³⁺, δ_Fe³⁺, ranges from about 0.25 to 0.49 mm/s (at 298 K relative to Fe metal) with the quadrupole splitting maximum, Δ_Fe³⁺, ranging from ∼1.2 to ∼1.6 mm/s. Both δ_Fe³⁺ and δ_Fe²⁺ are negatively correlated with total iron oxide content and Fe³⁺/ΣFe. The dominant oxygen coordination number Fe³⁺ changes from 4 to 6 with decreasing Fe³⁺/ΣFe. The distortion of the Fe³⁺-O polyhedra of the quenched melts (glasses) decreases as the Fe³⁺/ΣFe increases. These polyhedra do, however, coexist with lesser proportions of polyhedra with different oxygen coordination numbers. The δ_Fe²⁺ and Δ_Fe²⁺ distribution maxima at 298 K range from ∼0.95 to 1.15 mm/s and 1.9 to 2.0 mm/s, respectively, and decrease with increasing Fe³⁺/ΣFe. We suggest that these hyperfine parameter values for the most part are more consistent with Fe²⁺ in a range of coordination states from 4- to 6-fold. The lower δ_Fe²⁺-values for the most oxidized melts are consistent with a larger proportion of Fe²⁺ in 4-fold coordination compared with more reduced glasses and melts.     

Experimental data on the speciation of sulfur as a function of oxygen fugacity in basaltic melts

The speciation of sulfur as a function of oxygen fugacity was calculated in glasses of basaltic composition saturated experimentally with either sulfide or sulfate phases. The experiments were conducted on mixtures of synthetic and natural materials equilibrated at 1300 °C and 1 GPa in a piston-cylinder apparatus. Sulfur speciation was calculated by measuring the peak shift of the sulfur Kα radiation relative to a sulfide standard, whereas oxygen fugacity was calculated from the composition of olivine and spinel present in the assemblages. The results are consistent with sulfur being present as sulfite (S⁴⁺) in addition to sulfate (S⁶⁺) in oxidized melts. Therefore, sulfur speciation derived from SKα peak shifts should be seen as ”sulfate mole fraction equivalents“ (X(S⁶⁺)_eq.). Using the data available, an empiric function:

X(S6+)_eq.=0.86/(1+exp(2.89−2.23ΔFMQ))     

Relationship between Se/S and sulfur isotope ratios of hydrothermal sulfide minerals

The pH and f_{O 2} dependences of the [Se^2–]/[S^2–] ratio in chloride solutions at 100°, 200° and 300°C are predicted thermodynamically. Under the high f_{O 2} conditions where sulfate species are dominant in solution, the [Se^2–]/[S^2–] ratio always increases with increasing pH and/or f_{O 2}. Under the low f_{O 2} conditions where sulfide species are dominant in solution, the pH and f_{O 2} dependences of the [Se^2–]/[S^2–] ratio are seriously affected by the presence of native selenium. With native selenium present, the [Se^2–]/[S^2–] ratio decreases with increasing f_{O 2}, but almost independent of pH in geologically important pH regions. When native selenium is absent, the [Se^2–]/[S^2–] ratio is solely a function of pH and independent of f_{O 2}. Combining the above with the pH and f_{O 2} dependences of ³⁴S value of aqueous sulfur species, we discuss the possible influences of the pH and f_{O 2} of ore-forming solutions on the relationship between the Se/S ratio and ³⁴S value of hydrothermal sulfide minerals. The results are applied to some Japanese sulfide ore deposits.     

Vanadium K edge XANES of synthetic and natural basaltic glasses and application to microscale oxygen barometry

A new microscale oxybarometer for solar system basaltic glasses, based on vanadium K-edge X-ray absorption near edge structure (XANES) spectroscopy, is described. Vanadium is unique among abundant elements in siliceous materials in that it can potentially occur in nature in four valence states: V²⁺, V³⁺, V⁴⁺ and V⁵⁺. Consequently, the vanadium redox system is a robust oxybarometer covering at least six orders of magnitude in buffer-relative oxygen fugacity. The method was calibrated using synthetic glass standards produced under known fO₂ and temperature conditions. Correction for temperature differences among standards and unknowns was quantified using microXANES data for isobaric synthetic glass couples. Application of the method to lunar, martian, and terrestrial glasses yielded fO₂ estimates from 1.6 log units more reduced than the iron-wüstite (IW) buffer (IW-1.6) for lunar glasses, to IW + 4.0 for terrestrial glass inclusions. The martian and terrestrial results are in good agreement with previous estimates by other methods. The inferred fO₂ values for lunar pyroclastic glasses are ∼0.5 log unit more reduced than previous estimates, but the differences are comparable to analytical uncertainties. Micro-extended X-ray absorption fine structure spectra were consistent with the valence states determined by microXANES and provided additional constraints on vanadium site geometry. These results demonstrate the value of this new oxybarometer, which can be applied nondestructively to individual grains in conventional thin sections with ∼ micrometer resolution and ∼100 ppm elemental sensitivity.     

Entropy and structure of oxidized and reduced iron-bearing silicate glasses

The influence of ferrous and ferric iron on the low-temperature heat capacity and vibrational entropy of silicate glasses has been determined by adiabatic calorimetry. Two pairs of samples based on sodium disilicate and calcium Tschermak molecule compositions have been studied. Along with previous data for another Fe-bearing glass, these results have been used to complement the available set of composition independent partial molar relative entropies of oxides in silicate glasses with S₂₉₈ − S₀ values of 56.7 and 116 J/mol for FeO and Fe₂O₃, respectively. The calorimetric data indicate that the fraction of fivefold coordinated Al is significant in the CaO-“FeO”-Al₂O₃-SiO₂ system and that association of Ca²⁺ and Na⁺ with Fe³⁺ in tetrahedral coordination for charge compensation does not entail significant changes in coordination for these two cations. At very low temperatures, however, the heat capacity is no longer an additive function of composition because of unexpectedly high positive deviations from Debye laws. These anomalies are stronger for the reduced than the oxidized glasses and considerably larger than for iron-free glasses, but their origin cannot be established from the present measurements.     

Temperature dependence of sulfide and sulfate solubility in olivine-saturated basaltic magmas

The sulfur concentration at pyrrhotite- and anhydrite-saturation in primitive hydrous basaltic melt of the 2001-2002 eruption of Mt. Etna was determined at 200 MPa, T = 1050-1250 °C and at log fO₂ from FMQ to FMQ+2.2 (FMQ is Fayalite-Magnetite-Quartz oxygen buffer). At 1050 °C Au sample containers were used. A double-capsule technique, using a single crystal olivine sample container closed with an olivine piston, embedded in a sealed Au₈₀Pd₂₀ capsule, was developed to perform experiments in S-bearing hydrous basaltic systems at T > 1050 °C. Pyrrhotite is found to be a stable phase coexisting with melt at FMQ-FMQ+0.3, whereas anhydrite is stable at FMQ+1.4-FMQ+2.2. The S concentration in the melt increases almost linearly from 0.12 ± 0.01 to 0.39 ± 0.02 wt.% S at FeS-saturation and from 0.74 ± 0.01 to 1.08 ± 0.04 wt.% S at anhydrite-saturation with T ranging from 1050-1250 °C. The relationships between S concentration at pyrrhotite and/or anhydrite saturation, MgO content of the olivine-saturated melt, T, and log fO₂ observed in this study and from previous data are used to develop an empirical model for estimating the magmatic T and fO₂ from the S and MgO concentrations of H₂O-bearing olivine-saturated basaltic melts. The model can also be used to determine maximum S concentrations, if fO₂ and MgO content of the melt are known. The application of the model to compositions of melt inclusions in olivines from Mt. Etna indicates that the most primitive magmas trapped in inclusions might have been stored at log fO₂ slightly higher than FMQ+1 and at T = 1100-1150 °C, whereas more evolved melts could have been trapped at T ? 1100 °C. These values are in a good agreement with the estimates obtained by other independent methods reported in the literature.     

Martian regolith in Elephant Moraine 79001 shock melts? Evidence from major element composition and sulfur speciation

Shock veins and melt pockets in Lithology A of Martian meteorite Elephant Moraine (EETA) 79001 have been investigated using electron microprobe (EM) analysis, petrography and X-ray Absorption Near Edge Structure (XANES) spectroscopy to determine elemental abundances and sulfur speciation (S²⁻ versus S⁶⁺). The results constrain the materials that melted to form the shock glasses and identify the source of their high sulfur abundances. The XANES spectra for EETA79001 glasses show a sharp peak at 2.471 keV characteristic of crystalline sulfides and a broad peak centered at 2.477 keV similar to that obtained for sulfide-saturated glass standards analyzed in this study. Sulfate peaks at 2.482 keV were not observed. Bulk compositions of EETA79001 shock melts were estimated by averaging defocused EM analyses. Vein and melt pocket glasses are enriched in Al, Ca, Na and S, and depleted in Fe, Mg and Cr compared to the whole rock. Petrographic observations show preferential melting and mobilization of plagioclase and pyrrhotite associated with melt pocket and vein margins, contributing to the enrichments. Estimates of shock melt bulk compositions obtained from glass analyses are biased towards Fe- and Mg- depletions because, in general, basaltic melts produced from groundmass minerals (plagioclase and clinopyroxene) will quench to a glass, whereas ultramafic melts produced from olivine and low-Ca pyroxene megacrysts crystallize during the quench. We also note that the bulk composition of the shock melt pocket cannot be determined from the average composition of the glass but must also include the crystals that grew from the melt - pyroxene (En_72-75Fs_20-21Wo_5-7) and olivine (Fo_75-80). Reconstruction of glass + crystal analyses gives a bulk composition for the melt pocket that approaches that of lithology A of the meteorite, reflecting bulk melting of everything except xenolith chromite.Our results show that EETA79001 shock veins and melt pockets represent local mineral melts formed by shock impedance contrasts, which can account for the observed compositional anomalies compared to the whole rock sample. The observation that melts produced during shock commonly deviate from the bulk composition of the host rock has been well documented from chondrites, rocks from terrestrial impact structures and other Martian meteorites. The bulk composition of shock melts reflects the proportions of minerals melted; large melt pockets encompass more minerals and approach the whole rock whereas small melt pockets and thin veins reflect local mineralogy. In the latter, the modal abundance of sulfide globules may reach up to 15 vol%. We conclude the shock melt pockets in EETA79001 lithology A contain no significant proportion of Martian regolith.     

Partitioning of Se, As, Sb, Te and Bi between monosulfide solid solution and sulfide melt - Application to magmatic sulfide deposits

The chalcogenes (S, Se, Te), semimetals (As, Sb) and the metal Bi are important ligands for noble metals and form a wide range of compositionally diverse minerals with the platinum-group elements (PGE). With the exception of S, few experimental data exist to quantify the behavior of these elements in magmatic sulfide systems. Here we report experimental partition coefficients for Se, Te, As, Sb, and Bi between monosulfide solid solution (mss) and sulfide melt, determined at 950 °C at a range of sulfur fugacities (fS2) bracketed by the Fe-FeS (metal-troilite) and the Fe_1−×S-S_x (mss-sulfur) equilibria. Selenium is shown to partition in mss-saturated sulfide melt as an anion replacing S²⁻. Arsenic changes its oxidation state with fS₂ from predominantly anionic speciation at low fS₂, to cationic speciation at high fS₂. The elements Sb, Te, and Bi are so highly incompatible with mss that they can only be present in sulfide melt as cations and/or as neutral metallic species. The partition coefficients derived fall with increasing atomic radius of the element. They also reflect the positions of the respective elements in the Periodic Table: within a group (e.g., As, Sb, Bi) the partition coefficients fall with increasing atomic radius, and within a period the elements of the 15th group are more incompatible with mss than the neighboring elements of the 16th group.     

Controls on gold solubility in arc magmas: An experimental study at 1000 °C and 4 kbar

In order to (1) explain the worldwide association between epithermal gold-copper-molybdenum deposits and arc magmas and (2) test the hypothesis that adakitic magmas would be Au-specialized, we have determined the solubility of Au at 4 kbar and 1000 °C for three intermediate magmas (two adakites and one calc-alkaline composition) from the Philippines. The experiments were performed over a fO₂ range corresponding to reducing (∼NNO−1), moderately oxidizing (∼NNO+1.5) and strongly oxidizing (∼NNO+3) conditions as measured by solid Ni-Pd-O sensors. They were carried out in gold containers, the latter serving also as the source of gold, in presence of variable amounts of H₂O and, in a few additional experiments, of S. Concentrations of Au in glasses were determined by LA-ICPMS. Gold solubility in melt is very low (30-240 ppb) but increases with fO₂ in a way consistent with the dissolution of gold as both Au¹⁺ and Au³⁺ species. In the S-bearing experiments performed at ∼NNO−1, gold solubility reaches much higher values, from ∼1200 to 4300 ppb, and seems to correlate with melt S content. No systematic difference in gold solubility is observed between the adakitic and the non-adakitic compositions investigated. Oxygen fugacity and the sulfur concentration in melt are the main parameters controlling the incorporation and concentration of gold in magmas. Certain adakitic and non-adakitic magmas have high fO₂ and magmatic S concentrations favorable to the incorporation and transport of gold. Therefore, the cause of a particular association between some arc magmas and Au-Cu-Mo deposits needs to be searched in the origin of those specialized magmas by involvement of Au- and S-rich protoliths. The subducted slab, which contains metal-rich massive sulfides, may constitute a potentially favorable protolith for the genesis of magmas specialized with respect to gold.     

Interaction between sulphide and H2O in silicate melts

Reaction between dissolved water and sulphide was experimentally investigated in soda-lime-silicate (NCS) and sodium trisilicate (NS3) melts at temperatures from 1000 to 1200 °C and pressures of 100 or 200 MPa in internally heated gas pressure vessels. Diffusion couple experiments were conducted at water-undersaturated conditions with one half of the couple being doped with sulphide (added as FeS or Na₂S; 1500-2000 ppm S by weight) and the other with H₂O (∼3.0 wt.%). Additionally, two experiments were performed using a dry NCS glass cylinder and a free H₂O fluid. Here, the melt was water-saturated at least at the melt/fluid interface. Profiling by electron microprobe (sulphur) and infrared microscopy (H₂O) demonstrate that H₂O diffusion in the melts is faster by 1.5-2.3 orders of magnitude than sulphur diffusion and, hence, H₂O can be considered as a rapidly diffusing oxidant while sulphur is quasi immobile in these experiments.In Raman spectra a band at 2576 cm⁻¹ appears in the sulphide - H₂O transition zone which is attributed to fundamental S-H stretching vibrations. Formation of new IR absorption bands at 5025 cm⁻¹ (on expense of the combination band of molecular H₂O at 5225 cm⁻¹) and at 3400 cm⁻¹ was observed at the front of the in-diffusing water in the sulphide bearing melt. The appearance and intensity of these two IR bands is correlated with systematic changes in S K-edge XANES spectra. A pre-edge excitation at 2466.5 eV grows with increasing H₂O concentration while the sulphide peak at 2474.0 eV decreases in intensity relative to the peak at 2477.0 eV and the feature at 2472.3 eV becomes more pronounced (all energies are relative to the sulphate excitation, calibrated to 2482.5 eV). The observations by Raman, IR and XANES spectroscopy indicate a well coordinated S²⁻ - H₂O complex which was probably formed in the glasses during cooling at the glass transition. No oxidation of sulphide was observed in any of the diffusion couple experiments. On the contrary, XANES spectra from experiments conducted with a free H₂O fluid show complete transformation of sulphide to sulphate near the melt surface and coexistence of sulphate and sulphide in the center of the melt. This can be explained by a lower H₂O activity in the diffusion couple experiments or by the need of a sink for hydrogen (e.g., a fluid which can dissolve high concentration of hydrogen) to promote oxidation of sulphide by H₂O via the reaction S²⁻ + 4H₂O = SO₄²⁻ + 4H₂. Sulphite could not be detected in any of the XANES spectra implying that this species, if it exists in the melt, it is a subordinate or transient species only.     

Influence of glass polymerisation and oxidation on micro-Raman water analysis in alumino-silicate glasses

The development of an accurate analytical procedure for determination of dissolved water in complex alumino-silicate glasses via micro-Raman analysis requires the assessment of the spectra topology dependence on glass composition. We report here a detailed study of the respective influence of bulk composition, iron oxidation state and total water content on the absolute and relative intensities of the main Raman bands related to glass network vibrations (LF: ∼490 cm⁻¹; HF: ∼960 cm⁻¹) and total water stretching (H₂O_T: ∼3550 cm⁻¹) in natural glasses. The evolution of spectra topology was examined in (i) 33 anhydrous glasses produced by the re-melting of natural rock samples, which span a very large range of polymerisation degree (NBO/T from 0.00 to 1.16), (ii) 2 sets of synthetic anhydrous basaltic glasses with variable iron oxidation state (Fe³⁺/Fe_T from 0.05 to 0.87), and (iii) 6 sets of natural hydrous glasses (C_H2OT from 0.4 to 7.0 wt%) with NBO/T varying from 0.01 to 0.76.In the explored domain of water concentration, external calibration procedure based on the H₂O_T band height is matrix-independent but its accuracy relies on precise control of the focusing depth and beam energy on the sample. Matrix-dependence strongly affects the internal calibrations based on H₂O_T height scaled to that of LF or HF bands but its effect decreases from acid (low NBO/T, SM) to basic (high NBO/T, SM) glasses. Structural parameters such as NBO/T (non-bridging oxygen per tetrahedron) and SM (sum of structural modifiers) describe the matrix-dependence better than simple compositional parameters (e.g. SiO₂, Na₂O + K₂O). Iron oxidation state has only a minor influence on band topology in basalts and is thus not expected to significantly affect the Raman determinations of water in mafic (e.g. low SiO₂, iron-rich) glasses. Modelling the evolution of the relative band height with polymerisation degree allows us to propose a general equation to predict the dissolved water content in natural glasses:

     

Sulfur speciation in natural hydrothermal waters, Iceland

The speciation of aqueous dissolved sulfur was determined in hydrothermal waters in Iceland. The waters sampled included hot springs, acid-sulfate pools and mud pots, sub-boiling well discharges and two-phase wells. The water temperatures ranged from 4 to 210 °C, the pH_T was between 2.20 and 9.30 at the discharge temperature and the SO₄ and Cl concentrations were 0.020-52.7 and <0.01-10.0 mmol kg⁻¹, respectively. The analyses were carried out on-site within ∼10 min of sampling using ion chromatography (IC) for sulfate (SO₄²⁻), thiosulfate (S₂O₃²⁻) and polythionates (S_xO₆²⁻) and titration and/or colorimetry for total dissolved sulfide (S²⁻). Sulfite (SO₃²⁻) could also be determined in a few cases using IC. Alternatively, for few samples in remote locations the sulfur oxyanions were stabilized on a resin on site following elution and analysis by IC in the laboratory. Dissolved sulfate and with few exceptions also S²⁻ were detected in all samples with concentrations of 0.02-52.7 mmol kg⁻¹ and <1-4100 μmol kg⁻¹, respectively. Thiosulfate was detected in 49 samples of the 73 analyzed with concentrations in the range of <1-394 μmol kg⁻¹ (S-equivalents). Sulfite was detected in few samples with concentrations in the range of <1-3 μmol kg⁻¹. Thiosulfate and SO₃²⁻ were not detected in <100 °C well waters and S₂O₃²⁻ was observed only at low concentrations (<1-8 μmol kg⁻¹) in ∼200 °C well waters. In alkaline and neutral pH hot springs, S₂O₃²⁻ was present in significant concentrations sometimes corresponding to up to 23% of total dissolved sulfur (S_TOT). In steam-heated acid-sulfate waters, S₂O₃²⁻ was not a significant sulfur species. The results demonstrate that S₂O₃²⁻ and SO₃²⁻ do not occur in the deeper parts of <150 °C hydrothermal systems and only in trace concentrations in ∼200-300 °C systems. Upon ascent to the surface and mixing with oxygenated ground and surface waters and/or dissolution of atmospheric O₂, S²⁻ is degassed and oxidized to SO₃²⁻ and S₂O₃²⁻ and eventually to SO₄²⁻ at pH >8. In near-neutral hydrothermal waters the oxidation of S²⁻ and the interaction of S²⁻ and S⁰ resulting in the formation of S_x²⁻ are considered important. At lower pH values the reactions seemed to proceed relatively rapidly to SO₄²⁻ and the sulfur chemistry of acid-sulfate pools was dominated by SO₄²⁻, which corresponded to >99% of S_TOT. The results suggest that the aqueous speciation of sulfur in natural hydrothermal waters is dynamic and both kinetically and source-controlled and cannot be estimated from thermodynamic speciation calculations.     

The Fe/ΣFe ratios of MORB glasses and their implications for mantle melting

The Fe³⁺/ΣFe ratio of 104 MORB glasses from the Pacific, the Atlantic, the Indian, and the Red Sea spreading centers have been determined using wet chemical Fe²⁺ analyses and electron microprobe FeO_total measurements. The data provide a new estimate for the MORB oxygen fugacity (fO₂) of 0.41 ± 0.43 (1sigma, N = 100) log units below the fayalite-magnetite-quartz buffer (FMQ), equivalent to a Fe³⁺/ΣFe = 0.12 ± 0.02 (1sigma, N = 104). This new fO₂ estimate is 0.8 log units more oxidized than the average fO₂ proposed by Christie et al. (1986) (FMQ-1.20 ± 0.44; Fe³⁺/ΣFe = 0.07 ± 0.01; N = 87). This slight difference may be related in part to the 3.5% underestimation of the Fe²⁺ concentration determined by Christie et al. (1986) compared with this study. MORB oxygen fugacity does not display any significant difference between the three main oceanic domains, or between enriched and depleted MORB. Yet, the iron red-ox state ratio shows a broad increase during fractional crystallization. Detailed study of magmatic suites highlights the lack of systematic Fe³⁺/ΣFe ratio fractionation during differentiation. Despite the large variations of inferred partial melting degrees (from 5 to 20%), the present data set does not provide any evidence of Fe³⁺/ΣFe relationships with partial melting proxies such as Na_8.0.Based on the Fe³⁺ systematics during partial melting, it is suggested that the oxidation state of MORB reflects a “buffered mantle melting process” resulting in the apparent compatible behavior of Fe³⁺ during partial melting, and in the relatively constant Fe³⁺/ΣFe ratio irrespective of the extent of melting. This result implies that partial melting processes may be open relative to oxygen. We propose a model where the Fe³⁺/ΣFe ratio in the melt is buffered during partial melting. The MORB Fe₂O₃ systematics can be accounted for by using a fO₂ of FMQ-1 that is equivalent to the average fO₂ reported for abyssal peridotites.     

Redox equilibria of iron and silicate melt structure: Implications for olivine/melt element partitioning

Olivine/melt partitioning of ΣFe, Fe²⁺, Mg²⁺, Ca²⁺, Mn²⁺, Co²⁺, and Ni²⁺ has been determined in the systems CaO-MgO-FeO-Fe₂O₃-SiO₂ (FD) and CaO-MgO-FeO-Fe₂O₃-Al₂O₃-SiO₂ (FDA3) as a function of oxygen fugacity (f_O2) at 0.1 MPa pressure. Total iron oxide content of the starting materials was ∼20 wt%. The f_O2 was to used to control the Fe³⁺/ΣFe (ΣFe: total iron) of the melts. The Fe³⁺/ΣFe and structural roles of Fe²⁺ and Fe³⁺ were determined with ⁵⁷Fe resonant absorption Mössbauer spectroscopy. Changes in melt polymerization, NBO/T, as a function of f_O2 was estimated from the Mössbauer data and existing melt structure information. It varies by ∼100% in melts coexisting with olivine in the FDA3 system and by about 300% in the FD system in the Fe³⁺/ΣFe range of the experiments (0.805-0.092). The partition coefficients ( in olivine/wt% in melt) are systematic functions of f_O2 and, therefore, NBO/T of the melt. There is a -minimum in the FDA3 system at NBO/T-values corresponding to intermediate Fe³⁺/ΣFe (0.34-0.44). In the Al-free system, FD, where the NBO/T values of melts range between ∼1 and ∼2.9, the partition coefficients are positively correlated with NBO/T (decreasing Fe³⁺/ΣFe). These relationships are explained by consideration of solution behavior in the melts governed by Qⁿ-unit distribution and structural changes of the divalent cations in the melts (coordination number, complexing with Fe³⁺, and distortion of the polyhedra).     

Landscape variations in stream water SO4 and δSSO4 in a boreal stream network

Despite reduced anthropogenic deposition during the last decades, deposition sulphate may still play an important role in the biogeochemical cycles of S and many catchments may act as net sources of S that may remain for several decades. The aim of this study is to elucidate the temporal and spatial dynamics of both SO₄²⁻ and δ³⁴S_SO4 in stream water from catchments with varying percentage of wetland and forest coverage and to determine their relative importance for catchment losses of S. Stream water samples were collected from 15 subcatchments ranging in size from 3 to 6780 ha, in a boreal stream network, northern Sweden. In forested catchments (<2% wetland cover) S-SO₄²⁻ concentrations in stream water averaged 1.7 mg L⁻¹ whereas in wetland dominated catchments (>30% wetland cover) the concentrations averaged 0.3 mg L⁻¹. A significant negative relationship was observed between S-SO₄²⁻ and percentage wetland coverage (r² = 0.77, p < 0.001) and the annual export of stream water SO₄²⁻ and wetland coverage (r² = 0.76, p < 0.001). The percentage forest coverage was on the other hand positively related to stream water SO₄²⁻ concentrations and the annual export of stream water SO₄²⁻ (r² = 0.77 and r² = 0.79, respectively). The annual average δ³⁴S_SO4 value in wetland dominated streams was +7.6‰ and in streams of forested catchments +6.7‰. At spring flood the δ³⁴S_SO4 values decreased in all streams by 1‰ to 5‰. The δ³⁴S_SO4 values in all streams were higher than the δ³⁴S_SO4 value of +4.7‰ in precipitation (snow). The export of S ranged from 0.5 kg S ha⁻¹ yr⁻¹ (wetland headwater stream) to 3.8 kg S ha⁻¹ yr⁻¹ (forested headwater stream). With an average S deposition in open field of 1.3 kg S ha⁻¹ yr⁻¹ (2002-2006) the mass balance results in a net export of S from all catchments, except in catchments with >30% wetland. The high temporal and spatial resolution of this study demonstrates that the reducing environments of wetlands play a key role for the biogeochemistry of S in boreal landscapes and are net sinks of S. Forested areas, on the other hand were net sources of S.