
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

52135
10.1038/s41467-024-52135-6
Article
Antarctic krill sequester similar amounts of carbon to key coastal blue carbon habitats
http://orcid.org/0000-0003-1099-6705
Cavan E. L. e.cavan@imperial.ac.uk

1
http://orcid.org/0000-0002-9758-1096
Mackay N. 2
http://orcid.org/0000-0003-1441-8769
Hill S. L. 3
http://orcid.org/0000-0002-5931-4325
Atkinson A. 4
http://orcid.org/0000-0002-9583-5910
Belcher A. 35
http://orcid.org/0000-0002-1604-7263
Visser A. 6
1 https://ror.org/041kmwe10 grid.7445.2 0000 0001 2113 8111 Imperial College London, Ascot, Berkshire UK
2 https://ror.org/03yghzc09 grid.8391.3 0000 0004 1936 8024 University of Exeter, Exeter, Devon UK
3 https://ror.org/01rhff309 grid.478592.5 0000 0004 0598 3800 British Antarctic Survey, High Cross, Cambridge, UK
4 https://ror.org/05av9mn02 grid.22319.3b 0000 0001 2106 2153 Plymouth Marine Laboratory, Prospect Place, The Hoe, Plymouth, UK
5 https://ror.org/00pggkr55 grid.494924.6 UK Centre for Ecology and Hydrology, Bush Estate, Midlothian, UK
6 https://ror.org/04qtj9h94 grid.5170.3 0000 0001 2181 8870 VKR Centre for Ocean Life, Technical University of Denmark, 2800 Kongens, Lyngby, Denmark
8 9 2024
8 9 2024
2024
15 784228 9 2023
28 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
The carbon sequestration potential of open-ocean pelagic ecosystems is vastly under-reported compared to coastal vegetation ‘blue carbon’ systems. Here we show that just a single pelagic harvested species, Antarctic krill, sequesters a similar amount of carbon through its sinking faecal pellets as marshes, mangroves and seagrass. Due to their massive population biomass, fast-sinking faecal pellets and the modest depths that pellets need to reach to achieve sequestration (mean is 381 m), Antarctic krill faecal pellets sequester 20 MtC per productive season (spring to early Autumn). This is equates USD$ 4 − 46 billion depending on the price of carbon, with krill pellet carbon stored for at least 100 years and with some reaching as far as the North Pacific. Antarctic krill are being impacted by rapid polar climate change and an expanding fishery, thus krill populations and their habitat warrant protection to preserve this valuable carbon sink.

Antarctic krill are important facilitators of carbon sequestration, transferring as much carbon to the deep ocean through their faeces as is stored in seagrass beds and mangrove forests.

Subject terms

Carbon cycle
Marine biology
Ecosystem services
https://doi.org/10.13039/501100000761 Imperial College London issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Marine life has an important role in locking carbon away from the atmosphere in ocean systems, known as ‘blue carbon’. Coastal vegetation such as seagrass, mangroves and salt marshes dominate blue carbon research and policy1,2. Another important form of carbon sequestration which is ubiquitous in the global oceans is through the open-ocean biological carbon pump which stores carbon in the deep-sea3, and at present is not regularly framed in the blue carbon context. The biological carbon pump works by phytoplankton fixing dissolved CO2 into organic carbon during photosynthesis; when the cells die, they may sink to the deep ocean, along with the faecal pellets of phytoplankton grazers (zooplankton) and higher organisms such as fish4. If the sinking carbon is not remineralised, it may be sequestered for decades to centuries5, such that without this biological pump there would be 50 % more CO2 in the atmosphere6,7.

A key player in the biological carbon pump is the Antarctic krill, Euphausia superba, whose faecal pellets, moults and carcasses can completely dominate Southern Ocean carbon fluxes to the ocean’s interior8–14 in the Austral growth season (October through to April, which is our study period). Krill face human-induced threats; through climatic warming reducing sea-ice15 with record Antarctic sea-ice loss this year16, which is an important habitat for larval krill17. An additional issue is the ongoing expansion of the krill fishery in the absence of a management system that accounts for temporal and spatial variability in krill ecology18. Antarctic ecosystem services are globally important but these services are poorly understood in comparison to those elsewhere in the world19. Quantification of metrics that allow global comparisons is an important step towards closing this information gap and ensuring that decision-makers understand the importance of these services. Here we aim to address this by estimating carbon sequestration for the entire circumpolar Antarctic krill population. This advance will allow comparison with coastal vegetation blue carbon systems, showcasing the relevance of open-ocean animals in the same context. The approach used here can also be applied to other pelagic organisms (such as forage fish, jellyfish or gelatinous zooplankton), advancing our understanding of pelagic life in sequestering carbon.

The publication of a circumpolar Antarctic krill density database spanning 90 years (KRILLBASE20) facilitated the first estimates of krill faecal pellet carbon production (FPCprod) and subsequent export from the surface ocean on a large spatial scale. Belcher et al. 10 estimated that in just the marginal ice zone during the productive spring to early autumn months, krill FPCprod and export from the surface ocean (< 100 m) is ~40 MtC yr−1 with export greatest where krill numbers are highest (often in the Atlantic Southern Ocean). However, carbon export is not equivalent to carbon sequestration. In order to become sequestered, exported carbon must avoid remineralisation and sink into a water mass that will not be upwelled or ventilated elsewhere in the ocean for at least 100 years21.

In this study, we use krill density20 multiplied by our assessed faecal pellet egestion rate for Antarctic krill (see Supplementary Note S2 and Supplementary Tables S1–3) to estimate a faecal pellet carbon production flux (FPCprod) in the upper ocean (~20 m). We use this as our ‘export’ flux and attenuate this to the depth the carbon would need to reach to remain in the ocean for at least 100 years. We take sequestration depth data from an ocean circulation model22. This results in our krill faecal pellet carbon sequestration flux (FPCflux), defined as pellet carbon flux at an ocean depth that is sequestered for at least 100 years. To calculate the FPCflux, we attenuate the FPCprod flux following a power-law function known as ‘Martin’s b’, where the ‘b’ attenuation coefficient is −0.3 for krill, a lower attenuation than the global average due to the fast-sinking nature of krill pellets (see “Methods”). We then use the Social Cost of CO2 to convert our krill FPCflux estimates into a single comparable metric: USD ($), and compare it to other, better-known sequestration sources such as mangroves, seagrasses, and salt marshes. Although other life-history traits of Antarctic krill can be important in carbon fluxes, such as their moults and active CO2 transport from migrations, the data and knowledge available to investigate spatial and temporal patterns is greatest for krill faecal pellets, and hence this is our focus. We find Antarctic krill faecal pellets can sequester 20 MtC over the productive Austral spring and summer seasons. This single source of sequestration, from just a single species, is broadly similar to global totals for much more widely documented blue carbon habitats, such as mangroves, salt marshes and seagrasses.

Results

Circumpolar krill faecal pellet carbon sequestration flux

To estimate the circumpolar sequestration of krill faecal pellets, we combined four key elements and here map the data for the month of January, namely a climatological mean krill distribution (Fig. 1a); literature estimates of their faecal egestion to give mixed layer Krill FPCprod fluxes; the rate of attenuation of this flux with depth (here Martin’s b = − 0.3); and the depth needed to achieve sequestration (Fig. 1b). The latter ‘carbon sequestration depth’ is defined as the water column depth where the time-to-surface of the water is at least 100 years, calculated using an ocean circulation inverse model (OCIM)22. The OCIM is a global coarse-resolution (2°) model with 24 vertical levels, for which a steady-state physical circulation has been optimised for consistency with observed passive tracer distributions, sources and sinks, using an adjoint method23. The carbon sequestration depth ranges from 135 to 760 m and is deepest in the Indian Ocean. However, the spatial pattern of estimated carbon sequestration (Fig. 1c) was predominantly driven by the distribution of krill density, with the highest values in the SW Atlantic sector. The mean carbon sequestration depth for the whole data set was 381 m (Table 1), suggesting that particles only need to sink into the mid-mesopelagic zone of the Southern Ocean to remain sequestered for > 100 years. Therefore, some upper ocean observations of sinking particulate organic carbon flux from krill could have observed fluxes that were contributing to carbon sequestration. In all but two of the spatial cells (0.4 % of the total), the length of time carbon is stored increases monotonically with depth. In some of the krill hotspots, such as the Scotia Sea (Atlantic Sector), the carbon sequestration depth is fairly shallow, and a faecal pellet may only need to sink to ~200 m before becoming part of a water mass that would remain away from the surface ocean for 100 years. This highlights how important the choice of depth horizon is in flux observations and models, suggesting that arbitrary depth thresholds might bias estimates of carbon sequestration in the ocean24,25.Fig. 1 Circumpolar krill faecal pellet carbon sequestration flux.

The spatially varying parameters (a) krill density (ind m−2, January mean for data spanning 1926–2016) and (b) carbon sequestration depth, used along with egestion and attenuation rate to give (c) nominal krill faecal pellet carbon sequestration flux (FPCflux), attenuated to the carbon sequestration depth. The resolution is 2° latitude by 6° longitude. Note: the carbon sequestration depths are assumed constant in time (b), but data in (a, c) are for the month of January. See Supplementary Fig. S1 for the circumpolar densities per month for the entire time series. The solid black line represents the location of the Antarctic Polar Front.

Table 1 Mean and ranges of krill faecal pellet production (FPCprod) flux in the mixed layer and mean krill FPCflux at the carbon sequestration depth for all months, October to April, across the entire spatial dataset (All) and by Southern Ocean region (Fig. 1)

	Mean	Total	
Region	FPCprod

(mgC m−2 d−1)

	FPCflux

(mgC m−2 d−1)

	C sequest. depth (m)	C sequestered

(MtC) (billion USD $)

	
Atlantic	14.9 (0–253)	6.8 (0–120)	373 (188–758)	13.7	2.6–32.2	
Indian	7.1 (0–276)	2.9 (0–112)	449 (137–712)	4.1	0.8–9.5	
Pacific	3.6 (0–129)	1.8 (0–64)	313 (137–620)	1.9	0.3–4.4	
All	9.3 (0–276)	4.2 (0–120)	381 (137–758)	19.7	3.7–46.1	
Total carbon sequestered in MtC and USD is also given per region and summed over all regions. The range of carbon sequestered reported in USD ($) is given by applying the value of carbon as $51 per tonne and $640 per tonne (see Methods). The sequestered carbon flux (FPCflux) in each cell is a function of krill density, with density varying spatially and temporally (range = 0–600 ind. m−2) and driving the ranges in carbon fluxes.

Nominal krill FPCprod flux at 20 m depth ranged from 0 – 276 mgC m−2 d−1 with the highest mean flux in the Atlantic sector of the Southern Ocean (and Table 1), as this is where Antarctic krill are most densely populated (Fig. 1a). These krill pellet export fluxes are similar in magnitude to observed zooplankton and krill pellet export rates in the Southern Ocean11. We use a lower egestion rate of 0.46 mgC krill−1 d−1 compared to earlier studies such as Belcher et al. 10, who used a rate of 3.2 mgC d−1 and hence our range is an order of magnitude smaller (~300 compared to ~ 3000 mgC m−2 d−1). Our revised egestion rate estimate is derived from a number of independent methods (Supplementary Tables S1–S3) and is intended to reflect the range of feeding conditions experienced between spring to autumn and at the circumpolar scale, whereas the higher values in the literature represent short observation periods during a time of high food availability26.

The krill pellet export fluxes at 20 m were attenuated to simulate the degradation of pellets with depth due to consumption and bacterial remineralisation using a power-law coefficient27 based on observations from krill pellet fluxes presented in Belcher et al. 13 and including more recent krill pellet fluxes from Pauli et al. 28 (see Supplementary Table S4). This results in a krill pellet attenuation rate (Martin’s b) of −0.30, which shunts more carbon to depth than the ‘global’ value of b of −0.86 commonly applied to all sinking POC types. The lower attenuation of krill pellets results mainly from the rapid sinking of these pellets9,29. In addition, our attenuation rate is an average derived from observations which compare fluxes at mixed layer and mesopelagic depths and are likely affected by additional egestion between these two depths11. In the absence of detailed information on the depth distribution of egestion, this value appropriately compensates for the simplifying assumption in our calculation that faecal pellet carbon is only attenuated between 20 m and the carbon sequestration depth, rather than also being produced.

Sequestered krill faecal pellet carbon

Krill FPCflux at the sequestration depth ranged from 0 – 120 mgC m−2 d−1 (Fig. 1c), remaining highest on average in the Atlantic sector (Fig. 1c and Table 1). The total carbon that could be sequestered each Austral spring/summer was estimated by summing the pellet carbon fluxes (FPCflux, Fig. 1c) at the carbon sequestration depth over the area (m−2, Fig.1) and the number of days in the months (October to April) sampled. Total nominal krill FPCprod from within the mixed layer was 44 MtC yr−1, with 20 MtC yr−1 (45 %) sequestered (Fig. 2). This sequestered carbon equates to a range in monetary value of USD$4 – 46 billion (Table 1, Fig. 2), depending on the value of the Social Cost of CO2 per tonne of CO2 used.Fig. 2 Total carbon sequestered by krill faecal pellets.

Total carbon sequestered in Mt C (million tonnes) and billion USD per sector of the Southern Ocean. In this graphical representation, an economic value ($135 per tonne of CO2 sequestered) at the lower end of the considered range ($51–640) is used. The total across all three Southern Ocean regions is 20 Mt C and USD $4–46 billion using the full range of the Social Cost of CO2. See Table 1 for the ranges in fluxes that lead to the total carbon sequestered presented here and Methods for details on uncertainty and error.

Following the trend with the krill densities and pellet fluxes, most of the carbon is injected into the interior ocean in the Atlantic sector of the Southern Ocean, where the carbon sequestration depth is fairly shallow. Using spatially-resolved net primary productivity (NPP) values for the Southern Ocean30, we find that on average ~2.5 % of NPP is routed to krill pellet carbon sequestration (i.e., FPCflux/NPP) across the Southern Ocean, with values being up to 74 % in swarm locations (see Supplementary Table S5). The ‘transient time’,i.e., the time from photosynthesis by phytoplankton capturing dissolved inorganic carbon to krill pellet carbon being sequestered is likely less than a week, considering gut passage times of less than a day and pellet sinking rates of ~300 m d−129. This shows how, in particular, krill swarms are extremely efficient vectors of carbon sequestration, quickly grazing down phytoplankton and the subsequent mass production of fast-sinking faecal pellets10. The proportion of krill FPCprod that is sequestered (i.e., the transfer efficiency) is much higher (45 %) than the commonly reported 1–10 % in the biological pump because of the low attenuation rate of krill pellets and because pellets only need to sink a few hundred meters to achieve sequestration, compared to the 1000 m or 2000 m depths often used to define ‘transfer efficiency’22,31. This contributes to the growing body of work moving biological carbon pump research on from calculating carbon fluxes at arbitrary discrete depths, to factoring in carbon sequestration timescales25.

To determine the additional contribution of sequestered krill pellet carbon flux relative to plankton (phytoplankton and copepods) typically represented in biogeochemical and biological pump models, we use a published food-web model output32 (updated from ref. 22) to estimate the total carbon sequestered by plankton in our study region. The model by Nowicki et al. 32 estimates that circumpolar phytoplankton and copepods sequester 160 MtC compared to the 20 MtC we estimate for Antarctic krill, meaning that krill faecal pellets account for 12 % of the total ‘plankton’ carbon sequestration in the Southern Ocean. This 12 % likely represents carbon sequestration incorrectly allocated to phytoplankton detritus or copepod faecal pellets in models that do not explicitly represent such large swarming organisms as krill. Therefore, our type of model is important for conservation purposes as it can highlight the locations of high pulses of carbon to the ocean interior (see Supplementary Fig. S2a) from organisms of interest, here krill. When we include our best estimates of the additional carbon sequestration potential from krill moults (20 MtC) and active transport (26 MtC) (see Methods for details), then krill could sequester a total of 66 Mt Cy−1. The current generation of Earth System Models (ESMs) underestimate or entirely omit the contribution of these larger organisms to the ocean carbon cycle. Incorporating micronekton, such as krill, and behaviour, such as diel vertical migration, should be priorities for biological parameter development in ESMs.

Validating our approach

Our total pellet carbon sequestration estimates (Fig. 2) are based on detailed information on the distribution and density (number of krill individuals per unit area) of krill15, information on carbon sequestration depth and a revised set of parameters to represent egestion rate and the attenuation rate (Martin’s b) of sinking faecal pellets. Our values for each parameter are within the range of values reported in the literature (see “Methods” and Supplementary Tables S1–S4). Our revised egestion rate is more conservative (i.e., results in lower per capita krill FPCprod fluxes than previously used10), whereas our attenuation rate is based on krill pellet flux observations10, and because krill pellets sinking rapidly, the attenuation is more gradual than the value typically used to model global total POC flux (i.e., it preserves more of the carbon to sequestration depth). Our circumpolar krill density (5.7 × 1014 krill individuals in 16 million km2 of Southern Ocean) is similar to a previously reported circumpolar krill density of 5.4 × 1014 krill33. We compare our model with krill pellet carbon observations in our study area to assess whether our methodology is appropriate and representative of krill pellet fluxes.

Our model compares well with direct observations of krill faecal pellet carbon and total particulate organic carbon (POC) fluxes measured in the Southern Ocean. First we compare to South Georgia (53°S 38°W) sediment trap data at 300 m depth13. Extracting the krill faecal pellet flux in January from the annual sediment trap time series in Manno et al. 13 results in krill pellet flux at 300 m of 46 mgC m−2 d−1. We attenuated krill FPCprod fluxes from our model to 300 m at the location of the South Georgia sediment trap, which gave a similar krill FP flux of 40 mgC m−2 d−1. We also compare our model to observations of total POC fluxes further south at the Palmer Station (64°S 66°W) on the Western Antarctic Peninsula14. Here the mean total POC flux (pellets plus other sinking particles) in January from a 20-year time-series sediment trap at 170 m depth was 16 mgC m−2 d−1. Our modelled krill pellet flux at this location and depth is 6 mgC m−2 d−1, suggesting that 38 % of the total POC flux is made up of adult krill pellets near the sea ice (Supplementary Fig. S2b). Krill pellets tend to dominate the Palmer Station sediment trap in Austral summer, contributing up to 86 % of the total observed POC flux12. Our value of 38 % is lower than this because our model only accounts for adult krill > 40 mm, and close to the sea-ice and continent, there is a much larger population of juvenile krill, compared to, for instance, South Georgia15. These results give confidence that our large-scale estimates of krill carbon fluxes are reasonable and within the expected ranges based on observations.

Uncertainty in our estimates

We ran a sensitivity analysis, varying each parameter in turn by ± 10% (Fig. 3a, Supplementary Table S6) to identify the relative influence of each on our estimate of total carbon sequestered (as in Fig. 2). These parameters were krill density, egestion rate, attenuation rate and sequestration depth. Krill density and egestion rate have the joint largest influence. These two parameters together provide the krill FPCprod flux which sets the amount of carbon that can be injected into the deep ocean. As krill density is the main forcing data for our model and it has a large range over orders of magnitude (0–600 ind. m−2), krill distribution sets the spatial spread of the carbon fluxes (Fig. 1). The carbon sequestration depth is the other forcing data used in our model but had the least influence in our sensitivity analysis. Whilst it does vary regionally (Fig. 1b) the changes in sequestration depth are more gradual and less pronounced, hence in terms of forcing data krill density has the largest impact on carbon sequestration.Fig. 3 Sensitivity analysis.

a parameter sensitivity analysis changing each parameter by 10 %, and (b) uncertainty analysis where parameter values are changed based on observations reported in the literature. The circular gridlines give the percentage change in total carbon sequestered, with 100 % equivalent of 20 Mt C. Parameters which set the rate of export of pellet carbon out of the surface ocean, krill density and egestion rate, have the largest effect on overall carbon sequestered (Fig. 3a). However, the largest natural variability and/or variation in observations to date is associated with the attenuation of pellets (b) with depth (Fig. 3b), even though the 25th and 75th percentiles are represented for attenuation, whereas more extreme values (5th and 95th percentile) are used for egestion and sequestration depth (see “Methods”).

We conducted an additional analysis, varying each parameter based on the range of plausible values suggested by our analysis (for krill density, egestion rate and carbon sequestration depth) or observations reported in the literature (for Martin’s b). This allows us to identify where the scientific understanding is least certain and/or where natural variability is highest, and more research is needed to further constrain these parameter estimates. Here, the greatest influence on total carbon sequestered is the huge range of Martin’s b observations (−2.46 to 1.81) (Fig.  3band Supplementary Table S4). This range reflects both the spatial and temporal variability in attenuation of pellets, limitations with measurements which can arise from the time delay between pellets leaving the surface and reaching the mesopelagic34, and egestion beneath the assumed export depth8,35. We consider b values of − 0.61 and + 0.13, which are the 25th and 75th percentile of values from the literature (see Supplementary Table S4) and result in total carbon sequestered of 9 and 63 MtC, 44 % or 320 % of our best estimate (20 MtC) respectively (Fig. 3b and Supplementary Table S6). The total carbon sequestered with b set to the ‘global’ accepted value of − 0.86 is 4 Mt C. The next largest effect, increased egestion rates, results in a range of total carbon sequestration values from 24 to 267 % of our best estimate (Fig. 3). To improve modelling studies such as ours, we need to understand what factors control the pellet attenuation rate, with candidates including krill diet, composition of mesopelagic consumers and water temperatures. This will enable spatial and temporally varying attenuation rates to be applied.

Tracking krill carbon through the oceans

Carbon transferred to the deep by Antarctic krill does not necessarily stay in the Southern Ocean. As bacteria break down the carbon-rich pellets, they produce dissolved inorganic carbon, which can be transported around the global oceans. To track the fate of krill pellet carbon after it has been injected into the ocean’s interior, we ran an additional model analysis using an Ocean Circulation Inverse Model (OCIM; see “methods”)22. This allows us to determine the steady-state interior distribution of dissolved inorganic carbon originating from krill faecal pellets, i.e., the total amount of carbon stored from krill pellets in the ocean at any one time. Most of the krill carbon remains in the Atlantic sector (Fig. 4), where krill pellet carbon fluxes are highest (Figs. 1, 2); meanwhile, some are transported to the Pacific and Indian Oceans, including in the Northern Hemisphere with very little reaching the North Atlantic Ocean. Once krill pellet carbon leaves the surface ocean in the form of krill faecal pellet fluxes, it remains there with an average residence time of 219 years before it comes into contact with the surface again. Based on our krill pellet export fluxes and attenuation rate (b = − 0.30) the total amount of carbon derived from Antarctic krill pellets that the oceans can hold is estimated as 8.7 Gt C (Fig. 4). This attenuation rate results in ~ 20 % of the export carbon reaching the benthos as particulate pellets, where additional carbon sequestration mechanisms such as burial and incorporation into the bodies of long-living invertebrates can occur36. The total krill-carbon stored in the global water column reduces to 1.7 Gt C for only 58 years if the non-krill pellet attenuation rate of b = − 0.86 is applied instead (Supplementary Figs. S3, S4). This emphasises how important the remineralisation rate, or the attenuation rate, is in estimating long-term carbon storage by sinking particles. Given the evidence that for krill pellets, the carbon is shunted deeper into the ocean than other types of particles (Supplementary Table S4), the b parameter should vary with particle type in model simulations.Fig. 4 Global pelagic carbon footprint from Antarctic krill faecal pellets.

Equilibrium dissolved inorganic carbon (DIC) concentrations (gC m−2) from the remineralisation (attenuation) of Antarctic krill faecal pellets with depth, assuming a remineralisation coefficient of b = − 0.30. The highest water-column concentrations exist in the Atlantic Southern Ocean, where krill pellet fluxes are highest (Figs. 1, 2). This map does not include the remaining pellets that escape remineralisation and fuel Antarctic communities on the seafloor.

Discussion

Krill in the context of blue carbon

Defining sinking krill pellet fluxes in the context of the physical circulation allows us to quantify the carbon sequestered (defined here as locked away for >100 years), and thus discuss them in the context of other blue carbon stores. Our estimated sequestration store of Southern Ocean krill faecal pellets (20 MtC) over an Austral productive season is of the same magnitude as other global blue carbon stores, with saltmarsh, mangroves and seagrass estimated to sequester 13, 24 and 44 MtC yr−11, respectively (Fig. 5). The krill-pellet carbon injected into the oceans is three orders of magnitude lower in concentration (1.3 tC km−2 yr−1) than for coastal vegetative blue carbon stores that predominantly store carbon in highly concentrated soils and sediments (e.g., seagrass = 156 tC km2 yr−1)1 (Fig. 5). It is the vast Southern Ocean habitat that Antarctic krill occupy (16 million km2 – our study – to 19 million km233) compared to coastal vegetation (e.g., seagrass = 0.32 million km2 globally) that makes the total carbon sequestered by krill equal to that of the coastal blue carbon stores. Given our conservative approach to faecal pellet export estimation, the significant quantities of carbon in living krill37 and the additional flux that will result from krill moults (~20 MtC), winter feeding and respiration during migrations (~ 26 MtC)4,12,13 (Fig. 5) it is likely that Antarctic krill is amongst the world’s most important carbon-storing organisms. Given the comparability of these numbers associated with Antarctic krill to other carbon sequestration mechanisms, we suggest our approach can be used to quantify the carbon sequestration value of other pelagic marine life, such as copepods, gelatinous zooplankton and potentially low-mid trophic level nekton such as forage fish.Fig. 5 Comparison of krill pellet blue carbon with coastal vegetation blue carbon stores.

The krill pellet carbon sequestered in MtC per October-April growing season is given and converted to per km2. Krill sequester less per surface area of the ocean than coastal vegetation blue carbon, but the vast area krill inhabit puts the pellet carbon sequestered in the same order of magnitude as the coastal vegetation stores, according to data in Bertram et al. We also include best estimates of the contribution of sinking moults and daily migration (active transfer of CO2). If these are correct, krill could sequester a similar amount of carbon as seagrasses do globally. There is insufficient data to yet do this analysis for krill carcasses, but we include them graphically as they are an important carbon sink. Note, that other krill habitat size estimates find krill habitat to be 19 km2 using a different grid resolution.

Conservation perspectives

The krill-dominated habitats of the Southern Ocean are akin to the better-documented salt marsh, mangroves and seagrass habitats: not only do they support unique, valuable and iconic species, but they are also key conduits for carbon storage. Both of these facets warrant the utmost importance of conserving into the future. The value for krill is in the order of billions of dollars in terms of carbon storage that benefits society as an ecosystem service. Krill are under threat from warming resulting in loss of habitat, potential changes to food availability and competition, and from the krill fishery15,19. Whilst there are some fishery restrictions and protected areas in place in the Southern Ocean, more needs to be done to ensure this carbon sink can remain active into the future. It is not clear how a different upper ocean community (i.e., with less krill) may look, and the implications this has for carbon sequestration. Though other low-trophic level pelagic organisms are increasing in numbers in the Southern Ocean (e.g., salps)38, their efficiency in transferring carbon deep enough for sequestration is not always as high as krill. For instance, salp pellets do not always penetrate into the deeper ocean in such high numbers due to their fragility28,39, compared to dense krill faecal pellets. These subtleties illustrate the complexity of the biological carbon pump.

Antarctic krill are clearly important vectors of carbon sequestration. The large uncertainty in our understanding of future pelagic communities and carbon export under climate change means we cannot predict with certainty how future changes in the krill population will affect the magnitude of carbon sequestration in the Southern Ocean40. Our study highlights the importance of Antarctic krill, and likely other components of the biological pump, in sequestering carbon which is on par with coastal blue carbon stores. This study and others14,41,42 are increasingly showing the importance of krill and other zooplankton in structuring Southern Ocean food webs and activating a strong biological pump. Understanding the relative roles of the key phyto-, zooplankton and micronekton functional groups in carbon sequestration is therefore an urgent priority. Valuing pelagic vulnerable marine ecosystems in terms of carbon sequestration, like that of Antarctic krill, emphasises how crucial it is to meet climate goals and work towards including carbon policies in resource management.

Methods

Krill density

Krill density data were obtained from KRILLBASE20, which compiles postlarval krill numerical densities (number per unit sea surface area) based on net sampling around Antarctica spanning 1926–2016. These data were filtered following previously published procedures15, where only data collected from ‘hauls’ or ‘stratified pooled hauls’ where the top sampling depth was shallower than 20 m, and the bottom sampling depth deeper than 50 m were used. The data set contains a column of ‘standardised krill under 1m2’, which represents krill density normalised to the first of January, using efficient sampling gear during the night-time (when krill are likely to be in the surface layers). This standardisation uses empirical algorithms to extrapolate what the krill density would be on the 1st of January, depending on when the sample was collected (Table 4 of Atkinson et al. 20, and Supplementary Equations S1–S3). This allows for almost full circumpolar coverage of krill density for the month of January. The same model can, therefore, be used to calculate what the krill density would be on the first date of each month (Supplementary Fig. S5), resulting in krill density estimates for the full Austral spring-summer season, from October to April (Supplementary Fig. S1). See Supplementary Note S1 and Supplementary Equations S1–S3 for more information.

The data were then projected onto a 2° × 6° lat/lon grid, by averaging (mean) krill densities (N¯) in each month (t) over this spatial scale. This is more conservative than a 3 × 9° grid (Supplementary Fig. S6) and in line with a resolution also used by Atkinson et al. 33. It balances the risk of inappropriate extrapolation (when data are only available for a small fraction of a grided square), with the risk of underrepresenting the extent of krill habitat (as would occur with the 2° × 2° resolution of OCIM) (Supplementary Fig. S7). Where krill densities were > 600 ind. m−2 (n = 7, 0.3 % of data, ranged from 673 – 1681 ind. m−2) these were capped at 600 ind. m−2. This excludes bias from a few extremely high net catches sampling in a dense swarm but still allows for high krill densities to occur, representing regions where swarms are commonly found. The 7 data capped (one for each month) were in fact, from just one net haul, and removing data from this observation reduced the estimated total carbon sequestered from 21.5 MtC yr−1 to our value reported of 20 MtC yr−1.

Faecal pellet flux

For each month (t), krill faecal pellet carbon production flux (FPCprod) at 20 m for each 2 × 6° cell (k) was calculated by multiplying the monthly mean krill density (N¯, # m−2) by a krill faecal pellet egestion rate (E) (Eq. 1).1 FPCprodk,tmgCm−2d−1=N¯k,t⋅E

High daily rates determined from laboratory studies used in some previous studies (e.g., 4.03 mg C ind−1 d−1 (ref. 26) and 3.2 mg C ind−1 d−1 (ref. 10)) are unlikely to be sustained throughout the summer season or at the circumpolar scale. Assuming higher egestion rates, as in these previous studies, would have commensurate effects on the estimated magnitude and value of krill-mediated carbon sequestration, as shown in our sensitivity analysis. We, therefore, used three different methods to estimate egestion rates, and we assessed their validity by converting them to estimates of annual food consumption (in carbon units) by the circumpolar krill stock and comparing this to an indicative estimate of circumpolar primary production (1949 Mt C y−1)43. The resulting median egestion rate for an adult krill of 40 mm length was 0.46 mg C ind−1 d−1, with a 5-95th percentile range of 0.11 to 1.23 mg C ind−1 d−1. The results of this analysis and more detailed methods are available in Supplementary Note S2 and Supplementary Tables S1–S3.

The mesopelagic krill FP carbon flux (FPCflux) was calculated using Eq. 2, which attenuates the krill FPCprod between 20 m (z0) and any mesopelagic depth (z) using Martin’s b value of − 0.30, based on data collated by Belcher et al. 10 with the addition of Pauli et al. 28. As the attenuation rate is a crucial factor in determining pelagic sequestration, we also repeated the calculation using the globally applied estimate of − 0.86, increasing the remineralisation of the FPs as they sink to give a more conservative estimate. As Martin’s b describes the rate of transfer of POC flux to depth, it implicitly represents the sinking rate, remineralisation rate and zooplankton fragmentation rates of sinking faecal pellets.2 FPfluxk,t,zmgCm−2d−1=FPCprodk,t⋅zz0b

Carbon sequestration depth

The main purpose of this study was to calculate the amount of FP flux that would remain sequestered in the oceans for at least 100 years, and thus ocean circulation must be considered. For this, we used the Ocean Circulation Inverse Model (OCIM) model output product from DeVries & Weber22, which outputs the time-to-surface of a parcel of water throughout the water column, known as the First Passage Time (FPT). This allows us to find the depth for each cell at which the FPT would be at least 100 years, and to conclude that if the FP sank below this depth, then the carbon, either as particulate organic or dissolved inorganic carbon (POC or DIC) would remain away from the surface of the ocean for 100 years or more and is therefore sequestered.

We calculated the depth where FPT = 100 years or more, and this became our sequestration depth or FPT100. In the OCIM model, the sequestration times (years) are based in part on a winter mixed layer depth, as anything sinking shallower than that would be ventilated within a year. FPT in the ocean does not increase monotonically with depth due to overlapping water masses. However, for the OCIM outputs used in this study, FPT was found to increase monotonically with depth in all but two horizontal grid cells. This result is somewhat unexpected, given, for example, the structure of water masses in the Southern Ocean, where upwelling Circumpolar Deep Water sits underneath downwelling mode and intermediate waters, potentially highlighting weaknesses in the coarse resolution OCIM’s representation of the circulation. If the analysis were carried out using a finer resolution physical model, the final sequestration could decrease in some places or increase in others; unfortunately, no finer resolution version of the OCIM currently exists, and analysis of the sequestration of krill carbon using a forward model is outside the scope of this study.

To determine the krill FP carbon sequestration (FPCflux), we applied Martin’s b curve again, but this time only for the FPT100 depth, which varied in space (k) but not in time:3 FPCfluxk,tmgCm−2d−1=FPCprodk,t⋅FPT100kz0b

The final total circumpolar carbon sequestered by krill FPs (Table 1 and Fig. 2) was found by multiplying FPCflux by the area of each cell (Ak, m−2) and the number of days in the month (mk), and summed across the whole sampling area (k) and months analysed (t) (October to April):4 FPseqtotGtC=∑k,tFPCfluxk,t⋅Ak⋅mt⋅1e−18

Error bars are not presented in Fig. 2 as the final data here are sums, not means, and the main forcing data of krill densities are not normally distributed. Instead in Table 1 we report the ranges of fluxes present in our data, which lead to the total summed carbon stored. The statistical error can only be calculated based on krill density means and the associated error, such as standard deviation (20 ± 62 ind. m−2). The variation in the mean is large as the range in krill density is large (0–600 ind. m−2), with a median krill density of 1 ind. m−2. Calculating uncertainty in our carbon sequestration results based on krill densities leads to negative carbon sequestration values. We, therefore, report the ranges of fluxes in Table 1 and run sensitivity and uncertainty analyses (see below and Fig. 3).

Sensitivity analysis

The relative influence of each parameter used to estimate the final carbon sequestered (GtC) was investigated by sequentially increasing or decreasing each parameter by 10 % and investigating the impact on the final sequestered carbon at the sequestration (FPT100) depth for krill faecal pellets. A second analysis varied each parameter sequentially by either the 5th and 95th percentile values of those in our data or analysis (egestion rate and sequestration depth). For the attenuation rate (Martin’s b), we compiled observations of krill pellet fluxes from the literature but restricted the percentiles to the 25th and 75th in our uncertainty analysis as the effect of changing attenuation was large (Fig. 3b). Finally, for krill densities as the data were so skewed and over a large range, we instead used as the upper limit the arithmetic mean of all KRILLBASE records reported by Atkinson33 of 36 ind. m−2 and centred our mean (20 ind. m−2) so the lower limit was 4 (i.e., ± 16 ind. m−2).

Economic value

The total sequestered krill faecal pellet carbon (Eq. 4) was then converted from mass in Gt to its worth to society in monetary terms by multiplying by the Social Cost of CO2 (SCCO2, Eq. 5). At the time of writing the US was still using their interim SCCO2 of only USD$51 per tCO2, but a 2021 study found the global average SCCO2 was USD$640 per tCO21. However, a 2018 IPCC report predicts values could range from USD$135 – 5500 per tCO2 by 2030 to reach only 1.5 °C of warming44. Here we use a range of SSCO2 from USD$51 – 640 per tCO2. To convert from the monetary value of CO2 to the monetary value of organic carbon, the SSCO2 was multiplied by the ratio of total mass to carbon in a CO2 molecule:5 FPseq$=FPseqtot⋅SSCO2⋅4412

Comparison to other blue carbon mechanisms

We compare the sequestered carbon from krill faecal pellets to coastal vegetation blue carbon stores. As krill faecal pellets are only one element of the contribution of krill to carbon sequestration, we also considered whether we could estimate the contribution from krill moults, carcasses and migrations, both daily and seasonally. We are confident in our ability to make the best estimates on the magnitude of moults and migrations. Firstly, Manno et al. 13 found that moult carbon contribution at 300 m equalled that of krill pellets in the Southern Ocean. We parsimoniously assume the sequestration by moults is equal to that of pellets. For migration, given that ~20 % of krill live deeper than 400 m45, we estimate that the respired DIC released from these deep-dwelling krill, either from daily or seasonal migrations, is 26 MtC yr−1. Please refer to the Supplementary Note S3 for more details. The values for salt marshes, seagrass and mangroves in terms of Mt C stored each year were taken from Bertram et al. 1. Finally, the total biological carbon pump in our study region of the Southern Ocean was calculated using the OCIM model output22, by taking POC flux at the sequestration depth (FPT100) and summing over the area and time (October to April) of our sampling period.

Tracing respired krill pellet carbon globally

We used the OCIM transport matrix to track carbon respired from krill pellets as they sink through the ocean. Summed krill pellet export fluxes (FPCexp) over the spring to early autumn period and the krill attenuation rate of b = − 0.30 were used to find the steady-state distribution in the ocean of respired carbon originating from the pellets. Further, the global integral of this distribution as well as the average residence time (in years) was estimated. Note this does not include particulate carbon that escapes remineralisation and reaches the seafloor. This method has already been published in Boyd et al. 21 and Nowicki32. We also ran this analysis for a higher attenuation rate of b = − 0.86 (Supplementary Fig. S4). The matrix equation is:6 dC/dt=C(A-S)+Q=0

where A is the transport matrix, C(x,y,z) [gC m-3] the concentration of dissolved inorganic carbon (DIC) within the ocean, Q(x, y, z) [gC m−3 yr−1] is the source of DIC (derived from observed surface pellet fluxes and a Martin curve), and S [yr−1] represents the instantaneous removal of C from the surface. The steady-state solution: that is7 C*=−Q/(A-S)

is an estimate of the distribution of DIC emitted by Q in the world’s oceans before it comes into contact with the atmosphere again. Integrating this over the entire ocean gives the total sequestered DIC. Specifically, writing V(x, y, z) [m3] as the volume of each grid box in the transport matrix representation, Cseq = ∑ C*V. Dividing this by the net rate at which DIC is injected into the ocean (i.e., the volume integral of Q) gives the mean residence time for DIC injected via Q below the surface.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Supplementary information

Supplementary Information

Peer Review File

Reporting Summary

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52135-6.

Acknowledgements

E.L.C. was supported by an Imperial College Research Fellowship, a Natural Environment Research Council (NERC) grant NE/Y004515/1 and a WWF research grant (GB085708) along with A.A. and S.L.H. S.L.H. was additionally supported by the NERC British Antarctic Survey (BAS) ALI-Science Southern Ocean Ecosystems project. AA’s contribution was also supported by the NERC PICCOLO programme. N.M. was funded by NERC grant NE/W001543/1. A.W.V. is supported by the Centre for Ocean Life, a VKR Centre of excellence funded by the Villum Foundation, and by the Horizon2020 projects ECOTIP (grant agreement #869383) and OCEAN-ICU (grant agreement #101083922). Figure 5 was produced by Visual Knowledge.

Author contributions

E.L.C. devised the research and led the analysis and production of the figures and data presented in the main text. NM provided the carbon sequestration depth from the OCIM model. S.L.H. calculated the pellet egestion rate via the three different approaches. A.A. estimated the carbon sequestered via krill respiration. S.L.H. and A.A. provided expertise and advice on the circumpolar krill densities/abundance and the spatial area of krill habitat. A.B. and E.L.C. revised the attenuation rates. A.W.V. ran the transport matrix model to trace krill carbon in the global oceans. E.L.C. led the writing, and all authors contributed, with significant contributions from S.L.H.

Peer review

Peer review information

Nature Communications thanks Rebecca Trinh, Thomas Weber and the other anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

The two forcing data sets used in this study can be accessed online. KRILLBASE density data can be downloaded from the British Antarctic Survey at https://www.bas.ac.uk/project/krillbase/#data. The OCIM model and data can be downloaded at https://tdevries.eri.ucsb.edu/models-and-data-products/.

Code availability

The code used in this study for analysis and plotting is available at 10.6084/m9.figshare.26511022.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Bertram C The blue carbon wealth of nations Nat. Clim. Chang. 2021 11 704 709 10.1038/s41558-021-01089-4
Bertram, C. et al. The blue carbon wealth of nations. Nat. Clim. Chang. 11, 704–709 (2021).10.1038/s41558-021-01089-4
2. Macreadie PI The future of Blue Carbon science Nat. Commun. 2019 10 3998 10.1038/s41467-019-11693-w 31488846
Macreadie, P. I. et al. The future of Blue Carbon science. Nat. Commun. 10, 3998 (2019).31488846 10.1038/s41467-019-11693-w
3. Volk, T. & Hoffert, M. Ocean carbon pumps: Analysis of relative strengths and efficiencies in Ocean-driven atmospheric CO2 changes. in The Carbon Cycle and Atmospheric CO2: Natural Variations Archean to Present 99–110 (American Geophysical Union, 1985).
4. Turner JT Zooplankton fecal pellets, marine snow, phytodetritus and the ocean’s biological pump Prog. Oceanogr. 2015 130 205 248 10.1016/j.pocean.2014.08.005
Turner, J. T. Zooplankton fecal pellets, marine snow, phytodetritus and the ocean’s biological pump. Prog. Oceanogr. 130, 205–248 (2015).10.1016/j.pocean.2014.08.005
5. Cavan, E. L., Laurenceau-Cornec, E. C., Bressac, M. & Boyd, P. W. Exploring the ecology of the mesopelagic biological pump. Prog. Oceanogr. 176, 102125 (2019).
6. Kwon E Primeau F Sarmiento J The impact of remineralization depth on the air-sea carbon balance Nat. Geosci. 2009 2 630 635 10.1038/ngeo612
Kwon, E., Primeau, F. & Sarmiento, J. The impact of remineralization depth on the air-sea carbon balance. Nat. Geosci. 2, 630–635 (2009).10.1038/ngeo612
7. Parekh, P., Dutkiewicz, S., Follows, M. J. & Ito, T. Atmospheric carbon dioxide in a less dusty world. Geophys. Res. Lett. 33, 10.1029/2005GL025098 (2006).
8. Cavan EL The importance of Antarctic krill in biogeochemical cycles Nat. Commun. 2019 10 1 13 30602773
Cavan, E. L. et al. The importance of Antarctic krill in biogeochemical cycles. Nat. Commun. 10, 1–13 (2019).30602773
9. Belcher A The potential role of Antarctic krill faecal pellets in efficient carbon export at the marginal ice zone of the South Orkney Islands in spring Polar Biol. 2017 40 2001 2013 10.1007/s00300-017-2118-z 32009725
Belcher, A. et al. The potential role of Antarctic krill faecal pellets in efficient carbon export at the marginal ice zone of the South Orkney Islands in spring. Polar Biol. 40, 2001–2013 (2017).32009725 10.1007/s00300-017-2118-z
10. Belcher, A. et al. Krill faecal pellets drive hidden pulses of particulate organic carbon in the marginal ice zone. Nat. Commun. 10, 10.1038/2Fs41467-019-08847-1 (2019).
11. Cavan EL Attenuation of particulate organic carbon flux in the Scotia Sea, Southern Ocean, is controlled by zooplankton fecal pellets Geophys. Res. Lett. 2015 42 821 830 10.1002/2014GL062744
Cavan, E. L. et al. Attenuation of particulate organic carbon flux in the Scotia Sea, Southern Ocean, is controlled by zooplankton fecal pellets. Geophys. Res. Lett. 42, 821–830 (2015).10.1002/2014GL062744
12. Gleiber MR Steinberg DK Ducklow HW Time series of vertical flux of zooplankton fecal pellets on the continental shelf of the western Antarctic Peninsula Mar. Ecol. Prog. Ser. 2012 471 23 36 10.3354/meps10021
Gleiber, M. R., Steinberg, D. K. & Ducklow, H. W. Time series of vertical flux of zooplankton fecal pellets on the continental shelf of the western Antarctic Peninsula. Mar. Ecol. Prog. Ser. 471, 23–36 (2012).10.3354/meps10021
13. Manno, C., Fielding, S., Stowasser, G., Murphy, E. J. & Thorpe, S. E. Continuous moulting by Antarctic krill drives major. Nat. Commun. 11, 6051 (2020).
14. Trinh R Ducklow HW Steinberg DK Fraser WR Krill body size drives particulate organic carbon export in West Antarctica Nature 2023 618 526 530 10.1038/s41586-023-06041-4 37316721
Trinh, R., Ducklow, H. W., Steinberg, D. K. & Fraser, W. R. Krill body size drives particulate organic carbon export in West Antarctica. Nature 618, 526–530 (2023).37316721 10.1038/s41586-023-06041-4
15. Atkinson, A. et al. Krill (Euphausia superba) distribution contracts southward during rapid regional warming. Nat. Clim. Chang. 9, 142–147 (2019).
16. Purich A Doddridge EW Record low Antarctic sea ice coverage indicates a new sea ice state Commun. Earth Environ. 2023 4 314 10.1038/s43247-023-00961-9
Purich, A. & Doddridge, E. W. Record low Antarctic sea ice coverage indicates a new sea ice state. Commun. Earth Environ. 4, 314 (2023).10.1038/s43247-023-00961-9
17. Siegel V Loeb V Recruitment of Antarctic krill Euphausia superba and possible causes for its variability Mar. Ecol. Prog. Ser. 1995 123 45 56 10.3354/meps123045
Siegel, V. & Loeb, V. Recruitment of Antarctic krill Euphausia superba and possible causes for its variability. Mar. Ecol. Prog. Ser. 123, 45–56 (1995).10.3354/meps123045
18. Meyer B Successful ecosystem-based management of Antarctic krill should address uncertainties in krill recruitment, behaviour and ecological adaptation Commun. Earth Environ. 2020 1 1 12 10.1038/s43247-020-00026-1
Meyer, B. et al. Successful ecosystem-based management of Antarctic krill should address uncertainties in krill recruitment, behaviour and ecological adaptation. Commun. Earth Environ. 1, 1–12 (2020).10.1038/s43247-020-00026-1
19. Cavanagh, R. D. et al. Future risk for Southern Ocean ecosystem services under climate change. Front. Mar. Sci. 7, 10.3389/fmars.2020.615214 (2021).
20. Atkinson A KRILLBASE: A circumpolar database of Antarctic krill and salp numerical densities, 1926-2016 Earth Syst. Sci. Data 2017 9 193 210 10.5194/essd-9-193-2017
Atkinson, A. et al. KRILLBASE: A circumpolar database of Antarctic krill and salp numerical densities, 1926-2016. Earth Syst. Sci. Data 9, 193–210 (2017).10.5194/essd-9-193-2017
21. Boyd PW Claustre H Levy M Siegel DA Weber T Multi-faceted particle pumps drive carbon sequestration in the ocean Nature 2019 568 327 335 10.1038/s41586-019-1098-2 30996317
Boyd, P. W., Claustre, H., Levy, M., Siegel, D. A. & Weber, T. Multi-faceted particle pumps drive carbon sequestration in the ocean. Nature 568, 327–335 (2019).30996317 10.1038/s41586-019-1098-2
22. DeVries T Weber T The export and fate of organic matter in the ocean: New constraints from combining satellite and oceanographic tracer observations Glob. Biogeochem. Cycles 2017 31 535 555 10.1002/2016GB005551
DeVries, T. & Weber, T. The export and fate of organic matter in the ocean: New constraints from combining satellite and oceanographic tracer observations. Glob. Biogeochem. Cycles 31, 535–555 (2017).10.1002/2016GB005551
23. DeVries T Holzer M Radiocarbon and helium isotope constraints on deep Ocean ventilation and mantle-3He sources J. Geophys. Res. Ocean. 2019 124 3036 3057 10.1029/2018JC014716
DeVries, T. & Holzer, M. Radiocarbon and helium isotope constraints on deep Ocean ventilation and mantle-3He sources. J. Geophys. Res. Ocean. 124, 3036–3057 (2019).10.1029/2018JC014716
24. Palevsky HI Doney SC How choice of depth horizon influences the estimated spatial patterns and global magnitude of Ocean carbon export flux Geophys. Res. Lett. 2018 45 4171 4179 10.1029/2017GL076498
Palevsky, H. I. & Doney, S. C. How choice of depth horizon influences the estimated spatial patterns and global magnitude of Ocean carbon export flux. Geophys. Res. Lett. 45, 4171–4179 (2018).10.1029/2017GL076498
25. Baker CA Martin AP Yool A Popova E Biological carbon pump sequestration efficiency in the North Atlantic: A leaky or a long‐term sink? Global Biogeochem. Cycles 2022 36 e2021GB007286 10.1029/2021GB007286
Baker, C. A., Martin, A. P., Yool, A. & Popova, E. Biological carbon pump sequestration efficiency in the North Atlantic: A leaky or a long‐term sink? Global Biogeochem. Cycles 36, e2021GB007286 (2022).10.1029/2021GB007286
26. Clarke A Quetin LB Ross RM Laboratory and field estimates of the rate of faecal pellet production by Antarctic krill, Euphausia superba Mar. Biol. 1988 98 557 563 10.1007/BF00391547
Clarke, A., Quetin, L. B. & Ross, R. M. Laboratory and field estimates of the rate of faecal pellet production by Antarctic krill, Euphausia superba. Mar. Biol. 98, 557–563 (1988).10.1007/BF00391547
27. Martin JH Knauer GA Karl DM Broenkow WW VERTEX: carbon cycling in the northeast Pacific Deep Sea Res. Part A Oceanogr. Res. Pap. 1987 34 267 285 10.1016/0198-0149(87)90086-0
Martin, J. H., Knauer, G. A., Karl, D. M. & Broenkow, W. W. VERTEX: carbon cycling in the northeast Pacific. Deep Sea Res. Part A Oceanogr. Res. Pap. 34, 267–285 (1987).10.1016/0198-0149(87)90086-0
28. Pauli N-C Krill and salp faecal pellets contribute equally to the carbon flux at the Antarctic Peninsula Nat. Commun. 2021 12 7168 10.1038/s41467-021-27436-9 34887407
Pauli, N.-C. et al. Krill and salp faecal pellets contribute equally to the carbon flux at the Antarctic Peninsula. Nat. Commun. 12, 7168 (2021).34887407 10.1038/s41467-021-27436-9
29. Atkinson A Schmidt K Fielding S Kawaguchi S Geissler PA Variable food absorption by Antarctic krill: Relationships between diet, egestion rate and the composition and sinking rates of their fecal pellets Deep. Res. Part II Top. Stud. Oceanogr. 2012 59–60 147 158 10.1016/j.dsr2.2011.06.008
Atkinson, A., Schmidt, K., Fielding, S., Kawaguchi, S. & Geissler, P. A. Variable food absorption by Antarctic krill: Relationships between diet, egestion rate and the composition and sinking rates of their fecal pellets. Deep. Res. Part II Top. Stud. Oceanogr. 59–60, 147–158 (2012).10.1016/j.dsr2.2011.06.008
30. 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).
31. Henson S Lampitt R Johns D Variability in phytoplankton community structure in response to the North Atlantic Oscillation and implications for organic carbon flux Limnol. Oceanogr. 2012 57 1591 1601 10.4319/lo.2012.57.6.1591
Henson, S., Lampitt, R. & Johns, D. Variability in phytoplankton community structure in response to the North Atlantic Oscillation and implications for organic carbon flux. Limnol. Oceanogr. 57, 1591–1601 (2012).10.4319/lo.2012.57.6.1591
32. Nowicki M DeVries T Siegel DA Quantifying the carbon export and sequestration pathways of the Ocean’s biological carbon pump Glob. Biogeochem. Cycles 2022 36 e2021GB007083 10.1029/2021GB007083
Nowicki, M., DeVries, T. & Siegel, D. A. Quantifying the carbon export and sequestration pathways of the Ocean’s biological carbon pump. Glob. Biogeochem. Cycles 36, e2021GB007083 (2022).10.1029/2021GB007083
33. Atkinson A Siegel V Pakhomov EA Jessopp MJ Loeb V A re-appraisal of the total biomass and annual production of Antarctic krill Deep. Res. Part I Oceanogr. Res. Pap. 2009 56 727 740 10.1016/j.dsr.2008.12.007
Atkinson, A., Siegel, V., Pakhomov, E. A., Jessopp, M. J. & Loeb, V. A re-appraisal of the total biomass and annual production of Antarctic krill. Deep. Res. Part I Oceanogr. Res. Pap. 56, 727–740 (2009).10.1016/j.dsr.2008.12.007
34. Giering S Particle flux in the oceans: Challenging the steady state assumption Glob. Biogeochem. Cycles 2017 31 159 171 10.1002/2016GB005424
Giering, S. et al. Particle flux in the oceans: Challenging the steady state assumption. Glob. Biogeochem. Cycles 31, 159–171 (2017).10.1002/2016GB005424
35. Tarling GA Johnson ML Satiation gives krill that sinking feeling Curr. Biol. 2006 16 R83 R84 10.1016/j.cub.2006.01.044 16461267
Tarling, G. A. & Johnson, M. L. Satiation gives krill that sinking feeling. Curr. Biol. 16, R83–R84 (2006).16461267 10.1016/j.cub.2006.01.044
36. Bax N Perspective: Increasing blue carbon around Antarctica is an ecosystem service of considerable societal and economic value worth protecting Glob. Chang. Biol. 2020 27 5 12 10.1111/gcb.15392 33064891
Bax, N. et al. Perspective: Increasing blue carbon around Antarctica is an ecosystem service of considerable societal and economic value worth protecting. Glob. Chang. Biol. 27, 5–12 (2020).33064891 10.1111/gcb.15392
37. Bar-On Y Rob P Ron M The biomass distribution on Earth Proc. Natl. Acad. Sci. USA 2018 115 6506 6511 10.1073/pnas.1711842115 29784790
Bar-On, Y., Rob, P. & Ron, M. The biomass distribution on Earth. Proc. Natl. Acad. Sci. USA 115, 6506–6511 (2018).29784790 10.1073/pnas.1711842115
38. Atkinson A Siegel V Pakhomov EA Rothery P Long-term decline in krill stock and increase in salps within the Southern Ocean Nature 2004 432 100 103 10.1038/nature02996 15525989
Atkinson, A., Siegel, V., Pakhomov, E. A. & Rothery, P. Long-term decline in krill stock and increase in salps within the Southern Ocean. Nature 432, 100–103 (2004).15525989 10.1038/nature02996
39. Iversen MH Sinkers or floaters? Contribution from salp pellets to the export flux during a large bloom event in the Southern Ocean Deep. Res. Part II Top. Stud. Oceanogr. 2017 138 116 125 10.1016/j.dsr2.2016.12.004
Iversen, M. H. et al. Sinkers or floaters? Contribution from salp pellets to the export flux during a large bloom event in the Southern Ocean. Deep. Res. Part II Top. Stud. Oceanogr. 138, 116–125 (2017).10.1016/j.dsr2.2016.12.004
40. Henson, S. et al. Uncertain response of ocean biological carbon export in a changing world. Nat. Geosci. 15, 248–254 (2022).
41. Yang, G., Atkinson, A., Pakhomov, E. A., Hill, S. L. & Racault, M.-F. Massive circumpolar biomass of Southern Ocean zooplankton: Implications for food web structure, carbon export, and marine spatial planning. Limnol. Oceanogr. 67, 2516–2530 (2022).
42. Le Quéré C Role of zooplankton dynamics for Southern Ocean phytoplankton biomass and global biogeochemical cycles Biogeosciences 2015 12 11935 11985
Le Quéré, C. et al. Role of zooplankton dynamics for Southern Ocean phytoplankton biomass and global biogeochemical cycles. Biogeosciences 12, 11935–11985 (2015).
43. Arrigo, K. R., van Dijken, G. L. & Bushinsky, S. Primary production in the Southern Ocean, 1997–2006. J. Geophys. Res. Ocean. 113, 10.1029/2007JC004578 (2008).
44. Hourcade, J. et al. Strengthening and Implementing the Global Response. (2018).
45. Schmidt K Seabed foraging by Antarctic krill: Implications for stock assessment, bentho‐pelagic coupling, and the vertical transfer of iron Limnol. Oceanogr. 2011 56 1411 1428 10.4319/lo.2011.56.4.1411
Schmidt, K. et al. Seabed foraging by Antarctic krill: Implications for stock assessment, bentho‐pelagic coupling, and the vertical transfer of iron. Limnol. Oceanogr. 56, 1411–1428 (2011).10.4319/lo.2011.56.4.1411
