Foraging behaviour of black guillemots at three Norwegian sites during the breeding season

The intensification of coastal development poses potential threats for coastal seabirds, and understanding their habitat use is a key factor to guide conservation and management. In sub-arctic areas, black guillemots (Cepphus grylle) use coastal habitats year-round, which makes them vulnerable to the increasing human activities in these areas. In mainland Norway, one of the species’ strongholds, black guillemots are red-listed after substantial population declines. However, their fine-scale foraging behaviour has received little attention to date. We collected and analysed GPS tracking data from adult black guillemots at three sites located over a latitudinal gradient of 250 km in central/northern Norway. Maximum foraging ranges of 33 km at Sklinna (65°12′N) for incubating birds, and 18 km at both Vega (65°34′N) and Sklinna for chick-rearing birds, are among the longest reported for this species. At all three sites, foraging probability was highest in shallow waters (< 50 m depth) close to the colony and declined with increasing water depth and distance from colony. However, birds from Vega also foraged over deeper waters. Kelp presence was of high importance at Sklinna, but apparently less important at Røst (67°26’N) and Vega. We also found distinct differences in foraging activity across the day and with tidal height among the sites. Inter-site differences in habitat use and foraging activity may be explained by differences in the availability of habitats and suitable prey. Our study highlights the importance of shallow marine areas for black guillemots and shows that habitat use can vary substantially between sites.


Introduction
Seabirds are increasingly used as indicator species for ecosystem health and food availability (Piatt et al. 2008). This is because many of the species are highly sensitive to environmental changes, which also make them among the most threatened taxonomic groups of birds (Croxall et al. 2012).
Seabirds are facing a multitude of anthropogenic stressors including introduced predators, climate change, and bycatch (Dias et al. 2019). In addition to changes in the food web due to ocean warming, spread of invasive species, pollution and human exploitation of fish stocks (Sydeman et al. 2017), human activities in the coastal zone are increasing worldwide (Crain et al. 2008;Brown et al. 2018). This coastal development includes aquaculture (Quick et al. 2004), harvesting of kelp forests (Christensen-Dalsgaard et al. 2020) and wind power installations at sea (Benjamins et al. 2020;Peschko et al. 2020a;, as well as increasing ship and boat traffic (Fliessbach et al. 2019;Dehnhard et al. 2020b), all of which may cause disturbance and lead to displacement of seabirds.
The foraging habitats used by seabirds during the breeding season are particularly important and relevant to take into account for sustainable marine spatial planning and protection plans, since reduced prey availability around the breeding site can have severe consequences for breeding 1 3 87 Page 2 of 18 success Ponchon et al. 2014;Divoky et al. 2015). This is amplified by the fact that the spatial range in which seabirds may find food resources is limited particularly during the breeding season, when parental birds engage in central-place foraging to optimize time at their nest sites attending their eggs and chicks (Bell 1990). Especially in those species that carry single prey items back to the nest to feed their chicks ('single-prey loaders', e.g. guillemots and terns), prey availability in close proximity to the breeding site during the chick-rearing period is of high importance, such that chicks receive a sufficiently high number of suitably sized and energy-rich prey items throughout each day (Monaghan et al. 1994;McLeay et al. 2009;Gaglio et al. 2018).
Black guillemots (Cepphus grylle) breed in the temperate and arctic zones of the Atlantic and Pacific Ocean and associated seas. The species is a single-prey loader that raises up to two chicks (Cairns 1987b;Divoky et al. 2021). In the Arctic, the species is often associated with sea ice (Divoky et al. 2016). In sub-arctic and temperate seas, however, black guillemots stay within the coastal zone (here including the continental shelf) and usually within 100 km of the nearest land all year-round and typically undertake rather short migrations (< 500 km) during the non-breeding season (Ewins 1988;Ewins and Kirk 1988;Bakken et al. 2003;Baak et al. 2021). They are therefore particularly sensitive to regional or local changes in coastal habitats. Black guillemots are among the seabird species most vulnerable to bycatch in coastal gillnet fisheries (Fangel et al. 2015;Baerum et al. 2019;Christensen-Dalsgaard et al. 2019), and the introduced American mink (Neogale vison) affects breeding success and adult survival of black guillemots in many colonies in Norway (Nordström et al. 2003;Fjeld et al. 2020). Increasing boat traffic and the associated disturbance may also be a problem for many coastal species (e.g. Dehnhard et al. 2020b) including black guillemots (Ronconi and St. Clair 2002). In addition, climate change can lead to changes in their food base and is expected to negatively affect population trends at the southern boundary of the species' distribution range in the Baltic Sea (Buchadas and Hof 2017). Climate change is also considered a severe threat to kelp forest ecosystems (Smale 2020;Smith et al. 2023), which represent important foraging areas for black guillemots in Scotland (Owen 2015;Johnston et al. 2021). This also means that kelp harvesting poses a potential threat for black guillemots. Large-scale experiments have shown that kelp harvesting affects fish assemblages and is thought to decrease ecosystem resilience . Assessments of impacts of kelp harvesting on seabirds are scarce, and so far limited to European shags (Gulosus aristotelis) (Christensen-Dalsgaard et al. 2020) and great cormorants (Phalacrocorax carbo) (Lorentsen et al. 2010b), while the industry in Norway is expanding northwards (Fiskeridirektoratet 2022).
Black guillemots are coastal foragers with typically shallow-and often benthic-diving patterns in the range of about 15-20 m, although they can dive down to a depth of at least 43 m (Masden et al. 2013;Shoji et al. 2015;Divoky et al. 2021). Diving activity of black guillemots in Northern Ireland occurred only during daylight hours (Shoji et al. 2015). Studies from Iceland, Scotland and Northern Ireland showed that the foraging range of chick-rearing individuals is typically less than 10 km, although trips up to 24 km have been observed (Petersen 1981;Sawyer 1999;Owen 2015;Shoji et al. 2015;Johnston et al. 2018Johnston et al. , 2021. In Hudson Bay, Canada, foraging distances up to 15 km were observed during the breeding season (both incubation and chick-rearing), including observations of non-breeders and failed breeders (Cairns 1987a). In Scotland, the species forages particularly in kelp forests and often in environments with strong tidal currents (Masden et al. 2013;Owen 2015;Waggitt et al. 2017;Johnston et al. 2021), while the distribution of sea ice determines the habitat use in the high Arctic (Cairns 1987a;Divoky et al. 2021).
Matching their coastal diving foraging behaviour, black guillemot diet has been characterised by a number of benthic and kelp-forest associated fish species, varying with latitude, including sculpins (Cottidae), butterfish (Pholis gunnellus), gadoids (Gadiformes), sandeels (Ammodytidae) and squat lobsters (sub-order Anamura) in more temperate climates, and Arctic cod (Boreogadus saida) in the high Arctic (Slater and Slater 1972;Ewins 1990;Lønne and Gabrielsen 1992;Barrett and Anker-Nilssen 1997;Byers et al. 2010;Barrett et al. 2016; this study-see Supplement 1). Most diet data originate from observations of chick-feeding adults at the colonies, and black guillemots typically show a high chickfeeding rate in the early morning and during evening hours (Ewins 1990;Shoji et al. 2015).
The overarching goal of this study was to gain a better knowledge about the foraging behaviour of adult black guillemots in Norway during the breeding season. At least 10% of the global population of black guillemots breed in mainland Norway, with approximately 35,000 breeding pairs in 2005 (Barrett et al. 2006). However, the population is declining in many areas and black guillemots are now listed as "near threatened" in mainland Norway, with an estimated overall decline in the range of 15-30% over the last 30 years (Stokke et al. 2021). GPS tracking or seabird-at-sea observation data that could inform about the habitat use of black guillemots have essentially been lacking for the whole of mainland Norway. A better understanding of the foraging behaviour is much needed for conservation purposes and marine spatial planning, concerning especially the threats from bycatch and the growing kelp-harvesting industry in Norway.
Our study specifically aimed to document the fine-scale spatial habitat use by black guillemots and the environmental characteristics associated with these areas, and assess potential differences in foraging activity and habitat preferences between three different populations spread around the Arctic circle in northern and central Norway. We tested for variation in spatial habitat use, trip duration and foraging distance between the sexes and among sites, as well as for variation in environmental habitat preferences among the three populations. In addition, we examined how black guillemot foraging activity changed across the tidal cycle, as this may influence energetic costs of benthic diving, but also prey distribution (Wirjoatmodjo and Pitcher 1984;Hall et al. 1996;Enstipp et al. 2006) and the birds' diurnal cycle (Ewins 1990;Shoji et al. 2015). With almost permanent daylight at our study sites during summer, and twilight prevailing throughout the darkest hours of the night, we predicted that black guillemots would extend foraging into night-time hours.

Methods
Fieldwork was conducted on Heimøya, Sklinna (65°12′N, 10°59′E) in 2013, 2018 and 2019, on Hernyken, Røst (67°25′N, 11°53′E) in 2019 and 2021, and on Muddvaer, Vega (65°34′N, 11°41′E), in 2021. The population at Heimøya was stable with about 450 individual black guillemots during the study period, while the total number of birds in the entire Sklinna archipelago is likely somewhat higher. The population at Muddvaer held about 100 individuals in 2021 and has a stable population trend since 2013, whereas Røst has more than 1500 pairs (Anker-Nilssen 2009), with 40-50 pairs at Hernyken and nearby islets, a population that has increased by 5% per year over the last decade (Anker-Nilssen unpubl. data).
Both Sklinna and Røst are offshore archipelagos that are surrounded by deeper waters (Fig. 1). Sklinna measures about 4 by 4 km in size, with several islands and skerries (i.e. rocky islets). A few small skerries are located about 8 km south of Sklinna, while the nearest other land is 20 km to the south (outer islands of the larger Vikna archipelago, closer to the mainland). Twenty kilometres east of Sklinna is another extended shallow marine area with numerous skerries (Hortavaer). Røst, by contrast, is about 20 by 10 km in size, with several hundred islands and skerries. From Hernyken, the nearest land outside Røst (Vaerøy) is found 37 km to the northeast (including 20 km of open sea). Vega is a large island separated from the mainland of Norway by 12 km of open water. Muddvaer forms a small (3 by 2 km in size) archipelago, which is located 4 km south of Vega, and there are numerous skerries surrounding Muddvaer.
Black guillemots were tracked with 18 g ® Mataki GPS loggers (Debug Innovations, Cambridge, UK) in 2013 (Sklinna only) and 5-11 g ® PathTrack nanofix GPS loggers (Pathtrack Ltd, Otley, UK) in 2018, 2019 and 2021 (see Table 1). Both types of loggers remotely downloaded the recorded GPS data (at 5-min intervals for Mataki loggers, and either 30-min or 60-min intervals for PathTrack loggers) to one or two base stations placed in the colony. Mataki loggers were programmed to record data every 5 min, Path-Track loggers were programmed to record data at 4-min intervals at Sklinna and Muddvaer, and 1-or 5-min intervals at Røst. Both logger types were attached with 1-2 g ® Tesa tape on the lower back of the birds and thus added mass represented between 1.6 and 5.5% of the birds' body mass (Table 1). Birds were caught either by hand or with a crook on the nest, by mist net when approaching or leaving the colony, or with a leg-hold noose mat at a resting location at the edge of the colony. Birds were measured and weighed during handling, and for those that were not part of earlier studies, a small blood sample was taken from the foot web for molecular sexing. Handling time at deployment did not exceed 15 min. We attempted a second capture to recover loggers only in those cases where the nest site was known, but in most such cases, the logger had already fallen off.
Not all birds could be assigned to a nest, and the nest content was not known for all nests at Sklinna and Muddvaer, while all birds at Røst were chick-rearing (see Table 1). We therefore cannot rule out that some of the birds that were tracked were non-breeders or had failed their breeding attempt. Loggers collected on average data for 2 days (range 1-7). We did not conduct systematic nest checks after deployment. Due to the short tracking durations per individual, we assumed that those birds that were incubating were doing so during the entire duration during which GPS loggers were attached and collected data. At Røst and Muddvaer, chick-rearing individuals that were deployed with GPS loggers were observed to continue to feed their chicks. There were no indications that nests failed either due to handling of the birds or from other causes, while one of their parents was carrying a GPS logger. Since our study design was not balanced with respect to breeding stages (Table 1), we performed all statistical analyses both for the entire dataset (i.e. pooling all available data) and for chick-rearing individuals only (i.e. subsampling the data from Sklinna and Muddvaer).

Identification of foraging locations
We identified the most likely foraging locations based on expectation-maximization binary clustering (EMbC) (Garriga et al. 2016b). EMbC uses velocity and turning angle to classify movement data into four different clusters aligned with likely behavioural states: low velocities and low turns (LL; "resting"), low velocities and high turns (LH; "intensive search"), high velocities and low turns (HL; "travelling") and high velocities and high turns (HH; "extensive search") (Garriga et al. 2016b). EMbC has been shown to be useful across a broad range of species (e.g. Cecere et al. 2020;De Pascalis et al. 2020;Dehnhard et al. 2020a; and comes with the advantage of requiring less supervision, less a-priori assumptions and less computational power than other approaches (Garriga et al. 2016b).
GPS data were analysed in the EMbC R-package (Garriga et al. 2016a), using the stack clustering function stbc, which accounts for potential among-individual behavioural differences by annotating behavioural states for each individual. GPS data were not interpolated, and thus, sampling intervals differed within and between individuals. EMbC has proven to be robust towards variation in sampling intervals (Garriga et al. 2016b). The same approach (i.e. stack-clustering of GPS data with different sampling intervals) worked well for European shags, where sampling intervals differed substantially more than in our current dataset .
The stack clustering function was run on data that were pooled for all years and sites, in order to guarantee the same cut-offs among years and sites. The pre-and postsmoothing options were set to zero. Previous studies on European shags  and Antarctic fulmarine petrels (Dehnhard et al. 2020a) found a good match between EMbC stages LL and LH with simultaneously collected diving and wet-dry data, respectively. Lacking simultaneously collected diving data in this current study, we here adopted the same approach and defined all GPS time stamps with EMbC states LL and LH as foraging locations in further analyses. Fig. 1 Location of the three study sites (panel a, white stars) and GPS positions recorded at the different sites (white dots). Black guillemots were tracked during the chick-rearing period, but the datasets from Muddvaer and Sklinna included individuals with unknown breeding status (and thus potential non-breeders). Black dots mark foraging locations as identified through EMbC (see Methods). Except for panel a, land areas are shown in grey, depths in shades of blue, with grey lines marking the 50 m isobaths. Green areas indicate kelp forests. Note that the scale differs between the plots We excluded locations within 300 m of the centre of all colonies, since these locations are strongly associated with resting, preening and social display activities, besides some foraging that happens close to the colony. Distinguishing between foraging and non-foraging activities in this area is thus very challenging and would require the use of dive loggers in addition to GPS loggers. However, the proximity to the colony in itself means that this area is of high importance for the birds, while we lack information about the importance of areas further away. Similarly, any roosting places at or near foraging sites were removed, i.e. when GPS time stamps were located on land and not at sea. One apparent locational outlier was detected and removed.
When comparing foraging trip metrics between populations (Table 2), we excluded incomplete foraging trips. We defined foraging trips as movement paths ≥ 5 min ≥ 300 m away from the colony. Trips were considered incomplete when (1) locations at the colony were not available either immediately before or after the trip, (2) gaps of > 10 min existed between the last and/or first location in the trip and the next or previous location at the colony, and (3) gaps of > 30 min existed between locations during the trip. Due to the inability of GPS loggers to acquire locations when submerged (i.e. when the bird is diving), trips recorded by both logger types frequently included gaps between GPS locations. The 30-min cut-off to define incomplete trips Table 1 Overview of successful GPS deployments on black guillemots, including type, mass, settings and deployment period (start-end dates) for the loggers and information on sex, breeding stage, body mass (mean ± SE, range) and added mass in % of body mass When calculating added mass, the maximum weight of tape (2 g) was included. Information in parentheses refers to the number of individuals (for GPS interval was chosen as a conservative measure based on average trip duration (see Table 2) and the obtained GPS intervals. Since black guillemots breed in sheltered cavities under rocks, short visits to the nests to feed the chicks may be difficult to detect with the chosen settings of the GPS loggers and the above-named cut-offs. To make sure trips were identified correctly, we therefore also visually inspected distance from colony over time and could confirm that the abovelisted thresholds also correctly identified such trips where parental birds only returned very briefly to their nests to feed the chicks.

Environmental variables
To identify marine habitat preferences, we selected three environmental variables that characterise the coastal habitat black guillemots typically exploit (Johnston et al. 2018): sea depth, seabed slope and the presence/absence of kelp. We downloaded bathymetry data from GEBCO (https:// www. gebco. net/ data_ and_ produ cts/ gridd ed_ bathy metry_ data/), which comes as a spatial grid with a resolution of 15 × 15 arc seconds. Slope (in degrees) was calculated from the bathymetry data as the maximum change from the location grid cell to its eight closest neighbours using the R-package terra (Hijmans 2021). Kelp data as modelled by the Norwegian Institute for Water Research (NIVA) based on depth, wave activity, salinity and several other factors at 10 m spatial resolution (Bekkby et al. 2009;2020) were obtained from The Norwegian Environment Agency (https:// geoco rtex01. miljo direk torat et. no/ Html5 Viewe r/? viewer= natur base). Sea bottom substrate data from the European Marine Observation and Data Network (EMODnet) (http:// gis. ices. dk/ geone twork/ srv/ eng/ catal og. searc h#/ metad ata/ 01bf1 f24-fdcd-4ee7-af8b-e62cf 72fe2 f9) were considered, but since data coverage in the coastal areas was poor, we refrained from including this variable in our analyses.
To assess whether foraging activity differed with tidal height, we downloaded hourly water level observations for the three sites from the website of the Norwegian Mapping Authority (https:// www. kartv erket. no/ en/ at-sea/ se-havni va). We accessed daily data for time of sunrise and sunset at the breeding sites from www. timea nddate. no. These were not included in statistical analyses but used for visualisation.

Statistics
All statistical procedures were carried out in R (version R-4.2.2; Development Core Team 2022). Due to the low number of complete foraging trips in some of the years, we pooled all available data per site and refrained from investigating inter-annual effects. We investigated differences in foraging trip metrics between sites, breeding stages (where known) and sexes, using linear mixed models (LMMs) in the R-package lme4 (Bates et al. 2015). P-values were calculated with lmerTest (Kuznetsova et al. 2014). Models were run with either maximum distance from colony or trip duration as dependent variable, and bird ID as random factor. We ran two LMMs with site, breeding stage and sex (no interaction terms) as explanatory variables. Breeding stage could influence both maximum distance from colony and trip duration. Since chick-rearing was the only breeding stage during which birds from all three sites were tracked, we additionally ran two LMMs as described above using only data from this stage to check whether breeding stage had a strong influence on model results and thus inter-site differences.
To investigate marine habitat preferences, we used resource-selection functions (Boyce and McDonald 1999;Signer et al. 2019;Northrup et al. 2022), i.e. we compared covariates associated with the foraging locations of the birds with random locations within the spatially available habitat. We defined available habitat as the area within reach for birds around their colonies, thus creating a circular buffer around each colony. As in other studies comparing the habitat preferences of birds originating from different colonies (Christensen-Dalsgaard et al. 2017;Péron et al. 2018;Dehnhard et al. 2022), the radius was set as the maximum distance between a registered foraging location and the colony, separately per site (cf. Table 2). Since there was evidence for breeding stage to have an effect on maximum distance from colony, and particularly chick-rearing birds to travel shorter distances than incubating birds (see Results), we applied a different radius for those birds from Sklinna that were known to be chick-rearing. The radius was thus largest for incubating birds and those with unknown breeding status from Sklinna (33 km). Chick-rearing birds from Sklinna had a radius of 18 km, and thus the same as chick-rearing birds and those with unknown breeding status from Muddvaer (both 18 km). With 11 km, birds from Røst (all chick-rearing) had the smallest radius. To create a representative sample of the available habitats within these areas, five point locations were created randomly for each foraging location, separately per individual, within the defined available area, using the R-package sf (Pebesma 2018). Land areas within the circular buffers were removed before generating random locations. All foraging locations and random locations were intersected with the environmental layers using the R-packages sf and terra.
We then ran generalised additive mixed models (GAMMs) with a binomial distribution (1 = foraging locations, 0 = available, i.e. random locations). Like generalised additive models, GAMMs allow the fitting of nonlinear responses to predictor variables, which is a major advantage as animals rarely respond linearly to their environment (Aarts et al. 2008;Dehnhard et al. 2020a). The GAMMs were run using the R-package mgcv (version 1.8-38; Wood 2017) with a logit link function. Similarly to Christensen-Dalsgaard et al. (2017), models were fitted using thin-plate regression smoothing (Wood 2017). We initially set the maximum number of knots for smooth terms to 5 in order to avoid overfitting and used the gam.check function to check whether models with less knots had a better fit. GAMMs were run separately for each site. For Muddvaer and Sklinna, where the breeding stage was unknown for some individuals, we ran models twice: once with the full dataset and once only for chickrearing birds. We followed a forward-stepwise approach to add environmental covariates and modelled the environmental habitat preferences separately for each site. The initial models therefore contained one environmental covariate, either as fixed factor (kelp) or as smooth term (depth, slope, distance from colony) as well as Bird ID as random effect. After identifying the best-performing environmental variable, we assessed the inclusion of a second and third environmental variable. To avoid collinearity, we only included environmental covariates in the same model which had a mutual Spearman's rank correlation coefficient of ≤ 0.5 (for the Muddvaer dataset depth and slope were correlated, and for the Røst dataset depth and distance from colony were correlated). Model selection was based on the Akaike information criterion (AIC), and we did not attempt to fit more than three terms into the final model to avoid over-fitting.
Finally, we analysed whether time of day and waterlevel changes caused by tidal activity had any impact on foraging activity. For this purpose, we used the entire dataset, i.e. including also locations within 300 m of the centre of the colony. Foraging activity (as binomial response variable) was defined as time stamps associated with foraging locations (coded as 1, based on EMbC, as described above); all other time stamps (including all locations within 300 m of the colony) were defined as nonforaging activity (coded as 0). As for the other analyses, we pooled data from all years per site, and once more, for Muddvaer and Sklinna, where the breeding stage was unknown for some individuals, we ran models twice: once with the full dataset and once only for chick-rearing birds. The main GAMM thus contained foraging activity as binary response variable, and both hour of the day and water level as explanatory variables. We used thin-plate regression smooth terms for water level, and cyclic cubic regression splines for hour of the day (since hour of the day is circular). For the number of knots, we proceeded as described above. Bird ID was again included as a random effect.
Instead of referring to p-values as cut-offs to significance testing, we used the language of evidence throughout this paper, thus referring to strong or weak evidence for/against differences between groups (Muff et al. 2021).

Foraging trip metrics between sites, breeding stages and sexes
Using all available data from complete foraging trips, including those of potentially non-or failed breeders, there was no evidence that maximum distance from colony or trip duration in black guillemots differed among sites (LMMs, all F 2 ≤ 1.87, p ≥ 0.209) or sexes (LMMs, all F 1 ≤ 0.76, p ≥ 0.408). There was weak evidence that both maximum distance from colony and trip duration were longer during incubation compared to chick-rearing (LMMs, all F 2 ≥ 3.65, p ≤ 0.080), matching with the longest foraging trip being conducted by an incubating bird from Sklinna (33 km). During chick-rearing, birds foraged within 18 km around Vega and Sklinna, and 11 km around Røst, but typically closer to the colonies as reflected by the average trip distances (Table 2; Supplement 2). Using a subset of only chick-rearing birds, we found no evidence for maximum distance from colony or trip duration to differ among sites (LMMs; all F 2 ≤ 2.14, p ≥ 0.136) or between sexes (F 1 ≤ 0.55, p ≥ 0.472).

Marine habitat preferences
Distance from colony and depth were the two most important environmental variable to explain foraging habitat, whether the full dataset was used, or only data from chickrearing individuals. Distance from colony explained ≥ 43% of deviance in the models, and distance from colony was the single-most important variable at both Røst and Muddvaer (Table 3). At Sklinna, depth was the best predictor of foraging habitat (51% of deviance explained for all birds, and 68% of deviance explained for chick-rearing individuals), while depth was of secondary importance at Røst (27% of deviance explained) and Muddvaer (≤ 9% of deviance explained). Kelp presence was of high importance at Sklinna but was the least supported variable in models for both Røst and Muddvaer (  (Table 3). Since depth was correlated with distance from colony at Røst, the final composite model for Røst did not include depth, but slope and kelp. At Muddvaer, slope was correlated with depth, and therefore, the final composite model for Muddvaer did not include slope but kelp. Final GAMMs thus had a different structure for each of the three sites, reflecting somewhat different habitat preferences among the three sites (Table 3). For Muddvaer, model structure was identical for the full dataset and chick-rearing individuals. For Sklinna, distance from colony was of higher importance in the models for chick-rearing individuals compared to the full dataset (Table 3).
GAMM response curves for depth, distance from colony and slope showed some variation among sites for birds independent of their breeding status ( Fig. 2; Figure S3.2). While black guillemots at Røst and Sklinna showed a clear preference to forage in shallow waters, with foraging probability declining steeply with depth, this pattern was less pronounced for birds from Muddvaer, which showed a peak in foraging probability at a water depth of about 120 m (Fig. 2). Foraging probability decreased with distance from colony for birds from all three sites, but in all cases there was a second peak in foraging probability further away, the location of which differed between sites (~ 7 km for Røst, ~ 12 km for Muddvaer and ~ 16 km for Sklinna; Fig. 2). Finally, foraging probability in response to seabed slope varied strongly among sites. Birds from Røst foraged over areas that were essentially flat (< 2°), and seabed slope in the available habitat around the colony was never more pronounced than 5°, limiting the importance of this factor. The available habitat for birds from Sklinna and Muddvaer showed a larger variation, with slopes down to almost 20°, but foraging activity at both sites was limited to flat or moderately sloped seabed areas (Fig. 2). The response curves of models run on the full dataset for Sklinna and Muddvaer, i.e. including birds with unknown breeding status (both sites) and incubating birds (Sklinna only), were similar to those of the chick-rearing individuals ( Figure S3.2) with one exception: The second peak in foraging probability for birds from Sklinna was further away (~ 12-20 km) than for chick-rearing individuals, and there was another small increase in foraging probability at 30 km distance from colony ( Figure S3.2).

Effects of time of day and water-level changes
There was strong-to-very-strong evidence that water level and hour of the day affected foraging activity of black guillemots (GAMM all chi-square ≥ 15.71, all p ≤ 0.001; Fig. 3). Model results and model responses were very similar when considering only chick-rearing birds (Figure S3.3). The nature of the responses to water level and hour of the day differed, however, markedly among sites. At Røst and Muddvaer, foraging activity over the course of the day showed very pronounced but inverse patterns. At Muddvaer, foraging activity increased during the afternoon and was highest between 20:00 and 05:00 local time, and thus during several hours of civil twilight, whereas Røst birds were relatively inactive during these hours and foraged mostly during the middle of the day. At Sklinna, foraging activity showed a less pronounced pattern throughout the day than at the other sites, but was somewhat higher from the afternoon and throughout the For each site, the model structure of the six best-fitting models, including the best additive models (_best, _2nd) is given. For Sklinna and Muddvaer, models were run once with the full dataset, and once with a subset including only chick-rearing birds to improve comparability among colonies. Models are ranked by AIC, and we further present the adjusted R 2 values and the explained deviance. distcol distance from colony. Bird ID was included as random factor in all models *For the Muddvaer dataset, depth and slope were correlated, and for the Røst dataset depth and distcol were correlated. Correlated variables were not fitted in the same models to avoid collinearity issues  (Fig. 3).

Foraging range and trip duration
Foraging range and trip duration did not differ statistically among our three study sites. Irrespective of breeding status, most of the foraging activity took place close to the colonies (cf. Figure 1, Figure S3.1) and thus within the typical foraging range of < 10 km around the colonies during chick-rearing reported elsewhere in Europe (Petersen 1981;Sawyer 1999;Owen 2015;Shoji et al. 2015;Johnston et al. 2018;Owen et al. 2019). However, in contrast to these previous studies and the results from Røst, chick-rearing birds from Sklinna and Muddvaer travelled up to 18 km away from their colonies to find food. We also found evidence for breeding stage to affect foraging ranges, with chick-rearing birds travelling shorter distances than incubating birds. This has also been observed in Brünnich's guillemots (Uria lomvia; Ito et al. 2010), yet it does not seem to be a universal pattern in alcids, since it was not found in common guillemots (Uria aalge; Gulka and Davoren 2019) and little auks (Alle alle; Jakubas et al. 2014). A recent study from Caithness, Scotland, also found black guillemots to travel further (up to 26 km) from their colony during the incubation period than during chick-rearing (up to 12 km; Johnston et al. 2021). Maximum foraging ranges of black guillemots from Sklinna (33 km during incubation and 18 km during chick-rearing, respectively) were even longer and are among the longest reported so far for the species during the breeding period. A number of factors can explain variation in foraging range among breeding sites of black guillemots. Local distribution and patchiness of suitable foraging habitat as well as prey availability seem to be the most important factors (addressed below). In addition, breeding status is clearly important, possibly also chick age and whether parents are provisioning one or two chicks: Being able to carry only one fish at a time in the beak while raising two chicks requires a higher feeding frequency of both parents-and may force them to forage closer to the breeding site (Cairns 1987b;Divoky et al. 2021). Unfortunately, we only had exact data on chick age and brood size for 5 GPS birds at Røst and could not investigate this relationship any further.
GPS loggers in themselves may also affect the behaviour of the tracked birds, due to the added weight and drag Fig. 3 Response curves from the GAMM testing the effect of hour of the day (left) and water level (right) on foraging activity of all individuals (also those with unknown breeding status) at the three sites. For Røst and Sklinna, data from all years were pooled. The model included the smoothed effects for both water level and hour of the day, as well as the interaction of both with site. Chi-square and p-values refer to the summary statistics of the GAMM. Grey areas for the plots on the left-hand side indicate the maximum twilight hours at the sites when birds were tracked. Hour of the day refers to local time (i.e. GMT + 2 h) (Barron et al. 2010;Vandenabeele et al. 2011Vandenabeele et al. , 2012Evans et al. 2020). Black-legged kittiwakes (Rissa tridactyla) equipped with GPS-loggers performed longer-lasting trips compared to control birds (Heggøy et al. 2015). The same pattern was found in a meta-analysis studying potential impacts of biologging (Bodey et al. 2018). In our study, we used logger types of different sizes, and the added mass (including tape) ranged between 1.5% (Røst 2021) and 5.5% (Sklinna 2013) of the birds' body mass (Table 1), exceeding the commonly recommended threshold of 3% (Vandenabeele et al. 2012) at both Sklinna and Muddvaer. Altogether, we cannot rule out logger effects and would certainly recommend to always use the smallest loggers available (both in size and mass), and not to vary logger sizes across sites and years, when logistically and financially feasible. Finally, since black guillemots are very efficient in removing their loggers (in some cases within less than 24 h), we are confident that the logger attachment had no long-lasting negative impacts on the tracked birds.

Environmental habitat preferences
Matching previous studies on habitat use of black guillemots in the North Atlantic (reviewed in Johnston et al. 2018), we found that shallow (< 50 m) marine habitat near the colony was of high importance for this species also at the three Norwegian study sites. However, there were distinct differences in habitat preferences among sites. Firstly, distance from colony was more important at Røst and Muddvaer compared to Sklinna, where depth was the most important variable. Secondly, the presence of kelp forests was strongly associated with foraging habitat around Sklinna, but less so at Røst and Muddvaer. Thirdly, while models suggested foraging areas of all three populations were associated with shallow waters, birds from Muddvaer appeared to be foraging in areas with a depth of about 120 m. Finally, we found that environmental habitat preferences differed only slightly (and only at Sklinna) between the subset of chick-rearing birds compared to the full dataset including incubating birds and those with an unknown breeding status: The main difference was that distance from colony was more important than the presence of kelp forest for chick-rearing birds, but the final model response curves looked almost identical for both groups (cf. Figure 2 versus Figure S3.2). Our results match somewhat with those for little terns (Sternula albifrons), where failed breeders conducted longer-lasting foraging trips than successful breeders (Perrow et al. 2015). In lesser black-backed gulls (Larus fuscus), however, foraging trip metrics were similar between failed and successful breeders, but use of terrestrial versus marine habitats differed between the two groups (Camphuysen et al. 2015). Although we were unable to determine the breeding status of all the birds we tracked, our results reflect the overall importance of proximity to the breeding site, shallow marine habitat and kelp forests for black guillemots during the breeding period.
The differences in habitat preferences among sites may be explained by differences in the availability of shallow marine habitat, kelp forest distribution and differences in diet. Being a rather small archipelago (about 4 by 4 km) surrounded by deeper waters, Sklinna by itself does not hold much of the shallow marine area preferred by black guillemots (cf. Fig. 1). By contrast, the availability of suitable foraging habitats in the vicinity of the colonies at Røst (lying on a shallow plateau of about 20 by 10 km in size, also surrounded by deeper waters) and Muddvaer (surrounded by extensive shallow waters close to the Norwegian mainland) is much larger. These differences in habitat availability and the fact that birds from Sklinna often crossed deeper waters before reaching shallow foraging locations (see Fig. 1), likely explain why depth was a more important variable in the habitat models for Sklinna than distance from colony. Furthermore, the majority of benthic marine habitat of < 50 m depth around Sklinna is covered by kelp forests (cf. Fig. 1d). By contrast, kelp forests overall cover less of the shallow areas around Muddvaer and Røst (cf. Fig. 1b and 1c), which also have large patches of sandy bottom as well as beds of dead coralline algae, a substrate that kelp cannot attach to (Bekkby et al. 2009). In addition, not all rocky substrate may hold kelp forest in these areas either, since overgrazing by green sea urchins (Strongylocentrotus droebachiensis) has decimated kelp forests around Vega and other parts of Northern Norway, and kelp forests are only slowly recovering (Christie et al. 2019;Eikrem et al. 2019). Finally, the chick diet of black guillemots at Røst and Sklinna is generally similar, but butterfish is a much more common prey item at Røst (Barrett and Anker-Nilssen 1997;Lorentsen et al. 2010a). This was also the case in the years during which GPS-tracking took place (Fayet and Anker-Nilssen, unpublished data). However, butterfish is not exclusively associated with kelp forests (Koop and Gibson 1991;Shorty and Gannon 2013). This could explain why kelp forests were not an equally important foraging habitat at Røst as at Sklinna. It is further possible that the adults targeted predominately kelp forested areas close to their colony to find food for their chicks, but travelled further away and westwards (cf. Fig. 1b) for self-feeding. That adult diet differs from chick diet is not uncommon in alcids (Davoren and Burger 1999;Wilson 2004;Breed et al. 2009;Myksvoll et al. 2013), including black guillemots (Ewins 1990). Finally, black guillemots are known to show individual specialisation in both foraging area (Johnston 2019;Owen et al. 2019) and diet (Slater and Slater 1972). For example, some individuals may specialise to feed on butterfish, and others on sculpins. We did not record the detailed diet of those individuals equipped with GPS-loggers and thus do not know whether and how many of the individuals tracked specialised on particular diet items.
In contrast to Røst and Sklinna, chick diet at Muddvaer was dominated by lesser sandeel (Ammodytes marinus), followed by haddock (Melanogrammus aeglefinus), saithe and butterfish (Supplement 1). Sandeels are strongly associated with sandy bottoms, and due to their high energy-content, they are also considered to be excellent food for seabird chicks (Wanless et al. 2018). The importance of sandeel in the chick diet, and the association of sandeel with sandy substrates, may thus explain the lower importance of kelp forests in habitat models for Muddvaer. Sandeels show high annual variability in their presence (Frederiksen et al. 2005) and high sensitivity towards global warming (Lindegren et al. 2018). Unfortunately, fieldwork at Muddvaer took place only during one breeding season, and we therefore do not know whether sandeels are of equally high importance in other years and throughout the entire breeding season. Saithe and haddock were also important prey items at Muddvaer and may not have been associated exclusively with kelp forests: Based on their size, both haddock and saithe delivered to chicks at Muddvaer and Røst were 0-group fish (< 1 year old; Neuheimer et al. 2008;Hillersøy and Lorentsen 2012). Both of these fish species have a pelagic juvenile phase, before settling to the seabed (haddock) and into kelp forests (saithe), respectively (Olsen et al. 2009). As indicated by their silvery colour, at least some of the saithe were still in their pelagic life stage, thus contributing to an explanation for the lower importance of kelp forests as foraging habitat around Muddvaer and Røst.
The peak in foraging probability in areas with a depth of ~ 120 m may be explained by the prey selection of birds at Muddvaer. This is more than twice the maximum diving depth recorded for black guillemots (50 m;Cairns 1987a;Masden et al. 2013;Johnston et al. 2018), and it is thus highly unlikely they were performing benthic dives in such deep waters. Feeding on juvenile saithe and haddock in their pelagic life stages, however, may well happen over deeper waters, without the need for black guillemots to dive down to the bottom. Also, sandeels occur at water depths down to 120 m, although they prefer depths of 30-70 m (Wright et al. 1998;van der Kooij et al. 2008). In spring and summer, they do, however, exhibit a diurnal change between staying in the sand at night, and feeding in fast-swimming schools in the upper water column, often close to the surface, during daytime (Winslade 1974;van der Kooij et al. 2008). As the black guillemots' foraging activity over deeper waters (> 30 m, i.e. outside of kelp forests) occurred mostly between midnight and 07:00 local time (Supplement 4; see also below), it still remains unclear what prey were targeted there.
Highly relevant for conservation, the modelled habitat preference of black guillemots for foraging depth (Fig. 2) matches not only the shallow coastal areas, but looks almost identical to the variation in modelled bycatch rates of diving seabirds in response to fishing depth of coastal gillnetting in Norway (cf. Figure 3 in Baerum et al. 2019). In other words, their preference for shallow coastal areas makes black guillemots particularly vulnerable to bycatch. Nevertheless, Baerum et al. (2019) found other species to be even more frequently taken as bycatch in coastal gillnet fisheries.
Another concern is the importance of kelp forest in habitat models, specifically around Sklinna, and to a lesser extent also around Røst and Muddvaer. This raises the question to what extent kelp harvesting affects prey availability and, thereby, the breeding success of black guillemots. Kelp harvesting takes place in the same areas as those utilized by foraging black guillemots from Sklinna (c.f. maps in Christensen-Dalsgaard et al. 2020). In 2022, the Norwegian Fishery Directorate extended the harvesting areas northwards, with the northern border put immediately south of Muddvaer, without indicating if this will remain a longer-term limitation (Fiskeridirektoratet 2022). Two of the nine birds from Muddvaer tracked in this study entered areas that now have been trawled or will be in the coming years.
The issues of kelp harvesting and bycatch highlight the very common management dilemma that most often, even in the case of coastal species, only the seabirds' breeding sites are protected, but not the areas they rely on to find their food (Yorio 2009;Critchley et al. 2018). In the light of the many declining seabird populations in Norway (Fauchald et al. 2015;Anker-Nilssen et al. 2022), where almost twothirds of the seabird species are now on the national red list (Stokke et al. 2021), it is high time to implement effective conservation actions and protect important foraging habitats.

Differences in foraging activity with time of day and water level among populations
Rather surprisingly, we found that foraging activity in response to water level and time of day differed distinctly among the study populations. At Røst, foraging activity was highest at low water levels, while at Sklinna, it was highest at intermediate water levels, and lower (but also associated with higher model uncertainty) at very high and low water levels. Finally, birds at Muddvaer showed high foraging activity also at high water levels. Tidal range is on average between 1.6 and 1.7 m for all three sites and thus would by itself not lead to differential energetic costs for diving among sites. It therefore seems more likely that differences in habitat and prey availability explain the differential responses to water levels among the three sites.
Black guillemots are known to occur in habitats with strong tidal currents (Masden et al. 2013;Owen 2015;Waggitt et al. 2017). Tidal currents are strongest at intermediate water levels, which is when the foraging activity at Sklinna was highest. Possibly, certain fish species are easier to catch under strong current conditions, e.g. if the accessibility of lower trophic level prey is positively related to tidal currents and the fish change behaviour accordingly (e.g. Hall et al. 1996). Other benthic fish species may be more accessible during high or low tide (Wirjoatmodjo and Pitcher 1984). In addition, black guillemots from Muddvaer also foraged over deeper waters, where they may have fed on pelagic sandeels or young gadoids in the upper water column and therefore were less dependent on effects of tidal currents on prey availability in benthic/kelp forest habitats.
We found clear indications for black guillemots from both Sklinna and Muddvaer to forage in civil twilight conditions during the night, whereas birds from Røst avoided nighttime foraging. A previous study performed in Northern Ireland also reported that diving only occurred during daytime hours, and not between 22:00 and 03:00 GMT (Shoji et al. 2015). We cannot rule out that the EMbC misidentified resting on the water as foraging-testing this would require simultaneous GPS and diving data. Similarly, social displays, which often take place near the colony in the early morning hours, especially at high tides, may have been categorized as foraging. However, it seems unlikely that misidentified foraging behaviour would cause patterns to differ so strongly among populations. Moreover, our data confirm that birds from Sklinna and Muddvaer were regularly more than 5 km away from their colony at night-time hours, whereas birds from Røst stayed within a range of about 1 km of the colony at night (Supplement 4). In addition, night-time foraging of black guillemots from Muddvaer was identified in waters more than 75 m deep (Supplement 4). This suggests that black guillemots from Sklinna and Muddvaer were actively foraging at night, and certainly not engaging in social activities, which is only known to happen near their colonies (i.e. within the 300 m range where all behaviour was coded as "non-foraging"; see Methods). In contrast to the study by Shoji et al. (2015), our more northern study locations are not totally dark in summer. Civil twilight (i.e. when the sun is at an angle of 0-6° under the horizon) prevailed throughout all nights of our study, with exception for the last 5 nights of GPS tracking of one bird at Røst in first half of August 2021, when nautical twilight (sun angle 6-12° under the horizon) prevailed for up to 3 h, but it was still not completely dark. This late tracking event at Røst may thus have had a slight influence on the patterns we found. On the other hand, nocturnal foraging has been observed in both Brünnich's and common guillemots (Regular et al. 2011;Elliott and Gaston 2015;Dunn et al. 2020;Patterson et al. 2022). It thus seems very likely that black guillemots are able to forage under nautical twilight conditions, especially since the species also winters along the coastline of North Norway and along the ice edge in Svalbard, as well as in other high-Arctic environments (e.g. Divoky et al. 2016). Daylength during winter is short at these latitudes, and other diving seabirds wintering in the Arctic are also known to forage during civil and nautical twilight conditions (e.g. Moe et al. 2021).
Prey choice may also have had an influence on the feeding patterns. However, we do not know what type of prey the GPS birds were feeding on, nor if foraging at night was more associated with self-feeding than providing for chicks.
In the larger perspective, and independent of the birds' foraging strategies, the fact that black guillemots in northern and central Norway during the breeding season also forage at night-time may be important information for conservation actions. As discussed above, one of the main threats for black guillemots is bycatch in gillnet fisheries (Fangel et al. 2015;Baerum et al. 2019;Christensen-Dalsgaard et al. 2019). Diel restrictions on gillnet-fishing have been considered in some fisheries as a measure to reduce bycatch (Melvin et al. 1999). For example, a recent study from the UK suggested night-only gillnet fishing as a potential mitigation measure to reduce bycatch of razorbills (Alca torda), common guillemots and European shags (Cleasby et al. 2022). Considering our results, such a policy would likely have limited effect on black guillemot bycatch in central and northern parts of Norway.

Conclusions
Our study confirmed that shallow marine habitat within 10-30 km of the breeding site is the most important foraging habitat for black guillemots. The presence of kelp forests was also important, but not equally important at all sites. Differences among sites in foraging habitat preferences as well as diurnal foraging patterns, including also night-time foraging at Sklinna and Muddvaer, are likely associated with differences in habitat availability and bathymetry as well as different availability of local fish prey. Our study is the first to provide information about the fine-scale foraging behaviour of black guillemots in Norway, and thereby provides important knowledge to inform and improve conservation actions.
In 2015, the UN formulated the global goal to protect 30% of the marine environment, including coastal zones, by 2030 (Maestro et al. 2019). As of 2021, Norway had only protected 3.6% of its territorial waters (Statistics Norway 2022). Based on our results and parallel work on other species (e.g. Dehnhard et al. 2022), protection of shallow marine areas and kelp forests in the vicinity of important breeding sites of coastal, fish-eating seabirds appears to be a good strategy for conservation of marine biodiversity in a changing world with increasing anthropogenic pressures.