Seabird diving behaviour reveals the functional significance of shelf-sea fronts as foraging hotspots

Oceanic fronts are key habitats for a diverse range of marine predators, yet how they influence fine-scale foraging behaviour is poorly understood. Here, we investigated the dive behaviour of northern gannets Morus bassanus in relation to shelf-sea fronts. We GPS (global positioning system) tracked 53 breeding birds and examined the relationship between 1901 foraging dives (from time-depth recorders) and thermal fronts (identified via Earth Observation composite front mapping) in the Celtic Sea, Northeast Atlantic. We (i) used a habitat-use availability analysis to determine whether gannets preferentially dived at fronts, and (ii) compared dive characteristics in relation to fronts to investigate the functional significance of these oceanographic features. We found that relationships between gannet dive probabilities and fronts varied by frontal metric and sex. While both sexes were more likely to dive in the presence of seasonally persistent fronts, links to more ephemeral features were less clear. Here, males were positively correlated with distance to front and cross-front gradient strength, with the reverse for females. Both sexes performed two dive strategies: shallow V-shaped plunge dives with little or no active swim phase (92% of dives) and deeper U-shaped dives with an active pursuit phase of at least 3 s (8% of dives). When foraging around fronts, gannets were half as likely to engage in U-shaped dives compared with V-shaped dives, independent of sex. Moreover, V-shaped dive durations were significantly shortened around fronts. These behavioural responses support the assertion that fronts are important foraging habitats for marine predators, and suggest a possible mechanistic link between the two in terms of dive behaviour. This research also emphasizes the importance of cross-disciplinary research when attempting to understand marine ecosystems.

SLC, 0000-0002-5409-4458; KLS, 0000-0003-0843-0956; SCV, 0000-0002-0976-0167 Oceanic fronts are key habitats for a diverse range of marine predators, yet how they influence fine-scale foraging behaviour is poorly understood. Here, we investigated the dive behaviour of northern gannets Morus bassanus in relation to shelf-sea fronts. We GPS (global positioning system) tracked 53 breeding birds and examined the relationship between 1901 foraging dives (from time-depth recorders) and thermal fronts (identified via Earth Observation composite front mapping) in the Celtic Sea, Northeast Atlantic. We (i) used a habitat-use availability analysis to determine whether gannets preferentially dived at fronts, and (ii) compared dive characteristics in relation to fronts to investigate the functional significance of these oceanographic features. We found that relationships between gannet dive probabilities and fronts varied by frontal metric and sex. While both sexes were more likely to dive in the presence of seasonally persistent fronts, links to more ephemeral features were less clear. Here, males were positively correlated with distance to front and crossfront gradient strength, with the reverse for females. Both sexes performed two dive strategies: shallow V-shaped plunge dives with little or no active swim phase (92% of dives) and deeper U-shaped dives with an active pursuit phase of at least 3 s (8% of dives). When foraging around fronts, gannets were sprattus, and a range of demersal fishes, scavenged from fishing boats [39][40][41]. Previous work has shown that both in the Celtic Sea and the Benguela upwelling region off western South Africa, gannets increase foraging effort (estimated via two-dimensional movement data) within regions where frontal activity is increased [35,42]. By examining their diving behaviour, we attempt to establish why.

Device deployment
Fieldwork was conducted on Grassholm, Wales, UK (51°43 N, 5°28 W; figure 1), during the breeding season in July 2012 and 2013. Chick-rearing gannets were caught on the nest during changeover (to ensure chicks were not left unattended, and so deployments began immediately with a foraging trip), using a brass crook attached to the end of an approximately 5 m carbon fibre pole. Birds were selected opportunistically, away from the edge of the colony (under licence from Natural Resources Wales). Upon capture, birds were weighed (to the nearest 50 g) and 1-2 ml of blood taken via the tarsal vein (under licence from the UK Home Office), a small aliquot of which was later used for molecular sexing (commercially outsourced to AvianBiotech.com). Birds were then equipped, using Tesa ® tape, with (i) a 30 g GPS logger (i-gotU GT-120, Mobile Active Technology Inc.) attached to either the dorsal surface of the central pair of tail feathers (2012) or the central back feathers (2013) and (ii) a 5.7 g or 10.5 g TDR (CEFAS G5 or LOTEK LAT 1810, respectively) attached to the ventral surface of the central pair of tail feathers (under licence from the British Trust for Ornithology). Total handling time was around 12 min. The maximum combined weight of deployed loggers (40.5 g) was 1.37% the average bird body weight (2948.8 g ± 33.0 g). Deployment durations ranged from 1 to 7 days.

Particulars of device data acquisition
The GPS recorder logged location at 1 min intervals with an accuracy of ±4.4 m [43]. The G5 TDR logged pressure and temperature every 0.1 s (10 Hz) during dives, defined as wet periods (detected via a wet/dry sensor) below a depth of 1.5 m. Pressure resolution was 4 cm of water column with an accuracy of ±1 m. The LOTEK LAT 1810 TDR logged temperature and pressure continuously at 1 s intervals (1 Hz). Pressure resolution was 2.5 cm of water column with an accuracy of ±1 m.

Animal behaviour metrics
2.3.1. Global positioning system tracks GPS fixes at night (between the end of civil dusk to the beginning of civil dawn) were stripped from the dataset to eliminate periods when birds rest on the water [44,45]. In addition, all activity within 2 km of the breeding colony was removed to account for bathing and rafting [46]. Tracks were then split into individual foraging trips. In some instances, the device deployment period exceeded the battery life of the GPS logger. Resultant incomplete foraging trips were included in analyses, unless otherwise stated.

Dive events
A bespoke algorithm, written in Matlab, was used to identify dive events by birds equipped with a LOTEK TDR. To be consistent with the technicalities of the CEFAS logger, dives were defined as periods where the registered depth was at least 1.5 m. This also accounted for shallow subsurface activity associated with non-foraging behaviours such as sitting on the water and bathing [38,47].
To allocate a location for each dive, GPS tracks were first interpolated to a resolution of 1 s (in time), using a cubic spline interpolation. The timestamp at the beginning of a dive event was then used to assign a corresponding location from the high-resolution interpolation. Allocations mismatched by more than a second (owing to the GPS logger battery life being surpassed or because a dive occurred within the colony exclusion zone) were excluded from analyses. A small number of shallow dives (five in total) occurred approximately 40 min after the end of civil dusk. These were considered atypical and excluded from analyses [44,45]

Dive characteristics
Gannets predominantly employ two dive strategies (figure 1a,b). Short and typically shallow plunge dives (termed V-shaped; figure 1a) involve little or no active swim phase [38,48], and may be used to catch fast, responsive prey in the upper section of the water column [49][50][51]. Indeed, this strategy is often associated with the predation of near-surface pelagic fish such as herring and mackerel [49,50]. In contrast, longer deeper pursuit dives (U-shaped; figure 1b) involve a substantial active swim phase that likely incurs higher energetic costs [37,52]. As such, this foraging strategy may only be used following a failed plunge attempt [37,48], or when prey are distributed at deeper depths (e.g. shoals of capelin Mallotus villosus or sand eels Ammodytes sp. at depth and/or near the seabed; [49,51]). In some instances, these dives may also be used to catch multiple prey items [51], although this is likely to be dependent upon the responsiveness and swimming speeds of prey [49,50].  influence dive strategy (which may reflect changes in prey accessibility and catchability), each dive was classified as either U-or V-shaped (figure 1a,b), as determined by the length of the active swim phase (Ushaped dives were those with an active swim phase of at least 3 or at least 4 s dependent upon the logger sampling rate; 1 s and 10 Hz, respectively; [38]). The start and end of the active swim phase was defined, using gradients in the vertical change in depth (electronic supplementary material, S1 and figure S1; [37,53]). In addition, for each dive, the maximum dive depth (metres) and total dive duration (seconds) was calculated.

Front metrics
Fronts vary in their strength, persistence and predictability [8]. To be able to determine the relative importance of these characteristics on foraging behaviour, thermal front activity across the Celtic Sea was summarized into three metrics that reflected this variability: (i) cross-front gradient strength Gdens, (ii) distance to closest front Fdist, and (iii) seasonal front frequency Ffreq. These were produced over two temporal scales: (i) 7 day composites [31] and (ii) seasonal composites [32], details of which are provided below (see §2.4.1 Seven day composites and §2.4.2 Seasonal composites). As the occurrences of thermal and chlorophyll fronts are typically linked across shelf-seas (e.g. at tidal-mixing fronts; [34,54]), we did not include extra analyses, using chlorophyll a derived front metrics in this study.

Seven day composites
Gdens and Fdist were based on 7 day composite front maps centred to the date of a GPS/TDR fix (figure 2). First, raw (level 0) advanced very high-resolution radiometer (AVHRR) infrared data were converted to an index of SST (level 2). SST data were then mapped across the Celtic Sea with a spatial resolution of approximately 1.2 km 2 /pixel. Thermal fronts were detected over frames of 32 × 32 pixels, using single image edge detection (SIED; [55]) with a temperature difference threshold of 0.4°C across the front [31] comparable to [35]. The SIED front map generated from a single satellite image is unsuitable for the description of fronts owing to cloud cover in the study region. Therefore, all frontal segments obtained during the 7 day window were combined to obtain a more synoptic frontal picture [31]. If cloud persisted for the entire 7 days in certain regions, then these were marked as missing in the front metrics, and any corresponding bird tracks were excluded from analysis.
Gdens was then taken as the spatially smoothed average temperature gradient across all frontal pixels detected over a 7 day composite to give a continuous distribution of frontal intensity. A Gaussian filter with a width of five pixels was used for the spatial smoothing [34,35]. This metric indicates the intensity of contemporaneous frontal structures that may be either transient or persistent in occurrence.
Fdist was taken as the distance from any point to the closest simplified front. Simplified fronts were defined, using a clustering algorithm to identify continuous contours through the strongest frontal   Table 1. Overview of the three modelling approaches applied to determine the influence of frontal activity on (a) habitat usage and the probability of a dive event, (b) dive shape (U versus V) and (c) depth and duration of U-and V-shaped dives. Terms between the curly brackets are interchangeable to represent where models were fitted separately to avoid issues associated with multicollinearity. The initial fixed component with all potential explanatory variables (before model reduction/selection) is shown. The random component comprises a random intercept of BirdID and either (i) a nested spatial correlation structure CorStructSp (a) or (ii) a nested continuous temporal correlation structure   pixels on the spatially smoothed Gdens front map [35]. This metric quantifies the distance from each dive/pseudo-absence location (see §2.5 Statistical analysis) to the continuous and discrete surface signature of more defined frontal systems, and generally excludes the smaller ephemeral features that are often picked up by the Gdens metric.

Seasonal composites
Ffreq was based on all front maps generated between June and August inclusive for each year (2012 and 2013) separately (figure 1). Maps generated before spatial smoothing were used, and Ffreq was taken as the percentage of total detections in which a frontal temperature gradient greater than or equal to 0.04°C was observed (see Average front gradient in [31]). This threshold reduced noise in front detections associated with minor discrepancies in temperature observations [32]. This metric indicates areas where fronts frequently manifest across a season and so are generally persistent and highly predictable in occurrence.

Statistical analysis
Three modelling approaches were used to investigate the influence of our three front metrics (Gdens, Fdist and Ffreq), on gannet dive behaviour (table 1). First, habitat use versus availability was modelled against frontal activity. This was achieved using generalized linear mixed effects models (GLMMs) from the MASS package in R [56] with a binomial error structure and complementary log-log (cloglog) link function [57,58]. For each dive event, the locations of five pseudo-absences were randomly selected from within the bounds of the 95% utilization distribution of the population sample (see electronic supplementary material, S2 and figure S2; [59]). Second, dive shape was modelled against frontal activity. This approach also used GLMMs from the MASS package in R [56] with a binomial error structure and complementary log-log (cloglog) link function. Finally, the influence of frontal activity on dive depth and duration was modelled separately for U-and V-shaped dives using linear mixed effects models (LMMs) from the R package nlme [60]. Across the three analyses, the Gdens, Fdist and Ffreq front metrics were tested for in separate models to avoid issues pertaining from multicollinearity. Sex was included as a fixed factor and in a two-way interaction with each front metric to account for sexual segregation in gannet foraging behaviour [47,53,61]. Time of day (DecTime) was included in the dive shape, depth and duration models as a continuous quadratic function to allow for potential changes in the profile of a dive in the middle of the day, possibly owing to the diel migration of prey above and below the thermocline [38,62]. In all models, a random intercept of bird was included to avoid pseudo-replication and account for any individual differences in dive behaviour [6]. Where appropriate, a nested autocorrelation structure was also incorporated to allow for similarities between dives that occurred in temporal clusters. In the habitat-use availability analysis, this was fitted as a spatial correlation structure, using the coordinates of each dive or pseudo-absence location on a universal transverse Mercator projection, whilst for the dive shape, depth and duration analyses, a continuous time correlation structure was used (table 1).
As GLMMs from the MASS package are estimated using penalized quasi-likelihood, maximumlikelihood (ML) selection techniques (e.g. Akaike's information criteria, AIC) were not available, and so the best temporal correlation structure (e.g. exponential, rational quadratic, autoregressive) for the random component (table 1) of the model was selected through inspection of residual plots. Model reduction of the fixed component (table 1) was performed by removing variables with small parameter estimates, relatively large standard errors, confidence intervals that passed through zero and large p-values [63][64][65].
For each LMM, the most appropriate random structure (table 1) was determined via restricted maximum-likelihood (REML) estimation [57] and selected by a comparison of AIC values and residual plots. Model selection of the fixed effects (table 1) was conducted by backwards and forwards selection via ML estimation and the use of AIC and likelihood ratio tests. The most parsimonious model was then refitted using REML to obtain parameter estimates and associated p-values [57].
Models were evaluated by plotting Pearson (GLMM) or normalized (LMM) residuals against all potential explanatory variables, bird ID, distance to colony, latitude, longitude, tag type, time and year to check for any patterns indicative of a violation of model assumptions. Fitted versus predicted values were inspected to check for satisfactory model fit, and for the binomial GLMMs the area under the receiving operator characteristic curve (AUC; [66,67]) was calculated. Pseudo-R 2 -values were generated as an indication of variance explained [68], using the MuMIn package in R [69]. All analyses were performed in R v. 3.0.2 (R Development Core Team 2013) and Matlab R2011b.

Overview of gannet diving behaviour
The total number of dives made across all foraging trips was 1901 (figure 1); 712 and 1189 of these were made by males and females, respectively. The number of dives made per complete foraging trip averaged 22.2 ± 3.8, although 17.6% of these trips had no dives. All incomplete foraging trips included dive events. Six birds did not dive at all during deployments.
Dives tended to occur in short bursts rather than being spread out equally over the foraging trip. 25.1%, 39.9% and 58.2% of all dives were followed by another dive within 5, 10 and 20 min respectively, and only 22.9% of dives were spaced more than 1 h from the previous dive. Across complete foraging trips, the number of dives made per hour ranged from 0 to 5.5 with an average of 0.9 ± 0.13. The percentage of time spent underwater during a complete foraging trip was low and ranged from 0% to 0.7% with an average of 0.1 ± 0.02%.     were negatively correlated with distance to front (Fdist) and cross-front gradient strength (Gdens), and positively correlated with seasonal front frequency (Ffreq; figure 3 and table 2).

Dive shape
V-shaped (92.1% of dives) dives were more common than U-shaped dives (7.9% of dives; figure 4). All birds performed V-shaped dives and 69.6% performed U-shaped dives. Dive shape varied significantly with front frequency (Ffreq; figure 5 and table 3). When diving in areas of high front frequency, the probability that a bird performed a U-shaped dive compared with a V-shaped dive halved (from approx. 0.12 to 0.06). Dive shape did not change in response to DecTime, any of the other front metrics (Fdist and Gdens), sex or an interaction between sex and frontal activity.

Dive depth
Overall dive depths ranged from 1.6 to 14.9 m and were greater when birds performed U-shaped compared with V-shaped dives ( figure 4). V-and U-shaped dive depths were not significantly related to any of the front metrics. Females consistently dived deeper than males when performing V-shaped dives (4.6 ± 0.4 versus 3.4 ± 0.3 m respectively; p = 0.003). U-shaped dive depths averaged 6.3 ± 0.3 m and did not differ between sexes. DecTime had no influence on dive depth for either strategy, and there was no significant effect of an interaction between sex and frontal activity.

Dive duration
Overall dive durations ranged from 0.7 to 39.1 s and were greater when birds engaged in U-shaped compared with V-shaped dives ( figure 4). V-shaped dives were significantly shorter in proximity to fronts (Fdist; table 4 and figure 6). This response was more prominent in males (table 4 and figure 6), who had significantly shorter dives than females (table 4). V-shaped dive duration was not significantly related to any of the other front metrics (Ffreq and Gdens). V-shaped dives were longer in the middle of the day (DecTime; table 4 and electronic supplementary material, figure S3). U-shaped dive duration averaged    12.9 ± 0.8 s and did not vary in response to DecTime, any of the front metrics, sex or an interaction between sex and frontal activity.

Discussion
This study provides novel insights into the influence of physical oceanography on habitat use and dive behaviour by a medium-ranging piscivorous predator. We build upon prior observations that     gannets intensify restricted search behaviours in areas with seasonally persistent fronts [35] to show that these habitats are also favoured for diving, although there was variation between the sexes. In addition, we show that when diving around fronts, gannets are half as likely to engage in U-shaped compared with V-shaped dives and the average duration of V-shaped dives is significantly shortened, which was independent of sex. Combined, these findings are of particular significance because, as well as confirming the importance of shelf-sea fronts as foraging habitat, differences in diving behaviour provide a possible functional mechanism underpinning the links between large marine predators and these physical features.

Persistent fronts as predictable foraging habitats
Our analysis of gannet dive distributions showed females and, to a lesser extent, males, preferentially dived in regions of persistent frontal activity. This adds to a growing body of evidence documenting the importance of such features as foraging habitats for marine predators [21,34,35,[70][71][72]. However, patterns were far less clear around shorter-term, more ephemeral fronts. Here, male dive probabilities were positively correlated with distance to front (Fdist) and cross-front gradient strength (Gdens), with the reverse for females. Across shelf-seas, persistence is a key feature of the bioaggregating fronts [8] that are associated with high levels of primary productivity and biomass accumulation which sustains low-to mid-trophic-level enhancement and increases prey abundance (e.g. tidal-mixing fronts; [16,20,21,73]). Moreover, dependent upon spatio-temporal scale [74,75], these features may occur in a highly predictable manner [34], which likely aids individuals in efficiently locating their prey through learning, knowledge transfer and/or memory [4,[76][77][78][79][80]. The significance of small ephemeral features (Gdens) therefore requires further study.

Subsurface dive behaviour around fronts
The principal purpose of this study was to investigate the subsurface movements of gannets in relation to shelf-sea fronts, to better understand their functional significance. In addition to providing persistent and predictable foraging habitats, fronts are also thought to increase the catchability and accessibility of prey [20,22]. In gannets, foraging strategies are especially energetically expensive [48,52,81,82], and to maximize efficiency individuals adjust their underwater movements in response to the behaviours and depth distributions of their prey [49,51,83,84]. V-shaped dives dominated gannet foraging strategies across the Celtic Sea, which possibly suggests this method of prey capture is better suited than a U-shaped dive strategy to the types of prey naturally encountered in the region (e.g. fast swimming pelagic species; [39,49,50]). We hypothesize that the shorter duration and more frequent use of passive V-shaped dives (i.e. with little or no active swim phase) around fronts is related to improvements in the availability of this prey.
Strong biophysical coupling at fronts is thought to attract large numbers of the mid-trophic-level fish (e.g. garfish, herring and mackerel) that gannets feed on [39]. Resultant high prey densities [22] may increase encounter probabilities [85] which aid in capture, allowing for faster and shorter dives [37,48]. Moreover, pursuit dives (U-shaped) are thought to be less suitable for catching highly responsive shoaling fish [86], which instead are better ambushed during fast V-shaped dives [50]. Fronts may also concentrate fish close to the surface, making them easier for gannets to catch. This is due to increased primary productivity and biomass accumulation around the near-surface thermocline [16,20,73], and because some fish may actively avoid cool bottom-boundary layer waters (e.g. mackerel; [50,87] V-shaped dives were shallower than U-shaped dives (3.43/4.62 m for males and females respectively versus 6.32 m) and so an increase in their use around fronts may reflect this shallow distribution of fish [37,51]. In some instances, gannets may additionally cue in on other marine predators such as cetaceans [88][89][90] or other seabirds [90] that also forage around fronts [71,91,92]. Under such instances, mixedspecies foraging aggregations may increase foraging success through the disorganization of school cohesiveness [93], or by preventing prey escaping to deeper waters when under attack [94].

Identifying important habitat features for foraging marine predators in dynamic ecological systems
There was a low signal-to-noise ratio across all of our analyses, and as such, the resultant variances explained were low (less than 12%) whilst model fits were sometimes poor (AUC of use-availability models less than 0.6; [66]). This was likely to be a by-product of the highly dynamic nature of shelf-sea environments coupled to the multiple trophic connections required to link physical features to gannets. However, despite this complexity, we were able to identify biologically plausible mechanisms at play, consistent with our a priori predictions about the importance of fronts as foraging habitat. Composite front mapping techniques probably played an important role in this by objectively defining and identifying frontal features across multiple spatio-temporal scales, which allowed transient ephemeral features (Gdens) to be separated from the more persistent and predictable fronts used by gannets (Ffreq).
In some instances, suitable habitat locations may be unused owing to mechanisms other than those being investigated. For example, frequent frontal zones around the coast of Ireland (figure 1) were probably avoided as a result of competition with neighbouring colonies and resultant space segregation [95], whilst those occurring around the mouth of the Bristol Channel are possibly associated with waters whose turbidity hinders the foraging ability of a visual forager such as the gannet [96]. The use of a habitat-use availability analysis, that included only areas deemed accessible for foraging as defined by kernel density analyses of GPS tracks, minimized the influence of these potentially interfering processes.
Individual and sex-specific foraging specialization may further obfuscate relationships between gannet dive behaviours and physical oceanography [6,40,41,53,61]. Correlations between the dive distributions of males and frequent frontal zones (Ffreq) were weaker than those observed with females ( figure 3c,d). Moreover, there were contrasting negative and positive relationships with distance to front (Fdist; figure 3a,b). Sex-specific differences in habitat use by gannets have also been observed in the North Sea, where males preferentially forage in near-shore coastal regions, whilst females target offshore areas of intermediate SSTs where fronts would ordinarily manifest [53]. The underlying mechanisms driving sex-specific differences in front use described here are unclear, but could be related to contrasting parental roles [97], interference competition [98], habitat segregation [53] and/or differences in nutritional requirements/prey preference [61]. As such, we highlight the need for a comprehensive knowledge of inter-and intraspecies-specific behaviours when investigating the drivers of marine predator habitat selection [6,95,99].
As well as feeding naturally on pelagic fishes, gannets also feed on discards from commercial fisheries [40,41,61], and scavenging may erode the relationships between diving gannets and fronts. To test this, we re-modelled habitat-use availability excluding all dives associated with fishing vessels (length > 15 m) and their pseudo-absences, and re-analysed changes in dive duration and shape around fronts including presence/absence of fishing vessels as a two-level factor. Dives were assumed to be at a fishing boat when within 10 km and 1 h of a logged vessel location (17.6% and 20.9% of female and male dives, respectively, as indicated by the UK vessel monitoring system provided by the Centre for Environment, Fisheries and Aquaculture Science; [39,41,100]). We found that across these re-analyses, model outputs were consistent with previous investigations that did not account for/exclude dives associated with fishing boats (for dive duration and shape re-analyses, the factor variable of fishing boat presence was not retained following model reduction). As such, whilst scavenging waste from fisheries may increase noise in this system (the influence of boats under 15 m could not be ascertained), there is no evidence that this systematically altered our findings, which is perhaps an indication of the strength of association between gannets and fronts in the Celtic Sea.

Conclusion
It has been suggested that shelf-sea fronts constitute key components in the functioning of marine ecosystems by providing marine predators with persistent, predictable and productive foraging habitats  [8,20,21]. Our work reiterates this assertion and highlights the key role these features play in shaping both the distributions and foraging behaviours of gannets. Moreover, we suggest an increase in the use of short V-shaped dive strategies around fronts reflects improved prey accessibility and catchability. Our study highlights the complexities of interactions between marine vertebrate predators and their environment, and the subsequent importance of collaboration between disciplines (spatial ecology, oceanography and remote-sensing). We show that studies combining fine-scale foraging behaviours and remotely sensed measurements of physical oceanography can provide valuable insights towards the mechanisms that drive the at-sea distributions of marine predators. As such, there is a pressing need for cross-disciplinary research when attempting to understand marine vertebrate ecology and how marine ecosystems function.
Ethics. Gannets were caught on Grassholm Island, Wales under a Schedule 1 Licence from Natural Resources Wales.
The British Trust for Ornithology Special Methods Panel approved device attachment. Blood was taken under licence from the UK Home Office. All works were also approved by the University of Exeter's ethical review panel.
Data accessibility. Supporting datasets can be found in the electronic supplementary material.