Range contraction and increasing isolation of a polar bear subpopulation in an era of sea-ice loss

Impact of Holocene environmental change on the evolutionary ecology of an Arctic top predator

The Arctic is among the most climatically sensitive environments on Earth, and the disappearance of multiyear sea ice in the Arctic Ocean is predicted within decades. As apex predators, polar bears are sentinel species for addressing the impact of environmental variability on Arctic marine ecosystems. By integrating genomics, isotopic analysis, morphometrics, and ecological modeling, we investigate how Holocene environmental changes affected polar bears around Greenland. We uncover reductions in effective population size coinciding with increases in annual mean sea surface temperature, reduction in sea ice cover, declines in suitable habitat, and shifts in suitable habitat northward. Furthermore, we show that west and east Greenlandic polar bears are morphologically, and ecologically distinct, putatively driven by regional biotic and genetic differences. Together, we provide insights into the vulnerability of polar bears to environmental change and how the Arctic marine ecosystem plays a vital role in shaping the evolutionary and ecological trajectories of its inhabitants.

Drivers of polar bear behavior and the possible effects of prey availability on foraging strategy

BACKGROUND

Change in behavior is one of the earliest responses to variation in habitat suitability. It is therefore important to understand the conditions that promote different behaviors, particularly in areas undergoing environmental change. Animal movement is tightly linked to behavior and remote tracking can be used to study ethology when direct observation is not possible.

METHODS

We used movement data from 14 polar bears (Ursus maritimus) in Hudson Bay, Canada, during the foraging season (January-June), when bears inhabit the sea ice. We developed an error-tolerant method to correct for sea ice drift in tracking data. Next, we used hidden Markov models with movement and orientation relative to wind to study three behaviors (stationary, area-restricted search, and olfactory search) and examine effects of 11 covariates on behavior.

RESULTS

Polar bears spent approximately 47% of their time in the stationary drift state, 29% in olfactory search, and 24% in area-restricted search. High energy behaviors occurred later in the day (around 20:00) compared to other populations. Second, olfactory search increased as the season progressed, which may reflect a shift in foraging strategy from still-hunting to active search linked to a shift in seal availability (i.e., increase in haul-outs from winter to the spring pupping and molting seasons). Last, we found spatial patterns of distribution linked to season, ice concentration, and bear age that may be tied to habitat quality and competitive exclusion.

CONCLUSIONS

Our observations were generally consistent with predictions of the marginal value theorem, and differences between our findings and other populations could be explained by regional or temporal variation in resource availability. Our novel movement analyses and finding can help identify periods, regions, and conditions of critical habitat.

Movements of southern elephant seals (Mirounga leonina) from Davis Base, Antarctica: combining population genetics and tracking data

Marine animals such as the southern elephant seal (Mirounga leonina) rely on a productive marine environment and are vulnerable to oceanic changes that can affect their reproduction and survival rates. Davis Base, Antarctica, acts as a moulting site for southern elephant seals that forage in Prydz Bay, but the mitochondrial haplotype diversity and natal source populations of these seals have not been characterized. In this study, we combined genetic and animal tracking data on these moulting seals to identify levels of mitochondrial haplotype diversity, natal source population, and movement behaviours during foraging and haul-out periods. Using partial sequences of the mitochondrial control region, we identified two major breeding mitochondrial lineages of seals at Davis Base. We found that the majority of the seals originated from breeding stocks within the South Atlantic Ocean and South Indian Ocean. One seal was grouped with the Macquarie Island breeding stock (South Pacific Ocean). The Macquarie Island population, unlike the other two stocks, is decreasing in size. Tracking data revealed long-distance foraging activity of the Macquarie Island seal around Crozet Islands. We speculate that changes to the Antarctic marine environment can result in a shift in foraging and movement strategies, which subsequently affects seal population growth rates.

Migratory vertebrates shift migration timing and distributions in a warming Arctic

Climate warming in the Arctic has led to warmer and earlier springs, and as a result, many food resources for migratory animals become available earlier in the season, as well as become distributed further northwards. To optimally profit from these resources, migratory animals are expected to arrive earlier in the Arctic, as well as shift their own spatial distributions northwards. Here, we review literature to assess whether Arctic migratory birds and mammals already show shifts in migration timing or distribution in response to the warming climate. Distribution shifts were most prominent in marine mammals, as expected from observed northward shifts of their resources. At least for many bird species, the ability to shift distributions is likely constrained by available habitat further north. Shifts in timing have been shown in many species of terrestrial birds and ungulates, as well as for polar bears. Within species, we found strong variation in shifts in timing and distributions between populations. Ou r review thus shows that many migratory animals display shifts in migration timing and spatial distribution in reaction to a warming Arctic. Importantly, we identify large knowledge gaps especially concerning distribution shifts and timing of autumn migration, especially for marine mammals. Our understanding of how migratory animals respond to climate change appears to be mostly limited by the lack of long-term monitoring studies.

An on-ice aerial survey of the Kane Basin polar bear (Ursus maritimus) subpopulation

There is an imminent need to collect information on distribution and abundance of polar bears (Ursus maritimus) to understand how they are affected by the ongoing decrease in Arctic sea ice. The Kane Basin (KB) subpopulation is a group of high-latitude polar bears that ranges between High Arctic Canada and NW Greenland around and north of the North Water polynya (NOW). We conducted a line transect distance sampling aerial survey of KB polar bears during 28 April–12 May 2014. A total of 4160 linear kilometers were flown in a helicopter over fast ice in the fjords and over offshore pack ice between 76° 50′ and 80° N′. Using a mark-recapture distance sampling protocol, the estimated abundance was 190 bears (95% lognormal CI: 87–411; CV 39%). This estimate is likely negatively biased to an unknown degree because the offshore sectors of the NOW with much open water were not surveyed because of logistical and safety reasons. Our study demonstrated that aerial surveys may be a feasible method for obtaining abundance estimates for small subpopulations of polar bears.

A novel mark-recapture-recovery survey using genetic sampling for polar bears Ursus maritimus in Baffin Bay

Changes in sea-ice dynamics are affecting polar bears Ursus maritimus across their circumpolar range, which highlights the importance of periodic demographic assessments to inform management and conservation. We used genetic mark-recapture-recovery to derive estimates of abundance and survival for the Baffin Bay (BB) polar bear subpopulation—the first time this method has been used successfully for this species. Genetic data from tissue samples we collected via biopsy darting were combined with historical physical capture and harvest recovery data. The combined data set consisted of 1410 genetic samples (2011-2013), 914 physical captures (1993-1995, 1997), and 234 harvest returns of marked bears (1993-2013). The estimate of mean subpopulation abundance was 2826 (95% CI = 2284-3367) in 2012-2013. Estimates of annual survival (mean ± SE) were 0.90 ± 0.05 and 0.78 ± 0.06 for females and males age ≥2 yr, respectively. The proportion of total mortality of adult females and males that was attributed to legal harvest was 0.16 ± 0.05 and 0.26 ± 0.06, respectively. Remote sensing sea-ice data, telemetry data, and spatial distribution of onshore sampling indicated that polar bears were more likely to use offshore sea-ice habitat during the 1990s sampling period compared to the 2010s. Furthermore, in the 1990s, sampling of deep fjords and inland areas was limited, and no offshore sampling occurred in either time period, which precluded comparisons of abundance between the 1993-1997 and 2011-2013 study periods. Our findings demonstrate that genetic sampling can be a practical method for demographic assessment of polar bears over large spatial and temporal scales.

Sea ice reduction drives genetic differentiation among Barents Sea polar bears

Loss of Arctic sea ice owing to climate change is predicted to reduce both genetic diversity and gene flow in ice-dependent species, with potentially negative consequences for their long-term viability. Here, we tested for the population-genetic impacts of reduced sea ice cover on the polar bear (Ursus maritimus) sampled across two decades (1995-2016) from the Svalbard Archipelago, Norway, an area that is affected by rapid sea ice loss in the Arctic Barents Sea. We analysed genetic variation at 22 microsatellite loci for 626 polar bears from four sampling areas within the archipelago. Our results revealed a 3-10% loss of genetic diversity across the study period, accompanied by a near 200% increase in genetic differentiation across regions. These effects may best be explained by a decrease in gene flow caused by habitat fragmentation owing to the loss of sea ice coverage, resulting in increased inbreeding of local polar bears within the focal sampling areas in the Svalbard Archipelago. This study illustrates the importance of genetic monitoring for developing adaptive management strategies for polar bears and other ice-dependent species.

Dietary ecology of Alaskan polar bears (Ursus maritimus) through time and in response to Arctic climate change

Arctic climate change poses serious threats to polar bears (Ursus maritimus) as reduced sea ice makes seal prey inaccessible and marine ecosystems undergo bottom-up reorganization. Polar bears' elongated skulls and reduced molar dentition, as compared to their sister species the grizzly bear (Ursus arctos), are adaptations associated with hunting seals on sea ice and a soft, lipid-rich diet of blubber and meat. With significant declines in sea ice, it is unclear if and how polar bears may be altering their diets. Clarifying polar bear dietary responses to changing climates, both today and in the past, is critical to proper conservation and management of this apex predator. This is particularly important when a dietary strategy may be maladaptive. Here, we test the hypothesis that hard-food consumption (i.e., less preferred foods including bone), inferred from dental microwear texture analysis, increased with Arctic warming. We find that polar bears demonstrate a conserved absence of hard-object feeding in Alaska through time (including approximately 1000 years ago), until the 21st century, consistent with a highly conserved and specialized diet of soft blubber and flesh. Notably, our results also suggest that some 21st-century polar bears may be consuming harder foods (e.g., increased carcass utilization, terrestrial foods including garbage), despite having skulls and metabolisms poorly suited for such a diet. Prior to the 21st century, only polar bears with larger mandibles demonstrated increased hard-object feeding, though to a much lower degree than closely related grizzly bears which regularly consume mechanically challenging foods. Polar bears, being morphologically specialized, have biomechanical constraints which may limit their ability to consume mechanically challenging diets, with dietary shifts occurring only under the most extreme scenarios. Collectively, the highly specialized diets and cranial morphology of polar bears may severely limit their ability to adapt to a warming Arctic.

Marine mammal conservation: over the horizon

Marine mammals can play important ecological roles in aquatic ecosystems, and their presence can be key to community structure and function. Consequently, marine mammals are often considered indicators of ecosystem health and flagship species. Yet, historical population declines caused by exploitation, and additional current threats, such as climate change, fisheries bycatch, pollution and maritime development, continue to impact many marine mammal species, and at least 25% are classified as threatened (Critically Endangered, Endangered or Vulnerable) on the IUCN Red List. Conversely, some species have experienced population increases/recoveries in recent decades, reflecting management interventions, and are heralded as conservation successes. To continue these successes and reverse the downward trajectories of at-risk species, it is necessary to evaluate the threats faced by marine mammals and the conservation mechanisms available to address them. Additionally, there is a need to identify evidence-based priorities of both research and conservation needs across a range of settings and taxa. To that effect we: (1) outline the key threats to marine mammals and their impacts, identify the associated knowledge gaps and recommend actions needed; (2) discuss the merits and downfalls of established and emerging conservation mechanisms; (3) outline the application of research and monitoring techniques; and (4) highlight particular taxa/populations that are in urgent need of focus.

Transient benefits of climate change for a high-Arctic polar bear (Ursus maritimus) subpopulation

Kane Basin (KB) is one of the world's most northerly polar bear (Ursus maritimus) subpopulations, where bears have historically inhabited a mix of thick multiyear and annual sea ice year-round. Currently, KB is transitioning to a seasonally ice-free region because of climate change. This ecological shift has been hypothesized to benefit polar bears in the near-term due to thinner ice with increased biological production, although this has not been demonstrated empirically. We assess sea-ice changes in KB together with changes in polar bear movements, seasonal ranges, body condition, and reproductive metrics obtained from capture-recapture (physical and genetic) and satellite telemetry studies during two study periods (1993-1997 and 2012-2016). The annual cycle of sea-ice habitat in KB shifted from a year-round ice platform (~50% coverage in summer) in the 1990s to nearly complete melt-out in summer (<5% coverage) in the 2010s. The mean duration between sea-ice retreat and advance increased from 109 to 160 days (p = .004). Between the 1990s and 2010s, adult female (AF) seasonal ranges more than doubled in spring and summer and were significantly larger in all months. Body condition scores improved for all ages and both sexes. Mean litter sizes of cubs-of-the-year (C0s) and yearlings (C1s), and the number of C1s per AF, did not change between decades. The date of spring sea-ice retreat in the previous year was positively correlated with C1 litter size, suggesting smaller litters following years with earlier sea-ice breakup. Our study provides evidence for range expansion, improved body condition, and stable reproductive performance in the KB polar bear subpopulation. These changes, together with a likely increasing subpopulation abundance, may reflect the shift from thick, multiyear ice to thinner, seasonal ice with higher biological productivity. The duration of these benefits is unknown because, under unmitigated climate change, continued sea-ice loss is expected to eventually have negative demographic and ecological effects on all polar bears.

Comparing sea ice habitat fragmentation metrics using integrated step selection analysis

Habitat fragmentation occurs when continuous habitat gets broken up as a result of ecosystem change. While commonly studied in terrestrial ecosystems, Arctic sea ice ecosystems also experience fragmentation, but are rarely studied in this context. Most fragmentation analyses are conducted using patch-based metrics, which are potentially less suitable for sea ice that has gradual changes between sea ice cover, than distinct "long-term" patches. Using an integrated step selection analysis, we compared the descriptive power of a patch-based metric to a more novel metric, the variation in local spatial autocorrelation over time. We used satellite telemetry data from 39 adult female polar bears (Ursus maritimus) in Hudson Bay to examine their sea ice habitat using Advanced Microwave Scanning Radiometer 2 data during sea ice breakup in May through July from 2013-2018. Spatial autocorrelation resulted in better model fits across 64% of individuals, although both metrics were more effective in describing movement patterns than habitat selection. Variation in local spatial autocorrelation allows for the visualization of sea ice habitat at complex spatial and temporal scales, condensing a targeted time period of habitat that would otherwise have to be analyzed daily.

Interrelated ecological impacts of climate change on an apex predator

Climate change has broad ecological implications for species that rely on sensitive habitats. For some top predators, loss of habitat is expected to lead to cascading behavioral, nutritional, and reproductive changes that ultimately accelerate population declines. In the case of the polar bear (Ursus maritimus), declining Arctic sea ice reduces access to prey and lengthens seasonal fasting periods. We used a novel combination of physical capture, biopsy darting, and visual aerial observation data to project reproductive performance for polar bears by linking sea ice loss to changes in habitat use, body condition (i.e., fatness), and cub production. Satellite telemetry data from 43 (1991-1997) and 38 (2009-2015) adult female polar bears in the Baffin Bay subpopulation showed that bears now spend an additional 30 d on land (90 d in total) in the 2000s compared to the 1990s, a change closely correlated with changes in spring sea ice breakup and fall sea ice formation. Body condition declined for all sex, age, and reproductive classes and was positively correlated with sea ice availability in the current and previous year. Furthermore, cub litter size was positively correlated with maternal condition and spring breakup date (i.e., later breakup leading to larger litters), and negatively correlated with the duration of the ice-free period (i.e., longer ice-free periods leading to smaller litters). Based on these relationships, we projected reproductive performance three polar bear generations into the future (approximately 35 yr). Results indicate that two-cub litters, previously the norm, could largely disappear from Baffin Bay as sea ice loss continues. Our findings demonstrate how concurrent analysis of multiple data types collected over long periods from polar bears can provide a mechanistic understanding of the ecological implications of climate change. This information is needed for long-term conservation planning, which includes quantitative harvest risk assessments that incorporate estimated or assumed trends in future environmental carrying capacity.

Sea ice loss increases genetic isolation in a high Arctic ungulate metapopulation

Sea ice loss may have dramatic consequences for population connectivity, extinction-colonization dynamics, and even the persistence of Arctic species subject to climate change. This is of particular concern in face of additional anthropogenic stressors, such as overexploitation. In this study, we assess the population-genetic implications of diminishing sea ice cover in the endemic, high Arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) by analyzing the interactive effects of landscape barriers and reintroductions (following harvest-induced extirpations) on their metapopulation genetic structure. We genotyped 411 wild reindeer from 25 sampling sites throughout the entire subspecies' range at 19 microsatellite loci. Bayesian clustering analysis showed a genetic structure composed of eight populations, of which two were admixed. Overall population genetic differentiation was high (mean F_ST  = 0.21). Genetic diversity was low (allelic richness [AR] = 2.07-2.58; observed heterozygosity = 0.23-0.43) and declined toward the outer distribution range, where populations showed significant levels of inbreeding. Coalescent estimates of effective population sizes and migration rates revealed strong evolutionary source-sink dynamics with the central population as the main source. The population genetic structure was best explained by a landscape genetics model combining strong isolation by glaciers and open water, and high connectivity by dispersal across winter sea ice. However, the observed patterns of natural isolation were strongly modified by the signature of past harvest-induced extirpations, subsequent reintroductions, and recent lack of sea ice. These results suggest that past and current anthropogenic drivers of metapopulation dynamics may have interactive effects on large-scale ecological and evolutionary processes. Continued loss of sea ice as a dispersal corridor within and between island systems is expected to increase the genetic isolation of populations, and thus threaten the evolutionary potential and persistence of Arctic wildlife.

Canadian polar bear population structure using genome-wide markers

Predicting the consequences of environmental changes, including human-mediated climate change on species, requires that we quantify range-wide patterns of genetic diversity and identify the ecological, environmental, and historical factors that have contributed to it. Here, we generate baseline data on polar bear population structure across most Canadian subpopulations (n = 358) using 13,488 genome-wide single nucleotide polymorphisms (SNPs) identified with double-digest restriction site-associated DNA sequencing (ddRAD). Our ddRAD dataset showed three genetic clusters in the sampled Canadian range, congruent with previous studies based on microsatellites across the same regions; however, due to a lack of sampling in Norwegian Bay, we were unable to confirm the existence of a unique cluster in that subpopulation. These data on the genetic structure of polar bears using SNPs provide a detailed baseline against which future shifts in population structure can be assessed, and opportunities to develop new noninvasive tools for monitoring polar bears across their range.

Are polar bear habitat resource selection functions developed from 1985-1995 data still useful?

Greenhouse-gas-induced warming in the Arctic has caused declines in sea ice extent and changed its composition, raising concerns by all circumpolar nations for polar bear conservation.Negative impacts have been observed in three well-studied polar bear subpopulations. Most subpopulations, however, receive little or no direct monitoring, hence, resource selection functions (RSF) may provide a useful proxy of polar bear distributions. However, the efficacy of RSFs constructed from past data, that is, reference RSFs, may be degraded under contemporary conditions, especially in a rapidly changing environment.We assessed published Arctic-wide reference RSFs using tracking data from adult female polar bears captured in the Beaufort Sea. We compared telemetry-derived seasonal distributions of polar bears to RSF-defined optimal sea ice habitat during the period of RSF model development, 1985-1995, and two subsequent periods with diminished sea ice: 1996-2006 and 2007-2016. From these comparisons, we assessed the applicability of the reference RSFs for contemporary polar bear conservation.In the two decades following the 1985-1995 reference period, use and availability of optimal habitat by polar bears declined during the ice melt, ice minimum, and ice growth seasons. During the ice maximum season (i.e., winter), polar bears used the best habitat available, which changed relatively little across the three decades of study. During the ice melt, ice minimum, and ice growth seasons, optimal habitat in areas used by polar bears decreased and was displaced north and east of the Alaska Beaufort Sea coast. As optimal habitat diminished in these seasons, polar bears expanded their range and occupied greater areas of suboptimal habitat.Synthesis and applications: Sea ice declines due to climate change continue to challenge polar bears and their conservation. The distribution of Southern Beaufort Sea polar bears remained similar during the ice maximum season, so the reference RSFs developed from data collected >20 years ago continue to accurately model their winter distribution. In contrast, reference RSFs for the ice transitional and minimum seasons showed diminished predictive efficacy but were useful in revealing that contemporary polar bears have been increasingly forced to use suboptimal habitats during those seasons.

Geographical assignment of polar bears using multi-element isoscapes

Wide-ranging apex predators are among the most challenging of all fauna to conserve and manage. This is especially true of the polar bear (Ursus maritimus), an iconic predator that is hunted in Canada and threatened by global climate change. We used combinations of stable isotopes (¹³C,¹⁵N,²H,¹⁸O) in polar bear hair from > 1000 individuals, sampled from across much of the Canadian Arctic and sub-Arctic, to test the ability of stable isotopic profiles to ‘assign’ bears to (1) predefined managed subpopulations, (2) subpopulations defined by similarities in stable isotope values using quadratic discriminant analysis, and (3) spatially explicit, isotopically distinct clusters derived from interpolated (i.e. ‘kriged’) isotopic landscapes, or ‘isoscapes’, using the partitioning around medoids algorithm. A four-isotope solution provided the highest overall assignment accuracies (~80%) to pre-existing management subpopulations with accuracy rates ranging from ~30–99% (median = 64%). Assignment accuracies of bears to hierarchically clustered ecological groups based on isotopes ranged from ~64–99%. Multivariate assignment to isotopic clusters resulted in highest assignment accuracies of 68% (33–77%), 84% (47–96%) and 74% (53–85%) using two, three and four stable isotope groups, respectively. The resulting spatial structure inherent in the multiple stable isotopic compositions of polar bear tissues is a powerful forensic tool that will, in this case, contribute to the conservation and management of this species. Currently, it is unclear what is driving these robust isotopic patterns and future research is needed to evaluate the processes behind the pattern. Nonetheless, our isotopic approach can be further applied to other apex mammalian predators under threat, such as the large felids, providing that isotopic structure occurs throughout their range.

Spatial differences in genetic diversity and northward migration suggest genetic erosion along the boreal caribou southern range limit and continued range retraction

With increasing human activities and associated landscape changes, distributions of terrestrial mammals become fragmented. These changes in distribution are often associated with reduced population sizes and loss of genetic connectivity and diversity (i.e., genetic erosion) which may further diminish a species' ability to respond to changing environmental conditions and lead to local population extinctions. We studied threatened boreal caribou (Rangifer tarandus caribou) populations across their distribution in Ontario/Manitoba (Canada) to assess changes in genetic diversity and connectivity in areas of high and low anthropogenic activity. Using data from >1,000 caribou and nine microsatellite loci, we assessed population genetic structure, genetic diversity, and recent migration rates using a combination of network and population genetic analyses. We used Bayesian clustering analyses to identify population genetic structure and explored spatial and temporal variation in those patterns by assembling networks based on R ST and F ST as historical and contemporary genetic edge distances, respectively. The Bayesian clustering analyses identified broad-scale patterns of genetic structure and closely aligned with the R ST network. The F ST network revealed substantial contemporary genetic differentiation, particularly in areas presenting contemporary anthropogenic disturbances and habitat fragmentation. In general, relatively lower genetic diversity and greater genetic differentiation were detected along the southern range limit, differing from areas in the northern parts of the distribution. Moreover, estimation of migration rates suggested a northward movement of animals away from the southern range limit. The patterns of genetic erosion revealed in our study suggest ongoing range retraction of boreal caribou in central Canada.

Tracking arctic marine mammal resilience in an era of rapid ecosystem alteration

Global warming is significantly altering arctic marine ecosystems. Specifically, the precipitous loss of sea ice is creating a dichotomy between ice-dependent polar bears and pinnipeds that are losing habitat and some cetaceans that are gaining habitat. While final outcomes are hard to predict for the many and varied marine mammal populations that rely on arctic habitats, we suggest a simplified framework to assess status, based upon ranking a population's size, range, behavior, and health. This basic approach is proposed as a means to prioritize and expedite conservation and management efforts in an era of rapid ecosystem alteration.