- Open Access
Flexibility of foraging strategies of the great skua Stercorarius skua breeding in the largest colony in the Barents Sea region
Frontiers in Zoologyvolume 15, Article number: 9 (2018)
Foraging strategies of seabird species often vary considerably between and within individuals. This variability is influenced by a multitude of factors including age, sex, stage of annual life cycle, reproductive status, individual specialization and environmental conditions.
Using GPS-loggers, we investigated factors affecting foraging flight characteristics (total duration, maximal range, total distance covered) of great skuas Stercorarius skua of known sex breeding on Bjørnøya, Svalbard, the largest colony in the Barents Sea region. We examined influence of sex (females are larger than males), phase of breeding (incubation, chick-rearing), reproductive status (breeders, failed breeders) and bird ID (they are known for individual foraging specialization). Our analyses revealed that only bird ID affected foraging flight characteristics significantly, indicating a high degree of plasticity regardless of sex, reproductive status or phase of breeding. We recognized three main groups of individuals: 1) those preying mainly on other seabirds in the breeding colonies (6%), 2) those foraging at sea (76%) and kleptoparasiting other seabirds and/or foraging on fish and/or offal discarded by fishing vessels, and 3) those alternating between preying on other seabirds in breeding colonies and foraging at sea (18%). Despite marked size sexual dimorphism, we found no apparent sex differences in flight characteristics. Birds after egg- or chick-loss and thus not constrained as central foragers did not modify their foraging flights.
Great skuas breeding on Bjørnøya displayed a high degree of plasticity regardless of sex, reproductive status or phase of breeding. We recognized groups of individuals regularly preying in the seabird colonies, foraging at sea, and alternating between both strategies. This suggests foraging specialization of some individuals.
Foraging strategies of seabird species often vary considerably between and within individuals. This variability results from targeting usually widely-dispersed, spatiotemporally patchy prey and the influence of a multitude of other factors including age, sex, stage of annual life cycle, reproductive status, individual specialization and environmental conditions. These factors have major implications for our understanding of seabird ecology, because they affect the use of resources, level of intra-specific competition and niche partitioning .
Sexual differences in foraging behaviour are apparent in many seabirds. This may reflect habitat specialization and/or sex-specific nutrient requirements [2,3,4] and reduces competition between males and females. Sexual size dimorphism is thought to play a functional role in flight performance and is used to explain differences in the at-sea distribution of male and female seabirds . Partial or complete sexual foraging segregation has been reported in several seabird groups including both sexually dimorphic and monomorphic species (e.g. [6, 7]). Sexual segregation in foraging may be especially expressed in years of greater environmental stochasticity, when food availability is reduced .
Age may affect foraging strategy considerably, and immature birds avoid competition with adults by feeding in suboptimal habitats which increases their foraging time (e.g. [9,10,11]). Younger or less experienced breeders may forage less efficiently and/or at a lower trophic level, which may be reflected in reduced breeding performance compared to older birds [9, 12,13,14]. Immatures may be far more exploratory and lack route or foraging site fidelity compared to adult breeders . Inferior foraging success in younger individuals reflects their poorer skills in identifying, catching or handling prey or in selecting suitable locations, and weaker motor control or physiological fitness .
Within a single breeding season, seabirds may adopt a temporally flexible foraging strategy to satisfy different energy and time demands of incubation and chick-rearing and both their own energetic requirements and those of their offspring. The changing degree of central-place constraint in successive phases of the breeding and seasonal changes in prey availability may affect distribution, activity patterns or diet of seabirds [4, 15, 16]. Foraging strategies may also vary within a single breeding season as a response to fluctuations in prey availability that is driven by prey biology, environmental conditions and/or prey depletion near the colony due to intense foraging (“Storer-Ashmole’s halo” [17,18,19,20]). Moreover, some species (albatrosses, petrels, little auk Alle alle) adopt an unimodal foraging strategy during incubation (trip of similar duration) and bimodal strategy during chick-rearing, alternating short trips to nearby locations to collect food for chicks, and long trips to further location mainly for self-feeding (e.g. [21,22,23,24]).
One of the less studied factors affecting foraging strategy is reproductive status. Failed breeders often continue to associate with the colony, operating as central-place foragers but expand their foraging areas . They may be partially or completely segregated from breeders, probably to avoid competition .
Many marine birds and mammals show individual feeding specializations in terms of distribution, behaviour, diet or other aspects of resource acquisition that remain after accounting for the group effects of sex, age, breeding stage and status [26, 27]. Individual specialization may be expressed in fidelity to feeding sites, consistency in foraging trip characteristics, diving patterns and other at-sea activity, habitat use, specific prey items or trophic level in the short or long term . Foraging specializations are probably learned during individual exploratory behaviours early in life, which then become canalized with age and experience [28, 29].
In this study, we investigated foraging flights of great skuas Stercorarius skua breeding on Bjørnøya (Svalbard), one of the northernmost breeding areas within this species’ range . The great skua is a dietary generalist exploiting a wide range of prey. It searches for and catches food exclusively on the wing, mainly by splash-diving onto surface fish shoals. It may also kleptoparasitise and prey on seabirds or forage on discards from fishing vessels [31, 32]. The great skua is an important predator of seabirds during the breeding season taking eggs, chicks and adult birds. Dietary specialization of great skua pairs and colonies has been documented in several studies [5, 33,34,35,36]. Remains of seabirds (mainly black-legged kittiwakes Rissa tridactyla and northern fulmars Fulmarus glacialis) and fish were found in 98 and 38% of pellets of great skuas breeding on Bjørnøya in 2008–2009. Considerable fractions (33 and 61% in 2008 and 2009, respectively) of pairs breeding there showed dietary specialization on seabirds (≥70% pellets contained only seabird remains) . The great skua breeds mainly colonially on flat or gently sloping ground. The female lays 1–2 eggs that are incubated for 26–32 days by both sexes, although mainly by the female. The semi-altricial chicks stay in the nest area for 40–51 days after hatching. The chicks are guarded mainly by the female and fed by both parents . The great skua and other skuas (Stercorariinae), like other birds with a raptorial lifestyle, display sexual size dimorphism with females being larger than males .
Despite extensive studies on the great skua foraging and breeding ecology (e.g. [35, 38, 39]), foraging flights during the chick-rearing period have been studied so far only in Shetland  and St. Kilda, Outer Hebrides .
The aim of this study was to investigate factors affecting foraging flight characteristics (maximal range, total distance covered and total duration) of great skuas breeding on Bjørnøya. We considered sex, stage of breeding (incubation or chick-rearing), reproductive status (breeders vs. failed breeders) and bird ID. We expected that:
Because of various time budgets of incubating and chick-rearing individuals, birds during incubation would have longer foraging flights than during chick-rearing given the necessity of adults to feed chicks in regular intervals;
Due to various time budgets and duties, failed breeders free of central place forager constraint would have longer foraging flights to more distant areas enabling them to feed beyond the cost-effective flight distance to the colony;
Regardless of breeding status and phase of the breeding period and considering sexual size dimorphism, larger females will perform flights of longer duration and range compared to smaller males;
Individual great skuas that specialized in foraging on local resources (other seabirds on the same island) would perform shorter trips characterized by high repeatability in utilisation densities.
Bjørnøya (74°30′N, 19°01′E) is a 178 km2 island in the western Barents Sea at the south-western edge of the shallow Spitsbergen Bank (Fig. 1). This bank is, in summer, characterized by a mixture of cold Arctic water and melt-water from the Barents Sea ice. The waters around Bjørnøya are characterized by non-stratified Spitsbergen Bank water near the island. This water is surrounded to the south and west by a frontal zone where cold fresher Arctic water mixes with warm saline Atlantic water that is found in the deeper parts to the west, south and southeast of Bjørnøya. The zone where the two water masses meet is called the Polar Front. Its location is determined mainly by bathymetry and is oriented parallel to the slopes of the shallow Spitsbergen Bank with the surface expression following the 100–250 m isobath [41,42,43]. Bjørnøya accommodates large colonies of seabirds including black-legged kittiwakes, Brünnich’s guillemots Uria lomvia, common guillemots U. aalge and little auks . Great skuas were first observed breeding on Bjørnøya in 1970 [30, 44] but the population size has since increased rapidly, becoming the largest colony in the Barents Sea region [currently estimated at 750–1000 breeding pairs (H. Strøm unpubl. data) and one of the northernmost within the species’ range.
To characterize foraging flights of great skuas, we used global positioning system (GPS) loggers (Harrier, Skua and Uria models, Ecotone, Sopot, Poland) recording time, position and momentary speed. We deployed 21 loggers on one pair member from 21 nests. We captured the birds on nests during the incubation period in June 2014 in the colony in Flakmyrvatna area in NW part of the island. The study area contained 107 breeding pairs. The logger weight (including attachment, Harrier: 19.1–19.3 g; Skua: 26.3 g; Uria: 11.2 g) was equivalent to 1.2–1.5% (Harrier, N = 14 individuals), 1.7–2.1% (Skua, N = 4) and 0.7–0.9% (Uria, N = 3) of the bird’s body mass. The lightest logger type (Uria) was attached to the tail feathers. Other types of loggers were attached with full harness (12 Harriers, 4 Skuas) or with a loop around head and tail (2 Harriers). The GPS-loggers used a bidirectional radio link with base stations installed in the colony, allowing remote data download. To save battery life, the base station automatically switched off the loggers while they were within the download range of the base station. The loggers restarted to record positions when birds left the base station signal range. Sampling interval was set to 15 min; however, when the battery voltage was low, the interval increased to 60 min. The field-tested accuracy of the GPS receiver was ±10 m for 95% of positions. We analysed only records with sampling rates of up to 30 min. Within longer intervals between recorded positions, birds may have had time to return to the colony and start a new trip.
We sexed the studied birds molecularly based on DNA extracted from blood following a modified protocol of . We amplified CHD genes from extracted DNA by PCR using primers 2550F and 2718R . Analyses were carried out at the Norwegian University of Science and Technology, Trondheim.
We received data from 17 of 21 individuals. In the case of three loggers, communication with the base station failed and in the last case, the logger changed recording intervals, forcing us to exclude data from the analyses. In total, 402 trips of 17 birds were analysed with 4 to 87 trips recorded per individual (Table 1).
To monitor the status of GPS-logger equipped individuals and their nests, we controlled the nest content 4–10 d days after the logger deployment. We found that all studied great skuas failed breeding. Among the 21 nests, 53% failed at the incubation stage. In the remaining 47%, at least one chick hatched but all died after 1 to 17 days. Similar breeding failure was recorded in the whole colony area, where no birds were caught. No chicks reached the age of 20 days in 107 nests monitored. We do not know the precise reason of this failure, but we observed frequent conspecific predation which can be related to unfavourable feeding conditions (HS – unpubl. data).
Based on the geographical positions recorded by the GPS-loggers, we analysed the following foraging flights characteristics: (1) maximum range of flights – distance from the colony to the distal point reached on each foraging trip; (2) the total distance covered (km) as the sum of the distances (km) between all GPS positions along each individual’s track; in the case of incomplete trips (i.e. without the first or last position), the missing part was estimated based on distance of the first/the last position nearest to the colony; (3) total trip duration, defined as the time between departure and return to colony; in the case of the incomplete trips, the lacking part was estimated based on the momentary speed and distance of the first/the last position nearest to the colony. To characterize the whole area utilized by studied individuals, we calculated minimal convex polygon (MCP).
Considering the sample size for particular studied features (sex, reproductive status, phase of breeding), we were able to perform analyses including sex and ID of birds for the incubation periods for 17 individuals. For nine individuals, we compared flight characteristics during incubation before and after eggs loss. For one male and one female, we were able to compare flight characteristics during incubation, chick-rearing and after chick loss.
To characterize flights of males and females during the incubation period, we used Conditional Inference Tree (CIT). This is a non-parametric class of regression tree, examining the relationship between multiple explanatory variables and a single response variable using a recursive binary-partitioning process. Model outputs produce an ‘inverted tree’, in which the root at the top contains all observations, which is divided into two branches at the node. The aim of splitting the data at each step is to establish groups that had a between-variation as large, and within-variations as small, as possible. The node provides information about the explanatory variable name and its probability value. Branches are further split into two subsequent nodes and so on. . CIT uses a machine learning algorithm to determine when splitting is no longer valid using statistically-determined stopping criterion, an a priori p value . CIT is robust to typical regression problems such as over-fitting, collinearity, and bias with regard to the types of explanatory variables used [47, 48]. We conducted CIT analyses in R software  using party package . We checked significance of nodes in CIT analyses using structural change test implemented in strucchange package in R .
For data from the incubation period with satisfactory sample size, we calculated utilization distributions (UD, a probability distribution constructed from data providing the location of an individual in space at different points in time) for home ranges (95% kernel density) and core ranges (50% kernel density) of males and females in Geospatial Modelling Environment (GME) ver. 0.7.4.0 software (www.spatialecology.com/gme/) using kde function using plugin bandwidth selection. To investigate overlap between utilization distributions (i.e. 95 and 50% kernel density isopleth) of males and females during the incubation period before egg loss, individuals before and after egg loss, and individuals during incubation and chick-rearing periods, we used adehabitat package in R  with algorithm BA, i.e. the Bhattacharyya’s affinity, a statistical measure of affinity between two populations. It ranges from zero (no overlap) to 1 (identical UDs) . To measure individual consistency in UD between particular trips, we calculated the mean Bhattacharyya’s affinity of all pairwise combinations of the recorded trips and then calculated the mean value for particular individuals.
To estimate how many flights were to seabird colonies at the south end of Bjørnøya, we counted all flights within a 1 km buffer around the seabird cliff. We mapped data from the GPS loggers, performed spatial analyses and produced all figures with maps using ArcMap/ArcGIS 10.3.1 (Environmental Systems Research Institute, Redlands, CA, USA).
Great skuas breeding on Bjørnøya performed foraging flights lasting 0.8–224 h (median 5.2 h, N = 17 individuals), reaching areas at maximal distance from the colony ranging from 1.6 to 1003 km (median 88.1 km), covering total distances ranging from 3.6 to 2790 km (median 199 km) (Table 2, Fig. 1). Minimal convex polygon of all GPS position recorded ranged from 0.03 to 16,946 km2 (median 141.7 km2) (Table 2).
Factors affecting foraging flights characteristics during incubation
To investigate factors affecting foraging flight characteristics during incubation (we considered only the period before egg loss), we used bird ID and sex as explanatory variables. The Conditional Inference Tree (CIT) results indicate that bird ID best characterized total foraging trip duration variability. CIT recognized three groups of individuals with various mean foraging times: 29.2 h (node 5 represented by 4 individuals), 13.8 h (node 4 represented by 3 individuals) and 6.3 h (node 3 represented by 10 individuals) (Fig. 2a).
In the case of maximal range of foraging flights, bird ID characterized the best observed variability. Two main nodes were recognized: 106.6 km (node 2 represented by 12 individuals) and 212.4 km (node 3 represented by 5 individuals) (Fig. 2b).
CIT results indicated that bird ID best characterized total distance covered. Two main nodes: 246.8 km (node 2 represented by 12 individuals) and 613.8 km (node 3 represented by 5 individuals) (Fig. 2c) were recognized.
In the case of maximal convex polygon of all recorded GPS positions, bird ID characterized the best observed variability with two main nodes: 81.4 km2 (node 2 represented by 13 individuals) and 2085.1 km2 (node 3 represented by 4 individuals) (Fig. 2d).
In the case of individuals incubating eggs (11 females and 6 males), the 95% utilisation distributions showed Bhattacharyya’s affinity (BA) 0.80 between females and males (Fig. 3). In the case of core range, BA values were lower: 0.42 (Fig. 3b).
Factors affecting foraging flights of birds that lost eggs
To characterise foraging flight of great skuas during the incubation period (including periods before and after egg loss), we used sex, status (BEL – before egg loss, AEL – after egg loss), and bird ID as explanatory variables. CIT results indicated that bird ID and sex best characterized total foraging trip duration variability. CIT recognized two main groups – node 2 represented by 5 individuals with mean total trip time 22.8 h, and node 3. The latter node was split into two nodes: node 7 (represented by 2 individuals and mean trip duration 12.1 h) and node 4. The latter node was split between sexes. In that group containing 3 individuals, females performed flights of shorter duration (mean 5.0 h) than the male (8.0 h) (Fig. 4a).
In the case of maximal range of foraging flights, bird ID and sex characterized the best observed variability. Two main nodes were recognized: node 2 represented by 6 individuals and characterized by more distant flights (mean 174.3 km) and node 3. The latter node represented by 3 individuals was split by sex. Females were characterised by higher maximal range of foraging flights (mean 96.4 km) compared to the male (58.4 km) (Fig. 3c, 4b).
In the case of CIT total distance covered during foraging flights, bird ID and sex characterized the best observed variability. Two main nodes were recognized: 518.8 km (node 2 represented by 6 individuals) and node 3 represented by 3 individuals. The latter node was split by sex. Females were characterised by higher total distance covered during foraging flights (mean 230.8 km) compared to the male (154.1 km) (Fig. 4c).
In the case of maximal convex polygon of all recorded GPS positions, bird ID characterized the best observed variability. Two main nodes were recognized: 2207.7 km2 (node 2 represented by 5 individuals) and 262.9 km2 (node 3 represented by 4 individuals) (Fig. 4d).
In the case of nine incubating individuals that lost eggs, the 95% utilisation distributions showed a Bhattacharyya’s affinity (BA) of 0.74 between the phases before and after egg loss. In the case of core range, the BA value was lower - 0.36.
Factors affecting foraging fights of individuals that lost chicks
To characterise foraging flights of great skuas during the incubation and chick-rearing periods (including periods before and after chick death), we used status (INC – incubation, BCD – before chick death, ACD – period after chick death), and bird ID (two individuals: I – female and O - male) as explanatory variables. CIT results indicated that bird ID best characterized all characteristics of the flights (Fig. 3c, 5).
We found that in the case of two individuals tracked throughout incubation (INC) and chick-rearing (BCD – before chick death, ACD - after chick death), the 95% utilisation distributions showed Bhattacharyya’s affinities (BA) of 0.71 between INC and BCD, 0.56 between INC and ACD, and 0.36 between BCD and ACD. In the case of core range, the BA values were 0.30 between INC and BCD, 0.25 between INC and ACD, and 0.14 between BCD and ACD.
One individual, the female HAR 30, tracked throughout incubation, flew exclusively towards the seabird colonies at the south end of the island (Fig. 6a). The majority of her flights (6 of 7; 86%) targeted the seabird colonies (i.e. within the 1 km buffer zone around the colony). Her flights were characterized by the shortest total duration time, the lowest range, total distance covered and MCP area of all the studied individuals (Table 2). Three other individuals (females HAR 17, HAR 20 and male HAR 26) of 17 studied (17.6%), regularly visited the same colonies, but also flew to other areas. Flights within the 1 km buffer around the colony made up 13%, 29% and 11% of all recorded flights of HAR17, HAR20 and HAR26, respectively. Remaining individuals (76%) flew out to sea, mainly in a NW direction from the colony (Fig. 6).
The consistency of utilisation density of particular flights (i.e. mean BA values per individual) ranged from 0.36 to 0.79 in the case of home range (HR, 95% kernel density) and from 0.06 to 0.32 in the case of core range (CR, 50% kernel density) (Table 1). The trips of HAR 30 were characterized by the highest consistency (the highest values of mean Bhattacharyya’s affinity) of both HR and CR. Almost half (47%) of the studied great skuas had mean BA values for HR higher than average value for all individuals. In the case of CR, only 29% individuals were characterized by BA values higher than average value.
Our study revealed that individual specialization had the greatest effect on foraging flights of great skuas breeding on Bjørnøya. Three main strategies were identified: foraging on other seabirds, foraging at sea and a generalist strategy mixing the two. Similarly, great skuas from Unst, Shetland exhibited dietary specialization: a small proportion fed almost exclusively upon seabirds, a small proportion fed as generalists and most feed on fishery discards . Also the southern hemisphere species, the brown skua Stercorarius antarcticus lonnbergi breeding on King George Island (Maritime Antarctica) adopted three main strategies; foraging in a penguin colony, in a storm-petrel colony and at sea . Individual specialization may have a selective advantage where resources are to some extent predictable. The specialist foraging preference is probably passed to offspring through learning as they need to develop particular skills to successfully pursue different foraging modes including kleptoparasitism, predation of selected prey species or scavenging [1, 35, 54].
The study from Unst, Shetland revealed that great skuas specialized in seabird predation were able to spend less time foraging than individuals feeding predominantly on fish in the open sea . In concordance with this, the total trip duration of the female from Bjørnøya that specialized in foraging in the seabird colonies (median 2.1 h) was similar to the mean values 1.1–3.0 h reported for the radio-tracked seabird specialists from Shetland . The proportion of great skuas foraging mostly in seabird colonies in our study (5.8%) was lower than the frequency of seabird specialists calculated based on pellet composition analyses performed in the same colony, i.e. 33% in 2008 and 62% in 2009 . It is, however, possible that great skuas also killed seabirds at sea, not only in the colonies.
Great skuas breeding on Bjørnøya flew further (up to 1003 km) than birds from Foula, Shetland (up to 219 km) . On Foula, the frequency of nearby foraging flights (up to 40 km) was higher (~ 40%) compared to Bjørnøya (15%) (Table 3). This may be attributed to a lower proportion of terrestrial foraging areas on Bjørnøya.
The mean foraging trip duration in our study (5.2 h) was similar to 4.95 h reported for great skuas from Unst, Shetland that specialized in feeding on fish . However, the minimum convex polygons of great skuas breeding on Unst (mean values 1.0–14.4 km2) were considerably smaller than those in our study (median values 29.4–4629.7 km2). This may be explained by the higher proportion of individuals specialized as local seabird predators on Unst (20%) than on Bjørnøya (6%) ). Moreover, some pairs of great skuas on Unst defended feeding territories within a section of seabird colony  and thereby diminished considerably their minimum convex polygon.
The majority of the individuals studied on Bjørnøya flew out to sea, over the shallow shelf area and along the shelf break slope (Fig. 1). The Svalbard Bank, with a minimum depth of less than 40 m is considered as the most productive area in the Barents Sea with primary production about 2–3 times higher than in the adjacent, deeper waters . The shelf break zone is also a good foraging habitat for zooplanktivorous and piscivorous fish and birds including black-legged kittiwakes, Brünnich’s guillemots, common guillemots and little auks [43, 57,58,59]. It is thus possible that the studied great skuas may kleptoparasite or/and prey on seabirds foraging there or returning to the colony. A great skua dietary study from Bjørnøya from 2008 to 2009 revealed that most of fish otoliths found in pellets originated from Gadiformes fish . Most of the species within this group are mid-water or bottom dwelling species  and thus probably only available as discards from fishing vessels . This source of food is a locally important diet component of the great skua breeding in the Outer Hebrides  or Shetland . The fishing activity in the Barents Sea in 2014 included areas around Bjørnøya , but the level of discards is unknown . It has been reported, however, that discards of Atlantic cod Gadus morhua from the Norwegian shrimp fishery in the Barents Sea consists mainly of 1- and 2-year-olds  (i.e. with mean body length 10.4 and 18.4 cm, respectively ) and thus optimal food for great skuas. Great skuas from Bjørnøya may also prey directly on pelagic fish, such as the capelin Mallotus villosus during foraging flights at sea by dipping or surface-seizing. This is supported by an observation of adult and chick regurgitates consisting of capelin remains (H. Strøm unpubl. data), although it is also possible that these remains originated from kleptoparasitism. In Shetland, great skua diet composition varied among the colonies and was dependent on the colony size. The majority of individuals from large colonies fed on fish, including discards, and only a small proportion specialised in killing seabirds . With regard to the size and the rapid growth of the Bjørnøya population, it is likely that they have adapted to feed on widely abundant fish such as capelin or fish discarded from fisheries.
Almost half of the studied birds were consistent in their home range areas (mean values of BA higher than the average) suggesting frequent foraging in the same areas during consecutive flights. This would be expected of individuals repeatedly visiting seabird colonies or those exploring areas characterized by similar oceanographic features as fronts in the shelf break zone providing opportunity to find their own fish or kleptoparasitise other seabirds gathered there to forage. Individuals alternating between flights to sea and to the seabirds colonies had lower repeatability of utilisation densities values than the average for all individuals studied.
That neither sex nor phase of reproduction of the studied birds affected foraging flight characteristics significantly was unexpected. This may be interpreted as high inter-individual variability in foraging being more important than other sources of variability. The only sex difference in total trip duration found was in the small subgroup of birds (three individuals) during incubation but, in this case, sex effect was also a proxy for bird ID. Our result corroborates a study of great skuas and other skuas (Stercorariinae) that found no support for the theory that sexual dimorphism evolved as a result of specialized roles during breeding (smaller males with lower wing loading foraging more efficiently due to greater agility) .
Surprisingly, after losing eggs or chicks, individuals freed of the restraint of central-place foraging did not modify their flights characteristic as did failed breeding northern giant petrels Macronectes halli and southern giant petrels M. giganteus . Using similar home ranges before and after egg loss may indicate foraging site fidelity by the Bjørnøya birds. It suggests that feeding grounds utilized by the studied great skuas were at least good enough for self-maintenance.
We are aware of possible limitations of our study. Firstly, we did not know the age of the studied individuals. It has been reported that time spent foraging to provide food for chicks increases with age of great skua parents breeding in Shetland . Secondly, the 2014 breeding season was a complete failure for great skuas from Bjørnøya  with no chicks reaching an age beyond 20 days. Furthermore, there was a 26% reduction in the number of breeding pairs (to 107 nests in the whole monitoring area) compared to 2013 (HS, pers. obs.). We suppose that the breeding failure may have been the result of high conspecific predation rate driven by suboptimal feeding conditions on foraging grounds. In 2014, the stock of the only semi-pelagic or pelagic gadoid fish, the polar cod Boreogadus saida, which could have been preyed directly by great skuas or indirectly by kleptoparasiting other seabirds (this fish is frequently preyed by Brünnich’s guillemots and black-legged kittiwakes breeding on or sampled around Bjørnøya [57, 68]), was reduced and its distribution shifted northeastwards, beyond the Barents Sea shelf . The capelin stock at the time of the study was also significantly reduced [62, 69] thereby possibly influencing great skuas directly or/and indirectly. In years when capelin is scarce, the piscivorous guillemots breeding on Bjørnøya may rely more on euphausiids for food . In such years, the kleptoparasitic behaviour of great skuas may be less energetically profitable than in years with a higher contribution of fat-rich fish in their diets. In this context, foraging flights characteristics may also have been different than in other seasons. Finally, logger deployment might be considered as a cause of the breeding failure  but this was unlikely due to the overall breeding failure recorded throughout the colony. Moreover, other studies of great skuas fitted with GPS devices revealed that foraging flights characteristics and territory attendance rates were similar to those of control individuals .
The characteristics of great skua foraging flights on Bjørnøya during a poor breeding season were individual rather than a response to sex or breeding status (after or before egg or chick loss). Almost half of the studied birds were consistent in their home range areas during consecutive flights suggesting foraging specialization. Some individuals fed regularly in the seabird colonies at the southern part of the island, some foraged at sea in areas with oceanographic features as fronts providing opportunity to kleptoparasiting or even direct hunting for fish, and some birds used both foraging strategies.
Our study demonstrates the effectiveness of GPS-tracking as a tool for investigating foraging specialization of great skuas. Future studies could usefully apply the same approach, ideally supported by diet composition estimation (e.g. stable isotope ratios, pellet composition), to investigate consistency and flexibility in the foraging strategies.
Phillips RA, Lewis S, González-Solís J, Daunt F. Causes and consequences of individual variability and specialization in foraging and migration strategies of seabirds. Mar Ecol Progr Ser. 2017;578:117–50. https://doi.org/10.3354/meps12217.
Phillips RA, Silk JRD, Phalan B, Catry P, Croxall JP. Seasonal sexual segregation in two Thalassarche albatross species: competitive exclusion, reproductive role specialization or foraging niche divergence? Proc R Soc B Biol Sci. 2004;271:1283–91. https://doi.org/10.1098/rspb.2004.2718.
Phillips RA, McGill RAR, Dawson DA, Bearhop S. Sexual segregation in distribution, diet and trophic level of seabirds: insights from stable isotope analysis. Mar Biol. 2011;158:2199–208. https://doi.org/10.1007/s00227-011-1725-4.
Hedd A, Montevecchi WA, Phillips RA, Fifield DA, Garthe S. Seasonal sexual segregation by monomorphic sooty shearwaters Puffinus griseus reflects different reproductive roles during the pre-laying period. PLoS One. 2014;9:e85572. https://doi.org/10.1371/journal.pone.0085572.
Catry P, Phillips RA, Furness RW. Evolution of reversed sexual size dimorphism in Skuas and jaegers. Auk. 1999;116:158–68. https://doi.org/10.2307/4089462.
Ismar SMH, Raubenheimer D, Bury SJ, Millar CD, Hauber ME. Sex-specific foraging during parental care in a size-monomorphic seabird, the Australasian gannet (Morus serrator). Wilson J Ornithol. 2017;129:139–47. https://doi.org/10.1676/1559-4491-129.1.139.
Cleasby I, Wakefield E, Bodey T, Davies R, Patrick S, Newton J, et al. Sexual segregation in a wide-ranging marine predator is a consequence of habitat selection. Mar Ecol Prog Ser. 2015;518:1–12. https://doi.org/10.3354/meps11112.
Paiva VH, Pereira J, Ceia FR, Ramos JA. Environmentally driven sexual segregation in a marine top predator. Sci Rep. 2017;7:2590. https://doi.org/10.1038/s41598-017-02854-2.
Daunt F, Afanasyev V, Adam A, Croxall J, Wanless S. From cradle to early grave: juvenile mortality in European shags Phalacrocorax aristotelis results from inadequate development of foraging proficiency. Biol Lett. 2007;3(4):371. https://doi.org/10.1098/rsbl.2007.0157.
Fayet AL, Freeman R, Shoji A, Padget O, Perrins CM, Guilford T. Lower foraging efficiency in immatures drives spatial segregation with breeding adults in a long-lived pelagic seabird. Anim Behav. 2015;110:79–89. https://doi.org/10.1016/j.anbehav.2015.09.008.
Votier SC, Fayet AL, Bearhop S, Bodey TW, Clark BL, Grecian J, et al. Effects of age and reproductive status on individual foraging site fidelity in a long-lived marine predator. Proc R Soc B Biol Sci. 2017;284:20171068. https://doi.org/10.1098/rspb.2017.1068.
Elliott KH, Le Vaillant M, Kato A, Gaston AJ, Ropert-Coudert Y, Hare JF, et al. Age-related variation in energy expenditure in a long-lived bird within the envelope of an energy ceiling. J Anim Ecol. 2014;83:136–46.
Le Vaillant M, Le Bohec C, Prud’Homme O, Wienecke B, Le Maho Y, Kato A, et al. How age and sex drive the foraging behaviour in the king penguin. Mar Biol. 2013;160:1147–56. https://doi.org/10.1007/s00227-013-2167-y.
Le Vaillant M, Ropert-Coudert Y, Le Maho Y, Le Bohec C. Individual parameters shape foraging activity in breeding king penguins. Behav Ecol. 2016;27:352–62. https://doi.org/10.1093/beheco/arv146.
Ito M, Takahashi A, Kokubun N, Kitaysky A, Watanuki Y. Foraging behavior of incubating and chick-rearing thick-billed murres Uria lomvia. Aquat Biol. 2010;8:279–87. https://doi.org/10.3354/ab00229.
Navarro J, Louzao M, Igual JM, Oro D, Delgado A, Arcos JM, et al. Seasonal changes in the diet of a critically endangered seabird and the importance of trawling discards. Mar Biol. 2009;156:2571–8. https://doi.org/10.1007/s00227-009-1281-3.
Hipfner JM, McFarlane-Tranquilla L, Addison B, KA H. Trophic responses to the hatching of offspring in a central-place foraging seabird. J Ornithol. 2013;154:965–70. https://doi.org/10.1007/s10336-013-0962-3.
Elliott KH, Woo KJ, Gaston AJ, Benvenuti S, Dall’Antonia L, Davoren GK. Central-place foraging in an Arctic seabird provides evidence for Storer-Ashmole’s halo. Auk. 2009;126:613–25. https://doi.org/10.1525/auk.2009.08245.
Ashmole NP. The regulation of numbers of tropical oceanic birds. Ibis. 1963;103:458–73. https://doi.org/10.1111/j.1474-919X.1963.tb06766.x.
Deguchi T, Wada A, Watanuki Y, Osa Y. Seasonal changes of the at-sea distribution and food provisioning in rhinoceros auklets. Ecol Res. 2010;25:123–37. https://doi.org/10.1007/s11284-009-0639-9.
Jakubas D, Wojczulanis-Jakubas K, Iliszko L, Darecki M, Stempniewicz L. Foraging strategy of the little auk Alle alle throughout breeding season - switch from unimodal to bimodal pattern. J Avian Biol. 2014;45:551–60.
Jakubas D, Iliszko L, Wojczulanis-Jakubas K, Stempniewicz L. Foraging by little auks in the distant marginal sea ice zone during the chick-rearing period. Polar Biol. 2012;35:73–81.
Catard A, Weimerskirch H, Cherel Y. Exploitation of distant Antarctic waters and close shelf-break waters by white-chinned petrels rearing chicks. Mar Ecol Prog Ser. 2000;194:249–161.
Chaurand T, Weimerskirch H. The regular alternation of short and long foraging trips in the blue petrel Halobaena caerulea: a previously undescribed strategy of food provisioning in a pelagic seabird. J Anim Ecol. 1994;63:275–82. https://doi.org/10.2307/5546.
González-Solís J, Croxall JP, Afanasyev V. Offshore spatial segregation in giant petrels Macronectes spp.: Differences between species, sexes and seasons. Aquat Conserv. 2007;17(S1):S22–S36. doi:https://doi.org/10.1002/aqc.911.
Bolnick DI, Kirkpatrick M. The relationship between intraspecific assortative mating and reproductive isolation between divergent populations. Curr Zool. 2012;58:484–92. https://doi.org/10.1086/343878.
Dall SRX, Bell AM, Bolnick DI, Ratnieks FLW. An evolutionary ecology of individual differences. Ecol Lett. 2012;15:1189–98. https://doi.org/10.1111/j.1461-0248.2012.01846.x.
Guilford T, Freeman R, Boyle D, Dean B, Kirk H, Phillips R, et al. A dispersive migration in the Atlantic puffin and its implications for migratory navigation. PLoS One. 2011;6:e21336. https://doi.org/10.1371/journal.pone.0021336.
Wakefield ED, Cleasby IR, Bearhop S, Bodey TW, Davies RD, Miller PI, et al. Long-term individual foraging site fidelity-why some gannets don’t change their spots. Ecology. 2015;96:3058–74.
Cramp S. The complete birds of the western Palearctic on CD-ROM. In: CDROM. Oxford: Oxford University Press; 1998.
Furness RW. The Skuas. Calton: T & A D Poyser; 1987.
Phillips RA, Catry P, Thompson DR, Hamer KC, Furness RW. Inter-colony variation in diet and reproductive performance of great skuas Catharacta skua. Mar Ecol Prog Ser. 1997;152:285–93.
Phillips RA, Bearhop S, Hamer KC, Thompson DR. Rapid population growth of great Skuas Catharacta skua at St Kilda: implications for management and conservation. Bird Study. 1999;46:174–83.
Votier SC, Bearhop S, Ratcliffe N, Phillips RA, Furness RW. Predation by great skuas at a large Shetland seabird colony. J Appl Ecol. 2004;41:1117–28.
Votier SC, Bearhop S, Ratcliff N, Furness RW. Reproductive consequences for great Skuas specializing as seabird predators. Condor. 2004;106:275–87. https://doi.org/10.1650/7261.
Bayes JC, Dawson MJ, Potts GR. The food and feeding behaviour of the great Skua in the Faroes. Bird Study. 1964;11:272–9. https://doi.org/10.1080/00063656409476077.
Knutsen AE. Diet and breeding success of great skuas (Cathracta skua) on Bjørnøya, Norway, M.Sc. thesis: Norwegian University of Life Sciences (UMB); 2010.
Caldow RWG, Furness RW. The effect of food availability on the foraging behaviour of breeding great Skuas Catharacta skua and Arctic Skuas Stercorarius parasiticus. J Avian Biol. 2000;31:367–75. https://doi.org/10.1034/j.1600-048X.2000.310313.x.
Hamer KC, Furness RW, Caldow RWG. The effects of changes in food availability on the breeding ecology of great skuas Catharacta skua in Shetland. J Zool. 1991;223:175–88.
Votier SC, Crane JE, Bearhop S, de León A, McSorley CA, Mínguez E, et al. Nocturnal foraging by great skuas Stercorarius skua: implications for conservation of storm-petrel populations. J Ornithol. 2006;147:405–13.
Loeng H. Features of the physical oceanographic conditions of the Barents Sea. Polar Res. 1991;10:5–18. https://doi.org/10.1111/j.1751-8369.1991.tb00630.x.
Harris CL, Plueddemann AJ, Gawarkiewicz GG. Water mass distribution and polar front structure in the western Barents Sea. J Geophys Res. 1998;103:2905–17. https://doi.org/10.1029/97JC02790.
Węsławski JM, Stempniewicz L, Mehlum F, Kwaśniewski S. Summer feeding strategy of the little auk (Alle alle) from Bjørnøya, Barents Sea. Polar Biol. 1999;21:129–34.
Anker-Nilssen T, Barrett RT, Bustnes JO, Erikstad KE, Fauchald P, Lorentsen S-H, et al. SEAPOP studies in the Lofoten and Barents Sea area in 2006. NINA Report 249. 2007. https://munin.uit.no/bitstream/handle/10037/11174/article.pdf?sequence=2.
Griffiths R, Double MC, Orr K, Dawson RJG. A DNA test to sex most birds. Mol Ecol. 1998;7:1071–5. https://doi.org/10.1046/j.1365-294x.1998.00389.x.
Fridolfsson A-K, Ellegren H. A simple and universal method for molecular sexing of non-ratite birds. J Avian Biol. 1999;30:116. https://doi.org/10.2307/3677252.
Quinn GP, Keough MJ. Experimental design and data analysis for biologists. Cambridge: Cambridge University Press; 2002.
Hothorn T, Hornik K, Zeileis A. Unbiased recursive partitioning: a conditional inference framework. J Comput Graph Stat. 2006;15:651–74. https://doi.org/10.1198/106186006X133933.
R Development Core Team. R: a language and environment for statistical computing. R foundation for statistical computing. Austria: Vienna; 2015.
Zeileis A, Leisch F, Hornik K, Kleiber C. Strucchange: an R package for testing for structural change in linear regression models. J Stat Softw. 2002;7:1–38. http://epub.wu.ac.at/1124/
Calenge C. The package “adehabitat” for the R software: a tool for the analysis of space and habitat use by animals. Ecol Model. 2006;197:516–9. https://doi.org/10.1016/j.ecolmodel.2006.03.017.
Fieberg J, Kochanny CO. Quanitfying home-range overlap: the importance of the utilization distribution. J Wildl Manag. 2005;69:1346–59. https://doi.org/10.2193/0022-541X(2005)69.
Carneiro APB, Manica A, Trivelpiece WZ, Phillips RA. Flexibility in foraging strategies of Brown Skuas in response to local and seasonal dietary constraints. J Ornithol. 2015;156:625–33.
Phillips RA, Catry P, Silk JRD, Bearhop S, McGill R, Afanasyev V, et al. Movements, winter distribution and activity patterns of Falkland and brown skuas: insights from loggers and isotopes. Mar Ecol Prog Ser. 2007;345:281–91.
Thaxter CB, Ross-Smith VH, Clark JA, Clark NA, Conway GJ, Marsh M, et al. A trial of three harness attachment methods and their suitability for long-term use on lesser black-backed gulls and great Skuas. Ringing Migr. 2014;29:65–76.
Kędra M, Renaud PE, Andrade H, Goszczko I, Ambrose WG. Benthic community structure, diversity, and productivity in the shallow Barents Sea bank (Svalbard Bank). Mar Biol. 2013;160:805–19. https://doi.org/10.1007/s00227-012-2135-y.
Mehlum F, Nordlund N, Isaksen K. The importance of the “polar front” as a foraging habitat for guillemots Uria spp. breeding at Bjornoya, Barents Sea. J Mar Syst. 1998;14:27–43.
Jakubas D, Iliszko LM, Strøm H, Darecki M, Jerstad K, Stempniewicz L. Foraging behavior of a high-Arctic zooplanktivorous alcid, the little auk, at the southern edge of its breeding range. J Exp Mar Bio Ecol. 2016;475:89–99.
Gjøsæter H. Commercial fisheries (fish, seafood and marine mammals). In: Sakshaug E, Johnsen G, Kovacs K, editors. Ecosystem Barents sea. Trondheim: Tapir Academic Press; 2009. p. 373–414.
Furness RW, Hislop JRG. Diets and feeding ecology of great skuas Catharacta skua during the breeding season in Shetland. J Zool. 1981;195:1–23.
Phillips RA, Thompson DR, Hamer KC. The impact of great skua predation on seabird populations at St Kilda: a bioenergetics model. J Appl Ecol. 1999;36:218–32.
ICES. 31 ICES WGIBAR REPORT 2015. Annex 5: Barents Sea Ecosystem State 2014. 2015.
ICES. Report of the Working Group on the Integrated Assessments of the Barents Sea. WGIBAR 2017 Report 16-18 March 2017. Murmansk: ICES; 2017. ICES CM 2017/SSGIEA:04.
Mehl S, Aglen A, Alexandrov DI, Bogstad B, Dingsør GE, Gjøsæter H, et al. Fish investigations in the Barents Sea winter 2013-2014. IMR/PINRO joint Report series. 2014.
Votier SC, Bearhop S, Crane JE, Arcos JM, Furness W, Bearhop S, et al. Seabird predation by skuas Stercorarius skua - intra-specific competition for food? J Avian Biol. 2007;38:234–46.
Hamer KC, Furness RW. Age-specific breeding performance and reproductive effort in great Skuas Catharacta skua. J Anim Ecol. 1991;60:693–704.
Barrett R, Anker-Nilssen T, Bustnes J, Christensen-Dalsgaard S, Descamps S, Erikstad KE, et al. Key-site monitoring in Norway 2014, Including Svalbard and Jan Mayen. 2015. http://www.seapop.no/opencms/export/sites/SEAPOP/no/filer/short-reports/2015/SEAPOP-Short-Report-1-2015.pdf.
Mehlum F, Gabrielsen GW. The diet of high-arctic seabirds in coastal nad ice-covered pelagic areas near the Svalbard archipelago. Polar Res. 1993;12:1–20.
ICES. Final Report of the Working Group on the Integrated Assessments of the Barents Sea (WGIBAR), 22-26 February 201. Murmansk: ICES; 2016. ICES CM 2016/SSGIEA:04.
Furness RW, Crane JE, Bearhop S, Garthe S, Käkelä A, Käkelä R, et al. Techniques to link individual migration patterns of seabirds with diet specialization, condition and breeding performance. Ardea. 2006;94:631–8.
We thank Elínborg Sædís Pálsdóttir and Magdalena Hadwiczak for assistance in the field. We appreciate the improvements of our English as suggested by Rob Barrett.
The study was supported by funds from Poland through University of Gdańsk, and from Norway through the SEAPOP program (http://seapop.no/), and the Norwegian Polar Institute.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on request.
The study was conducted with permission of the Governor of Svalbard and Norwegian Animal Research Authority.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.