Water and the oxidation state of subduction zone 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.

Highly oxidized intraplate basalts and deep carbon storage

Deep carbon cycle is crucial for mantle dynamics and maintaining Earth's habitability. Recycled carbonates are a strong oxidant in mantle carbon-iron redox reactions, leading to the formation of highly oxidized mantle domains and deep carbon storage. Here we report high Fe3+/∑Fe values in Cenozoic intraplate basalts from eastern China, which are correlated with geochemical and isotopic compositions that point to a common role of carbonated melt with recycled carbonate signatures. We propose that the source of these highly oxidized basalts has been oxidized by carbonated melts derived from the stagnant subducted slab in the mantle transition zone. Diamonds formed during the carbon-iron redox reaction were separated from the melt due to density differences. This would leave a large amount of carbon (about four times of preindustrial atmospheric carbon budget) stored in the deep mantle and isolated from global carbon cycle. As such, the amounts of subducted slabs stagnated at mantle transition zone can be an important factor regulating the climate.

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.

Latitude-dependent oxygen fugacity in arc magmas

The redox state of arc mantle has been considered to be more oxidized and diverse than that of the mid-ocean ridge, but the cause of the variation is debated. We examine the redox state of the Cenozoic global arc mantle by compiling measured/calculated fO2 of olivine-hosted melt inclusions from arc magma and modeled fO2 based on V/Sc and Cu/Zr ratios of arc basaltic rocks. The results indicate that the redox state of Cenozoic arc mantle is latitude dependent, with less oxidized arc mantle in the low latitudes, contrasting with a near constant across-latitude trend in the mid-ocean ridges. We propose that such a latitude-dependent pattern in the arc mantle may be controlled by the variation in the redox state of subducted sediment, possibly related to a latitudinal variation in the primary production of phytoplankton, which results in more organic carbon and sulfide deposited on the low-latitude ocean floor. Our findings provide evidence for the impact of the surface environment on Earth's upper mantle.

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.

Zircons reveal the history of fluctuations in oxidation state of crustal magmatism and supercontinent cycle

We apply a zircon redox index to a global compilation of detrital zircons to track the variation of oxidation state, expressed as ΔFMQ, through Earth's history. Those from I-type rocks, which comprise mantle and crustal igneous protoliths, including tonalite-trondhjemite-granodiorites (TTGs), generally have a high oxidation state (ΔFMQ > 0). In contrast, zircons from igneous rocks derived from supracrustal source rocks (S-type) are commonly reduced (ΔFMQ < 0). With the probability density function derived from the Gaussian-Kernel-Density-Estimation, we use the maximum likelihood estimation (MLE) to distinguish S-type from I-type zircons through Earth's history using zircon redox. Voluminous S-type magma production shows a ca. 600 Ma cyclicity that is closely related to the supercontinent cycle. We link a cyclic drop in redox values after 2.6 Ga to periodic S-type magma generation associated with burial and melting of metasedimentary rocks during supercontinent assembly and amalgamation. The ΔFMQ of the detrital zircons rise at ∼3.5 Ga followed by a consistent average ΔFMQ > 0 over the last 3 Ga. Given that the redox state of magmas is independent of crustal thickness and silica variation, and elevated values are likely more closely related to tectonic setting, we suggest that the consistent average ΔFMQ > 0 from ca. 3.5 Ga onwards relates to recycling of oceanic lithosphere back into the mantle in what eventually became established as subduction zones. The more reduced magmas associated with sedimentary sources, became established at 2.6 Ga, presumably in response to continental rocks rising above sea-level, and follow peaks of productivity associated with the supercontinent cycle.

Earth's early continental crust formed from wet and oxidizing arc magmas

Formation of continental crust has shaped the surface and interior of our planet and generated the land and mineral resources on which we rely. However, how the early continental crust of Earth formed is still debated1-7. Modern continental crust is largely formed from wet and oxidizing arc magmas at subduction zones, in which oceanic lithosphere and water recycle into the mantle8-10. The magmatic H2O content and redox state of ancient rocks that constitute the early continental crust, however, are difficult to quantify owing to ubiquitous metamorphism. Here we combine two zircon oxybarometers11,12 to simultaneously determine magmatic oxygen fugacity (fO2) and H2O content of Archaean (4.0-2.5 billion years ago) granitoids that dominate the early continental crust. We show that most Archaean granitoid magmas were ≥1 log unit more oxidizing than Archaean ambient mantle-derived magmas13,14 and had high magmatic H2O contents (6-10 wt%) and high H2O/Ce ratios (>1,000), similar to modern arc magmas. We find that magmatic fO2, H2O contents and H2O/Ce ratios of Archaean granitoids positively correlate with depth of magma formation, requiring transport of large amounts of H2O to the lower crust and mantle. These observations can be readily explained by subduction but are difficult to reconcile with non-subduction models of crustal formation3-7. We note an increase in magmatic fO2 and H2O content between 4.0 and 3.6 billion years ago, probably indicating the onset of subduction during this period.

Copper isotope evidence for sulfide fractionation and lower crustal foundering in making continental crust

The continental crust is strongly depleted in copper compared with its building blocks-primary arc magmas-and this depletion is intrinsically associated with continental crust formation. However, the process by which Cu removal occurs remains enigmatic. Here we show, using Cu isotopes, that subduction-zone processes and mantle melting produce limited fractionation of Cu isotopes in arc magmas, and, instead, the heterogeneous Cu isotopic compositions of lower crustal rocks, which negatively correlate with Cu contents, suggest segregation or accumulation of isotopically light sulfides during intracrustal differentiation of arc magmas. This is supported by the extremely light Cu isotopic compositions of lower crustal mafic cumulates and heavy Cu isotopic compositions of differentiated magmas in thick continental arcs. Intracrustal differentiation of mantle-derived magmas and subsequent foundering of sulfide-rich mafic cumulates preferentially removes isotopically light Cu, leaving a Cu-depleted and isotopically heavy continental crust.

Garnet crystallization does not drive oxidation at arcs

Arc magmas, the building blocks of continental crust, are depleted in total iron (Fe), have higher ratios of oxidized Fe to total Fe (Fe³⁺/∑Fe), and record higher oxygen fugacities (fO₂'s) compared with magmas erupted at mid-ocean ridges. Garnet crystallization could explain these observations if garnet removes substantial amounts of Fe²⁺, but not Fe³⁺, from magma, yet this model for continental crust generation has never been tested experimentally. Analysis of garnets and melts in laboratory experiments show that the compatibilities of Fe²⁺ and Fe³⁺ in garnet are of similar magnitudes. Our results indicate that fractional crystallization of garnet-bearing cumulates will remove 20% of total Fe from primary arc basalts but negligibly alter the Fe³⁺/∑Fe ratio and fO₂ of the melt. Garnet crystallization is unlikely to be responsible for the relatively oxidized nature of basaltic arc magmas or the Fe-depletion trend observed in continental crust.

Petrogenesis of the Early Cretaceous Zhouguan Granodiorite in Jiaodong Peninsula: Evidence from Mineralogy, Geochemistry, Geochronology, and Sr-Nd Isotopes

The Jiaodong Peninsula is the most important gold mineralization area in China, and the formation of gold deposits is closely related to granitoids. The isotopic ages of the Early Cretaceous granodiorites in the northwestern Jiaodong Peninsula are concentrated in the range of 111~123 Ma, and are coeval with the formation of the gold deposits in the area. However, the studies on the geotectonic settings of the granodiorites, especially their petrogenesis and relationship with gold deposits in the northwestern Jiaodong Peninsula, are scarce. Based on field and petrographic observations, geochemistry, EPMA analysis, zircon U-Pb chronology, and Sr-Nd isotopes of the Early Cretaceous Zhouguan granodiorite in the Jiaodong area, the formation age of Zhouguan granodiorite is determined as 115 Ma ± 0.77 Ma; the analysis of EPMA shows that biotite is mainly composed of Fe-biotite and Mg-biotite, with its MgO content ranging from 9.797% to 11.635%. The crystallization temperature of biotite is in the range of 500 °C~625 °C and the emplacement depth of the rock mass is 3.98~8.71 km. The amphibole in the mass mainly includes magnesiohornblende, pargasite, and magnesiosadanagaite; among them, the former two are of crustal origin, while magnesiosadanagaite is of mantle origin. The crystallization pressure and depth of the former two are in the range of 0.75~3.02 kbar and 2.81~11.4 km, respectively, while the crystallization pressure and depth for the latter is 4.64 kbar and 17.53 km, respectively. The (87Sr/86Sr) values range from 0.710424 to 0.711074 and the (143Nd/144Nd) values range from 0.511530 to 0.511808. The parental magma of the Zhouguan granodiorite is highly oxidized with high-water content that is favorable for Au enrichment. Combined with the Nb-Y and Yb-Ta diagrams, a model describing the formation of Zhouguan granodiorite is proposed.

Volatile Characteristics of Apatite in Dacite from the Eastern Manus Basin and Their Geological Implications

As one of the youngest back-arc basins, the evolutionary behavior of magmatic volatiles in the Eastern Manus Basin has been poorly studied. Recently, apatite has received widespread attention for its powerful function in recording information on magmatic volatiles. In this paper, by determining the major element compositions and primary volatile abundances (F, Cl, SO3) of apatites in dacite, we analyze the compositions of volatiles before magma eruption in the Eastern Manus Basin, as well as their indications for magmatic oxygen fugacity and petrogenesis, so as to improve the study about the evolution of magmatic volatiles in this region. Experimental data indicate that apatite saturation temperatures range from 935 to 952 °C. All the apatites are magmatic apatites with F contents of 0.87−1.39 wt.%, Cl contents of 1.24−1.70 wt.%, and SO3 ≤ 0.06 wt.%. Analysis reveals that the apatites in this study crystallized from volatile-undersaturated melts, so their chemical compositions can be used as indicators of magmatic compositions. According to calculations, the minimum S concentrations of the host melts range from 2−65 ppm or 8−11 ppm. The crystallization and separation of magnetite caused the reduction state of melts, and the relatively low oxygen fugacity (ΔFMQ = −0.2 ± 0.9) caused low SO3 contents in apatites. In addition, F and Cl contents of the host melt were calculated to be 185−448 ppm and 1059−1588 ppm, corresponding to the H2O contents of 1.4−2.1% and 1.2−1.5% (error ± 30−40%), respectively. The high Cl/F ratio and H2O contents of samples indicate the addition of slab-derived fluids in the mantle source region of the Eastern Manus Basin. High F contents of the melts may be influenced by F-rich sediments, as well as the release of F from lawsonite and phengite decomposition. High Cl appears to originate from the dual influence of subduction-released fluids and Cl-rich seawater-derived components. Further, it is estimated that 14−21% of the total Cl concentrations in melts were added directly from subduction-released fluids, or higher.

Genesis and Metallogenic Characteristics of the Zhunsujihua Granitic Intrusions in Sonid Left Banner, Inner Mongolia, China

The Zhunsujihua porphyry Mo-Cu deposit is located in the Erenhot–East Ujimqin metallogenic belt in northeastern China. Granodioritic intrusions in the mining area are dominated by granodiorite and granodiorite porphyry, but the Mo mineralization is limited within the granodiorite. Zircon LA-ICP-MS U-Pb dating yields crystallization ages of 301.5 ± 3.0 Ma for granodiorite and 296.0 ± 3.0 Ma for granodiorite porphyry. These ages constrain the magmatic activity at the Zhunsujihua deposit that took place during the subduction of the Paleo-Asian oceanic plate. Whole-rock geochemical data suggest that the granodioritic intrusions belong to calc-alkaline and high-K calc-alkaline series, and are characterized by enrichment in K, Rb, U, Th, and Pb, and depletion in Nb, Ta, Ti, and P. The negative Eu, Ba, and Sr anomalies suggest that they have experienced extensive fractionation of plagioclase. Trace element compositions of zircons from the Zhunsujihua deposit provide constraints on the oxygen fugacity (ƒO2) of the magma, which is shown to high values with ∆FMQ = +0.5 to +5.6. The wide range of zircon εHf (t) (+1.3~+9.4) values, positive whole-rock εNd (t) (+2.5~+3.9) values, and relatively low initial (87Sr/86Sr)i (0.70367~0.70561) ratios indicate that the magmas mainly originated from a juvenile lower crust source derived from depleted mantle, but mixed with pre-existing crustal components. Moreover, the juvenile lower crust represents the main source of Mo for the Zhunsujihua deposit. A high magmatic oxygen fugacity and fractional crystallization played key roles in forming the Zhunsujihua deposit.

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.

Formation of Gold Alloys during Crustal Differentiation of Convergent Zone Magmas: Constraints from an AU-Rich Websterite in the Stanovoy Suture Zone (Russian Far East)

Gold is typically transported by mafic and evolved magmas into the upper crust to be deposited in shallow oxidized porphyry and epithermal environments. However, the magmatic behavior of gold is still poorly understood and warrants further attention. Additional insights into the magmatic evolution of gold and other noble metals can be provided by investigations of primitive convergent zone magmas and products of their differentiation that contain primary-textured Au-alloys. One of the best examples of such Au-rich ultramafic cumulates is the Triassic (232–233 Ma) Ildeus intrusion, which was emplaced within the Mesozoic Stanovoy subduction zone in the Russian Far East. Some websterites from the Ildeus intrusion, representing cumulates crystallized from a primitive convergent zone magma, are enriched in Au (up to 596 ppm) and contain abundant Cu-Ag-Au micro-particles. Most of these Au-alloy micro-particles display compositions similar to those previously found in explosive pyroclastic rocks in the Lesser Khingan iron district, mantle wedge peridotites in Kamchatka and Cretaceous adakites in the Stanovoy suture zone. Textural and compositional characteristics suggest that Cu-Ag-Au alloys precipitated from a primitive calc-alkaline melt during its crustal differentiation in a Mesozoic paleo-subduction zone. Some large Cu-Ag-Au grains display an internal honeycomb-like structure with alternating Cu-rich and Cu-poor zones. Heating experiments under atmospheric conditions recorded a substantial loss of Cu from primary magmatic Cu-Ag-Au alloys, which appears to be a process characteristic of oxidized hydrothermal ore systems. We suggest that the later-stage hydrothermal alteration of differentiated igneous conduits containing magmatic gold alloys results in the formation of Cu-free gold mineralization comparable to the upper crustal porphyry and epithermal environments.

Evidence for the oxidation of Earth's crust from the evolution of manganese minerals

Analysis of manganese mineral occurrences and valence states demonstrate oxidation of Earth's crust through time. Changes in crustal redox state are critical to Earth's evolution, but few methods exist for evaluating spatially averaged crustal redox state through time. Manganese (Mn) is a redox-sensitive metal whose variable oxidation states and abundance in crustal minerals make it a useful tracer of crustal oxidation. We find that the average oxidation state of crustal Mn occurrences has risen in the last 1 billion years in response to atmospheric oxygenation following a 66 ± 1 million-year time lag. We interpret this lag as the average time necessary to equilibrate the shallow crust to atmospheric oxygen fugacity. This study employs large mineralogical databases to evaluate geochemical conditions through Earth's history, and we propose that this and other mineral data sets form an important class of proxies that constrain the evolving redox state of various Earth reservoirs.

Differential Fluid Activity in a Single Exhumed Continental Subduction Unit from Local P-T-M(H2O) Records of Zoned Amphiboles (North Muya, Eastern Siberia)

The behavior of the continental lithosphere in the Alpine-type subduction zones, which primarily depends on its thickness, thermal regime of subduction and availability of fluids/melts, remains an important issue for both metamorphic petrology and geochemistry as well as for resolving the thermomechanical properties of subduction paleo-interfaces. Rehydrated (amphibole- and zoisite-bearing) eclogites from the Neoproterozoic North Muya high-pressure complex (northern Central Asian Orogenic belt, eastern Siberia) were studied in order to assess their peak burial depths, degree of prograde dehydration, and further retrograde hydration extent within a subducted and exhumed continental unit. Three medium-grained eclogites from different localities of HP complex show similarly dry peak assemblages of pyrope-almandine-grossular garnet (XGrs up to 0.30, XPrp up to 0.25) + Na-rich omphacite (up to 44 mol % of jadeite) + rutile + quartz, which are variably replaced by secondary plagioclase + clinopyroxene ± amphibole symplectitic aggregate. The eclogites were subjected to burial at similar peak depths (up to ~17–21 kbar) but different peak temperatures (within ~600–730 °C) with or without notable heating and re-equilibration due to crustal thickening. Variable degrees of exhumation-induced pervasive rehydration led to growth of individual zoned porphyroblastic barroisite-hornblende amphibole ([B]Na = 0.03–0.45) ± zoisite over the primary eclogitic assemblage or after notable thermally-driven development of symplectitic aggregate after omphacite. Amphibole compositions together with the zoisite presence/absence in different samples reflect continuous rehydration by addition of ~0.5–1.5 wt.% at different exhumation conditions, from nearly peak eclogitic P–T (~17–21 kbar) to granulite- and amphibolite-facies depths within the plagioclase stability field (<14 kbar). This diversity most likely required irregular distribution of internally sourced, low-volume, hydrous metamorphic fluid (i.e., from host felsic rocks or metasediments) acting at different depths of the subduction interface. From the performed PTX calculations, I suggest that nearly isochemical (i.e., without any significant modification of the bulk-rock composition other than incorporation of additional H2O), retrograde hydration by only at lower- to middle-crust conditions did not significantly influence the density and the rheology of the subducted continental slices due to both (1) a limited abundance of dense metabasic rocks, which are commonly more fluid-rich (e.g., due to chlorite or amphibole alteration), and (2) the initially dry nature of mafic and felsic continental rocks. The limited dehydration and rehydration scales exemplified by the North Muya eclogites and therefore low availability of hydrous metamorphic fluids may have accounted for the high buoyancy of the eclogitic crust and explained the absence of contemporaneous suprasubduction magmatism in the regional context at ca. ~630 Ma.

Serpentinite-derived slab fluids control the oxidation state of the subarc mantle

Recent geochemical evidence confirms the oxidized nature of arc magmas, but the underlying processes that regulate the redox state of the subarc mantle remain yet to be determined. We established a link between deep subduction-related fluids derived from dehydration of serpentinite ± altered oceanic crust (AOC) using B isotopes and B/Nb as fluid proxies, and the oxidized nature of arc magmas as indicated by Cu enrichment during magma evolution and V/Yb. Our results suggest that arc magmas derived from source regions influenced by a greater serpentinite (±AOC) fluid component record higher oxygen fugacity. The incorporation of this component into the subarc mantle is controlled by the subduction system’s thermodynamic conditions and geometry. Our results suggest that the redox state of the subarc mantle is not homogeneous globally: Primitive arc magmas associated with flat, warm subduction are less oxidized overall than those generated in steep, cold subduction zones.

Out-of-sequence skeletal growth causing oscillatory zoning in arc olivines

Primitive olivines from the monogenetic cones Los Hornitos, Central-South Andes, preserve dendritic, skeletal, and polyhedral growth textures. Consecutive stages of textural maturation occur along compositional gradients where high Fo–Ni cores of polyhedral olivines (Fo_92.5, Ni ~3500 ppm) contrast with the composition of dendritic olivines (Fo < 91.5, Ni < 3000 ppm), indicating sequential nucleation. Here we present a new growth model for oscillatory Fo–Ni olivine zoning that contrasts with the standard interpretation of continuous, sequential core-to-rim growth. Olivine grows rapidly via concentric addition of open-structured crystal frames, leaving behind compositional boundary layers that subsequently fill-in with Fo–Ni-depleted olivine, causing reversals. Elemental diffusion modeling reveals growth of individual crystal frames and eruption at the surface occurred over 3.5–40 days. Those timescales constrain magma ascent rates of 40–500 m/h (0.011 to 0.14 m/s) from the deep crust. Compared to ocean island basalts, where dendritic and skeletal olivines have been often described, magmas erupted at arc settings, experiencing storage and degassing, may lack such textures due to fundamentally different ascent histories.

Arc olivines are commonly explained through a paradigm of core-to-rim sequential growth and oscillatory zoning is interpreted to represent magma mixing. Here the authors show Fo–Ni–P oscillatory zoned olivines can grow as out-of-sequence crystal frames and complex zoning can occur in closed systems.

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.

Insights on the deep carbon cycle from the electrical conductivity of carbon-bearing aqueous fluids

The dehydration and decarbonation in the subducting slab are intricately related and the knowledge of the physical properties of the resulting C-H-O fluid is crucial to interpret the petrological, geochemical, and geophysical processes associated with subduction zones. In this study, we investigate the C-H-O fluid released during the progressive devolatilization of carbonate-bearing serpentine-polymorph chrysotile, with in situ electrical conductivity measurements at high pressures and temperatures. The C-H-O fluid produced by carbonated chrysotile exhibits high electrical conductivity compared to carbon-free aqueous fluids and can be an excellent indicator of the migration of carbon in subduction zones. The crystallization of diamond and graphite indicates that the oxidized C-H-O fluids are responsible for the recycling of carbon in the wedge mantle. The carbonate and chrysotile bearing assemblages stabilize dolomite during the devolatilization process. This unique dolomite forming mechanism in chrysotile in subduction slabs may facilitate the transport of carbon into the deep mantle.

Titanium isotopes constrain a magmatic transition at the Hadean-Archean boundary in the Acasta Gneiss Complex

Plate subduction greatly influences the physical and chemical characteristics of Earth's surface and deep interior, yet the timing of its initiation is debated because of the paucity of exposed rocks from Earth's early history. We show that the titanium isotopic composition of orthogneisses from the Acasta Gneiss Complex spanning the Hadean to Eoarchean transition falls on two distinct magmatic differentiation trends. Hadean tonalitic gneisses show titanium isotopic compositions comparable to modern evolved tholeiitic magmas, formed by differentiation of dry parental magmas in plume settings. Younger Eoarchean granitoid gneisses have titanium isotopic compositions comparable to modern calc-alkaline magmas produced in convergent arcs. Our data therefore document a shift from tholeiitic- to calc-alkaline-style magmatism between 4.02 and 3.75 billion years (Ga) in the Slave craton.

Crustal thickening and endogenic oxidation of magmatic sulfur

Porphyry ore deposits, Earth's most important resources of copper, molybdenum, and rhenium, are strongly associated with felsic magmas showing signs of high-pressure differentiation and are usually found in places with thickened crust (>45 kilometers). This pattern is well-known, but unexplained, and remains an outstanding problem in our understanding of porphyry ore deposit formation. We approach this problem by investigating the oxidation state of magmatic sulfur, which controls the behavior of ore-forming metals during magma differentiation and magmatic-hydrothermal transition. We use sulfur in apatite to reconstruct the sulfur oxidation state in the Gangdese batholith, southern Tibet. We find that magma sulfate content increased abruptly after India-Eurasia collision. Apatite sulfur content and the calculated magma S6+/ΣS ratio correlate with whole-rock dysprosium/ytterbium ratio, suggesting that residual garnet, favored in thickened crust, exerts a first-order control on sulfur oxidation in magmatic orogens. Our findings link sulfur oxidation to internal petrogenic processes and imply an intrinsic relationship of magma oxidation with synmagmatic crustal thickening.

Geochemical evidence for a widespread mantle re-enrichment 3.2 billion years ago: implications for global-scale plate tectonics

Progressive mantle melting during the Earth’s earliest evolution led to the formation of a depleted mantle and a continental crust enriched in highly incompatible elements. Re-enrichment of Earth’s mantle can occur when continental crustal materials begin to founder into the mantle by either subduction or, to a lesser degree, by delamination processes, profoundly affecting the mantle’s trace element and volatile compositions. Deciphering when mantle re-enrichment/refertilization became a global-scale process would reveal the onset of efficient mass transfer of crust to the mantle and potentially when plate tectonic processes became operative on a global-scale. Here we document the onset of mantle re-enrichment/refertilization by comparing the abundances of petrogenetically significant isotopic values and key ratios of highly incompatible elements compared to lithophile elements in Archean to Early-Proterozoic mantle-derived melts (i.e., basalts and komatiites). Basalts and komatiites both record a rapid-change in mantle chemistry around 3.2 billion years ago (Ga) signifying a fundamental change in Earth geodynamics. This rapid-change is recorded in Nd isotopes and in key trace element ratios that reflect a fundamental shift in the balance between fluid-mobile and incompatible elements (i.e., Ba/La, Ba/Nb, U/Nb, Pb/Nd and Pb/Ce) in basaltic and komatiitic rocks. These geochemical proxies display a significant increase in magnitude and variability after ~3.2 Ga. We hypothesize that rapid increases in mantle heterogeneity indicate the recycling of supracrustal materials back into Earth’s mantle via subduction. Our new observations thus point to a ≥ 3.2 Ga onset of global subduction processes via plate tectonics.

The Role of Magma Mixing in Generating Granodioritic Intrusions Related to Cu–W Mineralization: A Case Study from Qiaomaishan Deposit, Eastern China

The newly exploited Qiaomaishan Cu−W deposit, located in the Xuancheng ore district in the MLYRB, is a middle-sized Cu–W skarn-type polymetallic deposit. As Cu–W mineralization is a rare and uncommon type in the Middle-Lower Yangtze River Belt (MLYRB), few studies have been carried out, and the geochemical characteristics and petrogenesis of Qiaomaishan intrusive rocks related to Cu–W mineralization are not well documented. We studied two types of ore-bearing intrusive rocks in the Qiaomaishan region, i.e., pure granodiorite porphyry and granodiorite porphyry with mafic microgranular enclaves (MMEs). Age characterization using zircon LA–ICP–MS showed that they were formed almost simultaneously, around 134.9 to 135.1 Ma. Granodiorite porphyries are high Mg# adakites, characterized by high-K calc-alkaline and metaluminous features that are enriched in LILEs (e.g., Sr and Ba) and LREEs, but depleted in HFSEs (e.g., Nb, Ta, and Ti) and HREEs. Moreover, they have enriched Sr–Nd–Hf isotopic compositions (with whole-rock (87Sr/86Sr)i ratios (0.706666−0.706714), negative εNd(t) values of −9.1 to −8.6, negative zircon εHf(t) values of −12.2 to −6.7, and two-stage Hf model ages (TDM2) between 1.5 and 2.0 Ga). However, compared to host rocks, the granodiorite porphyry with MMEs shows variable geochemical compositions, e.g., high Mg#, Cr, Ni, and V contents and enriched with LILEs. In addition, they have more depleted ISr, εNd(t), and εHf(t) values (0.706025 to 0.706269, −6.4 to −7.4, and −10.6 to −5.7, respectively), overlapping with regions of Early Cretaceous mafic rocks derived from enriched lithospheric mantle in the MLYRB. Coupled with significant disequilibrium textures and geochemical features of host rocks and MMEs, we propose that those rocks have resulted from mixing the felsic lower crust-derived magma and the mafic magma generated from the enriched mantle. The mixed magmas subsequently rose to shallow crust to form the ore-bearing rocks and facilitate Cu–W mineralization.

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.

The role of buoyancy in the fate of ultra-high-pressure eclogite

Eclogite facies metamorphism of the lithosphere forms dense mineral assemblages at high- (1.6–2.4 GPa) to ultra-high-pressure (>2.4–12 GPa: UHP) conditions that drive slab-pull forces during its subduction to lower mantle conditions. The relative densities of mantle and lithospheric components places theoretical limits for the re-exposure, and peak conditions expected, of subducted lithosphere. Exposed eclogite terranes dominated by rock denser than the upper mantle are problematic, as are interpretations of UHP conditions in buoyant rock types. Their subduction and exposure require processes that overcame predicted buoyancy forces. Phase equilibria modelling indicates that depths of 50–60 km (P = 1.4–1.8 GPa) and 85–160 km (P = 2.6–5 GPa) present thresholds for pull force in end-member oceanic and continental lithosphere, respectively. The point of no-return for subducted silicic crustal rocks is between 160 and 260 km (P = 5.5–9 GPa), limiting the likelihood of stishovite–wadeite–K-hollandite-bearing assemblages being preserved in equilibrated assemblages. The subduction of buoyant continental crust requires its anchoring to denser mafic and ultramafic lithosphere in ratios below 1:3 for the continental crust to reach depths of UHP conditions (85–160 km), and above 2:3 for it to reach extreme depths (>160 km). The buoyant escape of continental crust following its detachment from an anchored situation could carry minor proportions of other rocks that are denser than the upper mantle. However, instances of rocks returned from well-beyond these limits require exceptional exhumation dynamics, plausibly coupled with the effects of incomplete metamorphism to retain less dense low-P phases.

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.

Global atmospheric oxygen variations recorded by Th/U systematics of igneous rocks

Scientists have been working on the dating of atmospheric oxygenation in Earth’s history for decades. However, most previous studies relied on evidence from sedimentary rocks. Here, we show that igneous rocks can also be linked with surface oxidation by a key geological process: plate subduction. We here make an attempt to apply the Th/U of worldwide arc igneous rocks as an indicator for the timing of atmospheric oxygenation over the Earth’s history. The results are coincident with the previously defined Great Oxidation Event and Neoproterozoic Oxygenation Event.

Atmospheric oxygen has evolved from negligible levels in the Archean to the current level of about 21% through 2 major step rises: The Great Oxidation Event (GOE) in the early Proterozoic and the Neoproterozoic Oxygenation Event (NOE) during the late Proterozoic. However, most previous methods for constraining the time of atmospheric oxygenation have relied on evidence from sedimentary rocks. Here, we investigate the temporal variations of the Th/U of arc igneous rocks since 3.0 billion y ago (Ga) and show that 2 major Th/U decreases are recorded at ca. 2.35 Ga and ca. 0.75 Ga, coincident with the beginning of the GOE and NOE. The decoupling of U from Th is predominantly caused by the significant rise of atmospheric oxygen. Under an increasingly oxidized atmosphere condition, more uranium in the surface environment became oxidized from the water-insoluble U⁴⁺ to the water-soluble U⁶⁺ valance and incorporated in the sea water and altered oceanic crust. Eventually, the subduction of this altered oceanic crust produced the low-Th/U signature of arc igneous rocks. Therefore, the sharp decrease of Th/U in global arc igneous rocks may provide strong evidence for the rise of atmospheric oxygen. We suggest that the secular Th/U evolution of arc igneous rocks could be an effective geochemical indicator recording the global-scale atmospheric oxygen variation.

The Role of N2 as a Geo-Biosignature for the Detection and Characterization of Earth-like Habitats

Since the Archean, N₂ has been a major atmospheric constituent in Earth's atmosphere. Nitrogen is an essential element in the building blocks of life; therefore, the geobiological nitrogen cycle is a fundamental factor in the long-term evolution of both Earth and Earth-like exoplanets. We discuss the development of Earth's N₂ atmosphere since the planet's formation and its relation with the geobiological cycle. Then we suggest atmospheric evolution scenarios and their possible interaction with life-forms: first for a stagnant-lid anoxic world, second for a tectonically active anoxic world, and third for an oxidized tectonically active world. Furthermore, we discuss a possible demise of present Earth's biosphere and its effects on the atmosphere. Since life-forms are the most efficient means for recycling deposited nitrogen back into the atmosphere at present, they sustain its surface partial pressure at high levels. Also, the simultaneous presence of significant N₂ and O₂ is chemically incompatible in an atmosphere over geological timescales. Thus, we argue that an N₂-dominated atmosphere in combination with O₂ on Earth-like planets within circumstellar habitable zones can be considered as a geo-biosignature. Terrestrial planets with such atmospheres will have an operating tectonic regime connected with an aerobic biosphere, whereas other scenarios in most cases end up with a CO₂-dominated atmosphere. We conclude with implications for the search for life on Earth-like exoplanets inside the habitable zones of M to K stars.

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.

Importance of Magmatic Water Content and Oxidation State for Porphyry-Style Au Mineralization: An Example from the Giant Beiya Au Deposit, SW China

The Beiya Au deposit is the largest Cenozoic Au deposit in the Jinshajiang-Ailaoshan porphyry metallogenic belt. Numerous studies document that high water content and fO2 are vital factors for the generation of Au mineralization. In this belt, only the Wandongshan and Hongnitang districts are considered to be of economic importance, while the other districts, such as Bailiancun, are barren. So in order to reveal the importance of water content and oxidation state for Beiya porphyry-style Au mineralization, the amphiboles and zircons compositions are used to evaluate the physicochemical conditions (e.g., pressure, temperature, fO2, and water content) of the Wandongshan ore-fertile porphyries and Bailiancun ore-barren porphyries observed in the Beiya Au deposit. The results show that the water content of the Wandongshan parent magma (≤4.11 ± 0.4 wt %) are slightly higher than those of the parent magma at Bailiancun (≤3.91 ± 0.4 wt %), while the emplacement pressure of the Wandongshan parent magma (31.5–68.6 MPa) is much lower than that of the parent magma at Bailiancun (142.3–192.8 MPa), indicating that the Wandongshan magma reached water saturation earlier. In addition, the Wandongshan porphyries crystallized from more oxidized magma (average of ΔFMQ = +3.5) with an average temperature of 778 °C compared to the Bailiancun porphyries (average of ΔFMQ = +1.5) with a mean magmatic temperature of 770 °C. The Ce4+/Ce3+ ratio of zircon in the Wandongshan ore-related intrusions (average Ce4+/Ce3+ of 62.00) is much higher than that of the Bailiancun barren porphyries (average Ce4+/Ce3+of 23.15), which further confirmed Wandongshan ore-related magma is more oxidized than the Bailiancun barren magma. Therefore, melts that are more enriched in water and with a high oxidation state will be more fertile to form an economic porphyry-style Au system.

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

Subduction zone magmas are more oxidised on eruption than those at mid-ocean ridges. This is attributed either to oxidising components, derived from subducted lithosphere (slab) and added to the mantle wedge, or to oxidation processes occurring during magma ascent via differentiation. Here we provide direct evidence for contributions of oxidising slab agents to melts trapped in the sub-arc mantle. Measurements of sulfur (S) valence state in sub-arc mantle peridotites identify sulfate, both as crystalline anhydrite (CaSO₄) and dissolved SO₄^2- in spinel-hosted glass (formerly melt) inclusions. Copper-rich sulfide precipitates in the inclusions and increased Fe³⁺/∑Fe in spinel record a S⁶⁺-Fe²⁺ redox coupling during melt percolation through the sub-arc mantle. Sulfate-rich glass inclusions exhibit high U/Th, Pb/Ce, Sr/Nd and δ³⁴S (+ 7 to + 11‰), indicating the involvement of dehydration products of serpentinised slab rocks in their parental melt sources. These observations provide a link between liberated slab components and oxidised arc magmas.

Multi-level magma plumbing at Agung and Batur volcanoes increases risk of hazardous eruptions

The island of Bali in Indonesia is home to two active stratovolcanoes, Agung and Batur, but relatively little is known of their underlying magma plumbing systems. Here we define magma storage depths and isotopic evolution of the 1963 and 1974 eruptions using mineral-melt equilibrium thermobarometry and oxygen and helium isotopes in mineral separates. Olivine crystallised from a primitive magma and has average δ¹⁸O values of 4.8‰. Clinopyroxene records magma storage at the crust-mantle boundary, and displays mantle-like isotope values for Helium (8.62 R_A) and δ¹⁸O (5.0–5.8‰). Plagioclase reveals crystallisation in upper crustal storage reservoirs and shows δ¹⁸O values of 5.5–6.4‰. Our new thermobarometry and isotope data thus corroborate earlier seismic and InSAR studies that inferred upper crustal magma storage in the region. This type of multi-level plumbing architecture could drive replenishing magma to rapid volatile saturation, thus increasing the likelihood of explosive eruptions and the consequent hazard potential for the population of Bali.

The redox "filter" beneath magmatic orogens and the formation of continental crust

The two most important magmatic differentiation series on Earth are the Fe-enriching tholeiitic series, which dominates the oceanic crust and island arcs, and the Fe-depleting calc-alkaline series, which dominates the continental crust and continental arcs. It is well known that calc-alkaline magmas are more oxidized when they erupt and are preferentially found in regions of thick crust, but why these quantities should be related remains unexplained. We use the redox-sensitive behavior of europium (Eu) in deep-seated, plagioclase-free arc cumulates to directly constrain the redox evolution of arc magmas at depth. Primitive arc cumulates have negative Eu anomalies, which, in the absence of plagioclase, can only be explained by Eu being partly reduced. We show that primitive arc magmas begin with low oxygen fugacities, similar to that of mid-ocean ridge basalts, but increase in oxygen fugacity by over two orders of magnitude during magmatic differentiation. This intracrustal oxidation is attended by Fe depletion coupled with fractionation of Fe-rich garnet. We conclude that garnet fractionation, owing to its preference for ferrous over ferric iron, results in simultaneous oxidation and Fe depletion of the magma. Favored at high pressure and water content, garnet fractionation explains the correlation between crustal thickness, oxygen fugacity, and the calc-alkaline character of arc magmas.

Immiscible hydrous Fe-Ca-P melt and the origin of iron oxide-apatite ore deposits

The origin of iron oxide-apatite deposits is controversial. Silicate liquid immiscibility and separation of an iron-rich melt has been invoked, but Fe–Ca–P-rich and Si-poor melts similar in composition to the ore have never been observed in natural or synthetic magmatic systems. Here we report experiments on intermediate magmas that develop liquid immiscibility at 100 MPa, 1000–1040 °C, and oxygen fugacity conditions (fO₂) of ∆FMQ = 0.5–3.3 (FMQ = fayalite-magnetite-quartz equilibrium). Some of the immiscible melts are highly enriched in iron and phosphorous ± calcium, and strongly depleted in silicon (<5 wt.% SiO₂). These Si-poor melts are in equilibrium with a rhyolitic conjugate and are produced under oxidized conditions (~FMQ + 3.3), high water activity (aH₂O ≥ 0.7), and in fluorine-bearing systems (1 wt.%). Our results show that increasing aH₂O and fO₂ enlarges the two-liquid field thus allowing the Fe–Ca–P melt to separate easily from host silicic magma and produce iron oxide-apatite ores.

The origin of iron oxide-apatite deposits remains enigmatic and controversial. Here, the authors perform experiments on intermediate magmas and show that increasing aH₂O and fO₂ enlarges the two-liquid field thus allowing the Fe–Ca–P melt to separate easily from host silicic magma and produce iron oxide-apatite ores.

Redox-influenced seismic properties of upper-mantle olivine

Lateral variations of seismic wave speeds and attenuation (dissipation of strain energy) in the Earth's upper mantle have the potential to map key characteristics such as temperature, major-element composition, melt fraction and water content. The inversion of these data into meaningful representations of physical properties requires a robust understanding of the micromechanical processes that affect the propagation of seismic waves. Structurally bound water (hydroxyl) is believed to affect seismic properties but this has yet to be experimentally quantified. Here we present a comprehensive low-frequency forced-oscillation assessment of the seismic properties of olivine as a function of water content within the under-saturated regime that is relevant to the Earth's interior. Our results demonstrate that wave speeds and attenuation are in fact strikingly insensitive to water content. Rather, the redox conditions imposed by the choice of metal sleeving, and the associated defect chemistry, appear to have a substantial influence on the seismic properties. These findings suggest that elevated water contents are not responsible for low-velocity or high-attenuation structures in the upper mantle. Instead, the high attenuation observed in hydrous and oxidized regions of the upper mantle (such as above subduction zones) may reflect the prevailing oxygen fugacity. In addition, these data provide no support for the hypothesis whereby a sharp lithosphere-asthenosphere boundary is explained by enhanced grain boundary sliding in the presence of water.

A record of deep-ocean dissolved O2 from the oxidation state of iron in submarine basalts

The oxygenation of the deep ocean in the geological past has been associated with a rise in the partial pressure of atmospheric molecular oxygen (O₂) to near-present levels and the emergence of modern marine biogeochemical cycles. It has also been linked to the origination and diversification of early animals. It is generally thought that the deep ocean was largely anoxic from about 2,500 to 800 million years ago, with estimates of the occurrence of deep-ocean oxygenation and the linked increase in the partial pressure of atmospheric oxygen to levels sufficient for this oxygenation ranging from about 800 to 400 million years ago. Deep-ocean dissolved oxygen concentrations over this interval are typically estimated using geochemical signatures preserved in ancient continental shelf or slope sediments, which only indirectly reflect the geochemical state of the deep ocean. Here we present a record that more directly reflects deep-ocean oxygen concentrations, based on the ratio of Fe³⁺ to total Fe in hydrothermally altered basalts formed in ocean basins. Our data allow for quantitative estimates of deep-ocean dissolved oxygen concentrations from 3.5 billion years ago to 14 million years ago and suggest that deep-ocean oxygenation occurred in the Phanerozoic (541 million years ago to the present) and potentially not until the late Palaeozoic (less than 420 million years ago).

Evaluating the importance of metamorphism in the foundering of continental crust

The metamorphic conditions and mechanisms required to induce foundering in deep arc crust are assessed using an example of representative lower crust in SW New Zealand. Composite plutons of Cretaceous monzodiorite and gabbro were emplaced at ~1.2 and 1.8 GPa are parts of the Western Fiordland Orthogneiss (WFO); examples of the plutons are tectonically juxtaposed along a structure that excised ~25 km of crust. The 1.8 GPa Breaksea Orthogneiss includes suitably dense minor components (e.g. eclogite) capable of foundering at peak conditions. As the eclogite facies boundary has a positive dP/dT, cooling from supra-solidus conditions (T > 950 ºC) at high-P should be accompanied by omphacite and garnet growth. However, a high monzodioritic proportion and inefficient metamorphism in the Breaksea Orthogneiss resulted in its positive buoyancy and preservation. Metamorphic inefficiency and compositional relationships in the 1.2 GPa Malaspina Pluton meant it was never likely to have developed densities sufficiently high to founder. These relationships suggest that the deep arc crust must have primarily involved significant igneous accumulation of garnet–clinopyroxene (in proportions >75%). Crustal dismemberment with or without the development of extensional shear zones is proposed to have induced foundering of excised cumulate material at P > 1.2 GPa.

Highly oxidising fluids generated during serpentinite breakdown in subduction zones

Subduction zones facilitate chemical exchanges between Earth’s deep interior and volcanism that affects habitability of the surface environment. Lavas erupted at subduction zones are oxidized and release volatile species. These features may reflect a modification of the oxidation state of the sub-arc mantle by hydrous, oxidizing sulfate and/or carbonate-bearing fluids derived from subducting slabs. But the reason that the fluids are oxidizing has been unclear. Here we use theoretical chemical mass transfer calculations to predict the redox state of fluids generated during serpentinite dehydration. Specifically, the breakdown of antigorite to olivine, enstatite, and chlorite generates fluids with high oxygen fugacities, close to the hematite-magnetite buffer, that can contain significant amounts of sulfate. The migration of these fluids from the slab to the mantle wedge could therefore provide the oxidized source for the genesis of primary arc magmas that release gases to the atmosphere during volcanism. Our results also show that the evolution of oxygen fugacity in serpentinite during subduction is sensitive to the amount of sulfides and potentially metal alloys in bulk rock, possibly producing redox heterogeneities in subducting slabs.

Identification of an Early-Middle Jurassic oxidized magmatic belt, south Gangdese, Tibet, and geological implications

The south Gangdese region is the site of subduction of the Neo-Tethys and subsequent continental collision. Compared with widespread Cretaceous and Cenozoic magmatism, Early-Middle Jurassic magmatic rocks and related deposits are rarely reported. Our work identified a >200km long felsic rock belt associated with Cu mineralization in the south Gangdese region. We report here zircon U-Pb ages, zircon Ce⁴⁺/Ce³⁺ values, and mineral assemblages of two Cu mineralized intrusions within the belt. A hornblende granite and a diorite porphyry were emplaced at 177.3Ma and 166.3Ma, respectively. Geological occurrence and magmatic hematite-magnetite-chalcopyrite intergrowths suggest that Cu mineralization formed coeval with Jurassic intrusions. Mineralized intrusions have high zircon Ce⁴⁺/Ce³⁺ and Eu_N/Eu_N^∗ ratios, and hematite-magnetite intergrowths, suggesting their parent magmas were highly oxidized. Hornblende is common and primary fluid inclusions are found in titanite and apatite, indicating their parent magmas were water-saturated and exsolved volatile phases at early stage of magmatic evolution. Those magma characters contribute to the formation of porphyry Cu deposits. Given that majority subduction-related porphyry Cu systems have been eroded following uplift and denudation, the well-preserved Early-Middle Jurassic Cu mineralized igneous rocks in south Gangdese are favorable prospecting targets for subduction-related porphyry Cu deposits.

The geobiological nitrogen cycle: From microbes to the mantle

Nitrogen forms an integral part of the main building blocks of life, including DNA, RNA, and proteins. N₂ is the dominant gas in Earth's atmosphere, and nitrogen is stored in all of Earth's geological reservoirs, including the crust, the mantle, and the core. As such, nitrogen geochemistry is fundamental to the evolution of planet Earth and the life it supports. Despite the importance of nitrogen in the Earth system, large gaps remain in our knowledge of how the surface and deep nitrogen cycles have evolved over geologic time. Here, we discuss the current understanding (or lack thereof) for how the unique interaction of biological innovation, geodynamics, and mantle petrology has acted to regulate Earth's nitrogen cycle over geologic timescales. In particular, we explore how temporal variations in the external (biosphere and atmosphere) and internal (crust and mantle) nitrogen cycles could have regulated atmospheric pN₂ . We consider three potential scenarios for the evolution of the geobiological nitrogen cycle over Earth's history: two in which atmospheric pN₂ has changed unidirectionally (increased or decreased) over geologic time and one in which pN₂ could have taken a dramatic deflection following the Great Oxidation Event. It is impossible to discriminate between these scenarios with the currently available models and datasets. However, we are optimistic that this problem can be solved, following a sustained, open-minded, and multidisciplinary effort between surface and deep Earth communities.

Observing eruptions of gas-rich compressible magmas from space

Observations of volcanoes from space are a critical component of volcano monitoring, but we lack quantitative integrated models to interpret them. The atmospheric sulfur yields of eruptions are variable and not well correlated with eruption magnitude and for many eruptions the volume of erupted material is much greater than the subsurface volume change inferred from ground displacements. Up to now, these observations have been treated independently, but they are fundamentally linked. If magmas are vapour-saturated before eruption, bubbles cause the magma to become more compressible, resulting in muted ground displacements. The bubbles contain the sulfur-bearing vapour injected into the atmosphere during eruptions. Here we present a model that allows the inferred volume change of the reservoir and the sulfur mass loading to be predicted as a function of reservoir depth and the magma's oxidation state and volatile content, which is consistent with the array of natural data.

Satellite observations are an important tool in volcano monitoring, but observations such as ground deformation and gas emissions are treated independently. Here, the authors present a model coupling them through their link to magma volatile contents and storage depths prior to eruption

Zinc isotope evidence for sulfate-rich fluid transfer across subduction zones

Subduction zones modulate the chemical evolution of the Earth's mantle. Water and volatile elements in the slab are released as fluids into the mantle wedge and this process is widely considered to result in the oxidation of the sub-arc mantle. However, the chemical composition and speciation of these fluids, which is critical for the mobility of economically important elements, remain poorly constrained. Sulfur has the potential to act both as oxidizing agent and transport medium. Here we use zinc stable isotopes (δ⁶⁶Zn) in subducted Alpine serpentinites to decipher the chemical properties of slab-derived fluids. We show that the progressive decrease in δ⁶⁶Zn with metamorphic grade is correlated with a decrease in sulfur content. As existing theoretical work predicts that Zn-SO₄^2- complexes preferentially incorporate heavy δ⁶⁶Zn, our results provide strong evidence for the release of oxidized, sulfate-rich, slab serpentinite-derived fluids to the mantle wedge.

Sulfur and metal fertilization of the lower continental crust

Mantle-derived melts and metasomatic fluids are considered to be important in the transport and distribution of trace elements in the subcontinental lithospheric mantle. However, the mechanisms that facilitate sulfur and metal transfer from the upper mantle into the lower continental crust are poorly constrained. This study addresses this knowledge gap by examining a series of sulfide- and hydrous mineral-rich alkaline mafic-ultramafic pipes that intruded the lower continental crust of the Ivrea-Verbano Zone in the Italian Western Alps. The pipes are relatively small (< 300 m diameter) and primarily composed of a matrix of subhedral to anhedral amphibole (pargasite), phlogopite and orthopyroxene that enclose sub-centimeter-sized grains of olivine. The 1 to 5 m wide rim portions of the pipes locally contain significant blebby and disseminated Fe-Ni-Cu-PGE sulfide mineralization. Stratigraphic relationships, mineral chemistry, geochemical modelling and phase equilibria suggest that the pipes represent open-ended conduits within a large magmatic plumbing system. The earliest formed pipe rocks were olivine-rich cumulates that reacted with hydrous melts to produce orthopyroxene, amphibole and phlogopite. Sulfides precipitated as immiscible liquid droplets that were retained within a matrix of silicate crystals and scavenged metals from the percolating hydrous melt, associated with partial melting of a metasomatized continental lithospheric mantle. New high-precision chemical abrasion TIMS U-Pb dating of zircons from one of the pipes indicates that these pipes were emplaced at 249.1 ± 0.2 Ma, following partial melting of lithospheric mantle pods that were metasomatized during the Eo-Variscan oceanic to continental subduction (~420-310 Ma). The thermal energy required to generate partial melting of the metasomatized mantle was most likely derived from crustal extension, lithospheric decompression and subsequent asthenospheric rise during the orogenic collapse of the Variscan belt (< 300 Ma). Unlike previous models, outcomes from this study suggest a significant temporal gap between the occurrence of mantle metasomatism, subsequent partial melting and emplacement of the pipes. We argue that this multi-stage process is a very effective mechanism to fertilize the commonly dry and refractory lower continental crust in metals and volatiles. During the four-dimensional evolution of the thermo-tectonic architecture of any given terrain, metals and volatiles stored in the lower continental crust may become available as sources for subsequent ore-forming processes, thus enhancing the prospectivity of continental block margins for a wide range of mineral systems.

Seismic evidence of effects of water on melt transport in the Lau back-arc mantle

Processes of melt generation and transport beneath back-arc spreading centres are controlled by two endmember mechanisms: decompression melting similar to that at mid-ocean ridges and flux melting resembling that beneath arcs. The Lau Basin, with an abundance of spreading ridges at different distances from the subduction zone, provides an opportunity to distinguish the effects of these two different melting processes on magma production and crust formation. Here we present constraints on the three-dimensional distribution of partial melt inferred from seismic velocities obtained from Rayleigh wave tomography using land and ocean-bottom seismographs. Low seismic velocities beneath the Central Lau Spreading Centre and the northern Eastern Lau Spreading Centre extend deeper and westwards into the back-arc, suggesting that these spreading centres are fed by melting along upwelling zones from the west, and helping to explain geochemical differences with the Valu Fa Ridge to the south, which has no distinct deep low-seismic-velocity anomalies. A region of low S-wave velocity, interpreted as resulting from high melt content, is imaged in the mantle wedge beneath the Central Lau Spreading Centre and the northeastern Lau Basin, even where no active spreading centre currently exists. This low-seismic-velocity anomaly becomes weaker with distance southward along the Eastern Lau Spreading Centre and the Valu Fa Ridge, in contrast to the inferred increase in magmatic productivity. We propose that the anomaly variations result from changes in the efficiency of melt extraction, with the decrease in melt to the south correlating with increased fractional melting and higher water content in the magma. Water released from the slab may greatly reduce the melt viscosity or increase grain size, or both, thereby facilitating melt transport.

An experimental study of the partitioning of trace elements between rutile and silicate melt as a function of oxygen fugacity

Subduction zone or arc magmas are known to display a characteristic depletion of High Field Strength Elements (HFSE) relative to other similarly incompatible elements, which can be attributed to the presence of the accessory mineral rutile (TiO2) in the residual slab. Here we show that the partitioning behavior of vanadium between rutile and silicate melt varies from incompatible (∼0.1) to compatible (∼18) as a function of oxygen fugacity. We also confirm that the HFSE are compatible in rutile, with D(Ta)> D(Nb)>> (D(Hf)>/∼ D(Zr), but that the level of compatibility is strongly dependent on melt composition, with partition coefficients increasing about one order of magnitude with increasing melt polymerization (or decreasing basicity). Our partitioning results also indicate that residual rutile may fractionate U from Th due to the contrasting (over 2 orders of magnitude) partitioning between these two elements. We confirm that, in addition to the HFSE, Cr, Cu, Zn and W are compatible in rutile at all oxygen fugacity conditions.

Effects of Eyjafjallajökull volcanic ash on innate immune system responses and bacterial growth in vitro

BACKGROUND

On 20 March 2010, the Icelandic volcano Eyjafjallajökull erupted for the first time in 190 years. Despite many epidemiological reports showing effects of volcanic ash on the respiratory system, there are limited data evaluating cellular mechanisms involved in the response to ash. Epidemiological studies have observed an increase in respiratory infections in subjects and populations exposed to volcanic eruptions.

METHODS

We physicochemically characterized volcanic ash, finding various sizes of particles, as well as the presence of several transition metals, including iron. We examined the effect of Eyjafjallajökull ash on primary rat alveolar epithelial cells and human airway epithelial cells (20-100 µg/cm(2)), primary rat and human alveolar macrophages (5-20 µg/cm(2)), and Pseudomonas aeruginosa (PAO1) growth (3 µg/104 bacteria).

RESULTS

Volcanic ash had minimal effect on alveolar and airway epithelial cell integrity. In alveolar macrophages, volcanic ash disrupted pathogen-killing and inflammatory responses. In in vitro bacterial growth models, volcanic ash increased bacterial replication and decreased bacterial killing by antimicrobial peptides.

CONCLUSIONS

These results provide potential biological plausibility for epidemiological data that show an association between air pollution exposure and the development of respiratory infections. These data suggest that volcanic ash exposure, while not seriously compromising lung cell function, may be able to impair innate immunity responses in exposed individuals.