Fixed-Nitrogen Loss Associated with Sinking Zooplankton Carcasses in a Coastal Oxygen Minimum Zone (Golfo Dulce, Costa Rica)

Oxygen minimum zones (OMZs) in the ocean are of key importance for pelagic fixed-nitrogen loss (N-loss) through microbial denitrification and anaerobic ammonium oxidation (anammox). Recent studies document that zooplankton is surprisingly abundant in and around OMZs and that the microbial community associated with carcasses of a large copepod species mediates denitrification. Here, we investigate the complex N-cycling associated with sinking zooplankton carcasses exposed to the steep O2 gradient in a coastal OMZ (Golfo Dulce, Costa Rica). 15N-stable-isotope enrichment experiments revealed that the carcasses of abundant copepods and ostracods provide anoxic microbial hotspots in the pelagic zone by hosting intense anaerobic N-cycle activities even in the presence of ambient O2. Carcass-associated anaerobic N-cycling was clearly dominated by dissimilatory nitrate reduction to ammonium (DNRA) at up to 30.8 nmol NH4+ individual-1 d-1, followed by denitrification (up to 10.8 nmol N2-N individual-1 d-1), anammox (up to 1.6 nmol N2-N individual-1 d-1), and N2O production (up to 1.2 nmol N2O-N individual-1 d-1). In contrast, anaerobic N-cycling mediated by free-living bacteria proceeded mainly through anammox and denitrification in the anoxic bottom water, which underpins the distinctive microbial metabolism associated with zooplankton carcasses. Pelagic N-loss is potentially enhanced by zooplankton carcasses both directly through N2 and N2O production, and indirectly through NH4+ production that may fuel free-living anammox bacteria. We estimate that in the hypoxic water layer of Golfo Dulce, carcass-associated N2 and N2O production enhance N-loss as much as 1.4-fold at a relative carcass abundance of 36%. In the anoxic bottom water, however, N-loss is likely enhanced only marginally due to high ambient rates and low zooplankton abundance. Thus, zooplankton carcasses may enhance N-loss mainly at the hypoxic boundaries of OMZs which are usually more extensive in open-ocean than in coastal settings. Notably, these contributions by zooplankton carcasses to pelagic N-loss remain undetected by conventional, incubation-based rate measurements.


INTRODUCTION
Marine mesozooplankton, commonly dominated by copepods, is increasingly recognized to host unique microenvironments in the pelagic macroenvironment (Tang et al., 2010). Copepod guts are oxygen-depleted and low in pH (Tang et al., 2011;Glud et al., 2015), which is remarkable considering their small size and the high ambient oxygen and pH levels. Anoxic guts have otherwise been reported for much larger benthic invertebrates that are exposed to low ambient oxygen levels in sediments (Plante and Jumars, 1992;Stief and Eller, 2006;Stief et al., 2009). Additionally, copepods and their carcasses constitute local sources of particulate and dissolved organic matter and nutrients, which attracts and feeds free-living and animalassociated bacteria (Steinberg et al., 2000;Tang et al., 2006;Saba et al., 2011). Thus, in pelagic ecosystems, copepods and other mesozooplankton represent vastly abundant "microbial hotspots" with unique biogeochemical features (Tang, 2005;Tang et al., 2010;Nuester et al., 2014). In this sense, zooplankters, and in particular their carcasses, resemble sinking phytodetritus aggregates that host dense microbial communities within a low-oxygen microenvironment (Ploug et al., 1997;Glud et al., 2015;Klawonn et al., 2015;Stief et al., 2016). The presence of zooplankton carcasses in the water column results from nonconsumptive mortality, caused by e.g., parasitism or turbulences, which allows extensive microbial degradation of carcasses during their descent to the seabed (Tang et al., 2011). Estimates of the relative abundance of zooplankton carcasses have lately been improved with new staining techniques Elliott and Tang, 2009), which revealed surprisingly high carcass abundances in natural settings, often in the order of 10-30% of the total zooplankton population (Tang et al., , 2014. For marine habitats of tropical and subtropical regions, relative carcass abundances as high as 20-70% ) and 10-90% (Elliott and Tang, 2009) are reported.
Denitrification activity in carcasses of the relatively large copepod Calanus finmarchicus is most pronounced at low ambient O 2 levels (Glud et al., 2015). Evidence is rapidly accumulating that copepods and other mesozooplankton are indeed abundant in and around oceanic oxygen minimum zones (OMZs) (Escribano et al., 2009;Teuber et al., 2013;Wishner et al., 2013;Hirche et al., 2014;Parris et al., 2014). Also the microbial communities associated with copepod carcasses, including those attached to the exoskeleton, are exposed to the low ambient O 2 levels and will likely respond with increased rates of anaerobic metabolism. This will reinforce the functioning of zooplankton carcasses as O 2 -depleted, pelagic hotspots with possible large-scale biogeochemical implications. Globally, OMZs are responsible for up to 50% of the oceanic fixednitrogen loss (N-loss) via denitrification and anammox (DeVries et al.,, 2013). Copepod carcasses may directly contribute to this N-loss via denitrification and possibly anammox (Glud et al., 2015). Model results further suggest that migrating zooplankton actively transports organically bound nitrogen to depth where it is partially released as NH + 4 and may fuel free-living anammox bacteria and thus contribute indirectly to pelagic N-loss (Bianchi et al., 2014). The quantitative importance of this pathway was recently challenged though because zooplankton may downregulate NH + 4 excretion upon exposure to extreme hypoxia or anoxia (Kiko et al., 2016).
Here, we measured diverse microbial N-cycle activities directly associated with the carcasses of abundant groups of zooplankton in a temperature-stratified coastal basin with anoxic bottom water (Golfo Dulce, Costa Rica). This coastal marine OMZ is an environmental showcase for pelagic Nloss through anammox and denitrification activity (Dalsgaard et al., 2003). Copepod and ostracod carcasses were incubated in seawater sampled from three water layers in which different temperature and O 2 conditions prevail. The seawater was amended with 15 NO − 2 to quantify both aerobic and anaerobic pathways of the microbial N-cycle. The basic hypothesis was that zooplankton carcasses represent anoxic microsites at which anaerobic N-cycling can take place. The ultimate goal was to assess the quantitative contribution by zooplankton carcasses to N-loss in a coastal OMZ, which would not be captured by conventional, incubation-based techniques to measure rates of pelagic N-cycling.

Study Site and Sample Collection
Zooplankton and water samples were collected at the northern end of Golfo Dulce (08 • 41 ′ 60 ′′ N, 83 • 23 ′ 24 ′′ W) which corresponds approximately to Station 1 in a previous study (Ferdelman et al., 2006; Figure 1A). This coastal marine basin is up to 215 m deep and characterized by a temperature-stratified water column and anoxic bottom water with high anammox activity (Thamdrup et al., 1996;Dalsgaard et al., 2003;Ferdelman et al., 2006;Morales-Ramírez et al., 2015). All samplings were made at the same time of day (11 am-1 pm) during a 1-week field campaign in January 2015.
Zooplankton for 15 N-incubation experiments was sampled on three occasions by hauling the 0-40 m depth interval with a conical net (0.5 m in diameter, 100 µm mesh size). The catch of 1-2 hauls was kept in 10 L of oxygenated seawater collected at 20 m depth with Niskin bottles and transported to the  Table S2 for detailed list of identified taxa.
laboratory within 1 h. Zooplankton for community analysis was collected on three occasions from three discrete depth intervals (0-40, 40-70, and 70-100 m) using the same net equipped with a remotely controlled closing device. These zooplankton samples were immediately preserved with formaldehyde (4% final concentration) for later taxonomic determination. Seawater for 15 N-incubation experiments was collected at approximately the center of the depth intervals for which the zooplankton community analyses were made (i.e., 20, 60, and 90 m). Multiple water column profiles of O 2 concentration and temperature were measured during the time period in which zooplankton and seawater were sampled, using a CTD (Sea & Sun Technology) equipped with an O 2 microsensor (Revsbech, 1989).

Zooplankton Community Analysis
Zooplankton samples from the three depth intervals were aliquoted with a Folsom plankton splitter, identified to genus or species level, and counted in a Bogorov chamber. For the different samples, the aliquot size was varied between 1/8 and 1/32 to arrive at a minimum count of 200 individuals. Based on these counts and the water volume cleared with the vertical hauls, the in situ abundances of the different zooplankton taxa were calculated.

N-Stable-Isotope Incubations
Microbial N-transformations associated with zooplankton carcasses were measured for copepods and ostracods as the two taxonomic groups of zooplankton that were abundant in all samples. Carcasses were incubated in seawater from three water layers adjusted to the respective in situ temperature and O 2 conditions [i.e., upper water layer: 26 • C/100% air saturation (AS), intermediate water layer: 21.5 • C/0-25% AS, and lower water layer: 17 • C/0% AS] ( Table S1). The resulting six combinations of zooplankton incubations were compared to Frontiers in Marine Science | www.frontiersin.org the respective seawater incubations. Each of these combinations was run in 3-5 replicates in 100-µm-filtered seawater (i.e., in the presence of free-living bacteria), while 1 additional replicate was run in 0.2-µm-filtered seawater (i.e., in the absence of free-living bacteria).
On the day of collection, the zooplankton samples were screened for dead copepods and dead ostracods. Depending on availability, 25-50 copepod and ostracod carcasses were carefully transferred into replicate 20-mL vials filled with 100-µmfiltered seawater from the respective water layer. The copepod community comprised various species (Table S2) and care was taken that every vial contained carcasses representing most of the copepod diversity. In contrast, the ostracod community comprised only one species (Table S2). Dry weight and carbon and nitrogen contents of copepod and ostracod carcasses were obtained for 4 separately collected batches of 25 individuals each. Carcasses were transferred onto pre-combusted glass fiber filters (GC50, Advantec), freeze-dried overnight, weighed, and analyzed on an elemental analyzer (Delta V Advantage, Thermo Scientific) against protein and acetanilide standards.
The aliquoted copepods and ostracods were pre-incubated overnight (12 h) in 100 µm-filtered seawater collected at the respective water depth adjusted to near-in situ temperature and O 2 conditions. This pre-incubation time was chosen to mimic the sinking of zooplankton carcasses through a 30-40 m depth interval (Kirillin et al., 2012) during which carcass degradation may be initiated. The 15 N-incubations were prepared on the next day using acid-washed, sample-rinsed glass bottles (25 mL) equipped with optode spots for contactless O 2 measurements (PyroScience, Germany) and sealed with butyl stoppers to enable sample extraction. Preparation, incubation, and sampling of the glass bottles were made as described by Stief et al. (2016). Briefly, the bottles were filled with the appropriate mixture of oxygen-and helium-flushed seawater (for 15 min each) to arrive at the desired air saturation level of O 2 . This degassing procedure increases the pH of seawater by <0.1 pH units, whereas it increases the pH of freshwater by ∼1 pH unit. 15 N-tracer was added as Na 15 NO 2 (98 atom% 15 N; Sigma-Aldrich) at a final concentration of ∼3 µmol L −1 , which allowed to trace both aerobic and anaerobic pathways of the microbial N-cycle. In situ NO − 2 concentrations in GD were up to 0.72 µmol L −1 at the time of the current study (Padilla et al., 2017). The pre-incubated carcasses were added using a plastic pipette. The bottles were quickly sealed without entrapping gas bubbles and mounted on a plankton wheel at the respective temperature and in darkness. The rotation of the plankton wheel assured that the zooplankton carcasses remained continuously suspended. Oxygen measurements and water sampling were carried out every 2 h for a total of 8 h. Water sampling using two syringes (one for sample extraction and one for replacement of water) served also to maintain the ambient O 2 level within a reasonable range of variation. On the basis of the preceding O 2 measurement, the O 2 concentration in the replacement water was adjusted to a level that would compensate any observed concentration decrease. Oxygen consumption rates were calculated according to Stief et al. (2016).

Nitrogen Analyses
Each 3-mL water sample was split into 1.5 mL for N 2 and N 2 O analyses and 1.5 mL for dissolved inorganic nitrogen (DIN) analyses. The N 2 /N 2 O sample was injected into a heliumflushed and half-evacuated 3-mL exetainer (Labco, U.K.) that contained 50 µL ZnCl 2 (50% w/v) for sample preservation. The DIN sample was immediately frozen at −20 • C. Isotopically labeled dinitrogen ( 15 N-N 2 ) was analyzed in the headspace of exetainer samples on a gas chromatography-isotopic ratio mass spectrometer (Delta V Plus, Thermo Scientific) (Dalsgaard et al., 2012) with the excess above natural abundance calculated according to Nielsen (1992). Many samples taken at the first time point had higher N 2 concentrations than the consecutive samples, probably due to incomplete equilibration between the water phase and the gas present in tiny gas bubbles or within the butyl stoppers, and were thus excluded from rate calculations. Total nitrous oxide was analyzed in the same exetainers on a gas chromatograph (GC 7890, Agilent Technologies). 15 N-labeled NO − 3 (McIlvin and Altabet, 2005), NO − 2 (Füssel et al., 2012), and NH + 4 (Warembourg, 1993) were analyzed with the cadmium reduction/sulfamic acid, sulfamic acid, and hypobromite assay, respectively, followed by 15 N-N 2 analysis on the gas chromatography-isotopic ratio mass spectrometer. The coefficients of variation of replicated 15 NO − 3 , 15 NO − 2 , and 15 NH + 4 standards (n = 7-10) amounted to 4, 5, and 6%, respectively. Total nitrate and nitrite were analyzed spectrophotometrically (Garcia-Robledo et al., 2014). Total ammonium was analyzed with the salicylate method (Bower and Holm-Hansen, 1980).

Rate Calculations
Turnover rates of N compounds ( 15 N-labeled and total N-pools) were calculated from linear concentration changes during the incubation period and corrected for the dilution due to repeated sampling . Rates determined in seawater incubations with seawater only were subtracted from those determined in zooplankton incubations to arrive at individualspecific process rates. The 15 N-turnover rates were used to calculate the process rates of NO − 2 oxidation, denitrification, and dissimilatory nitrate reduction to ammonium (DNRA): (a) NO − 2 oxidation rate equals 15 NO − 3 production rate divided by 15 N-labeling percentage (FN), (b) denitrification rate equals 30 N 2 production rate divided by FN 2 , and (c) DNRA rate equals 15 NH + 4 production rate divided by FN. For different processes and water depths, FN values of two different N-pools were considered: (a) for the NO − 2 oxidation rate and for the upper water layer, the measured 15 N-labeling percentage of the NO − 2 pool (FN NO2 ) was used and (b) for denitrification and DNRA and for the intermediate and lower water layers, the measured 15 Nlabeling percentage of the NO x (i.e., NO − 3 + NO − 2 ) pool (FN NOx ) was used. The latter option was chosen because in incubations with seawater from the intermediate and lower water layers, strong NO − 3 consumption accompanied 15 NO − 2 consumption, which diluted the 15 N-labeled NO x pool. Nitrite produced by NO − 3 consumption may bypass the ambient NO − 2 pool and may be directly turned over in the denitrification and DNRA reactions (De Brabandere et al., 2014).
Anammox rates (AMX total ) were calculated as: where P 29 and P 30 denote the production rates of singlelabeled and double-labeled N 2 , respectively (Thamdrup and Dalsgaard, 2002). This calculation neglects the contribution by anammox to 30 N 2 production resulting from the accumulation of 15 NH + 4 through DNRA, which is justified because the 15 Nlabeling percentage of the NH + 4 pool remained ≤4% . Previous studies found a tendency for anammox rates determined with 15 NO − 2 to be higher than rates determined with 15 NH + 4 (De Brabandere et al., 2014). In a comparison of anammox rates determined with either of the two labeled substrates during our expedition, the 15 NO − 2 -based rates at 90 m, the only depth where we detected anammox, were only slightly (14%) higher than those based on 15 NH + 4 (Laura A. Bristow, pers. comm.). Thus, the potential bias resulting from the use of 15 NO − 2 in our study should not have any substantial impact on the results and conclusions.
In incubations with seawater from the intermediate water layer, high consumption rates of NO − 3 total and NH + 4 total reduced the net production rates of 15 NO − 3 and 15 NH + 4 . The concentration time series of 15 NO − 3 and 15 NH + 4 were thus corrected to allow calculation of gross production rates: where C denotes concentrations of 15 N-labeled ( 15 N) or total N-pools (N total ) measured at time point t (i) or the following time point t (i+1) .

Statistics
Linear regression analysis was used to test, if individual concentration time series had a slope significantly different from zero. The process rate derived from a concentration time series that did not meet this criterion was taken as zero. One-sample ttests were used to test, if mean process rates were different from zero at a significance level of p = 0.05 (one-tailed). Grubbs' test for outliers was used to test, if incubations in 0.2-µm-filtered seawater had significantly lower process rates than incubations in 100-µm-filtered seawater. Student's t-tests were used to identify differences in individual-specific process rates between copepods and ostracods within a given water layer at a significance level of p = 0.05 (two-tailed). One-way ANOVA was used to identify differences in individual-specific process rates between the three water layers within a given taxonomic group of zooplankton at a significance level of p = 0.05 (two-tailed), followed by Holm-Sidak post-hoc tests. A non-parametric correlation analysis was conducted for N 2 O production rates vs. turnover rates of other N-cycle processes.

Environmental Settings in Golfo Dulce
At the time of the study, the water column at the targeted station in Golfo Dulce (GD) was temperature-stratified with the thermocline located at ca. 30-50 m water depth (Figure 2A). The water column was further divided into an upper oxic, an intermediate hypoxic/anoxic, and a lower anoxic layer ( Figure 2B). Oxygen concentrations were relatively stable directly at the surface and below 60 m water depth. In the depth interval 20-60 m, however, O 2 concentration showed large variations between the five sampling days. Additionally, rapid fluctuations between hypoxic and almost fully anoxic conditions were observed for individual sampling days on which O 2 profiles were repeated 2-3 times within a time period of 4 h.
Zooplankton abundance was high in the upper depth interval and progressively decreased to 10-fold lower levels in the lower depth interval (Figure 2C). The zooplankton community was dominated by small copepods making up ca. 50% of the total zooplankton abundance, whereas ostracod abundance was ca. 10-fold lower. The most common copepods were Oncaea sp., Oithona sp., and unidentified copepodites; the only ostracod identified was Cypridina sp. which had a similar body size as the most common copepods (Figures 1B,C; Table S2).

Oxygen Dynamics in Zooplankton and Seawater Incubations
Oxygen concentrations in zooplankton and seawater incubations were relatively stable for the oxic and anoxic water layer (Figures S1A,C). Strong fluctuations in O 2 concentration were observed in zooplankton incubations in water from the hypoxic/anoxic layer, with anoxic conditions occurring at three out of nine sampling time points ( Figure S1B), which resembled the O 2 dynamics observed in the intermediate water layer in GD ( Figure 2B).
Carcass-associated O 2 consumption rates (O 2 -Cons ZP ) were always significantly different from zero and were not significantly different between copepods and ostracods within any given water layer ( Figure S2a,e; Tables S3-S5). However, for both the copepod and the ostracod carcasses, O 2 -Cons ZP rates were significantly higher in the hypoxic/anoxic than in the oxic water layer (Table S6).

Nitrogen Cycling Associated with Zooplankton Carcasses
Generally, zooplankton incubations in 0.2-µm-filtered seawater (white circles in Figure 3, Figures S2, S3) did not result in lower-end, statistical outliers (p > 0.05, Grubbs' test for n = 72 individual measurements), which suggests that the observed N-cycling in the zooplankton incubations was due  Table S2 for detailed list of identified taxa. to carcass-associated bacteria rather than free-living bacteria in the ambient seawater. However, it cannot be ruled out that (i) carcass-associated bacteria were detaching during the 15 Nincubations and remained actively involved in N-cycling as free-living bacteria and (ii) zooplankton carcasses were leaking substrates into the surrounding water and thereby promoted Ncycling by free-living bacteria. In both cases though, significant contributions by free-living bacteria to anaerobic N-cycling would have only been possible under fully anoxic conditions in the incubation bottles. The rates of the diverse carcass-associated N-cycling pathways showed several clear patterns (Figure 3; Table S3): (i) Carcassassociated N-cycling rates were relatively low in the oxic water layer and significantly different from zero only for N 2 O production (N 2 O-Prod ZP ) and DNRA activity (DNRA ZP ) (Table S4), (ii) carcass-associated N-cycling rates were in most cases not different between copepods and ostracods within any given water layer (Table S5), and (iii) the rates of carcassassociated NO − 2 -oxidation (NO 2 -Ox ZP ), NO − 3 -reduction (NO 3 -Red ZP ), N 2 O-Prod ZP , denitrification (DEN ZP ), and anammox (AMX ZP ) were highest in the hypoxic/anoxic water layer, while only DNRA ZP rates increased consistently with water depth (Table S6).
In addition to the consumption of the added 15 NO − 2 tracer ( Figures S3b,f,j), high rates of NO 3 -Red ZP activity were observed in the hypoxic/anoxic and anoxic water layers (Figures S2f,j). In the hypoxic/anoxic water layer, NO 3 -Red ZP scaled with the net production of NO − 2 total (Figures S2f,g), which suggests high rates of carcass-associated dissimilatory nitrate reduction to nitrite (DNRN ZP ).
Carcass-associated N 2 O production (N 2 O-Prod ZP ) was significantly correlated with DNRA ZP in the oxic water layer, with NO 2 -Ox ZP and NO 3 -Red ZP in the hypoxic/anoxic water layer, and with DEN ZP in the anoxic water layer ( Figure S4, Table S7).
The partitioning between the diverse N-cycling pathways was strikingly different for zooplankton and seawater incubations (Figure 4). The relative share of DNRA was without exception larger in zooplankton than in seawater incubations, with the absolute rate and the relative share of DNRA ZP being particularly high in the anoxic water layer. Denitrification was most important in carcasses incubated in seawater from the hypoxic/anoxic layer. The relative share of N 2 O-Prod ZP was low in all zooplankton incubations, but N 2 O production was dominant in seawater incubations for the hypoxic/anoxic layer. Anammox activity was most important in seawater incubations for the anoxic layer and was not detected in the other two water layers. Only for copepod and ostracod carcasses incubated in seawater from the hypoxic/anoxic layer, a small relative share of AMX ZP was observed.

Direct and Indirect Contribution by Zooplankton Carcasses to Pelagic N-Loss
Zooplankton carcasses may contribute to pelagic N-loss directly by hosting bacteria that produce N 2 and N 2 O and indirectly by hosting bacteria that produce NH + 4 which, together with NH + 4 released from carcass degradation, may fuel anammox by free-living bacteria. The contribution by zooplankton carcasses to the combined rate of pelagic N 2 O production, denitrification, DNRA, and anammox was estimated from the individual-specific N-cycling rates (Figure 3) and the in situ abundance of copepods and ostracods in GD (Figure 2C), assuming a relative carcass abundance of 36%, which corresponds to the mean of minimum and maximum values reported in a recent meta-study (Tang et al., 2014). Due to the dominance of copepods in GD, significant contributions to pelagic N-cycling are almost exclusively linked to this taxonomic group of zooplankton (Figure 5). In the oxic water layer, it is in fact only the presence of copepod carcasses, which enables the occurrence of N 2 O production, denitrification, DNRA, and anammox, albeit at relatively low rates. In the hypoxic/anoxic water layer though, the estimated N-cycling rates associated with copepod carcasses rival or even exceed those measured in seawater only, while the contribution by ostracod carcasses is small. In contrast, in the anoxic water layer, zooplankton carcasses are estimated to only marginally enhance the rather intense pelagic N-cycling, except for a strong contribution by copepod carcasses to DNRA rates. However, the estimated contribution by zooplankton to pelagic N-loss may actually be higher than observed here, if also live specimens contribute to pelagic N-cycling, e.g., through microbial processes occurring in their anoxic guts.
These contributions by zooplankton carcasses to pelagic Ncycling correspond to or may translate into an enhancement of pelagic N-loss. N 2 O-Prod ZP , DEN ZP , and AMX ZP correspond to  Figure 3 were extrapolated to the in situ abundance of copepods and ostracods presented in Figure 2C and assuming a relative carcass abundance of 36% (Tang et al., 2014). N-cycling rates measured in seawater incubations are shown for comparison. Each panel is split into an upper and a lower part to account for the low rates measured in the oxic water layer. Means and standard error of 4-6 replicate measurements are shown. SW = 0 indicates rates in seawater incubations not significantly different from zero; * indicates rates in zooplankton incubations lower than 0.02 nmol N L −1 d −1 .
N-loss in all water layers, while DNRA ZP and carcass degradation might contribute to N-loss only in the anoxic water layer by fueling free-living anammox bacteria with NH + 4 . The most pronounced enhancement of pelagic N-loss is estimated to occur in the hypoxic/anoxic water layer where zooplankton abundance is still relatively high and pelagic N-loss proceeds solely through ambient N 2 O production, but not denitrification and anammox. Here, carcass-associated N-cycling is estimated to increase pelagic N-loss by a factor as high as 1.4 (Table 1). In the anoxic water layer, however, carcass-associated N-cycling is estimated to increase the intense pelagic N-loss by a factor of 1.02 only. For the oxic water layer, pelagic N-loss was not detected in seawater incubations and therefore it is not possible to calculate an enhancement factor for the anyway low carcass-mediated N-loss.
In the anoxic water layer, NH + 4 total production (NH 4 -Prod ZP ) due to carcass degradation and DNRA ZP was on average 1.9-fold higher than NH + 4 production by DNRA ZP alone (Figures S2l,  S5a). This multiplication factor further increases the estimates given in Figure 5C and is already accounted for in Table 1. The degradation-related origin of this NH + 4 total is confirmed by the significant correlation between NH 4 -Prod ZP and O 2 -Cons ZP , which suggests C/N ratios of organic matter degradation of 9.2 ( Figure S5b). This value is similar to the C/N ratios calculated for the biomass of the copepod and ostracod carcasses (see above). Notably, in the hypoxic/anoxic water layer, a net uptake of NH + 4 total by zooplankton carcasses was observed ( Figure  S2h), even though the high O 2 -Cons ZP rates ( Figure S2e) would suggest high NH 4 -Prod ZP rates due to carcass degradation. This net uptake of NH + 4 total by zooplankton carcasses is most likely due to intense N-assimilation by rapidly growing bacteria Depth-integrated rates of pelagic N-loss (mean ± standard error) were calculated as the sum of volumetric rates of N 2 O production, denitrification, and anammox measured in seawater incubations multiplied with the thickness of the respective depth interval. The contribution by zooplankton carcasses to pelagic N-loss was estimated for a relative carcass abundance of 36% (Tang et al., 2014), using data presented in Figure 5. associated with the carcasses, as indicated by a 2-to 5-fold increase of O 2 -Cons ZP rates during the 8-h incubation period, which was not observed in either of the other water layers (data not shown).

Carcass-Associated N-Loss in the Presence of Oxygen
The 15 N-incubation experiments conducted under different environmental conditions imply that the relative enhancement of pelagic N-loss by zooplankton carcasses is greatest in oxic Frontiers in Marine Science | www.frontiersin.org 8 May 2017 | Volume 4 | Article 152 FIGURE 6 | Conceptual scheme of zooplankton-carcass-associated N-cycling and potential effects on pelagic N-cycling. Microbial N-cycling may take place in the anoxic interior and on the exoskeleton of zooplankton carcasses. Anaerobic N-cycling associated with zooplankton carcasses contributes to pelagic N-loss both directly via N 2 and N 2 O production (through denitrification and anammox) and indirectly via NH + 4 production (through DNRA and mineralization) that may fuel free-living anammox bacteria in anoxic water layers. Carcass-derived NH + 4 may also fuel microbial nitrification in oxic and hypoxic water layers. and hypoxic water layers. In the presence of O 2 , the N-loss mediated by free-living bacteria is generally very limited. At ambient O 2 levels higher than 1-20 µmol L −1 , denitrification and anammox activities by free-living bacteria are usually inhibited (Kalvelage et al., 2011;Dalsgaard et al., 2014) and bacterial or archaeal nitrifiers contribute to N-loss only marginally through N 2 O production (Frame and Casciotti, 2010;Löscher et al., 2012). Thus, under oxic and hypoxic conditions, zooplankton carcasses must act as (partially) anoxic pelagic hotspots that host microbial N-cycling pathways which directly contribute to N-loss, namely N 2 O production, denitrification, and anammox ( Figure 6).
In the oxic water layer, N-loss via N 2 O and N 2 production was exclusively enabled by the presence of copepod and ostracod carcasses because the corresponding background rates in seawater were not significantly different from zero. The individual-specific N 2 O and N 2 production rates were very low though and added up to only ∼0.25 nmol N L −1 d −1 , despite the high abundance of copepods and ostracods in this water layer. In comparison, the N-loss rate measured in seawater incubations for the anoxic layer was >100 nmol N L −1 d −1 , which was similar to the N-loss rate measured in anoxic bottom water of GD in an earlier study (Dalsgaard et al., 2003).
In the hypoxic/anoxic water layer, individual-specific N 2 O and N 2 production rates were 1-2 orders of magnitude higher than in the oxic layer. These increased rates are likely explained by the expansion of the anoxic volume inside the carcasses due to the lowered ambient O 2 level as previously shown for sinking diatom and cyanobacterial aggregates (Klawonn et al., 2015;Stief et al., 2016;Lundgaard et al., 2017) and large copepod carcasses (Glud et al., 2015). Additionally, the occurrence of short anoxic phases during the incubation period may have boosted anaerobic N-cycling. It can also not be ruled out that the higher individual-specific N-cycling rates are linked to differences in the bacterial community composition between the oxic and the hypoxic/anoxic water layer.
Taking the relatively high copepod and ostracod abundances in the hypoxic/anoxic water layer into account, the extrapolated carcass-associated N-loss was here 14-fold higher than in the oxic water layer. Strikingly, the presence of carcasses is estimated to increase pelagic N-loss 1.4-fold in this water layer where the background rates in seawater already reached ∼10% of the Nloss in the anoxic water layer. In the hypoxic/anoxic water layer, pelagic N-loss proceeded exclusively through N 2 O rather than N 2 production, which is commonly observed in low-oxygen aquatic environments: N 2 O yields of nitrification and denitrification are often high at low ambient O 2 levels (Goreau et al., 1980;Bonin and Raymond, 1990), while N 2 production via denitrification and anammox may be completely inhibited above low-micromolar levels of O 2 (Kalvelage et al., 2011;Dalsgaard et al., 2014). It needs to be noted though that the seawater incubations did not reach anoxic conditions and thus the rates of anaerobic N-cycling may be underestimated relative to the zooplankton incubations.
In the presence of O 2 , carcass-associated N-loss proceeded mainly through DEN ZP and to a smaller extent through N 2 O-Prod ZP and AMX ZP . Elevated rates of DEN ZP at lower ambient O 2 levels have previously been reported for carcasses of the copepod C. finmarchicus (Glud et al., 2015) that has a much larger prosome length (∼2,600 µm) than Oncaea sp. and Oithona sp. (<500 µm). Apparently, reduced ambient O 2 levels turn the copepod carcasses into more and more O 2 -depleted microsites enhancing the potential for anaerobic N-cycling. In the current study, the DEN ZP rate measured for copepod carcasses under hypoxic/anoxic conditions (0-100 µmol O 2 L −1 ) and at 21.5 • C amounted to 3.3 nmol N individual −1 d −1 . The much larger C. finmarchicus carcasses reached the same rate at ∼10 µmol O 2 L −1 and at 7 • C (Glud et al., 2015). Moreover, the DEN ZP rate of C. finmarchicus carcasses increased exponentially to >50 nmol N individual −1 d −1 at 1 µmol O 2 L −1 (Glud et al., 2015).
Carcass-associated N-loss is also mediated by N 2 O-Prod ZP , albeit at lower rates than by DEN ZP . The origin of the N 2 O produced by zooplankton carcasses remains unresolved because only N 2 O total was measured. The most prominent sources of N 2 O are nitrification (strictly speaking: ammonia oxidation) and denitrification. Nitrous oxide production associated with larger invertebrates has been ascribed to both nitrification and denitrification (Svenningsen et al., 2012;Heisterkamp et al., 2013). As in their benthic counterparts, the live copepods and ostracods may host denitrification activity in their anoxic gut and nitrification activity on their exoskeleton. Both processes may also be involved in N 2 O production by carcasses, provided that these possess an anoxic interior and an oxic body surface. Indeed, N 2 O-Prod ZP was significantly correlated with both NO 2 -Ox ZP and NO 3 -Red ZP in the Frontiers in Marine Science | www.frontiersin.org 9 May 2017 | Volume 4 | Article 152 hypoxic/anoxic water layer, which is compatible with both NO − 3 -producing and NO − 3 -consuming pathways as possible sources of N 2 O.
Carcass-associated anammox activity was detected in the hypoxic/anoxic water layer, but not in the anoxic water layer. This is surprising given the pronounced O 2 sensitivity of anammox bacteria and also because a prevalence of anammox bacteria as colonizers of the carcasses would rather be expected in the anoxic water layer. However, anammox bacteria were also identified and anammox activity was detected in a sponge that hosts both oxic and anoxic microenvironments (Hoffmann et al., 2009).
In the hypoxic/anoxic water layer, high rates of DNRN ZP might be inferred from the co-occurring NO − 3 consumption and NO − 2 production at an almost equimolar ratio. Nitrite produced by DNRN ZP may theoretically fuel the anammox process (De Brabandere et al., 2014), which was not detected in seawater incubations for this water layer though.
The estimated 1.4-fold enhancement of N-loss by zooplankton carcasses in the hypoxic/anoxic water layer is based on a relative carcass abundance of 36% and on only about half of the total community of mesozooplankton in GD. Additionally, it can be assumed that also live specimens of copepods and ostracods host anaerobic N-cycling pathways, the rates of which still need to be quantified though. Thus, the 1.4-fold enhancement of N-loss represents a relatively conservative estimate. Irrespective of the exact extent of carcass-associated N-loss at high ambient O 2 levels, it is likely to be overlooked in studies that address N-loss in OMZs. First, pre-filtration of seawater for rate measurements will remove the carcasses; second, if non-filtered seawater is used, only a sparse number of incubation vials will contain a single carcass; third, oxic and hypoxic water layers are often excluded from sampling schemes because N-loss mediated by free-living bacteria is not expected to occur in the presence of O 2 .

Carcass-Associated N-Loss in the Absence of Oxygen
Under fully anoxic conditions, DNRA was clearly the dominant N-cycling pathway associated with zooplankton carcasses. Also under oxic and hypoxic conditions, DNRA ZP rates were high compared to DEN ZP , AMX ZP , and N 2 O-Prod ZP rates. The crucial difference between DNRA ZP proceeding in the presence or absence of O 2 is that exclusively in the anoxic bottom water of GD, the carcass-derived NH + 4 might fuel the activity of free-living anammox bacteria. It may thus be speculated that zooplankton carcasses sinking through the anoxic bottom water may contribute to N-loss to a large extent indirectly. This mechanism seems plausible for GD because NH + 4 concentrations were below 50 nmol L −1 in the anoxic, non-sulfidic water layer at the time of the current study (Laura A. Bristow, pers. comm.). Additionally, it has been shown that the in situ anammox activity in GD can be stimulated by experimental NH + 4 addition (Dalsgaard et al., 2003). Ammonium produced by DNRA ZP adds to the NH + 4 excreted by live zooplankton or produced by the degradation of carcasses and fecal pellets (Figure 6). Our results indicate that for the copepods and ostracods in GD, DNRA ZP approximately doubles the rate of NH + 4 production due to carcass degradation and thus significantly enforces the potential indirect contribution to pelagic N-loss.
It has previously been suggested that migrating zooplankton supplies NH + 4 to the anammox process in OMZs (Bianchi et al., 2014). The underlying mechanism is known as "active transport" of substrates and nutrients by migrating zooplankton (Steinberg et al., 2000;Schnetzer and Steinberg, 2002). Organic matter export to depth is thought to be more efficiently mediated by migrating zooplankton that congregate at the boundary or within OMZs during daytime than by sinking organic aggregates that are rapidly degraded already during their descent (Kalvelage et al., 2013;Bianchi et al., 2014). However, NH + 4 excretion by live zooplankton may be down-regulated due to anoxic conditions (Kiko et al., 2016). Our findings demonstrate though that NH + 4 release by zooplankton is also due to carcass degradation and carcass-associated DNRA activity, in particular if specimens die close to their daytime migration depth. Notably, DNRA ZP rates were highest under anoxic conditions where NH + 4 excretion by live specimens may be particularly low.
The carcasses might be preferentially colonized by DNRA bacteria that were residing in the gut of the live specimens. Anoxia and high availability of organic carbon relative to NO − 3 in the gut of invertebrates may favor DNRA over denitrification (Giblin et al., 2013). Interestingly, many of the bacterial groups recurring in the different copepod microbiomes (i.e., Bacteroidetes, Firmicutes, Actinobacteria, Vibrio sp.) comprise representatives that colonize the gut of other animals and also possess nrf, the key gene for the DNRA process (Mohan et al., 2004;Welsh et al., 2014;Decleyre et al., 2015).
An estimate for the anoxic bottom water in GD indicates that the extrapolated rate of NH 4 -Prod ZP (incl. DNRA ZP ) is ∼2-fold higher than the DNRA rate measured in seawater. This additional NH + 4 supplied by carcasses corresponds to 6.3% of the anammox rate measured in seawater (105 nmol N L −1 d −1 ). Rates of denitrification and N 2 O production in seawater and zooplankton incubations were insignificant compared to the DNRA rates and thus the zooplankton carcasses potentially enhance N-loss in the anoxic water layer only marginally.

Ecological Implications
Integrated over the top 100 m of the stratified water column of GD, zooplankton carcasses potentially enhance N-loss between 1.03-fold (considering direct contributions through N 2 and N 2 O production) and 1.05-fold (considering direct and indirect contributions through N 2 , N 2 O, and NH + 4 production). The largest fraction of direct contributions is to be expected in the hypoxic/anoxic water layer (1.4-fold increase), while indirect contributions come into play in the anoxic water layer only and are still minor compared to background rates in the ambient water (1.02-fold increase). Thus, zooplankton carcasses have mainly the potential to functionally extend the anoxic volume of OMZs at their hypoxic boundaries. The low levels of ambient O 2 in these zones allow relatively high zooplankton abundances and very high individual-specific N-cycling rates, but efficiently inhibit pelagic N-loss through anaerobic N-cycling processes mediated by free-living bacteria.
However, the vertical extent of the oxycline in GD is relatively small (∼20 m) compared to the thickness of the anoxic bottom water where the bulk of pelagic N-loss is mediated by free-living bacteria (∼100 m). Therefore, the significant contribution by zooplankton carcasses to pelagic N-loss in the hypoxic/anoxic water layer does not translate into a major enhancement of depthintegrated N-loss in GD. In contrast, open-ocean OMZs possess upper and lower hypoxic boundaries spanning depth intervals of several hundreds of meters and regularly host high zooplankton abundances (Wishner et al., 2013;Hirche et al., 2014). In such settings, the relative contribution by zooplankton carcasses to depth-integrated N-loss is likely much larger than observed in GD. The relative carcass abundance in zooplankton in and around OMZs has to our knowledge not been investigated so far. It may be assumed that relative carcass abundance increases with decreasing ambient O 2 level due to non-predatory mortality. In that case, the contributions by carcasses to pelagic N-cycling will gain additional importance in the hypoxic and anoxic zones of the OMZs.
The functional extension of the anoxic volume of OMZs has previously been proposed to be mediated by sinking diatom aggregates Stief et al., 2016;Lundgaard et al., 2017), but may as well be caused by other pelagic microbial hotspots with an anoxic center, such as sinking cyanobacterial aggregates, live specimens of zooplankton, and fecal pellets (Tang et al., 2011;Glud et al., 2015;Klawonn et al., 2015). Globally, OMZs are projected to expand significantly within the next century due to climate change (Stramma et al., 2008;Keeling et al., 2010). Thus, anaerobic N-cycling inside pelagic microbial hotspots will most likely intensify and increase the N-loss from the ocean, aside from any possible negative feedbacks (Kalvelage et al., 2011).
The contribution by zooplankton carcasses to N-loss through N 2 production might be counteracted by carcass-associated N 2 fixation activity. However, individual-specific N 2 production rates measured in the current study were at least three orders of magnitude higher than N 2 fixation rates measured for small copepods in another study (Scavotto et al., 2015). Despite the possible variability of N-cycling associated with mesozooplankton from different habitats, this large difference in rates suggests that carcass-associated microbial communities exhibit net consumption rather than net production of fixed nitrogen in the pelagic zone.
Copepods and ostracods are by far not the only groups of zooplankton that may host anaerobic N-cycle activities. In GD, other zooplankton species might also be of importance as colonization sites for N-cycle bacteria and as potentially anoxic microsites. The gastropod Cuvierina sp. reached very high abundances in GD during the time of the current study and might host microbial biofilms involved in N-cycling on their shell (Heisterkamp et al., 2013). In freshwater ecosystems, copepods and daphnids are abundant and may qualify as pelagic microbial hotspots ). The large Calanus species need more attention regarding the diverse N-cycle processes they may host besides denitrification (Glud et al., 2015). In the global ocean, krill and copepods compete for the status as the most abundant multicellular animals on Earth. Due to their large body size, krill most likely possess anoxic guts and might thus also represent anoxic pelagic microsites, despite high ambient O 2 levels in their natural habitat. Establishing the pathways and magnitude of N-cycle activities associated with these abundant groups of zooplankton will allow first estimates of their quantitative contribution to pelagic N-cycling in the global ocean.

AUTHOR CONTRIBUTIONS
PS and RG designed the study. PS, AL, AM, and BT carried out sample collection, incubation experiments, sample analysis, and data analysis. PS wrote the manuscript with input from all co-authors.