Dietary habits of polar bears in Foxe Basin, Canada: possible evidence of a trophic regime shift mediated by a new top predator

Dietary plasticity and broad North Atlantic origins inferred from bulk and amino acid-specific δ15N and δ13C favour killer whale range expansions into Arctic waters

Killer whales (Orcinus orca) occur seasonally in the eastern Canadian Arctic (ECA), where their range expansion associated with declining sea ice have raised questions about the impacts of increasing killer whale predation pressure on Arctic-endemic prey. We assessed diet and distribution of ECA killer whales using bulk and compound-specific stable isotope analysis (CSIA) of amino acids (AA) of 54 skin biopsies collected from 2009 to 2020 around Baffin Island, Canada. Bulk ECA killer whale skin δ15N and δ13C values did not overlap with potential Arctic prey after adjustment for trophic discrimination, and instead reflected foraging history in the North Atlantic prior to their arrival in the ECA. Adjusted killer whale stable isotope (SI) values primarily overlapped with several species of North Atlantic baleen whales or tuna. Amino acid (AA)-specific δ15N values indicated the ECA killer whales fed primarily on marine mammals, having similar glutamic acid δ15N-phenylalanine δ15N (δ15NGlx-Phe) and threonine δ15N (δ15NThr) as mammal-eating killer whales from the eastern North Pacific (ENP) that served as a comparative framework. However, one ECA whale grouped with the fish-eating ENP ecotype based δ15NThr. Distinctive essential AA δ13C of ECA killer whale groups, along with bulk SI similarity to killer whales from different regions of the North Atlantic, indicates different populations converge in Arctic waters from a broad source area. Generalist diet and long-distance dispersal capacity favour range expansions, and integration of these insights will be critical for assessing ecological impacts of increasing killer whale predation pressure on Arctic-endemic species.

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.

Measuring repeatability of compositional diet estimates: An example using quantitative fatty acid signature analysis

By measuring the temporal consistency, or repeatability, in the diets of predators, we can gain a better understanding of the degree of individual specialization in resource utilization and implications for predator-prey interactions, population dynamics, and food web structure. To measure repeatability, we require repeated diet estimates of individuals over time, such as those derived from quantitative fatty acid signature analysis (QFASA), a popular diet estimation technique. However, diet estimates are often lengthy compositional vectors with many zeros, as some prey will not be consumed by all individuals, precluding the use of previously proposed measures of repeatability. In this paper, we propose a novel approach for inferring repeatability for multivariate data and, in particular, compositional diet estimates. We extend the commonly used measure of repeatability for univariate data to the multivariate compositional setting by utilizing the mean squares obtained from a nonparametric multivariate analysis of variance, and an appropriate choice of statistical distance. Our measure and its extension are compatible with both balanced and unbalanced data sets. Associated confidence intervals via nonparametric bootstrapping are also developed for the case of QFASA diet estimates that incorporate both sampling error and measurement error, where the latter error arises because the diets of predators are estimated. Simulation study results suggest that for practical levels of repeatability, our methods yield confidence intervals with the desired coverage probability even when the sample size relative to the dimension of the data (i.e., number of prey species eaten) is small. We tested our methods using QFASA diet estimates for free-ranging Northwest Atlantic grey seals. Given the importance of understanding how predator diets vary over time and space, our method may find broad application to other compositional diet estimates, including those derived from the stomach or fecal contents, and stable isotope analyses.

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.

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.

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 < 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.

Set sustainable goals for the Arctic gateway coordinated international governance is required to resist yet another tipping point

Global warming is reducing the Arctic sea-ice and causing energetic stress to marine key predatory species such as polar bears and narwhals contributing to the ongoing pollution already threatening the biodiversity and indigenous people of the vulnerable region. Now, the opening of the Arctic gateway and in particular the increase in shipping activities causes further stress to marine mammals in the region. These shipping activities are foreseen to happen in the Northwest and Northeast Passage, Northern Sea Route and Transpolar Sea Route in the Arctic Ocean, which could be yet another step towards a crucial tipping point destabilizing global climate, including weathering systems and sea-level rise. This calls for international governance through the establishment of Arctic International National Parks and more Marine Protected Areas through the Arctic Council and UN's Law of the Sea to ensure sustainable use of the Arctic Ocean and adjacent waters.

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.

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.

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.

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.

Fatty acid-based diet estimates suggest ringed seal remain the main prey of southern Beaufort Sea polar bears despite recent use of onshore food resources

Polar bears (Ursus maritimus) from the southern Beaufort Sea (SB) subpopulation have traditionally fed predominantly upon ice-seals; however, as the proportion of the subpopulation using onshore habitat has recently increased, foraging on land-based resources, including remains of subsistence-harvested bowhead whales (Balaena mysticetus) and colonial nesting seabirds has been observed. Adipose tissue samples were collected from this subpopulation during the springs of 2013-2016 and analyzed for fatty acid signatures. Diet estimates were generated for the proportional consumption of ringed seal (Pusa hispida), bearded seal (Erignathus barbatus), and beluga whale (Delphinapterus leucas), relative to onshore foods, including bowhead whale remains and seabird, as represented by black guillemot (Cepphus grylle mandtii) nestlings and eggs. Quantitative fatty acid signature analysis (QFASA) estimated that the ice-obligate prey, ringed seal, remained the predominant prey species of SB polar bears (46.4 ± 1.8%), with much lower consumption of bearded seal (19.6 ± 2.0%), seabird (17.0 ± 1.2%), bowhead whale (15.0 ± 1.4%), and hardly any beluga whale (2.0 ± 0.5%). Adult and subadult females appeared to depend more on the traditional ringed seal prey than adult and subadult males. Diet estimates of SB polar bears showed significant interannual variability for all prey (F 12, 456 = 3.17, p < .001). Longer-term estimates suggested that both types of onshore prey, bowhead whale remains and seabird, have represented a moderate proportion of the food resources used by SB polar bears since at least the start of the 21st Century.

Preservation of fatty acid signatures in three vertebrate species after six months of storage at various temperatures

Fatty acid (FA) signatures (FAS) are important tools to assess the foraging ecology of wild animals. The present study was conducted to assess how well the general FAS and the proportions of individual FA are preserved in fat samples stored at different temperatures (–196, –80, –20, +4 and +20°C). Using three species (laboratory rat, American mink and rainbow trout), FAS were determined immediately upon sampling. Thereafter, eight subsamples per storage temperature from the inner part of the sample unaffected by oxygen and light were re-analyzed after 1, 2, 3, 7, 28, 84 and 168 days. Each time the remaining sample was sealed in its vial after replacing air with nitrogen gas. The results were tested with the mixed model and discriminant analyses. Generally, the FAS were well preserved regardless of storage temperature, and only a few major FA showed significant changes even after the 6-month period at room temperature. After an initial first-day change in proportions, presumably due to post-mortem enzymatic activities, the remaining minor changes could not be clearly attributed to either further autolysis, decomposition or autoxidation. In the discriminant analysis, the species-specific differences dominated and remained distinct even after 6 months. Furthermore, the analysis mostly classified the samples preserved at sub- and above-freezing temperatures separate from each other, and the general deviation from the initial analysis results was present as early as after 1 day. If FAS are to be analyzed in a very precise manner, the analysis should be performed immediately upon sampling. However, FAS remain adequately reliable for long periods of time even without preservation in deep freeze, widening the availability of potential samples for studies on foraging ecology and related disciplines.

Effects of biometrics, location and persistent organic pollutants on blood clinical-chemical parameters in polar bears (Ursus maritimus) from Svalbard, Norway

In the present study, blood clinical-chemical parameters (BCCPs) were analysed in 20 female and 18 male Svalbard polar bears (Ursus maritimus) captured in spring 2007. The aim was to study how age, body condition (BC), biometrics, plasma lipid content and geographical location may confound the relationship between persistent organic pollutants (POPs) including PCBs, HCB, chlordanes, DDTs, HCHs, mirex and OH-PCBs and the concentrations of 12 specific BCCPs (hematocrit [HCT], hemoglobin [HB], aspartate aminotransferase [ASAT], alanine aminotransferase [ALAT], γ-glutamyltransferase [GGT], creatine kinase [CK], triglycerides [TG], cholesterol [CHOL], high-density lipoprotein [HDL], creatinine (CREA], urea, potassium (K]), and to investigate if any of these BCCPs may be applied as potential biomarkers for POP exposure in polar bears. Initial PCA and O-PLS modelling showed that age, lipids, BC and geographical location (longitude and latitude) were important parameters explaining BCCPs in females. Following subsequent partial correlation analyses correcting for age and lipids, multiple POPs in females were still significantly correlated with HCT and HDL (all p < 0.05). In males, age, BM, BC and longitude were important parameters explaining BCCPs. Following partial correlation analyses correcting for age, biometrics, lipids and longitude in males, multiple POPs were significantly correlated with HCT, ASAT, GGT and CHOL (all p < 0.05). In conclusion, several confounding parameters has to be taken into account when studying the relations between BCCPs and POPs in polar bears. When correcting for these, in particular HCT may be used as a simple cost-efficient biomarker of POP exposure in polar bears. Furthermore, decreasing HDL concentrations and increasing CHOL concentration with increasing POP concentrations may indicate responses related to increased risk of cardiovascular disease. We therefore suggest to further study POP exposure and lipidome response to increase knowledge of the risk of cardiometabolic syndrome in polar bears.

High contributions of sea ice derived carbon in polar bear (Ursus maritimus) tissue

Polar bears (Ursus maritimus) rely upon Arctic sea ice as a physical habitat. Consequently, conservation assessments of polar bears identify the ongoing reduction in sea ice to represent a significant threat to their survival. However, the additional role of sea ice as a potential, indirect, source of energy to bears has been overlooked. Here we used the highly branched isoprenoid lipid biomarker-based index (H-Print) approach in combination with quantitative fatty acid signature analysis to show that sympagic (sea ice-associated), rather than pelagic, carbon contributions dominated the marine component of polar bear diet (72–100%; 99% CI, n = 55), irrespective of differences in diet composition. The lowest mean estimates of sympagic carbon were found in Baffin Bay bears, which were also exposed to the most rapidly increasing open water season. Therefore, our data illustrate that for future Arctic ecosystems that are likely to be characterised by reduced sea ice cover, polar bears will not only be impacted by a change in their physical habitat, but also potentially in the supply of energy to the ecosystems upon which they depend. This data represents the first quantifiable baseline that is critical for the assessment of likely ongoing changes in energy supply to Arctic predators as we move into an increasingly uncertain future for polar ecosystems.

Ecological Change Drives a Decline in Mercury Concentrations in Southern Beaufort Sea Polar Bears

We evaluated total mercury (THg) concentrations and trends in polar bears from the southern Beaufort Sea subpopulation from 2004 to 2011. Hair THg concentrations ranged widely among individuals from 0.6 to 13.3 μg g^-1 dry weight (mean: 3.5 ± 0.2 μg g^-1). Concentrations differed among sex and age classes: solitary adult females ≈ adult females with cubs ≈ subadults > adult males ≈ yearlings > cubs-of-the-year ≈ 2 year old dependent cubs. No variation was observed between spring and fall samples. For spring-sampled adults, THg concentrations declined by 13% per year, contrasting recent trends observed for other Western Hemispheric Arctic biota. Concentrations also declined by 15% per year considering adult males only, while a slower, nonsignificant decrease of 4.4% per year was found for adult females. Lower THg concentrations were associated with higher body mass index (BMI) and higher proportions of lower trophic position food resources consumed. Because BMI and diet were related, and the relationship to THg was strongest for BMI, trends were re-evaluated adjusting for BMI as the covariate. The adjusted annual decline was not significant. These findings indicate that changes in foraging ecology, not declining environmental concentrations of mercury, are driving short-term declines in THg concentrations in southern Beaufort Sea polar bears.

Simultaneous estimation of diet composition and calibration coefficients with fatty acid signature data

Knowledge of animal diets provides essential insights into their life history and ecology, although diet estimation is challenging and remains an active area of research. Quantitative fatty acid signature analysis (QFASA) has become a popular method of estimating diet composition, especially for marine species. A primary assumption of QFASA is that constants called calibration coefficients, which account for the differential metabolism of individual fatty acids, are known. In practice, however, calibration coefficients are not known, but rather have been estimated in feeding trials with captive animals of a limited number of model species. The impossibility of verifying the accuracy of feeding trial derived calibration coefficients to estimate the diets of wild animals is a foundational problem with QFASA that has generated considerable criticism. We present a new model that allows simultaneous estimation of diet composition and calibration coefficients based only on fatty acid signature samples from wild predators and potential prey. Our model performed almost flawlessly in four tests with constructed examples, estimating both diet proportions and calibration coefficients with essentially no error. We also applied the model to data from Chukchi Sea polar bears, obtaining diet estimates that were more diverse than estimates conditioned on feeding trial calibration coefficients. Our model avoids bias in diet estimates caused by conditioning on inaccurate calibration coefficients, invalidates the primary criticism of QFASA, eliminates the need to conduct feeding trials solely for diet estimation, and consequently expands the utility of fatty acid data to investigate aspects of ecology linked to animal diets.