Behavioural and energetic consequences of competition among three overwintering swan (Cygnus spp.) species

Winter numbers of the northwest European population of Bewick’s Swans (Cygnus columbianus bewickii) declined recently by c. 40%. During the same period, numbers of two sympatric and ecologically-similar congeners, the Mute Swan (Cygnus olor) and Whooper Swan (Cygnus cygnus) showed increases or stability. It has been suggested that these opposing population trends could have a causal relationship, as Mute and Whooper Swans are larger and competitively dominant to Bewick’s Swans in foraging situations. If so, effects of competition of Mute and Whooper Swans on Bewick’s Swans should be detectable as measurable impacts on behaviour and energetics. Here, we studied the diurnal behaviour and energetics of 1083 focal adults and first-winter juveniles (“cygnets”) of the three swan species on their winter grounds in eastern England. We analysed video recordings to derive time-activity budgets and these, together with estimates of energy gain and expenditure, were analysed to determine whether individual Bewick’s Swans altered the time spent on key behaviours when sharing feeding habitat with other swan species, and any consequences for their energy expenditure and net energy gain. All three swan species spent a small proportion of their total time (0.011) on aggressive interactions, and these were predominantly intraspecific (≥ 0.714). Mixed-effects models indicated that sharing feeding habitat with higher densities of Mute and Whooper Swans increased the likelihood of engaging in aggression for cygnet Bewick’s Swans, but not for adults. Higher levels of interspecific competition decreased the time spent by Bewick’s Swan cygnets on foraging, whilst adults showed the opposite pattern. When among low densities of conspecifics (< c. 200 individuals/km2), individual Bewick’s Swans spent more time on vigilance in the presence of higher densities of Mute and Whooper Swans, whilst individuals within higher density Bewick’s Swan flocks showed the opposite pattern. Crucially, we found no evidence that greater numbers of interspecific competitors affected the net energy gain of either adult or cygnet Bewick’s Swans. We found no evidence that Bewick’s Swan net energy gain was affected by sharing agricultural feeding habitat with larger congeners during winter. This was despite some impacts on the aggression, foraging and vigilance behaviours of Bewick’s Swans, especially among cygnets. It is unlikely therefore that competition between Bewick’s Swans and either Mute or Whooper Swans at arable sites in winter has contributed to the observed decline in Bewick’s Swan numbers. Further research is needed, however, to test for competition in other parts of the flyway, including migratory stopover sites and breeding areas.


Background
The conservation of avian abundance and diversity represents a major global endeavour, requiring accurate information on the causes of population declines (Kirby et al. 2008). In recent decades, the northwest European population of Bewick's Swans (Cygnus columbianus bewickii) has become a growing conservation concern. Repeated censuses have revealed that numbers wintering in northwest Europe have declined by 39.4% from a peak of 29,780 individuals in January 1995, to 18,057 in January 2010, with only a slight recovery to 20,149 in January 2015, and as a result the Bewick's Swan has been classified as Endangered in Europe (BirdLife International 2015; Beekman et al. 2019). To date, the demographic and environmental causes of this population decline have not been resolved fully (Wood et al. 2018;Nuijten et al. 2020a). To help inform conservation efforts, there is a need for further research to assess the impacts that threats to Bewick's Swans might have on swan numbers.
The observed decline in Bewick's Swan numbers was in contrast to many other populations of large herbivorous waterbirds in Europe, which have shown stable or increasing population trends since the Bewick's Swan population peaked in the mid-1990s (Rees et al. 2019). These include two congeners of the Bewick's Swan that are also native to Europe: the Mute Swan (Cygnus olor), which is largely sedentary throughout northwest Europe, and the two populations of migratory Whooper Swan (Cygnus cygnus), one of which breeds in Iceland and winters in the UK and Ireland, while the other breeds in northwest Russia and Fennoscandia, and winters in mainland northwest Europe (Rees et al. 2019). Mute Swan numbers across Europe have shown regional variation, but in almost all regions have either increased or have remained stable since the mid-1990s (Rees et al. 2019;Wood et al. 2019a). Whooper Swan numbers have more than doubled since 1995 (Hall et al. 2016;Laubek et al. 2019).
Both Whooper Swans and Mute Swans show considerable overlap with Bewick's Swans in their use of winter habitat in northwest Europe, with all three species feeding on agricultural crops and frequently associating in mixed species flocks (Rees et al. 1997). As the three swan species exploit the same food resources within the same habitats, the potential for competition amongst the three species exists. Observations have indicated that the larger Mute and Whooper Swans are dominant over the smaller Bewick's Swans in behavioural interactions (Black and Rees 1984;. In a study of all three swan species feeding on grain provided at a lake during winter, Black and Rees (1984) reported that Whooper Swans could displace both Bewick's and Mute Swans from feeding locations. Mute and Whooper Swans could therefore make winter feeding areas less attractive to Bewick's Swans through interference and depletion competition which reduced their rates of energy gain. There are potential fitness consequences to such competition, as for capital breeders such as swans, the energy and body reserves required to breed successfully is gained during the preceding winter (Drent et al. 2006) and at staging sites during migration (Nolet and Drent 1998). Competition at feeding sites therefore has been suggested as a possible contributory factor in the decline of the Bewick's Swan population (Nagy et al. 2012;Ponting 2014), but there have been no empirical tests of this hypothesis since the onset of the decline in Bewick's Swan numbers. Where the three swan species associate in mixed species flocks, any such effects of competition of Mute and Whooper Swans on Bewick's Swans should be detectable as measurable differences in behaviour (e.g. more time spent in aggressive interactions) and energetics (e.g. reduced net energy gain) (Amat 1990). Such an assessment would need to account for variables which might confound any effects of interspecific competition, such as the density of intraspecific competitors; a study of interference competition among Bewick's Swans by Gyimesi et al. (2010) found that individuals suffered progressively reduced food intake rates at greater competitor densities. It is not known, however, whether Bewick's Swans show similar behavioural and energetic responses to rising conspecific and heterospecific competitor densities, or whether the species identity of the surrounding individuals influences those responses.
In this study we assessed the behavioural and energetic consequences for the three swan species where they shared diurnal feeding habitat. Our study had two key objectives. Given that empirical data on interactions amongst overwintering swans are limited, our first objective was to quantify the time spent engaged in aggressive behaviours with conspecifics and heterospecifics. Our second objective was to test whether the behaviour or energetics of swans were altered by the presence of other swan species.
Keywords: Avian behaviour, Energy expenditure, Ethology, Interference competition, Interspecific interactions, Time activity budgets, Video observations

Study system
Agricultural land in lowland regions of northwest Europe represents important winter habitat for Bewick's, Mute, and Whooper Swans (Owen and Cadbury 1975;Rees et al. 1997;Davis et al. 2014). Our study focused on the area surrounding the Ouse Washes (52°31′ N, 0°16′ E) in eastern England; an expanse of former wetland, now drained for agriculture, which features a mosaic of predominantly arable fields and drainage ditches (Nisbet 1955;Thomas et al. 1981). This area supports internationally important numbers of Bewick's and Whooper Swans, as well as nationally important numbers of Mute Swans (Frost et al. 2019 (Fig. 1). Mute Swans have also shown a recent localised decline in numbers from 1151 in 2007/08 to 287 in 2017/18 ( Fig. 1), although their national numbers within Britain were stable over this period (Wood et al. 2019a).
Historically, swan species spent some time feeding on flooded grasslands within the Ouse Washes themselves (Nisbet 1955;Owen and Cadbury 1975), but in recent decades the swans have switched to feed on arable crops almost exclusively, with < 1% of individuals observed using pastures (Wood et al. 2019c). Overwintering swans feed on arable crops in the fields during the day, and return to the seasonally flooded Ouse Washes to roost (Ponting 2014;Wood et al. 2019c). All three swan species excavate and consume partially-buried post-harvest remains of root crops including Sugar Beet (Beta vulgaris) and Potatoes (Solanum tuberosum), along with some Maize (Zea mays), which are sown in spring for harvest during autumn and early winter (Rees 2006;Wood et al. 2019c). All three swan species graze on the above-ground biomass of cereals and brassicas such as Wheat (Triticum aestivum) or Oilseed Rape (Brassica napus), which are sown in the fields during autumn and early winter (Rees 2006;Wood et al. 2019c).

Data collection
Video footage of adults and cygnets of all three swan species was obtained during October-March inclusive between winters 2015/16 and 2017/18 (Additional file 1: Table S1). Each month, an area within 15 km of the Ouse Washes was surveyed for swans over a three-day period, based on the area known to have been used by the swans in previous winters (Ponting 2014;Wood et al. 2019c). The numbers of swans using the study area over a threeday period was considered to be relatively stable, especially given its position on the far edge of the flyways of both Bewick's and Whooper Swans (Brazil 2003;Rees 2006). Moreover, different swan species have been shown to exhibit the same seasonal pattern of habitat use, with post-harvest remains of root crops preferred in early winter, and early-growth cereals and brassicas preferred in late winter (Wood et al. 2019b). We attempted to achieve representative coverage of age classes, species, and swan densities found within the landscape (Additional file 1: Table S2). As all three swan species regularly intermix, all swans within the same field were considered to be sharing feeding habitat. The densities of (i) Bewick's Swans and (ii) their interspecific competitors in each field (individuals/km 2 ) were calculated as the total numbers of (i) Bewick's Swans and (ii) Mute and Whooper Swans combined, divided by the total field size. Where swans were encountered during surveys, a random subset of the individuals present were filmed using a tripod-mounted Canon SX50 HS video camera (Canon Inc., Japan) for up to 20 min, following the focal-bird sampling methods used in previous studies of swan behaviour (e.g. Keane and O'Halloran 1992;O'Hare et al. 2007;Tatu et al. 2007;Wood et al. 2019c). All individuals filmed within a given video were assigned an alphanumeric code (hereafter termed observation identity) that was unique to that Fig. 1 The long-term trends in the peak numbers of Bewick's, Mute, and Whooper Swans counted on the Ouse Washes as part of the Wetland Bird Survey (WeBS) monitoring programme (Frost et al. 2019) video, such that swans that were sharing feeding habitat at the same time shared the same observation identity. We measured the air temperature (± 0.5 °C) at the start of each video recording using an immersion thermometer (Breaksafe Thermometer, Brannan, UK). In this study we only considered diurnal behaviour and energy expenditure, as Bewick's Swans and the other two species typically roost in different areas within the Ouse Washes (Ponting 2014), and so competition was considered to only be likely during daylight hours.

Time-activity budgets
We used the footage obtained to construct time-activity budgets (Paulus 1988) for each individual using the focalindividual sampling method (Altmann 1974). Analysis of swan behaviours was performed using the Behavioral Observation Research Interactive Software (BORIS) developed by Friard and Gamba (2016), an event recorder that permits frame-by-frame viewing. For each video, all individuals within the video frame were watched separately, with the duration of each behaviour type noted for the duration that the focal individual remained visible. Eight behaviour categories, used in previous studies to synthesise the main activities of non-breeding swans (Nolet et al. 2002;Rees 2006;Włodarczyk 2017), were used to describe the swans' behaviour: resting, foraging, vigilance, preening (including other maintenance behaviours such as stretching), walking, flying, aggressive interactions, and non-aggressive social interactions. An aggressive interaction was defined as any interaction between the focal individual and one or more other individuals which involved giving or receiving physical or vocal threat displays or attacks . We also noted whether each aggressive interaction was intraspecific or interspecific (including which species the interspecific interactions were with).

Energetics
We calculated the rate of net energy gain of each observed individual by estimating their expenditure and gain of energy during their observation period. To estimate energy expenditure, we combined time-budgets with activity-specific multipliers of the basal metabolic rate (BMR) for that species and age class (Krivtsov and Mineyev 1991;Clausen et al. 2013). The basal metabolic rate for each age class in each species was calculated from the allometric equation presented by Hughes and Green (2005): where z represents the BMR value at 1 kg body mass (reported as 4.64 Watts for Anseriformes by Zar 1968) and M was the mean body mass of the swan in kilograms (averaged for males and females; Kear 2005) (Additional file 1: Table S3).
The baseline activity-specific energetic costs (Watts) for adult Bewick's Swans were those measured experimentally by Nolet et al. (2002). We used these values to determine the multipliers of basal metabolic rate (xBMR) for each behaviour for an adult Bewick's Swan, which could then be used to estimate the energetic costs of each behaviour for all of the other swan species and age classes, assuming that the costs were proportionally equivalent (Additional file 1: Table S4). Energy expenditure for each focal individual during an observation period could thus be calculated by summing the energy expenditure spent on each behavioural activity.
To estimate the gross energy gain of swans feeding on cereal and root crops, we used the equations derived by Wood et al. (2019b), as these accurately predicted resource selection by overwintering swans within our study area. The gross energy gain (E Gain ; Joules per observation) of each observed individual was estimated as: where I D was the species and age-class-specific instantaneous intake rate on a given density of either cereal or root crops (g Dry Matter/s; Wood et al. 2019b), Q was the mean energy density (J/DM g; Wood et al. 2019b) of the crop type, A was the proportional assimilation efficiency for a given crop type, whilst F Time was the total time (in seconds) that the focal swan was observed to spend on foraging. For the parameters I D and Q, mean monthly values for cereal and root biomass and energy densities were available from Wood et al. (2019b). For the parameter A, temporally-invariant proportional assimilation efficiency values of 0.63 and 0.84 were used for cereal and root crops, respectively (Nolet et al. 2002;Amano et al. 2004).
We also accounted for additional energetic costs of thermoregulation during cold or warm weather (Bech 1980;Jenssen et al. 1989). Bech (1980) found that Mute Swans had lower and upper critical temperatures of 0.8 °C and 15.0 °C, respectively, and hence we added an additional thermoregulatory cost to each observation made where the air temperature was outside of the range 0.8-15.0 °C. We applied this additional energetic cost by first estimating p BMR, which represented the BMR of the individual outside of its thermal neutral zone (BMR c ) as a proportion of its BMR within its thermal neutral zone: where V a was the instantaneous VO 2 per unit mass reported for Mute Swans within the thermal neutral zone (0.654 mL O 2 /(g·s)), m was the mean body mass of the swan in grams (averaged for males and females; Kear 2005) (Additional file 1: Table S3), and e was the energy yielded per ml of oxygen consumed (20 J/mL O 2 ; Nolet et al. 2002). Similarly, V c was the instantaneous VO 2 per unit mass reported for Mute Swans outside of its thermal neutral zone, which was reported by Bech (1980) for both the lower and upper critical temperatures of 0.8 °C and 15.0 °C, respectively, as: where t was the air temperature (°C). The air temperatures measured during each observation could therefore be used to adjust the baseline energetic costs of each behavioural activity.
We estimated the rate of net energy gain of each focal swan during an observation as the proportion of time the swan spent foraging multiplied by the expected instantaneous energy gain, minus energy expenditure. To facilitate comparisons across species and age classes, we converted our estimates of net energy gain from Watts to a multiple of BMR (hereafter: xBMR) for each individual. The estimates for Bewick's Swans were used in our statistical models (as described below), with the values derived for Mute and Whooper Swans reported for comparative purposes.

Statistical analyses
All statistical analyses were carried out using R version 3.6.3 (R Core Team 2020). To address our first objective regarding aggressive interactions among swans, we used two-sample binomial tests for equality of proportions (Crawley 2013) to assess the significance of the differences in the proportion of intraspecific and interspecific aggressive interactions recorded. Separate tests were conducted for each species and age class, with significant differences attributed where P < 0.05.
To address our second objective regarding impacts of competition on the behaviour and energetics of individual Bewick's Swans, we ran separate sets of models for each of our four response variables: (i) the probability of aggression observed during an observation; (ii) the proportion of time spent on foraging during an observation; (iii) the proportion of time spent on vigilance during an observation; and (iv) the rate of net energy gain during an observation. We ran linear mixed effects models, using the glmmTMB R package (Brooks et al. 2017). For each of our response variables, our set of candidate models consisted of all additive and two-way interactive combinations (that did not contain correlated variables, see below) of (i) the density of Bewick's Swans sharing feeding habitat with the focal individual (individuals/km 2 ); (ii) the combined density of Mute and Whooper Swans sharing feeding habitat with the focal individual (individuals/ km 2 ); and (iii) the age class of the focal individual (cygnet or adult). To account for non-independence of individuals within the same flock, we modelled observation identity as a random intercept in each model. Within each model, each observation was weighted by the duration of the observation period (min), to reflect that observation periods were not equal among focal individuals, and that longer observations potentially contained more information than shorter observation periods. For models of foraging, vigilance, and net energy gain, these linear mixed effects models (LMMs) featured Gaussian error structures. In all models of foraging and vigilance, the proportions of time spent on these behaviours were logit-transformed using the car package (Fox et al. 2018) to ensure that model assumptions were met (Warton and Hui 2011); similarly, net energy gain values (xBMR) were rescaled between 0 and 1 using the scales package (Wickham and Seidel 2020), and then logit-transformed using the same approach. In contrast, we used generalized linear mixed effects models (GLMMs), with binomial error structures, for the probability of aggression. Earlier attempts to use zero-inflated generalised linear mixed models (ZIGLMMS) for the time spent on aggressive interactions (based on the approach of ) resulted in models which failed to converge, so instead we opted to model the probability of occurrence for aggressive interactions. Model assumptions were checked using the performance package (Lüdecke et al. 2020). In addition, we calculated the variance inflation factor (VIF) associated with each explanatory variable, using the performance package (Lüdecke et al. 2020). We found that interactions between focal age class and Bewick's Swan density had VIF values of < 5.0 (Dormann et al. 2013), and so models containing this interaction were removed from our sets of candidate models for each of the four response variables, and were not considered further, so that multicollinearity did not influence the efficiency and reliability of parameter estimation for our models (Dormann et al. 2013).
For each of our response variables, we assessed the performance of each non-collinear candidate model using second-order Akaike's Information Criterion (AICc) values (Burnham et al. 2011), with any model with a ΔAICc value of ≤ 6.0 considered to have received support in the data (Richards et al. 2011). To be considered competitive, a model also had to perform better (AICc ≥ 6.0) than the null model, i.e. a model comprised of an intercept and no other parameters. In addition to the AICc threshold, we also took an additional step to minimize the effects of uninformative parameters (sensu; Arnold 2010) during model comparisons; models were judged competitive only if the AICc value associated with a more complex form of the model (i.e. one containing additional parameters) was lower than the simpler model (Harrison et al. 2018). We assessed the proportion of the variance in each response variable explained by a candidate model based on the two values of adjusted R 2 proposed for mixed effects models by Nakagawa et al. (2017), calculated using the performance package (Lüdecke et al. 2020). The marginal R 2 indicated the proportion of variance explained by the fixed effects alone, whilst the conditional R 2 indicated the proportion of variance explained by both the fixed and random effects.

Aggression
Across all species and age classes, the mean (± 95% CI) proportion of time spent by focal swans engaged in aggressive interactions was 0.011 ± 0.002 (Table 1). Intraspecific interactions accounted for 0.010 ± 0.002 of the total time budgets, whilst interspecific interactions represented 0.010 ± 0.001 of the total time for which swans were observed (Additional file 1: Table S5). Thus, across all species and age classes, intraspecific and interspecific aggression accounted for 0.872 and 0.128 of the total time spent on aggressive behaviours, respectively. The proportion of aggressive interactions that involved conspecifics was significantly greater than the proportion involving heterospecifics for all swans except Mute Swan cygnets (Table 2).

Foraging
Based on our focal observations, the mean (± 95% CI) proportion of observed time spent on foraging by Bewick's Swans was 0.414 (± 0.003) for adults and 0.420 (± 0.049) for cygnets (Table 1). Whilst the foraging effort observed for adult Bewick's Swans was the lowest observed for any species and age class, the greatest foraging effort was 0.573 (± 0.098) by cygnet Mute Swans (Table 1). These Table 3 A summary of our best-supported models (ΔAICc < 6.0) and null models of Bewick's Swan behavioural activities and net energy gain (NEG) For each behaviour the best-supported model is indicated in bold. k represents the number of parameters in the model. R 2 m and R 2 c indicate the proportions of variance explained by the fixed effects alone and both the fixed and random effects, respectively. RL, wi, and ER represent the relative likelihood, Akaike weight, and evidence ratio, respectively. Model parameters: intercept (1), density of Bewick's Swans sharing feeding habitat with the focal individual (B), density of Mute and Whooper Swans sharing feeding habitat with the focal individual (C), age class of the focal individual (A), and observation identity (U). The symbols + and: represented additive and interactive effects, respectively, within models A comparison of our candidate models indicated that the proportion of time spent by Bewick's Swans on foraging behaviour was best explained by a model comprised of the age class of the focal individual, the density of other Bewick's Swans sharing the feeding habitat, and the combined density of Mute and Whooper Swans sharing the feeding habitat with the focal individual. This bestsupported model accounted for 35.4% of the total Akaike weights, whilst the fixed effects alone and the fixed and random effects together were responsible for 13.1% and 59.0% of the between-individual variation in Bewick's Swan foraging effort (Table 3). A further five candidate models had associated AICc values within 6.0 of the best-supported models; however, each of these models accounted for a lower percentage of the total Akaike weights, had lower marginal and conditional R 2 values, and higher evidence ratios, than the best-supported model. Crucially, the variable describing the density of Mute and Whooper Swans was contained in all six of the top models, and in each case a positive effect size was detected. The best-supported model moreover performed substantially better than the null model (Table 3).

Table 4 Parameter estimates for the models of Bewick's Swan behaviour and energetics with the lowest AICc values
The variance and SD associated with the random effects, as well as the Variance Inflation Factor (VIF) associated with each parameter, are also indicated. Foraging and vigilance variables were logit-transformed response variable, whilst net energy gain was rescaled logit-transformed response variable. Model parameters: intercept (1), density of Bewick's Swan sharing feeding habitat with focal individual (B), density of Mute and Whooper Swans sharing feeding habitat with focal individual (C), age class of focal individual (A), and observation identity (U  . 2 The effect of combined Mute and Whooper Swan density within shared feeding habitat on the probability that adult and cygnet Bewick's Swans would be observed in an aggressive interaction, as predicted by our best-supported model (Table 3) The best-supported model for foraging indicated that the proportion of time spent by Bewick's Swans on foraging was positively related to the density of Bewick's Swans sharing the feeding habitat (Table 4; Fig. 3). However, a negative interaction between Mute and Whooper Swan density and age class was detected, such that the proportion of time spent by cygnets on foraging declined with Mute and Whooper Swan density, in contrast to the positive relationship found for adult Bewick's Swans (Table 4; Fig. 3). Fig. 3 The effects of Bewick's Swan density and combined Mute and Whooper Swan density on the proportion of time spent foraging by a adult and b cygnet Bewick's Swans, as predicted by our best-supported model (Table 3) Vigilance Based on our observations of focal swans, the mean (± 95% CI) proportion of time spent engaged in vigilance behaviour ranged from 0.177 (± 0.016) for adult Whooper Swans to 0.253 (± 0.032) for adult Mute Swans (Table 1). In comparison, adult and cygnet Bewick's Swans respectively spent 0.250 (± 0.022) and 0.233 (± 0.033) of their time on vigilance behaviour (Table 1).
Our model selection process identified eight noncollinear models that were potentially competitive, with ΔAICc values of < 6.0 ( Table 3). The model with the lowest AICc value was comprised of additive effects of age class, Bewick's Swan density, and the combined density of Mute and Whooper Swans, as well as an interaction between Bewick's Swan densities and combined Mute and Whooper Swan density; this model accounted for 37.0% of the total Akaike weights, whilst the fixed effects alone and the fixed and random effects together were responsible for 7.3% and 59.5% of the between-individual variation in Bewick's Swan vigilance (Table 3). The model with the next-lowest AICc value was identical to the first, except that it contained one additional parameter, an interaction between age class and Bewick's Swan density; however, despite the addition of this parameter the AICc value increased and so we considered the parameter to be uninformative and the model was not considered further (Table 3). The remaining six potentially competitive models each received substantially lower support in the data, with each model accounting for ≤ 7.6% of the total Akaike weights and their fixed effects accounting for ≤ 6.7% of the between-individual variation in vigilance. We therefore considered the lowest-AICc model to be our best-supported model of Bewick's Swan vigilance behaviour.
The best-supported model indicated that the proportion of time spent by Bewick's Swans on vigilance differed between adults and cygnets, with adults devoting more time to vigilance compared with cygnets (Table 4; Fig. 4). Individuals within low density Bewick's Swan flocks (< c. 200 individuals/km 2 ) spent more time on vigilance irrespective of age when amongst higher densities of Mute and Whooper Swans, whilst individuals within high density Bewick's Swan flocks showed the opposite pattern, as indicated by a negative interaction between these variables (Table 4; Fig. 4).

Net energy gain
The mean (± 95% CI) estimated rates of net energy gain values for Bewick's Swans were 2.0 (± 0.5) xBMR for adults and 1.4 (± 0.6) xBMR for cygnets (Table 1). The value for Bewick's Swan cygnets was the lowest estimated for any swan species and age class, whilst the maximum was 5.0 (± 0.7) xBMR for adult Mute Swans (Table 1).
Selection of the candidate models of Bewick's Swan net energy gain indicated that 13 non-collinear models had associated ΔAICc values of < 6.0 (Table 3). The implication of this is that no combination of explanatory variables received substantially more support in the data than the null model (ΔAICc = 0.72), and hence the evidence for any effects of these explanatory variables must be considered to be weak. In accordance with this interpretation, none of the candidate models achieved a marginal R 2 value of > 3.2%, suggesting that the explanatory variables could explain little of the between-individual variation in net energy gain (Table 3).

Discussion
In our study, we present the first detailed assessment of the potential effects of competition between the Bewick's Swan, and two larger competitors, the Mute Swan and Whooper Swan. Our investigation of individual behaviour and energetics indicated that, whilst Bewick's Swans exhibited some differences in behaviour in the presence of Mute and Whooper Swans, there were no detectable effects on estimated rates of net energy gain at our study site.
The overall level of aggression observed in our study population was low, with a mean (± SD) proportion of time spent on aggression by our three focal species of only 0.011 ± 0.035. In comparison, proportions of 0.033 ± 0.047 and 0.020 ± 0.036 were reported for swans and all waterbirds, respectively, in a recent meta-analysis of behavioural studies (Wood et al. 2017). Indeed, some previous studies of wintering waterbirds have reported up to 0.35 of time spent on aggressive interactions (e.g. Black et al. 1992). The low levels of aggression amongst our study population may reflect the open landscape and abundant food resources in the area surrounding the Ouse Washes, which allows the birds to space themselves out and thereby limit interference competition. Compared with natural plant communities, the coverage, biomass, and energetic density of agricultural crops is much more homogeneous (Robinson and Sutherland 2002). Hence, there is likely to be limited variation in feeding patch quality within a field, which will likely minimise both interference and depletion competition. In such a habitat there is little advantage to an individual in attempting to monopolise feeding areas, which reduces the likelihood of aggression over food resources.
Where aggressive interactions were observed, our data showed that these were typically directed towards individuals of the same species. Most aggression from Mute and Whooper Swans is therefore towards their conspecifics, rather than towards Bewick's Swans. Only for Mute Swan cygnets could a statistically significant bias towards conspecific aggression not be demonstrated, which likely reflects the small sample size obtained for this group. Our findings concur with earlier research on these three species in aquatic habitats, which found that proportionally between 0.59 and 0.80 of all observed aggression was intraspecific ). An earlier study of Mute Swans by Włodarczyk and Minias (2015) reported similarly that 0.80 of all aggressive behaviour involved conspecifics. Despite this predominance of intraspecific aggression, we did detect that Bewick's Swan cygnets had an increased likelihood of engaging in aggressive interactions when sharing feeding habitat with higher densities of interspecific competitors. In contrast, adult Bewick's Fig. 4 The effects of Bewick's Swan density and combined Mute and Whooper Swan density on the proportion of time spent on vigilance behaviour by a adult and b cygnet Bewick's Swans, as predicted by our best-supported model (Table 3) Swans showed no such density-dependent response. Our results may have been because cygnets are less experienced in social interactions and so are more likely to become involved in aggressive interactions when surrounded by higher numbers of competitors. It should be noted, however, that the proportion of between-individual variance accounted for by our best-supported model was relatively low (0.012; Table 3), which may reflect that the behaviours of individuals are influenced by a wide range of intrinsic and extrinsic variables (Rees et al. 2005;. Adult and cygnet Bewick's Swans showed opposing responses in their foraging behaviour to increasing interspecific competitor densities, with adults increasing their foraging effort whilst cygnets showed reduced foraging. The response of cygnets may be linked to our earlier finding regarding aggression, as individuals that engaged in aggressive interactions potentially had less time to spend on foraging. In contrast, the positive association between adult Bewick's Swan foraging effort and interspecific competitor densities may have been due to perceived depletion competition (Gyimesi et al. 2010), whereby individuals increased their foraging efforts to consume a greater share of a limited resource within a field. An alternative explanation for the positive association between adult Bewick's Swan foraging effort and interspecific competitor densities was that, because overwintering swans preferentially select the most profitable feeding habitats (Wood et al. 2019b), high density flocks indicated high quality feeding habitats, and so individuals would be expected to maximise foraging effort in such high quality habitats. Despite these effects, it is important to note that the mean proportions of time spent foraging by adult (0.414) and cygnet (0.420) Bewick's Swans were within the range of value reported previously for swans (e.g. Keane and O'Halloran 1992;O'Donoghue and O'Halloran 1994;O'Hare et al. 2007;Tatu et al. 2007;Nuijten et al. 2020c), which suggests that the Bewick's Swans in our study area were not suffering from limited access to food, either due to low availability or high interference competition.
There was also some evidence that sharing feeding habitat with interspecific competitors could affect the proportion of time that Bewick's Swans spent engaged in vigilance. However, the effect was complex, which may reflect the fact that for swans, both conspecific and heterospecific swans may provide vigilance against predators and humans, but also represent a potential threat due to the possibility of aggression and competition for shared food resources. When among low densities of conspecifics (< c. 200 individuals/km 2 ), individual Bewick's Swans spent more time on vigilance in the presence of higher densities of Mute and Whooper Swans, whilst individuals within higher density Bewick's Swan flocks showed the opposite pattern. This could indicate that, as winter Bewick's Swan numbers have decreased in recent decades, the remaining Bewick's Swans could have become more likely to find themselves within smaller Bewick's Swan flocks with higher vigilance costs. Individual Bewick's Swans among larger numbers of conspecifics may receive greater protection from interspecific aggression due a dilution effect (Delm 1990), i.e. whereby an individual Bewick's Swan has a lower chance of being targeted for interspecific aggression when it is within a larger group of Bewick's Swans. The contrasting conspecific and heterospecific effects on vigilance suggests that Bewick's Swans differentiate between other Bewick's Swans and their larger interspecific competitors within shared feeding habitat. Increased competition among higher densities of Bewick's Swans may also lead to increased foraging effort and hence reduced time available to spend on vigilance, given the inherent tradeoffs between different types of behaviour; if an individual spends more time on foraging, then the amount of time spent on another behaviour must be reduced to compensate (Metcalfe and Furness 1984). The extent to which foraging and vigilance behaviours are traded off against each other by individual swans is unclear, and future research could assess how swans trade-off expenditure of time on different behavioural activities.
Despite some evidence that Bewick's Swans altered their patterns of behaviour when sharing feeding habitat with interspecific competitors, their overall rates of net energy gain at the study site were not affected. Our results suggest that swans have at least some capacity to adjust the relative amounts of time devoted to key behaviours in order to maintain consistent levels of energy gain, as predicted by Wood et al. (2021). The mean estimates of net energy gain derived indicated that all species and age classes could achieve positive energy gains during winter. This finding concurs with earlier modelling work conducted for Bewick's and Whooper Swans in our study area (Wood et al. 2019b). The achievement of positive net energy gains during winter is particularly critical for migratory species such as Bewick's Swans, which must depart the winter grounds with sufficient energy reserves to reach a subsequent staging area or stopover site. For capital breeders such as swans, the energy and body reserves required to breed successfully is gained during the preceding winter (Drent et al. 2006) and at staging sites during migration (Nolet and Drent 1998). Our finding that Bewick's Swans wintering in eastern England have positive net energy gains suggests that competition with larger congeners during winter does not hinder their subsequent migration or breeding success.

Conclusions
In this study we found no evidence that Bewick's Swan net energy gain was affected when sharing agricultural feeding habitat with larger congeners during winter. These farmland habitats are the major feeding habitats of Bewick's Swans during winter (Beekman et al. 2019). We therefore consider it unlikely that competition between Bewick's Swans and either Mute or Whooper Swans in these agricultural habitats is a major cause of the observed decline in Bewick's Swan winter numbers. Further evidence against an impact of interspecific competition on Bewick's Swan numbers can be found in the spatial pattern of changes in numbers within the countries that comprise the winter grounds. Bewick's Swans continue to decline in the westerly countries of the range, such as the UK and the Netherlands, but are increasing in more easterly countries such as Germany (Augst et al. 2019;Beekman et al. 2019;Nuijten et al. 2020b). Whooper and Mute Swans, however, have both been increasing in numbers in both the western and eastern parts of the Bewick's Swans' winter grounds (Hall et al. 2015;Laubek et al. 2019;Wood and Włodarczyk 2020). If interspecific competition on the winter grounds was a major driver of Bewick's Swan decline through effects on behaviour and energetics, we would have expected to observe declines wherever the species overlapped, rather than only in the countries in the west of the range. The spatial pattern of declines and increases in Bewick's Swan numbers across the winter range are, however, consistent with short-stopping due to climate change, which has caused milder winters across Europe in recent years (Nuijten et al. 2020b). Further to the effects of climate, future research into fluctuations in Bewick's Swan numbers should investigate the role that known sources of mortality, such as illegal shooting (Newth et al. 2019), could have played in observed decline. In addition, other factors that could have led to deleterious environmental changes in the overwintering grounds, such as the expansion of energy infrastructure that poses collision risk to swans and other large birds, should be considered in future assessments.
It is important to note, however, that our findings do not preclude impacts of competition in other parts of the species' range, such as stopover sites and breeding areas. Some stopover sites have limited aquatic food resources that have been shown to be depleted by migrating Bewick's Swans during their stopover period (Nolet and Drent 1998). Increases of interspecific competitors at such sites could further limit the food supplies available to the sub-dominant Bewick's Swans. There is also the potential for increased competition to affect Bewick's Swan use of aquatic food resources such as macrophyte tubers, as studies have shown that greater summer grazing of above-ground macrophyte tissues by waterbirds such as Mute Swans can reduce the subsequent availability of tubers for Bewick's Swans (Hidding et al. 2009).