Reconstruction of the effect of biogeochemical evolution on BIF across geological timescale is extremely difficult due to the uncertainty of changes in key environmental factors. For example, oxygen level is a crucial factor determining whether Fe(II) oxidation is influenced by oxic or anoxic processes. However, it is difficult to determine the evolution of oxygen level in paleo-oceans with heterogeneous chemical structures, despite the widespread consensus of an overall rise of oxygen in the atmosphere of the early Earth. In comparison, temperature cooling is more predictable. From a long-term geological perspective, temperature cooling of the paleo-oceans is an indisputable fact, although the absolute magnitude of cooling is still controversial. Nevertheless, the cooling trend of temperature of the Archean-Proterozoic oceans would greatly reshape the structure and functional composition of microbial communities, thereby altering biogeochemical processes in the ocean, including the N and Fe cycling and BIF deposition. In the present work, by integrating all the results, we developed a four-stage evolution scenario for the coupled N-Fe redox processes, providing a possible explanation for the secular changes of iron mineral precipitation in BIFs in the Archean-Paleoproterozoic ocean (Fig. 4):
Small scale BIFs (i.e., with typical laminated Si-rich and Fe-rich layers) began to deposit as early as 3.8 Ga. ago.50 The ocean was hot (most probably higher than 55 oC, even over 75 oC)1, 47, 51, highly reducing, and almost devoid of any microbial activity (if any). Iron minerals in these oldest BIFs are overwhelmingly dominated by magnetite, accompanied by some Fe-silicates, such as grunerite, but few contain hematite52. During this prebiotic era, abiotic N2 fixation was the main process of introducing a low flux of fixed N, (initially NO2- ) (~ 106-1010 mol N∙year-1 or even higher), into the ferruginous ocean30, 53. The NO2- would be rapidly reduced through chemical oxidation of Fe(II) especially at high temperature, followed by precipitation of iron in the form of magnetite (2NO2- + 9Fe2+ + 16HCO3-→N2 + 3Fe3O4 + 8H2O + 16CO2, as shown in Fig. S6-c), and/or cronstedtite in the presence of Si. Therefore, although NO2- would not be accumulated to a considerable concentration, its abiotic supply was sufficient to precipitate a large proportion of the oldest BIFs (totally ~ 1011-1012 mol Fe in the Paleoarchean BIFs estimated by the BIF deposit size and average 20 wt% Fe content, Fig. 3a). Moreover, the primary production of Fe(II)-Fe(III)-mixed minerals (magnetite, silicates) eases the vexatious requirement of a reduction step of ferric iron precursors (e.g., ferrihydrites), as proposed previously54. This scenario offers a reasonable explanation for the higher abundance of magnetite in the early BIFs7.
After the Mesoarchean, deposition of BIFs increased slightly, until the Neoarchean. Magnetite is still the predominant iron mineral in most BIFs (particularly > 2.9 Ga)7. During the Mesoarchean, biological N2 fixation bloomed and greatly enlarged the NH4+ flux to the ocean9, 10. Microbial oxidation of NH4+ may also emerge as early as in the Mesoarchean55, as AOA could have oxidized ammonia at extremely low levels of O2 or even in the anoxic water 15, 16. Moreover, local “oxygen oasis”, due to the emerging oxygenic photosynthesis55, may have gradually stimulated NH4+ oxidation12. Therefore, microbial nitrification provided additional NOx- flow, mainly NO2-, as the temperature was still likely much higher than modern oceans (over 50 oC, as indicated in Figs. 2a and 4). With the increase of NO2- flux, chemical oxidation of Fe(II) by NO2-, which primarily forms magnetite and cronstedtite, was still the predominant process of coupled N-Fe processes. This scenario offers an explanation for the increase of magnetite-dominant BIFs volume (totally ~ 1014-1015 mol Fe in Mesoarchean BIFs as estimated by the BIFs deposit size and average 20 wt% Fe content, Fig. 3b)6.
BIF deposition increased significantly after the Neoarchean, reaching a peak from the mid-Neoarchean to the early Paleoproterozoic (2.7–2.4 Ga), right before GOE. Intriguingly, hematite content increased dramatically in these BIFs, especially after 2.6 Ga, suggesting that certain important transition of biogeochemical processes might have occurred and significantly affected the Fe redox cycle. Indeed, the oxygenation of surface water was more pervasive at this time, which greatly stimulated aerobic nitrification and further increased the NOx- flux. For example, the extremely positive excursions in δ15N was observed in the 2.7–2.5 Ga sediments, which was interpreted as a great rise in N loss via coupled nitrification-denitrification processes56, 57. Therefore, the abundant NOx- might be extensively reduced by Fe(II) oxidation considering the super-large BIF deposition during this period. However, as temperature continued to decrease (lower than 50 oC), microbial nitrification of NH4+ would produce a mixture of NO2- and NO3- with an increasing share of NO3-. With the increased NOx- pool/flux, Fe(II) oxidation coupled with microbial NO3- reduction and chemical NO2- reduction produced a mineral assemblage of magnetite, silicates, and Fe(III) hydroxides (i.e., goethite, lepidocrocite, Fig. 3c). Over time at certain temperature and pressure, these ferric hydroxides can transform to hematite58. Furthermore, the relative dominance of Fe(III) hydroxides produced by N-Fe redox coupling would gradually increase at lower temperatures (Table 1), which offers a plausible explanation for the transition of magnetite-dominant to hematite-dominant BIFs during this period.
A dramatic decrease or even stop of BIF deposition was observed after 1.8 Ga when the temperature further cooled down, suggesting a quiescence of Fe(II) oxidation at this time (Fig. 3d). Two possible reasons may account for this transition. First, canonical nitrification-denitrification cycles might outcompete the coupled N-Fe cycles, due to rapid and complete oxidation of NH4+ (into NO3-) and complete reduction of NO3- to N2O and N2 gases (Fig. S1d, S4c) at low temperatures. Second, the development of an euxinic zone (S2- rich layer, caused by increased sulfate input from oxidative weathering of continental sulfides) between the nitrogenous and ferruginous layers provided an important reductant (i.e., S2-) to consume NOx- and/or precipitate Fe(II)59, both of which would weaken the coupling of the N-Fe cycles. This scenario offers a possible explanation for the extremely low BIFs during the Meso- and Neoproterozoic (without considering the Rapitan IFs deposited during the Neoproterozoic).
In summary, our proposed evolution model of the coupled N-Fe cycling processes, as a response to temperature decrease, provides a reasonable explanation for the successive change of abundance and mineralogy of BIFs. An important reason for the overlooked role of the coupled N-Fe cycles in the BIF formation in most previous studies is the assumption of the low accumulation of oxidative NOx- species in the reducing Precambrian ocean. However, we argue that the increasing flux and changing speciation of NOx- are extremely important for triggering extensive Fe(II) oxidation and mineral formation. A recent study found that the networked N-Fe reactions would provide substantial N2O emission and N burial (by association with Fe minerals) in the Precambrian marine ecosystem60, supporting the pervasive N-Fe reactions. Here we further demonstrated that N and Fe cycles could also be directly coupled to affect Fe mineralogy. In addition, our observed iron mineral assemblages formed by the dynamically evolving N-Fe coupling also have some other under-appreciated geological and environmental consequences, such as the bioavailability of N (via change of N speciation) and P (via adsorption to iron minerals) overtime, thereby controlling primary productivity and emission of greenhouse gasses N2O and CO260. Finally, further work is required for additional assessment of the contribution of N-Fe coupling to BIF deposition, which include, but not limited to, the competition by different Fe(II) oxidation pathways, quantitative impacts of other evolutionary factors (e.g., pH, Ca/Mg), and the N and Fe isotopic footprints recorded in geological archives in response to temperature gradients.