Iron Isotopes Reveal a Benthic Iron Shuttle in the Palaeoproterozoic Zaonega Formation: Basinal Restriction, Euxinia, and the Effect on Global Palaeoredox Proxies

: The Zaonega Formation in northwest Russia (~2.0 billion years old) is amongst the most complete successions that record the middle of the Palaeoproterozoic era. As such, geochemical data from the formation have played a central role in framing the debate over redox dynamics in the aftermath of the Great Oxidation Event (GOE). However, uncertainty over local redox conditions and the degree of hydrographic restriction in the formation has led to contradictory interpretations regarding global oxygen (O 2 ) fugacity. Here, we provide new iron (Fe) isotope data together with major and trace element concentrations to constrain the local physiochemical conditions. The Zaonega Formation sediments show authigenic Fe accumulation (Fe/Al (cid:29) 1 wt.%/wt.%) and δ 56 Fe ranging from − 0.58‰ to +0.60‰. Many of the data fall on a negative Fe/Al versus δ 56 Fe trend, diagnostic of a benthic Fe shuttle, which implies that Zaonega Formation rocks formed in a redox-stratiﬁed and semi-restricted basin. However, basin restriction did not coincide with diminished trace metal enrichment, likely due to episodes of deep-water exchange with metal-rich oxygenated seawater, as evidenced by simultaneous authigenic Fe(III) precipitation. If so, the Onega Basin maintained a connection that allowed its sediments to record signals of global ocean chemistry despite signiﬁcant basinal effects.


Introduction
The Proterozoic eon (2500-541 Ma) hosted some of the most fundamental global environmental transitions in Earth's history, including the multi-step rise of oxygen in the atmosphere [1]. The most pronounced amongst these are the Great Oxidation Event (~2502-2320 Ma) [2][3][4], the Lomagundi Event, which has been proposed to coincide with an "O 2 overshoot" (~2200-2060 Ma) [5][6][7], and the terminal oxygenation of the atmosphere in the Neoproterozoic Oxidation Event~800-550 Ma [8]. Intermixed with these positive shifts in atmospheric O 2 fugacity is a retrograde redox shift-an "O 2 crash"-that follows the Lomagundi Event and the O 2 overshoot [9,10] (though see Ohmoto et al. [11] for a contradictory view). The timing, extent, and environmental significance of this O 2 crash are contentious, largely due to conflicting interpretations based on different geochemical redox proxies [9,10,[12][13][14][15][16][17]. Resolving the redox dynamics of the post-Lomagundi time period requires new geochemical proxy data from the sedimentary rock record to be placed into a framework that takes into account both global and local elemental cycling [18].
The Palaeoproterozoic Zaonega Formation comprises a sedimentary succession that records the end of the Lomagundi Event and the O 2 overshoot [19]. As such, it has been the target of a number of attempts to assess global redox conditions using a variety of geochemical proxies, such as carbon and nitrogen isotopes [10,12,20], rare earth element patterns [12,21], sulphur isotopes [17], selenium isotopes [22,23], and the abundances and isotopes of several redox-sensitive metals [13,17,[24][25][26][27]. Many of the interpretations in those previous studies rely on the assumption that the Zaonega Formation faithfully records chemical signals from the global ocean, so that inferences from these rocks may be extrapolated to assess global conditions. However, this assumption is being increasingly challenged on the basis that the Zaonega Formation was influenced by significant syn-depositional magmatic activity resulting in hydrothermal circulation, hydrocarbon generation and seepage, basinal reconfiguration and restriction, as well as local modulation of both the sulphur and carbon cycles [12,15,16,20,21,28,29].
The redox structure of the basin and its hydrographical openness are especially crucial in interpreting the redox-sensitive element record, since trace metal drawdown rates are dependent on water column oxygen fugacity or sulphide concentration, with restricted conditions in anoxic basins resulting in drawdown and isotope fractionation of the trace metal pool [30,31]. Yet, recently published interpretations of the Zaonega Formation do not agree in this regard. On one hand, sulphur isotope modeling results [15,16] and persistent positive europium anomalies [12,21] seem to suggest a basin with poor access to the open ocean, but where anoxia and sulphidic conditions were restricted mostly to pore waters. On the other hand, selenium isotopes and trace metal abundances in the Zaonega Formation seem to require more anoxic and less restricted conditions [13,17,23,24].
In this study, we report iron (Fe) elemental ratios and isotopes from a 102-metres-thick section of the upper Zaonega Formation in the northern Onega Basin-the same section previously investigated by Paiste et al. [15], Mänd et al. [13], and Kreitsmann et al. [12]-in order to resolve issues surrounding the interpretation of palaeoenvironmental depositional conditions. Due to its redox sensitivity and importance to microbial processes, Fe geochemistry is a powerful proxy for water column redox structure, sediment biogeochemistry, and basinal configuration [32][33][34], and provides key insights into the Palaeoproterozoic Zaonega Formation and its attendant controversies.

The Fe Isotope Palaeoredox Proxy
Iron is the most abundant redox-sensitive element in Earth's sedimentary environments and its geochemistry has been routinely used in assessments of palaeoredox conditions from ancient sediments [32,34]. To investigate Fe redox cycling, a variety of chemical indices have been proposed. For example, positive deviations in the ratio of bulk Fe to bulk aluminium (Al)-a redox-insensitive element hosted in the detrital fraction of sedimentsfrom the detrital mean (0.472 ± 0.302 wt.%/wt.%) [35] can indicate authigenic enrichment of Fe in sediments, whereas negative deviations may indicate microbial Fe(III) reduction and mobilisation [36,37]. The Fe/Al ratio is commonly elevated under anoxic water masses, where Fe(II) can be drawn down into sediments through capture in sulphide minerals, so that a trend towards higher Fe/Al ratios is observed moving from shallower oxic facies to deeper anoxic ones [38].
The most frequently utilised Fe-based redox proxy is the sequential Fe extraction method and the speciation data that it yields [34,39]. Based on calibrations in modern anoxic basins (i.e., the Black Sea), the ratio of highly reactive Fe to total Fe is indicative of redox conditions in the water column overlying the sediments. Under oxic conditions, the upper limit for this ratio is 0.22, and the lower limit for firmly anoxic conditions is 0.38 [40]. Values between these thresholds are considered equivocal, since water column iron enrichments could be masked by high sedimentation rates or transformation of unsulphidised highly reactive iron into sheet silicates [40][41][42]. Furthermore, if the fraction of pyrite-bound Fe within the total reactive Fe pool exceeds 0.7-0.8, the overlying water column is considered to reflect euxinic conditions [34,40].
However, these relatively simple frameworks are complicated by factors such as bathymetry, sedimentation rate, and local terrigenous and hydrothermal inputs that can significantly affect the distribution of Fe within sediments [43,44], but which are rarely resolvable for ancient sedimentary basins. In addition, studies of modern environments have shown that Fe distribution in sediments may be influenced by a strong redoxcline in the water column rather than reflecting redox conditions of the entire sedimentary basin [45]. The effects of post-depositional processes on Fe distribution patterns in sedimentary rocks are also important to consider, since diagenesis and metamorphism readily mobilise Fe and generate proxy signatures that are unrelated to environmental redox conditions [41,46]. Thus, pairing Fe speciation with other proxies (e.g., mineralogy, petrography, trace elements, and isotope characteristics) is necessary for a more robust interpretation of palaeoenvironmental redox conditions. Iron isotopes provide a particularly powerful tool in assessing Fe redox and precipitation processes [32,33]. The Fe isotope composition of a sample is commonly expressed as the permil deviation of the ratio 56 Fe over 54 Fe from that of the reference standard IRMM-014: Redox cycling produces the most substantial Fe isotope fractionations, since 56 Fe is preferentially enriched in Fe(III) over Fe(II). Equilibrium fractionation factors between aqueous Fe(II) and Fe(III)-bearing phases, such as ferrihydrite (Fe[OH] 3 ), goethite (FeO[OH]), or magnetite (Fe 3 O 4 ), have been estimated to be between~+3‰ and~+1‰ [47][48][49][50]. This means that sediments rich in authigenic Fe(III) commonly have high δ 56 Fe values-for example, ferrihydrite produced from aqueous Fe(II) by anoxygenic iron-based photosynthesis (also known as photoferrotrophy) is enriched in 56 Fe and results in a δ 56 Fe shift of~+1.5‰ [51]. Correspondingly, if a substantial fraction of dissolved Fe(II) is drawn down by oxidation and Fe(III) oxyhydroxide precipitation, the residual Fe(II) pool can shift towards negative δ 56 Fe values, both on a basin scale [52] and globally [53].
By contrast, redox cycling within sediments produces 56 Fe-depleted Fe(II)-most commonly this occurs through dissimilatory iron reduction (DIR), which couples the microbial reduction of Fe(III) with the oxidation of organic carbon. The generated dissolved Fe(II) carries δ 56 Fe values~−1.3‰ to −2.6‰ lower than the Fe(III) substrate and can subsequently precipitate as magnetite, siderite (FeCO 3 ), or pyrite (FeS 2 ), or escape into the water column [54][55][56]. In fact, DIR is such a dominant process in modern sediments that it supplies a substantial portion of 56 Fe-depleted dissolved Fe(II) into the oceans along with hydrothermal venting (δ 56 Fe typically~−0.5-0.1‰) [57,58].
Ferrous iron drawdown through sulphide precipitation in sulphidic waters or sediments with no redox change can also induce a fractionation. This process occurs through several steps involving aqueous FeS and greigite (Fe 3 S 4 ) and commonly results in the formation of pyrite [59,60]. The net fractionation factor for this process is typically considered to be negative, producing δ 56 Fe shifts from 0.5‰ to −2‰ [61,62]. Given that equilibrium and kinetic fractionations in pyrite precipitation produce opposing isotope effects, the sign and extent of fractionation depends on the extent and rate of Fe(II) precipitation: Late-stage pyrite precipitating slowly can evolve to positive values [63].
While using pyrite δ 56 Fe values recorded in ancient organic-rich mudstones (hereafter, mudstone refers to fine-grained siliciclastic-dominated rocks) to set constraints on local and global Fe cycling is possible [32], different influences on the Fe isotope budget can produce similar δ 56 Fe shifts, making the interpretation of this record at times challenging. For example, if pyrite precipitates in the water column, δ 56 Fe values can vary based on how much of the dissolved Fe(II) pool has previously been drawn down by Fe(II) oxidation, which preferentially removes 56 Fe [52]. For pyrite produced diagenetically in the sediments through DIR, a detrital Fe(III) source will produce more negative δ 56 Fe pyrite values than an authigenic Fe(III) source. Finally, since the resulting fractionation is governed by the rate of pyrite precipitation, both local sulphide and Fe(II) availability can modulate precipitation rates and, therefore, net fractionation [63].
A highly diagnostic pattern in pyrite δ 56 Fe is present in anoxic and sulphidic basins, such as the modern Black Sea. There, a "benthic Fe shuttle" delivers Fe(II) derived through DIR from suboxic sediments on basin margins towards deeper, more anoxic parts of the basins where it is deposited nearly quantitatively as pyrite [38,64]. Since biological reduction preferentially mobilises 54 Fe, this results in a gradient towards higher authigenic iron abundance (Fe/Al) and more negative δ 56 Fe from shelf to basin. For example, oxic shelf sediments in the Black Sea carry δ 56 Fe values averaging +0.16 ± 0.02‰, while sulphidic sediments further off-shore average −0.13 ± 0.04‰ [38]. This mechanism has similarly been applied to sediments deposited during ancient oceanic anoxic events [65]. Crucially, the presence of such a trend has been shown to depend on several parameters, including the redox structure of the basin and the shelf-to-basin ratio, both of which are closely related to the degree of hydrographic openness [43]. In more open settings than the Black Sea, even a positive sedimentary δ 56 Fe-Fe/Al trend is possible as smaller fractions of dissolved Fe(II) are oxidatively recaptured under suboxic conditions [66]. Using the existing body of work on Fe-cycling processes in modern anoxic settings as a guide, it is possible to reconstruct ancient Fe cycling regimes based on pyrite δ 56 Fe data.

Geological Background
The Zaonega Formation forms part of the sedimentary succession within the >1200-km 2 Palaeoproterozoic Onega Basin in northwestern Russia, on the northern shore of Lake Onega ( Figure 1). The Onega Basin succession begins in the early Palaeoproterozoic with the deposition of terrigenous clastic sediments interlayered with lava flows in a continental rift setting [67]. Stabilisation and shallowing of the basin led to the deposition of >800 m of halite-sulphate evaporites and evaporitic carbonates of the Tulomozero Formation [68]. Following a deepening of the basin, >1500 m of mudstones, dolostones, and mafic lavas of the Zaonega Formation were deposited during a period of high volcanic activity [19]. The basin was subsequently filled by the accumulation of hundreds of metres of mafic lavas and tuffs of the Suisari Formation [69]. Deposition terminated with the emplacement of lacustrine siliciclastics of the Kondopoga Formation [67]. Subsequently, during the 1890-1790 Ma Svecofennian orogeny, the basin was metamorphosed to greenschist facies and deformed into a series of northwest-southeast-oriented folds [67].
The lithology of the Zaonega Formation is dominated by interbedded organic-rich mudstones and dolostones-the latter forming a larger portion in the lower and upper part of the succession, the former dominant in the middle portion-as well as turbidites and a smaller number of tuffs [16]. Mudstones in the formation can be exceptionally rich in organic carbon (up to 70 wt.% in some layers) [19] and pyrite (up to 16.5%) [17,70,71]. The more organic-rich layers contain some of Earth's earliest phosphorus-rich sediments [72]. The lithological and geochemical composition of the formation attests to a biologically productive basinal ecosystem and dissimilatory metabolic processes that led to sharp and fluctuating redox gradients in the sediments and water column [72]. Additionally, >35% of the formation in the northern part of the basin, and ~70% in the southern part, is comprised of mafic volcanic rocks and intrusions, attesting to a strong volcanic influence on the sedimentary environment [67]. Some of the sills have peperitic contacts with the sedimentary layers, indicating syn-depositional emplacement [70,73]. The resulting thermal gradient led to syn-depositional oil and gas generation, in turn fueling basin-wide methanotrophy in the sediments or the water column [29,74], and resulting in oil spills [19]. The generation of hydrothermal fluids further led to emplacement of silicic veins that permeate the rock in the middle part of the Zaonega Formation [15] and Additionally, >35% of the formation in the northern part of the basin, and~70% in the southern part, is comprised of mafic volcanic rocks and intrusions, attesting to a strong volcanic influence on the sedimentary environment [67]. Some of the sills have peperitic contacts with the sedimentary layers, indicating syn-depositional emplacement [70,73]. The resulting thermal gradient led to syn-depositional oil and gas generation, in turn fueling basin-wide methanotrophy in the sediments or the water column [29,74], and resulting in oil spills [19]. The generation of hydrothermal fluids further led to emplacement of silicic veins that permeate the rock in the middle part of the Zaonega Formation [15] and have led to dedolomitisation and re-equilibration of the carbon and oxygen isotope compositions on carbonate beds, whereas internal parts of some beds have remained largely intact [12,20,28].
Age constraints for the Zaonega Formation are given by a 2090 ± 70 Ma dolomite Pb-Pb age in the underlying Tulomozero Formation [75], and by carbonate δ 13 C stratigraphy: The lower Zaonega Formation records the termination of the positive Lomagundi carbon isotope excursion [76], which has previously been dated to~2060 Ma in Fennoscandia [7]. The minimum depositional age is given by several cross-cutting dykes, some of which were syn-depositional, and have U-Pb ages of 1919 ± 18 Ma [77], 1956 ± 5 Ma [78], and 1961 ± 5.1 Ma [79]. Volcanics in the overlying Suisari Formation have been dated to 1969 ± 18 Ma (Re-Os) [69] and 1988 ± 34 Ma (Pb-Pb) [80]. Additionally, a single zircon from a tuff layer in the lower Zaonega Formation yielded a U-Pb age of 1982 ± 4.5 Ma [79], whereas preliminary Re-Os dates constrain the upper Zaonega Formation to~2050 Ma [81,82]. Ultimately, these age constraints place the deposition of the upper Zaonega Formation between~2050-1988 Ma.

Sample Material
The material for this study comes from drill cores OnZaP-1 and OnZaP-3, drilled 500 m apart in the northern Onega Basin near Shunga village (62.5870 N, 34.9310 E, and 62.5920 N, 34.9280 E, respectively). Using a P-rich dolomite marker horizon that occurs throughout the Onega Basin, Paiste et al. [15] correlated the two cores to form a 102 m thick combined section. The lower part of the section from the base at 102 m to 53 m is dominated by increasingly organic-and pyrite-rich mudstones, interbedded with some dolomite-calcite beds that are cross-cut by a suite of quartz, mica, and pyrobitumen veins. At 53 m, a layer of phosphorus-rich mudstone occurs and is overlain by a massive dolomite bed up to 44 m with thin mudstone interlayers and tens-of-centimetres to several-metres-thick black silica veins. Above 33 m, the section is dominated by relatively organic-and pyrite-poor mudstones and carbonates, which show rare cross-bedding features. See Paiste et al. [15] for a detailed lithological description.
Two sets of samples were chosen for chemical analyses: Set MSP0001 (n = 134) and set MSP0010 (n = 78). Part of the elemental concentration data and methodology reported here are adapted from Paiste et al. [15] and Mänd et al. [13].

Element Concentrations
Total organic carbon content for sample set MSP0001 was determined from powdered samples at the Geological Survey of Norway in Trondheim using a LECO SC-444 analyser. The detection limit was 0.1 wt.% and the precision was better than 10%. Similar data for sample set MSP0010 were measured by loss on ignition at the European Institute for Marine Studies (IUEM) in Brest, France. Approximately 0.5 g per each sample was dried and pulverised, then loaded into ceramic beakers and combusted at 500 • C for 24 h, after which the mass loss was measured. Repeat measurements of 7 samples determined that the error was typically <1 wt.%.
Elemental concentrations and Fe isotope compositions for sample set MSP0010 were measured at IUEM in Brest, France. Samples were powdered using a tungsten carbide crusher and an agate mill. To remove organic carbon,~0.5 g of each sample was combusted at 500 • C for 24 h. For major element analysis, samples were digested overnight at 80 • C in concentrated HF and HNO 3 , then neutralised with boric acid to retain Si. Concentrations were measured on a Jobin Yvon Horiba Ultima 2 ICP-OES. External accuracy was monitored using a suite of 8 reference materials; measured values were within 2% of recommended values in all cases and relative standard deviation between repeat measurements was <2% (Table S1).
Sample preparation for trace element and Fe isotope analysis of sample set MSP0010 was done in a class 1000 clean laboratory at IUEM using only purified acids from a TFE sub-boiling distillation system. All samples were handled in acid-washed TFE and PFA beakers and vials. First, 1.5 mL of concentrated HF and 1.5 mL of concentrated HNO 3 was added to the samples, heated to 80 • C overnight, then allowed to evaporate at 80 • C. Second, 3.3 mL of concentrated HCl and 0.7 mL of concentrated HNO 3 (i.e., aqua regia) was added, heated to 35 • C for 3 h, and allowed to evaporate at 80 • C once more. Finally, samples were taken up in 4 mL of 6N HCl and stored in PFA vials.
To measure elemental concentrations, 100 mg of each sample solution was diluted in 4 mL 2% HNO 3 and spiked with indium, which was monitored to correct for instrument drift. Concentrations were measured on a Thermo Scientific Element 2 sector field ICP-MS, and external precision was assessed against BHVO-2, ShTX-1, and ShWYO-1 geostandards. For the elements discussed here, reference material measurements were within 20% of known values and relative standard deviation between repeats was <12%.

Iron Isotopes
Iron isotope measurements were performed at IUEM and IFREMER (Brest, France) on aliquots of fully digested samples (described above) following the protocol developed by Rouxel et al. [83]. First, Fe in dissolved samples was separated from other elements using column chromatography: The samples were passed through polypropylene columns loaded with 2.0 mL (wet volume) AG1-X8 anion-exchange resins. Prior to sample loading, the columns were washed with~10 mL 3N HNO 3 ,~10 mL ultrapure (18 MΩ cm −1 ) water, and~10 mL 0.24N HCl, then conditioned with 4 mL of 6N HCl. An aliquot of each fully dissolved sample calculated to contain~40 µg of Fe was then loaded onto the columns, after which the matrix was eluted with the addition of 22 mL 6N HCl. Iron was then eluted with the addition of 16 mL of 0.24N HCl and captured in TFE beakers. The eluted Fe solution was evaporated at 90 • C and the Fe re-dissolved in 4 mL of 0.28N HNO 3 .
Iron isotope compositions were determined at IFREMER using a Thermo-Scientific Neptune multicollector ICP-MS in medium or high-resolution mode. Prior to analysis, each purified Fe sample was doped with 0.02 mL of a 100 µg g −1 Ni solution with a known isotope composition. During Fe isotope measurements, Ni isotopes were monitored to correct for matrix effects [84,85]. Furthermore, instrumental mass bias was monitored and corrected for using sample-standard-bracketing, wherein each sample analysis was bookended by measurement of an IRMM-014 standard solution [86]. Repeat Fe isotope measurements of the IRMM-014 (n = 107) reference material gave an internal reproducibility of 0.07‰ to 0.08‰ (2SD). Measurements of the geostandards BHVO-2 (n = 5) and IF-G (n = 3) gave average δ 56 Fe values of +0.12 ± 0.07‰ and +0.64 ± 0.08‰, respectively, which agree, within error, with previously reported values [87].

Principal Component Analysis
Principal component analysis (PCA) was performed using the factoextra package for the R programming language [88]. The analysis was performed using 19 chemical variables, including proxies for detrital material (Al, Ti, Si, Li, Fe), hydrothermal input (Si, (Eu/Eu*) SN ), carbonate content (Ca, Mg, Mn, Li), phosphate (Ca, P), anoxia (TOC, V EF , Mo EF , U EF , Re EF ), and iron cycling (Fe, Fe/Al, δ 56 Fe). All variables were centred and normalised by their standard deviation in order to standardise their variances. Six samples from the altered 77-70 m and 53-44 m intervals (see below) were excluded from PCA analyses. Six further outliers were iteratively removed from PCA analysis, for a final total of 24 samples used (Table S2). Data analysis and visualisation employed the pandas [89], pyrolite [90], matplotlib [91], mpltern [92], and seaborn [93] packages for the Python programming language.

Geochemical Data
The results of select geochemical analyses of the OnZaP section are displayed in Vanadium varies between 7 and 3018 µg g −1 and has a median value of 245 µg g −1 (n = 206). In order to correct for the fluctuating influence of the detrital sediments, authigenic V content can be represented through enrichment factors [95]. To calculate enrichment factors, the element in question is normalised to an immobile element (commonly Al), and the ratio is further normalised to a proxy composition for the average detrital sediment (here the upper continental crust) [96]: Iron isotope ratios in the OnZaP section range from −0.58 ± 0.07‰ to + 0.60 ± 0.07‰ δ 56 Fe (errors for individual measurements are given as 2σ) and display a median value of −0.11‰ (n = 38). Below 80 m and above 35 m, values are close to the crustal mean of 0‰, but in the middle of that range, δ 56 Fe becomes more variable and negative. Values significantly above the crustal mean are concentrated in a short interval from 60 to 53 m.
Positive europium (Eu) anomalies are associated with the input of high-temperature (>250 • C) hydrothermal fluids and can be recorded in authigenic phases in sedimentary rocks. Europium anomalies are commonly calculated based on concentrations normalised to Post-Archean Australian Shale (subscript "SN") [97], and use the normalised concentrations of neighbouring rare earth elements as the expected Eu abundance (Eu*): In unfiltered samples of the Zaonega Formation, (Eu/Eu*) SN ranges from 0.54 to 2.98, with a median of 1.16 (n = 72). Grey vertical spans denote average crustal composition [35,96], yellow horizontal spans denote zones of secondary silicification [15]. Note the logarithmic scales on the TOC and VEF axes.

Principal Component Analysis
Principal component analysis is commonly used for dimension reduction in highdimensionality datasets in order to find trends in the data and to explore relationships between variables [98]. In a multi-dimensional space containing the input data with each original dimension or axis representing one input variable, PCA places new axes (principal components; PC) in a way that maximises the proportion of total data variability explained by the first PC (PC1), then does the same for the residual variability (PC2) and so on with further PCs. Principal components can, therefore, be used to visualise total sample heterogeneity in a more concise manner than is possible using the original variables   [35,96], yellow horizontal spans denote zones of secondary silicification [15]. Note the logarithmic scales on the TOC and V EF axes.

Principal Component Analysis
Principal component analysis is commonly used for dimension reduction in highdimensionality datasets in order to find trends in the data and to explore relationships between variables [98]. In a multi-dimensional space containing the input data with each original dimension or axis representing one input variable, PCA places new axes (principal components; PC) in a way that maximises the proportion of total data variability explained by the first PC (PC1), then does the same for the residual variability (PC2) and so on with further PCs. Principal components can, therefore, be used to visualise total sample heterogeneity in a more concise manner than is possible using the original variables ( Figure 3A  The results for the first four PCs, which collectively account for 83% of variability, are shown on Figure 3 and in Table 1. Among the parameters studied here, the highest amount of variability (35%, PC1) can be ascribed to a contrast between Ti, P, Li, Al, Mn, Mg, Fe, and Ca on the one hand, and Si and (Eu/Eu*)SN on the other. In other words, this component describes contrast of detrital components and some authigenic phases, such as carbonate and phosphate, against Si-richness and hydrothermal influence. A few metal enrichment factors (VEF, ReEF) and δ 56 Fe display some covariance with Si. Principal component 2 (29% of variability) is strongly controlled by indicators of anoxia and reductive element drawdown, such as TOC, metal enrichment factors (MoEF, UEF, and VEF), and Fe accumulation (Fe, Fe/Al), with minor opposite contributions from Li, Si, and δ 56 Fe. The third principal component (13% of variability) includes major contributions from Al and Fe, and opposite ones from Mg, Mn, P, ReEF, VEF, P, (Eu/Eu*)SN, and Ca, possibly representing a contrast between clay-rich and carbonate lithologies. The fourth component (6% of variability) is strongly controlled by δ 56 Fe, together with Ca, Mn, Fe, and Fe/Al, which contrast modestly with Ti, Li, P, Al, and metal enrichment factors. The results for the first four PCs, which collectively account for 83% of variability, are shown on Figure 3 and in Table 1. Among the parameters studied here, the highest amount of variability (35%, PC1) can be ascribed to a contrast between Ti, P, Li, Al, Mn, Mg, Fe, and Ca on the one hand, and Si and (Eu/Eu*) SN on the other. In other words, this component describes contrast of detrital components and some authigenic phases, such as carbonate and phosphate, against Si-richness and hydrothermal influence. A few metal enrichment factors (V EF , Re EF ) and δ 56 Fe display some covariance with Si. Principal component 2 (29% of variability) is strongly controlled by indicators of anoxia and reductive element drawdown, such as TOC, metal enrichment factors (Mo EF , U EF , and V EF ), and Fe accumulation (Fe, Fe/Al), with minor opposite contributions from Li, Si, and δ 56 Fe. The third principal component (13% of variability) includes major contributions from Al and Fe, and opposite ones from Mg, Mn, P, Re EF , V EF , P, (Eu/Eu*) SN , and Ca, possibly representing a contrast between clay-rich and carbonate lithologies. The fourth component (6% of variability) is strongly controlled by δ 56 Fe, together with Ca, Mn, Fe, and Fe/Al, which contrast modestly with Ti, Li, P, Al, and metal enrichment factors.

. Post-Depositional Alteration
Pyrite-the dominant Fe phase in much of the OnZaP section [15]-is typically considered to be resistant to metamorphic resetting of δ 56 Fe [63]. Yet, previous studies of Fe geochemistry in the Zaonega Formation have revealed post-depositional alteration of pyrite and Fe-rich carbonates [15,24], likely related to thermal gradients and hydrothermal circulation driven by the emplacement of magma bodies into partially unlithified Zaonega Formation sediments [29]. Therefore, in order to glean palaeoredox insights from the data reported here, the effect of alteration on the preserved proxies must first be assessed.
Asael et al. [24] studied the Fe geochemistry of a section from the FAR-DEEP 13A drill core, collected~300 m from the OnZaP 1 and 3 drill sites. The FAR-DEEP 13A core has been correlated with the OnZaP section based on lithology and C isotope stratigraphy [16], so that the top of the section studied by Asael et al. [24] (130-175 m depth in the 13A core), corresponds to~40 m below the bottom of the OnZaP section. In that study, Fe speciation results indicated a significant amount of pyrrhotite (Fe 0.8-1 S), a phase commonly associated with pyrite metamorphism [99]. The relative pyrite to pyrrhotite content correlated with the S/Fe ratio, suggesting that pyrrhotite formed through thermal breakdown of pyrite and the loss of S, whereas the Fe pool remained immobile and no shifts in the bulk δ 56 Fe were expected [24].
In the OnZaP cores, as studied by Paiste et al. [15], very little pyrrhotite was detected by XRD, suggesting only minor S mobility, although some acid-volatile sulphide, commonly associated with Fe monosulphides (e.g., pyrrhotite and sphalerite), was liberated during HCl-digestion. Still, other evidence for local Fe remobilisation is present in calcitised dolostone beds as large (>100 µm) inclusion-rich euhedral pyrite, likely of late diagenetic or metamorphic origin, which may have formed from the recrystallisation of previously disseminated pyrite or on account of the expulsion of Fe from dolostones during calcitisation. Pyrite in mudstone beds, on the other hand, is commonly smaller in size (typically~10 µm), more disseminated, and likely early diagenetic in nature (though larger, late-stage pyrite crystals have also been observed). The strongest evidence for fluid alteration in the OnZaP section was seen in the visibly altered silica-rich 77-70 m interval, and in the 53-44 m dolostone interval, which hosts thick silica veins [15].
Overall, Fe alteration in the Zaonega Formation likely occurred in a closed system and resulted in the redistribution of Fe between different mineral phases with little net Fe gain or loss. This means that little change in bulk-sample Fe/Al and δ 56 Fe are expected, especially in mudstone intervals. However, iron speciation is known to have been significantly affected by Fe redistribution in the carbonate beds [15], and potential S mobility [24].
We performed a series of checks to reduce the potential effects that alteration may have on our analyses and interpretations. First, the MSP0010 sample set (from which all δ 56 Fe values were measured) was chosen specifically from mudstone intervals, avoiding conspicuous mono-mineral clusters and silicate veins. Second, we removed data points from the silicified 77-70 m and 53-44 m intervals from further palaeoenvironmental consideration. Although no anomalous Fe/Al ratios are observed in these silicified intervals, they host uniformly low δ 56 Fe values (<0.2‰), which is consistent with the possible input of 56 Fedepleted dissolved Fe(II) during hydrothermal alteration ( Figure 4A). Third, we focussed the discussion below primarily on bulk rock proxies, rather than mineral phase-specific indices (i.e., Fe speciation) since the former are less susceptible to diagenetic/metamorphic repartitioning of elements-for the same reason, we measured δ 56 Fe in whole-rock samples, not in pyrite extracts. This conservative approach allows us to conclude that secondary alteration is unlikely to have resulted in the wholesale overprint of primary δ 56 Fe signals in the OnZaP section, although it may have contributed to slightly increased scatter in the relationships explored below. However, we do note that the potential effects of hydrocarbon migration on Fe phases and δ 56 Fe signals in the Zaonega Formation remain unconstrained as we are not aware of any study that has investigated such effects in natural or experimental settings.
Overall, Fe alteration in the Zaonega Formation likely occurred in a closed system and resulted in the redistribution of Fe between different mineral phases with little net Fe gain or loss. This means that little change in bulk-sample Fe/Al and δ 56 Fe are expected, especially in mudstone intervals. However, iron speciation is known to have been significantly affected by Fe redistribution in the carbonate beds [15], and potential S mobility [24].
We performed a series of checks to reduce the potential effects that alteration may have on our analyses and interpretations. First, the MSP0010 sample set (from which all δ 56 Fe values were measured) was chosen specifically from mudstone intervals, avoiding conspicuous mono-mineral clusters and silicate veins. Second, we removed data points from the silicified 77-70 m and 53-44 m intervals from further palaeoenvironmental consideration. Although no anomalous Fe/Al ratios are observed in these silicified intervals, they host uniformly low δ 56 Fe values (<0.2‰), which is consistent with the possible input of 56 Fe-depleted dissolved Fe(II) during hydrothermal alteration ( Figure 4A). Third, we focussed the discussion below primarily on bulk rock proxies, rather than mineral phasespecific indices (i.e., Fe speciation) since the former are less susceptible to diagenetic/metamorphic repartitioning of elements-for the same reason, we measured δ 56 Fe in wholerock samples, not in pyrite extracts. This conservative approach allows us to conclude that secondary alteration is unlikely to have resulted in the wholesale overprint of primary δ 56 Fe signals in the OnZaP section, although it may have contributed to slightly increased scatter in the relationships explored below. However, we do note that the potential effects of hydrocarbon migration on Fe phases and δ 56 Fe signals in the Zaonega Formation remain unconstrained as we are not aware of any study that has investigated such effects in natural or experimental settings. . Panel (A) shows the data against typical crustal Fe/Al ratios (vertical grey span) [35], and δ 56 Fe range (grey horizontal span) [100], as well as the hydrothermal δ 56 Fe range (yellow horizontal span) [101,102]. Purple cross symbols are data from silicified intervals [15]. Panel (B) presents an interpretation in which the trends in δ 56 Fe values are explained through a benthic Fe shuttle (similar to the Black Sea) [38], and an open-marine authigenic Fe shuttle component (similar to the Peru margin) [66].

Principal Influences on Iron Geochemistry
Principal component analysis ( Figure 3, Table 1) reveals that the geochemistry of the Zaonega Formation sedimentary rocks is most significantly controlled (PC1) by variability in Si-richness, which seems to coincide with (Eu/Eu*)SN, an indicator of hydrothermal influence (see Section 6.3). The second component (PC2) is dominated by TOC, and authigenic Fe/trace metal enrichments. This correlation potentially describes the extent of anoxia, or conditions that are conducive to the development of authigenic metal enrichments. The third component (PC3) may describe an aspect of carbonate richness, whereas PC4, accounting for only 6% of variance, is most affected by δ 56 Fe. . Panel (A) shows the data against typical crustal Fe/Al ratios (vertical grey span) [35], and δ 56 Fe range (grey horizontal span) [100], as well as the hydrothermal δ 56 Fe range (yellow horizontal span) [101,102]. Purple cross symbols are data from silicified intervals [15]. Panel (B) presents an interpretation in which the trends in δ 56 Fe values are explained through a benthic Fe shuttle (similar to the Black Sea) [38], and an open-marine authigenic Fe shuttle component (similar to the Peru margin) [66].

Principal Influences on Iron Geochemistry
Principal component analysis ( Figure 3, Table 1) reveals that the geochemistry of the Zaonega Formation sedimentary rocks is most significantly controlled (PC1) by variability in Si-richness, which seems to coincide with (Eu/Eu*) SN , an indicator of hydrothermal influence (see Section 6.3). The second component (PC2) is dominated by TOC, and authigenic Fe/trace metal enrichments. This correlation potentially describes the extent of anoxia, or conditions that are conducive to the development of authigenic metal enrichments. The third component (PC3) may describe an aspect of carbonate richness, whereas PC4, accounting for only 6% of variance, is most affected by δ 56 Fe.
Iron geochemistry contributes most significantly to PCs 2 and 4, suggesting that redox conditions during deposition played a role in determining Fe distribution patterns. Curiously, δ 56 Fe has only weak contributions to the first two PCs, and a strong contribution to PC4, meaning that variations in δ 56 Fe do not coincide strongly with variations in other geochemical parameters. The first two PCs place a weak loading on δ 56 Fe, with the most conspicuous opposite loading placed on Fe and Fe/Al. This implies a potential anticorrelation between these variables, whereas unidirectional covariation between the same components is seen on PC4. This may also suggest that δ 56 Fe data are affected by several different mechanisms resulting in opposing trends-PCA generally fails to deconvolute such overlapping trends [98]. Given these interpretations, we approach our subsequent discussion on Fe geochemistry from the perspective of hydrothermal influence, detrital composition, and redox conditions.

Detrital Versus Authigenic Iron Sources
We investigated the influence of detrital components on δ 56 Fe in Figure 5. There is significant variability in the Si/Al ratio that seems to exert some influence on δ 56 Fe-only relatively Al-rich samples host fractionated δ 56 Fe ( Figure 5A,B). Possibly, the Si-rich and unfractionated end member represents a sandy or silty detrital component, as opposed to a more clay-rich component with higher Al content. Coarser grain size would imply a shallower palaeoenvironmental setting with higher rates of deposition and a lower likelihood of developing anoxia, given that equilibrium would have been maintained with a weak to mildly oxidising atmosphere, that could limit non-quantitative Fe redox cycling and the generation of fractionated δ 56 Fe. Clay deposition, on the other hand, occurs at lower rates and is associated with higher rates of organic matter loading, which may be characteristic of deposition beneath a chemocline, and in turn facilitate the development of anoxia and redox driven Fe cycling. However, grain size variation in the OnZaP section is not very large and what little variation there is may largely be controlled by turbidite deposition [15]. Alternatively, if Si-richness describes hydrothermal influence, the association of Si-rich samples with non-fractionated δ 56 Fe implies that hydrothermal fluids did not host fractionated Fe.
Samples with high Li contents denote another clearly distinguishable detrital component, which associates with Al-richness. This relationship may develop as a result of Li being concentrated in clay minerals [103]. Figure 5C reveals that the most Li-rich sediments tend to host negative δ 56 Fe values, and likely represent the weathering product of an uncommonly 56 Fe-depleted parent rock.

Hydrothermal Iron
Iron isotope variability is often studied in relation to authigenic Fe enrichments, proxied by Fe/Al, as the highest fractionations are associated with redox or precipitation processes resulting in authigenic Fe accumulation in sediments [32]. A link between these two variables is also implied by the PCA results. When viewed on a regression plot ( Figure 4A), there appears to be a general trend towards lower δ 56 Fe values as Fe/Al rises, though the correlation is weak (R 2 = 0.15). However, there is significant scatter and several samples that display an opposite trend towards positive δ 56 Fe at higher Fe/Al. One possible interpretation for the main negative trend relates to Fe(II) delivered to the Onega basin from hydrothermal vents. High-temperature dissolution of Fe typically results in δ 56 Fe compositions lower than the igneous background and the δ 56 Fe of such Fe(II), therefore, trends towards negative values (from −0.67‰ to −0.09‰, more commonly between −0.5‰ and −0.1‰) [101,102]. Since most δ 56 Fe values in the Zaonega Formation lie within the reported ranges of δ 56 Fe for crustal and hydrothermal Fe [101,102], it is possible that variations in δ 56 Fe follow a mixing relationship between detrital Fe dominated by crustal isotope compositions and Fe(II) derived from basinal hydrothermal input.
A commonly used proxy for tracking the influence of hydrothermal fluids is the europium anomaly ((Eu/Eu*) SN ; Figure 5D). Unlike neighbouring rare earth elements, trivalent Eu can be reduced to a divalent state in acidic, high-temperature fluids, in which case it complexes with chlorides and remains highly soluble compared to other rare earth elements [104]. Hydrothermal vents, therefore, commonly display relative enrichments in Eu, and while the anomaly is quickly diluted in modern open ocean settings away from the venting sites, it is a useful proxy for hydrothermal influence in more restricted basins [104]. Strong Eu anomalies have previously been reported from authigenic apatite and carbonates in the Zaonega Formation [12,21], whereas authigenic Fe(II) sulphides and Fe(III) oxyhydroxides, if precipitating from hydrothermally influenced seawater, are known to preserve significant Eu anomalies [104,105]. Therefore, if δ 56 Fe variations are driven mainly by varying proportions of hydrothermal Fe(II) input, a tendency towards higher Eu anomalies with lower δ 56 Fe values may possibly be expressed. A commonly used proxy for tracking the influence of hydrothermal fluids is the europium anomaly ((Eu/Eu*)SN; Figure 5D). Unlike neighbouring rare earth elements, trivalent Eu can be reduced to a divalent state in acidic, high-temperature fluids, in which case it complexes with chlorides and remains highly soluble compared to other rare earth elements [104]. Hydrothermal vents, therefore, commonly display relative enrichments in Eu, and while the anomaly is quickly diluted in modern open ocean settings away from the venting sites, it is a useful proxy for hydrothermal influence in more restricted basins [104]. Strong Eu anomalies have previously been reported from authigenic apatite and carbonates in the Zaonega Formation [12,21], whereas authigenic Fe(II) sulphides and Fe(III) oxyhydroxides, if precipitating from hydrothermally influenced seawater, are known to preserve significant Eu anomalies [104,105]. Therefore, if δ 56 Fe variations are driven mainly by varying proportions of hydrothermal Fe(II) input, a tendency towards higher Eu anomalies with lower δ 56 Fe values may possibly be expressed.
However, when we plotted (Eu/Eu*)SN against δ 56 Fe in the Zaonega Formation sedi- However, when we plotted (Eu/Eu*) SN against δ 56 Fe in the Zaonega Formation sediments, we found no such relationship-rather, δ 56 Fe fractionations seem to be associated specifically with samples that display very little to no Eu anomaly ( Figure 5D). Still, this in itself is not conclusive evidence against a hydrothermal influence, since in bulk samples, detrital siliciclastic material is likely to host the majority of the Eu budget with only minor contributions from authigenic phases (with the exception of apatite that efficiently scavenges rare earth elements; samples hosting this mineral account for many of the positive Eu anomalies seen of Figure 5D). In the absence of filtering to deconvolute detrital Eu from authigenic Eu [12], it can be assumed that (Eu/Eu*) SN variations in our samples can also reflect high-temperature Eu disproportionation between plagioclase and other igneous components in the detrital sediment source [106].
Therefore, a hydrothermal influence on the δ 56 Fe budget remains a possibility. Yet, it cannot explain the subset of samples displaying positive δ 56 Fe (Figure 4), nor the apparent covariation with redox-sensitive metal concentrations (Figure 2). Hence, there is most probably an alternative influence on δ 56 Fe.

Iron Redox-Shuttling
A common mechanism responsible for covariation between Fe/Al and δ 56 Fe in modern anoxic basins is the benthic Fe shuttle [38,64,107]. Dissimilatory iron reduction occurs widely in shelf sediments and is a major source for Fe(II) into the ocean as much of the resulting Fe(II) can escape from pore water into the bottom water [57,108]. In modern redox-stratified and (semi-)restricted basins-most notably the Black Sea [38]-there is a net transfer of this reactive Fe from the shelf into the basin [107]. Deep currents transport DIR-derived Fe as colloids or fine particulates to the deeper, anoxic parts of the basin [109], where it is quantitatively captured as pyrite upon reaction with dissolved sulphide [107]. The end result is that sediments in the anoxic basin become enriched in authigenic Fe relative to the continental crust (elevated Fe/Al ratios), whereas shelf sediments become comparatively depleted (depleted Fe/Al) [64]. In parallel, since DIR preferentially mobilises 54 Fe, the shelf sediments will tend to become more positive in δ 56 Fe. Partial reoxidation of dissolved Fe(II) during basin-wards transport will preferentially remove even more of the 56 Fe, leading to significant negative authigenic δ 56 Fe inputs into basinal sediments [38]. Anticorrelation between Fe/Al and δ 56 Fe has also been detected in Phanerozoic [65] and Neoproterozoic [110] drill core data, reflecting facies change from oxidised shelves to anoxic deep basins, settings conducive to benthic Fe shuttling.
On Figure 4B, data from the Zaonega Formation are plotted next to Black Sea data from Severmann et al. [38]. Broadly, these datasets lay out a similar trend, suggesting that a benthic Fe shuttle may be responsible for much of the variation in δ 56 Fe. A major difference with the Black Sea dataset is that Zaonega Formation samples display much higher Fe/Al ratios. This may be due to higher Fe(II) concentrations in the deep Proterozoic oceans, given the likelihood that they remained pervasively anoxic-although we have previously argued that redox-sensitive metal accumulations and isotope ratios in the Zaonega Formation suggest that large parts of the oceans, in this sense surface waters, were oxygenated at the time these rocks were deposited [13]. The shift in Fe/Al may alternatively imply an unusually high Fe/Al ratio of background detritus-average Fe/Al ratios reach~1 in the upper part of the OnZaP section ( Figure 2). Finally, the shift may be due to a more vigorous benthic Fe shuttle in the Zaonega Formation compared to the Black Sea.
A subset of relatively Fe/Al-rich samples display an opposite trend and constitute the highest δ 56 Fe component (Figure 4). While the small number of samples (~5) makes any conclusion tentative, the coherence of the trend and the position of these samples in a constrained stratigraphic horizon between 54 and 58 m (Figure 2) hints at a palaeoenvironmental significance. Mansor and Fantle [63] suggested that positive sediment δ 56 Fe fractionations can reflect late-stage pyrite formed at slow rates of precipitation, possibly due to low sulphate or iron concentrations, in which case the positive equilibrium isotope fractionation factor inherent in pyrite precipitation overwhelms the negative kinetic isotope effect. However, sulphur isotope systematics of the OnZaP section have previously been used to argue that pyrite precipitation occurred rapidly due to high rates of organic carbon remineralisation in the sediments [15,16].
Authigenic accumulation of 56 Fe-enriched Fe is perhaps more likely the result of partial Fe(II) oxidation and precipitation as Fe(III) phases, which commonly host positive δ 56 Fe signatures [47][48][49][50]. In fact, a positive Fe/Al-δ 56 Fe trend has been noted in sediments deposited on the lower boundary of the anoxic wedge on the hydrographically open Peru margin [66]. In the Peru margin oxygen minimum zone, DIR produces a reductive Fe(II) flux, which leaks from sediments; as the dissolved Fe(II) migrates towards the suboxic margins of the oxygen minimum zone, a fraction of that Fe(II) will be oxidised and precipitated, producing sediments with Fe/Al enrichments and positive δ 56 Fe values. The data from Scholz et al. [66] that define this mechanism have been added for comparison purposes to Figure 4B. Compared to the Peru margin data, the positive sub-trend in the Zaonega Formation is shifted towards markedly higher Fe/Al ratios, likely for the same reasons as the main negative trend.
Ferric oxyhydroxide precipitation can, in principle, account for both positively and negatively fractionated δ 56 Fe values through progressive isotope distillation (Rayleigh fractionation) of the Fe(II) due to preferential removal of positively fractionated Fe [53]. However, distillation implies drawdown of the Fe pool and lower authigenic Fe contents in the negative δ 56 Fe end member, inconsistent with the Zaonega Formation data, where both positively and negatively fractionated samples are associated with higher Fe/Al ratios. Instead, mixing of Fe from authigenic Fe(III) oxyhydroxides and benthic shuttle-derived authigenic sulphides can better explain the data, including the scatter observed between δ 56 Fe and Fe/Al, despite most Fe likely having been subsequently converted to pyrite during diagenesis. A simultaneous influence of both of these two mechanisms, depending on the location of the sampling site relative to the oxygen-minimum zone and the deep basin, was noted by Scholz et al. [44] in the modern Guaymas basin in the Gulf of California. Namely, sediments in the oxygen-minimum zone were depleted in Fe/Al and slightly enriched in 56 Fe, suggesting Fe(II) release during DIR and loss into the water column. Slight accumulation of Fe(III) oxyhydroxides in sites below the oxygen-minimum zone was noted through elevated Fe/Al and δ 56 Fe. Sites in the deeper basin, where most of the released Fe(II) ended up accumulating, displayed higher Fe/Al together with less positive δ 56 Fe. We suggest that, similarly to the Guaymas basin, the Fe cycle in the Zaonega Formation was influenced by both partial Fe(III) precipitation and Fe(II) scavenging.

Redox Development in the Upper Zaonega Formation
The OnZaP section records a depositional setting that likely passed through facies shifts and experienced changes in basin chemistry and hydrographical conditions as the sediments accumulated. The bottom of the section (up to~80 m) records δ 56 Fe values close to 0‰, as well as low Fe/Al ratios, (~0.5; Figure 2). However, both TOC and trace metal enrichment is very high in that interval, implying a predominantly anoxic environment and arguing against high detrital deposition rates that could diminish authigenic Fe accumulation. Furthermore, XRD-based Fe speciation results imply that most Fe in that interval is hosted within pyrite, arguing against sulphide limitation [15]. In a basin hosting a benthic Fe shuttle, these types of sediments would be expected to capture Fe relatively efficiently. However, the lack of a benthic Fe shuttle signal in this interval suggests that the basin experienced diminished suboxic cycling and benthic Fe flux at that time. Possibly, this interval marks a period of more open hydrographic conditions, such that the benthic Fe flux ensuing from shelf areas was diluted over the open ocean. Scholz [43] identified basinal restriction as one of the most important parameters for determining if sediments capture a benthic Fe shuttle. The lack of a clear benthic Fe shuttle trend was also noted in a section of the FAR-DEEP 13A drill core by Asael et al. [24], stratigraphically below the OnZaP section. This suggests that the conditions preventing a benthic Fe flux from being recorded persisted through much of the middle Zaonega Formation succession.
The shift to clearly more variable and negative δ 56 Fe above 80 m coincides with an increase in metal enrichment factors and TOC (Figure 2), as well as pyrite abundance [15]. In short, the depositional setting became more conducive to Fe drawdown via reaction with sulphide. Whether this sulphide was present in the open water column or confined to sediments is a matter of debate-high Mo concentrations (as reported in Ref. [13]) are usually associated with sulphidic conditions [109], whereas S isotope systematics in the Zaonega Formation are more easily explained through the confinement of sulphidic conditions to the pore waters [15,16].
The interval between 60-53 m hosts the highest values for many redox indicators, which corresponds to efficient Fe drawdown, and explains why the interval has some of the highest Fe/Al ratios and the most negative δ 56 Fe values ( Figure 2). However, a small number of samples in this interval display markedly positive δ 56 Fe values and may imply partial Fe(III) oxyhydroxide drawdown, which can only occur under an oxidising water column. This discrepancy may potentially result from highly variable redox conditions, where a sharp sulphidic-oxic redoxcline fluctuated between the water column and the sediments. This is consistent with previous interpretations by Lepland et al. [72] including phosphate-rich sediments and putative fossilised sulphur oxidising microbial communities in a nearby correlative outcrop that were interpreted to have formed in sediments with a rapidly fluctuating sulphidic-suboxic redox boundary. Likewise, dramatically elevated trace metal enrichments previously identified in this interval are consistent with a fluctuating redox setting that is especially conducive to trace metal drawdown [13], where oxidative pulses supply the trace elements and sulphidic episodes sequester those trace metals efficiently into the sediment pile.
Above 35 m, δ 56 Fe values rise again to near 0‰, concomitant with a drop in TOC and metal enrichment factors. In this interval, putative evaporite mineral pseudomorphs together with rare crossbedding have been interpreted to reflect a shallowing depositional setting, whereas the decrease in relative pyrite and highly reactive Fe abundance have been interpreted to reflect waning anoxic conditions [15]. The overall baseline Fe/Al still remains elevated at~1, on account of increasing Fe concentrations in dolomite and Fe/Mgmicas (possibly a sign of changing sediment provenance) and also the appearance of Fe-bearing carbonate minerals ankerite and siderite [12,15]. The latter minerals may imply high levels of dissolved Fe(II), but limited sulphide content (i.e., ferruginous conditions) at the depositional site [111], which is consistent with the lack of a benthic Fe shuttle isotope signature in this section. Further, deviations in δ 56 Fe from the crustal mean persist (−0.18‰ to +0.18‰), suggesting Fe isotope fractionation during partial Fe(II) drawdown. Alternatively, Fe-rich carbonates in the upper OnZaP section have been interpreted as having formed diagenetically [12], in which case the δ 56 Fe values may reflect porewaterspecific processes.

The Onega Basin Palaeoenvironment
The initiation of benthic Fe shuttling in the Zaonega Formation (in the interval between 80 and 30 m) has implications for the basinal redox structure and hydrographic setting that can, in turn, affect the interpretation of trace element proxies and global redox conditions. A benthic Fe shuttle comes with two principal requirements-a strong redox gradient ensuring heterogeneity in Fe release and capture efficiencies in different parts of the basin, and a degree of basinal restriction. On the redox side, a benthic Fe shuttle first requires the presence of suboxic settings in shallower parts of the basin that allows for the accumulation of biologically available Fe(III) phases in the sediments that can then subsequently be converted through DIR to Fe(II) [38]. It also requires that Fe(II) diffusing out of the sediments is not quantitatively oxidised back to Fe(III) oxyhydroxides upon reaching the water column [64,65]. Secondly, euxinia is usually invoked to explain how the resulting Fe(II) flux ends up being quantitatively trapped in the deeper basin [38,112]. Both of these settings would need to have been simultaneously present in the Zaonega Formation to explain the overall δ 56 Fe-Fe/Al relationship.
Arguably more important than the redox condition is the shelf-to-basin ratio, which is closely tied to the degree of hydrographic restriction [43]. Specifically, the larger the areal ratio of the Fe source (the suboxic shelf) to the Fe sink (the euxinic basin), the more pronounced the signature of the benthic Fe flux will be, relative to lithogenic detrital Fe. For example, the modern Black Sea and the Guaymas basin, both of which host a benthic Fe shuttle, have very large shelf-to-basin ratios, a result of their restricted settings [44,64]. Conversely, the Peru Margin that opens to the Pacific Ocean does not display substantial accumulations of 56 Fe-depleted Fe [66]. The δ 56 Fe data presented here suggest basinal restriction during deposition of the Zaonega Formation, in agreement with previous interpretations based on S isotope systematics [15], and tectonic reconstructions in which the Onega Basin is viewed as a rift basin with the potential for a large shelf-to-basin ratio [19].
It is difficult to understand why Fe/Al ratios in the middle of the OnZaP section are so highly elevated ( 1). This may point either to high detrital Fe/Al ratios (unlikely, as Fe/Al below 80 m is generally~0.5), a change in sediment provenance, or an especially effective benthic Fe shuttle. The latter may result from the unusually high bioproductivity [13,16,29,72] and organic carbon loading rates in the Zaonega Formation (up to 70 wt.% TOC). On the sink side, biomass can contribute to high Fe/Al by directly scavenging Fe [113], or indirectly by encouraging the establishment of anoxia in sediments and the conversion of sulphate to sulphide, leading to more efficient Fe capture. On the source side, organic carbon loading can also drive increased rates of DIR in shallower settings, as has been noted in the modern Black Sea [114].
In contrast to the prevailing negative δ 56 Fe trends, the positive δ 56 Fe values in the 60-53 m interval imply the episodic presence of suboxia, which can partially draw down the dissolved Fe(II) pool resulting in positive fractionations via partial Fe(II) oxidation [66]. This implies a highly variable redox structure possibly driven by varying organic matter loading, or periodic influxes of oxic water into the deeper basin ( Figure 6). While the low number of highly positive δ 56 Fe samples make this conclusion tentative, strong support is provided by numerous other geological and geochemical inferences: The accumulation of hundreds of metres of massive halites and anhydrites in the Tulomozero Formation, which underlies the Zaonega Formation in the Onega Basin, requires both restriction in order to promote evaporation, and episodic replenishment in order to sustain salt accumulation for hundreds of metres [115]. Similar cycles of restriction and openness are implied in the Zaonega Formation by the Se and N isotope record [23] and by extremely elevated authigenic trace metal concentrations [13]. A modern analogue for such cycles is seen in Baltic Sea anoxic deeps, which see bi-decadal influxes of oxic water, which replenish the local trace metal pool and also lead to especially pronounced sedimentary trace metal enrichments [43,116].
Overall, the Fe record in the Zaonega Formation is consistent with other records of the middle Palaeopoterozoic-a time of elevated oxygen levels compared to the preceding Archean eon and the subsequent Mesoproterozoic [1,5]. Pyrite δ 56 Fe in organic-rich mudstones of this time period has been shown to be generally elevated (−0.5‰ to +1‰), compared to those from the Archean and Mesoproterozoic [53,117]. Rouxel et al. [53] interpreted this as being the result of high rates of sulphate input into seawater [116,118], which drove extensive pyrite precipitation, leading to the preferential removal of 54 Fe from seawater. Zaonega Formation δ 56 Fe data fit squarely within that range. Significantly, thẽ 1‰ variation in δ 56 Fe within the Zaonega Formation suggests a relatively small dissolved Fe reservoir in which Fe cycling proceeded largely quantitatively, contrasting with δ 56 Fe isotope data from both the Archean [117] and Neoproterozoic [110].

Conclusions and Implications
In this study, we measured Fe isotopes throughout a section of the Zaonega Formation, which were found to vary between −0.58‰ and +0.60‰. While some of the observed variation was driven by changes in the composition of detrital material, the bulk reflects two competing processes that led to authigenic Fe accumulation. Most of the negatively fractionated Fe, prevailing between 80-30 m in the OnZaP section, can be ascribed to a benthic Fe shuttle, akin to that seen in the modern Black Sea [38]. This Fe shuttle requires a redox gradient from oxidising shallow to anoxic and sulphidic deep waters, consistent with the inference that the deep Onega Basin at the time represented a predominantly anoxic basin in a generally oxygenated world. This finding agrees best with the results of Paiste et al. [15], who stressed the importance of basinal restriction on S isotope

Conclusions and Implications
In this study, we measured Fe isotopes throughout a section of the Zaonega Formation, which were found to vary between −0.58‰ and +0.60‰. While some of the observed variation was driven by changes in the composition of detrital material, the bulk reflects two competing processes that led to authigenic Fe accumulation. Most of the negatively fractionated Fe, prevailing between 80-30 m in the OnZaP section, can be ascribed to a benthic Fe shuttle, akin to that seen in the modern Black Sea [38]. This Fe shuttle requires a redox gradient from oxidising shallow to anoxic and sulphidic deep waters, consistent with the inference that the deep Onega Basin at the time represented a predominantly anoxic basin in a generally oxygenated world. This finding agrees best with the results of Paiste et al. [15], who stressed the importance of basinal restriction on S isotope characteristics-only in restricted settings will a distinctive Fe shuttle signature be recorded [43].
However, positive δ 56 Fe values in the 60-53 m interval attest to the episodic delivery of authigenic Fe(III) to the sediments, requiring periodic influxes of oxidised waters into the predominantly anoxic basin. Hence, restriction was not total, and intermittent connection with the open ocean persisted. These types of redox fluctuation-driven drawdownreplenishment cycles are especially conducive to the generation of high sedimentary enrichments of certain trace metals, explaining their extreme accumulations reported previously in the Zaonega Formation [13,23].
In totality, the Onega Basin can be seen as a complex sedimentary system affected to a large degree by local processes and conditions, such as high rates of organic matter loading and basinal restriction. Nevertheless, this basin remained in communication with the open ocean, and was thus able to incorporate global signals of element cycling. Therefore, the geochemistry of the Zaonega Formation can be used to deduce global environmental conditions in the Palaeoproterozoic but doing so requires nuanced assessment of the coupled effects of local and global redox and hydrographic conditions. Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/min11040368/s1, Table S1: Geochemical data and standard measurements, Table S2: Compilation of Fe geochemical data from modern anoxic basins.