Increased Feeding and Nutrient Excretion of Adult Antarctic Krill, Euphausia superba, Exposed to Enhanced Carbon Dioxide (CO2)

Abstract

Ocean acidification has a wide-ranging potential for impacting the physiology and metabolism of zooplankton. Sufficiently elevated CO2 concentrations can alter internal acid-base balance, compromising homeostatic regulation and disrupting internal systems ranging from oxygen transport to ion balance. We assessed feeding and nutrient excretion rates in natural populations of the keystone species Euphausia superba (Antarctic krill) by conducting a CO2 perturbation experiment at ambient and elevated atmospheric CO2 levels in January 2011 along the West Antarctic Peninsula (WAP). Under elevated CO2 conditions (∼672 ppm), ingestion rates of krill averaged 78 µg C individual−1 d−1 and were 3.5 times higher than krill ingestion rates at ambient, present day CO2 concentrations. Additionally, rates of ammonium, phosphate, and dissolved organic carbon (DOC) excretion by krill were 1.5, 1.5, and 3.0 times higher, respectively, in the high CO2 treatment than at ambient CO2 concentrations. Excretion of urea, however, was ∼17% lower in the high CO2 treatment, suggesting differences in catabolic processes of krill between treatments. Activities of key metabolic enzymes, malate dehydrogenase (MDH) and lactate dehydrogenase (LDH), were consistently higher in the high CO2 treatment. The observed shifts in metabolism are consistent with increased physiological costs associated with regulating internal acid-base equilibria. This represents an additional stress that may hamper growth and reproduction, which would negatively impact an already declining krill population along the WAP.

About the Authors:

Grace K. Saba

* E-mail: [email protected]

Affiliation: Institute of Marine and Coastal Sciences, Rutgers University, New Brunswick, New Jersey, United States of America

Oscar Schofield

Affiliation: Institute of Marine and Coastal Sciences, Rutgers University, New Brunswick, New Jersey, United States of America

Joseph J. Torres

Affiliation: College of Marine Science, University of South Florida, St. Petersburg, Florida, United States of America

Erica H. Ombres

Affiliation: College of Marine Science, University of South Florida, St. Petersburg, Florida, United States of America

Deborah K. Steinberg

Affiliation: Virginia Institute of Marine Science, College of William & Mary, Gloucester Point, Virginia, United States of America

Introduction

The Antarctic krill, Euphausia superba, is a key species in Antarctic food webs [1], [2]. E. superba is a major consumer of phytoplankton [3]–[5] and a primary food source for many of the top predators in the Southern Ocean including baleen whales, seals, penguins, and flighted sea birds [6]. Through their feeding and excretion processes, krill are a major source of regenerated nutrients, which in turn support phytoplankton growth [7]–[9]. However, krill along the West Antarctic Peninsula (WAP) region, the northernmost part of the mainland of Antarctica extending into the Southern Ocean, have declined two-fold since the mid-1970s due to profound changes along the WAP in the past decades [2]. The west coast of the northern WAP is changing from a cold, dry polar climate to a warmer, humid subantarctic climate [10] as a result of a 6°C increase in mid-winter surface atmospheric temperatures (>than 5× the global average) in the past 50 years [11], [12]. This rapid warming is increasing the heat content of seawater over the shelf [13] and reducing both amount and duration of sea ice [14], [15]. Concurrent changes in the biomass and composition of krill food source have occurred [10], [16]. WAP summertime chlorophyll a (Chl a) has declined by 12% over the past 30 years [10], and, in the northern WAP region, there has been a non-uniform shift in the phytoplankton size fraction from large diatoms to small cells (<20 µm) [10], [17], on which krill can not efficiently feed [18]–[22]. Ocean acidification poses an additional threat to krill populations, as there are predictions that by the end of this century, the Southern Ocean will be the first region to be affected by seawater chemistry changes associated with enhanced carbon dioxide (CO2) [23], [24].

Human activities have driven the rapid 40% increase in atmospheric carbon dioxide CO2, from preindustrial levels of 280 ppm (parts per million) to current levels of nearly 397 ppm [25], [26]. Present-day atmospheric CO2 concentration is projected to double by the end of the 21st century [25], [27]. Nearly one-third of emitted anthropogenic CO2 is absorbed by the oceans [28], [29], resulting in reductions in seawater pH and alterations in carbonate chemistry (i.e., reductions in carbonate, CO32−, ions). The current rapid rate of oceanic CO2 uptake, one million metric tons of CO2 per hour [30], is paralleled by rates of acidification at least ten times faster than any change seen in the fossil record over the past 65 million years [31]. These rapid changes are expected to cause adverse ecosystem wide effects [32]–[35].

Not only are Antarctic krill exposed to seasonal fluctuations of seawater pCO2/pH [24], they can migrate between the surface and at depth during ontogenetic migration (700–1000 m; [36], [37]) and diurnal vertical migration (400+m; [38], [39]). As such they are currently exposed not only to a wide range of seawater pCO2/pH on short (∼daily) time scales, but also to hypercapnic water at depth [40]. Model projections using the IPCC IS92a scenario demonstrated that Southern Ocean seawater pCO2, within the depth range utilized by Antarctic krill, could rise up to ∼1400 ppm by the year 2100 [40].

Elevated seawater CO2 can impact marine organisms both via decreased carbonate saturation that affects calcification rates and via disturbance to acid-base (metabolic) physiology [23], [41], [42]. Organisms have different responses to hypercapnia, the CO2-induced acidification of body fluids. Many studies have demonstrated either no effect of CO2 on metabolism of organisms [43] or reduced metabolic activity under elevated CO2 [34], [41], [44]–[48]. Some organisms in environments predisposed to high fluctuations in pH and CO2 over short time scales, such as coastal upwelling regions, cannot compensate and suppress metabolism when they encounter low pH [34], [41], [44], [45]. Other organisms can fully compensate extracellular fluid pH; however, as a result acid-base and ion equilibria reach new steady state values [45], [49]. The result of this adaptation is that there are extra costs of compensation. For instance, organisms may have higher demands for acid-base regulator proteins (e.g., [50]) and would have to work harder to maintain or alter internal acid-base equilibria. Furthermore, their oxygen transport system may be compromised [41], [51], making them less effective at picking up oxygen (O2) and forcing them to process more water to extract the O2 they demand. For example, E. superba have a pH sensitive respiratory protein that could impair oxygen transport depending on blood buffering [52]. The total amount of oxygen that can be carried in the blood in support of routine activities is reduced by half in E. superba under an increase in pCO2 to 1000 ppm (drop in arterial pH from 8.1 to 7.9) [52]. The two strategies (incomplete or complete compensation) may affect the energetics and performance of an organism differently.

Most studies testing CO2/pH effects on crustacean zooplankton thus far have focused on growth, development, or mortality of various life stages [40], [53], [54], and a majority of them show decreased hatching success, irregular larval development, or decrease in larval size under conditions of high CO2 (elevated above ambient) or low pH (below ambient). The subtle effects on the physiology and metabolism of marine zooplankton due to ocean acidification are extremely understudied. For example, no studies have directly measured zooplankton feeding, nutrient release, or metabolism under CO2 levels predicted for the future ocean. In one study, however, amphipods exposed to high CO2 concentrations (low pH), exhibited an increase in the expression of the metabolic enzyme glyceraldehyde-3-phosphate dehydrogenase gene (gapdh gene), suggesting that metabolic changes occurred in response to acidification [55]. Increases in metabolic enzyme expression and ventilatory frequency and effort have also been demonstrated for some fish and elasmobranchs [56]–[58] and brittle stars [59]. Such increases may be attributed to enhanced metabolic costs. Compensation costs of enhanced CO2 may cause long-term shifts in respiration and metabolic equilibria and eventually hamper growth and reproduction of organisms [41], [42], and may eventually negatively impact an already declining krill population.

We conducted a field CO2 perturbation experiment along the WAP during the austral summer (January 2011) to determine krill feeding and nutrient excretion rates at ambient (∼390 ppm) and “high” (750 ppm) atmospheric CO2 concentrations. We selected the “high” CO2 level of 750 ppm because it represents the mean predicted atmospheric CO2 concentration by 2100. The experimental design strategy of selecting ambient and 750 ppm when using two CO2 treatments is also what is recommended in the Guide to Best Practices for Ocean Acidification Research and Data Reporting [60]. However, in the present/future oceans krill are/will be exposed to large fluctuations in pCO2 and pH due to strong deviations of seawater chemistry with atmospheric CO2 concentration [61]; thus, the results from this study likely reflect mean responses to enhanced CO2. We hypothesized that under conditions of high CO2 in perturbation experiments, krill grazing rate and metabolism (nutrient excretion, metabolic enzyme activity) would be increased, reflecting extra costs of compensation due to maintaining internal acid-base balance.

Materials and Methods

Ethics statement

All animal work has been conducted according to relevant national and international guidelines. No specific permits were required for the described field studies per articles of The Antarctic Treaty and the National Science Foundation, which is charged with enforcement of the Treaty in the USA and its possessions. The location (Antarctica) is protected under the provisions of The Antarctic treaty, but these provisions do not extend to valid scientific research supported by a recognized national research program of a signatory nation. No endangered or protected species were involved.

CO2 Perturbations

Seawater for sampling the initial phytoplankton composition in the water column was collected at various depths, ∼100 km offshore of Adelaide Island along the West Antarctic Peninsula in January 2011 (66.51°W, 69.87°S) aboard the R/V Laurence M. Gould using Niskin bottles affixed to a conductivity, temperature, depth (CTD) profiler. A known volume of collected seawater from each depth was filtered onto a GF/F filter and flash frozen for HPLC pigment analysis. The taxonomic composition of the phytoplankton assemblages was derived quantitatively from an analysis of HPLC pigment data using CHEMTAX (V195) [62], [63].

Seawater collected at 20 m depth from the CTD cast used for pigment analysis was also used for the experimental incubations. Seawater was gently siphoned from the Niskin bottles via silicon tubing affixed with 200 µm mesh screen caps, to exclude large zooplankton, into thirty-two, 2-liter round, clear, acid-cleaned polycarbonate bottles equipped with custom-made gas inflow and sampling ports. The seawater in 16 bottles was bubbled continuously with ambient air using aquarium pumps, and the seawater in another 16 bottles was bubbled continuously with commercial air/pCO2 gas mixtures (custom mixed by Linde/Spectra Gases) at the target “high” CO2 level of 750 ppm. All bottles were placed in a flow-through seawater tank on the ship’s deck, maintaining a temperature within±1°C of initial ambient sea surface temperature (0.65°C). One layer of neutral density screening was used to reduce light to ∼50% of surface irradiance. The seawater was bubbled with target pCO2 concentrations for 24 hours (Table 1). During this equilibration period, Antarctic krill were collected via near-surface net tows (2-m square-frame net with 700 µm mesh and a non-filtering cod end) near the same location as the seawater collection. Tows were sorted, and 40 adult Euphausia superba of similar size (20 non-gravid, and 20 gravid females) were gently placed into a 30-liter tub filled with filtered seawater. The non-gravid krill were selected solely based on size. Their sex was not determined, so both males and females were likely included in the treatments. The selected krill were allowed to empty their guts (∼12 hours) prior to the experimental incubations. After the 24-hour equilibration, 3 seawater bottles from each the ambient (390 ppm) and high (750 ppm) CO2 treatments were sacrificed for a suite of replicate measurements (see below), and served as the initial, or T0, time point. Additionally, 10 non-gravid and 10 gravid presorted krill were sacrificed for initial sampling of metabolic enzyme activity and biochemical composition (see below). From the remaining 13 seawater bottles for each CO2 treatment, 1 adult non-gravid krill was added to each of 5 bottles (1 krill/bottle), 1 adult gravid female was added to each of 5 bottles (1 krill/bottle), and 3 bottles remained without krill and served as controls. The bottles were incubated and continuously bubbled with target CO2 concentrations at a flow rate of 1 ml min−1 (∼1 bubble sec−1) as described above for an additional 24 hours. Another suite of seawater samples and the remaining krill from the incubation bottles were collected at this final time point, Tf.

[Figure omitted. See PDF.]

Table 1. Carbonate chemistry.

https://doi.org/10.1371/journal.pone.0052224.t001

Seawater Sample Collection and Analysis

The seawater samples collected at T0 and Tf include carbonate chemistry parameters (salinity, dissolved inorganic carbon [DIC], total alkalinity [AT], dissolved reactive silicate, phosphate [PO43−]); bacterial abundance; phytoplankton and microzooplankton abundance, size, and identification; chlorophyll, particulate carbon and nitrogen, dissolved organic carbon (DOC), urea, ammonium (NH4+), nitrate (NO3−), and total dissolved nitrogen (TDN).

Salinity was determined with a Guideline AutoSal salinometer from samples collected from each incubation bottle. BOD bottles (300 ml) for DIC/ AT analysis were allowed to overflow at least one to two volumes and then filled via siphoning directly from each of the incubation bottles. The samples were each fixed with 200 µl of saturated mercuric chloride, sealed with a pre-greased glass stopper followed by tape, and stored in a cool, dark location until analysis. DIC was determined using a coulometer to measure the amount of released CO2 after the sample is mixed with phosphoric acid [64]. AT was determined with an open-cell, potentiometric titration of seawater with 0.1 M HCl following the potential of a pH electrode [64]. Peak-area measurements from the DIC and AT analyses were calibrated using certified reference materials (CRMs) obtained from Andrew Dickson at UCSD Scripps Institute of Oceanography. Additionally, an aliquot of seawater from each incubation bottle was siphoned into a 60 ml syringe, filtered through a 0.8 µm syringe filter into a 50 ml polypropylene centrifuge tube, and frozen until later spectrophotometric determination of dissolved reactive silicate [65]. Final carbonate system parameters and pH were calculated using CO2 calc software [66] using a total pH scale (mol/kg-SW), K1 and K2 constants [67] with refits [68], and the acidity constant of KHSO4 in seawater [69].

Whole water samples for algal and protozoan cell counts were collected in duplicate from each experimental bottle, one preserved with acid Lugol’s solution (final concentration 5%) and one preserved with 37% buffered formaldehyde. Subsamples for algal identification of major taxa (cryptophytes and large [>20 µm] diatoms) and cell counts were settled in 1 ml Sedgewick rafters, and three replicate frames of at least 100 cells were counted with a Nikon DIAPHOT-TMD inverted microscope at 600X magnification. Subsamples (50–100 ml) for protozoans (ciliates, heterotrophic dinoflagellates) were settled for at least 24 hours in Utermöhl settling chambers, after which the entire contents were counted under an inverted microscope [70], [71]. An aliquot of seawater from each incubation bottle was also filtered onto a GF/F filter, which was wrapped in foil and frozen for fluorometric chlorophyll a (chl a) analysis [72]. Clearance and ingestion rates of Euphausia superba on prey were calculated according to the equations of Frost [73]. Chl a ingestion rates were converted to carbon (C) using the C:Chl ratio of 63∶1 determined previously in the study region [74] using the methods of Ducklow et al. [75]. Cell volumes of ciliates were calculated according to geometric shapes with linear measurements made via microscopy (minimum of n = 50 per cell type). Carbon contents for ingestion rate conversions were then estimated using cell C to volume conversions for ciliates from Putt and Stoecker [76].

A known volume of seawater from each incubation bottle was also filtered onto a combusted GF/F filter for particulate carbon and nitrogen analysis using a Carbon-Hydrogen-Nitrogen elemental analyzer (Carlo Erba Instruments, NA 1500 Series 2). The remaining filtrate for each replicate was analyzed for dissolved organic and inorganic nutrient concentrations. Concentration of DOC was measured with a Shimadzu TOC analyzer V (minimum detection limit, MDL = 0.5–1.0 µmol l−1; coefficient of variance, CV = 2–644%) after acidification and purging of dissolved inorganic carbon [77], [78]. Ammonium was measured with the phenol/hypochlorite Koroleff method with MDL = 0.05 µmol l−1 and CV = 2.5% [72], [79] and urea was measured with the diacetyl monoxime procedure with MDL = 0.05 µmol l−1 and CV = 2% [80]. Concentrations of PO43− (Koroleff method) (MDL = 0.05 µmol l−1; CV = 2–3%) were determined with a QuikChem 8500 AutoAnalyzer [79].

Krill nutrient release rates (ng individual−1 hour−1) were calculated according to the following equation:where ΔCt is the change in nutrient concentrations (ng l−1 day−1) in the treatment bottles and ΔCc is the average change in nutrient concentrations (ng l−1 day−1) in the control bottles; V is the incubation volume (l), N is the number of grazers in the treatment bottles, and T is incubation time (24 hours day−1). Nutrient uptake by algae and bacteria likely occurred during the incubations, and this nutrient decline is incorporated in this equation in the controls as ΔCc.

Metabolic Enzyme Activities and Chemical Composition of Krill

Krill collected at the start and end of the perturbation experiment were placed into individual 5 ml cryovials, flash frozen in liquid nitrogen, and stored at −80 °C until analysis for biochemical composition in the laboratory. A subset of replicates were analyzed for a suite of measurements for metabolic enzyme activities and biochemical composition. Two adult non-gravid krill and two adult gravid female krill for each of the two CO2 treatments at T0 and Tf were thawed, measured (total length), and weighed (wet weight) prior to placement in individual scintillation vials and placed in a 55°C drying oven for 72 hours. Individual krill were re-weighed (dry weight) then homogenized and subsampled (n = 2 per krill) for total particulate C and N and particulate organic carbon and nitrogen (POC and PON, respectively) after vapor phase acidification with concentrated hydrochloric acid on a Carbon-Hydrogen-Nitrogen elemental analyzer (Carlo Erba Instruments, NA 1500 Series 2) according to Hedges and Stern [81]. Particulate carbonates were calculated as the difference between total particulate C and POC [81]. The remaining three adult non-gravid and three gravid krill for each of the two CO2 treatments at T0 and Tf were homogenized and split for analysis of the following: malate dehydrogenase (MDH) and lactate dehydrogenase (LDH) activities [82] and protein content [83].

Statistical Analysis

Statistical comparisons of the effects of CO2 on feeding rates, nutrient release rates, metabolic enzyme activity, and biochemical composition were made by 1-way ANOVA, employing the p = 0.05 level of significance.

Data Management

Data reported here are available on the Palmer, Antarctica Long Term Ecological Research (PAL LTER) DataZoo website: http://pal.lternet.edu/data/.

Results

Carbonate Chemistry

At the start of the incubations with krill, T0, the pH and pCO2 of seawater in ambient bottles averaged 8.12 and 325 ppm, respectively, and the seawater in high CO2 bottles had an average pH and pCO2 of 7.91 and 554 ppm, respectively (Table 1). Tf seawater pH and pCO2 in ambient treatment bottles with krill averaged 7.99 and 452 ppm, respectively, while seawater pH and pCO2 in Tf high CO2 treatment bottles with krill averaged 7.84 and 672 ppm, respectively (Table 1). The differences in pH and pCO2 between ambient and high CO2 treatments were significantly different at T0 and at Tf (one-way ANOVA, p<0.05). Total alkalinity (AT) at Tf was on average 27 µmol kg−1 higher in the seawater treatments with krill added compared to seawater controls with no krill. Additionally, while seawater dissolved inorganic carbon (DIC) in the high CO2 treatment was similar between gravid and non-gravid krill, seawater AT was significantly higher (36 µmol kg−1) in the gravid krill treatment, which resulted in significantly lower calculated pCO2 ( = 563 ppm in gravid,  = 779 ppm in non-gravid; p<0.01) and caused high variability in AT and pCO2 when gravid and non-gravid krill were averaged together (Table 1).

Composition of Phytoplankton and Microzooplankton

The upper water column at the location where we collected seawater and krill for the experiment was well-mixed and uniform. In situ chl a biomass in the top 50 m was uniformly low (2.6 µg L−1) and then dropped to near undetectable levels below 50 m. Cryptophytes contributed on average 57% (± 1.2% SD) to total phytoplankton chl a and thus were the dominant phytoplankton group in our study area (Table 2). Contributions to total phytoplankton chl a by other taxa were much lower (diatoms, 35%; mixed flagellates, 4%; prasinophytes, 2%; and type 4 haptophytes, 1%; Table 2). Abundances of ciliates in our experimental bottles were very low (∼2–3 cells ml−1; Table 2), and heterotrophic dinoflagellates were not detected in our 100 ml seawater subsamples.

[Figure omitted. See PDF.]

Table 2. Initial composition and abundance of phytoplankton and microzooplankton.

https://doi.org/10.1371/journal.pone.0052224.t002

Feeding Rates

Chlorophyll carbon ingestion rates of krill exposed to high CO2 averaged 78 µg C individual−1 d−1 (0.05 % body C d−1) and were 3.5 times higher than krill ingestion rates at ambient CO2 (p<0.05 for all and gravid krill; Fig. 1). Ingestion rates of ciliates by krill were much lower compared to chlorophyll ingestion rates, averaging 4.7 and 6.3 µg C individual−1 d−1 in the ambient and high CO2 treatments for all krill, respectively (Fig. 1).

[Figure omitted. See PDF.]

Figure 1. Euphausia superba ingestion rates.

Chlorophyll a and ciliate carbon ingestion rates of krill exposed to ambient (Amb) and high ( = 672 ppm) concentrations of CO2. Chlorophyll a was converted to C using a C:chl ratio of 63∶1 estimated in our study area by Bernard et al. [74]. Ingestion rates of krill on ciliates were converted to C using cell C to volume conversions from Putt and Stoecker [76]. Mean of n = 5 for non-gravid and gravid krill and n = 10 for all krill; error bars = 2×standard error. x and y denote a significant difference in total C ingestion rates between treatments (x>y, p<0.05).

https://doi.org/10.1371/journal.pone.0052224.g001

Nutrient Release Rates

Release rates of dissolved organic carbon (DOC), ammonium (NH4+), and phosphate (PO43−) by krill were up to 3.0, 1.5, and 1.5 times higher, respectively, in the high CO2 treatment compared to ambient CO2 (Fig. 2). DOC, NH4+, and PO43− release averaged 0.2% body C d−1, 1.1% body N d−1, and 1.9% body P d−1 (based on a body N:P molar ratio of 20.4 in adult Euphausia superba, [84]), respectively, for krill in the high CO2 treatment. Release rates of urea, however, were consistently about 17% lower in the high CO2 treatment compared to ambient. This resulted in a lower proportion of urea release (% total measured nitrogen: urea+NH4+) by krill in the high CO2 treatment ( = 29%) compared to ambient ( = 38%).

[Figure omitted. See PDF.]

Figure 2. Euphausia superba organic and inorganic nutrient release rates.

Dissolved organic carbon (DOC; a), ammonium (NH4+; b), phosphate (PO43−; c), and urea (d) release rates of krill exposed to ambient and high ( = 672 ppm) concentrations of CO2. Mean of n = 5 for non-gravid and gravid krill and n = 10 for all krill; error bars = 2×standard error. For DOC release, x and y denote a significant difference between treatments (x>y, p<0.05).

https://doi.org/10.1371/journal.pone.0052224.g002

Krill Chemical Composition and Enzyme Activity

Total, organic, and inorganic particulate carbon (C) contents (%C, %POC, %PIC) of krill exposed to high CO2 were significantly higher than krill incubated in ambient CO2 conditions (p<0.05 for %C and %POC in non-gravid and gravid krill and for %PIC in all and non-gravid krill; Fig. 3). Nitrogen (N) contents of krill, however, were lower in the high CO2 treatment compared to ambient, and this was most pronounced in non-gravid krill (p<0.05 for all and non-gravid krill; Fig. 3). Higher C and lower N drove higher C:N ratios in krill in the high CO2 treatment (p<0.05 for all and non-gravid krill; Fig. 3). Compared to gravid krill, non-gravid krill had lower %C and %POC (p<0.05 for T0, Tf ambient, and Tf high CO2), %PIC (p<0.05 for T0 and Tf ambient), and molar C:N (p<0.05 for T0, Tf ambient, and Tf high CO2). However, non-gravid krill had higher %N (p<0.05 for T0 and Tf ambient) and protein content (Table 3; Fig. 3). The activity of enzymes malate dehydrogenase (MDH) and lactate dehydrogenase (LDH), proxies for overall metabolism and respiration [85]–[87], were consistently higher and protein content was consistently lower in the high CO2 treatment compared to ambient (Fig. 4). Additionally, MDH and LDH activities increased from T0 to Tf in the high CO2 treatment and remained unchanged from T0 to Tf ambient treatment, except for LDH in the Tf ambient NG treatment, which increased from T0. Protein content of krill, however, increased slightly from T0 to Tf in the ambient treatment, but decreased from T0 to Tf in the high CO2 treatment. However, there were no significant differences in MDH, LDH, citrate synthase (CS) activity, or protein content between krill type or CO2 treatments (p>0.05).

[Figure omitted. See PDF.]

Figure 3. Euphausia superba chemical composition.

Chemical composition (a: % carbon, %C; b: % nitrogen, %N; c: C:N molar ratio; d: % particulate organic carbon, % POC; and e: % particulate inorganic carbon, % PIC) of krill exposed to ambient and high ( = 672 ppm) concentrations of CO2. Mean of n = 4 for non-gravid and gravid krill and n = 8 for all krill; error bars = 2×standard error. v, w, x, and y denote significant differences between treatments (v>w>x>y, p<0.05).

https://doi.org/10.1371/journal.pone.0052224.g003

[Figure omitted. See PDF.]

Table 3. Physical and chemical characteristics of Euphausia superba.

https://doi.org/10.1371/journal.pone.0052224.t003

Discussion

Our study is the first to report physiological responses of krill to elevated CO2. Previous studies on organismal response to ocean acidification were limited to calcification rates, growth, and the development of early life stages, and studies focusing on zooplankton are scarce. This study demonstrates that E. superba respond to elevated CO2 by increasing ingestion rate, nutrient release rates, and metabolic activity.

Carbonate Chemistry

Total alkalinity (AT) at the end of the incubation was higher in the seawater treatments with krill added compared to seawater controls with no krill. This was most pronounced in the high CO2 gravid krill treatment, contributing to high variability in AT and pCO2 when calculated for all krill (Tf High+krill; Table 1). Higher AT in treatment bottles with krill is likely the result of increased dissolved organic matter (DOM), as evidenced by an accumulation of DOC in these treatments, either through krill excretion or leaching from their egested fecal pellets. The contribution of DOM to AT has previously assumed to be quantitatively insignificant and thus neglected from algorithms used to calculate AT. However, recent studies demonstrated that the accumulation of DOC (comprised of weak acids/bases) significantly increased AT [88–[90]. The contribution of dissolved organic nitrogen and phosphorus (DON and DOP, respectively) to AT has not yet been evaluated. Additionally, the magnitude of the contribution of DOM is dependent upon the composition of dissolved organic compounds [89]. DOC release rates (magnitude of DOC accumulation in the experimental bottles) in gravid krill were not significantly different from non-gravid krill, yet AT was higher in the gravid krill treatments. Alternate potential explanations for this discrepancy in AT between gravid and non-gravid krill include: differential composition of dissolved organic compounds released, unequal release of DON and DOP, or differences in fecal pellet production rate or chemical composition yielding differences in surrounding seawater chemistry. These parameters were not examined in our study; thus, the reason for the differences in AT between gravid and non-gravid krill remains unknown but worthy of additional study.

Seawater pCO2 in the high CO2 treatment did not reach target levels of 750 ppm after 48 hours of bubbling ( = 581 ppm in Tf High control and  = 672 ppm in Tf High+krill; Table 1), suggesting the seawater in the 2L bottles did not fully equilibrate during the course of the incubation. Nonetheless, pCO2 was significantly higher, and pH and ΩAr were significantly lower, in the high CO2 treatment compared to ambient, and there were significant differential responses in krill feeding, nutrient release, and chemical composition.

Feeding Rates

Ingestion rates of adult E. superba presented here were on the low end of those measured in previous studies. Ingestion rates in our study ranged from 23–99 µg C krill−1 d−1, which is equivalent to 0.01–0.07 µg chl ind−1 h−1, 0.005–0.03 µg chl g ww−1 h−1, and 0.01–0.07% body C d−1. These ingestion rates are within range of those found for E. superba in the fall (0.01 µg chl ind−1 h−1, [91]; near 0 µg chl g ww−1 h−1, [92]) and summer (50–445 µg C krill−1 d−1; [3]), but lower than others report for E. superba in the summer (129–447 µg C krill−1 d−1; [74]). Clearance rates of E. superba in our study ( = 8.1 and 10.3 ml krill−1 h−1 for non-gravid and gravid krill, respectively) were 1–3 orders of magnitude lower than those determined for E. superba by Price et al. ([93]; 100–1400 ml krill−1 h−1). A possible reason for reduced feeding activities of krill in our study is container size, which has been shown to affect clearance rates in krill in a previous study [93]. Clearance rates of E. superba were an order of magnitude higher in 50 L tubs [93] than those in 5 L [93], 4.2 L [19], and 2 L [94], [95] bottles. This response has been attributed to stress, abnormal swimming behavior [93], and insufficient food concentrations causing the “refiltration” of already cleared seawater [96] in small containers. However, our measured clearance rates were 1–2 orders of magnitude lower than previous studies with krill in similarly sized 2 L bottles (134 ml krill−1 h−1, [94]; 217 ml krill−1 h−1, [95]); thus, we suggest reduced feeding activity of krill in the present study was due to low available food concentration and phytoplankton composition, specifically the dominance of small cryptophytes.

The minimum carbon ingestion rate required by E. superba to meet minimal respiratory costs (MCU) and the minimum particulate carbon concentration in seawater needed to meet respiratory costs (Cm) can be estimated from these equations by Holm-Hansen and Huntley [97]: MCU (µg C krill−1 h−1) = 0.452W0.975, where W is krill dry weight in mg; and Cm (mg C L−1) = (MCU*RQ)/(a*F), where RQ is the respiratory quotient (0.8; [93], [98]), a is assimilation efficiency (0.8; [93], [98]), and F is filtration, or clearance rate (ml krill−1 h−1). Using measured krill dry weights ( = 241 mg, non-gravid;  = 477 mg, gravid) and clearance rates ( = 8.1 ml krill−1 h−1, non-gravid;  = 10.3 ml krill−1 h−1, gravid), MCU for non-gravid and gravid krill in our study is 95 and 185 µg C krill−1 h−1, respectively, and Cm for non-gravid and gravid krill in our study is 12 and 18 mg C L−1, respectively. However, the maximum C ingestion rate obtained by krill in our study (gravid, high CO2) was only 4.7 µg C krill−1 h−1, and particulate C concentrations measured at T0 prior to krill being added to the bottles averaged only 0.2 mg C L−1. Both parameters are 1–2 orders of magnitude lower than that required to meet minimal respiratory costs, and thus were essentially starved during the 24-hour incubation. Reduced feeding activity of krill in the present study was most likely due to a combination of low available food concentration and dominance of cryptophytes. Adult E. superba can feed efficiently on microzooplankton (ciliates, heterotrophic dinoflagellates), copepods, and large (>20 µm) phytoplankton (i.e., diatoms) [19], [93], [99]–[101]. Abundances of ciliates in our experimental bottles were very low (∼2–3 cells ml−1; Table 2), heterotrophic dinoflagellates were not detected in our 100 ml seawater subsamples, copepods were removed from the seawater with 200 µm mesh screen prior to the experiment, in situ chl a biomass in the top 50 m where we collected water and krill for the experiment was low (2.6 µg L−1), and cryptophytes contributed ∼57% to total chlorophyll (Table 2). Because of their small size (<10 µm), cryptophytes are not efficiently grazed on by krill [18]–[22]. Additionally, the presence of cryptophytes can inhibit krill grazing, as shown for E. superba in cryptophyte-dominated assemblages [22]. Thus, it is possible that krill were not actively feeding and were starving in situ before we collected them for the incubation.

We also consider potential direct effects of bubbling on feeding processes of krill in our study. No direct comparison of krill feeding rates in aerated vs. non-aerated bottles has been conducted; thus, there is no scientific evidence to support or reject any claims that bubbling caused low feeding rates of krill in our study. Although our estimated grazing rates for krill were low, they were still within range of those estimated in previous studies [3], [91], [92], in which bottles were not aerated during experimental incubations. Additionally, in a study conducted with copepods feeding on phytoplankton [102], aeration had no apparent direct affect on algae consumption. Moreover, Euphausia pacifica krill remained active while feeding in 1-gallon glass jars with gentle aeration for over two months in a study conducted by Yen et al. [103], and E. superba maintained in aerated buckets remained healthy and survived for a longer time compared to krill in non-aerated buckets at the same stock density [104]. Thus, any negative impacts of bubbling on krill feeding rates in our study were likely insignificant compared to impacts of low food availability and the dominance of cryptophytes.

Despite low feeding rates in krill during the incubation, there was a significant CO2−dependent response in krill ingestion rates. Ingestion rates in krill were higher in the high CO2 treatment compared to ambient. We hypothesize that increased feeding at high CO2 reflects the increased energetic cost of maintaining internal acid-base and ionic equilibria. These extra costs of compensation could include a higher demand for acid-base regulator proteins, which was demonstrated in juvenile cephalopods under short-term exposures to elevated CO2 [50]. Acid-base compensation under elevated CO2 may compromise the oxygen transport system in krill [52]. This will require the organism to process more water to extract the oxygen they demand [41], [51], likely increasing muscular activity and difficulty of feeding. The impacts of this response are reflected not only in krill ingestion rates, but also in krill nutrient release rates and chemical composition.

Nutrient Release Rates

Higher rates of DOC, NH4+, and PO4+release in krill in the high CO2 treatments was likely the direct result of higher ingestion rates by krill as shown for zooplankton in previous studies [105], [106]. Release rates of DOC by krill, to our knowledge have only been measured in one other study conducted by Ruiz-Halpern et al. [107]. Their rates of E. superba DOC release (ca. 80–202 µmol C g DW−1 h−1) were higher than those reported here (1.9–2.8 umol C g DW−1 h−1); however, they conducted their experiments immediately after collection (and likely feeding) and their incubations were short term (5–240 minutes compared to our 24 h incubation). DOC and nutrient release rates decline rapidly with time when zooplankton are incubated in filtered seawater (not feeding), and this was evident by the reduction of DOC release rates from 202 µmol C g DW−1 h−1 in a 15 min. incubation to 80 µmol C g DW−1 h−1 in a 240 min. incubation [107]. Ammonium (NH4+) release rates of adult E. superba presented here (4.3–24.3 µg N krill−1 h−1; 0.7–1.5% body N d−1;  = 12.6 and 15.3 µg N krill−1 h−1 for krill in ambient and high CO2 conditions, respectively) were higher compared to those measured in most previous studies (0.3–3.7 µg N krill−1 h−1, [108]; 0.6–1.3 µg N krill−1 h−1, [109];<0.5% body N d−1, [110]), but within range of those determined by Ikeda and Mitchell ([111]; 1.2–10.3 µg N krill−1 h−1, 0.7–1.0% body N d−1). Phosphate release rates of E. superba in our study (1.0–6.4 µg P krill−1 h−1; 1.3–2.4 % body P d−1; 0.006–0.009 µg P mg DW−1 h−1;  = 2.5 and 3.0 µg P krill−1 h−1 for krill in ambient and high CO2 conditions, respectively) were generally higher than or within range of than those measured for E. superba in the summer by Ikeda and Hing Fay ([109]; 0.3–1.6 µg P krill−1 h−1) and Ishii et al. ([84]; 3.2 µg P krill−1 h−1), but lower than those measured for krill feeding on copepods in the winter (0.026 µg P mg DW−1 h−1; [110]). Krill feeding carnivorously tend to have low atomic N:P release ratios (2.02; [110]), whereas starved krill or krill feeding on phytoplankton in the summer have higher N:P release ratios [110], [111]. The atomic N:P release ratio determined in our study averaged 12.9 and 18.4 when calculated for NH4+and NH4++urea, respectively, which is similar to the N:P ratio of 15.3 determined for starved krill in Ikeda and Mitchell [111].

Release rates of urea have not been reported previously for Antarctic krill. The proportion of urea release to total measured nitrogen (20.7–49.5%) were much higher than those found for a tropical euphausiid (1% of total N excreted; [112]). Although direct urea release was not determined in E. superba by Ruiz-Halpern et al. [107], dissolved organic N (DON) release (difference between measured total N and NH4+excretion rates), which would include urea, was about 52% of total measured nitrogen in their study. NH4+is generally the main nitrogenous excretory product of zooplankton [113]. However, urea or DON can be a significant proportion of total N released by zooplankton [114]–[118]. While NH4+release rates were consistently higher in krill in the high CO2 treatment, urea release rates in this treatment were consistently about 17% lower compared to ambient, suggesting potential differences in catabolic processes of krill between treatments. Due to the hypothesized increase in energetic costs under elevated CO2, krill in the high CO2 treatment were not only ingesting more food but were also metabolizing more N-rich protein (reflected in slightly lower protein contents of krill in the high CO2 treatment at Tf; Fig. 4), which could have led to differences in the release of N byproducts.

[Figure omitted. See PDF.]

Figure 4. Euphausia superba metabolic enzyme activity and protein content.

Malate dehydrogenase activity (MDH; a), lactate dehydrogenase activity (LDH; b), and protein (c) of krill exposed to ambient and high ( = 672 ppm) concentrations of CO2. Mean of n = 3 for non-gravid and gravid krill and n = 6 for all krill; error bars = 2 × standard error.

https://doi.org/10.1371/journal.pone.0052224.g004

Krill Chemical Composition and Enzyme Activity

Carbon, nitrogen, protein, LDH, and MDH contents determined for E. superba in our study were comparable to those found in previous studies [84], [87], [109], [119]. Significant differences in krill carbon content (%C, %POC, and %PIC), like nutrient excretion, may have resulted from different ingestion rates whereby higher ingestion rates increased carbon content. The larger size and expanded thoracic cavity of gravid krill, and the presence of eggs, likely resulted in the higher %C, %POC, %PIC (from the carbonate-containing exoskeleton), and C:N compared to non-gravid krill. Additionally, the non-gravid treatment potentially included both male and female krill. Male krill contain fewer lipids and cholesterol compared to female krill [120], and likely have fewer energetic requirements. Relative feeding rates between krill sexes are unknown; however, we expect them to be lower in males compared to females as demonstrated for Acartia tonsa and Centropages hamatus copepods [121]–[123].

Higher %PIC of all krill in the high CO2 treatment compared to ambient as well as the increase of %PIC of non-gravid krill from T0 to Tf in the high CO2 treatment provides evidence that calcification, or the mineralization of exoskeletal material, is not prohibited at elevated CO2 concentration and may in fact be enhanced as previously shown for calcification in brittle stars [59]. Slight, yet consistent, higher activities of MDH and LDH in krill exposed to elevated CO2 resulted from increased energetic demands in this treatment. Similarly, enhanced MDH activities were found in the Mediterranean fish Sparus aurata under low pH conditions [124].

If changes in krill chemical composition were caused only by differences in ingestion rates, then nitrogen content (%N) of krill would also increase with ingestion rate in the high CO2 treatment. However, krill nitrogen contents were lower (driving C:N higher) in the high CO2 treatment. This suggests that the compensation for higher energetic demand for acid-base regulation in krill exposed to elevated CO2 not only increases metabolic activity (significant increases in feeding and nutrient release; slight, yet consistent increases in enzyme activity), but also creates stoichiometric changes within the krill caused by differential partitioning of C and N and the utilization and catabolism of proteins (reducing body N, increasing N excretion).

Conclusions

Our results suggest that Euphausia superba respond to elevated CO2 by increasing ingestion rates, nutrient release rates, and metabolic activity, reflecting enhanced energetic requirements, but at what cost? Increased growth and metabolism was also observed in brittle stars at elevated CO2, but the cost - muscle wastage - was substantial [59]. Stoichiometric changes in krill caused by a decrease in %N and increases in %C, %PIC, and C:N as well as shifts in krill N excretory products at elevated CO2 could be indications of biochemical changes that we were unable to determine during this limited short-term field study. Subtle ocean-acidification induced shifts in physiological processes could affect growth and reproduction and accelerate population declines. Krill compensating for higher energy requirements at high CO2 will increase feeding and nutrient release rates, which, under favorable food conditions, may provide sufficient energy to maintain growth and reproduction. However, in time periods (i.e., winter) or locations (i.e., northern WAP) with lower food availability, and more importantly with continued reductions in phytoplankton biomass as a result of rapid climate change along the WAP [10], [125], krill may not be able to sustain increased energetic costs. Additionally, females may be less tolerant to these future changes compared to male krill due to their relatively higher energetic requirements. Rapid warming in the WAP region will intensify this response, as metabolic rates increase with temperature [126], suggesting that the combined effects of ocean warming and ocean acidification (enhanced energetic costs, decreased oxygen transport) will be detrimental to Antarctic krill. We do not yet know the response of Antarctic krill exposed to chronic elevated CO2 or whether or not krill have the capacity to fully compensate under elevated CO2. Future ocean acidification studies with Antarctic krill should focus on prolonged exposures, which will be necessary to pinpoint the underlying physiological responses to increase CO2, determine potential adaptive strategies of krill to high CO2, and to understand the associated feedbacks on the food web and biogeochemical cycles.

Acknowledgments

We thank Raytheon Polar Services, the captain and crew of the R.V. Laurence M. Gould, and the Palmer Station crew for field assistance. We also thank Hugh Ducklow, Matthew Erickson, Bethan Jones, Lihini Aluwihare, Kate Ruck, Kaycee Coleman, Emily Pirl, Travis Miles, Michael Garzio, and Charlotte Fuller for their assistance with field sampling and processing and analyzing samples. We also acknowledge the other Palmer LTER team members and Debora Iglesias-Rodriguez for experimental support (National Oceanography Center, UK). We thank Sarah Cooley and Brad Seibel for thoughtful discussions regarding carbonate chemistry and krill physiology, respectively. Additionally, we thank Simon Wright for providing the CHEMTAX software and occasional assistance in analysis of HPLC pigments.

Author Contributions

Conceived and designed the experiments: GS OS DS. Performed the experiments: GS DS. Analyzed the data: GS JT EO. Contributed reagents/materials/analysis tools: GS OS JT EO. Wrote the paper: GS. Critical revising of manuscript for important intellectual content: GS OS JT EO DS. Final approval of the version submitted to PLOS ONE: GS OS JT EO DS.

Citation: Saba GK, Schofield O, Torres JJ, Ombres EH, Steinberg DK (2012) Increased Feeding and Nutrient Excretion of Adult Antarctic Krill, Euphausia superba, Exposed to Enhanced Carbon Dioxide (CO2). PLoS ONE 7(12): e52224. https://doi.org/10.1371/journal.pone.0052224

References

1. Everson I (2000) Krill Biology, Ecology and Fisheries. Oxford: Blackwell Science. 372 p.

2. Atkinson A, Siegel V, Pakhomov E, Rothery P (2004) Long-term decline in krill stock and increase in salps within the Southern Ocean. Nature 432: 100–103.

3. Perissinotto R, Pakhomov EA, McQuaid CD, Froneman PW (1997) In situ grazing rates and daily ration of Antarctic krill Euphausia superba feeding on phytoplankton at the Antarctic Polar Front and the Marginal Ice Zone. Mar Ecol Prog Ser 160: 77–91.

4. Froneman PW, Pakhomov EA, Perissinotto R, McQuaid CD (2000) Zooplankton structure and grazing in the Atlantic sector of the Southern Ocean in late austral summer 1993. Part 2. Biochemical zonation. Deep Sea Res I 47: 1687–1702.

5. Pakhomov EA, Froneman PW (2004) Zooplankton dynamics in the eastern Atlantic sector of the Southern Ocean during the austral summer 1997/1998 – part 2: grazing impact. Deep Sea Res II 51: 2617–2631.

6. Laws RM (1985) The ecology of the Southern Ocean. Am Sci 73: 26–40.

7. Smetacek V, Assmy P, Henjes J (2004) The role of grazing in structuring Southern Ocean pelagic ecosystems and biogeochemical cycles. Antarct Sci 16: 541–558.

8. Tovar-Sanchez A, Duarte CM, Hernández-León S, Sañudo-Wilhelmy SA (2007) Krill as a central node for iron cycling in the Southern Ocean. Geophys Res Lett 34: L11601.

9. Whitehouse MJ, Atkinson A, Rees AP (2011) Close coupling between ammonium uptake by phytoplankton and excretion by Antarctic krill, Euphausia superba. Deep Sea Res I 58: 725–732.

10. Montes-Hugo M, Doney SC, Ducklow HW, Fraser W, Martinson D, et al. (2009) Recent Changes in Phytoplankton Communities Associated with Rapid Regional Climate Change Along the Western Antarctic Peninsula. Science 323: 1470–1473.

11. Skvarca P, Rack W, Rott H (1999) 34 year satellite time series to monitor characteristics, extent and dynamics of Larsen B, Antarctic Peninsula. Ann Glaciol 29: 255–260.

12. Vaughan DG, Marshall GJ, Connolley WM, Parkinson C, Mulvaney R, et al. (2003) Recent rapid regional climate warming on the Antarctic Peninsula. Climatic Change 60: 243–274.

13. Martinson DG, Stammerjohn SE, Smith RC, Iannuzzi RA (2008) Palmer, Antarctica, Long-Term Ecological Research program first 12 years: physical oceanography, spatio-temporal variability. Deep Sea Res II 55: 1964–1987.

14. Cook AJ, Fox AJ, Vaughan DG, Ferrigno JG (2005) Retreating Glacier Fronts on the Antarctic Peninsula over the Past Half-Century. Science 308: 541–544.

15. Stammerjohn SE, Martinson DG, Smith RC, Iannuzzi RA (2008) Sea Ice in the Western Antarctic Peninsula Region: Spatio-Temporal Variability from Ecological and Climate Change Perspectives. Deep Sea Res II 55: 2041–2058.

16. Hosie GW, Schultz MB, Kitchener JA, Cochran TG, Richards K (2000) Macrozooplankton community structure off East Antarctica (80–150°E) during the Austral summer of 1995/1996. Deep-Sea Res II 47: 2437–2463.

17. Moline MA, Claustre H, Frazer TK, Schofield O, Vernet M (2004) Alteration of the food web along the Antarctic Peninsula in response to a regional warming trend. Glob Change Biol 10: 1973–1980.

18. Meyer MA, El-Sayed SZ (1983) Grazing of Euphausia superba Dana on natural phytoplankton populations. Polar Biol 1: 193–197.

19. Boyd CM, Heyraud M, Boyd CN (1984) Feeding of the Antarctic krill, Euphausia superba. Crust Biol 4: 1289141.

20. Quetin LB, Ross RM (1985) Feeding by Antarctic krill, Euphausia superba: does size matter? In: Siegfried WR, Condy PR, Laws RM, editors. Antarctic Nutrient Cycles and Food Webs. Berlin: Springer-Verlag. 372–377.

21. Weber LH, El-Sayed SZ (1985) Spatial variability of phytoplankton and the distribution and abundance of krill in the Indian sector of the Southern Ocean. In: Siegfried WR, Condy PR, Laws RM, editors. Antarctic Nutrient Cycles and Food Webs. Berlin: Springer-Verlag. 283–293.

22. Haberman KL, Ross RM, Quetin LB (2003) Diet of Antarctic krill (Euphausia superba Dana): II. Selective grazing in mixed phytoplankton assemblages. J Exp Mar Biol Ecol 283: 97–113.

23. Orr JC, Fabry VJ, Aumont O, Bopp L, Doney SC, et al. (2005) Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms. Nature 437: 681–686.

24. McNeil BI, Matear RJ (2008) Southern Ocean acidification: A tipping point at 450–ppm atmospheric CO2. Proc Natl Acad Sci USA 105: 18860–18864.

25. Solomon S, Qin D, Manning M, Marquis M, Averyt K, et al.. (2007) Climate Change 2007: The Physical Science Basis: Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. New York: Cambridge University Press. 996 p.

26. Earth System Research Laboratory, National Oceanic & Atmospheric Administration. Available: http://www.esrl.noaa.gov/gmd/ccgg/trends/. Accessed 2012 Jun 16.

27. Houghton JT, Fiho LG, Meiro LG, Gallander BA, Harris N, et al.. (1996) Climate change: the science of climate change. Cambridge: Cambridge University Press. 584 p.

28. Sabine CL, Feely RA, Gruber N, Key RM, Lee K, et al. (2004) The oceanic sink for anthropogenic CO2. Science 305: 367–371.

29. Sabine CL, Feely RA (2007) The oceanic sink for carbon dioxide. In: Reay D, Hewitt N, Grace J, Smith K, editors. Greenhouse Gas Sinks. Oxfordshire: CABI Publishing. 31–49.

30. Brewer PG (2009) A changing ocean seen with clarity. Proc Natl Acad Sci USA 106: 12213–12214.

31. Ridgwell A, Schmidt DN (2010) Past constraints on the vulnerability of marine calcifiers to massive carbon dioxide release. Nat Geosci 3: 196–200.

32. Riebesell U (2004) Effects of CO2 enrichment on marine phytoplankton. J Oceanogr 60: 719–729.

33. Raven J, Caldeira K, Elderfield H, Hoegh-Guldberg O, Liss P, et al.. (2005) Ocean acidification due to increasing atmospheric carbon dioxide. Policy Document 12/05. London: The Royal Society. 68 p.

34. Fabry VJ, Seibel BA, Feely RA, Orr JC (2008) Impacts of ocean acidification on marine fauna and ecosystem processes. ICES J Mar Sci 65: 414–432.

35. Doney SC, Fabry VJ, Richard A, Feely RA, Kleypas JA (2009) Ocean acidification: the other CO2 problem. Ann Rev Mar Sci 1: 169–192.

36. Siegel V (1988) A concept of seasonal variation of krill (Euphausia superba) distribution and abundance west of the Antarctic Peninsula. In: Sahrhage D, editor. Antarctic Ocean and Resources Variability. Berlin: Springer-Verlag. 219–230.

37. Nicol S (2006) Krill, currents, and sea ice: Euphausia superba and its changing environment. BioSci 56: 111–120.

38. Godlewska M (1996) Vertical migrations of krill (Euphausia superba Dana). Pol. Arch. Hydrobiol. 14: 9–63.

39. Lascara CM, Hofmann EE, Ross RM, Quetin LB (1999) Seasonal variability in the distribution of Antarctic krill, Euphausia superba, west of the Antarctic Peninsula. Deep-Sea Res. I 46: 951–984.

40. Kawaguchi S, Kurihara H, King R, Hale L, Berli T, et al.. (2010) Will krill fare well under Southern Ocean acidification? Biol Lett. doi:10.1098/rsbl.2010.0777.

41. Pörtner H-O, Langenbuch M, Reipschläger A (2004) Biological impact of elevated ocean CO2 concentration: lessons from animal physiology and earth history. J Oceanogr 60: 705–718.

42. Pörtner H-O, Langenbuch M, Michaelidis B (2005) Synergistic effects of temperature extremes, hypoxia, and increases in CO2 on marine animals: From Earth history to global change. J Geophys Res 110: C09S10.

43. Gutowska MA, Pörtner H-O, Melzner F (2008) Growth and calcification in the cephalopod Sepia officinalis under elevated seawater pCO2. Mar Ecol Prog Ser 373: 303–309.

44. Somero GN (1985) Intracellular pH, buffering substances and proteins: imidazole protonation and the conservation of protein structure and function. In: Gilles R, Gilles-Baillien M, editors. Transport Processes, Iono- and Osmoregulation. Berlin: Springer-Verlag. 454–468.

45. Seibel BA, Walsh PJ (2003) Biological impacts of deep-sea carbon dioxide injection inferred from indices of physiological performance. J Exp Biol 206: 641–650.

46. Rosa R, Seibel BA (2008) Synergistic effect of climate-related variables suggests future physiological impairment in a top oceanic predator. Proc Natl Acad Sci USA 52: 20776–20780.

47. Maas AE, Elder LE, Dierssen HM, Seibel BA (2011) The metabolic response of Antarctic pteropods (Mollusca: Gastropoda) to regional productivity: implications for biogeochemical cycles. Mar Ecol Prog Ser 441: 129–131.

48. Seibel BA, Maas AE, Dierssen HM (2012) Energetic plasticity underlies a variable response to ocean acidification in the pteropod, Limacina helicina antarctica. PLoS One 7(4): e30464.

49. Seibel BA, Walsh PJ (2001) Potential impacts of CO2 injection on deep sea biota. Science 294: 319–320.

50. Hu MY, Tseng Y–C, Stumpp M, Gutowska MA, Kiko R, et al. (2011) Elevated seawater pCO2 differentially affects branchial acid-base transporters over the course of development in the cephalopod Sepa officinalis. Am J Physiol Regul Integr Comp Physiol 300(5): R1100–R1114.

51. Pörtner H-O (1990) An analysis of the effects of pH on oxygen binding by squid (Illex illecebrosus, Loligo pealei) haemocyanin. J Exp Biol 150: 407–424.

52. Bridges C, Savel A, Stöcker W, Markl J, Linzen B (1983) Structure and function of krill (Euphausia superba) haemocyanin–adaptation to life at low temperature. Life Chem Rep Suppl 1: 353–356.

53. Watanabe Y, Yamaguchi A, Ishida H, Harimoto T, Suzuki S, et al. (2006) Lethality of increasing CO2 levels on deep-sea copepods in the western North Pacific. J Oceanogr 62: 185–196.

54. Bechmann RK, Taban IC, Westerlund S, Godal BF, Arnberg M, et al. (2011) Effects of ocean acidification on early life stages of shrimp (Pandalus borealis) and mussel (Mytilus edulis).. J Toxicol Env Health A 74: 424–438.

55. Hauton C, Tyrrell T, Williams J (2009) The subtle effects of sea water acidification on the amphipod Gammarus locusta. Biogeosci 6(8): 1479–1489.

56. Burleson ML, Smatresk NJ (2000) Branchial chemoreceptors mediate ventilatory responses to hypercapnic acidosis in channel catfish. Comp Biochem Phys A Mol Integr Physiol 125: 403–414.

57. McKendry JE, Milsom WK, Perry SF (2001) Branchial CO2 receptors and cardiorespiratory adjustments during hypercarbia in Pacific spiny dogfish (Squalus acanthias). J Exp Biol 204: 1519–1527.

58. McKenzie DJ, Taylor EW, Dalla Valle AZ, Steffensen JF (2002) Tolerance of acute hypercapnic acidosis by the European eel (Anguilla anguilla). J Comp Physiol 172: 339–346.

59. Wood HL, Spicer JI, Widdicombe S (2008) Ocean acidification may increase calcification rates, but at a cost. Proc R Soc B 275: 1767–1773.

60. Barry JP, Tyrrell T, Hansson L, Plattner G-K, Gattuso J-P (2011) Atmospheric CO2 targets for ocean acidification perturbation experiments. In: Riebesell U, Fabry VJ, Hansson L, Gattuso J-P, editors. Guide to best practices for ocean acidification research and data reporting. Luxembourg: Publications Office of the European Union. 53–66.

61. Andersson AJ, Mackenzie FT (2012) Revisiting four scientific debates in ocean acidification research. Biogeosci 9: 893–905.

62. Wright SW, Jeffrey SW, Mantoura RFC, Llewellyn CA, Bjornland T, et al. (1991) Improved HPLC method for the analysis of chlorophylls and carotenoids from marine phytoplankton. Mar Ecol Prog Ser 77: 183–196.

63. Mackey MD, Mackey DJ, Higgins HW, Wright SW (1996) CHEMTAX- A program for estimating class abundances from chemical markers: application to HPLC measurements of phytoplankton pigments. Mar Ecol Prog Ser 144: 265–283.

64. Dickson AG, Sabine CL, Christian JR (2007) Guide to best practices for ocean CO2 measurements. PICES Special Publication 3: 1–191.

65. Strickland JDH, Parsons TR (1968) A Practical Handbook for Seawater Analysis. Bulletin 167. (Fisheries Research Board of Canada). 311 p.

66. Robbins LL, Hansen ME, Kleypas JA, Meylan SC (2010) CO2calc–A user-friendly seawater carbon calculator for Windows, Max OS X, and iOS (iPhone): U.S. Geological Survey Open-File Report 2010–1280. 17 p.

67. Mehrbach C, Culberson CH, Hawley JE, Pytkowicz RM (1973) Measurement of the apparent dissociation constants of carbonic acid in seawater at atmospheric pressure. Limnol Oceanogr 18: 897–907.

68. Dickson AG, Millero FJ (1987) A comparison of the equilibrium constants for the dissociation of carbonic acid in seawater media: Deep Sea Res A. 34(10): 1733–1743.

69. Dickson AG (1990) Standard potential of the reaction AgCl(s)+.5H2(g) = Ag(s)+HCl(aq) and the standard acidity constant of the ion HSO4– in synthetic sea water from 273.15 to 318.15 K. J Chem Thermodyn. 22(2): 113–127.

70. von Utermöhl H (1931) Neue Wege in der quantitativen Erfassung des Planktons. (mit besondere). Bericksichtigung des Ultraplanktons) Verh Int Verein Theor Angew Limnol 5: 567–595.

71. Hasle GR (1978) The inverted microscope method. In: Sournia A, editor. Phytoplankton Manual. Paris: UNESCO. 88–96.

72. Parsons TR, Maita Y, Lalli CM (1984) A Manual of Chemical and Biological Methods for Seawater Analysis. New York: Pergamon Press. 173 p.

73. Frost BW (1972) Effects of size and concentration of food particles on the feeding behavior of the marine planktonic copepod Calanus pacificus. Limnol Oceanogr 17: 805–815.

74. Bernard KS, Steinberg DK, Schofield OME (2012) Summertime grazing impact of the dominant macrozooplankton off the Western Antarctic Peninsula. Deep Sea Res I 62: 111–122.

75. Ducklow HW, Kirchman DL, Quinby HL, Carlson CA, Dam HG (1993) Stocks and dynamics of bacterioplankton carbon during the spring bloom in the eastern North Atlantic Ocean. Deep Sea Res II 40: 245–263.

76. Putt M, Stoecker DK (1989) An experimentally determined carbon:volume ratio for marine “oligotrichous” ciliates from estuarine and coastal waters. Limnol Oceanogr 34: 1097–1103.

77. Peltzer ET, Fry B, Doering PH, McKenna JH, Norrman B, et al. (1996) A comparison of methods for the measurement of dissolved organic carbon in natural waters. Mar Chem 54: 85–96.

78. Sharp JH, Carlson CA, Peltzer ET, Castle-Ward DM, Savidge KB, et al. (2002) Final dissolved organic carbon broad community intercalibration and preliminary use of DOC reference materials. Mar Chem 77: 239–253.

79. Grasshoff K, Kremling K, Ehrhardt M (1999) Methods of Seawater Analysis, 3rd edition. New York: Wiley. 600 p.

80. Price NM, Harrison PJ (1987) Comparison of methods for the analysis of dissolved urea in seawater. Mar Biol 94: 307–317.

81. Hedges JI, Stern JH (1984) Carbon and nitrogen determinations in carbonate-containing solids. Limnol Oceanogr 29: 20–34.

82. Walsh PJ, Bedolla C, Mommsen TP (1989) Scaling and sex-related differences in toadfish (Opsanus beta) sonic muscle enzyme activities. Bull Mar Sci 45: 68–75.

83. Lowry OH, Rosebrough NJ, Farr AL, Randall RJ (1951) Protein measurement with the Folin phenol reagent. J Biol Chem 193: 265–275.

84. Ishii H, Omori M, Maeda M, Watanabe Y (1987) Metabolic rate and elemental composition of the Antarctic krill, Euphausia superba Dana. Polar Biol 7: 379–382.

85. Lehninger AL (1982) Principles of biochemistry. New York: Worth. 1011 p.

86. Geiger SP, Kawall HG, Torres JJ (2001) The effect of the receding ice edge on the condition of copepods in the northwestern Weddell Sea: results from biochemical assays. Hydrobiol 453/454: 79–90.

87. Donnelly J, Kawall H, Geiger SP, Torres JJ (2004) Metabolism of Antarctic micronektonic crustacea across a summer ice-edge bloom: respiration, composition, and enzymatic activity. Deep Sea Res II 51: 2225–2245.

88. Kim H-C, Lee K, Wonyong C (2006) Contribution of phytoplankton and bacterial cells to the measured alkalinity of seawater. Limnol Oceanogr 51: 331–338.

89. Kim H-C, Lee K (2009) Significant contribution of dissolved org7anic matter to seawater alkalinity. Geophys Res Lett 36: L20603 .

90. Koeve W, Oschlies A (2012) Potential impact of DOM accumulation on fCO2 and carbonate ion computations in ocean acidification experiments. Biogeosci 9: 3787–3798.

91. Quetin LB, Ross RM, Amsler MO (1987) Field ingestion rates of Euphausia superba. EoS 68(50): 109.

92. Quetin LB, Ross RM (1991) Behavioral and physiological characteristics of the Antarctic krill, Euphausia superba. Am Zool 31: 49–63.

93. Price HJ, Boyd KR, Boyd CM (1988) Omnivorous feeding behavior of the Antarctic krill Euphausia superba. Mar Biol 97: 67–77.

94. Kato M, Segawa S, Tanoue E (1982) Filtering and ingestion rates of the Antarctic krill, Euphausia superba Dana. Trans Tokyo Univ Fish 5: 167–175.

95. Morris DJ (1984) Filtration rates of Euphausia superba Dana: under or overestimates? J Crust Biol 4(1): 185–197.

96. Clarke A, Morris DJ (1983) Towards an energy budget for krill: the physiology and biochemistry of Euphausia superba Dana. Polar Biol 2: 69–86.

97. Holm-Hansen O, Huntley M (1984) Feeding requirements of krill in relation to food sources. J Crust Biol 4(1): 156–173.

98. Ross RM (1982) Energetics of Euphausia pacifica I. Effects of body carbon and nitrogen and temperature on measured and predicted production. Mar Biol 68: 1–13.

99. Hopkins TL (1985) Food web of an Antarctic midwater ecosystem. Mar Biol 89: 197–212.

100. Pakhomov EA, Perissinotto R, Froneman PW, McQuaid CD (1997) Energetics and feeding dynamics of Euphausia superba in the South Georgia region during the summer of 1994. J Plankton Res 19: 399–423.

101. Perissinotto R, Gurney L, Pakhomov EA (2000) Contribution of heterotrophic material to diet and energy budget of Antarctic krill, Euphausia superba. Mar Biol 136: 129–135.

102. Jones RH, Flynn KJ, Anderson TR (2002) Effect of food quality on carbon and nitrogen growth efficiency in the copepod Acartia tonsa. Mar Ecol Prog Ser 235: 147–156.

103. Yen J, Brown J, Webster DR (2003) Analysis of the flow field of the krill, Euphausia pacifica. Mar Fresh Behav Physiol 36: 307–319.

104. King R, Nicol S, Cramp P, Swadling KM (2003) Krill maintenance and experimentation at the Australian Antarctic Division. Mar Fresh Behav Physiol 36: 271–283.

105. Corner EDS, Head RN, Kilvington CC, Pennycuick L (1976) On the nutrition and metabolism of zooplankton: X. Quantitative aspects of Calanus helgolandicus feeding as a carnivore. J Mar Biol Assoc UK 56: 345–358.

106. Kiørboe T, Mohlenberg F, Hamburger K (1985) Bioenergetics of the planktonic copepod Acartia tonsa: relation between feeding, egg production and respiration, and composition of specific dynamic action. Mar Ecol Prog Ser 26: 85–97.

107. Ruiz-Halpern S, Duarte CM, Tovar-Sanchez A, Pastor M, Horstkotte B, et al. (2011) Antarctic krill as a major source of dissolved organic carbon to the Antarctic ecosystem. Limnol Oceanogr 56: 521–528.

108. Meyer B, Auerswald L, Siegel V, Spahic S, Pape C, et al. (2010) Seasonal variation in body composition, metabolic activity, feeding, and growth of adult krill Euphausia superba in the Lazarev Sea. Mar Ecol Prog Ser 398: 1–18.

109. Ikeda T, Hing Fay E (1981) Metabolic activity of zooplankton from the Antarctic Ocean. Aust J Mar Fresh Res 32: 921–930.

110. Huntley ME, Nordhausen W, Lopez MDG (1994) Elemental composition, metabolic activity, and growth of Antarctic krill Euphausia superba during winter. Mar Ecol Prog Ser 107: 23–40.

111. Ikeda T, Mitchell AW (1982) Oxygen uptake, ammonium excretion and phosphate excretion by krill and other Antarctic zooplankton in relation to their body size and chemical composition. Mar Biol 71: 283–298.

112. Roger C (1982) Macroplancton et micronekton de l’Atlantique tropical. II. Cycles de l’azote et du phosphore. Remarques sur la mesure de la production. Océanographie tropicale 17(2): 177–185.

113. Regnault M (1987) Nitrogen excretion in marine and freshwater Crustacea. Biol Rev 62: 1–24.

114. Miller CA, Glibert PM (1998) Nitrogen excretion by the calanoid copepod Acartia tonsa: results of mesocosm experiments. J Plankton Res 20: 1767–1780.

115. Miller CA, Roman MR (2008) Effects of food nitrogen content and concentration on the forms of nitrogen excreted by the calanoid copepod, Acartia tonsa. J Exp Mar Biol Ecol 359: 11–17.

116. Steinberg DK, Saba GK (2008) Nitrogen consumption and metabolism in marine zooplankton. In: Capone DG, Bronk DA, Mulholland MR, Carpenter EJ, editors. Nitrogen in the Marine Environment, 2nd Edition. Boston: Academic Press. 1135–1196.

117. Saba GK, Steinberg DK, Bronk DA (2009) Effects of diet on release of dissolved organic and inorganic nutrients by the copepod Acartia tonsa. Mar Ecol Prog Ser 386: 147–161.

118. Saba GK, Steinberg DK, Bronk DA (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 404: 47–56.

119. Cullen M, Kaufmann RS, Lowery MS (2003) Seasonal variation in biochemical indicators of physiological status in Euphausia superba from Port Foster, Deception Island, Antarctica. Deep Sea Res II 50: 1878–1798.

120. Albessard E, Mayzaud P, Cuzin-Roudy J (2001) Variation of lipid classes among organs of the northern krill Meganyctophanes norvegica, with respect to reproduction. Comp Biochem Physiol A Mol Integr Physiol 129: 373–390.

121. Conover RJ (1956) Oceanography of Long Island Sound, 1952–1954. VI. Biology of Acartia clausi and A. tonsa. Bull Bingham Oceanogr Coll 15: 156–233.

122. Mauchline J (1998) The biology of calanoid copepods. In: Blaxter JHS, Southward AJ, Tyler PA, editors. Advances in Marine Biology, Vol. 33. San Diego, CA: Academic Press. 710 p.

123. Saage A, Vadstein O, Sommer U (2009) Feeding behavior of adult Centropages hamatus (Copepoda, Calanoida): functional response and selective feeding experiments. J Sea Res 62: 16–21.

124. Michaelidis B, Spring A, Pörtner H-O (2007) Effects of long-term acclimation to environmental hypercapnia on extracellular acid–base status and metabolic capacity in Mediterranean fish Sparus aurata. Mar Biol 150: 1417–1429.

125. Schofield O, Ducklow HW, Martinson DG, Meredith MP, Moline MA, et al. (2010) How do polar marine ecosystems respond to rapid climate change? Science 328: 1520.

126. Ikeda T (1985) Metabolic rates of epipelagic marine zooplankton as a function of body mass and temperature. Mar Biol 85: 1–11.