# Seasonality of downward carbon export in the Pacific Southern Ocean revealed by multi-year robotic observations

Léo Lacour, Joan Llort, Nathan Briggs, Peter G Strutton, Philip W Boyd

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

## To cite this version:

Léo Lacour, Joan Llort, Nathan Briggs, Peter G Strutton, Philip W Boyd. Seasonality of downward carbon export in the Pacific Southern Ocean revealed by multi-year robotic observations. Nature Communications, 2023, 14, pp.1278. 10.1038/s41467-023-36954-7. hal-04018961

HAL Id: hal-04018961

https://hal.science/hal-04018961

Submitted on 8 Mar 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Seasonality of downward carbon export in the Pacific Southern Ocean revealed by multi-year robotic observations

Received: 11 November 2021

Accepted: 27 February 2023

Published online: 08 March 2023

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

Check for updates

Léo Lacour<fnref n="1" />,²✉, Joan Llort ©³, Nathan Briggs⁴, Peter G. Strutton ©<fnref n="1" />,⁵ & Philip W. Boyd ©<fnref n="1" />

At high latitudes, the biological carbon pump, which exports organic matter from the surface ocean to the interior, has been attributed to the gravitational sinking of particulate organic carbon. Conspicuous deficits in ocean carbon budgets challenge this as a sole particle export pathway. Recent model estimates revealed that particle injection pumps have a comparable downward flux of particulate organic carbon to the biological gravitational pump, but with different seasonality. To date, logistical constraints have prevented concomitant and extensive observations of these mechanisms. Here, using year-round robotic observations and recent advances in bio-optical signal analysis, we concurrently investigated the functioning of two particle injection pumps, the mixed layer and eddy subduction pumps, and the gravitational pump in Southern Ocean waters. By comparing three annual cycles in contrasting physical and biogeochemical environments, we show how physical forcing, phytoplankton phenology and particle characteristics influence the magnitude and seasonality of these export pathways, with implications for carbon sequestration efficiency over the annual cycle.

The biological carbon pump (BCP) is considered a major contributor to the Southern Ocean (SO) carbon sink. It removes 2.7 (±0.6) PgC from the euphotic zone (upper -100 m) annually, representing ~30% of the global BCP<fnref n="1" />,². Since the BCP concept was initially proposed³, it has been widely assumed that it was mediated primarily by sinking particulate organic carbon (POC), now termed the Biological Gravitational Pump (BGP)<fnref n="1" />,⁵. However, present estimates of biological carbon demand in the mesopelagic (-100–1000 m depth) exceed carbon inputs attributed to sinking particles into this stratum by two- to threefold⁶,⁷. Such discrepancies highlight the need to reassess the pathways that contribute to downward carbon export.

Recently, other particle injection mechanisms (PIPs) have been invoked as additional pathways that help to balance the mesopelagic carbon budget⁴. Studies⁸⁻¹¹ have provided evidence that downward

export of organic matter also occurs through localized (1–10 km) eddy-driven subduction of POC. This process, known as the Eddy Subduction Pump (ESP), leads to episodic injection of POC-rich waters below the mixed layer. Similarly, the (sub)seasonal variability of the mixing layer depth leads to detrainment of organic matter to the dark mesopelagic when the mixing layer shoals and conversely entrainment to the surface layer when the mixing layer deepens. This physical mechanism is called the mixed layer pump (MLP)¹²⁻¹⁴. At basin scale, the ESP and MLP export is equivalent to ~25% of the BGP export⁴. However, at the (sub)mesoscale (<100 km) these mechanisms can be equal to, or even greater than the BGP¹⁵,¹⁶. This spatial mismatch illustrates the difficulty of carrying out comprehensive inter-comparison of these pumps in the field. Other PIPs such as the mesopelagic migrant pump¹⁷, seasonal lipid pump¹⁸ and large-scale

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

<footnote n="1">¹Institute for Marine and Antarctic Studies, University of Tasmania, Hobart, Australia. ²Sorbonne Université, CNRS, Laboratoire d'Océanographie de Villefranche, LOV, Villefranche-sur-Mer, France. ³Barcelona Supercomputing Center, Earth Sciences Dept., Barcelona, Spain. ⁴National Oceanography Centre, Southampton, UK. ⁵Australian Research Council Centre of Excellence for Climate Extremes, University of Tasmania, Hobart, Australia.</footnote>
subduction¹⁹ pump are also important contributors to the BGP, although assessing their relative contribution is even more challenging⁴ and beyond this study.

As opposed to the BGP, the physically driven PIPs can potentially transport all types of particles to depth, from small non-sinking particles, including healthy phytoplankton cells, to large fast-sinking aggregates such as fecal pellets⁴,⁸,¹³. The composition and size of exported particles will influence their fate in the mesopelagic with respect to depth of remineralisation, sequestration time scale, and consumption by midwater biota⁴,¹⁵,²⁰. To determine the relative importance of each of these pathways for carbon sequestration, it is essential to study these mechanisms in the context of the seasonality in characteristics of the upper ocean particle assemblage⁴. Although previous studies¹¹,²¹⁻²³ revealed that both the PIPs and BGP vary seasonally, concomitant and extensive in situ observations of these pathways, including upper ocean particle characteristics, are extremely rare. Therefore, little is known about their relative importance over complete annual cycles. Better characterisation of these mechanisms and their seasonality will fundamentally advance our understanding of ocean biological carbon export and sequestration, and help to close regional ocean carbon budgets.

Biogeochemical-Argo (BGC-Argo) floats with multi-year missions and high-frequency sampling offer a promising way to jointly investigate the PIPs and the BGP over a broad range of time and space scales. Such platforms have already been successfully used to characterise the MLP at sub-seasonal to seasonal scales¹²,¹³, the ESP at pan-Antarctic scale⁹, and the seasonality of the BGP²²,²⁴. Here, we jointly investigate these export pathways by refining and linking a range of previously developed techniques. The detection of subsurface anomalous features related to the MLP and ESP⁸,¹² uses concurrent measurements of chlorophyll *a* (Chl) fluorescence and particulate backscattering (proxies of phytoplankton concentration and POC), along with oxygen, temperature and salinity (see details in figure captions and Methods). For the BGP, we use two complementary approaches. Spikes in signals of backscattering and chlorophyll *a* fluorescence provide insights into the sinking of aggregates through the mesopelagic²⁵. Analysis of spikes were used in conjunction with short term (order of days) particle accumulation on transmissometer sensors, acting as optical sediment traps (OST; see details in figure captions and Methods)²⁶. This combination of approaches links a bespoke sensor constellation with recent advances in bio-optical signal analysis to compare the BGP with two PIPs, the MLP and the ESP, through three annual cycles. We focus on a single float, which sampled contrasting biogeochemical provinces across the Pacific sector of the SO, to infer the causal mechanisms that set the magnitude and seasonality of each of these pathways.

## Sampling across Southern Ocean biogeochemical provinces

We present in situ observations collected by a BGC-Argo float (WMO 7900791) that travelled >6000 km across the Pacific sector of the SO (Fig. 1a). This float was deployed in May 2016 in sub-polar waters and sampled the water column from 500 m to the surface every 1–5 days during three annual cycles. It moved south across the Polar Front in September 2016 as evidenced by the abrupt decrease in surface temperature (from 6 °C to 3 °C, Fig. 1c). The float remained in the vicinity of the Polar Front, a region of intense eddy kinetic energy (Fig. S12), for almost an annual cycle, until July 2017. Then, it entered polar waters where it measured near-zero surface temperatures. From June to August 2018, the float was close to the sea-ice edge (Fig. 1b) in an area characterized by very low salinity, and highly stratified surface waters (Fig. 1d and S1c) which prevented any deep mixing events during winter. Indeed, the mixing layer depth, defined here as the maximum vertical gradient in Chl (MLD<sub>bio</sub>, see Methods), remained shallow (<100 m) as opposed to the two previous winters where mixing reached 250–350 m (the black thick line in Fig. 1c, d). The float ended its mission in May 2019. This long trajectory across different oceanic

provinces enabled the characterization of the export pathways over a broad range of time and space scales in contrasting environments. The downside of this long sampling trajectory is the difficulty in studying the seasonality of these pathways when the float moves across different water masses. The trajectory was therefore segmented into three bloom cycles during which the contiguous nature of water masses was verified (Fig. S2 and Methods).

## Characteristics of the particle assemblage

Environmental forcing and ocean physics partly shape the seasonality of phytoplankton stocks and community composition, hereafter called phenology, and by extension the characteristics of the particle assemblage in the upper ocean²⁷. These characteristics will ultimately influence how particles are transported into the oceans' interior via gravitational settling or physical injection⁵,²⁰. We thus commence by discussing the differences in particle characteristics across the three annual cycles before analysing, in the next section, their ramifications for the seasonality of the BGP, the ESP and the MLP.

The contrasting environments between the three annual cycles corresponded to clear differences in bio-optical proxies for concentration, size and composition of phytoplankton community and associated particles (Fig. 2). The 2016–17 bloom, near the Polar Front, was the most intense as reflected in POC, derived from the particulate backscattering coefficient $b_{bp}$ (up to 175 mg C m⁻³, Fig. 2a and Methods), and the concentration of small (<100 µm) fluorescing particles $Chl_s$ (Fig. S9 and Methods). This bloom was characterised by very low $Chl_s/b_{bs}$ values (where $b_{bs}$ is the backscattering coefficient of <100 µm particles), a ratio which varies according to changes in phytoplankton community structure, photoacclimation and nutrient status²⁸⁻³⁰. Here, we attribute low $Chl_s/b_{bs}$ values to the growth of coccolithophores, calcifying phytoplankton which form liths, bio-mineral shells with a high refractive index and thus high backscattering signal³¹. Coccolithophore blooms are recurrent features in the SO, conspicuous from space³²,³³. Indeed, satellite records confirmed the presence of an intense coccolithophore bloom during this summer in the area where the float was profiling (Fig. S3). Bio-minerals act as ballast by increasing particle specific gravity and sinking speeds, and protecting POC from remineralisation³⁴, enhancing the transfer efficiency of particles through the mesopelagic. This likely explains the massive invasion of POC to 450 m following the coccolithophore bloom (Fig. 2a). Such particle characteristics should typically increase the relative importance of the BGP compared to physically driven pathways³⁵. Our POC estimates might be overestimated here, due to high backscattering by these bio-minerals.

The 2017–18 bloom was less intense with lower POC levels (<50 mg C m⁻³, Fig. 2a), but high $Chl_s/b_{bs}$ values at its peak (apex, Fig. 2a, b). Such surface high $Chl_s/b_{bs}$ ratios associated with high-latitude spring blooms have been attributed to diatom-dominated events²⁸,²⁹, which in this case would have been stimulated by dissolved iron supply from winter deep mixing and sea-ice melt (Fig. 1d). An iron stress index ($\alpha_{NPQ}$) was derived from the non-photochemical quenching signal in float fluorescence profiles³⁶ (see Methods). Higher $\alpha_{NPQ}$ values in Fig. 2c correspond to greater phytoplankton community iron stress. The decrease in POC (Fig. 2a) and mean size of particles (Fig. 2c and Methods), in the surface layer after the bloom apex in December 2017 (Fig. 2a, c), coincided with an increase in phytoplankton iron stress (Fig. 2c, grey bars), suggesting that the diatom bloom terminated due to iron limitation. Iron status is a key factor that influences particle characteristics. Indeed, the evolution from iron-replete to iron-limited phytoplankton as the bloom progresses is generally correlated with a change in the buoyancy of particles produced in the surface layer, from neutrally buoyant particles to fast-sinking aggregates³⁷.

The 2018–19 bloom was only partially sampled as data were missing for nearly two months early in the productive season, likely due to the presence of sea-ice which prevented any data transmission.
<figure><img src="image_6.png" /><figcaption>**Fig. 1 |** Contrasting environmental conditions in the Pacific sector of the **Southern Ocean**. **a** Float surfacing positions (every 1–5 days) during its 36-month mission. Background map is a climatology of Absolute Dynamic Topography (ADT). Light grey line indicates the climatological position of the Polar Front (ADT = −0.48 m) used in Ardyna et al.⁷² and derived from Swart et al.⁷³. **b** Minimum</figcaption></figure>

distance between float positions and the sea ice edge, defined as the 15% sea ice concentration limit. Vertical sections of c temperature and d salinity recorded by the float. Light grey lines are isopycnals and the dark line is the mixing layer depth as defined by biological criteria (MLD<sub>bio</sub>, see Methods). Periods with missing data are blank in panels c and d.

In January/February 2019, the presence of a subsurface chlorophyll maximum, characterized by high Chl<sub>s</sub>/b<sub>bs</sub> values at a depth of MLD<sub>bio</sub> due to photoacclimation (Fig. 2b), relatively small particles in the surface layer (Fig. 2c), and high iron limitation index (Fig. 2c) suggest that the bloom was at a late stage when sampled, typical of open-water conditions in the seasonal sea-ice zone³⁸. The second intense peak of POC recorded in April 2019 was observed in the Drake Passage, near the Antarctic Peninsula, a region known to be influenced by iron-enriched shelf waters³⁹. Accordingly, the iron limitation index dropped to low values similar to those observed during the first 2016–17 bloom (Fig. 2c).

## Seasonality of multiple carbon export pathways

Here, we explore the effect of phytoplankton phenology and seasonality in particle characteristics on the seasonality of the BGP, ESP and MLP.

The BGP was characterised by two independent methods, the optical spike and the optical sediment trap approaches. Optical spikes were used to quantify the concentration of large (>100 µm) fast-sinking fluorescing (Chl<sub>l</sub>) and backscattering (b<sub>bl</sub>) particles in the water column (Fig. S5), from which downward fluxes F<sub>Chl spike</sub> and F<sub>POC spike</sub> were derived by multiplying these concentrations with typical mesopelagic particle sinking speeds in high latitude blooms (see Methods). A constant bulk sinking speed is assumed here because information on time-varying sinking speeds was not available. However, uncertainty in bulk sinking speed derived from plume tracking was propagated to uncertainty in our F<sub>POC spike</sub> estimates (see Methods). During the 2016–17 bloom, both F<sub>Chl spike</sub> and F<sub>POC spike</sub> peaked prior to the maximum in depth-integrated POC (i.e., stock) over the mixing layer (early January, Fig. 3a) that was attributed to the coccolithophore bloom. F<sub>POC spike</sub> reached a maximum value of ~900 mg C m<sup>-2</sup> d<sup>-1</sup> in early December when the mean particle size and Chl<sub>s</sub> stock in the upper
<figure><img src="image_7.png" /></figure>

**Fig. 2 | Characteristics of the particle assemblage over three annual cycles, quantified by float-derived bio-optical proxies.** Timeseries of a particulate organic carbon (POC) derived from b<sub>bp</sub>, and b chlorophyll to backscattering ratio Chl<sub>s</sub> / b<sub>bs</sub>. Black lines show the mixing layer depth (MLD<sub>bio</sub>). c Mean particle size (diameter) in the upper 50 m estimated from high-frequency variations in b<sub>bp</sub> and the beam attenuation coefficient c<sub>p</sub> using the variance-to-mean ratio method<sup>52</sup> (see Methods). The blue line represents the median and vertical lines represent inter-

quartile ranges, over 10-d bins. Grey bars show the iron stress index α<sub>NPQ</sub> averaged over 30-d bins. α<sub>NPQ</sub> was derived from the non-photochemical signal in float fluorescence profiles (see Methods). As milestones, triangles on top of each panel mark the timing of onset (first positive net phytoplankton accumulation), climax (maximum accumulation) and apex (maximum phytoplankton carbon stock) for each bloom, following the method in Uchida et al. 2019<sup>68</sup> (see Methods).

layer were also maximum (Figs. 2c and 3a). The high F<sub>Chl spike</sub>/F<sub>POC spike</sub> ratio in the mesopelagic (Fig. S7) suggests that the exported material was dominated by fresh labile phytoplankton aggregates. Inversely, low F<sub>Chl spike</sub>/F<sub>POC spike</sub> following the coccolithophore bloom indicates an export of detrital matter with potentially high mineral content such as coccoliths associated with fecal pellets. Therefore, the seasonal succession in dominant phytoplankton groups, reflected in the different timing of Chl<sub>s</sub> and POC peaks, influences the lability of aggregates in the mesopelagic and potentially their penetration depth (Fig. S6). Despite the lowest surface POC and Chl<sub>s</sub> stocks of all three blooms, the 2017–18 bloom had relatively high F<sub>POC spike</sub> (~300 mg C m<sup>-2</sup> d<sup>-1</sup>, similar to the 2018–19 bloom) and F<sub>Chl spike</sub> values, and the highest F<sub>Chl spike</sub>/F<sub>POC spike</sub> ratio. Iron limitation during the 2017–18 diatom bloom may have led to high level of aggregation as indicated by the high mean particle size in the surface layer (Fig. 2c), also reported in another study in the subarctic Pacific<sup>37</sup>.

The OST flux measured at 300 m was divided into two components, a continuous flux of small slow-sinking particles and a pulsed flux of large fast-sinking particles (see Methods)<sup>26</sup>. The 2016–17 bloom illustrates the decoupling between the continuous and pulsed fluxes

(Fig. 3c). The latter reached ~1100 mg C m⁻² d⁻¹ in mid-January following the coccolithophore bloom, while the continuous flux peaked 20–30 days later in early February 2017 and was an order of magnitude smaller. The delay between the continuous and pulsed flux can be explained by the differential sinking speed between small and large particles, but may also reflect disaggregation processes in the mesopelagic where large particles from the pulsed flux fragment into smaller particles later detected in the continuous flux record<sup>25</sup>. Another important feature revealed by the OST is the relatively high continuous flux in late summer and fall 2018 when mixing layer POC (Fig. 3a) and particle size (Fig. 2c) were close to the annual minimum. This late long-lasting flux of slow-sinking small particles contributed ~60% of the annual continuous flux and was equivalent to ~25% of the total OST flux (continuous + pulsed; Fig. S13).

The eddy subduction pump contribution to the BCP was quantified by first identifying subducted water parcels along the float trajectory. Subducted water parcels were detected as subsurface anomalies in POC, spice (density-compensated changes in temperature and salinity) and oxygen profiles (Fig. S11 and Methods). Such features were detected only during 2016–17 in the vicinity of the Polar
<figure><img src="image_8.png" /><figcaption>**Fig. 3** | Seasonality of multiple carbon export pathways. **a** Stock of particulate organic carbon (POC) and small-particle chlorophyll (Chl<sub>s</sub>) in MLD<sub>bio</sub>. **b** Flux of large-particle POC (blue bars) and chlorophyll (green bars) in a 100-m bin below MLD<sub>bio</sub>, derived from b<sub>bl</sub> and Chl<sub>l</sub> respectively (see Methods). **c** Continuous flux of small slow-sinking particles (slowly accumulating on the transmissometer window; red bars) and pulsed flux of large fast-sinking particles (creating discontinuities in transmissometer records; blue bars) measured by the optical sediment trap (OST) at -300 m (see Methods). Panels **b** and **c** show two complementary approaches to</figcaption></figure>

characterise the BGP. **d** POC export driven by two particle injection pumps (PIPs): the mixed layer pump (MLP, brown bars) and the eddy subduction pump (ESP, yellow bars). In all panels, bars represent a 10-d median of the stock or flux. The length of the error bars represents the interquartile range (see Methods). Similarly to Fig. 2, triangles on top of panel **a** mark the timing of onset (first positive net phytoplankton accumulation), climax (maximum accumulation) and apex (maximum phytoplankton carbon stock) for each bloom, following the method in Uchida et al.<sup>68</sup> (see Methods).

Front, a highly energetic region prone to (sub)mesoscale circulation and ESP events⁶ (Fig. S12). No clear seasonal patterns were observed in ESP fluxes ($F_{ESP}$), with a maximum export of ~180 mg C m⁻² d⁻¹ on average in 10-d bins in June, September and December 2016, January and April 2017 (Fig. 3d). This temporal distribution suggests that a great diversity of particle types can be exported trough this pathway. Indeed, ESP events in September occurred early in the productive season prior to the bloom climax (Fig. S13). Thus, these events actively transported small and freshly produced organic material as deep as ~400 m (Figs. 2c, 3d and S13). Such labile particles would likely not have reached these depths via the BGP, which illustrates the biogeochemical significance of the ESP. A second train of ESP events occurred concurrently with the main peak in surface Chl<sub>s</sub> stock and mean particle size (December–January 2016) suggesting that the ESP also transported large fresh aggregates produced in the euphotic zone, in addition to small suspended particles and dissolved compounds¹³. We cannot rule out the hypothesis that subsurface anomalous water particles sampled by the float were potentially subducted weeks/months earlier and persisted at depth, as suggested in a study in the North Atlantic²¹. However, the strict detection thresholds that we applied for oxygen, spice and POC anomalies should select only recent ESP events.

The mixed layer pump flux ($F_{MLP}$) was quantified as the sum of the POC stock detrained from and entrained into the productive layer due to changes in the mixing layer depth during a given time period (see Methods)¹³,¹⁴. The mean flux over 10-d bins is reported in Fig. 3d (brown bars). Negative values indicate that the entrained POC was greater than the detrained POC during the 10-d period.

Conversely, positive values indicate a net POC export out of the productive layer. POC export occurred during the transition between winter deep mixing and spring stratification (2-3 months), as also reported in previous work¹²,¹⁶, with maximum values in late August 2016 (~140 mg C m⁻² d⁻¹). Lower seasonal variation in the mixing layer depth in 2018–19 due to the influence of melt water resulted in lower $F_{MLP}$ values (Fig. 3d). Time-integration of $F_{MLP}$ over complete annual cycles demonstrates that the MLP is an efficient pathway to transfer POC below the maximum winter mixing layer depth (Fig. S13). Indeed, the net export of POC was 5.7 g C m⁻² in 2016–17 and 2.3 g C m⁻² in 2017–18 (not calculated in 2018–19 due to missing data for nearly two months). The strong sub-seasonal variability of the mixing layer depth, with intermittent stratification and deep mixing events, also led to transient inputs of fresh organic material as deep as 350 m, as revealed by high Chl<sub>s</sub>/b<sub>BS</sub> ratio in the mesopelagic zone (Fig. S8). These transient inputs may sustain the mesopelagic biota before the main export pulse following the spring bloom¹²,¹³. The last export events, in early November 2016 and late October 2017, before the permanent summer stratification occurred at a time when particles in the surface started to increase in size (Fig. 2c), suggesting that relatively large particles can also be actively transferred to depth, similar to the ESP. The clear difference in the depth strata of particles injected by the two PIPs (Figs. S11 and S14) indicate that POC stocks exported by these physically mediated pathways were not double counted. However, these results illustrate the possible interplay between a PIP and the BGP, where large particles could sink further following a physical injection event.
## Implications for oceanic carbon storage

Multi-year and high-resolution float observations provide a new and more comprehensive view of export pathways that contribute to the BCP over the annual cycle. The main insight from our study is the seasonal succession of mechanisms contributing to a long and sustained transfer of organic carbon to the deep ocean. Export due to the ESP was observed early in the bloom, before the bloom climax, and again at the bloom apex, but only in one of three years. The MLP was most active during the seasonal transition from deep mixing to stratification (between bloom onset and climax), and again during sub-seasonal stratification events. There was large interannual variability in the magnitude of the MLP. Finally, the gravitational pulsed flux was most intense at the bloom apex, while the continuous flux persisted throughout the year, with maximum values delayed in time with respect to the bloom apex—in fall in the case of the 2017–18 bloom. The contribution of all these processes over a complete annual cycle should be accounted for when computing regional mesopelagic carbon budgets.

Total annual POC export from small slow-sinking particles (continuous OST flux) reached $11 \text{ g C m}^{-2}$ in 2016–17 and $8 \text{ g C m}^{-2}$ in 2017–18 (Fig. S13), of which 43 and 27% was exported prior to the bloom apex during the first and second annual cycles, respectively. As for large fast-sinking particles (pulsed OST flux), these relative contributions increase to 65 and 52%, with total annual export of $53 \text{ g C m}^{-2}$ and $5 \text{ g C m}^{-2}$ in 2016–17 and 2017–18, respectively (similar values for the optical spike approach). These results stress the importance of persistent observations of the BGP, not only following spring bloom collapse. In comparison, total annual export from the MLP was $5.7 \text{ g C m}^{-2}$ in 2016–17 and $2.3 \text{ g C m}^{-2}$ in 2017–18 representing 9 and 17% of total export from the BGP (pulsed + continuous OST flux). This contribution falls within the range estimated in the North Atlantic with a similar approach (5–25%)¹⁴, and is slightly lower than the average contribution at high latitudes (23%)¹³. These relative contributions increase to 50 and 165% at the time of the bloom climax in 2016–17 and 2017–18, respectively (Fig. S13), indicating that the MLP is the dominant carbon export pathway during the pre-bloom period when strong variability in the mixing layer depth occurs. Evaluating the total annual export from the ESP is more challenging due to the difficulty to integrate fluxes from episodic and local ESP events²³. Our method also does not estimate net ESP export as upward fluxes are not considered¹⁶. We thus compare the mean gross ESP export with total BGP export solely at the time of the bloom apex. In 2016–17, the ESP contribution to total BGP export was about 19%, a value similar to the North Atlantic spring bloom (–25%)⁸ and at the peak bloom in the Drake Passage (–18%)²³. The absence of ESP events during the second and third annual cycles in less energetic regions (Fig. S12) suggests that, in agreement with Llort et al.⁹, this pathway might be less widespread in the Southern Ocean than previously reported in the model-based study of Omand et al.⁸. On the other hand, episodic and local ESP events might not be efficiently captured by a single BGC-Argo float, in contrast to a glider, which can sample regions of interest such as fronts at high resolution. However, the emerging global float network could overcome this limitation and provide a valuable statistical view of the importance of the ESP at the basin scale.

A final important insight from our study is the identification of the timing and depth strata for the potential interplay between the PIPs and BGP⁴. Float-derived optical proxies revealed that PIPs potentially inject large fresh aggregates as deep as 400 m, in addition to suspended particles and dissolved organic matter. These fresh aggregates can then sink further through gravitational sinking³⁵. With vertical velocities of order of 100–1000 m d⁻¹⁴⁰,⁴¹, physically mediated processes, such as the MLP and ESP, considerably accelerate the transit of particles through the main remineralisation horizon. The interplay between the PIPs and BGP makes it challenging to assess their relative contribution to carbon sequestration, but the overlap that we find in

their timing clearly suggest that their joint contributions can boost the overall efficiency of the BCP.

## Methods

### Float data processing

The BGC-Argo float used is a Teledyne Webb Research APEX float, equipped with a Sea-Bird SBE41 CTD sensor, ECO fluorometer and scattering sensors measuring chlorophyll *a* fluorescence and the volume scattering function (at -124°, 470 and 700 nm wavelengths), OC4 radiometer measuring downwelling irradiance at 412, 443, 490 and 550 nm, C-Rover transmissometer measuring the beam attenuation coefficient at 660 nm ($c_p$, in m⁻¹), and an Aanderaa Oxygen Optode.

The float mission included CTD and bio-optical day and night profiles every 1.5–6 days, from 500 m depth to the surface, with a parking depth of 300 m. The vertical sampling resolution ranged from 3 to 10 m depending on the float ascent speed. All profile data were interpolated at 1 m resolution.

The CTD and trajectory data were quality-controlled using the standard Argo protocol⁴². ECO raw signals were converted to chlorophyll *a* concentration (Chl, in mg Chl m⁻³) and particulate backscattering coefficient ($b_{bp}$, in m⁻¹) following the BGC-Argo procedures⁴³,⁴⁴. Bio-optical data were quality-controlled following the BGC-Argo quality control manual⁴⁵. In addition, Chl was corrected for non-photochemical quenching (NPQ) following Xing et al.⁴⁶. Briefly, for each profile where the sun elevation angle was >5°, the maximum Chl value above the mixed layer depth, defined here as a density difference of 0.01 kg m⁻³ with a reference value at 5 m depth, was extrapolated toward the surface. As an additional condition, the depth of the extrapolated Chl value had to be shallower than the 15 µmol photons m⁻² s⁻¹ isolune⁴⁶. As the instantaneous Photosynthetically Available Radiation (iPAR) was not directly measured, we estimated iPAR profiles from the measured downwelling irradiance at 4 wavelengths. At each depth, a spline interpolation was used to compute the irradiance spectra from 400 to 1000 nm at 1 nm resolution, with irradiance value set to zero at 1000 nm. iPAR was then calculated by integrating the interpolated irradiance spectra from 400 to 700 nm. Oxygen data were calibrated using air measurements following Johnson et al.⁴⁷, with a mean gain factor of 1.0557. The $b_{bp}$ at 700 nm was converted into particulate organic carbon (POC, in mg C m⁻³) using the SO-specific relationship in Johnson et al.⁴⁸ POC = $3.12*10⁴$ (± $2.47*10³$) * $b_{bp700}$ + $3.04$ (± $6.78$). Similarly, the $b_{bp}$ at 470 nm was converted into phytoplankton carbon ($c_{phyto}$, in mg C m⁻³) following Graff et al.⁴⁹: $c_{phyto}$ = $12.128*10³$ * $b_{bp470}$ + $0.59$. Note that these carbon values might be over-estimated during the 2016–17 coccolithophore bloom due to the high refractive index of coccoliths.

In this study, we distinguish between the mixed layer, the zone of relatively homogeneous water formed by the history of mixing, and the mixing layer, the zone in which mixing is currently active. The mixing layer depth was computed as the maximum vertical gradient of Chl (MLD$_{bio}$) following Lacour et al.¹². The underlying concept is that Chl, as a proxy of phytoplankton concentration, is homogeneous over the whole mixing layer if turbulent mixing overcomes vertical variations in phytoplankton net growth rate³⁰. MLD$_{bio}$ should capture the high-frequency variability of the mixing layer at time scales typical of phytoplankton growth⁵¹. Note that in summer stratified waters, phytoplankton can thrive below the mixing layer. In this situation, MLD$_{bio}$ corresponds more likely to the base of the euphotic zone.

### Mean particle size

Mean particle diameter in the upper 50 m, weighted by particle cross-sectional area, was estimated from high-frequency variations in $b_{bp}$ and $c_p$ using the variance-to-mean ratio method⁵², adapted for use on profiling floats following Rembauville et al.²⁹. This method extracts particle size information from the “spikiness” of optical profiles where
particle concentrations are too high to separate individual spikes from a small-particle baseline. Briefly, mean diameter was estimated from c_p via Eq. (1):

$$ \bar{A}_{cp} = \frac{\text{var}(c_{pdetrended})}{\text{mean}(c_p)} \frac{V}{Q_c} \frac{1}{\alpha(\tau)} \quad (1) $$

$$ \alpha(\tau) = \begin{cases} 1 - (3\tau)^{-1}, & \text{if } \tau \ge 1 \\ \tau - \frac{\tau^2}{3}, & \text{if } \tau \le 1 \end{cases} $$ $$ \tau = \left( \frac{t_{res}}{t_{samp}} \right) $$

$$ \bar{D}_{cp} = 2\sqrt{\bar{A}_{cp}\pi^{-1}} $$

where $t_{res}$ is the residence time of a particle in the sample volume (0.1 s), $t_{samp}$ is the duration of a single measurement (1 s), V is the transmissometer sample volume (12.5 ml), $Q_c$ is the optical attenuation efficiency of the particles (assumed to be 2 following Bohren and Huffman$^{53}$), var(c_p detrended) is the variance of c_p after detrending using an 5-point running median, and $D_{cp}$ is the area-weighted mean particle diameter. The calculation was then repeated for b\_bp at both 470 and 700 nm, by replacing c_p with b\_bp in Eq. 1, and using V = 0.62 ml (Briggs et al.$^{52}$), $t_{samp}$ = 1 s, $t_{res}$ = 0.06 s, and $Q_{bb}$ = 0.02$^{52}$. In order to reduce noise in the individual size estimates, the three estimates were then combined together into 10 d bins, whose medians and 25th and 75th percentiles are reported here.

## Optical spikes in the mesopelagic

Chl and b\_bp were partitioned into six components following Briggs et al.$^{54}$: deep sensor blanks, including a background of small refractory particles (Chl_r and b\_br); small, labile fluorescing (Chl_s) and backscattering (b\_bs) particles; and large, fast-sinking fluorescing (Chl_l) and backscattering (b\_bl) particles. The division between small and large corresponds approximately to a particle chlorophyll content of 60 pg for Chl_s versus Chl_l and a particle diameter of 100 µm for b\_bs versus b\_bl. Timeseries of Chl_l and b\_bl at 470 and 700 nm are shown in Supplementary Fig. S5. Chl_l corresponds to fresh phytoplankton aggregates while b\_bl additionally includes fecal and detrital matter$^{25}$. The vertical fluxes of large particles were computed as follows:

$$ F_{Chl spike} = \langle Chl_l \rangle^* w_{Chl_l} $$

$$ F_{POC spike} = \langle b_{bl700} \rangle^* Q_{b_{bl}/POC}^* w_{b_{bl}} $$

where $w_{b_{bl}} = 74(\pm 45)md^{-1}$ and $w_{Chl_l} = 98(\pm 64)md^{-1}$ are the mean (and standard deviation) sinking speeds of large backscattering and fluorescing particles, respectively, as estimated in Briggs et al.$^{54}$. These estimates are based on 18 particle export pulse events for b\_bl and 25 events for Chl_l, including 22 events in total in the Southern Ocean. $Q_{b_{bl}/POC} = 3.12*10^4(\pm 2.47*10^3)$ is the slope of the relationship between b\_bp and POC developed in Johnson et al.$^{48}$. Given that b\_bl is a fraction of b\_bp, the intercept of the relationship was not used. Brackets represent a 10-d bin average of all Chl_l and b\_bl values in a 100-m bin below the mixing layer depth. The sinking speed w is assumed to be constant over time. The probability distribution of the flux F was derived from the mean and standard deviation of each term in the calculation, using a Monte Carlo simulation, with 5000 iterations, and assuming each term is normally distributed. The mean and standard deviation of $\langle Chl_l \rangle$ and $\langle b_{bl700} \rangle$ represent the uncertainty introduced by the spatio-temporal variability. Vertical bars in Fig. 3b show the median of the

distribution of F while the length of error bars represent its interquartile range. Similarly, Fig. 3a shows the median and interquartile range of POC and Chl stocks in the mixing layer, with uncertainties introduced by the b\_bp to POC conversion and the spatio-temporal variability in 10-d bins.

## Optical sediment trap methodology

The Optical Sediment Trap (OST) method uses the rate of change of particle attenuance (ATN = -ln(transmittance), in m$^2$ m$^{-2}$) measured by the transmissometer during the float drifting period at the parking depth of 300 m$^{24,26,55}$. The accumulation of small slow-sinking particles onto the upward-facing window of the transmissometer drives a smoothly increasing ATN while the accumulation of large fast-sinking particles produces discontinuities in transmissometer records. Both signals were converted into a continuous and a pulsed ATN flux (in units of m$^2$ m$^{-2}$ d$^{-1}$), respectively, following procedures described in Estapa et al.$^{26,55}$ with few modifications to take into account the use of a different platform with different parking behaviour and sampling frequency (every 1 h in this study). First, parking phases of duration <24 h were not considered to ensure enough data points. The first 3 data points (3 h) were removed from the analysis as the float takes time to stabilize at the target depth. Optical spikes, defined as an increase in ATN of 0.002 m$^2$ m$^{-2}$ h$^{-1}$ followed by a decrease within 3 h, were also removed and a 5-point median filter was applied. Then, the remaining data points were divided into linear segments interspersed with discontinuities, or 'jumps', following the procedure in Estapa et al.$^{55}$, except that discontinuities were not identified a priori with a threshold criteria but resulted from the subdivision of segments not meeting the fitting criteria. For each parking phase, the continuous flux of slow-sinking particles was computed as the mean slope of all linear segments weighted by their length. The pulsed flux was computed as the sum of all positive discontinuities normalized by the duration of the parking phase. Finally, both continuous and pulsed ATN fluxes were converted to POC flux (in mg C m$^{-2}$ d$^{-1}$) following Estapa et al. (2019):$^{22}$

$$ F_{OST} = \langle F_{ATN} \rangle^* 2133(\pm 173) + 1.4(\pm 5.6) $$

where $\langle F_{ATN} \rangle$ is the mean of either the continuous or the pulsed attenuance flux in 10-d bins. Values in parentheses are the standard deviations of the parameters of the linear fit. Note that this relationship is based on measurements made in the North Atlantic and may be inappropriate in other basins with very different particle assemblages. A SO-specific relationship does not currently exist. As for the $F_{spike}$, the probability distribution of $F_{OST}$ was computed using a Monte Carlo approach, with uncertainties introduced by the ATN to POC conversion and the spatio-temporal variability in 10-d bins.

Figure 3b, c allows us to compare the two flux estimates of the BGP, namely $F_{OST}$ and $F_{POC spike}$. For the three blooms, the pulsed $F_{OST}$ and $F_{POC spike}$ were of similar order of magnitude, but the temporal dynamics were different, especially during the 2016–17 bloom when $F_{POC spike}$ peaked in early December, ~30 days before the pulsed $F_{OST}$. This difference is likely due to the difference in the depth of the measured flux (100-m bin below the mixing layer versus the parking depth of 300 m). Indeed, the flux of fresh aggregates recorded in $F_{Chl spike}$ and $F_{POC spike}$ was quickly attenuated with depth (Fig. S6), and was therefore less intense at 300 m compared to the flux following the coccolithophore bloom in early January.

## Mixed layer pump

The methodology to estimate the magnitude of the mixed layer pump (MLP) was inspired by Dall’Olmo et al.$^{13}$ and Bol et al.$^{14}$. The stock of entrained/detrained POC due to the change in the mixing layer depth
between two subsequent time steps was calculated as follows:

$$\int POC(t) = \begin{cases} + \int_{\max(MLD_{bio}(t+1),100 m)}^{MLD_{bio}(t)} POC(z,t) dz & \text{if } MLD_{bio}(t)>MLD_{bio}(t+1) \\ - \int_{\max(MLD_{bio}(t+1),100 m)}^{MLD_{bio}(t+1)} POC(z,t) dz & \text{if } MLD_{bio}(t+1)>MLD_{bio}(t) \end{cases}$$

In this calculation, $\int POC$ is positive when the mixing layer shoals ($MLD_{bio}(t) > MLD_{bio}(t+1)$) and POC stock is detrained from the mixing layer. The layer delimited by $MLD_{bio}(t)$ at the bottom and $MLD_{bio}(t+1)$ at the top is defined as the remnant layer (Fig. S11). Conversely, $\int POC$ is negative when the mixing layer deepens and POC stock is entrained from deeper waters. Following Dall'Olmo et al.¹³, POC stocks from the surface to 100 m were not included in the calculation to avoid accounting for POC produced below the mixing layer, typically in summertime stratified waters. The nominal depth of 100 m corresponds approximately to the base of the euphotic layer in high latitude regions and is also consistent with other studies estimating carbon export⁵. The sum of $\int POC$ over a given period, divided by the duration of the period, provides an estimate of the net carbon export out of the productive layer $F_{MLP}$, expressed in mg C m⁻² d⁻¹. This approach assumes a spatial homogeneity in physical and biogeochemical properties over the area covered by the float displacement during the given period⁵⁶. To smooth out any (sub)mesoscale spatial heterogeneities, $F_{MLP}$ was computed using a 20-d centred moving sum, and then averaged over 10-d bins to be consistent with flux estimates from the other pathways. Following a Monte Carlo approach, this process was repeated 5000 times, with uncertainties introduced by the $b_{bp}$ to POC conversion. Error bars in Fig. 3 highlight the uncertainties associated to both the conversion to POC and the spatio-temporal variability in 10-d bins.

## Eddy subduction pump

To detect subsurface features resulting from ESP events, we used an updated version of the method described in Llort et al.⁹ and Omand et al.⁸. This method relies on the observation that the vertical extent of submesoscale features that drive ESP is of the order of -10 m and can be smoothed out by averaging the float's vertical variability over larger (i.e., -100 m) vertical scales. We applied a 20-bin centred running median over single vertical profiles interpolated at 5 m vertical resolution. As a result, we obtained submesoscale-free vertical variability that can be compared to the original profiles to identify anomalous features. We focused only on anomalies found between the bottom of the mixing layer depth and 500 m. As anomalies can be related to other mechanisms than submesoscale vertical circulation, individual profiles were classified as an ESP event if negative anomalies of Apparent Oxygen Utilisation (AOU') and spice ($\pi'$) coincided in depth with a positive POC anomaly (POC'). If more than one anomaly was detected in the same profile, we only considered the shallowest one. The running median was computed with a centred filter, i.e., +/-10 bin median for each bin depth. To define the relevant anomalies, we applied detection thresholds to only consider anomalies with AOU' < -5 µmol kg⁻¹, POC' > 0 and $\pi'$ < -0.05. Tests with positive spiciness anomalies (i.e., $\pi'$ > 0.05) resulted in no new events detected. This result is likely because, in the Southern Ocean, spiciness is negative at the surface and decreases with depth⁵⁷, unlike other regions such the North Atlantic where Omand et al.⁸ found positive $\pi'$ anomalies associated with ESP events. We have also modified the method to better constrain the thickness of the anomalous features, an important metric to estimate the amount of POC exported. Here, we detected the top and bottom depths of both AOU' and $\pi'$ as the first and last bin depths where the anomalies were still negative. We then compared the detected depths on AOU' and $\pi'$ to define the deepest (shallowest) of the two as the anomaly bottom (top) depth.

For each detected feature, the associated vertical C export was computed as follows:

$$F_{ESP} = POC_{surf} * \langle w_{sub} \rangle$$

where $\langle w_{sub} \rangle$ is the mean vertical velocity of the injected water parcel, defined as:

$$\langle w_{sub} \rangle = (MLD_{bio} - Z_{ev}) / \Delta t$$

where

$$\Delta t = (POC_{surf} - POC_{ev}) / R$$

$R$ is the respiration rate of organic matter in the mesopelagic and is equal to 0.33(±0.16)mg Cm⁻²d⁻¹, the mean (and standard deviation) respiration rate between 100 and 400 m estimated from BGC-Argo floats in Southern Ocean open waters (Hennon et al.⁵⁸, their Figure 4). $POC_{surf}$ is the concentration of POC in the mixing layer averaged over all float profiles in the 10-d bin containing the detected subsurface feature. The standard deviation of POC concentration was also used to assess the spatio-temporal variability of $POC_{surf}$. Likewise, the mean and standard deviation of $MLD_{bio}$ were calculated in the same 10-d bin. $Z_{ev}$ and $POC_{ev}$ are the median depth and POC concentration of the subsurface feature, respectively. Following a similar approach as for the other export pathways, the probability distribution of $F_{ESP}$ was computed using a Monte Carlo simulation, with uncertainties introduced by the respiration rate $R$, the $b_{bp}$ to POC conversion and the spatio-temporal variability of $POC_{surf}$ and $MLD_{bio}$. This process was repeated for each subsurface feature detected in 10-d bins. Figure 3 shows the resulting median and interquartile range of $F_{ESP}$ for each 10-d bin. Note that $F_{ESP}$ is a gross export and does not include potential upward fluxes associated with ESP events. This approach relies on 4 assumptions:

1. The probability distribution of $POC_{surf}$ was representative of the $POC_{ev}$ concentration in the water parcel before being subducted.

2. The difference between $POC_{surf}$ and $POC_{ev}$ was only due to respiration of organic matter since the water parcel left the surface mixed layer.

3. The average respiration rate in the water parcel was constant in the water parcel.

4. The distance travelled by the water parcel since leaving the mixing layer until being detected was equivalent to the difference between the depth of detection $Z_{ev}$ and $MLD_{bio}$.

Note that all ESP features were detected much deeper than MLP remnant layers resulting in no double counting of physically driven POC export by the MLP and ESP (Fig. S14).

## Iron stress index

The phytoplankton iron stress index was computed following the method in Ryan-Keogh & Thomalla³⁶. The concept underpinning this method is that NPQ variability is linked to both iron and light availability and has the potential to provide important diagnostic information on phytoplankton physiology³⁹. To remove the effect of in situ light availability on NPQ variability, Ryan-Keogh & Thomalla³⁶ proposed to compute $\alpha_{NPQ}$ the initial slope of the NPQ-PAR curve. Thereby, $\alpha_{NPQ}$ could be used as a proxy for iron limitation, with higher values being associated with greater phytoplankton iron stress. In our study, NPQ as a function of depth was quantified as the difference between the quenching corrected fluorescence profile and the quenched one, normalized by the latter. For each profile, we plotted our iPAR estimates against NPQ values. We then applied a linear fit to the NPQ-iPAR curve in the region of low iPAR values (between 15 and 75 µmol photons m⁻² s⁻¹), where the slope of the linear fit gives $\alpha_{NPQ}$. Linear fits with $R^2 < 0.8$
were rejected (1% of the profiles). We did not fit a Platt-like model as in Ryan-Keogh & Thomalla³⁶ because our iPAR values were too low to induce a saturation plateau and because we were interested only in the initial slope of the NPQ-iPAR curve. Note that αNPQ, as a proxy for iron limitation, has to be interpreted with caution as shifts in phytoplankton community composition and changes in the light regime and thus phytoplankton photoacclimation status impact the variability of αNPQ⁶⁰,⁶¹. For that reason, we also computed the median light level within the mixing layer ($I_{ML} = PAR_{SAT}e^{-0.5K_{d(PAR)}MLD_{bio})$ which is often used in photoacclimation models⁶². Here, $PAR_{SAT}$ is the daily mean PAR from MODIS Aqua (4 km) and $K_{d(PAR)}$ is the diffuse attenuation coefficient of PAR (in units of m⁻¹). $K_{d}$ was first computed at 490 nm by fitting a fourth-degree polynomial function to the logarithm of the downwelling irradiance $E_{d(490)}$ as a function of depth, measured by the float, and then calculating the mean slope over the first 50 m⁶³,⁶⁴. $K_{d(490)}$ was then converted to $K_{d(PAR)}$ following equation 9 in Morel et al. (2007)⁶⁵. We did not find a clear relationship between $I_{ML}$ and αNPQ (Fig. S4). The influence of shifts in phytoplankton community composition on αNPQ variability cannot be ruled out, but phytoplankton community composition also shifts in response to iron availability⁶⁶.

## Phenology metrics

Phytoplankton phenology can be characterised by three metrics:⁶⁷ (i) the bloom onset when phytoplankton carbon stock starts accumulating (i.e. when the accumulation rate r changes from negative to positive), (ii) the bloom climax when r reaches its annual maximum value, (iii) the bloom apex when phytoplankton carbon stock P reaches its annual maximum value. Following Uchida et al.⁶⁸, P was estimated by vertically integrating $C_{phyto}$, derived from the backscattering at 470 nm, over the whole water column. To avoid including non-phytoplankton particulate matter (e.g. refractory material) in the calculation, $C_{phyto}$ was masked out at depths where $Chl_s = Chl - Chl_l - Chl_r \le 0$, where $Chl_r = 0.3 mgm^{-3}$ (Fig. S9). The timeseries of $C_{phyto}$ was first linearly interpolated on equally spaced 5-day timeseries, and then smoothed with a 30-day (7 points) running average to filter out short-term fluctuations and focus on seasonality. The accumulation rate was then calculated as $= \frac{1}{P} \frac{6P}{6t}$. The bloom onset and apex were used as milestones to compare the timing of the 2016–17 and 2017–18 blooms with the seasonality of the carbon export pathways (Figs. S10 and S13). The time axis in Fig. S13 was rescaled by the onset and the apex of the bloom, so that 0 corresponds to the onset and 1 to the apex.

## Quasi-Lagrangian framework of the float trajectory

During its 36-month mission, the float visited different oceanic provinces and crossed strong water mass boundaries, such as in September 2016 when the float crossed the Polar Front. Thus, observed changes in biogeochemical properties were not solely due to temporal changes, confounding the study of the seasonality of these properties. We therefore divided the timeseries into three periods in which the contiguous nature of the water masses was verified based on temperature and salinity properties (24-Sep-2016 to 22-May-2017; 10-Oct-217 to 7-Jun-2018; 5-Nov-2018 to 1-May-2019; Fig. S2). The absence of strong water mass contrasts allows us to assume a quasi-Lagrangian framework, where changes in biogeochemical properties can be interpreted as temporal changes. This approach is commonly used in float studies⁶⁹⁻⁷¹. For completeness, the Figures show the full float timeseries. However, we focus the Discussion on the three quasi-Lagrangian periods where we can confidently interpret the relationship between phytoplankton phenology in the surface, and the export and fate of POC in the underlying mesopelagic zone. All carbon flux estimates from the different export pathways remain robust for the whole timeseries.

## Data availability

BGC-Argo data used in this study can be downloaded from the Argo Global Data Assembly Center (ftp://ftp.ifremer.fr/ifremer/argo/). These data were collected and made freely available by the International Argo Program and the national programs that contribute to it: (http://www.argo.ucsd.edu, https://www.ocean-ops.org). The Argo Program is part of the Global Ocean Observing System. Ssalto/Duacs altimeter product was produced and distributed by the Copernicus Marine and Environment Monitoring Service (CMEMS, https://marine.copernicus.eu/). ERA5 wind speed and net heat flux reanalysis products were downloaded from https://cds.climate.copernicus.eu/cdsapp#/home. Monthly GlobColour products were downloaded from CMEMS and MODIS products from the NASA Ocean Color website (https://oceancolor.gsfc.nasa.gov). The sea-ice product was downloaded from https://nsidc.org/data/nsidc-0081.

## References

1. Arteaga, L., Haëntjens, N., Boss, E., Johnson, K. S. & Sarmiento, J. L. Assessment of export efficiency equations in the Southern Ocean applied to satellite-based net primary production. *J. Geophys. Res. Ocean.* **123**, 2945–2964 (2018).

2. Schlitzer, R. Carbon export fluxes in the Southern Ocean: results from inverse modeling and comparison with satellite-based estimates. *Deep. Res. Part II Top. Stud. Oceanogr.* **49**, 1623–1644 (2002).

3. Volk, T. & Hoffert, M. I. Ocean carbon pumps: analysis of relative strengths and efficiencies in ocean-driven atmospheric CO₂ changes. in *The Carbon Cycle and Atmospheric CO 99–110 (American Geophysical Union, 1985)*. https://doi.org/10.1029/gm032p0099

4. Boyd, P. W., Claustre, H., Levy, M., Siegel, D. A. & Weber, T. Multifaceted particle pumps drive carbon sequestration in the ocean. *Nature* **568**, 327–335 (2019).

5. Siegel, D. A., Devries, T. & Cetinić, I. Quantifying the Ocean's biological pump and its carbon cycle impacts on global scales. *Ann. Rev. Mar. Sci.* **15**, 329–356 (2022).

6. Burd, A. B. et al. Assessing the apparent imbalance between geochemical and biochemical indicators of meso- and bathypelagic biological activity: What the @#! is wrong with present calculations of carbon budgets? *Deep. Res. Part II Top. Stud. Oceanogr.* **57**, 1557–1571 (2010).

7. Emerson, S. Annual net community production and the biological carbon flux in the ocean. *Global Biogeochem. Cycles* **28**, 14–28 (2014).

8. Omand, M. M. et al. Eddy-driven subduction exports particulate organic carbon from the spring bloom. *Science* **348**, 222–225 (2015).

9. Llort, J. et al. Evaluating Southern Ocean carbon eddy-pump from biogeochemical-argo floats. *J. Geophys. Res. Ocean.* **123**, 971–984 (2018).

10. Stukel, M. R. & Ducklow, H. W. Stirring up the biological pump: vertical mixing and carbon export in the southern ocean. *Global Biogeochem. Cycles* **31**, 1420–1434 (2017).

11. Erickson, Z. K. & Thompson, A. F. The seasonality of physically driven export at submesoscales in the Northeast Atlantic Ocean. *Global Biogeochem. Cycles* **32**, 1144–1162 (2018).

12. Lacour, L., Briggs, N., Claustre, H., Ardyna, M. & Dall'Olmo, G. The intraseasonal dynamics of the mixed layer pump in the subpolar North Atlantic Ocean: a biogeochemical-argo float approach. *Global Biogeochem. Cycles* **33**, 266–281 (2019).

13. Dall'Olmo, G., Dingle, J., Polimene, L., Brewin, R. J. W. & Claustre, H. Substantial energy input to the mesopelagic ecosystem from the seasonal mixed-layer pump. *Nat. Geosci.* **9**, 820–823 (2016).

14. Bol, R., Henson, S. A., Rumyantseva, A. & Briggs, N. High-frequency variability of small-particle carbon export flux in the Northeast Atlantic. *Global Biogeochem. Cycles* **32**, 1803–1814 (2018).
15. Stukel, M. R., Song, H., Goericke, R. & Miller, A. J. The role of subduction and gravitational sinking in particle export, carbon sequestration, and the remineralization length scale in the California Current Ecosystem. *Limnol. Oceanogr.* **63**, 363–383 (2018).

16. Resplandy, L., Lévy, M. & McGillicuddy, D. J. Effects of eddy-driven subduction on ocean biological carbon pump. *Global Biogeochem. Cycles* **33**, 1071–1084 (2019).

17. Irigoien, X. et al. Large mesopelagic fishes biomass and trophic efficiency in the open ocean. *Nat. Commun.* **5**, 3271 (2014).

18. Jónasdóttir, S. H., Visser, A. W., Richardson, K. & Heath, M. R. Seasonal copepod lipid pump promotes carbon sequestration in the deep North Atlantic. *Proc. Natl. Acad. Sci. USA* **112**, 12122–12126 (2015).

19. Levy, M. et al. Physical pathways for carbon transfers between the surface mixed layer and the ocean interior. *Global Biogeochem. Cycles* **27**, 1001–1012 (2013).

20. Durkin, C. A. et al. A visual tour of carbon export by sinking particles. *Global Biogeochem. Cycles* **35**, 1–17 (2021).

21. Johnson, A. R. & Omand, M. M. Evolution of a subducted carbon-rich filament on the edge of the North Atlantic Gyre. *J. Geophys. Res. Ocean.* **126**, e2020JC016685 (2021).

22. Estapa, M. L., Feen, M. L. & Breves, E. Direct observations of biological carbon export from profiling floats in the Subtropical North Atlantic. *Global Biogeochem. Cycles* **33**, 282–300 (2019).

23. Davies, A. R., Veron, F. & Oliver, M. J. Biofloat observations of a phytoplankton bloom and carbon export in the Drake Passage. *Deep. Res. Part I Oceanogr. Res. Pap.* **146**, 91–102 (2019).

24. Bishop, J. K. B. & Wood, T. J. Year-round observations of carbon biomass and flux variability in the Southern Ocean. *Global Biogeochem. Cycles* **23**, 1–12 (2009).

25. Briggs, N. et al. High-resolution observations of aggregate flux during a sub-polar North Atlantic spring bloom. *Deep. Res. Part I Oceanogr. Res. Pap.* **58**, 1031–1039 (2011).

26. Estapa, M., Durkin, C., Buesseler, K., Johnson, R. & Feen, M. Carbon flux from bio-optical profiling floats: calibrating transmissometers for use as optical sediment traps. *Deep. Res. Part I Oceanogr. Res. Pap.* **120**, 100–111 (2017).

27. Trudnowska, E. et al. Marine snow morphology illuminates the evolution of phytoplankton blooms and determines their subsequent vertical export. *Nat. Commun.* **12**, 2816 (2021).

28. Cetinić, I. et al. A simple optical index shows spatial and temporal heterogeneity in phytoplankton community composition during the 2008 North Atlantic Bloom Experiment. *Biogeosciences* **12**, 2179–2194 (2015).

29. Rembauville, M. et al. Plankton assemblage estimated with BGC-argo floats in the Southern Ocean: implications for seasonal suc-cessions and particle export. *J. Geophys. Res. Ocean.* **122**, 8278–8292 (2017).

30. Geider, R. J. Light and temperature dependence of the carbon to chlorophyll a ratio in microalgae and cyanobacteria: implications for physiology and growth of phytoplankton. *New Phytol.* **106**, 1–34 (1987).

31. Terrats, L., Claustre, H., Cornec, M., Mangin, A. & Neukermans, G. Detection of coccolithophore blooms with biogeochemical-argo floats. *Geophys. Res. Lett.* **47**, e2020GL090559 (2020).

32. Balch, W. M. et al. The contribution of coccolithophores to the optical and inorganic carbon budgets during the Southern Ocean Gas Exchange Experiment: new evidence in support of the Great Calcite Belt hypothesis. *J. Geophys. Res. Ocean.* **116**, 1–14 (2011).

33. Balch, W. M., Gordon, H. R., Bowler, B. C., Drapeau, D. T. & Booth, E. S. Calcium carbonate measurements in the surface global ocean based on Moderate-Resolution Imaging Spectroradiometer data. *J. Geophys. Res. C Ocean.* **110**, 1–21 (2005).

34. Klaas, C. & Archer, D. E. Association of sinking organic matter with various types of mineral ballast in the deep sea: Implications for the rain ratio. *Global Biogeochem. Cycles* **16**, 63-1-63-14 (2002).

35. Dever, M., Nicholson, D., Omand, M. M. & Mahadevan, A. Size-differentiated export flux in different dynamical regimes in the Ocean. *Global Biogeochem. Cycles* **35**, 1-20 (2021).

36. Ryan-Keogh, T. J. & Thomalla, S. J. Deriving a proxy for iron limitation from chlorophyll fluorescence on buoyancy gliders. *Front. Mar. Sci.* **7**, 1-13 (2020).

37. Boyd, P. W. et al. The evolution and termination of an iron-induced mesoscale bloom in the northeast subarctic Pacific. *Limnol. Oceanogr.* **50**, 1872-1886 (2005).

38. Sedwick, P. N., Di Tullio, G. R. & Mackey, D. J. Iron and manganese in the Ross Sea, Seasonal iron limitation in Antarctic. *J. Geophys. Res. Ocean.* **105**, 11321-11336 (2000).

39. Hopkinson, B. M. et al. Iron limitation across chlorophyll gradients in the southern Drake Passage: phytoplankton responses to iron addition and photosynthetic indicators of iron stress. *Limnol. Oceanogr.* **52**, 2540-2554 (2007).

40. D'Asaro, E. A. Convection and the seeding of the North Atlantic bloom. *J. Mar. Syst.* **69**, 233-237 (2008).

41. Spall, M. A. Frontogenesis, subduction, and cross-front exchange at upper ocean fronts. *J. Geophys. Res.* **100**, 2543-2557 (1995).

42. Wong, A., Keeley, R., Carval, T. & Argo Data Management Team. *Argo Quality Control Manual For CTD and Trajectory Data*. 1-58 (Ifremer, 2020).

43. Schmechtig, C. et al. *Processing Bio-Argo particle backscattering at the DAC level Version. Argo Data Management* (Ifremer, 2016). https://doi.org/10.13155/39459

44. Johnson, K. et al. *Processing Bio-Argo nitrate concentration at the DAC Level. Argo Data Management* (Ifremer, 2018). https://doi.org/10.13155/39468

45. Schmechtig, C., Thierry, V. & Team, B. A. Argo Quality Control Manual For Biogeochemical Data. 1-54 (Bio-Argo group, 2016).

46. Xing, X., Briggs, N., Boss, E. & Claustre, H. Improved correction for non-photochemical quenching of in situ chlorophyll fluorescence based on a synchronous irradiance profile. *Opt. Express* **26**, 24734 (2018).

47. Johnson, K. S., Plant, J. N., Riser, S. C. & Gilbert, D. Air oxygen calibration of oxygen optodes on a profiling float array. *J. Atmos. Ocean. Technol.* **32**, 2160-2172 (2015).

48. Johnson, K. S. et al. Biogeochemical sensor performance in the SOCCOM profiling float array. *J. Geophys. Res. Ocean.* **122**, 6416-6436 (2017).

49. Graff, J. R. et al. Analytical phytoplankton carbon measurements spanning diverse ecosystems. *Deep. Res. Part I Oceanogr. Res. Pap.* **102**, 16-25 (2015).

50. Huisman, J., Van Oostveen, P. & Weissing, F. J. Critical depth and critical turbulence: two different mechanisms for the development of phytoplankton blooms. *Limnol. Oceanogr.* **44**, 1781-1787 (1999).

51. Lacour, L. et al. Unexpected winter phytoplankton blooms in the North Atlantic subpolar gyre. *Nat. Geosci.* **10**, 836-839 (2017).

52. Briggs, N. T., Slade, W. H., Boss, E. & Perry, M. J. Method for estimating mean particle size from high-frequency fluctuations in beam attenuation or scattering measurements. *Appl. Opt.* **52**, 6710-6725 (2013).

53. Bohren, C. F. Absorption and Scattering of Light by Small Particles (eds Bohren, C. F. & Huffman, D. R.) (1983). https://doi.org/10.1088/0031-9112/35/3/025

54. Briggs, N., Dall'Olmo, G. & Claustre, H. Major role of particle fragmentation in regulating biological sequestration of CO₂ by the oceans. *Science* **367**, 791-793 (2020).

55. Estapa, M. L., Buesseler, K., Boss, E. & Gerbi, G. Autonomous, high-resolution observations of particle flux in the oligotrophic ocean. *Biogeosciences* **10**, 5517-5531 (2013).
56. Dall'Olmo, G. & Mork, K. A. Carbon export by small particles in the Norwegian Sea. *Geophys. Res. Lett.* **41**, 2921–2927 (2014).

57. Tailleux, R. Spiciness theory revisited, with new views on neutral density, orthogonality, and passiveness. *Ocean Sci.* **17**, 203–219 (2021).

58. Hennon, T. D., Riser, S. C. & Mecking, S. Profiling float-based observations of net respiration beneath the mixed layer. *Global Biogeochem. Cycles* **30**, 920–932 (2016).

59. Schallenberg, C. et al. Diel quenching of Southern Ocean phytoplankton fluorescence is related to iron limitation. *Biogeosciences* **17**, 793–812 (2020).

60. Arrigo, K. R. et al. Photophysiology in two major southern ocean phytoplankton taxa: photosynthesis and growth of phaeocystis antarctica and fragilariopsis cylindrus under different irradiance levels. *Integr. Comp. Biol.* **50**, 950–966 (2010).

61. Serôdio, J. & Lavaud, J. A model for describing the light response of the nonphotochemical quenching of chlorophyll fluorescence. *Photosynth. Res.* **108**, 61–76 (2011).

62. Behrenfeld, M. J. et al. Revaluating ocean warming impacts on global phytoplankton. *Nat. Clim. Chang.* **6**, 323–330 (2016).

63. Organelli, E. et al. A novel near-real-time quality-control procedure for radiometric profiles measured by bio-argo floats: protocols and performances. *J. Atmos. Ocean. Technol.* **33**, 937–951 (2016).

64. Lacour, L., Larouche, R. & Babin, M. In situ evaluation of spaceborne CALIOP lidar measurements of the upper-ocean particle backscattering coefficient. *Opt. Express* **28**, 26989 (2020).

65. Morel, A. et al. Examining the consistency of products derived from various ocean color sensors in open ocean (Case 1) waters in the perspective of a multi-sensor approach. *Remote Sens. Environ.* **111**, 69–88 (2007).

66. Caputi, L. et al. Community-level responses to iron availability in open ocean plankton ecosystems. *Global Biogeochem. Cycles* **33**, 391–419 (2019).

67. Llort, J., Lévy, M., Sallée, J. B. & Tagliabue, A. Onset, intensification, and decline of phytoplankton blooms in the Southern Ocean. *ICES J. Mar. Sci.* **72**, 1971–1984 (2015).

68. Uchida, T. et al. Southern Ocean phytoplankton blooms observed by biogeochemical floats. *J. Geophys. Res. Ocean.* **124**, 7328–7343 (2019).

69. Moreau, S., Boyd, P. W. & Strutton, P. G. Remote assessment of the fate of phytoplankton in the Southern Ocean sea-ice zone. *Nat. Commun.* **11**, 3108 (2020).

70. Mignot, A., Ferrari, R. & Claustre, H. Floats with bio-optical sensors reveal what processes trigger the North Atlantic bloom. *Nat. Commun.* **9**, 190 (2018).

71. Johnson, K. S., Plant, J. N., Dunne, J. P., Talley, L. D. & Sarmiento, J. L. Annual nitrate drawdown observed by SOCCOM profiling floats and the relationship to annual net community production. *J. Geo-phys. Res. Ocean.* **122**, 6668–6683 (2017).

72. Ardyna, M. et al. Delineating environmental control of phytoplankton biomass and phenology in the Southern Ocean. *Geophys. Res. Lett.* **44**, 5016–5024 (2017).

73. Swart, S., Speich, S., Ansorge, I. J. & Lutjeharms, J. R. E. An altimetry-based gravest empirical mode south of Africa: 1. Development and validation. *J. Geophys. Res. Ocean.* **115**, https://doi.org/10.1029/2009JC005299 (2010).

## Acknowledgements

L.L. was supported by a European Union's Horizon 2020 Marie Skłodowska-Curie grant (no. 892653). Visits by L.L. to IMAS/UTas were supported by the Australian Research Council Special Research

Initiative for Antarctic Gateway Partnership (Project ID SR140300001), the Australian Research Council Centre of Excellence for Climate System Science (Project ID B0018492), and the UTas visiting scholar program. J.L. was supported by the European Union's Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No. 754433. N.B. was supported by a European Research Council Consolidator Grant (GOCART, agreement number 724416). P.S. was supported by the Australian Research Council Centre of Excellence for Climate Extremes (Project ID CE170100023). P.W.B. was funded by the Australian Research Council by Laureate Fellowship FL160100131. We thank the Argo Data Management team (ADMT) and the BGC-Argo Data Management team (BGC ADMT). The specific float used in this study was funded by the ARC SRI for Antarctic Gateway Partnership (Project ID SR140300001). We thank the captain and crew of RV Investigator voyage IN2016_V03 and chief scientists Bernadette Sloyan and Susan Wijffels for deploying the float.

## Author contributions

L.L. and J.L. designed the study and conducted the data analysis. L.L. wrote the manuscript. J.L., P.S. and P.W.B. designed the float mission and contributed to manuscript review and editing. N.B. conducted the mean particle diameter analysis and contributed to manuscript review and editing. P.S. and P.W.B. helped to design the study.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-023-36954-7.

**Correspondence** and requests for materials should be addressed to Léo. Lacour.

**Peer review information** *Nature Communications* thanks Lionel Arteaga, Amala Mahadevan, Laure Resplandy and the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer review reports are available.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023