Next Article in Journal
The Performance and Metabolism of Dairy Cows Receiving an Ultra-Diluted Complex in the Diet during the Transition Period and Early Lactation
Next Article in Special Issue
The Science Underlying Giant Panda Conservation Translocations
Previous Article in Journal
Associations between Surface and Rectal Temperature Profiles of Low-Birth-Weight Piglets
Previous Article in Special Issue
Initial Post-Release Performance of Cultured Cyprinus chilia Juveniles in a Shallow Lake in Southwestern China
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 13
Issue 20
10.3390/ani13203260
animals-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating the Spatial Distribution and Future Conservation Requirements of the Spotted Seal in the North Pacific

by
Leyu Yang
1,2,
Hongfei Zhuang
2,
Shenghao Liu
2,
Bailin Cong
2,
Wenhao Huang
2,3,
Tingting Li
2,
Kaiyu Liu
1,2 and
Linlin Zhao
2,*
1
School of Advanced Manufacturing, Fuzhou University, Jinjiang 362200, China
2
Key Laboratory of Marine Eco-Environmental Science and Technology, First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China
3
College of Environmental Science and Engineering, Ocean University of China, Qingdao 266100, China
*
Author to whom correspondence should be addressed.
Animals 2023, 13(20), 3260; https://doi.org/10.3390/ani13203260
Submission received: 6 September 2023 / Revised: 4 October 2023 / Accepted: 7 October 2023 / Published: 19 October 2023
(This article belongs to the Special Issue Protecting Endangered Species)
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

To understand the impact of climate change on marine mammals, we focused on the spotted seal population in the North Pacific. This ice-breeding species exhibits distinct variations across different regions. Our study aimed to quantify their ecological niches and conduct a conservation gap analysis. We found clear niche divergence among three populations and observed habitat contraction driven by climate change, potentially leading to breeding habitat loss in certain areas. Unfortunately, existing marine protected areas do not adequately cover most spotted seal habitats. By incorporating local adaptation into species distribution modeling, our research provides valuable insights for designing effective conservation policies to protect the different geographical populations of spotted seals in the face of climate change. This study highlights the importance of considering local adaptation in conservation and management strategies for marine mammal species.

Abstract

Local adaptation has been increasingly involved in the designation of species conservation strategies to response to climate change. Marine mammals, as apex predators, are climatechange sensitive, and their spatial distribution and conservation requirements are critically significant for designing protection strategies. In this study, we focused on an ice-breeding marine mammal, the spotted seal (Phoca largha), which exhibits distinct morphological and genetic variations across its range. Our objectives were to quantify the ecological niches of three spotted seal populations, construct the species-level model and population-level models that represent different regions in the Bering population (BDPS), Okhotsk population (ODPS) and southern population (SDPS), and conduct a conservation gap analysis. Our findings unequivocally demonstrated a clear niche divergence among the three populations. We predicted habitat contraction for the BDPS and ODPS driven by climate change; in particular, the spotted seals inhabiting Liaodong Bay may face breeding habitat loss. However, most spotted seal habitats are not represented in existing marine protected areas. Drawing upon these outcomes, we propose appropriate conservation policies to effectively protect the habitat of the different geographical populations of spotted seals. Our research addresses the importance of incorporating local adaptation into species distribution modeling to inform conservation and management strategies.

1. Introduction

Climate change poses a significant peril to global biodiversity in the 21st century, unleashing a series of profound and unpredictable changes on our planet [1]. In comparison to terrestrial communities, marine communities stand as more vulnerable sentinels, acutely attuned to the environmental changes wrought by this climatic change [2]. The effects of climate change on the marine environment have changed the life history and distribution landscape of marine species, with the potential to cause habitat destruction and even species extinction [3]. Indeed, mounting evidence attests to the rapid reconfiguration of species distributions along depth or latitudinal gradients in response to changing climates [4,5]. Marine mammals, serving as highly susceptible apex predators within marine ecosystems, assume a pivotal role in maintaining biodiversity and regulating ecosystem processes [6]. Consequently, obtaining comprehensive and precise knowledge of the current distributions of marine mammals, as well as reliable projections under future climate change scenarios, is of paramount importance in developing effective strategies for resource management and conservation.
Species distribution models (SDMs) are useful tools in this endeavor, as they can predict potential changes in species habitats by determining statistical relationships between species occurrence data and environmental predictors, and can also be used to forecast how suitable areas may vary under different climate change scenarios [7]. Traditionally, SDMs have been constructed at the species level based on the “niche conservatism” hypothesis, which suggests that individuals of the same species have similar niche spaces and exhibit consistent responses to climate change across their range [8,9]. Nevertheless, over an extensive evolutionary timeframe, species with a wide range may gradually adapt to local conditions, leading to niche divergence [10]. Recent research highlights a growing recognition of the importance of integrating local adaptation or intraspecific variation into climate responses, with an increasing number of studies emphasizing this crucial aspect [11,12]. By modeling habitat suitability below the species level, more accurate estimations of species ranges and climate change projections can be attained. Notably, within the realm of species conservation, the consideration of climate change responses within phylogeographic lineages has gained currency for certain taxonomic groups in terrestrial communities [13,14], and marine communities [15,16].
Here, we examined how predictions of climate change responses can differ when considering local adaptation for spotted seals. The International Union for Conservation of Nature’s Red List of Threatened Species ranked spotted seals (Phoca largha) as “least concern”, indicating a low risk of extinction [17]. However, in China, spotted seals have the highest protection level among rare and endangered species (class-I protection) due to increasing threats from climate change and habitat loss [18]. North Pacific spotted seals have eight specialized breeding habitats and limited mobility, which contribute to their vulnerability to climate change [19]. Based on morphological and genetic differences, spotted seals were divided into three distinct geographical populations [20]: the Bering population (BDPS), the Okhotsk population (ODPS) and the Southern population (SDPS). Throughout their extensive evolutionary narratives, the three populations have taken up residence within divergent ecological domains, thereby hinting at the plausible presence of localized adaptation. While the precise degree of distributional overlap and genetic interchange amongst these populations remains shrouded in obscurity, the likelihood of their existence cannot be discounted [21,22]. Neglecting to account for such local adaptations when employing SDMs in conservation or management decisions may result in erroneous characterizations of species’ responses to environmental changes throughout their ranges, thus misdirecting conservation efforts [23]. Therefore, to further estimate the impacts of climate change on spotted seals at a higher taxonomic resolution below the species level has become an urgent issue in their long-term conservation and management.
In this study, we quantified realized niches (i.e., the portion of the fundamental niche currently used by the species), developed SDMs and made future predictions to examine how climate change might influence spotted seals by constructing species-level versus population-level models. We sought to address the following hypotheses: (1) conspicuous disparities exist in spatial niches among the three distinct populations; (2) the three populations respond differently to climate change, and population-level SDMs are more reliable than species-level results; (3) the potential distribution of spotted seals under future climatic conditions will shift to higher latitudes; and (4) there are large gaps in spotted seal conservation outside protected areas that are not effectively protected. Our results emphasize the value of developing SDMs below the species level and serve as a useful guide for designing climate-adapted conservation and management strategies for spotted seals within more precise taxonomic units.

2. Materials and Methods

2.1. Study Area and Spotted Seal Occurrence Data

This study focused on the distribution range of spotted seals in the North Pacific, with the main study area located in the temperate and cold temperate coastal and littoral regions (90° E—240° W, 0° N—80° N; Figure 1). The spotted seal occurrence data were collected from the literature, the Global Biodiversity Information Facility (GBIF; http://www.gbif.org/, accessed on 21 December 2022) [24], and the Ocean Biogeographic Information System (https://obis.org/, accessed on 4 December 2022) [25] (Supplementary Table S1). To minimize sampling bias, we spatially thinned the occurrence data to match the resolution of the environmental data (5′ × 5′, approximately 9.2 km × 9.2 km) using the R package “spThin” [26], and only one random spotted seal distribution record was used in each raster. Following this data-filtering procedure, 1990 records were kept in order to construct the SDM at the species level (hereafter “species model”). Of these, 680 records belonged to the BDPS that was used to construct the SDM at the population level (hereafter “BDPS model”), 853 records belonged to the ODPS that was used to construct the ODPS model, and 457 records belonged to the SDPS that was used to construct the SDPS model.

2.2. Environmental Predictor Variables

Habitat surroundings have a significant influence on the distribution of spotted seals, and considering a combination of bioenvironmental relevance and data availability, 10 environmental variables that may influence the distribution of spotted seals were selected for this study (Table 1). Current and future environmental data were downloaded from online datasets: the water depth and distance to shore were downloaded from the Global Marine Environment Datasets (https://gmed.auckland.ac.nz/, accessed on 15 April 2023) [27], and the remaining predictors were downloaded from the Ocean Raster for Analysis of Climate and Environment (https://bio-oracle.org/, accessed on 6 April 2023) [28]. Considering the correlation between predictors, we completed Pearson’s correlation factor analysis between environmental layers using the R package usdm [29], retaining only environmental variables with correlation factor values < |0.7| [30] (Supplementary Figure S1). Finally, seven predictors, including water depth (Dep), distance to shore (DTS), chlorophyll concentration (Chl), current velocity (CV), ice thickness (IT), salinity, and water temperature (Tmean), were retained for modeling analysis.
To project the future habitat suitability of spotted seals, we considered four representative concentration pathway (RCP) scenarios (i.e., RCP 2.6, RCP 4.5, RCP 6.0 and RCP 8.5), and two periods (i.e., 2050s: the average for 2040s–2050s, and 2100s: the average for 2090–2100). We obtained the corresponding projections of future marine environmental layers from Bio-ORACLE v2.0. This study assumed no change in water depth and distance to shore [31].

2.3. Estimates of Niche Divergence

To assess whether the three spotted seal populations occupy different niche spaces, we used n-dimensional hypervolume to characterize their realized niche [32]. For this, we first performed principal component analysis (PCA) on seven selected environmental variables and retained the top four principal components, which cumulatively explained 83.6% of the total variance (Supplementary Figure S2). Then, we calculated the principal component retention values corresponding to each population using a Gaussian approach with the R package hypervolume [33]. Finally, the niche divergence between populations can be assessed by overlapping the hypervolume of each population using the R package BAT [34]. Total niche divergence (βTotal) was divided into the following two processes: niche shift (spatial replacement between hypervolumes) and niche contraction/expansion (net difference between hypervolumes). The βTotal ranged from 0 to 1, indicating the two hypervolumes of complete overlap to complete separation [11].

2.4. SDMs Establishment and Projection

We conducted SDM analysis based on the “biomod 2” package in the R platform (version 4.2.3) [35]. This package contained 10 modeling algorithms: generalized additive mode, generalized linear model (GLM), generalized boosting model (GBM), random fores, surface range envelope (SRE), artificial neural network (ANN), flexible discriminant analysis (FDA), classification tree analysis (CTA), multiple adaptive regression splines (MARS), and maximum entropy (Maxent). Since true absence data were lacking, we randomly simulated the same number of pseudo-absence records as that of presence records in the environmental conditions [36]. The dataset was divided into five groups during the modeling process, with an equal number of records in each group, four of which were used for model training and the remaining one for model testing. To evaluate the predictive performance of each model, the 5-fold cross-validation process was repeated 10 times. We used the TSS (the true skill statistic) and AUC (the area under the ROC curve) values to assess the accuracy of the models.
This study selected TSS > 0.8 and AUC > 0.9 as model selection standards [37,38] and used a weighted-average algorithm to build an integrated model for reducing the uncertainty of individual models. To better explain habitat suitability, we transformed continuous habitat suitability predictions into a binary map by maximizing the probability threshold of the TSS [39]. We applied a randomized method to measure Pearson correlations between all predictor and assessment variables [40] to assess the relative importance of each variable in predicting species distributions. Finally, we built two levels of species and population ensemble models to predict the potential distribution of habitat for the whole species and three geographical populations (the BDPS, ODPS and SDPS) under current and future (2050s, 2100s) climate scenarios under RCP 2.6 and RCP 8.5.

2.5. Protection Gap Analysis

The Global Marine Protected Areas layers were sourced from the World Database on Protected Areas (https://www.protectedplanet.net/, accessed on 4 May 2023) [41], while data on protected areas in the Yellow and Bohai Seas of China were sourced from a published article [42]. First, we overlaid the layers of existing protected areas and the range of spotted seals to analyze the proportion of the existing protected area covered in the spotted seal distribution range and the uncovered spatial area in QGIS 3.28.6 software (https://www.qgis.org/en/site/, accessed on 5 April 2023). We then conducted a conservation gap analysis of the uncovered spatial areas to identify uncovered habitat areas. Finally, by integrating current and future climate change scenarios, we projected trends in the range of spotted seals and further identified conservation gaps for these species under future climate scenarios.

3. Results

3.1. Niche Divergence among the Three Populations

According to the results of niche divergence studies of different geographical populations, the BDPS has the widest ecological range. The four-dimensional hypervolume for the BDPS, ODPS, and SDPS, respectively, was 1642.20, 353.46, and 145.16. The niche divergence between two populations was very high, with values of 0.81 (BDPS:ODPS), 0.92 (BDPS:SDPS), and 0.86 (ODPS:SDPS) shown by the paired comparison of hypervolumes. Contraction/expansion accounted for more than 85% of the niche divergence between the BDPS and the ODPS or SDPS, whereas niche transitions had a much smaller role (15%). The main cause of the niche divergence between the ODPS and SDPS was contraction/expansion (>65%), with niche shift accounting for a little part (35%) (Table 2).
When the three populations’ general niches were compared, it was clear that PCA1 was the primary focus of niche divergence (Figure 2a), which was mostly explained by water depth, distance to shore, and chlorophyll concentration (Supplementary Figure S3). The paired-niche comparison revealed that the PCA1 was primarily responsible for the niche divergence between the BDPS and ODPS (Figure 2b), which was primarily explained by water depth, ice thickness, and current velocity (Supplementary Figure S3); the PCA1 was responsible for the niche divergence between the BDPS and SDPS (Figure 2c), which was primarily explained by water temperature, chlorophyll concentration, and salinity (Supplementary Figure S3); and the PCA1 was also responsible for the niche divergence between the BDPS and SDPS (Figure 2d), which was primarily explained by ice thickness, chlorophyll concentration and water temperature (Supplementary Figure S3).

3.2. Current SDMs Projections

Based on the TSS and AUC values in the model results, eight models were selected to build the weighted species-level ensemble model after removing MaxEnt and SRE, and nine models were selected to build the weighted population-level ensemble model after removing SRE from the ten single models. The higher values of AUC and TSS for all four ensemble models indicated high predictive performance (Table 3).
The species-level model showed that water temperature and depth considerably contributed to the distribution of spotted seals, while current velocity and chlorophyll concentration contributed little to the model (Figure 3a). According to this model, spotted seals preferred to live in waters with a temperature range of 0 °C to 15 °C and a depth of 0 m to 1000 m (Supplementary Figure S4). The population-level model showed that the main factors influencing the potential distribution of spotted seals differ between populations. Specifically, the BDPS distribution was most influenced by water temperature and ice thickness; this population preferred to inhabit areas with water temperature ranging from 0 °C to 8 °C (Supplementary Figure S5a) and had a higher probability of occurrence at ice thicknesses of 0–1 m (Supplementary Figure S5b). The ODPS distribution was most influenced by water temperature and depth; this population preferred to inhabit areas with water temperature ranging from 2 °C to 12 °C (Supplementary Figure S5c) and water depth ranging from 0 m to 2000 m (Supplementary Figure S5d). The SDPS distribution was mainly influenced by water temperature and distance to shore (Figure 3b); this population preferred to live in a water temperature range of 2 °C to 15 °C (Supplementary Figure S5e) and had a higher probability of living in near-shore waters not more than 800 meters from the coast (Supplementary Figure S5f).
According to the results of the model, the distribution and suitable habitat of the three spotted seal populations under current environmental conditions are as follows: the BDPS is mainly distributed in the Chukchi Sea, the Bering Sea, the coast of Kamchatka Island and the northern part of the Sea of Okhotsk, with the largest area of suitable habitat (6.24 × 106 km2) (Figure 4c,f). The ODPS is mainly distributed in the Sea of Okhotsk and near Sakhalin Island, extending northward to the Aleutian Peninsula, with the second largest area of suitable habitat (3.54 × 106 km2) (Figure 4d,g). The SDPS is mainly distributed on Hokkaido Island, the Sea of Japan and the Yellow Sea of China, reaching as far as the East China Sea, with the smallest area of suitable habitat (1.08 × 106 km2) (Figure 4e,h).
The species model results are shown in Figure 4a,b. Both the continuous and binary projections showed that the predicted suitable area from the species model is similar to that from the BDPS model in the Chukchi Sea and the Bering Sea, is similar to that from the ODPS model in the Okhotsk Sea and is smaller than that from the SDPS model in the southernmost region.

3.3. Habitat Suitability under Future Climate Scenarios

This study built SDMs at both the species and population levels to predict climate change impacts on potentially suitable habitats for spotted seals. Habitat-scale changes are influenced by climate change scenarios, particularly under the pessimistic scenario of uncontrolled greenhouse gas emissions (RCP 8.5), and the suitable habitat is predicted to vary considerably within the appropriate range (Table 4). For all climate change scenarios, both the BDPS and ODPS show a decreasing trend in the extent of suitable habitat, with the greatest decrease under RCP 8.5 in the 2050s, while the SDPS shows an increasing trend in the extent of suitable habitat. For the BDPS, suitable habitat for spotted seals is mainly stable in the Chukchi Sea and around the Bering Strait, with the northern coasts of the Chukchi and Taymyr Peninsulas also serving as potential habitats (Figure 5c,f). For the ODPS, suitable habitat is mainly stable along the northwestern coast of the Kamchatka Peninsula, while most of the suitable habitat in the Sea of Okhotsk will be lost (Figure 5d,g). For the SDPS, although large areas of suitable habitat are found around the Aleutian Islands and the Gulf of Alaska, the southernmost spotted seal colony in Liaodong Bay will be lost (Figure 5e,h).
At the species level, the southernmost Liaodong Bay spotted seal breeding area is gradually being lost; as the range of the Arctic Ocean seas north of the Chukchi Peninsula expands toward higher latitudes, the distribution of spotted seals may occur at higher latitudes in the Thamel Peninsula (Figure 5a,b). Range changes in the size of the predicted species tended to decrease under RCP 2.6 and increase under RCP 8.5, with the change reaching 36.16% under RCP 8.5 in the 2100s.

3.4. Spotted Seal Conservation Gap Analysis

According to the results of the protection gap analysis, 278,617 km2 was found to be protected, representing only 5.65% of the spotted seal range, i.e., more than 94% of the area is not covered by protected areas (Figure 6a). In the overlaid map, the Eastern Bering Sea, the Aleutian Islands and the Sea of Okhotsk overlap significantly with the range of spotted seals. Although the number of protected areas around the Yellow and Bohai Seas in China is high, they are very small. Further analysis suggests that outside these existing marine protected areas (MPAs), there are both important nonbreeding areas and breeding habitat for spotted seals, such as the western coast of the Kamchatka Peninsula, the western Bering Sea and Peter the Great Bay, which remain unprotected.
In this study, the potential range of spotted seals under future climate scenarios was overlaid with existing MPAs, and it was found that spotted seals were protected in only 6.04–6.22% of their range under the RCP 2.6 and RCP 8.5 scenarios in the 2050s and 2100s (Table 5). In summary, the overall change in the protected gap area is insignificant in the four future climate scenarios (Figure 6b,c).

4. Discussion

4.1. Consideration of Local Adaptation

The study quantified the realized niches of three spotted seal populations and found that the BDPS exhibited a substantially broader niche space than the ODPS and SDPS, primarily due to niche contraction/expansion. These findings underscore that geographically separated spotted seal populations inhabit dissimilar ecological niches; therefore, niche conservatism does not apply to this species. Then, we constructed SDMs for the BDPS, ODPS, and SDPS individually, which accounted for local adaptation, and these models revealed that the three populations showed differing responses to climate change predictors. In summary, this study emphasizes the importance of considering local adaptation in projecting the potential distribution of species to inform conservation and management decisions in a climate change scenario.
Local adaptation and intraspecific variation were incorporated into SDMs based on the recognition that populations of species inhabiting widely different habitats over significant time scales will often show adaptations to their respective local conditions, resulting in intraspecific niche variation. This local adaptation may be reflected by morphological and genetic differences [43]. For instance, studies have revealed significant differences in the nonmetric features of the skull between spotted seals from the central and eastern Bering Sea [44] and a phylogeographic break between spotted seals breeding in the Yellow Sea–Japan Sea region and those breeding in the Sea of Okhotsk, Bering Sea and Chukchi Sea. The above studies, in conjunction with the observed niche divergence among the BDPS, ODPS, and SDPS established in this study, stress the importance of building SDMs at the population level to account for local adaptation.
In this study, despite the fact that the species and population level SDMs predicted a similar change trend in the species range, with a reduction in suitable habitat for the BDPS and ODPS and an expansion of suitable habitat for the SDPS under future climate change scenarios, the magnitudes of the range change predicted by the two types of models varied. The population level model produced more encouraging findings for the BDPS and ODPS with less loss of appropriate habitat. Due to the inclusion of potential local responses in population-level models, our climate change estimates were in fact less pessimistic. Our results are in line with a number of published studies that suggest adaptive genetic variation within a species can reduce the species’ susceptibility to climate change [45,46,47].

4.2. Impacts of Climate Change on Spotted Seals

The niche divergence of the BDPS, ODPS and SDPS was mainly due to niche contraction/expansion, to which water temperature, ice thickness and chlorophyll concentration contributed most. Spotted seals are cool-temperature marine mammals and water temperature is an important influence on both the physiology and behavior of spotted seals, affecting them indirectly by altering the distribution of prey, predators and disease-causing vectors [48]. Also, spotted seals are often dependent on sea ice for breeding and foraging, and thicker ice may provide better breeding and foraging conditions, while thinner ice may limit these activities. Moreover, chlorophyll concentration variations may reflect the productivity of marine ecosystems and the base of the food chain [49]. Differences in nutrient conditions and chlorophyll concentration in different regions may lead to different availability of food resources, and spotted seal populations may choose to adapt to different food resources according to chlorophyll concentration in different regions, resulting in population differentiation. Climate-induced changes in these factors cause species to alter their current distribution patterns to track ecological niches. In summary, these three climatic factors directly or indirectly affect the physiology, behavior and availability of food resources of spotted seals, thus leading to the divergence of adaptation strategies among different populations.
Different populations of spotted seals may have different levels of adaptiveness and vulnerability to climate change. Our population-level SDMs predicted that the BDPS will colonize the northern coast of the Chukchi, and as far west as the Taymyr Peninsulas, while retaining most of the current suitable areas, indicating the resilience of the population to climate change. In contrast, the ODPS preserves most of the currently suitable areas, and parts of the suitable habitat in the Sea of Okhotsk will be lost, indicating that the population is less resilient to climate change. The SDPS will have large areas of suitable habitat near the Aleutian Islands and the Gulf of Alaska, but spotted seals of the SDPS would not be capable of shifting their range to the northeast to reach this area and the southernmost spotted seal breeding area in Liaodong Bay will be lost, indicating the vulnerability of this population to climate change.
Overall, the thinning and breaking up of ice caused by climate change will expose vast regions of the northern Bering Sea and Chukchi Sea, which is likely to increase suitable habitat for northern spotted seals, thus offsetting the loss of habitat in the south. Similar findings were reported in studies of climate change-induced geographic translocation of species [15,50,51,52].

4.3. Model Predictive Accuracy

Integrated habitat suitability models built with a weighted integration technique, especially for rare species, can increase model prediction accuracy and avoid overfitting problems without compromising explanatory power [53]. The TSS values of the four ensemble models were all 0.86, their AUC values were all 0.95, and their errors for the environmental importance results were all less than 0.2, showing that the predictions were very accurate representations of the current and future distributions of spotted seals under various climate scenarios. An integrated model was used in this work to estimate the whole distributions of spotted seals, and the results were mostly in line with the known range of spotted seals (temperate and cold temperate coastal and shore).
Although the results were as expected, there were limitations to the method of model species distribution to predict the range of spotted seals. We used future environmental data that did not occur objectively but were predictions based on atmosphere-ocean general circulation models (AOGCMs), as well as the effects of many factors such as changes in food availability, changes in ocean currents that persist between breeding sites, and the ability of spotted seals to swim on their own. When predicting the suitable distribution of spotted seals, the results of the study may have overestimated the range of the spotted seal. Predictions made in the future that take these elements into account in the model will be closer to the species’ actual distribution. Moreover, to improve or validate the predictive power of SDMs, independent geographically or temporally separated data should be collected [54]. Considering the difficulty and cost of field surveys, emerging environmental DNA methods could be used to determine the presence of potential distribution areas for spotted seals as predicted in this study, and these data could be used in future work on SDMs.

4.4. Management and Conservation of Spotted Seals

Spotted seals have received little attention from conservationists or managers despite significant human and climate change challenges. Spotted seals and other marine mammals are vulnerable in the face of global warming, and their potential extinction could have far-reaching consequences for the functioning of global marine ecosystems in the future [55]. Water temperature and ice thickness are two important environmental factors that might affect the distribution of spotted seal populations geographically, yet they are also strongly related to climate change. To develop climate-adapted conservation and management methods, it is imperative to assess how the changing climate is affecting the appropriateness of spotted seal habitat. In general, all three populations exhibit a propensity to migrate toward the poles under a warming climate. However, the SDPS are more vulnerable to climate than the BDPS and ODPS, and their population sizes have already been significantly reduced from historical levels and may be at risk of population genetic extinction. Therefore, the SDPS deserves more attention and protection in the face of climate stress.
Marine protected areas (MPAs) have proven to be an effective tool for protecting endangered species and maintaining ecosystem services [56]. The European Union target of ‘30 by 30’ that is 30% of the ocean protected (as MPAs) by 2030. For the conservation of spotted seals and their habitat, many countries and regions have established MPAs [57,58]. These marine protected areas were crucial for the conservation of spotted seals and their habitats. However, when we overlaid the spotted seal range, we discovered that only 5.65–6.22% of the range was protected, meaning that more than 94% of the area was unprotected. Therefore, protected areas (marine reserves, nature reserves and national parks) need to be expanded, and the establishment of protected areas across international borders should be considered to better protect spotted seals and their habitats [59].
The current and potential future distribution of spotted seals is mainly in coastal waters; as a result, protecting the species from anthropogenic environmental contamination and hunting pressure in these seas is essential. Given the different vulnerabilities of the three populations to climate change, we need to develop local adaptive conservation and management measures for the populations in different areas. Since the SDPS is most sensitive to climate change, we recommend that stricter hunting restrictions, such as bans on poaching for genitalia and culling by fishers, be imposed on this population. For ODPS, although there is a degree of adaptability to climate change and little risk of population extinction, direct or indirect commercial fisheries interactions may have a significant cumulative effect. Therefore, we should increase management efforts to control marine development activities such as sand mining, oil and gas exploration and water pipelines to reduce damage to spotted seals (survival environment). Although the BDPS has the ability to adapt to climate change, we should also pay attention to the timing and routes of the breeding migrations of spotted seals, avoid fishing and shipping operations during migrations, and raise awareness of conservation.
Finally, we must emphasize that MPAs and other suggested conservation measures will only guarantee that adequate habitats for spotted seals are safeguarded from human impact now and in the future. However, as greenhouse gas concentrations rise, the SDPS will lose its ideal habitat. This population will progressively decline towards extinction if the issue is not managed for the long-term. Therefore, reducing human-caused greenhouse gas emissions is the ultimate solution for the sustainability of these populations.

5. Conclusions

In conclusion, this study represents the first step in estimating climate impacts on the potential distribution of spotted seals in the North Pacific considering local adaptation. Population-level SDMs are more reliable than species-level SDMs because of the different responses of the three spotted seal populations to environmental predictor variables. Additionally, conservation efforts should be dedicated to the establishment of MPAs, first in the stable spots predicted to remain climatically suitable for the species, and second in the currently suitable areas. The comparison of current and predicted habitat suitability maps presented in our study serves as a crucial tool allowing us to delineate the most promising regions for establishing both types of measures. In future studies, other analytical methods and multiple data sources should be incorporated to improve our ability to predict the potential distribution of spotted seals and deliver more accurate information for related conservation and management.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani13203260/s1, Figure S1: Results of Pearson correlation analysis for 10 environmental variables; Figure S2: The results of principal component analyses based on all the environmental variables involved in the model using the R package FactoMineR; Figure S3: The degree of contribution of environmental variables to the principal component that has the greatest influence on ecological niche differentiation between populations; Figure S4: The response curves of spotted seals probability of occurrence against water temperature and water depth based on the species-level model; Figure S5: The response curves of population probability of occurrence against the two most important driving factors based on the population-level model. Table S1: Reference list of occurrence data for spotted seal (Phoca largha) in this research [21,60,61,62,63,64,65,66].

Author Contributions

Conceptualization, L.Y. and L.Z.; methodology, L.Y.; software, L.Y.; validation, S.L., B.C. and W.H.; formal analysis, T.L.; investigation, K.L.; resources, L.Z.; data curation, S.L.; writing—original draft preparation, L.Y.; writing—review and editing, H.Z.; visualization, B.C.; supervision, L.Z.; project administration, L.Z.; funding acquisition, L.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (2022YFF0802204) and the National Natural Science Foundation of China (grant no. 32201433).

Institutional Review Board Statement

Ethical review and approval were waived for this study due to all data sources for spotted seals were downloaded online.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

We thank the modelling help from Linjie Li and Bin Sun in Ocean University of China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dawson, T.P.; Jackson, S.T.; House, J.I. Beyond Predictions: Biodiversity Conservation in a Changing Climate. Science 2011, 332, 53–58. [Google Scholar] [CrossRef] [PubMed]
  2. Sorte, C.J.B.; Williams, S.L.; Carlton, J.T. Marine range shifts and species introductions: Comparative spread rates and community impacts. Glob. Ecol. Biogeogr. 2010, 19, 303–316. [Google Scholar] [CrossRef]
  3. Cheung, W.W.L.; Sarmiento, J.L.; Dunne, J.; Frölicher, T.L.; Lam, V.W.Y.; Deng Palomares, M.L.; Watson, R.; Pauly, D. Shrinking of fishes exacerbates impacts of global ocean changes on marine ecosystems. Nat. Clim. Chang. 2012, 3, 254–258. [Google Scholar] [CrossRef]
  4. Hastings, R.A.; Rutterford, L.A.; Freer, J.J.; Collins, R.A.; Simpson, S.D.; Genner, M.J. Climate Change Drives Poleward Increases and Equatorward Declines in Marine Species. Curr. Biol. 2020, 30, 1572–1579. [Google Scholar] [CrossRef]
  5. Fossheim, M.; Primicerio, R.; Johannesen, E.; Ingvaldsen, R.B.; Aschan, M.M.; Dolgov, A.V. Recent warming leads to a rapid borealization of fish communities in the Arctic. Nat. Clim. Chang. 2015, 5, 673–677. [Google Scholar] [CrossRef]
  6. Tian, J.; Lu, Z.; Sanganyado, E.; Gan, Z.; Wang, Z.; Kong, Z.; Wu, J.; Liu, W. Tissue distribution and trophic magnification of trace elements in typical marine mammals in Bohai and north Yellow Seas. Mar. Pollut. Bull. 2022, 184, 114135. [Google Scholar] [CrossRef] [PubMed]
  7. Pacifici, M.; Foden, W.B.; Visconti, P. Assessing species vulnerability to climate change. Nat. Clim. Chang. 2015, 5, 215–224. [Google Scholar] [CrossRef]
  8. Smith, A.B.; Godsoe, W.; Rodriguez-Sanchez, F.; Wang, H.H.; Warren, D. Niche Estimation Above and Below the Species Level. Trends. Ecol. Evol. 2019, 34, 260–273. [Google Scholar] [CrossRef] [PubMed]
  9. Pearson, R.G.; Dawson, T.P. Predicting the impacts of climate change on the distribution of species: Are bioclimate envelope models useful? Glob. Ecol. Biogeogr. 2003, 12, 361–371. [Google Scholar] [CrossRef]
  10. Lee-Yaw, J.A.; Kharouba, H.M.; Bontrager, M.; Mahony, C.; Csergo, A.M.; Noreen, A.M.; Li, Q.; Schuster, R.; Angert, A.L. A synthesis of transplant experiments and ecological niche models suggests that range limits are often niche limits. Ecol. Lett. 2016, 19, 710–722. [Google Scholar] [CrossRef]
  11. Carvalho, J.C.; Cardoso, P. Decomposing the Causes for Niche Differentiation Between Species Using Hypervolumes. Front. Ecol. Evol. 2020, 8, 00243. [Google Scholar] [CrossRef]
  12. Peterson, M.L.; Doak, D.F.; Morris, W.F. Incorporating local adaptation into forecasts of species’ distribution and abundance under climate change. Glob. Chang. Biol. 2019, 25, 775–793. [Google Scholar]
  13. Zhao, X.; Li, X.; Zhang, Z.; Garber, P.A.; Yu, M.; Qiao, H.; Li, M. Differential response to climate change and human activities in three lineages of Sichuan snub-nosed monkeys (Rhinopithecus roxellana). Divers. Distrib. 2022, 28, 2416–2428. [Google Scholar] [CrossRef]
  14. Zhang, P.; Grenouillet, G.; Dong, X.; Zheng, Y.; Lek, S.; Chang, J. Capturing response differences of species distribution to climate and human pressures by incorporating local adaptation: Implications for the conservation of a critically endangered species. J. Environ. Manag. 2021, 284, 111998. [Google Scholar] [CrossRef] [PubMed]
  15. Li, L.; Zhao, L.; Fu, J.; Sun, B.; Liu, C. Predicting the habitat suitability for populations of Pacific cod under different climate change scenarios considering intraspecific genetic variation. Ecol. Indic. 2022, 142, 109248. [Google Scholar] [CrossRef]
  16. Hu, Z.M.; Zhang, Q.S.; Zhang, J.; Kass, J.M.; Mammola, S.; Fresia, P.; Draisma, S.G.A.; Assis, J.; Jueterbock, A.; Yokota, M.; et al. Intraspecific genetic variation matters when predicting seagrass distribution under climate change. Mol. Ecol. 2021, 30, 3840–3855. [Google Scholar] [CrossRef] [PubMed]
  17. Boveng, P.L. Phoca Largha. The IUCN Red List of Threatened Species 2016: e.T17023A45229806. 2016. Available online: https://www.iucnredlist.org/species/17023/45229806 (accessed on 6 September 2023).
  18. Wang, F.; Yu, Y.; Wu, H.; Wu, W.; Wang, L.; An, L.; Cai, W. Microplastics in spotted seal cubs (Phoca largha): Digestion after ingestion? Sci. Total Environ. 2021, 785, 147426. [Google Scholar] [CrossRef]
  19. Shaughnessy, P.D.; Fay, F.H. A review of the taxonomy and nomenclature of North Pacific Harbour seals. J. Zool. 1977, 182, 385–419. [Google Scholar] [CrossRef]
  20. Boveng, P.L.; Bengtson, J.L.; Buckley, T.W.; Cameron, M.F.; Dahle, S.P.; Kelly, B.P.; Megrey, B.A.; Overland, J.E.; Williamson, N.J. Status Review of the Spotted Seal (Phoca largha); Department of Commerce, National Oceanic and Atmospheric Administration, National Marine Fisheries Service: Silver Spring, MD, USA, 2009; 153p, NMFS-AFSC-200. [Google Scholar]
  21. Trukhin, A.M. Changes in the Breeding Range of Spotted Seals, Phoca largha (Pallas, 1811) (Carnivora, Pinnipedia), in the Western Sea of Japan: Causes and Effects. Russ. J. Mar. Biol. 2022, 48, 101–107. [Google Scholar] [CrossRef]
  22. Podlesnykh, A.V.; Katin, I.O. Genetic variation in the spotted seal (Phoca largha Pallas, 1811) from the Rimsky-Korsakov Archipelago (Peter the Great Bay, western sea of Japan) as inferred from mitochondrial DNA control region sequences. Zool. Anz. 2022, 301, 174–178. [Google Scholar] [CrossRef]
  23. Hamann, A.; Aitken, S.N. Conservation planning under climate change: Accounting for adaptive potential and migration capacity in species distribution models. Divers. Distrib. 2013, 19, 268–280. [Google Scholar] [CrossRef]
  24. The Global Biodiversity Information Facility (GBIF): An On-Line, Free and Open Access to Biodiversity Data. Available online: https://www.gbif.org/occurrence/download/0229487-220831081235567 (accessed on 21 December 2022).
  25. The Ocean Biogeographic Information System (OBIS): An On-Line, Worldwide Atlas for Accessing, Modeling and Mapping Marine Biological Data in a Multidimensional Geographic Context. Available online: https://datasets.obis.org/downloads/540d223c-5ef1-458e-b271-6f786e8d1d8b.zip (accessed on 4 December 2022).
  26. Aiello-Lammens, M.E.; Boria, R.A.; Radosavljevic, A.; Vilela, B.; Anderson, R.P. spThin: An R package for spatial thinning of species occurrence records for use in ecological niche models. Ecography 2015, 38, 541–545. [Google Scholar] [CrossRef]
  27. Basher, Z.; Bowden, D.; Costello, M.J. Global Marine Environment Datasets (GMED). World Wide Web Electronic Publication. Version 2.0 (Rev.02.2018). Available online: http://gmed.auckland.ac.nz (accessed on 15 April 2023).
  28. Assis, J.; Tyberghein, L.; Bosch, S.; Verbruggen, H.; Serrão, E.A.; De Clerck, O.; Tittensor, D. Bio-ORACLE v2.0: Extending Marine Data Layers for Bioclimatic Modelling. Available online: https://bio-oracle.org/ (accessed on 6 April 2023).
  29. Naimi, B.; Hamm, N.A.; Groen, T.A.; Skidmore, A.K.; Toxopeus, A.G. Where is positional uncertainty a problem for species distribution modelling? Ecography 2014, 37, 191–203. [Google Scholar] [CrossRef]
  30. Fournier, A.; Barbet-Massin, M.; Rome, Q.; Courchamp, F. Predicting species distribution combining multi-scale drivers. Glob. Ecol. Conserv. 2017, 12, 215–226. [Google Scholar] [CrossRef]
  31. Zhang, Z.; Kass, J.M.; Mammola, S. Lineage-level distribution models lead to more realistic climate change. Divers. Distrib. 2021, 27, 684–695. [Google Scholar] [CrossRef]
  32. Hutchinson, G.E. Concluding Remarks. Cold Spring Harb. Symp. Quant. Biol. 1957, 22, 415–427. [Google Scholar] [CrossRef]
  33. Blonder, B.; Morrow, C.B.; Maitner, B.; Harris, D.J.; Lamanna, C.; Violle, C.; Enquist, B.J.; Kerkhoff, A.J. New approaches for delineating n-dimensional hypervolumes. Methods Ecol. Evol. 2018, 9, 305–319. [Google Scholar] [CrossRef]
  34. Cardoso, P.; Rigal, F.; Carvalho, J.C. BAT–Biodiversity Assessment Tools, an R package for the measurement and estimation of alpha and beta taxon, phylogenetic and functional diversity. Methods Ecol. Evol. 2015, 6, 232–236. [Google Scholar] [CrossRef]
  35. Thuiller, W.; Georges, D.; Engler, R.; Breiner, F.; Georges, M.D.; Thuiller, C.W. biomod2: Ensemble Platform for Species Dis-tribution Modeling. Available online: https://cran.r-project.org/web/packages/biomod2/index.html (accessed on 6 April 2023).
  36. Barbet-Massin, M.; Jiguet, F.; Albert, C.H.; Thuiller, W. Selecting pseudo-absences for species distribution models: How, where and how many? Methods Ecol. Evol. 2012, 3, 327–338. [Google Scholar] [CrossRef]
  37. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  38. Swets, J.A. Measuring the Accuracy of Diagnostic Systems. Science 1988, 240, 1285–1293. [Google Scholar] [CrossRef]
  39. Liu, C.; White, M.; Newell, G.; Pearson, R. Selecting thresholds for the prediction of species occurrence with presence-only data. J. Biogeogr. 2013, 40, 778–789. [Google Scholar] [CrossRef]
  40. Zanardo, N.; Parra, G.J.; Passadore, C.; Möller, L.M. Ensemble modelling of southern Australian bottlenose dolphin Tursiops sp. distribution reveals important habitats and their potential ecological function. Mar. Ecol. Prog. Ser. 2017, 569, 253–266. [Google Scholar] [CrossRef]
  41. UNEP-WCMC and IUCN (2023), Protected Planet: The World Database on Protected Areas (WDPA) and World Database on Other Effective Area-Based Conservation Measures (WD-OECM). Available online: https://www.protectedplanet.net (accessed on 4 May 2023).
  42. Zhuang, H.; Shao, F.; Zhang, C.; Xia, W.; Wang, S.; Qu, F.; Wang, Z.; Lu, Z.; Zhao, L.; Zhang, Z. Spatial-temporal shifting patterns and in situ conservation of spotted seal (Phoca largha) populations in the Yellow Sea ecoregion. Integr. Zool. 2023. [CrossRef] [PubMed]
  43. Rice, D.W. A List of the Marine Mammals of the World; Department of Commerce, National Oceanic and Atmospheric Admin-istration, National Marine Fisheries Service: Silver Spring, MD, USA, 1977; 13p, NMFS-SSRF-711. [Google Scholar]
  44. Fedoseev, G. Use of Nonmetrical Characters of Skulls of Bering Sea Seals in a Study of the Phenotypic Structure of Their Populations; Department of Commerce, National Oceanic and Atmospheric Administration, National Marine Fisheries Service: Silver Spring, MD, USA, 1984; pp. 49–54, NMFS 12. [Google Scholar]
  45. Razgour, O.; Forester, B.; Taggart, J.B.; Bekaert, M.; Juste, J.; Ibanez, C.; Puechmaille, S.J.; Novella-Fernandez, R.; Alberdi, A.; Manel, S. Considering adaptive genetic variation in climate change vulnerability assessment reduces species range loss projections. Proc. Natl. Acad. Sci. USA 2019, 116, 10418–10423. [Google Scholar] [CrossRef] [PubMed]
  46. Ikeda, D.H.; Max, T.L.; Allan, G.J.; Lau, M.K.; Shuster, S.M.; Whitham, T.G. Genetically informed ecological niche models improve climate change predictions. Glob. Chang. Biol. 2017, 23, 164–176. [Google Scholar] [CrossRef]
  47. Oney, B.; Reineking, B.; O’Neill, G.; Kreyling, J. Intraspecific variation buffers projected climate change impacts on Pinus contorta. Ecol. Evol. 2013, 3, 437–449. [Google Scholar] [CrossRef]
  48. Grebmeier, J.M.; Overland, J.E.; Moore, S.E.; Farley, E.V.; Carmack, E.C.; Cooper, L.W.; Frey, K.E.; Helle, J.H.; McLaughlin, F.A.; McNutt, S.L. A major ecosystem shift in the northern Bering Sea. Science 2006, 311, 1461–1464. [Google Scholar] [CrossRef]
  49. Polovina, J.J.; Howell, E.A.; Kobayashi, D.R.; Seki, M.P. The transition zone chlorophyll front updated: Advances from a decade of research. Prog. Oceanogr. 2017, 150, 79–85. [Google Scholar] [CrossRef]
  50. Yang, T.; Liu, X.; Han, Z. Predicting the Effects of Climate Change on the Suitable Habitat of Japanese Spanish Mackerel (Scomberomorus niphonius) Based on the Species Distribution Model. Front. Mar. Sci. 2022, 9, 927790. [Google Scholar] [CrossRef]
  51. Sun, B.; Zhao, L.; Shao, F.; Lu, Z.; Tian, J.; Liu, C. Estimating the impacts of climate change on the habitat suitability of common minke whales integrating local adaptation. Front. Mar. Sci. 2022, 9, 923205. [Google Scholar] [CrossRef]
  52. Chen, Y.; Shan, X.; Ovando, D.; Yang, T.; Dai, F.; Jin, X. Predicting current and future global distribution of black rockfish (Sebastes schlegelii) under changing climate. Ecol. Indic. 2021, 128, 107799. [Google Scholar] [CrossRef]
  53. Breiner, F.T.; Guisan, A.; Bergamini, A.; Nobis, M.P.; Anderson, B. Overcoming limitations of modelling rare species by using ensembles of small models. Methods Ecol. Evol. 2015, 6, 1210–1218. [Google Scholar] [CrossRef]
  54. Bahn, V.; McGill, B.J. Testing the predictive performance of distribution models. Oikos 2013, 122, 321–331. [Google Scholar] [CrossRef]
  55. Albouy, C.; Delattre, V.; Donati, G.; Frölicher, T.L.; Albouy-Boyer, S.; Rufino, M.; Pellissier, L.; Mouillot, D.; Leprieur, F. Global vulnerability of marine mammals to global warming. Sci. Rep. 2020, 10, 548. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, W.; Li, J. conservation of biodiversity in China: Advances and prospects. Biodivers. Sci. 2021, 29, 133–149. [Google Scholar] [CrossRef]
  57. Ma, Z.; Chen, Y.; Melville, D.S.; Fan, J.; Liu, J.; Dong, J.; Tan, K.; Cheng, X.; Fuller, R.A.; Xiao, X.; et al. Changes in area and number of nature reserves in China. Conserv. Biol. 2019, 33, 1066–1075. [Google Scholar] [CrossRef] [PubMed]
  58. Witherell, D.; Woodby, D. Application of Marine Protected Areas for Sustainable Production and Marine. Mar. Fish. Rev. 2005, 67, 1–28. [Google Scholar]
  59. Hamilton, C.D.; Lydersen, C.; Aars, J.; Acquarone, M.; Atwood, T.; Baylis, A.; Biuw, M.; Boltunov, A.N.; Born, E.W.; Boveng, P.; et al. Marine mammal hotspots across the circumpolar Arctic. Divers. Distrib. 2022, 28, 2729–2753. [Google Scholar] [CrossRef]
  60. Trukhin, A.M.; Permyakov, P.A.; Ryazanov, S.D.; Lobanov, V.B.; Kim, H.W.; Choi, Y.M.; Sohn, H. Migrations of young spotted seals (Phoca largha) from Peter the Great Bay, Sea of Japan/East Sea, and the pattern of their use of seasonal habitats. PLoS ONE 2021, 16, e0244232. [Google Scholar]
  61. Solovyeva, M.A.; Kuznetsova, D.M.; Glazov, D.M.; Burkanov, V.N.; Rozhnov, V.V. The First Data on the Use of the Pacific Ocean Area by the Okhotsk Population of the Spotted Seal Phoca largha. Biol. Bull. 2021, 48 (Suppl. S1), S219–S224. [Google Scholar] [CrossRef]
  62. Solov’eva, M.A.; Glazov, D.M.; Kuznetsova, D.M.; Rozhnov, V.V. Movements of spotted seals (Phoca largha) in the Sea of Okhotsk according to satellite tagging data. Russ. J. Ecol. 2016, 47, 405–411. [Google Scholar] [CrossRef]
  63. Han, J.; Lu, Z.; Tian, J.; Ma, Z.; Wang, Z.; Yang, Y.; Wang, Q.; Song, X.; Peng, Z. Release studies on spotted seals (Phoca largha) using satellite telemetry tracking technique. Acta Theriol. Sin. 2013, 33, 300. [Google Scholar]
  64. Mizuno, A.W.; Wada, A.; Ishinazaka, T.; Hattori, K.; Watanabe, Y.; Ohtaishi, N. Distribution and abundance of spotted seals Phoca largha and ribbon seals Phoca fasciata in the southern Sea of Okhotsk. Ecol. Res. 2002, 17, 79–96. [Google Scholar] [CrossRef]
  65. Lowry, L.F.; Burkanov, V.N.; Frost, K.J.; Simpkins, M.A.; Davis, R.; DeMaster, D.P.; Suydam, R.; Springer, A. Habitat use and habitat selection by spotted seals (Phoca largha) in the Bering Sea. Can. J. Zool. 2000, 78, 1959–1971. [Google Scholar] [CrossRef]
  66. Lowry, L.F.; Frost, K.J.; Davis, R.; DeMaster, D.P.; Suydam, R.S. Movements and behavior of satellite-tagged spotted seals (Phoca largha) in the Bering and Chukchi Seas. Polar Biol. 1998, 19, 221–230. [Google Scholar] [CrossRef]
Animals 13 03260 g001
Figure 1. Map of the study area and occurrence records of spotted seal. Red dots represent the occurrence records of the Bering population (BDPS), green dots represent those of the Okhotsk population (ODPS), and blue dots those of the southern population (SDPS). The dotted black lines are drawn along 43° N latitude and 157° E longitude, which are considered the boundaries between the SDPS and ODPS and between the ODPS and BDPS.
Figure 1. Map of the study area and occurrence records of spotted seal. Red dots represent the occurrence records of the Bering population (BDPS), green dots represent those of the Okhotsk population (ODPS), and blue dots those of the southern population (SDPS). The dotted black lines are drawn along 43° N latitude and 157° E longitude, which are considered the boundaries between the SDPS and ODPS and between the ODPS and BDPS.
Animals 13 03260 g001
Animals 13 03260 g002
Figure 2. Three spotted seal populations with realized niches quantified by four-dimensional hypervolume. The larger blue, orange and green dots indicate the mean niche position (niche centroid) of the Bering population (BDPS), Okhotsk population (ODPS), and southern population (SDPS), respectively. The overlapping hypervolumes for the three populations (a), BDPS and ODPS (b), BDPS and SDPS (c) and ODPS and SDPS (d).
Figure 2. Three spotted seal populations with realized niches quantified by four-dimensional hypervolume. The larger blue, orange and green dots indicate the mean niche position (niche centroid) of the Bering population (BDPS), Okhotsk population (ODPS), and southern population (SDPS), respectively. The overlapping hypervolumes for the three populations (a), BDPS and ODPS (b), BDPS and SDPS (c) and ODPS and SDPS (d).
Animals 13 03260 g002
Animals 13 03260 g003
Figure 3. Importance of environmental variables in driving whole-species distribution based on the species-level model (a) and population distributions based on the population-level model (b). Dep, water depth; DTS, distance to shore; Tmean, water temperature; Sal, salinity; IT, ice thickness; CV, current velocity.
Figure 3. Importance of environmental variables in driving whole-species distribution based on the species-level model (a) and population distributions based on the population-level model (b). Dep, water depth; DTS, distance to shore; Tmean, water temperature; Sal, salinity; IT, ice thickness; CV, current velocity.
Animals 13 03260 g003
Animals 13 03260 g004
Figure 4. Habitat suitability map of spotted seals projected by the ensemble models under current climate scenarios. (a,b) Respective binary and continuous plots of the species; (c,f) respective binary and continuous plots for BDPS; (d,g) corresponding plots for ODPS; (e,h) corresponding plots for SDPS.
Figure 4. Habitat suitability map of spotted seals projected by the ensemble models under current climate scenarios. (a,b) Respective binary and continuous plots of the species; (c,f) respective binary and continuous plots for BDPS; (d,g) corresponding plots for ODPS; (e,h) corresponding plots for SDPS.
Animals 13 03260 g004
Animals 13 03260 g005
Figure 5. Predicted change in the suitable habitat based on the ensemble models in the 2050s under RCP 8.5 (5085) and in the 2100s under RCP 8.5 (0085) for the species (a,b), BDPS (c,f), ODPS (d,g) and SDPS (e,h).
Figure 5. Predicted change in the suitable habitat based on the ensemble models in the 2050s under RCP 8.5 (5085) and in the 2100s under RCP 8.5 (0085) for the species (a,b), BDPS (c,f), ODPS (d,g) and SDPS (e,h).
Animals 13 03260 g005
Animals 13 03260 g006
Figure 6. Analysis of the conservation gap for spotted seals under current and future climate scenarios. (a) Analysis of the conservation gap under the climate scenario; (b) analysis of the conservation gap in the 2050s under RCP 8.5 (5085); (c) analysis of the conservation gap in the 2100s under RCP 8.5 (0085).
Figure 6. Analysis of the conservation gap for spotted seals under current and future climate scenarios. (a) Analysis of the conservation gap under the climate scenario; (b) analysis of the conservation gap in the 2050s under RCP 8.5 (5085); (c) analysis of the conservation gap in the 2100s under RCP 8.5 (0085).
Animals 13 03260 g006
Table 1. Ten environmental variables initially selected for this study and their units, spatial resolution and sources.
Table 1. Ten environmental variables initially selected for this study and their units, spatial resolution and sources.
Environment VariableUnitSpatial ResolutionSource
water depthm5 arc minuteshttps://gmed.aucklandac.nz/, accessed on 15 April 2023
distance to shorekm5 arc minuteshttps://gmed.auckland.ac.nz/, accessed on 15 April 2023
calcitemol·m−35 arc minuteshttps://bio-oracle.org/, accessed on 6 April 2023
chlorophyll concentrationmg·m−35 arc minuteshttps://bio-oracle.org/, accessed on 6 April 2023
currents velocitym·s−15 arc minuteshttps://bio-oracle.org/, accessed on 6 April 2023
dissolved oxygenmol·m−35 arc minuteshttps://bio-oracle.org/, accessed on 6 April 2023
sea ice concentrationfraction5 arc minuteshttps://bio-oracle.org/, accessed on 6 April 2023
ice thicknessm5 arc minuteshttps://bio-oracle.org/, accessed on 6 April 2023
salinityPSS5 arc minuteshttps://bio-oracle.org/, accessed on 6 April 2023
water temperature°C5 arc minuteshttps://bio-oracle.org/, accessed on 6 April 2023
Table 2. Total niche differentiation (βTotal) between populations and the proportion of niche shift and niche contraction/expansion. BDPS, Bering distinct population segment; ODPS, Okhotsk distinct population segment; SDPS, Southern distinct population segment.
Table 2. Total niche differentiation (βTotal) between populations and the proportion of niche shift and niche contraction/expansion. BDPS, Bering distinct population segment; ODPS, Okhotsk distinct population segment; SDPS, Southern distinct population segment.
Populations PairβTotalNiche ShiftNiche Contraction/
Expansion
BDPS-ODPS0.810.12(15%)0.69(85%)
BDPS-SDPS0.920.04(4%)0.88(96%)
ODPS-SDPS0.860.27(32%)0.59(68%)
Table 3. Number of models and evaluating indicators for the ensemble models built at the species level and population level. NME, number of models used in ensemble modeling; TSS, the true skill statistics; AUC, the area under the receiver operating characteristic curve.
Table 3. Number of models and evaluating indicators for the ensemble models built at the species level and population level. NME, number of models used in ensemble modeling; TSS, the true skill statistics; AUC, the area under the receiver operating characteristic curve.
EnsembleTSSAUCNME
Species model0.8570.9538
BDPS model0.8610.9499
ODPS model0.9320.9759
SDPS model0.9500.9789
Table 4. Predicted size of changes [9] in species range based on the species-level and population-level models under future climate scenarios. RCP 2.6 (8.5), the representative concentration pathway 2.6 (8.5); 2050s (2100s), at the middle (end) of the 21st century.
Table 4. Predicted size of changes [9] in species range based on the species-level and population-level models under future climate scenarios. RCP 2.6 (8.5), the representative concentration pathway 2.6 (8.5); 2050s (2100s), at the middle (end) of the 21st century.
RCPBDPSODPSSDPSSpecies
2050s2100s2050s2100s2050s2100s2050s2100s
RCP 2.6−32.48−38.34−32.80−47.6113.1517.86−1.43−1.73
RCP 8.5−38.24−63.94−44.25−66.519.5762.912.2736.16
Table 5. Spotted seal areas (km2) and corresponding percentages protected under the current and future climate scenarios.
Table 5. Spotted seal areas (km2) and corresponding percentages protected under the current and future climate scenarios.
Climate ScenarioArea Protected (km2)Percentage of Protection [9]
current278,6175.65
2050s RCP 2.6278,0626.17
2050s RCP 8.5277,7766.04
2100s RCP 2.6279,4036.22
2100s RCP 8.5277,4495.56
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, L.; Zhuang, H.; Liu, S.; Cong, B.; Huang, W.; Li, T.; Liu, K.; Zhao, L. Estimating the Spatial Distribution and Future Conservation Requirements of the Spotted Seal in the North Pacific. Animals 2023, 13, 3260. https://doi.org/10.3390/ani13203260

AMA Style

Yang L, Zhuang H, Liu S, Cong B, Huang W, Li T, Liu K, Zhao L. Estimating the Spatial Distribution and Future Conservation Requirements of the Spotted Seal in the North Pacific. Animals. 2023; 13(20):3260. https://doi.org/10.3390/ani13203260

Chicago/Turabian Style

Yang, Leyu, Hongfei Zhuang, Shenghao Liu, Bailin Cong, Wenhao Huang, Tingting Li, Kaiyu Liu, and Linlin Zhao. 2023. "Estimating the Spatial Distribution and Future Conservation Requirements of the Spotted Seal in the North Pacific" Animals 13, no. 20: 3260. https://doi.org/10.3390/ani13203260

APA Style

Yang, L., Zhuang, H., Liu, S., Cong, B., Huang, W., Li, T., Liu, K., & Zhao, L. (2023). Estimating the Spatial Distribution and Future Conservation Requirements of the Spotted Seal in the North Pacific. Animals, 13(20), 3260. https://doi.org/10.3390/ani13203260

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop