New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Neoproterozoic to early Phanerozoic rise in island arc redox state due to deep ocean oxygenation and increased marine sulfate levels
Edited by Mark H. Thiemens, University of California, San Diego, La Jolla, CA, and approved March 13, 2019 (received for review December 21, 2018)

Significance
Igneous island arc rocks are more oxidized than midocean-ridge basalts. Whether this originates from an oxidized mantle source or during differentiation or eruption is debated. Newly compiled Fe3+/ΣFe and V/Sc ratios presented here indicate that island arc rocks became more oxidized 800–400 Ma. We attribute this increase in the redox state of island arc rocks to the oxidation of their mantle source from subduction of oxidized oceanic crust following deep-ocean oxygenation and increased marine sulfate concentrations over the same time interval. This provides evidence that the oxidized nature of modern island arc magmas is due to an oxidized mantle source and represents a rare example of a change in the surface biogeochemistry influencing the igneous rock record.
Abstract
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.
Sometime between 800 and 400 Ma, atmospheric O2 concentrations are thought to have risen to sufficient concentrations (∼10 to 50% of modern levels) to oxygenate the majority of the deep oceans (e.g., refs. 1⇓⇓⇓⇓⇓⇓⇓⇓–10). This transition has been of long-standing interest as it marks a fundamental shift in the biogeochemical cycles on Earth’s surface to those similar to what is seen today and has been linked to the proliferation of animals (e.g., refs. 1, 2, 5, 6, 11, and 12). Additionally, over this same time interval marine sulfate concentrations are thought to have increased by as much as ∼10-fold and reached near-modern levels (e.g., refs. 13⇓⇓⇓–17). Increases in marine sulfate and deep-ocean O2 concentrations at this time have been proposed to have led to an increase in the redox state of oceanic crust via oxidation during hydrothermal circulation of oxygenated water (10) and/or circulation of seawater-derived sulfate (18, 19) through this crust. Before the oxygenation of the deep ocean and the rise in seawater sulfate levels these oxidative processes are likely to have been significantly diminished.
Igneous rocks from modern arc environments are well known to be oxidized relative to midocean-ridge basalts (MORB) (e.g., refs. 20⇓⇓⇓–24). Understanding the origin of this difference has been of enduring interest as the elevated redox state, often reported in terms of oxygen fugacity (fO2), of arc basalts relative to MORB has been proposed to control their differing mineral assemblages and thus differing trends in chemical composition (e.g., ref. 20), as well as the formation of arc-associated economic gold and copper deposits (e.g., ref. 25). One leading hypothesis for the origin of the elevated fO2 values of arc rocks vs. MORB is that oxidized material derived from subducting slabs oxidizes the mantle source of arc magmas (i.e., the subarc mantle) (18, 21, 24, 26⇓⇓⇓–30). Based on this, it has been hypothesized that an increase in the redox state of subducting slabs due to the oxygenation of the deep ocean and increased marine sulfate concentrations could be reflected by an elevation in the fO2 of igneous arc rocks in the Phanerozoic relative to in the Precambrian (10, 18, 19). However, whether the fO2 of arc rocks has or has not changed with time is neither known nor has it been directly evaluated.
Note that from this point on we use “redox state” to describe the relative balance of oxidized vs. reduced elements (e.g., Fe3+ vs. Fe2+ or S6+ vs. S2−) in any rock, including both pristine igneous rocks and those that may have experienced secondary (i.e., subsolidus) alteration (e.g., oxidation by oxygenated waters). We use oxygen fugacity (fO2), which is the fugacity of O2 in a magma or rock assemblage at a specified temperature and pressure, only for the description of the redox state of what are considered primary (i.e., unaltered) igneous mineral assemblages and melts. Finally, when fO2 values of differing rock types (e.g., island arc rocks vs. MORB) are compared this is always done implicitly relative to an igneous redox buffer [e.g., quartz–fayalite–magnetite (QFM)].
To test this hypothesis, we compiled two geochemical parameters for island arc rocks through time that have been widely used to estimate the fO2 of igneous rocks and their mantle sources, specifically Fe3+/ΣFe (i.e., molar Fe3+/[Fe2++Fe3+]) (e.g., refs. 24 and 31) and V/Sc ratios (32⇓⇓–35). Both ratios positively covary with source rock fO2. V/Sc ratios were chosen over other proposed trace-element-ratio oxybarometers such Zn/Fe (36), V/Ga (34), and V/Yb (37) as V/Sc is the only one that has been widely applied to quantitatively reconstruct the fO2 of modern basalts (including island arcs and MORB) and ancient (i.e., Precambrian) igneous rocks (32⇓⇓–35, 38). Island arc rocks were chosen for this study over continental arc rocks as their parental magmas neither transit through nor crystallize in older continental crust, which, if assimilated, could change their redox state.
Our proposed test is simple: (i) If the oxygenation of the deep ocean and concomitant increase in marine sulfate concentration caused a significant increase in the redox state of oceanic crustal material subducted into the mantle and (ii) this oxidized material influences the fO2 of the mantle source of island arc magmas, then it follows that both Fe3+/ΣFe and V/Sc ratios of igneous rocks from island arc environments should increase around the time of the Neoproterozoic–Phanerozoic boundary (∼540 Ma).
This test may fail for a variety of reasons. First, Fe3+/ΣFe ratios are susceptible to modification during melt differentiation, degassing, metamorphism, and exposure to the atmosphere, and such processes could obscure any primary information on the redox state of the source of island arcs. Second, serpentinization of ultramafic rocks also leads to the oxidation of oceanic crust and the upper lithosphere that is transported to the mantle during subduction and, therefore, could also lead to the oxidation of the subarc mantle (e.g., ref. 39). Although the degree of serpentinization of oceanic lithosphere is uncertain (one order of magnitude uncertainty; ref. 40), this process could have led to a sufficient flux of subducted oxidized material to oxidize island arc mantle sources since the initiation of subduction. Importantly, this would allow for the oxidation of the subarc mantle without invoking the involvement of oceanic crust oxidized by marine dissolved O2 or sulfate. Third, it has been proposed that the fO2 of the mantle sources of island arc rocks and MORB are similar (33, 34, 36, 41) or even that some primary arc melts can be more reduced than MORB (42). If this is correct, then the processes that cause the elevated redox state of island arc volcanic rocks must occur instead during ascent, melt differentiation, and/or eruption (33, 34, 36, 41, 43). If island arc rocks are more oxidized vs. MORB due to the subduction of oxidized serpentinite or due to melt differentiation or eruption, then no temporal change in the Fe3+/ΣFe or V/Sc ratios of island arc rocks would be predicted—for example, being uniformly elevated in island arc igneous rocks relative to MORB at all ages in the first case or uniform and similar in primitive, undegassed island arc igneous rocks (or even reduced) relative to primitive MORB at all ages in the second case.
Materials and Methods
We compiled Fe3+/ΣFe and V/Sc ratios of both extrusive and intrusive rocks previously interpreted to have formed in island arc environments from 3.1 billion years ago to the present. Implicitly, we have made the common assumption that plate tectonics and thus subduction began in the Archean, although the timing continues to be debated (see discussion in refs. 44 and 45). Island arc localities were identified based on trace-element arguments (e.g., refs. 46⇓–48), field relationships and lithologic associations (e.g., volcaniclastic sedimentary rocks and paired metamorphic belts), presence of boninites (e.g., refs. 49 and 50), geochronology supporting a tectonic history as an island arc, and prior compilations (e.g., refs. 51 and 52)—all localities have previously and independently interpreted as ancient island arcs. Details of each locality are given in Dataset S1, including which of the above criteria were fulfilled for each locality. In general, most localities had at least two of the above criteria met; however, identification of Paleoproterozoic and Archean localities was primarily based on trace element arguments. We note that Dataset S1 is not meant to represent a comprehensive list of all island arc localities in the geologic record, but rather those for which whole-rock Fe3+/ΣFe and/or V/Sc ratios were available.
Insofar as it was possible, cumulates were avoided for both records. Metamorphic grade for the localities, when reported, varied from unmetamorphosed to amphibolite facies and did not vary systematically with age (Dataset S1). Fe3+/ΣFe values are based on whole-rock wet chemical extractions (all data in Dataset S2 and locality averages in Dataset S3). For Fe3+/ΣFe ratios, subaerially erupted rocks were excluded (again insofar as it was possible) to avoid samples potentially oxidized by atmospheric O2 during eruption following the initial oxygenation of the atmosphere ∼2.4 to 2.3 billion years ago (e.g., refs. 5 and 11).
For the V/Sc compilation (all data in Dataset S4 and locality averages in Dataset S5), rocks were avoided that showed obvious evidence of clinopyroxene crystallization, which modifies V/Sc ratios (e.g., refs. 32 and 53). Following previous workers (32⇓⇓–35, 53), we used a MgO weight percent cutoff to select samples derived from melts that did not undergo clinopyroxene crystallization. We chose a value of 6.5 wt % MgO (with samples needing to have >6.5 wt % MgO) based on an examination of where CaO begins to decrease when CaO vs. MgO is plotted for precompiled island arc igneous rocks from the GEOROC database (georoc.mpch-mainz.gwdg.de/georoc/Entry.html). This break in slope marks the onset of clinopyroxene crystallization and typically occurs in this database between 5 and 6.5 wt % MgO. Our MgO cutoff is lower than that used previously for MORB (8 wt % MgO; ref. 32) as the presence of water in island arc parental melts expands the stability field of olivine relative to pyroxene (e.g., ref. 54) but is consistent with previous MgO cutoffs used in studies of V and Sc in island arc rocks (5 to 6 wt %; refs. 34 and 53). We also excluded rocks with >15 wt % MgO to avoid incorporation of samples derived from exceptionally high degrees of melt or that represent a mixture of melt and accumulated olivine (e.g., ref. 32). As discussed below, the precise choice of the MgO cutoff does not affect any of our interpretations. We also present previously compiled MORB V/Sc ratios from Gale et al. (55) for comparison as well as compilations of V/Sc for ancient basalts used previously to reconstruct the fO2 of the upper mantle (32, 35). For these, as was the case for our treatment of MORB, we again selected samples with MgO contents between 8 and 15 wt %.
Results and Discussion
Fe3+/ΣFe Ratios of Island Arc Igneous Rocks Through Time.
Fe3+/ΣFe ratios of igneous rocks and melt inclusions have been used extensively to reconstruct the fO2 of the parental melts of young igneous rocks, for example MORB dredged from the seafloor or samples taken from active island arcs (e.g., refs. 31, 22, 56, 24, and 57⇓⇓–60). We use this approach on whole rocks to attempt to reconstruct the redox state of island arc rocks through time. Glasses are preferred for Fe3+/ΣFe–based fO2 reconstructions as they are quenched melts. In using whole rocks, we follow previous workers in assuming that examined whole rocks are representative of the bulk composition of the melt from which they crystallized (e.g., refs. 56 and 60).
Use of whole-rock Fe3+/ΣFe ratios from ancient (i.e., Precambrian) samples for the purpose of reconstructing the original redox state of a sample is uncommon due to the typical assumption that Fe3+/ΣFe ratios of ancient samples may be compromised by metamorphism or hydrothermal alteration (e.g., ref. 61). As such, it is critical to evaluate the possible effects of alteration and metamorphism on the compiled whole-rock Fe3+/ΣFe ratios, which we do below. We first present the compiled Fe3+/ΣFe ratios to see if there are any differences in Fe3+/ΣFe ratios of igneous island arc rocks as a function of time. Then we evaluate, through a variety of tests, whether the Fe3+/ΣFe ratios are controlled by primary or secondary processes (e.g., oxidation by atmospheric O2 before collection in the field, metamorphic/hydrothermal alteration, degassing, or differentiation). Following this, we use a completely different oxybarometer, V/Sc, to independently test our interpretations.
Fig. 1 shows the average Fe3+/ΣFe ratios for each island arc locality vs. age. An increase in Fe3+/ΣFe ratios around the Precambrian–Cambrian boundary (541 Ma) is visually apparent. Data histograms are given in SI Appendix, Fig. S1. To test whether there is a significant difference in average Fe3+/ΣFe ratios of igneous island arc rocks from the Phanerozoic vs. the Precambrian, we calculated the average Fe3+/ΣFe ratios over a given age range based on the average Fe3+/ΣFe ratio of each locality in that age range. For this test, we chose the following age ranges: the Late Paleozoic to Cenozoic (420 to 0 Ma), the Early Paleozoic (541 to 420 Ma), and the Precambrian (>541 Ma). These age ranges allow us to explicitly test the hypothesis posed above as to whether island arc rocks are more oxidized in the Phanerozoic vs. the Precambrian following the oxygenation of the deep ocean and increased marine sulfate concentrations ∼800 to 400 Ma. The Early Paleozoic is separated as its own group because some recent studies indicate the large-scale oxygenation of the deep ocean occurred during or after this time interval (2, 6, 9, 10). The Neoproterozoic (1,000 to 541 Ma), which has also been proposed to be a time of increasing atmospheric O2 (e.g., refs. 1, 3, 5, 7, and 12), is not treated separately here due to the limited number of arc localities available from this time frame (n = 5). However, treating the Neoproterozoic as a separate group for Fe3+/ΣFe ratios does not change our conclusions.
Average Fe3+/ΣFe ratios of island arcs rocks through time. Each data point is the average value of a specific island arc locality. Data include both intrusive and extrusive igneous rocks. Solid horizontal lines are averages of the locality means over the given time period: 0.265 ± 0.018 (1 SE) for the Precambrian, 0.307 ± 0.028 for the Early Paleozoic, and 0.382 ± 0.021 for Late Paleozoic–Cenozoic. Shading represents 2 SE. Vertical error bars for individual points are 1 SE. Horizontal age uncertainties vary depending on the age constraints available for each island arc (Dataset S1).
Mean Fe3+/ΣFe ratios for the chosen age ranges are 0.265 ± 0.018 (1 SE; 15 localities, 367 analyses) for the Precambrian, 0.307 ± 0.028 (12 localities, 317 analyses) for the Early Paleozoic, and 0.382 ± 0.021 (17 localities, 528 analyses) for Late Paleozoic–Cenozoic (<420 Ma) (Fig. 1). Our average Fe3+/ΣFe ratio for Late Paleozoic–Cenozoic island arc rocks are similar to mean values of similarly aged island arcs based on data from precompiled igneous island arc rocks from GEOROC, 0.376 ± 0.013 (1 SE) (Fig. 1), and a value of 0.38 as given by Lécuyer and Ricard (23). The mean value for the Late Paleozoic–Cenozoic is significantly different (P < 0.05) from the Precambrian mean value (P = 0.0005; Tukey’s honestly significant difference test), while the Early Paleozoic mean is neither statistically distinct (P > 0.05) from the Late Paleozoic–Cenozoic mean (P = 0.07) nor the Precambrian mean (P = 0.52). This indicates that mean Fe3+/ΣFe ratios of igneous island arc rocks from Phanerozoic are indeed statistically more elevated than those from the Precambrian. Finally, we examined the fidelity of our choice of age bins by fitting a moving average to the data (SI Appendix, Fig. S2). This approach also shows an increase in Fe3+/ΣFe after 541 Ma and is thus consistent with our analysis above.
The average Late Paleozoic–Cenozoic island arc rock Fe3+/ΣFe ratio of ∼0.38 is elevated compared with typical measurements of Fe3+/ΣFe based on Fe K-edge X-ray absorption near-edge structure (XANES) measurements of modern island arc glasses from, for example, the Marianas [0.239 ± 0.003 (1 SE) (58, 59, 62)] and Vanuatu [0.283 ± 0.019 (63)]. To better understand this difference, we compared measurements of Fe3+/ΣFe ratios based on wet chemical extractions from GEOROC vs. XANES measurements on igneous rocks from these two island arc localities as well as Hawaii (SI Appendix, Fig. S3 and ref. 64). Whole-rock Fe3+/ΣFe ratios from GEOROC are offset to higher values compared with XANES-based estimates of glasses by, on average, 0.074 ± 0.024 (1σ). The distributions also differ in that the whole-rock samples show a tail to elevated Fe3+/ΣFe that is not present in the XANES data (SI Appendix, Figs. S1 and S3). When distributions are unimodal (e.g., for Marianas and Hawaii) the modes of the smoothed distributions of the whole-rock and XANES data are similar (within 0.015; SI Appendix, Fig. S3). A similar difference is also observed when comparing Fe3+/ΣFe values of ancient submarine basalts preserved on the continents (i.e., ophiolites) vs. XANES on MORB. That is, average Fe3+/ΣFe ratios for Precambrian submarine basalts from ophiolites range from 0.20 to 0.26 (10), 0.06 to 0.12 higher than that determined by XANES for modern MORB (0.14; ref. 65). Precambrian whole-rock samples were used for this comparison to avoid issues associated with oxidative alteration of MORB during hydrothermal circulation that is hypothesized to have occurred in Phanerozoic samples (10).
We interpret the elevated Fe3+/ΣFe ratios of the whole-rock vs. XANES data (SI Appendix, Figs. S1 and S3) as being due to a small amount of oxidation of samples on Earth’s surface. This oxidation appears to bias average whole-rock Fe3+/ΣFe ratios to slightly above (∼0.074) eruptive values but is not an issue for XANES measurements as these are made on pristine, unaltered glasses. As our compilation is based on submarine rocks, oxidation by atmospheric O2 during eruption is not the cause of elevated whole-rock vs. XANES Fe3+/ΣFe ratios. The GEOROC data include samples from both subaerial and submarine settings. It has a mean Fe3+/ΣFe ratio similar (within 0.01) to our compilation of submarine rocks. This likely indicates that oxidative alteration by atmospheric O2 during eruption does not strongly bias the GEOROC data for subaerially erupted samples. Rather, we propose that the elevated whole-rock vs. XANES Fe3+/ΣFe ratios are most plausibly due to recent subaerial exposure on the Earth’s surface following crystallization and exposure as otherwise the Precambrian submarine basalts would not also show this elevation. This is because these older samples formed in environments (the deep ocean) that are thought to have been essentially devoid of O2 (e.g., refs. 5 and 11). We expect that any oxidation that has occurred recently during exposure would affect the mean Fe3+/ΣFe ratios of all of the localities approximately equally regardless of age as we are unaware of any reason that ancient samples exposed at Earth’s surface would oxidize significantly more slowly in the presence of modern atmospheric O2 than younger samples. We discuss the effects of oxidation by dissolved O2 due to temporal changes in marine dissolved O2 concentrations in the next section.
Finally, we tested whether the presence of a small number of highly oxidized samples (presumably via oxidation by modern atmospheric O2) biases our record such that more highly oxidized and altered samples are Phanerozoic in age. To conduct this test, we filtered our dataset such that samples with Fe3+/ΣFe ratios >0.8, 0.7, 0.6, or 0.5 were excluded (SI Appendix, Table S1). Differences in mean Fe3+/ΣFe ratios between the Precambrian vs. Late Paleozoic samples change from a maximum difference of 0.116 (±0.028, 1 SE) without any filter to a minimum difference of 0.089 (±0.022, 1 SE) with removal of all samples with Fe3+/ΣFe >0.5. This shows that the differences we observe between Precambrian and Late Paleozoic–Cenozoic samples are not driven by a small number of samples with exceptionally high Fe3+/ΣFe ratios.
Based on these tests, we proceed based on the working hypothesis that the averages of the compiled Fe3+/ΣFe ratios in each age range can be used to understand relative changes in Fe3+/ΣFe values with time. However, the observed systematic (though small) elevation in whole-rock Fe3+/ΣFe ratios relative to XANES measurements, (which we interpreted above to indicate that whole-rock samples are commonly slightly oxidized postcrystallization) makes it challenging to use the Fe3+/ΣFe values to calculate the fO2 of the magmas accurately; therefore, we do not attempt this here.
Causes of the Elevated Phanerozoic vs. Precambrian Fe3+/ΣFe Ratios of Island Arc Rocks.
We now explore several possible explanations for the origin of the increase in the average island arc rock Fe3+/ΣFe values from ∼0.27 to 0.38 from the Precambrian to the late Paleozoic. This increase with time could result from a variety of secular changes, including (i) temporal changes in the nature or degree of alteration or metamorphism, (ii) changes in the degree of differentiation or degassing of island arc magmas with time, or (iii) changes in the fO2 of island arc mantle sources with time.
Temporal changes in the nature or degree of alteration or metamorphism.
One possibility is that the observed changes in Fe3+/ΣFe ratios with time could result from postdepositional alteration of island arc igneous rocks caused either by a general reductive metamorphism of ancient samples (with Precambrian samples originally being more oxidized) relative to younger samples or to increased seafloor oxidation of younger arc samples (with Phanerozoic samples being originally more reduced than their current levels). As iron is the main redox active element in igneous rocks, closed-system metamorphism will only redistribute Fe3+ and Fe2+ among different mineral assemblages but will not change whole-rock Fe3+/ΣFe ratios (66). In contrast, open-system metamorphism can introduce reductants or oxidants to the system and thus can alter whole-rock Fe3+/ΣFe ratios (66).
The lower Fe3+/ΣFe ratios for Precambrian vs. Late Paleozoic and younger rocks could be unrelated to original rock redox states if metamorphism is typically reducing given that older rocks have a greater likelihood to be metamorphosed than younger rocks. Stolper and Keller (10) presented arguments for why metamorphism does not appear to universally lower Fe3+/ΣFe ratios of ancient igneous rocks. In brief, they showed that in contrast to Archean samples, numerous Meso-Paleoproterozoic (1 to 2.5 billion years ago) subaerially erupted igneous rocks exist that are highly oxidized (Fe3+/ΣFe > 0.4) and attributed this oxidation to the presence of O2 in the Proterozoic atmosphere, which began to accumulate starting ∼2.5 to 2.3 billion years ago (e.g., refs. 5 and 11). As a specific example, greenschist metamorphic facies subaerial volcanic rocks from the 2,058-Ma Kuetsjärvi volcanic formation reach Fe3+/ΣFe > 0.75 (67). This shows that metamorphism does not universally lower Fe3+/ΣFe ratios of ancient igneous rocks.
As our samples are submarine in origin, oxidation by atmospheric O2 during eruption for extrusive samples is not an issue for our dataset, unlike for subaerial volcanic rocks. Stolper and Keller (10) also showed that the Fe3+/ΣFe ratios of submarine basalts increase beginning between 541 to 420 Ma. They interpreted this as occurring due to the Phanerozoic oxygenation of the deep ocean carrying dissolved O2 into hydrothermal systems circulating through oceanic crust—that is, due to oxidative alteration on the seafloor. Although large-scale hydrothermal circulation systems are best known from midocean-ridge and back-arc spreading centers, island arcs also create local hydrothermal cells (e.g., ref. 68). Consequently, the observed Phanerozoic increase in average submarine island arc rock Fe3+/ΣFe ratios could be due to oxidation of submarine island arc igneous rocks by O2 during hydrothermal circulation following the oxygenation of the deep ocean sometime between 800 and 400 Ma (e.g., refs. 1, 2, 4⇓–6, and 10). We note, though, that the increase in Late Paleozoic–Cenozoic island arc rock Fe3+/ΣFe ratios is approximately one-third of that seen in submarine basalts from ophiolites (10), indicating that hydrothermal circulation is unlikely to be the primary control on Fe3+/ΣFe ratios of island arc igneous rocks.
Regardless, to test for time-variant subsolidus oxidation of the compiled samples during hydrothermal circulation, we compared Fe3+/ΣFe ratios of intrusive vs. extrusive island arc igneous rocks from our compilation as a function of age (SI Appendix, Table S2). Modern extrusive rocks generally experience more oxidative alteration during hydrothermal circulation on the seafloor than intrusive rocks due to both their emplacement at Earth’s surface and higher permeabilities; for example, hydrothermally altered basalts from drilled oceanic crust and exposed ophiolite sequences and are typically elevated in Fe3+/ΣFe by >0.2 compared with sheeted dykes and gabbros (e.g., refs. 69 and 70). Thus, extrusive rocks are expected to be more prone to hydrothermal oxidation on the seafloor than are intrusive rocks. Average Fe3+/ΣFe for the Precambrian, Early Paleozoic, and Late Paleozoic–Cenozoic extrusive island arc rocks are 0.272 ± 0.022 (1 SE), 0.314 ± 0.030, and 0.366 ± 0.023, respectively. Average Fe3+/ΣFe for intrusive island arc rocks for the same age ranges are 0.266 ± 0.028, 0.245 ± 0.017, and 0.425 ± 0.039 (SI Appendix, Fig. S4). Thus, for both the intrusive and extrusive groups, the Late Paleozoic–Cenozoic mean Fe3+/ΣFe ratios for island arc rocks are greater (beyond 2 SE) than for the Precambrian samples. Additionally, the mean values for extrusive and intrusive rocks overlap at the 2 SE level for the specific age bins. Likewise, for the GEOROC island arc database (all Late Paleozoic or younger) the average Fe3+/ΣFe ratios for extrusive and intrusive island arc rocks are 0.389 ± 0.015 (1 SE) and 0.361 ± 0.021 and thus within 2 SE of each other. Consequently, there appears to be no significant difference in average Fe3+/ΣFe ratios between intrusive and extrusive island arc rocks through time for the given age ranges, and the observed increase in Fe3+/ΣFe ratios with time is seen in both rock types. Finally, there is no apparent correlation between Fe3+/ΣFe and the water content of the rocks (as given by the loss on ignition weight percent), indicating that hydrothermal alteration does not inherently promote oxidation or reduction of iron (SI Appendix, Fig. S5). Based on this, we infer that subsolidus oxidation during hydrothermal circulation of island arc rocks is not the cause of the observed difference in average Fe3+/ΣFe of late Paleozoic–Cenozoic vs. Precambrian island arc rocks.
Changes in the degree of differentiation or degassing of island arc magmas with time.
Differentiation of a magma through fractional crystallization of minerals that preferentially incorporate Fe2+ over Fe3+ (e.g., olivine) will increase the Fe3+/ΣFe value of a melt if a magma behaves as a closed system in regard to fO2 (71⇓⇓–74). In contrast, crystallization and removal of Fe3+-rich minerals (e.g., magnetite) from a melt will decrease the melt Fe3+/ΣFe ratio (62, 75⇓–77). We tested whether differences in the degree of fractionation has any control on the difference in mean Fe3+/ΣFe ratios between the various age ranges by filtering samples based on their SiO2 or MgO contents, where higher SiO2 and lower MgO contents are indicative of greater differentiation. Even after removal of samples with >70, 65, and 60 wt % SiO2 or MgO <1, 2, 3, 4, and 5 wt %, mean Fe3+/ΣFe values for the Precambrian vs. Late Paleozoic–Cenozoic do not overlap at the 2 SE level and change in value by at most 0.036 after filtering (SI Appendix, Table S1). This indicates that differences in degree of fractionation with time do not significantly influence the observed changes in Fe3+/ΣFe ratios with time.
Degassing of a mixed H–C–O–S fluid can cause Fe3+/ΣFe ratios of a melt to increase, decrease, or potentially not vary at all depending on which volatile species are lost from the melt (e.g., refs. 57, 58, 62, and 78⇓⇓⇓–82). For example, H2 loss due to H2O dissociation can increase melt fO2 and thus Fe3+/ΣFe ratios (82) and has been called on as a potential explanation for the difference in Fe3+/ΣFe ratios of island arc basalts vs. MORB (33, 36). If a change in the amount of degassing of H2O (or another volatile that causes the oxidation of residual melt) is the cause of the increase of Fe3+/ΣFe ratios of island arcs rocks in our compilation, then we would expect that this difference would be more pronounced in extrusive rocks, which experience more degassing than intrusive rocks. However, as already discussed, Fe3+/ΣFe ratios of intrusive vs. extrusive rocks for a given age range are similar and show a similar difference in Fe3+/ΣFe for Phanerozoic vs. Precambrian rocks. This indicates that the more oxidized nature of Phanerozoic vs. Precambrian island arc rocks is not due to a global change in the style of degassing in Phanerozoic vs. Precambrian island arc rocks.
Changes in the redox state of island arc mantle sources with time.
Given the above considerations, we propose that the increase in Fe3+/ΣFe ratios of island arc igneous rocks from the Precambrian to the Late Paleozoic (541 to 420 Ma) is primary. We further hypothesize that this increase results from changes in the fO2 of the mantle sources of island arc magmas and consequently of primitive island arc rocks. This necessarily requires that the Fe3+/ΣFe ratios of island arc rocks primarily reflect the fO2 of their mantle sources (e.g., refs. 22 and 24) as opposed to being controlled by oxidation/reduction during crystallization, ascent, and/or eruption (e.g., refs. 33 and 36). Additionally, this indicates that oxidizing fluids derived from subducted serpentinized ultramafic rocks do not fully control the redox budget of the subarc mantle as otherwise Fe3+/ΣFe ratios of island arc rocks of all ages should be elevated. This presumes there were not large-scale changes around the Precambrian–Cambrian boundary in rates of serpentinization or how serpentinite-derived oxidants are transferred to the subarc mantle—we are unaware of any evidence for such a change.
Below we propose that the Phanerozoic increase in Fe3+/ΣFe ratios of island arc igneous rocks is due to increased dissolved O2 and sulfate concentrations in the deep ocean 800 to 400 Ma. The higher O2 and sulfate concentrations increased the redox state of sediments and hydrothermally altered oceanic crust that is eventually subducted to the subarc mantle, oxidizing it and, as a result, increasing the fO2 of the parental melts of island arc rocks (10, 18, 19). We additionally comment below on the observed scatter in average locality Fe3+/ΣFe ratios observed in light of this proposal. Before discussing this further, we first explore an obvious prediction of our proposal: If the increase in Fe3+/ΣFe does indeed represent an increase in the fO2 of the subarc mantle, then other redox-sensitive parameters that can track the fO2 of island arc igneous rocks should also increase over this same time period. In the next section we evaluate whether this is the case for the widely used V/Sc oxybarometer before exploring the geological implications of the compiled records further below.
V/Sc Ratios of Island Arc Rocks Through Time.
An alternative parameter to Fe3+/ΣFe ratios commonly used to track the fO2 of a mantle source and its partial melts is the whole-rock V/Sc ratio (32, 33). Specifically, as the fO2 of a peridotite source increases, the V/Sc ratio increases in a partial melt due to the increased bulk incompatibility of V as its oxidation state increases (34, 83). In contrast to Fe3+/ΣFe ratios, V/Sc ratios of igneous rocks have been argued to be largely impervious to low-temperature alteration, metamorphism, and the effects of differentiation and degassing for high-MgO rocks (32, 33). Although V/Sc ratios of basalts are also sensitive to degree of partial melting and the initial V/Sc ratio of their mantle source (32, 33, 84), if these parameters are roughly constant through time for island arc magmatism, then V/Sc ratios of island arc rocks are expected to be sensitive to variations in the fO2 of a peridotitic source. Thus, to complement and further test the significance of the Fe3+/ΣFe-ratio record, we also compiled V/Sc ratios of island arc igneous rocks from the past 3.1 billion years.
In Fig. 2, we provide the mean V/Sc value of each island arc locality as a function of age as well as the average value for of all localities from a given age range. Data histograms are given in SI Appendix, Fig. S6. Here we have combined our compilation of Late Paleozoic–Cenozoic island arc samples with those of the GEROC database (which cover the same time frame) to increase data coverage from the Late Paleozoic to Cenozoic. We also provide average V/Sc ratios of modern MORB based on mean values at a given plate boundary, 6.83 ± 0.10 (1 SE; data from ref. 55), which is in agreement with that given by Li and Lee (32) of 6.74. Additionally, we present measurements of V/Sc ratios of Precambrian basalts previously used to reconstruct the fO2 of the upper mantle (Fig. 2 and refs. 32 and 35). These yield a mean V/Sc ratio of 6.38 ± 0.24 (1 SE), within 2 SE of modern MORB. This lack of a discernible change in V/Sc in nonsubduction related basalts is consistent with previous studies that found no clear change in the fO2 of the magma-producing upper mantle potentially from the Archean or at least Proterozoic (2.5 billion years ago) to today (32, 35, 85, 86). Note that we have not included V/Sc ratios of eclogites (35) or komatiites (85) used to reconstruct past mantle fO2 values given their different formational environments, degrees of melting, and metamorphic histories compared with typical basalts. Such rocks potentially indicate lower (∼1 log unit) mantle fO2 values compared with modern before 2.5 billion years ago (35), but similar mantle fO2 values from the Proterozoic on (2.5 to 0 billion years ago; refs. 35 and 86).
Average V/Sc ratios of island arc rocks through time. Each data point is the average value of a specific island arc locality. Data include both intrusive and extrusive igneous rocks. Solid horizontal lines are averages of the formational means over the given time period: 6.49 ± 0.19 (1 SE) for the Precambrian, 6.74 ± 0.35 for the Early Paleozoic, and 7.57 ± 0.17 for the Late Paleozoic–Cenozoic. Shading represents 2 SE. Vertical error bars for individual points are 1 SE. Horizontal age uncertainties vary depending on the age constraints available for each island arc (Dataset S1). The gray horizontal bar is the V/Sc range of the averages of modern MORB for specific plate boundaries (from ref. 55) and ancient basalts used previously to reconstruct past upper mantle fO2 (data from refs. 32 and 35): average of 6.69 and range of 5.59 to 7.59.
In contrast, mean V/Sc ratios of the island arc localities using the same age ranges as for the Fe3+/ΣFe ratios are 6.49 ± 0.19 (1 SE; 27 locations, 307 analyses) for the Precambrian, 6.74 ± 0.35 (9 locations, 129 analyses) for the Early Paleozoic, and 7.57 ± 0.17 (11 locations compiled by us and 24 from GEOROC, 98 analyses compiled by us and 2,855 from GEOROC) for the Late Paleozoic–Cenozoic. The average V/Sc ratio for the Late Paleozoic–Cenozoic differs significantly (P < 0.05) from that of the Precambrian (P = 0.0008; Tukey’s honestly significant difference test), while the Early Paleozoic mean is neither statistically distinct (P > 0.05) from the Late Paleozoic–Cenozoic mean (P = 0.10) nor the Precambrian mean (P = 0.95). Use of a 6.5 wt % MgO filter does not control the result: Late Paleozoic–Cenozoic island arc rock V/Sc ratios are elevated relative to Precambrian island arc rocks by 1.08 using an MgO filter of 6.5 wt % and by 1.20 using a filter of 8.0 wt %. Finally, our average V/Sc ratio for the Late Paleozoic–Cenozoic time period (7.57 ± 0.17) is within 0.12 of the value for similarly aged island arc volcanic rocks independently compiled in Turner and Langmuir (87) of 7.69 ± 0.28 (1 SE) when an MgO filter of 6.5 to 15 wt % is applied to that data. Again, we examined the fidelity of choice of our age bins by fitting a moving average data to the data (SI Appendix, Fig. S7). As was the case for the Fe3+/ΣFe ratio, this approach shows an increase in V/Sc ratios after 541 Ma, and is thus consistent with our analysis above.
An increase of ∼1 in the V/Sc ratio corresponds to an increase in ∼0.5 log units of fO2 at 10 to 20% partial melting of a spinel lherzolite with V and Sc contents the same as that for the source of MORB (32, 33). Some Cenozoic localities, however, record average V/Sc values of 9 to 10, which could correspond to an increase in fO2 of a log unit or more at similar melting conditions. Thus, like the Fe3+/ΣFe ratios, the V/Sc ratios are consistent with a statistically significant increase in the redox state of island arc mantle sources and primary melts from the Neoproterozoic to Phanerozoic. Importantly, this increase does not appear to occur in nonisland arc basalts and is instead specific to island arc rocks (Fig. 2). Additionally, we note that our compilation and that of ref. 87 yield average V/Sc ratios of modern island arc rocks that are elevated vs. MORB by ∼0.8 (and differ beyond 2 SE)—this is in contrast to previous work that argued for no significant difference between the two (33, 34). If the elevation in V/Sc observed here between modern island arcs and MORB is interpreted following typical frameworks (32⇓–34), then it would indicate the source of island arc rocks is more oxidized (higher fO2) relative to the source of MORB. If correct, then the V/Sc and Fe3+/ΣFe ratios of island arc rocks and MORB are actually in qualitative agreement on this issue (22, 24, 58, 59).
In SI Appendix, Fig. S8 we directly compare V/Sc vs. Fe3+/ΣFe ratios of time period means and locality averages. Linear regressions of both yield a positive correlation, as would be expected. This direct comparison of V/Sc and Fe3+/ΣFe ratios is complicated by two issues. First, only one-third of the non-GEOROC localities have both Fe3+/ΣFe and V/Sc ratios available, and most (92%) of the V/Sc and Fe3+/ΣFe measurements in our compilation are not on the same samples. Thus, a clear area of future work is a direct comparison of Fe3+/ΣFe vs. V/Sc measured on the same samples for Precambrian vs. Phanerozoic aged rocks.
Second, although both V/Sc and Fe3+/ΣFe ratios in partial melts of the mantle are sensitive to source fO2, the exact dependence of V/Sc on fO2 is a function of the degree of melting (33, 84), modal mineralogy of the source (33), and the initial V/Sc ratio of the source (84). Lower initial V/Sc ratios of the subarc mantle due to prior melt extraction (e.g., from back-arc spreading centers) has recently been proposed to strongly influence the relationship between V/Sc ratios of arc basalts and the fO2 of the source region (84). We explored this issue quantitatively by modeling V and Sc partitioning during partial melting of the mantle (SI Appendix). We examined how V/Sc vs. Fe3+/ΣFe ratios vary for a primitive basalt given plausible ranges in source fO2, source V/Sc, modal mineralogy, and melt fraction (SI Appendix, Fig. S9). In SI Appendix, Fig. S10 we overlay the potential modeled ranges of V/Sc vs. Fe3+/ΣFe on our compiled data corrected for an increase in 0.074 in Fe3+/ΣFe ratios due to later oxidation (as based on the XANES comparison above). All data fall within a modeled fO2 range of −1 to 4 log units below or above QFM, a typical range for MORB and island arc sources (21, 22, 29, 58, 59). This shows that our measured V/Sc and Fe3+/ΣFe are reasonable and in agreement. Our model reinforces the conclusion of ref. 84 that the precise relationship between V/Sc and melt fO2 is model-dependent.
We expect that the generation of island arcs has been similar through time (similar degrees of melt and similar amount of prior melt depletion), and as such that V/Sc ratios can be used to track relative changes in the fO2 of island arc sources with time. Based on this, we propose the observed increase in V/Sc ratios to be most simply interpreted to reflect an increase in the fO2 of the subarc mantle sometime between the end of the Neoproterozoic through to the early Paleozoic. Most importantly, this is in agreement with and provides strong independent evidence for our interpretation of the Fe3+/ΣFe data as reflecting an increase in the redox state of the subarc mantle between 800 and 400 Ma.
Finally, we also examined V/Yb ratios for our compilation as this ratio has independently been proposed to be a oxybarometer (37). We observe a hyperbolic relationship between V/Yb vs. Yb (SI Appendix, Fig. S11) such that at lower (<2 ppm) Yb concentrations V/Yb ratios are clearly elevated relative to samples with higher (>2 ppm) Yb concentrations. We expect that this hyperbolic relationship derives, in part, from imprecisions and inaccuracies in Yb measurements such that at the already low Yb concentrations [most (62%) samples have <2 ppm Yb] small (<0.5 ppm) measurement errors generate the hyperbolic relationship. Given this complexity (which is not present in the V/Sc dataset; SI Appendix, Fig. S11) we do not apply the V/Yb oxybarometer to our dataset.
Porphyry Copper Deposits.
It has been proposed that oxidized arc melts (relative to MORB) are a necessary condition for the generation of porphyry copper deposits commonly associated with island and especially continental arcs (25, 88), although this is debated (e.g., ref. 41). Porphyry copper deposits in the Precambrian are rare and a rise in their occurrence begins at the start of the Phanerozoic (89) coinciding with the increase in Fe3+/ΣFe and V/Sc ratios of the island arc compilation presented here (Fig. 3 A–C). Some have ascribed this increase to a bias in the rock record in which preferential erosional loss of Precambrian deposits has occurred (90). Alternatively both Evans and Tomkins (18) and Richards and Mumin (19) proposed that the Phanerozoic increase is a primary feature of the rock record and that the increase is due to the oxygenation of the deep ocean, which increased marine sulfate concentrations, which in turn led to an increase in both the sulfur content and redox state of arc parental melts via subduction and incorporation of this sulfate into the subarc mantle. The key point is that although the formation (91) and preservation of porphyry copper deposits is potentially complex their temporal record is consistent with and has previously and independently been linked to the idea of a Phanerozoic increase in the redox state of the mantle source of island arc magmas as we have proposed here.
Parameters associated with changes in the redox state of island arc rocks with time. (A) Fe3+/ΣFe ratios of island arc rocks through time as in Fig. 1. (B) V/Sc ratios of island arc rocks through time as in Fig. 2. (C) Abundance of porphyry copper deposits with time (data from ref. 89). (D) Fe3+/ΣFe ratios of submarine basalts from ophiolites through time indicative of the redox state of altered oceanic crust (data from ref. 10).
Effects of Deep-Ocean O2 and Marine Sulfate Concentrations on the fO2 of Island Arc Mantle Sources.
Based on the island arc Fe3+/ΣFe and V/Sc records and the temporal distribution of porphyry copper deposits, we propose an increase in the redox state of the mantle source of island arcs occurred sometime between the end of the Neoproterozoic and the Early Paleozoic. The question is, what geological process (or processes) drove this change? As discussed, a leading hypothesis for why modern igneous island arc rocks are more oxidized vs. MORB is that oxidized material from the subducted slab is transferred to and thus oxidizes the subarc mantle, which melts to form island arc rocks. Thus, an increase in the redox state of subducted oceanic crust from the Neoproterozoic to Early Phanerozoic provides a plausible explanation for the observed increase in the redox state of island arcs over this time frame.
Over the same time interval that our compilation indicates that island arc magmas became more oxidized (800 to 400 Ma), the atmosphere is thought to have accumulated sufficient O2 (10 to 50% of modern) to oxygenate the majority of the deep oceans (e.g., refs. 1⇓⇓⇓⇓⇓⇓⇓⇓–10). Additionally, marine sulfate concentrations are thought to have also increased to near modern levels over this same time frame (e.g., refs. 13⇓⇓⇓–17). Based on this, we propose that this oxygenation and elevation in marine sulfate concentrations increased the redox state of altered oceanic crust that, when eventually subducted to the mantle, oxidized the mantle source of island arc rocks.
There are at least three ways the oxygenation of the deep ocean and an increase in marine sulfate concentrations could have increased the redox state of subducted oceanic crust. (i) Deep ocean oxygenation would have increased the redox state of deep-sea sediments by increasing sedimentary Fe3+/ΣFe ratios and potentially lowering organic carbon contents (burial of organic carbon is enhanced in anoxic vs. oxic settings; e.g., ref. 92). (ii) Circulation of O2 derived from deep ocean waters would have oxidized iron and sulfur in oceanic crust—indeed, Fe3+/ΣFe ratios of submarine basalts from ophiolites show a threefold increase from the Neoproterozoic to Late Paleozoic–Cenozoic (Fig. 3D and ref. 10). This independently demonstrates that there was an increase in the redox state of altered oceanic basalts over the same time frame the redox state of island arc rocks increases. (iii) Increased dissolved sulfate levels in the ocean could have increased the redox state of subducted oceanic crust at this time (18, 19) either during high-temperature (>150 °C) hydrothermal circulation with anhydrite precipitation (if anhydrite does not redissolve before subduction; refs. 93 and 94) and/or oxidation of Fe2+-bearing minerals by sulfate, which occurs in experiments at elevated (>200 °C) temperatures (95). Which of these processes is most important in causing the observed increase in island arc redox state is not clear, and it is possible, and perhaps likely, that all played a role—in other words, our record does not identify the oxidant (Fe3+ vs. sulfate) or source of the oxidant (sediments vs. altered igneous rocks) of the subarc mantle. A schematic cartoon of how changes in deep-ocean O2 and sulfate concentrations with time could have influenced the redox state of island arc rocks is given in Fig. 4. We note that we have not considered changes in carbonate subduction fluxes with time as it has been argued that subducted carbonate is not a source of oxidants for island-arc melts (18) because at the depths that these melts are derived (⪅120 km), carbon exists as C4+ in carbonate or CO2 (96). However, if these melts formed at greater depths where graphite is stable, subducted carbonate could oxidize the subarc mantle (40).
(A and B) Schematic cartoon of how proposed changes in the redox state of the deep ocean alters the redox state of igneous island arc rocks.
For our given age bins, there is significant scatter in Fe3+/ΣFe and V/Sc values (Figs. 1 and 2). For example, although most Late Paleozoic–Cenozoic samples are more oxidized than Precambrian equivalents, some Late Paleozoic–Cenozoic localities are lower in Fe3+/ΣFe and V/Sc than the average value for the Precambrian. Some of this scatter likely results from alteration of some samples or methodological differences for the measurement of FeO and Fe2O3 wt % or V and Sc abundances between studies. However, it not expected that there should by a single value for Fe3+/ΣFe or V/Sc ratios of island arc rocks for a given age range once subduction of oxidized material begins. This is because the extent of subarc oxidation during subduction has been proposed to be controlled by the age of the arc, the total amount of oceanic crust subducted, and convergence rate, all of which influence the integrated flux of oxidants to a subduction zone and the ability of those oxidants to be released at depth and infiltrate the subarc mantle vs. escape up the slab–mantle interface (29). Thus, based on previous work (29), records of Fe3+/ΣFe and V/Sc ratios from similarly aged island arc localities are expected to show a range of values.
The Fe3+/ΣFe and V/Sc records are unlikely to be able to provide a definitive estimate of when precisely the deep ocean became oxygenated or experienced increased sulfate concentrations given both the scatter in the records and the inherent lag of potentially up to ∼108 years between when oceanic crust becomes oxidized on the seafloor, when it is subducted and releases these oxidants to the subarc mantle, when the subarc melts, and when those melts erupt and crystalize to form island arc rocks. For example, Evans and Tomkins (18) estimated that for an oxygenated deep ocean 2.4 × 1012 mol/y of oxidants are subducted. This combined with assumptions about how material is transferred from the slab to the subarc mantle and mantle mixing rates, they (18) modeled that it takes between 104 and 107 years to noticeably (>0.3 log unit change in fO2) oxidize the subarc mantle and that it takes 106 to 108 years for the subarc mantle to reach a steady state fO2 of 1 to 3 log units above the MORB source. This does not account for the ∼107 years needed to oxidize the oceanic crust on the seafloor.
Additionally, although a rise in both marine seawater sulfate (e.g., refs. 13⇓⇓⇓–17) and deep-ocean O2 concentrations (e.g., refs. 1, 2, 4⇓–6, and 10) are thought to have occurred toward the end of the Neoproterozoic through the early Phanerozoic, whether these occurred simultaneously is not known. Thus, an increase in marine sulfate levels could predate and cause the oxidation of the subarc mantle by sulfate (in part or fully) before the oxygenation of the deep ocean. Regardless, the presented island arc record is broadly consistent with independent work that indicates the deep oceans were anoxic and relatively sulfate-poor (compared with modern) before the Neoproterozoic (>1,000 Ma), and that permanent deep-ocean oxygenation and sulfate accumulation occurred sometime during or after the Neoproterozoic (1, 2, 4⇓–6, 8, 10, 12⇓⇓⇓⇓–17, 97). We note that atmospheric O2 and marine sulfate levels are proposed to have temporarily reached near modern levels ∼2 billion years ago (in the aftermath of the initial oxygenation of the atmosphere) before subsiding back to lower values within ∼108 years (98, 99). Our data compilation lacks the temporal resolution needed to test whether any oxidation of the subarc mantle occurred due to subduction of oxidized material generated during this brief excursion in O2 and sulfate levels.
So far, we have focused on how the oxygenation of the deep ocean and the associated increase in marine sulfate concentrations could have led to an increase in the fO2 of mantle sources of island arc magmas. However, it is worth speculating on how the redox state of subducted oceanic crust could have differed in the Precambrian when the ocean is thought to have been anoxic with either Fe2+ or sulfide as the main redox-active dissolved species (5, 11, 97, 100). Anoxic deep oceans could have led to an increase in the relative amount of organic carbon reaching sediments vs. that sinking out of the mixed layer of the ocean (e.g., ref. 92) as well as a decrease in the Fe3+/ΣFe ratio of Fe in deep-sea sediments. Additionally, dissolved Fe2+ or sulfide in the deep ocean could have been fixed in altered oceanic basalts during hydrothermal alteration as iron-rich clays or sulfide minerals (e.g., pyrite). If sufficient organic carbon and/or Fe2+ and reduced sulfur-bearing minerals were sequestered in sediments or altered oceanic crust, subduction of oceanic crust could have actually reduced the subarc mantle in the Precambrian before the oxygenation of the deep ocean.
Although speculative, this idea is consistent with the data presented here. The V/Sc ratios of some Precambrian (and particularly Mesoproterozoic and older) island arc rocks are lower than those for modern MORB (Fig. 2), potentially indicative of a more reduced source for island arc parental melts. We note that there is another potential explanation for the lower V/Sc ratios of some Mesoproterozoic and older island arc rocks vs. modern MORB. Specifically, a depleted mantle source for island arc parental melts has been proposed to lower the initial V/Sc of the mantle source for modern island arcs (84). Without the compensatory effect of oxidized material in the subducted slabs to increase the subarc mantle fO2 (and thus raise the V/Sc ratio of the island arc parental melts), this prior melt depletion could have acted to lower the V/Sc ratios of Mesoproterozoic and older island arc rocks. Whether the subduction of reduced altered oceanic crust vs. source depletion controls the V/Sc ratios of Mesoproterozoic and older island arc rocks is unclear. Importantly, both are consistent with an increase in the flux of oxidized material subducted to the subarc mantle causing an increase in V/Sc ratios of island arc rocks as we have proposed here.
Implications for the Formation of Calc-Alkaline Rocks.
As parental melts of MORB crystallize, total iron contents increase forming the “tholeiitic” or “iron-enrichment” trend (e.g., refs. 73, 101, and 102), whereas in some island, and especially continental, arc magmas, total iron contents remain constant or decrease as melts differentiate, following the “iron-depletion” or “calc-alkaline” trend (e.g., refs. 102 and 103). One explanation for these different trends is that the elevated fO2 of arc magmas results in spinels rich in the Fe3+-bearing magnetite component to saturate early during crystallization, extracting significant iron from the melt and preventing iron enrichment (e.g., refs. 20 and 71). Importantly, Archean and early Proterozoic igneous rocks exist that follow the calc-alkaline trend (e.g., refs. 104 and 105). Additionally, it has recently been argued that the ratio of tholeiitic vs. calc-alkaline rocks preserved in continental crust has not changed over the past ∼3.8 billion years (45). As the Fe3+/ΣFe and V/Sc compilations suggest that the fO2 of Precambrian island arcs and MORB mantle sources were similar (Fig. 2), it must follow, according to the hypothesis developed here, that an elevated fO2 of island arc sources (and thus continental arc sources as well) relative to MORB sources is not an essential ingredient in the generation of calc-alkaline rocks characterized by suppression of iron enrichment during crystallization. Instead, other explanations for the calc-alkaline trend, such as the elevated water contents of the mantle sources of primitive island arc magmas (e.g., refs. 106 and 107), must also be important in creating iron-depletion trends in Precambrian arc rocks.
Conclusions
In previous discussions of the evolution of the redox state of the mantle and ocean/atmosphere, a focus has been on how changes in the fO2 of the mantle could (or could not) have caused the initial rise in atmospheric O2 levels ∼2.4 to 2.3 billion years ago via the control of mantle fO2 on the redox state of volcanic gasses emitted to the atmosphere (e.g., refs. 32, 35, 61, 64, 108, and 109). We interpret the results presented here to indicate that the opposite is the case for the later rise of atmospheric O2 sometime between 800 to 400 Ma to near-modern (10 to 50%) levels and the associated rise in marine sulfate concentrations. Over this time period, rather than changes in mantle fO2 controlling the redox state of the atmosphere and ocean, we propose that changes in the redox state of the atmosphere and deep ocean influenced the fO2 of the subarc mantle (i.e., the source region of island arc melts) through subduction of oceanic crust that was oxidized via interaction with an oxygenated deep ocean with elevated sulfate concentrations. This in turn led to an increase in the redox state of island arc magmas, which are the building blocks of new continental crust. This represents a rare example of changes in the biogeochemistry of Earth’s surface leading to significant changes in the geochemical evolution of the mantle and the chemical composition of an important class of igneous rocks. Finally, by using ancient island arc rocks and taking advantage of the change in the boundary conditions of the redox state of the deep oceans through time, this work presents independent and distinct evidence that the more oxidized nature of modern igneous island arc rocks vs. MORB is derived from the transfer of oxidized material in the subducted slab from sediments or altered mafic crust to the subarc mantle as opposed to serpentinization of subducted slabs or processes associated with transit of magmas through the crust and eruption.
Acknowledgments
We thank E. Cottrell, M. Hirshmann, and E. Stolper for helpful discussions; Z. Brook for help compiling the data; J. Blundy and C. Langmuir for feedback on an earlier draft; and M. Brounce and L. Kump for constructive reviews.
Footnotes
- ↵1To whom correspondence may be addressed. Email: dstolper{at}berkeley.edu or cbucholz{at}caltech.edu.
Author contributions: D.A.S. and C.E.B. designed research, performed research, analyzed data, and wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1821847116/-/DCSupplemental.
Published under the PNAS license.
References
- ↵
- Canfield DE,
- Poulton SW,
- Narbonne GM
- ↵
- Dahl TW, et al.
- ↵
- ↵
- ↵
- ↵
- ↵
- Blamey NJ, et al.
- ↵
- Lenton TM, et al.
- ↵
- Wallace MW,
- Shuster A,
- Greig A,
- Planavsky NJ,
- Reed CP
- ↵
- ↵
- ↵
- Planavsky NJ, et al.
- ↵
- ↵
- Canfield DE,
- Farquhar J
- ↵
- Lyons TW,
- Gill BC
- ↵
- Algeo TJ,
- Luo GM,
- Song HY,
- Lyons TW,
- Canfield DE
- ↵
- Turner EC,
- Bekker A
- ↵
- Evans K-A,
- Tomkins A-G
- ↵
- Richards JP,
- Mumin AH
- ↵
- Gill JB
- ↵
- ↵
- ↵
- Lécuyer C,
- Ricard Y
- ↵
- Kelley KA,
- Cottrell E
- ↵
- Mungall JE
- ↵
- ↵
- Wood BJ,
- Bryndzia LT,
- Johnson KE
- ↵
- ↵
- Evans K,
- Elburg M,
- Kamenetsky V
- ↵
- Bénard A, et al.
- ↵
- Christie DM,
- Carmichael IS,
- Langmuir CH
- ↵
- Li Z-XA,
- Lee C-TA
- ↵
- ↵
- ↵
- Aulbach S,
- Stagno V
- ↵
- ↵
- Laubier M,
- Grove TL,
- Langmuir CH
- ↵
- ↵
- Alt JC, et al.
- ↵
- Evans KA
- ↵
- Lee C-TA, et al.
- ↵
- Dauphas N, et al.
- ↵
- Ulmer P,
- Kaegi R,
- Müntener O
- ↵
- Hawkesworth C,
- Cawood PA,
- Dhuime B
- ↵
- Keller B,
- Schoene B
- ↵
- Pearce JA,
- Cann JR
- ↵
- ↵
- ↵
- ↵
- ↵
- Furnes H,
- Dilek Y,
- De Wit M
- ↵
- Safonova I,
- Kotlyarov A,
- Krivonogov S,
- Xiao W
- ↵
- Woodhead J,
- Eggins S,
- Gamble J
- ↵
- Kushiro I
- ↵
- Gale A,
- Dalton CA,
- Langmuir CH,
- Su Y,
- Schilling J
- ↵
- ↵
- ↵
- ↵
- Brounce M,
- Kelley KA,
- Cottrell E,
- Reagan MK
- ↵
- Grocke SB,
- Cottrell E,
- de Silva S,
- Kelley KA
- ↵
- ↵
- Kelley KA,
- Cottrell E
- ↵
- Gentes Z
- ↵
- Brounce M,
- Stolper E,
- Eiler J
- ↵
- Zhang HL,
- Cottrell E,
- Solheid PA,
- Kelley KA,
- Hirschmann MM
- ↵
- ↵
- Rybacki K,
- Kump L,
- Hanski E,
- Melezhik V
- ↵
- de Ronde CE, et al.
- ↵
- Gillis KM,
- Robinson PT
- ↵
- ↵
- Osborn EF
- ↵
- ↵
- Juster TC,
- Grove TL,
- Perfit MR
- ↵
- Tang M,
- Erdman M,
- Eldridge G,
- Lee C-TA
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- Carmichael IS,
- Ghiorso MS
- ↵
- ↵
- ↵
- ↵
- Prytulak J, et al.
- ↵
- ↵
- Nicklas RW,
- Puchtel IS,
- Ash RD
- ↵
- Turner SJ,
- Langmuir CH
- ↵
- Richards JP
- ↵
- Singer DA,
- Berger VI,
- Moring BC,
- Groat CG
- ↵
- Kesler SE,
- Wilkinson BH
- ↵
- Blundy J,
- Mavrogenes J,
- Tattitch B,
- Sparks S,
- Gilmer A
- ↵
- ↵
- Alt JC
- ↵
- Alt JC,
- Davidson GJ,
- Teagle DA,
- Karson JA
- ↵
- ↵
- Dasgupta R
- ↵
- ↵
- Bachan A,
- Kump LR
- ↵
- Blättler CL, et al.
- ↵
- ↵
- ↵
- ↵
- ↵
- Brown GC
- ↵
- ↵
- ↵
- Grove TL,
- Brown SM
- ↵
- Kump LR,
- Kasting JF,
- Barley ME
- ↵
Citation Manager Formats
Sign up for Article Alerts
Article Classifications
- Physical Sciences
- Earth, Atmospheric, and Planetary Sciences