Oxidising agents in sub-arc mantle melts link slab devolatilisation and arc magmas

Mantle oxidation by sulfur drives the formation of giant gold deposits in subduction zones

Oxidation of the sub-arc mantle driven by slab-derived fluids has been hypothesized to contribute to the formation of gold deposits in magmatic arc environments that host the majority of metal resources on Earth. However, the mechanism by which the infiltration of slab-derived fluids into the mantle wedge changes its oxidation state and affects Au enrichment remains poorly understood. Here, we present the results of a numerical model that demonstrates that slab-derived fluids introduce large amounts of sulfate (S6+) into the overlying mantle wedge that increase its oxygen fugacity by up to 3 to 4 log units relative to the pristine mantle. Our model predicts that as much as 1 wt.% of the total dissolved sulfur in slab-derived fluids reacting with mantle rocks is present as the trisulfur radical ion, S3-. This sulfur ligand stabilizes the aqueous Au(HS)S3- complex, which can transport Au concentrations of several grams per cubic meter of fluid. Such concentrations are more than three orders of magnitude higher than the average abundance of Au in the mantle. Our data thus demonstrate that an aqueous fluid phase can extract 10 to 100 times more Au than in a fluid-absent rock-melt system during mantle partial melting at redox conditions close to the sulfide-sulfate boundary. We conclude that oxidation by slab-derived fluids is the primary cause of Au mobility and enrichment in the mantle wedge and that aqueous fluid-assisted mantle melting is a prerequisite for formation of Au-rich magmatic hydrothermal and orogenic gold systems in subduction zone settings.

Sulfur species and gold transport in arc magmatic fluids

The sulfur species present in magmatic fluids affect the global redox cycle, the Earth's climate and the formation of some of the largest and most economic ore deposits of critical metals. However, the speciation of sulfur under conditions that are relevant for upper crustal magma reservoirs is unclear. Here we combine a prototype pressure vessel apparatus and Raman spectroscopy to determine sulfur speciation in arc magmatic fluid analogues in situ over a range of geologically relevant pressure-temperature-redox conditions. We find that HS-, H2S and SO2 are the main sulfur species in the experimental fluids, while the concentrations of S6+ species and S2 - and S3 - sulfur radical ions are negligible, in contrast to previous experimental work. The measured gold solubilities in the experimental fluids are highest when sulfur is dominantly present as HS- and H2S species and greatly exceed thermodynamic predictions. Our results imply that HS-, rather than sulfur radicals, accounts for the high solubilities of gold in magmatic-hydrothermal fluids. We also find that magmatic sulfur degassing is a redox process under oxidizing conditions and may lead to additional magma oxidation beyond that imparted by slab-derived fluxes and crystallization.

Mobilization of isotopically heavy sulfur during serpentinite subduction

Primitive arc magmas are more oxidized and enriched in sulfur-34 (34S) compared to mid-ocean ridge basalts. These findings have been linked to the addition of slab-derived volatiles, particularly sulfate, to arc magmas. However, the oxidation state of sulfur in slab fluids and the mechanisms of sulfur transfer in the slab remain inconclusive. Juxtaposed serpentinite and eclogitic metagabbro from the Voltri Massif (Italy) provide evidence for sulfur mobilization and associated redox processes during infiltration of fluids. Using bulk rock and in situ δ34S measurements, combined with thermodynamic calculations, we document the transfer of bisulfide-dominated, 34S-enriched fluids in equilibrium with serpentinite into adjacent metagabbro. We argue that the process documented in this study is pervasive along the subduction interface and infer that subsequent melting of these reacted slab-mantle interface rocks could produce melts that display the characteristic oxygen fugacity and sulfur isotope signatures of arc magmas worldwide.

Incipient carbonate melting drives metal and sulfur mobilization in the mantle

We present results from high-pressure, high-temperature experiments that generate incipient carbonate melts at mantle conditions (~90 kilometers depth and temperatures between 750° and 1050°C). We show that these primitive carbonate melts can sequester sulfur in its oxidized form of sulfate, as well as base and precious metals from mantle lithologies of peridotite and pyroxenite. It is proposed that these carbonate sulfur-rich melts may be more widespread than previously thought and that they may play a first-order role in the metallogenic enhancement of localized lithospheric domains. They act as effective agents to dissolve, redistribute, and concentrate metals within discrete domains of the mantle and into shallower regions within Earth, where dynamic physicochemical processes can lead to ore genesis at various crustal depths.

Sulfur disproportionation in deep COHS slab fluids drives mantle wedge oxidation

Sulfur degassed at volcanic arcs calls for dissolved sulfate ions (S6+) released by subduction-zone fluids, oxidizing (in association with carbon) the subarc mantle, but sulfur speciation in subduction fluids at subarc depths remains unclear. We apply electrolytic fluid thermodynamics to model the dissolution behavior of pyrite in metacarbonate sediments as a function of P, T and rock redox state up to 4.3 gigapascals and 730°C. At subarc depth and the redox conditions of the fayalite-magnetite-quartz oxygen buffer, pyrite dissolution releases oxidized sulfur in fluids by disproportionation into sulfate, bisulfite, and sulfide species. These findings indicate that oxidized, sulfur-rich carbon-oxygen-hydrogen-sulfur (COHS) fluids form within subducting slabs at depths greater than 100 kilometers independent from slab redox state and that sulfur can be more effective than the concomitantly dissolved carbon at oxidizing the mantle wedge, especially when carbonates are stable.

Sulfur and chlorine budgets control the ore fertility of arc magmas

Continental arc magmas supply the ore-forming element budget of most globally important porphyry-type ore deposits. However, the processes enabling certain arc segments to preferentially generate giant porphyry deposits remain highly debated. Here we evaluate the large-scale covariation of key ore-forming constituents in this setting by studying silicate melt inclusions in volcanic rocks from a fertile-to-barren segment of the Andean Southern Volcanic Zone (33–40 °S). We show that the north-to-south, fertile-to-barren gradient is characterized by a northward increase in S and Cl concentrations and a simultaneous decrease in Cu. Consequently, we suggest that the concentration of S and Cl rather than the concentration of ore metals regulates magmatic-hydrothermal ore fertility, and that the loss of volatiles prior to arrival in the upper crust impacts ore-forming potential more than magmatic sulfide saturation-related ore metal scavenging.

Earth’s largest copper deposits form in continental arcs, yet it is not well understood what determines whether a magmatic system generates economic mineralization or not. Here the authors show that the abundance of chlorine and sulfur, rather than the abundance of ore metals controls magmatic ore fertility.

Oxidized primary arc magmas: Constraints from Cu/Zr systematics in global arc volcanics

Arc volcanics are more oxidized than mid-ocean ridge basalts (MORB), but it is debated whether this is a mantle feature or a result of magmatic evolution. Copper, a sulfur-loving element, has been used to trace the behavior of redox-sensitive sulfur during mantle melting and infer similar redox states of sub-arc and sub-ridge mantle. Previous studies, however, neglected elevated sulfur contents in the sub-arc mantle, leading to underestimation of oxygen fugacities, and did not recognize systematic Cu variations in arc volcanics. Here, we show that the Cu/Zr ratio is a sensitive indicator that responds to sulfur content, oxygen fugacity, and extent of melting of the mantle. Because of higher mantle S contents, Cu systematics of arc magmas require one log unit higher oxygen fugacities of sub-arc than sub-ridge mantle. Low Cu contents of thick-crusted arc volcanics result from low extents of melting of sulfur-rich mantle, obviating the need for deep crustal sulfide fractionation, with substantial implications for the origin of porphyry-Cu deposits.

Potassium distribution and isotope composition in the lithospheric mantle in relation to global Earth's reservoirs

Recent analytical advances have provided means to measure potassium (K) isotopes in various terrestrial materials, but little is known about K distribution and stable isotope composition in the lithospheric mantle because of: (a) common low K abundances, (b) potential contamination and alteration, (c) diversity of mantle rocks and minerals hosting K in different tectonic settings. We report K abundances and δ41K values for well-studied whole-rock (WR) mantle xenoliths (spinel and garnet peridotites and pyroxenites) from mobile belts, a craton, a subduction zone, as well as for K-rich phases (mica, amphibole, silicate glass) and xenolith-bearing volcanic materials (67 samples). The xenolith materials show extremely broad ranges of K content (7 μg/g to 6.6 wt.%) and δ41K (-2.77‰ to 0.62‰). They contrast with the narrow δ41K range for host volcanic materials (-0.53‰ to -0.27‰) and literature data on oceanic basalts (melting products of upwelling asthenosphere: -0.43 ± 0.17‰, 2sd). Amphibole-bearing subduction zone peridotites show the highest WR δ41K values (0.40 to 0.62‰) likely inherited from fluids released into the mantle wedge from subducted oceanic crust. All other WR samples yield negative δ41K: -0.06‰ to -2.77‰ for peridotites and -0.17‰ to -0.52‰ for pyroxenites. The δ41K in K-rich mantle phases range from positive values (0.16 to 0.57‰) for phlogopite in strongly metasomatized peridotites to negative values (-0.27 to -0.94‰) for phlogopite, amphibole and glass pockets from other samples, which cannot be explained by equilibrium inter-mineral fractionation inferred from ab initio calculations. We attribute the broad δ41K variations to (a) isotope fractionation during fluid-rock interaction in the mantle, and (b) distinct sources of K-bearing fluids. Kinetic isotope fractionation during fluid percolation and diffusion is inferred for composite xenoliths (phlogopite and pyroxenite veins in peridotites) that have lower δ41K in the hosts than in the veins due to slower migration of 41K than 39K from the veins (former fluid channels) to host mantle. Overall, the K isotope fractionation in the lithospheric mantle appears to be greater than for magmatic fractionation in the crust. The average δ41K of normal off-craton continental lithospheric mantle calculated from the least altered fertile and lightly metasomatized lherzolites is -0.57 ± 0.28‰ (2sd). This value is within error (though a little lower) of estimates for both continental crust and MORB and OIB mantle sources indicating that these major silicate Earth reservoirs have similar bulk δ41K values, apparently due to low or negligible K isotopic fractionation during their formation by magmatic differentiation and melting. By contrast, K isotopes in modern and fossil subduction zones are fractionated via fluid-related equilibrium and kinetic processes.

Oxidized sulfur-rich arc magmas formed porphyry Cu deposits by 1.88 Ga

Most known porphyry Cu deposits formed in the Phanerozoic and are exclusively associated with moderately oxidized, sulfur-rich, hydrous arc-related magmas derived from partial melting of the asthenospheric mantle metasomatized by slab-derived fluids. Yet, whether similar metallogenic processes also operated in the Precambrian remains obscure. Here we address the issue by investigating the origin, fO₂, and S contents of calc-alkaline plutonic rocks associated with the Haib porphyry Cu deposit in the Paleoproterozoic Richtersveld Magmatic Arc (southern Namibia), an interpreted mature island-arc setting. We show that the ca. 1886–1881 Ma ore-forming magmas, originated from a mantle-dominated source with minor crustal contributions, were relatively oxidized (1‒2 log units above the fayalite-magnetite-quartz redox buffer) and sulfur-rich. These results indicate that moderately oxidized, sulfur-rich arc magma associated with porphyry Cu mineralization already existed in the late Paleoproterozoic, probably as a result of recycling of sulfate-rich seawater or sediments from the subducted oceanic lithosphere at that time.

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

Uncovering and quantifying the subduction zone sulfur cycle from the slab perspective

Sulfur belongs among H₂O, CO₂, and Cl as one of the key volatiles in Earth’s chemical cycles. High oxygen fugacity, sulfur concentration, and δ³⁴S values in volcanic arc rocks have been attributed to significant sulfate addition by slab fluids. However, sulfur speciation, flux, and isotope composition in slab-dehydrated fluids remain unclear. Here, we use high-pressure rocks and enclosed veins to provide direct constraints on subduction zone sulfur recycling for a typical oceanic lithosphere. Textural and thermodynamic evidence indicates the predominance of reduced sulfur species in slab fluids; those derived from metasediments, altered oceanic crust, and serpentinite have δ³⁴S values of approximately −8‰, −1‰, and +8‰, respectively. Mass-balance calculations demonstrate that 6.4% (up to 20% maximum) of total subducted sulfur is released between 30–230 km depth, and the predominant sulfur loss takes place at 70–100 km with a net δ³⁴S composition of −2.5 ± 3‰. We conclude that modest slab-to-wedge sulfur transport occurs, but that slab-derived fluids provide negligible sulfate to oxidize the sub-arc mantle and cannot deliver ³⁴S-enriched sulfur to produce the positive δ³⁴S signature in arc settings. Most sulfur has negative δ³⁴S and is subducted into the deep mantle, which could cause a long-term increase in the δ³⁴S of Earth surface reservoirs.

Sulfur is one of the key volatiles in Earth’s chemical cycles; however, sulfur speciation, isotopic composition, and flux during the subduction cycle remain unclear. Here, the authors provide direct constraints on subduction zone sulfur recycling from high-pressure rocks and explore implications for arc magmatism.

Subducting serpentinites release reduced, not oxidized, aqueous fluids

The observation that primitive arc magmas are more oxidized than mid-ocean-ridge basalts has led to the paradigm that slab-derived fluids carry SO₂ and CO₂ that metasomatize and oxidize the sub-arc mantle wedge. We combine petrography and thermodynamic modelling to quantify the oxygen fugacity (fO₂) and speciation of the fluids generated by serpentinite dehydration during subduction. Silicate-magnetite assemblages maintain fO₂ conditions similar to the quartz-fayalite-magnetite (QFM) buffer at fore-arc conditions. Sulphides are stable under such conditions and aqueous fluids contain minor S. At sub-arc depth, dehydration occurs under more reducing conditions producing aqueous fluids carrying H₂S. This finding brings into question current models in which serpentinite-derived fluids are the cause of oxidized arc magmatism and has major implications for the global volatile cycle, as well as for redox processes controlling subduction zone geodynamics.

Molybdenum systematics of subducted crust record reactive fluid flow from underlying slab serpentine dehydration

Fluids liberated from subducting slabs are critical in global geochemical cycles. We investigate the behaviour of Mo during slab dehydration using two suites of exhumed fragments of subducted, oceanic lithosphere. Our samples display a positive correlation of δ^98/95Mo_{NIST 3134} with Mo/Ce, from compositions close to typical mantle (−0.2‰ and 0.03, respectively) to very low values of both δ^98/95Mo_{NIST 3134} (−1‰) and Mo/Ce (0.002). Together with new, experimental data, we show that molybdenum isotopic fractionation is driven by preference of heavier Mo isotopes for a fluid phase over rutile, the dominant mineral host of Mo in eclogites. Moreover, the strongly perturbed δ^98/95Mo_{NIST 3134} and Mo/Ce of our samples requires that they experienced a large flux of oxidised fluid. This is consistent with channelised, reactive fluid flow through the subducted crust, following dehydration of the underlying, serpentinised slab mantle. The high δ^98/95Mo_{NIST 3134} of some arc lavas is the complement to this process.

Fluid liberation and migration from subducted oceanic slabs play a critical role in arc magmatism but the volume and origin of the released fluids is unclear and difficult to trace. Here, the authors use Molybdenum isotope ratios to tackle these problems.

A metasomatized lithospheric mantle control on the metallogenic signature of post-subduction magmatism

Ore deposits are loci on Earth where energy and mass flux are greatly enhanced and focussed, acting as magnifying lenses into metal transport, fractionation and concentration mechanisms through the lithosphere. Here we show that the metallogenic architecture of the lithosphere is illuminated by the geochemical signatures of metasomatised mantle rocks and post-subduction magmatic-hydrothermal mineral systems. Our data reveal that anomalously gold and tellurium rich magmatic sulfides in mantle-derived magmas emplaced in the lower crust share a common metallogenic signature with upper crustal porphyry-epithermal ore systems. We propose that a trans-lithospheric continuum exists whereby post-subduction magmas transporting metal-rich sulfide cargoes play a fundamental role in fluxing metals into the crust from metasomatised lithospheric mantle. Therefore, ore deposits are not merely associated with isolated zones where serendipitous happenstance has produced mineralisation. Rather, they are depositional points along the mantle-to-upper crust pathway of magmas and hydrothermal fluids, synthesising the concentrated metallogenic budget available.

Arc magmas oxidised by water dissociation and hydrogen incorporation in orthopyroxene

Elevated H₂O concentrations and oxygen fugacities are two fundamental properties that distinguish magmas formed in subduction zones from new crust generated at mid ocean ridges. However, the mechanism of magma oxidation, and how it relates to the increase in H₂O remains unclear. In this study, we use infrared spectroscopy of mantle wedge orthopyroxene to trace the temporal and spatial evolution of oxygen fugacity during transport of hydrous arc melts towards the crust. A positive correlation between equilibrium oxygen fugacity and orthopyroxene H₂O concentrations for the peridotite samples studied allowed the assignment of specific, commonly-observed absorption bands to redox-sensitive crystallographic defects. H₂O content associated with these redox-sensitive defects increases in concentration across individual crystals, uniquely preserving the time-dependent transition from reduced to oxidised conditions during the migration of hydrous melts through the mantle wedge. A separate, but related process of reaction with H₂ occurring primarily during the earliest stages of melt-mantle reaction may be fundamental in generating the oxidised nature of hydrous melts. Our study proposes that the oxidised nature of arc magmas may not be a primary feature, but is instead acquired progressively as hydrous primary melts react with the surrounding mantle.

Neoproterozoic to early Phanerozoic rise in island arc redox state due to deep ocean oxygenation and increased marine sulfate levels

A rise in atmospheric O2 levels between 800 and 400 Ma is thought to have oxygenated the deep oceans, ushered in modern biogeochemical cycles, and led to the diversification of animals. Over the same time interval, marine sulfate concentrations are also thought to have increased to near-modern levels. We present compiled data that indicate Phanerozoic island arc igneous rocks are more oxidized (Fe3+/ΣFe ratios are elevated by 0.12) vs. Precambrian equivalents. We propose this elevation is due to increases in deep-ocean O2 and marine sulfate concentrations between 800 and 400 Ma, which oxidized oceanic crust on the seafloor. Once subducted, this material oxidized the subarc mantle, increasing the redox state of island arc parental melts, and thus igneous island arc rocks. We test this using independently compiled V/Sc ratios, which are also an igneous oxybarometer. Average V/Sc ratios of Phanerozoic island arc rocks are elevated (by +1.1) compared with Precambrian equivalents, consistent with our proposal for an increase in the redox state of the subarc mantle between 800 and 400 Ma based on Fe3+/ΣFe ratios. This work provides evidence that the more oxidized nature of island arc vs. midocean-ridge basalts is related to the subduction of material oxidized at the Earth's surface to the subarc mantle. It also indicates that the rise of atmospheric O2 and marine sulfate to near-modern levels by the late Paleozoic influenced not only surface biogeochemical cycles and animal diversification but also influenced the redox state of island arc rocks, which are building blocks of continental crust.