Abstract
Zooplankton are important consumers of sinking particles in the ocean’s twilight zone. However, the impact of different taxa depends on their feeding mode. In contrast to typical suspension-feeding zooplankton, flux-feeding taxa preferentially consume rapidly sinking particles that would otherwise penetrate into the deep ocean. To quantify the potential impact of two flux-feeding zooplankton taxa [Aulosphaeridae (Rhizaria), and Limacina helicina (euthecosome pteropod)] and the total suspension-feeding zooplankton community, we measured depth-stratified abundances of these organisms during six cruises in the California Current Ecosystem. Using allometric–scaling relationships, we computed the percentage of carbon flux intercepted by flux feeders and suspension feeders. These estimates were compared to direct measurements of carbon flux attenuation (CFA) made using drifting sediment traps and 238U–234Th disequilibrium. We found that CFA in the shallow twilight zone typically ranged from 500 to 1000 μmol organic C flux remineralized per 10-m vertical depth bin. This equated to approximately 6–10% of carbon flux remineralized/10 m. The two flux-feeding taxa considered in this study could account for a substantial proportion of this flux near the base of the euphotic zone. The mean flux attenuation attributable to Aulosphaeridae was 0.69%/10 m (median = 0.21%/10 m, interquartile range = 0.04–0.81%) at their depth of maximum abundance (∼100 m), which would equate to ∼10% of total flux attenuation in this depth range. The maximum flux attenuation attributable to Aulosphaeridae reached 4.2%/10 m when these protists were most abundant. L. helicina, meanwhile, could intercept 0.45–1.6% of carbon flux/10 m, which was slightly greater (on average) than the Aulosphaeridae. In contrast, suspension-feeding zooplankton in the mesopelagic (including copepods, euphausiids, appendicularians, and ostracods) had combined clearance rates of 2–81 L m-3 day-1 (mean of 19.6 L m-3 day-1). This implies a substantial impact on slowly sinking particles, but a negligible impact on the presumably rapidly sinking fecal pellets that comprised the majority of the material collected in sediment traps. Our results highlight the need for a greater research focus on the many taxa that potentially act as flux feeders in the oceanic twilight zone.
Introduction
Zooplankton play diverse roles in the cycling of many elements in the ocean including iron, zinc, sulfur, and mercury (; ; ; Schmidt et al., 2016; Gorokhova et al., 2018). However, their greatest importance to global biogeochemistry is likely derived from their roles in the biological carbon pump (BCP; ; Turner, 2015; Steinberg and Landry, 2017). The BCP refers to the processes that transport organic carbon fixed by phytoplankton in the euphotic zone into the mesopelagic realm (Silver and Gowing, 1991; ; Siegel et al., 2016). The BCP leads to net transport of CO2 from the surface ocean into the deep ocean where it can be sequestered for periods of decades to millennia (). Estimates of the present magnitude of the BCP range from 5 to 12 Pg C year-1 (Henson et al., 2011; Laws et al., 2011; Siegel et al., 2014; ), however, the responses of mesozooplankton and the BCP to future climate change remain unknown.
Research on the role of mesozooplankton in the BCP has focused primarily on the epipelagic zone where the relationship between zooplankton and the BCP can change based on the community composition of phytoplankton and zooplankton. Zooplankton can play an important role in combining smaller particles into large, rapidly sinking fecal pellets (; Komar et al., 1981; Turner, 2002; Wilson et al., 2008) and promote aggregation and sinking through discarded mucous feeding webs (; Hansen et al., 1996; Robison et al., 2005). Zooplankton can also decrease the magnitude of the biological pump when their grazing activities exert top-down control on phytoplankton production (; Goericke, 2002) or fragment larger particles into smaller ones (), and also play an important role in nutrient regeneration in the euphotic zone (; ; Saba et al., 2011). Active transport by diel vertically migrating zooplankton is also an important (potentially dominant) component of the biological pump in many marine ecosystems (Steinberg et al., 2000; Hannides et al., 2009; ; Stukel et al., 2018b; ; Hernández-Leon et al., unpublished; Kelly et al., unpublished; Kiko et al., unpublished).
Although fewer studies have quantified zooplankton impacts in the mesopelagic, zooplankton play a substantial role in consuming, disaggregating, and transforming sinking particles (Steinberg and Landry, 2017). These organisms may play a substantial role in modulating marine snow flux in the mesopelagic (Lampitt et al., 1993) and have been hypothesized to play roles as “gatekeepers” that modulate carbon transfer from the euphotic zone to the mesopelagic (Jackson and Checkley, 2011). However, zooplankton are phylogenetically and functionally diverse with a wide array of feeding strategies (Kiørboe, 2011). It is thus important to consider the relative importance of zooplankton with different feeding modes. For instance, suspension-feeding salps and crustaceans may consume particles in relative proportion to their abundance in the water column, although they may show selectivity based on particle size or other characteristics (). In contrast, a flux-feeding pteropod will intercept particles in proportion to the speed with which they sink through the water column (Jackson et al., 1993), while cruise-feeding zooplankton may search for and colonize large aggregates (Kiørboe and Thygesen, 2001). Yet other taxa may break aggregates apart due to their swimming and feeding behaviors or partially consume fecal pellets leading to breakage and decreased settling velocities (Goldthwait et al., 2004; Iversen and Poulsen, 2007). Understanding these interactions is important to understanding observed decadal scale changes in carbon flux attenuation (CFA) in the mesopelagic (Lomas et al., 2010).
In this study we investigate the relative importance of representatives of two groups of flux-feeding zooplankton (phaeodarians and euthecosome pteropods) with respect to sinking particle flux attenuation. Phaeodarians are a group of siliceous protists from the supergroup Rhizaria that typically thrive in the deep ocean (Nakamura and Suzuki, 2015). We focus in this manuscript on one group of large Phaeodaria (Aulosphaeridae) that have a typical diameter of ∼2-mm and are common in the shallow twilight zone in the CCE (Ohman et al., 2012; , ). These mesopelagic phaeodarians are likely flux feeders that have relatively slow (for protists) growth rates and rely on the rain of sinking particles from above for their food (Gowing, 1986, 1989; Gowing and Bentham, 1994; Stukel et al., 2018a). Euthecosome pteropods are pelagic molluscs that produce large mucous feeding webs to trap food, including swimming zooplankton and sinking particles (; Lalli and Gilmer, 1989). We focus on the species Limacina helicina, which is common in cold waters worldwide from the Antarctic to the Arctic including the California Current (Hunt et al., 2010; ). L. helicina produces a mucous feeding web with a typical diameter of 40–55 mm (). We draw our data from a decade’s worth of field campaigns with differing objectives and hence changing methodology. We thus do not intend this to be a definitive assessment of the role of flux feeders in the mesopelagic ecosystem. Rather we intend it as an initial quantitative investigation of the differing roles of suspension feeders and two specific taxa of flux feeders. We find that each of these flux feeders has the potential to mediate a substantial portion of the CFA in the shallow twilight zone, although their impact is greatly reduced in the deeper mesopelagic. In contrast, suspension-feeding zooplankton are abundant through the mesopelagic, but likely only play a substantial role in the attenuation of the flux of slowly sinking particles. There is, however, substantial uncertainty around our core conclusions, and we hope that this uncertainty will spur future targeted investigations.
Materials and Methods
Field Sampling
In situ measurements were made on six process cruises of the CCE LTER Program (P0704, April 2007; P0810, October 2008; P1106, June 2011; P1208, August 2012; P1408, August 2014; and P1604 April 2016, Figure 1). On these cruises, we used a quasi-Lagrangian sampling scheme to track water parcels for a period of 2–5 days while quantifying biotic and abiotic standing stocks and rates (Landry et al., 2009, 2012). After preliminary site surveys with a free-fall Moving Vessel Profiler (Ohman et al., 2012), quasi-Lagrangian experiments (hereafter referred to as “cycles”) were initiated with the deployment of a surface-tethered drifting sediment trap with a 3 × 1-m drogue centered at 15-m depth to track the mixed layer (Stukel et al., 2013). A second, identically drogued, experimental array was deployed and recovered daily while being used as a platform for in situ incubations (Landry et al., 2009). During each cycle, paired day–night Multiple Opening and Closing Net and Environmental Sensing System (MOCNESS) net tows were used to determine vertical patterns of mesozooplankton abundance (Powell and Ohman, 2015). On approximately 10 CTD-Niskin rosette casts per cycle, we used an Underwater Vision Profiler (UVP5) to determine vertical profiles of rhizarians (Ohman et al., 2012; ). We also measured the vertical flux (and flux attenuation) of sinking particles using the aforementioned drifting sediment traps and measurements of water column 238U–234Th deficiency (Stukel et al., 2015). We divided our Lagrangian experiments into oligotrophic cycles (<0.5 μg Chl a L-1) or high biomass cycles (>0.5 μg Chl a L-1). For details on cruise conditions, see Supplementary Materials.
FIGURE 1
Zooplankton Collection and Enumeration
Mesozooplankton depth-stratified abundances were determined from day–night MOCNESS tows (1 m2 net, 202-μm mesh) on the P0704 and P0810 cruises (typically two pairs of tows per cycle). Each net collected organisms over a ∼50-m depth interval with consecutive nets tripped from a depth of 450 m to the surface. Samples were preserved in 1.8% formaldehyde and analyzed using a ZooScan digital scanner with ZooProcess software (Gorsky et al., 2010). Zooplankton were sorted into broad taxonomic groups (e.g., copepods, euphausiids, doliolids, etc.) using machine learning algorithms, the taxonomic assignments of 100% of the vignettes validated manually, then automatically sized (as Feret diameter). Organismal length was then used with allometric equations (see below) to quantify organism biomass and clearance rate. For additional details on ZooScan processing, see Stukel et al. (2013), Powell and Ohman (2015), and Ohman and Romagnan (2016). On two cycles from the P1208 cruise, these samples were also sorted to enumerate the abundance of pteropods, specifically L. helicina ().
Because they are not well preserved in net tows, large rhizarians (>600-μm) were quantified using an UVP5 (Picheral et al., 2010; ). The UVP5 is an in situ imaging camera that was mounted downward facing on the bottom of the ship’s CTD-Niskin rosette and deployed an average of 10 times per Lagrangian cycle. It reliably images organisms that are >600-μm in diameter, although avoidance issues should be expected for strongly swimming taxa. The UVP5 images a volume of ∼1 L per image at ∼6 Hz. Data were analyzed as described in , . Briefly, images were automatically analyzed to separate out vignettes representing organisms or marine snow aggregates. ZooProcess software was utilized to generate morphometric information (e.g., diameter) and classify the organisms into broad taxonomic groups. Classifications were then 100% manually validated. Phaeodarian taxonomic groups included Aulosphaeridae, which was the dominant rhizarian present and is the only group of organisms enumerated by UVP5 that is utilized in this study.
Sediment Trap Deployments
VERTEX-style particle interceptor tube (PIT) sediment traps were deployed at the beginning and recovered at the end of each cycle (Knauer et al., 1979; Stukel et al., 2013). PITs consisted of a polycarbonate tube with 7-cm inner diameter and 8:1 aspect ratio with a baffle on top comprised of 13 smaller, beveled tubes. On the P0704 cruise, PITs were deployed at a depth of 100 m. On P0810, P1106, and P1208, PITs were deployed at a depth of 100 m and near the base of the euphotic zone (as estimated based on fluorescence profiles from MVP transects) if the base of the euphotic zone was shallower than 75 m. On the P1408 and P1604 cruises, PITs were deployed at the base of the euphotic zone, 100 m, and 150 m.
PITs were deployed with a dense formaldehyde-filtered seawater brine. Deployments lasted from 2.25 to 4.25 days. After recovery, the interface separating brine water from overlying ambient seawater was identified and the overlying seawater was immediately removed. Samples were then gravity filtered through a 200-μm mesh filter, the filter examined under a stereomicroscope, and swimming zooplankton removed. Three to five replicates per depth were filtered through pre-combusted GF/F filters and used for particulate organic carbon analyses (¼ to ½ of a tube). An additional three samples were filtered through pre-combusted quartz (QMA) filters and used for C:234Th ratio measurements (Stukel et al., 2019). On the P0704, P0810, and P1604 cruises, two replicates (one half tube each) per depth were saved in formaldehyde and analyzed under a stereomicroscope to quantify fecal pellet abundance (Stukel et al., 2013; Morrow et al., 2018). For detailed methods and information on additional analyses made from these traps, see Gutierrez-Rodriguez et al. (2018), Morrow et al. (2018), and Stukel et al. (2019).
234Th Analyses
Water column 234Th activity was measured (typically two profiles per cycle and 10–12 depths per profile, spanning the upper 200 m of the water column) using standard small volume techniques (; Pike et al., 2005). Briefly, 4-L samples were spiked with tracer 230Th and thorium was co-precipitated with manganese oxide. Samples were beta counted on a RISO low-level background beta counter and re-counted >6 half-lives later. Samples were dissolved and spiked with 229Th. The 229:230Th ratio was determined by inductively coupled plasma mass spectrometry to determine the yield of the initial thorium filtration. For additional details, see Stukel et al. (2019). 238U–234Th deficiency was quantified after determining 238U activity from relationships with salinity published in Owens et al. (2011). 238U–234Th deficiency was combined with sediment trap organic carbon flux measurements to estimate twilight zone CFA as outlined below.
Carbon Flux Attenuation
Within the CCE, which has high mesoscale variability and pronounced horizontal currents, we consider our drifting sediment traps to provide a more accurate estimate of carbon flux than 234Th (see Supplementary Materials). Our supposition that the sediment traps have no substantial over- or under-collection bias is supported by a total of 56 paired sediment trap and 238U–234Th deficiency measurements showing good agreement (see “Results” section). Consequently, we use sediment trap values of carbon flux at deployment depths (typically near the base of the euphotic zone and at 100 m) and utilize 238U–234Th measurements to generate smooth profiles of CFA above, between, and below sediment trap deployment depths. Specifically, carbon flux at a depth horizon (D) can be quantified from 238U–234Th deficiency using a steady-state without advection equation:
where CTh(D) is the C:234Th ratio of sinking particles at the depth horizon of interest, λ234 is the 234Th decay constant, and Def(z) is the 238U–234Th deficiency at depth z, which is equal to the activity of 238U minus the activity of 234Th. CFA can thus be calculated from the first derivative of Eq. 1:
Stukel et al. (2019) found a strong relationship between the C:234Th ratio of sinking particles and the ratio of vertically integrated POC to vertically integrated total water column 234Th (vC:234Thtot). Using this equation allows us to determine CTh(D) as a smoothly varying function of D. For additional details see Supplementary Appendix S1.
Particle Sinking Speed
We estimate sinking speeds from microscopic analyses of fecal pellets collected in the sediment traps (Stukel et al., 2013; Morrow et al., 2018) and a relationship between fecal pellet size and sinking rate. A strong relationship between size and sinking speed is a consistent finding of studies that have quantified fecal pellet sinking rates (Small et al., 1979; ; Turner, 2015). Using fecal pellet sinking rate as a function of equivalent spherical diameter (ESD) data reviewed in Stukel et al. (2014), a power-law relationship suggests that sinking speed can be predicted from ESD:
where sinking speed is in units of m day-1 and ESD is in units of mm. This relationship is derived from multiple studies spanning a range of taxonomic groups including appendicularians (Ploug et al., 2008), copepods and euphausiids (Smayda, 1971; Turner, 1977; Small et al., 1979; Yoon et al., 2001; Ploug et al., 2008), chaetognaths (), and thaliaceans and pteropods (; Madin, 1982; Yoon et al., 2001). We thus expect it to be broadly representative of sinking rates of fecal pellets produced by the mixed zooplankton assemblages encountered in the southern CCE.
Flux-Feeding Calculations
The impact of flux-feeding zooplankters (such as a thecosome pteropod or phaeodarian) on particle flux can be quantified based on the effective cross-sectional area over which particles are collected (Jackson et al., 1993):
where F is the flux of sinking particles and Nx is the numerical concentration of suspension feeders (NAulo is the abundance of Aulosphaeridae and Nptero is the abundance of pteropods; all a function of depth). σx is the cross-sectional area over which the organisms are intercepting particles: σx = π/4 × ESDeff2. Following Stukel et al. (2018a), ESDeff is the effective diameter over which the organisms collect sinking particles and is dependent on both the diameter of the collection apparatus and the average cross-sectional area of sinking particles. For the CCE, the effective diameter of sinking fecal pellets was calculated as Dpar = 403 μm. We then calculate:
where ESDAulo is the equivalent spherical diameter of Aulosphaeridae cells (measured separately by UVP5 for each cell) and ESDweb is the equivalent spherical diameter of a pteropod feeding web [assumed to be 45 mm based on summary in ]. The factor of 1.25 in Eq. 6 represents the ratio of cell diameter (reported by UVP5) to diameter including radial spines.
Suspension-Feeding Calculations
The impact of suspension-feeding zooplankton (within which we include both true filter-feeders, e.g., appendicularians, and other organisms such as some herbivorous copepods that may use feeding currents that can be modeled as suspension-feeding) on particle flux attenuation can be calculated as:
where Nx is the number concentration of suspension feeders with a clearance rate of CR0, and S is the sinking speed of the particles. The impact of suspension-feeding zooplankton thus depends on the range of sinking speeds of sinking particles. We calculated clearance rates using the allometric scaling relationship determined in Kiørboe (2011):
where CR0 is the clearance rate (mL day-1), BC is the individual carbon content (g), a = 7.31 (±0.27), b = 1.01 (±0.05), and the Q10 used to account for temperature effects was assumed (following Hansen et al., 1997) to be 2.8.
We utilized zooplankton data from nighttime MOCNESS tows sorted into broad taxonomic and size groups using ZooScan. The choice to use only nighttime biomass likely leads to a conservative estimate of flux attenuation, because it excludes the activity of diel vertical migrants. This decision was made because these vertical migrants likely feed primarily in the surface layers and are not actively feeding at depth. The carbon biomass of individual organisms was determined using length:carbon relationships for different taxa as outlined in Table 1 of Stukel et al. (2013). We calculated clearance rates for four groups of potentially suspension-feeding mesozooplankton: copepods, euphausiids, appendicularians, and ostracods, while recognizing that some taxa include omnivores or predators. Nauplii and doliolids were also sorted in the samples, but their biomasses were low and hence are not considered further. Other taxa, including chaetognaths, were abundant, but are not suspension feeders.
Results
Carbon Flux and Carbon Flux Attenuation
We used two independent estimates of particle flux (sediment traps and 238U–234Th) to test whether or not our estimates of particle flux (and particle flux attenuation) suffered from methodological biases (Figure 2). Comparisons between 234Th flux collected in sediment traps and 234Th flux estimated from 238U:234Th disequilibrium [computed using a one-dimensional, steady-state model without upwelling, Savoye et al. (2006)] showed strong agreement, with occasional outliers. The median ratio of sediment trap-derived flux to deficiency-derived flux was 0.998 suggesting near perfect agreement on a typical deployment. However, the overall mean of the sediment trap dataset was 6% higher than that of the Th deficiency dataset, suggesting that when there was a substantial disagreement, the sediment trap results were likely to be higher. This is not surprising in a dynamic system, where bloom decay can lead to spikes in particle export on faster time-scales than the approximately monthly temporal integration time-scale of 238U–234Th deficiency approaches. The overall agreement between sediment trap and thorium-based approaches suggests that our sediment traps had neither an over- nor an under-collection bias. Furthermore, our results suggest that although 234Th-derived flux may not perfectly match with contemporaneous processes occurring in the surface layer, there is no reason to suspect that flux attenuation calculations based on 234Th measurements will have a systematic bias.
FIGURE 2
Across the 30 Lagrangian cycles, POC flux (measured by sediment trap) at 100 m depth ranged from 2.6 to 24.9 mmol C m-2 day-1. Export was higher in coastal areas than in the oligotrophic, offshore domain and there was a strong correlation between POC flux and both primary productivity and surface chl. When independent sediment trap export measurements were made at the base of the euphotic zone and at 100 m depth, the correlation (Spearman’s ρ) between flux at these two depths was 0.72 (p = 0.002). This across-depth correlation was substantially weaker when only data from frontal regions (P1106 and P1208 Cruises) were analyzed (ρ = 0.49, p = 0.36) than when non-front data were considered (ρ = 0.92, p = 4.5 × 10-4). For all paired samples, the median ratio of sediment trap flux at 100 m to sediment trap flux near the base of the euphotic zone (shallow trap depths varied from 47 to 70 m) was 0.72, suggesting that approximately 28% of sinking POC could be expected to be remineralized before a depth of 100 m on these cycles.
To determine continuous profiles of carbon flux and CFA, we merged sediment trap data with 234Th data. We restricted this flux attenuation analysis to the depth range from the base of the euphotic zone to a depth of 150 m, because at deeper depths we cannot constrain the C:234Th ratio with certainty. The depth range of our flux attenuation calculations thus corresponds with the upper twilight zone, where zooplankton roles in flux attenuation have been hypothesized to be particularly important (Jackson and Checkley, 2011). CFA in the upper twilight zone was highly variable (Figure 3A). When comparing across all cycles, flux attenuation decreased with depth from a mean (across all cycles) flux attenuation of 913 μmol C m-2 day-1 decrease in flux over a 10 m depth range between 80 and 90 m (interquartile range was 57–874 μmol C m-2 day-1/10 m) to a mean of 495 μmol C m-2 day-1/10 m (interquartile = 328–692 μmol C m-2 day-1/10 m) between 140 and 150 m depth. Additional patterns can be seen when considering the oligotrophic and high biomass cycles separately. For the oligotrophic cycles, the base of the euphotic zone (1% light level) was typically in the range of 60–80 m. We thus considered flux attenuation starting at the deeper limit of this range. Flux attenuation increased with depth from a mean of 231 μmol C m-2 day-1/10 m (interquartile = -82 to 418.5) in the 80–90 m depth range to a mean of 503 μmol C m-2 day-1/10 m (interquartile = 233–799) in the 100–110 m depth range. Beneath this depth flux attenuation declined gradually (Figure 3B). For the high biomass cycles a similar pattern was seen. Flux attenuation over the 50–60 m depth range in the high biomass cycles averaged 595 μmol C m-2 day-1/10 m (interquartile = -338 to 1221) and increased to 1503 μmol C m-2 day-1/10 m (interquartile = 514–1.981) in the 80–90 m depth range.
FIGURE 3
In contrast to absolute flux attenuation, the percentage of flux remineralized over a 10-m depth range did not decrease beneath 100 m (Figure 3D–F). Instead, relative flux attenuation increased somewhat from near the base of the euphotic zone to 150 m. This pattern was relatively consistent when restricting analyses to the oligotrophic cycles (Figure 3E), but was less distinct in the high biomass cycles (Figure 3F). Nevertheless, across depth and regions, relative CFA was typically in the range of 6–10% of carbon flux remineralized over a 10-m depth range. The median relative flux attenuation in the shallow twilight zone varied from 6.0 to 9.0% of flux/10 m (mean ranged from 6.5 to 13%/10 m). This sets a reasonable expectation for the amount of flux attenuation that must be mediated by the combined mesopelagic microbial, mesozooplankton, and nekton communities.
Suspension-Feeding Zooplankton and Particle Flux Attenuation
Suspension-feeding crustaceans are the most abundant (by carbon biomass) ecological category of zooplankton in the CCE (Lavaniegos and Ohman, 2007). To quantify the potential role of these suspension feeders, we utilized data from nighttime MOCNESS tows in the CCE. Suspension-feeder biomass in the mesopelagic was typically dominated by copepods and euphausiids (Supplementary Figure S2). Between 100 and 450 m depth, copepod biomass was typically in the range of 1–3 mg C m-3. Euphausiid biomass was almost always <1 mg C m-3 at depths deeper than 200 m (and typically <0.4 mg C m-3) but was greater in the shallow mesopelagic. Ostracods had substantially lower biomass, never exceeding 0.6 mg C m-3 in the mesopelagic. Appendicularians, while occasionally abundant in the euphotic zone, were also relatively minor contributors to total suspension-feeder biomass in the mesopelagic.
Community clearance rates ranged from 2 to 81 L m-3 day-1 in the mesopelagic, which corresponds to mesozooplankton clearing between 0.2 and 8.1% of the water each day (Figure 4). Clearance rates were dominated by copepods and euphausiids, with copepods dominating beneath 150 m and both having substantial contributions at shallower depths.
FIGURE 4
The impact of suspension feeders on particle flux attenuation depends on the average sinking rates of particles. We chose P0704-1 as a cycle with typical community clearance rates and calculated flux attenuation for particles with a full range of sinking speeds (Figure 5). For particles with a sinking speed of 1 m day-1, nearly all flux would be consumed by suspension feeders within 50 m of the depth of particle creation. By contrast, for particles sinking at a speed of 10 m day-1, only ∼70% of flux would be expected to be consumed before particles reach a depth of 450 m. For particles sinking at a rate of >50 m day-1, the impact of suspension feeders becomes negligible. Comparing data from all cycles, we see a similar pattern. If particles sink at 1 m day-1, anywhere from ∼5 to 50% of sinking particles would be consumed over a 10 m depth range (Figure 6A). For particles sinking with a speed of 10 m day-1, particle flux attenuation was typically <4%/10 m and for particles sinking at a speed of 100 m day-1, particle flux attenuation was always <1%/10 m and typically <0.4%/10 m.
FIGURE 5
FIGURE 6
To investigate typical particle sinking speeds, we measured the abundance and size of recognizable fecal pellets collected in sediment traps. These pellets were almost always the dominant visually identifiable component of the sinking material, although in the oligotrophic regions the identifiable pellets typically comprised less than half of total carbon flux. We applied allometric–scaling relationships between sinking speed and fecal pellet size to estimate fecal pellet settling velocities. Across all samples, half of the fecal pellet carbon flux was mediated by pellets sinking at a speed slower than 141 m day-1 (Figure 7). Only 10% of the flux was due to pellets sinking slower than 78 m day-1 and only 1% was derived from pellets sinking slower than 46 m day-1. Even for the cycle most dominated by small fecal pellets, half of the pellet flux was due to pellets sinking faster than 85 m day-1. When considering these results in light of the calculated impact of suspension-feeding zooplankton on flux attenuation of particles with average sinking speeds in the range of 100 m day-1, it becomes clear that suspension-feeding mesozooplankton are not likely to be playing a dominant role in CFA in the mesopelagic. This does not mean, however, that they play no role in flux attenuation. Rather, the abundance and activity of suspension feeders in the epipelagic may be the reason that so few slowly sinking particles were collected in sediment traps beneath the euphotic zone.
FIGURE 7
Flux-Feeding Zooplankton and Particle Flux Attenuation
To investigate the potential role of flux- feeding zooplankton in CFA, we quantified the abundances of two prominent taxa of flux feeders in the CCE: Aulosphaeridae (a phaeodarian) and L. helicina (a thecosome pteropod). These are certainly not the only flux-feeding zooplankton in the CCE, hence our calculations of the total contribution of flux feeders should be considered quite conservative.
Aulosphaeridae abundance consistently peaked in the shallow twilight zone between 50 and 150 m depth (Figure 8A–E). Abundances were substantially higher on the P0810 and P1106 cruises (peak cycle average abundances reaching >500 cells m-3) than on the warm period cruises (P1408 and P1604, peak abundances <50 cells m-3). Unlike most other plankton community and biogeochemical measurements, there was not a strong correlation between Aulosphaeridae abundance and primary productivity (cf. ). Using Eqs 4 and 5, we quantified the potential impact of Aulosphaeridae on CFA (Figure 8F–H). Despite high variability, Aulosphaeridae often had an important role in flux attenuation. Near the depth of peak Aulosphaeridae abundance (∼100 m), the mean (across all cycles) flux attenuation attributable to Aulosphaeridae was 0.69%/10 m (median = 0.21%/10 m, interquartile range = 0.04–0.81%). The maximum Aulosphaeridae flux attenuation (Cycle P0810-1, 100–110 m depth range) was 4.2%/10 m. For comparison, the total flux attenuation mediated by all abiotic and biotic factors (as assessed using sediment traps and 234Th) was typically in the range of 6–10%/10 m. Thus, the Aulosphaeridae (a single family of phaeodarians) can play a substantial role in CFA, although at most times its contribution is <5% of total flux attenuation.
FIGURE 8
Limacina helicina was only quantified in MOCNESS samples from two cycles (P1208-3, which had a 31-m deep euphotic zone, and P1208-4, which had a 70-m euphotic zone). L. helicina abundances typically declined from peak abundances of ∼0.5–2 individuals m-3 in the upper 100 m of the water column to <0.2 individuals m-3 beneath 100 m depth (Figure 9A). We computed (using Eqs 4 and 6) the potential role of L. helicina in intercepting sinking particles and found that they could potentially intercept between 4 and 10% of sinking particles between the base of the euphotic zone and a depth of 400 m (Figure 9B). However, their impact was concentrated just below the base of the euphotic zone, where on Cycle P1208-3 they consumed an average of 1.2–1.6% of carbon flux/10 m and on P1208-4 they consumed an average of 0.45–0.81%/10 m (Figure 9C,D). As for the Aulosphaeridae, the evidence suggests that this species can intercept a substantial portion of the sinking particles, but their impact is concentrated on the region immediately beneath the euphotic zone.
FIGURE 9
Discussion
Zooplankton as Gatekeepers to the Mesopelagic
Zooplankton play diverse roles in the epi- and mesopelagic (Steinberg and Landry, 2017). Our results show that suspension-feeding mesozooplankton (especially copepods and euphausiids) are abundant in the twilight zone and can have substantial clearance rates on the magnitude of 50 L m-3 day-1 (Figure 4 and Supplementary Figure S2). However, these clearance rates are not sufficient to give them a meaningful impact on CFA of the particles sinking at ∼100 m day-1, which dominate flux in the region (Figure 5, 6). Instead, their activity leads to near complete consumption of slowly sinking particles near the base of the euphotic zone. This is likely a general result, because studies focused on in situ measurement of particle sinking speeds often find typical velocities on the order of 100 m day-1 beneath the euphotic zone (
We also note that our calculations were based on nighttime net tows, thus largely excluding any potential impacts of diel vertically migrating suspension feeders. This decision was based on the assumption that diel vertical migrants feed primarily in the surface layers, not at depth. If we calculate clearance rates based on daytime net tows, we find increased mesopelagic clearance rates, particularly at depths between 150 and 250 m. However, this increased clearance rate (daytime values were typically 50% higher than nighttime values) still resulted in a low impact of suspension-feeders on rapidly sinking (100 m day-1) particles (median across all cycles was <0.4%/10 at all depths), except for Cycle P0810-3. During this cycle, abundant vertically migrating euphausiids were present and, if they were feeding at depth, the total suspension-feeding community could have been responsible for attenuation of 4% of C flux/10 m during daytime hours (in the depth range 200–250 m). Notably, Euphausia pacifica (one of the dominant euphausiids in the CCE) has been shown to play an important role in disaggregation of marine snow in the euphotic zone (
Despite their comparatively weak capacity for intercepting rapidly sinking particles, suspension-feeding zooplankton in the mesopelagic have substantial carbon demands in the CCE (Kelly et al., unpublished) and many other ecosystems (Hernández-Leon and Ikeda et al., 2005; Steinberg et al., 2008b;
In contrast, flux-feeding zooplankton can be substantial loss terms for rapidly sinking particles. The mean percentage flux attenuation for Aulosphaeridae in the shallow mesopelagic was 0.69%/10-m depth horizon corresponding to the depth of maximum Aulosphaeridae abundance (Figure 8). For comparison, in the same depth range total flux averaged 7.1%/10 m (Figure 3). This implies that a single family of giant Rhizaria (∼2-mm) may be responsible for nearly 10% of the flux attenuation in the layer immediately beneath the euphotic zone. Similarly, the pteropod L. helicina is abundant beneath the euphotic zone and may be responsible for ∼10–20% of total flux attenuation in the depth ranges where it is most common. These are certainly not the only flux feeders in the CCE. Many other pteropod and rhizarian taxa occur in the CCE (Kling and Boltovskoy, 1995;
Also of import, these flux feeders not only intercept and consume sinking particles, but also produce sinking fecal pellets. L. helicina fecal pellets have been shown to contribute 19% of the POC flux in a coastal bay near Antarctica (Manno et al., 2009), while “mini-pellets” produced by phaeodarians were abundant contributors to sinking flux in the Eastern Tropical Pacific and Northeast Atlantic (Gowing and Silver, 1985; Lampitt et al., 2009). Furthermore, due to their dense aragonite shells and siliceous tests, respectively, dead pteropods and phaeodarians likely contribute substantially to sinking flux and are often found in sediment trap material (Takahashi and Honjo, 1981;
Microbes and Zooplankton in the Mesopelagic
The abundance (and impact on flux attenuation) of suspension-feeding and flux-feeding zooplankton generally decreases with depth in the mesopelagic. In particular, the flux-feeding zooplankton that we believe play a disproportionately strong role in twilight zone flux attenuation both decreased sharply in abundance beneath a depth of ∼100 m. However, the percentage (but not absolute magnitude) of CFA increased slightly with depth between 80 and 150 m depth (Figure 3). While the substantial variability in flux attenuation profiles prevents us from making definitive conclusions about these patterns, this mismatch between flux attenuation patterns determined from bulk flux estimates (sediment traps and 234Th) and the abundance of flux feeders points to the importance of other sources of carbon remineralization. Stemmann et al. (2004a) suggested based on modeling results for the Mediterranean Sea that flux feeders may be particularly important in the shallow twilight zone, while microbial degradation becomes increasingly important at the lower particle concentrations experienced in deeper waters. Such a result seems plausible for multiple reasons. First, free-living zooplankton will have substantially greater difficulty obtaining sufficient energy for basal metabolic needs and respiration in the deeper ocean, while particle-attached microbes should have ample food at all depths. Second, the abundance of particle-attached microbes on a sinking particle may be expected to increase with age (and depth) of the particle, because it may be continually colonized by new microbes while it sinks, and these microbes are likely to multiply on the sinking particle.
Results from other studies in the CCE can highlight the potential impact of microbes. Simon et al. (1990) estimated that marine snow turnover times with respect to bacterial remineralization ranged from 20 to 100 d. This is similar to the turnover times expected to be mediated by suspension-feeding zooplankton, and suggests that bacteria have only a minor impact on the flux of rapidly sinking particles. However, Samo et al. (2012) measured bacterial carbon production rates (including particle-attached and free-living bacteria) at a depth of 100 m that ranged from ∼10 to 100 μmol C m-3 day-1. For comparison, we quantified that the total flux attenuation at this depth was 75 μmol C m-3 day-1 (median = 66 μmol C m-3 day-1; 95% confidence interval of up to 313 μmol C m-3 day-1). It is also unclear whether decomposition rates for marine snow aggregates are comparable to those of the fecal pellets that dominate sediment trap material in the CCE. Gowing and Silver (1983) suggested that fecal pellets collected near our study region had substantial contributions of interior bacteria that likely originated as enteric or ingested, but digestion-resistant bacteria. These communities would likely be substantially different, and with different biogeochemical impacts, than those found on aggregates.
Heterotrophic protists can also play dominant roles in the microbial communities consuming fecal pellets (Poulsen and Iversen, 2008). Gutierrez-Rodriguez et al. (2018) used 18S sequencing to investigate protistan sequences in sediment traps on our P1408 cruise. They found a consistent and substantial increase in the relative contribution of dinophytes in un-preserved trap samples relative to formaldehyde-preserved trap samples. Dinophytes typically increased from <10% of total protistan reads to >50%, despite deployment times that were only 3.25 days. Stramenopiles (mainly heterotrophic nanoflagellate taxa) also increased in the unpreserved samples. This suggests potentially rapid growth rates for these protists on sinking particles, with a commensurate role in feeding on either detritus contained in the sinking particles or other microbial taxa transported within the particles.
Multiple studies from other regions have investigated the differing roles of zooplankton and microbes in the mesopelagic (Simon et al., 2002; Robinson et al., 2010; Steinberg and Landry, 2017). Steinberg et al. (2008b) found that bacterial carbon demand exceeded zooplankton demand in the North Pacific subtropical gyre by a factor of 3, although the two were approximately equal in the subarctic gyre. The combined respiration of these two groups was found to be larger than that entering the mesopelagic through sinking particles, a result that has been consistently found in other regions including the subtropical and north Atlantic (
Epipelagic-Mesopelagic Coupling in the CCE
Our result that suspension-feeding zooplankton are abundant in the mesopelagic, but play only minor roles in consuming sinking particles, raises important questions about how they satisfy their metabolic demands. Based on results from our P0704 and P0810 cruises, Kelly et al. (unpublished) estimated that total mesozooplankton (mostly suspension-feeders) respiration in the mesopelagic ranged from 3.2 to 18 mg C m-2 day-1 across these cruises. These values substantially exceed the total carbon that we calculate suspension-feeders would be able to consume from rapidly sinking particles settling through the mesopelagic. However, the Kelly et al. (unpublished) results also suggest a resolution to this apparent imbalance. Diel vertically migrating taxa were abundant in the region and their active transport provides substantial energy subsidies to mesopelagic food webs. Mesopelagic resident organisms derive energy directly (through predation) and indirectly (through their release of dissolved organic carbon that supports microbial communities) from vertical migrants. Indeed, results showed that nearly half of the organic matter consumed by mesopelagic resident mesozooplankton was derived (directly or indirectly) from foodweb pathways that originated with active transport by vertical migrants rather than pathways originating from sinking particle flux (Kelly et al., unpublished).
This highlights an important reality of mesopelagic ecosystems. Although we most often think of these ecosystems as being supported by sinking particles (and to a lesser extent diel vertical migration), mesopelagic communities feature complex ecological relationships between particle-attached and free-living bacteria, protistan bacterivores, and zooplankton (Robinson et al., 2010). These zooplankton communities are diverse and have many feeding modes including suspension-feeding (likely on protists and suspended organic matter), flux-feeding, predation, and parasitism (Steinberg et al., 2008a; Wilson et al., 2010;
Conclusion
Our results offer compelling evidence that two flux-feeding zooplankton taxa (the phaeodarian Aulophaeridae and the pterpopod L. helicina) can exert substantial impact on sinking particles in the shallow twilight zone. Both taxa have the capacity to consume greater than 1% of sinking particles per 10-m vertical depth range at the depths at which they are most abundant, which equates to ∼10–20% of total CFA across those same depth ranges (typically the first 50–100 m beneath the euphotic zone). These are unlikely to be the only flux feeders in the mesopelagic. Indeed many other rhizarians and pteropods are found in the CCE and many other taxa are suspected of full or part-time flux-feeding behavior. Flux-feeders are likely particularly important in consuming rapidly sinking particles that would otherwise penetrate into the deep ocean. The abundance of suspension-feeding zooplankton in the ocean suggests that slowly sinking particles will have a relatively minor contribution to sinking carbon flux in the mesopelagic because they will be rapidly consumed. While our results suggest that feeding mode leads to very different biogeochemical importance for these functional groups of zooplankton, substantial work is needed to directly quantify the impact of these organisms in situ and to investigate the spatial and temporal distributions of each feeding trait amongst the diverse zooplankton communities found in the ocean (
Statements
Data availability statement
The datasets generated for this study are available on request to the corresponding author.
Author contributions
MS and TK were responsible for the sediment trap and thorium measurements. MO was responsible for the MOCNESS zooplankton data. MO and TB were responsible for the UVP deployments. TB was responsible for the UVP rhizarian analyses. MS wrote the manuscript. All authors edited the final version of the manuscript.
Funding
This work was funded by the NSF Bio OCE grants to the CCE LTER Program: OCE-0417616, OCE-1026607, OCE-1637632, and OCE-1614359. A portion of this work was performed at the National High Magnetic Field Laboratory, which is supported by the National Science Foundation Cooperative Agreement No. DMR-1644779 and the State of Florida.
Acknowledgments
We thank the captains and crews of the R. V. Melville, R. V. Thompson, and R. V. Sikuliaq as well as our many collaborators in the CCE LTER Program, without them this study would not have been possible. We also thank N. Bednaršek for the pteropod enumerations. Data used in this manuscript can be found on the CCE LTER DataZoo website: https://oceaninformatics.ucsd.edu/datazoo/catalogs/ccelter/datasets.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2019.00397/full#supplementary-material
References
1
AlcarazM.AlmedaR.CalbetA.SaizE.DuarteC. M.LasternasS.et al (2010). The role of arctic zooplankton in biogeochemical cycles: respiration and excretion of ammonia and phosphate during summer.Polar Biol.331719–1731. 10.1007/s00300-010-0789-9
2
AlldredgeA. L. (1976). Discarded appendicularian houses as sources of food, surface habitats, and particulate organic matter in planktonic environments.Limnol. Oceanogr.2114–23.
3
AlldredgeA. L.GotschalkC. (1988). In situ settling behavior of marine snow.Limnol. Oceanogr.33339–351. 10.1016/j.scitotenv.2016.09.115
4
ArchibaldK. M.SiegelD. A.DoneyS. C. (2019). Modeling the impact of zooplankton diel vertical migration on the carbon export flux of the biological pump.Global Biogeochem. Cycles33181–199. 10.1029/2018gb005983
5
ArmstrongR. A.PetersonM. L.LeeC.WakehamS. G. (2009). Settling velocity spectra and the ballast ratio hypothesis.Deep Sea Res. II561470–1478. 10.1016/j.dsr2.2008.11.032
6
ArnostiC. (2011). “Microbial extracellular enzymes and the marine carbon cycle,” in Annual Review of Marine ScienceVol. 3edsCarlsonC. A.GiovannoniS. J. (Palo Alto: Annual Reviews), 401–425. 10.1146/annurev-marine-120709-142731
7
AsherE.DaceyJ.StukelM.LongM.TortellP. (2016). Processes driving seasonal variability in DMS, DMSP, and DMSO concentrations and turnover in coastal Antarctic waters.Limnol. Oceanogr.62104–124. 10.1002/lno.10379
8
BainesS. B.ChenX.TwiningB. S.FisherN. S.LandryM. R. (2016). Factors affecting Fe and Zn contents of mesozooplankton from the Costa Rica Dome.J. Plankton Res.38331–347. 10.1093/plankt/fbv098
9
BaltarF.ArísteguiJ.GasolJ. M.SintesE.HerndlG. J. (2009). Evidence of prokaryotic metabolism on suspended particulate organic matter in the dark waters of the subtropical North Atlantic.Limnol. Oceanogr.54182–193. 10.4319/lo.2009.54.1.0182
10
BartonA. D.PershingA. J.LitchmanE.RecordN. R.EdwardsK. F.FinkelZ. V.et al (2013). The biogeography of marine plankton traits.Ecol. Lett.16522–534. 10.1111/ele.12063
11
BathmannU. V.NojiT. T.Von BodungenB. (1991). Sedimentation of pteropods in the Norwegian Sea in autumn.Deep Sea Res. Part A Oceanogr. Res. Pap.381341–1360. 10.1016/0198-0149(91)90031-a
12
BednaršekN.MožinaJ.VogtM.O’brienC.TarlingG. (2012). The global distribution of pteropods and their contribution to carbonate and carbon biomass in the modern ocean.Earth Syst. Sci. Data4167–186. 10.5194/essd-4-167-2012
13
BednaršekN.OhmanM. D. (2015). Changes in pteropod distributions and shell dissolution across a frontal system in the California current system.Mar. Ecol. Progr. Ser.52393–103. 10.3354/meps11199
14
Benitez-NelsonC. R.BuesselerK. O.Van Der LoeffM. R.AndrewsJ.BallL.CrossinG.et al (2001). Testing a new small-volume technique for determining Th-234 in seawater.J. Radioanal. Nuclear Chem.248795–799.
15
BianchiD.GalbraithE. D.CarozzaD. A.MislanK. A. S.StockC. A. (2013). Intensification of open-ocean oxygen depletion by vertically migrating animals.Nat. Geosci.6545–548. 10.1038/ngeo1837
16
BiardT.KrauseJ. W.StukelM. R.OhmanM. D. (2018). The significance of giant phaeodarians (Rhizaria) to biogenic silica export in the california current ecosystem.Global Biogeochem. Cycles32987–1004. 10.1029/2018gb005877
17
BiardT.OhmanM. D. (2019). Vertical niche definition of test-bearing protists (Rhizaria) into the twilight zone revealed by in situ imaging. bioRxiv
18
BiardT.StemmannL.PicheralM.MayotN.VandrommeP.HaussH.et al (2016). In situ imaging reveals the biomass of giant protists in the global ocean.Nature532504–507. 10.1038/nature17652
19
BodeM.HagenW.SchukatA.TeuberL.Fonseca-BatistaD.DehairsF.et al (2015). Feeding strategies of tropical and subtropical calanoid copepods throughout the eastern Atlantic Ocean - Latitudinal and bathymetric aspects.Progr. Oceanogr.138268–282. 10.1016/j.pocean.2015.10.002
20
BoydP. W.SherryN. D.BergesJ. A.BishopJ. K. B.CalvertS. E.CharetteM. A.et al (1999). Transformations of biogenic particulates from the pelagic to the deep ocean realm.Deep Sea Res. Part II Top. Stud. Oceanogr.462761–2792. 10.1016/s0967-0645(99)00083-1
21
BrulandK. W.SilverM. W. (1981). Sinking rates of fecal pellets from gelatinous zooplankton (salps, pteropods, doliolids).Mar. Biol.63295–300. 10.1007/bf00395999
22
BuitenhuisE.Le QuereC.AumontO.BeaugrandG.BunkerA.HirstA.et al (2006). Biogeochemical fluxes through mesozooplankton.Global Biogeochem. Cycles20:18.
23
BurdA. B.HansellD. A.SteinbergD. K.AndersonT. R.AristeguiJ.BaltarF.et al (2010). Assessing the apparent imbalance between geochemical and biochemical indicators of meso- and bathypelagic biological activity: what the @$#! is wrong with present calculations of carbon budgets?Deep Sea Res. II571557–1571. 10.1016/j.dsr2.2010.02.022
24
BurdA. B.JacksonG. A. (2009). Particle Aggregation.Annu. Rev. Mar. Sci.165–90.
25
ChristiansenS.HovingH.-J.SchütteF.HaussH.KarstensenJ.KörtzingerA.et al (2018). Particulate matter flux interception in oceanic mesoscale eddies by the polychaete Poeobius sp.Limnol. Oceanogr.632093–2109. 10.1002/lno.10926
26
CrippsG.LindequeP.FlynnK. J. (2014). Have we been underestimating the effects of ocean acidification in zooplankton?Glob. Chang. Biol.203377–3385. 10.1111/gcb.12582
27
DaggM. (1993). Sinking particles as a possible source of nutrition for the large calanoid copepod Neocalanus cristatus in the subarctic Pacific Ocean.Deep Sea Res. Part I401431–1445. 10.1016/0967-0637(93)90121-i
28
DeVriesT.PrimeauF.DeutschC. (2012). The sequestration efficiency of the biological pump.Geophys. Res. Lett.39:L13601. 10.1029/2012GL051963
29
DeVriesT.WeberT. (2017). The export and fate of organic matter in the ocean: new constraints from combining satellite and oceanographic tracer observations.Global Biogeochem. Cycles31535–555. 10.1002/2016gb005551
30
DillingL.AlldredgeA. L. (2000). Fragmentation of marine snow by swimming macrozooplankton: a new process impacting carbon cycling in the sea.Deep Sea Res. I471227–1245. 10.1016/s0967-0637(99)00105-3
31
DucklowH. W.SteinbergD. K.BuesselerK. O. (2001). Upper ocean carbon export and the biological pump.Oceanography1450–58. 10.5670/oceanog.2001.06
32
FabryV.DeuserW. (1992). Seasonal changes in the isotopic compositions and sinking fluxes of euthecosomatous pteropod shells in the Sargasso Sea.Paleoceanography7195–213. 10.1029/91pa03138
33
FowlerS. W. (1977). Trace elements in zooplankton particulate products.Nature269:51. 10.1038/269051a0
34
FrangoulisC.ChristouE. D.HecqJ. H. (2005). “Comparison of marine copepod outfluxes: nature, rate, fate and role in the carbon and nitrogen cycles,” in Advances in Marine BiologyVol. 47ed.SmithD. (London: Academic Press Ltd.), 253–309. 10.1016/s0065-2881(04)47004-7
35
FuchsH. L.FranksP. J. S. (2010). Plankton community properties determined by nutrients and size-selective feeding.Mar. Ecol. Progr. Ser.4131–15. 10.3354/meps08716
36
GieseckeR.GonzalezH. E.BathmannU. (2010). The role of the chaetognath Sagitta gazellae in the vertical carbon flux of the Southern Ocean.Polar Biol.33293–304. 10.1007/s00300-009-0704-4
37
GilmerR.HarbisonG. (1986). Morphology and field behavior of pteropod molluscs: feeding methods in the families Cavoliniidae, Limacinidae and Peraclididae (Gastropoda: Thecosomata).Mar. Biol.9147–57. 10.1007/bf00397570
38
GlibertP. M. (1998). Interactions of top-down and bottom-up control in planktonic nitrogen cycling.Hydrobiologia3631–12. 10.1073/pnas.0809671106
39
GoerickeR. (2002). Top-down control of phytoplankton biomass and community structure in the monsoonal Arabian Sea.Limnol. Oceanogr.471307–1323. 10.4319/lo.2002.47.5.1307
40
GoldthwaitS.YenJ.BrownJ.AlldredgeA. (2004). Quantification of marine snow fragmentation by swimming euphausiids.Limnol. Oceanogr.49940–952. 10.4319/lo.2004.49.4.0940
41
GorokhovaE.SoerensenA. L.MotwaniN. H. (2018). Mercury-methylating bacteria are associated with zooplankton: a proof-of-principle survey in the Baltic Sea.bioRxiv
42
GorskyG.OhmanM. D.PicheralM.GaspariniS.StemmannL.RomagnanJ. B.et al (2010). Digital zooplankton image analysis using the ZooScan integrated system.J. Plankton Res.32285–303. 10.1093/plankt/fbp124
43
GowingM. (1989). Abundance and feeding ecology of Antarctic phaeodarian radiolarians.Mar. Biol.103107–118. 10.1007/bf00391069
44
GowingM. M. (1986). Trophic biology of phaeodarian radiolarians and flux of living radiolarians in the upper 2000 m of the North Pacific central gyre.Deep Sea Res.33655–674. 10.1016/0198-0149(86)90059-2
45
GowingM. M.BenthamW. N. (1994). Feeding ecology of phaeodarian radiolarians at the VERTEX North Pacific time series site.J. Plankton Res.16707–719. 10.1093/plankt/16.6.707
46
GowingM. M.SilverM. W. (1983). Origins and microenvironments of bacteria mediating fecal pellet decomposition in the sea.Mar. Biol.737–16. 10.1007/bf00396280
47
GowingM. M.SilverM. W. (1985). Minipellets - a new and abundant size class of marine fecal pellets.J. Mar. Res.43395–418. 10.1357/002224085788438676
48
Gutierrez-RodriguezA.StukelM. R.Lopes Dos SantosA.BiardT.ScharekR.VaulotD.et al (2018). High contribution of Rhizaria (Radiolaria) to vertical export in the California Current Ecosystem revealed by DNA metabarcoding.ISME J.13964–976. 10.1038/s41396-018-0322-7
49
HamnerW.MadinL.AlldredgeA.GilmerR.HamnerP. (1975). Underwater observations of gelatinous zooplankton: Sampling problems, feeding biology, and behavior1.Limnol. Oceanogr.20907–917. 10.4319/lo.1975.20.6.0907
50
HannidesC. C. S.LandryM. R.Benitez-NelsonC. R.StylesR. M.MontoyaJ. P.KarlD. M. (2009). Export stoichiometry and migrant-mediated flux of phosphorus in the North Pacific Subtropical Gyre.Deep Sea Res. I5673–88. 10.1016/j.dsr.2008.08.003
51
HansenJ. L. S.KiorboeT.AlldredgeA. L. (1996). Marine snow derived from abandoned larvacean houses: sinking rates, particle content and mechanisms of aggregate formation.Mar. Ecol. Progr. Ser.141205–215. 10.3354/meps141205
52
HansenP. J.BjornsenP. K.HansenB. W. (1997). Zooplankton grazing and growth: scaling within the 2-2,000-mu m body size range.Limnol. Oceanogr.42687–704. 10.4319/lo.1997.42.4.0687
53
HaussH.ChristiansenS.SchütteF.KikoR.Edvam LimaM.RodriguesE.et al (2016). Dead zone or oasis in the open ocean? Zooplankton distribution and migration in low-oxygen modewater eddies.Biogeosciences131977–1989. 10.5194/bg-13-1977-2016
54
HensonS. A.SandersR.MadsenE.MorrisP. J.Le MoigneF.QuartlyG. D. (2011). A reduced estimate of the strength of the ocean’s biological carbon pump.Geophys. Res. Lett.38:L04606.
55
Hernandez-LeonS.IkedaT. (2005). A global assessment of mesozooplankton respiration in the ocean.J. Plankton Res.27153–158. 10.1093/plankt/fbh166
56
HuntB.StrugnellJ.BednarsekN.LinseK.NelsonR. J.PakhomovE.et al (2010). Poles apart: the “bipolar” pteropod species Limacina helicina is genetically distinct between the Arctic and Antarctic oceans.PLoS One5:e9835. 10.1371/journal.pone.0009835
57
IversenM. H.PoulsenL. K. (2007). Coprorhexy, coprophagy, and coprochaly in the copepods Calanus helgolandicus, Pseudocalanus elongatus, and Oithona similis.Mar. Ecol. Progr. Ser.35079–89. 10.1021/acs.est.5b05905
58
JacksonG. A.CheckleyD. M.Jr.DaggM. (2015). Settling of particles in the upper 100 m of the ocean detected with autonomous profiling floats off California.Deep Sea Res. I9975–86. 10.1016/j.dsr.2015.02.001
59
JacksonG. A.CheckleyD. M. (2011). Particle size distributions in the upper 100 m water column and their implications for animal feeding in the plankton.Deep Sea Res. I58283–297. 10.1016/j.dsr.2010.12.008
60
JacksonG. A.NajjarR. G.ToggweilerJ. R. (1993). Flux feeding as a mechanism for zooplankton grazing and its implications for vertical particulate flux.Limnol. Oceanogr.381328–1331. 10.4319/lo.1993.38.6.1328
61
KikoR.HaussH.BuchholzF.MelznerF. (2016). Ammonium excretion and oxygen respiration of tropical copepods and euphausiids exposed to oxygen minimum zone conditions.Biogeosciences132241–2255. 10.5194/bg-13-2241-2016
62
KiørboeT. (2011). How zooplankton feed: mechanisms, traits and trade-offs.Biol. Rev.86311–339. 10.1111/j.1469-185x.2010.00148.x
63
KiørboeT.ThygesenU. H. (2001). Fluid motion and solute distribution around sinking aggregates. II. Implications for remote detection by colonizing zooplankters.Mar. Ecol. Progr. Ser.21115–25. 10.3354/meps211015
64
KlingS. A.BoltovskoyD. (1995). Radiolarian vertical distribution patterns across the southern California current.Deep Sea Res. I42191–231. 10.1016/0967-0637(94)00038-t
65
KnauerG. A.MartinJ. H.BrulandK. W. (1979). Fluxes of particulate carbon, nitrogen, and phosphorus in the upper water column of the Northeast Pacific.Deep Sea Res.2697–108. 10.1016/0198-0149(79)90089-x
66
KomarP. D.MorseA. P.SmallL. F.FowlerS. W. (1981). An analysis of sinking rates of natural copepod and euphausiid fecal pellets.Limnol. Oceanogr.26172–180. 10.4319/lo.1981.26.1.0172
67
KosobokovaK.HircheH.-J.ScherzingerT. (2002). Feeding ecology of Spinocalanus antarcticus, a mesopelagic copepod with a looped gut.Mar. Biol.141503–511. 10.1007/s00227-002-0848-z
68
LalliC. M.GilmerR. W. (1989). Pelagic Snails: The Biology of Holoplanktonic Gastropod Mollusks.Palo Alto, CA: Stanford University Press.
69
LampittR.SalterI.JohnsD. (2009). Radiolaria: major exporters of organic carbon to the deep ocean.Global Biogeochem. Cycles23:GB1010.
70
LampittR.WishnerK.TurleyC.AngelM. (1993). Marine snow studies in the Northeast Atlantic Ocean: distribution, composition and role as a food source for migrating plankton.Mar. Biol.116689–702. 10.1007/bf00355486
71
LandryM. R.OhmanM. D.GoerickeR.StukelM. R.BarbeauK. A.BundyR.et al (2012). Pelagic community responses to a deep-water front in the California current ecosystem: overview of the A-front study.J. Plankton Res.34739–748. 10.1093/plankt/fbs025
72
LandryM. R.OhmanM. D.GoerickeR.StukelM. R.TsyrklevichK. (2009). Lagrangian studies of phytoplankton growth and grazing relationships in a coastal upwelling ecosystem off Southern California.Prog. Oceanogr.83208–216. 10.1016/j.pocean.2009.07.026
73
LavaniegosB. E.OhmanM. D. (2007). Coherence of long-term variations of zooplankton in two sectors of the California current system.Prog. Oceanogr.7542–69. 10.1016/j.pocean.2007.07.002
74
LawsE. A.D’saE.NaikP. (2011). Simple equations to estimate ratios of new or export production to total production from satellite-derived estimates of sea surface temperature and primary production.Limnol. Oceanogr. Methods9593–601. 10.4319/lom.2011.9.593
75
LomasM.SteinbergD. K.DickeyT.CarlsonC.NelsonN.CondonR. H.et al (2010). Increased ocean carbon export in the sargasso sea linked to climate variability is countered by its enhanced mesopelagic attenuation.Biogeosciences757–70. 10.5194/bg-7-57-2010
76
MadinL. P. (1982). Production, composition, and sedimentation of salp fecal pellets in oceanic waters.Mar. Biol.6739–45. 10.1007/bf00397092
77
MannoC.TirelliV.AccorneroA.Fonda UmaniS. (2009). Importance of the contribution of Limacina helicina faecal pellets to the carbon pump in Terra Nova Bay (Antarctica).J. Plankton Res.32145–152. 10.1093/plankt/fbp108
78
MayorD. J.SandersR.GieringS. L.AndersonT. R. (2014). Microbial gardening in the ocean’s twilight zone: detritivorous metazoans benefit from fragmenting, rather than ingesting, sinking detritus: fragmentation of refractory detritus by zooplankton beneath the euphotic zone stimulates the harvestable production of labile and nutritious microbial biomass.Bioessays361132–1137. 10.1002/bies.201400100
79
McDonnellA. M. P.BuesselerK. O. (2010). Variability in the average sinking velocity of marine particles.Limnol. Oceanogr.552085–2096. 10.4319/lo.2010.55.5.2085
80
MichaelsA. F.CaronD. A.SwanbergN. R.HowseF. A.MichaelsC. M. (1995). Planktonic sarcodines (Acantharia, Radiolaria, Foraminifera) in surface waters near Bermuda: abundance, biomass and vertical flux.J. Plankton Res.17131–163. 10.1093/plankt/17.1.131
81
MorrowR. M.OhmanM. D.GoerickeR.KellyT. B.StephensB. M.StukelM. R. (2018). Primary productivity, mesozooplankton grazing, and the biological pump in the california current ecosystem: variability and response to El niño.Deep Sea Res. I14052–62. 10.1016/j.dsr.2018.07.012
82
NakamuraY.SuzukiN. (2015). “Phaeodaria: diverse marine cercozoans of world-wide distribution,” in Marine Protists, edsOhtsukaS.SuzakiT.HoriguchiT.SuzukiN.NotF. (Berlin: Springer), 223–249. 10.1007/978-4-431-55130-0_9
83
OhmanM. D.PowellJ. R.PicheralM.JensenD. W. (2012). Mesozooplankton and particulate matter responses to a deep-water frontal system in the southern California current system.J. Plankton Res.34815–827. 10.1093/plankt/fbs028
84
OhmanM. D.RomagnanJ. B. (2016). Nonlinear effects of body size and optical attenuation on diel vertical migration by zooplankton.Limnol. Oceanogr.61765–770. 10.1002/lno.10251
85
OwensS. A.BuesselerK. O.SimsK. W. W. (2011). Re-evaluating the 238U-salinity relationship in seawater: implications for the 238U-234Th disequilibrium method.Mar. Chem.12731–39. 10.1016/j.marchem.2011.07.005
86
PicheralM.GuidiL.StemmannL.KarlD. M.IddaoudG.GorskyG. (2010). The underwater vision profiler 5: an advanced instrument for high spatial resolution studies of particle size spectra and zooplankton.Limnol. Oceanogr. Methods8462–473. 10.4319/lom.2010.8.462
87
PikeS. M.BuesselerK. O.AndrewsJ.SavoyeN. (2005). Quantification of 234Th recovery in small volume sea water samples by inductively coupled plasma-mass spectrometry.J. Radioanal. Nucl. Chem.263355–360. 10.1007/s10967-005-0594-z
88
PlougH.IversenM. H.FischerG. (2008). Ballast, sinking velocity, and apparent diffusivity within marine snow and zooplankton fecal pellets: implications for substrate turnover by attached bacteria.Limnol. Oceanogr.531878–1886. 10.4319/lo.2008.53.5.1878
89
PoulsenL. K.IversenM. H. (2008). Degradation of copepod fecal pellets: key role of protozooplankton.Mar. Ecol. Prog. Ser.3671–13. 10.3354/meps07611
90
PoulsenL. K.KiorboeT. (2005). Coprophagy and coprorhexy in the copepods Acartia tonsa and Temora longicornis: clearance rates and feeding behaviour.Mar. Ecol. Prog. Ser.299217–227. 10.3354/meps299217
91
PowellJ. R.OhmanM. D. (2015). Changes in zooplankton habitat, behavior, and acoustic scattering characteristics across glider-resolved fronts in the Southern California current system.Prog. Oceanogr.13477–92. 10.1016/j.pocean.2014.12.011
92
ProudR.CoxM. J.BrierleyA. S. (2017). Biogeography of the global ocean’s mesopelagic zone.Curr. Biol.27113–119. 10.1016/j.cub.2016.11.003
93
ReinthalerT.Van AkenH.VethC.ArísteguiJ.RobinsonC.WilliamsP. J. L. B.et al (2006). Prokaryotic respiration and production in the meso-and bathypelagic realm of the eastern and western North Atlantic basin.Limnol. Oceanogr.511262–1273. 10.4319/lo.2006.51.3.1262
94
RichardsonA. J. (2008). In hot water: zooplankton and climate change.Ices J. Mar. Sci.65279–295. 10.1093/icesjms/fsn028
95
RobinsonC.SteinbergD. K.AndersonT. R.ArísteguiJ.CarlsonC. A.FrostJ. R.et al (2010). Mesopelagic zone ecology and biogeochemistry – a synthesis.Deep Sea Res. II571504–1518. 10.1016/j.dsr2.2010.02.018
96
RobisonB. H.ReisenbichlerK. R.SherlockR. E. (2005). Giant larvacean houses: rapid carbon transport to the deep sea floor.Science3081609–1611. 10.1126/science.1109104
97
SabaG. K.SteinbergD. K.BronkD. A. (2011). The relative importance of sloppy feeding, excretion, and fecal pellet leaching in the release of dissolved carbon and nitrogen by Acartia tonsa copepods.J. Exp. Mar. Biol. Ecol.40447–56. 10.1016/j.jembe.2011.04.013
98
SamoT. J.PedlerB. E.BallG. I.PasulkaA. L.TaylorA. G.AluwihareL. I.et al (2012). Microbial distribution and activity across a water mass frontal zone in the California current ecosystem.J. Plankton Res.34802–814. 10.1093/plankt/fbs048
99
SavoyeN.Benitez-NelsonC.BurdA. B.CochranJ. K.CharetteM.BuesselerK. O.et al (2006). 234Th sorption and export models in the water column: a review.Mar. Chem.100234–249. 10.1016/j.marchem.2005.10.014
100
SchmidtK.SchlosserC.AtkinsonA.FieldingS.VenablesH. J.WaludaC. M.et al (2016). Zooplankton gut passage mobilizes lithogenic iron for ocean productivity.Curr. Biol.262667–2673. 10.1016/j.cub.2016.07.058
101
SchukatA.BodeM.AuelH.CarballoR.MartinB.KoppelmannR.et al (2013). Pelagic decapods in the northern benguela upwelling system: distribution, ecophysiology and contribution to active carbon flux.Deep Sea Res. Part I Oceanogr. Res. Pap.75146–156. 10.1016/j.dsr.2013.02.003
102
SiegelD. A.BuesselerK. O.BehrenfeldM. J.Benitez-NelsonC. R.BossE.BrzezinskiM. A.et al (2016). Prediction of the export and fate of global ocean net primary production: the EXPORTS science plan.Front. Mar. Sci.3:22. 10.3389/fmars.2016.00022
103
SiegelD. A.BuesselerK. O.DoneyS. C.SailleyS. F.BehrenfeldM. J.BoydP. W. (2014). Global assessment of ocean carbon export by combining satellite obervations and food-web models.Global Biogeochem. Cycles28181–196. 10.1002/2013gb004743
104
SilverM. W.GowingM. M. (1991). The “particle” flux: origins and biological components.Prog. Oceanogr.2675–113. 10.1016/0079-6611(91)90007-9
105
SimonM.AlldredgeA. L.AzamF. (1990). Bacterial carbon dynamics on marine snow.Mar. Ecol. Prog. Ser.65205–211. 10.3354/meps065205
106
SimonM.GrossartH. P.SchweitzerB.PlougH. (2002). Microbial ecology of organic aggregates in aquatic ecosystems.Aquat. Microb. Ecol.28175–211. 10.3354/ame028175
107
SmallL. F.FowlerS. W.UnluM. Y. (1979). Sinking rates of natural copepod fecal pellets.Mar. Biol.51233–241. 10.1021/acs.est.5b05905
108
SmaydaT. J. (1971). Normal and accelerated sinking of phytoplankton in sea.Mar. Geol.11105–122. 10.1016/0025-3227(71)90070-3
109
SteinbergD. K.CarlsonC. A.BatesN. R.GoldthwaitS. A.MadinL. P.MichaelsA. F. (2000). Zooplankton vertical migration and the active transport of dissolved organic and inorganic carbon in the Sargasso Sea.Deep Sea Res. I47137–158. 10.1016/s0967-0637(99)00052-7
110
SteinbergD. K.CopeJ. S.WilsonS. E.KobariT. (2008a). A comparison of mesopelagic mesozooplankton community structure in the subtropical and subarctic North Pacific Ocean.Deep Sea Res. II551615–1635. 10.1016/j.dsr2.2008.04.025
111
SteinbergD. K.Van MooyB. A. S.BuesselerK. O.BoydP. W.KobariT.KarlD. M. (2008b). Bacterial vs. zooplankton control of sinking particle flux in the ocean’s twilight zone.Limnol. Oceanogr.531327–1338. 10.1002/bies.201400100
112
SteinbergD. K.LandryM. R. (2017). Zooplankton and the ocean carbon cycle.Annu. Rev. Mar. Sci.9413–444. 10.1146/annurev-marine-010814-015924
113
StemmannL.JacksonG. A.GorskyG. (2004a). A vertical model of particle size distributions and fluxes in the midwater column that includes biological and physical processes-part II: application to a three year survey in the NW mediterranean sea.Deep Sea Res. I51885–908. 10.1016/j.dsr.2004.03.002
114
StemmannL.JacksonG. A.IansonD. (2004b). A vertical model of particle size distributions and fluxes in the midwater column that includes biological and physical processes-part I: model formulation.Deep Sea Res. I51865–884. 10.1016/j.dsr.2004.03.001
115
StukelM. R.KahruM.Benitez-NelsonC. R.DecimaM.GoerickeR.LandryM. R.et al (2015). Using lagrangian-based process studies to test satellite algorithms of vertical carbon flux in the eastern North Pacific Ocean.J. Geophys. Res. Oceans1207208–7222. 10.1002/2015jc011264
116
StukelM. R.BiardT.KrauseJ.OhmanM. D. (2018a). Large phaeodaria in the twilight zone: their role in the carbon cycle.Limnol. Oceanogr.632579–2594. 10.1002/lno.10961
117
StukelM. R.DécimaM.LandryM. R.SelphK. E. (2018b). Nitrogen and isotope flows through the costa rica dome upwelling ecosystem: the crucial mesozooplankton role in export flux.Global Biogeochem. Cycles321815–1832. 10.1029/2018gb005968
118
StukelM. R.KellyT. B.AluwihareL. I.BarbeauK. A.GoerickeR.KrauseJ. W.et al (2019). The carbon:234thorium ratios of sinking particles in the California current ecosystem 1: relationships with plankton ecosystem dynamics.Mar. Chem.2121–15. 10.1016/j.marchem.2019.01.003
119
StukelM. R.MislanK. A. S.DécimaM.HmeloL. (2014). “Detritus in the pelagic ocean,” in Eco-DAS IX Symposium Proceedings, ed.KempP. F. (Waco, TX: Association for the Sciences of Limnology and Oceanography), 49–76.
120
StukelM. R.OhmanM. D.Benitez-NelsonC. R.LandryM. R. (2013). Contributions of mesozooplankton to vertical carbon export in a coastal upwelling system.Mar. Ecol. Prog. Ser.49147–65. 10.3354/meps10453
121
TakahashiK.HonjoS. (1981). vertical flux of radiolaria: a taxon-quantitative sediment trap study from the Western Tropical Atlantic.Micropaleontology27140–190.
122
TrullT. W.BrayS. G.BuesselerK. O.LamborgC. H.ManganiniS.MoyC.et al (2008). In situ measurement of mesopelagic particle sinking rates and the control of carbon transfer to the ocean interior during the vertical flux in the global ocean (VERTIGO) voyages in the North Pacific.Deep Sea Res. II551684–1695. 10.1016/j.dsr2.2008.04.021
123
TurnerJ. T. (1977). Sinking rates of fecal pellets from marine copepod Pontella meadii.Mar. Biol.40249–259. 10.1007/bf00390880
124
TurnerJ. T. (2002). Zooplankton fecal pellets, marine snow and sinking phytoplankton blooms.Aquat. Microb. Ecol.2757–102. 10.3354/ame027057
125
TurnerJ. T. (2015). Zooplankton fecal pellets, marine snow, phytodetritus and the ocean’s biological pump.Prog. Oceanogr.130205–248. 10.1016/j.pocean.2014.08.005
126
UttalL.BuckK. (1996). Dietary study of the midwater polychaete poeobius meseres in monterey bay, California.Mar. Biol.125333–343. 10.1007/bf00346314
127
WakehamS. G.LeeC.HedgesJ. I.HernesP. J.PetersonM. J. (1997). Molecular indicators of diagenetic status in marine organic matter.Geochim. Cosmochim. Acta615363–5369. 10.1016/s0016-7037(97)00312-8
128
WilsonS. E.SteinbergD. K.BuesselerK. O. (2008). Changes in fecal pellet characteristics with depth as indicators of zooplankton repackaging of particles in the mesopelagic zone of the subtropical and subarctic North Pacific Ocean.Deep Sea Res. II551636–1647. 10.1016/j.dsr2.2008.04.019
129
WilsonS. E.SteinbergD. K.ChuF. L. E.BishopJ. K. B. (2010). Feeding ecology of mesopelagic zooplankton of the subtropical and subarctic North Pacific Ocean determined with fatty acid biomarkers.Deep Sea Res. I571278–1294. 10.1016/j.dsr.2010.07.005
130
WishnerK. F.OutramD. M.SeibelB. A.DalyK. L.WilliamsR. L. (2013). Zooplankton in the eastern tropical north pacific: boundary effects of oxygen minimum zone expansion.Deep Sea Res. Part I Oceanogr. Res. Pap.79122–140. 10.1016/j.dsr.2013.05.012
131
YoonW. D.KimS. K.HanK. N. (2001). Morphology and sinking velocities of fecal pellets of copepod, molluscan, euphausiid, and salp taxa in the northeastern tropical Atlantic.Mar. Biol.139923–928. 10.1007/s002270100630
Summary
Keywords
biological pump, carbon export, remineralization length scale, mesozooplankton ecology, pteropods, marine biogeochemistry, sinking particles, marine snow
Citation
Stukel MR, Ohman MD, Kelly TB and Biard T (2019) 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. doi: 10.3389/fmars.2019.00397
Received
26 March 2019
Accepted
26 June 2019
Published
12 July 2019
Volume
6 - 2019
Edited by
Rainer Kiko, GEOMAR Helmholtz Centre for Ocean Research Kiel, Germany
Reviewed by
John Patrick Dunne, Geophysical Fluid Dynamics Laboratory (GFDL), United States; Hiroaki Saito, The University of Tokyo, Japan; Santiago Hernández-León, University of Las Palmas de Gran Canaria, Spain
Updates

Check for updates
Copyright
© 2019 Stukel, Ohman, Kelly and Biard.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael R. Stukel, mstukel@fsu.edu
This article was submitted to Marine Biogeochemistry, a section of the journal Frontiers in Marine Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.