Inter-individual variation in provisioning rate, prey size and number, and links to total prey biomass delivered to nestlings in the Collared Flycatcher (Ficedula albicollis)

Background: In bird species where offspring growth and survival rely on parents’ food provisioning, parents can maximise their fitness by increasing the quantity and/or the quality of preys delivered to their offspring. Many studies have focused on inter-individual variation in feeding rate, yet this measure may not accurately reflect the total amount of food (i.e. energy) provided by parents if there is large variation in the quantity and quality of preys at each feeding. Here, we explored the relative role of individual (sex, age, body condition), breeding (hatching date, brood size) and environmental (temperature) factors on feeding rate, prey number, size and quality, and their contribution to total prey biomass delivered to the nestlings of 164 Collared Flycatcher (Ficedula albicollis) parents in 98 nests. Results: Preys delivered to the nest were mainly larvae (53.6%) and flying insects (45.6%). Feeding rate increased with brood size and age, and was higher in males than females. Mean prey number decreased, but mean prey size increased, as the season progressed and parents feeding their brood with primary larvae brought more preys per visit. Relationships between feeding rate, mean prey number and size remained when taking into account the provisioning quality: parents brought either a large number of small prey or a small number of larger items, and the force of the trade-offs between feeding rate and mean prey number and size depended on the quality of the provisioning of the parents. Whatever the percentage of larvae among preys in the provisioning, the variance in total prey biomass was foremost explained by feeding rate (65.1% to 76.6%) compared to mean prey number (16.4% to 26%) and prey size (2.7% to 4%). Conclusions: Our study shows that variation in feeding rate, prey number, size, but not quality (i.e. percentage of larvae), were influenced by individual factors (sex and age) and breeding decisions (brood size and timing of breeding) and that, whatever the provisioning strategy adopted, feeding rate was the best proxy of the total biomass delivered


Background
In birds, provisioning behaviour is strongly associated with reproductive success due to the high-energy demands by young during the period of fast growth (Daunt et al. 2007;Austin et al. 2020). To maximise the energy delivered to their young, parents can use different nonexclusive provisioning behaviour strategies: they can (a) deliver preys at a fast rate, (b) bring many prey items per foraging trip, (c) bring large items, and/or (d) deliver preys with high energy and nutrient contents (Ramsay and Houston 2003;Sexton et al. 2017). Most of the studies on provisioning behaviour so far have focused on the factors accounting for variations in feeding rate (i.e. strategy a) (Christe et al. 1996;Møller and Jennions 2001). The sources of variations in prey number, size and quality (strategies b to d) and their importance in shaping the total amount of energy delivered to the nestlings remain less explored (Wright et al. 1998;Mägi et al. 2009;Bowers et al. 2014;Funghi et al;2019). For example, higher breeding success could be counterintuitively associated to lower feeding rate when chicks are provisioned with higher quality preys and parents need less visits to reach offspring energy demands (Mägi et al. 2009). In this case, provisioning rate alone would not be a good proxy for breeding success. Hence, studies exploring detailed provisioning strategies beyond measures of feeding rate are needed in order to identify the main sources of variation in provisioning behaviour and test the relative contribution of the different provisioning behavioural parameters to the total biomass delivered to the nest. Furthermore, as provisioning strategies and energy fluxes can be influenced by the trade-off between prey quality and quantity (Magrath et al. 2004;Mägi et al. 2009;Espíndola-Hernández et al. 2017), studies investigating the influence of prey quality on quantitative metrics of provisioning rate (such as feeding rate, prey number and size) may help to gain new insights on individual variation in provisioning strategies.
Inter-individual variation in provisioning has frequently been reported in birds. The sources of variation can be broadly divided in at least three categories. First, variation in provisioning can be related to individual factors such as sex, body condition, age and experience. For example, males are often found to feed their young at lower rates than females, which can be driven by their uncertainty in paternity (Sheldon and Ellegren 1998;Gao et al. 2020). Parental effort may also be traded against body condition maintenance due to the trade-off in the allocation of limited energy between current and future reproduction (Verhulst and Nilsson 2008). Individuals may also increase reproductive effort when reaching old ages, to maximise their last reproductive output (Pärt et al. 1992;Froy et al. 2013), or when inexperienced, to reach the chicks' energy demands (Cichoñ 2003). Interindividual differences in cognitive abilities might also shape variation in provisioning behaviour because skills such as searching for, catching and handling preys might have to be learned so that wiser birds may demonstrate a higher foraging efficiency (Rutz et al. 2006;Cauchard et al. 2017;McLeay et al. 2017;Franks and Thorogood 2018). Second, variation in provisioning behaviour can also arise from breeding decisions such as the timing of breeding and the number of chicks to feed. Indeed, individuals starting breeding earlier might secure higher quality breeding sites/territories (Kokko 1999), although breeding in synchrony (i.e. not too early or too late) with the peak of food resources determines the access to food resources (Both et al. 2006;Verhulst and Nilsson 2008). Furthermore, parents of large broods could maximise the energy delivered to the brood by concentrating on caterpillars (Moreno et al. 1995) rather than nutritious (as they contain more of the amino acid taurine important of nestling development; Ramsay and Houston 2003) but low energy contents spiders (Wiebe and Slagsvold 2015;Espíndola-Hernández et al. 2017). Finally, environmental factors such as temperature or rainfall are likely to influence provisioning by affecting both parent's ability to forage (Sanz and Tinbergen 1999;McKechnie and Wolf 2010;Wiley and Ridley 2016) and prey availability (Funghi et al. 2019;Vincens and Bosch 2000).
In this study, we explored the relative role of individual factors (sex, age, body condition), breeding decisions (timing of the season, brood size) and environmental factors (temperature) on feeding rate, number, size and type of preys delivered to the nestlings in a wild population of Collared Flycatchers (Ficedula albicollis). Collared Flycatchers typically prey on flying insects on the wing, but during the rearing period, they also forage on the foliage and on the ground for caterpillars and other crawling insects (Gustafsson 1989), which both have a higher energy content than flying insects and require less energy to capture (Barba et al. 1996). Because provisioning parameters such as feeding rate, number and size of preys can differ across individuals based on the quality of their provisioning (e.g. percentage of larvae among preys, Espíndola-Hernández et al. 2017), we first investigated the individual, breeding and environmental factors influencing the percentage of larvae among preys. Then, we controlled for this percentage of larvae when analysing individual variation in feeding rate, number and size of preys depending on the same individual, breeding and environmental, factors. This approach allowed us to examine to what extent feeding rate alone, which is widely used in the literature as a non-invasive proxy of food intake, provides a reliable measure of the quality and quantity of the prey biomass delivered to the nestlings.

Study site and population monitoring
The study was conducted on 164 (95 females and 69 males, corresponding to 98 nests) Collared Flycatchers, breeding on the island of Gotland, Sweden (57° 07′ N, 18° 20′ E) between April and July 2015. This patchy population of small, migratory hole-nesting passerine birds has been monitored each year since 1980. Flycatchers readily accept to breed in artificial nest boxes, providing an easy access to recording breeding data and behaviour at the nest (Doligez et al. 2004). Collared Flycatchers display biparental care and feed on average 3 to 5 young from late May to early July during 16 days until fledging.
Nests were visited at least every four days to collect data on laying and hatching dates, clutch and brood size (at day 8, 2-3 days after provisioning videos), and final breeding success. When chicks were 8 to 14 days old, adults were caught and identified using individually numbered rings. They were sexed (based on plumage sexual dimorphism (i.e. males are black and white while females are brownish) and aged (yearling vs. 2 years old or older, based on plumage characteristics, Svensson 1992), and measured (tarsus length, to the nearest 0.1 mm, and body mass, to the nearest 0.1 g). Body condition was computed as the body mass divided by tarsus length (Doligez et al. 2004).
Data on daily temperature were obtained from the weather centre 02679 in Hoburgen, Gotland (56° 55′ N, 18° 09′ E; 36 m), which is located around 20 km from our field site. We used the mean daily temperature in Celsius degrees.

Food provisioning
Food provisioning behaviour was recorded during 1h20 between 06:00 and 13:00 (previous study in the same population showed no variation in the provisioning rate using the same temporal window; Sheldon et al. 1997) when nestlings were 5-6 days old, i.e. during the peak of daily parental activity and nestling growth (Krist et al. 2004;Rosivall et al. 2009). An IR camera was installed in the nest box on the day before the recording to let birds habituate to the device. Videos to be analysed were randomly attributed between two observers (EIM and RL) blind to the individual and environmental variables subsequently used in statistical models. Using the motion-detection software programme Motion Meerkat (Weinstein 2015), we recorded for each parental visit the following parameters: (1) entrance and exit time, (2) number of preys delivered, (3) prey categories divided into 4 groups: caterpillar, flying insect, spider or undetermined, (4) prey size on a scale of 1-4, with 1 for preys smaller than beak length, 2 for preys of about beak length, 3 for preys twice longer than beak length and 4 for preys three times longer than beak length or bigger.
For each individual, we extracted five descriptors of provisioning behaviour: (a) mean number of prey per visit, (b) mean prey size per visit, (c) feeding rate per hour, computed as the total number of feeding visits that occurred from the time the bird first appeared on the video until the end of the video divided by the time interval in hours, (d) total prey biomass delivered to the nestlings per hour, computed as mean number of prey per visit × mean size of prey per visit × feeding rate per hour, and e) the percentage of larvae over the total number of preys delivered to the nestlings per hour.
Preliminary analyses of 147 visits by the two observers showed that they differ in how they scored prey number and size per visit: one observer counted more preys (paired t-test: t 146 = -4.00, P < 0.001; mean [95%CI] difference = 0.35 [0.18; 0.53]) but sized preys as smaller (t 146 = 8.69, P < 0.001; -0.50 [-0.61; -0.38]) than the other observer. There were no differences in the prey category. Hence, we standardized mean prey size and mean prey number using mean differences between observers as correction factors for subsequent statistical analyses.

Statistical analyses
The quality of the food delivered to the nestlings by a parent can influence his strategy (i.e. provisioning frequency and number/size/type of items per visit). We thus first examined the variation of the proportion of larvae delivered to the nestlings and controlled for it when further investigating the individual and environmental factors on provisioning parameters.
We investigated the influence of individual and environmental factors on provisioning behaviour by running three sets of generalised linear mixed effect models with (1) mean prey number per visit, (2) mean prey size per visit and (3) feeding rate. We included as random effects the forest patch to account for potential spatial variation in prey availability and nest box to account for the non-independence of both parents in a pair. We included sex, age, body condition, hatching date, brood size at day 8, mean daily temperature and the percentage of larvae as fixed variables, plus all possible two-way interactions between sex or age and the other fixed effects. Interactions were removed from the models if non-significant.
We then characterized food provisioning strategies by assessing partial correlations between mean prey number and size per visit and provisioning rate. To outline differences when looking for provisioning quantity or quality, we divided parents into three categories depending on whether they were bringing few larvae in the preys, a mix of larvae and other insects, or mostly larvae. We used tertiles in the percentage of larvae in the preys to build up our three categories.
Finally, we assessed the relative contributions of mean prey number and size per visit and feeding rate to the total prey biomass delivered to the nestlings per hour using used multiple linear regression models. We split the analyses in parents provisioning their nest with larvae in the preys, a mix of larvae and other insects, or mostly larvae, to investigate whether the type of preys brought influenced the relationships between provisioning variables and total prey biomass delivered to the nestlings.
Data analyses were carried out in the version 3.5.0 of the statistical freeware R cran (R Core Team 2016). Continuous explanatory variables were scaled before analysis to improve the interpretability of our model estimates (Schielzeth 2010). Mean estimates and 95% confidence intervals are reported in the results; intervals that do not overlap with 0 are considered as significant.

Results
Provisioning behaviour was recorded during 80 min on 98 nests, corresponding to 164 parents actively provisioning their chicks (95 females and 69 males). The main prey items delivered to the nestlings were larvae (53.6%; 2033 out of 3795 preys) followed by flying insects (45.6%; 1731 preys). Spiders and preys that could not be classified represented less than 1% of the preys (23 spiders and 8 non-classified).

Factors influencing the proportion of larvae delivered to the nestlings
Although the proportion of larvae delivered to the nestlings greatly varied across individuals (Fig. 1), none of the individual, breeding and environmental factors considered in this study significantly explained variation in this trait (all P-values > 0.1; Table 1).

Relationships between provisioning variables
We examined the relationships between mean prey size, number and feeding rate while taking into account for the quality of prey items provisioned by the parents (i.e. percentage of larvae). The first and second tertile in the distribution of the percentage of larvae brough to the nest were at 44.4% and 66.7% (Fig. 1). Thus, we compared the relationships between provisioning variables in parents provisioning their chicks with < 44.4% of larvae in the preys (lower tertile; n = 54 parents) versus parents provisioning their chicks between 44.4 and 66.7% of larvae in the preys (middle tertile; n = 53 parents) versus parents provisioning their chicks with > 66.7% of larvae in the preys (upper tertile; n = 57 parents).
Results show a strong trade-off between the mean prey number and mean prey size per visit for a given parent independently of the percentage of larvae in the preys (Table 2). Thus, for each visit, parents provisioned either large numbers of small prey items or small numbers of large items (Fig. 5a). There was evidence for a significant trade-off between feeding rate and mean number of preys brought per visit in parents bringing a mixed of larvae and other insects to their offspring (Table 2), those parents making either fewer visits with many prey items per visit or many visits with few prey items per visit (Fig. 5b). Finally, there was also evidence for a trade-off between feeding rate and mean prey size foremost apparent in parents bringing mostly larvae among the preys (Table 2), those parents making either few visits with large prey items per visit or many visits with small prey items per visit (Fig. 5c).

Provisioning variables explaining the total prey biomass delivered to the nestlings
Whatever the percentage of larvae delivered to the nestlings by a parent, the variance in total prey biomass per hour delivered to the nestlings was mainly explained by variation in feeding rate (65.1% to 67.6%; Table 3): parents provisioning their chicks with few larvae, a mixed percentage of larvae, or many larvae, increased their total prey biomass delivered to the nestlings while increasing their feeding rate ( . The variance in total prey biomass was also significantly explained by mean prey number and size delivered to the nestlings (Table 3), but to a lower extent, with prey number explaining 16.4% to 26% of the variance and prey size explaining 2.7% to 4.0% of the variance. Parents provisioning their chicks with few larvae, a mixed percentage of larvae, or many larvae, increased their total prey biomass delivered to the nestlings while increasing the number preys brought per feeding visit

Discussion
In this study, we examined to what extent individual, breeding and environmental factors influence parental provisioning strategies, and investigated the relative contribution of feeding rate, prey number and size in influencing the total prey biomass delivered to the nestlings. We found collared flycatcher parents to provision their chicks mainly with larvae (53.6%) and flying insects (45.6%). We found that neither sex, age, body condition, hatching date, brood size or temperature influenced the percentage of larvae among the preys, even if we found a great inter-individual variation in this percentage (Fig. 1). This result could seem surprising, since previous studies found flycatcher's diet to vary with time or sex. Parent Collared Flycatchers were reported to bring mostly caterpillars over other preys when caterpillar abundance is high around the second half of May (Török and Tȯth 1988). In the closely related Pied Flycatcher (Ficedula hypoleuca), one study reported that females foraged on caterpillars in the tree canopy and males on flying insects caught the wings (Mänd et al. 2013), whereas no differences between sexes in prey types were found in two other studies (Moreno et al. 1995;Siikamäki et al. 1998). Interestingly, the sex difference in foraging niche reported by Mänd et al. (2013) was reduced when environmental conditions became worse (i.e. experimentally increased level of hunger of nestlings). Since the environmental conditions experienced by our population in spring 2015 were relatively unfavourable (i.e. reproductive success in 2015 was lower than in other years, pers. obs.), we cannot exclude that we could have found sex and timing of breeding differences if the study was performed under different conditions. Individuals might also learn specific prey choices throughout their life, so that older birds may select more profitable prey types for nestlings (Forslund and Pärt 1995;Wiebe and Slagsvold 2015). However, again, we did not observe any difference with age in our study although a large sample size on very old birds (4 years old and older) would be needed to look at senescence in provisioning parameters in this species.
As previously found in the literature, feeding rate increased mainly with brood size: the more chicks there are the more food parents have to bring back to the nest (Sanz and Tinbergen 1999), including for experimentally enlarged broods (Sendecka et al. 2007). We also found feeding rate to increase with age, in line with previous studies in Collared Flycatchers and other bird species (Pärt et al. 1992;Cichoñ 2003;Froy et al. 2013). This is either because birds may increase reproductive effort when reaching old ages to maximise their last reproductive output (Froy et al. 2013), or because young birds may show a reduced experience in provisioning (Cichoñ 2003) or in securing high quality territories/sites with high quality prey (Slotow and Rothstein 1991;Franks and Thorogood 2018). Feeding rate also varied with sex, but we found males to provision their young more often than females. This result can be surprising, as the reverse is usually found in the literature in species with biparental care (review in Lewis et al. 2002 for seabirds and Kilner 2002 for passerine); however, no difference is often reported too (Moreno et al. 1995;Peck and Congdon 2006). Differences in feeding rate between sexes can arise if males engage in activities such as nest/territory defence against competitors at the expense of parental care (Qvarnström 1997) or males reduce their parental effort when being uncertain about their paternity in species with extra-pair mating, such as the Collared Flycatcher (Gao et al. 2020). Our result is likely due to the timing of our video recordings: we recorded provisioning early in the morning during the nestling rearing period when females, on top of provisioning, still partly brood nestlings that do not fully thermoregulate yet, while males only provision nestlings (Schwagmeyer et al. 1999). Accordingly, on our videos, males were observed from time to time delivering the prey items that they brought back to the brooding female inside the nest box, and the female then delivered the food to the young. This behaviour both increased the feeding rate of males that were saving time on each visit and decreased the feeding rate of females that were spending more time in the nest box. In our study, the feeding rate was not influenced by the timing of the breeding (i.e. hatching date) nor by the percentage of larvae among the preys. However, our results showed that mean prey size and number per visit varied with the timing of the breeding, as previously found in the literature in various species too (Bowers et al. 2014), and the percentage of larvae among the preys: parents that fed their young mostly with larvae (i.e. > 77%) brought more preys per visit. Although, we cannot exclude that using only prey length, and thus ignoring their width, height as well as their mass, may have led to a low accurate estimate of prey size.
We then compared the relationship between feeding rate, mean prey size and number per visit, taking into account the percentage of larvae among the preys brought by the parent. Parents with low, medium and high percentage of larvae in their provisioning showed negative correlation between the mean prey number and size per visit, as previously found in the literature (Siikamäki et al. 1998;Wright et al. 1998). However, the feeding rate itself increased when mean prey number decreased, a trade-off that was foremost apparent in parents with a mixed percentage of larvae during provisioning: these parents made either few visits with many prey items per visit or many visits with few prey items per visit. In parents provisioning their brood mainly with larvae or flying insects, the relationships between feeding rate and prey number were negative, but not significantly so. We also found feeding rate to increase when mean prey size decrease, but mainly in parents with a high percentage of larvae during provisioning: these parents made either few visits with large prey items per visit or many visits with small prey items per visit. These results support the hypothesis that quantitative measures of provisioning strategies (i.e. feeding rate, prey number and size) can be influenced by variation in prey quality, as previously reported in the three other insectivorous birds: the Brown Songlark (Cincloramphus cruralis) (Magrath et al. 2004), the Great Tit (Parus major) (Mägi et al. 2009) and the Thorn-tailed Rayadito (Aphrastura spinicauda) (Espíndola-Hernández et al. 2017). In the present study, prey quality had nonetheless only minor consequences on variation in provisioning strategies. This might be partly explained by the crude measure of quality (larvae versus flying) that we used, and ideally information per type of prey on energy content as well as search and handling time should be used to build an index of prey quality (Barba et al. 1996;Mägi et al. 2009;Espíndola-Hernández et al. 2017).
Finally, whatever the percentage of larvae among the preys delivered to the nestlings, the total prey biomass was mainly influenced by the feeding rate: 70.1%, 76.6% and 65.1% of the total prey biomass variance was explained by the feeding rate of parents provisioning their chicks with law, mixed or many larvae respectively. Mean prey number explained a lesser part of this variation, explaining 17.4%, 16.4% and 26.0% respectively. As for mean prey size, it only explained a negligible part of the total biomass, i.e. 2.7%, 4.0% and 3.5% of the variance. Although nestlings' differences in food requirements and in begging behaviour due to sexual differences might influence inter-individual chicks' difference in food intake (Espíndola-Hernández et al. 2017), these results support the use of feeding rate as a good non-invasive proxy of nestlings total food intake during rearing. Our estimate of biomass is based however on estimates of size from pictures, and future studies are needed to validate this crude estimate as an accurate measure of biomass in g per se.
Besides individual and breeding factors, environmental factors such as meteorological conditions should influence provisioning performance by affecting both parent's ability to forage and prey availability (Funghi et al. 2019). For example, poor environmental conditions, as reflected by low temperature, might affect birds' abilities to locate preferred prey types and/or increase the energetic cost of looking for prey items. Indeed, low temperatures reduce ectothermic preys' metabolism and thus growth and activity for caterpillars, flying insects and spiders (Vincens and Bosch 2000). Moreover, foraging in colder weather is more energy demanding for parents for maintaining both body temperature and activity level (Sanz and Tinbergen 1999), so that birds may not be able to invest as much energy in provisioning compared to warmer days. In the same way, high temperature can cause a stress to breeders that rapidly exceed their physiological tolerance limits (McKechnie and Wolf, 2010) and affect their investment in young (Wiley and Ridley 2016). In our study, environmental factors showed only a small, non-significant, contribution to variation in provisioning traits: only the mean number of preys tended to be affected by temperature, with parents bringing more prey items per visit on warmer days. During our study, daily temperature ranged from 11.6 °C to 14.8 °C, and it is thus possible that variance was too small in our sample to detect an effect. We were not able to quantify the influence of rainfall because provisioning behaviour was not recorded during heavy rain days, to avoid the risk of further disturbance for the parents known to be in difficult conditions already (our personal observations showed that provisioning was strongly reduced in periods of heavy rain).

Conclusions
We investigated simultaneously individual, breeding and environmental factors to better understand their effect on Collared Flycatchers provisioning performance. While brood size, sex and age explained variation in feeding rate, only hatching date and the quality of the provisioning explained variation in mean prey number and size per visit. Despite the trade-off in provisioning strategies between prey number and size (i.e. parents make either few visits with many or large prey items per visit, or many visits with few or small prey items per visit), and between feeding rate and mean prey number/size under certain conditions, the total prey biomass brought per hour to the nest was mainly influenced by feeding rate. This result supports the use of feeding rate as an important measure of provisioning in the avian literature.