A Comparison between Four Analytical Methods for the Measurement of Fe(II) at Nanomolar Concentrations in Coastal Seawater

Dissolved Fe(II) in seawater is deemed an important micronutrient for microbial organisms, but its analysis is challenging due to its transient nature. We conducted a series of Fe(II) method comparison experiments, where spikes of 5 to 31 nM Fe(II) were added to manipulated seawaters with varying dissolved oxygen (37 to 156 μM) concentrations. The observed Fe(II) concentrations from four analytical methods were compared: spectrophotometry with ferrozine, stripping voltammetry, and flow injection analysis using luminol (with, and without, a pre-concentration column). Direct comparisons between the different methods were undertaken from the derived apparent Fe(II) oxidation rate constant (k1). Whilst the two luminol based methods produced the most similar concentrations throughout the experiments, k1 was still subject to a 20–30% discrepancy between them. Contributing factors may have included uncertainty in the calibration curves, and different responses to interferences from Co(II) and humic/fulvic organic material. The difference in measured Fe(II) concentrations between the luminol and ferrozine methods, from 10 min–2 h after the Fe(II) spikes were added, was always relatively large in absolute terms (&gt;4 nM) and relative to the spike added (&gt;20% of the initial Fe(II) concentration). k1 derived from ferrozine observed Fe(II) concentrations was 3–80%, and 4–16%, of that derived from luminol observed Fe(II) with, and without, pre-concentration respectively. The poorest comparability of k1 was found after humic/fulvic material was added to raise dissolved organic carbon to 120 μM. A luminol method without pre-concentration then observed Fe(II) to fall below the detection limit (&lt;0.49 nM) within 10 min of a 17 nM Fe(II) spike addition, yet other methods still observed Fe(II) concentrations of 2.7 to 3.7 nM 30 min later. k1 also diverged accordingly with the ferrozine derived value 4% of that derived from luminol without pre-concentration. These apparent inconsistencies suggest that some inter-dataset differences in measured Fe(II) oxidation rates in natural waters may be attributable to differences in the analytical methods used rather than arising solely from substantial shifts in Fe(II) speciation.

section_content: Standardization of trace metal clean sampling techniques, inter-comparison exercises and the widespread use of consensus seawater reference samples have led to reproducible measurements of dissolved (<0. 20 µm) Fe concentrations in full depth profiles of all major ocean basins (Mawji et al., 2015), vastly improving our understanding of Fe biogeochemistry in the ocean. Such extensive exercises, however, have not yet been widely conducted with transient species such as Fe(II) which, due to their short-lived and highly reactive nature, present many challenges to the analyst. The rapid oxidation rate of Fe(II) in oxic seawater (King et al., 1995; Santana-Casiano et al., 2005; Sarthou et al., 2011), and the sensitivity of Fe(II) concentration and oxidation/reduction rates to multiple variables including pH, temperature, light intensity, dissolved oxygen (O 2 ), reactive oxygen species (ROS) and natural organic material (NOM; Davison and Seed, 1983; Millero et al., 1987; Rose and Waite, 2002) make sample collection logistically challenging and mean that the accuracy of Fe(II) measurements is difficult to verify. 

Fe redox chemistry is an important feature of the marine Fe cycle. In surface waters, photochemical processes form dissolved Fe(II) from Fe(III) either directly by photoreduction of Fe(III) species, or indirectly via ROS (O'Sullivan et al., 1991; Barbeau, 2006; Croot et al., 2008). Dissolved Fe(II), whilst short-lived in oxic seawater, accounts for a large fraction of dissolved Fe (DFe) entering the ocean from shelf sediments in oxygen minimum zones (OMZs; Lohan and Bruland, 2008; Vedamati et al., 2014), from hydrothermal vents (Statham et al., 2005; Sedwick et al., 2015), from rainwater (Zhuang et al., 1995; Willey et al., 2008) and possibly also from ice melt (Lin and Twining, 2012). From an energetic perspective, Fe(II) is expected to be more bioavailable than Fe(III) (Sunda et al., 2001). Furthermore, regardless of whether dissolved Fe(II) is actually more bioavailable than dissolved Fe(III) to specific cellular uptake systems, redox cycling in oxic waters is always expected to have the net effect of maintaining Fe in more labile phases (Croot et al., 2001; Emmenegger et al., 2001) due to the enhanced solubility of Fe(II) relative to Fe(III). 

The impracticality of creating stable transient-species reference materials, means that uncertainties remain about how comparable and how accurate existing Fe(II) data are. This is especially the case at nanomolar Fe(II) concentrations (which can occur, for example, in coastal waters, estuaries and OMZs) where multiple methods are available, each with different potential interferences and analytical constraints. In broad terms three different analytical approaches have been taken to measure Fe(II) concentrations in natural waters: luminol based chemiluminescence (Seitz and Hercules, 1972; O'Sullivan et al., 1995; Bowie et al., 2002), spectrophotometric methods using ferrozine (Stookey, 1970; Waterbury et al., 1997) or other synthetic ligands (Smith et al., 1952; Hennessy et al., 1984), and voltammetry (Gledhill and Van Den Berg, 1995). Luminol based flow injection techniques are by far the most commonly used to measure Fe(II) concentrations at sea, largely because of the sub-nanomolar detection limit that can be achieved (Bowie et al., 2002; Hansard and Landing, 2009), but other methods are still actively used in several different environments and applications. Ferrozine, for example, remains the most common analytical approach used in both autonomous Fe(II) sensors (Chin et al., 1994; Sarradin et al., 2005) and also in freshwater and estuarine environments (Kieber et al., 2001; Statham et al., 2012). 

Here, we compare four analytical methods for the determination of nanomolar Fe(II) concentrations in coastal seawater. We aim to assess the uncertainty associated with measuring and interpreting Fe(II) concentrations in natural waters via different analytical approaches. The four approaches selected for comparison were: ferrozine (Waterbury et al., 1997), Cathodic Stripping Voltammetry (CSV; Gledhill and Van Den Berg, 1995), and flow injection analysis (FIA) using luminol chemiluminescence both with (Bowie et al., 2002) and without (Croot and Laan, 2002) a pre-concentration step. In this study, in order to assess whether Fe(II) methods produce comparable Fe(II) concentrations and oxidation rates under environmentally relevant conditions, the experimental conditions for a series of method comparison experiments were selected to test the hypothesis that Fe(II) stability in coastal seawater is substantially affected by small changes in NOM composition and initial Fe(II) concentration. Varying spikes of Fe(II) were added to batches of Mediterranean seawater under pre-bloom conditions, post-bloom conditions (after a mesocosm experiment) and after addition of a humic/fulvic spike to represent high NOM estuarine waters.
section_name: INTRODUCTION

section_content: A series of eight method comparison experiments (Table 1 ) were conducted in Crete (June 2016). Coastal seawater was collected using acid cleaned (0. 3 M HCl soak for 3 days followed by 3 deionized water and 3 seawater rinses) tubing and high density polyethylene (HDPE) containers. Seawater was used either as collected or, in order to adjust the sample matrix with respect to NOM that may adversely affect Fe redox chemistry in coastal seawater (Santana-Casiano et al., 2000; Rose and Waite, 2003), under post-bloom conditions after a mesocosm experiment. The 12 day mesocosm was conducted as part of the Ocean Certain project in 1,500 L HDPE containers (cleaned as above) where macronutrients had been added daily [6 nM phosphate, 48 nM nitrate, 48 nM ammonium and silicic acid as necessary to maintain an excess according to the Redfield ratio (Redfield, 1934) ] to stimulate a phytoplankton bloom. In both cases, seawater was stored under reduced lighting conditions for >3 days prior to the Fe(II) method comparison experiments, in order to lower the concentration of ROS. 

In addition to the comparison of luminol (with, and without, pre-concentration), ferrozine and voltammetry methods for Fe(II) analysis in Mediterranean seawater, potential artifacts associated with spectrophotometric analysis of Fe(II) concentrations were investigated using three different commercially available synthetic Fe(II) ligands: ferrozine (as per the 4-method comparison in coastal seawater), ferene and 2,4,6tripyridyl-2-triazine (TPTZ). As spectrophotometric methods are most commonly used in freshwater, this comparison was undertaken using aged River Itchen (Southampton, UK) water which was filtered (0. 2 µm) into 1 L low density polyethylene (LDPE) bottles and then stored at 6 • C in the dark. The River Itchen was selected for its relatively high pH (typically ∼8) and modest DFe (100-500 nM) content (Statham et al., 2012; Hopwood et al., 2015).
section_name: MATERIALS AND METHODS

section_content: Eight method comparison experiments were conducted in Mediterranean seawater using a combination of the selected methods to measure Fe(II) concentrations in either 20 or 2 L containers over time periods of 0. 5-2 h after 5 to 31 nM spikes of dissolved Fe(II) were added (1-8, Table 1 ). Both 20 L and 2 L containers were pre-cleaned with mucasol detergent (Sigma-Aldrich) for 1 day, followed by 1 week in 1. 2 M HCl with 3 deionized water rinses after each stage. Experiments were carried out in a temperature controlled room (measured range 19. 9-21. 5 • C), where coastal seawater was stored prior to commencing experiments, to ensure a near constant seawater temperature. The two FIA instruments (A-no pre-concentration, B-with a preconcentration column) were placed adjacent to the experiment container. Analysis using voltammetry and spectroscopy was undertaken in a nearby temporary trace metal clean room where all surfaces were covered with plastic sheeting, the air was continuously filtered through polyester dust filters and clean suits were worn by the analysts. Stock solutions, reagents, manipulated seawater solutions and sample vials were prepared and handled in laminar flow hoods within the temporary trace metal clean room. An Fe(II) stock solution was prepared from ammonium Fe(II) sulfate hexahydrate (99. 997% Sigma Aldrich) in de-ionized water acidified with 1 µL HCl (trace metal grade, Fisher) per mL solution. A 1 µM Fe(II) stock solution was made daily and diluted to make standard additions and spikes. Prior to commencing experiments, N 2 (99. 999% purity) was gently bubbled through aged seawater for 3-12 h to lower the dissolved O 2 concentration and thus reduce the anticipated oxidation rate of added Fe(II) spikes (Santana-Casiano et al., 2005). 20 L collapsible low density polyethylene (LDPE) containers, filled with coastal seawater, were subject to external pressure to minimize the air headspace and gas exchange during the experiments. Experiments comparing only the two luminol based methods (numbers 6-8, Table 1 ) were undertaken in 2 L opaque HDPE containers with a constant, gentle stream of N 2 gas bubbles throughout the experiment. 

For all method comparison experiments, an Fe(II) spike was added at time zero and followed by immediate mixing (inversion of the containers by hand) for ∼30 s. For one experiment (number 4, Table 1 ) a spike of 10 mL humic/fulvic stock solution (unfiltered) was added before the seawater was bubbled with N 2. The two luminol-based Fe(II) FIA methods had polytetrafluoroethylene (PTFE) sampling lines running straight into the experimental containers, i. e., no sub-sampling of discrete aliquots was required from either 20 or 2 L containers. For the voltammetry and ferrozine methods, sub-samples were collected using a plastic syringe through a sample line also fitted into the 20 L experimental container. The sample line was first rinsed using the experimental seawater. 50 mL aliquots of seawater were then collected and immediately emptied into trace metal clean, opaque 125 mL HDPE sample bottles which were pre-laced with reagents. The time at which these sub-samples were collected was recorded (specifically the time at which the 50 mL sample was emptied into the pre-prepared vials). Full details of apparatus and analytical procedures for the four Fe(II) methods is included in Supplementary Material. 

As much standardization as possible was introduced across the four techniques for Fe(II) analysis in Mediterranean seawater. The same 1 µM Fe(II) stock solutions were used by all analysts, vials and reagent bottles were all pre-prepared according to the same cleaning procedure and, where chemicals were used by multiple methods, the same reagent supplies were used. For two reasons it was not considered necessary, or desirable, for analysts to be unaware of the spiked concentration of Fe(II). First, all methods had to be calibrated over the expected range of Fe(II) concentrations (something which could normally be estimated roughly based on ambient water conditions). Second, there was already considerable uncertainty concerning the starting dissolved O 2 concentration, and therefore also the Fe(II) concentration at the first measured time point (which was always >1 min after the Fe(II) spike was made at t = 0 due to the requirement for mixing). The dissolved O 2 concentration was only known after the completion of each method comparison experiment. 

For each spiked method comparison experiment (1-8, Table 1 ) the calculated Fe(II) oxidation rate constant (assumed to be first order with respect to [Fe(II)]) was derived from Millero et al. (1987) using measured salinity, temperature, dissolved O 2, pH (Table 2 ) and the corresponding values for K w and I. The expected Fe(II) oxidation rate was also estimated from Roy and Wells (2011) using the two surface (5 m depth) oxidation rate constants reported (at 25 • C and pH 8. 0) from Ocean Station Papa (North Pacific). These apparent rate constants were adjusted for dissolved O 2 as per equations 1 and 2 (originally assumed to be saturated, now presented for measured dissolved O 2 concentrations, Table 2 ); but not for the discrepancies in salinity, temperature or pH as these were small and the relationships with k not derived in the original dataset. 

-

Total dissolvable Fe (TdFe) samples (unfiltered) were collected at the last time point for each of the method comparison experiments 1-8 (Table 1 ) in trace metal clean 125 mL LDPE (Nalgene) bottles. TdFe samples were then acidified to pH <2. 0 by the addition of HCl (150 µL, UPA grade, Romil) and stored upright for 6 months prior to analysis. Samples were then diluted using 1 M distilled HNO 3 (Spa grade, Romil, distilled using a sub-boiling PFA distillation system, DST-1000, Savillex), and subsequently analyzed by high resolution inductively coupled plasma-mass spectrometry (HR-ICP-MS, ELEMENT II XR, ThermoFisherScientific) with calibration by standard addition. H 2 O 2 (unfiltered) samples were also collected at the last time point in opaque HDPE (Nalgene) bottles and analyzed by luminol chemiluminescence (Yuan and Shiller, 1999), ∼2 h after collection (stored in the dark at 20 • C) to allow time for Fe(II) to decay to background levels. Total organic carbon (TOC, unfiltered) samples were collected at the same time in precombusted 50 mL glass vials, preserved by the addition of 100 µL HCl (trace metal grade, Fisher) and subsequently analyzed by high-temperature combustion analysis (Farmer and Hansell, 2007) (Shimadzu TOC-L Total Organic Carbon Analyzer). 

For method comparison experiments 1-8 (Table 1 ), water samples for carbonate parameters and dissolved O 2 analysis were collected 5 min after commencing each experiment (approximately corresponding to the first voltammetry/ferrozine time point). Dissolved O 2 was measured using an Oxyminisensor (World Precision Instruments). Total alkalinity and dissolved inorganic carbon measurements were made on a VINDTA 3C system (Marianda, Kiel, Germany), then used to calculate pH (reported at the experimental temperature on the free pH scale) as per Tynan et al. (2016). Conversion of pH from carbonate parameters determined at 25 • C, to pH free values at the experimental temperatures was done using CO2SYS in MATLAB (van Heuven et al., 2011). 

A humic/fulvic organic stock solution was prepared by adding 50 mg Suwannee River humic acid standard I and 50 mg Suwannee River fulvic acid standard II (International Humic Substances Society, IHSS, used as received) to 20 mL de-ionized water (Milli-Q,18. 2 M •cm) and then allowed to stand for 2 days at 6 • C in the dark without filtration. Co(II) and V(IV) interferences, with luminol A and B, were tested by spiking Co(II) and V(IV) into the seawater matrix used for method comparison experiments 7 and 8 (Table 1 ) after Fe(II) concentration fell below the detection limit of both methods. Calculated spikes of 0. 10 nM Co(II), 1. 0 nM Co(II), 2. 0 nM Co(II), 0. 95 nM V(IV), 4. 7 nM V(IV) and 47 nM V(IV) were added sequentially.
section_name: Method Comparison Experiments in Spiked Seawater

section_content: A comparison of different Fe(II) ligands used for spectrophotometric analysis was made using River Itchen (Southampton, UK) water (experiment 9, Table 1 ). The ligands ferene (3-(2-pyridyl)-5,6-bis(2-(5-furyl sulfonic acid))-1,2,4-triazine (Hennessy et al., 1984) and TPTZ (Collins and Diehl, 1960) were used to determine Fe(II) concentrations spectrophotometrically with the LWCC apparatus as described for ferrozine (Supplementary Material), but with the first monitored wavelength for peak absorbance changed from 562 to 593 nm. 10 mM stock solutions of ferene were made in de-ionized water, whereas 10 mM stock solutions of TPTZ were made in de-ionized water acidified to pH 3 by addition of HCl (trace metal grade, Fisher). Prior to analysis, filtered (0. 2 µm) River Itchen water was stored in the dark at 6 • C for 1 week. 500 µL ammonium acetate buffer (pH 8) was then added to 10 mL aged river water aliquots followed by a Fe(II) ligand (ferene/TPTZ/ferrozine) to a final ligand concentration of 40 µM. After a mixing time of between 2 min and 20 h in the dark at room temperature, the sample was loaded into a LWCC and absorbance measured. Further experiments were conducted with varying concentrations (40-1,000 µM) of ferrozine.
section_name: Method Comparison of Ferrozine, Ferene, and TPTZ in River Water

section_name: RESULTS

section_content: Measurement of TOC (74 µM in coastal seawater, 80 µM in post-mesocosm water), and pH (8. 22 to 8. 24) in the experimental seawaters verified that the composition was relatively homogeneous across the different water batches (Table 2 ). TOC determined on Deep Sea Water reference material (Batch 8-2008, Hansell Laboratory) was 45 ± 0. 9 µM (consensus value 41-43 µM). After 3-7 days stored under reduced lighting conditions, H 2 O 2 was significantly lower (12-27 nM) than the original concentration measured when the seawater was collected (90 nM for coastal seawater and 70 nM for post-mesocosm water). Combined with the applied deoxygenation, the reduced H 2 O 2 concentrations should have slowed the oxidation rate of added Fe(II) spikes (Millero and Sotolongo, 1989) Experiment number (Exp., same as per Table 1 ), salinity (S), total organic carbon (TOC, unfiltered), total dissolvable Fe (TdFe, unfiltered), voltammetry measured reactive iron (Fe R ), detection limit (DL), not determined (ND). 

of Fe(II) oxidation with time, it should not have affected the comparison of methods. Analysis of the Certified Reference Materials NASS-7 and CASS-6 via ICP-MS yielded Fe concentrations of 6. 21 ± 0. 62 nM (NASS-7, certified 6. 29 ± 0. 47) and 26. 6 ± 0. 71 nM (CASS-6, certified 27. 9 ± 2. 1). TdFe (collected at the end of experiments 1-8), was generally within the range expected for coastal seawaters (Table 2 ). For the experiments where large Fe(II) spikes were added (up to 31 nM), the added Fe(II) spike was no longer present in solution (as DFe), or in suspension (as TdFe), at the end of the experiments (∼1 h after spike addition). This is consistent with the low solubility of Fe(III) in seawater (Millero et al., 1995; Millero, 1998) and the added Fe(II) spike adsorbing and/or precipitating on the container walls after oxidation. One exception was the humic/fulvic spiked experiment (Table 2 ), where 32 nM TdFe remained in suspension at the end of the experiment which most likely indicated the presence of more stable Fe(III)-organic complexes or colloids (Kuma et al., 1998) following Fe(II) oxidation. All the experiments used relatively high spikes of Fe(II) (5. 2 to 31 nM) which was necessary to keep the Fe(II) concentration above the detection limit (Table 3 ) of all methods for multiple datapoints. Whilst not representative of open ocean or oxic shelf waters, where reported Fe(II) concentrations are typically picomolar (Ussher et al., 2007; Hansard et al., 2009; Sarthou et al., 2011), our experimental conditions (Table 2 ) form a reasonable representation of conditions within tropical OMZs where Fe(II) concentrations on the order of 10-100 nM can be found in coastal waters (Hong and Kester, 1986; Lohan and Bruland, 2008; Vedamati et al., 2014). 

Detection limits for Fe(II) were defined as 3 standard deviations of observed concentrations in seawater stored in the dark at 21 • C (determined after the addition of all components to seawater other than the Fe(II) spike) for luminol A/voltammetry/ferrozine, and as 3 standard deviations of observed concentrations for an analytical cycle without sample flow for luminol B. In all cases detection limits were sufficiently low to measure the initial decline in Fe(II) concentration. Mean detection limits (Table 3 ) exclude experiment 4 where the addition of riverine NOM caused much higher detection limits for all methods except voltammetry.
section_name: Summary of Seawater Method Comparison Experiments

section_content: Luminol A

In all Fe(II) method comparison experiments in seawater (1-8, Table 1 ), the luminol A approach showed a decline in Fe(II) concentrations from 52 to 110% of the initial Fe(II) concentration added to below the detection limit of the method (Figure 1 ). The propagated uncertainty in the concentration of Fe(II) added to the seawater matrix at time zero ([Fe(II)] t0 ) was ± 9%. The initial Fe(II) measurement, which for luminol A occurred 72-247 s after the spike addition, was always expected to be less than [Fe(II)] t0. 

The most likely explanation for the slight overestimate of Fe(II) concentration (which was always <101% [Fe(II)] t0, except for experiment 8) was the modest uncertainty associated with the initial spiked concentration.
section_name: Comparison of Fe(II) Results from Spiked Seawater Experiments

section_content: The luminol B approach produced data with a lower temporal resolution due to the longer analytical cycle required for preconcentration (3. 5 min compared to 1. 0 min for luminol A) and 1 and 2. 

was only used in experiments 3-8 (Table 1 ). As a result of a poor calibration for experiment 5, data for this experiment was flagged (the suspected reason was the incomplete removal of organic material on the 8-HQ column from the previous experiment). For the five other method comparison experiments, the first measured Fe(II) concentration (152-242 s after t = 0) was 68 to 89% of [Fe(II)] t0. The trend in all experiments determined using the luminol B approach was generally similar to luminol A, with Fe(II) falling below the detection limit of the method in all experiments (Figure 1 ). Experiment 4, which included an addition of riverine humic/fulvic material, was a notable exception where a residual, and apparently stable, Fe(II) concentration of 4 nM was still measured over 30 min after Fe(II) concentration declined below the detection limit of luminol A (0. 49 nM) (Figure 1D ). Experiments 6, 7, and 8 were conducted on smaller water volumes comparing only the two luminol methods (luminol A and B) at intermediate O 2 concentrations, and the two methods generally produced similar Fe(II) concentrations (Figures 1F-H ). To test the effect of known inorganic interferences on both methods simultaneously, spikes of V(IV) and Co(II) were added to the seawater matrix of experiments 7 and 8 after Fe(II) concentration fell below the detection limit of both luminol A and B. The equivalent Fe(II) signal was then derived from the resulting chemiluminescence peaks (Figure 2 ). The apparent Fe(II) signal arising from V(IV) was similar for both methods (Figure 2 ). However, the Co(II) interference was curiously greatest with the luminol B method, despite the addition of DMG to the luminol B reagent mixture to suppress Co(II) interference (Ussher et al., 2009). With the luminol A approach, the propagated standard deviation in the apparent Fe(II) concentration arising from all Co(II) spikes was large compared to the "equivalent Fe(II)" concentration measured (0. 2 ± 0. 1, 0. 1 ± 0. 1, and 0. 1 ± 0. 1 nM, respectively).
section_name: Luminol B

section_content: During the method comparison experiments, within-sample replicates from the voltammetry method were quite variable (e. g., in experiment 5 up to 20%, Table 2 ). This was likely a result of the high apparent reactive Fe(III) concentration observed after the addition of the Fe(II) spike. Furthermore, it is difficult to explain the large increase in Fe R observed post Fe(II) spike, but these Fe R concentrations resulted in problems of carryover and reproducibility for the voltammetry method. Fe(II) concentrations estimated by voltammetry decreased after the first time point to below detection, with the exception of experiment 4 (Figure 1D ). The voltammetry method thus agreed with the results obtained from the luminol methods, except when riverine NOM was present. The depleted O 2 content of the experimental seawater was likely a contributing factor to the instability of Fe R, as the dissolved O 2 content of sampled water likely increased between sample collection and Fe R measurement, and was not necessarily constant between replicate sub-samples for Fe R.
section_name: Voltammetry

section_content: The ferrozine method produced data with a temporal resolution similar to that of the luminol B method. Approximately 5 min was required to complete analysis of a sample and flush the LWCC prior to loading the next sample. For the method comparison experiments (Table 1 ), initial Fe(II) concentrations of 31-151% [Fe(II)] t0 were obtained using the ferrozine method. The lowest recovery of the calculated spike added at the first sampled time point (31%) was for experiment 4 (Figure 1D ), where a fulvic/humic spike was added, and this recovery was indeed the lowest for any method in any experiment (the recovery range without this experiment was 72 to 153%). Observed Fe(II) concentrations using the ferrozine method never declined to below the detection limit (0. 3-0. 4 nM, Table 3 ) in any experiment. An apparently stable nanomolar Fe(II) concentration was found in experiments using coastal seawater (Figures 1A-C ), in experiments using post-mesocosm water (Figures 1D-H ) and also in experiments using aged river water (Figure 3 ).
section_name: Ferrozine

section_content: A comparison of different Fe(II) binding ligands used for spectrophotometric Fe(II) analysis suggested that all three ligand-Fe(II) methods reported a low nanomolar, stable Fe(II) concentration in aged river water (Figure 3 ). The observed Fe(II) concentration using the ferrozine method increased both with ferrozine concentration, and with increased time between ferrozine addition and absorbance measurement (Figure 3 ).
section_name: Comparison of Synthetic Fe(II) Ligands in Riverwater

section_name: DISCUSSION

section_content: The results presented here show that Fe(II) is an extremely transient species which is analytically challenging to quantify in seawater. During method comparison experiments in coastal seawater, where the concentration of Fe(II) declined due to rapid oxidation (Figure 1 ), it was difficult to directly compare individual datapoints for each method as they were subject to small offsets in measurement and/or sample timing. Even for the two FIA luminol methods, used here with automated timing, there was a non-negligible time difference comparing direct injection (luminol A) with pre-concentration (luminol B), which had implications for how the measured Fe(II) concentrations were interpreted. A qualitative comparison of datapoints from the different analytical techniques could be made visually (Figure 1 ). A more robust quantitative comparison, which should not be affected by a small measurement time offset, was made by deriving the apparent oxidation rate constant of Fe(II) oxidation (Equation 1 ). 

The apparent oxidation rate of nanomolar Fe(II) concentrations in seawater (Equations 1 and 2) is assumed to be first order with respect to the concentration of Fe(II) (Millero et al., 1987; Santana-Casiano et al., 2005), and thus the apparent rate constant k 1 (s -1 ) was derived from the gradient of a plot of ln[(Fe(II)] against time (e. g., Figure 4 ). The measured decline in Fe(II) concentration (Figure 1 ) following spikes of Fe(II) to seawater should approximate the oxidation rate, because oxidation is anticipated to be the dominant Fe(II) removal process under these experimental conditions (Table 2 ). 

The voltammetry determined Fe(II) concentrations were subject to a large uncertainty because of the high (5. 2 to 31 nM) Fe(II) spikes added. Voltammetry is not a technique ideally suited to analyses of high, or highly variable, DFe concentrations because it relies on subtracting measured Fe(III) from measured DFe [thus the propagated error for Fe(II) concentration is always larger than for a direct measurement of DFe or Fe(III)]. Also, voltammetry required an analysis time of 30 min per sample and was thereby slow compared to other available methods. Given the low temporal resolution of this data we therefore opted not to calculate k 1 for voltammetry data. 

The calculated oxidation rate of Fe(II) with O 2 in experiments 1-8 (Table 4 ) was derived using reported rate constants determined under laboratory conditions (Millero et al., 1987) and from measurements using a luminol chemiluminescence method in surface North Pacific seawater (Roy and Wells, 2011). As discussed further in prior work (Hansard et al., 2009; Sarthou et al., 2011), there are multiple reasons why these calculated rate constants are for indicative purposes only (e. g., O 2 may have changed over the duration of the experiments, the effect of DOC is not accounted for in this simple rate equation, the contribution of H 2 O 2 to the oxidation rate was likely nonnegligible, and the values derived from Roy and Wells (2011) are not corrected for small discrepancies in salinity or temperature). Yet, the differences between the luminol and ferrozine based methods are striking nonetheless (Figure 5 ). 

Luminol A produced rate constants comparable to those measured (after adjusting the dissolved O 2 concentration) in surface Pacific seawater with a luminol based method (Roy and Wells, 2011) and universally smaller than those derived from Millero et al. (1987), although we note that derived constants are very sensitive to how K w is parameterized. Luminol B was also generally comparable to the range derived from Roy and Wells (2011), except for experiments 4 and 5 (where the linear fit was notably poorer: R 2 0. 63 and 0. 62, respectively). For the two luminol methods (A and B), although the experimentally derived apparent rate constants were never identical (even considering the modest uncertainty in k 1 ), k 1 was similar in experiments 3, 6, 7, and 8 (Table 4 ). Whereas, the discrepancy in k 1 between luminol A and ferrozine measurements was always marked. The ferrozine derived Fe(II) time series exhibited a poor linear fit and thus the relative uncertainty was always high compared to luminol A and B. Curiously though, in the humic/fulvic spiked experiment (experiment 4) the pre-concentration method (luminol B) also indicated a residual stable Fe(II) concentration, TABLE 4 | A comparison of different methods for determining Fe(II) concentrations, which should not be affected by the measurement time offset between different methods, was made by deriving the apparent Fe(II) oxidation constant (k 1 ). SE, the standard error, (and R 2 ) is also shown for the linear regression. similar to the ferrozine technique. The ferrozine/luminol B derived values of k 1 were also similar for this experiment (Table 4 ). Between experiments 3 and 4 all conditions other than the humic/fulvic spike were similar (Table 2 ).
section_name: Derivation of Apparent Oxidation Rate Constants

section_name: Method

section_content: Lower than expected values of k 1, and correspondingly long halflives with respect to oxidation, have been reported in various natural waters using the ferrozine method (Kieber et al., 2001; Statham et al., 2005; Klar et al., in press ) with the foremost hypothesis to explain this phenomenon being the existence of stable Fe(II)-organic species. Whilst all methods here, consistent with data by Roy and Wells (2011), produced half-lives slower than expected based on relatively simple calculations from rate constants, ferrozine in all cases produced the longest half-lives (Figure 5 ). In no experiments were these half-lives similar to those derived from luminol A data. 

Various drawbacks have been identified with all available Fe(II) methods. It is relatively straightforward, for example, to demonstrate that the unmodified ferrozine method appears to over-estimate nanomolar Fe(II) concentrations in some natural waters (Figure 3 ; Murray and Gill, 1978; Box, 1984; Cowart et al., 1993), which could be attributed to reduction of labile Fe(III) phases in the presence of NOM after the ferrozine ligand is added. Short mixing times between ferrozine addition and absorbance measurement, combined with handling samples in the dark and buffering to near-neutral pH, is reported to minimize, although not eliminate, the interference due to labile Fe(III) reduction (Box, 1984; Pullin and Cabaniss, 2001). This interference may potentially be more important when a LWCC is used at low (nM) Fe(II) concentrations because of the more intense light source required as part of this apparatus. The contribution of light-exposure during analysis to a supposed overestimation of Fe(II) concentration cannot alone however explain the increase in measured Fe(II) concentrations in aged river water when the time interval between ferrozine addition and absorbance measurement was increased (Figure 3 ). (Roy and Wells, 2011), from method comparison experiments reported here, and calculated from rate constants (Millero et al., 1987). 

Conversely, the binding of Fe(II) to NOM could lead to an underestimation of Fe(II) concentration by the ferrozine method when a short mixing time is used due to slow complexation of organically associated Fe(II) (Kieber et al., 2005; Hopwood et al., 2014). This issue appears however to be confined exclusively to freshwater environments as there is presently no conclusive evidence that sufficiently strong Fe(II)-NOM interactions exist in seawater to impede complexation of Fe(II) by addition of a synthetic ligand (Croot and Hunter, 2000). 

Ferrozine derived Fe(II) results in the method comparison experiments (Figure 1 ) were therefore generally consistent with the known limitations of the unmodified ferrozine method in natural waters. The first datapoints were generally in agreement with the initial spike of Fe(II) added, confirming that the later over-estimation of Fe(II) could not consistently be attributed to an erroneous baseline. It remains possible that the residual Fe(II) concentration observed at the end of all experiments by ferrozine, and occasionally luminol B, (Figure 1 ) represents an Fe(II)species which is resistant to oxidation, and indeed the presence of organic compounds capable of retarding Fe(II) oxidation in seawater has been widely discussed in prior work (e. g., González et al., 2014; Santana-Casiano et al., 2014). However, it is not expected that such high residual Fe(II) concentrations (2. 7-9. 1 nM) would remain in aged Mediterranean seawater. Furthermore, such stable Fe(II) species should, in principle, be detected by all analytical approaches. 

Non-linearity in plots of ln[Fe(II)] against time, as observed for ferrozine Fe(II) data here (e. g., Figure 4 ), is often observed at low nanomolar Fe(II) concentrations, as under these conditions the assumption of a first order rate with respect to Fe(II) concentration does not appear to be universally correct (King et al., 1995; Santana-Casiano et al., 2006). However, here the non-linearity was much more pronounced using the ferrozine method. Yet it is not clear why such a large difference arises in both Fe(II) concentrations and k 1 between the ferrozine and luminol A method. Comparison of ferrozine with other synthetic ligands added to aged alkaline river water suggested that an increase in measured Fe(II) is to be expected with time whenever a strong Fe(II) ligand is added as part of a spectrophotometric method. It should however be noted that the rates of increase in measured Fe(II) concentration reported here (Figure 3 ) are likely not representative of all experimental conditions or all natural waters. The change in Fe(II) concentration after the addition of a synthetic Fe(II) ligand to natural water may be very sensitive to the type(s) of NOM present in solution, and thus experiments with river water and IHSS humic/fulvic additions are not directly comparable to experiments with additions of marine derived NOM. No attempt was made to optimize the methods for ferene and TPTZ to minimize the observed increase in Fe(II) concentration [the methods were used simply as previously described (Collins and Diehl, 1960; Hennessy et al., 1984) ] and thus nor should the comparison (Figure 3 ) be interpreted as demonstrating that the overestimation of Fe(II) concentration is consistently greater when using ferene or TPTZ instead of ferrozine. 

Results from Fe(II) method comparison experiment 4 (Figure 1D ) suggested that terrestrially derived NOM can impact the determination of Fe(II) concentrations in a number of different ways. Terrestrially derived NOM is known to contain Fe(III) chelators (Powell and Wilson-Finelli, 2003; Krachler et al., 2005), and suggested to contain Fe(II) chelators (Kieber et al., 2005; Statham et al., 2012), that result in the formation of relatively stable complexes and/or colloids. Dissolved organic Fe(II) species, with retarded Fe(II) oxidation rates relative to inorganic Fe(II), could therefore be formed by the addition of humic/fulvic compounds to seawater. Yet the addition of strong Fe(III) chelators into solution would be expected to facilitate Fe(II) oxidation (Boukhalfa and Crumbliss, 2002; Ussher et al., 2005) and indeed the acceleration of Fe(II) oxidation in seawater by NOM is used to explain faster than expected rates of decay observed in seawater. The formation of Fe(II)-organic species could possibly change the recovery of Fe(II) by the different methods employed, for example by altering the retention efficiency of Fe(II) on the 8-HQ column used in luminol B (Ussher et al., 2005) and this could produce bias in the derived oxidation rate. Also, as evidenced by a change in detection limit (Table 3 ), the chemiluminescence produced by luminol is generally suppressed by humic/fulvic material (Rose and Waite, 2001). 

In experiment 4, luminol A was the only method to determine that Fe(II) concentrations fell below the detection limit (0. 49 nM for luminol A) with an apparent rate constant greater than in any other experiment (the standard error on k 1 was also disproportionately high). This could, in isolation, be interpreted as an acceleration of the Fe(II) oxidation rate in the presence of NOM. However, we must also consider the consistent order of magnitude difference in k 1 derived from luminol A and ferrozine. In method comparison experiment 4, luminol B, voltammetry and ferrozine methods consistently observed a stable, >4 nM concentration of Fe(II) in solution for over 30 min after the addition of an Fe(II) spike. The only clear distinction between these three methods collectively and luminol A, is that luminol B, ferrozine and voltammetry all work via the complexation of Fe(II) to a ligand (8-HQ, ferrozine and Dp, respectively) whereas luminol A relies on only the catalytic effect of Fe(II) reacting with luminol. It remains a possibility that, on the timescale of the comparison experiments (0. 5-1 h), the development of Fe(II)-organic species contributes to a retarded oxidation rate. However, because the effect of organic material is method-dependent, it is also apparent that there is at least one method interference related to the presence of the added NOM. Whilst Fe(III)-organic speciation in seawater is at least moderately well-defined on a global scale (see, for example, Gledhill and Buck, 2012), Fe(II)-organic speciation in natural waters is almost entirely uncharacterized and, as we demonstrate here, is difficult to probe with existing methods which can exhibit conflicting results even in quite simple experiments. Thus, it is very difficult to evaluate if Fe(II)-organic species are present as a significant fraction of measured Fe(II) concentrations in natural waters and to assess the potential artifacts of any such Fe(II)-organic species on the accuracy of Fe(II) measurements. More generally, it is difficult to unambiguously confirm the classic interpretation for discrepancies in Fe(II) oxidation rates in natural waters; that faster than expected Fe(II) oxidation rates are driven by fractions of NOM (e. g., Roy and Wells, 2011; Lee et al., 2017), whereas slower than expected Fe(II) oxidation rates are driven by a fraction of Fe(II) present as stable Fe(II)-organic species (e. g., Willey et al., 2005; Klar et al., in press) ; when contradictory conclusions are drawn from two independent methods in the same series of experiments.
section_name: Potential Explanations for Inter-Method Differences in Fe(II) Concentration and Oxidation Rate

section_content: With FIA for the determination of total DFe in seawater (an analytical method for which it is possible to verify the accuracy of measurements) a mathematical approach shows that the determination of the correct gradient for a calibration curve is the single greatest error in the determination of DFe concentrations (Floor et al., 2015). For Fe(II) measurement via FIA using standard addition of Fe(II), this is therefore also very likely to be the case. Here each Fe(II) method was independently calibrated for each experiment and thus differences may have arisen because of random errors when comparing k 1 for luminol A and B. The relatively high detection limit for luminol B (Table 3 ) compared to luminol A partially arose because the pre-concentration system (luminol B) was optimized for low, or sub-, nanomolar concentrations. However, it is curious that whilst the interference for V(IV) was similar for both methods, the interference from Co(II) was more pronounced for the luminol B method (Figure 2 ). This was unexpected as the added DMG is known to suppress this interference (Ussher et al., 2009). It was suspected that the DMG decomposed during transport as several months previously the same batch was used successfully. 

A specific difficulty arises when calibrating the Fe(II) FIA methods at intermediate O 2 concentrations where Fe(II) is oxidizing as it is analyzed. Unless a standard addition of Fe(II) is made internally, exactly as the luminol reagent and seawater mix, the concentration of Fe(II) entering the FIA sample line should always be greater than the concentration of Fe(II) that enters the flow cell as luminescence is produced. Calibrating the method by standard addition into the same matrix as the samples via the same loading procedure thereby means that standard peaks decline as spiked solutions are run. For standard additions, the rate of decline of an Fe(II) spike should be identical to the rate of decline of Fe(II) in a sample-assuming that the sample matrix is non-changeable and ignoring possible complications due to hypothesized Fe(II)-organic species. Yet in reality the sample matrix cannot be maintained at its exact initial condition with respect to dissolved O 2 unless it is either practically anoxic (in which case it can theoretically be maintained in its ambient state with a flow of inert gas), or at equilibrium with atmospheric O 2. At intermediate O 2 concentrations, the dissolved O 2 content of a sample is likely to increase as it is analyzed increasing the Fe(II) oxidation rate. The alternative approach is running standard solutions with a stable concentration of Fe(II) (for example, as per the luminol B method, using a low pH standard matrix) (Hansard and Landing, 2009). In this case, the rate of decay of Fe(II) in the sample as it is analyzed is not accounted for. Yet, the calibration slope is likely subject to less uncertainty. It is not clear which is preferable from an analytical perspective.
section_name: Limitations to Analysis of Fe(II) Concentration in Seawater

section_content: Inter-study differences in how Fe(II) methods are calibrated appears to remain a minor issue compared to the apparent difference reported here between the ferrozine and luminol observed Fe(II) concentrations and k 1. It is yet to be established why, and under what circumstances, significant inter-method differences arise between absolute measurements of Fe(II) concentration and apparent oxidation rates. Whilst the accuracy of DFe measurements can be verified by analysis of a reference material, the transient nature of Fe(II) means that such a simple and robust verification of dissolved Fe(II) concentrations is not possible and thus greater uncertainty is generally acknowledged in Fe(II) concentrations. Yet changes in Fe(II) oxidation rate in natural waters are widely assumed to be a robust indicator of changes in Fe(II) speciation. Here, we suggest that Fe(II) oxidation rates in natural waters can be method-dependent and thus caution should be used when comparing apparent rate constants derived from different methods, especially across gradients in NOM and more generally between different sample matrixes. Process studies in natural waters using k 1 /half-life as an indicator of Fe(II) speciation should ideally demonstrate that the analytical methods used produce the expected decay rate of Fe(II) spikes to an aged sample matrix with environmentally relevant concentrations of Fe(II) added. As is already often the case (e. g., Sarthou et al., 2011), Fe(II) methods should also be frequently recalibrated using the most appropriate sample matrix in order to account, as best as possible, for the effects of any shifts in the sample matrix composition on measured Fe(II) concentration and oxidation rate.
section_name: Implications for Best Practice in Fe(II) Analysis

section_content: The almost universal adaptation of a luminol based method for measuring Fe(II) in seawater, largely because of its sensitivity and the possibility to achieve low picomolar detection limits, means that it should be logistically feasible to standardize a precise definition of "Fe(II)". Thus, even if Fe(II) was measured with unknown accuracy and not all interfering species or method artifacts were identified, Fe(II) data in specified environments would at least be internally consistent and precision could be improved. Such standardization has however yet to occur as at present there are methodological differences reported in sample handling, preservation and analysis techniques. Fortunately, whilst some uncertainties remain with respect to interfering species, it seems that different luminol FIA techniques produce similar Fe(II) concentrations and apparent oxidation rates in a simple seawater matrix. 

Yet here we found a consistent discrepancy between the results of the luminol and unmodified ferrozine methods both in observed Fe(II) concentrations and in apparent oxidation rate constants. This highlights the particular difficulty in measuring Fe(II) concentrations, and suggests that the uncertainty in the accuracy of Fe(II) oxidation rates with any presently available method is considerably larger than the standard error of the apparent oxidation rate constant. A safe conclusion is that deviations from expected Fe(II) oxidation rates in natural waters should be interpreted with caution as the artifacts of available Fe(II) methods are not fully understood.
section_name: CONCLUSIONS

section_content: ACKNOWLEDGMENTS D. Polyviou ( University of Southampton ) is thanked for H 2 O 2 measurements, E. Skjoldal ( University of Bergen )

for O 2 measurements, E. Reggiani ( NIVA ) for carbonate measurements and J. Gallego-Urrea ( University of Gothenburg ) for assistance with CO2SYS. M. Rijkenberg ( NIOZ ) is thanked for recommendations on calibration methods for Fe(II) FIA. Labview software for operating the luminol A FIA system was designed by P. Croot, M. Heller, C. Neill, and W. King. Linear regression (and standard errors) were calculated in SigmaPlot.

section_name: FUNDING

section_content: Analytical work in Crete was funded by the European Commission (OCEAN-CERTAIN, FP7-ENV-2013-6. 1-1 ;no: 603773 ). MH, MG, and EA gratefully acknowledge funding from the Collaborative Research Centre 754 ( SFB 754 ) " Climate-Biogeochemistry Interactions in the Tropical Ocean " funded by the German Research Foundation (DFG).
section_name: FUNDING

section_name: AUTHOR CONTRIBUTIONS

section_content: All authors contributed to the design of the study and interpreted the data. MH, AB, and MG conducted laboratory analysis. MH and AB wrote the initial draft of the paper and then all authors contributed to its revision.
section_name: AUTHOR CONTRIBUTIONS

section_content: The Supplementary Material for this article can be found online at: http://journal. frontiersin. org/article/10. 3389/fmars. 2017. 00192/full#supplementary-material
section_name: SUPPLEMENTARY MATERIAL

section_content: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. 

Copyright © 2017 Hopwood, Birchill, Gledhill, Achterberg, Klar and Milne. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
section_name: Conflict of Interest Statement: