Abstract
Iron formations (IF) preserve a history of Precambrian oceanic elemental abundance that can be exploited to examine nutrient limitations on early biological productivity. However, in order for IF to be employed as paleomarine proxies, lumped-process distribution coefficients for the element of interest must be experimentally determined or assumed. This necessitates consideration of bulk ocean chemistry and which authigenic ferric iron minerals controlled the sorption reactions. It also requires an assessment of metal mobilization reactions that might have occurred in the water column during particle descent and during post-depositional burial. Here, we summarize recent developments pertaining to the interpretation and fidelity of the IF record in reconstructions of oceanic trace element evolution. Using an updated compilation, we reexamine and validate temporal trends previously reported for the nickel content in IF (see Konhauser et al., 2009). Finally, we reevaluate the consequences of methanogen Ni starvation in the context of evolving views of the Archean ocean-climate system and how the Ni famine may have ultimately facilitated the rise in atmospheric oxygen. Key Words: Methane—Biogeochemistry—Chemical evolution—Iron formation—Archean. Astrobiology 15, 804–815.
1. Introduction
I
Within the past 5 years, that hypothesis has been questioned on two grounds. The first is that the partitioning coefficients used to extrapolate dissolved Ni concentrations in the Precambrian oceans may have been underestimated because they did not consider that the initial ferric iron phase could have been green rust (Zegeye et al., 2012). This impacts extrapolated paleoseawater Ni concentrations, as green rust sorbs Ni more efficiently than the precursor ferric hydroxide phase, such as ferrihydrite (Fe[OH]3), necessitating lower equilibrium dissolved Ni concentrations to explain the IF data. The second is that Ni may have been remobilized from ferrihydrite during the particles' descent through the water column and upon burial (Frierdich et al., 2011). If Ni was lost from the particles, then what is recorded in IF is only a fraction of what might initially have been deposited. These two critiques operate in opposite directions, implicating either lower Ni (in the case of a green rust sorbent) or higher Ni (in the case of diagenetic Ni loss) in Archean seawater than that estimated by Konhauser et al. (2009).
In this work, we had three objectives. First, we provide an updated Ni/Fe database expanded by 676 filtered data points, nearly two-thirds greater than what was available at the time of the original study (1085), to determine whether the initial trends in Ni/Fe ratios through time are conserved in the expanded data set. Second, we discuss how recent experimental studies on IF mineral precipitation and diagenesis bear on the utility of the IF trace metal record as a proxy for evolving Precambrian seawater composition. Third, we incorporate evolving views of the Archean ocean-climate system to paint a more complete picture of oxygen's rise during the GOE.
2. Nickel Concentrations in Precambrian IF and Archean Seawater
After updating the IF compilation of Ni/Fe by nearly 700 points (Fig. 1), the trends shown here demonstrate a remarkable fidelity to those first described by Konhauser et al. (2009). A dramatic decline in Ni between 2.7 and 2.5 Ga remains clear. When combined with the distribution coefficients obtained from Konhauser et al. (2009), this updated record would still indicate a paleomarine Ni concentration of up to ∼400 nM in the Archean. Following the precipitous drop in Ni, concentrations would approach ∼200 nM by 2.5 Ga, with modern levels established by the Neoproterozoic. It is encouraging that the first-order, secular trends in paleomarine Ni concentrations preserved in the IF record are also recorded by other proxy records, such as sedimentary to early diagenetic pyrite (Large et al., 2014).

Compilation of Ni/Fe ratios in IF through time; this reflects an updated version of Fig. 1 in the work of Konhauser et al. (2009) with an additional 676 new analyses. Crosses indicate laser ablation analyses of individual grains, and red circles bulk-digest data. (Color graphics available at
An alternative way to evaluate paleomarine Ni concentrations from IF data without using experimentally derived partition coefficients is by examining the scaling between Ni and Fe in the IF record itself (Fig. 2). In doing so, it becomes apparent that there is a correlation between Fe and Ni deposition; that is, as Fe increases so too does Ni, as expected for partition coefficient-type behavior of Ni sorption to Fe minerals. On such a cross plot, we may also superimpose the scaling expected between a trace metal and Fe based on hypothesized partitioning scenarios (cf. Robbins et al., 2013). The advantage of this approach is that hypothesized partitioning scenarios can be overlain on the IF record data and subsequently evaluated. If the hypothesized partitioning is reflective of the actual deposition of IF minerals, then the proposed lines should constrain the majority of IF data. Indeed, when the KD values predicted by Konhauser et al. (2009) are superimposed on the cross plot (Fig. 2), almost the entire range of Ni/Fe values in the IF record fall within the predicted ranges. This concordance does not attest to the accuracy or reality of experimentally determined KD values, as the slope in solid-phase Ni versus Fe space is calculated as KD*[Ni]sw and thus depends equally on [Ni]sw. In the case of the Konhauser et al. (2009) partitioning scenarios, [Ni]sw values were originally extrapolated from the rock record data using experimentally determined KD values, such that the presented KD*[Ni]sw scenarios (slopes in Ni vs. Fe space) are necessarily tuned to this data set, and apparent correspondence may be considered circular. However, what Fig. 2 effectively demonstrates is (1) that partitioning scenarios can be readily evaluated against rock record data, and should be in order to be considered tenable, and (2) that looking forward, supposedly free parameters (KD, [Ni]sw) could conceivably be tuned against rock record data to explore plausible scenarios in KD-[Ni]sw space, as was done by Robbins et al. (2013) for Zn/Fe in the absence of experimentally determined partition coefficients.

Ni-Fe cross plot showing partitioning of Ni as recorded by IF. Blue lines indicate the anticipated scaling of Ni with Fe based on the partitioning scenarios of Konhauser et al. (2009). Also shown as a red line is the anticipated scaling based on the partitioning scenario determined for green rust in Lake Matano (Zegeye et al., 2012). (Color graphics available at
An important consideration in using the IF record to extrapolate the paleomarine concentrations of trace elements is the influence of organic carbon. The reason being that it is becoming increasingly accepted that microorganisms were in some way involved in the primary oxidation of dissolved Fe(II) to the precursor ferrihydrite particles (Konhauser et al., 2007a; Pecoits et al., 2015). Two possible roles for bacteria are envisioned. The original model of biological IF precipitation assumed the oxidation of hydrothermal Fe(II), via abiotic oxidation by any available O2 generated through oxygenic photosynthesis, that is, cyanobacteria (Cloud, 1972; Klein and Beukes, 1989), and/or via biotic oxidation by chemolithotrophic bacteria, that is, microaerophilic Fe(II)-oxidizers, such as the genera Gallionella (Holm, 1989). Phototrophic Fe(II)-oxidizing bacteria (the so-called photoferrotrophs) were discovered two decades ago (Widdel et al., 1993) and represent an attractive pathway for explaining early IF precipitation in the absence of O2 (Ehrenreich and Widdel, 1994; Heising et al., 1999; Posth et al., 2008, 2013a; Köhler et al., 2010; Pecoits et al., 2015). These bacteria fix CO2 phototrophically using Fe2+ as an electron donor, potentially generating large quantities of Fe(III) in the process. The attractiveness of this concept is that it explains IF deposition in the absence of molecular oxygen using dissolved Fe2+ that would have been abundant in seawater at that time, and they would have been capable of oxidizing enough Fe(II) to explain even the largest IF deposits (Konhauser et al., 2002). In fact, calculations based on experimentally determined Fe(II)-oxidation rates by these organisms under light regimes representative of ocean water at depths of a few hundred meters suggest that, even in the presence of cyanobacteria, anoxygenic phototrophs living beneath a wind-mixed surface layer would have had the first opportunity to utilize Fe2+ (Kappler et al., 2005).
If the precursor ferrihydrite was biologically generated, then it stands to reason that organic carbon may have influenced Ni sorption to the precursor ferrihydrite particles. For instance, it has been shown that biogenic ferrihydrite, consisting of cell-mineral aggregates, differs from abiogenic minerals in terms of chemical composition, particle size, and density (Posth et al., 2010), as well as surface sorptive properties (Moon and Peacock, 2012, 2013; Muehe et al., 2013). To this end, Eickhoff et al. (2014) examined the partitioning of Ni to both abiogenic and biogenic ferrihydrite. In these experiments, the biogenic ferrihydrite consisted of a cell-mineral aggregate precipitated by either a freshwater or marine photoferrotroph. When normalized to specific surface area, biogenic ferrihydrite sorbed less Ni, apparently due to microbially derived organics or whole cells binding to the mineral surface and thereby competing with Ni for sorption sites. Eickhoff et al. (2014) suggested that, given the potential for co-precipitating organic matter to depress partitioning of Ni, (i) previous estimates for Ni concentrations in Archean seawater may be too low, and (ii) the decline of Ni in Precambrian ocean may have occurred closer to the GOE, at 2.45 Ga. However, comparing the data from Eickhoff et al. (2014) to those presented by Konhauser et al. (2009) requires caution. While Konhauser et al. (2009) used a concentration of Fe (0.179 mM) plausible for Archean oceans (Fe: 0.02–0.5 mM, Holland, 1973; Morris, 1993), the concentrations used in the Eickhoff study were constrained by experimental limitations; a concentration of 2 mM Fe(II) was used for microbial Fe(II) oxidation to grow enough cells for a visualization by scanning transmission X-ray microscopy (STXM), and 10 mg/L (170 μM) Ni was required for STXM analysis (Hitchcock et al., 2009), a factor of 400 greater than the values extrapolated by Konhauser et al. (2009) and ∼104 greater than modern seawater. Using these particular concentrations excludes a Si/Fe ratio relevant for Archean oceans, which makes it difficult to draw direct conclusions about partitioning of Ni considering the strong inhibitory role Si appears to play during trace element sorption to ferrihydrite (Konhauser et al., 2007b, 2009).
In addition to considerations based on the presence of silica and organic matter, others have suggested that IF may have been dominated by precursor mineral phases other than ferrihydrite with significantly different Ni sorption behavior. Based on field observations from the stratified and ferruginous water column in Lake Matano in Indonesia, Zegeye et al. (2012) suggested that carbonated green rust constituted an important shuttle for water column Ni to lakebed sediments. From their data, they suggested that green rust is three times more efficient at sorbing Ni compared to ferrihydrite. When considering this possibility, and its applicability to the interpretation of the IF record, several factors must be considered. First, green rust is not an abundant ferric iron precipitate today. At modern marine hydrothermal vents (e.g., Loihi) and in acid mine drainage settings, where Fe(II) oxidation is widespread, the dominant ferric oxyhydroxide mineral precipitated is ferrihydrite, with lesser amounts of goethite (e.g., Toner et al., 2012, and Clark et al., 1997, respectively). Green rust occurs only as a transient phase during the experimental formation and reduction of ferric oxyhydroxides (e.g., Cornell and Schwertmann, 2003, and Parmar et al., 2001, respectively) and Fe(II) oxidation by some nitrate-reducing bacteria (Pantke et al., 2012), as well as in some soils (e.g., Trolard et al., 1997), and at the transitional boundary in Lake Matano (Zegeye et al., 2012), but it is certainly not abundant in nature. Second, there is no direct evidence to support the formation of green rust in the Precambrian oceans, other than potentially as a transient phase, despite the ocean composition being fundamentally different than today. Recent work on the IF record has supported the assertion that the primary precipitates forming IF were likely ferric oxyhydroxides, in particular ferrihydrite (Ahn and Buseck, 1990; Bekker et al., 2010; Sun et al., 2015), although a recent study by Rasmussen et al. (2015) suggests that the primary precipitates comprised a ferrous-ferric silicate phase, such as greenalite. The Fe-isotope signatures of some Archean and Proterozoic IF suggest a biogenic or abiogenic Fe-oxide or oxyhydroxide mineral phase, rather than green rust, ankerite, or siderite (Planavsky et al., 2012). Third, both freshwater and marine photoferrotrophs that oxidize Fe(II) during anoxygenic photosynthesis and may have been ultimately responsible for iron formation generally produce poorly crystalline ferrihydrite, even under fully anoxic conditions (e.g., Kappler and Newman, 2004; Eickhoff et al., 2014; Wu et al., 2014). Fourth, when the partitioning coefficients of Zegeye et al. (2012) are evaluated in the face of Ni/Fe ratios preserved in the IF record, there is an obvious and fundamental disconnect with respect to the partitioning scenario occurring in Lake Matano (Fig. 2). Dissolved Ni concentrations at the studied depths in Lake Matano (31–50 nM) are broadly comparable to modern oceanic values (∼10 nM). However, if the partitioning scenarios observed in Lake Matano were applied to the IF record, Archean Ni concentrations would have to have been some 100–1000 times lower than modern oceanic values. This scenario is difficult to reconcile with the general secular evolution of both Earth's crust and sediments from Ni-rich to Ni-poor over geological time (Konhauser et al., 2009).
The authors attribute this disconnect in part to the deficiency of silica in Lake Matano (300–420 μM; Crowe et al., 2008) relative to Precambrian oceans, whose presence results in depressed Ni sorption (Konhauser et al., 2009). The experimental partitioning data of Konhauser et al. (2009), performed for several different silica concentrations (670 μM, saturation with respect to cristobalite, and 2200 μM, saturation with respect to amorphous silica), indicate that dissolved silica at levels approaching saturation with respect to amorphous silica suppresses Ni sorption by a factor of approximately 10 in terms of the KD value relative to the silica-free condition (0.00123 μM −1 at 2.20 mM Si versus 0.01341 μM −1 at 0 mM Si; Konhauser et al., 2009). This highlights the necessity for considering elevated silica concentrations when examining the partitioning of Ni into Fe minerals. Furthermore, Jones et al. (2015) estimated Precambrian silica concentrations on the order of 600–1500 μM, the lower end being not too dissimilar to concentrations currently present in Lake Matano. Thus, the presence or absence of silica alone does not appear to account for the drastic difference in Ni sorption efficiency between green rust in Lake Matano and synthetic ferrihydrite. Further experiments are warranted, but the simplest explanation may be that green rust played little role in trace element shuttling during Archean IF deposition.
3. Nickel Mobility during Diagenesis
One potentially troubling aspect regarding the interpretation of the IF record of evolving paleomarine Ni is the secondary mobilization of Ni from the ferrihydrite-rich sediments during burial. In this regard, Frierdich et al. (2011) suggested that IF were unsuitable proxies for Ni concentrations through time as interactions with Fe(II)-rich fluids could have stimulated the release of Ni; this is a function of redox-driven recrystallization of Fe minerals exposed to ferrous iron. For instance, Frierdich and Catalano (2012) experimentally showed a release of 2.0% and 11.2% of Ni from hematite and goethite, respectively, when Ni-substituted minerals reached equilibrium with Fe(II) solutions. A greater magnitude of release was observed when Ni-substituted hematite and goethite were sequentially exposed to fresh Fe(II) solutions every 2 weeks for up to seven cycles. Under this scenario, 6.6% and 27.6% of Ni were released from hematite and goethite, respectively. When IF are considered open systems, as in the work of Frierdich et al. (2011) and Frierdich and Catalano (2012), then it is certainly possible that the IF record only a fraction of the initially sorbed Ni. Importantly, this implies that the resultant paleomarine concentrations extrapolated from the KD values determined by Konhauser et al. (2009) would be too low and that decreases in the marine Ni reservoir over time were even more dramatic than as originally proposed by Konhauser et al. (2009).
If we consider the upper and lower limits of Ni release observed, 27.6% and 2%, respectively, we can apply a correction factor to the Ni/Fe ratios in the IF record, in order to account for this loss. Considering these two extremes, a factor of between 1.38 and 1.02 can be surmised. In the case of 27.6% of Ni being mobilized, this would be the equivalent to dividing either the Ni or Ni/Fe value from the respective analyses by 0.724 (i.e., 1/0.724 = 1.38) or multiplying by the corresponding correction factor (i.e., 1.38). Conceptually, this would treat the IF record as representing only 72.4% of the initial Ni, with the other 27.6% being mobilized, prior to the application of a correction factor. Accordingly, this factor becomes lower when one considers that less Ni is mobilized, converging to 1.00 as the amount mobilized decreases. Given that the Ni record in IF shows ∼40× decrease in paleomarine-dissolved Ni concentrations from the Precambrian to the modern, a correction factor between 1.38 and 1.02 imparts relatively little difference. For example, at 400 nM Ni, the highest concentrations predicted for the Archean by Konhauser et al. (2009), the corrected concentration becomes 552 nM. If anything this would imply that the previous estimate of Konhauser et al. (2009) and those presented in this article are more conservative than previously argued, a possibility in line with the work of Eickhoff et al. (2014) discussed above.
It is crucial to point out that experiments reacting Ni-substituted Fe-minerals with Fe(II) solutions (e.g., Frierdich et al., 2011; Frierdich and Catalano, 2012) were performed in trace element–free solutions. This is a situation far removed from ferrihydrite forming in equilibrium with seawater or subsequently transforming in contact with sediment pore-waters. When precipitating particles sorb Ni in a similar manner to KD experiments, particle composition (and correspondingly, Ni/Fe ratios) are determined by equilibrium with surrounding fluids. Placing Ni-substituted hematite or goethite into a Ni-free Fe(II) solution represents a system at disequilibrium and would clearly promote mobilization of Ni by leaching. The importance of this can be viewed by comparing the mobilization of Ni by Fe(II) solutions when the fluid is periodically replaced and when the solution remains unchanged over the length of the experiment (Frierdich and Catalano, 2012). This comparison suggests that the magnitude of Ni mobilized is significantly reduced when some Ni is present. Once Ni is released from either goethite or hematite, the Fe(II) solutions reach equilibrium with the mineral phases, and no further release of Ni is observed. Under such a case only ∼2–11% of Ni is mobilized and, viewed in light of the IF record, has little effect on extrapolated paleomarine Ni concentrations. Furthermore, this redox-driven recrystallization becomes more complicated when impurities in the minerals are considered. Frierdich et al. (2012) showed that the mobilization of Ni is attenuated when the mineral structure contains other impurities, such as Sn, Al, or Cr; this scenario is more realistic in terms of Fe particles precipitating in the Precambrian oceans and sinking to the bottom sediments to ultimately become IF.
The remobilization of Ni from IF sediment has been recently examined in a set of experiments designed to consider diagenetic mineral transformation. As discussed above, given the likely catalytic role that microbes played in Fe(II) oxidation during IF precipitation, the organic remains and ferric oxyhydroxide particles would together have settled through the water column and become deposited at the seafloor (e.g., Li et al., 2011). Given that the bulk water column was anoxic, perhaps with the exception of an upper layer of oxygenated waters above the chemocline of a stratified ocean (e.g., Klein and Beukes, 1989), ferrihydrite would have represented a favorable electron acceptor for the oxidation of the cellular remains. Coupling the reduction of Fe(III) minerals to the oxidation of organic matter not only explains the low content of organic carbon in the IF (<0.5%; Gole and Klein, 1981), but it also explains the prevalence of light carbon isotope compositions in associated carbonate minerals (Perry et al., 1973; Walker, 1984; Baur et al., 1985), the light iron isotope compositions of magnetite and siderite (e.g., Johnson et al., 2008; Craddock and Dauphas, 2011), and the textures of the reduced iron phases (e.g., Kolo et al., 2009; Köhler et al., 2013; Y.-L. Li et al., 2013). In this regard, Posth et al. (2013b) and Köhler et al. (2013) both demonstrated that by incubating ferrihydrite at 1.2 kbar and 170°C for 14 days the mineralogical transition from ferrihydrite to hematite is complete. The addition of organic carbon to these diagenetic experiments led to the formation of reduced Fe minerals in the post-diagenetic capsules, mainly a mixture of hematite, magnetite, and siderite, a composition broadly consistent with the modern-day mineralogy of IF.
Robbins et al. (2015) took this one step further, demonstrating that, during these diagenetic mineral transformations, Ni remains largely immobile. In fact, more than 91% of Ni was retained following the mineral transformations characteristic of these diagenetic experiments. Post-diagenetic analysis of the particles indicated a mineralogy dominated by hematite with Ni strongly retained in the particles. Critically, those experiments captured the dehydration of ferrihydrite to hematite, which can occur at temperatures below 140°C (Becker and Clayton, 1976). Ni was retained on both biogenic ferrihydrite and abiogenic ferrihydrite, in the presence and absence of glucose. Both biogenic ferrihydrite and abiogenic ferrihydrite, in the presence of glucose, should provide ideal conditions for Fe(III) reduction and mobilization of Ni via redox-driven reactions. When considered alongside other evidence for authigenic seawater signatures in IF, such as the preservation of marine-like rare earth element (REE) patterns (e.g., Bau and Dulski, 1996; Bolhar et al., 2004; Hofmann, 2005; Pecoits et al., 2009; Haugaard et al., 2013), these findings suggest IF should be a robust indicator of paleomarine concentrations, and by extension Ni/Fe ratios, even after diagenetic mineral transformations. Moreover, δ18O isotope heterogeneity in BIF hematite grains is thought to reflect near-primary depositional features, and clearly contrasts with the limited range of O isotope compositions in quartz, which are best explained by homogenization through recrystallization and exchange with secondary fluids (W. Li et al., 2013). Similarly, Fe isotope heterogeneity between hematite and magnetite layers is interpreted to reflect variable extents of oxidation and precipitation from aqueous Fe(II) of variable isotopic composition, depending on the proportion of hydrothermal versus diagenetic input, and as such, are considered not to have been overprinted by post-depositional processes (Frost et al., 2006; Steinhoefel et al., 2010; Y.-L. Li et al., 2013).
Finally, it is important to stress that the observations made by Konhauser et al. (2009) were based on temporal trends, that is, changing dissolved Ni concentrations in seawater from the Archean through to the present. Irrespective of whether extrapolated paleomarine Ni concentrations were overestimated (e.g., Zegeye et al., 2012) or underestimated (Frierdich et al., 2011; Frierdich and Catalano, 2012; Eickhoff et al., 2014), the key observation is the drop in dissolved Ni concentrations between 2.7 and 2.5 Ga and the fact that Ni concentrations in IF prior to 2.7 Ga were never reached again. The IF record is clear on this, suggesting that the authigenic marine nature of this signal remains undisputed. Suppressed methanogenesis remains a likely outcome of diminishing Ni concentrations in seawater (see section below).
4. Archean Methanogen Famine Revisited
It has been postulated that, as long as biogenic methane production was high, oxygen could not accumulate in the atmosphere because of the high chemical reactivity of the two gases (e.g., Zahnle et al., 2006). In their model, oxidative weathering of terrestrial sulfide minerals increased the supply of dissolved sulfate
Without question, much of the nickel famine hypothesis rests on the assumption that methane levels were high in the Archean (e.g., Kasting and Siefert, 2002) and in some way a determinant for the timing of the rise of atmospheric oxygen. Although abiogenic supplies of methane were likely minimal during the Hadean due to a relatively high early mantle oxygen fugacity (e.g., Wade and Wood, 2005), when methanogenic bacteria evolved, perhaps as early as 4.1 Ga (Battistuzzi et al., 2004), they provided a flux of methane sufficient to sustain an atmospheric CH4 mixing ratio of >10−3, or some 1000 times that of today (Pavlov et al., 2000). Evidence for biogenic methane accumulation exists in Paleoarchean chert-hosted fluid inclusions (Ueno et al., 2006) and in marine sediments, where δ13Corg values of Neoarchean kerogen are often lower than −40‰ and can reach −60‰ (Eigenbrode and Freeman, 2006). Such isotopically light marine organic matter either reflects the primary productivity of methanogens (e.g., Londry et al., 2008), incorporation of methane into biomass of methanotrophic bacteria (Hayes et al., 1983), or the formation of a 12C-enriched organic haze by photolysis of methane (Pavlov et al., 2001).
Our initial hypothesis also necessitates that seawater Ni concentrations did indeed decline after 2.7 Ga. This view has merit because a hotter Archean mantle (e.g., Berry et al., 2008) led to more mantle source melting, which produced abundant Ni-rich komatiite and olivine-rich basalt as new seafloor, seamounts, and continental crust (e.g., Kamber et al., 2005). Second, Archean tholeiites, which make up most of the oceanic crust, were richer in Ni than modern tholeiites (∼140 vs. ∼90 ppm; Arndt, 1991). Third, a major peak in preserved komatiite abundance occurred during the most intense period of mantle plume magmatism and continental crustal growth recorded in Earth's history, between 2.72 and 2.66 Ga (Condie, 1998; Barley et al., 2005), after which their abundance decreased rapidly. Ultimately, the hotter Archean mantle produced more Ni-rich ultramafic rocks, resulting in a much greater supply of Ni to the oceans than any time thereafter (Smithies et al., 2005).
Perhaps it is no coincidence that the demise of methanogens and the rise of cyanobacteria are linked. In a competition for resources, if one microorganism does well, the others are held in check through competitive exclusion (e.g., Hibbing et al., 2010; Weber et al., 2014). Microbes can further enhance their competitiveness through the production of antimicrobial compounds or by having evolved a better resistance to some external chemical. In the case of methanogens, although no studies to date have shown any pathogenic species, they do have toxin/antitoxin systems (Cavicchioli et al., 2003), while a study of acidic peat bogs showed that the addition of the antibiotic rifampicin inhibited growth of acetogens without affecting the growth of methanogens (Bräuer et al., 2004). In either case, it is possible that methanogens pre-2.7 Ga were widely dispersed throughout the oceans, from below the photic zone of the water column to the bottom sediments, with photoautotrophs dominating the surface. Based on the predicted iron:phosphorous (Fe:P) ratios in evolving Precambrian upwelling waters, Jones et al. (2015) suggested that photoferrotrophs would likely have faired better than the cyanobacteria as long as the oceans were highly ferruginous (with Fe:P >424, as based on the Redfield ratio), but as soon as dissolved Fe:P fell below 424, iron would be exhausted, leaving P to fuel other modes of photoautotrophy, such as oxygenic photosynthesis. Similarly, Swanner et al. (2015) experimentally demonstrated that highly ferruginous waters were toxic to cyanobacteria because high Fe2+ concentrations caused intracellular reactive oxygen species to accumulate within the cell, leading to decreased photosynthetic efficiency and lower growth rates. Rampant methanogenesis, whether supported by higher H2 fluxes or organic matter produced by anoxygenic phototrophy, would have contributed to the maintenance of the largely anoxic conditions. In turn, this permitted widespread Fe(II) dispersal that would have prevented cyanobacterial proliferation should the Fe(II)-related brakes identified by Jones et al. (2015) and Swanner et al. (2015) prove critical.
In addition to the drop in seawater Ni concentrations after 2.7 Ga, other events around this time would have converged to cause a contraction in methanogens while simultaneously facilitating the rise of cyanobacteria. First, the increasing supply of solutes and sediments from weathering of the newly emergent continents (e.g., Lalonde and Konhauser, 2015) led to the formation of stable continental platforms. This, in turn, facilitated the diversification of stromatolites that utilized these shallow-water settings (McNamara and Awramik, 1992) and the evolution of new types of cyanobacteria with greater motility and faster growth rates, which could better compensate for continuous accretion (Knoll, 1984). With increased consumption of CO2, indiscriminate sedimentation of massive amounts of calcium carbonate took place in situ on the seafloor, leading to decimeter- to meter-thick beds that extended over thousands of square kilometers (Grotzinger and Knoll, 1999). Moreover, increased nutrient supply would also have had the effect of increasing the productivity of planktonic cyanobacteria, as envisioned for increased phosphorus delivery to the ocean following the snowball Earth glaciations in the Neoproterozoic (Planavsky et al., 2010). The end result was that cyanobacteria now expanded throughout the shallow littoral zones, and the oxygen they produced would have marginalized the methanogens that were already in decline.
Second, there should have been waning fluxes of H2 (methanogen's preferred reductant) to the bulk ocean waters with decreasing serpentinization of the seafloor (i.e., 2Mg1.8Fe0.2SiO4 + 2.933H2O → Mg2.7Fe0.3Si2O5(OH)4 + 0.9Mg(OH)2 + 0.033Fe3O4 + 0.033H2), accompanying the decrease in production of olivine-rich oceanic crust at the end of the 2.7 Ga mantle plume event and the concomitant growth and stabilization of continental cratons (Barley et al., 2005). Kump and Barley (2007) argued that the shift from predominantly submarine volcanism to more subaerial volcanism further decreased abiotic H2 fluxes since the latter yield more oxidized gases (H2O, CO2, and SO2). A loss of H2 would likely have contributed to a declining biogenic methane supply (Kharecha et al., 2005) even if methanogens continued to grow directly on seafloor crust—as they do today—where the source of both Ni and H2 remained locally abundant (e.g., Brazelton et al., 2006; Amend et al., 2011).
Third, seawater concentrations of Fe2+ declined episodically on the shallow regions of the continental shelves due to the cessation of plume-driven submarine volcanism and a lowering of sea level. On the semi-restricted carbonate platforms, such as that represented by the ca. 2.65–2.46 Ga Campbellrand Subgroup in South Africa, lower concentrations of Fe2+ allowed for “oxygen oases” to develop. In fact, stromatolites persist on kilometer-scale outcrops, suggestive of an active marine biosphere (Sumner, 1997). The presence of dissolved oxygen in seawater along the slope is inferred based on shifts toward heavier δ15N values in shale organic carbon that are consistent with the onset of an oxic N cycle (Godfrey and Falkowski, 2009), varying abundances in rhenium (Re) and molybdenum (Mo) in shales that reflect cycling of these elements in oxygenated seawater (Kendall et al., 2010), and coupled Fe and Mo isotopes that indicate the formation of ferric oxyhydroxides in the upper water column that were likely linked to cyanobacteria rather than photoferrotrophs (Czaja et al., 2012). The platform had rimmed margins (Sumner and Beukes, 2006), which could have sheltered cyanobacteria in lagoons from the delivery of Fe(II)-rich water that was sourced from depth, based on the deposition of Fe-rich shales along the slope and IF deeper in the basin (Klein and Beukes, 1989). Geochemical indicators of oxygen appear during a sea-level regression, but disappear during transgressions (Godfrey and Falkowski, 2009), and are temporally limited, lasting perhaps tens of millions of years (Sumner and Beukes, 2006).
The transition from a methane-dominated ocean to a cyanobacteria-dominated ocean is perhaps best evidenced in the 150-million-year section (2.72–2.57 Ga) of late Archean shallow- and deep-water sediments of the Hamersley Basin in Western Australia. Eigenbrode and Freeman (2006) reported δ13Corg values for 175 kerogen samples, with a range of values between −57‰ and −28‰. The most isotopically depleted kerogens are from the oldest samples, and the magnitude of their depletion is best accounted for by the biological processes of methane assimilation and oxidation. Methane oxidation requires electron acceptors, such as O2, nitrate
As the GOE commenced ca. 2.45 Ga, aerobic chemolithoautotrophic weathering of pyrite on land produced significant acidity, which would have increased sulfate and nutrient fluxes to seawater (Reinhard et al., 2009; Konhauser et al., 2011). With increased sulfate supply to the oceans, marine SRB would have begun to dominate the anoxic regions of the water column, leading to the further marginalization of methanogens to the sulfate-poor bottom sediments, as is the case today. The end result of the progressive decline in methane production (a potent greenhouse gas) was a cooling climate, eventually culminating in the onset of major Paleoproterozoic glaciations (Zahnle et al., 2006).
Kasting (2013) suggested that a decline in methanogenesis could not have led to the progressive oxidation of Earth's surface. According to his model, if the methane that was being produced originated from atmospheric H2 (Reaction 4), then a decrease in methanogenesis should have decreased atmospheric total hydrogen mixing ratio, thereby decreasing both H2 loss and O2 production. If the methane instead originated from fermentation and methanogenesis within seafloor sediments, the net reaction would have been as in Reaction 5. Each mole of CH4 that was produced should have led to 2 mol of escaping H2 (Reaction 6; see Catling et al., 2001), but the 2 mol of CH2O that would have been buried had the methanogens not been there would have been the equivalent of 4 mol of H2 escaping (see Kasting, 2013, Reaction 10); thus O2 production is again decreased. So, according to Kasting (2013), this mechanism cannot have contributed to the rise of atmospheric O2.
This model makes two assumptions whose validity is worth examining. First, in the case of H2-driven methanogenesis, it is assumed that atmospheric H2 is consumed that would have otherwise contributed to the total atmospheric mixing ratio of hydrogen. If instead volcanic or mantle H2 tended to be fixed in organic matter (for example, by photoferrotrophy; Croal et al., 2009) and buried without methanogenesis, the net effect would be to prevent the transfer of H2 and CH4 to the atmosphere, such that methanogenesis becomes an important driver of hydrogen escape, not a brake (Catling et al., 2001). The second assumption, for the case of CH2O methanogenesis, is related in that the reducing power that generated the organic matter undergoing methanogenesis comes at the expense of atmospheric H2. If methanogenesis releases the equivalent of 2 mol H2 to the atmosphere to participate in hydrogen escape, but its absence results in the burial of 4 mol H2 equivalents that would not otherwise contribute to hydrogen escape (for example, burial of biomass produced by non-H2-based phototrophy), then burial is the more efficient pathway for generating oxidizing equivalents at Earth's surface, and any decrease in hydrogen escape arising from depressed methanogenesis becomes subordinate.
The landscape of possible metabolic and reactive transport scenarios is admittedly complex, and the ensemble of consequences of depressed methanogenesis in Earth's redox balance and the GOE itself warrants further scrutiny. To complicate matters, increases in primary productivity (e.g., the first widespread oxygenic photosynthesis in open oceans) would stimulate methanogenesis and atmospheric methane mixing ratios, plausibly to the point of organic haze formation (Domagal-Goldman et al., 2008; Zerkle et al., 2012; Kurzweil et al., 2013). Such a haze would be bi-stable and self-supporting by positive feedbacks, requiring a significant drop in methane fluxes to transition away from a thick haze (Zerkle et al., 2012; Kurzweil et al., 2013). The evolving nickel scarcity provides a directional, compelling hypothesis for why methane fluxes might drop in the face of the ramping biological productivity often evoked for the GOE.
5. Conclusion
An expanded database of Ni/Fe ratios in iron formation through geological time reaffirms a profound and unidirectional shift in nickel availability during the late Archean. Extrapolation from rock record data to paleomarine dissolved Ni concentrations in ancient seawater, for example, using partition coefficient approaches, is not necessarily straightforward; and there exist a variety of caveats, mostly pertaining to the system specificity of partition coefficient approaches (e.g., mineralogy, the role of biomass and competing ions during sorption) as well as potential diagenetic effects. In reviewing these caveats, we conclude that some are not necessarily applicable to Precambrian oceans, others may have faults stemming from experimental design, while others may be of little impact to reconstructed Ni concentrations. That said, caution is welcome when extrapolating elemental concentrations from the rock record to ancient seawater, and we detail approaches for evaluating and potentially tuning trace element partitioning scenarios in light of the rock record data itself. While the biospheric effects of the nickel famine are likely multiple and more nuanced than the scenario originally envisioned by Konhauser et al. (2009), especially with respect to global redox balance and the possibility of organic haze in the atmosphere, all experimental and rock record evidence to date reaffirm that marine Ni concentrations dropped precipitously and unidirectionally between 2.7 and 2.5 Ga, never to recover.
Footnotes
Acknowledgments
Funding was provided by the National Sciences and Engineering Research Council of Canada (NSERC) to Kurt O. Konhauser; a Vanier Canada Graduate Scholarship to Leslie J. Robbins; the UnivEarthS Labex program at Sorbonne Paris Cité to Ernesto Pecoits; LabexMER to Stefan V. Lalonde; the European Research Council to Andreas Kappler; and NERC (Natural Environment Research Council) and STFC (Science & Technology Facilities Council) to Caroline Peacock.
