Determining the species assemblage and habitat use of cetaceans in the Svalbard Archipelago, based on observations from 2002 to 2014
Published: June 1, 2018
Latest article update: Aug. 7, 2023
This study used 13 years of cetacean sighting data (2002–2014) from waters around the Svalbard Archipelago to determine key habitats for year-round resident species as well as seasonally resident species, and to explore spatial overlap between these groups via a combination of kernel density estimation and Maxent modelling. The data set consists of observations made by research vessels conducting various marine studies, coast guard ships and marine-cruise tourist operators. Data are reported from the seasonal period in which there is daylight (March-November), though 95% of the observations occurred June- September. Changes over the study period were investigated, within the limits of the data, to explore whether range shifts may be occurring. Fifteen cetacean species were reported. Among the resident ice-associated cetaceans, only white whales were reported frequently; they were seen exclusively in coastal habitats, in accordance with their known use of tidal glacier fronts for feeding in this region. Narwhal and bowhead whales were rare. Seasonally resident minke whales, fin whales, humpback whales, blue whales and sperm whales as well as small dolphins were seen frequently, in broad and somewhat overlapping habitats. Other less common seasonal residents included killer whales, northern bottlenose whales and sei whales; harbour porpoises and long-finned pilot whales were also reported, but rarely. Shifts over the study period towards higher latitudes, and into coastal environments, were observed for several seasonally resident species. These expansions are likely linked to warming ocean temperatures and a precipitous decline in sea-ice cover in the area.
Keywords
Climate change, citizen science, Arctic, whales, marine mammals, Maxent modelling
Understanding both the geographic range and the habitat use of a species is fundamental to characterizing its ecology and enabling predictions regarding how it will react to environmental changes (Gaston & Fuller 2009). Large, highly mobile marine predators that do not easily carry tracking devices present a challenge in this regard, as is the case for all cetaceans. However, understanding the distribution of these animals is particularly important because they exert strong top-down control on their prey (Ainley et al. 2007; Bowen 1997; Coyle et al. 2007; Surma & Pitcher 2015) and are a major source of interference competition with other predators (Fraser et al. 1992; Ainley & Hyrenbach 2010), and thus have far-reaching effects on the ecosystems they occupy. It has recently become clear that cetaceans are also important in stimulating productivity, particularly in iron- limited regions (e.g., Lavery et al. 2014). The dramatic changes in distribution and abundance of the “Great Whales” due to extensive hunting in previous centuries have had significant impacts on the environments they occupy (Hacquebord 1999; Rugh et al. 2003). Many cetacean species are now recovering, but this is occurring during a time of rapid environmental change, particularly in polar regions, so it is likely that new distributional patterns will become established (Smetacek & Nicol 2005; MacLeod 2009; Scheffers et al. 2016; Thomas et al. 2016).
The Svalbard Archipelago, located in the High Arctic within the North-east Atlantic (Fig. 1), has an ecological and cultural history that is intimately associated with cetaceans (e.g., Arlov 1994). The first large-scale commercial harvesting of whales, which commenced in the early 1600s, was based in and around this group of islands (Hacquebord 1999). Vast numbers of whales were taken in the centuries that followed, involving a wide array of species (Burnett 2012). The bowhead whale (Balaena mysticetus) population in this region was driven to the brink of extinction before becoming protected in 1939. Other whale species were also depleted, and iteratively protected. Only minke whales (Balaenoptera acutorostrata) are currently harvested in the region; regular surveys ensure that catches are sustainable (Bogstad et al. 2015). However, very little is known about current trends in abundance or distribution and habitat use of Arctic resident or other seasonally occurring cetacean species in this Arctic gateway, which is a known hotspot for cetaceans (Kovacs & Lydersen 2006). Updated knowledge regarding the cetacean community around the Svalbard Archipelago is of particular importance because the region is undergoing profound physical and biological changes associated with warming ocean temperatures. Air and water temperatures are increasing (see Descamps et al. 2017 for a review), while sea-ice coverage is concomitantly decreasing (Onarheim et al. 2014; Pavlova et al. 2014; Laidre et al. 2015). These changes to the physical environment are dramatically altering biota in the region, with boreal invertebrate and fish species replacing traditionally Arctic species (e.g., Fossheim et al. 2015; Dalpadado et al. 2016; Gluchowska et al. 2016). The changing ice conditions in particular are also resulting in substantial increases in ship traffic and industrial interest in the region, which represent potential risks to cetaceans (e.g., Reeves et al. 2014).
Despite the near extinction of bowhead whales due to overharvesting, all three ice-associated endemic Arctic whales - bowhead whales, white whales (Delphinapterus leucas) and narwhal (Monodon mono- ceros) - still occupy the waters around Svalbard (Kovacs et al. 2009). These year-round residents of the High Arctic are all well adapted to living in ice-covered seas, where the sea ice shelters their young from storms, provides protection from killer whales (Orcinus orca) and offers an environment with little competition for available prey (Kovacs et al. 2011). These species have evolved life-history traits that buffer the normal, high levels of inter-annual fluctuation in climate in the Arctic (Harington 2008). However, the current rapid declines of their sea-ice habitats are a serious concern (Kovacs et al. 2011; Laidre et al. 2015), and the first responses to these declines are likely to be changes in distribution (Tynan & DeMaster 1997; Gilg et al. 2012).
Svalbard waters also provide important summer foraging habitats for seasonally resident, migratory baleen whales, including the minke whale, fin whale (Balaenoptera physalus), humpback whale (Megaptera novaeangliae) and blue whale (B. musculus) (Kovacs & Lydersen 2006). White-beaked dolphins (Lagenorhyn- chus albirostris), killer whales, long-finned pilot whales (Globicephala melas), sperm whales (Physeter macro- cephalus) and northern bottlenose whales (Hyperoodon ampullatus) are also known to occur in this region (Kovacs et al. 2009; Skem-Mauritzen et al. 2011).
Monitoring cetacean distributions over relatively large spatial scales in the Arctic is logistically and financially challenging. However, data can be gathered on these highly visible animals in some areas by engaging “citizen science”: using volunteers under the guidance of professional scientists and naturalists to collect robust data in order to answer scientific questions (Silvertown 2009; Thiel et al. 2014). The use of citizen science to assess habitat use by cetaceans in the Svalbard area is feasible because of the attractiveness of these animals to the flourishing marine cruise tourist industry. This study used 13 years of cetacean sighting data (2002-2014) from waters around the Svalbard Archipelago (Fig. 1) to provide a baseline description of the assemblage of cetaceans in this region, with the aim of identifying key habitats for year- round resident Arctic species and seasonally resident species and to explore spatial overlap between these groups. Temporal changes over the period of study were also investigated, within the limits of the data, to explore whether range changes may be occurring.
Cetacean sightings data have been collected in the waters surrounding the Svalbard Archipelago in the Norwegian Arctic (74-84°N, 0-35°E; Fig. 1) since 2002. These data are stored in the NPI’s MMSDB (http://www.npolar.no/en/services/mms).
Contributors to this database include Norwegian Coast Guard vessels, research vessels operating in Svalbard waters from the Institute of Marine Research, the Norwegian Polar Institute and other organizations, the Governor of Svalbard’s field inspectors, and marine cruise tourist operators (hereafter called EVs, most notably vessels associated with the AECO. The latter group accounts for approximately half of the data reported.
Most data are reported to the MMSDB via a custom-designed data-inputting software package (available at http://www.npolar.no/en/services/mms). These records include species identification, the number of individuals (pod size), the latitude and longitude of the observation, date, ship name, contact person and any other comments that the observers deem worth reporting related to the animals’ behaviour, sex, age classes, etc. EVs all have at least two observers around the clock (24-hour light in summer), minimally two crew members stationed on the bridge and sometimes also several dedicated cetacean experts (as well as pinniped, polar bear and bird experts) actively looking for animals.
In July 2015, interviews were conducted with staff aboard some of the EVs that contribute regularly to the MMSDB, including: MV Sea Adventurer and MV Ocean Nova (Quark Expeditions); MS Quest (Polar Quest); MV Ortelius (Oceanwide Expeditions); MS Malmo and MS Origo (Expeditions Origo), as well as MS Havsel (a vessel used for tourists and film crews in Svalbard waters). Staff confirmed that all observers were able to identify all cetacean species seen in the Svalbard region confidently, and that they recorded observations as “unidentified whale” or “unidentified dolphin” if there was any uncertainty. Staff also confirmed that all cetacean observations are routinely recorded, as well as the number of animals seen; estimates are made for particularly large groups. NPI marine mammal scientists oversee all data entering the MMSDB and check unexpected or rare observations with persons reporting them; rare sightings are usually confirmed via photographic evidence.
The data set used herein includes all of the cetacean observations during the months with daylight in the Svalbard region (March through November) as research vessels have activity throughout most of the period with light at these high latitudes. EVs operate principally from June through September. Only the Norwegian Coast Guard operates year-round, and the few observations reported from the dark period (December-February) were not included in this study because of the risk of misidentification. EVs typically travel from Tromso in northern Norway to Longyearbyen in Svalbard, before cruising within the Svalbard Archipelago. The favoured route circumnavigates Spitsbergen, the largest island in the archipelago, although many boats also go along the continental slope to the west of Spitsbergen, and around the islands in the north-east when conditions permit. The Norwegian Coast Guard and research vessels have much greater spatial coverage, moving throughout most of the study area, including heavily ice-covered waters.
In this study, the terms “observation”, “sighting” and “record” are used to denote a single recorded observation, which can consist of an individual animal of one species, or a group of animals of one species seen together. Dolphins (of the genus Lagenorhynchus) were treated as a single group for analyses (hereafter referred to as “small dolphins”)
because of the large number of reports of undifferentiated dolphins; white-beaked dolphins are known to be the most common dolphin in the Svalbard region (Kovacs et al. 2009), but white-sided dolphins (I. acutus) also occur.
In order to place cetacean observations within the bathymetric and physical oceanographic context in which they occurred, a grid with cell size of 12.5 km2 (selected because this is the minimum resolution of the environmental data used) was created using the Create Fishnet tool within ArcMap (version 10.4.1), from coordinates 74°N, 0°E in the south-west to 84°N, 35°E in the northeast and selected environmental parameters were mapped onto this study area. The International Bathymetric Chart of the Arctic Ocean version 3.0 (Jakobsson et al. 2012) was used for the bathymetry of the study region. Areas where glaciers had retreated since the production of this chart were assigned a depth of 10 m. To ensure that the bathymetric and other oceanographic parameters did not include values from land, land areas within the grid cells were eliminated by combining the Svalbard land polygon and the grid polygon, and deleting the land layer. Monthly composites of the MODIS SST 11 (i daytime, at 9 km resolution, were used for the SST of the region for each month (March-November inclusive) within each year of the study period (NASA GSFC, Ocean Biology Processing Group 2014). Mid-monthly five-day median sea-ice concentration data sets from the ASI algorithm SSMI/S at 12.5 km resolution were used for sea-ice data for each month (March- November inclusive) within each year of the study period (see Kaleschke et al. 2001; Spreen et al. 2008; Kern et al. 2010).
ArcMap software was used to assign environmental parameters to cetacean sighting locations. Mean values for water depth, slope (gradient of the slope in degrees from the horizontal plane), aspect (horizontal direction of slope face in degrees from north) and distance from coast (distance from the nearest point of land in the Svalbard Archipelago) were calculated for each grid cell. The Spatial Analyst Tool was used to create slope and aspect rasters from the bathymetric map. The mean depth, slope and aspect were then calculated for each grid cell using the Zonal Statistics as Table tool. The distance to the shore was calculated by converting each grid cell into a point feature at the centre of each respective cell using the Feature to Point tool. The distance from the Svalbard land polygon was then calculated for each of these points, and the points were rejoined to the grid. Mean values for sea-ice concentration and SST were calculated for each grid cell for every month in which observations occurred. Because of variation in the spatial satellite coverage, areas with missing values for these two parameters were assigned the value of the grid cell to the north of them. Each cetacean sighting was assigned all of the physical parameters of the grid cell in which it fell; these data were then edited so that each observation had the appropriate SST and sea-ice data for the month in which it occurred. A few months had poor satellite coverage of SST in the Svalbard region, so interpolation was not possible; therefore, some of the records (n = 32) did not have SST values assigned to them, and these observations were excluded from the modelling.
KDEs were calculated using the Kernel Density Analysis tool within ArcMap, with a search radius of 50 km and a cell size of 1 km, for each species with over 100 observations to display where each of the species were observed in this study (see Silverman 1986; Powell 2000). The 50 km search radius was used as a general daily displacement distance for cetaceans using foraging grounds (e.g., Hobbs et al. 2005; Zerbini et al. 2006; Dalia Rosa et al. 2008; Heide-Jorgensen et al. 2012; Kennedy et al. 2014; Heide-Jorgensen et al. 2015). The resulting KDEs were displayed using 10 quantiles. Images were resampled during display using bilinear interpolation to facilitate visualization of contours of density. KDEs were also produced in two temporal categories: “early” (2002-08, seven years) and “late” (2009- 2014, six years), with the two periods having approximately equal numbers of observations, (n = 3005 and n = 2967 respectively) for selected seasonally resident species (minke whale, fin whale, blue whale and humpback whale) to explore potential changes in habitat use. These KDEs were subtracted from one another for each species, to visualize change in the density of sightings of individual species between the two periods. It is important to note that the Kernel Density Analysis tool within ArcMap does not distinguish between land and water, so it smoothes the aquatic distributions of cetaceans across land points. This is not ideal in grid cells with a lot of land and little water (such as in narrow straits).
Modelling of cetacean habitat suitability was conducted using Maxent software (version 3.3.3k; see Phillips et al. 2006). This software is appropriate for the MMSDB data because, unlike generalized linear models and generalized additive models, Maxent does not require animal absence to be recorded (Kanaji & Okazaki 2014). The set of observations, with all environmental parameters attached, including the correct SST and sea-ice concentration data for the appropriate month of each individual observation, was converted into a CSV file, and this was added into “Samples” in the Maxent software. Models were displayed relative to a grid of background data from August 2014. This was a CSV file composed of the 12.5 km2 resolution grid, with
each cell containing the mean values for each of the bathymetric parameters over that area, as well as the mean SST and ice concentration from August 2014, selected because August is the month of greatest observation frequency (Supplementary Fig. Sla) and 2014 was the latest year of the study period, and presumed to be the most representative of current conditions.
Models were run using 70% of the observations for model training, and the remaining 30% for model testing. All other model settings were left at standard values, which assumes species to be evenly distributed across the entire environment (see Merow et al. 2013). The modelled habitat suitability CSV files were converted into point features, and then into rasters using the Feature to Raster tool in ArcMap, with an output cell size of 12.5 km. All modelled outputs were resampled during display using bilinear interpolation.
AUC and percentage contribution of each variable were calculated for all species models. AUC was used to measure model performance. With presence-only data the maximum achievable AUC is less than 1 and a random model prediction gives an AUC of 0.5 (Phillips et al. 2006). The percentage contribution is calculated by randomly altering the variable values amongst the training points and measuring the subsequent decrease in training AUC, with large decreases indicating a heavy dependence by the model on that variable (Phillips 2012).
Comparisons of modelled habitats were undertaken using ENM Tools software (version 1.4.4; see Warren et al. 2010). Niche similarity was measured using Schoener’s D and a measure I, derived from Hellinger distance. Both measures range from 0 to 1, with 0 denoting no overlap in the modelled grids of environmental suitability and 1 denoting 100% overlap, i.e., with all grid cells being equally suitable for both species (Warren et al. 2010). Schoener’s D treats the modelled suitability scores as proportional to species abundance, whereas I treats them as probability distributions (Warren et al. 2010). The I statistic was used to classify degrees of overlap, with I > 0.9 denoting a high level of overlap and I > 0.8 denoting a moderate level of overlap.
A few caveats must be considered when interpreting this data set. Vessel activity is greatest around the coastal regions, particularly on the west coast of Spitsbergen. Few vessels venture into areas with heavy ice cover, although EV visits to the ice edge are quite common. The opportunistic nature of the observations means there is no way to account for multiple sightings of the same individual, although tourist ships have a system whereby they avoid being in sight of each other as much as is possible, which minimizes resighting of the same individual on short- time frames. Ship traffic is greatest from June through September so effort is concentrated in these months (Supplementary Fig. Sla). Additionally, the increase in EV activity through the study period (from approximately 50 000 cruise tourists in 2002 to 60 000 in 2014; MOSJ 2017) may have increased effort towards the later years of this study. Finally, the tourist visitation season has grown longer, and tourism has expanded northward, as sea-ice cover has declined (markedly from 2006 onward).
The MMSDB included 6676 cetacean sighting records in the period 2002-2014; 613 (9%) of these were not identified to species (or group in the case of dolphins), and 91 (1%) were recorded in the dark winter months (December-February) and hence are not included in the analyses herein. MMSDB cetacean records included 15 species (Table 1) during the study period. Minke whales were the most frequently observed whale species over the study period (Table 1, Fig 2a); however, the total number of individual minke whales was low compared to some other species because of the solitary behavioural tendencies in this species (Table 1, Fig. 2b). Fin whales, small dolphins and humpback whales were the next most frequently observed taxa (Table 1, Fig. 2a), with small dolphins accounting for the greatest number of animals amongst all of the cetaceans (Table 1, Fig. 2b). White whales were observed frequently, and accounted for the second greatest total number of animals because of their large pod sizes (Table 1, Fig. 2). The two other resident, ice- associated whales, the bowhead whale and the narwhal, were observed during most years of the study, though quite infrequently, with 30 and 18 observations, respectively (Table 1). Blue whales and sperm whales each had over 100 observations in the database. Several of the other seasonally resident species - killer whales, northern bottlenose whales and sei whales (B. borealis) - were observed relatively infrequently, with <41 observations. Two species - harbour porpoises (Phocoena phocoena) and long-finned pilot whales - were rare, with <19 observations (Table 1).
Maxent habitat models produced high training and test AUC values (>0.86) for all species that had numerous sightings (>100) (Supplementary Table SI). High scores were also achieved in models inferring habitat suitability for killer whales, northern bottlenose whales and sei whales. Model outputs are presented for all of these species. Distance to land and ice cover were important parameters for many species, while slope, SST, depth and aspect were critical to only a few (see Supplementary Table SI for details). The remaining species were deemed unsuitable for modelling because of either few observations (long-finned pilot whales and harbour porpoises) or lack of coverage of core habitat (bowhead whales and narwhals).
White whales were the most frequently observed resident cetacean, with 286 observations and pod sizes of over 50 individuals seen regularly, although average pod size was 10 (Table 1, Supplementary Fig. S2). This species was found close to the coast (Fig. 3a) (median distance = 1.47 km; Fig. 4d), in shallow water (median = 57 m; Fig. 4c). White whales were observed frequently in areas with considerable amounts of sea-ice coverage, with a quarter of the recorded observations coming from areas with at least 30% ice cover (Fig. 4a). However, half of the observations of this species were made in localities with less than 2% ice cover. White whales were observed in areas with somewhat lower SSTs than the seasonally resident species, with ca. 25% of the records being made in waters where SST was <2°C (Fig. 4b). Although white whales occurred in most fjords in Svalbard, they were most commonly seen along the west coast of Spitsbergen (Fig. 3a, Supplementary Fig. S2). Habitat modelling suggested that white whales occupy the narrowest range of habitats of all the species investigated (Fig. 3b); their habitat was concentrated in fjords and straits, very close to the coast, with habitat suitability dropping to <5% at a distance of 20 km from land (Fig. 3b).
Bowhead whales were reported infrequently (Fig. 5). Pod sizes ranged from one to seven animals in the 30 observations in the MMSDB (Table 1). Most of the observations of this species occurred along the slope of the Yermak Plateau to the north-west of Spitsbergen, in an area known as Whalers Bay, although some sightings did occur in the north-east parts of the archipelago and along the west coast.
Similar to bowhead whales, the other resident Arctic species, the narwhal, was observed infrequently (Table 1, Fig. 5). Typically, observations of narwhals were made at high latitudes, in areas that had relatively dense ice cover, such as in W ahlenbergfj or den on the west side of Nordaustlandet, or associated with pack ice north of 80°N. Narwhals were observed both in shallow coastal environments and in deep waters beyond the continental slope. Narwhals were found at higher latitudes than other species (Fig. 5).
Minke whales were reported throughout much of the study area (Fig. 6a), to as far north as 81.6°N (Supplementary Fig. S3). The greatest densities of sightings occurred around Bjsrnsya, along the continental shelf to the west and north of Spitsbergen, and at the 100 m slope edge east ofEdgeoya. This species occupied areas with depths from <20 m to >4000 m, however they were typically found in areas with intermediate depths (median depth = 188 m; Fig. 4c) over the continental slope and therefore fairly close to land (median distance = 38.2 km; Fig. 4d). Most observations came from open-water areas, although a few sightings came from areas with some drift ice (Fig. 4a). Habitat modelling suggested that minke whales occupy a broad array of open-water and coastal habitats (Fig. 6b), all the way up to the ice edge. The habitat model for this species suggested that areas over deep water beyond the continental slope and areas with high ice cover were unsuitable habitats (<20%; Fig. 6b).
Similar to minke whales, fin whales were observed throughout most of the study area (Fig. 7a), as far north as 81.5°N (Supplementary Fig. S4). Greatest densities of sightings for this species came from around Bjornoya and along the continental slope to the west and north of Spitsbergen (Fig. 7a). Fin whales were observed over a broad array of depths (median depth = 250 m; Fig. 7c); however, they tended to be found in slightly deeper areas than minke whales (Fig. 7a). Fin whales occupied areas with open water; however, a few observations came from areas with loose drift ice (Fig. 4a). Fin whales were less prevalent to the east of Edgeoya compared to minke whales, but were more common in Storfjorden in south-eastern Svalbard. Many fin whale observations were made close to the Knipovich Ridge west of Spitsbergen. Modelling suggested that fin whales occupy a broad array of habitats (Fig. 7b), with the most suitable habitats occurring along the steeper areas of the continental slope around the 500 m depth contour (Fig. 7b). Areas with particularly shallow and flat bathymetry were modelled to be unsuitable habitats (<20%), including Spitsbergenbanken and around Edgeoya. Additionally, deep areas beyond the slope, shallow areas of the Barents Sea and areas with dense ice cover were also found to be unsuitable habitats (<10%).
Humpback whales had a patchy distribution, which extended as far north as 81.1°N (Fig. 8a, Supplementary Fig. S5). This species was observed frequently along the continental slope and in shallow areas west of Spitsbergen. However, contrary to the distribution of the other rorquals in this study, humpback whales were also found in large numbers to the south of Nordaustlandet. This species tended to occupy shallower depths, within a narrower range than the other rorquals, with 75% of the records coming from areas with depths of <263 m (Fig. 4c, Supplementary Fig. S5). This species was not observed in areas with sea ice (Fig. 4a). The Maxent model demonstrated that humpback whale habitat occurred along the continental shelf and slope, but not beyond the slope into deep water (Fig. 8b). The major fjord systems around Spitsbergen were fairly suitable habitat (>30%), as were some areas southeast of Spitsbergen (suitability >40%).
Blue whales were observed west and north of Spitsbergen (Fig. 9a). They were most frequently seen in coastal habitats (median distance = 11.5 km; Fig. 4d) and occasionally along the continental slope (Fig. 9a). The northernmost record was 81.5°N (Supplementary Fig. S6). Blue whales were not observed in the east of the Svalbard Archipelago, nor were they seen frequently beyond the continental slope to the west (Fig. S6). They occurred in waters of moderate depths (median depth = 214 m; Fig. 4c), sometimes in close proximity to loose sea ice; a quarter of the observations were from areas with at least 16% ice cover (Fig. 4a). The most suitable habitats identified from the Maxent model were steep coastal areas (Fig. 9b), reflecting the contribution of the distance to land and slope parameters to the model (Supplementary Table SI). Hinlopenstretet was particularly suitable (>75%), as was Isfjorden (40-70% suitability) (Fig. 9b). The continental slope along the 500 m depth isobath was also a fairly suitable habitat (>40%). However, suitability dropped rapidly to less than 10% at greater depths, and was less than 1% beyond the slope. Habitat tended to be unsuitable for blue whales to the east of Spitsbergen (Fig. 9b).
Small dolphins were observed across much of the study area (Fig. 10a), as far north as 80.8°N (Supplementary Fig. S7). They were encountered most frequently along the shelf-break to the south and west of Svalbard, and in deep-water areas (median depth = 355 m; Fig. 4c) affiliated with the Knipovich Ridge (Fig. 10a). Small dolphin observations occurred most commonly in open water, far from land (median distance = 80.2 km; Fig. 4d), with only 37 out of 1031 observations being within 10 km of the coastline. Small dolphins tended to be observed in waters with slightly higher temperatures than most other seasonally resident species (median = 5.7° C; Fig. 4b), and were not seen in areas with sea ice (Fig. 4a). Modelling suggested that the continental slope from a depth of 500-1000 m (Fig. 10b) was the best habitat for small dolphins. Areas with sea ice cover were modelled to be unsuitable (<10%).
Sperm whales were observed most frequently in the south-west of the archipelago, west of Bjornoya (Fig. Ila), primarily in deep water (median depth = 1594 m; Fig. 4c) along and beyond the continental slope, far from land (median distance = 119.2 km; Fig. 4d). Only 34 out of 130 observations came from north of 76°N; however, several observations were made close to 80°N (Supplementary Fig. S8). Sperm whales were associated with the highest SSTs (median temperature = 7.5°C; Fig. 4(b). Suitable habitats occurred from the continental slope out into Fram Strait (Fig. 11b), primarily along the slope between the 500 m and 2000 m isobaths (suitability 60-85%; Fig. 11b).
Most areas shallower than the 500 m isobath were modelled to be unsuitable (<25%) for sperm whales, apart from some areas around Kong Karls Land and off north-west Spitsbergen. Areas with sea ice were modelled to be unsuitable habitats (<5%; Fig. 1 lb).
Some few seasonally resident species had few observations, so their habitat models should be viewed with caution, but the locations of their occurrences were deemed noteworthy particularly in the context of expected expansion into the Svalbard area with climate change induced alterations of habitat and prey (habitat models for species with 31-41 observations are shown in Supplementary Fig. S9). Killer whales were observed over a broad area and across a broad range of habitats north to 80°N (Fig. 5). Sei whales were the least frequently observed seasonally resident baleen whale. They were seen as far north as the northern tip of Spitsbergen (79.87° N, 14.88° E; Fig. 5). Northern bottlenose whales were primarily seen in deep offshore waters south-east of Spitsbergen, close to the Knipovich Ridge (Fig. 5). Only eight of 37 observations took place in waters of <1800 m depth, all of these just west of the entrance of IsHorden. Harbour porpoises and long-finned pilot whales were both observed only infrequently around Svalbard, in relatively shallow areas mostly in the southern parts of the study area (Fig. 5).
Niche similarity analyses suggested that no two species were suited to identical niches (see Table 2). Habitat models suggest that white whales occupied a unique niche. However, some species did show overlap. Minke whales, fin whales, humpback whales and small dolphins had high or moderate degrees of habitat overlap (Table 2). The habitat of blue whales overlapped with minke whales, fin whales and humpback whales.
Killer whale habitat models need to be viewed with caution because of low numbers of sightings of this species in the MMSDB. However, the areas suggested to contain suitable habitat for them overlapped the ranges of many other cetaceans including small dolphins, minke whales, fin whales, humpback whales and, to a degree, sperm whales (see Fig. 5 Supplementary S9a and species-specific figures) as well as sei whales. Regions with heavy ice cover, where bowhead whale and narwhal were observed, abutted the extremes of the range in which killer whales were observed.
Minke whales, fin whales, blue whales, humpback whales and sperm whales were all observed more frequently at higher latitudes towards the end of the study period, whereas other species such as the small dolphins did not exhibit such a latitudinal change (Fig. 12). Blue whales exhibited the steepest gradient for this trend. Five sightings of blue whales have been made in areas above 81°N since 2012; all three of these observations were groups of three individuals (Supplementary Fig. 86). Blue whales have also contributed to a greater proportion of the observations in recent years (Fig. 2a); 29 observations were made in the early period of this study, whereas 172 observations of this species were made in the late period (see Fig. 2a).
Blue whales, fin whales, minke whales and humpback whales have all been observed more frequently in near-shore areas within the archipelago in late years compared to early years. All four species were also sighted at more northerly latitudes in recent years (Fig. 13). Similarly, white whales exhibited an increase in the number of observations and total animals throughout the study period (Fig. 2). Only 14 recorded observations of this species were made from 2002 to 2005; 91 observations were made between 2006 and 2009; and 181 observations were made between 2010 and 2014. White whales were also unusual in being reported most frequently in June and July (the months with fourth and third highest total reporting, respectively), whereas all other species with >100 observations were observed most frequently in August, when the monthly peak in sightings occurred (Supplementary Fig. 810). The long-term patterns were seen despite no systematic increase in the numbers of sightings reported to the database through the study period. The number of sightings was on average 466 ± 262 annually, but varied from 100 to 900 over the study period with no specific time trend (see Supplementary Fig. Sib).
The logistical and financial challenges of monitoring cetaceans over large temporal and spatial scales make databases such as the MMSDB, which use opportunistic sightings, a cost-effective way of gathering data on cetacean habitat use. Previous cetacean studies that utilize citizen science have provided useful insight into the abundance and distribution of cetaceans elsewhere (e.g., Tonachella et al. 2012; Cheney et al. 2013; Beck et al. 2014; Bruce et al. 2014). This project aimed to provide a baseline description of the assemblage and habitat use by cetaceans in the Svalbard region, and provide insight into any evidence of changes that may be occurring. However, it is important to note that the sightings in this study were not part of systematically designed surveys, and we have no accurate measures of sighting effort, though the routes covered by the EVs have remained very similar in most years in the study period. Some considerations must be taken into account when interpreting the results because biases do exist within the data. Most of the Spitsbergen (the largest island in the archipelago) coast is heavily trafficked by EVs, which consistently travel to locations including specific glacier fronts, walrus haul-out sites and bird colonies; they also go to areas where they have heard of recent animal sightings (especially bears and whales). Years of very dense pack ice to the north and east naturally restrict visitation to these areas, whilst overall observer effort has expanded into newly openwater areas to the north where the sea ice has receded. Boats tend to spend little time in Storfjorden; exemplified by the few observations of white whales there, despite this area being known from tracking studies to be an area heavily used by this species (e.g., Lydersen et al. 2001). There is also relatively little effort in the deep waters beyond the shelf break, such as those associated with the Knipovich Ridge; species occupying these areas will be under-represented. Visual shipbased detection of cetaceans also includes inherent bias. Different species vary in their detectability, which can further be affected by sea state (e.g., Barlow 2015), and in their behaviour in relation to boat activity (e.g., Richardson & Wursig 1997).
The NPI's MMSDB records confirm that the waters around the Svalbard Archipelago comprise important summer foraging habitat for a wide variety of seasonally resident cetacean species, as well as being core habitat for white whales. Bowhead whales and narwhals are seen in the area, though infrequently in open-water areas. This study also clearly demonstrates that the IUCN and other organizations presenting overviews of the distribution of cetacean species need to update northerly range boundaries for several cetaceans (Kovacs & Lydersen 2008) within the North Atlantic Arctic including sperm whales (Taylor et al. 2008a), northern bottlenose whales (Taylor et al. 2008b) and sei whales (Reilly et al. 2008a).
White whales have received more research attention than any other cetacean species in Svalbard waters. Distribution has been studied using satellite-tracking studies (Lydersen et al. 2001; Lydersen et al. 2002; Lydersen et al. 2014), Marginal Ice Zone aerial surveying (Vacquie-Garcia et al. 2017) and passive acoustic monitoring (Castellote et al. 2013). Additionally, dietary studies and toxicological studies have been undertaken (e.g., Dahl et al. 2000; Villanger et al. 2011). The KDEs in the current study confirm the tightly coastal distribution of this species in the Svalbard area, consistent with earlier satellite tracking records, as well as the preferential use of Hords that terminate in tidal glaciers. However, the habitat modelling undertaken herein for this species did not identify areas in the northern parts of Storfjorden as good habitats, although this area is known to be heavily utilized (Lydersen et al. 2001), and it also overemphasized the attractiveness of Hirdopenstretet as habitat compared to what both sightings and tracking suggest. These differences likely arise because of the extraordinarily coastal distribution of this white whale population - most locations from satellite tags are within tens of metres from land, which makes the estimation of ice cover and other variables challenging at small spatial scales (environmental grid resolution of 12.5 km). This, in combination with the fact that the actual locations of glacier fronts are changing rapidly and the fact that depth, and SST values are few and crude this close to shore undoubtedly contribute to poor model performance for this species.
White whale observations occurred mostly in June and July, whereas background data for all species were mapped for August 2014, which might also influence the modelling results for this species. This study found that white whales were frequently observed in areas with moderate ice cover, even after the collapse of annually-formed sea ice in western Svalbard in 2006, which is related to their use of tidal glacier fronts, where pieces of glacier ice are often found in significant concentrations. These areas provide refugia for Arctic fish and invertebrate species in Svalbard because of the presence of ice and the cold temperatures of glacial runoff waters. They also have high prey biomass for mammals and birds because circulation patterns induced by the freshwater outflows advect large amounts of production towards the glaciers (Lydersen et al. 2014). The preferred prey of white whales, polar cod (Boreogadus saida) (Dahl et al. 2000), are increasingly restricted to areas in front of glaciers because of warm Atlantic water inflows into Hords on the west and north coasts of Spitsbergen (see Hop & Gjosaeter 2013; Dalpadado et al. 2016; for details). The tightly coastal distribution of Svalbard white whales, without any offshore seasonal migration, is in contrast to white whales elsewhere in their range, which routinely exhibit clear offshore seasonal migration patterns (e.g., Goering & McRoy 1974; Davis et al. 1980; Heide- Jorgensen et al. 2010; Asselin et al. 2012).
The global white whale population was recently listed as Near Threatened on the IUCN Red List (Jefferson et al. 2012), largely on the basis of expected negative impacts of climate change on some populations. Actual trends are known for only a quarter of the recognized populations (five out of 18, three of which are declining [Laidre et al. 2015]). Nothing is known about abundance or trends for the Svalbard population. This species does have an extensive range, south to the St Lawrence Estuary in North America (Lesage & Kingsley 1998), and its flexible dietary habits should serve to reduce the impacts of climate-induced habitat changes (see Laidre et al. 2008). However, their extensive use of the lipid-rich polar cod as prey, and the observed declines in this Arctic fish species in the Barents region (Hop & Gjosaeter 2013), could pose a problem for white whales in Svalbard if they do not readily prey shift. Moreover, even if they do shift to more open-water temperate prey, they are likely to face increased competition from other cetacean species. Increased sighting frequency through this study could be a positive sign of increased abundance, but unstructured observational data alone are insufficient to conclude numerical trends. This species should be assessed, and monitored in the future in Svalbard.
Bowhead whales have been reported to the MMSDB only infrequently, which is not surprising given that the Spitsbergen bowhead whale population is Critically Endangered (Reilly et al. 2012). However, it should be noted that a recent survey in the Marginal Ice Zone north of Svalbard using both ship and helicopter transects found that this species occupied areas with ice cover from 50 to 80%, with sightings only occurring from the air (Vacquie-Garcia et al. 2017). This survey suggests that several hundred bowhead whales occur in the surveyed area, rather than the tens of animals estimated for this population some decades ago. These findings, in combination with surprisingly high occurrences of bowhead whale calls on passive acoustic devices located in Eram Strait gives cause for optimism regarding an increasing trend in this population (Wiig et al. 2007; Wiig et al. 2010; Moore et al. 2012; Stafford et al. 2012). Elsewhere in their range, this species has been shown to be remarkably flexible in their selection of prey and feeding areas and some populations are increasing (Zeh & Punt 2005; Moore et al. 2010; Nielsen et al. 2015). However, Svalbard bowhead whales appear to exhibit a very strong preference for sea-ice habitats where high concentrations of their Arctic calanoid prey occur (Lydersen et al. 2012; Falk-Petersen et al. 2015) and are therefore likely currently experiencing range contraction.
Similar to the situation for bowhead whales, narwhal habitat was impossible to infer in this study because of the small number of sightings in the open-water areas covered by the MMSDB. This species is tightly associated with sea ice, using leads and cracks deep into heavily ice-covered areas (e.g., Laidre et al. 2008; Reeves et al. 2014). There have been very few studies on the distribution, habitat use and general ecology of narwhals in the Svalbard region (Gjertz 1991; Lydersen et al. 2007) compared to the Canadian Arctic and West Greenland (e.g., Finley & Renaud 1980; Heide-Jorgensen et al. 2002; Laidre et al. 2004; Laidre & Heide-Jorgensen 2005; Dietz et al. 2008; Richard et al. 2010; Reeves et al. 2014). Satellite tracking of a few young narwhals in Svalbard revealed that they do sometimes occupy coastal ice-free areas in late summer, diving into deep trenches to the north-east of Nordaustlandet (Lydersen et al. 2007), which provide some of the only deep-water areas in the region. This species is known to dive to depths of 1000 m when hunting their preferred prey, Greenland halibut (Reinhardtius hippoglossoides) in other Arctic areas (Heide- Jorgensen et al. 2015). A ship and helicopter survey conducted in 2015 in the Marginal Ice Zone north of Svalbard estimated close to 1000 narwhal in the surveyed area (Vacquie-Garcia et al. 2017), which represents less than half of the habitat available between Svalbard and the eastern side of Franz Josef Land. This species was found deep into the ice, extending as far as the survey lines were flown (185 km), so their complete distribution was certainly not covered. Their extreme affiliation with heavy sea ice explains why few open-water observations are reported to the MMSDB.
Minke whales, fin whales, humpback whales and small dolphins were the most frequently observed taxa in Svalbard waters, and all exhibited broad and overlapping habitat use. However, some key differences were also observed. Minke whales covered the broadest area from coastal regions to beyond the continental shelf. Humpback whales tended to occur in shallower waters over the continental shelf, and fin whales generally occupied steeper and deeper regions. These three baleen whales co-occur throughout much of their global range by various forms of niche partitioning, such as preferred depths for foraging and specific food preferences (e.g., Goldbogen et al. 2011; Friedlaender et al. 2015). The diet of these three cetaceans in the Barents Sea is known to consist primarily of krill, although they also target capelin (Mallotus villosus), herring (Clupea harengus), polar cod and amphipods to varying degrees (Lydersen et al. 1991; Skern-Mauritzen et al. 2011). Recent increases in the temperature and salinity of Icelandic waters, and subsequent shifts in the fish assemblages have resulted in minke whales shifting their diet from high quantities of krill and capelin to larger benthic gadoids, primarily Atlantic cod (Gadhus morhua) and haddock (Melanogrammus aeglefinus) (Vikingsson et al. 2014). Small dolphins showed a particular preference for the continental slope and did not occur often in coastal environments. They are known to target capelin in this region, and shift their distribution to track this fish species (Skern-Mauritzen et al. 2011). However, they also eat a wide range of other small fish, including herring and some gadoids (Fall & Skern-Mauritzen 2014). The increased abundance of boreal fish species in Svalbard (Berge et al. 2015; Fossheim et al. 2015) and plasticity in prey choice by these seasonally resident balaenopterids and dolphins are likely to lead to an increased presence of these cetaceans in Svalbard waters.
Blue whales across the world are most frequently sighted in association with continental slopes. This was the case early (2002-08) in the study period reported herein. Upwelling in these regions stimulates high productivity and supports large densities of krill (Fiedler et al. 1998; Gill 2002; Wall et al. 2009). However, in the latter years (2009-2014) of the current study, blue whales used coastal and Hord environments more often. This shift in location is almost certainly driven by incursions of warm Atlantic Water and associated krill into Svalbard’s Hords in recent years, concomitant with decreased formation of coastal sea ice in winter (see Eriksen & Dalpadado 2011; Buchholz et al. 2012; Dalpadado et al. 2016).
Sperm whales and northern bottlenose whales also occurred in Svalbard in modest numbers. The habitat model for sperm whales suggested suitable areas along 1000 m isobaths and in deeper waters along the continental slope, consistent with their known use of deep water areas associated with their cephalopod prey throughout their range (e.g., Weir et al. 2001). The northern bottlenose whale distribution in Svalbard waters was consistent with studies of this species elsewhere in the North Atlantic, with good habitat identified beyond the 1000 m isobath (Reeves et al. 1997; Weir et al. 2001; Whitehead & Hooker 2012). Both sperm whales and northern bottlenose whales rely heavily on the squid species Gonatus fabricii in Norwegian waters, adults of which have been recorded as deep as 2700 m (Bjorke 2001).
Sei whales were sighted only infrequently in this study. However, they were seen at the northern tip of Spitsbergen, their most northerly recorded sighting (see Reilly et al. 2008a), suggesting a possible range expansion might be taking place. Although sample size for sei whales in the MMSDB deems it appropriate to be cautious, the habitat suitability model for sei whales suggested that this species occupies a habitat similar to minke, fin and humpback whales; sei whales tend to prefer waters along the continental slope, concordant with the known habitat use of this species in the North Atlantic (Prieto et al. 2014). Sei whales consume copepods and krill (Jonsgard & Darling 1977); the increase in krill abundance and diversity in Svalbard waters (Buchholz et al. 2012) could provide a food supply for this species to exploit in the future.
Killer whales were observed throughout the study area right up to the ice edge, although they were only seen infrequently. Killer whales were hunted along the Norwegian coast until 1981 (0ien 1988). This species is known to prey on all of the ice-associated whales in other Arctic locales (Ferguson et al. 2010; Higdon et al. 2012), and white whales in Svalbard have been suggested to exhibit behaviours that could be anti-predatory responses to killer whales, including spending most of their time very close to shore and not vocalizing often (Lydersen et al. 2001; Karlsen et al. 2002). In the north-eastern Canadian Arctic, killer whale sightings have increased exponentially since 1900, associated with the decreasing sea-ice extent opening up new hunting grounds (Higdon & Ferguson 2009). Given their occasional presence throughout Svalbard waters, as well as the crudely modelled suitability of a broad array of habitats in the region, especially dose to the ice edge, there is potential for killer whales to become a significant predator in this ecosystem in the future (e.g., Baylis et al. 2015).
Long-finned pilot whales, harbour porpoises and white-sided dolphins were all reported to the MMSDB during the study period, but must currently be considered vagrants to the area. Long-finned pilot whales have previously been recorded from Bjornoya, and a solitary harbour porpoise was sighted in association with white whales on the north coast of Spitsbergen between 2005 and 2007 (Kovacs et al. 2009). These three species are not currently regular residents of Svalbard, however given the ongoing increases in their preferred prey species in these waters (Gannon et al. 1997; Eriksen & Dalpadado 2011; Sveegaard et al. 2012; Golikov et al. 2013; Berge et al. 2015), future changes in the abundance and northerly distribution of these cetaceans are likely in the coming decades.
This study provided some indication that latitudinal shifts may be occurring in the ranges of several cetacean species in this region. Such shifts follow a broader pattern of cetacean range changes associated with increasing water temperatures across the North Atlantic (e.g., MacLeod et al. 2005; MacLeod 2009; Lambert et al. 2014). This study provides the first indication for an increased presence of blue whales in Svalbard waters. The northward expansion of this species is consistent with observations from the central north-east Atlantic
between 1987 and 2001 (Pike et al. 2009), and in Icelandic waters from 1996 to 2004 (Vikingsson et al. 2015). Northward shifts were also observed in minke whales, fin whales, humpback whales (all of which were also increasingly sighted in coastal environments in recent years in this study), sei whales and sperm whales. These observed shifts coincide with recent increases in the Atlantic Water temperatures reaching this High-Arctic region, and concomitant newfound abundance of many of the whales’ prey species in the waters around and within the Hords of Svalbard. Boreal fish species including Atlantic cod, haddock, herring and mackerel (Scomber scombrus), none of which was found regularly in Svalbard waters in previous decades, have all been increasing in abundance since the 2000s and are now regularly caught in the Hords of west Spitsbergen (Berge et al. 2015). Krill is increasing in abundance and has recently experienced reproductive success in Svalbard waters (Eriksen & Dalpadado 2011; Buchholz et al. 2012). Several squid species have also expanded their ranges into the Arctic (see Golikov et al. 2013). Cetacean species in the North Atlantic are likely to continue undergoing range changes associated with their prey’s response to increasing water temperature, and the study herein provides an indication that this is already occurring around Svalbard. The extent of these distribution changes will be determined somewhat by the thermal preferences of the species; but likely more importantly by competition with ecologically similar species from slightly different thermal regimes, and presence of their temperature-dependent prey (see MacLeod 2009 and references therein).
The large numbers of seasonally resident species seen in the Svalbard region and the observed expansion into higher latitudes could also reflect population increases within these species, as has been suggested for cetaceans in the Southern Chukchi Sea (Clarke et al. 2013). Historical whaling records indicate the presence of several seasonally resident cetaceans in Svalbard waters during the early 20th century. There are records of blue whales being caught along the west coast of Spitsbergen, including within IsHorden and north of 80°N in the 1930s (Smith no date). Humpback whales and fin whales were also caught in this region, and there is at least one record of a sperm whale having been caught at 80.68°N, 12.33°E in 1938 (Smith no date) The implementation of the International Whaling Commission’s “moratorium” on commercial whaling in 1986 led to cessation of the vast majority of whaling across the globe, allowing for recovery of many populations (e.g., Sears & Calambokidis 2002; Magera et al. 2013; Rocha et al. 2015; Thomas et al. 2016). For example, North Atlantic blue whale and humpback whale populations have been estimated to be increasing at 5% and 3% per year, respectively, in recent decades (Sears & Calambokidis 2002; Reilly et al. 2008b). It is therefore difficult to differentiate between range changes and population recovery.
This study provides the first comprehensive baseline description of the cetacean assemblage occupying Svalbard waters from March through November. The year-round Arctic resident white whales were observed exclusively in coastal habitats, in accordance with their previously documented use of tidal glacier fronts as feeding grounds, and their use of shallow near shore corridors during transits between glaciers. Bowhead whales and narwhals were observed very infrequently, likely on account of their extreme affinity for heavy sea-ice cover in this region, beyond the reach of most reporting vessels. Seasonally resident minke whales, fin whales, humpback whales and small dolphins were all observed often around Svalbard. The continental slope was an important habitat for these species, however they were also observed across the whole study areas in somewhat overlapping habitats, which is likely explained by their flexibility in prey choice and niche partitioning.
Many of the seasonally resident cetaceans showed possible signs of range expansions during the decadal+ study period. This was most pronounced for blue whales, the relative abundance of which increased dramatically in Svalbard waters over the study period. Minke whales, fin whales and humpback whales also used increasingly higher latitudes (and coastal habitats) as the study period progressed. In this study, sperm whales, northern bottlenose whales and sei whales were all seen at latitudes further north than their reported known ranges; definitions of northerly range boundaries for these cetacean species therefore need to be updated by the IUCN and other organizations presenting global or regional overviews. Such range expansions are likely linked to the increasing water temperatures and the concomitant increase in temperate prey species, including krill and various fish, in these high latitude coastal habitats. The continued use of the MMSDB and similar databases can provide a cost-effective means for collecting observational data to track major shifts in distribution and abundance in locations where there are sufficient observational platforms, and help identify areas where research effort should be focussed.
We thank the AECO, the Marine Research Institute, the Norwegian Coast Guard and all other contributors to the MMSDB. SST data herein is courtesy of NASA GSFC, Ocean Ecology Laboratory, Ocean Biology Processing Group. ASI algorithm SSMI/S sea-ice concentration data, originally computed at and provided by the Institut Framjais de Recherche pour 1’Exploitation de la Mer, Brest, France, were obtained as five-day median- filtered and gap-filled product for 2002-2014 from the Integrated Climate Data Center (icdc.cen.uni-hamburg. de), University of Hamburg. Siri Uldal created the database software frames. Comments from Dr Janet Clarke and an anonymous reviewer helped us improve the manuscript.
No potential conflict of interest was reported by the authors.
This work was supported by the Norwegian Research Council ICE Whales project (no. 244488/E10), and the Norwegian Polar Institute and the World Wide Fund for Nature (WWF).