Effects of ocean acidification on respiration, growth and C:N ratio of arctic Hyas araneus larvae

Early life stages of marine crustaceans respond sensitively to elevated seawater PCO2. However, the underlying physiological mechanisms have not been studied well. We therefore investigated the effects of elevated seawater PCO2 on oxygen consumption, dry weight, elemental composition, median developmental time (MDT) and mortality in zoea I larvae of the spider crab Hyas araneus (Svalbard 79°N/11°E; collection, May 2009; hatch, December 2009). At the time of moulting, oxygen consumption rate had reached a steady state level under control conditions. In contrast, elevated seawater PCO2 caused the metabolic rate to rise continuously leading to a maximum 1.5-fold increase beyond control level a few days before moulting into the second stage (zoea II), followed by a pronounced decrease. Dry weight of larvae reared under high CO2 conditions was lower than in control larvae at the beginning of the moult cycle, yet this difference had disappeared at the time of moulting. MDT of zoea I varied between 45 ± 1 days under control conditions and 42 ± 2 days under the highest seawater CO2 concentration. The present study indicates that larval development under elevated seawater PCO2 levels results in higher metabolic costs during premoulting events in zoea I. However, H. araneus zoea I larvae seem to be able to compensate for higher metabolic costs as larval MDT and survival was not affected by elevated PCO2 levels.

High-resolution in situ oxygen microprofiles, porewater and solid phase geochemistry from the Crimean shelf (Black Sea) from Maria S. Merian cruise MSM15/1

The outer western Crimean shelf of the Black Sea is a natural laboratory to investigate effects of stable oxic versus varying hypoxic conditions on seafloor biogeochemical processes and benthic community structure. Bottom-water oxygen concentrations ranged from normoxic (175 µmol O2/L) and hypoxic (< 63 µmol O2/L) or even anoxic/sulfidic conditions within a few kilometers' distance. Variations in oxygen concentrations between 160 and 10 µmol/L even occurred within hours close to the chemocline at 134 m water depth. Total oxygen uptake, including diffusive as well as fauna-mediated oxygen consumption, decreased from 15 mmol/m**2/d on average in the oxic zone, to 7 mmol/m**2/d on average in the hypoxic zone, correlating with changes in macrobenthos composition. Benthic diffusive oxygen uptake rates, comprising respiration of microorganisms and small meiofauna, were similar in oxic and hypoxic zones (on average 4.5 mmol/m**2/d), but declined to 1.3 mmol/m**2/d in bottom waters with oxygen concentrations below 20 µmol/L. Measurements and modeling of porewater profiles indicated that reoxidation of reduced compounds played only a minor role in diffusive oxygen uptake under the different oxygen conditions, leaving the major fraction to aerobic degradation of organic carbon. Remineralization efficiency decreased from nearly 100 % in the oxic zone, to 50 % in the oxic-hypoxic zone, to 10 % in the hypoxic-anoxic zone. Overall, the faunal remineralization rate was more important, but also more influenced by fluctuating oxygen concentrations, than microbial and geochemical oxidation processes.

Fatty acid composition, and bacteria and meiofauna abundance in sediment cores PS71/013 and PS71/085 during POLARSTERN cruise ANT-XXIV/2

During the RV Polarstern ANT XXIV-2 cruise to the Southern Ocean and the Weddell Sea in 2007/2008, sediment samples were taken during and after a phytoplankton bloom at 52°S 0°E. The station, located at 2960 m water depth, was sampled for the first time at the beginning of December 2007 and revisited at the end of January 2008. Fresh phytodetritus originating from the phytoplankton bloom first observed in the water column had reached the sea floor by the time of the second visit. Absolute abundances of bacteria and most major meiofauna taxa did not change between the two sampling dates. In the copepods, the second most abundant meiofauna taxon after the nematodes, the enhanced input of organic material did not lead to an observable increase of reproductive effort. However, significantly higher relative abundances of meiofauna could be observed at the sediment surface after the remains of the phytoplankton bloom reached the sea floor. Vertical shifts in meiofauna distribution between December and January may be related to changing pore-water oxygen concentration, total sediment fatty acid content, and pigment profiles measured during our study. Higher oxygen consumption after the phytoplankton bloom may have resulted from an enhanced respiratory activity of the living benthic component, as neither meiofauna nor bacteria reacted with an increase in individual numbers to the food input from the water column. Based on our results, we infer that low temperatures and ecological strategies are the underlying factors for the delayed response of benthic deep-sea copepods, in terms of egg and larval production, to the modified environmental situation.

Seawater carbonate chemistry, thickness and carbonate elemental composition of the test of juvenile sea urchins in a laboratory experiment

Ocean surface CO2 levels are increasing in line with rising atmospheric CO2 and could exceed 900 µatm by year 2100, with extremes above 2000 µatm in some coastal habitats. The imminent increase in ocean pCO2 is predicted to have negative consequences for marine fishes, including reduced aerobic performance, but variability among species could be expected. Understanding interspecific responses to ocean acidification is important for predicting the consequences of ocean acidification on communities and ecosystems. In the present study, the effects of exposure to near-future seawater CO2 (860 µatm) on resting (M O2rest) and maximum (M O2max) oxygen consumption rates were determined for three tropical coral reef fish species interlinked through predator-prey relationships: juvenile Pomacentrus moluccensis and Pomacentrus amboinensis, and one of their predators: adult Pseudochromis fuscus. Contrary to predictions, one of the prey species, P. amboinensis, displayed a 28-39% increase in M O2max after both an acute and four-day exposure to near-future CO2 seawater, while maintaining M O2rest. By contrast, the same treatment had no significant effects on M O2rest or M O2max of the other two species. However, acute exposure of P. amboinensis to 1400 and 2400 µatm CO2 resulted in M O2max returning to control values. Overall, the findings suggest that: (1) the metabolic costs of living in a near-future CO2 seawater environment were insignificant for the species examined at rest; (2) the M O2max response of tropical reef species to near-future CO2 seawater can be dependent on the severity of external hypercapnia; and (3) near-future ocean pCO2 may not be detrimental to aerobic scope of all fish species and it may even augment aerobic scope of some species. The present results also highlight that close phylogenetic relatedness and living in the same environment, does not necessarily imply similar physiological responses to near-future CO2.

Physiology, growth and development of larval krill ''Euphausia superba'' in autumn and winter in the Lazarev Sea, Antarctica

The physiological condition of larval Antarctic krill was investigated during austral autumn 2004 and winter 2006 in the Lazarev Sea, to provide better understanding of a critical period of their life cycle. The condition of larvae was quantified in both seasons by determining their body length (BL), dry mass (DM), elemental- and biochemical composition, as well as stomach content analysis, and rates of metabolism and growth. Overall the larvae in autumn were in better condition under the ice than in open water, and for those under the ice there was a decrease in condition from autumn to winter. Thus growth rates of furcilia larvae in open water in autumn were similar to winter values under the ice (mean 0.008 mm/d), whereas autumn, under ice values were higher: 0.015 mm/d. Equivalent larval stages had up to 30% lower BL and 70% lower DM in winter compared to autumn, with mean oxygen consumption 44% lower (0.54 µl O2 DM/h). However, their ammonium excretion rates doubled (from 0.03-0.06 µg NH4 DM/h) so their mean O:N ratio was 46 in autumn and 15 in winter. Thus differing metabolic substrates were used between autumn and winter, suggesting a flexible overwintering strategy, as suggested for adults. The larvae were eating small copepods (Oithona spp.) and/or protozoans as well as autotrophic food under the ice. However, pelagic Chlorophyll a (Chl a) was a good predictor for growth in both seasons. The physics (current speed/ice topography) probably has a critical part to play in whether larval krill can exploit the food that may be associated with sea ice or be advected away from such suitable feeding habitat.

(Table 2) Stable carbon and oxygen isotope record of benthic foraminifera of eastern Mediterranean Sea sediment cores

The paleoproductivity, paleo-oxygenation, and paleohydrographic configuration of the southeastern Mediterranean during the late Holocene was reconstructed on the basis of the isotopic composition of the epibenthic Heterolepa floridana, shallow-endobenthic Uvigerina mediterranea, and the deeper endobenthic Bulimina inflata from two high-resolution cores GA-112 (470 m) and GA-110 (670 m). The Delta d13C between H. floridana and U. mediterranea reveals four intervals of enhanced productivity, from 3.3-2.6, 2.3-1.9, 1.5-1.1, and 0.8-0.4 kyr BP, coinciding with increased nutrient supply by the Nile River. The entire basin was well aerated, with oxygen consumption varying between 1.0 and 3.5 mL O2/L. Oxygen consumption increases toward present day, probably because of higher accumulation of total organic carbon at 1.7 kyr BP, coinciding with the appearance of the mesotropic benthic species. The hydrographic configuration of the basin has changed during the course of the last 3.75 kyr. The Levantine Intermediate Water (LIW) deepens below 470 m between 3.3 and 2.0 kyr, and especially between 2.5 and 2.0 kyr. During the last 1.5 kyr, the LIW becomes shallower than 470 m, similar to the present day. The change in the hydrographic configuration reflects changes in evaporation/precipitation ratio and in temperature.

Euphausiid respiration model revamped, link to model results

Euphausiids constitute major biomass component in shelf ecosystems and play a fundamental role in the rapid vertical transport of carbon from the ocean surface to the deeper layers during their daily vertical migration (DVM). DVM depth and migration patterns depend on oceanographic conditions with respect to temperature, light and oxygen availability at depth, factors that are highly dependent on season in most marine regions. Changes in the abiotic conditions also shape Euphausiid metabolism including aerobic and anaerobic energy production. Here we introduce a global krill respiration model which includes the effect of latitude (LAT), the day of the year of interest (DoY), and the number of daylight hours on the day of interest (DLh), in addition to the basal variables that determine ectothermal oxygen consumption (temperature, body mass and depth) in the ANN model (Artificial Neural Networks). The newly implemented parameters link space and time in terms of season and photoperiod to krill respiration. The ANN model showed a better fit (r**2=0.780) when DLh and LAT were included, indicating a decrease in respiration with increasing LAT and decreasing DLh. We therefore propose DLh as a potential variable to consider when building physiological models for both hemispheres. We also tested for seasonality the standard respiration rate of the most common species that were investigated until now in a large range of DLh and DoY with Multiple Linear Regression (MLR) or General Additive model (GAM). GAM successfully integrated DLh (r**2= 0.563) and DoY (r**2= 0.572) effects on respiration rates of the Antarctic krill, Euphausia superba, yielding the minimum metabolic activity in mid-June and the maximum at the end of December. Neither the MLR nor the GAM approach worked for the North Pacific krill Euphausia pacifica, and MLR for the North Atlantic krill Meganyctiphanes norvegica remained inconclusive because of insufficient seasonal data coverage. We strongly encourage comparative respiration measurements of worldwide Euphausiid key species at different seasons to improve accuracy in ecosystem modelling.

Benthic oxygen flixes on the Washington shelf and slope

Benthic oxygen fluxes calculated from in situ microelectrode profiles arc compared with benthic flux chamber O2 uptake measurements on a transect of eight stations across the continental shelf and three stations on the slope of Washington State. Station depths ranged from 40 to 630 m and bottom-water oxygen concentrations were 127-38 µM. The fluxes measured by the two methods were similar on the slope, but on the shelf, the chamber flux exceeded the microelectrode flux by as much as a factor of 3-4. We attribute this difference to pore-water irrigation, a process which apparently accounts for the oxidation of a significant amount of organic C in the continental shelf sediments. Combining our diffusive flux data with other data demonstrates clearly that the bottomwater oxygen concentration must play some significant role in determining the sedimentary oxygen consumption rate. Numerical simulation of the microelectrode 0, profiles suggests that roughly half the diffusive 0, flux must be consumed within - 1 mm of the sediment surface. If this conclusion is correct, then the magnitude of the diffusive flux depends both on the bottom-water oxygen concentration and on the supply rate of labile C to the sediment surf'ace.

Part of the global DOC versus AOU (dissolved organic carbon/apparent oxygen utilization) data compilation, WOCE14

Total organic carbon (TOC) samples were collected at 6 stations spaced ~800 km apart in the eastern South Atlantic, from the Equator to 45°S along 9°W. Analyses were performed by high temperature catalytic oxidation (HTCO) in the base laboratory. Despite the complex advection and mixing patterns of North Atlantic and Antarctic waters with extremely different degrees of ventilation, TOC levels below 500 m are quasi-constant at 55±3 µmol C/l, pointing to the refractory nature of deep-water TOC. On the other hand, a TOC excess from 25 to 38 g C/m**2 is observed in the upper 100 m of the permanently stratified nutrient-depleted Equatorial, Subequatorial and Subtropical upper ocean, where vertical turbulent diffusion is largely prevented. Conversely, TOC levels in the nutrient-rich upper layer of the Subantarctic Front only exceeds 9 g C/m**2 the deep-water baseline. As much as 70% of the TOC variability in the upper 500 m is due to simple mixing of reactive TOC formed in the surface layer and refractory TOC in deep ocean waters, with a minor contribution (13%) to oxygen consumption in the prominent subsurface AOU maximum at 200-400 m depth.

Metabolism of mesozooplankton in the North-Eastern Aegean Sea during SES_GR2 in October 2008

The dataset is based on samples taken during October 2008 in the North-Eastern Aegean Sea. NH4 excretion rate: Mesozooplankton is collected by vertical tows within the Black sea water body mass layer in the NE Aegean, using a WP-2 200 µm net equipped with a large non-filtering cod-end (10 l). Macrozooplankton organisms are removed using a 2000 µm net. A few unsorted animals (approximately 100) are placed inside 8 bottles of 350 or 650 ml filled with GF/F or 0.2 µm Nucleopore filtered seawater and then on a wheell at dim light and maintaining the in situ temperature. 4 bottles without animals are used as control. After 24hours bottles are opened and water samples taken for NH4 chemical analysis. Then the bottle content is filtered on pre-combusted preweighted CF/F filters, which are then dried at 60 C and weighted. Calculations are made as described by Ikeda et al. (2000). Samples for the NH4 determination were collected in pre-cleaned 50 ml Duran bottles and analysed onboard immediately after collection. Ammonium concentration was measured on a Perkin Elmer Lambda 25 UV/VIS Spectrometer according to the method of Koroleff (1970). PO4 excretion rate: Mesozooplankton is collected by vertical tows within the Black sea water body mass layer in the NE Aegean, using a WP-2 200 µm net equipped with a large non-filtering cod-end (10 l). Macrozooplankton organisms are removed using a 2000 µm net. A few unsorted animals (approximately 100) are placed inside 8 bottles of 350 or 650 ml filled with GF/F or 0.2 µm Nucleopore filtered seawater and then on a wheell at dim light and maintaining the in situ temperature. 4 bottles without animals are used as control. After 24hours bottles are opened and water samples taken for PO4 chemical analysis. Then the bottle content is filtered on pre-combusted preweighted CF/F filters, which are then dried at 60 C and weighted. Calculations are made as described by Ikeda et al. (2000). Samples for the determination of PO4 were collected in pre-cleaned 50 ml polyethylene volumetric tubes and analysed on board immediately after collection. PO4 concentration was measured on a Perkin Elmer Lambda 25 UV/VIS Spectrometer following the protocol of Murphy and Riley (1962). O2 consumption rate: Mesozooplankton is collected by vertical tows within the Black sea water body mass layer in the NE Aegean, using a WP-2 200 µm net equipped with a large non-filtering cod-end (10 l). Macrozooplankton organisms are removed using a 2000 µm net. A few unsorted animals (approximately 100) are placed inside 8 bottles of 350 or 650 ml filled with GF/F or 0.2 µm Nucleopore filtered seawater and then on a wheell at dim light and maintaining the in situ temperature. 4 bottles without animals are used as control. After 24hours bottles are opened and water samples taken for O2 chemical analysis. Then the bottle content is filtered on pre-combusted preweighted CF/F filters, which are then dried at 60 C and weighted. Calculations are made as described by Ikeda et al. (2000). For the dissolved O2 determination, the samples were fixed immediately after collection and analysed with the Winkler method as modified by Carpenter (1965a and 1965b). Carbon specific CO2 respiration rate: O2 consumption rate was converted to CO2 production using a RQ value of 0.87 (Mayzaud et al. 2005). Conversion of mesozooplankton dry weight to carbon was done using the % of carbon content measured in the same station from the SESAME dataset of zooplankton biomass. Carbon specific NH4 excretion rate: Conversion of mesozooplankton dry weight to carbon was done using the % of carbon content measured in the same station from the SESAME dataset of zooplankton biomass. Carbon specific PO4 excretion rate: Conversion of mesozooplankton dry weight to carbon was done using the % of carbon content measured in the same station from the SESAME dataset of zooplankton biomass.

Part of the global DOC versus AOU (dissolved organic carbon/apparent oxygen utilization) data compilation, Southern Ocean

Recent evidence that dissolved organic carbon (DOC) is a significant component of the organic carbon flux below the photic layer of the ocean (1), together with verification of high respiration rates in the dark ocean (2), suggests that the downward flux of DOC may play a major role in supporting respiration there. Here we show, on the basis of examination of the relation between DOC and apparent oxygen utilization (AOU), that the DOC flux supports ~10% of the respiration in the dark ocean. The contribution of DOC to pelagic respiration below the surface mixed layer can be inferred from the relation between DOC and apparent oxygen utilization (AOU, µM O2), a variable quantifying the cumulative oxygen consumption since a water parcel was last in contact with the atmosphere. However, assessments of DOC/AOU relations have been limited to specific regions of the ocean (3, 4) and have not considered the global ocean. We assembled a large data set (N = 9824) of concurrent DOC and AOU observations collected in cruises conducted throughout the world's oceans (fig. S1, table S1) to examine the relative contribution of DOC to AOU and, therefore, respiration in the dark ocean. AOU increased from an average (±SE) 96.3 ± 2.0 µM at the base of the surface mixed layer (100 m) to 165.5 ± 4.3 µM at the bottom of the main thermocline (1000 m), with a parallel decline in the average DOC from 53.5 ± 0.2 to 43.4 ± 0.3 µM C (Fig. 1). In contrast, there is no significant decline in DOC with increasing depth beyond 1000 m depth (Fig. 1), indicating that DOC exported with overturning circulation plays a minor role in supporting respiration in the ocean interior (5). Assuming a molar respiratory quotient of 0.69, the decline in DOC accounts for 19.6 ± 0.4% of the AOU within the top 1000 m (Fig. 1). This estimate represents, however, an upper limit, because the correlation between DOC and AOU is partly due to mixing of DOC-rich warm surface waters with DOC-poor cold thermocline waters (6). Removal of this effect by regressing DOC against AOU and water temperature indicates that DOC supports only 8.4 ± 0.3% of the respiration in the mesopelagic waters.

Respiration rates and electron transport system activities of copepod species obtained during Maria S. Merian cruise MSM17/3

Respiration rates and electron transport system (ETS) activities were measured in dominant copepod species from the northern Benguela upwelling system in January-February 2011 to assess the accuracy of the ETS assay in predicting in vivo respiration rates. Individual respiration rates varied from 0.06 to 1.60 µL O2/h/ind, while ETS activities converted to oxygen consumption ranged from 0.14 to 4.46 µL O2/h/ind. ETS activities were significantly correlated with respiration rates (r**2 = 0.79, p = 0.0001). R:ETS ratios were lowest in slow-moving Eucalanidae (0.11) and highest in diapausing Calanoides carinatus copepodids CV (0.76) while fast-moving copepods showed intermediate R:ETS (0.23-0.37). 82% of the variance of respiration rates could be explained by differences in dry mass, temperature and the activity level of different copepod species. Three regression equations were derived to calculate respiration rates for diapausing, slow- and fast-moving copepods, respectively, based on parameters such as body mass and temperature. Thus, knowledge about the activity level and behavioral characteristics of copepod species can significantly increase the predictive accuracy of metabolic models, which will help to better understand and quantify the impact of copepods on nutrient and carbon fluxes in marine ecosystems.

Artificial oxygen fluxes measured by the eddy correlation method using stirring-sensitive oxygen microsensor and oxygen optodes in a flume experiment

In the last decade, the aquatic eddy correlation (EC) technique has proven to be a powerful approach for non-invasive measurements of oxygen fluxes across the sediment water interface. Fundamental to the EC approach is the correlation of turbulent velocity and oxygen concentration fluctuations measured with high frequencies in the same sampling volume. Oxygen concentrations are commonly measured with fast responding electrochemical microsensors. However, due to their own oxygen consumption, electrochemical microsensors are sensitive to changes of the diffusive boundary layer surrounding the probe and thus to changes in the ambient flow velocity. The so-called stirring sensitivity of microsensors constitutes an inherent correlation of flow velocity and oxygen sensing and thus an artificial flux which can confound the benthic flux determination. To assess the artificial flux we measured the correlation between the turbulent flow velocity and the signal of oxygen microsensors in a sealed annular flume without any oxygen sinks and sources. Experiments revealed significant correlations, even for sensors designed to have low stirring sensitivities of ~0.7%. The artificial fluxes depended on ambient flow conditions and, counter intuitively, increased at higher velocities because of the nonlinear contribution of turbulent velocity fluctuations. The measured artificial fluxes ranged from 2 - 70 mmol m**-2 d**-1 for weak and very strong turbulent flow, respectively. Further, the stirring sensitivity depended on the sensor orientation towards the flow. Optical microsensors (optodes) that should not exhibit a stirring sensitivity were tested in parallel and did not show any significant correlation between O2 signals and turbulent flow. In conclusion, EC data obtained with electrochemical sensors can be affected by artificial flux and we recommend using optical microsensors in future EC-studies. Flume experiments were conducted in February 2013 at the Institute for Environmental Sciences, University of Koblenz-Landau Landau. Experiments were performed in a closed oval-shaped acrylic glass flume with cross-sectional width of 4 cm and height of 10 cm and total length of 54 cm. The fluid flow was induced by a propeller driven by a motor and mean flow velocities of up to 20 cm s-1 were generated by applying voltages between 0 V and 4 V DC. The flume was completely sealed with an acrylic glass cover. Oxygen sensors were inserted through rubber seal fittings and allowed positioning the sensors with inclinations to the main flow direction of ~60°, ~95° and ~135°. A Clark type electrochemical O2 microsensor with a low stirring sensitivity (0.7%) was tested and a fast-responding needle-type O2 optode (PyroScience GmbH, Germany) was used as reference as optodes should not be stirring sensitive. Instantaneous three-dimensional flow velocities were measured at 7.4 Hz using stereoscopic particle image velocimetry (PIV). The velocity at the sensor tip was extracted. The correlation of the fluctuating O2 sensor signals and the fluctuating velocities was quantified with a cross-correlation analysis. A significant cross-correlation is equivalent to a significant artificial flux. For a total of 18 experiments the flow velocity was adjusted between 1.7 and 19.2 cm s**-1, and 3 different orientations of the electrochemical sensor were tested with inclination angles of ~60°, ~95° and ~135° with respect to the main flow direction. In experiments 16-18, wavelike flow was induced, whereas in all other experiments the motor was driven by constant voltages. In 7 experiments, O2 was additionally measured by optodes. Although performed simultaneously with the electrochemical sensor, optode measurements are listed as separate experiments (denoted by the attached 'op' in the filename), because the velocity time series was extracted at the optode tip, located at a different position in the flume.

Lake Baikal amphipods under climate change: Thermal constraints and ecological consequences, links to supplementary material

Lake Baikal, the world's most voluminous freshwater lake, has experienced unprecedented warming during the last decades. A uniquely diverse amphipod fauna inhabits the littoral zone and can serve as a model system to identify the role of thermal tolerance under climate change. This study aimed to identify sublethal thermal constraints in two of the most abundant endemic Baikal amphipods, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, and Gammarus lacustris, a ubiquitous gammarid of the Holarctic. As the latter is only found in some shallow isolated bays of the lake, we further addressed the question whether rising temperatures could promote the widespread invasion of this non-endemic species into the littoral zone. Animals were exposed to gradual temperature increases (4 week, 0.8 °C/d; 24 h, 1 °C/h) starting from the reported annual mean temperature of the Baikal littoral (6 °C). Within the framework of oxygen- and capacity-limited thermal tolerance (OCLTT), we used a nonlinear regression approach to determine the points at which the changing temperature-dependence of relevant physiological processes indicates the onset of limitation. Limitations in ventilation representing the first limits of thermal tolerance (pejus (= "getting worse") temperatures (Tp)) were recorded at 10.6 (95% confidence interval; 9.5, 11.7), 19.1 (17.9, 20.2), and 21.1 (19.8, 22.4) °C in E. verrucosus, E. cyaneus, and G. lacustris, respectively. Field observations revealed that E. verrucosus retreated from the upper littoral to deeper and cooler waters once its Tp was surpassed, identifying Tp as the ecological thermal boundary. Constraints in oxygen consumption at higher than critical temperatures (Tc) led to an exponential increase in mortality in all species. Exposure to short-term warming resulted in higher threshold values, consistent with a time dependence of thermal tolerance. In conclusion, species-specific limits to oxygen supply capacity are likely key in the onset of constraining (beyond pejus) and then life-threatening (beyond critical) conditions. Ecological consequences of these limits are mediated through behavioral plasticity in E. verrucosus. However, similar upper thermal limits in E. cyaneus (endemic, Baikal) and G. lacustris (ubiquitous, Holarctic) indicate that the potential invader G. lacustris would not necessarily benefit from rising temperatures. Secondary effects of increasing temperatures remain to be investigated.