Glacial ice supports a distinct and undocumented polar bear subpopulation persisting in late 21st-century sea-ice conditions.

Polar bears are susceptible to climate warming because of their dependence on sea ice, which is declining rapidly. We present the first evidence for a genetically distinct and functionally isolated group of polar bears in Southeast Greenland. These bears occupy sea-ice conditions resembling those projected for the High Arctic in the late 21st century, with an annual ice-free period that is >100 days longer than the estimated fasting threshold for the species. Whereas polar bears in most of the Arctic depend on annual sea ice to catch seals, Southeast Greenland bears have a year-round hunting platform in the form of freshwater glacial mélange. This suggests that marine-terminating glaciers, although of limited availability, may serve as previously unrecognized climate refugia. Conservation of Southeast Greenland polar bears, which meet criteria for recognition as the world's 20th polar bear subpopulation, is necessary to preserve the genetic diversity and evolutionary potential of the species.

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.

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.

Vulnerability of Arctic marine mammals to vessel traffic in the increasingly ice-free Northwest Passage and Northern Sea Route.

The fabled Northwest Passage and Northern Sea Route that were once the quests of early Western explorers are now increasingly sea ice-free, with routine vessel transits expected by midcentury. The potential impacts of this novel vessel traffic on endemic Arctic marine mammal (AMM) species are unknown despite their critical social and ecological roles in the ecosystem and widely recognized susceptibility to ice loss. We developed a vulnerability assessment of 80 subpopulations of seven AMM species to vessel traffic during the ice-free season. Vulnerability scores were based on the combined influence of spatially explicit exposure to the sea routes and a suite of sensitivity variables. More than half of AMM subpopulations (42/80) are exposed to open-water vessel transits in the Arctic sea routes. Narwhals (Monodon monoceros) were estimated to be most vulnerable to vessel impacts, given their high exposure and sensitivity, and polar bears (Ursus maritimus) were estimated to be the least vulnerable because of their low exposure and sensitivity. Regions with geographic bottlenecks, such as the Bering Strait and eastern Canadian Arctic, were characterized by two to three times higher vulnerability than more remote regions. These pinch points are obligatory pathways for both vessels and migratory AMMs, and so represent potentially high conflict areas but also opportunities for conservation-informed planning. Some of the species and regions identified as least vulnerable were also characterized by high uncertainty, highlighting additional data and monitoring needs. Our quantification of the heterogeneity of risk across AMM species provides a necessary first step toward developing best practices for maritime industries poised to advance into this rapidly changing seascape.

Harvesting wildlife affected by climate change: a modelling and management approach for polar bears.

The conservation of many wildlife species requires understanding the demographic effects of climate change, including interactions between climate change and harvest, which can provide cultural, nutritional or economic value to humans.We present a demographic model that is based on the polar bear Ursus maritimus life cycle and includes density-dependent relationships linking vital rates to environmental carrying capacity (K). Using this model, we develop a state-dependent management framework to calculate a harvest level that (i) maintains a population above its maximum net productivity level (MNPL; the population size that produces the greatest net increment in abundance) relative to a changing K, and (ii) has a limited negative effect on population persistence.Our density-dependent relationships suggest that MNPL for polar bears occurs at approximately 0·69 (95% CI = 0·63-0·74) of K. Population growth rate at MNPL was approximately 0·82 (95% CI = 0·79-0·84) of the maximum intrinsic growth rate, suggesting relatively strong compensation for human-caused mortality.Our findings indicate that it is possible to minimize the demographic risks of harvest under climate change, including the risk that harvest will accelerate population declines driven by loss of the polar bear's sea-ice habitat. This requires that (i) the harvest rate - which could be 0 in some situations - accounts for a population's intrinsic growth rate, (ii) the harvest rate accounts for the quality of population data (e.g. lower harvest when uncertainty is large), and (iii) the harvest level is obtained by multiplying the harvest rate by an updated estimate of population size. Environmental variability, the sex and age of removed animals and risk tolerance can also affect the harvest rate. Synthesis and applications. We present a coupled modelling and management approach for wildlife that accounts for climate change and can be used to balance trade-offs among multiple conservation goals. In our example application to polar bears experiencing sea-ice loss, the goals are to maintain population viability while providing continued opportunities for subsistence harvest. Our approach may be relevant to other species for which near-term management is focused on human factors that directly influence population dynamics within the broader context of climate-induced habitat degradation.

Habitat selection by two beluga whale populations in the Chukchi and Beaufort seas.

There has been extensive sea ice loss in the Chukchi and Beaufort seas where two beluga whale (Delphinapterus leucas) populations occur between July-November. Our goal was to develop population-specific beluga habitat selection models that quantify relative use of sea ice and bathymetric features related to oceanographic processes, which can provide context to the importance of changing sea ice conditions. We established habitat selection models that incorporated daily sea ice measures (sea ice concentration, proximity to ice edge and dense ice) and bathymetric features (slope, depth, proximity to the continental slope, Barrow Canyon, and shore) to establish quantitative estimates of habitat use for the Eastern Chukchi Sea ('Chukchi') and Eastern Beaufort Sea ('Beaufort') populations. We applied 'used v. available' resource selection functions to locations of 65 whales tagged from 1993-2012, revealing large variations in seasonal habitat selection that were distinct between sex and population groups. Chukchi whales of both sexes were predicted to use areas in close proximity to Barrow Canyon (typically <200 km) as well as the continental slope in summer, although deeper water and denser ice were stronger predictors for males than females. Habitat selection differed more between sexes for Beaufort belugas. Beaufort males selected higher ice concentrations (≥40%) than females (0-40%) in July-August. Proximity to shore (<200 km) strongly predicted summer habitat of Beaufort females, while distance to the ice edge was important for male habitat selection, especially during westward migration in September. Overall, our results indicate that sea ice variables were rarely the primary drivers of beluga summer-fall habitat selection. While diminished sea ice may indirectly affect belugas through changes in the ecosystem, associations with bathymetric features that affect prey availability seemed key to habitat selection during summer and fall. These results provide a benchmark by which to assess future changes in beluga habitat use of the Pacific Arctic.

Decadal shifts in autumn migration timing by Pacific Arctic beluga whales are related to delayed annual sea ice formation.

Migrations are often influenced by seasonal environmental gradients that are increasingly being altered by climate change. The consequences of rapid changes in Arctic sea ice have the potential to affect migrations of a number of marine species whose timing is temporally matched to seasonal sea ice cover. This topic has not been investigated for Pacific Arctic beluga whales (Delphinapterus leucas) that follow matrilineally maintained autumn migrations in the waters around Alaska and Russia. For the sympatric Eastern Chukchi Sea ('Chukchi') and Eastern Beaufort Sea ('Beaufort') beluga populations, we examined changes in autumn migration timing as related to delayed regional sea ice freeze-up since the 1990s, using two independent data sources (satellite telemetry data and passive acoustics) for both populations. We compared dates of migration between 'early' (1993-2002) and 'late' (2004-2012) tagging periods. During the late tagging period, Chukchi belugas had significantly delayed migrations (by 2 to >4 weeks, depending on location) from the Beaufort and Chukchi seas. Spatial analyses also revealed that departure from Beaufort Sea foraging regions by Chukchi whales was postponed in the late period. Chukchi beluga autumn migration timing occurred significantly later as regional sea ice freeze-up timing became later in the Beaufort, Chukchi, and Bering seas. In contrast, Beaufort belugas did not shift migration timing between periods, nor was migration timing related to freeze-up timing, other than for southward migration at the Bering Strait. Passive acoustic data from 2008 to 2014 provided independent and supplementary support for delayed migration from the Beaufort Sea (4 day yr-1 ) by Chukchi belugas. Here, we report the first phenological study examining beluga whale migrations within the context of their rapidly transforming Pacific Arctic ecosystem, suggesting flexible responses that may enable their persistence yet also complicate predictions of how belugas may fare in the future.