# Small sinking particles control anammox rates in the Peruvian oxygen minimum zone

Clarissa Karthäuser, Soeren Ahmerkamp, Hannah K. Marchant, Laura A.

Bristow, Helena Hauss, Morten H. Iversen, Rainer Kiko, Joeran Maerz, Gaute

Lavik, Marcel M. M. Kuypers

► To cite this version:

Clarissa Karthäuser, Soeren Ahmerkamp, Hannah K. Marchant, Laura A. Bristow, Helena Hauss, et al.. Small sinking particles control anammox rates in the Peruvian oxygen minimum zone. Nature Communications, 2021, 12, 10.1038/s41467-021-23340-4. insu-03661720

HAL Id: insu-03661720

https://insu.hal.science/insu-03661720

Submitted on 7 May 2022

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

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

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

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

ARTICLE

https://doi.org/10.1038/s41467-021-23340-4

OPEN

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

Check for updates

# Small sinking particles control anammox rates in the Peruvian oxygen minimum zone

Clarissa Karthäuser ID 1,8, Soeren Ahmerkamp 1,2,8✉, Hannah K. Marchant ID 1,2✉, Laura A. Bristow ID 1,6, Helena Hauss ID 3, Morten H. Iversen ID 2,4, Rainer Kiko 3,7, Joeran Maerz ID 5, Gaute Lavik 1 & Marcel M. M. Kuypers 1

Anaerobic oxidation of ammonium (anammox) in oxygen minimum zones (OMZs) is a major pathway of oceanic nitrogen loss. Ammonium released from sinking particles has been suggested to fuel this process. During cruises to the Peruvian OMZ in April-June 2017 we found that anammox rates are strongly correlated with the volume of small particles (128-512 µm), even though anammox bacteria were not directly associated with particles. This suggests that the relationship between anammox rates and particles is related to the ammonium released from particles by remineralization. To investigate this, ammonium release from particles was modelled and theoretical encounters of free-living anammox bacteria with ammonium in the particle boundary layer were calculated. These results indicated that small sinking particles could be responsible for ~75% of ammonium release in anoxic waters and that free-living anammox bacteria frequently encounter ammonium in the vicinity of smaller particles. This indicates a so far underestimated role of abundant, slow-sinking small particles in controlling oceanic nutrient budgets, and furthermore implies that observations of the volume of small particles could be used to estimate N-loss across large areas.

<footnote n="1">¹Max Planck Institute for Marine Microbiology, Bremen, Germany. ²MARUM—Center for Marine Environmental Sciences, University of Bremen, Bremen, Germany. ³GEOMAR Helmholtz Center for Ocean Research Kiel, Kiel, Germany. ⁴Alfred Wegener Institute Helmholtz Center for Polar and Marine Research, Bremerhaven, Germany. ⁵Max Planck Institute for Meteorology, Hamburg, Germany. ⁶Department of Biology, University of Southern Denmark, Odense, Denmark. ⁷Laboratoire d'Océanographie de Villefranche-sur-Mer, Villefranche-sur-Mer, France. ⁸These authors contributed equally: Clarissa Karthäuser, Soeren Ahmerkamp. ✉email: sahmerka@mpi-bremen.de; hmarchan@mpi-bremen.de</footnote>
The upwelling of cold, nutrient-rich water drives high pri-
mary production rates in the surface waters of (sub)tro-
pical eastern ocean boundaries. The resulting export of
organic matter in the form of aggregates and particles fuels high
microbial respiration rates which, combined with sluggish ven-
tilation, leads to the formation of oxygen minimum zones
(OMZs). OMZs (when defined by $O_2 \le 20 \mu mol L^{-1}$) make up
<1% of the ocean volume, however 20–40% of global oceanic
fixed nitrogen (N)-loss occurs within them²⁻⁴ due to the micro-
bially mediated processes of anammox (anaerobic oxidation of
ammonium¹,⁵⁻⁷) and denitrification⁸⁻¹⁰. Hence, despite their
small volume, OMZs have a profound impact on biogeochemical
cycles.

In the past decades, OMZs have expanded substantially¹¹,¹²,
and climate models predict that global warming will accelerate
deoxygenation¹³. However, the current and future extent of
OMZs, and N-loss patterns within them, remain difficult to
reproduce and predict in current Earth System Models¹⁴⁻¹⁶. To
solve this issue, the drivers and impacts of physical and biological
processes including ventilation, net primary production, respira-
tion and carbon export need to be better represented in
models¹⁴,¹⁷⁻¹⁹. For example, the representation of OMZs, as well
as estimates of N-loss in low oxygen waters has been improved by
the inclusion of aggregate and particle dynamics into biogeo-
chemical models²⁰,²¹. However, the role that sinking particles
play in regulating N-loss within anoxic waters remains poorly
understood.

It has been shown previously that rates of anammox correlate strongly with satellite-derived estimates of organic nitrogen export out of the photic zone by sinking aggregates and particles⁵. Generally, this correlation has been attributed to the release of ammonium associated with the remineralization of sinking organic matter⁵,⁹. Ammonium is a key substrate for anammox, but is present at very low concentrations in OMZ waters. Although anammox bacteria colonize resuspended organic particles within some OMZs²², it is currently unclear how the microorganisms that carry out anammox interact with sinking organic material and gain access to particle-derived ammonium. Whilst incubation experiments using size-fractionated water samples have indicated that anammox rates are higher in the presence of particles, metagenomic and metatranscriptomic studies have indicated that anammox bacteria reside in the water column, rather than on the particles themselves²³⁻²⁶.

The influence of differently sized particles is likely a key factor required to understand the mechanistic link between organic matter which sinks into anoxic OMZ waters, the apparently free-living anammox community and associated N-loss rates. In general, if composed of the same material, smaller particles sink slower than larger particles, thus smaller particles have longer residence times and are older by the time they reach deeper waters²⁷⁻²⁹. Furthermore, organic carbon appears to be respired more quickly on smaller and suspended particles³⁰,³¹. Together these factors mean that smaller particles undergo more mineralization than larger particles as they sink through the water column²⁷. Whilst this means that smaller particles have low transfer efficiencies in terms of carbon export from the subsurface to deep waters²⁷,³²,³³, the remineralization of smaller particles has the potential to play an important role in releasing ammonium within the upper anoxic water column, where anammox rates are often highest⁸,³⁴.

In this study, we aimed to constrain the relationship between
anammox rates, the export of organic nitrogen out of the photic
zone, and the abundance and total volume of sinking particles
within the anoxic part of the Peruvian Upwelling System OMZ.
The Peruvian Upwelling System forms part of the boundary

current system of the Eastern Tropical South Pacific and is one of
the most biologically productive regions in the world³⁵. During
austral autumn (April and June 2017), we carried out anammox
rate measurements in anoxic bulk and size-fractionated waters
and correlated them to imaged particle abundance and volume at
high resolution over depth profiles. Thereafter, to gain mechan-
istic insights into the relationship between particles and ana-
mmox rates, we used the in situ particle profiles combined with
derived settling velocities to model export of organic nitrogen by
sinking particles, N-release associated with different particle size
classes within the OMZ, and encounter rates of anammox bac-
teria with particles. We also used the measured anammox rates
and particle abundances to derive a simple relationship in order
to extrapolate anammox rates across the anoxic waters of the
Peruvian OMZ at high spatial resolution.

Results and discussion

**Physico-chemical conditions in the Peruvian Upwelling System.** The influence of sinking particles and aggregates on anammox rates was investigated on two cruises to the Peruvian Upwelling System in the Eastern Tropical South Pacific in austral autumn (April and June 2017). These cruises took place during the declining phase of a Coastal Warming Event, which was characterized by an unusual development of warm water off the coast of Peru while neutral sea surface temperatures persisted in the equatorial Pacific³⁶⁻³⁹. As a result, increased surface water temperatures of up to 24 °C were observed at some of the sampled stations (Supplementary Fig. 1). Despite this unusual event, most other conditions in the OMZ were typical for austral autumn⁵,⁷,⁸,³⁴,⁴⁰,⁴¹. The oxic–anoxic interface (throughout defined as where $O_2$ dropped below 1.5 µM, which was the detection limit of the sensor) was generally around 50–100 m water depth (Fig. 1 and Supplementary Fig. 1), although, as is typical for this region, at some shallow onshore stations, the OMZ was observed to expand up to 10 m below the surface (Supplementary Fig. 1, Supplementary Tables 1 and 2). As observed previously, surface chlorophyll *a* concentrations were variable, both spatially and temporally (Supplementary Video 1, Supplementary Fig. 2c), nevertheless surface chlorophyll *a* was generally higher at the onshore stations (defined throughout as water depths <600 m), than at offshore stations (water depths >600 m). These depth cut-offs are based on ref. ⁵, where 600 m was the average depth of the lower OMZ boundary, roughly equivalent to the edge of the Peruvian shelf, and are henceforth used throughout to differentiate between onshore and offshore OMZ stations. Chlorophyll *a* concentrations were less variable over depth, and dropped below 0.2 mg L⁻¹ by a depth of 50 m (Supplementary Fig. 1).

**Abundance and characterization of sinking organic material.**

To characterize the export of the particulate organic matter formed in the OMZ surface waters, the size and abundance of particles and aggregates with equivalent spherical diameters (ESD) between 128 and 2048 µm were determined from 5 m below the surface down to the seafloor at onshore stations and down to at least 500 m at offshore stations using an Underwater Vision Profiler (UVP; an underwater camera system that can be mounted to a rosette sampler⁴²). Particles were subsequently binned into four size classes based on ESD (128–256, 256–512, 512–1024, 1024–2048 µm). The size classes were chosen to ensure that there was a similar total volume of particles within each size class, as particle abundance and volume generally decrease with increasing size⁴³,⁴⁴ (Supplementary Fig. 3). The average particle and aggregate volume was higher for the onshore stations than for the offshore stations at all depths (Fig. 2a and b). The total volume of the imaged particles generally peaked at the base of the
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Physico-chemical parameters, particle abundance and anammox rates from an exemplary onshore and offshore station in the Peruvian upwelling system in April 2017. In situ particle abundances and anammox rates from a onshore station 423 (total water depth 241 m) and c offshore station 549 (total water depth 4350 m). b and d Depth profiles of oxygen (purple line), ammonium (open circles) and nitrite (open triangles) concentrations. Particles were sorted into four size classes with equivalent spherical diameters (ESD) of 128–256, 256–512, 512–1024, and 1024–2048 µm. Particle abundances in the smallest size class reached up to 398 particles L⁻¹ in the upper 20 m of the onshore station. Anammox rates were determined in time series incubations at six discrete depths. Significant rates are shown with closed circles, rates that are not significantly different from zero are shown with open circles (see Supplementary Table 1). Error bars represent the standard error of the slope. The base of the euphotic zone (where photosynthetically active radiation (PAR) dropped to <1%, based on Aqua-MODIS satellite data) and the oxic-anoxic interface (where O₂ dropped to <1.5 µM) are indicated.</figcaption></figure>

euphotic zone at the onshore stations (up to 3.6 ± 2.8 mm³ L⁻¹) and towards the centre of the euphotic zone at the offshore stations (up to 3.2 ± 2.4 mm³ L⁻¹), after which they declined to minima of 0.5 mm³ L⁻¹ in both the onshore and offshore OMZ (Fig. 2, Supplementary Figs. 1 and 2a and b). At some shallow onshore stations, particle and aggregate volume increased again towards the seafloor, likely due to sediment resuspension.

We assume that all of the imaged aggregates and particles were sinking, similar to previous studies based on in situ particle profiling with a UVP system<sup>45,46</sup>. While suspended particles may comprise a substantial amount of total POC in the ocean<sup>27,30,47</sup>, such particles are likely smaller than those imaged by the UVP (which were >128 µm). For example, marine microgels are generally <100 µm (ref. 48) and a recent study which examined in situ settling velocities with a camera system, observed that more than 95% of particles with ESDs larger than 80 µm were sinking with velocities >10 m day⁻¹ in the mesopelagic zone<sup>49</sup>. Our assumption that the UVP imaged particles were sinking is also supported by the presence of particles in the 128–2048 µm size range within the tray of a Marine Snow Catcher, which is designed to only collect sinking particles.

Categorization of 910 particles and aggregates collected using a Marine Snow Catcher revealed that they were generally detrital material that appeared to be composed of similar matter independent of size class or packing density. Although the matter consisted of aggregated material rather than single particles, we hereafter refer to particles for simplicity. Particles in the two smaller size classes (128–256 and 256–512 µm) consisted mainly of loosely and densely packed detrital material (fluffy detritus, 44% and compact detritus, 34%), as well as some faecal pellets (6%) and various other types of particles such as dead zooplankton, foraminifera and other material that could not clearly be identified (15% other material). The two larger size classes (512–1024, 1024–2048 µm) consisted of loosely packed detrital material (fluffy detritus, 51%), loosely packed detrital material that contained faecal pellets (25%), densely packed detrital material (compact detritus, 11%) and other material (11%). Only a small proportion of the larger size classes consisted of individual faecal pellets (3%, Supplementary Table 3). Overall, these categorizations were consistent with those from in situ observations made using the UVP and are representative of the sinking, but not the suspended material. The UVP profiles indicated that particles within the smallest size class (128–256 µm) were the most abundant at all depths (Supplementary Table 1) and overall abundance decreased with increasing particle size. The observed distribution of particles over depth and between size classes is consistent with studies in other ocean regions, where a higher export of smaller particles leads to their comparatively greater abundance below the euphotic zone<sup>45,47,50,51</sup>. Additionally, the higher abundance of smaller particles in the OMZ core might result from the break-up of larger aggregates through fragmentation (Supplementary Discussion)<sup>20,45,47,50,51</sup> or active transport by migrating zooplankton<sup>52,53</sup>, although this latter option is unlikely as few faecal pellets were observed at this depth.

**Export of particles as a source of N for anammox in the Peruvian OMZ.** To determine the links between sinking organic matter and N-loss rates in the Peruvian OMZ, we determined anammox and denitrification rates in bulk water collected from the anoxic region of the OMZ and compared them to particle abundances and volume. On average, anammox rates were 33.5 ± 26.1 nmol N₂ L⁻¹ day⁻¹ at onshore stations and 8.3 ± 15.4 at offshore stations (Table 1, Supplementary Table 1, Supplementary Fig. 1). Rates were generally highest 20–60 m below the oxic–anoxic interface after which they decreased with increasing depth. Both the offshore rates and comparatively high rates at the productive onshore stations (which reached up to 89.6 nmol N₂ L⁻¹ day⁻¹), as well as the rate distribution over depth are consistent with those measured before in the Eastern Tropical South Pacific<sup>5,7,8,34,40,41,54–57</sup>. In comparison, denitrification rates were on average 13.5 and 4.5 times lower than anammox rates at onshore
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Particle volume and export production in the Peruvian upwelling system during April 2017. Average particle volume profiles in (a) offshore and (b) onshore stations determined by the Underwater Vision Profiler camera system (UVP). Particles were sorted into four size classes with equivalent spherical diameters (ESD) of 128–256, 256–512, 512–1024, and 1024–2048 µm. The shaded envelopes correspond to standard deviation in abundances from 56 offshore and 37 onshore stations (cropped at 2 mm³ L⁻¹ for visualization purposes). The base of the euphotic zone (where photosynthetically active radiation (PAR) dropped to <1%, based on Aqua-MODIS satellite data) and the oxic–anoxic interface (where O₂ dropped to <1.5 µM) are indicated. c Export production estimated from satellite products (pseudocolor map, see the “Methods” section) and from particle abundances at the base of the euphotic zone (circles). The white line is the 600 m isobath on which basis the onshore and offshore stations were separated. Stations where anammox incubations were performed are circled in red, the elongated red circle encompasses four stations.</figcaption></figure>

<table><caption>**Table 1** Averaged and integrated anammox and denitrification rates, and export of organic nitrogen out of the euphotic zone into the oxygen minimum zone (OMZ) from all sampled onshore and offshore stations (for more information refer to Supplementary Tables 1a and 2).</caption><thead><tr><th></th><th></th><th>Unit</th><th>Onshore (mean ± SD)</th><th>Offshore (mean ± SD)</th><th>Combined (mean ± SD)</th></tr></thead><tbody><tr><td rowspan="2">OMZ average<sup>a</sup></td><td>Anammox</td><td>nmol N<sub>2</sub> L<sup>-1</sup> day<sup>-1</sup></td><td>33.5 ± 26.1</td><td>8.3 ± 15.4</td><td>20.1 ± 24.4</td></tr><tr><td>Denitrification</td><td></td><td>2.5 ± 4.0</td><td>1.8 ± 2.4</td><td>2.1 ± 3.2</td></tr><tr><td rowspan="2">Integrated rate over depth<sup>a</sup></td><td>Anammox</td><td>mmol N<sub>2</sub> m<sup>-2</sup> day<sup>-1</sup></td><td>3.1 ± 3.3</td><td>1.2 ± 1.7</td><td>2.2 ± 2.8</td></tr><tr><td>Denitrification</td><td></td><td>0.2 ± 0.5</td><td>0.5 ± 0.5</td><td>0.4 ± 0.5</td></tr><tr><td rowspan="4">Organic nitrogen export</td><td>128–256 µm</td><td>mmol N m<sup>-2</sup> day<sup>-1</sup></td><td>8.0 ± 7.9</td><td>6.0 ± 2.1</td><td>7.1 ± 5.8</td></tr><tr><td>256–512 µm</td><td></td><td>6.7 ± 5.1</td><td>4.5 ± 0.9</td><td>5.7 ± 3.8</td></tr><tr><td>512–1024 µm</td><td></td><td>5.3 ± 4.5</td><td>3.2 ± 1.5</td><td>4.4 ± 3.5</td></tr><tr><td>1024–2048 µm</td><td></td><td>2.1 ± 1.9</td><td>1.2 ± 1.0</td><td>1.7 ± 1.5</td></tr></tbody></table>

and offshore stations, respectively (Table 1, Supplementary Table 1 and Supplementary Fig. 4). Based on these results, denitrification was a minor N-loss process in the Peruvian OMZ during austral autumn, similar to previous observations⁵,⁵⁴.

The export of organic nitrogen out of the photic zone was calculated using two methods. One method was based on satellite measurements of chlorophyll *a* and temperature and an assumed particle export to primary production ratio (similar to ref. ⁵). The other method was based on combining the in situ measurements of particle abundance and size (using the UVP camera system) with measured particle carbon contents, previously determined C: N ratios and calculated settling velocities (UVP method). Both methods resulted in similar estimates of organic nitrogen export, although those from the satellite were slightly lower (Supplementary Discussion and Supplementary Fig. 5). As the UVP-based estimates were carried out simultaneously to the water

sampling, they capture the local variability of organic nitrogen export more precisely, and therefore we focus on the UVP dataset for comparison to the anammox rates. Based on the UVP method, organic nitrogen export from the photic zone ranged between 6.5 and 61.7 mmol N m⁻² day⁻¹ onshore, and 6.9 and 45.7 mmol N m⁻² day⁻¹ offshore (Fig. 2c, Supplementary Fig. 2d), and was always at least 1.8-fold higher than depth-integrated anammox rates (Supplementary Table 2 and Supplementary Fig. 6). As such, the export of organic nitrogen was sufficient to support the measured N-loss, which is in accordance with previous results from the Peruvian OMZ⁵.

**Anammox bacteria are free-living.** To further investigate the relationship between anammox rates and export of organic N via sinking particles, we determined where anammox activity was located within the anoxic waters of the OMZ. Previous studies

<footnote><sup>a</sup>In situ oxygen <1.5 µM.</footnote>
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Anammox rates determined from unfiltered water (bulk) and water filtered through either a 10 or 1.6 µm filter in the anoxic waters of the Peruvian upwelling system during April and June 2017. Boxplots depict the 25-75% quantile range, with the centre line depicting the median (50% quantile); whiskers encompass data points within 1.5 times the interquartile range. A paired t-test (Supplementary Table 4) showed that there was no significant difference in rates between any of the treatments. Anammox rates from one depth are outside the range of the graph and had anammox rates of 14.1 nmol N₂ L⁻¹ day⁻¹ (bulk water), 41.1 nmol N₂ L⁻¹ day⁻¹ (10 µm filtered) and 8.7 nmol N₂ L⁻¹ day⁻¹ (1.6 µm filtered). These values were included in the statistical analysis. At an additional 23 depths anammox rates were below detection limit in all size fractions and these were not included in the analysis. See Supplementary Fig. 7 for individual data points and Supplementary Table 1b for the standard error of the slope of each individual rate measurement.</figcaption></figure>

based on size-fractionated water samples have indicated that anammox bacteria are more likely to be part of the free-living pelagic community than the particle-associated community in various OMZs, including the Eastern Tropical North Pacific, Eastern Tropical South Pacific, Black Sea and Cariaco Basin<sup>23-25,58</sup>. However, FISH-based studies have indicated that anammox bacteria can be associated with particles in the Namibian OMZ where sediment resuspension is high<sup>22</sup>.

The bulk water incubations that we carried out already indicated that anammox activity was unlikely to be occurring directly on particles (see Supplementary Discussion). To confirm that anammox activity was not directly associated with the particles, we carried out size-fractionation experiments at 11 stations to compare anammox rates in water that was either unfiltered, or in water where particles had been removed by filtration through a 1.6 µm filter or a 10 µm filter. Anammox rates in the <10 µm fraction, the <1.6 µm fraction and bulk water samples (i.e. non-fractionated) were not significantly different from each other (Fig. 3, Supplementary Fig. 7 and Supplementary Table 4). This indicates that anammox bacteria were not removed along with the particles by the filtration steps, and are therefore most likely free-living within the anoxic OMZ waters rather than attached to particles.

**Correlation of anammox rates to total particle volume.** To constrain the relationship between the apparently free-living anammox bacteria and particles within the OMZ, anammox rates from the bulk water incubations were compared to the total volume of particles within different size classes determined at the same time and depth with the UVP (Fig. 4). We found significant positive correlations between anammox rates from anoxic waters and the corresponding volume of particles within each of the particle size classes (Spearman's rank correlation; *p* < 0.05, see Fig. 4, Supplementary Fig. 8 and Supplementary Table 5). These correlations were strongest for the two smaller size classes. The observed stronger correlation between anammox rates and smaller particles was further supported by a multiple linear regression analysis, which revealed that the volume of the smallest particles was the most important variable needed to predict anammox rates

<figure><img src="image_7.png" /><figcaption>**Fig. 4** Relationship between volumetric anammox rates and particle volumes from different size classes in the Peruvian upwelling system during April 2017. Particles throughout the water column were quantified with an Underwater Vision Profiler camera system (UVP) and binned into four size classes with equivalent spherical diameters (ESD) of 128-256, 256-512, 512-1024, and 1024-2048 µm. Anammox rates were determined from the slope of ²⁹N₂ production over time in anoxic incubations after addition of ¹⁵NO₂⁻ (taking into consideration any contribution to ²⁹N₂ production from denitrification if a denitrification rate with *p* < 0.05 could be detected). For all particle size classes, there was a significant positive correlation (Spearman's rank correlation; *p* = <0.05) between particle volume and anammox rates from the anoxic part of the oxygen minimum zone (O₂ < 1.5 µM). The line is the linear regression from which the slope and *R*² was calculated (see Supplementary Table 5 for all relevant statistics). Please note the different scales for the particle volume **a-d**. One outlier was removed from the figures and correlation as it was more than 6 standard deviations from the mean in the smallest size fraction and it was the only sample from the euphotic zone (24.3 m depth). In **a-d**, open and filled symbols depict offshore stations and onshore stations, respectively. In panel **a**, the correlation with the smallest size class (140-270 µm) during a cruise in February 2013 (M93) is included for comparison and depicted with grey x symbols. This dataset was not included in the linear regression but is included in Supplementary Table 1c. Correlations with the particle abundance in number per liter yielded similar results (see Supplementary Table 5).</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Fig. 5** Transport of nitrogen (N) from the euphotic zone through the water column. **a** Modelled export of organic nitrogen from the euphotic zone into the oxygen minimum zone, **c** modelled nitrogen remineralization associated with particles in the anoxic water column (**e**) theoretical number of encounters between anammox bacteria and the particle boundary layer (see text for further information). **b, d** and **f** indicate the ratio of processes between the smaller (128–512 µm) and larger particles (512–2048 µm) from (**a, c** and **e**,) respectively. **Ind:** individuals. The outlines of the violin plots depict the kernel density estimation of the data points shown. The solid black line and solid red line indicate mean and median, respectively.</figcaption></figure>

(see Supplementary Table 6 for more details). The slope of the correlation between anammox rates and particles in the smallest size class (128–256 µm) was 390 nmol N₂ mm⁻³ day⁻¹.

Intriguingly, we also found a significant linear correlation between anammox rates and the volume of the smallest particle size class, during a previous cruise (M93) to the same region in February 2013 (Fig. 4a, Supplementary Fig. 9), although the slope was higher (~590). It should be noted that in contrast to the dataset presented here, the M93 cruise occurred in austral summer when productivity is relatively high and the water column was sulfidic at some stations. We excluded the data from those stations, as sulfidic events have been shown to alter N-loss patterns in the Peruvian OMZ³⁴,⁴⁰,⁵⁹,⁶⁰. Nonetheless, the similarity in the correlation indicates that the positive relationship between measured anammox rates and measured volume of smaller particles is a consistent feature within the Peruvian OMZ. To explore the underlying mechanisms of this relationship, we compared the organic N-export flux associated with each particle

size class, modelled the potential ammonium release by the particles and estimated the likelihood of anammox bacteria to encounter the released ammonium.

**Modelling of potential nitrogen release from different particle size classes.** The export of organic nitrogen into the anoxic waters of the OMZ was modelled for each particle size class based on measured particle abundance at the base of the euphotic zone, measured size-specific carbon content, calculated sinking velocities and a C:N value of 7.2 as was determined previously for the region (Eq. (1)). The model assumes that all of the imaged particles were sinking (see above). Overall, when all the stations where the UVP was deployed were taken together, the smaller particles (<512 µm) were estimated to be responsible for the export of 1.8 times more organic nitrogen into the anoxic zone than the larger size classes (Fig. 5a, b). This relationship between the size classes persisted when a different function was used to calculate sinking velocities (Supplementary Fig. 10) and when using size-dependent C:N ratios (Supplementary Fig. 11 and Supplementary Discussion) to represent the preferential remineralization of organic N that often occurs during organic matter degradation (see Supplementary Discussion and Supplementary Figs. 9–13)²⁸,²⁹,⁶¹–⁶³.

The amount of organic nitrogen bound in particles that enter the OMZ is only the first step in substrate supply to anammox bacteria. The particulate organic nitrogen must subsequently be remineralized to inorganic nitrogen (i.e. ammonium). The continual decline in total particle numbers and volume with depth implied that remineralization of organic carbon and nitrogen within particles occurred throughout the OMZ. The highest carbon flux attenuation was observed within the upper oxic waters of the OMZ, indicating that this was where remineralization rates were highest (Supplementary Fig. 15). In this zone, carbon-specific degradation rates were estimated to be between 0.33 and 1.6 day⁻¹ (see Supplementary Discussion), which is in line with previous rates from the Eastern Tropical North Pacific determined from measurements of oxygen respiration (0.13–5 day⁻¹, ref. ³⁰). Based on a median particle number of 66 particles per liter this is equivalent to remineralization of between 3.2 and 6.0 µmol C L⁻¹ day⁻¹ (Supplementary Table 8), or, assuming a C:N ratio of 7.2, N-remineralization of between 0.44 and 0.83 µmol N L⁻¹ day⁻¹.

By comparing the time scale of diffusion with the settling velocity, we also estimated the distance that particles could sink before remineralized nitrogen would mostly diffuse into the water column. This so-called diffusion length scale ranged between 0.2 mm for the smallest particles and 20 cm for the largest particles. These results indicate that reduced nitrogen formed by remineralization in the oxic, euphotic water column mainly diffuses out of the particles within that same zone. It is therefore unlikely that particles transport ammonium generated by aerobic remineralization in the upper water column into the anoxic zone. Therefore, while the diffusive release of reduced nitrogen is likely a key factor that drives the enhanced N-cycling activity that is frequently observed at the oxic-anoxic interface within OMZs⁵,⁷,⁴¹,⁵⁵,⁶⁴, it would not support anammox activity below the oxic-anoxic interface.

Anammox activity within the rest of the OMZ therefore seems to be supported by anaerobic remineralization and subsequent nitrogen release from particles within the anoxic part of the OMZ. Carbon-specific degradation rates (derived from the UVP measurements and measured particle carbon content, see Supplementary Fig. 16) from 50 m below the oxic-anoxic interface were estimated to be on average 0.22 day⁻¹ onshore and 0.15 day⁻¹ offshore. This would result in a remineralization
of 0.3 µmol C L⁻¹ day⁻¹ at onshore stations and 0.16 µmol C L⁻¹ day⁻¹ at offshore stations. Assuming a C:N stoichiometry of 7.2 for the sinking organic matter (see ref. ⁵), this would be equivalent to the release of 0.04 and 0.02 µmol N L⁻¹ day⁻¹. Of this, 71% and 76%, (on average 74%) could be attributed to the smaller particle size classes at the onshore and offshore stations, respectively (Eq. (7), Fig. 5c, d). The relatively higher importance of smaller particles persisted even when varying the parametrization of the settling velocity and C:N ratios in the model and when considering the impact of particle fragmentation (see Supplementary Discussion).

Intriguingly, the modelled estimates of nitrogen release from the particles were higher than the organic matter remineralization (i.e. ammonium release) that is estimated to occur in the Peruvian OMZ via heterotrophic nitrate-based and nitrite-based respiratory processes (i.e nitrate reduction to nitrite, denitrification and dissimilatory nitrate/nitrite reduction to ammonium, DNRA). Nitrate reduction to nitrite, which is suggested to be the main remineralization process, has typical onshore rates of 0.2 µmol N L⁻¹ day⁻¹ (refs. ⁵,⁵⁵). Remineralization associated with nitrate reduction to nitrite would therefore lead to the release of around 0.015 µmol NH₄⁺ L⁻¹ day⁻¹, around 40% of our estimated ammonium release (0.04 µmol N L⁻¹ day⁻¹), which is likely underestimated. Preferential N-mineralization within particles and particle fragmentation (which cannot be fully resolved via modelling with this dataset) could both theoretically lead to higher ammonium release (see Supplementary Discussion). Furthermore, our modelling approach only takes sinking particles into account, whereas it is possible that some fraction of the smaller particles are suspended in the water column²⁷. If this were the case, we would underestimate N-release, as suspended particles would have a longer residence time in the OMZ (see Supplementary Discussion). Our results therefore imply that non-nitrate-based and nitrite-based respiratory processes lead to organic nitrogen remineralization within the particles. Such processes might be microaerobic respiration by prokaryotes⁶⁴⁻⁶⁶, and microzooplankton transported within particles, or fermentation, which has recently been shown to be a significant source of ammonium within the Eastern Tropical South Pacific OMZ⁶⁷.

**Encounters between anammox bacteria and remineralized nitrogen released by particles.** Remineralized organic nitrogen released from sinking particles does not immediately diffuse evenly throughout the water column. Instead, ammonium and other solutes diffuse from particles into their boundary layer, leaving a plume in the wake of sinking particles (see Supplementary Fig. 12, refs. ⁶⁸⁻⁷⁰). The concentrations of nutrients in the diffusive boundary layer and plume tails can be two to three orders of magnitude higher than in the surrounding seawater⁷⁰⁻⁷³. To estimate encounter rates between free-living anammox bacteria and the diffusive boundary layer and plumes of particles, we combined particle abundance data with sinking velocities, boundary layer theory (Eq. (8), Fig. 5e) and previously determined anammox bacterial cell numbers⁴¹. The resulting encounter model showed an average of 125 encounters per hour per mL (regardless of size class), which equates to 0.3% of the anammox bacteria population encountering a particle-associated nutrient hotspot every hour. When combined with typical cell-specific anammox rates of ~1–10 fmol cell⁻¹ h⁻¹ (refs. ⁷⁴,⁷⁵), these encounter rates would support in situ rates between 3 and 30 nmol L⁻¹ day⁻¹, which is in line with the measured rates (Table 1). It is likely that this is a conservative estimate as the encounter model assumes that anammox bacteria are not motile⁷⁶,⁷⁷, which limits encounters with sinking particles and resulting nutrient hotspots. So far, *Scalindua* (which are the most abundant anammox bacteria within the Peruvian OMZ⁷⁸) are

not known to be motile, however it should be noted that for motile bacteria at the same cell numbers, encounter rates could be increased by 179–400 times (Supplementary Fig. 14; Supplementary Discussion)⁷⁹.

The encounter model revealed that on average, non-motile anammox bacteria were between 1.5 and 2.6 times more likely to encounter the ammonium diffusing from particles within the two smaller size classes than that from the larger size classes (Fig. 5f, Supplementary Fig. 10). This was mainly due to the higher abundance and increased boundary layer thickness to size ratio of the smaller particles (Supplementary Fig. 12; Supplementary Discussion). For motile bacteria, this estimate would increase to 3.8 times. This further highlights the important role that smaller particles seem to play in regulating anammox rates and regional nitrogen loss.

**Extrapolation of anammox to regional scales based on particle abundances.** The simple linear relationship between the measured volume of small (128–256 µm) particles and measured anammox rates within the anoxic part of the OMZ (Fig. 4a) suggests that the distribution of small particles can be used to estimate nitrogen loss rates over large areas. The slope of the linear regression between particle volume in the smallest size class (128–256 µm) and anammox rates (390· V₁₂₆⁻₂₅₆⁻¹₅; Fig. 4a), implies an N₂ production per particle volume of 390 nmol N₂ mm⁻³ day⁻¹. Furthermore, the negative intercept implies that a minimum volume of small particles (0.04 mm³ L⁻¹) is required to fuel a significant anammox rate. This is supported by the observation that anammox is undetectable in the OMZ waters containing very few particles. Based on our mechanistic modelling, the cessation in anammox rates at low particle volumes likely results from ammonium release and encounter rates becoming too low (see Supplementary Fig. 17). The relationship from the linear regression was applied to a UVP dataset which covered a large area of the Peruvian OMZ and was collected during March–July 2017 (cruises M135–M138). UVP measurements were separated for the “inner” and “outer” OMZ, based on distance from the coast (175 km, as in refs. ⁵,⁴¹,⁸⁰), which led to an estimated average nitrogen loss within the anoxic waters of around 6.1 and 1.5 mmol N m² day⁻¹ for the inner and outer OMZ, respectively (Fig. 6). Subsequently, we extrapolated these rates over the area of the inner Peruvian OMZ (3.2 × 10⁵ km², as in refs. ⁴¹,⁸⁰) and the outer Peruvian OMZ (12 × 10⁵ km² as in ref. ⁵). This resulted in an estimated monthly N-loss of 1.0 Tg N, and 0.6 Tg N for the inner and outer OMZ, respectively, during the period covered by this dataset.

If we extrapolate the monthly rates across an entire year, they would lead to an N-loss estimate of 19 Tg N yr⁻¹ for the anoxic part of the Peruvian OMZ. This annual N-loss estimate is very similar to previous estimates which range from 10 to 25 Tg N yr⁻¹ (refs. ⁵,⁴¹,⁸⁰⁻⁸²). This annual extrapolation would be improved by additional data from the early winter to early summer seasons, sulfidic events and also by the inclusion of anammox activity from the oxycline, which we did not investigate in this study. However, it demonstrates the potential of using particle size distribution and abundances to extrapolate measured anammox rates. Quantification of nitrogen loss in OMZs is often based on a few measurements extrapolated to a large region, potentially leading to uncertainties related to spatial heterogeneities. In contrast, profiling of aggregate and particle abundance is becoming a standard oceanographic tool, with camera systems such as the UVP often installed on rosette-CTD samplers⁸³,⁸⁴. Therefore, the application of this approach in combination with rate measurements in other OMZs and seasons would allow data to be collected at a high resolution over large areas, which can then be integrated into Earth System Models.
<figure><img src="image_9.png" /><figcaption>**Fig. 6** Regional estimate of nitrogen (N)-loss via anammox based on the volume of small particles (128–256 µm). Anammox rates in the anoxic water column of Peruvian oxygen minimum zone (O₂ < 1.5 µM) derived from UVP measurements carried out from March to July 2017 (black circles depict UVP casts during cruises M135–M138). Extrapolating these rates over a whole year would lead to a nitrogen loss of around 19 Tg. The red circles depict stations where the anammox rates were measured. The white line is the 600 m isobath.</figcaption></figure>

Our results show that anammox rates in the OMZ of the Peruvian Upwelling System are correlated to the abundance and volume of sinking particles. Furthermore, we show that highly abundant smaller, slower sinking particles may play a greater role in controlling anammox rates than rarer, faster sinking larger particles. Modelling approaches showed that the more abundant smaller particles transport between 1.2 and 1.8-fold more organic nitrogen into the OMZ than larger particles. Smaller particles also release around 75% more remineralized nitrogen within the OMZ than larger particles, due to their higher abundance and slower sinking rates. Furthermore, anammox bacteria are twice as likely to encounter remineralized nitrogen released from smaller particles (Fig. 7). Generally, smaller sinking particles are considered to have less of an impact on biogeochemical cycles than larger particles, as they have a lower carbon transfer efficiency from the upper water column to the deep sea. However, our results highlight an important role of smaller, slower sinking particles in controlling nitrogen loss, and show that their abundance and distribution could be used to estimate nitrogen loss across large areas. Taking the role of smaller particles into account in Earth System Models may have the potential to more accurately predict global nitrogen loss patterns in the present and future ocean.

## Methods

**Water sampling and physico-chemical measurements.**

A conductivity-temperature-depth (CTD)—rosette equipped with 24, 10 L Niskin bottles was used to collect water samples. Duplicate sensors were employed to record temperature and salinity (CTD type SBE 9-plus), oxygen (SBE 43 type oxygen sensor, Seabird) and fluorescence (WET Labs Fluorometer, USA). The precision for the salinity sensor was 0.002 g kg⁻¹ and for the temperature sensor 0.002 °C. The detection limit of the oxygen sensor was 1.5 µmol L⁻¹ (ref. ⁸⁵) and it was calibrated regularly during the cruises against Winkler titration⁸⁶. Chlorophyll *a* concentrations were derived from the fluorescence using the factory calibration⁸⁷. High resolution in situ particle size spectra were determined using an UVP with an average pixel size of 151 µm (UVP5, ref. ⁴²), which was attached to the CTD rosette. Particle sizes in the range 128–151 µm were determined using the well-established and standard regression approach⁴², therefore we were able to detect particles with ESD between 128 and 4096 µm (the upper and lower limits of the UVP system). The lowering speed of the CTD ranged between 0.5 and 1 m s⁻¹ and data were averaged in 5 m depth bins. Particles within each depth bin were binned

<figure><img src="image_10.png" /><figcaption>**Fig. 7** Schematic illustration showing interactions between anammox bacteria and particles within the Peruvian oxygen minimum zone (OMZ).</figcaption></figure>

Aggregated organic matter formed by primary production in euphotic surface waters sinks into the anoxic OMZ (export production). Subsequently, remineralization of the particles leads to the formation of ammonium, the limiting substrate for anammox. The ammonium diffuses into the particle boundary layer, leaving a plume in the wake of the sinking particles. Free-living anammox bacteria encounter the ammonium and utilize it as a substrate to carry out anammox, leading to the formation of dinitrogen gas ($N_2$).

into size classes of 128–256, 256–512, 512–1024, 1024–2048 and 2048–4096 µm based on the pixel size of the UVP and to account for the decrease in particle abundance with increasing size. Particles in the 2048–4096 µm size bin were incredibly rare (on average 0.022 L⁻¹) so they were excluded from further analysis.

Nutrient samples were collected during the CTD rosette upcast and measured immediately on board. Nitrate and nitrite measurements were carried out spectrophotometrically on a QuAAtro autoanalyzer (Seal Analytical GmbH, Germany) with a precision of 0.1 µmol L⁻¹ (using a standard method as described by Hydes et al.⁸⁸). Ammonium concentrations were measured according to the standard fluorometric method⁸⁹.

**Anammox and denitrification rate measurements and depth-integration.** Bulk incubation experiments to determine anammox rates were conducted at 13 stations, with 6 depths per station in most cases, which were distributed between the upper boundary of the OMZ ($O_2 < 10 \mu mol L^{-1}$) down to a depth of 464 m in the OMZ core (Supplementary Tables 1 and 2). For each incubation, water was filled into 250 mL serum bottles from the Niskin bottles using gas tight iso-versinic tubing. Bottles were allowed to gently overflow two to three times their own volume and then stoppered bubble-free with 20 mm butyl rubber septa (SUPELCO) that had previously been stored under a helium atmosphere. Subsequently ¹⁵N-labelled incubations were carried out as described by Dalsgaard et al.⁹⁰. Briefly, serum bottles were bubbled with helium for 15 min and amended with 5 µM ¹⁵NO₂ after 10 min, after which they were bubbled with helium for another 5 min. Water was aliquoted through gas tight tubing into 12 mL glass vials (Exetainers, Labco Ltd), which were allowed to gently overflow with twice their volume before being capped headspace free. Caps were stored under a helium atmosphere for several months prior to use to reduce oxygen contamination⁹¹. Helium headspaces of 2 mL were introduced into each exetainer and flushed twice with helium⁶⁵. Oxygen was monitored in a subset of the incubations, three glass vials per station were prepared with an optode spot (Pyroscience, TROXSP5 Trace Oxygen Sensor Spot) and filled at random during aliquoting, in addition to the incubation vials. Oxygen concentrations were measured using a fibre-optic cable connected to an amplifier (Pyroscience, Firesting, FSO2-4). The oxygen concentrations in the samples were always below the detection limit (0.06 µmol L⁻¹). Exetainer incubation experiments were incubated in the dark at 14 ± 0.5 °C, which reflected the in situ temperature. Incubations were terminated
after approximately 0, 6, 12 18 and 24 h by addition of 100 µL of a saturated mercuric chloride solution.

At 11 additional stations (7 offshore, 4 onshore; see Supplementary Table 1b) anammox rates were determined in size-fractionated samples. Water was collected from the CTD—rotsette into 2 L glass bottles through gas tight tubing and allowed to overflow to avoid oxygen contamination. The bottle was sealed bubble free using a deoxygenated butyl rubber septa that had previously been stored under a helium atmosphere, using a needle to avoid the formation of a bubble. Subsequently, the water was degassed with helium for 15 min. The water was dispensed into 160 mL glass serum bottles using a helium overpressure through an inline 47 mm diameter filter holder (PALL, USA) which was flushed with helium to avoid oxygen contamination. For bulk water samples, no filter was added to the holder, for a second serum bottle, the filter holder contained a 10 µm pore size nylon filter (Millipore, Bellerica, MA, USA), and for the third set a 1.6 µm pore size filter was used (GF/A filter, Whatman, GE Healthcare Biosciences, Pittsburgh, PA, USA). Serum bottles were bubbled with helium for 10 min and the tracer was added (5 µM $^{15}$NO$_{2}$) and then bubbled for a further 5 min before being dispensed into 12 mL glass vials (Exetainers, LabCo Ltd.), as described above. As a precaution against external microbial contamination, all bottles and hoses were washed with HCl before use, and sterilized stoppers and needles were used for each experiment. The aliquoted samples for the size fractionation incubations were then treated the same as the incubations described above.

Isotopic ratios of $^{14}$N$^{14}$N, $^{14}$N$^{15}$N and $^{15}$N$^{15}$N dinitrogen gas were measured in the headspace of the incubation experiments with a GC-IRMS upon return to the laboratory in Bremen, Germany (customized TraceGas coupled to a multicollector IsoPrime100, Manchester, UK). The $^{15}$NO$_{2}^{-}$ labelling percentage was determined by measuring nitrite concentrations before and after label addition using the Griess method$^{86}$. $N_2$ production from denitrification and anammox was calculated from Eqs. (2) and (3), respectively, in Thamdrup et al.$^{7}$ and rates were subsequently calculated from the slope of the linear regression of the $N_2$ production as a function of time. If the slope parameter was not significantly different from zero (Student's t-test, $p < 0.05$), then rates were considered as undetectable (i.e. zero).

Areal denitrification and areal anammox rates were calculated by integrating the respective $N_2$ production over depth. Only sampling depths with oxygen concentrations below 1.5 µmol L$^{-1}$ were considered. The lowest depth of the integration was either the lower edge of the OMZ or the seafloor. In most cases $N_2$ production rates dropped to zero before that depth.

**Statistical correlation tests and linear regressions.** All statistical tests were carried out using the R software version 3.6.1$^{92}$. The Shapiro.test function was used to test the variables for normality. Spearman's rank correlation coefficients were calculated to either for correlations, using the cor.test function. Populations were compared either with Student's t-test (function t.test) or Wilcoxon signed rank test (function wilcox.test) depending on normality of the data. Linear regressions were carried out to determine the correlation of anammox rates with particle abundance (particle volume per liter) for the four size classes 128–256, 256–512, 512–1024, 1024–2048 µm. The linear model (lm) function of R was used for individual and multiple linear regressions. Figures were made using the ggplot2$^{93}$ package and Adobe Illustrator.

**Collection of sinking organic material.** Sinking organic material was collected non-destructively using a Marine Snow Catcher (OSIL, UK) lowered at 1 m s$^{-1}$. The snow catcher was closed in the oxycline or upper part of the OMZ (from oxygen concentrations of 6 µmol L$^{-1}$ down to 150 m below the oxic–anoxic interface, at water depths between 50 and 230 m, dependent on the station; see Supplementary Table 2). The waters of the upper OMZ were previously shown to have the highest N cycling activity and abundance of sinking organic material within the OMZ$^{5,64}$. After retrieval, the snow catcher was placed in an upright position in a shaded area on deck for 1–2 h to allow the sinking organic material to settle to the bottom compartment of the snow catcher and be collected in the snow catcher tray. Subsequently, the overlying water was gently removed at a flow rate of 600 mL min$^{-1}$ and discarded. After this, the snow catcher tray was transferred to the lab and individual particles and aggregates were picked by eye using a wide bore pipette and taking utmost care. The particle or aggregate was transferred to a small glass Petri dish filled with anoxic sterile filtered seawater from the same station. A scaled image was taken of each particle or aggregate for size determination (see below).

**Size determination and categorization of sinking organic material.** In order to determine the volume and type of the sinking organic material collected from the snow catcher, scaled images were taken of each particle or aggregate using a SLR-Camera (Nikon D-90) fitted with an AF-S Nikkor zoom lens ($f=16–88$ mm). Area was determined by manually identifying the edges of the organic material and integrating the inner pixels using a MATLAB script. The volume of each particle was estimated from the equivalent sphere diameter (ESD), assuming a spherical shape. Based on their visual properties, the sinking organic material was categorized manually, following ecotaxa categories where possible (https://ecotaxaobs.vlfr.fr; ref. $^{94}$). The main categories were fluffy detritus (loosely packed material), compact detritus (densely packed material), faecal pellets and other sinking

material. In addition to the ecotaxa categories, a category for loosely packed detritus containing faecal pellets was created for this dataset.

**Total organic carbon (TOC) content determination of single particles and aggregates.** TOC was determined using a combined elemental-analyser/isotopic-ratio spectrometer (EA-IRMS, Delta Plus XP coupled to a Flash Element Analyzer 112). Sixty-nine individual particles or aggregates were picked, transferred into sterile-filtered seawater for photography (as described above) and to remove surrounding material, frozen in a droplet of seawater and transported to the laboratory in Eppendorf tubes. After thawing, Eppendorf tubes were centrifuged at 100 rcf for 10 min, and most of the overlying water was removed. The pellet was resuspended in 10–20 µL of the remaining seawater and transferred into a cylindrical 5 × 9 mm tin capsule (HEKAtech GmbH, Germany). Subsequently, the remaining material was collected by rinsing the tube with 10 µL of MilliQ water twice and transferring it to the tin capsule. The $^{12}$C concentrations of the materials were usually close to the detection limit of the EA-IRMS. In order to increase the signal-intensity, we amended the sample with 4 µg of $^{15}$N$^{13}$C-labelled urea (99 atom % $^{13}$C, 98 atom % $^{15}$N, Sigma Aldrich). The samples were dried in a 60 °C oven and the tin capsules were folded into pellets. Following the same protocol, a standard series was prepared from caffeine with known $^{12}$C contents ranging between 0 and 12 µg C. Subsequently, the carbon content of the samples and standard series was measured on the EA-IRMS.

**Export production.** Export production was estimated based on both satellite imagery and UVP data. Briefly, for the satellite-based estimate chlorophyll a, temperature and PAR data for the region were downloaded from Giovanni, NASA (https://giovanni.gsfc.nasa.gov/giovanni/, last visit: 11.05.2020, ref. $^{95}$) for the time of the cruise. Primary production was calculated based on chlorophyll a concentrations using a temperature and photosynthetically active radiation (PAR)-based transfer function$^{96}$. We derived carbon-based export production using a particle export to primary production ratio (pe-ratio) determined via two methods: (1) using the global empirical relationship described in ref. $^{77}$ (based on temperature and chlorophyll a) and (2) using a pe-ratio of 0.4 as described in a previous report from the Peruvian OMZ$^{97}$. The UVP-based estimate of export of organic nitrogen out of the euphotic zone was derived for particle size classes between 128 and 2048 µm and was based on the particle flux (J) at the base of the euphotic zone and the upper boundary of the OMZ (similar to refs. $^{46,98}$):

$$J = c_f \int_{r_0}^{r_f} n(r)m(r)U(r)dr \quad (1)$$

where $n(r)$ is the particle concentration, $m(r)$ the size-specific carbon content of the particles, $U(r)$ the particle settling velocity, $r$ is the sphere-equivalent radius of the particles, $r_0$ and $r_f$ are the lower and upper limits for integration, respectively, and $c_f$ is the conversion factor that relates carbon to nitrogen export: $c_f = 1$ for carbon export production and $c_f = 0.14$ (based on a previously defined C:N ratio of 7.2 in sinking organic matter$^{5}$) for nitrogen export. The measured carbon content (Supplementary Fig. 14) of individual particles was used to rescale the power-law function described by Alldredge 1998 (ref. $^{99}$):

$$m(r) = 8.9V(r)^{0.52} \quad (2)$$

where $m(r)$ is the carbon content (µg agg$^{-1}$), $V(r)$ is the total volume of the particle (mm$^3$) which was estimated assuming spherical shapes of the particles: $V(r) = 4/3\pi r^3$. The particle settling velocity was calculated based on Stoke's law assuming that in average the particles are well represented by a spherical shape:

$$U(r) = \frac{2r^2g \Delta \rho}{9\nu \rho} \quad (3)$$

where $g$ is the acceleration by gravity (9.81 m s$^{-2}$), $\nu$ the kinematic viscosity of seawater and $\rho$ the density of sea-water. $\nu$ and $\rho$ were estimated based on the CTD profiles at the specific depths of the particles. The excess density of the particles was calculated as

$$\Delta \rho = (1 - \theta)\rho_s + \rho(\theta - 1) \quad (4)$$

where $\rho_s$ is the solid-hydrated density (1150 kg m$^{-3}$ for diatom aggregates$^{100}$). The fractal characteristic of marine aggregates leads to a size-dependent porosity$^{100}$ ($r$ in mm):

$$\theta(r) = -2.6 \cdot 10^{-3} r^{-1.6} + 1 \quad (5)$$

**Nitrogen release and encounter rates.** In order to investigate the relationship between particle abundance and anammox rates, the UVP data was used to estimate the release of N from particles and the encounter rates of anammox bacteria with the released N. During the descent through the oxic and anoxic water column, ammonium is assumed to be released by remineralization of the organics bound in the particles. We assume that the particle turnover follows a first-order kinetic: $dm/dt = -R_{rem}m$, where $R_{rem}$ is the carbon-specific degradation rate (day$^{-1}$) and $t$ the time variable. Assuming a constant settling velocity of the particle through the water-column, $R_{rem}$ can be transformed into a length-scale: $L_{rem} = U/R_{rem}$, the so-called remineralization length-scale. When production, remineralization as well as
fragmentation are in balance, the flux attenuation follows an exponential decay:

$$J(z) = J_0 \cdot e^{\left(\frac{z-z_0}{L_{rem}}\right)} \quad (6)$$

where z is the depth and z₀ is a reference depth. This equation illustrates the important meaning of L<sub>rem</sub>. Substituting z = L<sub>rem</sub> immediately yields J(L<sub>rem</sub>)/J₀ = 0.37, implying that a particle decayed to 37% of its initial value after travelling a distance of L<sub>rem</sub>. De Soto et al. (2018) ref. 101 derived an extended model for the flux attenuation by introducing the parameter z<sub>β</sub> which has the units of depth (m) and can be interpreted as a depth where the degradation of particles decays:

$$J(z) = J_0 \cdot e^{\left[ -\frac{z_0}{L_{rem}} \left( e^{-\frac{z_0}{J_0}} - \left( e^{-\frac{z_0}{J_0}} \right)^{\frac{J_0}{z_0}} \right) \right]} \quad (7)$$

We estimated the remineralization length scale by fitting Eqs. (6) and (7) to the carbon flux profiles (Supplementary Fig. 13). Values for L<sub>rem</sub> in the core of the OMZ were found to range between 139 and 305 m with an average of 187 m (Supplementary Table 7) based on Eq. (6) and in between 57 and 195 m with an average of 107 m based on Eq. (7) (Supplementary Table 8). Applying Eq. (3) for the settling velocity, the carbon-specific degradation rates were estimated to range between 0.13 and 0.28 day<sup>-1</sup> with an average of 0.18 day<sup>-1</sup> for Eq. (6) and between 0.1 and 0.54 day<sup>-1</sup> with an average of 0.32 day<sup>-1</sup> for Eq. (7). Overall, L<sub>rem</sub> as well as R<sub>rem</sub> were in the same range and the models predicted similar trends. Henceforth, we will focus on the estimates based on Eq. (6) as those represent a more conservative range. Based on L<sub>rem</sub>, the N released by particles while settling through the water column was estimated as

$$NR = c_f \int_{r_0}^{r_1} n(r)m(r)R_{\text{rem}}(r)dr \quad (8)$$

Similar to the remineralization length scale, the diffusion length scale indicates the travelled distance of a particle before solutes are produced in the interior of particles, are released to the water column: L<sub>diff</sub> = U/R<sub>diff</sub>, where R<sub>diff</sub> is the diffusion timescale estimated as R<sub>diff</sub> = D/(r + δ)<sup>2</sup> with D the diffusion coefficient.

Encounters between the free-living anammox bacteria and the nitrogen that is released into the diffusive boundary layer of the particle were calculated. Therefore, we integrated the diffusive boundary layer area perpendicular to the settling vector over the distance travelled at a given time. Assuming that anammox bacteria are passive tracers in the water column and that the bacteria only need to be near the boundary layer to take up the released N, the encounter rate can be calculated as

$$ER = c_s \int_{r_0}^{r_1} n(r)\pi(2r\delta + \delta^2)U(r)dr \quad (9)$$

where c<sub>s</sub> is the concentration of anammox bacteria (4.3 × 10<sup>4</sup> cells mL<sup>-1</sup>, ref. 41) in a given water volume. δ = r/Sh is the thickness of the diffusive boundary layer which is estimated based on the Sherwood number Sh, Sh = 1 + 0.62·Re<sup>0.41</sup>Sc<sup>0.33</sup>, where Re = Ur/v is the Reynolds number and Sc = v/D the Schmidt number following (refs. 70,102). All integrations were performed numerically in Matlab (2018b Mathworks).

## Data availability

Nutrient data from M136 are deposited on Pangaea (https://doi.org/10.1594/PANGAEA.904404) in supplement to ref. 85. Nutrient data from M138 are deposited on Pangaea (https://doi.org/10.1594/PANGAEA.914948). UVP data are available upon request to R. Kiko and will be made available in a data collection to be found at https://doi.org/10.1594/PANGAEA.924375. The UVP data are also freely available via the particle module of ecotaxa at https://ecotaxaobs-vlfri.fr/part/ (particle projects 25, 95, 96, 97 and 98; please note that usage of ecotaxa requires user registration). The satellite imagery data used in this study was downloaded through the Giovanni online data system, developed and maintained by the NASA GES DISC. All other data are provided within the text or the supplementary material.

Received: 17 August 2020; Accepted: 26 April 2021;
Published online: 28 May 2021

## References

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

2. Codispoti, L. A. et al. The oceanic fixed nitrogen and nitrous oxide budgets: moving targets as we enter the anthropocene? Sci. Mar. 65, 85–105 (2001).

3. Gruber, N. & Sarmiento, J. L. Global patterns of marine nitrogen fixation and denitrification. Glob. Biogeochem. Cycles 11, 235–266 (1997).

4. Gruber, N. In The Ocean Carbon Cycle and Climate (eds Follows, M. & Oguz, T.) 97–148 (Springer Netherlands, 2004).

5. Kalvelage, T. et al. Nitrogen cycling driven by organic matter export in the South Pacific oxygen minimum zone. Nat. Geosci. 6, 228–234 (2013).

6. Peters, B. D. et al. Vertical modeling of the nitrogen cycle in the eastern tropical South Pacific oxygen deficient zone using high-resolution

concentration and isotope measurements. Glob. Biogeochem. Cycles 30, 1661–1681 (2016).

7. Thamdrup, B. et al. Anaerobic ammonium oxidation in the oxygen-deficient waters off northern Chile. Limnol. Oceanogr. 51, 2145–2156 (2006).

8. Dalsgaard, T., Thamdrup, B., Farías, L. & Revsbech, N. P. Anammox and denitrification in the oxygen minimum zone of the eastern South Pacific. Limnol. Oceanogr. 57, 1331–1346 (2012).

9. Babbin, A. R., Keil, R. G., Devol, A. H. & Ward, B. B. Organic matter stoichiometry, flux, and oxygen control nitrogen loss in the ocean. Science 344, 406–408 (2014).

10. Ward, B. B. et al. Denitrification as the dominant nitrogen loss process in the Arabian Sea. Nature 461, 78–81 (2009).

11. Schmidtko, S., Stramma, L. & Visbeck, M. Decline in global oceanic oxygen content during the past five decades. Nature 542, 335–339 (2017).

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

13. Oschlies, A., Brandt, P., Stramma, L. & Schmidtko, S. Drivers and mechanisms of ocean deoxygenation. Nat. Geosci. 11, 467–473 (2018).

14. Breitburg, D. et al. Declining oxygen in the global ocean and coastal waters. Science 359, 6371 (2018).

15. Cocco, V. et al. Oxygen and indicators of stress for marine life in multi-model global warming projections. Biogeosciences 10, 1849–1868 (2013).

16. Stramma, L. et al. Expansion of oxygen minimum zones may reduce available habitat for tropical pelagic fishes. Nat. Clim. Chang. 2, 33–37 (2012).

17. Frölicher, T. L., Rodgers, K. B., Stock, C. A. & Cheung, W. W. L. Sources of uncertainties in 21st century projections of potential ocean ecosystem stressors. Glob. Biogeochem. Cycles 30, 1224–1243 (2016).

18. Bopp, L. et al. Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models. Biogeosciences 10, 6225–6245 (2013).

19. Capone, D. G. & Hutchins, D. A. Microbial biogeochemistry of coastal upwelling regimes in a changing ocean. Nat. Geosci. 6, 711–717 (2013).

20. Bianchi, D., Weber, T. S., Kiko, R. & Deutsch, C. Global niche of marine anaerobic metabolisms expanded by particle microenvironments. Nat. Geosci. 11, 263–268 (2018).

21. Niemeyer, D., Kriest, I. & Oschlies, A. The effect of marine aggregate parameterisations on nutrients and oxygen minimum zones in a global biogeochemical model. Biogeosciences 16, 3095–3111 (2019).

22. Woebken, D., Fuchs, B. M., Kuypers, M. M. M. & Amann, R. Potential interactions of particle-associated anammox bacteria with bacterial and archaeal partners in the Namibian upwelling system. Appl. Environ. Microbiol. 73, 4648–4657 (2007).

23. Ganesh, S. et al. Size-fraction partitioning of community gene transcription and nitrogen metabolism in a marine oxygen minimum zone. ISME J. 9, 2682–2696 (2015).

24. Suter, E. A., Pachiadaki, M., Taylor, G. T., Astor, Y. & Edgcomb, V. P. Free-living chemoautotrophic and particle-attached heterotrophic prokaryotes dominate microbial assemblages along a pelagic redox gradient. Environ. Microbiol. 20, 693–712 (2018).

25. Fuchsman, C. A., Devol, A. H., Saunders, J. K., McKay, C. & Rocap, G. Niche partitioning of the N cycling microbial community of an offshore oxygen deficient zone. Front. Microbiol. 8, 2384 (2017).

26. Ganesh, S., Parris, D. J., DeLong, E. F. & Stewart, F. J. Metagenomic analysis of size-fractionated picoplankton in a marine oxygen minimum zone. ISME J. 8, 187–211 (2014).

27. Riley, J. S. et al. The relative contribution of fast and slow sinking particles to ocean carbon export. Glob. Biogeochem. Cycles 26, 1 (2012).

28. Walker, B. D., Beupré, S. R., Guilderson, T. P., McCarthy, M. D. & Druffel, E. R. M. Pacific carbon cycling constrained by organic matter size, age and composition relationships. Nat. Geosci. 9, 888–891 (2016).

29. Benner, R. & Amon, R. M. W. The size-reactivity continuum of major bioelements in the ocean. Annu. Rev. Mar. Sci. 7, 185–205 (2015).

30. 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).

31. Iversen, M. H. & Ploug, H. Ballast minerals and the sinking carbon flux in the ocean: carbon-specific respiration rates and sinking velocity of marine snow aggregates. Biogeosciences 7, 2613–2624 (2010).

32. Weber, T., Cram, J. A., Leung, S. W., DeVries, T. & Deutsch, C. Deep ocean nutrients imply large latitudinal variation in particle transfer efficiency. Proc. Natl Acad. Sci. USA 113, 8606–8611 (2016).

33. Maerz, J., Six, K. D., Stemmler, I., Ahmerkamp, S. & Ilyina, T. Microstructure and composition of marine aggregates as co-determinants for vertical particulate organic carbon transfer in the global ocean. Biogeosciences 17, 1765–1803 (2020).

34. Callbeck, C. M., Lavik, G., Stramma, L., Kuypers, M. M. M. & Bristow, L. A. Enhanced nitrogen loss by eddy-induced vertical transport in the offshore Peruvian oxygen minimum zone. PLoS ONE 12, e0170059 (2017).
35. Pennington, J. T. et al. Primary production in the eastern tropical Pacific: a review. *Prog. Oceanogr.* **69**, 285–317 (2006).

36. Rodríguez-Morata, C., Diaz, H. F., Ballesteros-Canovas, J. A., Rohrer, M. & Stoffel, M. The anomalous 2017 coastal El Niño event in Peru. *Clim. Dyn.* **52**, 5605–5622 (2019).

37. Garreaud, R. D. A plausible atmospheric trigger for the 2017 coastal El Niño. *Int. J. Climatol.* **38**, e1296–e1302 (2018).

38. Fi, O.-G. Coastal el niño 2017 or simply: the carnival coastal warming event? *MOJES* **2**, 8 (2017).

39. Echevin, V. et al. Forcings and evolution of the 2017 coastal El Niño off northern Peru and Ecuador. *Front. Mar. Sci.* **5**, 367 (2018).

40. Callbeck, C. M. et al. Oxygen minimum zone cryptic sulfur cycling sustained by offshore transport of key sulfur oxidizing bacteria. *Nat. Commun.* **9**, 1729 (2018).

41. Hamersley, M. R. et al. Anaerobic ammonium oxidation in the Peruvian oxygen minimum zone. *Limnol. Oceanogr.* **52**, 923–933 (2007).

42. 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).

43. Sheldon, R. W., Prakash, A. & Sutcliffe, W. H. The size distribution of particles in the ocean1. *Limnol. Oceanogr.* **17**, 327–340 (1972).

44. Stemmann, L. & Boss, E. Plankton and particle size and packaging: from determining optical properties to driving the biological pump. *Annu. Rev. Mar. Sci.* **4**, 263–290 (2012).

45. Kiko, R. et al. Biological and physical influences on marine snowfall at the equator. *Nat. Geosci.* **10**, 852–858 (2017).

46. Guidi, L. et al. Relationship between particle size distribution and flux in the mesopelagic zone. *Deep Sea Res. Part I* **55**, 1364–1374 (2008).

47. Baker, C. A. et al. Slow-sinking particulate organic carbon in the Atlantic Ocean: magnitude, flux, and potential controls. *Glob. Biogeochem. Cycles* **31**, 1051–1065 (2017).

48. Engel, A., Endres, S., Galgani, L. & Schartau, M. Marvelous marine microalg: on the distribution and impact of gel-like particles in the oceanic water-column. *Front. Mar. Sci.* **7**, 405 (2020).

49. Iversen, M. H. & Lampitt, R. S. Size does not matter after all: no evidence for a size-sinking relationship for marine snow. *Prog. Oceanogr.* **102445** (2020).

50. Cavan, E. L., Giering, S. L. C., Wolff, G. A., Trimmer, M. & Sanders, R. Alternative particle formation pathways in the eastern tropical north pacific's biological carbon pump. *J. Geophys. Res. Biogeosci.* **123**, 2198–2211 (2018).

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

52. Kiko, R. et al. Zooplankton-mediated fluxes in the Eastern Tropical North Atlantic. *Front. Mar. Sci.* **7**, 353 (2020).

53. Bianchi, D., Babbin, A. R. & Galbraith, E. D. Enhancement of anammox by the excretion of diel vertical migrators. *Proc. Natl Acad. Sci. USA* **111**, 15653–15658 (2014).

54. Galán, A., Faúndez, J., Thamdrup, B., Santibáñez, J. F. & Farias, L. Temporal dynamics of nitrogen loss in the coastal upwelling ecosystem off central Chile: evidence of autotrophic denitrification through sulfide oxidation. *Limnol. Oceanogr.* **59**, 1865–1878 (2014).

55. Lam, P. et al. Revising the nitrogen cycle in the Peruvian oxygen minimum zone. *Proc. Natl Acad. Sci. USA* **106**, 4752–4757 (2009).

56. Galán, A. et al. Anammox bacteria and the anaerobic oxidation of ammonium in the oxygen minimum zone off northern Chile. *Deep Sea Res. Part II* **56**, 1021–1031 (2009).

57. Ward, B. B. et al. Organic carbon, and not copper, controls denitrification in oxygen minimum zones of the ocean. *Deep Sea Res. Part I* **55**, 1672–1683 (2008).

58. Fuchsman, C. A., Staley, J. T., Oakley, B. B., Kirkpatrick, J. B. & Murray, J. W. Free-living and aggregate-associated Planctomycetes in the Black Sea. *FEMS Microbiol. Ecol.* **80**, 402–416 (2012).

59. Lavik, G. et al. Detoxification of sulphidic African shelf waters by blooming chemolithotrophs. *Nature* **457**, 581–584 (2009).

60. Callbeck, C. M. et al. *Arcobacter peruensis* sp. nov., a chemolithoheterotroph isolated from sulfide- and organic-rich coastal waters off Peru. *Appl. Environ. Microbiol.* **85**, 24 (2019).

61. Hach, P. F. et al. Rapid microbial diversification of dissolved organic matter in oceanic surface waters leads to carbon sequestration. *Sci. Rep.* **10**, 13025 (2020).

62. Martin, J. H., Knauer, G. A., Karl, D. M. & Broenkow, W. W. VERTEX: carbon cycling in the northeast Pacific. *Deep Sea Res. Part A* **34**, 267–285 (1987).

63. Van Mooy, B. A., 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).

64. Kalvelage, T. et al. Aerobic microbial respiration in oceanic oxygen minimum zones. *PLoS ONE* **10**, e0133526 (2015).

65. De Brabandere, L. et al. Vertical partitioning of nitrogen-loss processes across the oxic-anoxic interface of an oceanic oxygen minimum zone. *Environ. Microbiol.* **16**, 3041–3054 (2014).

66. Garcia-Robledo, E. et al. Cryptic oxygen cycling in anoxic marine zones. *Proc. Natl Acad. Sci. USA* **114**, 8319–8324 (2017).

67. Srain, B. M. et al. Fermentation and anaerobic oxidation of organic carbon in the oxygen minimum zone of the upwelling ecosystem off concepción, in Central Chile. *Front. Mar. Sci.* **7**, 533 (2020).

68. Jackson, G. A. Seascapes: the world of aquatic organisms as determined by their particulate natures. *J. Exp. Biol.* **215**, 1017–1030 (2012).

69. Zetsche, E., Larsson, A. I., Iversen, M. H. & Ploug, H. Flow and diffusion around and within diatom aggregates: effects of aggregate composition and shape. *Limnol. Oceanogr.* **65**, 8 (2020).

70. Kiørboe, T. Formation and fate of marine snow: small-scale processes with large- scale implications. *Sci. Mar.* **65**, 57–71 (2001).

71. Kiørboe, T. & Jackson, G. A. Marine snow, organic solute plumes, and optimal chemosensory behavior of bacteria. *Limnol. Oceanogr.* **46**, 1309–1318 (2001).

72. Stocker, R., Seymour, J. R., Samadani, A., Hunt, D. E. & Polz, M. F. Rapid chemotactic response enables marine bacteria to exploit ephemeral microscale nutrient patches. *Proc. Natl Acad. Sci. USA* **105**, 4209–4214 (2008).

73. Stocker, R. Marine microbes see a sea of gradients. *Science* **338**, 628–633 (2012).

74. Strous, M. et al. Missing lithotroph identified as new planctomycete. *Nature* **400**, 446–449 (1999).

75. Kuypers, M. M. M. et al. Anaerobic ammonium oxidation by anammox bacteria in the Black Sea. *Nature* **422**, 608–611 (2003).

76. van Niftrik, L. & Jetten, M. S. M. Anaerobic ammonium-oxidizing bacteria: unique microorganisms with exceptional properties. *Microbiol. Mol. Biol. Rev.* **76**, 585–596 (2012).

77. Zhao, R. et al. Geochemical transition zone powering microbial growth in subsurface sediments. *Proc. Natl. Acad. Sci.* **117**, 51 (2020).

78. Woebken, D. et al. A microdiversity study of anammox bacteria reveals a novel Candidatus Scalindua phylotype in marine oxygen minimum zones. *Environ. Microbiol.* **10**, 3106–3119 (2008).

79. Lambert, B. S., Fernandez, V. I. & Stocker, R. Motility drives bacterial encounter with particles responsible for carbon export throughout the ocean. *Limnol. Oceanogr.* **4**, 113–118 (2019).

80. Codispoti, L. A. & Packard, T. T. Denitrification rates in the eastern tropical South Pacific. *J. Mar. Res.* **38**, 3 (1980).

81. Deutsch, C., Gruber, N., Key, R. M., Sarmiento, J. L. & Ganachaud, A. Denitrification and N₂ fixation in the Pacific Ocean. *Glob. Biogeochem. Cycles* **15**, 483–506 (2001).

82. DeVries, T., Deutsch, C., Primeau, F., Chang, B. & Devol, A. Global rates of water-column denitrification derived from nitrogen gas measurements. *Nat. Geosci.* **5**, 547–550 (2012).

83. Lombard, F. et al. Globally consistent quantitative observations of planktonic ecosystems. *Front. Mar. Sci.* **6**, 196 (2019).

84. Giering, S. L. C. et al. Sinking organic particles in the ocean—flux estimates from in situ optical devices. *Front. Mar. Sci.* **6**, 834 (2020).

85. Lüdke, J. et al. *Influence of Intraseasonal Eastern Boundary Circulation Variability on Hydrography and Biogeochemistry off Peru* (PANGAEA, 2019).

86. Grasshoff, K., Kremling, K. & Ehrhardt, M. *Methods of Seawater Analysis* (Wiley, 1999).

87. Loginova, A. N., Thomsen, S. & Engel, A. Chromophoric and fluorescent dissolved organic matter in and above the oxygen minimum zone off Peru. *J. Geophys. Res. Oceans* **121**, 7973–7990 (2016).

88. Hydes, D. et al. in *The GO-SHIP Repeat Hydrography Manual: A Collection of Expert Reports and Guidelines*. IOCCP Report No. 14, ICPO Publication Series No. 134, Version 1, 2010 (eds Hood, E. M., C. L. Sabine & B. M. Sloyan) (UNESCO/IOC, 2010).

89. Holmes, R. M., Aminot, A., Kérouel, R., Hooker, B. A. & Peterson, B. J. A simple and precise method for measuring ammonium in marine and freshwater ecosystems. *Can. J. Fish. Aquat. Sci.* **56**, 1801–1808 (1999).

90. Dalsgaard, T., Canfield, D. E., Petersen, J., Thamdrup, B. & Acuña-González, J. N2 production by the anammox reaction in the anoxic water column of Golfo Dulce, Costa Rica. *Nature* **422**, 606–608 (2003).

91. De Brabandere, L., Thamdrup, B., Revsbech, N. P. & Foadi, R. A critical assessment of the occurrence and extend of oxygen contamination during anaerobic incubations utilizing commercially available vials. *J. Microbiol. Methods* **88**, 147–154 (2012).

92. Team, R. C. & DC, R. *A Language and Environment For Statistical Computing* (R Foundation for Statistical Computing, 2019).

93. Wickham, H. *ggplot2—Elegant Graphics for Data Analysis* (Springer International Publishing, 2016).

94. Picheral, M., Colin, S. & Irisson, J. O. *EcoTaxa, A Tool for the Taxonomic Classification of Images* http://ecotaxa OBS-vlfvr.fr (2017).
95. Acker, J. G. & Leptoukh, G. Online analysis enhances use of NASA earth science data. *EOS Trans. AGU* **88**, 14 (2007).

96. Behrenfeld, M. J. & Falkowski, P. G. Photosynthetic rates derived from satellite-based chlorophyll concentration. *Limnol. Oceanogr.* **42**, 1–20 (1997).

97. Dunne, J. P., Armstrong, R. A., Gnanadesikan, A. & Sarmiento, J. L. Empirical and mechanistic models for the particle export ratio. *Glob. Biogeochem. Cycles* **19**, 4 (2005).

98. Iversen, M. H., Nowald, N., Ploug, H., Jackson, G. A. & Fischer, G. High resolution profiles of vertical particulate organic matter export off Cape Blanc, Mauritania: degradation processes and ballasting effects. *Deep Sea Res. Part I* **57**, 771–784 (2010).

99. Alldredge, A. The carbon, nitrogen and mass content of marine snow as a function of aggregate size. *Deep Sea Res. Part I* **45**, 529–541 (1998).

100. Alldredge, A. L. & Gotschalk, C. In situ settling behavior of marine snow. *Limnol. Oceanogr.* **33**, 339–351 (1988).

101. De Soto, F., Ceballos-Romero, E. & Villa-Alfageme, M. A microscopic simulation of particle flux in ocean waters: application to radioactive pair disequilibrium. *Geochim. Cosmochim. Acta* **239**, 136–158 (2018).

102. Ploug, H., Hietanen, S. & Kuparinen, J. Diffusion and advection within and around sinking, porous diatom aggregates. *Limnol. Oceanogr.* **47**, 1129–1136 (2002).

## Acknowledgements

We thank Gabriele Klockgether and Philipp Hach for mass spectrometry and elemental analysis, Matija Lagator for help with C content measurements, Angelina Klett for help with particle size determination, Gerd Krahmann and Hermann Bange for CTD and nutrient data, Markus Dengler and Hermann Bange for their support during M136 and M138 respectively, and Simon Ramondenc for discussion of the multiple linear regression. We thank the captain and crew of R/V Meteor for their support during cruises M136 and M138. Further, we thank Bennett Lambert for discussions on encounter rates. This work was supported by the Max Planck Society and the DFG-funded collaborative research center 754 "Climate-biogeochemistry interactions in the tropical Ocean" (Work Packages B4, B3 and B8). R.K. and H.H. were furthermore supported by the "CUSCO—Coastal Upwelling System in a Changing Ocean" project (WP5) funded by the Federal Ministry of Education and Research (Germany). R.K. also acknowledges support via a "Make Our Planet Great Again" grant of the French National Research Agency within the "Programme d'Investissements d'Avenir"; reference "ANR-19-MPGA-0012". S.A. and J.M. acknowledge funding from the Max Planck Society (MPG) for the "Multiscale Approach on the Role of Marine Aggregates" (MARMA) project. H.K.M. and M.H. I. received funding from the DFG under Germany's Excellence Strategy (no. EXC-2077-390741603). M.H.I. was supported by the Alfred Wegener Institute Helmholtz Center for Polar and Marine Research and the HGF Young Investigator Group SeaPump "Seasonal and regional food web interactions with the biological pump": VH-NG-1000.

## Author contributions

C.K., L.A.B., S.A. and M.M.M.K. designed study, C.K., L.A.B. and G.L., collected samples and conducted experiments, C.K., S.A. and H.K.M. analysed and visualized data, R.K. and H.H. collected and processed UVP data, J.M. and M.H.I. contributed to data interpretation and organic carbon measurements, S.A. developed model, processed UVP, CTD and satellite data and extrapolated data, C.K., S.A., H.K.M. and M.M.M.K. wrote the manuscript with contributions from all co-authors.

## 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/s41467-021-23340-4.

**Correspondence** and requests for materials should be addressed to S.A. or H.K.M.

**Peer review information** *Nature Communications* thanks Andrew Babbin, Jeremy Rich, and other, anonymous, reviewers for their contributions to the peer review of this work.

**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_11.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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021