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

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

Population genetic structure and evolutionary history of North Atlantic beluga whales (Delphinapterus leucas) from West Greenland, Svalbard and the White Sea

Population structure in many Arctic marine mammal species reflects a dynamic interplay between physical isolating mechanisms and the extent to which dispersal opportunities are met. We examined variation within mtDNA and eight microsatellite markers to investigate population structure and demographic history in beluga whales in the North Atlantic. Genetic heterogeneity was observed between Svalbard and West Greenland that reveals limited gene flow over ecological time scales. Differentiation was also recorded between Atlantic belugas and two previously studied populations in the North Pacific, the Beaufort Sea and Gulf of Alaska. However, Bayesian cluster analysis of the nDNA data identified two population clusters that did not correspond to the respective ocean basins, as predicted, but to: (1) Arctic (Svalbard–White Sea–Greenland–Beaufort Sea) and (2) Subarctic (Gulf of Alaska) regions. Similarly, the deepest phylogeographic signal was between the Arctic populations and the Gulf of Alaska. Fitting an isolation-with-migration model yielded genetic abundance estimates that match census estimates and revealed that Svalbard and the Beaufort Sea likely diverged 7,600–35,400 years ago but have experienced recurrent periods with gene flow since then, most likely via the Russian Arctic during subsequent warm periods. Considering current projections of continued sea ice losses in the Arctic, this study identified likely routes of future contact among extant beluga populations, and other mobile marine species, which have implications for genetic introgression, health, ecology and behavior.     

Ecosystem regime shifts have not affected growth and survivorship of eastern Beaufort Sea belugas

Large-scale ocean-atmosphere physical dynamics can have profound impacts on the structure and organization of marine ecosystems. These changes have been termed “regime shifts”, and five different episodes have been detected in the North Pacific Ocean, with concurrent changes also occurring in the Bering and Beaufort Seas. Belugas from the Eastern Beaufort Sea (EBS) use the Bering Sea during winter and the Beaufort Sea during summer, yet the potential effects of regime shifts on belugas have not been assessed. We investigated whether body size and survivorship of EBS belugas harvested in the Mackenzie River delta region between 1993 and 2003 have been affected by previous purported regime shifts in the North Pacific. Residuals from the relationship between body length and age were calculated and compared among belugas born between 1932 and 1989. Residual body size was not significantly related to birth year for any regime, nor to the age group individuals belonged to during any regime. The percentage deviation in number of belugas born in any given year that survived to be included in the hunt (survivorship) did not show any significant trend within or between regimes. Accounting for lags of 1–5 years did not reveal any evidence of delayed effects. Furthermore, neither population index was significantly related to changes in major climatic variables that precede regime shifts. Our results suggest that EBS beluga body size and survivorship have not been affected by the major regime shifts of the North Pacific and the adjacent Bering and Beaufort Seas. EBS belugas may have been able to modify their diet without compromising their growth and survivorship. Diet and reproductive analyses over large and small time scales can help understand the mechanisms enabling belugas to avoid significant growth and reproductive effects of past regime shifts. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Allelic and haplotype variation of major histocompatibility complex class II DRB1 and DQB loci in the St Lawrence beluga (Delphinapterus leucas)

In order to assess levels of major histocompatibility complex (Mhc) variation within the St Lawrence beluga (Delphinapterus leucas) the variation at the beluga Mhc DRB1 class II locus was assessed by single-strand conformation polymorphism (SSCP) analysis of the peptide-binding region for 313 whales collected from 13 sampling locations across North America. In addition, samples from west Greenland and the St Lawrence were also typed at the DQB locus, allowing comparison to a previous study and assessment of linkage disequilibrium of alleles at the two loci. Comparisons of DRB1 and DQB allele frequencies among all sampling locations indicated genetic structure (α < 0.005). Most of this structure resulted from differences between the different wintering groups. Significant genetic structure (α = 0.05) exists among each pair of the following groups at both the DRB1 and DQB loci; St Lawrence, Hudson Strait, Bering Sea, Cunningham Inlet, and Davis Strait (minus Cunningham Inlet), except the St Lawrence and Hudson Strait for the DQB locus. In the St Lawrence population, six of the eight DRB1 alleles are present representing all five known allelic lineages. Evidence of linkage disequilibrium between the DRB1 and DQB is present in two sampling locations, the St Lawrence and Nuussuaq (α = 0.05). Analysis of probable DRB1–DQB haplotypes among groups of beluga suggests a haplotype reduction in the St Lawrence.     

Sequence variation at the major histocompatibility complex locus DQ beta in beluga whales (Delphinapterus leucas) [published erratum appears in Mol Biol Evol 1995 Nov;12(6):1174]

Genetic variation at the Major Histocompatibility Complex locus DQ beta was analyzed in 233 beluga whales (Delphinapterus leucas) from seven populations: St. Lawrence Estuary, eastern Beaufort Sea, eastern Chukchi Sea, western Hudson Bay, eastern Hudson Bay, southeastern Baffin Island, and High Arctic and in 12 narwhals (Monodon monoceros) sympatric with the High Arctic beluga population. Variation was assessed by amplification of the exon coding for the peptide binding region via the polymerase chain reaction, followed by either cloning and DNA sequencing or single-stranded conformation polymorphism analysis. Five alleles were found across the beluga populations and one in the narwhal. Pairwise comparisons of these alleles showed a 5:1 ratio of nonsynonymous to synonymous substitutions per site leading to eight amino acid differences, five of which were nonconservative substitutions, centered around positions previously shown to be important for peptide binding. Although the amount of allelic variation is low when compared with terrestrial mammals, the nature of the substitutions in the peptide binding sites indicates an important role for the DQ beta locus in the cellular immune response of beluga whales. Comparisons of allele frequencies among populations show the High Arctic population to be different (P < or = .005) from the other beluga populations surveyed. In these other populations an allele, Dele-DQ beta*0101-2, was found in 98% of the animals, while in the High Arctic it was found in only 52% of the animals. Two other alleles were found at high frequencies in the High Arctic population, one being very similar to the single allele found in narwhal.      

Population structure of North American beluga whales (Delphinapterus leucas) based on nuclear DNA microsatellite variation and contrasted with the population structure revealed by mitochondrial DNA variation

Beluga whales ( Delphinapterus leucas ) in North American waters migrate seasonally between wintering areas in broken pack ice and summering locations in estuaries and other open water areas in the Arctic and sub-Arctic. Results from our previous investigation of beluga whale mitochondrial DNA (mtDNA) revealed genetic heterogeneity among beluga from different summering locations that was interpreted as representing a high degree of summering site philopatry. However, mtDNA is maternally inherited and does not reflect mating that may occur among beluga from different summering locations in wintering areas or during annual migrations. To test the possibility that breeding occurs among beluga from different summering locations, genetic variability at five nuclear DNA (nDNA) microsatellite loci was examined in the same animals tested in the mtDNA study. Beluga samples ( n = 640) were collected between 1984 and 1994 from 24 sites across North America, mostly during the summer. Whales from the various sites were categorized into eight summering locations as identified by mtDNA analysis, as well as four hypothesized wintering areas: Bering Sea, Hudson Strait (Hudson Strait, Labrador Sea, southwest Davis Strait), Baffin Bay (North Water, east Davis Strait), and St Lawrence River. Microsatellite allele frequencies indicated genetic homogeneity among animals from summering sites believed to winter together but differentiation among whales from some of the wintering areas. In particular, beluga from western North America (Bering Sea) were clearly distinguished from beluga from eastern North America (Hudson Strait, Baffin Bay, and St Lawrence River). Based upon the combined data set, the population of North American beluga whales was divided into two evolutionarily significant units. However, the population may be further subdivided into management units to reflect distinct groups of beluga at summering locations.     

Population structure of North American beluga whales (Delphinapterus leucas) based on nuclear DNA microsatellite variation and contrasted with the population structure revealed by mitochondrial DNA variation

Beluga whales ( Delphinapterus leucas ) in North American waters migrate seasonally between wintering areas in broken pack ice and summering locations in estuaries and other open water areas in the Arctic and sub-Arctic. Results from our previous investigation of beluga whale mitochondrial DNA (mtDNA) revealed genetic heterogeneity among beluga from different summering locations that was interpreted as representing a high degree of summering site philopatry. However, mtDNA is maternally inherited and does not reflect mating that may occur among beluga from different summering locations in wintering areas or during annual migrations. To test the possibility that breeding occurs among beluga from different summering locations, genetic variability at five nuclear DNA (nDNA) microsatellite loci was examined in the same animals tested in the mtDNA study. Beluga samples ( n = 640) were collected between 1984 and 1994 from 24 sites across North America, mostly during the summer. Whales from the various sites were categorized into eight summering locations as identified by mtDNA analysis, as well as four hypothesized wintering areas: Bering Sea, Hudson Strait (Hudson Strait, Labrador Sea, southwest Davis Strait), Baffin Bay (North Water, east Davis Strait), and St Lawrence River. Microsatellite allele frequencies indicated genetic homogeneity among animals from summering sites believed to winter together but differentiation among whales from some of the wintering areas. In particular, beluga from western North America (Bering Sea) were clearly distinguished from beluga from eastern North America (Hudson Strait, Baffin Bay, and St Lawrence River). Based upon the combined data set, the population of North American beluga whales was divided into two evolutionarily significant units. However, the population may be further subdivided into management units to reflect distinct groups of beluga at summering locations.     

Bone fluoride concentrations in beluga whales from Canada

Beluga whales (Delphinapterus leucas) from the St. Lawrence Estuary have been reported to have dental and bone abnormalities. To determine whether these lesions could be caused by high exposure to fluorides, we measured bone fluoride levels in eight beluga whales stranded on the shores of the St. Lawrence Estuary (Quebec, Canada), and in nine beluga whales killed by Inuit hunters in the Hudson Bay (North Western Territories, Canada). In both groups, fluoride concentrations were higher than those found in terrestrial mammals intoxicated by fluorides. Unexpectedly, fluoride concentration was significantly higher in beluga whales from the Hudson Bay (mean +/- SD: 10.365 +/- 1.098 ppm) than in beluga whales from the St. Lawrence Estuary (4.539 +/- 875 ppm) and was positively correlated with age in the latter population. Differences in diet might explain the differences in fluoride concentrations found between these two populations.     

Satellite telemetry reveals population specific winter ranges of beluga whales in the Bering Sea

At least five populations (stocks) of beluga whales (Delphinapterus leucas) are thought to winter in the Being Sea, including the Bristol Bay, Eastern Bering Sea (Norton Sound), Anadyr, Eastern Chukchi Sea, and Eastern Beaufort Sea (Mackenzie) populations. Belugas from each population have been tagged with satellite‐linked transmitters, allowing us to describe their winter (January–March) distribution. The objectives of this paper were to determine: (1) If each population winters in the Bering Sea, and if so, where? (2) Do populations return to the same area each year (i.e., are wintering areas traditional)? (3) To what extent do the winter ranges of different populations overlap? Tagged belugas from all five populations either remained in, or moved into, the Bering Sea and spent the winter there. Each population wintered in a different part of the Bering Sea and populations with multiple years of data (four of five) returned to the same regions in multiple years. When data were available from two populations that overlapped in the same year, they did not occupy the shared area at the same time. Although our sample sizes were small, the evidence suggests belugas from different populations have traditional winter ranges that are mostly exclusive to each population.     

Microsatellite analysis of population structure in Canadian polar bears

Attempts to study the genetic population structure of large mammals are often hampered by the low levels of genetic variation observed in these species. Polar bears have particularly low levels of genetic variation with the result that their genetic population structure has been intractable. We describe the use of eight hypervariable microsatellite loci to study the genetic relationships between four Canadian polar bear populations: the northern Beaufort Sea, southern Beaufort Sea, western Hudson Bay, and Davis Strait - Labrador Sea. These markers detected considerable genetic variation, with average heterozygosity near 60% within each population. Interpopulation differences in allele frequency distribution were significant between all pairs of populations, including two adjacent populations in the Beaufort Sea. Measures of genetic distance reflect the geographic distribution of populations, but also suggest patterns of gene flow which are not obvious from geography and may reflect movement patterns of these animals. Distribution of variation is sufficiently different between the Beaufort Sea populations and the two more eastern ones that the region of origin for a given sample can be predicted based on its expected genotype frequency using an assignment test. These data indicate that gene flow between local populations is restricted despite the long-distance seasonal movements undertaken by polar bears.     

Summer diet of beluga whales inferred by fatty acid analysis of the eastern Beaufort Sea food web

Beluga whales (Delphinapterus leucas) are the most abundant odontocetes in Arctic waters and are thus thought to influence food web structure and function. The diet of the Beaufort Sea beluga population is not well known, partly due to the inherent difficulty of observing feeding behaviour in Arctic marine cetaceans. To determine which prey items are critical to the Beaufort Sea beluga diet we first examine and describe the Mackenzie Delta and Beaufort Sea food web using fatty acid analyses. Fatty acid profiles effectively partitioned prey items into groups associated with their habitat and feeding ecology. Next, the relative contribution of various prey items to beluga diet was investigated using fatty acids. Finally, beluga diet variability was examined as a function of body size, a known correlate of habitat use. Beluga appeared to feed predominantly on Arctic cod (Boreogadus saida) collected from near shore and offshore regions. Size related dietary differences suggested larger sized beluga preferred offshore Arctic cod given the shared high levels of long chain monounsaturates, whereas smaller sized beluga appeared to feed on prey in their near shore habitats that included near shore Arctic cod. The presence of Arctic cod groups in shallow near shore and deep offshore habitats may facilitate the behavioural segregation of beluga habitat use as it relates to their size and resource requirements. Given Arctic cod are a sea ice associated fish combined with the accelerated sea ice loss in this region, beluga whales may need to adapt to new dietary regimes.     

Assessing the abundance of Bristol Bay belugas with genetic mark‐recapture methods

The Bristol Bay stock of beluga whales (Delphinapterus leucas) is genetically distinct and resides in Bristol Bay year‐round. We estimated the abundance of this population using genetic mark‐recapture, whereby genetic markers from skin biopsies, collected between 2002 and 2011, were used to identify individuals. We identified 516 individual belugas in two inner bays, 468 from Kvichak Bay and 48 from Nushagak Bay, and recaptured 75 belugas in separate years. Using a POPAN Jolly‐Seber model, abundance was estimated at 1,928 belugas (95% CI = 1,611–2,337), not including calves, which were not sampled. Most belugas were sampled in Kvichak Bay at a time when belugas are also known to occur in Nushagak Bay. The pattern of genetic recaptures and data from belugas with satellite transmitters suggested that belugas in the two bays regularly mix. Hence, the estimate of abundance likely applies to all belugas within Bristol Bay. Simulations suggested that POPAN estimates of abundance are robust to most forms of emigration, but that emigration causes negative bias in both capture and survival probabilities. Because it is likely that some belugas do not enter the sampling area during sampling, our estimate of abundance is best considered a minimum population size.     

Panmictic population structure in the hooded seal (Cystophora cristata)

Two putative populations of hooded seals (Cystophora cristata) occur in the North Atlantic. The Greenland Sea population pup and breed on the pack ice near Jan Mayen ('West Ice') while the Northwest Atlantic population is thought to pup in the Davis Strait, in the Gulf of St. Lawrence (the 'Gulf'), and off southern Labrador or northeast Newfoundland (the 'Front'). We used microsatellite profiling of 300 individuals at 13 loci and mitochondrial DNA sequencing of the control region of 123 individuals to test for genetic differentiation between these four breeding herds. We found no significant genetic differences between breeding areas, nor evidence for cryptic nor higher level genetic structure in this species. The Greenland Sea breeding herd was genetically most distant from the Northwest Atlantic breeding areas; however, the differences were statistically nonsignificant. Our data therefore suggest that the world's hooded seals comprise a single panmictic genetic population.     

Intraspecific variability in the “vowel”‐like sounds of beluga whales (Delphinapterus leucas): Intra‐ and interpopulation comparisons

The beluga whale (Delphinapterus leucas) has a rich and complicated vocal repertoire. However, different populations use similar and common types of signals. We studied physical features of one of these types, “vowels,” in three Russian populations: the White Sea population (European North), the Chukotka population (the Bering Sea, Chukotka), and the Okhotsk Sea population (Russian Far East) as well as in four summer aggregations of the White Sea belugas over several years in duration. The pulse repetition rate (PRR) at half of the duration of the signal was measured. We found that the PRR of “vowels” collected in the same summer aggregation during different years is stable in time but varies between locations. The degree of variation corresponds with the geographic distance between different locations. Significant differences were discovered between populations separated by thousands of kilometers, and to a lesser extent, between summer aggregations inhabiting different bays of the White Sea. The variation in PRR between the locations can be caused by the divergence of signals owing to the accumulation of random errors during transmission of these signals from generation to generation, which progressed independently in different summer aggregations and populations.     

Hematology and serum chemistry values in the beluga (Delphinapterus leucas)

Normal values and ranges for 31 clinical hematology and serum chemistry tests are reported for the beluga or white whale (Delphinapterus leucas). The values were collected over a 6-yr period from eight belugas maintained for display at Sea World (San Diego, California, USA) facilities and represent long-term evaluations for each animal in a controlled environment. They represent the first report for a number of serum chemistry values for the beluga. Normal values such as these provide an important data base from which to detect diagnostically important changes in health status for belugas in a zoological setting. They also establish a baseline from which to evaluate differences in normal values in free-ranging belugas and from which to diagnose disease problems in wild populations.     

POPULATION INDEX ESTIMATES FOR THE ST. LAWRENCE BELUGAS, 1973–1995

The belugas ( Delpbinapterus leucas ) inhabiting the St. Lawrence estuary were freely hunted until 1979 and thereby much reduced in numbers. In the late 1970s concerns for their future, and for the effect of contaminants, habitat modification, and marine traffic, became acute, and they were declared "Endangered" in 1983.
The numerical progress of the population has been monitored since 1973 by continued sample surveys. Nine of 10 surveys were aerial, and of those 9, 7 used high-altitude aerial photography. Published indices of population size showed an average increase of about 17 belugas per year (SE = 4) and a 1995 smoothed value of about 650 (SE = 40) belugas. This apparent temporal trend might be open to question owing to slight variation in survey methods and in the area covered by different surveys, so factors were developed to correct for differences in survey coverage and in treatment of data and for the effect of hunting on population trajectory. Revised assessments of the population trajectory and of the present smoothed index value changed little as a result. A series of total-count surveys in 1987-1992 estimated a decline of about 6 belugas per year, and when they were combined with sample surveys flown in 1984-1995, an increase of 10 belugas per year was estimated and a 1995 standardized index of about 640 (SE = 43). It appeared that hunting in the 1970s could have consumed all of, perhaps even slightly more than, the potential for natural increase of this population.     

Key Factors Determining the Distribution of Beluga Whales, <Emphasis Type="Italic">Delphinapterus leucas</Emphasis> (Pallas, 1776), in River Estuaries of Western Kamchatka

The results of observations on the distribution of beluga whales, Delphinapterus leucas (Pallas, 1776), in three large rivers of western Kamchatka in the summer and autumn seasons are discussed. In the summer, the number of beluga whales in the Khairyuzova, Belogolovaya, and Moroshechnaya rivers reaches 111–250 individuals. Most of the belugas enter the rivers during the flood tidal phase: the number of animals in the estuaries increases along with the rising water level to the maximum value at spring tide. The belugas do not move upstream out of the estuaries and tend to remain in the zone of mixing riverine and marine waters, where 20 species of fish and three species of invertebrates have been identified. At ebb tide, the belugas leave for the sea; however, during a large run of salmon some individuals remain in the estuaries and continue hunting in deep-water areas. The main issue that causes beluga whales to form summer aggregations in Kamchatkan rivers is the hunt for salmon. The distribution of beluga whales in river estuaries is defined by the dynamics and intensity of salmon spawning runs. The preference of beluga whales for these rivers can be explained by the channel type of their estuaries.     

Beluga (Delphinapterus leucas) habitat selection in the eastern Beaufort Sea in spring, 1975�C1979

An understanding of the adaptability of belugas (Delphinapterus leucas) to changing ice conditions is required to interpret and predict possible changes in habitat selection in response to projected loss of sea ice throughout the circumpolar Arctic. We analyzed beluga observations made during spring aerial surveys for ringed seals conducted from 1975 to 1979 in the eastern Beaufort Sea. Despite inter-annual variability in the extent and distribution of sea ice, belugas consistently selected areas with water depths of 200–500 m and heavy ice concentrations (8/10 to 10/10) while areas of open water to light ice concentrations (0/10 to 1/10) were not selected. Belugas were also found in proximity to regions with ≥0.5 degrees seafloor slope which include the continental slope and other areas with the potential for oceanographic upwellings. In most years (4 of 5), fast-ice edges and coastal areas were not selected. In the lightest ice year analyzed, belugas showed less specificity in habitat selection as their distribution expanded and shifted shoreward to fast-ice edges. The observed distribution is discussed in terms of predator–prey relationships particularly with reference to beluga feeding on polar cod (Boreogadus saida). More research is required to examine and compare possible changes in distribution since the late 1970s and to investigate the factors driving the patterns described.     

Population-specific home ranges and migration timing of Pacific Arctic beluga whales (Delphinapterus leucas)

Two populations of beluga whales (Delphinapterus leucas), the Eastern Beaufort Sea (BS) and Eastern Chukchi Sea (ECS), make extensive seasonal migrations into the Pacific Arctic. However, the extent to which these populations overlap in time and space is not known. We quantified distribution and migration patterns for BS and ECS belugas using daily locations from whales tracked with satellite-linked transmitters. Home ranges and core areas in summer (July and August) and in each month (July–November), daily displacement, dispersal from core areas, and autumn migration timing were estimated. Distinct summer and fall distribution patterns and staggered autumn migration timing were identified for BS and ECS whales. Summer home ranges for each population had less than 10 % overlap. Monthly home ranges were also relatively distinct between populations except in September (up to 88 % home range overlap). A distinct east–west shift in focal area use occurred in September that persisted into October, with the two populations essentially switching longitudinal positions. Highest daily displacements occurred during the migratory period in September for BS whales and October for ECS whales, further indicating westward fall migration was offset between populations. Sexual segregation of males and females within a population also varied monthly. Autumn migration timing as well as differences in spatial and temporal segregation between BS and ECS beluga populations may be a result of maternally driven philopatry and population-specific adaptations to dynamically available resources. Our results contribute to the management of these populations by identifying seasonal area use and differences in migration patterns.