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 U-238-Th-234 disequilibrium. We found that CFA in the shallow twilight zone typically ranged from 500 to 1000 m 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 (similar to 100 m), which would equate to similar 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.

Abstract: Estimating interannual variability in carbon export is a key goal of many marine biogeochemical studies. However, due to variations in export mechanisms between regions, generalized models used to estimate global patterns in export often fail when used for intra-regional analysis. We present here a region-specific model of export production for the California Current Ecosystem (CCE) parameterized using intensive Lagrangian process studies conducted during El Niño-Southern Oscillation (ENSO) warm and neutral phases by the CCE Long-Term Ecological Research (LTER) program. We find that, contrary to expectations from prominent global algorithms, export efficiency (e-ratio = export / primary productivity) is positively correlated with temperature and negatively correlated with net primary productivity (NPP). We attribute these results to the substantial horizontal advection found within the region, and verify this assumption by using a Lagrangian particle tracking model to estimate water mass age. We further suggest that sinking particles in the CCE are comprised of a recently-produced, rapidly-sinking component (likely mesozooplankton fecal pellets) and a longer-lived, slowly-sinking component that is likely advected long distances prior to export. We determine a new algorithm for estimating particle export in the CCE from NPP (Export = 0.08 · NPP + 72 mg C m-2 d-1). We apply this algorithm to a two-decade long time series of NPP in the CCE to estimate spatial and interannual variability across multiple ENSO phases. Reduced export during the warm anomaly of 2014-2015 and El Niño 2015-2016 resulted primarily from decreased export in the coastal upwelling region of the CCE; the oligotrophic offshore region exhibited comparatively low seasonal and interannual variability in flux. The model resolves intra-regional patterns of in situ export measurements, and provides a valuable contrast to global export models.