- Open Access
Inter-individual differences in foraging tactics of a colonial raptor: consistency, weather effects, and fitness correlates
Movement Ecology volume 8, Article number: 28 (2020)
Consistent inter-individual differences in behavioural phenotypes may entail differences in energy efficiency and expenditure, with different fitness payoffs. In colonial-breeding species, inter-individual differences in foraging behaviour may evolve to reduce resource use overlap among conspecifics exploiting shared foraging areas. Furthermore, individual differences in foraging behaviour may covary with individual characteristics, such as sex or physiological conditions.
We investigated individual differences in foraging tactics of a colonial raptor, the lesser kestrel (Falco naumanni). We tracked foraging trips of breeding individuals using miniaturized biologgers. We classified behaviours from GPS data and identified tactics at the foraging trip level by cluster analysis. We then estimated energy expenditure associated to each tactic from tri-axial accelerometer data.
We obtained 489 foraging trips by 36 individuals. Two clusters of trips were identified, one (SF) characterized by more static foraging behaviour and the other (DF) by more dynamic foraging behaviour, with a higher proportion of flying activity and a higher energy expenditure compared to SF. Lesser kestrels showed consistent inter-individual differences in foraging tactics across weather condition gradients, favouring DF trips as solar radiation and crosswind intensity increased. DF trips were more frequent during the nestling-rearing than during the egg incubation stage. Nestlings whose tracked parent was more prone to perform DF trips experienced higher daily mass increase, irrespective of nestling feeding rates.
Our study provided evidence that breeding lesser kestrels flexibly adopted different foraging tactics according to contingent weather landscapes, with birds showing consistent inter-individual differences in the tendency to adopt a given tactic. The positive correlation between the tendency to perform more energy-demanding DF trips and nestling growth suggests that individual differences in foraging behaviour may play a role in maintaining key life-history trade-offs between reproduction and self-maintenance.
Inter-individual differences in behavioural phenotypes, that are consistent over time and across environmental contexts, have been frequently documented in animals . Individuals can also consistently differ in how they modulate specific behaviours in accordance to spatial and temporal environmental gradients, the so-called contextual plasticity . Ultimately, inter-individual differences in behavioural phenotypes and in contextual plasticity may be associated with fitness differences among individuals , with far-reaching ecological and evolutionary consequences . For instance, individual differences in foraging behaviour may favour foraging specialization, resulting in a reduction of intraspecific competition by limiting resource use overlap [3, 8, 46].
In colonial species, where conspecifics gather in groups to reproduce close to each other and exploit shared foraging areas [11, 38], resource depletion around colony sites commonly occurs [4, 22]. Such depletion may favour the evolution and maintenance of individual foraging differences, which may be important in compensating the negative fitness effects of intraspecific competition. Indeed, individual differences in dietary preferences have been documented in several colonial vertebrates, including birds (e.g. ), pinnipeds (e.g. ) and terrestrial mammals (e.g. ). Individual foraging differences can also result in inter-individual differences in foraging tactics [5, 9, 17, 29], which we define as a distinguishable combination of behavioural patterns (i.e. multiple behaviours) shared by different individuals to search for food (similarly to [17, 39]). In northern gannets Morus bassanus, the analysis of both food boluses and blood isotopes, combined with the analysis of at-sea foraging behaviour, has uncovered individual differences in foraging tactics, with some birds exploiting consistently and more frequently than others fishing vessels discards . In the same species, the analysis of foraging trips of birds breeding in two large colonies revealed consistent individual preferences in foraging areas during the nestling-rearing stage, but also large inter-individual differences in prey searching behaviour along environmental gradients . Within-population differences in foraging behaviour can also be unrelated to individual preferences and rather arise from individual characteristics, such as sex , age , reproductive stage , personality  or physiological status differences .
The lesser kestrel Falco naumanni is a small (ca. 120 g) colonial diurnal raptor, which mostly breeds in holes and crevices of buildings in towns and cities, and forages in farmland areas surrounding breeding sites . The species shows flexible foraging behaviour, whereby both flight and hunting mode vary in accordance to weather conditions: energy-saving soaring-gliding flight is more frequently adopted than energy-expensive flapping flight when solar radiation is high, and perch-hunting is more frequently used than flight-hunting when both wind speed and solar radiation are low . However, it is as yet unknown whether individuals consistently differ in their foraging tactic across weather condition gradients, or whether individuals differ in their behavioural response to weather conditions.
In this study, we characterized the foraging tactics adopted by breeding lesser kestrels that were tracked by miniaturized biologgers (including both a GPS and a tri-axial accelerometer) across multiple foraging trips. We classified behaviours from movement data and identified tactics at the foraging trip level based on the combination of different behaviours occurring within each trip, while energy expenditure associated to each tactic was estimated by accelerometer data. According to previous knowledge about lesser kestrels foraging ecology , we expected birds to adopt two main foraging tactics: a more energy-demanding tactic whereby birds mainly search for prey while flying within a foraging trip (dynamic foraging, DF), and a less energy-demanding tactic, resulting from trips with prolonged perching while waiting for prey detection (static foraging, SF).
We then investigated 1) whether there were inter-individual differences in the tendency to adopt a given foraging tactic across weather condition gradients (solar radiation, rain, wind), expecting birds mainly to adopt the DF tactic whenever conditions are favourable for soaring-gliding, i.e. with high solar radiation  and wind assistance (tailwind or crosswind) at trip departure (e.g. ). Furthermore, we explored 2) whether the tendency to adopt a specific tactic was explained by individual characteristics, such as sex and breeding stage (incubation or nestling-rearing). Foraging behaviour of lesser kestrels can indeed vary markedly between males and females and between the incubation and nestling-rearing stages . Finally, we investigated 3) the association between the tendency to perform a specific tactic and fitness-related traits. We expected the more energy-demanding DF tactic to be adopted preferentially by individuals in better body condition and to be associated with improved fitness, as estimated by higher breeding success and larger nestling body mass increase.
Target species, study area and general field methods
The lesser kestrel is a sexually dimorphic species, females being ca. 15% heavier than males . Females lay up to 5 eggs that both parents incubate for ca. 30 days. After hatching, both parents feed the nestlings until fledging, which takes place at 35–40 days of age. The study was carried out in the city of Matera (southern Italy; 40°39′ N, 16°36′ E), hosting a colony of ca. 1000 breeding pairs . We relied on nest-boxes placed on terraces of buildings in the old town, which were monitored 2–3 times per week to obtain detailed data about reproductive stage (laying date, incubation, hatching, nestling body mass at ca. 7 and ca. 14 days after hatching of the first egg) [53, 54]. Breeding individuals were captured by hand within nestboxes during the late incubation and early nestling-rearing stage. Upon capture, birds were individually marked, and body mass (using an electronic scale, accuracy 0.1 g) and keel length (using a dial calliper, accuracy 0.1 mm) were recorded.
GPS deployment and identification of foraging trips
We equipped 36 breeding lesser kestrels (13 females and 23 males) with Axy-Trek biologgers, including a GPS and a tri-axial accelerometer (TechnoSmArt Europe S.r.l., Rome, Italy), using a backpack Teflon harness. Devices (including the harness) weighed between 5.9 and 7.2 g, on average accounting for 4.5% of body mass (range: 3.8–5.5%) (hereafter, relative load of device). The accelerometer was set to record data at 25 Hz and the GPS to record one position per minute from 05:00 to 21:00 local time (i.e. ca. 20 min before sunrise and ca. 30 after sunset; devices were switched off during night-time to preserve battery power). Birds were tagged in the morning and devices were set to start the following day in order to collect data when the tagged birds were likely inured to the device. After 2–5 days, birds were recaptured and the device was removed. Movement data were collected during June 1–20 of the breeding seasons 2016, 2017 and 2018, when pairs were in the late incubation or early nestling-rearing stage. Only one member of each pair was tracked and none of the birds was tagged more than once.
Foraging trips were identified as those tracks starting and ending within a 50-m buffer around the nest or the roosting site, and heading to the rural surroundings, by means of ESRI ArcMap 10.2.1. Since the devices switched on at 05:00 local time, in some cases the first position of the first foraging trip of the day was already in the countryside surrounding the town; when the distance between such position and the nesting site was > 2 km, the foraging trip was discarded and not included in any analysis. We did not consider as foraging trips all those excursions which covered the urban area only, identified by means of the 2012 CORINE Land Cover (CLC) map (codes 111 and 112, respectively continuous and discontinuous urban habitat), because birds generally do not forage in the urban area (authors’ pers. obs.). Each trip was classified as occurring during incubation if only eggs were present in the nest of the target individual on the date when the foraging trip was performed, or as occurring during the nestling-rearing stage if at least one nestling was present in the nest on the date when the trip was performed.
Identification and characterization of foraging tactics
To identify behaviours adopted by lesser kestrels during foraging trips, we applied the Expectation Minimization binary Clustering (EMbC) algorithm to GPS data by means of the R package “EMbC” . The EMbC is a classification algorithm based on maximum likelihood which assigns a behavioural mode to each GPS position according to instantaneous velocity and turning angle between successive positions. The algorithm assigns positions to four behavioural modes (see Fig. S1): 1) low velocity and low turns, which we interpreted as ‘perching’ behaviour (see Fig. S2); low velocity and high turns, representing ‘intensive search’; 3) high velocity and low turns, representing ‘relocation’; 4) high velocity and high turns, representing ‘extensive search’ [27, 39]. As the EMbC algorithm disregards the temporal information, we accounted for the possible incorrect labelling of positions when a long-term predominant behavioural mode occurred by applying a post-processing smoothing using the smth() function (with default parameters) (see ; and Fig. S1).
To identify the foraging tactics, we applied a cluster analysis to the percentage of the four behaviours occurring in each trip . Cluster analysis was performed with a K-means procedure by means of the R package “stats” . The optimal number of clusters was assessed by means of the NbClust procedure from the R package “NbClust” , which computes 30 indexes for determining the optimal number of clusters, including the option of no clustering (one cluster only). It then suggests the best number of clusters based on the majority consensus rule. As the NbClust procedure identified two clusters as the best number, we applied the K-means algorithm with K = 2 over 10,000 iterations.
We then assessed whether the two identified trip clusters (trip types, hereafter) affected variation in spatio-temporal trip descriptors (trip duration, trip length, maximum distance from the nest site and tortuosity, i.e. ratio between total trip length and the maximum distance from the nest site ;) by means of linear mixed models (LMMs), including individual identity as a random intercept effect to account for non-independence of trips performed by the same individual.
Based on tri-axis accelerometer data, we calculated the overall dynamic body acceleration (ODBA) for each foraging trip, smoothing total acceleration over 1 s [59, 68]. ODBA is considered a proxy of energy expenditure in birds [23, 68]. It positively correlates with O2 consumption rates and CO2 production in great cormorants (Phalacrocorax carbo)  and with heart rate in two vulture species (Gyps fulvus and G. himalayensis) . We investigated whether energy expenditure was affected by trip type by means of a LMM including individual identity as a random intercept effect. Accelerometer data were available for 34 out of 36 tracked birds.
LMMs were fitted using the lmer function of the R package “lme4” . Residuals did not significantly deviate from a normal distribution.
Environmental factors affecting foraging tactics
To investigate the effect of environmental conditions on the tendency to perform different trip types, we first associated to each trip the following variables: 1) solar radiation (W/m2) at departure, which seems to be determinant for performing soaring-gliding flight ; 2) presence of rain during the trip (hereafter, ‘presence of rain’; 0 = rain absent; 1 = rain present), which we hypothesized may negatively affect the likelihood of performing foraging in flight. Rain was not considered as a continuous variable since it occurred in 14% of foraging trips only, and considering it as continuous would have resulted in a very skewed variable with an excess of zeroes; 3) tail-wind (TWC) and 4) cross-wind components (CWC), both of which are known to affect movement activity in soaring-gliding raptors . To control for potential differences in foraging behaviour between foraging habitats, we also computed 5) the percentage of positions in arable lands (the main habitat used for foraging, see below) for each trip (hereafter, ‘time in arable lands’). Nestbox identity was not included in the models since all tagged birds belonged to different nests, with the exception of two nestboxes which were sampled twice in different years but they were occupied by different individuals.
Solar radiation and rain data were recorded at a weather station located at 8 km from the nest sites (Matera, Contrada Matinelle, 40°41′ N; 16°31′ E). Wind data (speed and direction) were recorded at a different weather station, located at 15 km from the nest sites (Grottole 40°37′ N; 16°26′ E). All weather data were recorded at hourly intervals, and were associated to the GPS position that was closest in time.
TWC and CWC were calculated for each trip based on the mean value of wind speed and direction (WS and WD respectively) at time of departure and at time of returning, and the direction of the trip (TD), as follows:
TD, which we assumed to reflect the direction of the goal area, was calculated as the angle between the N-S axis (directed northwards) and the position of the farthest point of the trip from the nesting site. Positive TWC values imply that a bird flew globally with tail-wind on its way out of the colony towards the foraging grounds, whereas negative TWC values indicate the opposite (outgoing flights with headwinds). Large CWC values mean that a bird flew on average with high side-wind during the foraging trip.
To calculate the proportion of time spent in arable lands during each foraging trip, we assigned all GPS positions, excluding those identified as relocation by the EMbC, to the corresponding habitat type from CLC by means of ESRI ArcMap 10.2.1. We pooled together those CLC habitat types that were similar in habitat and structure, obtaining 6 habitat classes: artificial landscape (continuous and discontinuous urban fabric, infrastructures, industrial areas), arable lands, permanent crops (tree plantations, olive groves, vineyards), grasslands (pastures and natural grasslands), heterogeneous agricultural areas (annual crops associated with permanent crops, complex cultivation patterns, agro-forestry areas), and wooded areas (forests and bushes). Each trip was then characterized by the percentage of positions occurring in each habitat class (time spent in each habitat). While foraging, birds spent most of the time in arable lands (median 70.6%, 25th – 75th percentiles: 23.3 – 94.1%), and time in arable lands was negatively correlated with time in grasslands (r = − 0.69, n = 489 trips), the second most frequently used habitat (median 0.0%, 25th – 75th percentiles: 0.0 – 25.0%).
The probability to perform a given trip type (0 = SF, 1 = DF) was modelled by means of binomial generalised linear mixed models (GLMMs), with solar radiation, TWC, CWC, and presence of rain as fixed predictors, controlling for time in arable lands, breeding stage, sex and sampling year. Individual identity was included as a random intercept effect to control for non-independence of prey searching behaviour performed by the same individuals. All predictors were standardized (mean = 0 and SD = 1). Because of size and morphological differences, sexes may differ in their behavioural response to environmental variables; we hence included in the initial model all two-way interactions between sex and each weather variable (solar radiation, TWC, CWC, and presence of rain). The final binomial GLMM was obtained after removing weak (95% CI of parameter estimates intersecting 0) interactions in a single step. GLMMs were not overdispersed (ϕ always < 1.05).
Individual differences in foraging tactics and their correlates
The random intercept effect of the final binomial GLMM describes the extent to which individuals preferentially perform different trip types (i.e. whether foraging tactics can be regarded as an individual-specific trait). Intra-idividual consistency of the probability of performing DF trips was estimated as the proportion of variance explained by the random intercept effect, accounting for variance explained by fixed effects (adjusted repeatability, Radj). Radj was computed using the observation-level variance obtained via the delta method  and significance was assessed by a likelihood ratio test .
To investigate individual differences in the behavioural response to environmental gradients (behavioural reaction norms), which represent the degree of contextual plasticity (see ), we ran four binomial GLMMs with trip type as the binary dependent variable and sex, breeding stage, sampling year and one weather variable (solar radiation, presence of rain, TWC, or CWC) at a time as predictors, while including an individual-level random slope effect for that weather variable. Random slope models were fitted by including in binomial GLMMs one weather variable at a time to avoid model overparametrization and lack of convergence . Significance of the random slope effect was tested by a likelihood ratio test . Models were not overdispersed (ϕ always < 1.05). All binomial GLMMs were fitted using the glmer function of the R package “lme4” .
The individual tendency to perform DF trips was expressed as the individual-level random intercept estimate (hereafter, individual intercept) from the final binomial GLMM, higher values implying a stronger tendency to perform DF trips. Individual intercepts were computed as the conditional modes of the random effect evaluated at the parameter estimates (a.k.a. Best Linear Unbiased Predictors in LMMs). Uncertainty of individual intercepts was estimated by a simulation approach (n = 10,000 simulations; see ), and expressed as the SD of the simulations, using the REsim function of the R package “merTools” .
We tested the associations between the individual tendency to perform DF trips by the tracked parent and several fitness proxies, namely body condition of the tracked individual, its breeding success, and the daily body mass increase (DBMI) of its nestlings during the early nestling-rearing stage. We also investigated the association of the tendency to adopt the DF tactic with feeding frequency (number of foraging trips/h during the tracking period) and with the relative load of device (mass of biologging device relative to body mass). The latter association was tested (separately for males and females) to assess the possible effects of device load on foraging behaviour. All these associations were tested by computing weighted correlation coefficients (rw), where the weighting variable was the inverse of the SD of the individual tendency. Such weighting should at least partly account for uncertainty of conditional modes, as advocated by Houslay and Wilson . Weighted correlation coefficients were computed using the weightedCorr function of the R package “wCorr” . Significance of rw was tested by randomization, randomly shuffling variables 9999 times, and computing the probability of observing a more extreme value than the observed one .
To estimate body condition of tracked birds, we computed the scaled mass index (hereafter, SMI) (body mass scaled by a skeletal trait, in our case keel length ). SMI was calculated as in Podofillini et al. . To remove heterogeneity in SMI related to sex (see ), we computed the residuals of a linear model of SMI with sex as a predictor (hereafter, residual SMI).
The breeding success of the tracked birds was estimated as the fate of the brood, which was either coded as failed (0 = no nestlings alive at 14 days) or successful (1 = at least one nestling alive at 14 days after hatching of the first egg).
Nestling DBMI was computed as the mean of the daily relative body mass increase among all nestlings of a brood k between ca. 7 (mean ± SD = 6.6 ± 1.8 days) and ca. 14 days (mean ± SD = 14.7 ± 1.7 days) from hatching, as follows:
where BM is body mass of nestling j, n is the day post-hatching at which the first record of BM was taken and i is the number of days elapsed between the first and second measure of BM of a nestling, m is the brood size at day(n+i) (range: 1–4). We assumed that nestling DBMI denoted the ability of parents to foster the growth of their offspring, and that large value of DBMI could be considered as a proxy for high parental investment and high nestling fitness. To remove heterogeneity among sampling years and different brood sizes on nestling DBMI, we computed the residuals from a linear model of nestling DBMI including year and brood size at day n+i as predictors (hereafter, residual nestling DBMI). The correlation test was based on data from 22 parents whose eggs hatched and whose offspring were alive at 14 days from hatching of the first egg.
Feeding frequency was computed only for birds which were tracked during the nestling-rearing stage (n = 14). To remove heterogeneity among sampling years, sexes, and variation in brood size on feeding frequency (see [19, 31]), we computed residuals from a linear model of feeding frequency including year, sex, and brood size during the tracking days as predictors (hereafter, residual feeding frequency).
All analyses were ran on R ver 3.6.2 .
Identification and characterization of foraging tactics
We obtained 489 foraging trips from 36 breeding birds, the mean value being 14 trips (± 11 SD) per individual, ranging from a minimum of 1 to a maximum of 45 trips per individual (see also Table 1). The cluster analysis identified two clusters of trips, which we interpreted as two main foraging tactics (Fig. 1). The first cluster included trips characterized by high frequency of perching (mean ± SD proportion over all GPS positions of a trip: 0.53 ± 0.17), low frequency of intensive (0.15 ± 0.01) and extensive search (0.04 ± 0.05), and lower relocation positions (0.28 ± 0.12) compared to the other cluster. This cluster of trips was likely reflecting a relatively more SF tactic (Table 2). The second cluster was characterized by trips with a more dynamic and exploratory behaviour, with birds mostly searching for food while flying (perching: 0.07 ± 0.09; intensive search: 0.26 ± 0.16; extensive search: 0.07 ± 0.07; relocation: 0.60 ± 0.01), likely reflecting a relatively more DF tactic (Table 2). On average, SF trips lasted longer and were associated to lower values of ODBA compared to DF ones, whereas all other trip descriptors were not markedly different (Table 2). Figure 2 shows representative examples of DF and SF foraging trips performed by a single individual during both the incubation and the nestling rearing stages. Differences in the temporal sequence of behaviours clearly highlight that DF trips were characterized by more time actively spent searching for food, whereas SF trips showed prolonged perching periods (Fig. 2). There was no apparent spatial differentiation in exploited areas between SF and DF trips (Fig. 3).
Environmental factors affecting foraging tactics and individual differences in foraging behaviour
Tracked birds preferentially performed DF trips (65% of all trips; intercept-only binomial GLMM, estimate [95% CI] = 0.62 [0.18, 1.06], Z = 2.76, P = 0.006). The probability of performing DF trips varied among years and was moderately positively affected by solar radiation (effect size r = 0.28), with CWC and breeding stage having somewhat weaker effects (r = 0.13 and 0.20, respectively) (Table 3). With low CWC and solar radiation, birds mostly adopted SF, whereas the probability of adopting the DF tactic increased as CWC and solar radiation increased (Fig. 4). Birds were more likely to perform DF trips during the nestling-rearing compared to the incubation stage (Table 3). All other predictors had a negligible effect on the probability to perform DF trips (95% CI including 0, all r < 0.06). Males and females did not markedly differ in the probability to perform DF trips according to weather condition gradients (two-way interactions between sex and each weather variable, r always < 0.05, P always > 0.29).
The proportion of variance explained by individual identity was low (Radj = 0.13, χ2 = 15.17, df = 1, P < 0.001), indicating that individuals mostly adopted a flexible foraging behaviour. However, variation among individuals in the tendency to perform DF trips was substantial, with values ranging between −1.15 and 1.65 (logit scale; Fig. S3). Individual differences in the behavioural response to weather conditions gradients were negligible in most cases (random slope effects, solar radiation: χ2 = 0.09, df = 2, P = 0.96; presence of rain: χ2 = 2.02, df = 2, P = 0.36; TWC: χ2 = 5.60, df = 2, P = 0.06; CWC: χ2 = 1.10, df = 2, P = 0.58). Hence, individuals consistently differed in foraging tactic across solar radiation and CWC gradients (Fig. 4).
Correlates of individual variation in foraging tactics
The individual tendency to perform DF trips was very weakly correlated with residual SMI (rw = −0.03, n = 35, Prand = 0.87) and breeding success (rw = 0.15, n = 35, Prand = 0.41), but it was moderately positively correlated with the residual nestling DBMI (rw = 0.45, n = 22, Prand = 0.038) (Fig. 5). The body mass of nestlings increased on average by 11.4% (± 4.7 SD) between the two measurements, ranging between 4 and 29%. The latter correlation was not explained by a higher feeding frequency of birds performing more DF trips, as feeding frequency was weakly correlated with the tendency to perform DF trips (rw = 0.16, n = 14, Prand = 0.58) (Fig. 5).
Finally, the relative load of devices was very weakly associated with the individual tendency to perform DF trips in both sexes (males: rw = 0.08, n = 23, Prand = 0.71; females: rw = − 0.07, n = 12, Prand = 0.83).
By GPS-tracking individuals over multiple foraging trips, we investigated the degree of individual specialization in foraging tactics of breeding lesser kestrels. Foraging tactics adopted by birds during foraging trips varied between two extremes. On the one side, birds performed more static foraging (SF) trips, characterized by high frequency of perching and low proportion of searching behaviour, which were also long-lasting. On the other side, high frequency of both relocation and intensive search and low frequency of perching resulted in more dynamic foraging (DF) trips, that were short-lasting. Unsurprisingly, DF trips were associated with 1.7-fold higher ODBA compared to SF trips. Although birds performed more DF trips overall, the probability of performing DF trips increased with increasing solar radiation and crosswind, and was higher during the nestling-rearing than during the egg incubation stage. On top of this, individuals significantly differed in their tendency to perform different trip types, with no evidence of contextual plasticity (i.e. all individuals modulated their foraging tactic in a similar way in relation to weather condition gradients). The tendency to perform energy-demanding DF trips by the attending tracked parent was associated with higher nestling body mass increase, but not with a higher nestling feeding frequency.
The occurrence of both wide-ranging and relatively more static foraging tactics has been documented in many predators, including reptiles (e.g. ), fish (e.g. ), birds (e.g. ) and mammals . However, to our knowledge, the alternation of foraging tactics has seldom been analysed at the individual level. It is being increasingly appreciated that animal movement patterns and behaviour are shaped by the so-called ‘energy landscape’, i.e. the variation in the cost of transport across time and space, determined by the interaction between static landscape features and dynamic environmental conditions [1, 23, 61]. Lesser kestrels heavily rely on thermal soaring and gliding for foraging, especially when solar radiation is high . We showed that foraging lesser kestrels mainly performed SF trips when weather conditions were not ideal for soaring-gliding, i.e. with low solar radiation  and weak crosswinds, which are known to affect movement patterns in soaring raptors . As expected according to previous studies [43, 66], performing SF trips was an energy-sparing tactic compared to performing DF trips. Birds may thus use DF or SF depending on the relative energy advantage, determined by the wind and solar radiation landscape surrounding the breeding colony. Birds were mostly adopting the more energy-demanding tactic (i.e., DF) only when the energy landscape surrounding breeding sites allowed it.
Inter-individual differences in foraging tactics of colonial vertebrates may originate from divergent selection to mitigate intraspecific competition, by limiting resource use overlap among individuals sharing foraging areas (e.g., ). Under this scenario, we would not expect any significant fitness difference between individuals adopting preferentially one or the other foraging tactic. Although we did not measure other fitness components (i.e., parental survival), the higher body mass increase of nestlings whose parent was performing more DF trips might suggest that the tendency to adopt a given foraging tactic is related to individual characteristics, such as age/experience or physiological status [2, 26] rather than to the mitigation of intraspecific competition. Admittedly, the higher body mass increase of nestlings whose parent was preferentially performing DF trips should be viewed with caution because we could assess the behaviour of a single parent only. Notwithstanding, it suggests that, even when considering the uniparental contribution to nestling growth, an increased energy expenditure during foraging could result in faster nestling growth and thus better fitness prospects.
The higher body mass increase of nestlings of tracked parents was not a by-product of higher nestling feeding frequencies of birds preferentially performing DF trips. This raises the question of why nestlings reared by parents more prone to perform DF trips grew more. We may speculate that parents preferentially performing DF trips could have delivered more energy-rich prey to their nestlings than those performing more SF trips, resulting in faster mass growth. For instance, parents performing more DF trips may have been mostly targeting large crickets, that are the preferred lesser kestrel prey and have a higher fat content compared to e.g. vertebrate preys , whereas those performing more SF trips may have been targeting larger (but less energetic) prey items, such as lizards and mammals. Indeed, sit-and-wait predators generally catch larger prey compared to those taken by active predators .
The shorter duration of DF compared to SF trips could be related to the DF tactic being associated with group foraging by means of local enhancement processes for food finding . Such processes imply that individuals searching for food are attracted by feeding aggregations of other individuals and do not need to spend time searching for productive food patches [33, 40]. Social foraging should increase individual foraging efficiency when exploiting ephemeral and unpredictable resources . In the study area, we indeed regularly observed aggregations of foraging lesser kestrels performing DF to catch large orthopterans flushed during harvesting operations (see also ), while birds perching on wires or poles were generally observed alone.
We provided evidence for both individual foraging specialization and high flexibility in foraging tactics, with individuals consistently modulating their foraging tactic according to the concomitant weather landscape. The two foraging tactics were not equivalent in term of energy expenditure and consequences for fitness. Parents preferentially performing DF trips may have exploited group foraging to target more profitable, energy-rich prey in a shorter amount of time, resulting in increased nestling growth, though at the cost of a higher energy expenditure for transport. We may speculate that parents mainly performing DF trips may favour offspring growth over self-maintenance, whereas those mainly performing SF trips may do the opposite. Our results therefore suggest that inter-individual differences in foraging tactics may play a role in maintaining variation within populations in key intergenerational life-history trade-offs, such as those between parental reproductive effort and offspring survival, or between offspring growth and parental self-maintenance.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on request.
CORINE Land Cover
Daily body mass increase
Expectation Minimization binary Clustering
Generalised linear mixed model
Linear mixed model
Overall dynamic body acceleration
Scaled mass index
Direction of the trip
Amélineau F, Péron C, Lescroël A, Authier M, Provost P, Grémillet D. Windscape and tortuosity shape the flight costs of northern gannets. J Exp Biol. 2014;217(6):876–85.
Angelier F, Shaffer SA, Weimerskirch H, Trouvé C, Chastel O. Corticosterone and foraging behavior in a pelagic seabird. Physiol Biochem Zool. 2007;80(3):283–92.
Araújo MS, Bolnick DI, Layman CA. The ecological causes of individual specialisation. Ecol Lett. 2011;14(9):948–58.
Ashmole NP. The regulation of numbers of tropical oceanic birds. Ibis. 1963;103:458–73.
Austin R, De Pascalis F, JPY A, Haakonsson J, Votier SC, Ebanks-Petrie G, et al. A sex-influenced flexible foraging strategy in a tropical seabird, the magnificent frigatebird. Mar Ecol Prog Ser. 2019;611:203–14.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
Benhamou S. How to reliably estimate the tortuosity of an animal's path: straightness, sinuosity, or fractal dimension? J Theor Biol. 2004;229:209–20.
Bolnick DI, Svanbäck R, Fordyce JA, Yang LH, Davis JM, Hulsey CD, Forister ML. The ecology of individuals: incidence and implications of individual specialization. Am Nat. 2003;161(1):1–28.
Boyd C, Punt AE, Weimerskirch H, Bertrand S. Movement models provide insights into variation in the foraging effort of central place foragers. Ecol Model. 2014;286:13–25.
Catry I, et al. Easy but ephemeral food: exploring the trade-offs of agricultural practices in the foraging decisions of lesser kestrels on farmland. Bird Study. 2014;61:447–56.
Cecere JG, Bondì S, Podofillini S, Imperio S, Griggio M, Fulco E, Curcio A, Ménard D, Mellone U, Saino N, Serra L, Sarà M, Rubolini D. Spatial segregation of home ranges between neighbouring colonies in a diurnal raptor. Sci Rep. 2018;8(1):11762.
Ceia FR, Ramos JA. Individual specialization in the foraging and feeding strategies of seabirds: a review. Mar Biol. 2015;162:1923–38.
Charrad M, Ghazzali N, Boiteau V, Niknafs A. NbClust: Determining the Best Number of Clusters in a Data Set. 2015. Available from: https://cran.r-project.org/package=NbClust.
Cherel Y, Kernaléguen L, Richard P, Guinet C. Whisker isotopic signature depicts migration patterns and multi-year intra-and inter-individual foraging strategies in fur seals. Biol Lett. 2009;5(6):830–2.
Cryan PM, Stricker CA, Wunder MB. Evidence of cryptic individual specialization in an opportunistic insectivorous bat. J Mammal. 2012;93(2):381–9.
Dall SR, Bell AM, Bolnick DI, Ratnieks FL. An evolutionary ecology of individual differences. Ecol Lett. 2012;15(10):1189–98.
De Pascalis F, Imperio S, Benvenuti A, Catoni C, Rubolini D, Cecere JG. Sex-specific foraging behaviour is affected by wind conditions in a sexually-size dimorphic seabird. Anim Behav. 2020.
Dingemanse NJ, Kazem AJ, Réale D, Wright J. Behavioural reaction norms: animal personality meets individual plasticity. Trends Ecol Evol. 2010;25(2):81–9.
Donázar JA, Negro JJ, Hiraldo F. Functional analysis of mate-feeding in the lesser kestrel Falco naumanni. Ornis Scand. 1992;23:190–4.
Duriez O, Kato A, Tromp C, Dell'Omo G, Vyssotski AL, Sarrazin F, Ropert-Coudert Y. How cheap is soaring flight in raptors? A preliminary investigation in freely-flying vultures. PLoS One. 2014;9(1):e484887.
East P, Magnan P. Some factors regulating piscivory of brook trout, Salvelinus fontinalis, in lakes of the Laurentian shield. Can J Fish Aquat Sci. 1991;48(9):1735–43.
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(3):613–25.
Elliott KH, Chivers LS, Bessey L, Gaston AJ, Hatch SA, Kato A, et al. Windscapes shape seabird instantaneous energy costs but adult behavior buffers impact on offspring. Mov Ecol. 2014;2(1):17.
Emad A, Bailey P (2017). wCorr: weighted correlations. R package version 1.9.1. https://CRAN.R-project.org/package=wCorr.
Focardi S, Cecere JG. The Lévy flight foraging hypothesis in a pelagic seabird. J Anim Ecol. 2014;83(2):353–64.
Frankish CK, Manica A, Phillips RA. Effects of age on foraging behavior in two closely related albatross species. Mov Ecol. 2020;8(1):7.
Garriga J, Palmer JRB, Oltra A, Bartumeus F. Expectation-maximization binary clustering for Behavioural annotation. PLoS One. 2016;11(3):e0151984.
Greene CH. Patterns of prey selection: implications of predator foraging tactics. Am Nat. 1986;128:824–39.
Heithaus MR, Dill LM. Feeding strategies and tactics. In: Encyclopedia of Marine Mammals (eds W.F. Perrin, B. Würsig & J.G.M. Thewissen): Academic Press, San Diego; 2009. p. 414–23.
Hernández-Pliego J, Rodríguez C, Dell’Omo G, Bustamante J. Combined use of tri-axial accelerometers and GPS reveals the flexible foraging strategy of a bird in relation to weather conditions. PLoS One. 2017a;12(6):e0177892.
Hernández-Pliego J, Rodríguez C, Bustamante J. A few long versus many short foraging trips: different foraging strategies of lesser kestrel sexes during breeding. Mov Ecol. 2017b;5(1):8.
Houslay TM, Wilson AJ. Avoiding the misuse of BLUP in behavioural ecology. Behav Ecol. 2017;28(4):948–52.
Jackson AL, Ruxton GD, Houston DC. The effect of social facilitation on foraging success in vultures: a modelling study. Biol Lett. 2008;4(3):311–3.
Jaeger BC, Edwards LJ, Das K, Sen PK. An R 2 statistic for fixed effects in the generalized linear mixed model. J Appl Stat. 2017;44(6):1086–105.
Klaassen RH, Hake M, Strandberg R, Alerstam T. Geographical and temporal flexibility in the response to crosswinds by migrating raptors. Proc Biol Sci. 2010;278(1710):1339-46.
Knowles J.E. and C. Frederick (2019). merTools: Tools for Analyzing Mixed Effect Regression Models. R package version 0.5 https://CRAN.R-project.org/package=merTools.
La Gioia G., Melega L. & Fornasari L. 2017. Piano d’Azione Nazionale per il grillaio (Falco naumanni). Quad. Cons. Natura, 41, MATTM - Istituto Superiore per la Protezione e la Ricerca Ambientale (ISPRA).
Lascelles BG, Taylor PR, Miller MGR, Dias MP, Oppel S, Torres L, et al. Applying global criteria to tracking data to define important areas for marine conservation. Divers Distrib. 2016;22(4):422–31.
Louzao M, Wiegand T, Bartumeus F, Weimerskirch H. Coupling instantaneous energy-budget models and behavioural mode analysis to estimate optimal foraging strategy: an example with wandering albatrosses. Mov Ecol. 2014;2:8.
Machovsky-Capuska GE, Hauber ME, Libby E, Amiot C, Raubenheimer D. The contribution of private and public information in foraging by Australasian gannets. Anim Cogn. 2014;17:849–58.
Manly BFJ. Randomization and Monte Carlo methods in biology London. UK: Chapman and Hall; 1991.
McLaughlin RL. Search modes of birds and lizards: evidence for alternative movement patterns. Am Nat. 1989;133(5):654–70.
Nadjafzadeh M, Hofer H, Krone O. Sit-and-wait for large prey: foraging strategy and prey choice of White-tailed eagles. J Ornithol. 2016;157(1):165–78.
Nakagawa S, Schielzeth H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol Evol. 2013;4(2):133–42.
Nakagawa S, Johnson PCD, Schielzeth H. The coefficient of determination R2 and intra-class correlation coefficient from generalized linear mixed-effects models revisited and expanded. J R Soc Interface. 2017;14(134):20170213.
Navarro J, Grémillet D, Ramirez FJ, Afán I, Bouten W, Forero MG. Shifting individual habitat specialization of a successful predator living in anthropogenic landscapes. Mar Ecol Prog Ser. 2017;578:243-51.
Overington SE, Dubois F, Lefebvre L. Food unpredictability drives both generalism and social foraging: a game theoretical model. Behav Ecol. 2008;19(4):836–41.
Patrick SC, Bearhop S, Grémillet D, Lescroël A, Grecian WJ, Bodey TW, et al. Individual differences in searching behaviour and spatial foraging consistency in a central place marine predator. Oikos. 2014;123(1):33–40.
Patrick SC, Weimerskirch H. Consistency pays: sex differences and fitness consequences of behavioural specialization in a wide-ranging seabird. Biol Lett. 2014a;10(10):20140630.
Patrick SC, Weimerskirch H. Personality, foraging and fitness consequences in a long lived seabird. PLoS One. 2014b;9(2):e87269.
Peig J, Green AJ. New perspectives for estimating body condition from mass/length data: the scaled mass index as an alternative method. Oikos. 2009;118:1883–91.
Perry G. The evolution of search modes: ecological versus phylogenetic perspectives. Am Nat. 1999;153:98–109.
Podofillini S, Cecere JG, Griggio M, Curcio A, De Capua EL, Fulco E, Pirrello S, Saino N, Serra L, Visceglia M, Rubolini D. Home, dirty home: effect of old nest material on nest-site selection and breeding performance in a cavity-nesting raptor. Curr Zool. 2018;6:693–702.
Podofillini S, Cecere JG, Griggio M, Corti M, De Capua EL, Parolini M, Saino N, Serra L, Rubolini D. Benefits of extra food to reproduction depend on maternal condition. Oikos. 2019;128:943–59.
Pöysä H. Group foraging in patchy environments: the importance of coarse-level local enhancement. Ornis Scand. 1992;23:159–66.
R Core Team. R: A language and environment for statistica computing. R Fundation for Statistical Computing, Vienna, Austria. 2020. https://www.R-project.org/.
Réale D, Dingemanse NJ, Kazem AJ, Wright J. Evolutionary and ecological approaches to the study of personality. Philos Trans R Soc B. 2010;365:3937–46.
Rodríguez C, Tapia L, Kieny F, Bustamante J. Temporal changes in lesser kestrel (Falco naumanni) diet during the breeding season in southern Spain. J Raptor Res. 2010;44(2):120–9.
Ropert-Coudert Y, Wilson RP, Gremillet D, Kato A, Lewis S, Ryan PG. Electrocardiogram recordings in free-ranging gannets reveal minimum difference in heart rate during flapping versus gliding flight. Marine Ecol Prog Ser. 2007;331:309.
Schielzeth H, Forstmeier W. Conclusions beyond support: overconfident estimates in mixed models. Behav Ecol. 2009;20(2):416–20.
Shepard EL, Wilson RP, Rees WG, Grundy E, Lambertucci SA, Vosper SB. Energy landscapes shape animal movement ecology. Am Nat. 2013;182(3):298–312.
Sidiropoulos N, Sohi SH, Pedersen TL, Porse BT, Winther O, Rapin N, Bagger FO. SinaPlot: an enhanced chart for simple and truthful representation of single observations over multiple classes. J Comput Graph Stat. 2018;27(3):673–6.
Smith BR, Blumstein DT. Fitness consequences of personality: a meta-analysis. Behav Ecol. 2008;19(2):448–55.
Stamps J, Groothuis TGG. The development of animal personality: relevance, concepts and perspectives. Biol Rev. 2010;85:301–25.
Votier SC, Bearhop S, Witt MJ, Inger R, Thompson D, Newton J. Individual responses of seabirds to commercial fisheries revealed using GPS tracking, stable isotopes and vessel monitoring systems. J Appl Ecol. 2010;47(2):487–97.
Wakeley JS. Activity budgets, energy expenditures, and energy intakes of nesting ferruginous hawks. Auk. 1978;95(4):667–76.
Williams TM, Wolfe L, Davis T, Kendall T, Richter B, Wang Y, Bryce C, Elkaim GH, Wilmers CC. Instantaneous energetics of puma kills reveal advantage of felid sneak attacks. Science. 2014;346(6205):81–5.
Wilson RP, White CR, Quintana F, Halsey LG, Liebsch N, Martin GR, Butler PJ. Moving towards acceleration for estimates of activity-specific metabolic rate in free-living animals: the case of the cormorant. J Anim Ecol. 2006;75(5):1081–90.
Zuur A, Ieno EN, Walker N, Saveliev AA, Smith GM. Mixed effects models and extensions in ecology with R: Springer Science & Business Media, LLC; 2009.
We thank S. Podofillini, the Provincia di Matera (E. De Capua, N. Braia, A. Pierro), the Provveditorato Interregionale per le Opere Pubbliche (M. Cristallo, A. Cecca), M. Visceglia and E. Fulco for support during field activities. We thank anonymous reviewers and the Associate Editor for their constructive comments that helped us to considerably improve previous manuscript drafts. Meteorological data were obtained from ALSIA (courtesy of E. Scalcione).
Ethics approval and consent to participate
Capture, handling and tagging procedures were carried out by the Italian Institute for Environmental Protection and Research (ISPRA), under the authorization of Law 157/1992 [Art.4(1) and Art 7(5)], which regulates research on wild bird species. Birds were caught by hand at the nexboxes and released back into them within 15 min to minimize stress. They were handled by experienced staff only and no bird was injured by the capturing/handling procedure.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1. Contains the scatterplot of GPS positions in relation to flight velocity and turning angle, highlighting the four behavioural modes assigned by the EMbC algorithm (Figure S1.), the frequency histograms of ODBA values associated to GPS positions for each behavioural mode obtained by the EMbC algorithm, showing that behavioural modes largely differ in ODBA values (Figure S2.), and the ‘caterpillar plot’ illustrating the variation of the random effect estimates obtained by simulations from the final binomial GLMM (Figure S3.).
About this article
Cite this article
Cecere, J.G., De Pascalis, F., Imperio, S. et al. Inter-individual differences in foraging tactics of a colonial raptor: consistency, weather effects, and fitness correlates. Mov Ecol 8, 28 (2020). https://doi.org/10.1186/s40462-020-00206-w
- Dynamic foraging
- Behavioural plasticity
- Birds of prey
- Foraging in flight
- Foraging strategy