Thursday , November 21 2019
Home / argentina / Northern cod species face spawning habitat losses if global warming exceeds 1.5 ° C

Northern cod species face spawning habitat losses if global warming exceeds 1.5 ° C


The Atlantic Ocean and Arctic Rapid Climate Change poses a threat to some of the world's largest fish populations. The impact of warming and acidification may become accessible through a mechanism-based risk assessment and projections of future habitat suitability. We show that ocean acidification causes a narrowing of embryonic thermal ranges, which identifies the suitability of habitats as a critical life-history bottleneck for two abundant cod species. Embryonic tolerance ranges linked to climate simulations reveal ever-increasing CO2 emissions [Representative Concentration Pathway (RCP) 8.5] will deteriorate suitability of present spawning habitat for both Atlantic cod (Gadus morhua) and Polar cod (Boreogadus saida) by 2100. Moderate warming (RCP4.5) may avert the dangerous climate on the Atlantic Polar cod, which also loses the benefits of ice-covered ocean. The RCP2.6 emissions following, however, supports large universal habitat suitability for both species, suggesting that it is "below 2 ° C, if not 1.5 ° C," as pledged by the Paris Agreement.


Ocean warming and acidification (OWA), driven by unabated CO2 emissions, are expected to be constrained to the survival and reproduction of many marine organisms (1) Existing knowledge that is a physiological limit for early life stages define the vulnerability of species to OWA (2) Studies of worst-case impact risk management and gain societal acceptance for mitigation policy (3) However, even more important is the identification of emissions that are required to minimize the impact risks and to locate the potential refuge, which should receive priority in conservation (13) Yet, mechanism-based risk assessments that integrate vulnerable life stages and their specific habitat needs into a scenario context are barely available, especially for marine species inhabiting Arctic regions (4, 5)

The Subarctic and Arctic seas around Northern Europe (i.e., Iceland Sea, Norwegian Sea, East Greenland Sea, and Barents Sea) have higher rates of ocean warming, acidification and sea ice loss than most other marine areas on Earth (6) These ocean regions – previously termed the Seas of Norden (7) —Are inhabited by highly productive fish populations, most of the annual migrations to specific spawning locations (4) The biophysical features of suitable spawning habit support the survival of early life as well as their dispersal towards suitable nursery areas (8) Given that fish embryos are often more sensitive to environmental change than subsequent life stages (2), embryonic tolerance may act as a fundamental constraint on spawning habitat suitability. For example, thermal tolerance ranges may represent a biogeographic constraint (8) and are likely to be explained by incomplete development of cardiovascular and other homeostatic systems (9) Ocean acidification (OA) caused by elevated aquatic CO2 levels may exacerbate the disturbance of homeostasis (10), narrowing the thermal range (2, 11) and possibly reducing spawning habitat suitability by impairing egg survival.

Both Atlantic cod and Polar codings are key members of the Northern high-latitude fish fauna, but they differ in terms of thermal affinity and spawning preference4, 5) Atlantic cod is a "thermal generalist" that occupies temperate to Arctic waters between .51.5 ° and 20 ° C (12) In contrast, Polar cod is a "thermal specialist," endemic to the High Arctic and is rarely found at temperatures above 3 ° C (13) Because of overlapping temperature ranges of juvenile and adult life stages, both species coexist during their feeding migrations summer (14) During winter and spring, however, there are separate locations with different water temperatures and sea conditions (Fig. 1). Given that Atlantic cod prefers warmer waters (3 ° to 7 ° C) than Polar cod (−1 ° to 2 ° C), the latter species are considered to be particularly vulnerable to climate change (5, 14) For more information, please contact us for the article about reproduction, which is a nursery habitat for larvae and juveniles during spring and summer (5)

Fig. 1 Distribution patterns of Atlantic cod and Polar cod in the Seas of Norden.

(A) Atlantic cod; (B) Polar cod. Population of both reproductive species during winter and spring (Atlantic cod: March to May; Polar cod: December to March) at species-specific locations (ie, spawning habitats, blue-shaded areas) with characteristic temperature and sea conditions (Atlantic cod: 3 ° to 7 ° C, open water; Polar cod: −1 ° to 2 ° C, closed sea-ice cover). Green arrows indicate dispersal driven egg and larval by prevailing surface currents. During summer, the feeding grounds (green-shaded areas) of both species, partly overlap, for example, around Svalbard, which marks the northernmost distribution limit of Atlantic cod. Red symbols denote the origin of animals (spawning adults) used in this study. Distribution maps were redrawn after (4, 13, 33) NEW, the Northeast Water Poly; FJL, Franz-Joseph-Land; NZ, Novaya Zemlya.

The spawning aggregations of Atlantic and Polar Cod — often comprising millions of individuals — are important resources for humans and other marine predators. For instances, the Norwegian Atlantic cod fisheries alone generates annual revenue of US $ 800 million (15), while Polar cod is an essential food item for many marine birds and mammals (5) Estimating changes in spawning habitat suitability for these focal species therefore has high socioecological relevance (4) Functional responses from embryos to OWA incorporated into habitat models may help to identify the effects of varying emission scenarios, including the goal of limiting global warming to 1.5 ° C above preindustrial levels (16)

Here, we assessed the embryonic ranges of thermal tolerance under OA in Atlantic cod and Polar cod. Oxygen consumption rates (MO2) of eyed-stage embryos and larval morphometrics at hatch provide insight into energetic imposed constraints by OWA. Spawning habitat suitability was mapped across the Seas of Norden under different Representative Concentration Pathways (RCPs) by linking egg survival data with climate simulations of the Coupled Intercomparison Project Phase 5 (CMIP5). RCPs assume either "no greenhouse gas mitigation" (RCP8.5), "intermediate mitigation" (RCP4.5), or "strong mitigation" (RCP2.6). The increase in the global mean surface is relatively to the reference period 1850-1900 and is estimated to be the first estimate for the consequences of keeping global warming to "well below 2 ° C, if not 1.5 ° C," stated in the Paris Agreement (16)


Embryonic oxygen consumption (MO2Increased temperature but increased level of warmest temperatures (Atlantic cod: ≥9 ° C; Polar cod: .54.5 ° C; Fig. 2, A and B), which is, in combination with increased mortality under these conditions (Fig. 3), the indicative of severe heat stress. Embryos acclimated to lower temperatures (<9 ° / 4.5 ° C) and elevated Pco2 (partial pressure of CO2) consumed ~ 10% more oxygen compared to those reared under control Pco2. This trend is reversed upon warming, indicating that additional oxygen and associated energy demands under non-critical OA conditions are high temperatures, causing the upper thermal limit of metabolic maintenance to decline. Higher energy requirements under elevated Pco2 may result from increased acid-base regulation, protein turnover, and damage repair (9, 10) The allocation of energy to life – sustaining functions should receive priority over growth (17), as was evidenced by CO2– and warming-induced reductions in the size at at hatch (Fig. 2, C to F, and Fig. S2). The relative decrease in larval yolk-free body area due to elevated Pco2 averaged 10% for Atlantic cod (P < 0.001) and 13% for Polar cod (P < 0.001), with the smallest larvae hatching at the warmest temperature (Fig. 2, C and D, and table S1). Reductions in larval body size and dry weight (Fig. 2, E and F, and table S1) are in line with the CO2-induced reallocation of energy away from growth also seen in other fish species (18)

Fig. 2 Effects of elevated Pco2 on temperature-dependent oxygen consumption rates (MO2) and growth of Atlantic cod embryos and Polar cod embryos (right).

(A and B) MO2 was measured in eyed-stage embryos (image). Symbols are means (± SEM depicted as bars, n = 6 or 4). Performance curves (lines) are based on n = 28 data points. Dark and light shadings indicate 90 and 95% Bayesian credible confidence intervals, respectively. (C and D) Larval yolk-free body area at hatch was assessed as an indicator of the somatic growth and resource (yolk) utilization. The box plots overlaid with individual values ​​for show the 25th, 50th, and 75th percentile; whiskers mark 95% confidence intervals. (D) Sufficient sample sizes are not available at 6 ° C because most individuals edited or malformed. (E and F) Offsets between regression lines (with 95% confidence intervals) indicate CO2-related differences in size-weight relationships of newly hatched larvae (image). Individuals were pooled across temperature treatments (E: 0 ° to 12 ° C, F: 0 ° to 3 ° C). (A to F) Significant main effects of temperature, Pco2, or their interaction (T * Pco2) are indicated by black ★, whereas orange ★ denote significant CO2 effects within temperature treatments (Tukey post hoc test, n = 6 or 4 per treatment). See table S1 for details on statistical tests. N.a., not available.

Fig. 3 Effects of elevated Pco2 on temperature-dependent egg survival in Atlantic cod and Polar cod.

(A) Atlantic cod; (B) Polar cod. Symbols represent means (± SEM depicted as bars, n = 6). Thermal performance curves (TPCs, lines) of each species are based on n = 36 data points. Dark and light shadings indicate 90 and 95% Bayesian credible confidence intervals, respectively. TPCs were extrapolated into subzero temperatures by incorporating freezing tolerance thresholds from the literature (Materials and Methods). Significant main effects of temperature, Pco2, or their interaction (T * Pco2) are indicated by black ★, whereas orange ★ denote significant CO2 effects within temperature treatments (Tukey post hoc test, n = 6 or 4 per treatment). See table S1 for details on statistical tests.

Egg survival (≤0 ° and ≥9 ° C) and Polar cod (≥3 ° C) decreased the preferred preferred spawning temperatures, particularly under the influence of elevated Pco2 (Fig. 3 and table S1). Accordingly, our results confirm that embryonic tolerance ranges represent tight constraints on the thermal spacing of the Atlantic cod and Polar cod. CO2-induced mortality at their optimum temperature was less pronounced for Atlantic cod (6 ° C, Fig. 3A) than for Polar cod (0 ° to 1.5 ° C, Fig. 3B). This observation corresponds to the variation in CO2 the sensitivity of the optimum temperature conditions (OWA effects)18) However, both species experienced a similar CO2-Related decline in egg survival at their respective warmer threshold (−48% at 9 ° C for Atlantic cod and −67% at 3 ° C for Polar cod). Increased thermal sensitivity of embryos under projected Pco2 levels of implies narrowing of their thermal tolerance range and form of the species reproductive niche (2) As a consequence, the spatial extent of thermally suitable spawning habitat for Atlantic cod and Polar may not only shift to higher latitudes but also contract due to OWA.

Compared with contemporary (known) spawning sites of Atlantic cod and Polar cod in the study area (blue areas in Fig. 1; yellow dashed areas in Fig. 4), our baseline simulations (1985–2004) suggest that spawning occurs exclusively within the thermal optimum range of embryo development [>90% potential egg survival (PES), Fig. 4]. However, the area of ​​thermally suitable spawning habitat (PES> 90%) is larger than the area where spawning actually occurs. For example, despite suitable temperatures, no spawning of Atlantic cod is currently observed in the northeastern Barents Sea (19), indicating that spawning habitat suitability also depends on factors other than temperature. Mechanisms that preclude certain areas suitable for spawning may include a dispersal of eggs and larvae, unfavorable feeding conditions, and predation pressure (8, 19)

Fig. 4 Current (baseline) spawning habitat suitability for Atlantic cod and Polar cod in the Seas of Norden.

(A) Atlantic cod; (B) Polar cod. Spawning habitat suitability is expressed as PES (% PES, color coded) by combining experimental survival data (Fig. 3) with WOA13 temperature fields (1 ° × 1 °, upper 50 m of shelf seas) for the baseline period 1984–2005. Values ​​are averaged over spawning seasons (Atlantic cod: March to May; Polar cod: December to March) and referenced against locations where spawning has been documented[yellowdashedareas([yellowdashedareas(13, 33)]. The spatial extent of thermally suitable spawning habitat (PES> 90%) is typically larger than the "realized spawning habitat" because other limiting factors are not considered. Dotted magenta lines indicate the seasonal sea-ice edge positions (70%), which vary slightly between species due to varying species-specific spawning seasons.

By 2100, unabated OWA (RCP8.5) is a substantial decline in PES at the major sites of both species (Fig. 5, A to C). For Atlantic cod, the PES is around Iceland (−10 to −40%) and the Faroe Islands (−20 to −60%) and along the entire Norwegian coast (−20 to −60%), including the most important spawning sites at the Lofoten archipelago (at 68 ° N, Fig. 5A). In turn, extensive shelf regions off The Poles, across the northeastern Barents Sea will become more suitable (PES, +10 to + 60%) due to decreasing sea-ice cover. However, the habitat potential gains in the cold tolerance of Atlantic cod embryos under OA conditions and, possibly, unknown constraining factors (see above). Under RCP4.5, decrements in the PES of the Atlantic code for some southern spawning sites (eg, Faroe Islands: −10 to −40%) are mostly outweighed by thermal benefits (PES, +20 to + 60%) in the northeastern Barents Sea (between Svalbard, Franz Josef Land, and Novaya Zemlya; Fig. 5, D and F).

Fig. 5 Change in thermally suitable spawning habitat of Atlantic cod (left) and Polar cod (right) in the Seas of Norden under RCPs.

(A to C) RCP8.5: Unabated OWA. (D to F) RCP4.5: Intermediate warming (no acidification considered). (G to I) RCP2.6: Less than 2 ° C global warming (no acidification considered). Maps show the PES shift between the baseline period (1985-2004; spawning season of Atlantic cod: March to May; spawning season of Polar cod: December to March; see Fig. 3) and the median of CMIP5 multimodel-based projections ( seasonal sea surface temperature, 0 to 50 m; see Materials and Methods) for this century's end (2081-2100). Black shading indicates areas (cells, 1 ° × 1 °) with high uncertainty (that is, the shift in PES within that cell is smaller than the CMIP5 ensemble spread; see Materials and Methods). Dotted magenta lines represent sea-ice edge positions of each species-specific spawning season (defined as areas with ice concentrations> 70%). (C, F, and I) For each map, values ​​(change in PES) of individual cells are summarized by kernel density estimations, with the corresponding width to the relative occurrence of values. Box plots show the 25th, 50th, and 75th percentile; the end of the whiskers mark the 95% intervals.

Polar cod will likely experience the most dramatic losses of habitat south of Svalbard and Novaya Zemlya (PES, −40 to −80%; RCP8.5; Fig. 5B). Moreover, Polar cod is going to lose most of the under-ice habit except for a small refuge on the East Greenland shelf (Fig. 5B). Even warming without OA effects (RCP4.5; Fig. 5, E and F) will be substantially reduced from Svalbard (PES, −20 to −60%) and Novaya Zemlya (PES, −10 to −40%). The widespread loss under RCP8.5 and RCP4.5 scenarios may indirectly affect the reproductive success of polar cod, because the ice protects habitat feeding for early-life stages (5) Limiting global warming to about 1.5 ° C above the preindustrial levels (that is, the median temperature of RCP2.6) may not only minimize reductions in PES at the present core to less than 10% (Fig. 5, G) to I) but also maintain some sea-ice cover.


Our project suggest that OWA-driven impacts on egg survival and consequent changes in spatial habitat suitability may be the primary determinant of climate-dependent constraints on Atlantic cod and Polar cod. Present findings are in line with the hypothesis that tolerance thermal ranges and embryo are compressed by progressive OWA (2) Our results also corroborate the idea that there is an existential threat climate change for cold adapted species such as Polar cod (20), although we identified some cold refugia for this species in the High Arctic. Atlantic cod may follow the poleward of the displacement of its thermal optimum, which is currently dominated by Polar cod. The parallel decreases in habitat for the suitability of Iceland and the Norwegian coast (under RCP8.5) implies that, by 2100, spawning south of the Arctic Circle (e.g., south of the Lofoten) may not be more possible for Atlantic cod. The potential displacement of commercial important fish stocks across management boundaries and exclusive economic zones poses a not only to national fishermen and conservationist challenges (5) but also to international bodies and regulations, which intend to avoid overexploitation, resource conflicts, and the degradation of pristine ecosystems in the Arctic (4, 21)

However, if global warming is limited to 1.5 ° C above preindustrial levels, then changes in thermal suitability of present behavior are unlikely to exceed the critical thresholds of Atlantic cod and Polar cod. Residual risks may be reduced further as both species can be adapted to climatic changes, by responding either (i) via shifts in timing and / or location of spawning within present regions (22) or (ii) through transgenerational processes that enhance physiological tolerance (23) Uncertainties in our results also relate to (iii) the reliability and resolution of the CMIP5 climate projections (24)

First, the temporal window for spawning in the winter is limited to late winter extreme seasonality of light and associated primary production at high latitudes (> 60 ° N) (22) Significant changes in the spatial frequency are unlikely to occur in this region. Instead, north expansions of spawning during historic and ongoing warming are well documented, especially for Atlantic cod, which extended its activity to West Svalbard in the 1930s (25) However, core spawning areas (e.g., the Lofoten archipelago for the Barents Sea population) have always been occupied during the past century, possibly due to favorable combinations of biotic and abiotic factors that maximize recruitment success (8, 22) After spawning, dispersal of eggs and larvae to suitable nursery areas – sometimes over hundreds of kilometers – plays an important role in terms of life-cycle connectivity and population replenishment (8) Spawning at alternative locales (as required under RCP8.5 for both species and under RCP4.5 for Polar cod) can affect connectivity and increase the risk of advocacy losses and recruitment failure (8) Accordingly, the successful establishment of new spawning habitats will be largely dependent on factors in addition to egg survival (i.e., availability, predation pressure, and connectivity), all of which are currently difficult to predict (2, 22)

Second, our results assume that embryonic tolerance ranges are constant across different populations and generations (i.e., no evolutionary change within this century). These assumptions are supported by experimental data[egsimilartemperaturetemperatureforeggdevelopmentamongdifferentAtlanticcodpopulations([egsimilartemperatureoptimaforeggdevelopmentamongdifferentAtlanticcodpopulations(26); see also fig. S1], as well as by field observations[egconsistentnorthwardshiftofcodspawningactivityresponsetoprevious/ongoingwarming([egconsistentnorthwardshiftofcodspawningactivityinresponsetoprevious/ongoingwarming(17)]and phylogenetic analysis of thermal tolerance evolution in marine fish[eg<01°Cchangeinthermaltoleranceper1millionyears([eg<01°Cchangeinthermaltoleranceper1millionyears(27)]. Transgenerational plasticity (TGP) may promote short-term adaptation to environmental change via nongenetic inheritance (e.g., maternal transmission) (23) However, in the case of the theory of TGP, the Atlantic cod suggest that egg viability is impaired during similar warming if it is heat exposed during gonad maturation (28) This example of negative TGP corresponds to the majority (57%) of studies on TGP in fish that are observed either neutral (33%) or negative (24%) responses (29) Given the limited capacity for short-term adaptation, it is most likely that species have to abandon their traditional habitats as soon as physiological limits are exceeded (2) Accordingly, our refuge habitats that should receive priority to the implementation of marine reserves.

Third, CMIP5 climate projections include uncertainties (24) To some extent, these uncertainties can be reduced and assessed by multimodel results (see Materials and Methods). Global climate models (24) The model of the impact of climate change projects for the most improved areas of future studies, most elegant by means of multiresolution ocean models with unstructured meshes (30)

In light of embryonic intolerance to OWA, we show that with unabated greenhouse gas emissions, large areas used for spawning will become less suitable for the recruitment of Atlantic cod and Polar cod, possibly leading to cascading impacts on food and web and associated ecosystem services (4, 5) However, our results also emphasize mitigation measures, as pledged by the Paris Agreement, can mitigate climate change effects on both species. Given that current CO2 emissions trajectories yield a 1% chance of limiting global warming to 1.5 ° C above preindustrial levels (31), our results call for immediate emissions following scenarios compatible with warming by 1.5 ° C to convert irreversible ecosystem damage in the Arctic and elsewhere.



Atlantic cod were caught by longlining in the southern Barents Sea (Tromsoflaket: 70 ° 28′00 ″ N, 18 ° 00′00 ″ E) in March 2014. Nofima AS, Tromso, Norway) and held in a flow-through tank (25 m3) under ambient light, salinity [34 practical salinity units (PSU)], and temperature conditions (5 ° ± 0.5 ° C). Polar cod was caught in Kongsfjorden (West Svalbard: 78 ° 95′02 ″ N, 11 ° 99′84 ″ E) by trawling in January 2014. Selected fish were held in flow-through tanks (0.5 m)3) and transferred to the Aquaculture Research Station in Karvikå (NOFIMA, The Arctic University of Norway UiT, Tromso). At the station, the fish were kept in a flow-through tank (2 m3) at 3 ° ± 0.3 ° C water temperature (34 PSU) and complete darkness. In both experiments, the gametes used for in vitro fertilization were obtained by a spawning strip from n = 13 (Polar cod: 12) males and n = 6 females (table S2).

Fertilization protocol

All fertilizations were conducted within 30 minutes after stripping. Each egg batch was split in half and fertilized using ultraviolet (UV) -sterilized seawater (34 PSU), adjusted for the broodstock holding temperature (Atlantic cod: 5 ° C; Polar cod: 3 ° C) and two different Pco2 conditions[control[controlPco2: 400 μatm, pH(Free-Scale) 8.15; high Pco2: 1100 μatm, pHF 7.77]. A standardized dry fertilization protocol with milt aliquots from n = 3 males was used to maximize fertilization success (32)

Fertilization success

Fertilization success was assessed in subsamples (3 × 100 eggs per batch and Pco2 treatment), which were incubated within sealed petri dishes until the 8/16-cell stage (Atlantic cod: 12 hours, 5 ° C; Polar cod: 24 hours, 3 ° C) and photographed under a stereomicroscope for subsequent evaluation (S3 table ) These images were also used to determine the mean egg diameter of an egg batch (30 eggs per batch, table S3).

Incubation setup

According to different spawning seasons, both experiments could be conducted consecutively with the same experimental setup in 2014 (Polar cod: February to April; Atlantic cod: April to May). Eggs previously fertilized at either control or high Pco2 were maintained at the respective CO2 treatment and incubated until different temperatures (Atlantic cod: 0 °, 3 °, 6 °, 9 °, and 12 ° C; Polar cod: 0 °, 1.5 °, 3 °, 4.5 °, and 6 ° C) . Temperature preferences for Atlantic cod (3 ° to 7 ° C) (33) and Polar cod (≤2 ° C) (13) and projected warming scenarios for the respective region. Each treatment group of egg batches were subdivided into two stagnant incubators (20 incubators per female, 120 in each experiment). To not be survival estimates, while using egg survival (while larval morphometrics at hatch), while subsamples required for embryonic MO2 measurements were taken from the second incubator.

Initially, all incubators (volumes, 1000 ml) were filled with filtered (0.2 μm) and UV-sterilized seawater (34 PSU) and were adjusted to be positively buoyant eggs. With regard to oxygen supply in a stagnant incubator, it has a sufficient space to arrange themselves in a single layer beneath the water surface. We therefore adjusted the amount of eggs per incubator (Atlantic cod: ~ 300 to 500; Polar cod: ~ 200 to 300) and Polar cod (~ 1.45 mm) and Polar cod (~ 1.65 mm). Loaded incubators were then located in differently thermostatted shower baths (volume, 400 liters) to ensure smooth temperature change inside the incubator. The transparent, bottom-tapered incubators were sealed with a Styrofoam cover to prevent CO2 outgassing and temperature fluctuations. According to natural light regimes, Atlantic egg received dim light with a daily rhythm of 8 hours light / 16 hours darkness, and Polar cod eggs were kept in darkness except for light exposure during handling. Every 24 hours, 90% of the water volume of each incubator was replaced by filtered (0.2 μm) and UV sterilized to avoid oxygen depletion. An outlet valve was mounted at the bottom of the incubators to drain the seawater with dead eggs, which lose buoyancy and descend to the bottom. Each seawater bath contained two 60-liter reservoir tanks, which were used to adjust the condition to the corresponding temperature and Pco2 conditions. Water temperature inside the water baths are controlled every 15 minutes (± 0.1 ° C) via a multichannel aquarium computer (IKS-Aquastar, IKS Systems, Germany). Future Pco2 conditions were established by injection of pure CO2 submerged gas 60 liters reservoir tanks at each temperature. A multichannel feedback system (IKS-Aquastar), connected to individual pH probes (IKS-Aquastar) and solenoid valves, was used to control water pH and Pco2 values. The Pco2 the water tank with infrared Pco2 probe (Vaisala GMP 343, manual temperature compensation, ± 5 μm accuracy; Vaisala, Finland). The probe was equipped with a MI70 pump and aspiration pump, which was connected to a degassing membrane (G541, Liqui-Cel, 3M, USA) to measure Pco2 in air equilibrated to dissolved water gases (34) Factory calibration was confirmed by measurements of a bubbled with a technical gas mixture (1000 μatm CO2 in air, Air Liquide, Germany). Before the daily water exchange, the reservoir values ​​were measured in decimal places of laboratory electrode pH three (Mettler Toledo InLab Routine Pt 1000 with temperature compensation, Mettler Toledo, Switzerland), which was connected to a WTW 3310 pH meter. A National Bureau of Standards, a two-point calibration with buffers was performed on a daily base. To convert NBS to free proton concentration scale for seawater pH (35), the electrode was calibrated with tris-HCl seawater buffers (36), which was acclimated to the corresponding incubation temperature before each measurement. Seawater pH values ​​refer to the free pH scale (pHF) throughout this manuscript. Seawater parameters are summarized in fig. S3.

Egg survival

Egg mortality was recorded on a 24-hour basis for individuals within an incubator had either edited or hatched (Fig. S4). Once hatching commenced, free-swimming larvae were collected in the morning, euthanized with an overdose of tricaine methanesulfonate (MS-222), and counted after visual examination for morphological deformities under a stereomicroscope. The incidence of valid deformities was quantified as the percentage of exhibiting severe deformations of the yolk sac, cranium, or vertebral column. Egg survival is defined as a percentage of nonmalformed, viable larvae that are initial number of fertilized eggs (Fig. S5). The eggs in an incubator was estimated from the mean fertilization success of the respective egg batch (table S3).


Oxygen consumption rates (MO2) of the eyed-stage embryos (at 50% eye pigmentation, Fig. S4) measured in closed, temperature-controlled respiration chambers (OXY0 41 A, Collotec Meßtechnik GmbH, Germany). The double-walled chambers were connected to a flow-through thermostat to adjust the chamber to the corresponding incubation temperature of the eggs. Measurements were conducted in triplicate with 10 to 20 eggs of each female and combination treatment. Eggs were chambered with a volume of 1 ml of sterilized seawater adjusted to the corresponding Pco2 treatment. A magnetic microstructure (3 mm) was placed against the floating eggs to avoid oxygen stratification within the respiration chamber. The change in oxygen saturation was detected by micro-optodes (fiber-optic microsensor, flat broken tip, diameter: 140 μm, PreSens GmbH, Germany) connected to a Microx TX3 (PreSens GmbH, Germany). Recordings were stopped as soon as the oxygen saturation declined below 80% air saturation. Subsequently, the water volume of the respiration chamber and wet weight of the measured eggs (gww) were determined by weighing (±1 mg). Oxygen consumption was expressed as[nmolO[nmolO2 (gww * min)−1]and corrected for bacterial oxygen consumption (<5%) and optode drift, which was determined by blank measurements before and after three successive egg respiration measurements.

Larval morphometrics

Subsamples of 10 to 30 nonmalformed larvae from each female and treatment combination were photographed for subsequent measurements of larval morphometrics (standard length, yolk-free body area, total body area, and yolk sac area) using Olympus image analysis software (Stream Essentials, Olympus, Tokyo, Japan). Only samples obtained from the same daily cohort (during peak hatch at each temperature treatment) were used for statistical comparison. After being photographed, 10 to 20 larvae were freeze dried to determine individual dry weights (±0.1 μg, XP6U Micro Comparator, Mettler Toledo, Columbus, OH, USA). Replicates with less than 10 nonmalformed larvae were precluded from statistical analyses.

Statistical analysis

Statistics were conducted with the open source software R, version 3.3.3 ( Linear mixed effect models[package“lme4”([package“lme4”(37)]were used to analyze data on egg survival and MO2. In each case, we treated different levels of temperature and Pco2 as fixed factors and included “female” (egg batch) as a random effect. Differences in larval morphometrics (yolk-free body area, total body area, dry weight, standard length, and yolk sac area) were determined by multifactorial analysis of covariance. These models were run with temperature and Pco2 as fixed factors and egg diameter as a covariate. Levene’s and Shapiro-Wilk methods confirmed normality and homoscedasticity, respectively. The package “lsmeans” (38) was used for pairwise comparisons (P values were adjusted according to Tukey’s post hoc test method). All data are presented as means (± SEM) and statistical tests with P < 0.05 were considered significant. Results are summarized in table S1.

Curve fitting

Generalized additive models[package“mgcv”([package“mgcv”(39)]were used to fit temperature-dependent curves of successful development building on egg survival and MO2. This method has the benefit of avoiding a priori assumptions about the shape of the performance curve, which is crucial in assessing the impact of elevated Pco2 on thermal sensitivity. “Betar” and “Gaussian” error distributions were used for egg survival and MO2 data, respectively. To avoid overfitting, the complexity of the curve (i.e., the number of degrees of freedom) was determined by penalized regression splines and generalized cross-validation (39) Models of egg survival were constrained at thermal minima because eggs of cold-water fish can survive subzero temperatures far below any applicable in rearing practice. Following Niehaus et al. (40), we forced each model with artificial zero values (n = 6) based on absolute cold limits from the literature. These limits were set to −4°C for Atlantic cod (41) and −9°C for Polar cod assuming similar freezing resistance, as reported for another ice-associated fish species from Antarctica (42)

Spawning habitat maps

Fitted treatment effects on normalized egg survival data (fig. S6A; raw data are shown in Fig. 3) were linked to climate projections for the Seas of Norden to infer spatially explicit changes in the maximum PES under different RCPs. That is, the treatment fits were evaluated for gridded upper-ocean water temperatures (monthly averages) bilinearly interpolated to a horizontal resolution of 1° × 1° and a vertical resolution of 10 m. To account for species-specific reproduction behavior, we first constrained each map according to spawning seasonality and depth preferences reported for Atlantic cod[MarchtoMay50to400m([MarchtoMay50to400m(33)]and Polar cod[DecembertoMarch5to400m([DecembertoMarch5to400m(13)]. As both species produce pelagic eggs that immediately ascend into the upper mixed layer if spawned at greater depths (13, 33), we further limited the eligible depth range to the upper 50 m. PES at a given latitude and longitude was then estimated from the calculations by selecting the value at the depth of maximum egg survival (at 0 to 50 m depth). Egg dispersal was not considered since the major bulk of temperature- and acidification-related mortality occurs during the first week of development (fig. S4).

Oceanic conditions were expressed as climatological averages of water temperatures, sea-ice concentrations, and the pH of surface water. Our observational baseline is represented by monthly water temperatures[WOA13([WOA13(43)]and sea-ice concentrations[HadISST([HadISST(44)], averaged from 1985 to 2004, and by pH values averaged over the period 1972–2013[GLODAPv2([GLODAPv2(45, 46)]. Simulated ocean climate conditions were expressed as 20-year averages of monthly seawater temperatures and sea-ice concentrations and of 20-year averages of annual pH values of surface water. End-of-century projections were derived from climate simulations for 2081–2100 carried out in CMIP5 (45) We considered only those 10 ensemble members (see table S4) that provide data on each of the relevant parameters (water temperature, sea ice, and pH) under RCP8.5, RCP4.5, and RCP2.6 (47) Projected pH values and temperatures are shown in fig. S6 (E to L). To account for potential model biases, we diagnosed for each of the 10 CMIP5 models the differences between simulations and observations for the baseline period and subtracted these anomalies from the CMIP5-RCP results for 2081–2100. For 2081–2100, we considered the CMIP5-RCPs ensemble median of maximum PES and assessed the uncertainty of PES at a given location by defining a signal-to-noise ratio that relates the temporal change in PES between 2081–2100 and 1985–2004 (ΔPES) to the median absolute deviation (MAD) of results for 2081–2100. Model results are not robust where the temporal change in PES is smaller than the ensemble spread, i.e., ΔPES/MAD < 1. PES calculations for scenarios RCP2.6 and RCP4.5 were carried out for Pco2 = 400 μatm. The effect of elevated Pco2 (1100 μatm) on PES was only considered under scenario RCP8.5.


Supplementary material for this article is available at

Fig. S1 Thermal niches of adult Atlantic cod and Polar cod.

Fig. S2 Treatment effects on larval morphometrics at hatch.

Fig. S3. Water quality measurements.

Fig. S4. Effects of temperature and Pco2 on daily mortality rates of Atlantic cod and Polar cod.

Fig. S5. Effects of temperature and Pco2 on embryonic development of Atlantic cod and Polar cod.

Fig. S6. Spawning habitat maps for Atlantic cod and Polar cod are based on experimental egg survival data and climate projections under different emission scenarios.

Table S1. Summary table for statistical analyses conducted on data presented in Figs. 2 and 3 of the main text and in figs. S1 and S5.

Table S2. Length and weight of female and male Atlantic cod and Polar cod used for strip spawning and artificial fertilization.

Table S3. Mean egg diameter and fertilization success of egg batches (±SD, n = 3) produced by different females (n = 6).

Table S4. List of CMIP5 models that met the requirements for this study (for details, see the “Spawning habitat maps” section in the main text).

References (4855)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


  1. H.-O. Pörtner et al., in Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (Cambridge Univ. Press, 2014), pp. 411–484.

  2. O. Hoegh-Guldberg, R. Cai, E. S. Poloczanska, P. G. Brewer, S. Sundby, K. Hilmi, V. J. Fabry, S. Jung, The Ocean, in Climate Change 2014: Impacts, Adaptation, and Vulnerability. Part B: Regional Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel of Climate Change, V. R. Barros, C. B. Field, D. J. Dokken, M. D. Mastrandrea, K. J. Mach, T. E. Bilir, M. Chatterjee, K. L. Ebi, Y. O. Estrada, R. C. Genova, B. Girma, E. S. Kissel, A. N. Levy, S. MacCracken, P. R. Mastrandrea, L. L. White, Eds. (Cambridge Univ. Press, 2014), chap. 30, pp. 1655–1731.

  3. J. Blindheim, The seas of Norden, in Norden: Man and Environment, U. Varjo, W. Tietze, Eds. (Gebrüder Borntraeger, 1987), pp. 20–32.

  4. A. M. Ajiad, H. Gjøsæter, in The Barents Sea. Ecosystem, Resources, Management. Half a Century of Russian-Norwegian Cooperation, T. Jakopsen, V. K. Ozhigin, Eds. (Tapir Academic Press, 2011), pp. 315–328.

  5. FAO, The State of World Fisheries and Aquaculture (SOFIA) (FAO Fisheries and Aquaculture Department, 2018).

  6. UNFCCC, Adoption of The Paris Agreement FCCC/CP/2015/L.9/Rev.1 (2015).

  7. K. Brander, Spawning and life history information for North Atlantic cod stocks, ICES Cooperative Research Report (2005).

  8. A. G. Dickson, C. L. Sabine, J. R. Christian, Guide to Best Practices for Ocean CO2Measurements (North Pacific Marine Science Organization, 2007).

  9. S. Wood, M. S. Wood, Package ‘mgcv’. R package version, 1.7-29 (2017).

  10. R. A. Locarnini, A. V. Mishonov, J. I. Antonov, T. P. Boyer, H. E. Garcia, O. K. Baranova, M. M. Zweng, C. R. Paver, J. R. Reagan, D. R. Johnson, M. Hamilton, D. Seidov, World Ocean Atlas 2013 (NOAA, 2013), vol. 1, pp. 73–44.

  11. R. M. Key, A. Olsen, S. van Heuven, S. K. Lauvset, A. Velo, X. Lin, C. Schirnick, A. Kozyr, T. Tanhua, M. Hoppema, S. Jutterström, R. Steinfeldt, E. Jeansson, M. Ishi, F. F. Perez, T. Suzuki, Global Ocean Data Analysis Project, Version 2 (GLODAPv2), ORNL/CDIAC-162, NDP-P093 (Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, US Department of Energy, 2015).

Acknowledgments: We acknowledge the support of S. Hardenberg, E. Leo, M. Stiasny, C. Clemmensen, G. Göttler, F. Mark, and C. Bridges. Special thanks are dedicated to the staff of the Tromsø Aquaculture Research Station and the Centre for Marine Aquaculture. Funding: Funding was received from the research program BIOACID [Biological Impacts of Ocean Acidification by the German Federal Ministry of Education and Research (BMBF), FKZ 03F0655B to H.-O.P. and FKZ 03F0728B to D.S.]. Funding was also received from AQUAculture infrastructures for EXCELlence in European fish research (AQUAEXCEL, TNA 0092/06/08/21 to D.S.). F.T.D., M.B., H.-O.P., and D.S. were supported by the PACES (Polar Regions and Coasts in a Changing Earth System) program of the Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research (AWI). Previous and additional support from grants POLARIZATION (Norwegian Research Council grant no. 214184 to J.N.) and METAFISCH (BMBF grant no. FZK01LS1604A to H.-O.P. and F.T.D.) are also acknowledged. Author contributions: F.T.D. and D.S. devised the study and designed the experiments. F.T.D. conducted the experiments. J.N., V.P., and A.M. provided equipment and facility infrastructure. F.T.D. analyzed the experimental data. M.B. analyzed climate data and generated habitat maps. F.T.D. drafted the manuscript. F.T.D., D.S., M.B., and H.-O.P. wrote the manuscript. J.N., V.P., and A.M. edited the manuscript. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors. The experimental data supporting the findings of this study are available from PANGEA (, a member of the ICSU World Data System.

Source link