Demography of an apex predator at the edge of its range: impacts of changing sea ice on polar bears in Hudson Bay

An interspecific foraging association with polar bears increases foraging opportunities for avian predators in a declining Arctic seabird colony

Interspecific foraging associations (IFAs) are biological interactions where two or more species forage in association with each other. Climate-induced reductions in Arctic sea ice have increased polar bear (Ursus maritimus) foraging in seabird colonies, which creates foraging opportunities for avian predators. We used drone video of bears foraging within a common eider (Somateria mollissima) colony on East Bay Island (Nunavut, Canada) in 2017 to investigate herring gull (Larus argentatus) foraging in association with bears. We recorded nest visitation by gulls following n = 193 eider flushing events from nests during incubation. The probability of gulls visiting eider nests increased with higher number of gulls present (β = 0.14 ± 0.03 [SE], p < .001) and for nests previously visited by a bear (β = 1.14 ± 0.49 [SE], p < .02). In our model examining the probability of gulls consuming eggs from nests, we failed to detect statistically significant effects for the number of gulls present (β = 0.09 ± 0.05 [SE], p < .07) or for nests previously visited by a bear (β = -0.92 ± 0.71 [SE], p < .19). Gulls preferred to visit nests behind bears (χ2 = 18, df = 1, p < .0001), indicating gulls are risk averse in the presence of polar bears. Our study provides novel insights on an Arctic IFA, and we present evidence that gulls capitalize on nests made available due to disturbance associated with foraging bears, as eiders disturbed off their nest allow gulls easier access to eggs. We suggest the IFA between gulls and polar bears is parasitic, as gulls are consuming terrestrial resources which would have eventually been consumed by bears. This finding has implications for estimating the energetic contribution of bird eggs to polar bear summer diets in that the total number of available clutches to consume may be reduced due to avian predators.

Polar bear energetic and behavioral strategies on land with implications for surviving the ice-free period

Declining Arctic sea ice is increasing polar bear land use. Polar bears on land are thought to minimize activity to conserve energy. Here, we measure the daily energy expenditure (DEE), diet, behavior, movement, and body composition changes of 20 different polar bears on land over 19-23 days from August to September (2019-2022) in Manitoba, Canada. Polar bears on land exhibited a 5.2-fold range in DEE and 19-fold range in activity, from hibernation-like DEEs to levels approaching active bears on the sea ice, including three individuals that made energetically demanding swims totaling 54-175 km. Bears consumed berries, vegetation, birds, bones, antlers, seal, and beluga. Beyond compensating for elevated DEE, there was little benefit from terrestrial foraging toward prolonging the predicted time to starvation, as 19 of 20 bears lost mass (0.4-1.7 kg•day-1). Although polar bears on land exhibit remarkable behavioral plasticity, our findings reinforce the risk of starvation, particularly in subadults, with forecasted increases in the onshore period.

Survival improvements of marine mammals in zoological institutions mirror historical advances in human longevity

An intense public debate has fuelled governmental bans on marine mammals held in zoological institutions. The debate rests on the assumption that survival in zoological institutions has been and remains lower than in the wild, albeit the scientific evidence in support of this notion is equivocal. Here, we used statistical methods previously applied to assess historical improvements in human lifespan and data on 8864 individuals of four marine mammal species (harbour seal, Phoca vitulina; California sea lion, Zalophus californianus; polar bear, Ursus maritimus; common bottlenose dolphin, Tursiops truncatus) held in zoos from 1829 to 2020. We found that life expectancy increased up to 3.40 times, and first-year mortality declined up to 31%, during the last century in zoos. Moreover, the life expectancy of animals in zoos is currently 1.65-3.55 times longer than their wild counterparts. Like humans, these improvements have occurred concurrently with advances in management practices, crucial for population welfare. Science-based decisions will help effective legislative changes and ensure better implementation of animal care.

Arctic and red fox population responses to climate and cryosphere changes at the Arctic's edge

Responses of one species to climate change may influence the population dynamics of others, particularly in the Arctic where food webs are strongly linked. Specifically, changes to the cryosphere may limit prey availability for predators. We examined Arctic (Vulpes lagopus) and red fox (V. vulpes) population dynamics near the southern edge of the Arctic fox distribution using fur harvest records from Churchill, Manitoba, Canada between 1955 and 2012. Arctic foxes showed a declining population trend over time (inferred from harvest records corrected for trapping effort), whereas the red fox population trend was relatively stable. The positive relationship between the annual Arctic and red fox harvests suggested interspecific competition did not promote the Arctic fox decline. To investigate alternative mechanisms, we evaluated the relative influence of sea-ice phenology, snow depth, snow duration, winter thaws, and summer temperature on the harvest dynamics of both species in the most recent 32 years (1980-2012; n = 29) of our data. Arctic fox harvests were negatively related to the length of time Hudson Bay was free of sea ice. Shorter sea ice duration may reduce access to seal carrion as an alternative winter food source when lemming densities decline. Contrary to our prediction, red fox harvest was not related to summer temperature but was positively related to snow depth, suggesting winter prey availability may limit red fox population growth. Predators have an important ecological role, so understanding the influence of changes in the cryosphere on predator-prey interactions may better illuminate the broader influence of climate change on food-web dynamics.

Widespread Exposure to Mosquitoborne California Serogroup Viruses in Caribou, Arctic Fox, Red Fox, and Polar Bears, Canada

Northern Canada is warming at 3 times the global rate. Thus, changing diversity and distribution of vectors and pathogens is an increasing health concern. California serogroup (CSG) viruses are mosquitoborne arboviruses; wildlife reservoirs in northern ecosystems have not been identified. We detected CSG virus antibodies in 63% (95% CI 58%-67%) of caribou (n = 517), 4% (95% CI 2%-7%) of Arctic foxes (n = 297), 12% (95% CI 6%-21%) of red foxes (n = 77), and 28% (95% CI 24%-33%) of polar bears (n = 377). Sex, age, and summer temperatures were positively associated with polar bear exposure; location, year, and ecotype were associated with caribou exposure. Exposure was highest in boreal caribou and increased from baseline in polar bears after warmer summers. CSG virus exposure of wildlife is linked to climate change in northern Canada and sustained surveillance could be used to measure human health risks.

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.

Reproductive senescence in polar bears in a variable environment

Reproductive senescence is ubiquitous in mammals. However, patterns of senescence vary across reproductive traits, even within populations, perhaps because of differences in selection pressures, physiological constraints, and responses to environmental conditions. We investigated reproductive senescence in wild female polar bears (Ursus maritimus), using 31 years of capture-recapture data from the Svalbard area. We studied the influence of environmental conditions on age-specific litter production and litter size using generalized linear mixed models. Further, using a capture-recapture model that handles the dependency between vital rates of individuals belonging to the same family unit, we assessed maternal-age-related changes in first year cub and litter survival. We provide clear evidence for reproductive senescence in female polar bears. Litter production and litter size peaked in middle-aged females and declined sharply afterward. By contrast cub and litter survival did not decline after prime age. We found no evidence of terminal investment. The reproductive output of all females was affected by sea-ice conditions during the previous year and the Arctic Oscillation, with some effects differing greatly between age groups. Old females were affected the most by environmental conditions. Our results suggest that the decline in reproductive output is a combination of fertility and body-condition senescence, with a weak contribution of maternal-effect senescence, possibly due to benefits of experience. Further, as predicted by evolutionary theory, senescence appears to be a consequence of failures in early stages of the reproductive cycle rather than in late stages, and environmental variation affected old females more than prime-aged females. Our study emphasizes the need to study several reproductive traits and account for environmental variation when investigating reproductive senescence. Differences in senescence patterns across reproductive traits should be interpreted in light of evolutionary theory and while considering underlying physiological drivers.

Modeling spatiotemporal abundance and movement dynamics using an integrated spatial capture-recapture movement model

Animal movement is a fundamental ecological process affecting the survival and reproduction of individuals, the structure of populations, and the dynamics of communities. Methods to quantify animal movement and spatiotemporal abundances, however, are generally separate and therefore omit linkages between individual-level and population-level processes. We describe an integrated spatial capture-recapture (SCR) movement model to jointly estimate (1) the number and distribution of individuals in a defined spatial region and (2) movement of those individuals through time. We applied our model to a study of polar bears (Ursus maritimus) in a 28,125 km² survey area of the eastern Chukchi Sea, USA in 2015 that incorporated capture-recapture and telemetry data. In simulation studies, the model provided unbiased estimates of movement, abundance, and detection parameters using a bivariate normal random walk and correlated random walk movement process. Our case study provided detailed evidence of directional movement persistence for both male and female bears, where individuals regularly traversed areas larger than the survey area during the 36-day study period. Scaling from individual- to population-level inferences, we found that densities varied from <0.75 bears/625 km² grid cell/day in nearshore cells to 1.6-2.5 bears/grid cell/day for cells surrounded by sea ice. Daily abundance estimates ranged from 53 to 69 bears, with no trend across days. The cumulative number of unique bears that used the survey area increased through time due to movements into and out of the area, resulting in an estimated 171 individuals using the survey area during the study (95% credible interval 124-250). Abundance estimates were similar to a previous multiyear integrated population model using capture-recapture and telemetry data (2008-2016; Regehr et al., Scientific Reports 8:16780, 2018). Overall, the SCR-movement model successfully quantified both individual- and population-level space use, including the effects of landscape characteristics on movement, abundance, and detection, while linking the movement and abundance processes to directly estimate density within a prescribed spatial region and temporal period. Integrated SCR-movement models provide a generalizable approach to incorporate greater movement realism into population dynamics and link movement to emergent properties including spatiotemporal densities and abundances.

Anthropogenic food: an emerging threat to polar bears

Supplemental food from anthropogenic sources is a source of conflict with humans for many wildlife species. Food-seeking behaviours by black bears Ursus americanus and brown bears Ursus arctos can lead to property damage, human injury and mortality of the offending bears. Such conflicts are a well-known conservation management issue wherever people live in bear habitats. In contrast, the use of anthropogenic foods by the polar bear Ursus maritimus is less common historically but is a growing conservation and management issue across the Arctic. Here we present six case studies that illustrate how negative food-related interactions between humans and polar bears can become either chronic or ephemeral and unpredictable. Our examination suggests that attractants are an increasing problem, exacerbated by climate change-driven sea-ice losses that cause increased use of terrestrial habitats by bears. Growing human populations and increased human visitation increase the likelihood of human–polar bear conflict. Efforts to reduce food conditioning in polar bears include attractant management, proactive planning and adequate resources for northern communities to reduce conflicts and improve human safety. Permanent removal of unsecured sources of nutrition, to reduce food conditioning, should begin immediately at the local level as this will help to reduce polar bear mortality.

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.

Polar Bear Harvest Patterns Across the Circumpolar Arctic

Wildlife harvest remains a conservation concern for many species and assessing patterns of harvest can provide insights on sustainability and inform management. Polar bears (Ursus maritimus) are harvested over a large part of their range by local people. The species has a history of unsustainable harvest that was largely rectified by an international agreement that required science-based management. The objective of our study was to examine the temporal patterns in the number of polar bears harvested, harvest sex ratios, and harvest rates from 1970 to 2018. We analyzed data from 39,049 harvested polar bears (annual mean 797 bears) collected from 1970 to 2018. Harvest varied across populations and times that reflect varying management objectives, episodic events, and changes based on new population estimates. More males than females were harvested with an overall M:F sex ratio of 1.84. Harvest varied by jurisdiction with 68.0% of bears harvested in Canada, 18.0% in Greenland, 11.8% in the USA, and 2.2% in Norway. Harvest rate was often near the 4.5% target rate. Where data allowed harvest rate estimation, the target rate was exceeded in 11 of 13 populations with 1–5 populations per year above the target since 1978. Harvest rates at times were up to 15.9% of the estimated population size suggesting rare episodes of severe over-harvest. Harvest rate was unrelated to a proxy for ecosystem productivity (area of continental shelf within each population) but was correlated with prey diversity. In the last 5–10 years, monitored populations all had harvest rates near sustainable limits, suggesting improvements in management. Polar bear harvest management has reduced the threat it once posed to the species. However, infrequent estimates of abundance, new management objectives, and climate change have raised new concerns about the effects of harvest.

The Adrenal Cortisol Response to Increasing Ambient Temperature in Polar Bears (Ursus maritimus)

Our objective was to identify the upper ambient temperature threshold that triggers an increase in cortisol in response to increased thermoregulatory demands in polar bears. The results reported here include endocrine data collected over two years from 25 polar bears housed in 11 accredited zoological institutions across North America. The effects of ambient temperature, sex, age group (juvenile, adult, elderly), breeding season and humidity on fecal cortisol metabolite (FCM) concentrations (N = 8439 samples) were evaluated using linear mixed models. Ambient temperatures were placed into five different categories: <5 °C, 6−10 °C, 11−15 °C, 16−20 °C, and >20 °C. Ambient temperature and humidity had a significant (p < 0.05) effect on FCM concentrations with significant (p < 0.05) interactions of sex, age and breeding season. Once biotic factors were accounted for, there was a significant (p < 0.05) increase in FCM concentrations associated with ambient temperatures above 20 °C in adult polar bears. The implications of these findings for the management of both zoo and wild polar bears are discussed.

Demographic risk assessment for a harvested species threatened by climate change: polar bears in the Chukchi Sea

Climate change threatens global biodiversity. Many species vulnerable to climate change are important to humans for nutritional, cultural, and economic reasons. Polar bears Ursus maritimus are threatened by sea-ice loss and represent a subsistence resource for Indigenous people. We applied a novel population modeling-management framework that is based on species life history and accounts for habitat loss to evaluate subsistence harvest for the Chukchi Sea (CS) polar bear subpopulation. Harvest strategies followed a state-dependent approach under which new data were used to update the harvest on a predetermined management interval. We found that a harvest strategy with a starting total harvest rate of 2.7% (˜85 bears/yr at current abundance), a 2:1 male-to-female ratio, and a 10-yr management interval would likely maintain subpopulation abundance above maximum net productivity level for the next 35 yr (approximately three polar bear generations), our primary criterion for sustainability. Plausible bounds on starting total harvest rate were 1.7-3.9%, where the range reflects uncertainty due to sampling variation, environmental variation, model selection, and differing levels of risk tolerance. The risk of undesired demographic outcomes (e.g., overharvest) was positively related to harvest rate, management interval, and projected declines in environmental carrying capacity; and negatively related to precision in population data. Results reflect several lines of evidence that the CS subpopulation has been productive in recent years, although it is uncertain how long this will last as sea-ice loss continues. Our methods provide a template for balancing trade-offs among protection, use, research investment, and other factors. Demographic risk assessment and state-dependent management will become increasingly important for harvested species, like polar bears, that exhibit spatiotemporal variation in their response to climate change.

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.

Climate change impacts on seabirds and marine mammals: The importance of study duration, thermal tolerance and generation time

Understanding climate change impacts on top predators is fundamental to marine biodiversity conservation, due to their increasingly threatened populations and their importance in marine ecosystems. We conducted a systematic review of the effects of climate change (prolonged, directional change) and climate variability on seabirds and marine mammals. We extracted data from 484 studies (4808 published studies were reviewed), comprising 2215 observations on demography, phenology, distribution, diet, behaviour, body condition and physiology. The likelihood of concluding that climate change had an impact increased with study duration. However, the temporal thresholds for the effects of climate change to be discernibly varied from 10 to 29 years depending on the species, the biological response and the oceanic study region. Species with narrow thermal ranges and relatively long generation times were more often reported to be affected by climate change. This provides an important framework for future assessments, with guidance on response- and region-specific temporal dimensions that need to be considered when reporting effects of climate change. Finally, we found that tropical regions and non-breeding life stages were poorly covered in the literature, a concern that should be addressed to enable a better understanding of the vulnerability of marine predators to climate change.

Fatty acid profiles of feeding and fasting bears: estimating calibration coefficients, the timeframe of diet estimates, and selective mobilization during hibernation

Accurate information on diet composition is central to understanding and conserving carnivore populations. Quantitative fatty acid signature analysis (QFASA) has emerged as a powerful tool for estimating the diets of predators, but ambiguities remain about the timeframe of QFASA estimates and the need to account for species-specific patterns of metabolism. We conducted a series of feeding experiments with four juvenile male brown bears (Ursus arctos) to (1) track the timing of changes in adipose tissue composition and QFASA diet estimates in response to a change in diet and (2) quantify the relationship between consumer and diet FA composition (i.e., determine "calibration coefficients"). Bears were fed three compositionally distinct diets for 90-120 days each. Two marine-based diets were intended to approximate the lipid content and composition of the wild diet of polar bears (U. maritimus). Bear adipose tissue composition changed quickly in the direction of the diet and showed evidence of stabilization after 60 days. During hibernation, FA profiles were initially stable but diet estimates after 10 weeks were sensitive to calibration coefficients. Calibration coefficients derived from the marine-based diets were broadly similar to each other and to published values from marine-fed mink (Mustela vison), which have been used as a model for free-ranging polar bears. For growing bears on a high-fat diet, the temporal window for QFASA estimates was 30-90 days. Although our results reinforce the importance of accurate calibration, the similarities across taxa and diets suggest it may be feasible to develop a generalized QFASA approach for mammalian carnivores.

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.

Long-term increases in pathogen seroprevalence in polar bears (Ursus maritimus) influenced by climate change

The influence of climate change on wildlife disease dynamics is a burgeoning conservation and human health issue, but few long-term studies empirically link climate to pathogen prevalence. Polar bears (Ursus maritimus) are vulnerable to the negative impacts of sea ice loss as a result of accelerated Arctic warming. While studies have associated changes in polar bear body condition, reproductive output, survival, and abundance to reductions in sea ice, no long-term studies have documented the impact of climate change on pathogen exposure. We examined 425 serum samples from 381 adult polar bears, collected in western Hudson Bay (WH), Canada, for antibodies to selected pathogens across three time periods: 1986-1989 (n = 157), 1995-1998 (n = 159) and 2015-2017 (n = 109). We ran serological assays for antibodies to seven pathogens: Toxoplasma gondii, Neospora caninum, Trichinella spp., Francisella tularensis, Bordetella bronchiseptica, canine morbillivirus (CDV) and canine parvovirus (CPV). Seroprevalence of zoonotic parasites (T. gondii, Trichinella spp.) and bacterial pathogens (F. tularensis, B. bronchiseptica) increased significantly between 1986-1989 and 1995-1998, ranging from +6.2% to +20.8%, with T. gondii continuing to increase into 2015-2017 (+25.8% overall). Seroprevalence of viral pathogens (CDV, CPV) and N. caninum did not change with time. Toxoplasma gondii seroprevalence was higher following wetter summers, while seroprevalences of Trichinella spp. and B. bronchiseptica were positively correlated with hotter summers. Seroprevalence of antibodies to F. tularensis increased following years polar bears spent more days on land, and polar bears previously captured in human settlements were more likely to be seropositive for Trichinella spp. As the Arctic has warmed due to climate change, zoonotic pathogen exposure in WH polar bears has increased, driven by numerous altered ecosystem pathways.

Survival and abundance of polar bears in Alaska's Beaufort Sea, 2001-2016

The Arctic Ocean is undergoing rapid transformation toward a seasonally ice-free ecosystem. As ice-adapted apex predators, polar bears (Ursus maritimus) are challenged to cope with ongoing habitat degradation and changes in their prey base driven by food-web response to climate warming. Knowledge of polar bear response to environmental change is necessary to understand ecosystem dynamics and inform conservation decisions. In the southern Beaufort Sea (SBS) of Alaska and western Canada, sea ice extent has declined since satellite observations began in 1979 and available evidence suggests that the carrying capacity of the SBS for polar bears has trended lower for nearly two decades. In this study, we investigated the population dynamics of polar bears in Alaska's SBS from 2001 to 2016 using a multistate Cormack-Jolly-Seber mark-recapture model. States were defined as geographic regions, and we used location data from mark-recapture observations and satellite-telemetered bears to model transitions between states and thereby explain heterogeneity in recapture probabilities. Our results corroborate prior findings that the SBS subpopulation experienced low survival from 2003 to 2006. Survival improved modestly from 2006 to 2008 and afterward rebounded to comparatively high levels for the remainder of the study, except in 2012. Abundance moved in concert with survival throughout the study period, declining substantially from 2003 and 2006 and afterward fluctuating with lower variation around an average of 565 bears (95% Bayesian credible interval [340, 920]) through 2015. Even though abundance was comparatively stable and without sustained trend from 2006 to 2015, polar bears in the Alaska SBS were less abundant over that period than at any time since passage of the U.S. Marine Mammal Protection Act. The potential for recovery is likely limited by the degree of habitat degradation the subpopulation has experienced, and future reductions in carrying capacity are expected given current projections for continued climate warming.

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.

Are tissue samples obtained via remote biopsy useful for fatty acid-based diet analyses in a free-ranging carnivore?

Fundamental knowledge on free-ranging animals has been obtained through capture-based studies; however, these may be logistically intensive, financially expensive, and potentially inconsistent with local cultural values. Genetic mark–recapture using remote tissue sampling has emerged as a less invasive alternative to capture-based population surveys but provides fewer opportunities to collect samples and measurements for broader ecological studies. We compared lipid content, fatty acid (FA) composition, and diet estimates from adipose tissue of polar bears (Ursus maritimus) obtained from two collection methods: remote biopsies (n = 138) sampled from helicopters and hunter-collected tissue (n = 499) from bears harvested in Davis Strait and Gulf of Boothia, Nunavut, 2010 – 2018. Lipid content of adipose tissue was lower in remote biopsies than harvest samples likely because remote biopsies removed only the outermost layer of subcutaneous tissue, rather than the more metabolically dynamic innermost tissue obtained from harvest samples. In contrast, FA composition was similar between the two collection methods with relatively small proportional differences in individual FAs. For diet estimates in Davis Strait, collection method was not a predictor of prey contribution to diet. In Gulf of Boothia, collection method was a predictor for some prey types, but the differences were relatively minor; the rank order of prey types was similar (e.g., ringed seal; Pusa hispida was consistently the primary prey in diets) and prey proportions differed by &lt; 6% between the collection methods. Results from both methods showed that diets varied by geographic area, season, year, age class, and sex. Our study demonstrates that adipose tissue from remote biopsy provides reliable estimates of polar bear diet based on FA analysis and can be used to monitor underlying ecological changes in Arctic marine food webs.

Bridging gaps in demographic analysis with phylogenetic imputation

Phylogenetically informed imputation methods have rarely been applied to estimate missing values in demographic data but may be a powerful tool for reconstructing vital rates of survival, maturation, and fecundity for species of conservation concern. Imputed vital rates could be used to parameterize demographic models to explore how populations respond when vital rates are perturbed. We used standardized vital rate estimates for 50 bird species to assess the use of phylogenetic imputation to fill gaps in demographic data. We calculated imputation accuracy for vital rates of focal species excluded from the data set either singly or in combination and with and without phylogeny, body mass, and life-history trait data. We used imputed vital rates to calculate demographic metrics, including generation time, to validate the use of imputation in demographic analyses. Covariance among vital rates and other trait data provided a strong basis to guide imputation of missing vital rates in birds, even in the absence of phylogenetic information. Mean NRMSE for null and phylogenetic models differed by <0.01 except when no vital rates were available or for vital rates with high phylogenetic signal (Pagel's λ > 0.8). In these cases, including body mass and life-history trait data compensated for lack of phylogenetic information: mean normalized root mean square error (NRMSE) for null and phylogenetic models differed by <0.01 for adult survival and <0.04 for maturation rate. Estimates of demographic metrics were sensitive to the accuracy of imputed vital rates. For example, mean error in generation time doubled in response to inaccurate estimates of maturation time. Accurate demographic data and metrics, such as generation time, are needed to inform conservation planning processes, for example through International Union for Conservation of Nature Red List assessments and population viability analysis. Imputed vital rates could be useful in this context but, as for any estimated model parameters, awareness of the sensitivities of demographic model outputs to the imputed vital rates is essential.

Age-structured Jolly-Seber model expands inference and improves parameter estimation from capture-recapture data

Understanding the influence of individual attributes on demographic processes is a key objective of wildlife population studies. Capture-recapture and age data are commonly collected to investigate hypotheses about survival, reproduction, and viability. We present a novel age-structured Jolly-Seber model that incorporates age and capture-recapture data to provide comprehensive information on population dynamics, including abundance, age-dependent survival, recruitment, age structure, and population growth rates. We applied our model to a multi-year capture-recapture study of polar bears (Ursus maritimus) in western Hudson Bay, Canada (2012–2018), where management and conservation require a detailed understanding of how polar bears respond to climate change and other factors. In simulation studies, the age-structured Jolly-Seber model improved precision of survival, recruitment, and annual abundance estimates relative to standard Jolly-Seber models that omit age information. Furthermore, incorporating age information improved precision of population growth rates, increased power to detect trends in abundance, and allowed direct estimation of age-dependent survival and changes in annual age structure. Our case study provided detailed evidence for senescence in polar bear survival. Median survival estimates were lower (<0.95) for individuals aged <5 years, remained high (>0.95) for individuals aged 7–22 years, and subsequently declined to near zero for individuals >30 years. We also detected cascading effects of large recruitment classes on population age structure, which created major shifts in age structure when these classes entered the population and then again when they reached prime breeding ages (10–15 years old). Overall, age-structured Jolly-Seber models provide a flexible means to investigate ecological and evolutionary processes that shape populations (e.g., via senescence, life expectancy, and lifetime reproductive success) while improving our ability to investigate population dynamics and forecast population changes from capture-recapture data.

Migration dynamics of polar bears (Ursus maritimus) in western Hudson Bay

Migration is predicted to change both spatially and temporally as climate change alters seasonal resource availability. Species in extreme environments are especially susceptible to climate change; hence, it is important to determine environmental and biological variables that influence their migration. Polar bears (Ursus maritimus) are an Arctic apex carnivore whose migration phenology has been affected by climate change and is vulnerable to future changes. Here, we used satellite-linked telemetry collar data from adult female polar bears in western Hudson Bay from 2004 to 2016 and multivariate response regression models to demonstrate that 1) spatial and temporal migration metrics are correlated, 2) ice concentration and wind are important environmental variables that influence polar bear migration in seasonal ice areas, and 3) migration did not vary across the years of our study, highlighting the importance of continued monitoring. Specifically, we found that ice concentration, wind speed, and wind direction affected polar bear migration onto ice during freeze-up and ice concentration and wind direction affected migration onto land during breakup. Bears departed from land earlier with increased wind speed and the effect of wind direction on migration may be linked to prey searching and ice drift. Low ice concentration was associated with higher movement during freeze-up and breakup. Our findings suggest that migration movement may increase in response to climate change as ice concentration and access to prey declines, potentially increasing nutritional stress on bears.

Seal body condition and atmospheric circulation patterns influence polar bear body condition, recruitment, and feeding ecology in the Chukchi Sea

Polar bears (Ursus maritimus) are experiencing loss of sea ice habitats used to access their marine mammal prey. Simultaneously, ocean warming is changing ecosystems that support marine mammal populations. The interactive effects of sea ice and prey are not well understood yet may explain spatial-temporal variation in the response of polar bears to sea ice loss. Here, we examined the potential combined effects of sea ice, seal body condition, and atmospheric circulation patterns on the body condition, recruitment, diet, and feeding probability of 469 polar bears captured in the Chukchi Sea, 2008-2017. The body condition of ringed seals (Pusa hispida), the primary prey of females and subadults, was related to dietary proportions of ringed seal, feeding probability, and the body condition of females and cubs. In contrast, adult males consumed more bearded seals (Erignathus barbatus) and exhibited better condition when bearded seal body condition was higher. The litter size, number of yearlings per adult female, and the condition of dependent young were higher following winters characterized by low Arctic Oscillation conditions, consistent with a growing number of studies. Body condition, recruitment, and feeding probability were either not associated or negatively associated with sea ice conditions, suggesting that, unlike some subpopulations, Chukchi Sea bears are not currently limited by sea ice availability. However, spring sea ice cover declined 2% per year during our study reaching levels not previously observed in the satellite record and resulting in the loss of polar bear hunting and seal pupping habitat. Our study suggests that the status of ice seal populations is likely an important factor that can either compound or mitigate the response of polar bears to sea ice loss over the short term. In the long term, neither polar bears nor their prey are likely robust to limitless loss of their sea ice habitat.

Aerial survey estimates of polar bears and their tracks in the Chukchi Sea

Polar bears are of international conservation concern due to climate change but are difficult to study because of low densities and an expansive, circumpolar distribution. In a collaborative U.S.-Russian effort in spring of 2016, we used aerial surveys to detect and estimate the abundance of polar bears on sea ice in the Chukchi Sea. Our surveys used a combination of thermal imagery, digital photography, and human observations. Using spatio-temporal statistical models that related bear and track densities to physiographic and biological covariates (e.g., sea ice extent, resource selection functions derived from satellite tags), we predicted abundance and spatial distribution throughout our study area. Estimates of 2016 abundance ([Formula: see text]) ranged from 3,435 (95% CI: 2,300-5,131) to 5,444 (95% CI: 3,636-8,152) depending on the proportion of bears assumed to be missed on the transect line during Russian surveys (g(0)). Our point estimates are larger than, but of similar magnitude to, a recent estimate for the period 2008-2016 ([Formula: see text]; 95% CI 1,522-5,944) derived from an integrated population model applied to a slightly smaller area. Although a number of factors (e.g., equipment issues, differing platforms, low sample sizes, size of the study area relative to sampling effort) required us to make a number of assumptions to generate estimates, it establishes a useful lower bound for abundance, and suggests high spring polar bear densities on sea ice in Russian waters south of Wrangell Island. With future improvements, we suggest that springtime aerial surveys may represent a plausible avenue for studying abundance and distribution of polar bears and their prey over large, remote areas.

Modeling the demography of species providing extended parental care: A capture-recapture multievent model with a case study on polar bears (Ursus maritimus)

In species providing extended parental care, one or both parents care for altricial young over a period including more than one breeding season. We expect large parental investment and long-term dependency within family units to cause high variability in life trajectories among individuals with complex consequences at the population level. So far, models for estimating demographic parameters in free-ranging animal populations mostly ignore extended parental care, thereby limiting our understanding of its consequences on parents and offspring life histories.We designed a capture-recapture multievent model for studying the demography of species providing extended parental care. It handles statistical multiple-year dependency among individual demographic parameters grouped within family units, variable litter size, and uncertainty on the timing at offspring independence. It allows for the evaluation of trade-offs among demographic parameters, the influence of past reproductive history on the caring parent's survival status, breeding probability, and litter size probability, while accounting for imperfect detection of family units. We assess the model performance using simulated data and illustrate its use with a long-term dataset collected on the Svalbard polar bears (Ursus maritimus).Our model performed well in terms of bias and mean square error and in estimating demographic parameters in all simulated scenarios, both when offspring departure probability from the family unit occurred at a constant rate or varied during the field season depending on the date of capture. For the polar bear case study, we provide estimates of adult and dependent offspring survival rates, breeding probability, and litter size probability. Results showed that the outcome of the previous reproduction influenced breeding probability.Overall, our results show the importance of accounting for i) the multiple-year statistical dependency within family units, ii) uncertainty on the timing at offspring independence, and iii) past reproductive history of the caring parent. If ignored, estimates obtained for breeding probability, litter size, and survival can be biased. This is of interest in terms of conservation because species providing extended parental care are often long-living mammals vulnerable or threatened with extinction.

Physiological consequences of Arctic sea ice loss on large marine carnivores: unique responses by polar bears and narwhals

Rapid environmental changes in the Arctic are threatening the survival of marine species that rely on the predictable presence of the sea ice. Two Arctic marine mammal specialists, the polar bear (Ursus maritimus) and narwhal (Monodon monoceros), appear especially vulnerable to the speed and capriciousness of sea ice deterioration as a consequence of their unique hunting behaviors and diet, as well as their physiological adaptations for slow-aerobic exercise. These intrinsic characteristics limit the ability of these species to respond to extrinsic threats associated with environmental change and increased industrial activity in a warming Arctic. In assessing how sea ice loss may differentially affect polar bears that hunt on the ice surface and narwhals that hunt at extreme depths below, we found that major ice loss translated into elevated locomotor costs that range from 3- to 4-fold greater than expected for both species. For polar bears this instigates an energy imbalance from the combined effects of reduced caloric intake and increased energy expenditure. For narwhals, high locomotor costs during diving increase the risk of ice entrapment due to the unreliability of breathing holes. These species-specific physiological constraints and extreme reliance on the polar sea ice conspire to make these two marine mammal specialists sentinels of climate change within the Arctic marine ecosystem that may foreshadow rapid changes to the marine ecosystem.

A review of the combined effects of climate change and other local human stressors on the marine environment

Climate change (CC) is a key, global driver of change of marine ecosystems. At local and regional scales, other local human stressors (LS) can interact with CC and modify its effects on marine ecosystems. Understanding the response of the marine environment to the combined effects of CC and LS is crucial to inform marine ecosystem-based management and planning, yet our knowledge of the potential effects of such interactions is fragmented. At a global scale, we explored how cumulative effect assessments (CEAs) have addressed CC in the marine realm and discuss progress and shortcomings of current approaches. For this we conducted a systematic review on how CEAs investigated at different levels of biological organization ecological responses, functional aspects, and the combined effect of CC and HS. Globally, the effects of 52 LS and of 27 CC-related stressors on the marine environment have been studied in combination, such as industrial fisheries with change in temperature, or sea level rise with artisanal fisheries, marine litter, change in sediment load and introduced alien species. CC generally intensified the effects of LS at species level. At trophic groups and ecosystem levels, the effects of CC either intensified or mitigated the effects of other HS depending on the trophic groups or the environmental conditions involved, thus suggesting that the combined effects of CC and LS are context-dependent and vary among and within ecosystems. Our results highlight that large-scale assessments on the spatial interaction and combined effects of CC and LS remain limited. More importantly, our results strengthen the urgent need of CEAs to capture local-scale effects of stressors that can exacerbate climate-induced changes. Ultimately, this will allow identifying management measures that aid counteracting CC effects at relevant scales.

Simulating changes in polar bear subpopulation growth rate due to legacy persistent organic pollutants - Temporal and spatial trends

Although atmospheric concentrations of many conventional persistent organic pollutants (POPs) have decreased in the Arctic over the past few decades, levels of most POPs and mercury remain high since the 1990s or start to increase again in Arctic areas, especially polar bears. So far, studies generally focused on individual effects of POPs, and do not directly link POP concentrations in prey species to population-specific parameters. In this study we therefore aimed to estimate the effect of legacy POPs and mercury on population growth rate of nineteen polar bear subpopulations. We modelled population development in three scenarios, based on species sensitivity distributions (SSDs) derived for POPs based on ecotoxicity data for endothermic species. In the first scenario, ecotoxicity data for polar bears were based on the HC₅₀ (the concentration at which 50% of the species is affected). The other two scenarios were based on the HC₅ and HC₉₅. Considerable variation in effects of POPs could be observed among the scenarios. In our intermediate scenario, we predicted subpopulation decline for ten out of 15 polar bear subpopulations. The estimated population growth rate was least reduced in Gulf of Boothia and Foxe Basin. On average, PCB concentrations in prey (in μg/g toxic equivalency (TEQ)) posed the largest threat to polar bear subpopulations, with negative modelled population growth rates for the majority of subpopulations. We did not find a correlation between modelled population changes and monitored population trends for the majority of chemical-subpopulation combinations. Modelled population growth rates increased over time, implying a decreasing effect of PCBs, DDTs, and mercury. Polar bear subpopulations are reportedly still declining in four out of the seven subpopulations for which sufficient long-term monitoring data is available, as reported by the IUCN-PBSG. This implies that other emerging pollutants or other anthropogenic stressors may affect polar bear subpopulations.

Grizzly and polar bears as nonconsumptive cultural keystone species

Grizzly bears and polar bears often serve as ecological “flagship species” in conservation efforts, but although consumptively used in some areas and cultures they can also be important cultural keystone species even where not hunted. We extend the application of established criteria for defining cultural keystone species to also encompass species with which cultures have a primarily nonconsumptive relationship but that are nonetheless disproportionately important to well-being and identity. Grizzly bears in coastal British Columbia are closely linked to many Indigenous Peoples (including the Haíɫzaqv (Heiltsuk), Kitasoo/Xai’xais, and Nuxalk First Nations), where they are central to the identity, culture, and livelihoods of individuals, families, Chiefs, and Nations. Polar bears in Churchill, Manitoba, provide another example as a cultural keystone species for a mixed Indigenous and non-Indigenous community in which many of the livelihood benefits from the species are mediated by economic transactions in a globalized tourism market. We discuss context specificity and questions of equity in sharing of benefits from cultural keystone species. Our expanded definition of cultural keystone species gives broader recognition of the beyond-ecological importance of these species to Indigenous Peoples, which highlights the societal and ecological importance of Indigenous sovereignty and could facilitate the increased cross-cultural understanding critical to reconciliation.

Influence of sea ice dynamics on population energetics of Western Hudson Bay polar bears

The Arctic marine ecosystem has experienced extensive changes in sea ice dynamics, with significant effects on ice-dependent species such as polar bears (Ursus maritimus). We used annual estimates of the numbers of bears onshore in the core summering area, age/sex structure and body condition data to estimate population energy density and storage energy in Western Hudson Bay polar bears from 1985 to 2018. We examined intra-population variation in energetic patterns, temporal energetic trends and the relationship between population energetics and sea ice conditions. Energy metrics for most demographic classes declined over time in relation to earlier sea ice breakup, most significantly for solitary adult females and subadult males, suggesting their greater vulnerability to nutritional stress than other age/sex classes. Temporal declines in population energy metrics were related to earlier breakup and longer lagged open-water periods, suggesting multi-year effects of sea ice decline. The length of the open-water period ranged from 102 to 166 days and increased significantly by 9.9 days/decade over the study period. Total population energy density and storage energy were significantly lower when sea ice breakup occurred earlier and the lagged open-water period was longer. At the earliest breakup and a lagged open-water period of 180 days, population energy density was predicted to be 33% lower than our minimum estimated energy density and population storage energy was predicted to be 40% lower than the minimum estimated storage energy. Consequently, over the study, the total population energy density declined by 53% (mean: 3668 ± 386 MJ kg-1/decade) and total population storage energy declined by 56% (mean: 435900 ± 46770 MJ/decade). This study provides insights into ecological mechanisms linking population responses to sea ice decline and highlights the significance of maintaining long-term research programs.

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.

Drivers and consequences of apex predator diet composition in the Canadian Beaufort Sea

Polar bears (Ursus maritimus) rely on annual sea ice as their primary habitat for hunting marine mammal prey. Given their long lifespan, wide geographic distribution, and position at the top of the Arctic marine food web, the diet composition of polar bears can provide insights into temporal and spatial ecosystem dynamics related to climate-mediated sea ice loss. Polar bears with the greatest ecological constraints on diet composition may be most vulnerable to climate-related changes in ice conditions and prey availability. We used quantitative fatty acid signature analysis (QFASA) to estimate the diets of polar bears (n = 419) in two western Canadian Arctic subpopulations (Northern Beaufort Sea and Southern Beaufort Sea) from 1999 to 2015. Polar bear diets were dominated by ringed seal (Pusa hispida), with interannual, seasonal, age- and sex-specific variation. Foraging area and sea ice conditions also affected polar bear diet composition. Most variation in bear diet was explained by longitude, reflecting spatial variation in prey availability. Sea ice conditions (extent, thickness, and seasonal duration) declined throughout the study period, and date of sea ice break-up in the preceding spring was positively correlated with female body condition and consumption of beluga whale (Delphinapterus leucas), suggesting that bears foraged on beluga whales during entrapment events. Female body condition was positively correlated with ringed seal consumption, and negatively correlated with bearded seal consumption. This study provides insights into the complex relationships between declining sea ice habitat and the diet composition and foraging success of a wide-ranging apex predator.

The stress of Arctic warming on polar bears

Arctic ecosystems are changing rapidly in response to climate warming. While Arctic mammals are highly evolved to these extreme environments, particularly with respect to their stress axis, some species may have limited capacity to adapt to this change. We examined changes in key components of the stress axis (cortisol and its carrier protein-corticosteroid binding globulin [CBG]) in polar bears (Ursus maritimus) from western Hudson Bay (N = 300) over a 33 year period (1983-2015) during which time the ice-free period was increasing. Changing sea ice phenology limits spring hunting opportunities and extends the period of onshore fasting. We assessed the response of polar bears to a standardized stressor (helicopter pursuit, darting, and immobilization) during their onshore fasting period (late summer-autumn) and quantified the serum levels of the maximum corticosteroid binding capacity (MCBC) of CBG, the serum protein that binds cortisol strongly, and free cortisol (FC). We quantified bear condition (age, sex, female with cubs or not, fat condition), sea ice (breakup in spring-summer, 1 year lagged freeze-up in autumn), and duration of fasting until sample collection as well as cumulative impacts of the latter environmental traits from the previous year. Data were separated into "good" years (1983-1990) when conditions were thought to be optimal and "poor" years (1991-2015) when sea ice conditions deteriorated and fasting on land was extended. MCBC explained 39.4% of the variation in the good years, but only 28.1% in the poor ones, using both biological and environmental variables. MCBC levels decreased with age. Changes in FC were complex, but more poorly explained. Counterintuitively, MCBC levels increased with increased time onshore, 1 year lag effects, and in poor ice years. We conclude that MCBC is a biomarker of stress in polar bears and that the changes we document are a consequence of climate warming.

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.

Contrasting Temporal Patterns of Mercury, Niche Dynamics, and Body Fat Indices of Polar Bears and Ringed Seals in a Melting Icescape

Polar bears (Ursus maritimus) and ringed seals (Pusa hispida) have a strong predator-prey relationship and are facing climate-associated Arctic habitat loss and harmful dietary exposure to total mercury (THg) and other pollutants. However, little is known about whether both species inhabiting the same area exhibit similar temporal patterns in Hg concentration, niche dynamics, and body fat indices. We used THg, δ¹³C, and δ¹⁵N values of western Hudson Bay polar bear hair (2004-2016) and ringed seal muscle samples (2003-2015) to investigate temporal trends of these variables and multidimensional niche metrics, as well as body fat indices for both species. We found a decline in THg concentration (by 3.8% per year) and δ¹³C (by 1.5‰) in ringed seals suggesting a change in feeding habits and carbon source use over time, whereas no significant changes occurred in polar bears. In contrast, the polar bear 3-dimensional niche size decreased by nearly half with no change in ringed seal niche size. The δ¹³C spacing between both species increased by approximately 1.5× suggesting different responses to annual changes in sympagic-pelagic carbon source production. Ringed seal body fat index was higher in years of earlier sea ice breakup with no change occurring in polar bears. These findings indicate that both species are responding differently to a changing environment suggesting a possible weakening of their predator-prey relationship in western Hudson Bay.

Correlates of seasonal change in the body condition of an Arctic top predator

Climate-driven sea ice loss has led to changes in the timing of key biological events in the Arctic, however, the consequences and rate of these changes remain largely unknown. Polar bears (Ursus maritimus) undergo seasonal changes in energy stores in relation to foraging opportunities and habitat conditions. Declining sea ice has been linked to reduced body condition in some subpopulations, however, the specific timing and duration of the feeding period when bears acquire most of their energy stores and its relationship to the timing of ice break-up is poorly understood. We used community-based sampling to investigate seasonality in body condition (energy stores) of polar bears in Nunavut, Canada, and examined the influence of sea ice variables. We used adipose tissue lipid content as an index of body condition for 1,206 polar bears harvested from 2010-2017 across five subpopulations with varying seasonal ice conditions: Baffin Bay (October-August), Davis Strait and Foxe Basin (year-round), Gulf of Boothia and Lancaster Sound (August-May). Similar seasonal patterns were found in body condition across subpopulations with bears at their nadir of condition in the spring, followed by fat accumulation past break-up date and subsequent peak body condition in autumn, indicating that bears are actively foraging in late spring and early summer. Late season feeding implies that even minor advances in the timing of break-up may have detrimental effects on foraging opportunities, body condition, and subsequent reproduction and survival. The magnitude of seasonal changes in body condition varied across the study area, presumably driven by local environmental conditions. Our results demonstrate how community-based monitoring of polar bears can reveal population-level responses to climate warming in advance of detectable demographic change. Our data on the seasonal timing of polar bear foraging and energy storage should inform predictive models of the effects of climate-mediated sea ice loss.

Modeling optimal responses and fitness consequences in a changing Arctic

Animals must balance a series of costs and benefits while trying to maximize their fitness. For example, an individual may need to choose how much energy to allocate to reproduction versus growth, or how much time to spend on vigilance versus foraging. Their decisions depend on complex interactions between environmental conditions, behavioral plasticity, reproductive biology, and energetic demands. As animals respond to novel environmental conditions caused by climate change, the optimal decisions may shift. Stochastic dynamic programming provides a flexible modeling framework with which to explore these trade-offs, but this method has not yet been used to study possible changes in optimal trade-offs caused by climate change. We created a stochastic dynamic programming model capturing trade-off decisions required by an individual adult female polar bear (Ursus maritimus) as well as the fitness consequences of her decisions. We predicted optimal foraging decisions throughout her lifetime as well as the energetic thresholds below which it is optimal for her to abandon a reproductive attempt. To explore the effects of climate change, we shortened the spring feeding period by up to 3 weeks, which led to predictions of riskier foraging behavior and higher reproductive thresholds. The resulting changes in fitness may be interpreted as a best-case scenario, where bears adapt instantaneously and optimally to new environmental conditions. If the spring feeding period was reduced by 1 week, her expected fitness declined by 15%, and if reduced by 3 weeks, expected fitness declined by 68%. This demonstrates an effective way to explore a species' optimal response to a changing landscape of costs and benefits and highlights the fact that small annual effects can result in large cumulative changes in expected lifetime fitness.

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.

State of knowledge on current exposure, fate and potential health effects of contaminants in polar bears from the circumpolar Arctic

The polar bear (Ursus maritimus) is among the Arctic species exposed to the highest concentrations of long-range transported bioaccumulative contaminants, such as halogenated organic compounds and mercury. Contaminant exposure is considered to be one of the largest threats to polar bears after the loss of their Arctic sea ice habitat due to climate change. The aim of this review is to provide a comprehensive summary of current exposure, fate, and potential health effects of contaminants in polar bears from the circumpolar Arctic required by the Circumpolar Action Plan for polar bear conservation. Overall results suggest that legacy persistent organic pollutants (POPs) including polychlorinated biphenyls, chlordanes and perfluorooctane sulfonic acid (PFOS), followed by other perfluoroalkyl compounds (e.g. carboxylic acids, PFCAs) and brominated flame retardants, are still the main compounds in polar bears. Concentrations of several legacy POPs that have been banned for decades in most parts of the world have generally declined in polar bears. Current spatial trends of contaminants vary widely between compounds and recent studies suggest increased concentrations of both POPs and PFCAs in certain subpopulations. Correlative field studies, supported by in vitro studies, suggest that contaminant exposure disrupts circulating levels of thyroid hormones and lipid metabolism, and alters neurochemistry in polar bears. Additionally, field and in vitro studies and risk assessments indicate the potential for adverse impacts to polar bear immune functions from exposure to certain contaminants.