# Global vulnerability of marine mammals to global warming

Camille Albouy, Valentine Delattre, Giulia Donati, Thomas Frölicher, Severine Albouy-Boyer, Marta Rufino, Loïc Pellissier, David Mouillot, Fabien Leprieur

<figure><img src="image_1.png" /></figure>

► To cite this version:

Camille Albouy, Valentine Delattre, Giulia Donati, Thomas Frölicher, Severine Albouy-Boyer, et al. Global vulnerability of marine mammals to global warming. Scientific Reports, 2020, 10 (1), 10.1038/s41598-019-57280-3. http://hal-03403188

HAL Id: hal-03403188

https://hal.umontpellier.fr/hal-03403188

Submitted on 25 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Corrected: Author Correction

OPEN

# Global vulnerability of marine mammals to global warming

Camille Albouy<fnref n="1" />*, Valentine Delattre², Giulia Donati³,⁴, Thomas L. Frölicher⁵,⁶, Severine Albouy-Boyer⁷, Marta Rufino⁸,⁹, Loïc Pellissier³,⁴, David Mouillot² & Fabien Leprieur²,¹⁰

Although extinctions due to climate change are still uncommon, they might surpass those caused by habitat loss or overexploitation over the next few decades. Among marine megafauna, mammals fulfill key and irreplaceable ecological roles in the ocean, and the collapse of their populations may therefore have irreversible consequences for ecosystem functioning and services. Using a trait-based approach, we assessed the vulnerability of all marine mammals to global warming under high and low greenhouse gas emission scenarios for the middle and the end of the 21st century. We showed that the North Pacific Ocean, the Greenland Sea and the Barents Sea host the species that are most vulnerable to global warming. Future conservation plans should therefore focus on these regions, where there are long histories of overexploitation and there are high levels of current threats to marine mammals. Among the most vulnerable marine mammals were several threatened species, such as the North Pacific right whale (Eubalaena japonica) and the dugong (Dugong dugon), that displayed unique combinations of functional traits. Beyond species loss, we showed that the potential extinctions of the marine mammals that were most vulnerable to global warming might induce a disproportionate loss of functional diversity, which may have profound impacts on the future functioning of marine ecosystems worldwide.

During the past few decades, the Earth has entered a new era of rapid and potentially irreversible climate warming due to the positive radiative imbalances triggered by greenhouse gas emissions from human activities<fnref n="1" />,². The oceans have taken up 93% of the extra energy<fnref n="1" /> that has been accumulated in the Earth system in recent decades, and its temperature has increased much faster since 1991 than has been recorded previously². In addition, changes in ocean temperature also affect the sea level, sea ice extent and salinity (through changes in precipitation and evaporation). All these changes were found to have negative impacts on marine biota³,⁴ and especially on marine mammals (e.g.⁵⁻¹⁰).

One of the most common responses of marine mammals to temperature changes is shifts in their spatial distributions, which could result in modifications of the ranges of the species (e.g.¹¹⁻¹³). For example, Bryde's whales (Balaenoptera brydei), a widely spread subtropical and tropical species, was increasingly detected in the cooler waters off southern California during the period from 2000–2010¹⁴. In contrast, the white-beaked dolphin (Lagenorhynchus albirostris), a cold-water species, has reduced its range and is declining in abundance¹². As mentioned by Elliott & Simmonds¹², geography could constrain species with low dispersal ability to a particular area. For example, the distribution of an endemic species of porpoise, the vaquita (Phocoena sinus), is limited to the northern end of the Gulf of California. Changes in water temperatures could alter the life cycles of the prey of marine mammal and provoke mismatches between the abundance of prey and those of marine mammals. This situation could be particularly critical for migratory species that travel long distances between feeding and breeding areas¹² and species that depend on prey as a source of protein for lactation or for weaned calves¹²,¹⁵. The increase in ocean temperatures could also affect the reproductive success of marine mammals, as has been reported for female sperm whales, which have lower conception rates after long exposure to higher sea surface temperatures (SSTs) than usual. Finally, increases in ocean temperatures can have direct impacts on the survival

<footnote n="1">¹IFREMER, unité Ecologie et Modèles pour l'Halieutique, rue de l'Ile d'Yeu, BP21105, 44311, Nantes, cedex 3, France.
²MARBEC, Univ Montpellier, CNRS, Ifremer, IRD, Montpellier, France.
³Landscape Ecology, Institute of Terrestrial Ecosystems, ETH Zürich, 8092, Zürich, Switzerland.
⁴Swiss Federal Research Institute WSL, 8903, Birmensdorf, Switzerland.
⁵Climate and Environmental Physics, Physics Institute, University of Bern, Bern, Switzerland.
⁶Oeschger Centre for Climate Change Research, University of Bern, Bern, Switzerland.
⁷Cadres en Mission Nantes, Nantes, France.
⁸MARE - Marine and Environmental Sciences Centre, Faculty of Sciences, University of Lisbon, Campo Grande, 1749-016, Lisboa, Portugal.
⁹CCMAR, The Centre of Marine Sciences, Universidade do Algarve, Campus de Gambelas, 8005-139, Faro, Portugal.
¹⁰Institut Universitaire de France, Paris, France.
*email: albouycamille@gmail.com</footnote>
rates of marine mammals (e.g.<sup>16,17</sup>) by increasing the stress of organisms, fostering the development of pathogens<sup>13</sup> and increasing the propagation of pathogens to new species by causing species to experience range shifts<sup>13</sup>. Despite the large potential impacts of ocean warming on marine mammals, the global vulnerability of marine mammals to global warming is poorly understood. This is an important knowledge gap because 37% of marine mammals are currently considered endangered by the IUCN (3 species are labeled as critically endangered, 13 as endangered and 12 as vulnerable<sup>18,19</sup>). In addition, marine mammals are important drivers of ecosystem functioning<sup>13</sup> and fulfill key ecological roles worldwide, mostly via trophic dynamics due to their role as consumers at most trophic levels and their role in nutrient cycling (e.g.<sup>20-24</sup>).

Given the substantial observed and projected decline in Arctic sea ice cover<sup>1</sup>, it has been generally assumed that Arctic marine mammals were the most vulnerable in the face of climate change (e.g.<sup>20,25,26</sup>). However, marine mammals living in temperate regions could also be negatively affected as warmer sea temperatures can modify the phytoplankton blooms in spring and deregulate the whole food chain by reducing the amount of input energy in the ecosystem<sup>20</sup>. Identifying the most vulnerable species worldwide under projected global warming and the regions where these species are concentrated could be useful for the prioritization of species-level conservation plans (e.g.<sup>27,28</sup>). Furthermore, exploring the potential consequences of climate change at the assemblage level would enable the assessment of the potential losses of species diversity and allow planning of better conservation strategies. To that end, this level of investigation often requires going beyond species richness as this metric fails to consider that species may have different ecological roles and evolutionary heritages (see<sup>29</sup>). Thus, focusing on the functional and phylogenetic components of biodiversity should be preferred to provide an integrated approach for biodiversity conservation<sup>30-33</sup>. Phylogenetic diversity reflects the evolutionary history of a given assemblage and is important for explaining the role of species interactions and biogeographic histories in structuring communities<sup>34</sup>. Functional diversity enables the quantification of the value and range of species traits that influence the performances of certain species and thus the ecosystem functioning<sup>35</sup>. Hence, measuring phylogenetic and functional diversity together as complementary components of biodiversity is essential for understanding the complete structures, compositions and dynamics of natural communities<sup>36-39</sup>.

In this study, we assessed the vulnerability of marine mammals under future global warming using a trait-based approach<sup>25,40,41</sup> (see Supplementary Methods 1 and 2). We first built a global index of the intrinsic sensitivity of all marine mammals to global warming, i.e., the magnitude of the species-specific responses to ocean warming from the inherent biological and ecological features. We then proposed an exposure index that would quantify the magnitude of ocean warming within a the geographical range of a species using the projections from 11 different global Earth system models under two future greenhouse gas emission scenarios: a low greenhouse gas emission scenario (Representative Concentration Pathway (RCP) 2.6 with a radiative forcing of 2.6 in the year 2100), and a high greenhouse gas emission scenario (RCP8.5). We used sea surface temperature (SST) as a proxy for ocean warming as this variable was found to be a good predictor of marine mammal distributions<sup>42-44</sup>. Finally, we derived a vulnerability index as the product of the intrinsic sensitivity and the exposure to ocean warming<sup>33</sup>. From this index, we also identified the global geographic hotspots of marine mammal vulnerability to global warming. Finally, we assessed the effects of the potential extinctions caused by climate change on the phylogenetic and functional diversity of marine mammals. For this, we compared several scenarios by removing species either at random or according to their vulnerability to global warming or their IUCN Red List extinction risk.

## Methods

**Intrinsic sensitivity to global warming.** We evaluated the intrinsic sensitivity of marine mammal species to ocean warming using a trait-based approach<sup>40,45-48</sup>. The species sensitivity to ocean warming was defined as the inability of the species to persist in its habitat<sup>40,48</sup>. Based on an extensive review of the literature, we compiled information on fifteen traits in five main categories (feeding, habitat, reproduction, social behavior and biology) for 123 marine mammal species (see full description in Supplementary Method 1)<sup>28</sup>. First, we considered the traits that have been commonly used in climate change vulnerability assessments<sup>40</sup> and for a wide range of taxa (see<sup>40,45-48</sup> i.e., the three dimensions of rarity<sup>49</sup>: geographical range-related variable<sup>50</sup>, diet specialization and habitat specialization, which reflect the plasticity of a species<sup>51</sup>). Second, we accounted for specific traits that are known to mediate the sensitivity of marine mammals to climate change, such as traits related to reliance on pack ice<sup>25</sup>. Several traits specific to habitat requirements were obtained through analysis of the IUCN geographical range data (http://iucnredlist.org) and the Bio-ORACLE raster environmental spatial data (http://www.oracle.ugent.be; see Supplementary Method 1 for more details).

To provide a quantitative ranking of species sensitivity to global warming, we extended the methodology proposed by<sup>25</sup>, in which each considered trait was associated with specific quantitative ranking criteria that were evaluated on a three-point scale (with 0 being the least sensitive and 2 being the most sensitive). The ranks of the quantitative traits were evaluated by dividing the ranges of the values into three equal parts. For those quantitative traits, we performed a statistical sensitivity test on the variation of ranks from 1% to 33% (see Supplementary Method 2). When considering the qualitative traits, the modalities were ranked according to the relative influence of each modality on sensitivity to ocean warming. For example, epipelagic feeders are expected to be more sensitive to ocean warming than benthic or mesopelagic feeders since the upper layer of oceans are more exposed to warming than the deep layers<sup>1</sup> (Table S1). Hence, the epipelagic modality of the 'habitat vertical specialization' trait was labeled as 2 (most sensitive), while the benthic and mesopelagic modalities were both labeled as 1, and the habitat-generalist species were 0 (least sensitive). The sum of all trait values (between 0 and 2) resulted in an overall species-specific sensitivity ranking that allowed direct comparison between species. The resulting values were then divided by the maximum sensitivity value to set the sensitivity index between 0 (least sensitive species) and 1 (most sensitive species).
**Exposure and vulnerability to climate change.** The vulnerability of the species to global warming (V) was defined as the product of the intrinsic sensitivity (S) multiplied by the exposure (P)$^{40}$:

$$V = P \times S$$

Through this method, a species that was intrinsically sensitive to global warming would be considered vulnerable only if it was actually exposed to any change in temperature within its geographical range. Exposure was defined as the magnitude of the changes in SST within the geographical range of each species between the baseline in the present-day period (1971–2000) and the values in a near-future period (2030–2059) and a period at the end of 21st century (2070–2099)$^{52}$. To assess the uncertainties in the models and scenarios, we used 11 different Earth system models that participated in the Coupled Model Intercomparison Project phase 5 (CMIP5)$^{53}$ and that were run for the 1861–2100 period under the historical, RCP8.5 and RCP2.6 scenarios$^{54}$. We regridded the model output data onto a regular 1° × 1° grid (~10,000 km²) by using the Krige function of the R package gstat. We selected the following CMIP5 ESMs: MRI-CGCM3, IPSL-CM5A-LR, GFDL-ESM2G, GFDL-ESM2M, IPSL-CM5A-MR, MIROC-ESM, MPI-ESM-LR, GFDL-CM3, CSIRO, and CanESM2. We selected high and low greenhouse gas emission scenarios among the set of RCPs that were run with the CMIP5 models to explore the impacts of different socioeconomic climate policies on the climatic projections$^{55}$. RCP8.5 is a high-carbon emission, business-as-usual scenario with a net radiative forcing of 8.5 W/m² by 2100 and a simulated global atmospheric surface temperature change of approximately 3.7 ± 0.7 °C relative to the temperature from 1986–2005. RCP2.6 is a strong emission mitigation scenario with a global mean net radiative forcing of 2.6 W/m² by 2100 and a simulated global atmospheric surface temperature change of approximately 1.0 ± 0.4 °C$^{56}$ relative to 1986–2005. The Earth system model output was regridded onto a regular 1° × 1° grid (~10,000 km²) using the Krige function of the R package gstat, and we extracted the average value of the SST for both the baseline records and the projected SST from each model and for the two emissions scenarios.

We then calculated the changes in the SST at the global scale within the geographical range of each species from the baseline to the two future periods for both the RCP2.6 and RCP8.5 scenarios. We downloaded the geographical range maps from the IUCN (http://www.iucnredlist.org) for the 123 studied species and then derived a presence/absence matrix by overlapping the geographical ranges and counting how many species occurred in each grid cell (1° × 1° grid cells). To quantify the degree of exposure to ocean warming at the species level, the distributions of the absolute values of δ SST were divided into three equal parts according to statistical considerations. The lowest category was coded as 0, the middle category as 1 and the upper as 2. For each species geographical range that was divided into grid cells (1° × 1°), we counted the relative number of cells occurring in each category and deducted the exposure value according to the following equation:

$$P = a \times 0 + b \times 1 + c \times 2$$

where a, b and c are the percentage of cells in the lower, middle and upper categories, respectively. This index accounts for the heterogeneity in δ SST throughout the geographical range of a species and is independent of the geographical range size. Consequently, a species that will face great temperature changes (either increases or decreases) in a large proportion of its range was categorized as highly exposed to climate change, while a species that will likely face lower δ SST within its whole range was categorized as less exposed.

**Diversity erosion scenarios.** We investigated the potential erosion of marine mammal diversity under global warming through two complementary biodiversity components: the phylogenetic and functional diversities. Phylogenetic diversity was measured with the PD index$^{57}$ as the minimum total length of all the phylogenetic branches required to represent a given set of taxa on the phylogenetic tree$^{57}$. The PD calculation was applied to a pruned supertree of mammals$^{58}$. However, this supertree contained a large number of polytomies, which implied potentially poor branch length estimates in parts of the tree. To account for this phylogenetic uncertainty in our analyses, we also quantified PD using 100 trees that were chosen at random from the Bayesian posterior distribution of the fully resolved trees that were generated by a polytomy resolver$^{59}$.

The functional diversity was quantified using an index of functional richness (FRic$^{60}$). The FRic index relies on a multidimensional Euclidean space where the axes are functional traits (or factorial axes from a principal coordinates analysis (PCoA) computed using these traits) along which species are placed according to their trait values$^{60}$. This index measures the volume of the functional space that is occupied by a given species assemblage by calculating the convex hull volume, which is defined by the species at the vertices of the functional space that encompasses the entire trait space filled by all species in the assemblage$^{60}$. To calculate the FRic index, we used 13 functional traits, i.e., main diet, foraging depth range, foraging location, fasting strategy, female sexual maturity, weaning, gestation length, interval between litters, breeding sites, social group size, social behavior, adult maximum body mass and sexual dimorphism (see$^{28}$ for full details). These traits covered five functions of marine mammals (i.e., feeding, reproduction, habitat, social behavior and intrinsic biology) and reflected the trade-offs in resource allocation$^{20}$. The pairwise functional distances between the species were computed using Gower's distance, which allows the mixing of different types of variables while weighting the functional traits in order to give equal weight to each function (i.e., reproduction, social behavior, biology, feeding and habitat). Then, a principal coordinates analysis (PCoA) was performed using this functional distance matrix$^{61}$. According to the metric proposed by$^{62}$ to evaluate the quality of a functional space (i.e., the extent to which it is a faithful representation of the initial functional trait values), we retained the first four principal axes of the PCoA to build the multidimensional functional space. We preferred to use the FRic index instead of the FD index, which is based on a functional dendrogram, as a recent study showed that the FD index may lead to biased assessments of functional diversity and inaccurate ecological conclusions$^{62}$. In addition, additional analyses showed that the distances derived from a UPGMA dendrogram did not faithfully represent the initial Gower's distances.
For each diversity index (PD and FRic), the diversity was initially measured for the global pool of marine mammals. Extinctions were simulated by removing individual species from the global pool and then recalculating the diversity of the remaining pool. The order of the species to be removed from the global pool varied among four types of scenario. First, we developed a vulnerability scenario in which the order of the simulated extinctions was given by the ranking of the vulnerability of a species to climate change, where the most vulnerable species were assumed to go extinct first. Second, we constructed 999 random scenarios in which the order of the species to be removed from the pool was chosen randomly. We computed the mean and SD values of these scenarios. Third, we implemented a scenario that coupled the IUCN classification and the ranking of vulnerability to climate change, and we simulated extinctions by ordering the groups of species according to the IUCN categories (i.e., critically endangered (CR), endangered (EN), vulnerable (VU), near threatened (NT), least concern (LC) and data deficient (DD)). Within each group of species in the same IUCN category, the species were ranked according to their vulnerabilities to climate change. Hence, the species that was simulated to go extinct first was the most vulnerable to climate change among those classified as critically endangered. Finally, we implemented 999 random IUCN scenarios where extinction was predicted by ordering the groups of species according to the IUCN categories. We randomized the order of species removal within each group. We used the mean and SD for the graphical representation.

**Functional and phylogenetic originality.** To explore whether the most vulnerable species to climate change were the most "original" from a functional or phylogenetic point of view, we calculated two distinct metrics. First, we quantified the phylogenetic originality (or evolutionary distinctiveness) of each species using a metric that represented the relative contributions of the species to the phylogenetic diversity<sup>63</sup>. We calculated the mean evolutionary distinctiveness per species using the 100 resolved phylogenetic trees. Second, we derived the functional originality of each species from Gower's distance matrix<sup>36</sup>. This metric was calculated as the distance between each species and its nearest neighbor in the functional space and hence represented the degree of the isolation of a species in the functional space occupied by a given assemblage<sup>36</sup>.

## Results

**Projected climate change.** The CMIP5 Earth system models project an increase in global SST of 0.7 ± 0.51 °C (RCP2.6) and 2 ± 1.3 °C (RCP8.5), on average, between the baseline period (1971–2000) and the end of the 21st century (2070–2099). However, the projected increase in SST will not be homogeneous throughout the world ocean (see Figs. S1–S3 for the period 2030–2059). For example, by the end of the century and according to the RCP85 scenario, the Barents Sea, the Bering Sea and the Tasmanian Sea are projected to experience increases in SST of up to 4.5 °C, while other regions such as the Labrador Sea or the Southern Ocean are projected to show low levels of warming or even decreasing SSTs (Figs. S1, S2).

**Sensitivity, exposure and vulnerability to climate change.** The index of sensitivity to ocean warming, which ranged ranging from 0 to 1, was based on fifteen species traits reflecting the biological and demographic responses of marine mammals to climate change (see Methods and Supplementary Method 1). The sensitivity of the species to ocean warming (S) ranged from 0.35 to 1 (see Fig. S4; Table S2), with a mean value of 0.69 (sd = 0.12). The walrus (*Odobenus rosmarus*) and the North Pacific right whale (*Eubalaena japonica*) were the most sensitive species (S = 1), followed by the gray whale (*Eschrichtius robustus*; S = 0.92), the crabeater seal (*Lobodon Carcinophaga*; S = 0.92), the North Atlantic right whale (*Eubalena glacialis*, S = 0.88), the dugong (*Dugong dugon*, S = 0.88) and the narwhal (*Monodon Monoceros*; S = 0.88; Table S2).

At the end of the 21st century, according to the RCP2.6 and RCP8.5 scenarios, respectively, the exposure of the species to climate change ranged from 0 (not exposed to SST variation throughout the geographic range of the species between the baseline period and the end of the century) to 0.92 and 1.2 (most exposed to ocean warming), with mean values of 0.34 and 0.62, respectively (see Fig. S4; Tables S2, S3).

Combining the sensitivity and exposure indices allowed us to assess the vulnerability of marine mammals to global warming. This was best exemplified by the whole species belonging to the Eubalaena genus. For example, the North Pacific right whale (*Eubalaena japonica*) and the southern right whale (*Eubalaena australis*) were both highly sensitive to climate change according to their intrinsic biological and ecological characteristics but showed contrasting exposure to climate change within their current geographical ranges (Fig. 1). In fact, the North Pacific right whale would be more exposed to climate change at the end of the 21st century than the southern right whale according to both RCP scenarios, which makes the former species more vulnerable to global warming (Fig. 1). Overall, at the end of the century, the vulnerability index ranged from 0 to 0.74 and 1.06, according to the RCP2.6 and RCP8.5 scenarios, respectively, with mean values of 0.14 (±0.12 sd) and 0.42 (±0.22 sd) (see Fig. S4; Tables S2, S3).

At the end of the 21st century, the ranking of the 20 most vulnerable species based on the RCP8.5 scenario reveals the overrepresentation of IUCN-recognized endangered or vulnerable species (endangered = 40%, vulnerable = 15%) compared to the RCP2.6 scenario (RCP2.6: endangered = 5%, vulnerable = 5%; Fig. 2a). Despite these differences, the North Pacific right whale (*Eubalaena japonica*) and the gray whale (*Eschrichtius robustus*) were found to be the most vulnerable species to ocean warming according to both RCP scenarios (Fig. 2a). The vulnerability of a species to ocean warming had no phylogenetic signal (*Pagel* λ < 0.1 regardless of the RCP scenario and the time period considered; values obtained from over 100 resolved phylogenetic trees), indicating that future climatic impacts were not clustered within particular lineages of marine mammals.

At the end of the 21st century and for both RCP scenarios, no clear relationship was found between the vulnerabilities of different species to ocean warming and their levels of evolutionary distinctiveness, which represented their relative contributions to the evolutionary history (or phylogenetic diversity) in the global marine mammal fauna (Spearman rank correlation test, RCP2.6: *p* = 0.045, *p* = 0.63; RCP8.5: *p* = 0.025, *p* = 0.78; see Figs. 2b and
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Projected climate change within the geographical species ranges. Projected climate change within the geographical range of the North Pacific right whale (a) and the Southern right whale. (b) Colors show the magnitudes of the difference in the mean sea surface temperature between the baseline period (1971–2000) and the end of the century (2070–2099) according to the RCP8.5 scenario. The geographical range is represented in bright colors that contrast with the dark background. These maps and barplots allow visualization of the exposure, i.e., the degree of climate change (either warming for positive values of SST or cooling for negative values) that a species is likely to face within its range. Maps were created using R 3.6.0 software (https://www.r-project.org/).</figcaption></figure>

S5). Nonetheless, a few species were simultaneously highly vulnerable to climate change and evolutionarily distinct (Fig. 2b), namely, the dugong (*Dugong dugon*), the African manatee (*Trichechus senegalensis*) and the Florida manatee (*Trichechus manatus*). Moreover, these three species are classified as vulnerable according to the IUCN Red List.

Similarly, no clear relationship was found between the vulnerabilities of different species to ocean warming and their levels of functional originality (i.e., the degree of uniqueness of the species traits in the global pool of marine mammals) for both RCP scenarios at the end of the 21st century (Spearman rank correlation test, RCP2.6: ρ = 0.17, p = 0.07; RCP8.5: ρ = 0.011, p = 0.21; Fig. S5). However, contrary to the results for evolutionary distinctiveness, 9 and 17 species were both vulnerable to climate change and functionally distinct according to the RCP2.6 and RCP8.5 scenarios, respectively, at the end of the 21st century. In the RCP8.5 scenario, 3 of the 17 species were classified as endangered, 10 as least concern, 1 as near threatened and 3 as vulnerable. For example, the North Pacific right whale (*Eubalaena japonica*), the gray whale (*Eschrichtius robustus*), the dugong (*Dugong dugon*), the northern elephant seal (*Mirounga angustirostris*) and the walrus (*Odobenus rosmarus*) were both vulnerable to climate change and functionally unique (Figs. 2c and S5). Two of these species were considered to be endangered or vulnerable according to the IUCN Red List, namely, the North Pacific right whale and the dugong, respectively.

When considering the vulnerability of marine mammals to global warming at the assemblage level, our results showed a latitudinal gradient for both RCP scenarios at the end of the 21st century (Fig. 3c,d). The mean vulnerability to climate change was found to be higher in the northern hemisphere, peaking at the 57°5′ and 56°5′ latitudes for RCP2.6 and RCP8.5, respectively. The areas with the highest mean vulnerabilities to climate change
<figure><img src="image_4.png" /><figcaption>RCP 8.5</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_7.png" /></figure>

**Figure 2.** Vulnerability of marine mammals to climate change. (a) Ranking of the 20 most vulnerable marine mammals to climate change according to two RCP scenarios: RCP2.6 (left) and RCP8.5 (right). The current IUCN status is indicated by a color range from red (critically endangered) to light blue (least concern) and the data deficient status is shown in violet. (b) Relationship between the intrinsic sensibility of each marine mammal to climate change and its exposure to projected climate change within its geographical range according to the RCP8.5 scenario. Each bubble is proportional to the level of functional originality (i.e., the measure of the extent to which a given species is functionally unique among the global marine mammal fauna). In the top right corner, the drawings of the most vulnerable species also display the highest levels of evolutionary distinctiveness, namely, the dugong, the African manatee and the Florida manatee. (c) Relationship between the intrinsic sensitivity of each marine mammal to climate change and its exposure to projected climate change within its geographical range according to the RCP8.5 scenario. Each bubble is proportional to the level of evolutionary distinctiveness (i.e., which represents the relative contribution of a species to the evolutionary history of the global marine mammal fauna). In the top right corner, the drawings of the most vulnerable species also display the species with highest levels of functional originality, namely, the North Pacific right whale, the gray whale, the dugong, the Northern elephant seal and the walrus.
<figure><img src="image_8.png" /><figcaption>**Figure 3.** Projected changes in sea surface temperatures and the associated vulnerabilities of marine mammals. Projected changes in the sea surface temperature between the baseline period (1971–2000) and the end of the century (2070–2099) following the RCP2.6 (a) and RCP8.5 (b) scenarios as simulated by 11 different CMIP5 Earth system models (MRI-CGCM3, IPSL-CM5A-LR, GFDL-ESM2G, GFDL-ESM2M, IPSL-CM5A-MR, MIROC-ESM, MPI-ESM-LR, GFDL-CM3, CSIRO, CanESM2). Local assemblage-level vulnerability of marine mammals to climate change for the period (2070–2099) based on (c) the RCP26 scenario and (d) the RCP8.5 scenario. To evaluate the vulnerability of marine mammals to climate change at the assemblage level, we averaged the vulnerability of each species in each grid cell (1° × 1° grid ells, ~10,000 km²). Maps were created using R 3.6.0 software (https://www.r-project.org/).</figcaption></figure>

were the Baltic, Black and Caspian Seas; the Bering Sea; the area along the Japanese coast; and the northern part of the Pacific Ocean along the Canadian and American coasts (Figs. 3c,d, and S3 for the period from 2030–2059).

**Scenarios of functional and phylogenetic diversity erosion.** Regardless of the considered RCP, the vulnerability scenarios that consists of removing species from the global pool of marine mammals according to their level of vulnerability to climate change (i.e., from the most to the least vulnerable species), showed a faster decrease in functional diversity (i.e., the amount of functional trait space occupied by the species, as measured by the FRic index⁵⁰) than would be expected by chance by the end of the 21st century (Figs. 4a,b and S6 for the period 2030–2059). For example, the extinction of the top 20 most vulnerable marine mammals as a result of climate change, according to the RCP8.5 scenario, would induce a functional diversity loss of approximately 20% (Fig. 4a) and a decrease of 13% compared with a random loss.

The vulnerability-IUCN scenario, which consisted of removing species according to their vulnerabilities to climate change by considering the IUCN status, showed that the extinctions of species within the least concern category would induce a faster loss of functional diversity than expected by a random scenario of species removal within each IUCN category (Figs. 4b and S6). Finally, the removal of the most vulnerable species, based on both RCP scenarios for the end of the 21st century, showed that the loss of phylogenetic diversity would not be different than expected with random extinctions (Figs. 4c and S6). Similar results were found when considering the IUCN-vulnerability scenario (Figs. 4d and S6).

Discussion

Marine mammals are threatened both directly and indirectly by different human activities, such as fishing, whale watching, vessel collisions, acoustic disturbance, pollution, and modification or loss of available habitats⁶⁴. Consequently, 37% of marine mammals are already included in the IUCN Red List (3 species are labeled as critically endangered, 13 as endangered and 12 as vulnerable⁶⁵). However, the marine mammal vulnerability to global warming has never been assessed at the species and global levels. Here, using a species-level trait-based approach, we showed that many marine mammal species distributed across the northern hemisphere and belonging to different taxonomic groups (e.g., whales, dolphins and seals) were highly vulnerable to global warming (Fig. 3), even under a strong mitigation scenario (i.e., RCP2.6), which, according to the CMIP5 simulations, gives a two in three chance of limiting global warming to below 2 °C. A key finding from our study was that the North Pacific, which
<figure><img src="image_9.png" /><figcaption>**Figure 4.** Scenario of the potential loss of functional richness and phylogenetic diversity in the global marine mammal fauna. Potential loss of functional richness (**a,b**) and phylogenetic diversity (**c,d**) in the global marine mammal fauna according to the RCP8.5 scenario for the end of the century (2070–2099). The upper graphs (**a,b**) show two erosion scenarios; one was repeated 999 times with a randomly selected order of species extinctions (dark continuous line) and is shown with the corresponding standard deviation, while the second had the species extinctions occur in order of their vulnerabilities to climate change (green continuous line). The lower graphs (**c,d**) show two random erosion scenarios. One was repeated 999 times with a randomly selected order of species extinctions (dark continuous line) and a second where the species were grouped according to their IUCN categories (critically endangered (CR), endangered (EN), vulnerable (VU), near threatened (NT), least concern (LC) and data deficient (DD)) prior to randomization. Finally, the lower graph presents a scenario (“Vulnerability & IUCN”, multicolor line) with the species grouped by IUCN categories and placed in decreasing order within each group according to their vulnerabilities to climate change. The ranking of the species differs between each pair of vulnerability scenarios, preventing direct comparison.</figcaption></figure>

has already been identified as a hotspot of human threats for marine mammals¹⁹,⁶⁴, is also a hotpot of vulnerability to global warming for this group. This implies that marine mammals in this region face double jeopardy from both human activities (e.g., marine traffic, pollution and offshore oil and gas development) and global warming, with potential additive or synergetic effects and as a result, these ecosystems face irreversible consequences for marine ecosystem functioning.

The marine mammals that are the most sensitive to climate change generally show marked feeding and habitat specialization, as well as reduced or fragmented geographical ranges (e.g., the dugong *Dugong dugon* and the walrus *Odobenus rosmarus*), which is consistent with the results of earlier studies that focused on diverse terrestrial and marine organisms²⁵,⁴⁰,⁶⁶⁻⁶⁸. Dependence on sea ice also seems to be a common denominator for many of the most sensitive marine mammals that are most sensitive to climate change²⁵ (e.g., the polar bear *Ursus maritimus* and the beluga whale *Delphinapterus leucas*).

By definition, the vulnerability index was higher for species that were both sensitive and exposed to global warming. The most vulnerable species according to both future RCP scenarios was the North Pacific right whale (*Eubalaena japonica*). This species has already been identified by the IUCN as an endangered species because of the absence of evidence of a recovering trend and the extremely low estimated number of individuals (~400 in the Okhotsk Sea and ~100 in the rest of the North Pacific)⁶⁹. While historically common in many areas of the North Pacific and Bering Sea, the North Pacific right whale was strongly depleted by intensive whaling in the mid-19th century⁷⁰. The second most vulnerable species was the gray whale (*Eschrichtius robustus*), which is classified as least concern by the IUCN, as its population size has been estimated to be above the threshold for the threatened category and the eastern subpopulation has increased over the last three generations⁷¹. However, the gray whale subpopulations in the western North Pacific are listed as critically endangered by the IUCN. Indeed, decades of commercial whaling led to the complete collapse of this subpopulation in much of its historical range. The Pacific right whale and the gray whale were also found to be functionally unique, i.e., with unique combinations
of functional traits. Their potential extinction could have large consequences on marine ecosystem functioning⁷². For example, the gray whale can resuspend large amounts of sediment and nutrients in the water column, which in turn enhances nutrient cycling and brings some benthic crustaceans that serve as food for seabirds to the ocean surface⁷³. For this species, a recent study showed that dispersal occurred between the eastern North Pacific and the North Atlantic during the Pleistocene and Holocene warming periods⁶³, following the opening of the Bering Strait. While this study suggests that climate warming may create favorable thermal conditions in the North Atlantic, where the gray whale was historically present, it remains uncertain if the eastern Pacific population will be able to establish a persistent population in the North Atlantic over a long time period. As suggested by⁷⁴, the northern limit of the range of this species is only occupied during the summer months, and individuals following the traditional migratory routes may have limited movement through the Bering Strait. Then, the ongoing opening of the Northwest Passage, which is a consequence of climate warming, has created new opportunities for marine traffic and gas and oil activities, which may negatively impact gray whale migrations⁷³,⁷⁴. Overall, the gray whale and the Pacific right whale should be of particular concern for conservation prioritization given their high levels of vulnerability to climate change, their high functional originality and the current threats that they are facing.

For some species, such as the dugong or the walrus, the level of vulnerability differed between the RCP2.6 and the RCP8.5 scenarios. Both species were shown to be highly sensitive to climate change but not highly exposed, according to the RCP2.6 scenarios (see Appendices 6 and 7). In contrast, under the RCP8.5 high emission scenario, these two species would be intensively exposed to global warming and hence would be highly vulnerable (see Appendices 6 and 7). This result highlights that economic and political decisions towards the reduction of CO₂ emissions and the mitigation of global warming can have serious consequences for the level of threat a species will face. This is particularly important because these species, while not classified in the top 20 most vulnerable species, were shown to be highly distinct in their evolutionary histories, in the case of the dugong (Fig. 2b), and in their traits, in the case of the walrus (Fig. 2c). Thus, the extinction of these species could provoke the loss of unique and important evolutionary lineages as well as a disruption in ecosystem functioning.

Among the top 20 most vulnerable species to climate change, some are still abundant within their geographical ranges and have shown no evidence of population declines, according to the IUCN. This is particularly the case for the spotted seal (Phoca largha) or the Pacific white-sided dolphin (Lagenorhynchus obliquidens), which are classified as least concern (Fig. 2a). Several species are also classified as data deficient, such as Baird's beaked whale (Berardius bairdii). For these species, we recommend considering their potential vulnerabilities to climate change when setting their IUCN statuses. This is even more critical than the potential extinction of the most vulnerable species within the least concern IUCN category, which may lead to a disproportional loss of functional diversity (Fig. 4b) and could ultimately disrupt marine ecosystem functioning⁷⁵. More generally, it would be particularly useful to consider the level of vulnerability to climate change when evaluating the IUCN statuses of marine mammals to implement mitigation actions on species that are not currently threatened (or have insufficient data) but that could be threatened in the near future⁷⁶.

Scaling up our vulnerability analysis to the species assemblage level showed that, regardless of the RCP scenario, the North Pacific Ocean, Greenland Sea and Barents Sea (Fig. 3c,d) hosted the marine mammals that were most vulnerable to global warming. Indeed, these regions might face the strongest effects from global warming under both emissions scenarios (Fig. 3a,b) and have already undergone temperature increases 2–3 times higher than the changes to the global mean surface temperature over the past 150 years¹. These regions should therefore be of particular concern for spatial monitoring for the conservation of marine mammals. Indeed, there are multiple threats to the marine mammals⁷⁷ in the North Pacific Ocean, the Greenland Sea and the Barents Sea, such as marine traffic and offshore oil and gas exploitation that are known to impact cetaceans, such as through ocean noise pollution. These current threats could have additive or synergetic effects with climate change, which may therefore increase the overall vulnerability of marine mammals. In that context, future studies should focus on evaluating whether the combined effects of habitat loss, marine pollution and climate change are greater than the effects of each threat individually, which is challenging but crucial for providing effective conservation actions for marine mammals. In addition, these studies should consider the use of a finer spatial resolution to account for the species-habitat relationships (e.g.,⁷⁸), which is required to evaluate the potential impacts of climate change on the regional spatial distributions of marine mammals.

Our results also suggested that the potential extinction of the marine mammals that are most vulnerable to global warming may lead to a disproportionate loss of functional diversity in the global marine mammal fauna (Fig. 4a); many of the most vulnerable species displaying a high level of functional originality (see Fig. S6). This projected loss of functional diversity may ultimately threaten marine ecosystem stability and service provisioning³³,⁷⁵. Similar findings have been reported for birds⁷⁷,⁷⁹ and marine fishes⁸⁰ under climate change scenarios. These findings have a particular resonance in the context of the recent geological past. During the Plio-Pleistocene, large climatic oscillations and sea level changes caused numerous extinctions among the global marine megafauna, which directly induced an important loss of functional diversity that was not fully compensated by the evolution of new genera⁸¹. Marine megafauna have therefore been more sensitive to past environmental changes than previously assumed⁸², which implies that future climate change could pose a great challenge for large marine animals, especially mammals. However, our findings should be extended by downscaling our approach to a finer spatial resolution to assess the potential effects of climate change on the functional trait compositions of the local assemblages (for marine fishes⁸⁰ and for birds⁷⁹), which could have useful applications for prioritization in spatial conservation planning.

While the trait-based approach used in this study has the advantage of coupling the magnitude of climate change that a species would experience (exposure) with its intrinsic ability to tolerate changing climate (sensitivity), it nevertheless has some limitations. First, we did not consider the adaptive capacity of species to avoid the negative impacts of ocean warming through dispersal or microevolutionary changes⁴⁰,⁴¹. However, the intrinsic
and extrinsic abilities of marine mammals to colonize new favorable habitats and to avoid the impacts of ocean warming¹¹ are poorly documented (but see⁷³). For example, the population size that is often used as an indicator of adaptive capacity (i.e., species that are numerous should have more options for adapting or reestablishing themselves in local or new areas²⁵) is not uniform among all marine mammals. We therefore decided to focus on the sensitivity and exposure dimensions to avoid biasing our assessment towards the most studied marine mammals. Nevertheless, we recognize that adaptive capacity is an important component of climate change vulnerability⁴⁰,⁴¹, and we encourage future studies on marine mammals to focus on this issue. Second, when estimating exposure to global warming, we did not consider all components of the oceanic conditions under projected climate change scenarios, such as the loss of ice sea cover or changes in prey availability, which also influence the distributions of marine mammals⁸². We did not include the changes in Arctic sea ice cover because the CMIP5 models do not adequately represent the observed loss of Arctic sea ice cover over the last few decades, i.e., they largely underestimate this loss⁸³. Our results should also be interpreted in the lith of CMIP5 climate models that do not include the effects of meltwater from the ice sheets and ice shelves of Antarctica⁸⁴. In addition, the projections of prey availability under climate change are currently not available at the global scale. We therefore encourage the development of climatic-biological coupled models, such as the “Model of Intermediate Complexity for Ecosystem Assessments” (MICE), which links krill and whale population dynamics with climate change drivers (i.e., changes in sea surface temperature, primary productivity and sea ice). For example, a MICE for the Southern Ocean was developed⁸² and showed that the projected declines in several whale species by 2100 could be the consequences of reduced prey from warming and increased interspecific competition between whale species. This modelling approach would complement to our trait-based approach by providing a mechanistic understanding of the vulnerabilities of various species to climate change.

## Data availability

All traits data used to evaluate the species sensibility to climate change along with presence absence matrix and the phylogenetical trees are available: https://figshare.com/articles/Input_data_for_Global_vulnerability_of_marine_mammals_to_global_warming_11323304.

Received: 29 August 2019; Accepted: 20 December 2019;

Published online: 17 January 2020

## References

1. IPCC, Summary for Policymakers. In: IPCC Special Report on the Ocean and Cryosphere in a Changing Climate [H.-O. Pörtner, D. C. Roberts, V. Masson-Delmotte, P. Zhai, M. Tignor, E. Poloczanska, K. Mintenbeck, M. Nicolai, A. Okem, J. Petzold, B. Rama, N. Weyer (eds.)]. In press (2019).

2. Cheng, L., Abraham, J., Hausfather, Z. & Trenberth, K. E. How fast are the oceans warming? *Science*, **363**, 128–129 (2019).

3. Lewandowska, A. M. et al. Effects of sea surface warming on marine plankton. *Ecol. Lett.* **17**, 614–623 (2014).

4. Faillettaz, R., Beaugrand, G., Goberville, E. & Kirby, R. R. Atlantic Multidecadal Oscillations drive the basin-scale distribution of Atlantic bluefin tuna. *Sci. Adv.* **5** (2019).

5. Serpetti, N. et al. Impact of ocean warming on sustainable fisheries management informs the Ecosystem Approach to Fisheries. *Sci. Rep.* **7**, 13438 (2017).

6. Chambault, P. et al. Sea surface temperature predicts the movements of an Arctic cetacean: the bowhead whale. *Sci. Rep.* **8**, 9658 (2018).

7. Derville, S. et al. Whales in warming water: Assessing breeding habitat diversity and adaptability in Oceania's changing climate. *Glob. Chang. Biol.* **25**, 1466–1481 (2019).

8. Poff, N. L. et al. Functional trait niches of North American lotic insects: traits-based ecological applications in light of phylogenetic relationships. *J. North Am. Benthol. Soc.* **25**, 730–755 (2006).

9. Simmonds, M. P. & Isaac, S. J. The impacts of climate change on marine mammals: early signs of significant problems. *Oryx* **41**, 19 (2007).

10. Wild, S. et al. Long-term decline in survival and reproduction of dolphins following a marine heatwave. *Curr. Biol.* **29**, R239–R240 (2019).

11. Kaschner, K., Tittensor, D. P., Ready, J., Gerrodette, T. & Worm, B. Current and Future Patterns of Global Marine Mammal Biodiversity. *PLoS One* **6**, e19653 (2011).

12. Elliott, W. & Simmonds, M. Whales in Hot Water? The Impact of a Changing Climate on Whales, Dolphins and Porpoises: A call for action. *WWF-International, Gland Switz./WDCS, Chippenham, UK* (2007).

13. Evans, P. G. H. & Bjørge, A. Impacts of climate change on marine biodiversity. *MCCIP Sci. Rev.* 134–148, https://doi.org/10.14465/2013.arc15.134-148 (2013).

14. Kerosky, S. M. et al. Bryde's whale seasonal range expansion and increasing presence in the Southern California Bight from 2000 to 2010. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **65**, 125–132 (2012).

15. Schumann, N., Gales, N. J., Harcourt, R. G. & Arnould, J. P. Y. Impacts of climate change on Australian marine mammals. *Aust. J. Zool.* **61**, 146 (2013).

16. Harvell, C. D. et al. Climate warming and disease risks for terrestrial and marine biota. *Science* **296**, 2158–62 (2002).

17. Burek, K. A., Gulland, F. M. D. & O'Hara, T. M. Effects of climate change on arctic marine mammal health. *Ecol. Appl.* **18**, S126–S134 (2008).

18. Schipper, J. et al. The Status of the World's Land and Marine Mammals: Diversity, Threat, and Knowledge. *Science*. **322**, 225–230 (2008).

19. Davidson, A. D. et al. Drivers and hotspots of extinction risk in marine mammals. *Proc. Natl. Acad. Sci. USA* **109**, 3395–400 (2012).

20. Perrin, W. F., Würsig, B. G. & Thewissen, J. G. M. *Encyclopedia of marine mammals*. (Elsevier/Academic Press, 2009).

21. Bowen, W. D. Role of marine mammals in aquatic ecosystems. *Mar. Ecol. Prog. Ser.* **158**, 267–274 (1997).

22. Roman, J. & McCarthy, J. J. The Whale Pump: Marine Mammals Enhance Primary Productivity in a Coastal Basin. *PLoS One* **5**, e13255 (2010).

23. Berta, A. *Return to the Sea: The Life and Evolutionary Times of Marine Mammals*. (University of California Press, 2012).

24. Estes, J. A. et al. Trophic downgrading of planet Earth. *Science*. **333**, 301–306 (2011).

25. Laidre, K. L. et al. Quantifying the sensitivity of Arctic marine mammals to climate-induced habitat change. *Ecol. Appl.* **18**, S97–125 (2008).
26. Moore, S. E. & Huntington, H. P. Arctic marine mammals and climate change: Impacts and resilience. *Ecol. Appl.* **18**, S157-S165 (2008).

27. Pompa, S., Ehrlich, P. R. & Ceballos, G. Global distribution and conservation of marine mammals. *Proc. Natl. Acad. Sci.* **108**, 13600-13605 (2011).

28. Albouy, C., Delattre, V. L., Mérigot, B., Meynard, C. N. & Leprieur, F. Multifaceted biodiversity hotspots of marine mammals for conservation priorities. *Divers. Distrib.* **23**, 615–626 (2017).

29. Dalongeville, A., Andrello, M., Mouillot, D., Albouy, C. & Manel, S. Ecological traits shape genetic diversity patterns across the Mediterranean Sea: a quantitative review on fishes. *J. Biogeogr.* **43**, 845–857 (2015).

30. Mouillot, D. et al. Protected and threatened components of fish biodiversity in the mediterranean sea. *Curr. Biol.* **21**, 1044–1050 (2011).

31. Devictor, V. et al. Spatial mismatch and congruence between taxonomic, phylogenetic and functional diversity: the need for integrative conservation strategies in a changing world. *Ecol. Lett.* **13**, 1030–1040 (2010).

32. Guilhaumon, F. et al. Representing taxonomic, phylogenetic and functional diversity: New challenges for Mediterranean marine-protected areas. *Divers. Distrib.* **21** (2015).

33. Cadotte, M. W., Carscadden, K. & Mirotchnick, N. Beyond species: functional diversity and the maintenance of ecological processes and services. *J. Appl. Ecol.* **48**, 1079–1087 (2011).

34. Webb, C. O., Ackerly, D. D., McPeek, M. A. & Donoghue, M. J. Annual Review of Ecology and SystematicsPhylogenies and community ecology. **33**, 475–505 (2002).

35. Diaz, S. & Cabido, M. Vive la difference: plant functional diversity matters to ecosystem processes. *Trends Ecol. Evol.* **16**, 646–655 (2001).

36. Mouillot, D., Graham, N. A. J., Villéger, S., Mason, N. W. H. & Bellwood, D. R. A functional approach reveals community responses to disturbances. *Trends Ecol. Evol.* **28**, 167–177 (2013).

37. Webb, C. O., Ackerly, D. D., McPeek, M. A. & Donoghue, M. J. Phylogenies and community ecology. *Annu. Rev. Ecol. Syst.* **33**, 475–505 (2002).

38. Maherali, H. & Klironomos, J. N. Influence of Phylogeny on fungal community assembly and ecosystem functioning. *Science.* **316**, 1746–1748 (2007).

39. Mazel, F. et al. Prioritizing phylogenetic diversity captures functional diversity unreliably. *Nat. Commun.* **9**, 2888 (2018).

40. Foden, W. B. et al. Identifying the World's Most Climate Change Vulnerable Species: A Systematic Trait-Based Assessment of all Birds, Amphibians and Corals. *PLoS One* **8**, e65427 (2013).

41. Foden, W. B. et al. Climate change vulnerability assessment of species. *Wiley Interdiscip. Rev. Clim. Chang.* **10**, e551 (2019).

42. Kaschner, K., Watson, R., Trites, A. & Pauly, D. Mapping world-wide distributions of marine mammal species using a relative environmental suitability (RES) model. *Mar. Ecol. Prog. Ser.* **316**, 285–310 (2006).

43. MacLeod, C. Global climate change, range changes and potential implications for the conservation of marine cetaceans: a review and synthesis. *Endanger. Species Res.* **7**, 125–136 (2009).

44. Fullard, K. J. et al. Population structure of long-finned pilot whales in the North Atlantic: a correlation with sea surface temperature? *Mol. Ecol.* **9**, 949–958 (2000).

45. Dickinson, M. G., Orme, C. D. L., Suttle, K. B. & Mace, G. M. Separating sensitivity from exposure in assessing extinction risk from climate change. *Sci. Rep.* **4**, 6898 (2015).

46. Böhm, M. et al. Hot and bothered: Using trait-based approaches to assess climate change vulnerability in reptiles. *Biol. Conserv.* **204**, 32–41 (2016).

47. Conti, L., Schmidt-Kloiber, A., Grenouillet, G. & Graf, W. A trait-based approach to assess the vulnerability of European aquatic insects to climate change. *Hydrobiologia* **721**, 297–315 (2014).

48. Hossain, M. A. et al. Assessing the vulnerability of freshwater crayfish to climate change. *Divers. Distrib.* **24**, 1830–1843 (2018).

49. H Rabinowitz, R. A. The Biological Aspects of Rare Plant Conservation Edited by Hugh Syng Seven forms of rarity. In *The Biological Aspects of Rare Plant Conservation* (ed. Hugh Syng) 205–217 (John Wiley & Sons, Ltd, 1981).

50. Harnik, P. G., Simpson, C. & Payne, J. L. Long-term differences in extinction risk among the seven forms of rarity. *Proc. R. Soc. B Biol. Sci.* **279**, 4969–4976 (2012).

51. Slatyer, R. A., Hirst, M. & Sexton, J. P. Niche breadth predicts geographical range size: a general ecological pattern. *Ecol. Lett.* **16**, 1104–1114 (2013).

52. Thuiller, W. Patterns and uncertainties of species' range shifts under climate change. *Glob. Chang. Biol.* **10**, 2020–2027 (2004).

53. Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An overview of CMIP5 and the experiment Design. *Bull. Am. Meteorol. Soc.* **93**, 485–498 (2012).

54. van Vuuren, D. P. et al. The representative concentration pathways: an overview. *Clim. Change* **109**, 5–31 (2011).

55. Moss, R. H. et al. The next generation of scenarios for climate change research and assessment. *Nature* **463**, 747–756 (2010).

56. Collins, M. et al. Long-term climate change: Projections, commitments and irreversibility, in Climate Change The Physical Science Basis. In *Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climat* (ed. Stocker et al. T. F.) (2013).

57. Faith, D. P. J. Conservation evaluation and phylogenetic diversity. *Biol. Conserv.* **61**, 1–10 (1992).

58. Fritz, S. A., Bininda-Emonds, O. R. P. & Purvis, A. Geographical variation in predictors of mammalian extinction risk: big is bad, but only in the tropics. *Ecol. Lett.* **12**, 538–549 (2009).

59. Kissling, W. D. et al. Towards novel approaches to modelling biotic interactions in multispecies assemblages at large spatial extents. *J. Biogeogr.* **39**, 2163–2178 (2012).

60. Villeger, S., Mason, N. W. H. & Mouillot, D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. *Ecology* **89**, 2290–2301 (2008).

61. Legendre, P. & Legendre, L. *Numerical Ecology* (1998).

62. Maire, E., Grenouillet, G., Brosse, S. & Villéger, S. How many dimensions are needed to accurately assess functional diversity? A pragmatic approach for assessing the quality of functional spaces. *Glob. Ecol. Biogeogr.* **24**, 728–740 (2015).

63. Isaac, N. J. B., Turvey, S. T., Collen, B., Waterman, C. & Baillie, J. E. M. Mammals on the EDGE: Conservation Priorities Based on Threat and Phylogeny. *PLoS One* **2**, e296 (2007).

64. Avila, I. C., Kaschner, K. & Dormann, C. F. Current global risks to marine mammals: Taking stock of the threats. *Biol. Conserv.* **221**, 44–58 (2018).

65. IUCN. The IUCN Red List of Threatened Species. Version 2018-2 (2018).

66. González-Suárez, M., Gómez, A. & Revilla, E. Which intrinsic traits predict vulnerability to extinction depends on the actual threatening processes. *Ecosphere* **4**, art76 (2013).

67. Tingley, R., Hitchmough, R. A. & Chapple, D. G. Life-history traits and extrinsic threats determine extinction risk in New Zealand lizards. *Biol. Conserv.* **165**, 62–68 (2013).

68. Vasquez, R. & Simonetti, J. A. Life history traits and sensitivity to landscape change: the case of birds and mammals of mediterranean Chile. *Rev. Chil. Hist. Nat.* **72**, 517–525 (1999).

69. Ivashchenko, Y. & Clapham, P. Soviet catches of right whales Eubalaena japonica and bowhead whales Balaena mysticetus in the North Pacific Ocean and the Okhotsk Sea. *Endanger. Species Res.* **18**, 201–217 (2012).

70. Josephson, E., Smith, T. D. & Reeves, R. R. Historical distribution of right whales in the North Pacific. *Fish Fish.* **9**, 155–168 (2008).
71. Rugh, D. J., Hobbs, R. C., Lerczak, J. A. & Breiwick, J. Estimates of abundance of the eastern North Pacific stock of gray whales (Eschrichtius robustus) 1997–2002. *J. Cetacean Res. Manag.* 7, 1–12 (2005).

72. Roman, J. et al. Whales as marine ecosystem engineers. *Front. Ecol. Environ.* **12**, 377–385 (2014).

73. Alter, S. E. et al. Climate impacts on transocean dispersal and habitat in gray whales from the Pleistocene to 2100. *Mol. Ecol.* **24**, 1510–1522 (2015).

74. McLeod, E., Salm, R., Green, A. & Almany, J. Designing marine protected area networks to address the impacts of climate change. *Front. Ecol. Environ.* **7**, 362–370 (2009).

75. Frainer, A. et al. Climate-driven changes in functional biogeography of Arctic marine fish communities. *Proc. Natl. Acad. Sci. USA* **114**, 12202–12207 (2017).

76. Rosset, V. & Oertli, B. Freshwater biodiversity under climate warming pressure: Identifying the winners and losers in temperate standing waterbodies. *Biol. Conserv.* **144**, 2311–2319 (2011).

77. Thuiller, W. et al. A road map for integrating eco-evolutionary processes into biodiversity models. *Ecol. Lett.* **16**, 94–105 (2013).

78. Gregr, E. J. Insights into North Pacific right whale Eubalaena Japonica habitat from historic whaling records. *Endanger. Species Res.* **15**, 223–239 (2011).

79. Barbet-Massin, M. & Jetz, W. The effect of range changes on the functional turnover, structure and diversity of bird assemblages under future climate scenarios. *Glob. Chang. Biol.* **21**, 2917–2928 (2015).

80. Albouy, C. et al. Projected impacts of climate warming on the functional and phylogenetic components of coastal Mediterranean fish biodiversity. *Ecography.* **38**, 681–689 (2015).

81. Pimiento, C. et al. The Pliocene marine megafauna extinction and its impact on functional diversity. *Nat. Ecol. Evol.* **1**, 1100–1106 (2017).

82. Tulloch, V. J. D., Plagányi, É. E., Brown, C., Richardson, A. J. & Matear, R. Future recovery of baleen whales is imperiled by climate change. *Glob. Chang. Biol.* **25**, 1263–1281 (2019).

83. Maslowski, W., Clement Kinney, J., Higgins, M. & Roberts, A. The Future of Arctic Sea Ice. *Annu. Rev. Earth Planet. Sci.* **40**, 625–654 (2012).

84. Bronselaer, B. et al. Change in future climate due to Antarctic meltwater. *Nature* **564**, 53–58 (2018).

## Acknowledgements

This work was funded by the French National Research Agency (ANR MORSE Project 11-CEPL-006). We thank Bastien Mérigot and Christine Meynard for the helpful comments during preliminary analyses. TLF acknowledges funding from the Swiss National Science Foundation under Grant PP00P2_170687 and the European Union's Horizon 2020 research and innovation programme under grant agreement No 820989.

## Author contributions

C.A., V.D. and F.L. conceived the idea for the study. C.A., V.D., S.A.B. and T.F. collected the data. C.A. and V.D. performed the analyses with S.A.B., D.M., L.P., G.D. and M.R. providing input. C.A. and F.L. led the manuscript writing, and V.D., S.A.B., D.M., L.P., G.D., M.R. and T.F. edited the paper.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-019-57280-3.

Correspondence and requests for materials should be addressed to C.A.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020