Kevin A. Wood, Phoebe Ham, Jake Scales, Eleanor Wyeth, Paul E. Rose. 2020: Aggressive behavioural interactions between swans (Cygnus spp.) and other waterbirds during winter: a webcam-based study. Avian Research, 11(1): 30. DOI: 10.1186/s40657-020-00216-7
Citation:
Kevin A. Wood, Phoebe Ham, Jake Scales, Eleanor Wyeth, Paul E. Rose. 2020: Aggressive behavioural interactions between swans (Cygnus spp.) and other waterbirds during winter: a webcam-based study. Avian Research, 11(1): 30. DOI: 10.1186/s40657-020-00216-7
Kevin A. Wood, Phoebe Ham, Jake Scales, Eleanor Wyeth, Paul E. Rose. 2020: Aggressive behavioural interactions between swans (Cygnus spp.) and other waterbirds during winter: a webcam-based study. Avian Research, 11(1): 30. DOI: 10.1186/s40657-020-00216-7
Citation:
Kevin A. Wood, Phoebe Ham, Jake Scales, Eleanor Wyeth, Paul E. Rose. 2020: Aggressive behavioural interactions between swans (Cygnus spp.) and other waterbirds during winter: a webcam-based study. Avian Research, 11(1): 30. DOI: 10.1186/s40657-020-00216-7
Our understanding of any impacts of swans on other waterbirds (including other swans), and potential effects on waterbird community structure, remain limited by a paucity of fundamental behavioural and ecological data, including which species swans interact aggressively with and how frequently such interactions occur.
Methods
Behavioural observations of aggression by swans and other waterbirds in winters 2018/2019 and 2019/2020, were carried out via live-streaming webcams at two wintering sites in the UK. All occurrence sampling was used to identify all aggressive interactions between conspecific or heterospecifics individuals, whilst focal observations were used to record the total time spent by swans on aggressive interactions with other swans. Binomial tests were then used to assess whether the proportion of intraspecific aggressive interactions of each species differed from 0.5 (which would indicate equal numbers of intraspecific and interspecific interactions). Zero-inflated generalized linear mixed effects models (ZIGLMMs) were used to assess between-individual variation in the total time spent by swans on aggressive interactions with other swans.
Results
All three swan species were most frequently aggressive towards, and received most aggression from, their conspecifics. Our 10-min focal observations showed that Whooper (Cygnus cygnus) and Bewick's Swans (C. columbianus bewickii) spent 13.8± 4.7 s (means± 95% CI) and 1.4± 0.3 s, respectively, on aggression with other swans. These durations were equivalent to 2.3% and 0.2% of the Whooper and Bewick's Swan time-activity budgets, respectively. Model selection indicated that the time spent in aggressive interactions with other swans was best-explained by the number of other swans present for Whooper Swans, and an interactive effect of time of day and winter of observation for Bewick's Swans. However, the relationship between swan numbers and Whooper Swan aggression times was not strong (R2= 19.3%).
Conclusions
Whilst swans do exhibit some aggression towards smaller waterbirds, the majority of aggression by swans is directed towards other swans. Aggression focused on conspecifics likely reflects greater overlap in resource use, and hence higher potential for competition, between individuals of the same species. Our study provides an example of how questions relating to avian behaviour can be addressed using methods of remote data collection such as live-streaming webcams.
Reintroduction aims to establish a viable population of the focal species within its indigenous range (IUCN/SSC 1998, 2013). However, any reintroduction projects may fail to establish populations owing to either external or intrinsic factors that influence wildlife demography and long-term viability in different establishment phases (Armstrong and Seddon 2008; Robert et al. 2015). In particular, small and isolated populations are more at risk of extinction by the additional threats of chance events, such as demographic, environmental, genetic stochasticities, and catastrophic events. In these situations, population viability analysis (PVA), as an effective tool to assess the fates of populations under various scenarios of the deterministic forces and chance events, has been widely used in conservation biology for many endangered mammals and birds (Song 1996; Brook et al. 2000; Perkins et al. 2008; Hernandez-Camacho and Trites 2018; Kimanzi 2018; Brook et al. 2019; Blazquez et al. 2020; Fryxell et al. 2020; He et al. 2020). Specifically, PVA can provide some information about population dynamics, such as population size, extinction probability and genetic diversity change. For example, the viability of a reintroduced Bearded Vulture (Gypaetus barbatus) population was mainly affected by mortality and artificial supplementary release (Bustamante 1996). The current population of the Oriental White Stork (Ciconia boyciana) has a 12.5% extinction probability in the next 100 years. Additionally, nestling mortality and the proportion of adult females involved in reproduction greatly impact population growth (Sung et al. 2012). The PVA results of the Western Lowland Gorillas (Gorilla gorilla gorilla) suggested that the reintroduced gorilla populations have a reasonable chance of persistence over 200 years, but the reinforcement of the populations could significantly improve probabilities of population persistence and retention of genetic diversity (King et al. 2014). Simulation results of PVA of the Sharp-tailed Grouse (Tympanuchus phasianellus) showed that, in a fixed scenario, improvement of habitat quality could facilitate the establishment of a minimum viable population by increasing fecundity and wintering survival (Milligan et al. 2018).
In the case of the Crested Ibis (Nipponia nippon) whose population once was thought to be extirpated (Archibald et al. 1980), its wild population has recovered to approximately 4400 (Shaanxi Hanzhong Crested Ibis National Nature Reserve, unpublished data) and several reintroduction programs have been successively carried out in China (Yu et al. 2015), Japan (Nagata and Yamagishi 2016) and Korea (Yoon and Choi 2018). However, up to date, population viability analysis on the endangered species has been limited to the wild population, showing a probability of extinction of 19.7% by 2097; long-term population growth is restricted by the carrying capacity and inbreeding, as it is more sensitive to catastrophes and habitat conditions (Li and Li 1998). A reconstructed PVA model predicted that the wild population would not go extinct in the next 50 years and revealed that population growth is significantly affected by the carrying capacity (Li 2013). Differences in the parameters used in the two PVA models led to variable prediction results.
Since the first reintroduction of the Crested Ibis in Ningshan, Shaanxi, China (Yu et al. 2009), new pair formations have continued to occur, and the number of fledglings and the distribution range has continued to increase. More recently, the population has survived the establishment stage, with high post-release dispersal and low survival, and it is now considered to be at the persistence stage (Huo et al. 2014; Yu et al. 2015; Dong et al. 2018). A recent study on the survival showed that the number and the range of the reintroduced population gradually increased without further releases (Li et al. 2018). Moreover, stochastic estimations of the population growth rate and population size simulated from an age-classified Leslie matrix indicated that the reintroduced population is likely to be at regulation phase over the next 25 years (Wang et al. 2017). However, gaps in knowledge of the viability of this population and its drivers in the release environment impede the identification of the most effective population-level priorities for aiding in species recovery. Here, we used a VORTEX model to simulate how the future dynamics of the population in the next 50 years to determine the key factors affecting its survival through sensitivity analysis. The process we used can serve as an example for developing a framework of conservation priorities on the population level.
Methods
Study area
The fieldwork was conducted in Ningshan County (33°07′09″‒33°50′32″N, 108°02′15″‒108°56′44″E) of Shaanxi Province, which is located on the southern slope of the Qinling Mountains, China. The county is an important transportation hub connecting the Weihe Valley with the Sichuan-Chongqing region. Natural conditions of the study area have been reported in detail elsewhere (Yu et al. 2015). Most individuals of the reintroduced population were found in the middle reaches of the Changan River, one of the tributaries of the Hanjiang River (Li et al. 2018).
Study population
A total of 56 captive-bred individuals (27 females, 29 males) from two captive Crested Ibis populations were released to the wild from 2007 to 2011 (Yu et al. 2015). In addition, five individuals from the wild population in Yang County immigrated to Ningshan County through natural dispersal. By the end of the breeding season in 2020, a total of 220 offspring have fledged in Ningshan County. All released individuals and more than 95% of the wild-born fledglings were marked with colored alphanumeric plastic bands (provided by the National Bird Banding Center of China) for field identification of their age and sex (Yu et al. 2015). Prior to release, eleven individuals were equipped with lightweight radio transmitters (model RI-2D; frequency 216.368‒216.691 MHz; battery life ~ 18 months; Holohil Systems, Carp, Ontario, Canada) and two ibises with GPS trackers (Hunan Global Messenger Technology Co., Ltd., HQBG2512S) to monitor their subsequent fates.
Population modeling
A VORTEX model (Version 10.3.5.0) for PVA, specific to vertebrates with a long lifespan and low fecundity (Miller and Lacy 2005), was used to predict dynamics of small populations as well as the future probability of population extinction. The PVA model has the potential to simulate the fates of each individual and gene flows among generations based on various random variables such as demographic, environmental, catastrophic and stochastic genetic events (Lacy 1993).
The demographic parameters in the VORTEX model were estimated based on the monitoring data of reintroduced Crested Ibis population from 2007 to 2018 (Table 1), and the model was used to simulate population dynamics over the next 50 years after 2018 to estimate the influencing parameters by the standard sensitivity index (Sx) (Lacy 1993; Mccarthy 1995; Reed et al. 1998; Akcakaya and Sjögren-Gulve 2000; Morris and Doak 2002).
Table
1.
Parameters of the VORTEX model used to simulate Crested Ibis population growth
Variable
Explanation or notes
Default values
# of iterations
# of times the simulation is run
500
# of years
# of years the population is modeled
50 (2019‒2068)
Extinction definition
How extinction is defined in the model
Only one sex remains
# of populations
# of populations to be modeled
1
Inbreeding depression
The level and effect of inbreeding depression
Yes
Lethal equivalents
The severity of inbreeding depression
0
EV concordance of reproduction and survival
The concordance between reproduction and survival in different populations
The current populations of the Crested Ibis including wild, captive and reintroduced populations were derived from seven individuals approximately 40 years ago, so each of them is an inbred population (Li et al. 2014; Feng et al. 2019; Fu et al. 2019). However, the Crested Ibis has survived the negative effect of inbreeding depression, and its distribution range and population size progressively expanded. It may be that the protracted inbreeding history has facilitated the purging of deleterious recessive alleles (Zhu et al. 2018). And gene mutation and natural selection can also eliminate the harmful recessive genes (Byers and Waller 1999; Keller and Waller 2002). The previous study on the wild Crested Ibis population showed that inbreeding has no significant impact on population growth (Li 2013). Therefore we assumed that the lethal equivalent value to zero for each diploid in the target population. The study population had a balanced sex ratio (1:1) in the initial population completely composed of captive-bred individuals, as well as the current population including released ibises and wild-born individuals. Population extinction was defined as all individuals of one sex are disappeared.
We assumed that Crested Ibis would be density-dependent, given that it is a territorial bird with limited nesting sites and food sources. Thus, the formula used in the VORTEX model was:
in which P(N) is the percentage of adult female breeders when the population size is N; P(K) is the percentage of adult female breeders at carrying capacity (K); A is Allee parameter, which accounts for the decrease in the proportion of female breeders at low densities because of the increased difficulty of finding a mate at low population densities; And P(0) is the percentage of adult female breeders at low densities when there is no Allee effect. The previous study revealed that there has a weak Allee effect (Chen 2018) and a positive annual growth rate in the target population (Wang et al. 2017), therefore we assumed the parameter A to be zero. From 2007 to 2012, the focal ibis population, with a small size, was at the establishment stage (Additional file 1: Table S1), therefore we set P(0) to 35.70. For the territorial Crested Ibis with limited nesting sites and food sources, there will be intense intraspecific competition when the population size is close to the carrying capacity which may affect the number of female breeders in the population due to the limited availability resources (Henle et al. 2004). Thus, we assumed an ideal value of P(K) is 50%. The parameter B was set to 2, which means that P(N) is a quadratic function of N. This function was commonly used to simulate the growth of a population with density dependence (Ralls et al. 1988).
Two types of catastrophes faced by our study population are floods and spring cold weather. Floods will inundate paddy fields, which will seriously affect Crested Ibis foraging, and the spring cold weather will threaten the chicks and reduce the food abundance. Based on meteorological data and reports of Ningshan County, the frequency of flood and spring cold weather are approximately 14% and 15%, respectively. The effect of catastrophes on the reproduction and survival of Crested Ibis in Ningshan is unclear, so we assumed that 90% of the individuals survived a catastrophe and of those that survived, reproduction was reduced by 10% in the year of the catastrophes. Due to the wide range of Crested Ibis activities and understrength monitoring personnel, there is no detailed mortality data of Crested Ibis in Ningshan County at present. We used the survival of adults, subadults, and chicks to calculate the maximum mortality, and survival was 0.941, 0.849, and 0.549, respectively (Wang 2015). Therefore, the maximum mortality of adults, subadults, and chicks was roughly estimated to be 5.9%, 15.1% and 45.1%, respectively. To measure the impact of catastrophes on the population, we increased or reduced the frequency of catastrophes by 10% and 50% from the original baseline value. In the same way, the mortality rate of each age class was increased or decreased by 10%, to measure the impact of the mortality on the population.
In a previous study, the total area of suitable habitats in Ningshan County was estimated to be 167.36 km2 using MaxEnt model (Wei 2020). The wild Crested Ibis population growth in the core nesting range may have been experiencing the suppression from the environment capacity at present, with the nesting density in the core range of 0.73 nests/km2 (Wang et al. 2020), namely average home-range size for each breeding pair was approximately 1.37 km2 in the wild population. The home-range size was restricted by the food, space, and other habitat environmental factors. Yangxian County can provide sufficient food resources for the Crested Ibis through dense water systems, paddy fields, and diverse vegetation cover types. However, the breeding and foraging grounds used by ibises in our focal population are located in intermontane basins that are often developed for use as paddy fields and interspersed with hills (Yu et al. 2015). The food abundance and availability of the Crested Ibis in Ningshan County is lower than that of Yang County. According to the average home-range size of the wild population and the difference of food resources between the two habitats, the home-range size of Crested Ibis in the study area is bigger than that in Yang County. Therefore, we estimated the average territory area of a breeding pair is estimated to be 2 km2. The average number of fledgelings per nest is roughly estimated to be 3. Thus, we estimated the carrying capacity (K) of Ningshan County to be 420 ± 80 individuals, including 84 breeding pairs, with approximately three young birds fledged per nest.
Results
In our simulation, the reintroduced ibis population is unlikely to go extinct but exhibits a stead growing trend within the next 50 years. Under current natural conditions, the intrinsic rate of increase (r) of the population was estimated to be 0.084 with an SD (r) of 0.071. The number of individuals will increase to 367, and the genetic diversity of the population will drop to 0.97 by the end of the next 50 years.
The sensitivity index (Sx) of the population clearly showed that carrying capacity, sex ratio, frequency of catastrophes, and mortality were the four factors accounting for the population parameters (Table 2).
Table
2.
Results of different model parameters during the next 50 years and sensitivity analysis
The carrying capacity and sex ratio have great effects on population size and extinction probability. If the carrying capacity was set to double or halve the current carrying capacity, the simulation results showed that the carrying capacity is the main factor restricting the population size. The smaller the carrying capacity hold, the shorter the duration required for the population to reach the K value (Fig. 1a–c). The carrying capacity also affects the genetic diversity of the population. When the carrying capacity is reduced to half of the baseline value, the genetic diversity of the population is significantly different from that of the initial population (one-way ANOVA: F2, 147 = 28.614, P < 0.05) (Fig. 2). We set six conditions of sex ratio imbalance to simulate population dynamics, and the results showed that if the sex ratio becomes extreme imbalanced, the population size will increase slowly over the next 50 years or even decline to the brink of extinction (Table 2, Fig. 1d–i). When the sex ratio has shifted to females or males to the same degree (such as 25% male vs. 75% male), the population shows an opposite growth trend (Table 2, Fig. 1f, g).
Figure
1.
Impact of different model parameters on the reintroduced population size. The data graph shows a present conditions, bK = 840, (c) K = 210, d females: males = 1:2, e females: males = 2:1, f females: males = 1:3, g females: males = 3:1, h females: males = 1:4, and i females: males = 4:1
Figure
2.
Impact of carrying capacity on genetic diversity during the next 50 years. Three scenarios are imagined: the carrying capacity is 420 (blue), 840 (orange), and 210 (green), respectively
The results showed that although the frequency of catastrophes will not have a significant impact on population size, it will prevent the population from reaching the carrying capacity. That is the higher the frequency of catastrophes, the longer the population takes to reach the expected K value (Fig. 3). Sensitivity analysis of mortality rates for each age class showed that changes in infant mortality had a greater impact on future population growth than changes in other age groups (Table 2).
Figure
3.
Impact of different catastrophe frequencies on the population size. The data graph shows a baseline value of frequency of catastrophes, b frequency of catastrophes increased by 10%, c frequency of catastrophes decreased by 10%, d frequency of catastrophes increased by 50%, e frequency of catastrophes decreased by 50%
Since the release of 56 captive-reared Crested Ibises in 2007, long-term monitoring has been conducted to collect various types of field data. However, it is difficult to exactly quantify the mortality and a total number of individuals due to the large range of individual activities. In addition, we estimated the current population size to be approximately half of the carrying capacity. As a result, the parameters in the model, such as K and P(K), need to be roughly estimated. In conclusion, a certain amount of bias in simulation results may be produced because of insufficient field data.
According to the simulation results of the PVA model, the population size will constantly increase in the next 20 years, which means that the population will have reached the maximum carrying capacity. The population is resistant to slight fluctuations of catastrophe frequency, mortality, and sex ratio, indicating that it can also become a self-sustaining population without further releases. Despite being consistent with the results of population sustainability analysis, the population growth rate may still be affected by the following factors.
Impact of inbreeding depression
Inbreeding depression, commonly occurred in small populations, may lead to the decrease of population genetic diversity by increasing the frequency of homozygosity of harmful alleles. As a result, the reduced adaptability of offspring is likely to increase the risk of population extinction (Ellstrand and Elam 1993; Fischer et al. 2003).
However, under current circumstances, some studies indicated that the Crested Ibis has evolved a series of mechanisms to avoid inbreeding, such as dispersal and kin recognition (Macdonald and Johnson 2001; Liu and Zhang 2008). Besides, genetic diversity in small populations is also affected by genetic drift that can reduce the frequency of harmful genes (Bortoluzzi et al. 2020). Based on the current study and population status, we set the model parameters—lethal gene equivalent to zero in the present study. Therefore, the predicted population size might slightly overestimate the actual status of the population in the future. Combination of genetic diversity study with PVA should be made to estimate the exact lethal gene equivalent by investigating the relationship between observed heterozygosity and inbreeding coefficient, allowing more accurate prediction of the future dynamics for the ibis populations.
Impact of carrying capacity
In the sensitivity analysis, the higher sensitivity index occurs when the carrying capacity changes (Table 2), indicating that the carrying capacity has a great impact on the population size. From the results, we inferred there is a great potential for an expansion of the population, when the carrying capacity of the Crested Ibis is great. So carrying capacity is one of the primary factors restricting population size. The genetic diversity of the population can also be influenced by the carrying capacity because the mutation rate of small populations is much lower than its loss rate (Lande and Barrowclough 1987). When the mutation rate and loss rate of one population are at the same level under a certain scale of population size, the population can be adaptive to the change of environment and survive for a long time (Fan and Jiang 2007). Therefore, increasing the carrying capacity can increase the genetic heterozygosity of a population, which is beneficial to the long-term population growth.
Because of different emphases of the spatial carrying capacity and nutrient carrying capacity, the estimation results may vary to a great extent (Mentis and Duke 1976; Wallmo et al. 1977; Hobbs et al. 1982). For example, when estimating carrying capacity from the perspective of territory area, habitat fragmentation should be taken into account (Haila 2002; Wu and Li 2003). Carrying capacity may become smaller with the increase of habitat fragmentation, thus much attention should be paid to the edge effect caused by habitat fragmentation (Tallmon and Mills 2004). As for the present study, the carrying capacity of our focal population still remains unclear. Hence, more accurate estimation on the carrying capacity of the ibis population can be achieved by constructing the mathematical models under varied landscape scale in combined with individual nutrient requirements and habitat nutrient supply.
Effect of sex ratio
At present, the age structure of the Crested Ibis population in Ningshan County seems to be reasonable, with the sex ratio roughly at 1:1 (Li et al. 2018). The sensitivity analysis showed that the population size is highly sensitive to the sex ratio, which means that a certain degree of sex ratio deviation will exert a great impact on the subsequent fate of the population.
There is a varied impact of the sex ratio on population size of the current and wild population (Li 2013). A sex ratio imbalance of the same proportion (such as 25% male) will lead to a decline in the natality of Ningshan ibis population and then population extinction. However, the size of wild population in Yang County will only increase slowly, which may be due to the fact that the current population size of Ningshan County is much smaller than that of Yang County. Our results also showed that when the population sex ratio skews to males or females to the same extent (such as 25% male vs. 75% male), the population exhibits the opposite trend of increase or extinction, respectively. Under the influence of model parameter "mating monopoly", when there are fewer males in the population, the mating probability of the adult males is higher than that of the original model parameters, which may lead to an underestimation of simulation.
Conservation and management recommendations
Despite the initial successful reintroductions of the Crested Ibis in China, no predictive framework is available to quantify which general plan is vital to the recovery of the Crested Ibis within their former range. The present population is considered to be at the persistence stage, but threats at the population and metapopulation level remain. Some conservation management strategies that need to be addressed are listed as follows.
Nature reserves should be established to improve habitat quality in order to increase the carrying capacity. Additionally, local farmers should be encouraged to cultivate larger areas of paddy fields to ensure adequate food supply for the Crested Ibis. At the metapopulation level, habitat restoration and expansion are pressing issues for facilitating dispersal to promote individual exchanges between the wild and reintroduced populations. Priority should be given to determining the effects of predators on population dynamics and unknown reasons for individual mortality. Laws and regulations associated with management should be publicized, and punishment for illegal activities should be strictly enforced. Public education should be implemented to raise awareness of wildlife conservation.
Additional file 1: Table S1. The number of adult females and males that bred from 2008 to 2018. Table S2. Clutch size of Crested Ibis from 2008 to 2018.
Acknowledgements
We thank the staff of the Reintroduction Base of Crested Ibis in Ningshan County for providing logistical help and research aid. And special thanks to Robert Lacy and the IUCN/SSC for providing us the software for the VORTEX model.
Authors' contributions
YZ and ML conceived and designed the study. YZ, FW, ML and XL conducted the field work. YZ and FW carried out the analyses. YZ prepared the draft of the manuscript ZC, XYu and XYe revised it. All authors read and approved the final manuscript.
Availability of data and materials
The datasets used in the present study are available from the corresponding author on reasonable request.
Declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Altmann J. Observational study of behavior: sampling methods. Behaviour. 1974;49:227–66.
Amat JA. Food usurpation by waterfowl and waders. Wildfowl. 1990;41:107–16.
Anderson MJ, Urbine JL, Wilson C, Callabro L. Employment of web-based images and a live web cam in the examination of lateral neck-resting preferences in the American flamingo (Phoenicopterus ruber). J Caribb Ornithol. 2011;24:41–7.
Arnold TW. Uninformative parameters and model selection using Akaike's Information Criterion. J Wildl Manage. 2010;74:1175–8.
Bailey RO, Batt BDJ. Hierarchy of waterfowl feeding with Whistling Swans. Auk. 1974;91:488–93.
Barton K. MuMIn: Multi-Model Inference. Version 1.43.15. 2019. .
Beekman J, Koffijberg K, Wahl J, Kowallik C, Hall C, Devos K, et al. Long-term population trends and shifts in distribution for Bewick's Swans Cygnus columbianus bewickii wintering in northwest Europe. Wildfowl. 2019; Special Issue 5: 73–102.
Beven G. Coot feeding on weed disturbed by Mute Swans. Brit Birds. 1980;73:219–20.
Black JM, Rees EC. The structure and behaviour of the Whooper Swan population wintering at Caerlaverock, Dumfries and Galloway, Scotland: an introductory study. Wildfowl. 1984;35:21–36.
Bowler JM. Feeding strategies of Bewick's Swans (Cygnus columbianus bewickii) in winter. PhD Thesis. Bristol: University of Bristol; 1996.
Brazil MA. A case of unusual aggression by a Whooper Swan. Tori. 1983;32:155.
Brides K, Wood KA, Hearn RD, Fijen TPM. Changes in the sex ratio of the Common Pochard Aythya ferina in Europe and North Africa. Wildfowl. 2017;67:100–12.
Brooks ME, Kristensen K, van Benthem KJ, Magnusson A, Berg CW, Nielsen A, et al. glmmTMB balances speed and flexibility among packages for zeroinflated generalized linear mixed modeling. R J. 2017;9:378–400.
Burgess RM, Stickney AA. Interspecific aggression by Tundra Swans towards Snow Geese on the Sagavanirktok River Delta, Alaska. Auk. 1994;111:204–7.
Burnham KP, Anderson DR, Huyvaert KP. AIC model selection and multimodel inference in behavioral ecology: some background, observations, and comparisons. Behav Ecol Sociobiol. 2011;65:23–35.
Clopper CJ, Pearson ES. The use of confidence or fiducial limits illustrated in the case of the binomial. Biometrika. 1934;26:404–13.
Conover MR, Kania GS. Impact of interspecific aggression and herbivory by mute swans on native waterfowl and aquatic vegetation in New England. Auk. 1994;111:744–8.
Crawley MJ. The R Book. 2nd ed Chichester: Wiley; 2013.
Davis JB, Guillemain M, Kaminski RM, Arzel C, Eadie JM, Rees EC. Habitat and resource use by waterfowl in the northern hemisphere in autumn and winter. Wildfowl. 2014;4:17–69.
Delacour J. Waterfowl in large mixed collections. Int Zoo Yearb. 1973;13:15–9.
Dormann CF, Elith J, Bacher S, Buchmann C, Carl G, Carré G, et al. Collinearity: a review of methods to deal with it and a simulation study evaluating their performance. Ecography. 2013;36:27–46.
Eichorst B. Internet webcams provide opportunities for college student research on animal behavior and ecology: an example with birds. Am Biol Teacher. 2018;80:680–5.
Ely CR, Budeau DA, Swain UG. Aggressive encounters between Tundra Swans and Greater White-fronted Geese during brood rearing. Condor. 1987;89:420–2.
Frost TM, Calbrade NA, Birtles GA, Mellan HJ, Hall C, Robinson AE, et al. Waterbirds in the UK 2018/2019: The Wetland Bird Survey. Thetford: BTO/RSPB/ JNCC; 2020.
Gayet G, Guillemain M, Mesléard F, Fritz H, Vaux V, Broyer J. Are Mute Swans (Cygnus olor) really limiting fishpond use by waterbirds in the Dombes, Eastern France? J Ornithol. 2011;152:45–53.
Gayet G, Calenge C, Broyer J, Mesléard F, Vaux V, Fritz H, et al. Analysis of spatial point pattern shows no desertion of breeding Mute Swan areas by the other waterbirds within fishpond. Acta Ornithol. 2016;5:151–63.
Gillham ME. Feeding habits and seasonal movements of mute swans on two south Devon estuaries. Bird Study. 1956;3:205–12.
Gurtovaya EN. Aggressive interactions between Bewick's Swans and other Anseriformes in the breeding period. Casarca. 2000;6:167–76.
Gyimesi A, Stillman RA, Nolet BA. Cryptic interference competition in swans foraging on cryptic prey. Anim Behav. 2010;80:791–7.
Gyimesi A, van Lith B, Nolet BA. Commensal foraging with Bewick's Swans Cygnus bewickii doubles instantaneous intake rate of Common Pochard Aythya ferina. Ardea. 2012;100:55–62.
Holm S. A simple sequentially rejective multiple test procedure. Scand J Stat. 1979;6:65–70.
Johnsgard PA. Handbook of waterfowl behavior. New York: Cornell University Press; 1965.
Källander H. Commensal association of waterfowl with feeding swans. Waterbirds. 2005;28:326–31.
King JA. The ecology of aggressive behavior. Annu Rev Ecol Syst. 1973;4:117–38.
Lenth R. Emmeans: Estimated Marginal Means, aka Least-Squares Means. Version 1.4.5. 2020. .
Lind H. The rotation display of the Mute Swan Cygnus olor: synchronised neighbour responses as instrument in the territorial defence strategy. Ornis Scand. 1984;15:98–104.
Lüdecke D. Sjstats: Statistical Functions for Regression Models. Version 0.17.9. 2020. .
Lumsden HG. Trumpeter Swans and Mute Swans compete for space in Ontario. Ontario Birds. 2016;34:14–23.
Mac Nally R, Duncan RP, Thomson JR, Yen JD. Model selection using information criteria, but is the "best" model any good? J Appl Ecol. 2018;55:1441–4.
Marchowski D, Neubauer G. Kleptoparasitic strategies of Mallards towards conspecifics and Eurasian Coots. Ardea. 2019;107:110–4.
Metcalfe NB, Furness RW. Aggression in shorebirds in relation to flock density and composition. Ibis. 1987;129:553–63.
Nakagawa S, Johnson PC, Schielzeth H. The coefficient of determination R2 and intra-class correlation coefficient from generalized linear mixedeffects models revisited and expanded. J R Soc Int. 2017;14:20170213.
Newth JL, McDonald RA, Wood KA, Rees EC, Semenov I, Chistyakov A, et al. Predicting intention to hunt protected wildlife: a case study of Bewick's swans in the European Russian Arctic. Oryx. (in press).
O'Hare MT, Stillman RA, McDonnell JO, Wood LR. Effects of mute swan grazing on a keystone macrophyte. Freshwater Biol. 2007;52:2463–75.
Peiman K, Robinson B. Ecology and evolution of resource-related heterospecific aggression. Quart Rev Biol. 2010;85:133–58.
Pelligrini AD. The roles of aggressive and affiliative behaviors in resource control: a behavioral ecological perspective. Dev Rev. 2008;28:461–87.
Peluso AI, Royer EA, Wall MJ, Anderson MJ. The relationship between environmental factors and flamingo aggression examined with internet resources. Avian Biol Res. 2013;6:215–20.
Pöysä H. Resource utilization pattern and guild structure in a waterfowl community. Oikos. 1983;40:295–307.
R Core Team. R: A language and environment for statistical computing. [3.6.3]. Vienna, Austria: R Foundation for Statistical Computing; 2020.
Rees EC. Bewick's Swan. London: T & AD Poyser; 2006.
Rees EC, Cao L, Clausen P, Coleman JT, Cornely J, Einarsson O, et al. Conservation status of the world's swan populations, Cygnus sp. and Coscoroba sp.: a review of current trends and gaps in knowledge. Wildfowl. 2019; Special Issue 5: 35–72.
Richards SA. Dealing with overdispersed count data in applied ecology. J Appl Ecol. 2008;45:218–27.
RSPB. Birds A- Z: Bird Guides. Royal Society for the Protection of Birds. . Accessed October 2018.
Schulwitz SE, Spurling DP, Davis TS, McClure CJW. Webcams as an untapped opportunity to conduct citizen science: siz years of the American Kestrel Partnership's KestrelCam. Glob Ecol Conserv. 2018;15:e00434.
Scott DK. Social behaviour of wintering Cygnus columbianus bewickii. In: Matthews GVT, Smart M, editors. Proceedings of the Second International IWRB Swan Symposium, Sapporo, Japan, 1980. Slimbridge: International Waterfowl Research Bureau; 1981. p. 211–25.
Shimada T. Ducks foraging on swan faeces. Wildfowl. 2012;62:224–7.
Sladen WJL. Swans should not be hunted. In: Sears J, Bacon PJ, editors. Proceedings of the Third International IWRB Swan Symposium, Oxford, 1989. Slimbridge: International Waterfowl Research Bureau; 1991. p. 368–75.
Stone WB, Marsters AD. Aggression among captive Mute Swans. New York Fish Game J. 1970;17:51–3.
Tatu KS, Anderson JT, Hindman LJ, Seidel G. Diurnal foraging activities of mute swans in Chesapeake Bay, Maryland. Waterbirds. 2007;30:121–9.
Therres GD, Brinkler DF. Mute Swan interactions with other birds in Chesapeake Bay. In: Perry MC, editor. Mute Swans and Their Chesapeake Bay Habitats: Proceedings of a Symposium. Virginia: US Geological Survey; 2004. p. 43–6.
Tingay A. Aggression in the Black Swan. Emu. 1974;74:35–8.
Vogrin M. A Coot Fulica atra eating waterfowl droppings. Butll GCA. 1997;14:63–4.
Włodarczyk R, Minias P. Division of parental duties confirms a need for biparental care in a precocial bird, the mute swan Cygnus olor. Anim Biol. 2015;65:163–76.
Wood KA, Stillman RA, Goss-Custard JD. The effect of kleptoparasite and host numbers on the risk of food-stealing in an avian assemblage. J Avian Biol. 2015;46:589–96.
Wood KA, Ponting J, D'Costa N, Newth JL, Rose PE, Glazov P, et al. Understanding intrinsic and extrinsic drivers of aggressive behaviour in waterbird assemblages: a meta-analysis. Anim Behav. 2017;126:209–16.
Wood KA, Cao L, Clausen P, Ely CR, Luigujõe L, Rees EC, et al. Current trends and future directions in swan research: insights from the 6th International Swan Symposium. Wildfowl. 2019a; Special Issue 5:1–34.
Wood KA, Hilton GM, Newth JL, Rees EC. Seasonal variation in energy gain explains patterns of resource use by avian herbivores in an agricultural landscape: insights from a mechanistic model. Ecol Model. 2019b;409:108762.
Zuur AF, Ieno EN, Elphick CS. A protocol for data exploration to avoid common statistical problems. Methods Ecol Evol. 2010;1:3–14.
Jiabin Liu, Jiaojiao Yu, Wenlei Bi, et al. Genetic status assessment and future development forecast for an isolated giant panda population. Global Ecology and Conservation, 2025, 57: e03423.
DOI:10.1016/j.gecco.2025.e03423
2.
Soodong Lee, Chunghyeon Oh, Bonggyo Cho, et al. Understanding the home range characteristics of the first naturally bred pair of crested ibis(Nipponia nippon) released into the natural habitat. BMC Zoology, 2024, 9(1)
DOI:10.1186/s40850-024-00220-0
3.
Yilamujiang Tuohetahong, Ruyue Lu, Ruiyan Guo, et al. Climate and land use/land cover changes increasing habitat overlap among endangered crested ibis and sympatric egret/heron species. Scientific Reports, 2024, 14(1)
DOI:10.1038/s41598-024-71782-9
4.
Jes Hooper. Thinking with Civets: The Role of Zoos in the Decolonisation of Animal Tourism. Animals, 2023, 13(11): 1739.
DOI:10.3390/ani13111739
5.
Youwen Wang, Maria Vasilyeva, Sergei Stepanov, et al. Numerical Investigation and Factor Analysis of Two-Species Spatial-Temporal Competition System after Catastrophic Events. WSEAS TRANSACTIONS ON SYSTEMS, 2023, 22: 423.
DOI:10.37394/23202.2023.22.45
6.
Yulong Zhou, Xian Jiang, Zhanlong Chen. Spatial- and Temporal-Trajectory Analysis of the Crested Ibis (Nipponia nippon) by Fusing Multiple Sources of Data. Animals, 2023, 13(2): 237.
DOI:10.3390/ani13020237
7.
Yuji Okahisa, Hisashi Nagata. Evaluation of ongoing Crested Ibis (Nipponia nippon) reintroduction using an integrated population model and Bayesian population viability analysis. Ibis, 2022, 164(4): 1104.
DOI:10.1111/ibi.13071
8.
Min Li, Rong Dong, Yilamujiang Tuohetahong, et al. Impact of Allee effects on the establishment of reintroduction populations of endangered species: The case of the Crested Ibis. Global Ecology and Conservation, 2022, 35: e02103.
DOI:10.1016/j.gecco.2022.e02103
9.
Xian Jiang, Tingdong Yang, Dongping Liu, et al. An Automatic Identification Method of Crested Ibis (Nipponia nippon) Habitat Based on Spatiotemporal Density Detection. Animals, 2022, 12(17): 2220.
DOI:10.3390/ani12172220
10.
Min Li, Yilamujiang Tuohetahong, Feng Lin, et al. Factors affecting post-release survival and dispersal of reintroduced Crested Ibis (Nipponia nippon) in Tongchuan City, China. Avian Research, 2022, 13: 100054.
DOI:10.1016/j.avrs.2022.100054
11.
Xinhai Li, Jiayu Fu, Tianqing Zhai, et al. Understanding Recovery Is as Important as Understanding Decline: The Case of the Crested Ibis in China. Land, 2022, 11(10): 1817.
DOI:10.3390/land11101817
12.
Liming Ma, Xinhai Li, Tianqing Zhai, et al. Changes in the Habitat Preference of Crested Ibis (Nipponia nippon) during a Period of Rapid Population Increase. Animals, 2021, 11(9): 2626.
DOI:10.3390/ani11092626