# Exploring the Effects of Organic Matter Characteristics on Fe(II) Oxidation Kinetics in Coastal Seawater

J Magdalena Santana-Casiano, David González-Santana, Quentin Devresse, Helmke Hepach, Carolina Santana-González, Birgit Quack, Anja Engel, Melchor González-Dávila

► To cite this version:

J Magdalena Santana-Casiano, David González-Santana, Quentin Devresse, Helmke Hepach, Carolina Santana-González, et al.. Exploring the Effects of Organic Matter Characteristics on Fe(II) Oxidation Kinetics in Coastal Seawater. Environmental Science and Technology, 2022, 56, pp.2718 - 2728. 10.1021/acs.est.1c04512 . hal-05134803

HAL Id: hal-05134803

https://hal.univ-brest.fr/hal-05134803v1

Submitted on 29 Jun 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Exploring the Effects of Organic Matter Characteristics on Fe(II) Oxidation Kinetics in Coastal Seawater

J. Magdalena Santana-Casiano,* David González-Santana, Quentin Devresse, Helmke Hepach, Carolina Santana-González, Birgit Quack, Anja Engel, and Melchor González-Dávila

<figure><img src="image_2.png" /></figure>

Cite This: *Environ. Sci. Technol.* 2022, 56, 2718–2728

<figure><img src="image_3.png" /></figure>

Read Online

## ACCESS

<figure><img src="image_4.png" /></figure>

Metrics & More

<figure><img src="image_5.png" /></figure>

Article Recommendations

<figure><img src="image_6.png" /></figure>

Supporting Information

**ABSTRACT:** The iron(II) oxidation kinetic process was studied at 25 stations in coastal seawater of the Macaronesia region (9 around Cape Verde, 11 around the Canary Islands, and 5 around Madeira). In a physicochemical context, experiments were carried out to study the pseudo-first-order oxidation rate constant (k', min⁻¹) over a range of pH (7.8, 7.9, 8.0, and 8.1) and temperature (10, 15, 20, and 25 °C). Deviations from the calculated k'cal at the same T, pH, and S were observed for most of the stations. The measured t₁/₂ (ln 2/k', min) values at the 25 stations ranged from 1.82 to 3.47 min (mean 1.93 ± 0.76 min) and for all but two stations were lower than the calculated t₁/₂ of 3.21 ± 0.2 min. In a biogeochemical context, nutrients and variables associated with the organic matter spectral properties (CDOM and FDOM) were

analyzed to explain the observed deviations. The application of a multilinear regression model indicated that k' can be described (R = 0.921 and SEE = 0.064 for pH = 8 and T = 25 °C) from a linear combination of three organic variables, k'OM = k'cal -0.11* TDN + 29.9*bDOM + 33.4*C1humic where TDN is the total dissolved nitrogen, bDOM is the spectral peak obtained from colored dissolved organic matter (DOM) analysis when protein-like or tyrosine-like components are present, and C1humic is the component associated with humic-like compounds obtained from the parallel factor analysis of the fluorescent DOM. Results show that compounds with N in their structures mainly explain the observed k' increase for most of the samples, although other components could also play a relevant role. Experimentally, k' provides the net result between the compounds that accelerate the process and those that slow it down.

**KEYWORDS:** iron(II), oxidation kinetics, Macaronesia, coastal seawater, CDOM, FDOM

<figure><img src="image_7.png" /></figure>

## 1. INTRODUCTION

Iron (Fe) is a trace element whose speciation and concentration in the marine environment are affected by multiple abiotic and biotic processes.¹ Redox and complexation reactions control its solubility and therefore the fraction of dissolved and bioavailable Fe.² Processes such as photo-oxidation are a key factor in the control of Fe speciation in surface waters, particularly in areas subject to high irradiation and organic matter content.³ Iron is essential for organisms and plays an important role in the functioning of marine ecosystems.¹ Marine organisms have developed different strategies to assimilate Fe, ranging from the reduction of metal on the cell surface to the release of ligands that complex Fe.⁴ These ligands become part of the organic matter pool. The amount and type of organic matter present in the medium control the speciation of Fe.⁵ Dissolved organic matter (DOM) affects the Fe organic complexation, with the colloidal and particulate organic matter (COM and POM) acting as net dissolved Fe sinks.⁶,⁷

In the ocean, the thermodynamically stable form of Fe is the oxidized ferric iron (Fe(III)), with the soluble fraction in its majority complexed with organic ligands.⁸ However, ferrous iron (Fe(II)) can also be found in these conditions.⁹ In surface waters, many authors attribute it to photo-reduction processes in which complexed Fe(III), Fe(III)L, is involved.¹⁰ Abiotic redox reactions can also occur in the absence of light, causing the reduction of Fe(III) induced by oxidation of organic matter.¹¹ In anoxic areas and hydrothermal vents, Fe(II) is the predominant form and as soon as these waters are mixed with oxygen-rich waters, Fe(II) tends to oxidize.¹²,¹³ The persistence of Fe(II) coming from hydrothermal vents or anoxic sediments is attributed to the complexation with

Received: July 6, 2021

Revised: January 10, 2022

Accepted: January 11, 2022

Published: January 25, 2022

<figure><img src="image_8.png" /></figure>
organic ligands and the formation of iron sulfide (FeS).<sup>12,14</sup> In the water column, in situ processes such as the mineralization of POM during denitrification are possible sources of Fe(II).<sup>15</sup> Particle-bound Fe(III) can be reduced to Fe(II) by reduced sulfur compounds produced within reducing micro-environments of particles.<sup>16</sup> In the continental margin, the effect of allochthonous matter is important in the Fe(II) behavior controlling its oxidation.<sup>17,18</sup>

Knowing how long Fe(II) can remain in solution is one of the challenges in marine chemistry. This knowledge is required to understand the chemical behavior of this trace metal in the marine environment especially when considering that Fe(II) partakes and has played an essential role in the development of marine organisms since primordial times.<sup>19</sup> Information on the rate of oxidation of Fe(II) and the half-life time in the ocean is therefore essential for global biogeochemical models.<sup>20</sup>

This work aimed to study the Fe(II) oxidation kinetic constants in coastal seawater and to investigate the possible factors that could contribute to explain the different persistence of Fe(II) in the marine environment. These studies provide a rate constant (log k') and a half-life time (t₁/₂) data set under different pH and temperature conditions, which can be incorporated into ocean biogeochemical Fe models. The k' and t₁/₂ values at a fixed pH and T were also determined for every sample. This allows for the identification of non pH and T effects on the oxidation rate constants. Moreover, spectral parameters and ratios obtained from the colored and fluorescent DOM (CDOM and FDOM) were used to provide an equation that accounts for the deviation between measured k' and calculated k'cal. While the focus of this study is the identification of which variables associated with the organic matter spectral properties (CDOM and FDOM) can explain the deviations of the measured k' from the calculated k'cal at the same T and pH, additional attention should be given in future studies to the effects of other parameters such as Fe(II) and organic ligand concentration ratios and the effect of photo-generated compounds.

## 2. MATERIALS AND METHODS

**2.1. Sampling Protocol.** Samples were collected onboard the RV POSEIDON during the POS533 cruise. The cruise started on the 28th of February 2019 in Mindelo, Sao Vicente (Cape Verde), and ended on the island of Gran Canaria (Canary Islands) on the 22nd of March 2019. The cruise was divided into three sets of stations (Figure 1). The first set was located in the Cape Verde archipelago, the second set was positioned around the Canary archipelago, and the third set encompassed the Madeira area. The number and sampling locations are indicated in Table 1.

At each station, trace metal samples were collected at 20 m depth using a Teflon pump (Furon) with a 40 m Teflon tube connected with an AcroPak 1500 capsule w/Supor Memb 0.8/ 0.2 µm filter. The pump was left for 15 min in flowing seawater to rinse the inner hose, and the 1 L sample was collected in LDPE bottles (Nalgene) and stored at -20 °C until the land-based laboratory analysis. The material was previously cleaned following the trace metal GEOTRACES protocol,<sup>21</sup> and the experiments were performed in ISO Class 5 laminar flow hoods within an ISO Class 6 clean lab (QUIMA-IOCAG TM lab).

A lowered SeaBird SBE 9-plus CTD system equipped with two sets of pumped sensors measuring conductivity, temperature, and oxygen at 24 Hz was used. The CTD was mounted

<figure><img src="image_9.png" /><figcaption>Figure 1. Map of the archipelagos. Cape Verde Island, Canary Islands, and Madeira region. Station 4S* (Selvages) belongs to Madeira.</figcaption></figure>

<table><caption>Table 1. Station Number, Name, Location, and Bottom Depth of the sites where Fe(II) oxidation kinetics samples were collected during the POS533 Cruise<sup>a</sup></caption><thead><tr><th>Station number</th><th>Station name</th><th>Latitude (°N)</th><th>Longitude (°E)</th><th>Bottom Depth (m)</th><th>Station for OM</th></tr></thead><tbody><tr><td>1</td><td>Santo Antao</td><td>16.911</td><td>-25.166</td><td>100</td><td></td></tr><tr><td>4</td><td>Sao Nicolau</td><td>16.556</td><td>-24.386</td><td>100</td><td></td></tr><tr><td>6</td><td>Fogo</td><td>14.843</td><td>-24.263</td><td>500</td><td></td></tr><tr><td>9</td><td>Santiago SE</td><td>14.891</td><td>-23.507</td><td>100</td><td>9</td></tr><tr><td>11</td><td>Maio</td><td>15.093</td><td>-23.189</td><td>100</td><td>12</td></tr><tr><td>13</td><td>Boa Vista W</td><td>16.141</td><td>-22.973</td><td>100</td><td>13</td></tr><tr><td>15</td><td>Sal NW</td><td>16.680</td><td>-23.028</td><td>100</td><td></td></tr><tr><td>16</td><td>CVAO</td><td>16.892</td><td>-24.852</td><td>100</td><td></td></tr><tr><td>17</td><td>CVOO</td><td>17.583</td><td>-24.282</td><td>3200</td><td>18</td></tr><tr><td>21</td><td>Hierro SE</td><td>27.644</td><td>-17.958</td><td>200</td><td>22</td></tr><tr><td>23</td><td>H-G</td><td>27.959</td><td>-17.653</td><td>3000</td><td></td></tr><tr><td>25</td><td>Gomera SW</td><td>28.038</td><td>-17.316</td><td>100</td><td>25</td></tr><tr><td>31</td><td>Tenerife SW</td><td>28.028</td><td>-16.728</td><td>100</td><td>31</td></tr><tr><td>32</td><td>Tenerife S</td><td>28.009</td><td>-16.603</td><td>100</td><td></td></tr><tr><td>34</td><td>Tenerife E</td><td>28.317</td><td>-16.339</td><td>100</td><td>33</td></tr><tr><td>37</td><td>Tf-GC</td><td>28.320</td><td>-15.978</td><td>3000</td><td>37</td></tr><tr><td>39</td><td>Canaria W</td><td>27.811</td><td>-15.832</td><td>100</td><td>39</td></tr><tr><td>40</td><td>Canaria S</td><td>27.687</td><td>-15.704</td><td>100</td><td></td></tr><tr><td>42</td><td>Canaria NE</td><td>28.188</td><td>-15.391</td><td>100</td><td>41</td></tr><tr><td>44</td><td>ESTOC</td><td>29.167</td><td>-15.500</td><td>3500</td><td></td></tr><tr><td>45</td><td>Selvagens</td><td>30.018</td><td>-16.079</td><td>100</td><td>45</td></tr><tr><td>52</td><td>Desertas</td><td>32.404</td><td>-16.527</td><td>600</td><td>53</td></tr><tr><td>54</td><td>Ridge M</td><td>32.647</td><td>-16.590</td><td>600</td><td>55</td></tr><tr><td>61</td><td>Madeira c2</td><td>32.719</td><td>-17.238</td><td>100</td><td></td></tr><tr><td>63</td><td>Madeira c3</td><td>32.651</td><td>-17.114</td><td>100</td><td>63</td></tr></tbody></table>

on a SeaBird rosette with 12 bottles (10 L) which were used to take discrete water samples for sensor calibration and biogeochemical parameters.<sup>22</sup> The biogeochemical parameters sampled were nitrate + nitrite (NO₃ + NO₂), phosphate (PO₄), silicates, dissolved and particulate organic carbon

<footnote><sup>a</sup>Stations for OM are the closest stations where organic matter variables were collected.</footnote>
(DOC, POC), total dissolved and particulate nitrogen (TDN, PN), CDOM, and FDOM. The mixed layer depth was deeper than 20 m.²² Information about sampling and analysis is included in the Supporting Information.

**2.2. UV Irradiation.** An Ace Glass-incorporated ultra-violet photo-oxidation unit which operates at 120 V, 60 Hz was used to UV irradiate the samples. The quartz UV lamp works at 430 mA, 254 nm, and 3.5 W. For the irradiation procedure, each sample was placed in a 100 mL quartz tube and 12 tubes were irradiated at the same time. The samples were irradiated for 4 h. They were then kept in the dark and let to equilibrate with atmospheric air in a laminar flow for more than 72 h before the analysis.

**2.3. Oxidation Kinetic Experiments.** The oxidation kinetic experiments were carried out in a 60 mL LDPE container within a double-wall glass thermo-regulated cell, connected to a thermostatic bath (Julabo), as described in previous studies.²³,²⁴ For each oxidation kinetic study, 25 mL of the seawater sample was used. The initial concentration of added Fe(II) was 0.97 nM. Previous studies have shown that the inorganic Fe(II) oxidation rate constant is independent of the added Fe(II) concentration.²⁵,²⁶ The experiments were carried out at a fixed free-scale pH<sub>F</sub> = 8 with temperatures of 10, 15, 20, and 25 °C to get the temperature dependence and similarly, at a fixed temperature of 25 °C with pH<sub>F</sub> of 7.8, 7.9, 8.0, and 8.1 to get the pH dependence. In each oxidation kinetic experiment, the pH was controlled using a Titrino 719 (Metrohm) which automatically added 0.01 M hydrochloric acid (HCl, Panreac Hiperpur-plus). The pH electrodes were calibrated using a Tris buffer solution.²⁷ All the experiments were carried out under O₂ saturated conditions to avoid oxygen limitation. Experiments were performed in the dark to avoid the generation of reactive oxygen species (ROS) from photo-oxidation processes and the photo-transformation of Fe.

The concentration of total dissolved Fe(II), TdFe(II), in the samples for the kinetics experiments, was determined using the flow injection analysis by chemiluminescence (FIA-CL) technique with a FeLume system (Waterville Analytical)²⁸ as described in previous studies.¹²,²³ The FIA-CL technique uses luminol as the reagent. 5 L of the luminol reagent was prepared using 2.71 × 10⁻⁴ M of 5-amino-2,3-dihydro-1,4-phthalazine-dione (Sigma), 4.93 × 10⁻² M of Na₂CO₃ (Sigma-Aldrich), and 0.4 M of previously distilled 25% NH₃ (Panreac). The final pH was adjusted to 10.4 by adding 0.06 M Q-HCl. The luminol solution was stored in the dark due to its light sensitivity.

The kinetic studies were carried out in a thermo-regulated cell connected to a thermostatic bath (Julabo) with control to ±0.01 °C. For each study, the seawater sample was acclimated to the desired temperature (in situ and 25 °C). When the temperature was stable, the pH<sub>F</sub> for the sample was measured.²⁷ For each oxidation kinetic analysis, 50 mL of seawater was used. The seawater was placed in the thermostated cell, and the magnetic stirrer was switched on for 1 h to attain oxygen concentration equilibrium. When the solution stabilized at the desired temperature and pH, the sample line was introduced into the reaction cell. The samples were automatically mixed with the buffer just before being introduced into the detector. The sample plus luminol (1 mL min⁻¹ flux) and the ammonium acetate buffer (0.04 M ammonium acetate adjusted to pH 5.5 with acetic acid at 0.125 mL min⁻¹ flux) were introduced into the detector with a peristaltic pump. This modification provided continuous

registration of the measure. This was followed by the Fe(II) addition, and the time was registered.

Iron(II) oxidation kinetic studies in different aqueous media and conditions²⁵,²⁸⁻³⁷ have shown that the rate of oxidation with O₂ can be expressed as an apparent rate constant (k<sub>app</sub>), regardless of the mechanism that describes the process (eq 1).

$$ \frac{d[\text{Fe(II)}]}{dt} = -k_{\text{app}}[\text{Fe(II)}][\text{O}_2] \quad (1) $$

where k<sub>app</sub> = k[OH⁻]². The brackets indicate the total molar concentration.

The inorganic (i) and organic complexation (L) of Fe influences the kinetics rate constant,²⁹,³⁴ and the apparent rate constant includes the contribution of the inorganic and organic species of Fe(II)

$$ k_{\text{app}} = \sum_{i} k_{i} \alpha_{i} + \sum_{L} k_{L} \alpha_{L} \quad (2) $$

When the reaction is studied in excess O₂, the reaction can be considered pseudo-first-order (eq 3).

$$ \frac{d[\text{Fe(II)}]}{dt} = -k'[\text{Fe(II)}] \quad (3) $$

where k' = k<sub>app</sub>[O₂] in s⁻¹ and k<sub>app</sub> is expressed in M⁻¹ s⁻¹.

The slope corresponding to the plot of Ln [Fe(II)] versus time determines the k' (min⁻¹). In the Supporting Information (Figure S1), the linear relationship is shown for St1.¹²,²³,²⁸

The calculated k<sub>cal</sub> value was obtained from an empirical equation³⁵ for known conditions of temperature, pH, and salinity (eq 4). This equation is valid from 0.5²⁶ to 200 nM.³⁵ In this study, the experimental conditions were pH = 8, T = 25 °C, and S = 35.

$$ \log(k'_{\text{cal}}) = 35.407 - 6.7109 \times \text{pH} + 0.5342 \times \text{pH}^2 - \frac{5362.6}{T} - 0.04406 \times \sqrt{S} - 0.002847 \times S \quad (4) $$

The extended equation¹² calculated for a higher range of low temperatures applicable to deep waters (until 2 °C) can also be used.

**2.4. CDOM and FDOM Study.** For the CDOM study, the a₂₅₄ and a₃₂₅ absorption coefficients at the wavelengths 254 and 325 nm, respectively,³⁸ were selected together with the slope ratio (S<sub>R</sub>) and the E2/E3 ratios.³⁸,³⁹ The dimensionless S<sub>R</sub> was calculated from the ratio of the slope of the shorter wavelength region (275–295 nm) to that of the longer wavelength region (350–400 nm). The E2/E3 was calculated as the ratio of absorption at 250 to 365 nm.³⁹⁻⁴¹ From the FDOM study, the fluorescence peaks,⁴² the humification index (HIX), and the biological index (BIX) were determined⁴³,⁴⁴ and the parallel factor analysis (PARAFAC) was applied. The variable b<sub>DOM</sub> was obtained from the fluorescence peak analysis. b<sub>DOM</sub> is a peak that represents tyrosine-like and protein-like components (with fluorescence Ex/Em = 275/305), and it has been interpreted as a material of autochthonous origin.⁴⁵

The corrected excitation emission matrices (EEMs), the HIX⁴³ and BIX,⁴⁴ were used to determine the relative degree of humification and autotrophic productivity of fluorescent CDOM, respectively. HIX is the ratio of the fluorescence over Em 434–480 nm to that over 300–346 nm (at Ex 255 nm). BIX is the ratio of the fluorescence at Em 380 nm to that at 430 nm (at Ex 310 nm). HIX is the estimate of the degree of
DOM maturation with an increase in red-shifted emission which presumably arises from increasing conjugation (and possibly aromatization) of DOM.⁴³ BIX is an index for DOM sources. A BIX < 0.7 represents important terrestrial contributions, while a BIX > 1 is characteristic of DOM with a biological/aquatic bacterial origin.⁴⁴

The PARAFAC analysis for this study contained five components. Components C1_humic (Ex/Em: 235(365)/484) and C3_humic (Ex/Em: 230/398) had broad excitation and emission spectra, with excitation and emission maxima in the ultraviolet and visible region, respectively. They are traditionally referred to as humic-like components ⁴⁰. Component C2_autoDOM (Ex/Em: 240(300)/342) had a fluorescence signal which is similar to free and protein-bound amino acids and has been ascribed to autochthonous DOM ⁴⁵. Components C4 (Ex/Em: 230/306) and C5 (Ex/Em: 230(275)/340) were similar to tyrosine-like and tryptophan-like compounds.⁴⁵ All optical analyses including PARAFAC were conducted using the software R (v4.0.2.) in Rstudio⁴⁶ with the package staRdom (“spectroscopic analysis of DOM in R”).⁴⁷ All corrected EEMs of seawater in the POS533 cruise (n = 299) were used for modeling.

The identified fluorescent components were compared with the published data on an open-access spectral database⁴⁸ (OpenFluor, https://openfluor.lablicate.com). The database compares fluorescence data sets and determines Tucker congruence (set at 95%) as the similarity criteria between pairs of excitation and emission spectra. It provides a rapid way to test new PARAFAC models against those in the literature.

**2.5. Statistical Analysis.** The Pearson product-moment correlation and Spearman rank-order correlation statistical analysis tests between the oxidation rate constant ($k'$ at pH = 8 and T = 25 °C) and the measured biochemical variables were performed. A multiple linear regression (MLR) model was also applied to predict if the variability observed in the oxidation rate constants at a fixed temperature, pH, and oxygen saturated condition could be described by two or more organic spectral properties.

## 3. RESULTS AND DISCUSSION

The Fe(II) oxidation kinetics process was studied in coastal seawater considering both a physicochemical and a biogeochemical framework. In the physicochemical context, experiments were carried out to study $k'$ over a range of pH and T. The goal was to verify if there were differences in the calculated $k'_{cal}$ value from eq 4.¹²,³⁵ In the biogeochemical context, variables associated with the presence of OM were analyzed. The goal was to explain the deviations observed for $k'$ in the physicochemical context. These deviations are a consequence of the interaction between Fe and OM, expressed as ligands, L

$$ \mathrm{Fe}(\mathrm{II}) + \mathrm{L} \rightleftharpoons \mathrm{Fe}(\mathrm{II})\mathrm{L} \quad (5) $$

$$ \mathrm{Fe}(\mathrm{II})\mathrm{L} + \mathrm{O}_2 + \mathrm{H}^+ \rightarrow \mathrm{Fe}(\mathrm{III})\mathrm{L} + \mathrm{H}_2\mathrm{O} \quad (6) $$

$$ \mathrm{Fe}(\mathrm{III}) + \mathrm{L} \rightleftharpoons \mathrm{Fe}(\mathrm{III})\mathrm{L} \quad (7) $$

$$ \mathrm{Fe}(\mathrm{III})\mathrm{L} \rightleftharpoons \mathrm{Fe}(\mathrm{II}) + \mathrm{L}' \quad (8) $$

**3.1. pH Effect.** The dependence of log $k'$ on the pH obtained for each station was plotted (Supporting Information Figure S2) together with the expected calculated log $k'_{cal}$ under different pH conditions. In all the experiments, linear

dependences of log $k'$ with pH were observed. Deviations from log $k'_{cal}$ as a function of pH were observed for some stations.

The slope of the dependence of log $k'$ with the pH is shown in Figure 2. In Cape Verde, only St4 (Sao Nicolau, 1.68 ±

<figure><img src="image_10.png" /><figcaption>**Figure 2.** Representation of the slope of log $k'-$pH for the stations in Cape Verde, the Canary Islands, and the Madeira region. The slope was obtained from the plot of log $k'$ (min⁻¹) vs pH (7.8, 7.9, 8, and 8.1) (Supporting Information Figure S2). The slope for $k'_{cal}$ is included as St70.</figcaption></figure>

0.03) presented a value similar to the theoretical one (1.78 ± 0.03). The other stations presented lower slopes ranging between 0.60 ± 0.03 (St1, Santo Antao) and 1.50 ± 0.05 (St6, CVAO).

In the Canary Islands, St25 (Gomera), St32 (Tenerife-S), St39 (Gran Canaria-W), and St40 (Gran Canaria-S) presented a log $k'$-pH slope value close to the theoretical value. The other stations presented lower slopes ranging between 0.54 ± 0.03 (St37, between Tenerife and Gran Canaria) and 1.39 ± 0.08 (St34, Tenerife-E). In the north of the Canary Islands, the ESTOC site presented a slope of 0.97 ± 0.06.

Five stations were studied in the Madeira region. Station 45 (Selvagens) presented the lowest log $k'$-pH slope of 0.53 ± 0.03. At St54 (Ridge), the value was 1.56 ± 0.09, and St63 presented a slope (1.85 ± 0.03) similar to the theoretical one, and St61 presented the highest slope of the cruise (2.39 ± 0.06).

Two factors must be considered when interpreting the results. First, pH can change the speciation of both Fe(II) and organic ligands. Second, pH can affect the Fe(II)L complexation. The oxidation of Fe(II) in the absence of ligands is affected by pH as a result of speciation changes as demonstrated in previous research.³⁴ In seawater, degraded and recently produced organic compounds are present and can also interact with Fe(II). The log $k'$-pH slope shows the net effect of organic ligands on $k'$. If there is an interaction between the ligands and Fe(II), forming complexes with different strengths, there will be a displacement of the log $k'$-pH slope with respect to the calculated value. The slope will shift up or down if the net result produces an acceleration or a slow-down of the oxidation kinetics rate. Moreover, if the pH changes the speciation of these complexes, a change in slope will occur. This change is due to the net effect of the organic matter matrix. However, the scientific community is currently
not able to define the specific functional groups involved. In this sense, previous studies¹⁸,⁴⁹ found differences between the effects of allochthonous and autochthonous DOM on the Fe(II) oxidation rate constant.

**3.2. Temperature Effect.** The k' dependence with temperature was studied, and the activation energy (E_a) was calculated using the Arrhenius equation (eq 9).

$$ \ln k' = \ln A - E_a / RT \quad (9) $$

The ln k' versus 1/T for ESTOC, CVOO, and Selvagens is shown in the Supporting Information (Figure S3). The E_a for each studied station was plotted in Figure 3 and compared to

<figure><img src="image_11.png" /><figcaption>**Figure 3.** Activation energy for each station. The calculated value is also included as St70.</figcaption></figure>

the calculated E_a, cal (103 ± 3 kJ mol⁻¹) obtained from eq 4. For Cape Verde, the E_a ranged from 20.7 kJ mol⁻¹ in St1 (Santo Antao) to 106.4 kJ mol⁻¹ in St13 (Boa Vista W), the latter having an E_a within the theoretically calculated E_a. In the Canary Islands, St21 (El Hierro-SE) and St44 (ESTOC) presented an E_a (124.7 and 124.0 kJ mol⁻¹, respectively) higher than the theoretical E_a. The other stations presented low E_a values ranging from 20.8 kJ mol⁻¹ in St31 (Tenerife-W) to 97.5 kJ mol⁻¹ in St34 (Tenerife-E). In the Madeira area, St54 (Ridge) presented a high E_a value (117.7 kJ mol⁻¹), with coastal stations presenting 79.8 and 96.2 kJ mol⁻¹ (St63 and St61, respectively).

If the rate-controlling process is always the same, then the E_a obtained from the Arrhenius equation should not change for the same experimental conditions. Different E_a may involve different mechanisms or at least, different Fe(II)-organic matter species involved in the oxidation process.

For artificial seawater, without organic ligands, or seawater in which k' is not affected by the organic matter present,³⁵ the average slope is -5362 ± 162 and the activation energy is 103 ± 3 kJ mol⁻¹. The same results (-5434 ± 183) 104 ± 3 kJ were obtained in non-hydrothermally affected stations within the Mid-North Atlantic Ridge.¹² According to these studies, changes in the E_a are probably caused by the interaction of organic compounds with the Fe(II) species which affected the limiting Fe(II) oxidation step. Therefore, it may have a different oxidation reaction mechanism.

**3.3. The t₁/₂ at a Fixed pH and Temperature.** Previous studies showed log k' differences between stations when experiments were carried out under the same pH and T

conditions. From the pseudo-first-order rate constant (k'), the t₁/₂ was calculated as t₁/₂ = ln 2/k'. This variable represents the persistence time of Fe(II) in each station under the studied conditions. The t₁/₂ was calculated for pH = 8 and T = 25 °C conditions (Figure 4). The mean t₁/₂ for the three archipelagos

<figure><img src="image_12.png" /><figcaption>**Figure 4.** The t₁/₂ (min) for the stations in Cape Verde, the Canary Islands, and the Madeira region. The t₁/₂ theoretical value is included as St70.</figcaption></figure>

was 1.93 ± 0.76 min for Cape Verde, 1.82 ± 0.45 min for the Canary Islands, and 2.86 ± 0.66 min for Madeira, and the theoretical t₁/₂ was 3.21 ± 0.2 min. Overall, all stations presented values below the theoretical t₁/₂ except for St11 (3.09 min) in Cape Verde and St61 (3.10 min) in Madeira which presented values similar to the theoretical t₁/₂. Only two stations presented t₁/₂ values higher than the theoretical: St46 (3.47 min) and St52 (3.41 min).

Although factors such as O₂, pH, and temperature are the main variables that control the oxidation kinetics of Fe(II), photo-generated compounds (organic radicals and ROS)⁵⁰,⁵¹ and other biochemical variables can play a relevant role. When these other variables affect the oxidation process, either by accelerating or delaying it, the value obtained for k' is different from the calculated k'cal. Photochemical processes should be important in samples from the top few centimeters of the seawater column, while thermal processes control Fe transformation in deeper waters.⁵⁰

Variations in nutrient or DOC concentrations, the nature of the OM, and colloids present in the solution can lead to changes in the Fe(II) oxidation rate.¹²,²³,⁵²,⁵³ These factors were analyzed for their relationship with the observed changes in t₁/₂ in different stations at fixed pH and T.

**3.4. Nutrients.** The stations in the Cape Verde archipelago presented the highest and most variable nutrient concentrations, ranging from 0.03 to 1.16 µM for total inorganic nitrogen (nitrate and nitrite), 0.07 to 0.18 µM for phosphates, and 0.04 to 0.63 µM for silicates. The Canary Islands presented the lowest total inorganic nitrogen (0.08 ± 0.11 µM) while Madeira was characterized by the lowest phosphate (0.04 ± 0.01 µM) and silicate (0.03 ± 0.04 µM) concentrations. Concentrations for each station are shown in Supporting Information (Figure S4). The effects of nutrient concentrations, in particular N and Si, in the t₁/₂ are generally observed when nutrient concentrations are high, for example, in phytoplankton growing media or eutrophic environ-
<figure><img src="image_13.png" /><figcaption>Figure 5. Time evolution of k' for non-UV-irradiated seawater and UV-irradiated seawater for an ESTOC sample with 78 µM of DOC [(Fe(II))₀ = 0.97 nM, pH = 8, T = 20 °C, S = 36]. The continuous line defines the behavior for a sample that contains a DOC of 78 µM as in the ESTOC. The dashed red lines represent the calculated variation that would be obtained if the DOC content changes. The longest dashed line in green corresponds to a sample that has 91 µM DOC, as in the samples from Cape Verde. The dashed line in red corresponds to a sample that has 73 µM DOC, as in the samples from Madeira.</figcaption></figure>

ments.⁵²,⁵³ The concentrations found in these areas were low and did not exert an appreciable effect on k'.

**3.5. Total Dissolved Organic C and N.** The DOC and TDN concentrations varied by archipelago and oceanic stations (Supporting Information Figure S5). The mean DOC and TDN concentrations were, respectively, 91.3 ± 3.5 and 6.98 ± 0.7 µM in Cape Verde, 78.2 ± 0.4 and 5.2 ± 0.1 µM in the Canary Islands, and 73.8 ± 0.2 and 5.1 ± 0.1 µM in Madeira. A slight gradient was observed which increased from south to north between the archipelagos.

**3.6. Particulate Organic C and N.** Unlike DOC and TDN, the particulate organic C (POC) and N (PN) showed significant variations between stations of the same region (Supporting Information Figure S5). The south-to-north gradient was also observed for POC and PN with the lowest concentrations measured in Madeira, followed by the Canary Islands. The highest POC and PN concentrations were measured in Cape Verde. The differences observed between archipelagos in the particulate material suggested that variations in the content and distribution of the colloidal material could also occur.⁵⁴ However, the samples for Fe(II) were filtered through a 0.2 µm pore size filter, where the particulate material should not affect the obtained t₁/₂.¹² The colloidal material which has been demonstrated to influence the oxidation process¹² was not analyzed due to time constraints but should be taken into account in future work.

**3.7. Bulk of DOC and the UV-Irradiation Effect.** A 20 m sea-surface water sample of ESTOC was used to carry out studies with a known initial DOC content of 78 µM. The sample was divided into two sub-samples: one was irradiated, while the other was kept in the dark (non-irradiated). The evolution of k' over time was studied for 100 days at a fixed pH = 8 and T = 20 °C (Figure 5). For the non-irradiated seawater, k' was similar during all the experiments with a value of 0.092 ± 0.003 min⁻¹ and the t₁/₂(non-IR) = 7.5 ± 0.2 min. For the irradiated samples, high values were obtained the first few days

after the irradiation, presenting an exponential decay described by eq 10. k' changed from 3.59 ± 0.10 min⁻¹ on day 0 to 0.92 ± 0.02 min⁻¹ on day 30 at a rate of 0.088 ± 0.01 min⁻¹, after which it reached a plateau. The reduction of the initial k' by 50% was reached on day 7.

$$k' = 0.72(\pm 0.09) + 2.87(\pm 0.21) e^{-0.088(\pm 0.01)t} \quad (10)$$

std dev = ±0.14 min⁻¹.

The study carried out with a UV-irradiated and a non-irradiated sample presented differences between the two conditions. When the organic matter undergoes photo-oxidation, ROS intermediates, ROS (H₂O₂, O₂⁻, and OH·) are generated and these affect the Fe(II) oxidation kinetics rate constant.⁵⁵,⁵⁶ Organic radicals such as semiquinone radicals and/or peroxyl radicals can also be generated, becoming even more important than ROS under certain conditions.⁵⁰ Consequently, UV-irradiated samples should be left in the dark for at least 30 days to reach the plateau and measure a consistent k'.

The difference of k'*(day 30) - k'*(day 0) = 2.674 min⁻¹ is a measure of the photo-generated compound effect due to the irradiation process. The effect of OM present in the sample can be calculated from the difference between the non-irradiated and irradiated k'. The k' value changed by 0.82 min⁻¹, which accounted for an increase in the t₁/₂ of 6.78 min at pH = 8 and T = 20 °C conditions. The effect of three DOC concentrations (73, 78, and 91 µM for Madeira, ESTOC, and Cape Verde, respectively) is represented in Figure 5 assuming a change only in the amount but not in general composition.

The competitive effect between O₂ and H₂O₂ was studied in previous work,³⁶,⁵⁷ concluding that the oxidation of Fe(II) with H₂O₂ plays a relatively minor role in most natural waters. At the pH of seawater, O₂ is the most important oxidant when [H₂O₂] is below 200 nM and [Fe(II)] is at nanomolar
concentrations. Previous studies in the area showed that H₂O₂ concentrations are below 100 nM.⁵⁸

The nature of OM should also be considered. Culture studies carried out in laboratories show a high dependence of k' on the content of DOC.⁵⁹,⁶⁰ However, the same does not occur with oceanic samples¹²,²³ presumably because the dilution factor in the ocean may be important.

The properties of DOM are diverse and depend on its source (terrestrial or aquatic) and diagenetic state. The colorimetric properties of DOM give information about its origin and diagenetic status.³⁸,⁴⁴ The analysis carried out considered the fraction of OM that absorbs ultraviolet and visible light, through the CDOM, and considered the fluorescent properties of the FDOM. The CDOM and FDOM sampling did not always coincide with the kinetic study sampling. As a result, the analysis only compares common stations or geographically close stations (Table 1).

### 3.8. CDOM and FDOM.
Using spectral parameters and ratios obtained from CDOM and FDOM, different information about the characteristics of DOM was obtained (Supporting Information Figures S6–S11). The absorbances a₂₅₄ and a₃₂₅ are respectively proportional to the abundance of conjugated carbon double bonds⁶¹ and aromatic substances.³⁸,⁴⁰ The comparison of a₂₅₄ and a₃₂₅ measured for different samples did not show significant variations (Supporting Information Figure S6). For the entire study region, the a₂₅₄ average value was 1.07 ± 0.23 m⁻¹, with minimum and maximum values of 0.62–1.25 m⁻¹. The a₂₅₄ absorbances were in agreement with those obtained for temperate Atlantic waters⁶¹ where the higher a₂₅₄ absorbances were found in surface waters (1.43 ± 0.18 m⁻¹) decreasing to 0.87 ± 0.78 m⁻¹ in deep waters. In this study, the a₃₂₅ average was 0.11 ± 0.02 m⁻¹ and ranged between 0.06 and 0.15 m⁻¹. The determined coastal zone absorbances were slightly higher than those observed in oceanic waters of the North Atlantic subtropical gyre.³⁸ This may be due to the increased production in coastal waters.

The Sₘ and E2/E3 ratios are independent of the CDOM concentration and provide information on the average characteristics (chemistry, source, and diagenesis) of CDOM.³⁹ The average value for Sₘ (Supporting Information Figure S7) was 2.21 ± 1.29, ranging from 1.18 to 6.61, and was within the marine origin DOM Sₘ range. In open ocean waters, Sₘ varies from 1.5 to 4, while it reaches values lower than one when it has a terrestrial origin DOM.³⁹,⁶²

The E2/E3 ratio is used to track changes in the relative size of DOM molecules (Supporting Information Figure S8). As the molecular size increases, E2/E3 decreases due to stronger light absorption by high-molecular-weight CDOM at longer wavelengths.⁴¹ The E2/E3 ratio changed from 16 to 40 with an average of 24.39 ± 6.73. The highest ratios were obtained at St41, located close to St42 (Gran Canaria-E), while the lowest ratio was calculated at the oceanic St18.

The diagenetic state of DOM can be deduced from the BIX and HIX indexes.⁴⁴ The BIX index varied between 0.8 and 8.1 with an average of 1.80 ± 1.81 for the studied regions. In the Canary Islands, St31 (Tenerife-W) had the highest BIX index (8.1 ± 0.75). Stations 9 (Santiago) in Cape Verde and 21 (El Hierro-SE) in the Canary Islands presented values of 0.85 ± 0.39 and 0.96 ± 0.1, respectively (Supporting Information Figure S9). However, most of the stations measured presented BIX indices greater than 1. Consequently, DOM predominantly had an autochthonous origin. Increases in the BIX

index indicated a recently reworked one by bacteria DOM.⁴⁴ The HIX index ranged between 0.09 and 0.62 with an average of 0.48 ± 0.12. In this study, HIX indexes were always below 4, indicating autochthonous DOM from a biological origin.⁴⁴

The PARAFAC analysis characterized five components (Supporting Information Figures S10 and S11).⁶³ C1 varied between 0.008 and 0.015 RU. C1 was significantly higher in the Cape Verde region (ANOVA, Tukey, *p* < 0.05). C3 ranged from 0.004 to 0.025 RU. The Cape Verde region was significantly higher than the Canary Islands.⁴⁵,⁶⁴ C2autoDOM ranged from 0.004 to 0.018 RU. Two exceptions were found in the Canary Island region, at St31 and St33 (Tenerife W and E), where concentrations reached 0.54 and 0.06 RU, respectively. C4 was not significantly different between the regions (ANOVA *p* < 0.05) with an average of 0.018 ± 0.005 RU. C5 fluctuated from 0.012 to 0.12 RU with higher values at St31 (Tenerife-W) and St63 (Madeira-W) of 0.12 and 0.10 RU, respectively.

### 3.9. Statistical Analysis.
The measured biochemical variables (nutrients, DOC, TDN, CDOM-a₂₅₄, a₃₂₅, E2/E3, Sₘ, b-t-a-m-c peaks, FDOM-BIX, HIX, and C1-C5 from PARAFAC) at 13 stations were tested to identify any correlation with the oxidation rate constant k' at pH = 8 and T = 25 °C. In the previous work¹²,³⁵ the main controlling factors in the determination of k'cal were pH, temperature, salinity, and oxygen. The statistic results indicated a lack of significant correlations (*p* < 0.05) between k' and the biogeochemical variables. From the MLR model, the organic variable TDN and the spectral organic variables bDOM and C1humic were able to predict k' (eq 11) with an R-value of 0.921 and a standard error of estimate for k' of 0.064 min⁻¹.

$$k'_{\text{OM}} (\text{min}^{-1}) = k'_{\text{cal}} - 0.11 * \text{TDN} + 29.89 * b_{\text{DOM}} + 33.43 * \text{C1}_{\text{humic}} \quad (11)$$

where TDN is the total dissolved nitrogen (µM) and bDOM is the absorbance peak that appears when protein-like or tyrosine-like components are present,⁴² and C1humic is associated with humic-like components⁶³ in RU (Raman units). Equation 11 was able to explain 84 ± 10% of k' for the Macaronesia region (data available in Supporting Information Table S1). At fixed pH = 8, T = 25 °C, and k' = 0.218 min⁻¹, the statistical *p*-values were *p* < 0.001 for TDN, *p* = 0.015 for bDOM and *p* = 0.009 for C1humic.

The equation indicated that within the highly variable DOM, compounds containing nitrogen in their structure or nitrogen functional groups could exert an important effect on k'. The overall effect was that increasing concentrations of tyrosine-, protein- or humic-like compounds resulted in a higher observed k' than the k'cal. As a result, Fe(II) would have a lower t₁/₂ which would affect the permanence of Fe(II) in the ocean. This does not mean that only nitrogen-containing compounds have an impact on the Fe(II) oxidation process. Other functional groups have also been described in the literature.³⁴ It may be that in the study area, they are the most active.

In the presence of organic ligands L with N in the structure, eq 11 can be considered and explained by eqs 5 and 6

$$\text{Fe}(\text{II}) + \text{L} \rightleftharpoons \text{Fe}(\text{II})\text{L}$$

$$\text{Fe}(\text{II})\text{L} + \text{O}_2 + \text{H}^+ \rightarrow \text{Fe}(\text{III})\text{L} + \text{H}_2\text{O}$$
When the Fe(III)L (ferric chelate) stability constant is higher than that of Fe(II)L (ferrous chelate), oxidation of Fe(II)L to Fe(III)L by O₂ is a highly favored reaction.⁶⁵,⁶⁶ Moreover, for ligands that have a weaker nitrogen donor, the stability of Fe(II)L increases with basicity, but not as much as that of Fe(III)L.⁶⁵ Iron(II) adsorbed onto mineral surfaces and soluble Fe(II) chelates are important natural reductants. Studies with both Fe-goethite and Fe(II)-tiron as models of Fe(II) adsorbed and soluble Fe(II) chelates with different N–O containing compounds indicated that both the amino-functional groups and the pyridine ring are involved in complexation. Ring-N is more strongly involved than ring-O.⁶⁶

Specific components of the DON pool in the ocean include urea, dissolved combined and free amino acids, proteins, nucleic acids, amino sugars, and humic substances.⁶⁷ The chemical identity of most of these compounds and the mechanisms by which they are cycled are unknown.⁶⁸ Furthermore, DON originates from both allochthonous and autochthonous sources.⁶⁷ In this study, DON had a higher autochthonous origin (deduced from S<sub>R</sub>). Previous studies indicated that a substantial fraction of DON in the ocean has bacterial origin such as glutamine and glutamate, with the largest proportion released by planktonic marine cyanobacteria.⁶⁹,⁷⁰

Other compounds may also interact with Fe, such as
bacterial siderophores and planktonic exudates such as
polysaccharides and transparent exopolymers.⁵ In the same
way that microorganisms follow a seasonal cycle in surface
waters, the composition and concentration of autochthonous
ligands will be conditioned by seasonal variability and the
characteristic organism of each marine environment.

4. ENVIRONMENTAL IMPLICATIONS

The sources and molecular identities of DOM in the ocean are not yet fully understood, and the parameterizations of organic ligands in ocean biogeochemistry models still have significant uncertainties.⁷ It is known that the presence of DOM in seawater can affect the Fe(II) oxidation rate constant.<sup>17,34,37</sup> Previous laboratory studies<sup>35,71</sup> show that organic compounds can accelerate, reduce, or have no effect on Fe(II) oxidation kinetics. The *k'* variability in the presence of organic ligands suggests that the effect of DOM on Fe(II) oxidation is dependent on the molecules and their properties. This implies that specific ligands or groups of organic compounds have to be known to understand the effect that those ligands produce on Fe(II). Studies using environmental samples show that DOM from different sources presents varying effects on Fe(II) oxidation kinetics.<sup>17,18,37</sup> Oxidation rates for freshwater (e.g., river water and wastewater effluent) are generally higher than those for coastal waters.<sup>18</sup> The humic-type DOM (allochthonous origin) was defined as the key factor that accelerates the Fe(II) oxidation in freshwater samples. The lower oxidation rates of coastal seawater compared with those of freshwater and organic ligand-free seawater were thought to be associated with microbially derived autochthonous DOM. Few studies have considered the redox activity of compounds that makes up part of the DOM.<sup>72,73</sup> DOM is capable of acting as both an electron donor and an acceptor, keeping Fe in a redox cycle.<sup>2</sup> In any case, the results show the net effect of the different organic compounds that may be present in seawater.

This study remarks the important role of DOM in the Fe(II)
oxidation kinetic process and the consequences of Fe(II)
persistence in the marine environment. Although the

physicochemical variables pH and temperature control the Fe(II) oxidation rate in a non-oxygen limited medium, the biogeochemical context is important. The *k'* deviation from *k*'<sub>cal</sub> was explained through the spectral characterization of the organic matter. The observed variability of *k'* was correlated with the TDN and two spectral variables, *b*<sub>DOM</sub> and C1<sub>humic</sub>. However, it is necessary to indicate that the ratio of the organic ligand and Fe(II) will influence the fraction of Fe(II) complexed by organics and thereby the Fe(II) oxidation kinetics. Although we have obtained a correlation and it is an important advancement, studies related to the concentration ratio between Fe (II) and organics are necessary to the extent of the validity of the empirical equation to a range of Fe(II) concentration.

The nature of DOM in the medium may control the redox
cycle of Fe. In the continental margin, Fe(II) is influenced by
allochthonous contributions (i.e., rivers, marshes, and
estuaries). However, around the volcanic islands, DOM has
an autochthonous origin. The studied archipelagos presented
common characteristics: they have a volcanic origin and are
generally arid. There are no fluvial contributions in these
islands, but sporadic contributions through the ravines. Rainfall
is quite scarce in Cape Verde and the Canary Islands. This
produces a predominantly autochthonous DOM. Little is
known about the compounds or structures that make up
autochthonous CDOM in the ocean. This study demonstrated
that a higher degree of specificity in the OM characterization is
required if we want to determine the role that organic
compounds play in the persistence of Fe(II) in seawater and
the biogeochemical cycle of Fe..

■ ASSOCIATED CONTENT

● Supporting Information

The Supporting Information is available free of charge at
https://pubs.acs.org/doi/10.1021/acs.est.1c04512.

[Fe(II)] and ln [Fe(II)] versus time under different pH and T conditions; dependence of log *k'* with the pH; ln *k'* with 1/T; detailed sampling protocol analysis; figures for nutrients, DOC, POC, TDN, PON, CDOM, and FDOM; and TDN (μM), b (RU), and C1 (RU) data (PDF)

■ AUTHOR INFORMATION

Corresponding Author

J. Magdalena Santana-Casiano — Instituto de Oceanografia y Cambio Global, Universidad de Las Palmas de Gran Canaria, 35017 Las Palmas, Spain; ● orcid.org/0000-0002-7930-7683; Phone: +34 928 454448; Email: magdalena.santana@ulpgc.es

Authors

David González-Santana — Instituto de Oceanografia y Cambio Global, Universidad de Las Palmas de Gran Canaria, 35017 Las Palmas, Spain; Université de Brest, CNRS, IRD, Ifremer, LEMAR, F-29280 Plouzane, France

Quentin Devresse — GEOMAR—Helmholtz Centre for Ocean Research Kiel, 24105 Kiel, Germany

Helmke Hepach — GEOMAR—Helmholtz Centre for Ocean Research Kiel, 24105 Kiel, Germany

Carolina Santana-González — Instituto de Oceanografía y Cambio Global, Universidad de Las Palmas de Gran Canaria, 35017 Las Palmas, Spain
Birgit Quack — GEOMAR—Helmholtz Centre for Ocean Research Kiel, 24105 Kiel, Germany

Anja Engel — GEOMAR—Helmholtz Centre for Ocean Research Kiel, 24105 Kiel, Germany

Melchor González-Dávila — Instituto de Oceanografía y Cambio Global, Universidad de Las Palmas de Gran Canaria, 35017 Las Palmas, Spain

Complete contact information is available at:
https://pubs.acs.org/10.1021/acs.est.1c04512

## Funding

This work was financially supported by the ATOPFe project (CTM2017-83476-P) from the Ministerio de Ciencia e Innovación (Spain). This study also received funding from the European Union's Horizon 2020 research and innovation program under grant agreement no. 820989 (project COMFORT, *Our common future ocean in the Earth system-quantifying coupled cycles of carbon, oxygen, and nutrients for determining and achieving safe operating spaces with respect to tipping points*). The work reflects only the authors' view; the European Commission and their executive agency are not responsible for any use that may be made of the information the work contains.

## Notes

The authors declare no competing financial interest.

The Fe(II) database is available in the PANGAEA database as J.M.S.-C; B.Q. (2021): Fe(II) kinetic data from water samples during POSEIDON cruise POS533 (AIMAC). PANGAEA, https://doi.pangaea.de/10.1594/PANGAEA.934051. Online database OpenFluor is available at https://openfluor.lablicate.com. The PARAFAC model for this study is named "Macaronesia_POS533".

## ACKNOWLEDGMENTS

We want to express our gratitude to B.Q. from GEOMAR for inviting us to participate in the AIMAC project. Thanks to the master and crew of the *R/V Poseidon* for the support during the cruise. Special thanks go to Rui Caldeira, Cátia Azevedo, Claudio Cardoso, Ricardo Faria, and Jesus Reis from the Oceanic Observatory of Madeira for the CTD deployments and data. Dr. Gerd Krahmann (GEOMAR) for the evaluation of CTD-data. Kastriot Qelaj (GEOMAR) for analyzing the nutrient samples and Jon Roa, Tania Klüver, and Vivien Floren (all GEOMAR) for sample analysis of DOM and POM components.

## REFERENCES

(1) Boyd, P. W.; Ellwood, M. J.; Tagliabue, A.; Twining, B. S. Biotic and Abiotic Retention, Recycling and Remineralization of Metals in the Ocean. *Nat. Geosci.* **2017**, *10*, 167–173.

(2) Daugherty, E. E.; Gilbert, B.; Nico, P. S.; Borch, T. Complexation and Redox Buffering of Iron(II) by Dissolved Organic Matter. *Environ. Sci. Technol.* **2017**, *51*, 11096–11104.

(3) Voelker, B. M.; Morel, F. M. M.; Sulzberger, B. Iron Redox Cycling in Surface Waters: Effects of Humic Substances and Light. *Environ. Sci. Technol.* **1997**, *31*, 1004–1011.

(4) Shaked, Y.; Kustka, A. B.; Morel, F. M. M. A General Kinetic Model for Iron Acquisition by Eukaryotic Phytoplankton. *Limnol. Oceanogr.* **2005**, *50*, 872–882.

(5) Hassler, C. S.; van den Berg, C. M. G.; Boyd, P. W. Toward a Regional Classification to Provide a More Inclusive Examination of the Ocean Biogeochemistry of Iron-Binding Ligands. *Front. Mar. Sci.* **2017**, *4*, 19.

(6) Lønborg, C.; Carreira, C.; Jickells, T.; Álvarez-Salgado, X. A. Impacts of Global Change on Ocean Dissolved Organic Carbon (DOC) Cycling. *Front. Mar. Sci.* **2020**, *7*, 466.

(7) Pham, A. L. D.; Ito, T. Formation and Maintenance of the GEOTRACES Subsurface-Dissolved Iron Maxima in an Ocean Biogeochemistry Model. *Global Biogeochem. Cycles* **2018**, *32*, 932–953.

(8) Gledhill, M. The Organic Complexation of Iron in the Marine Environment: A Review. *Front. Microbiol.* **2012**, *3*, 69.

(9) Boye, M.; Aldrich, A.; van den Berg, C. M. G.; de Jong, J. T. M.; Nirmaier, H.; Veldhuis, M.; Timmermans, K. R.; de Baar, H. J. W. The Chemical Speciation of Iron in the North-East Atlantic Ocean. *Deep Sea Res., Part I* **2006**, *53*, 667–683.

(10) Fan, S.-M. Photochemical and Biochemical Controls on Reactive Oxygen and Iron Speciation in the Pelagic Surface Ocean. *Mar. Chem.* **2008**, *109*, 152–164.

(11) Santana-Casiano, J. M.; González-Dávila, M.; González, A. G.; Millero, F. J. Fe(III) Reduction in the Presence of Catechol in Seawater. *Aquat. Geochem.* **2010**, *16*, 467–482.

(12) González-Santana, D.; González-Dávila, M.; Lohan, M. C.; Artigue, L.; Planquette, H.; Sarthou, G.; Tagliabue, A.; Santana-Casiano, J. M. Variability in Iron (II) Oxidation Kinetics across Diverse Hydrothermal Sites on the Northern Mid Atlantic Ridge. *Geochim. Cosmochim. Acta* **2021**, *297*, 143–157.

(13) González-Santana, D.; Planquette, H.; Cheize, M.; Whitby, H.; Gourain, A.; Holmes, T.; Guyader, V.; Cathalot, C.; Pelleter, E.; Fouquet, Y.; Sarthou, G. Processes Driving Iron and Manganese Dispersal From the TAG Hydrothermal Plume (Mid-Atlantic Ridge): Results From a GEOTRACES Process Study. *Front. Mar. Sci.* **2020**, *7*, 568.

(14) Yücel, M.; Gartman, A.; Chan, C. S.; Luther, G. W. Hydrothermal Vents as a Kinetically Stable Source of Iron-Sulphide-Bearing Nanoparticles to the Ocean. *Nat. Geosci.* **2011**, *4*, 367–371.

(15) Moffett, J. W. Iron(II) in the world's oxygen deficient zones. *Chem. Geol.* **2021**, *580*, 120314.

(16) Canfield, D. E.; Stewart, F. J.; Thamdrup, B.; De Brabandere, L.; Dalsgaard, T.; Delong, E. F.; Revsbech, N. P.; Ulloa, O. A Cryptic Sulfur Cycle in Oxygen-Minimum-Zone Waters off the Chilean Coast. *Science* **2010**, *330*, 1375–1378.

(17) Rose, A. L.; Waite, T. D. Kinetics of Iron Complexation by Dissolved Natural Organic Matter in Coastal Waters. *Mar. Chem.* **2003**, *84*, 85–103.

(18) Lee, Y. P.; Fujii, M.; Kikuchi, T.; Terao, K.; Yoshimura, C. Variation of Iron Redox Kinetics and Its Relation with Molecular Composition of Standard Humic Substances at Circumneutral pH. *PLoS One* **2017**, *12*, No. e0176484.

(19) Ilbert, M.; Bonnefoy, V. Insight into the Evolution of the Iron Oxidation Pathways. *Biochim. Biophys. Acta, Bioenerg.* **2013**, *1827*, 161–175.

(20) Tagliabue, A.; Völker, C. Towards Accounting for Dissolved Iron Speciation in Global Ocean Models. *Biogeosciences* **2011**, *8*, 3025–3039.

(21) Cutter, G. A.; Andersson, P.; Codispoti, L.; Croot, P. L.; Place, P.; Hoe, T.; Kingdom, U.; Francois, R.; Sciences, O.; Lohan, M. C.; Circus, D.; Obata, H. *Sampling and Sample-Handling Protocols for GEOTRACES Cruises*, 2010.

(22) Mehlmann, M.; Quack, B.; Atlas, E.; Hepach, H.; Tegtmeier, S. Natural and Anthropogenic Sources of Bromoform and Dibromo-methane in the Oceanographic and Biogeochemical Regime of the Subtropical North East Atlantic. *Environ. Sci.: Processes Impacts* **2020**, *22*, 679–707.

(23) Santana-González, C.; González-Dávila, M.; Santana-Casiano, J. M.; Gladyshev, S.; Sokov, A. Organic Matter Effect on Fe(II) Oxidation Kinetics in the Labrador Sea. *Chem. Geol.* **2019**, *511*, 238–255.

(24) Santana-González, C.; Santana-Casiano, J. M.; González-Dávila, M.; Santana-del Pino, A.; Gladyshev, S.; Sokov, A. Fe(II) oxidation
kinetics in the North Atlantic along the 59.5° N during 2016. *Mar. Chem.* **2018**, *203*, 64–77.

(25) Millero, F. J.; Sotolongo, S.; Izaguirre, M. The Oxidation Kinetics of Fe(II) in Seawater. *Geochim. Cosmochim. Acta* **1987**, *51*, 793–801.

(26) Santana-González, C.; Santana-Casiano, J. M.; González-Dávila, M.; Fraile-Nuez, E. Emissions of Fe(II) and Its Kinetic of Oxidation at Tagoro Submarine Volcano, El Hierro. *Mar. Chem.* **2017**, *195*, 129–137.

(27) Millero, F. J. The pH of Estuarine Waters. *Limnol. Oceanogr.* **1986**, *31*, 839–847.

(28) King, D. W.; Lounsbury, H. A.; Millero, F. J. Rates and Mechanism of Fe(II) Oxidation at Nanomolar Total Iron Concentrations. *Environ. Sci. Technol.* **1995**, *29*, 818–824.

(29) Santana-Casiano, J. M.; González-Dávila, M.; Millero, F. J. The Oxidation of Fe (II) in NaCl–HCO3– and Seawater Solutions in the Presence of Phthalate and Salicylate Ions: A Kinetic Model. *Mar. Chem.* **2004**, *85*, 27–40.

(30) Stumm, W.; Lee, G. F. Oxygenation of Ferrous Iron. *Ind. Eng. Chem.* **1961**, *53*, 143–146.

(31) Kester, D. R.; Byrne, R. H., Jr.; Liang, Y.-J. Redox Reactions and Solution Complexes of Iron in Marine Systems; ACS Publications, 1975.

(32) Tamura, H.; Goto, K.; Yotsuyanagi, T.; Nagayama, M. Spectrophotometric determination of iron(II) with 1,10-phenanthroline in the presence of large amounts of iron(III). *Talanta* **1974**, *21*, 314–318.

(33) Emmenegger, L.; King, D. W.; Sigg, L.; Sulzberger, B. Oxidation Kinetics of Fe(II) in a Eutrophic Swiss Lake. *Environ. Sci. Technol.* **1998**, *32*, 2990–2996.

(34) Santana-Casiano, J. M.; González-Dávila, M.; Rodríguez, M. J.; Millero, F. J. The Effect of Organic Compounds in the Oxidation Kinetics of Fe(II). *Mar. Chem.* **2000**, *70*, 211–222.

(35) Santana-Casiano, J. M.; González-Dávila, M.; Millero, F. J. Oxidation of Nanomolar Levels of Fe(II) with Oxygen in Natural Waters. *Environ. Sci. Technol.* **2005**, *39*, 2073–2079.

(36) González-Dávila, M.; Santana-Casiano, J. M.; Millero, F. J. Competition between O2 and H2O2 in the Oxidation of Fe (II) in Natural Waters. *J. Solution Chem.* **2006**, *35*, 95–111.

(37) Rose, A. L.; Waite, T. D. Kinetic Model for Fe(II) Oxidation in Seawater in the Absence and Presence of Natural Organic Matter. *Environ. Sci. Technol.* **2002**, *36*, 433–444.

(38) Nelson, N. B.; Siegel, D. A. The Global Distribution and Dynamics of Chromophoric Dissolved Organic Matter. *Ann. Rev. Mar. Sci* **2013**, *5*, 447–476.

(39) Helms, J. R.; Stubbins, A.; Ritchie, J. D.; Minor, E. C.; Kieber, D. J.; Mopper, K. Absorption Spectral Slopes and Slope Ratios as Indicators of Molecular Weight, Source, and Photobleaching of Chromophoric Dissolved Organic Matter. *Limnol. Oceanogr.* **2008**, *53*, 955–969.

(40) Catalá, T. S.; Martínez-Pérez, A. M.; Nieto-Cid, M.; Álvarez, M.; Otero, J.; Emelianov, M.; Reche, I.; Arístegui, J.; Álvarez-Salgado, X. A. Dissolved Organic Matter (DOM) in the open Mediterranean Sea. *I. Basin-wide distribution and drivers of chromophoric DOM*. *Prog. Oceanogr.* **2018**, *165*, 35–51.

(41) De Haan, H.; De Boer, T. Applicability of Light Absorbance and Fluorescence as Measures of Concentration and Molecular Size of Dissolved Organic Carbon in Humic Lake Tjeukemeer. *Water Res.* **1987**, *21*, 731–734.

(42) Coble, P. G. Marine Optical Biogeochemistry: The Chemistry of Ocean Color. *Chem. Rev.* **2007**, *107*, 402–418.

(43) Zsolnay, A.; Baigar, E.; Jimenez, M.; Steinweg, B.; Saccomandi, F. Differentiating with Fluorescence Spectroscopy the Sources of Dissolved Organic Matter in Soils Subjected to Drying. *Chemosphere* **1999**, *38*, 45–50.

(44) Huguet, A.; Vacher, L.; Relexans, S.; Saubusse, S.; Froidefond, J. M.; Parlanti, E. Properties of Fluorescent Dissolved Organic Matter in the Gironde Estuary. *Org. Geochem.* **2009**, *40*, 706–719.

(45) Coble, P. G. Characterization of Marine and Terrestrial DOM in Seawater Using Excitation-Emission Matrix Spectroscopy. *Mar. Chem.* **1996**, *51*, 325–346.

(46) Ihaka, R.; Gentleman, R. R: A Language for Data Analysis and Graphics. *J. Comput. Graph Stat.* **1996**, *5*, 299–314.

(47) Pucher, M.; Wünsch, U.; Weigelhofer, G.; Murphy, K.; Hein, T.; Graeber, D. StaRdom: Versatile Software for Analyzing Spectroscopic Data of Dissolved Organic Matter in R. Water **2019**, *11*, 2366.

(48) Murphy, K. R.; Stedmon, C. A.; Wenig, P.; Bro, R. OpenFluor- An Online Spectral Library of Auto-Fluorescence by Organic Compounds in the Environment. *Anal. Methods* **2014**, *6*, 658–661.

(49) Lee, Y. P.; Fujii, M.; Kikuchi, T.; Natsuike, M.; Ito, H.; Watanabe, T.; Yoshimura, C. Importance of Allochthonous and Autochthonous Dissolved Organic Matter in Fe(II) Oxidation: A Case Study in Shizugawa Bay Watershed, Japan. *Chemosphere* **2017**, *180*, 221–228.

(50) Garg, S.; Jiang, C.; Waite, T. D. Impact of pH on Iron Redox Transformations in Simulated Freshwaters Containing Natural Organic Matter. *Environ. Sci. Technol.* **2018**, *52*, 13184–13194.

(51) Garg, S.; Rose, A. L.; Waite, T. D. Photochemical Production of Superoxide and Hydrogen Peroxide from Natural Organic Matter. *Geochim. Cosmochim. Acta* **2011**, *75*, 4310–4320.

(52) Samperio-Ramos, G.; Santana-Casiano, J. M.; González-Dávila, M. Effect of Ocean Warming and Acidification on the Fe(II) Oxidation Rate in Oligotrophic and Eutrophic Natural Waters. *Biogeochemistry* **2016**, *128*, 19–34.

(53) González, A. G.; Santana-Casiano, J. M.; Pérez, N.; González-Dávila, M. Oxidation of Fe(II) in Natural Waters at High Nutrient Concentrations. *Environ. Sci. Technol.* **2010**, *44*, 8095–8101.

(54) Boye, M.; Nishioka, J.; Croot, P.; Laan, P.; Timmermans, K. R.; Strass, V. H.; Takeda, S.; de Baar, H. J. W. Significant Portion of Dissolved Organic Fe Complexes in Fact Is Fe Colloids. *Mar. Chem.* **2010**, *122*, 20–27.

(55) González-Davila, M.; Santana-Casiano, J. M.; Millero, F. J. Oxidation of Iron (II) Nanomolar with H₂O₂ in Seawater. *Geochim. Cosmochim. Acta* **2005**, *69*, 83–93.

(56) Millero, F. J.; Sotolongo, S. The Oxidation of Fe(II) with H₂O₂ in Seawater. *Geochim. Cosmochim. Acta* **1989**, *53*, 1867–1873.

(57) Santana-Casiano, J. M.; González-Dávila, M.; Millero, F. J. The Role of Fe(II) Species on the Oxidation of Fe(II) in Natural Waters in the Presence of O₂ and H₂O₂. *Mar. Chem.* **2006**, *99*, 70–82.

(58) Heller, M. I.; Gaiero, D. M.; Croot, P. L. Basin Scale Survey of Marine Humic Fluorescence in the Atlantic: Relationship to Iron Solubility and H2O2. *Global Biogeochem. Cycles* **2013**, *27*, 88–100.

(59) González, A. G.; Santana-Casiano, J. M.; González-Dávila, M.; Pérez-Almeida, N.; Suárez de Tangil, M. Effect of Dunaliella Tertiolecta Organic Exudates on the Fe(II) Oxidation Kinetics in Seawater. *Environ. Sci. Technol.* **2014**, *48*, 7933–7941.

(60) Samperio-Ramos, G.; Santana-Casiano, J. M.; González-Dávila, M. Variability in the Production of Organic Ligands, by Synechococcus PCC 7002, under Different Iron Scenarios. *J. Oceanogr.* **2017**, *74*, 277–286.

(61) Lønborg, C.; Álvarez-Salgado, X. A. Tracing Dissolved Organic Matter Cycling in the Eastern Boundary of the Temperate North Atlantic Using Absorption and Fluorescence Spectroscopy. *Deep Sea Res., Part 1* **2014**, *85*, 35–46.

(62) Loginova, A. N.; Thomsen, S.; Engel, A. Chromophoric and Fluorescent Dissolved Organic Matter in and above the Oxygen Minimum Zone off Peru. *J. Geophys. Res.: Oceans* **2016**, *121*, 7973–7990.

(63) Catalá, T. S.; Reche, I.; Fuentes-Lema, A.; Romera-Castillo, C.; Nieto-Cid, M.; Ortega-Retuerta, E.; Calvo, E.; Álvarez, M.; Marrasé, C.; Stedmon, C. A.; Álvarez-Salgado, X. A. Turnover Time of Fluorescent Dissolved Organic Matter in the Dark Global Ocean. *Nat. Commun.* **2015**, *6*, 5986.

(64) Stedmon, C. A.; Thomas, D. N.; Granskog, M.; Kaartokallio, H.; Papadimitriou, S.; Kuosa, H. Characteristics of Dissolved Organic
Matter in Baltic Coastal Sea Ice: Allochthonous or Autochthonous
Origins? *Environ. Sci. Technol.* **2007**, *41*, 7273–7279.

(65) Martell, A. E.; Motekaitis, R. J.; Chen, D.; Hancock, R. D.; McManus, D. Selection of New Fe(III)/Fe(II) Chelating Agents as Catalysts for the Oxidation of Hydrogen Sulfide to Sulfur by Air. *Can. J. Chem.* **1996**, *74*, 1872–1879.

(66) Li, X.; Chen, Y.; Zhang, H. Reduction of Nitrogen-Oxygen Containing Compounds (NOCs) by Surface-Associated Fe(II) and Comparison with Soluble Fe(II) Complexes. *Chem. Eng. J.* **2019**, *370*, 782–791.

(67) Berman, T.; Bronk, D. Dissolved Organic Nitrogen: A Dynamic Participant in Aquatic Ecosystems. *Aquat. Microb. Ecol.* **2003**, *31*, 279–305.

(68) McCarthy, M. D.; Hedges, J. I.; Benner, R. Major Bacterial Contribution to Marine Dissolved Organic Nitrogen. *Science* **1998**, *281*, 231–234.

(69) Capone, D. G.; Ferrier, M. D.; Carpenter, E. J. Amino Acid Cycling in Colonies of the Planktonic Marine Cyanobacterium Trichodesmium Thiebautii. *Appl. Environ. Microbiol.* **1994**, *60*, 3989–3995.

(70) Glibert, P. M.; Bronk, D. A. Release of Dissolved Organic Nitrogen by Marine Diazotrophic Cyanobacteria, Trichodesmium Spp. *Appl. Environ. Microbiol.* **1994**, *60*, 3996–4000.

(71) Theis, T. L.; Singer, P. C. Complexation of Iron(II) by Organic Matter and Its Effect on Iron(II) Oxygenation. *Environ. Sci. Technol.* **1974**, *8*, 569–573.

(72) Garg, S.; Jiang, C.; David Waite, T. Mechanistic Insights into Iron Redox Transformations in the Presence of Natural Organic Matter: Impact of PH and Light. *Geochim. Cosmochim. Acta* **2015**, *165*, 14–34.

(73) Jiang, C.; Garg, S.; Waite, T. D. Hydroquinone-Mediated Redox Cycling of Iron and Concomitant Oxidation of Hydroquinone in Oxic Waters under Acidic Conditions: Comparison with Iron-Natural Organic Matter Interactions. *Environ. Sci. Technol.* **2015**, *49*, 14076–14084.

<figure><img src="image_14.png" /></figure>

Recommended by ACS

Oxidative Precipitation of Fe(II) in Porous Media:
Laboratory Experiment and Numerical Simulation

Zicheng Zhao, Ling Li, et al.

MARCH 14, 2023

ACS ES&T WATER

READ

<figure><img src="image_15.png" /></figure>

Effect of Stoichiometry on Nanomagnetite Sulfidation

Mingjun Nie, Juan Liu, et al.

FEBRUARY 06, 2023

ENVIRONMENTAL SCIENCE & TECHNOLOGY

READ

<figure><img src="image_16.png" /></figure>

Dissolved Organic Matter Photoreactivity Is Determined by
Its Optical Properties, Redox Activity, and Molecular
Composition

Stephanie M. Berg, Christina K. Remucal, et al.

APRIL 11, 2023

ENVIRONMENTAL SCIENCE & TECHNOLOGY

READ

<figure><img src="image_17.png" /></figure>

Vanadate Retention by Iron and Manganese Oxides

Macon J. Abernathy, Samantha C. Ying, et al.

AUGUST 05, 2022

ACS EARTH AND SPACE CHEMISTRY

READ

<figure><img src="image_18.png" /></figure>

Get More Suggestions >