Elsevier

Earth and Planetary Science Letters

Volume 432, 15 December 2015, Pages 142-151
Earth and Planetary Science Letters

Redox-variability and controls in subduction zones from an iron-isotope perspective

https://doi.org/10.1016/j.epsl.2015.09.036Get rights and content

Highlights

  • Subducting sediment has no influence on Fe isotope compositions of arc lavas.

  • Fe isotopes are decoupled from mantle redox but record mantle depletion.

  • Depleted arc sources are isotopically lighter than fertile mantle.

  • An oxidant other than Fe is required for oxidised arc lavas.

  • Fe isotopes in arc lavas are controlled by source depletion and crystal fractionation.

Abstract

An ongoing controversy in mantle geochemistry concerns the oxidation state of the sources of island arc lavas (IAL). Three key factors control oxidation–reduction (redox) of IAL sources: (i) metasomatism of the mantle wedge by fluids and/or melts, liberated from the underlying subducted slab; (ii) the oxidation state of the wedge prior to melting and metasomatism; and (iii) the loss of melt from IAL sources. Subsequently, magmatic differentiation by fractional crystallisation, possible crustal contamination and degassing of melts en route to and at the surface can further modify the redox states of IAL. The remote nature of sub-arc processes and the complex interplay between them render direct investigations difficult. However, a possible gauge for redox-controlled, high-temperature pre-eruptive differentiation conditions is variations in stable Fe isotope compositions (expressed here as δ57Fe) in erupting IAL because Fe isotopes can preserve a record of sub-surface mass transfer reactions involving the major element Fe. Here we report Fe isotope compositions of bulk IAL along the active Banda arc, Indonesia, which is well known for a prominent subducted sediment input. In conjunction with other arc rocks, δ57Fe in erupted Banda IAL indicates that fractional crystallisation and possibly crustal contamination primarily control their Fe isotope signatures. When corrected for fractional crystallisation and filtered for contamination, arc magmas that had variable sediment melt contributions in their sources show no resolvable co-variation of δ57Fe with radiogenic isotope tracers. This indicates that crustal recycling in the form of subducted sediment does not alter the Fe isotope character of arc lavas, in agreement with mass balance estimates. Primitive sources of IAL, however, are clearly isotopically lighter than those sourced beneath mid-ocean ridges, indicating either preferential Fe3+-depletion in the mantle wedge by prior, δ57Fe-heavy melt extraction, and/or addition of an isotopically-light slab-derived agent. Based on our findings and previous models of arc redox conditions, we propose a three-stage process to explain the Fe isotope composition of IAL: (i) prior melt depletion lowers Fe3+/ΣFe (Fe3+ over total Fe) in the residues, leaving refractory, δ57Fe-light and possibly reduced mantle wedge protoliths beneath arcs. The oxygen fugacity (fO2) of these refractory protoliths may be up to −2 log10 units reduced relative to the fayalite–magnetite–quartz synthetic oxygen buffer (ΔFMQ); (ii) oxidised, slab-derived fluids, Fe-poor but possibly rich in sulphate (i.e., S6+), trigger melting of depleted protoliths with minimal effect on δ57Fe. Melts derived from this fluid-modified wedge source have high Fe3+/ΣFe, oxidised by the reduction of S6+, but importantly retain the light δ57Fe from their mantle wedge source; (iii) after melt liberation from the mantle wedge, arc magmas initially become progressively oxidised and isotopically heavier in Fe through fractional crystallisation of ferromagnesian silicates. In summary, reduction consequent to Fe3+-rich melt extraction and subsequent oxidation, likely by S6+-rich fluids, results in a “redox yo-yo” in IAL sources. Fractional crystallisation will further oxidise and elevate δ57Fe in erupting IAL. Iron isotope signatures in IAL record a complex magmatic history with no simple relation between δ57Fe and calculated fO2 in erupted lavas. Records of higher fO2 in subduction zones compared to MORB sources are inherited from the subduction component.

Introduction

Lavas at convergent oceanic margins (hereafter referred to as island arc lavas; IAL) are chemically distinct from mid-ocean ridge basalts (MORB) particularly with respect to lithophile trace element abundances. Most IAL are sourced from the mantle by fluxed melting in the wedge (i.e., sub-arc asthenosphere) between the down-going and the overriding tectonic plates. In addition to the type of dry, decompressional melting occurring beneath MOR, a major fraction of melting in the mantle wedge is triggered by liquids sourced from the subducting plate. Fluids that cause melting in the wedge are liberated by dehydration reactions in the slab with the breakdown of hydrous phases hosted by altered oceanic crust. These fluids carry the so-called “subduction component” into the mantle wedge and upwards to the surface in arc melts (Arculus and Powell, 1986, Elliot et al., 1997, Gamble et al., 1996, Stolper and Newman, 1994, Thirlwall et al., 1996, Woodhead et al., 1998). The compositions and physical characteristics of IAL are consequent to the contribution of subduction component(s) and partial melting processes within the mantle wedge.

Element partitioning between residual minerals and the liberated slab-derived liquids at variable pressure–temperature (P–T) conditions in the downgoing slab define the compositional range of the subduction component (Pearce and Parkinson, 1993). For elements with multiple valence states, this partitioning is strongly dependent on the redox state of the system, reflected in the precisely calculable parameter of oxygen fugacity (fO2, Canil, 2002, Lee et al., 2005, Mallmann and O'Neill, 2007, McCanta et al., 2004). Determination of fO2 in the mantle wedge through studies of IAL have proved difficult due to weakly constrained processes of partial melting, fluid–rock interactions, and magmatic differentiation in the mantle wedge, together with compositional variations in subduction components and/or the local mantle wedge. However, despite these factors, it is well known that, on average, primitive IAL are more oxidised than MORB (see a recent review by Evans, 2012), yet the sources and locus of this oxidation process remain elusive (e.g., Lee et al., 2005).

Key to the redox history of igneous systems is the oxidation state of Fe (commonly expressed as the ratio of Fe3+ over total Fe, or ΣFe, i.e., Fe3+ + Fe2+), which is typically the predominant variable valence element in basaltic–andesitic igneous systems. Determinations of Fe3+/ΣFe in primitive IAL reveal an elevated oxidation state (expressed in log10fO2 units relative to equilibrium in the synthetic oxygen buffer system of fayalite–magnetite–quartz, i.e., ΔFMQ) of IAL, in the order of 1–4 log10 units above FMQ (Arculus, 1994, Carmichael, 1991, Osborn, 1959), compared to MORB with ΔFMQ+0.5 (Cottrell and Kelley, 2011).

In understanding the oxidation state of the sub-arc mantle wedge, three opposing scenarios have been suggested. Parkinson and Arculus (1999) argued for a more oxidised mantle wedge in comparison to the mantle sources of MORB. Their calculated fO2 for the mantle wedge from olivine–orthopyroxene-spinel thermobarometry of peridotites sourced from active island arcs, correspond to ΔFMQ=+0.5 to +1.7. Equally high estimates (ΔFMQ=+1 to +4) are deduced from spinel compositions in primitive IAL (Evans et al., 2012). Frost and Ballhaus (1998) and Parkinson and Arculus (1999) argued H2O (as part of the subduction component) is unlikely to be a mantle wedge oxidiser, the latter nevertheless suggesting a “solute-rich hydrous fluid” is most likely involved. By means of mass balance calculations of potential oxidising species in subduction zones, Evans (2012) suggested S6+, Fe3+ or CO2−3 could be involved. Although the evidence from arc peridotites suggests an oxidised mantle wedge, the similarity of Fe/Zn and V/Sc between IAL and MORB, noting the higher compatibility of oxidised Fe and V in a melt phase compared with residual solids relative to Zn and Sc respectively, argues for the same oxidation state of respective IAL and MORB sources (Lee et al., 2010, Lee et al., 2005, Mallmann and O'Neill, 2009).

Recent analyses of Fe isotope compositions in IAL are potentially key in this debate. Iron isotopes are considered to fractionate as a function of Fe redox state during partitioning of Fe between liquid and crystalline residues, e.g., during partial melting or by fractional crystallisation (Polyakov and Mineev, 2000, Shahar et al., 2008, Young et al., 2015), generally with the heavier Fe isotope following Fe3+. However, primitive arc melts are characterised by light Fe isotope compositions (Dauphas et al., 2009, Nebel et al., 2013) and may thus, counterintuitively, indicate a Fe3+-depleted (by prior melt extraction), or a Fe2+-enriched (by metasomatism) wedge source. Hence, our present collective interpretation of the sum of geochemical proxies allows for a more oxidised, equally oxidised or reduced IAL source mantle compared with MORB.

Iron isotopic characteristics of IAL are potentially important for resolving several problems regarding primitive arc magma genesis. These are the extents of prior melt depletion in the arc source(s) compared with those of MOR; whether transfer of oxidised Fe from slab-derived subducted sediments to the wedge is a significant redox process; and the influence of processes such as assimilation–fractional crystallisation (AFC), in modifying Fe isotope characteristics during magmatic evolution. In addressing these issues, we present Fe isotope compositions for a series of Holocene IAL and subducting sediments from the active Banda arc in Indonesia. Banda is a prime example wherein subducted sediment has clearly contributed towards the trace element and isotopic characteristics of the arc magmas (Nebel et al., 2011, Vroon et al., 2001, Vroon et al., 1995). We aim to simultaneously investigate the role of Fe in subducting sediment as an oxidiser in the mantle wedge, and to test whether Fe isotope compositions of IAL are reflective of the redox-character of their mantle wedge sources.

Section snippets

Geologic framework

The Banda subduction zone extends northeast of Timor, and exhibits an almost 180° curvature from south to north enclosing the Banda Sea (Fig. 1) with a radius of curvature of ca. 200 km (Sandiford, 2010, Spakman and Hall, 2010). The arc formed by the collision of the Australian/Indian plate with the Banda Sea microplate (Bird, 2003) and subduction of the former. Along-arc elemental and radiogenic isotope systematics reveal a contribution from subducted sediments to the arc lavas, with

Analytical methods

Sample preparation, including crushing, milling, and de-carbonatisation of sediments has been described in previous publications in detail (Vroon et al., 1995, Vroon et al., 1993). For Fe isotope analyses, approximately 25 mg of sample powder was dissolved in 5 ml Savillex vials in a mixture of concentrated HF–HNO3 for 48 h on a hotplate at 120 °C. After evaporation, all samples were equilibrated with HCl and re-dissolved in 1 ml of 9 M HCl. Chemical purification was performed with a DOWEX

Results

The Fe isotope compositions of the Banda arc rocks and sediments are reported in Table 1 and illustrated in Fig. 2. Banda arc rock isotope compositions vary from δ57Fe = +0.09 to δ57Fe = +0.23, overlapping in part the values observed for other island arcs (Dauphas et al., 2009) or MORB (Teng et al., 2013). Sample RO2C6 (Island of Romang) is an exception with an extremely heavy Fe isotope composition with δ57Fe = +1.00.

The Fe isotope compositions of Banda arc sediments show a range from δ57Fe =

The effects of fractional crystallisation

With one exception, the SiO2 and MgO contents of the analysed Banda arc samples range from 55.3 to 59.9 wt.%, and from 5.5 to 2.6 wt.% respectively; the exception is RO2C6 with SiO2 = 71.6 wt.%, and 0.31 wt.% MgO. All samples have experienced olivine and pyroxene fractional crystallisation from a lower SiO2–higher MgO basaltic parent prior to eruption. To account for the effects of Fe isotope fractionation caused by fractional crystallisation (Schuessler et al., 2007, Sossi et al., 2012, Teng

Conclusions

The Fe isotope compositions of Banda and other arc lavas imply their mantle sources are isotopically lighter than those sourcing mid-ocean ridges. The light Fe isotopic signature of primitive arc melts, irrespective of their actual oxidation state, is a consequence of prior episodes of extensive melt-depletion of their protoliths in the sub-arc mantle wedge. These protoliths were possibly further isotopically affected by metasomatic fluids. Subducting sediments are shown to leave no detectable

Acknowledgments

The project was funded by the Australian Research Council by a DECRA fellowship to ON (DE120100513), who also thanks Steve Eggins for access to clean room facilities and mass spectrometers. One anonymous reviewer and Helen Williams are acknowledged for their very constructive and insightful comments. Tamsin Mather is thanked for her additional comments and efficient editorial handling. A special thank you from O.N. to the Canberra T.I.E. team.

References (87)

  • K.A. Evans et al.

    Insights into subduction zone sulfur recycling from isotopic analysis of eclogite-hosted sulfides

    Chem. Geol.

    (2014)
  • A.J. Frierdich et al.

    Iron isotope fractionation between aqueous Fe(II) and goethite revisited: new insights based on a multi-direction approach to equilibrium and isotopic exchange rate modification

    Geochim. Cosmochim. Acta

    (2014)
  • J. Hermann et al.

    Experimental constraints on high pressure melting in subducted crust

    Earth Planet. Sci. Lett.

    (2001)
  • O. Ishizuka

    Early stages in the evolution of Izu–Bonin arc volcanism: new age, chemical, and isotopic constraints

    Earth Planet. Sci. Lett.

    (2006)
  • K. Klimm et al.

    The dissolution mechanism of sulphur in hydrous silicate melts. II: Solubility and speciation of sulphur in hydrous silicate melts as a function of fO(2)

    Chem. Geol.

    (2012)
  • C.A. Macris et al.

    Crystal chemical constraints on inter-mineral Fe isotope fractionation and implications for Fe isotope disequilibrium in San Carlos mantle xenoliths

    Geochim. Cosmochim. Acta

    (2015)
  • G. Mallmann et al.

    The effect of oxygen fugacity on the partitioning of Re between crystals and silicate melt during mantle melting

    Geochim. Cosmochim. Acta

    (2007)
  • M.T. McCulloch et al.

    Geochemical and geodynamic constraints on subduction zone magmatism

    Earth Planet. Sci. Lett.

    (1991)
  • O. Nebel et al.

    Hafnium and iron isotopes in early Archean komatiites record a plume-driven convection cycle in the Hadean Earth

    Earth Planet. Sci. Lett.

    (2014)
  • O. Nebel et al.

    Tungsten isotopes as tracers of core–mantle interactions: the influence of subducted sediments

    Geochim. Cosmochim. Acta

    (2010)
  • I.J. Parkinson et al.

    The redox state of subduction zones: insights from arc-peridotites

    Chem. Geol.

    (1999)
  • V.B. Polyakov et al.

    The use of Mossbauer spectroscopy in stable isotope geochemistry

    Geochim. Cosmochim. Acta

    (2000)
  • M.W. Schmidt

    Melting of pelitic sediments at subarc depths: 2. Melt chemistry, viscosities and a parameterization of melt composition

    Chem. Geol.

    (2015)
  • M.W. Schmidt et al.

    Experimentally based water budgets for dehydrating slabs and consequences for arc magma generation

    Earth Planet. Sci. Lett.

    (1998)
  • A. Shahar et al.

    Equilibrium high-temperature Fe isotope fractionation between fayalite and magnetite: an experimental calibration

    Earth Planet. Sci. Lett.

    (2008)
  • Z.D. Sharp et al.

    Water-soluble chlorides in massive seafloor serpentinites: a source of chloride in subduction zones

    Earth Planet. Sci. Lett.

    (2004)
  • D.M. Shaw

    Trace element fractionation during anatexis

    Geochim. Cosmochim. Acta

    (1970)
  • C. Spandler et al.

    Geochemical heterogeneity and element mobility in deeply subducted oceanic crust; insights from high-pressure mafic rocks from New Caledonia

    Chem. Geol.

    (2004)
  • A. Spicak et al.

    Seismic response to recent tectonic processes in the Banda Arc region

    J. Asian Earth Sci.

    (2013)
  • E. Stolper et al.

    The role of water in the petrogenesis of Mariana trough magmas

    Earth Planet. Sci. Lett.

    (1994)
  • B. Taylor et al.

    Back-arc basin basalt systematics

    Earth Planet. Sci. Lett.

    (2003)
  • F.-Z. Teng et al.

    Iron isotopic systematics of oceanic basalts

    Geochim. Cosmochim. Acta

    (2013)
  • M.F. Thirlwall et al.

    Resolution of the effects of crustal assimilation, sediment subduction, and fluid transport in island arc magmas: Pb–Sr–Nd–O isotope geochemistry of Grenada, Lesser Antilles

    Geochim. Cosmochim. Acta

    (1996)
  • A.G. Tomkins et al.

    Separate zones of sulfate and sulfide release from subducted mafic oceanic crust

    Earth Planet. Sci. Lett.

    (2015)
  • P.Z. Vroon et al.

    Oxygen isotope systematics of the Banda Arc: low delta O-18 despite involvement of subducted continental material in magma genesis

    Geochim. Cosmochim. Acta

    (2001)
  • P.Z. Vroon et al.

    Strontium, neodymium, and lead isotopic and trace-element signatures of the East Indonesian sediments – provenance and implications for Banda Arc magma genesis

    Geochim. Cosmochim. Acta

    (1995)
  • S. Weyer et al.

    Partial melting and melt percolation in the mantle: the message from Fe isotopes

    Earth Planet. Sci. Lett.

    (2007)
  • S. Weyer

    Fe-isotope fractionation during partial melting on Earth and the current view on the Fe-isotope budgets of the planets – reply to the comment of F. Poitrasson and to the comment of B.L. Beard and C.M. Johnson on “Iron isotope fractionation during planetary differentiation” by S. Weyer, A.D. Anbar, G.P. Brey, C. Munker, K. Mezger and A.B. Woodland

    Earth Planet. Sci. Lett.

    (2007)
  • H.M. Williams et al.

    Iron isotope tracing of mantle heterogeneity within the source regions of oceanic basalts

    Earth Planet. Sci. Lett.

    (2014)
  • H.M. Williams

    Systematic iron isotope variations in mantle rocks and minerals: the effects of partial melting and oxygen fugacity

    Earth Planet. Sci. Lett.

    (2005)
  • H.M. Williams

    Fractionation of oxygen and iron isotopes by partial melting processes: implications for the interpretation of stable isotope signatures in mafic rocks

    Earth Planet. Sci. Lett.

    (2009)
  • J. Woodhead et al.

    High-field strength element systematics in island arc basalts – evidence for multiphase melt extraction and a depleted mantle wedge Earth and

    Planet. Sci. Lett.

    (1993)
  • R.K. Workman et al.

    Major and trace element composition of the depleted MORB mantle (DMM)

    Earth Planet. Sci. Lett.

    (2005)
  • Cited by (82)

    View all citing articles on Scopus
    1

    Present address: Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, CNRS, F-75005 Paris, France.

    View full text