# Hypoxia-tolerant zooplankton may reduce biological carbon pump efficiency in the Humboldt current system off Peru

Anja Engel, Carolina Cisternas-Novoa, Helena Hauss, Rainer Kiko, Frédéric Le Moigne

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

## To cite this version:

Anja Engel, Carolina Cisternas-Novoa, Helena Hauss, Rainer Kiko, Frédéric Le Moigne. Hypoxia-tolerant zooplankton may reduce biological carbon pump efficiency in the Humboldt current system off Peru. Communications Earth & Environment, 2023, 4 (1), pp.458. 10.1038/s43247-023-01140-6. hal-04341664

HAL Id: hal-04341664

https://hal.science/hal-04341664v1

Submitted on 14 Dec 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 CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

# Hypoxia-tolerant zooplankton may reduce biological carbon pump efficiency in the Humboldt current system off Peru

Anja Engel✉️<sup>1,2</sup>, Carolina Cisternas-Novoa<sup>1,2</sup>, Helena Hauss<sup>1,3</sup>, Rainer Kiko<sup>1,4</sup> & Frédéric A. C. Le Moigne<sup>1,5</sup>

In the ocean, downward flux of particles produced in sunlit surface waters is the major component of the biological carbon pump, which sequesters atmospheric carbon dioxide and fuels deep-sea ecosystems. The efficiency of downward carbon transfer is expected to be particularly high in tropical upwelling systems where hypoxia occurring beneath the productive surface waters is thought to hamper particle consumption. However, observations of both particle feeders and carbon export in low-oxygen waters are scarce. Here, we provide evidence that hypoxia-tolerant zooplankton feed on sinking particles in the extensive Oxygen Minimum Zone (OMZ) off Peru. Using several arrays of drifting sediment traps and in situ imaging, we show geochemical and morphological transformations of sinking particles and substantial control of carbon export by zooplankton. Our findings challenge the assumption of a consistently efficient biological carbon pump in OMZs and further demonstrate the need to consider mesopelagic organisms when studying oceanic carbon sequestration.

<footnote n="1">¹GEOMAR, Helmholtz Centre for Ocean Research Kiel, Kiel, Germany. ²Memorial University, St John's, Newfoundland and Labrador, Canada. ³NORCE Norwegian Research Centre, Division Climate & Environment, Bergen, Norway. ⁴Sorbonne Université, Laboratoire d'Océanographie de Villefranche, Villefranche-sur-Mer, France. ⁵Univ. Brest, CNRS, IRD, Ifremer, LEMAR, France. ✉email: aengel@geomar.de</footnote>
The biological carbon pump (BCP) sequesters CO₂ from the ocean's surface to the deep sea¹, mitigating increasing atmospheric CO₂ concentrations². The storage time of carbon entering the dark ocean via the BCP depends on the depth at which organic matter is remineralized back to CO₂. The shallower the remineralization depth, the faster the CO₂ will eventually escape into the atmosphere³. Climate change induced by variability in atmospheric CO₂ concentration also relates to changes in the efficiency of the BCP (BCP<sub>eff</sub>)⁴, being a function of both the carbon export efficiency (E<sub>eff</sub>) and the carbon transfer efficiency (T<sub>eff</sub>):

$$BCP_{\text{eff}} = E_{\text{eff}} \times T_{\text{eff}} \qquad (1)$$

$$E_{\text{eff}} = \frac{F_{z0}}{\text{PP}} \qquad (2)$$

$$T_{\text{eff}} = \frac{F_z}{F_{z0}} \qquad (3)$$

where F<sub>z0</sub> is the particulate organic carbon (POC) export flux out of the surface ocean, PP is the amount of CO₂ fixed by primary production and (F<sub>z</sub>) is the POC flux at a particular depth z. T<sub>eff</sub> is related to the attenuation of particle flux, often characterized as the exponent b of the *Martin* power-law function⁵:

$$F_z = F_{100} \times \left( \frac{Z}{100} \right)^{-b} \qquad (4)$$

with 100 being the standard reference depth (m) for F<sub>z0</sub>. A higher *Martin* b coefficient means stronger flux attenuation and less efficient downward transport of carbon. BCP<sub>eff</sub> is controlled by multiple factors influencing E<sub>eff</sub>, T<sub>eff</sub>, or both, such as plankton growth, species diversity, particle sinking velocity, zooplankton feeding, and microbial remineralization⁶,⁷.

As a consequence of climate change, the ocean is losing oxygen, leading to an expansion and intensification of OMZs with severe consequences for marine biota and biogeochemical cycles⁸. The largest present-day OMZs are associated with highly productive Eastern Boundary Upwelling Systems, which are responsible for ~25% of carbon export globally⁹. Deoxygenation reduces heterotrophic activity and results in the selection of organisms with little or no metabolic demand for oxygen¹⁰. To thrive under low oxygen concentration (hypoxia: <60 µmol O₂ L⁻¹) or even

complete absence of oxygen (anoxia), microorganisms such as bacteria and archaea use anaerobic respiration pathways¹¹. Most zooplankton and nekton species avoid low-oxygen environments, these include many taxa known to feed on or break up sinking particles¹². Therefore, a widespread assumption is that reduced zooplankton activity in OMZs results in higher T<sub>eff</sub> and, hence, in higher efficiency of the BCP¹³⁻¹⁵. One of the ocean's most extensive and intense permanent OMZ is associated with the Humboldt Current System off Peru and Chile¹⁶, occupying ~2 × 10⁶ km³. In this region, several zooplankton species migrate into the OMZ core daily, surviving prolonged oxygen depletion¹⁷⁻¹⁹. The impact of hypoxia-tolerant particle feeders on the BCP remains unknown. Yet, their interactions with particle flux might challenge the current view²⁰ that (1) productive regions with a severe OMZ have a high T<sub>eff</sub>, and (2) future expansion of OMZs may enhance CO₂ sequestration overall. Here, we determined T<sub>eff</sub> for the hypoxic waters along the continental margin off Peru. Using surface-tethered sediment traps at 5–8 depths, with the shallowest trap placed below the surface mixed layer (ML) and the euphotic zone (EZ), i.e., 40–60 m, and the deepest trap at ≤600 m, we obtained a high-resolution dataset of export fluxes for the upper mesopelagic zone, where most of the flux attenuation occurs. We combined sediment trap sampling with zooplankton surveys using in situ imaging and Multinet vertical hauls. Based on our results, we argue for a revised view of BCP<sub>eff</sub> in OMZs that explicitly includes zooplankton-induced particle transformation.

## Results and discussion

**Carbon export fluxes and flux attenuation in the OMZ.** We deployed particle interceptor traps six times (T1–T4, T6, T7) at different locations off Peru (12.0°S–14.6°S and 77.3°W–78.6°W; Supplementary Fig. 1) during the RV METEOR cruises M136 and M138; four times in April and two times in June 2017 (see Methods; Supplementary Table 1). Phytoplankton abundance in the upper mixed layer varied considerably between the cruises and deployment sites as indicated by chlorophyll-a (Chla) concentrations²¹ (Supplementary Fig. 1). Chla concentrations were higher in April (avg.: 3.1 ± 2.6, max: 11 µg L⁻¹) than in June (avg.: 2.8 ± 1.6, max: 4 µg L⁻¹). Oxygen concentrations declined rapidly below 20–50 m water depth and increased again below 500 m (Fig. 1a–f), indicating severe hypoxic conditions

<figure><img src="image_4.png" /><figcaption>**Fig. 1** Carbon export fluxes off Peru. a–f Particulate organic carbon (POC) export fluxes (grey circles) were measured using sediment traps; the solid black lines and crosses correspond to regression curves (power law function)⁵ fitted to the data with the depth of the shallowest trap as reference depth for fluxes out of the mixed layer (F<sub>z0</sub>). Also shown are temperature (blue) and oxygen concentration (red) profiles, the depth of the mixed layer (MLD; dashed line) and the extent of the OMZ (purple shading; <15 µmol O₂ kg⁻¹).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Potential drivers of export flux attenuation. **a** Relationship between carbon export flux and satellite-derived primary production<sup>23</sup> (PP, during the period of deployment); **b**, **c**, **e**, **f** Export flux attenuation coefficient (Martin-b) and biochemical parameters indicative of particle quality as obtained in the deepest trap of each deployment. Ratios are [g]:[g] in (**b**, **c**, **f**), and [mol]:[mol] in (**e**); and **d** export flux attenuation coefficients and temperature as suggested by Marsay et al.<sup>29</sup> (dotted line), blue symbols: Martin-b values from well-oxygenated waters (>200 µmol L<sup>-1</sup>), red symbols: OMZ (<70 µmol L<sup>-1</sup>) with circles: previously published data (Supplementary Table 3), and triangles: data from this study. All b-values in **d** were calculated with a reference depth of 100 m for comparison.</figcaption></figure>

(< 15.0 µmol O₂ L⁻¹)²² in the upper mesopelagic zone. Temperature declined at shallow depth (20–30 m), indicating a ML-depth (MLD) of <50 m, slightly shallower than the EZ depth, which ranged between 25 and 48 m (Supplementary Table 2). Below the MLD, the temperature decreased steadily and within similar ranges at all deployment sites (Fig. 1a–f). All vertical profiles of downward fluxes of POC closely followed the Martin-type power-law function (Fig. 1a–f, Supplementary Table 2). Absolute carbon fluxes and flux attenuation in the OMZ differed clearly between the deployment sites. The highest flux at 100 m was recorded on the shelf in April, T2 (F₁₀₀ = 412 mg C m⁻² d⁻¹), and the lowest for T6 in June (F₁₀₀ = 51 mg C m⁻² d⁻¹). Overall, carbon export flux (F₁₀₀) was closely related to primary production (PP) (r = 0.94, *p* < 0.01), derived from remote sensing data using the vertically generalized productivity model (standard VGPM)²³ (Fig. 2a). Hence, despite large spatial and temporal differences in PP, export efficiency, Eₐff = Δ[F₁₀₀]/Δ[PP], was similar for all deployments, indicating that approximately 10% of the PP was exported below the mixed layer.

Since both the ML and EZ were shallow (<50 m) in the study region, we calculated transfer efficiencies (T<sub>eff</sub>) and attenuation coefficients (b-values) using the shallowest trap as reference depth for fluxes out of the surface ocean (F<sub>z0</sub>)⁴ (Fig. 1a–f). Attenuation

of carbon flux was highest at T7, showing a *Martin-b* of 1.04 and thus the least efficient carbon transfer to depth (T<sub>eff</sub> = 11.5%). In contrast, T<sub>eff</sub> at T4 was almost four times higher (44.6%), with a *Martin-b* of 0.33. Regarding all deployments, no significant relationship was observed between F<sub>z0</sub> and *Martin-b* values, indicating that the amount of POC exported did not control the flux attenuation within the mesopelagic OMZ. *Martin-b* values at similar deployment sites, i.e., T4 and T6, as well as T1 and T7, were higher in June than in April.

Inferred from the biogeochemical composition of material captured in the deepest traps, we considered several possible controls on flux attenuation and so T<sub>eff</sub>, such as: (a) the freshness of organic matter, indicated by the relative contribution of chlorophyll *a* to total particulate organic carbon ([POC]:[Chla]), (b) the effect of ballasting by biogenic silica ([POC]:[BSi]), (c) the lability of the particulate matter ([POC]:[PN]), and (d) the extent of zooplankton grazing on particles, indicated by the amount of pheopigments²⁴ (sum of pheophorbide-*a* and pyropheophorbide-*a* relative to Chla; [Pheo]:[Chla]) (Fig. 2). Between trap deployments, [POC]:[Chla] and [POC]:[BSi] varied considerably but did not control *Martin-b* coefficients. Molar [POC]:[PN] ratios in the deepest traps ranged from 7.7 to 10.3 but could not explain differences in export flux profiles either. Pheophorbide-*a* and pyropheophorbide-*a* are produced during zooplankton
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Zooplankton distribution in the OMZ. Relative frequency distribution of copepoda (top row), solitary collodaria with a black capsule (SolB, middle row), and euphausiids (bottom row) in relation to depth (a, c, e) and oxygen distribution (b, d, f) during the day (grey bars) and night (black bars) based upon in situ images collected with an Underwater Vision Profiler (UVP5). The scale bar is always 5 mm.</figcaption></figure>

degradation of *Chla* and are quantitatively the most important pigments in fecal pellets²⁵. Pheophorbides are widely used as biomarkers for zooplankton feeding activity on sinking particles²⁶,²⁷. The [Pheo]:[Chla] ratios of the particulate material collected during this study varied between 7 and 25 [g]:[g], indicating a substantial contribution of herbivore fecal material²⁶. Overall, carbon flux attenuation (*Martin*-b) was significantly correlated to the [Pheo]:[Chla] of the deepest traps (*p* < 0.01; *r* = 0.81; *n* = 6, Fig. 2f), which indicates that flux attenuation could be linked to zooplankton grazing. The [Pheo]:[Chla] ratio of material collected in traps above the OMZ averaged 6.9 ± 7.9 and was lower than within the OMZ (17.3 ± 9.3) (Supplementary Fig. 2) indicating a source of feces in the OMZ rather than settling of fecal pellets out of the oxygenated surface waters. Highest [Pheo]:[Chla] ratio was often observed close to the upper or lower boundary of the OMZ. In the deepest traps, we measured the highest [Pheo]:[Chla] ratio for T7 and the lowest for T4, suggesting that particles sinking out of the OMZ had undergone the highest and lowest zooplankton reworking at these sites, respectively. (Fig. 2f, Supplementary Fig. 2). Overall, [Pheo]:[-Chla] ratios were higher in June (21 ± 11) than in April (15 ± 7), pointing to higher grazing activity later during the season²⁸.

Previous observations from oxygenated waters have suggested that flux attenuation correlates with the median temperature of the upper mesopelagic waters²⁹. We therefore also compared our and previously published *b*-values using the 100 m reference depth (*F*₁₀₀) and differentiated between low-oxygen (<70 µmol O₂ L⁻¹) and high-oxygen (>70 µmol O₂ L⁻¹) systems (Fig. 2d, Supplementary Table 3). Flux-attenuation in high-oxygen systems matches the prediction over a wide range of temperatures fairly well. Low-oxygen systems, however, show

large variations of *b*-values despite similar temperatures. Although *Martin*-b values at some low-oxygen sites agreed with the suggested temperature dependence, it is evident that temperature is a poor predictor of POC flux attenuation values in OMZs. One explanation for the lower-than-expected *b*-values is that hypoxia reduces particle decomposition rates as demonstrated in incubation experiments³⁰,³¹, which purposefully excluded zooplankton. In the absence of larger heterotrophs, microbial decomposition of organic matter is the predominant process responsible for flux attenuation²⁰. Indeed, *b*-values at sites with limited indications of zooplankton grazing in our study (e.g., T4) and other suboxic regions (Supplementary Table 3) are below the temperature prediction line (Fig. 2d), while flux attenuation at T7 was higher than predicted. Variations in microbial decomposition were unlikely to explain the observed large differences in carbon flux attenuation during this study since the activity of microbial hydrolytic enzymes and the bacterial biomass production in the OMZ were similar between stations²¹. Instead, a coherent interpretation of our biogeochemical data is that feeding activity by zooplankton tolerating low oxygen leads to higher flux attenuation.

**Zooplankton feeding on and disrupting sinking particles.**

Pheopigment biomarkers for zooplankton grazing during this study were consistent with a microscopic analysis of sinking material caught in transparent gel-filled traps. Large phytoplankton aggregates were more frequent at T4 than at T7, while T7 showed higher abundance of large fecal pellets (Supplementary Fig. 3). Fecal pellet morphology at T7 indicated the presence of squat lobsters and euphausiids (Supplementary Fig. 4); their
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Carbon export flux attenuation in OMZs. Conceptual model of POC export flux attenuation in the water column of OMZs. Left: Low POC flux attenuation is observed when organic matter degradation is solely driven by anaerobic microbial communities that are adapted to hypoxia. Right: In addition to anaerobic microorganisms, zooplankton that have adapted to hypoxic conditions feed on, and break up sinking particles, thereby enhancing particle flux attenuation through the OMZ.</figcaption></figure>

presence in the OMZ was also confirmed by vertically stratified
plankton net samples collected at the recovery sites of the traps
(Supplementary Fig. 5). Both the squat lobster *Pleuroncodes*
*monodon* and the euphausiid *Euphausia mucronata* are endemic
to the Humboldt Current and perform diel vertical migrations
into the OMZ core (down to 250-350 m depth approximately)
during the daytime, thus tolerating hypoxia and even anoxia for
prolonged periods³²,³³. To further resolve the fine-scale dis-
tribution patterns of zooplankton, we deployed an Underwater
Vision Profiler 5 (UVP)³⁴ on the CTD-Rosette (equipped with an
O₂-sensor) during both cruises (M136, and M138) as well as
during an earlier cruise (M135) off Peru. 77 UVP profiles down to
at least 500 m depth were analyzed. Copepod abundance showed
no differences in vertical distribution between day and night, but
this might be related to the taxonomic resolution attainable with
the UVP, which does not allow resolution of species-specific diel
vertical migration patterns (e.g., as found for *Eucalanus inermis*
in the upper 200 m off Chile¹⁹). We found 88% of all copepods in
the oxygenated (>15 µmol kg⁻¹ O₂) surface layer (Fig. 3a, b), 0.1%
were found below the OMZ at O₂ concentrations >15 µmol kg⁻¹
O₂, 7% were found at the upper and lower oxycline at oxygen
concentrations between 4 and 15 µmol kg⁻¹ and 5% were found
in the OMZ at O₂ concentrations below 4 µmol kg⁻¹. Overall,
copepods seem to avoid the OMZ, but some species such as
*Eucalanus inermis* may use it as a diapause refuge, others as
feeding ground³⁵. In situ imaging also indicated the presence of

several giant rhizaria taxa in the OMZ off Peru (Fig. 3c, d; Sup-
plementary Fig. 6). Importantly, large parts of the rhizaria
populations were found at oxygen concentrations of
<4 µmol kg⁻¹ during both day- and night-time (Fig. 3d; Supple-
mentary Fig. 6b, d, f). In total, 63% of the solitary black collo-
darians (SolB) were found in the core of the OMZ at
<4 µmol O₂ kg⁻¹, 7% between 4 and 15 µmol O₂ kg⁻¹ and 30% at
>15 µmol O₂ kg⁻¹ (Fig. 3c). The solitary globular collodarians
(SolG) and *Aulosphaeridae* were more evenly distributed (Sup-
plementary Fig. 6c) and the *Aulacantha* showed a slight pre-
ference for the upper and lower oxycline (Supplementary Fig. 6e).
However, also in these three groups 64, 60 and 41% of all imaged
organisms, respectively, were found at <4 µmol O₂ kg⁻¹.
Euphausiids were observed in the anoxic core (200–400 m) of the
OMZ during the daytime and in the upper water column during
the night (Fig. 3a, b; Supplementary Figs. 4 and 5). Although we
cannot directly correlate the presence of rhizaria or euphausiids
to our flux attenuation values, as the data do not stem from
parallel deployments, the overall high presence of flux-feeding
rhizaria³⁶ within the OMZ and the presence of krill at the upper
oxycline strongly suggest zooplankton interacting with and
grazing on sinking particles at extremely low oxygen levels off
Peru. Crustaceans, especially large krill and *Pleuroncodes* may
contribute to flux attenuation via direct grazing on sinking phy-
toplankton cells and marine snow, and also via disruption of large
aggregates³⁷. Off California, *Euphausia pacifica* abundances of
$5-10 \text{ ind. m}^{-3}$ resulted in a net increase of aggregate numbers, likely via the fragmentation of larger aggregates into smaller ones³⁸. At T7 we found krill abundances of $6 \text{ ind. m}^{-3}$ at the lower oxycline during the night (Supplementary Fig. 5). In the Humboldt current system, dense krill swarms can host up to $50 \text{ ind. m}^{-3}$ off Chile¹⁸, and we found up to $90 \text{ ind. m}^{-3}$ off Peru. Dense zooplankton aggregations may also disrupt sinking particles during their diel migrations to and from their daytime residence depth of about at 200-400 m.

**Potential limitations of zooplankton effects on BCP under hypoxia.** Our results highlight the importance of hypoxia-adapted zooplankton in explaining variations of the BCP in OMZs underlying highly productive waters such as the Humboldt current system. Here, pelagic crustaceans and other zooplankton species might exploit and transform particles rapidly, in contrast to microbes that grow slowly under anoxic conditions²⁰,²¹ and cannot efficiently attenuate the flux of rapidly sinking particles (Fig. 4). To date, little is known about the foraging activities of zooplankton in the OMZ¹⁸,³⁹. An important conclusion of our study is that behavioral and physiological adaptations of zooplankton to environmental extremes are critical, but largely neglected, traits for understanding the BCP. The Humboldt coastal upwelling system is generated by ocean circulation and has exhibited permanent hypoxia since >14.7 ka with little fluctuation⁴⁰. Sustained by abundant nutrients, the production and supply of organic matter to the deeper hypoxic waters have always been high, creating an environment rich in food and low in competition, thus particularly rewarding for heterotrophs that can tolerate low-oxygen conditions. In this sense, 'old' OMZ ecosystems like the Humboldt Current may be unique in that they have allowed the proliferation of zooplankton specifically adapted to hypoxia. In systems that have recently suffered from oxygen deficiency due to eutrophication and climate warming, zooplankton tend to avoid hypoxic layers⁴¹. On the other hand, opportunistic zooplankton species may be able to exploit even transient hypoxia, as recently demonstrated by the mass occurrence of the gelatinous, flux-feeding polychaete *Poeobius sp.*, in 'dead zone' eddies of the Tropical North Atlantic⁴². Predicting the impacts of zooplankton on CO₂ sequestration in the ocean is challenging due to the high temporal and spatial heterogeneity of organismal abundance and distribution. Even where zooplankton have adapted to hypoxia, the distribution of these metazoans is not uniform. Zooplankton aggregate in patches is often associated with submesoscale (~10³–10⁴ m) and mesoscale (~10⁴–10⁵ m) hydrographical features⁴³. Patchiness of zooplankton abundance can also be an explanation for the observed range of [Pheo]:[–Chla] in the trap material and the corresponding variability in POC flux attenuation during this study, where meso- and sub-mesoscale physical variability due to filaments and eddies is typically high⁴⁴. Moreover, species-specific traits related to metabolic adaptations and acclimations, life cycle dynamics, and behavior are further components of the complex ecology of mesopelagic organisms. Overall, the patchy distribution and the various roles of mesopelagic organisms on flux attenuation and generation³⁹ limit, if not prohibit, simplified empirical approaches to estimate carbon export fluxes from physical and chemical environmental data alone⁴⁵,⁴⁶. Based on several independent methods, our results provide consistent evidence that zooplankton interact with and substantially reduce particle export fluxes in low-oxygen systems. Our study supports the emerging perception of zooplankton as gatekeepers³⁹ of the biological carbon pump, even within the OMZ. Given the potentially significant influence of planktonic fauna on the BCP, we argue that observations and knowledge of mesopelagic ecosystem dynamics

must be greatly expanded to make reliable predictions about the marine carbon cycle.

**Methods**

**Study site and sample collection.** Samples were collected in April and June 2017 during the RV *Meteor* cruises M136 and M138 off Peru. Details on the hydrography and O₂ concentrations in the study region during sampling are given elsewhere²¹. We collected sinking particles using surface-tethered sediment traps⁴⁷,⁴⁸ at six locations (Supplementary Fig. 1, Supplementary Table 1). Sediment traps were deployed for two to six days and collected particles at eight depths (Supplementary Table 1) between 50 and 600 m in the open ocean (T1, T3, T4, T6, and T7), and 5 depths between 40 and 300 m on the shelf (T2). One trap array (T5) was lost during the deployment. The first trap depth was always chosen below the base of the mixed layer based on CTD profiles. At each depth, twelve acrylic particle interceptor tubes (PITs) mounted in a PVC cross-frame were deployed. Each PIT was equipped with an acrylic baffle at the top to minimize the collection of swimmers⁴⁷,⁴⁸. The PITs were 7 cm in diameter and 53 cm in height with an aspect ratio of 7.5 and a collection area of 0.0038 m². The cross-frame and PITs were attached to a line with a bottom weight and a set of surface and subsurface floats. The procedures for PIT preparation and sample recovery followed Engel et al.⁴⁸. Samples were split into aliquots that were processed for the different biogeochemical analyses. At selected traps, we collected intact individual particles in situ using gel cups, i.e., cylindrical collectors placed at the bottom of the PIT tubes (collection area of 0.0020 m²) and filled with polyacrylamide gels⁴⁹. For T4 and T7, those gel cups were placed at 200, 300, 400, 500, and 600 m depth. Following recovery, gel cups were placed under a Zeiss® Stemi 508 stereomicroscope coupled to a Zeiss® Axiocam ERc 5 s digital camera. Gels were then photographed at a magnification of 6.5 against a laser-etched glass grid (each cell of 14 by 12.5 mm)⁴⁹, using the Zen lite software provided by Zeiss®. A specifically designed light system providing a homogeneous source of light from below the gels was used to avoid any shading on the photographs. We assessed the flux abundance of various types of particles in both traps 4 and 7 following methods presented in Ebersbach and Trull⁵⁰. We grouped observed particles into distinct categories: phytoplankton aggregates (PAGG), fecal aggregates (FAGG), and various fecal pellets (FP), including cylindrical, oval, and spherical fecal pellets. The flux abundance spectra of each type of particle at both traps 4 and 7 were then binned into five size classes in equivalent spherical diameter (ESD): 0.071–0.145, 0.145–0.296, 0.296–0.603, 0.603–1.228, >1.228 mm ESD, shown in Supplementary Fig. 3.

An underwater vision profiler (UVP5) was mounted in the CTD rosette frame throughout the cruises to obtain high-resolution particle abundance and size spectra data as well as zooplankton images³⁴. At the sites of traps 6 and 7, depth-stratified zooplankton samples were collected by day/night vertical haul pairs ($1 \text{ m s}^{-1}$) with a Hydrobios Multinet Midi (0.25 m² mouth opening, 5 nets, 200 µm mesh). Samples were rinsed into Kautex® jars and fixed in 4% borax-buffered formaldehyde in seawater solution. Samples were scanned on an Epson Perfection V750 Pro scanner in a modification of the ZooScan method⁵¹. Both the UVP5 and the Zooscan image dataset were uploaded to Ecotaxa³⁴ (http://ecotaxa OBS-vlfr.fr/) and predicted into different classes using the in-built random forest algorithm with deep learning features. Automated image sorting was then manually validated by experts. Biomass was estimated according to Lehette and Hernández-León⁵². Rhizaria categories were distinguished according to Biard et al.⁵³. The solitary black collodians (SolB, Fig. 2) are large spherical
organisms with a dense central part (black/dark). The solitary globular collodians (SolGlob, Supplementary Fig. 6) are large, spherical to-oval organisms with a homogenous grey (or dark- grey) surface. The central sphere is surrounded by a blurry halo generated by a network of pseudopodial extensions.

**Biogeochemical analysis.** Particulate organic carbon (POC), nitrogen (PN), and pigments were determined by filtering aliquots, of 100 to 200 mL of the trapped material onto pre-combusted GF/F filters (8 h at 500 °C) at low vacuum (< 200 mbar). POC and PN samples were filtered in triplicate, and pigments in duplicate. Filters for POC-PN were stored dry, acid-fumed (37% hydrochloric acid) to remove carbonates, and dried for 12 h at 60 °C before analysis. POC and PN were determined with an elemental analyzer (Euro EA, Hechatech). Biogenic silica (BSi) was determined in aliquots of 50 to 100 mL, filtered in duplicate onto 0.4 µm cellulose acetate filters. Samples were stored at -20 °C until analysis. For the measurements, filters were digested in NaOH at 85 °C for 135 min; the pH was adjusted to 8 with HCl. Silicate was measured spectrophotometrically⁵⁴. Pigments were extracted with 90% acetone and homogenized with glass beads. The tube containing the filter was centrifuged for 10 min at 5000 rpm, and the supernatant was collected. Extracts were filtered through 0.2 µm Teflon filters and stored at -80 °C until analysis. Chloropigments (chlorophyll a, pheophytin a, pheophorbide a, and pyropheophorbide a) were determined using reverse-phase High-Performance Liquid Chromatography (HPLC, DionexUlti Mate R 3000LC system equipped with an autosampler, a photodiode array and a fluorescence detector, ThermoScientific). Pigments were identified by using a comparison of their retention times and spectral properties with standards (DHI Water & Environment, Denmark). Data fits and statistical tests were performed with the software packages Microsoft Office Excel 2019 and Sigma Plot 14.0 (Systat).

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

**Data availability**

All data included in this study are publicly available. Engel, Anja; Cisternas-Novoa, Carolina; Hauss, Helena; Kiko, Rainer; Le Moigne, Frédéric A C: Biological carbon pump efficiency in the Humboldt current system off Peru: sediment traps during METEOR cruises M136 and M138. PANGAEA, https://doi.pangaea.de/10.1594/PANGAEA. 963289.

Received: 24 October 2023; Accepted: 22 November 2023;
Published online: 06 December 2023

**References**

1. Volk, T. & Hoffert, M. Ocean carbon pumps: analysis of relative strengths and efficiencies in ocean-driven atmospheric CO₂ changes. (eds E. Sundquist and W. S. Broecker) *The Carbon Cycle And Atmospheric CO₂: Natural Variations Archean To Present* (AGU, Washington, DC, 1985), 99–110.

2. Parekh, P., Dutkiewicz, S., Follows, M. J. & Ito, T. Atmospheric carbon dioxide in a less dusty world. *Geophys. Res. Lett.* **33**, L03610 (2006).

3. Kwon, E., Primeau, F. & Sarmiento, J. The impact of remineralisation depth on the air-sea carbon balance. *Nat. Geosci.* **2**, 630–635 (2009).

4. Buesseler, K. O., Boyd, P. W., Black, E. E. & Siegel, D. A. Metrics that matter for assessing the ocean biological carbon pump. *Proc. Natl Acad. Sci. USA* **117**, 9679–9687 (2020).

5. Martin, J. H., Knauer, G. A., Karl, D. M. & Broenkow, W. W. Vertex—carbon cycling in the Pacific. *Deep Sea Res. A* **34**, 267–285 (1987).

6. De la Rocha, C. & Passow, U. Factors influencing the sinking of POC and the efficiency of the biological carbon pump. *Deep Sea Res. II* **54**, 639–658 (2007).

7. Henson, S. A., Le Moigne, F. A. C. & Giering, S. L. C. Drivers of carbon export efficiency in the global ocean. *Glob. Biogeochem. Cycles* **33**, 891–903 (2019).

8. Stramma, L., Johnson, G. C., Sprintall, J. & Mohrholz, V. Expanding oxygen-minimum zones in the tropical ocean. *Science* **320**, 655–658 (2008).

9. Jahnke, R. A. A global synthesis. Carbon and nutrient fluxes in continental margins: a global synthesis. In: (eds Liu, K. K., Atkinson, L., Quinones, R., Talaue-McManus, L.) *Global Change: The IGBP Series*. (Springer-Verlag, 2010).

10. Lam, P. & Kuypers, M. M. M. Microbial nitrogen cycling processes in oxygen minimum zones. *Annu. Rev. Mar. Sci.* **3**, 317–348 (2011).

11. Wright, J. J., Konwar, K. M. & Hallam, S. J. Microbial ecology of expanding oxygen minimum zones. *Nat. Rev. Microbiol.* **10**, 381–394 (2012).

12. Engel, A., Kiko, R. & Dengler, M. Organic matter supply and utilization in oxygen minimum zones. *Ann. Rev. Mar. Sci.* **14**, 355–378 (2022).

13. Bretagnon, M., Paulmier, A., Garçon, V. & Dewitte, B. Modulation of the vertical particle transfer efficiency in the oxygen minimum zone off Peru. *Biogeosciences* **15**, 5093–5111 (2018).

14. Devol, A. H. & Harnett, H. E. Role of the oxygen minimum zone in transfer of organic carbon to the deep ocean. *Limnol. Oceanogr.* **25**, 1684–1690 (2001).

15. Tutasi, P. & Escribano, R. Zooplankton diel vertical migration and downward C flux into the oxygen minimum zone in the highly productive upwelling region off northern Chile. *Biogeosciences* **17**, 455–473 (2020).

16. Fuenzalida, R., Schneider, W., Garcés-Vargas, J., Bravo, L. & Lange, C. Vertical and horizontal extension of the oxygen minimum zone in the eastern South Pacific Ocean. *Deep Sea Res. II* **56**, 992–1003 (2009).

17. Wishner, K. F. et al. Ocean deoxygenation and zooplankton: very small oxygen differences matter. *Sci. Adv.* **4**, eaua5180 (2018).

18. Antezana, T. Adaptive behavior of *Euphausia mucronata* in relation to the oxygen minimum layer of the Humboldt Current. *Oceanogr. East. Pac.* **2**, 29–40 (2002).

19. Escribano, R., Hidalgo, P. & Krautz, C. Zooplankton associated with the oxygen minimum zone system in the northern upwelling region of Chile during March 2000. *Deep Sea Res. II* **56**, 1083–1094 (2009).

20. Cavan, E. L., Trimmer, M., Shelley, F. & Sanders, R. Remineralization of particulate organic carbon in an ocean oxygen minimum zone. *Nat. Commun.* **8**, 14847 (2017).

21. Maßmig, M., Lüdke, J., Krahmann, G. & Engel, A. Bacterial degradation activity in the eastern tropical South Pacific oxygen minimum zone. *Biogeosciences* **17**, 215–230 (2020).

22. Kavelage, T. et al. Aerobic microbial respiration in oceanic oxygen minimum zones. *PLoS ONE* **10**, e0133526 (2016).

23. Behrenfeld, M. J. & Falkowski, P. G. Photosynthetic rates derived from satellite-based chlorophyll concentration. *Proc. Oceanogr. C.2*, 1–20 (1997).

24. Welschmeyer, N. A. & Lorenzen, C. J. Chlorophyll budgets: Zooplankton grazing and phytoplankton growth in a temperate fjord and the Central Pacific 'Gyres'. *Limnol. Oceanogr.* **30**, 1–21 (1985).

25. Goericke, R., Strom, S. L. & Bell, M. A. Distribution and sources of cyclic pheophorbides in the marine environment. *Limnol. Oceanogr.* **45**, 200–211 (2000).

26. Stukel, M. R. et al. Sinking carbon, nitrogen, and pigment flux within and beneath the euphotic zone in the oligotrophic, open-ocean Gulf of Mexico. *J. Plankton Res.* **44**, 711–727 (2022).

27. Wakeham, S. G. et al. Organic biomarkers in the twilight zone—time series and settling velocity sediment traps during MedFlux. *Deep Sea Res. II* **56**, 1437–1453 (2009).

28. Taylor, G. T. et al. Ecosystem responses in the southern Caribbean Sea to global climate change. *Proc. Natl Acad. Sci. USA* **109**, 19315–19320 (2012).

29. Marsay, C. M. et al. Attenuation of sinking particulate organic carbon flux through the mesopelagic ocean. *Proc. Natl Acad. Sci. USA* **112**, 1089–1094 (2015).

30. Le Moigne, F. A. C. et al. On the effect of low oxygen concentrations on bacterial degradation of sinking particles. *Sci. Rep.* **7**, 16722 (2017). Art. No.

31. Van Mooy, B. A. S., Keil, R. G. & Devol, A. H. Impact of suboxia on sinking particulate organic carbon: enhanced carbon flux and preferential degradation of amino acids via denitrification. *Geochim. Cosmochim. Acta* **66**, 457–465 (2002).

32. Kiko, R. et al. The squat lobster *Pleuroncodes monodon* tolerates anoxic “dead zone” conditions off Peru. *Mar. Biol.* **162**, 1913–1921 (2015).

33. Kiko, R. & Hauss, H. On the estimation of zooplankton-mediated active fluxes in Oxygen Minimum Zone regions. *Front. Mar. Sci.* **6**, 741 (2019).

34. Picheral, M. et al. The Underwater Vision Profiler 5: an advanced instrument for high spatial resolution studies of particle size spectra and zooplankton. *Limnol. Oceanogr. Methods* **8**, 462–473 (2010).

35. Wishner, K. F., Seibel, B. & Outram, D. Ocean deoxygenation and copepods: coping with oxygen minimum zone variability. *Biogeosciences* **17**, 2315–2339 (2020).

36. Biard, T. & Ohman, M. D. Vertical niche definition of test-bearing protists (*Rhizaria*) into the twilight zone revealed by in situ imaging. *Limnol. Oceanogr.* **65**, 2583–2602 (2020).
37. 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).

38. Dilling, L. & Alldredge, A. L. Fragmentation of marine snow by swimming macrozooplankton: a new process impacting carbon cycling in the sea. *Deep Sea Res. I* **47**, 1227–1245 (2000).

39. Stukel, M. R., Ohman, M. D., Kelly, T. B. & Biard, T. The roles of suspension-feeding and flux-feeding zooplankton as gatekeepers of particle flux into the mesopelagic ocean in the northeast pacific. *Front. Mar. Sci.* **6**, 397 (2019).

40. Moffitt, S. E. et al. Paleooceanographic Insights on recent oxygen minimum zone expansion: lessons for modern oceanography. *PLoS ONE* **10**, e0115246 (2015).

41. Roman, M. R. et al. Impacts of hypoxia on zooplankton spatial distributions in the northern Gulf of Mexico. *Estuar. Coasts* **35**, 1261–1269 (2012).

42. Christiansen, S. et al. Particulate matter flux interception in oceanic mesoscale eddies by the polychaete *Poeobius sp. Limnol. Oceanogr.* **63**, 2093–2109 (2018).

43. Bertrand, A. et al. Broad impacts of fine-scale dynamics on seascAPE structure from zooplankton to seabirds. *Nat. Commun.* **5**, 5239 (2014).

44. Thomsen, S. et al. Do submesoscale frontal processes ventilate the oxygen minimum zone off Peru? *Geophys. Res. Lett.* **43**, S367–S386 (2016).

45. Laufkötter, C. et al. Temperature and oxygen dependence of the remineralization of organic matter. *Glob. Biogeochem. Cycl.* **31**, 1038–1050 (2017).

46. Henson, S. A. et al. A reduced estimate of the strength of the ‘ocean’s biological carbon pump. *Geophys. Res. Lett.* **38**, L04606 (2011).

47. Knauer, G. A., Martin, J. H. & Bruland, K. W. Fluxes of particulate carbon, nitrogen, and phosphorus in the upper water column of the northeast Pacific. *Deep Sea Res.* **26**, 97–108 (1979).

48. Engel, A., Wagner, H., Le Moigne, F. A. C. & Wilson, S. T. Particle export fluxes to the oxygen minimum zone of the Eastern Tropical North Atlantic. *Biogeosciences* **14**, 1825–1838 (2017).

49. Lundsgaard, C. Use of a high viscosity medium in studies of aggregates. In: Sediment Trap Studies in the Nordic Countries. 3. *Proceeding of the Symposium on Seasonal Dynamics of Planktonic Ecosystems and Sedimentation in Coastal Nordic Waters*. Finnish Environment Agency, 141–152 (1995).

50. Ebersbach, F. & Trull, T. W. Sinking particle properties from polyacrylamide gels during the Kerguelen Ocean and Plateau compared study (KEOPS): zooplankton control of carbon export in an area of persistent natural iron inputs in the Southern Ocean. *Limnol. Oceanogr.* **53**, 212–224 (2008).

51. Gorsky, G. et al. Digital zooplankton image analysis using the ZooScan integrated system. *J. Plankton Res.* **32**, 285–303 (2010).

52. Lehette, P. & Hernández-León, S. Zooplankton biomass estimation from digitized images: a comparison between subtropical and Antarctic organisms. *Limnol. Oceanogr. Methods* **7**, 304–308 (2009).

53. Biard, T. et al. In situ imaging reveals the biomass of giant protists in the global ocean. *Nature* **532**, 504–507 (2016).

54. Hansen, H. P. & Koroleff F. Determination of nutrients. In *Methods Seawater Anal* 3rd edn, 159–228 (1999).

### Acknowledgements

Jon Roa, Tania Klüver, Marie Massmig, Ruth Flerus, Kerstin Nachtigall, Claudia Sforna, and Christian Begler supported trap preparation, deployments, sampling, or sample analysis. This research was supported by the Helmholtz Association, the DFG Collaborative Research Center 754' Climate-Biogeochemistry Interactions in the Tropical

'Ocean' (B9; grant/award no. 27542298), the collaborative project CUSCO—Coastal Upwelling System in a Changing Ocean (grant no. 03F0813D) of the German Federal Ministry of Education and Research (BMBF) and by a Fellowship of the Excellence Cluster 'The Future 'Ocean' (CP1403 to F.A.C.L.M.) funded by the DFG. RK acknowledges support via a "Make Our Planet Great Again" grant of the French National Research Agency within the "Program d'Investissements d'Avenir"; reference "ANR-19-MPGA-0012". We thank the crew, officers, and captains of the F.S. Meteor as well as the scientific parties of cruises M135, M136, and M138. H.H. acknowledges financial support by the Bundesministerium für Bildung und Forschung (BMBF) through the project CO2Meso (03F0876A).

### Author contributions

A.E., C.C.N., and F.A.C.L.M. designed the study. C.C.N. and F.A.C.L.M. performed trap deployments at sea, and biogeochemical, and microscopic analysis. H.H. and R.K. were in charge of the UVP5 operation and zooplankton sampling. All authors analyzed the data. A.E. wrote the initial draft of the paper and C.C.N, F.A.C.L.M., H.H., and R.K. edited subsequent drafts.

### Funding

Open Access funding enabled and organized by Projekt DEAL.

### 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/s43247-023-01140-6.

**Correspondence** and requests for materials should be addressed to Anja Engel.

**Peer review information** *Communications Earth & Environment* thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editors: Clare Davis. A peer review file is available.

**Reprints and permission 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.

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

**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 Author(s) 2023