Holocene deglaciation drove rapid genetic diversification of Atlantic walrus.

Rapid global warming is severely impacting Arctic ecosystems and is predicted to transform the abundance, distribution and genetic diversity of Arctic species, though these linkages are poorly understood. We address this gap in knowledge using palaeogenomics to examine how earlier periods of global warming influenced the genetic diversity of Atlantic walrus (Odobenus rosmarus rosmarus), a species closely associated with sea ice and shallow-water habitats. We analysed 82 ancient and historical Atlantic walrus mitochondrial genomes (mitogenomes), including now-extinct populations in Iceland and the Canadian Maritimes, to reconstruct the Atlantic walrus' response to Arctic deglaciation. Our results demonstrate that the phylogeography and genetic diversity of Atlantic walrus populations was initially shaped by the last glacial maximum (LGM), surviving in distinct glacial refugia, and subsequently expanding rapidly in multiple migration waves during the late Pleistocene and early Holocene. The timing of diversification and establishment of distinct populations corresponds closely with the chronology of the glacial retreat, pointing to a strong link between walrus phylogeography and sea ice. Our results indicate that accelerated ice loss in the modern Arctic may trigger further dispersal events, likely increasing the connectivity of northern stocks while isolating more southerly stocks putatively caught in small pockets of suitable habitat.

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.

A risk assessment review of mercury exposure in Arctic marine and terrestrial mammals.

There has been a considerable number of reports on Hg concentrations in Arctic mammals since the last Arctic Monitoring and Assessment Programme (AMAP) effort to review biological effects of the exposure to mercury (Hg) in Arctic biota in 2010 and 2018. Here, we provide an update on the state of the knowledge of health risk associated with Hg concentrations in Arctic marine and terrestrial mammal species. Using available population-specific data post-2000, our ultimate goal is to provide an updated evidence-based estimate of the risk for adverse health effects from Hg exposure in Arctic mammal species at the individual and population level. Tissue residues of Hg in 13 species across the Arctic were classified into five risk categories (from No risk to Severe risk) based on critical tissue concentrations derived from experimental studies on harp seals and mink. Exposure to Hg lead to low or no risk for health effects in most populations of marine and terrestrial mammals, however, subpopulations of polar bears, pilot whales, narwhals, beluga and hooded seals are highly exposed in geographic hotspots raising concern for Hg-induced toxicological effects. About 6% of a total of 3500 individuals, across different marine mammal species, age groups and regions, are at high or severe risk of health effects from Hg exposure. The corresponding figure for the 12 terrestrial species, regions and age groups was as low as 0.3% of a total of 731 individuals analyzed for their Hg loads. Temporal analyses indicated that the proportion of polar bears at low or moderate risk has increased in East/West Greenland and Western Hudson Bay, respectively. However, there remain numerous knowledge gaps to improve risk assessments of Hg exposure in Arctic mammalian species, including the establishment of improved concentration thresholds and upscaling to the assessment of population-level effects.

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.

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

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

Interrelated ecological impacts of climate change on an apex predator.

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

Immunologic, reproductive, and carcinogenic risk assessment from POP exposure in East Greenland polar bears (Ursus maritimus) during 1983-2013.

Polar bears (Ursus maritimus) are among the world's highest trophic level marine predators and as such have some of the highest tissue concentrations of organohalogen contaminants (OHCs) among Arctic biota. In this paper we present the results of a three decade (1983-2013) risk assessment of OHC exposure and effects on reproduction, immunity, and cancer (genotoxicity) in polar bears from Central East Greenland. Risk of adverse effects are evaluated using a risk quotient (RQ) approach with derivation from measured OHC concentrations in polar bear tissue and critical body residues (CBR) extrapolated for polar bears using physiologically-based pharmacokinetic modelling (PBPK). The additive RQs for all OHCs in polar bears were above the threshold for all effect categories (RQ > 1) in every year, suggesting this population has been at significant and continuous risk of contaminant-mediated effects for over three decades. RQs peaked in 1983 (RQ > 58) and again in 2013 (RQ > 50) after a period of decline. These trends follow ΣPCB levels during that time, and contributed almost all of the risk to immune, reproductive, and carcinogenic effects (71-99% of total RQ). The recent spike in RQs suggests a major shift in polar bear contaminant exposure from climate related changes in food composition and hereby the increased risk of adverse health effects. In the context of lifetime exposure ΣPCB and PFOS levels showed the interactive importance of year of birth, age, and emission history. In conclusion, the results indicate that East Greenland polar bears have been exposed to OHC levels over the period of 1983-2013 that potentially and continuously affected individual and theoretically also population health, with a peaking risk in the more recent years.

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

Climate change is expected to result in range shifts and habitat fragmentation for many species. In the Arctic, loss of sea ice will reduce barriers to dispersal or eliminate movement corridors, resulting in increased connectivity or geographic isolation with sweeping implications for conservation. We used satellite telemetry, data from individually marked animals (research and harvest), and microsatellite genetic data to examine changes in geographic range, emigration, and interpopulation connectivity of the Baffin Bay (BB) polar bear (Ursus maritimus) subpopulation over a 25-year period of sea-ice loss. Satellite telemetry collected from n = 43 (1991-1995) and 38 (2009-2015) adult females revealed a significant contraction in subpopulation range size (95% bivariate normal kernel range) in most months and seasons, with the most marked reduction being a 70% decline in summer from 716,000 km2 (SE 58,000) to 211,000 km2 (SE 23,000) (p < .001). Between the 1990s and 2000s, there was a significant shift northward during the on-ice seasons (2.6° shift in winter median latitude, 1.1° shift in spring median latitude) and a significant range contraction in the ice-free summers. Bears in the 2000s were less likely to leave BB, with significant reductions in the numbers of bears moving into Davis Strait (DS) in winter and Lancaster Sound (LS) in summer. Harvest recoveries suggested both short and long-term fidelity to BB remained high over both periods (83-99% of marked bears remained in BB). Genetic analyses using eight polymorphic microsatellites confirmed a previously documented differentiation between BB, DS, and LS; yet weakly differentiated BB from Kane Basin (KB) for the first time. Our results provide the first multiple lines of evidence for an increasingly geographically and functionally isolated subpopulation of polar bears in the context of long-term sea-ice loss. This may be indicative of future patterns for other polar bear subpopulations under climate change.

Arctic marine mammal population status, sea ice habitat loss, and conservation recommendations for the 21st century.

Arctic marine mammals (AMMs) are icons of climate change, largely because of their close association with sea ice. However, neither a circumpolar assessment of AMM status nor a standardized metric of sea ice habitat change is available. We summarized available data on abundance and trend for each AMM species and recognized subpopulation. We also examined species diversity, the extent of human use, and temporal trends in sea ice habitat for 12 regions of the Arctic by calculating the dates of spring sea ice retreat and fall sea ice advance from satellite data (1979-2013). Estimates of AMM abundance varied greatly in quality, and few studies were long enough for trend analysis. Of the AMM subpopulations, 78% (61 of 78) are legally harvested for subsistence purposes. Changes in sea ice phenology have been profound. In all regions except the Bering Sea, the duration of the summer (i.e., reduced ice) period increased by 5-10 weeks and by >20 weeks in the Barents Sea between 1979 and 2013. In light of generally poor data, the importance of human use, and forecasted environmental changes in the 21st century, we recommend the following for effective AMM conservation: maintain and improve comanagement by local, federal, and international partners; recognize spatial and temporal variability in AMM subpopulation response to climate change; implement monitoring programs with clear goals; mitigate cumulative impacts of increased human activity; and recognize the limits of current protected species legislation.

Implications of the circumpolar genetic structure of polar bears for their conservation in a rapidly warming Arctic.

We provide an expansive analysis of polar bear (Ursus maritimus) circumpolar genetic variation during the last two decades of decline in their sea-ice habitat. We sought to evaluate whether their genetic diversity and structure have changed over this period of habitat decline, how their current genetic patterns compare with past patterns, and how genetic demography changed with ancient fluctuations in climate. Characterizing their circumpolar genetic structure using microsatellite data, we defined four clusters that largely correspond to current ecological and oceanographic factors: Eastern Polar Basin, Western Polar Basin, Canadian Archipelago and Southern Canada. We document evidence for recent (ca. last 1-3 generations) directional gene flow from Southern Canada and the Eastern Polar Basin towards the Canadian Archipelago, an area hypothesized to be a future refugium for polar bears as climate-induced habitat decline continues. Our data provide empirical evidence in support of this hypothesis. The direction of current gene flow differs from earlier patterns of gene flow in the Holocene. From analyses of mitochondrial DNA, the Canadian Archipelago cluster and the Barents Sea subpopulation within the Eastern Polar Basin cluster did not show signals of population expansion, suggesting these areas may have served also as past interglacial refugia. Mismatch analyses of mitochondrial DNA data from polar and the paraphyletic brown bear (U. arctos) uncovered offset signals in timing of population expansion between the two species, that are attributed to differential demographic responses to past climate cycling. Mitogenomic structure of polar bears was shallow and developed recently, in contrast to the multiple clades of brown bears. We found no genetic signatures of recent hybridization between the species in our large, circumpolar sample, suggesting that recently observed hybrids represent localized events. Documenting changes in subpopulation connectivity will allow polar nations to proactively adjust conservation actions to continuing decline in sea-ice habitat.

Population genomics reveal recent speciation and rapid evolutionary adaptation in polar bears.

Polar bears are uniquely adapted to life in the High Arctic and have undergone drastic physiological changes in response to Arctic climates and a hyper-lipid diet of primarily marine mammal prey. We analyzed 89 complete genomes of polar bear and brown bear using population genomic modeling and show that the species diverged only 479-343 thousand years BP. We find that genes on the polar bear lineage have been under stronger positive selection than in brown bears; nine of the top 16 genes under strong positive selection are associated with cardiomyopathy and vascular disease, implying important reorganization of the cardiovascular system. One of the genes showing the strongest evidence of selection, APOB, encodes the primary lipoprotein component of low-density lipoprotein (LDL); functional mutations in APOB may explain how polar bears are able to cope with life-long elevated LDL levels that are associated with high risk of heart disease in humans.

Global change effects on the long-term feeding ecology and contaminant exposures of East Greenland polar bears.

Rapid climate changes are occurring in the Arctic, with substantial repercussions for arctic ecosystems. It is challenging to assess ecosystem changes in remote polar environments, but one successful approach has entailed monitoring the diets of upper trophic level consumers. Quantitative fatty acid signature analysis (QFASA) and fatty acid carbon isotope (δ(13) C-FA) patterns were used to assess diets of East Greenland (EG) polar bears (Ursus maritimus) (n = 310) over the past three decades. QFASA-generated diet estimates indicated that, on average, EG bears mainly consumed arctic ringed seals (47.5 ± 2.1%), migratory subarctic harp (30.6 ± 1.5%) and hooded (16.7 ± 1.3%) seals and rarely, if ever, consumed bearded seals, narwhals or walruses. Ringed seal consumption declined by 14%/decade over 28 years (90.1 ± 2.5% in 1984 to 33.9 ± 11.1% in 2011). Hooded seal consumption increased by 9.5%/decade (0.0 ± 0.0% in 1984 to 25.9 ± 9.1% in 2011). This increase may include harp seal, since hooded and harp seal FA signatures were not as well differentiated relative to other prey species. Declining δ(13) C-FA ratios supported shifts from more nearshore/benthic/ice-associated prey to more offshore/pelagic/open-water-associated prey, consistent with diet estimates. Increased hooded seal and decreased ringed seal consumption occurred during years when the North Atlantic Oscillation (NAO) was lower. Thus, periods with warmer temperatures and less sea ice were associated with more subarctic and less arctic seal species consumption. These changes in the relative abundance, accessibility, or distribution of arctic and subarctic marine mammals may have health consequences for EG polar bears. For example, the diet change resulted in consistently slower temporal declines in adipose levels of legacy persistent organic pollutants, as the subarctic seals have higher contaminant burdens than arctic seals. Overall, considerable changes are occurring in the EG marine ecosystem, with consequences for contaminant dynamics.

Three decades (1983-2010) of contaminant trends in East Greenland polar bears (Ursus maritimus). Part 1: legacy organochlorine contaminants.

Legacy organochlorine contaminants were determined in adipose tissues from 294 polar bears (Ursus maritimus) sampled in East Greenland in 23 of the 28years between 1983 and 2010. Of 19 major legacy contaminants and congeners (ΣPCB, 4 PCB congeners (CB153, 180, 170/190), ΣDDT, p,p'-DDE, p,p' -DDD and p,p'-DDT, α- and ß-hexachlorocyclohexane (HCH), HCB, octachlorostyrene, dieldrin, oxychlordane, cis- and trans-chlordane, cis- and trans-nonachlor, heptachlor epoxide and BB-153), 18 showed statistically significant average yearly declines of -4.4% (range: -2.0 to -10.8%/year) among subadult polar bears (i.e. females<5years, males<6years). For example, the ∑PCB concentrations declined 2.7 fold from 22730ng/g lw (95% C.I.: 12470-32990ng/g lw) in 1983-1986 to 8473ng/g lw (95% C.I.: 6369-9776ng/g lw) in 2006-2010. Similar but fewer statistically significant trends were found for adult females and adult males likely due to smaller sample size and years. Despite declines as a result of international regulations, relatively high levels of these historic pollutants persist in East Greenland polar bear tissues.

Three decades (1983-2010) of contaminant trends in East Greenland polar bears (Ursus maritimus). Part 2: brominated flame retardants.

Brominated flame retardants were determined in adipose tissues from 294 polar bears (Ursus maritimus) sampled in East Greenland in 23 of the 28years between 1983 and 2010. Significant linear increases were found for sum polybrominated diphenyl ether (ΣPBDE), BDE100, BDE153, and hexabromocyclododecane (HBCD). Average increases of 5.0% per year (range: 2.9-7.6%/year) were found for the subadult polar bears. BDE47 and BDE99 concentrations did not show a significant linear trend over time, but rather a significant non-linear trend peaking between 2000 and 2004. The average ΣPBDE concentrations increased 2.3 fold from 25.0ng/g lw (95% C.I.: 15.3-34.7ng/g lw) in 1983-1986 to 58.5ng/g lw (95% C.I.: 43.6-73.4ng/g lw) in 2006-2010. Similar but fewer statistically significant trends were found for adult females and adult males likely due to smaller sample size and years. Analyses of δ(15)N and δ(13)C stable isotopes in hair revealed no clear linear temporal trends in trophic level or carbon source, respectively, and non-linear trends differed among sex and age groups. These increasing concentrations of organobromine contaminants contribute to complex organohalogen mixture, already causing health effects to the East Greenland polar bears.

What are the toxicological effects of mercury in Arctic biota?

This review critically evaluates the available mercury (Hg) data in Arctic marine biota and the Inuit population against toxicity threshold values. In particular marine top predators exhibit concentrations of mercury in their tissues and organs that are believed to exceed thresholds for biological effects. Species whose concentrations exceed threshold values include the polar bears (Ursus maritimus), beluga whale (Delphinapterus leucas), pilot whale (Globicephala melas), hooded seal (Cystophora cristata), a few seabird species, and landlocked Arctic char (Salvelinus alpinus). Toothed whales appear to be one of the most vulnerable groups, with high concentrations of mercury recorded in brain tissue with associated signs of neurochemical effects. Evidence of increasing concentrations in mercury in some biota in Arctic Canada and Greenland is therefore a concern with respect to ecosystem health.

Females roam while males patrol: divergence in breeding season movements of pack-ice polar bears (Ursus maritimus).

Intraspecific differences in movement behaviour reflect different tactics used by individuals or sexes to favour strategies that maximize fitness. We report movement data collected from n = 23 adult male polar bears with novel ear-attached transmitters in two separate pack ice subpopulations over five breeding seasons. We compared movements with n = 26 concurrently tagged adult females, and analysed velocities, movement tortuosity, range sizes and habitat selection with respect to sex, reproductive status and body mass. There were no differences in 4-day displacements or sea ice habitat selection for sex or population. By contrast, adult females in all years and both populations had significantly more linear movements and significantly larger breeding range sizes than males. We hypothesized that differences were related to encounter rates, and used observed movement metrics to parametrize a simulation model of male-male and male-female encounter. The simulation showed that the more tortuous movement of males leads to significantly longer times to male-male encounter, while having little impact on male-female encounter. By contrast, linear movements of females are consistent with a prioritized search for sparsely distributed prey. These results suggest a possible mechanism for explaining the smaller breeding range sizes of some solitary male carnivores compared to females.

Tissue-specific concentrations and patterns of perfluoroalkyl carboxylates and sulfonates in East Greenland polar bears.

Several perfluoroalkyl carboxylates (PFCAs) and perfluoroalkyl sulfonates (PFSAs) of varying chain length are bioaccumulative in biota. However, wildlife reports have focused on liver and with very little examination of other tissues, and thus there is a limited understanding of their distribution and potential effects in the mammalian body. In the present study, the comparative accumulation of C(6) to C(15) PFCAs, C(4), C(6), C(8) and C(10) PFSAs, and select precursors were examined in the liver, blood, muscle, adipose, and brain of 20 polar bears (Ursus maritimus) from Scoresby Sound, Central East Greenland. Overall, PFSA and PFCA concentrations were highest in liver followed by blood > brain > muscle ≈ adipose. Liver and blood samples contained proportionally more of the shorter/medium chain length (C(6) to C(11)) PFCAs, whereas adipose and brain samples were dominated by longer chain (C(13) to C(15)) PFCAs. PFCAs with lower lipophilicities accumulated more in the liver, whereas the brain accumulated PFCAs with higher lipophilicities. The concentration ratios (±SE) between perfluorooctane sulfonate and its precursor perfluorooctane sulfonamide varied among tissues from 9 (±1):1 (muscle) to 36 (±7):1 (liver). PFCA and PFSA patterns in polar bears indicate that the pharmacokinetics of these compounds are to some extent tissue-specific, and are the result of several factors that may include differing protein interactions throughout the body.

Influence of carbon and lipid sources on variation of mercury and other trace elements in polar bears (Ursus maritimus).

In the present study, the authors investigated the influence of carbon and lipid sources on regional differences in liver trace element (As, Cd, Cu, total Hg, Mn, Pb, Rb, Se, and Zn) concentrations measured in polar bears (Ursus maritimus) (n = 121) from 10 Alaskan, Canadian Arctic, and East Greenland subpopulations. Carbon and lipid sources were assessed using δ(13) C in muscle tissue and fatty acid (FA) profiles in subcutaneous adipose tissue as chemical tracers. A negative relationship between total Hg and δ(13) C suggested that polar bears feeding in areas with higher riverine inputs of terrestrial carbon accumulate more Hg than bears feeding in areas with lower freshwater input. Mercury concentrations were also positively related to the FA 20:1n-9, which is biosynthesized in large amounts in Calanus copepods. This result raises the hypothesis that Calanus glacialis are an important link in the uptake of Hg in the marine food web and ultimately in polar bears. Unadjusted total Hg, Se, and As concentrations showed greater geographical variation among polar bear subpopulations compared with concentrations adjusted for carbon and lipid sources. The Hg concentrations adjusted for carbon and lipid sources in Bering-Chukchi Sea polar bear liver tissue remained the lowest among subpopulations. Based on these findings, the authors suggest that carbon and lipid sources for polar bears should be taken into account when one is assessing spatial and temporal trends of long-range transported trace elements.

Spatial and temporal trends of selected trace elements in liver tissue from polar bears (Ursus maritimus) from Alaska, Canada and Greenland.

Spatial trends and comparative changes in time of selected trace elements were studied in liver tissue from polar bears from ten different subpopulation locations in Alaska, Canadian Arctic and East Greenland. For nine of the trace elements (As, Cd, Cu, Hg, Mn, Pb, Rb, Se and Zn) spatial trends were investigated in 136 specimens sampled during 2005-2008 from bears from these ten subpopulations. Concentrations of Hg, Se and As were highest in the (northern and southern) Beaufort Sea area and lowest in (western and southern) Hudson Bay area and Chukchi/Bering Sea. In contrast, concentrations of Cd showed an increasing trend from east to west. Minor or no spatial trends were observed for Cu, Mn, Rb and Zn. Spatial trends were in agreement with previous studies, possibly explained by natural phenomena. To assess temporal changes of Cd, Hg, Se and Zn concentrations during the last decades, we compared our results to previously published data. These time comparisons suggested recent Hg increase in East Greenland polar bears. This may be related to Hg emissions and/or climate-induced changes in Hg cycles or changes in the polar bear food web related to global warming. Also, Hg:Se molar ratio has increased in East Greenland polar bears, which suggests there may be an increased risk for Hg(2+)-mediated toxicity. Since the underlying reasons for spatial trends or changes in time of trace elements in the Arctic are still largely unknown, future studies should focus on the role of changing climate and trace metal emissions on geographical and temporal trends of trace elements.