# Winter Biogeochemical Cycling of Dissolved and Particulate Cadmium in the Indian Sector of the Southern Ocean (GEOTRACES GIpr07 Transect)

Ryan Cloete, Jean Loock, Natasha van Horsten, Susanne Fietz, Thato Mtshali, Hélène Planquette, Alakendra Roychoudhury

► To cite this version:

Ryan Cloete, Jean Loock, Natasha van Horsten, Susanne Fietz, Thato Mtshali, et al.. Winter Biogeochemical Cycling of Dissolved and Particulate Cadmium in the Indian Sector of the Southern Ocean (GEOTRACES GIpr07 Transect). Frontiers in Marine Science, 2021, 8, 10.3389/fmars.2021.656321. [https://doi.org/10.3366/fms.2021.656321]

HAL Id: hal-03366004

https://hal.science/hal-03366004

Submitted on 21 Oct 2021

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
<figure><img src="image_2.png" /></figure>

Check for
updates

# Winter Biogeochemical Cycling of Dissolved and Particulate Cadmium in the Indian Sector of the Southern Ocean (GEOTRACES Glpr07 Transect)

Ryan Cloete<fnref n="1" />*, Jean C. Loock<fnref n="1" />, Natasha R. van Horsten<fnref n="1" />,²,³, Susanne Fietz<fnref n="1" />, Thato N. Mtshali⁴, Hélène Planquette³ and Alakendra N. Roychoudhury<fnref n="1" />

OPEN ACCESS

**Edited by:**
Sunil Kumar Singh,
Physical Research Laboratory, India

**Reviewed by:**
Gregory De Souza,
ETH Zürich, Switzerland
Tarun K. Dalai,
Indian Institute of Science Education
and Research Kolkata, India

*Correspondence:
Ryan Cloete
15994619@sun.ac.za

**Specialty section:**
This article was submitted to
Marine Biogeochemistry,
a section of the journal
Frontiers in Marine Science

**Received:** 20 January 2021

**Accepted:** 05 July 2021

**Published:** 26 July 2021

**Citation:**
Cloete R, Loock JC, van Horsten NR, Fietz S, Mtshali TN, Planquette H and Roychoudhury AN (2021) Winter Biogeochemical Cycling of Dissolved and Particulate Cadmium in the Indian Sector of the Southern Ocean (GEOTRACES Glpr07 Transect). Front. Mar. Sci. 8:656321. doi: 10.3389/fmars.2021.656321

Winter distributions of dissolved cadmium (dCd) and particulate cadmium (pCd) were measured for the first time in the Indian sector of the Southern Ocean thereby contributing a unique spatial and seasonal dataset. Seven depth profiles, between 41°S and 58°S, were collected along the 30°E longitude during the 2017 austral winter to investigate the biogeochemical cycling of cadmium during a period characterized by contrasting upper water column dynamics compared to summer. Our results support an important role for biological uptake during winter months albeit weaker compared to summer. Distinct, biologically driven changes in cadmium cycling across the transect were observed. For example, surface ratios of pCd to phosphorus (P; pCd:P) increased from 0.37 to 1.07 mmol mol⁻<fnref n="1" /> between the subtropical zone (STZ) and the Antarctic zone (AAZ) reflecting increased Cd requirements for diatoms at higher latitudes which, in turn, was driven by a complex relationship between the availability of dCd and dissolved iron (dFe), zinc (dZn) and manganese (dMn). Vertical profiles of pCd:P displayed near-surface maxima consistent with (1) P occurring in two phases with different liabilities and the liability of Cd being somewhere in-between and (2) increasing dCd to phosphate (PO₄; dCd:PO₄) ratios with depth at each station. North of the Antarctic Polar Front (APF), a secondary, deeper pCd:P maximum may reflect an advective signal associated with northward subducting Antarctic Intermediate Water (AAIW). The strong southward increase in surface dCd and dCd:PO₄, from approximately 10–700 pmol kg⁻<fnref n="1" /> and 40–400 µmol mol⁻<fnref n="1" />, respectively, reflected the net effect of preferential uptake and regeneration of diatoms with high Cd content and the upwelling of Cd enriched water masses in the AAZ. Furthermore, distinct dCd versus PO₄ relationships were observed in each of the intermediate and deep water masses suggesting that dCd and PO₄ distributions at depth are largely the result of physical water mass mixing.

**Keywords:** trace metal, phytoplankton, GEOTRACES, Southern Ocean, dissolved cadmium, particulate cadmium, diatoms, winter

<footnote n="1">¹Centre for Trace Metal and Experimental Biogeochemistry (TracEx), Department of Earth Sciences, Stellenbosch University, Stellenbosch, South Africa, ²Southern Ocean Carbon & Climate Observatory, Council for Scientific and Industrial Research (CSIR), Stellenbosch, South Africa, ³CNRS, IRD, Ifremer, LEMAR, Univ Brest, Plouzané, France, ⁴Oceans and Coast, Department of Environment, Forestry and Fisheries, Cape Town, South Africa</footnote>
# INTRODUCTION

It is well established that numerous trace metals are required
nutrients for phytoplankton and that trace metal availability
in the world's oceans can influence phytoplankton growth and
community structure (Twining and Baines, 2013) with potential
consequences for the global carbon cycle. Data pertaining to
the trace metal cadmium (Cd) are scarce, particularly in the
Southern Ocean (Schlitzer et al., 2018), however, it has allowed
for a basic model of Cd cycling to be developed. The dominant
external sources of Cd to the surface ocean are riverine discharge
(Lambelet et al., 2013), natural and anthropogenic derived
atmospheric deposition (Duce et al., 1991) and seasonal ice
melt in the polar oceans (Hölemann et al., 1999). The internal
cycling of Cd in the Southern Ocean is further controlled by
complexation with natural organic ligands in surface waters
(Ellwood, 2004), year-round mixing between water masses
with different preformed Cd concentrations and upwelling at
higher latitudes (Abouchami et al., 2014; Sieber et al., 2019).
Additionally, deep winter mixing (or entrainment) is a seasonally
constrained physical supply mechanism whereby deeper mixing
in winter compared to summer allows access to nutrient rich
subsurface water masses found below the summer stratified layer
(Ellwood, 2008; Tagliabue et al., 2014).

In addition to its sources, Cd cycling can be extended
to differentiate between various source pools of trace metals.
The dissolved Cd pool (<0.2 µm pore size filtered seawater),
denoted as “d,” e.g., dCd hereafter, is recognized as the nutrient
reservoir available for phytoplankton uptake. The particulate
Cd pool (>0.45 µm pore size filter), denoted as “p,” e.g.,
pCd hereafter, is not considered bioavailable to phytoplankton
yet contributes significantly to the total Cd inventory in the
ocean (Turekian, 1977; Lam and Marchal, 2015). A host of
biophysicochemical processes control the cycling and exchange
between dissolved and particulate cadmium in the ocean. In the
surface ocean, uptake of dCd by phytoplankton and recycling
upon cell death are the primary processes that facilitate the
exchange between inorganic dCd and organic pCd. Subsequent
processes, such as zooplankton grazing, result in the formation
of larger particles which either sink and contribute to the
vertical export of pCd, or disaggregate and return Cd to
the dissolved pool through remineralisation, desorption
and dissolution, either in the euphotic or mesopelagic zone (Lam
and Marchal, 2015). Scavenging (or adsorption) of dCd onto
organic or hydrothermally sourced particles may enhance export
of Cd through formation of more rapidly sinking particulate
phases (Conway and John, 2015). Similarly, precipitation of solid
Cd sulfides in euxinic microenvironments (e.g., within particles)
increases pCd export in localized regions (Janssen et al., 2014).
As a result of the interplay of all these processes, increases in
water column dCd typically co-occur with decreases in pCd and
vice versa (Quay et al., 2015). However, dissolved and particulate
phases are not exact mirror images of each other considering
that water column processes lead to differences in their respective
oceanic residence times, e.g., years to thousands of years for
dissolved phases and days to months for particulate phases
(Lam and Marchal, 2015).

The first reliable open ocean dCd measurements (Boyle et al., 1976; Bruland et al., 1978) were well correlated with the major nutrient phosphate (PO₄) and showed a surface depletion (<2 pmol kg⁻¹) relative to the deeper waters (>1,000 pmol kg⁻¹), as is typical of nutrient-type behavior. Indeed, a diatom specific carbonic anhydrase (CA) enzyme containing Cd (CDCA) has been identified, supporting Cd uptake (Lane and Morel, 2000; Lane et al., 2005). Furthermore, Cd may stimulate the growth of zinc (Zn)-limited diatoms through a substitutive mechanism (Price and Morel, 1990). The substitution of Cd for Zn under low Zn conditions has also been observed in other phytoplankton groups which do not contain the CDCA enzyme, notably coccolithophores and chlorophytes (Lee and Morel, 1995; Xu et al., 2007). These findings provide a potential mechanism responsible for the observed coupling of Cd and PO₄. However, there are divergences in the global oceanic Cd-PO₄ relationship, with a pronounced break ('kink') in the slope at PO₄ ~1.3 µmol kg⁻¹, indicating a net preferential removal of Cd from surface waters (Elderfield and Rickaby, 2000). Many studies have investigated the assimilation of PO₄ and Cd in phytoplankton and revealed different factors influencing the uptake of Cd relative to PO₄. For example, a diatom dominated phytoplankton assemblage showed decreases in cellular Cd in response to the addition of the key micronutrient iron (Fe) suggesting Cd uptake is directly linked to Fe availability (Cullen et al., 2003). Subsequent experiments have shown that increased cellular Cd in response to low Fe availability may follow one of two potential mechanisms. First, biodilution occurs in Fe-limited phytoplankton where reduced growth rates result in slower accumulation of P while Cd uptake remains relatively constant ultimately causing elevated pCd:P signatures (Cullen, 2006). Second, increased dCd uptake can occur through non-specific Fe transporters when phytoplankton are exposed to low levels of Fe (Lane et al., 2008, 2009). In addition to Fe, increased cellular Cd is also inversely related to the availability of other essential micronutrients such as Zn and manganese (Mn; Sunda and Huntsman, 2000) and may be further modulated by dissolved carbon dioxide (CO₂) concentrations in surface waters (Cullen et al., 1999; Cullen and Sherrell, 2005). Here, and in addition to the previously mentioned biological uses, Cd may be taken up passively with these divalent metals and bound within the cell to prevent toxicity (Horner et al., 2013). In sum, the role of Cd as a micronutrient is dependent on a number of oceanographic properties including the availability of other bio-essential trace metal concentrations, notably Fe, Zn, and Mn.

Understanding the biogeochemical cycling of Cd in the
Southern Ocean is particularly important given that water masses
formed in this region are the primary communication channels
to the low latitude oceans (Broecker, 1991). The far reaching
influence of the Southern Ocean has been shown in terms of
macronutrient distributions (Frew and Hunter, 1992; Sarmiento
et al., 2004) as well as for Cd (Baars et al., 2014; Xie et al.,
2015; Middag et al., 2018; Sieber et al., 2019) highlighting
the combination of biological uptake in the surface Southern
Ocean, lateral transport and mixing between water masses with
different preformed nutrient signatures. While significant strides
have been made in identifying and understanding processes
controlling oceanic Cd distributions, spatial and seasonal gaps in the global database still exist. For example, the Indian sector of the Southern Ocean is particularly data scarce (Janssen et al., 2020). Additionally, most observations to date represent spring/summer seasons with few winter Cd measurements available (Ellwood, 2008). We therefore have little understanding of Cd cycling over seasonal cycles and of the conditions that characterize the Southern Ocean during winter, a period which may be critical in supplying Cd and other nutrients (e.g., Fe) to sustain phytoplankton growth over the subsequent spring and summer seasons (Tagliabue et al., 2014). Under the GEOTRACES framework, this study aims to improve our understanding of Cd biogeochemical cycling in the ocean by contributing the first winter dissolved and particulate measurements of Cd from the western Indian sector (WOCE meridional section IO6S; GEOTRACES section G1pr07) of the Southern Ocean.

# METHODOLOGY

## Sample Collection

Seawater samples were collected on-board the South African (SA) Agulhas II polar research vessel during the 2017 Winter Cruise (WC-17; 28/06/2017–13/07/2017). The cruise transect followed

the World Ocean Circulation Experiment (WOCE) IO6 transect along the 30°E meridian in the Indian sector of the Southern Ocean (Figure 1A). Seven stations, consisting of four deep and three shallow, were sampled between 41°00′S and 58°30′S allowing observations over important frontal systems and in different water masses. “Deep” stations were sampled from near the surface (25 m) to near sea-floor (<4,500 m) while “shallow” stations were sampled between near surface and 1,500 m depth.

Seawater samples were collected following a strict 'trace metal clean' GEOTRACES protocol (Cutter et al., 2017). A vertical profile sampling method was employed, using 24 trace clean, internally Teflon-coated polyvinyl chloride (PVC) 12 L GO-FLO sampling bottles mounted on an epoxy-coated aluminum frame conductivity-temperature-depth (CTD) rosette. The rosette also housed sensors used to obtain vertical profiles of temperature, salinity, fluorescence and photosynthetically active radiation (PAR). A Kevlar hydrowire with internal signal cables allowed for the transfer of data between the CTD and the on-board control room. The GO-FLO bottles were triggered at predetermined depths during the up cast. Directly upon recovery of the rosette, the GO-FLO bottles were covered in a PVC plastic wrap in addition to their ends being covered in plastic (PVC) shower caps, and were transported into a class 100 clean container lab for sub-sampling. Samples for dissolved trace metal determination

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

**FIGURE 1** | (A) Sample station locations during the 2017 Winter Cruise along the GEOTRACES G1pr07 transect, equivalent to World Ocean Circulation Experiment (WOCE) IO6 transect. Positions of the frontal systems intersected, and corresponding oceanic zones, are shown as well as the main surface currents. Acronyms in alphabetical order are as follows; AAZ, Antarctic Zone; AC, Agulhas Current; ACC, Antarctic Circumpolar Current; APF, Antarctic Polar Front; MIZ, Marginal Ice Zone; PFZ, Polar Frontal Zone; SACCF, Southern Antarctic Circumpolar Current Front; SAF, Subantarctic Front SAZ: Subantarctic Zone; SBdy, Southern Boundary; STF, Subtropical Front; STZ, Subtropical Zone. (B) Scatter plot of potential temperature (°C) against practical salinity used to characterize the water mass regime. Acronyms in alphabetical order are as follows; AABW, Antarctic Bottom Water; AAIW, Antarctic Intermediate Water; AASW, Antarctic Surface Water; LCDW, Lower Circumpolar Deep Water; NADW, North Atlantic Deep Water; SAMW, Subantarctic Mode Water; SASW, Subantarctic Surface Water; SICW, South Indian Central Water; STSW, Subtropical Surface Water; UCDW, Upper Circumpolar Deep Water. Figures constructed using Ocean Data View (ODV; Schlitzer, 2020).
were collected into 125 mL acid-cleaned low density polyethylene (LDPE) bottles after online filtration through 0.2 µm pore size capsule filters (Sartobran) and under slight N₂ gas (0.5 bar, 99.9999% N₂, BIP technology) overpressure. Samples were acidified (pH~1.7) under laminar flow using hydrochloric acid (30% Ultrapur® HCl, Merck) and stored at room temperature for later analysis. Thereafter, previously cleaned filters (25mm diameter, Supor, 0.45 µm pore size) were mounted on acid-washed Swinnex (Millipore) filter holders which were then attached to the GO-FLO spigots and used for filtration (5–10 L) of samples for particulate trace metal determination (Planquette and Sherrell, 2012). After filtration, filter holders were removed and disassembled before the sampled filter was transferred to an acid-washed polystyrene Petri-dish using clean plastic forceps. Sampled filters were subsequently frozen at –20°C for transport back to the shore lab. All handling of filters, filter holders and samples was done under a class 100 laminar flow hood.

## Dissolved Cadmium (dCd) Determination
### seaFAST and ICP-MS

Samples were analyzed by Inductively Coupled Plasma Mass Spectrometry (ICP-MS, Agilent 7900), after offline preconcentration by Solid Phase Extraction (SPE) on a seaFAST SC-4 DX module (Cloete et al., 2019). Briefly, 10 mL of seawater was taken up by the seaFAST (Elemental Scientific Inc.) module, buffered with an ammonium acetate buffer to a pH of 6.0 ± 0.2, and loaded onto a high affinity metal chelating resin column (Nobias PA1). Here the metal ions including Cd were bound to the resin and separated from the seawater matrix elements, e.g., Na, Mg, and Cl which passed through the column. The metal ions were subsequently eluted (in a 2% HNO₃ solution; Ultrapur®, Merck) from the resin column in low volumes (250 µL) resulting in a preconcentration factor of 40. On the ICP-MS, samples were introduced using a low self-aspirating perfluoroalkoxy (PFA) nebulizer with a flow rate of 0.2 mL min⁻¹. Cd was measured using the instruments Octopole Reaction System (ORS) in helium (He) collision mode to reduce plasma and matrix based interferences, e.g., molybdenum oxide (MoO⁺) interference on Cd, although the molybdenum concentration was extensively reduced by the seaFAST matrix removal system. In order to quantify any remaining MoO⁺ interference on Cd, a Cd-free Mo standard at typical seawater concentrations (10–100 nmol kg⁻¹) was analyzed. Instrument parameters were optimized to reduce the interference to a maximum of 1 ppt Cd (~9 pmol kg⁻¹). The Cd-free Mo standard was included throughout the analysis to monitor MoO⁺ production and correct the Cd results. Online internal standard addition for drift correction was not possible using the self-aspirating nebulizer. Instead, instrument drift was monitored by running a multi-element calibration standard (MES; supplied by Inorganic Ventures) every six samples. Where drift exceeded 5% relative to the starting concentration of the MES for a specific element, a drift correction was applied using the formula: 2*ConcMES_start/(ConcMES_a + ConcMES_b)*ConcSample, where a and b are the MES before and after each set of six samples.

## Accuracy and Precision

The accuracy of the method was verified by the analysis of numerous community reference seawater samples for dCd (Table 1). For the GEOTRACES SAFe D2, GSC and GSP reference seawater, our dCd concentrations were within analytical uncertainty estimates of consensus values confirming the method's accuracy over a wide concentration range. Our mean dCd value from the analysis of the NASS-7 certified reference material (CRM) was within analytical uncertainty of the certified dCd value. To monitor ICP-MS precision, the Winter Indian Southern Ocean Seawater (WISOS) internal control seawater was placed in the analysis sequence and results compared against each other and the calibrated mean (Table 1). The WISOS internal control was collected in bulk (20 L) from surface water at 55°S; 28°E during the 2017 Winter Cruise after filtration through a 0.2 µm filter into an acid cleaned LDPE container. Sub-samples (1 L) of the bulk seawater were then acidified (pH = 1.7) using hydrochloric acid (30% HCl, Ultrapur®, Merck). The calibrated dCd mean (888 ± 21 pmol kg⁻¹) concentration was established by replicate analysis (n = 10) of the WISOS internal control seawater in conjunction with the SAFe and NASS-7 reference materials. Results from the analysis of the in-run WISOS control sample were 899 ± 27 pmol kg⁻¹ (n > 30) confirming the method's precision over multiple analytical sessions.

<table><caption>**TABLE 1** | Results for dCd from the ICP-MS analysis of GEOTRACES SAFE D2, GSC, and GSP seawater reference materials (https://www.geotraces.org/standards-and-reference-materials/), NASS-7 certified seawater reference material and our own WISOS internal reference seawater are compared with respective consensus values.</caption><thead><tr><th></th><th>dCd<br>pmol kg<sup>-1</sup></th></tr></thead><tbody><tr><td><strong>SAFe D2</strong></td><td></td></tr><tr><td>Consensus</td><td>986 ± 23</td></tr><tr><td>Measured (n = 4)</td><td>982 ± 10</td></tr><tr><td><strong>GSC</strong></td><td></td></tr><tr><td>Consensus</td><td>355 ± 21</td></tr><tr><td>Measured (n = 5)</td><td>345 ± 21</td></tr><tr><td><strong>GSP</strong></td><td></td></tr><tr><td>Consensus</td><td>1.95 ± 1,95</td></tr><tr><td>Measured (n = 5)</td><td>4.07 ± 0.48</td></tr><tr><td><strong>NASS-7</strong></td><td></td></tr><tr><td>Certified</td><td>141.50 ± 14.40</td></tr><tr><td>Measured (n = 5)</td><td>132.69 ± 3.08</td></tr><tr><td><strong>WISOS reference material</strong></td><td></td></tr><tr><td>TracEx calibration (n = 10)</td><td>888.02 ± 21.22</td></tr><tr><td>Measured (n = 5)</td><td>899.41 ± 27.55</td></tr><tr><td><strong>Blanks</strong></td><td></td></tr><tr><td>Method (n = 5)</td><td>1.22 ± 0.30</td></tr><tr><td>Limit of detection (n = 5)</td><td>0.89</td></tr></tbody></table>

Consensus values for SAFe D2 as of 2013 and GSC and GSP as of 2019. Calibrated mean values for WISOS calculated by repeat analysis (n > 10) of a large volume surface seawater sample collected from 55°S; 28°E. Instrument and method blanks as well as ICP-MS limit of detection are also shown.
## Blanks and Limits of Detection

The method blank (seaFAST and ICP-MS) was determined by preconcentrating and analyzing a solution of 2% HNO₃ (Ultrapur®, Merck) in ultrapure deionized water (the same composition as the eluent used in sample pre-concentration). The method blank was subjected to the same pre-concentration procedure as the seawater samples. The mean method blank was 1.22 ± 0.30 pmol kg⁻¹ (n = 5). The limit of detection (LOD) of the instrument, calculated as 3 times the standard deviation of the preconcentrated blanks, was 0.89 pmol kg⁻¹ (n = 5).

## Data Processing

Samples for dCd were analyzed in duplicate. The final values represent the mean of duplicate measurements. In the two cases where the percentage relative standard deviation (%RSD) between the duplicate measurements was >10%, one of the values was deemed a suspected outlier and not used further. This was determined by curve fitting the data points based on the values above (shallower depth) and below (deeper depth) the suspect value as well as by comparison with other parameters (salinity, temperature, and nutrients) measured from the same GO-FLO bottle.

## Phosphate (PO₄) Determination

Seawater samples for macronutrient analysis were collected from GO-FLO bottles, filtered immediately after collection using a 0.2 µm pore size syringe filter (Anotop) into 50 mL Falcon® tubes and frozen at -20°C. Post cruise, PO₄ was determined manually by colorimetric method (Grasshoff, 1983) at the Marine Biogeochemistry Laboratory at the University of Cape Town (MBL-UCT; Weir et al., 2020). The analytical error of the PO₄ measurements was ±0.06 µmol kg⁻¹ (Weir et al., 2020).

## Particulate Cadmium (pCd) and Phosphorus (P) Determination

The pCd and P samples were processed and analyzed 6 months after sample collection according to established methods (Planquette and Sherrell, 2012). Filters containing the samples were acid reflux digested at 130°C in acid-cleaned Savillex screw-cap vial. Archive solutions were stored in 3 mL of 0.12 M HNO₃ (Ultrapur®, Merck), of which 250 µL was diluted up to 2 mL for analysis by sector field inductively coupled plasma mass spectrometry (SF-ICP-MS, Element XR Thermo Scientific). Samples were spiked with 1 µg L⁻¹ indium (In) as an internal standard to correct for instrument drift. Three certified reference materials (PACS 3, MESS 4 and BCR 414) were processed as samples and analyzed for pCd and P to assess the accuracy of the methodology (Table 2). The certified reference materials yielded mean percentage recoveries of 108, 188, and 115% for pCd and 100, 106, and 95% for P (PACS 3, MESS 4, and BCR 414, respectively). Blanks were determined by digesting and analyzing unused acid washed filters (n = 5). The detection limits, defined as three times the standard deviation of the blanks, for pCd and P were 0.07 pmol kg⁻¹ and 0.03 nmol kg⁻¹, respectively (n = 5).

<table><caption>**TABLE 2** | Results for pCd and P from the ICP-MS analysis of PACS 3, MESS 4, and BCR 414 certified references materials.</caption><thead><tr><th></th><th>pCd<br>µmol kg<sup>-1</sup></th><th>P<br>mmol kg<sup>-1</sup></th></tr></thead><tbody><tr><td><strong>PACS 3</strong></td><td></td><td></td></tr><tr><td>Certified</td><td>19.84 ± 1.16</td><td>29.38 ± 2.91</td></tr><tr><td>Measured</td><td>21.32 ± 0.44</td><td>29.41 ± 0.31</td></tr><tr><td>% recovery</td><td>108</td><td>100</td></tr><tr><td><strong>MESS 4</strong></td><td></td><td></td></tr><tr><td>Certified</td><td>2.49 ± 0.36</td><td>33.58 ± 5.17</td></tr><tr><td>Measured</td><td>4.69 ± 1.21</td><td>35.58 ± 0.51</td></tr><tr><td>% recovery</td><td>188</td><td>106</td></tr><tr><td><strong>BCR 414</strong></td><td></td><td></td></tr><tr><td>Certified</td><td>3.41 ± 0.12</td><td>428.84 ± 51.20</td></tr><tr><td>Measured</td><td>3.93 ± 0.60</td><td>409.34 ± 52.14</td></tr><tr><td>% recovery</td><td>115</td><td>95</td></tr><tr><td>Filter blank (n = 5)</td><td>0.09 pmol kg<sup>-1</sup></td><td>0.08 nmol kg<sup>-1</sup></td></tr><tr><td>Limit of detection</td><td>0.07 pmol kg<sup>-1</sup></td><td>0.03 nmol kg<sup>-1</sup></td></tr></tbody></table>

Filter blanks as well as ICP-MS limit of detection are also shown.

# RESULTS

## Hydrographic Setting

### Frontal Positions

The transect intersected several frontal systems which, in turn, separated the study area into physically and biogeochemically distinct regions (Figure 1A). Frontal positions were determined previously (Weir et al., 2020) following temperature, salinity and oxygen criteria outlined in the literature (Orsi et al., 1995; Belkin and Gordon, 1996; Pollard et al., 2002). The Subtropical Front (STF) located at 42.4°S, defined the boundary between the Subtropical Zone (STZ) to the north and the Subantarctic Zone (SAZ) to the south. The Subantarctic Front (SAF) was identified at 46.2°S and separated the SAZ in the north from the Polar Frontal Zone (PFZ) to the south which was bounded in the south by the Antarctic Polar front (APF) at 49.3°S. Further south, the Antarctic Zone (AAZ) extended until the Southern Boundary (SBdy) at 58.5°S, incorporating the southern Antarctic Circumpolar Current Front (sACCF) at 56.5°S. The Marginal Ice Zone (MIZ), defined by 30% ice cover, was encountered at 61.7°S (de Jong et al., 2018). The location of the sample stations meant that, from north to south, station IO08 was found in the STZ, stations IO07 and IO06 in the SAZ, station IO05 in the PFZ and stations IO04, IO02, and IO01 were within the AAZ.

## Water Mass Classification

Water masses (Figure 1B) were identified using potential temperature, salinity and oxygen criteria outlined in the literature (Orsi et al., 1995; Anilkumar et al., 2006). The high temperature and salinity signature of surface waters at the northernmost station reflected Subtropical Surface Water (STSW; >12°C, S > 35.3). Moving southward, the decrease in surface water temperature marked the transition from STSW to Subantarctic Surface Water (SASW; 7-12°C) and Antarctic Surface Water (AASW; -1-2.5°C). Below the surface
waters in the STZ, a layer of South Indian Central Water (SICW) was sampled. This water mass was underlain by the northwards subducting Subantarctic Mode Water (SAMW; 12–15°C; 26.2–26.6 kg m⁻³) which was sampled between ~300–500 m at 41°S. Similarly, Antarctic Intermediate Water (AAIW; 2.5–5.0°C, 26.7–27.4 kg m⁻³) was found subducting northwards below SAMW. In contrast to the intermediate waters, sampled deep waters consisted predominantly of southward upwelling Upper Circumpolar Deep Water (UCDW) and the deeper lying North Atlantic Deep Water (NADW) and Lower Circumpolar Deepwater (LCDW). The UCDW (27.2–27.8 kg m⁻³) was further characterized by oxygen minima (155–180 µmol kg⁻¹) and PO₄ maxima (2.1–2.6 µmol kg⁻¹). North of the APF, NADW (27.0–27.9 kg m⁻³) was identified by high salinity (S > 34.7) and deep PO₄ minima (1.9–2.3 µmol kg⁻¹) while LCDW (27.8–27.9 kg m⁻³), sampled below UCDW and NADW, was characterized by lower temperatures

(~1°C). Finally, Antarctic Bottom Waters (AABW; 27.8–27.9 kg m⁻³), distinguished by characteristic local temperature minimums (–0.5–1.0°C) were sampled at depth in the AAZ (IO02 and IO04).

## Surface Mixed Layer (SML) and Euphotic Depth

The SML at each station was defined as the depth interval between the sea surface and the Mixed Layer Depth (MLD) which was calculated as the depth (between 10 and 400 m) at which the Brunt Väisälä frequency (BVF) squared reaches a maximum while the euphotic zone depth was approximated as the depth at which 1% of surface PAR was measured (Weir et al., 2020). The MLD was shallowest at 122 m in the STZ, deepest in the SAZ at 190–220 m and intermediate in the PFZ (134 m) and AAZ (125–161 m). The euphotic zone depth was always shallower than the MLD and generally shoaled northwards from >56 m in the AAZ to 18 m in the STZ.

<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Section plots of (A) dCd, (B) PO₄ overlain by salinity, (C) pCd and (D) P. Each figure is separated into upper 1,000 m (**upper panel**) and below 1,000 m (**bottom panel**). Frontal positions, corresponding zones and water masses are shown (see **Figure 1** for acronyms). Figure constructed using Ocean Data View (ODV; Schlitzer, 2020). PO₄ data from Weir et al. (2020).</figcaption></figure>
# Distribution of Dissolved Elements

## Dissolved Cadmium

Depth profiles of dCd concentrations along the transect are shown in **Figure 2A** and **Supplementary Figure 1**. At each station, dCd profiles displayed typical nutrient-type behavior consistent with low surface concentrations and increased concentrations at depth. Across the transect, dCd was fairly uniform in the SML and was characterized by a strong north-south gradient which was associated with the frontal positions. Within the SML, dCd concentrations were $10 \pm 2$ pmol kg$^{-1}$ in the STZ and increased rapidly across the APF to $749 \pm 16$ pmol kg$^{-1}$ in the southern AAZ. The difference between surface and deep water concentrations was most apparent at the northern stations and decreased southwards. For example, the mean dCd in STSW was $33 \pm 21$ pmol kg$^{-1}$, ~20 times lower than the concentration of $688$ pmol kg$^{-1}$ ($n = 1$) in UCDW. At the southernmost station, this factor was ~1.4 between AASW ($713 \pm 30$ pmol kg$^{-1}$) and UCDW ($1018 \pm 65$ pmol kg$^{-1}$). Mid-depth maxima in dCd varied in concentration and with depth and coincided with local maxima in PO$_4$ and minima in O$_2$. The mid-depth dCd maximum ranged between $700$ pmol kg$^{-1}$ at 1,250 m depth in the STZ and $1,150$ pmol kg$^{-1}$ at 200 m depth in the AAZ. Below the mid-depth peak, concentrations were typically lower in NADW ($598 \pm 47$ pmol kg$^{-1}$) and showed increases in the LCDW ($800 \pm 120$ pmol kg$^{-1}$) and AABW ($856 \pm 88$ pmol kg$^{-1}$).

## Dissolved Phosphate

Major nutrient concentrations have been reported previously (Weir et al., 2020). Here we briefly describe the distribution of PO$_4$ (**Figure 2B** and **Supplementary Figure 2**) because of the extensive reference made to the dCd-PO$_4$ relationship. There was a strong southward increase in SML PO$_4$ concentrations, from $0.28 \pm 0.08$ µmol kg$^{-1}$ in STSW to $1.85 \pm 0.09$ µmol kg$^{-1}$ in AASW. Below the SML, PO$_4$ increased to maximum concentrations in UCDW at each station. The intermediate depth PO$_4$ maximum was most evident around the APF where PO$_4$ reached $2.43 \pm 0.16$ µmol kg$^{-1}$. The underlying NADW ($2.10 \pm 0.16$ µmol kg$^{-1}$) and LCDW ($2.14 \pm 0.14$ µmol kg$^{-1}$) were characterized by local PO$_4$ minima before slight increases in PO$_4$ in AABW ($2.20 \pm 0.02$ µmol kg$^{-1}$).

# Distribution of Particulate Elements

## Particulate Cadmium

The depth profiles of pCd (**Figure 2C** and **Supplementary Figure 1**) generally showed the opposite vertical distribution in comparison to dCd. Highest pCd concentrations were typically found within the surface water masses at each station, while strong negative gradients through subsurface water masses resulted in pCd concentrations having uniformly low ($<1$ pmol kg$^{-1}$) concentrations below 1,000 m at all stations. UCDW, LCDW, and AABW had mean pCd concentrations of $1.34 \pm 1.08$, $0.35 \pm 0.21$, and $0.39 \pm 0.12$ pmol kg$^{-1}$, respectively. In the SML, mean pCd was highest near the APF (St. 48°S; $18.83 \pm 3.76$ pmol kg$^{-1}$), and generally decreased with distance north and south of this maximal zone. To the north, mean pCd in the SML of the

STZ was $6.18 \pm 3.49$ pmol kg$^{-1}$ and $9.44 \pm 2.90$ pmol kg$^{-1}$ in the SAZ. To the south, mean pCd in the SML of the AAZ was $14.72 \pm 4.04$ pmol kg$^{-1}$.

## Particulate Phosphorus

The depth profiles of P (**Figure 2D** and **Supplementary Figure 2**) resembled those of pCd whereby surface waters were characterized by high P concentrations which decreased rapidly in subsurface water masses to extremely low values at depth. Like pCd, highest SML P concentrations ($14.29 \pm 1.78$ nmol kg$^{-1}$) were measured near the APF. North of the APF, concentrations in the SML ranged between $9.77 \pm 3.26$ and $12.84 \pm 4.49$ nmol kg$^{-1}$ in the SAZ and STZ, respectively. South of the APF, concentrations in the SML were between $8.54 \pm 1.19$ and $12.98 \pm 0.53$ nmol kg$^{-1}$. At greater depths, mean P concentrations in UCDW, LCDW, and AABW were $2.47 \pm 1.40$, $1.93 \pm 0.87$, and $1.48 \pm 0.25$ nmol kg$^{-1}$, respectively.

# DISCUSSION

## Comparison With Previously Reported dCd and pCd

### Global Comparison

The data presented here contributes the first measurements of dCd and pCd from the 30°E longitude in the Indian sector of the Southern Ocean and one of only a handful of winter trace metal datasets produced globally, thus closing key spatial and seasonal gaps in the global database. For comparative purposes, a global compilation of dCd and pCd measurements (**Figures 3A–D**) has been created from all available data included in the GEOTRACES Intermediate Data Product (IDP) 2017 (Schlitzer et al., 2018; see figure caption for full list of references). For dCd, our Southern Ocean data fell within the range of concentrations found in other ocean basins. Globally, the majority of vertical ocean profiles of dCd showed surface minima, from $0.04$ pmol kg$^{-1}$ in the Mediterranean Sea (Gerringa et al., unpublished) to $36$ pmol kg$^{-1}$ in the Atlantic Southern Ocean (Baars et al., 2014). Below the surface, most vertical ocean profiles displayed intermediate to deep dCd maxima, from $368$ pmol kg$^{-1}$ in the Mediterranean Sea (Gerringa et al., unpublished) to $1,246$ pmol kg$^{-1}$ in the North Pacific Ocean (Minami et al., 2015). For the full depth profiles, our Southern Ocean dCd data had a range of $9–1,150$ pmol kg$^{-1}$ which compared most similarly to that of the Atlantic Southern Ocean ($36–970$ pmol kg$^{-1}$; Baars et al., 2014) and Indian Ocean ($10–1,100$ pmol kg$^{-1}$; Thi Dieu Vu and Sohrin, 2013). For pCd, our Southern Ocean data fell within the range of concentrations found in other ocean basins and most closely resembled those observed in the South Atlantic, particularly in the upper 1,000 m. Highest pCd concentrations were found in surface waters of all datasets while extremely low values, between $0.02$ pmol kg$^{-1}$ (South Atlantic; Henderson et al., unpublished) and $0.06$ pmol kg$^{-1}$ (this study) characterized deep ocean distributions. For the upper 1,000 m, mean pCd concentrations were between $1.13$ and $5.04$ pmol kg$^{-1}$ in the North Atlantic (Hayes et al., 2018) and Equatorial Pacific (Lee et al., 2018), respectively, while the mean value for our Southern Ocean data was $5.92$ pmol kg$^{-1}$.
A

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

B

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

C

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

D

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

**FIGURE 3** | Box and whisker plots of (A) dCd for full depth profiles, (B) pCd for full depth profiles, (C) dCd for profiles to 1,000 m and (D) pCd for profiles to 1,000 m. Data compiled from all available data from the GEOTRACES intermediate data product (IDP) 2017 (Schlitzer et al., 2018). Plots in bold are data from this study. From bottom to top, plots are structured as follows: minimum, quartile 1, median, quartile 3, maximum. The cross (x) represents the mean value. Note the changes in the vertical scale. For dCd data references are as follows, Indian Southern Ocean data [this study (n = 132 full depth; 103 < 1,000 m)], Atlantic Southern Ocean [transect GIPY05 (n = 177 full depth; 116 < 1,000 m) (Baars et al., 2014)], Pacific Southern Ocean [transect GIPY06 (n = 82 full depth; 223 < 1,000 m) (Butler et al., 2013), transect GIPY02 (n = 0 full depth; 109 < 1,000 m) (Butler et al., 2013), transect GPpr01 (n = 0 full depth; 34 < 1,000 m) (Hassler and Ellwood, unpublished), transect GPpr02 (n = 77 full depth; 168 < 1,000 m (Ellwood, 2008)], North Atlantic Ocean [transect GA02 (n = 407 full depth; 227 < 1,000 m) (Middag et al., 2018), transect GA03 (n = 588 full depth; 303 < 1,000 m) (Hayes et al., 2018)], Tropical Atlantic Ocean [transect GA02 (n = 479 full depth; 259 < 1,000 m) (Middag et al., 2018), transect GA11 (n = 57 full depth; 27 < 1,000 m) (Xie et al., 2018)], South Atlantic Ocean [transect GA02 (n = 482 full depth; 261 < 1,000 m) (Xie et al., 2015; Middag et al., 2018), transect GA10 (n = 340 full depth; 159 < 1,000 m (Achterberg et al., unpublished), transect GIPY04 (n = 86 full depth; 103 < 1,000 m) (Boye et al., 2012)], Indian Ocean [transect GI04 (n = 167 full depth; 94 < 1,000 m) (Thi Dieu Vu and Sohrin, 2013)], North Pacific Ocean [transect GP02 (n = 256 full depth; 137 < 1,000 m) (Minami et al., 2015), transect GP18 (n = 94 full depth; 49 < 1,000 m) (Minami et al., 2015), transect GPpr07 (n = 25 full depth; 68 < 1,000 m) (Janssen et al., 2017)], Tropical Pacific Ocean [transect GP16 (n = 667 full depth; 350 < 1,000 m) (John et al., 2017)], South Pacific Ocean [transect GP13 (n = 342 full depth; 544 < 1,000 m) (Bowie et al., unpublished)] and Mediterranean Sea [transect GA04 (n = 1376 full depth; 1,057 < 1,000 m) (Gerringa et al., unpublished)]. For pCd, references are as follows, Indian Southern Ocean data [this study (n = 122 full depth; 94 < 1,000 m)], North Atlantic Ocean [transect GA03 (n = 0 full depth; 34 < 1,000 m) (Hayes et al., 2018)], Tropical Atlantic Ocean [transect GA06 (n = 103 full depth; 82 < 1,000 m) (Achterberg, unpublished)], South Atlantic Ocean [transect GA10 (n = 231 full depth; 134 < 1,000 m) (Henderson et al., unpublished)] and Tropical Pacific Ocean [transect GP16 (n = 703 full depth; 397 < 1,000 m) (Lee et al., 2018)].

## Southern Ocean dCd Comparison

Considering only surface (~25 m) dCd measurements allowed for some interesting observations to be made regarding inter-basin and seasonal differences within the Southern Ocean (Figure 4). At the lower latitudes (north of ~45°S), surface

dCd was low across all datasets (Figure 4A) and our winter concentrations (9–34 pmol kg⁻¹) were most similar to those measured in the Pacific sector of the Southern Ocean (17–23 pmol kg⁻¹). At the higher latitudes (south of ~45°S), surface dCd increased rapidly across all datasets.
A

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

B

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

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

**FIGURE 4** | Southern Ocean surface (~25 m) comparisons of (A) dCd against latitude and (B) dCd:PO₄ spot ratios against latitude from the various sectors of the Southern Ocean and during both winter and summer seasons. References as follows: Atlantic sector (Boye et al., 2012; Baars et al., 2014); Pacific/Indian sectors (Janssen et al., 2020); Pacific sector (summer: Butler et al., 2013; winter: Ellwood, 2008); Indian sector (This study).

While the sharpest dCd increase was observed in the winter Indian sector of the Southern Ocean, dCd:PO₄ spot ratios calculated for all datasets were more comparable (Figure 4B), particularly for the Pacific sector (winter and summer) and Indian sector (winter) of the Southern Ocean, suggesting dCd increases were proportional to increases in PO₄. Between 55°S and 60°S there was a decoupling between winter (only Indian sector dataset considered as the Pacific sector measurements stop at 52°S) and summer dCd:PO₄ whereby winter ratios were considerably higher. These seasonal differences are likely caused by two primary mechanisms. First, intense summer productivity at the higher latitudes is driven by diatoms with high intracellular pCd:P (Baars et al., 2014) causing reduced surface dCd and dCd:PO₄ compared to winter values. Second, deep winter mixing provides access to subsurface water masses with enriched dCd:PO₄ signatures (from diatom remineralisation and regional upwelling of UCDW with high dCd:PO₄; see section "Physical Resupply of dCd Through Upwelling and Deep Winter Mixing") thereby preferentially enriching surface waters with additional dCd. While these comparisons do present interesting insights into the seasonal cycling of dCd, interpretations should be made

with caution considering that the data are not direct seasonal
crossover comparisons.

Distribution of dCd and pCd in the Upper
Water Column

Relation With PO₄ and P

Seawater samples were collected in early to mid-July and therefore represent early winter conditions. Owing to increased surface mixing and decreased light exposure (Dong et al., 2008; Swart et al., 2015), the winter months are generally considered biologically dormant in comparison to the spring and summer months. However, estimates of relative productivity (inferred from measurements of Chl-*a*; Supplementary Figure 3), were comparable to spring measurements and were at the lower end of summer measurements (Weir et al., 2020). For example, north of the APF, our early winter (July) Chl-*a* concentrations averaged over the SML (0.31–0.54 µg kg⁻¹) were similar to Chl-*a* concentrations measured to the east during spring (0.34 µg kg⁻¹; Becquevort et al., 2000) and summer (0.39–1.23 µg kg⁻¹; Rembauville et al., 2016). In the AAZ, our SML Chl-*a* (0.20–0.30 µg kg⁻¹) was within the range of previous spring
concentrations (0.10–0.39 µg kg⁻¹; Becquevort et al., 2000) and below previous summer concentrations (0.49–0.98 µg kg⁻¹; Rembauville et al., 2016). This suggests that the winter euphotic layer was not biologically dormant (Weir et al., 2020) and that biological uptake by phytoplankton could potentially play an important role in controlling dCd and pCd distributions during winter.

The distribution of Cd and P were closely correlated across the transect, particularly in the upper water column (above the depth of the O₂ minimum) where a strong linear relationship between dCd vs. PO₄ [dCd (pmol kg⁻¹) = 431 ± 44 PO₄ (µmol kg⁻¹) - 198; r² = 0.89; n = 93] and pCd vs. P [pCd (pmol kg⁻¹) = 0.65 ± 0.11 P (nmol kg⁻¹) + 0.40; r² = 0.86; n = 85] was observed. Assuming that all P is of biogenic origin, and considering the close coupling of pCd and P, it is reasonable to assume that pCd is majorly biogenic in composition. Low dCd and PO₄ concentrations near the ocean surface and enrichments in subsurface waters are widely accepted to result from biological uptake by phytoplankton in the euphotic zone and subsequent remineralisation at depth (Abouchami et al., 2014; Baars et al., 2014; Sieber et al., 2019; Janssen et al., 2020). Biological uptake and assimilation of dCd and PO₄ drives the production of organic pCd and P as observed by the high surface pCd and P concentrations at each station (Figure 2 and Supplementary Figures 1, 2). Furthermore, vertical profiles of pCd co-varied with those of Chl-*a* in the upper 150 m and pCd maxima were either at or very close to the depth of the Chl-*a* maxima (Supplementary Figure 3). In the subsurface, decreasing pCd and P concentrations correspond to generally increasing dCd and PO₄ concentrations partly reflecting the aerobic remineralisation of sinking phytoplankton cells by the microbial community. The extremely low (<1 pmol kg⁻¹) concentrations of pCd, and relatively high concentrations of dCd (>500 pmol kg⁻¹), below 1,000 m suggests that variations in Cd distributions at depth were primarily controlled by the advection of water masses (with preformed Cd signatures) into the sampling region and were therefore not related to local particle transformations.

## Cd Uptake Influenced by the Availability of Dissolved Cd and Other Metals (Fe, Zn, and Mn)

The availability of Cd and the relative availability of several essential trace metals (e.g., Fe, Zn, and Mn) has been shown to impact Cd uptake and therefore directly influence phytoplankton pCd:P stoichiometries (Price and Morel, 1990; Lee and Morel, 1995; Sunda and Huntsman, 2000; Cullen and Sherrell, 2005; Middag et al., 2018). In an effort to better understand the link between Cd uptake and the availabilities of Cd, Fe, Zn and Mn, we compare surface (~25 m) pCd:P with dCd and dissolved ratios of Cd to Fe (dCd:dFe), Cd to Zn (dCd:dZn) and Cd to Mn (dCd:dMn; Figure 5; see Supplementary Tables 1, 2 for validation of dFe, dZn, and dMn). Ratios of pCd:P were low in the STZ and SAZ (0.37–0.52 mmol mol⁻¹), increased in the PFZ (0.86 mmol mol⁻¹) and were highest in the AAZ (0.91–1.07 mmol mol⁻¹). In the STZ, dCd was 9 pmol kg⁻¹ and increased roughly 80-fold to 768 pmol kg⁻¹ in the AAZ. All dCd:metal ratios were lowest in the STZ and increased at various rates through the SAZ and PFZ. In the PFZ, dCd:dZn was highest

(0.38 mol mol⁻¹) and decreased in the AAZ (0.18–0.23 mol mol⁻¹) while dCd:dFe (13.78–15.36 mol mol⁻¹) and dCd:dMn (1.62–2.29 mol mol⁻¹) were highest in the AAZ.

The uptake of Cd by phytoplankton is directly related to the concentration of free (bioavailable) dCd and inversely related the concentration of free dFe, dZn, and dMn. Since Mn forms only weak organic complexes, the bioavailable Mn concentration is assumed to be approximated by dMn (Sunda and Huntsman, 2000). For Cd and Zn, organic complexation may potentially impact biological uptake in the STZ and SAZ where dCd and dZn are inherently low, however, in the PFZ and AAZ, bioavailable dCd and dZn are high owing to increased concentrations which saturate metal binding ligands (Baars and Croot, 2011; Baars et al., 2014). For Fe, organic complexation dominates Fe speciation in the Southern Ocean (Croot et al., 2004) and, coupled with the low surface dFe (0.03–0.12 nmol kg⁻¹) measured throughout the transect, suggests that dFe was potentially at or near limiting concentrations. It follows that in regions with higher dCd and dCd:metal ratios, phytoplankton may assimilate more Cd relative to P providing at least one potential driver behind the elevated pCd:P in the PFZ and AAZ. While our data suggest dCd availability is the primary factor determining Cd uptake, as suggested previously (Middag et al., 2018), assessing the relative importance of the various dCd:metal ratios in influencing pCd:P is more complex. Lowest dFe was measured in the AAZ (0.03–0.05 nmol kg⁻¹) potentially limiting phytoplankton growth in this region resulting in increased cellular pCd:P via growth biodilution or increased uptake of Cd (relative to P) through non-specific divalent metal transporters, a mechanism which is intensified under low Fe (Lane et al., 2009). In the PFZ, pCd:P was relatively high (comparable to the lower end of pCd:P in the AAZ) yet dCd:dFe was more than 4 fold lower than in the AAZ. Based on dCd:dFe alone, we would expect pCd:P in the PFZ to be much lower in comparison to the AAZ. However, dCd:dZn in the PFZ was almost two-fold higher than the AAZ and dCd:dMn was high (although not as high as values calculated for the AAZ) possibly indicating increased competition success for dCd over dZn (and possibly dMn) or passive uptake of Cd to prevent toxicity (Horner et al., 2013) thereby increasing pCd:P ratios. In the northern AAZ, pCd:P was higher, and dCd lower, at 50°S compared to 56°S, which contradicts the primary role for dCd availability in governing pCd:P. At both stations, dCd:dFe was equally high, however, dFe at 50°S (0.03 nmol kg⁻¹; the lowest surface value measured along the transect) was lower than at 56°S (0.05 nmol kg⁻¹) potentially limiting phytoplankton growth (and increasing pCd:P) at 50°S to a greater degree than 56°S. Further south, dCd, pCd:P, and dCd:dFe decreased between 56°S and 58°S and dCd:dZn and dCd:dMn increased suggesting that, of the Cd:metal antagonisms considered as secondary controls (after the primary control by dCd), the dCd:dFe antagonism had the largest effect on pCd:P variations in this region.

## Relating pCd:P to Phytoplankton Stoichiometry

The literature regarding plankton metal stoichiometries in natural phytoplankton assemblages is scarce
<figure><img src="image_12.png" /><figcaption>**FIGURE 5** | Plot of dCd and pCd:P, dCd:dFe; dCd:dZn, and dCd:dMn spot ratios for surface waters (~25 m) across the transect. Note: where necessary units were converted to plot on a single vertical axis for direct comparison.</figcaption></figure>

(Twining and Baines, 2013). Nevertheless, phytoplankton cellular Cd quotas are reported to range by a factor of 100 (0.01–1 mmol mol⁻¹ P) compared to a factor of 20 for other bioactive metals such as Fe, Mn, Zn, and Cu (Ho et al., 2003; Twining and Baines, 2013). We compare our surface water (~25 m) pCd:P ratios with plankton derived Cd:P ratios to assess whether or not the surface water particulate signature is a good representation of the observed phytoplankton community based on observed phytoplankton community data from this transect (Weir et al., 2020, **Supplementary Figure 4**). We acknowledge the limitations of such a comparison which include differences in oceanographic setting (e.g., open ocean vs. coastal upwelling) and inherent environmental conditions (e.g., oligotrophic vs. nutrient replete conditions, differences in resident phytoplankton community assemblage and grazing pressures by local zooplankton) which may influence Cd uptake systematics. Regardless, we find that pCd:P do provide a reasonable estimation of the resident phytoplankton community. To illustrate, a flagellate-dominated assemblage from the North Atlantic had a pCd:P ratio of 0.51 ± 0.09 mmol mol⁻¹ (Kuss and Kremling, 1999), almost exactly matching our pCd:P ratio from the flagellate-dominated SAZ stations (0.52 mmol mol⁻¹ for both SAZ stations). However, these ratios were higher than the STZ (0.37 mmol mol⁻¹) station and lower than the PFZ station (0.86 mmol mol⁻¹) despite significant contributions from flagellates to the phytoplankton community at both stations. There was, however, a much larger contribution of diatoms relative to flagellates at the PFZ station compared to the STZ providing a potential driver of variations in pCd:P whereby diatom growth rates have been shown to directly influence pCd:P ratios while flagellates show no functional response in pCd:P (Finkel et al., 2007). In the PFZ, changes in diatom growth rate, potentially due to the availability of other trace metals [see section “Cd Uptake Influenced by the Availability of Dissolved Cd and Other Metals (Fe, Zn, and Mn)”], would therefore have

a greater bearing on the total phytoplankton community pCd:P ratio compared to the STZ. In the AAZ, the high pCd:P ratios (0.91–1.07 mmol mol⁻¹) compared most similarly to a ratio from a diatom-dominated assemblage from the high latitude Southern Ocean (1.29 ± 0.07 mmol mol⁻¹; Cullen et al., 2003) in agreement with our phytoplankton community assemblage data. The higher pCd:P ratios in the Cullen study could be the result of lower summer dFe and dMn concentrations (0.03 and 0.08 mol kg⁻¹, respectively) which may decrease diatom growth rates (and therefore increase pCd:P) to a greater extent compared to winter.

## Remineralisation Length Scales of Cd and P

In order to investigate remineralisation rates of pCd and P at each station, we plot pCd and P concentrations individually (as a percentage of their respective maximum concentrations at each station) with depth as well as corresponding pCd:P ratios with depth (**Figures 6A–D**). We show one representative profile for each Southern Ocean Zone (see **Supplementary Figure 5** for all other stations and **Supplementary Figure 6** for full depth profiles). Maximum pCd and P concentrations generally occurred at the same depth and were most commonly located between the euphotic depth and the base of the mixed layer (with an exception in the PFZ where maxima were co-located above the euphotic depth). We observe different trends in the vertical pCd:P profiles across the transect. In the AAZ (**Figure 6A**), pCd:P increased to a subsurface maximum (1.2 mmol mol⁻¹) at 100 m, consistent with a preferential remineralisation of P over this depth range. Below the pCd:P maximum, ratios decreased rapidly and the relative concentrations of pCd and P suggest a preferential remineralisation of Cd at depth. These observations are consistent with previous findings in the AAZ of the Southern Ocean, and in high-nutrient, low-chlorophyll areas throughout the global ocean, where it was hypothesized that a labile P fraction (associated with ATP and DNA in the cell) is preferentially remineralised over pCd near the surface resulting in an increase
<figure><img src="image_13.png" /></figure>

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

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

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

**FIGURE 6** | Particulate cadmium (black bars) and P (gray bars) normalized to profile maximum values (bottom axis) for the upper 1,000 m in (A) the AAZ, (B) the PFZ, (C) the SAZ, and (D) the STZ. The red dots are the absolute pCd:P ratios at each depth (top axis). Dashed horizontal lines represents the bottom of the euphotic zone at each station. Solid horizontal line represents the bottom of the MLD at each station. No euphotic zone data for 56°S (AAZ) due to PAR sensor failure.

in pCd:P (Bourne et al., 2018). Below the pCd:P maximum, pCd is preferentially remineralised over a second, more labile, P fraction (associated with phospholipid membranes) causing a decrease in pCd:P through the mesopelagic zone (Bourne et al., 2018).

A major difference between vertical profiles of pCd:P in the AAZ and those to the north was the rate at which pCd:P decreased with depth below the subsurface maximum. For example, pCd:P in the AAZ decreased by a factor of 7.5 from a subsurface maximum of 1.2 mmol mol⁻¹ to 0.16 mmol mol⁻¹ at 1,000 m. This factor was 3.6 in the PFZ (Figure 6B), 1.7 in the SAZ (Figure 6C) and 0.9 in the STZ (Figure 6D; there was no distinct pCd:P maximum here although a local subsurface maximum at 100 m was evident). Since all vertical pCd:P profiles ultimately converge on relatively consistent deep water (> 1,500 m) ratios, stations with higher ratios in surface waters (i.e., PFZ and AAZ) must decrease more rapidly, relative to depth, compared to stations in the SAZ and STZ. Vertical and latitudinal

differences in pCd:P profile shapes may simply be related to intracellular Cd concentrations, which vary as a function of the availability of Cd and other metals [see section “Cd Uptake Influenced by the Availability of Dissolved Cd and Other Metals (Fe, Zn, and Mn)”]. In addition, variations in pCd:P profiles are potentially driven by the proportion of diatoms in the resident phytoplankton community. In the AAZ (diatom dominated), preferential remineralisation of pCd below the pCd:P maximum will occur to a greater extent because of the elevated cellular pCd:P in diatoms. In the PFZ (diatom and flagellate dominated), SAZ (flagellate dominated), and STZ (flagellate dominated), the increasing contributions of flagellates, likely with lower cellular pCd:P compared to diatoms, provides a potential mechanism behind the slower rates of pCd remineralisation below the pCd:P maximum. Data pertaining to taxonomic differences in pCd:P are scarce and have not been equivocally demonstrated in the field. For other metals, however, notably for Fe, Mn, and Zn,
metal:P ratios of diatoms are higher than flagellates from the Southern Ocean (Twining and Baines, 2013) providing at least some support for a taxonomic effect on pCd:P. In addition to the near-surface pCd:P maximum, there is a secondary, deeper, maximum visible at the northern stations. This maximum deepens from 200 m at 48°S to 400 m at 45°S and 800 m at 41°S (Figures 6B–D) and appears to be spatially associated with

isopycnals corresponding to AAIW (Figure 7A). This feature may represent the advection of particles with relatively high pCd:P from the upper water column of the PFZ along the northward subducting flow path of AAIW. However, not all profiles fit the described vertical behavior of pCd:P. For example, in the SAZ (45°S), pCd:P decreased in the upper 150 m, a direct contrast to the other stations. Here, and likely other stations

<figure><img src="image_17.png" /><figcaption>**FIGURE 7** | Distribution plots of (A) pCd:P spot ratios for the upper 1,500 m and (B) dCd:PO₄ spot ratios separated into upper 500 m (upper panel) and below 500 m (bottom panel). Frontal position and corresponding zones are shown. In (A), the solid white line is the median isopycnal defining Antarctic Intermediate Water (AAIW) and the dashed white lines are the upper and lower isopycnal limits for AAIW. In (B), dotted lines represent the average range of the summer mixed layer depth [data from WOCE Global Data Resource (2018)]. Solid line represents the calculated winter mixed layer depth during this study. Figure constructed using Ocean Data View (ODV; Schlitzer, 2020).</figcaption></figure>
where deviations from general trends are evident, pCd:P may include signatures of dead organic matter, bacteria and possibly small zooplankton which may affect pCd:P ratios. Ultimately, the near-surface pCd:P maxima were consistent with increasing dCd:PO₄ ratios with depth at each station and in agreement with previous suggestions of pCd and P remineralisation length scales impacting ratios of dCd to PO₄ in surface and subsurface waters (Quay et al., 2015; Bourne et al., 2018).

## Physical Resupply of dCd Through Upwelling and Deep Winter Mixing

Considering biological uptake as the only control of surface dCd (and PO₄), and the strong preferential uptake of dCd (relative to PO₄) in the PFZ and AAZ compared to the north, we would expect surface dCd to remain low throughout the southern section of the transect. However, the mean SML dCd concentration in the PFZ and AAZ was between 318 and 714 pmol kg⁻¹ and significantly higher compared to stations to the north (Figure 2A). The strong surface dCd gradient was likely the result of the upwelling of underlying water masses, primarily UCDW, which supply dCd enriched waters to the surface layer. Upwelling of UCDW occurs in the AAZ, where UCDW was sampled between 300–600 m. This water is then transported northwards into the upper water column of the lower latitudes (Orsi et al., 1995). Importantly, UCDW carried a high mean dCd:PO₄ spot ratio of 367 ± 70 µmol mol⁻¹ (Figure 7B), similar to the previous summer spot ratios (365 µmol mol⁻¹; Baars et al., 2014). The dCd:PO₄ spot ratios in UCDW increased southwards and reached a maximum of ~500 µmol mol⁻¹ below the SML in the AAZ. The source of the ‘excess’ dCd relative to PO₄ was predominantly the remineralisation of particles with high pCd:P (Figure 7A) which was set during biological uptake within the SML, as observed elsewhere in the Southern Ocean (Abouchami et al., 2014; Baars et al., 2014; Middag et al., 2018; Sieber et al., 2019).

In addition to the strong latitudinal trends in dCd:PO₄ spot ratios, we observed a significant increase with depth (Figure 7B). The high dCd:PO₄ spot ratios found beneath the summer mixed layer support the notion of a physical resupply of excess dCd (i.e., preferential resupply of dCd relative to PO₄) to the SML through deep winter mixing. While our sampling took place before the period of maximum nutrient recharge, between August and November (Johnson et al., 2017), the early winter MLDs in the SAZ were similar to the maximum values (~200 m) expected in west Indian sector (Sallée et al., 2010) suggesting deep winter mixing and nutrient recharge was ongoing. Deep winter mixing is an important nutrient resupply mechanism in the Southern Ocean as shown previously for Fe (Tagliabue et al., 2014; Mtshali et al., 2019) and other trace metals (Cloete et al., 2019). Here, deeper mixing during the winter months allows access to subsurface waters with dCd enrichments relative to PO₄ (i.e., higher dCd:PO₄ spot ratios) found below the depleted summer mixed layer, potentially resupplying dissolved nutrients which may aid and sustain local phytoplankton production through the spring and summer months. For example, in the STZ and SAZ, dCd:PO₄ increased by two-threefold from 49 to

85 µmol mol⁻¹ at the depth of the average summer mixed layer in this region (60–80 m; WOCE Global Data Resource, 2018), to values of 162–178 µmol mol⁻¹ at the base of the winter mixed layer (~200 m; Figure 7B). Winter deep mixing may be of particular importance in the STZ and SAZ where summer surface dCd can be depleted to <10 pmol kg⁻¹ (Boye et al., 2012).

## dCd in Intermediate and Deep Waters and the Importance of Water Mass Mixing

While dCd and PO₄ are well correlated, our data show a change in the steepness of the regression between dCd and PO₄ (dCd/PO₄) at PO₄ > 1.0 µmol kg⁻¹ (Figure 8), a feature commonly referred to as the ‘kink,’ and which is evident in the global dCd/dPO₄ relationship (Middag et al., 2018), derived from the GEOTRACES 2017 Intermediate Data Product (Schlitzer et al., 2018). The samples with low dCd and PO₄ represent surface waters from the STZ and SAZ and plot on a relatively shallow linear slope (61 ± 9 µmol mol⁻¹) with a slightly negative intercept indicating a slight preferential uptake of Cd relative to PO₄. Our winter slope value is higher compared to summer STZ and SAZ surface waters (31 ± 6 µmol mol⁻¹; Boye et al., 2012; Baars et al., 2014) likely reflecting preferential dCd recharge (relative to PO₄) as a result of deep winter mixing as discussed previously (see section “Physical Resupply of dCd Through Upwelling and Deep Winter Mixing”). The higher dCd/PO₄ signatures of intermediate and deep water masses are collectively represented by the steep linear slope (Frew and Hunter, 1992). We calculate a dCd/PO₄ slope value of 453 ± 68 µmol mol⁻¹ (n = 73) for these waters which is similar to slopes observed in previous Southern Ocean winter measurements (444 ± 11 µmol mol⁻¹; Ellwood, 2008), and lower than slope values in the Atlantic sector of the Southern Ocean (590 ± 12 µmol mol⁻¹; Boye et al., 2012; Baars et al., 2014). The high slope values and negative intercepts in our winter data suggest a combination of preferential uptake (and subsurface regeneration) of dCd relative to PO₄ but also strong upwelling of dCd enriched water masses in the AAZ (Baars et al., 2014). Furthermore, we identify that in fact multiple linear regressions describe data from the steep slope (Figure 8) as shown previously in the Atlantic (Middag et al., 2018). Each slope represents a different water mass and it is assumed that mixing between water masses with different origins and hence, different preformed nutrient components, results in a new regression where the slope value falls between the slope values of the water mass end-members being mixed. For example, NADW has a dCd/dPO₄ slope value of 250 µmol mol⁻¹ at its source in the high latitude North Atlantic (Bruland and Franks, 1983) compared to a slope value of 311 ± 13 µmol mol⁻¹ (n = 8) measured here. As NADW flows southward across the equator, the increase in slope value arises from mixing with northward flowing AAIW (341 ± 31 µmol mol⁻¹; n = 30) at intermediate depths and at deeper depths with CDW (344 ± 77 µmol mol⁻¹; n = 38) and AABW (2934 ± 44 µmol mol⁻¹; n = 4). These Antarctic derived water masses have higher dCd/dPO₄ slope values such that the slope for NADW observed here lies between its source waters and the waters it interacted with. Our observations clearly
<figure><img src="image_18.png" /><figcaption>**FIGURE 8** | Plot of dCd/dPO₄ ratios, defined as the slope of regression between dCd and PO₄, for various water masses.</figcaption></figure>

show the strong influence of water mass end-member mixing on deep dCd distributions in this region and is consistent with data from the Atlantic and Pacific sectors of the Southern Ocean (Baars et al., 2014, Sieber et al., 2019; Janssen et al., 2020).

## CONCLUSION

This study presents the first winter measurements of dissolved cadmium (dCd) and particulate cadmium (pCd) along the 30°E longitude in the Indian Sector of the Southern Ocean and provides a unique opportunity to investigate Cd cycling during a period characterized by contrasting upper water column dynamics compared to summer. Our data suggest that biological uptake of trace metals is an important factor to consider during the winter season despite unfavorable growth conditions for phytoplankton. At all stations, dCd was lowest in the SML and increased with depth. Conversely, pCd was highest in the SML, but not necessarily in the shallowest euphotic zone sample, and decreased rapidly with depth. The vertical behaviors of both dCd and pCd are consistent with biological uptake in surface waters and remineralisation in the subsurface. Concentrations of dCd (relative to PO₄) near the surface were higher than comparable summer measurements, most evident in the high latitude AAZ, likely reflecting a combination of weaker biological uptake and preferential dCd recharge through deep winter mixing. We observed distinct regional trends in Cd uptake and remineralisation reflecting a combination of changes in phytoplankton community composition and the availability of both dCd other trace metals such as Fe, Zn and Mn. Our data also provided evidence for differential remineralisation rates of pCd vs. P with depth at each station with implications for ratios of dCd to PO₄ in the upper water column. The combined effect of diatom uptake and regeneration characteristics and the upwelling of

water masses in the PFZ and AAZ resulted in a strong southward increase in surface dCd. In the intermediate and deep waters, where dCd dominated the total Cd inventory, there were distinct water mass-dependent relationships between linear regressions of dCd versus dPO₄ which suggested that dCd distributions are largely the result of mixing between northward traveling waters of Antarctic origin with high dCd/dPO₄ signatures and southward traveling NADW with an inherently low dCd/dPO₄ signature.

## DATA AVAILABILITY STATEMENT

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

## AUTHOR CONTRIBUTIONS

RC, JL, NH, and TM collected field samples. RC and JL analyzed for dissolved Cd. NH and HP analyzed for particulate Cd. RC wrote the manuscript with contributions from SF and AR. AR provided funds for the study. All authors commented on previous drafts of the manuscript, contributed to the final submission draft and approved the submitted version.

## FUNDING

RC was supported through the National Research Foundation (NRF) Innovation Ph.D. studentship. This research was supported by NRF grants (UID# 93069, 105826, and 110715) to AR. SF acknowledges funding from NRF (UID 111210 and 110731).
ACKNOWLEDGMENTS

We would like to thank the South African National Antarctic
Programme (SANAP) as well as Captain Knowledge Bengu and
the crew of the R/V SA Agulhas II for their professionalism
and support during the winter 2017 voyage. We acknowledge
Chief Scientist Marcello Vichi and all the participants involved
in the expedition. We are grateful to the "Iron" team for
their help in collecting trace-clean seawater samples. We
thank the two reviewers for their efforts which resulted in
an improved manuscript. We would like to thank Gideon

Henderson, Eric Achterberg, Angela Milne, Maeve Lohan,
Andrew Bowie, Yoshiki Sohrin, and Michael Ellwood for
allowing permission to use data that are not published outside
of the IDP2017.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found
online at: https://www.frontiersin.org/articles/10.3389/fmars.
2021.656321/full#supplementary-material

REFERENCES

Abouchami, W., Galer, S. J. G., De Baar, H. J. W., Middag, R., Vance, D., Zhao, Y., et al. (2014). Biogeochemical cycling of cadmium isotopes in the Southern Ocean along the Zero Meridian. *Geochim. Cosmochim. Acta* **127**, 348–367.
doi: 10.1016/j.gca.2013.10.022

Anilkumar, N., Luis, A. J., Somayajulu, Y. K., Ramesh Babu, V., Dash, M. K., Pednekar, S. M., et al. (2006). Fronts, water masses and heat content variability in the Western Indian sector of the Southern Ocean during austral summer 2004. *J. Mar. Syst.* **63**, 20–34. doi: 10.1016/j.jmarsys.2006.04.009

Baars, O., Abouchami, W., Galer, S., Boye, M., and Croot, P. (2014). Dissolved cadmium in the Southern Ocean: Distribution, speciation, and relation to phosphate. *Limnol. Ocean.* **59**, 385–399. doi: 10.4319/lo.2014.59.2.0385

Baars, O., and Croot, P. L. (2011). The speciation of dissolved zinc in the Atlantic sector of the Southern Ocean. *Deep. Res. Part II Top. Stud. Oceanogr.* **58**, 2720–2732. doi: 10.1016/j.dsr2.2011.02.003

Becquevort, S., Menon, P., and Lancelot, C. (2000). Differences of the protozoan biomass and grazing during spring and summer in the Indian sector of the Southern Ocean. *Polar Biol.* **23**, 309–320. doi: 10.1007/s003000050450

Belkin, I. M., and Gordon, A. L. (1996). Southern Ocean fronts from the Greewich meridian to Tasmania. *J. Geophys. Res.* **101**, 3675–3696. doi: 10.1029/95JC02750

Bourne, H. L., Bishop, J. K. B., Lam, P. J., and Ohnemus, D. C. (2018). Global Spatial and Temporal Variation of Cd:P in Euphotic Zone Particulates. *Global Biogeochem. Cycles* **32**, 1123–1141. doi: 10.1029/2017GB005842

Boye, M., Wake, B. D., Lopez Garcia, P., Bown, J., Baker, A. R., and Achterberg, E. P. (2012). Distributions of dissolved trace metals (Cd, Cu, Mn, Pb, Ag) in the southeastern Atlantic and the Southern Ocean. *Biogeosciences* **9**, 3231–3246. doi: 10.5194/bg-9-3231-2012

Boyle, E. A., Sclater, F., and Edmond, J. M. (1976). On the marine geochemistry of cadmium. *Nature* **263**, 42–44. doi: 10.1038/263042a0

Broecker, W. S. (1991). The Great Ocean Conveyor. *Oceanography* **4**, 79–89. doi: 10.5670/oceanog.1991.07https

Bruland, K. W., and Franks, R. P. (1983). Mn, Ni, Cu, Zn and Cd in the Western North Atlantic. *Trace Met. Sea Water*. **9**, 395–396. doi: 10.1007/978-1-4757-6864-0_23

Bruland, K. W., Knauer, G. A., and Martin, J. H. (1978). Cadmium in Northeast Pacific Waters. *Limnol. Oceanogr.* **23**, 618–625. doi: 10.4319/lo.1978.23.4.0618

Butler, E. C. V., O'Sullivan, J. E., Watson, R. J., Bowie, A. R., Remenyi, T. A., and Lannuzel, D. (2013). Trace metals Cd, Co, Cu, Ni, and Zn in waters of the subantarctic and Polar Frontal Zones south of Tasmania during the “SAZ-Sense” project. *Mar. Chem.* **148**, 63–76. doi: 10.1016/j.marchem.2012.1.005

Cloete, R., Loock, J. C., Mtshali, T., Fietz, S., and Roychoudhury, A. N. (2019). Winter and summer distributions of Copper, Zinc and Nickel along the International GEOTRACES Section GIPY05: Insights into deep winter mixing. *Chem. Geol.* **511**, 342–357. doi: 10.1016/j.chemgeo.2018.10.023

Conway, T. M., and John, S. G. (2015). Biogeochemical cycling of cadmium isotopes along a high-resolution section through the North Atlantic Ocean. *Geochim. Cosmochim. Acta* **148**, 269–283. doi: 10.1016/j.gca.2014.09.032

Croot, P. L., Andersson, K., Öztürk, M., and Turner, D. R. (2004). The distribution and speciation of iron along 6°E in the Southern Ocean. *Deep. Res. Part II Top. Stud. Oceanogr.* **51**, 2857–2879. doi: 10.1016/j.dsr2.2003.10.012

Cullen, J. T. (2006). On the nonlinear relationship between dissolved cadmium and phosphate in the modern global ocean : Could chronic iron limitation of phytoplankton growth cause the kink ? *Limnol. Oceanogr.* **51**, 1369–1380. doi: 10.4319/lo.2006.51.3.1369

Cullen, J. T., and Sherrell, R. M. (2005). Effects of dissolved carbon dioxide, zinc, and manganese on the cadmium to phosphorus ratio in natural phytoplankton assemblages. *Limnol. Oceanogr.* **50**, 1193–1204. doi: 10.4319/lo.2005.50.4.1193

Cullen, J. T., Chase, Z., Coale, K. H., Fitzwater, S. E., and Sherrell, R. M. (2003). Effect of iron limitation on the cadmium to phosphorous ratio of natural phytoplankton assemblages from the Southern Ocean. *Limnol. Oceanogr.* **48**, 1079–1087. doi: 10.4319/lo.2003.48.3.1079

Cullen, J. T., Lane, T. W., Morel, F. M. M., and Sherrell, R. M. (1999). Modulation of cadmium uptake in phytoplankton by seawater CO2 concentration. *Nature* **402**, 165–167. doi: 10.1038/46007

Cutter, G., Casciotti, K., Croot, P., Geibert, W., Heimbürger, L.-E., Lohan, M., et al. (2017). *Sampling and Sample-handling Protocols for GEOTRACES Cruises. Version 3*. Toulouse: GEOTRACES International Project Office.

de Jong, E., Vichi, M., Mehlmann, C. B., Eayrs, C., De Kock, W., Moldenhauer, M., et al. (2018). *Sea Ice conditions within the Antarctic Marginal Ice Zone in winter 2017, onboard the SA Agulhas II, in: PANGAEA*. Cape Town: University of Cape Town.

Dong, S., Sprintall, J., Gille, S. T., and Talley, L. (2008). Southern ocean mixed-layer depth from Argo float profiles. *J. Geophys. Res. Ocean.* **113**, 1–12. doi: 10.1029/2006JC004051

Duce, R. A., Liss, P. S., Merrill, J. T., Atlas, E. L., Hicks, B. B., Millertl, J. M., et al. (1991). The atmospheric input of trace species to the world ocean. *Global Biogeochem. Cycles* **5**, 193–259. doi: 10.1029/91gb01778

Elderfield, H., and Rickaby, R. E. M. (2000). Oceanic Cd/P ratio and nutrient utilization in the glacial Southern Ocean. *Nature* **405**, 305–310. doi: 10.1038/35012507

Ellwood, M. J. (2004). Zinc and cadmium speciation in subantarctic waters east of New Zealand. *Mar. Chem.* **87**, 37–58. doi: 10.1016/j.marchem.2004.01.005

Ellwood, M. J. (2008). Wintertime trace metal (Zn, Cu, Ni, Cd, Pb and Co) and nutrient distributions in the Subantarctic Zone between 40-52°S; 155-160°E. *Mar. Chem.* **112**, 107–117. doi: 10.1016/j.marchem.2008.07.008

Finkel, Z. V., Quigg, A. S., Chiampi, R. K., Schofield, O. E., and Falkowski, P. G. (2007). Phylogenetic diversity in cadmium: Phosphorus ratio regulation by marine phytoplankton. *Limnol. Oceanogr.* **52**, 1131–1138. doi: 10.4319/lo.2007.52.3.1131

Frew, R. D., and Hunter, K. A. (1992). Influence of Southern Ocean waters on the cadmium-phosphate properties of the global ocean. *Nature* **360**, 144–146. doi: 10.1038/360144a0

Grasshoff, K. (1983). "Automated chemical analysis," in *Methods of Seawater Analysis*, eds K. Grasshoff, M. Ernhardt, and K. Kremling (New York, NY: Verlag Chemie Weinhein), 263–289. doi: 10.1002/9783527613984

Hayes, C. T., Black, E. E., Anderson, R. F., Baskaran, M., Buesseler, K. O., Charette, M. A., et al. (2018). Flux of Particulate Elements in the North Atlantic Ocean Constrained by Multiple Radionuclides. *Glob. Biogeochem. Cycles* **32**, 1738–1758. doi: 10.1029/2018GB005994
Ho, T.-Y., Quigg, A., Finkel, Z. V., Milligan, A. J., Wyman, K., Falkowski, P. G., et al. (2003). The elemental composition of some marine phytoplankton. *J. Phycol.* 39, 1145–1159. doi: 10.1111/j.0022-3646.2003.03-090.x

Hölemann, J. A., Schirmacher, M., and Prange, A. (1999). *Dissolved and Particulate Major and Trace Elements in Newly Formed Ice from the Laptev Sea (Transdrift III, October 1995)*, *Land-Ocean Systems in the Siberian Arctic*. Berlin: Springer-Verlag, doi: 10.1007/978-3-642-60134-7_11

Horner, T. J., Lee, R. B. Y., Henderson, G. M., and Rickaby, R. E. M. (2013). Nonspecific uptake and homeostasis drive the oceanic cadmium cycle. *Proc. Natl. Acad. Sci. U S A.* 110, 2500–2505. doi: 10.1073/pnas.1213857110/
/DCSupplemental.www.pnas.org/cgi//10.1073/pnas.12138
57110

Janssen, D. J., Abouchami, W., Galer, S. J. G., and Cullen, J. T. (2017). Fine-scale spatial and interannual cadmium isotope variability in the subarctic northeast Pacific. *Earth Planet. Sci. Lett.* 472, 241–252. doi: 10.1016/j.epsl.2017.0
4.048

Janssen, D. J., Conway, T. M., John, S. G., Christian, J. R., Kramer, D. I., Pedersen, T. F., et al. (2014). Undocumented water column sink for cadmium in open ocean oxygen-deficient zones. *Proc. Natl. Acad. Sci. U S A.* 111, 6888–6893. doi: 10.1073/pnas.1402388111

Janssen, D. J., Sieber, M., Ellwood, M. J., Conway, T. M., Barrett, P. M., Chen, X., et al. (2020). Trace metal and nutrient dynamics across broad biogeochemical gradients in the Indian and Pacific sectors of the Southern Ocean. *Mar. Chem.* 221:103773. doi: 10.1016/j.marchem.2020.103773

John, S. G., Helgoe, J., and Townsend, E. (2017). Biogeochemical cycling of Zn and Cd and their stable isotopes in the Eastern Tropical South Pacific. *Mar. Chem.* 201, 256–262. doi: 10.1016/j.marchem.2017.06.001

Johnson, K. S., Plant, J. N., Dunne, J. P., Talley, L. D., and Sarmiento, J. L. (2017). Annual nitrate drawdown observed by SOCCOM profiling floats and the relationship to annual net community production. *J. Geophys. Res. Ocean.* 122, 6668–6683. doi: 10.1002/2017JC012839

Kuss, J., and Kremling, K. (1999). Spatial variability of particle associated trace elements in near-surface waters of the North Atlantic (30 °N/60 °W to 60 °N/2 °W), derived by large volume sampling. *Mar. Chem.* 68, 71–86. doi: 10.1016/S0304-4203(99)00066-3

Lam, P. J., and Marchal, O. (2015). Insights into Particle Cycling from Thorium and Particle Data. *Ann. Rev. Mar. Sci.* 7, 159–184. doi: 10.1146/annurev-marine-010814-015623

Lambelet, M., Rehkämper, M., van de Flierdt, T., Xue, Z., Kreissig, K., Coles, B., et al. (2013). Isotopic analysis of Cd in the mixing zone of Siberian rivers with the Arctic Ocean-New constraints on marine Cd cycling and the isotope composition of riverine Cd. *Earth Planet. Sci. Lett.* 361, 64–73. doi: 10.1016/j.epsl.2012.11.034

Lane, E. S., Jang, K., Cullen, J. T., and Maldonado, M. T. (2008). The interaction between inorganic iron and cadmium uptake in the marine diatom Thalassiosira oceanica. *Limnol. Oceanogr.* 53, 1784–1789. doi: 10.4319/lo.2008.53.5.1784

Lane, E. S., Semeniuk, D. M., Strzepek, R. F., Cullen, J. T., and Maldonado, M. T. (2009). Effects of iron limitation on intracellular cadmium of cultured phytoplankton: Implications for surface dissolved cadmium to phosphate ratios. *Mar. Chem.* 115, 155–162. doi: 10.1016/j.marchem.2009.0
7.008

Lane, T. W., and Morel, F. M. M. (2000). A biological function for cadmium in marine diatoms. *Proc. Natl. Acad. Sci. U S A.* 97, 4627–4631. doi: 10.1073/pnas.
090091397

Lane, T. W., Saito, M. A., George, G. N., Prince, R. C., Pickering, I. J., and Morel, F. M. M. (2005). A cadmium enzyme from a marine diatom. *Nature* 435, 42–42. doi: 10.1038/435042a

Lee, J. G., and Morel, F. M. M. (1995). Replacement of zinc by cadmium in marine phytoplankton. *Mar. Ecol. Prog. Ser.* 127, 305–309. doi: 10.3354/meps127

Lee, J. M., Heller, M. I., and Lam, P. J. (2018). Size distribution of particulate trace elements in the U.S. GEOTRACES Eastern Pacific Zonal Transect (GP16). *Mar. Chem.* 201, 108–123. doi: 10.1016/j.marchem.2017.09.006

Middag, R., van Heuven, S. M. A. C., Bruland, K. W., and de Baar, H. J. W. (2018). The relationship between cadmium and phosphate in the Atlantic

Ocean unravelled. *Earth Planet. Sci. Lett.* 492, 79–88. doi: 10.1016/j.epsl.2018.0
3.046

Minami, T., Konagaya, W., Zheng, L., Takano, S., Sasaki, M., Murata, R., et al. (2015). An off-line automated preconcentration system with ethylenediaminetriacetate chelating resin for the determination of trace metals in seawater by high-resolution inductively coupled plasma mass spectrometry. *Anal. Chim. Acta* 854, 183–190. doi: 10.1016/j.aca.2014.1.1016

Mtshali, T. N., Horsten, N. R., Thomalla, S. J., Ryan-Keogh, T. J., Nicholson, S. A., Roychoudhury, A. N., et al. (2019). Seasonal depletion of the dissolved iron reservoirs in the sub-Antarctic zone of the Southern Atlantic Ocean. *Geophys. Res. Lett.* 2019:2018GL081355. doi: 10.1029/2018GL081355

Orsi, A. H., Whitworth, T., and Nowlin, W. D. (1995). On the meridional extent and fronts of the Antarctic Circumpolar Current. *Deep. Res. Part I* 42, 641–673. doi: 10.1016/0967-0637(95)00021-W

Planquette, H., and Sherrell, R. M. (2012). Sampling for particulate trace element determination using water sampling bottles: methodology and comparison to in situ pumps. *Limnol. Oceanogr. Methods* 10, 367–388. doi: 10.4319/lom.2012.10.367

Pollard, R. T., Lucas, M. I., and Read, J. F. (2002). Physical controls on biogeochemical zonation in the Southern Ocean. *Deep. Res. Part II Top. Stud. Oceanogr.* 49, 3289–3305. doi: 10.1016/S0967-0645(02)00084-X

Price, N. M., and Morel, F. M. M. (1990). Cadmium and Cobalt substitution for Zn in a marine diatom. *Nature* 344, 658–660. doi: 10.1038/344658a0

Quay, P., Cullen, J., Landing, W., and Morton, P. (2015). Processes controlling the distributions of Cd and PO4 in the ocean. *Global Biogeochem. Cycles* 29, 830–841. doi: 10.1002/2014GB004998

Rembauville, M., Blain, S., Caparros, J., and Salter, I. (2016). Particulate matter stoichiometry driven by microplankton community structure in summer in the Indian sector of the Southern Ocean. *Limnol. Oceanogr.* 61, 1301–1321. doi: 10.1002/lno.10291

Sallée, J. B., Speer, K., Rintoul, S., and Wijffels, S. (2010). Southern Ocean Thermocline Ventilation. *J. Phys. Oceanogr.* 40, 509–529. doi: 10.1175/2009JPO4291.1

Sarmiento, J. L., Gruber, N., Brzezinski, M. A., and Dunne, J. P. (2004). High-latitude controls of thermocline nutrients and low latitude biological productivity. *Nature* 427, 56–60. doi: 10.1038/nature10605

Schlitzer, R. (2020). *Ocean Data View*. Cambridge: Cambridge University Press, doi: 10.1017/CBO9781107415324.004

Schlitzer, R., Anderson, R. F., Dodas, E. M., Lohan, M., Geibert, W., Tagliabue, A., et al. (2018). The GEOTRACES Intermediate Data Product 2017. *Chem. Geol.* 493, 210–223. doi: 10.1016/j.chemgeo.2018.05.040

Sieber, M., Conway, T. M., Souza, G. F., De, Hassler, C. S., Ellwood, M. J., et al. (2019). High-resolution Cd isotope systematics in multiple zones of the Southern Ocean from the Antarctic Circumnavigation Expedition. *Earth Planet. Sci. Lett.* 527:115799. doi: 10.1016/j.epsl.2019.115799

Sunda, W. G., and Huntsman, S. A. (2000). Effect of Zn, Mn, and Fe on Cd accumulation in phytoplankton: Implications for oceanic Cd cycling. *Limnol. Oceanogr.* 45, 1501–1516. doi: 10.4319/lo.2000.45.7.1501

Swart, S., Thomalla, S. J., and Monteiro, P. M. S. (2015). The seasonal cycle of mixed layer dynamics and phytoplankton biomass in the Sub-Antarctic Zone: A high-resolution glider experiment. *J. Mar. Syst.* 147, 103–115. doi: 10.1016/j.jmarsys.2014.06.002

Tagliabue, A., Sallée, J. B., Bowie, A. R., Lévy, M., Swart, S., and Boyd, P. W. (2014). Surface-water iron supplies in the Southern Ocean sustained by deep winter mixing. *Nat. Geosci.* 7, 314–320. doi: 10.1038/ngeo2101

Thi Dieu Vu, H., and Sohrin, Y. (2013). Diverse stoichiometry of dissolved trace metals in the Indian Ocean. *Sci. Rep.* 3, 1–5. doi: 10.1038/srep01745

Turekian, K. K. (1977). The fate of metals in the oceans. *Geochim. Cosmochim. Acta* 41, 1139–1144. doi: 10.1016/0016-7037(77)90109-0

Twining, B. S., and Baines, S. B. (2013). The trace metal composition of marine phytoplankton. *Ann. Rev. Mar. Sci.* 5, 191–215. doi: 10.1146/annurev-marine-121211-172322
Weir, I., Fawcett, S., Smith, S., Walker, D., Bornman, T., and Fietz, S. (2020). Winter biogenic silica and diatom distributions in the Indian sector of the Southern Ocean. *Deep Sea Res. Part I Oceanogr. Res. Pap.* 166:103421. doi: 10.1016/j.dsr.2020.103421

WOCE Global Data Resource (2018). *NODC Standard Product: World Ocean Circulation Experiment (WOCE) Global Data Resource (GDR), versions 1-3, on CD-ROM and DVD*. Washington, D.C: NOAA.

Xie, R. C., Galer, S. J. G., Abouchami, W., and Frank, M. (2018). Limited impact of eolian and riverine sources on the biogeochemical cycling of Cd in the tropical Atlantic. *Chem. Geol.* 511, 371–379. doi: 10.1016/j.chemgeo.2018.10.018

Xie, R. C., Galer, S. J. G., Abouchami, W., Rijkenberg, M. J. A., De Jong, J., De Baar, H. J. W., et al. (2015). The cadmium-phosphate relationship in the western South Atlantic - The importance of mode and intermediate waters on the global systematics. *Mar. Chem.* 177, 110–123. doi: 10.1016/j.marchem.2015.06.011

Xu, Y., Tang, D., Shaked, Y., and Morel, F. M. M. (2007). Zinc, cadmium, and cobalt interreplacement and relative use efficiencies in the coccolithophore Emiliania huxleyi. *Limnol. Oceanogr.* 52, 2294–2305. doi: 10.2307/4502377

**Conflict of Interest:** 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.

The reviewer GD declared a past co-authorship with one of the authors HP to the handling editor.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2021 Cloete, Loock, van Horsten, Fietz, Mtshali, Planquette and Roychoudhury. 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) and the copyright owner(s) 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.