Apparent annual survival of staging ruffs during a period of population decline: insights from sex and site-use related differences

The ruff Philomachus pugnax, a lekking shorebird wintering in Africa and breeding across northern Eurasia, declined severely in its western range. Based on a capture-mark-resighting programme (2004–2011) in the westernmost staging area in Friesland (the Netherlands), we investigated changes in apparent annual survival in relation to age and sex to explore potential causes of decline. We also related temporal variation in apparent survival to environmental factors. We used the Capture-Mark-Recapture multievent statistical framework to overcome biases in survival estimates after testing for hidden heterogeneity of detection. This enabled the estimation of the probability to belong to high or low detectability classes. Apparent survival varied between years but was not related to weather patterns along the flyway, or to flood levels in the Sahel. Over time, a decline in apparent survival is suggested. Due to a short data series and flag loss in the last period this cannot be verified. Nevertheless, the patterns in sex-specific detectability and survival lead to new biological insights. Among highly detectable birds, supposedly most reliant on Friesland, males survived better than females ( $$\bar{\phi }_{\text{HD males}}$$ ϕ ¯ HD males  = 0.74, range 0.51–0.93; $$\bar{\phi }_{\text{HD females}}$$ ϕ ¯ HD females  = 0.51, range 0.24–0.81). Among low detectable birds, the pattern is reversed ( $$\bar{\phi }_{\text{LD males}}$$ ϕ ¯ LD males  = 0.64, range 0.37–0.89; $$\bar{\phi }_{\text{LD females}}$$ ϕ ¯ LD females  = 0.73, range 0.48–0.93). Probably the staging population contains a mixture of sex-specific migration strategies. A loss of staging females could greatly affect the dynamics of the western ruff population. Further unravelling of these population processes requires geographically extended demographic monitoring and the use of tracking devices.


Introduction
Increasing human pressures on habitats are driving many migratory species in decline (Berger 2004;Sanderson et al. 2006;Wilcove and Wikelski 2008;Limburg and Waldman 2009). Despite their great ability to move, due to specific habitat requirements, tight time schedules and density dependent processes, long-distance migrants are vulnerable to global change (Alerstam et al. 2003;Piersma 2007;Taylor and Norris 2007). An understanding of the resilience of migratory species, and the design of effective conservation measures for their protection, requires the identification of the reasons for population change. An essential step is to address the variation in demographic processes (survival, recruitment and movements) determining population growth rate (Caswell, 2001). The challenge to do this properly is considerable when studying long-distance migrants because: (i) the variation in demographic processes that may be related to immediate or delayed effect(s) of event(s) encountered in separated areas and times during the species migratory cycle (Harrison et al. 2011;Juillet et al. 2012) and (ii) the degree of migratory connectivity will influence the population dynamics (Webster et al. 2002;Taylor and Norris 2010). Capture-Mark-Recapture (CMR) methods based on the long-term monitoring of marked individuals constitute a versatile approach to assess population change. Taking into account imperfect detection of marked individuals, it allows the estimation of demographic rates in a robust way and provides a basis for inference on potential drivers (Lebreton et al. 1992;Gimenez et al. 2008).
The Ruff (Philomachus pugnax) is a long-distance migrating shorebird found across the African and Eurasian continents (van Rhijn 1991). Since the late 20th century, Ruffs showed a strong population decline over their westernmost migration route, the route that connects the inland wetlands of West-Africa to the wet grasslands of northwestern Europe and the tundras north and east of it (Zöckler 2002a;Rakhimberdiev et al. 2010;Verkuil et al. 2012). Breeding Ruffs first vanished from temperate latitudes partly because of the loss of wet grasslands habitat following European agriculture intensification (Thorup 2006), but their decline also spread into the sub-arctic region of Scandinavia (Väisänen et al. 2005;Øien and Aarvak 2010;Lindström and Green 2013). The latter suggests that global factors, including climate warming (Zöckler 2002b;Ims et al. 2008; Virkkala and Rajasärkkä 2011) may play a role. At the main staging area of 'western' Ruffs, in the Netherlands, peak roost counts of the spring passage population have declined from > 20000 in the late 1990s to fewer than 5000 in 2010 (Verkuil et al. 2012, Fig. 1b). In concert with a decline in numbers, the daily body mass increments of staging Ruffs decreased between 2001 and 2008, suggesting a decrease in staging habitat quality as a driver of this loss (Verkuil et al. 2012).
Changes occurring on the African nonbreeding grounds, including the repeated droughts, high hunting pressure and increasing human pressure in the Sahel region could be involved as well (Zwarts et al. 2009). Recent studies suggested the occurrence of a massive redistribution towards a more eastern flyway which would have contributed to the decline in numbers of western Ruffs (Rakhimberdiev et al. 2010;Verkuil et al. 2012).
The demographic processes underlying Ruff population changes remain to be investigated, and given the strong sexual dimorphism (males are much larger, and more ornamented, than females) this requires a sex-specific approach. The Ruff is a lekking species (Hogan-Warburg 1966;Widemo 1998) in which males and females differ not only by size but also by their ecology and migratory behaviour (van Rhijn 1991). The demographic processes and density-dependent processes relative to the spatio-temporal distributions of the two sexes are important determinants of their population dynamics (van Rhijn 1991;Höglund 1996;Anthony and Blumstein 2000;Gascoigne et al. 2009).
As a first step in the investigation of the demography of the western migrant Ruff population, we focused on survival probability. For shorebirds species, which are relatively long-lived with an early maturity but a highly variable recruitment due to the unpredictability of their breeding environments especially at high lati tude (see e.g., Aharon-Rotman et al. 2015), adult survival is likely an important determinant of population growth rate (Saether and Bakke 2002). Based on a capture-mark-resighting (CMR) dataset collected during eight successive spring migration seasons (March-April 2004-2011 in the staging area of the Nether lands in the province of Friesland, we used the CMR statistical framework to investigate annual variation in Ruff "apparent" survival probability (i.e., product of survival and fidelity) as a function of sex and age. Further, we explored relationships between annual variation in apparent survival and environmental conditions encountered during the entire annual cycle along the East-Atlantic migratory flyway.

Field methods and data collection
Our study area in southwest Friesland, the Netherlands, is situated along the northeast shores of Lake IJsselmeer and covers over 100 km 2 of agricultural habitat, mainly intensively managed grasslands for dairy farming along with traditionally managed wet herb-rich grasslands (Verkuil and de Goeij 2003;Groen et al. 2012, Fig.1a). The passage of northward migrating Ruffs starts in mid-March and lasts into May. From the 2nd week of March numbers build up to peak in the first half of April. Males stay for about three weeks (Verkuil et al. 2010), during which they undertake nuptial moult and gain body mass for the flight towards their breeding grounds (Jukema and Piersma 2000;Verkuil et al. 2010). Migration of subadults (first year of age) starts in early April and females (called 'reeves') only appear in mid-April (Wymenga 1999). Females are also much less numerous than males and due to insufficient resightings data, their staging duration could not be determined (see Verkuil et al. 2010). By the second week of May all Ruffs have left the staging area.
Ruffs were captured with traditional, partially wind-driven, "wilsternets" (equivalent to a 20 m long and 3 m high clap net; Piersma et al. 2005). Age was determined from leg colour [grey/green in subadults (first year bird); orange/ pink in adults] and the presence of brown (buff) middle cover tertiary feathers retained from juvenile plumage (Glutz von Blotzheim et al. 1975). As females are a third smaller than males, most birds could easily be assigned a sex, while the faeders, female-mimicking males, were discriminate using wing length (Jukema and Piersma 2006). At last, each captured Ruff was marked with a metal ring and a unique combination of 4 colour-rings and a flag before release.
Observations of colour-ringed individuals were made with telescopes (20-60×65 mm and 20-60×80 mm telescopes) from cars parked along country roads between agricultural fields or in nature reserves. Each spring, two teams of two persons were looking for colour-ringed individuals for up to 10 h a day, six days a week. The southern, central and northern parts of our study site were visited alternatively in order to achieve complete coverage every two days.
Our observation effort spread during the entire passage period (i.e., 15 March -15 May). In this way we guaranteed a good and consistent coverage of the passage population throughout the years of study. Staging Ruffs were not evenly spread across the study area, and some subareas were used more than others, in particular coastal polders close to roost sites (see Fig. 3.1A). Theses sub-areas necessarily required more time. However, as Ruffs are very mobile, turnover of individuals is usually high. Also, all observations for which the colour-ring combination was not identified with assurance were excluded from the analyses.
Between 2004 and 2010 we marked 5026 Ruffs. Among them were 4210 males of which 659 (15 %) were marked as subadult and 3551 as adult; and 816 females of which 157 (18 %) were subadult and 659 adult. From 2004From -2011 % of the marked females (n = 187) and 42 % of the marked males (n = 1769) were seen at least once after release. Tag loss (loss of the flag) was very limited until 2009. In 2010, the loss of flag became more common and we started recording it carefully. The 14 cases individuals without a flag were identified by reading their metal ring with a telescope indicated that flags were not lost immediately after ringing, but suggest a gradual loss. Degradation or heterogeneous quality of the plastics may be involved. Loss of colour-rings hardly ever occurred. Observations of birds which had lost their coloured flag comprised about 5 % (100 observations of marked Ruffs with no flag on a total of 2138 observations) and 10 % (206 observations of marked Ruffs with no flag on a total of 2146 observations) of our sightings in the spring seasons of 2010 and 2011, respectively. However as a same bird could be observed more than one time, the proportion of birds that have lost their flag must be lower. In the absence of physical recaptures to individually verify the loss of flags using information on the metal rings, tag-loss cannot be easily handled (see Juillet et al. 2010). We must therefore acknowledge that our estimates of apparent survival are likely negatively biased especially at the end of the study (see Arnason and Mills 1981).

Statistical modelling
We defined four groups: males and females, either marked as adults or marked as subadults. The encounter history of each individual was coded '1' when the individual was either captured or subsequently observed during a given encounter occasion (i.e., here a spring migration period); '0' when the individual was not observed. Significant trap-dependence, in this case 'trap-happiness' (i.e., individuals recently captured or observed were more likely to be re-observed subsequently -overall Test 2.CT, P < 0.001), was detected particularly among males (Test 2.CTmales marked as subadult, P = 0.018; marked as adult, P = 0.001). The test for transience was significant for adult males (i.e., presence of transients: individuals marked and never re-observed -test 3.SR, P = 0.016). In females, all test statistics were non-significant, but the sparse amount of data could have limited the power to detect such biases (see Pradel 1997). Results of GOF tests indicated overall non-homogeneous detection (Crespin et al. 2008). Such hidden heterogeneity in detection among marked Ruffs may be related to our sampling or linked to the behaviour of the birds.
To account for uncertainty on heterogeneity of detection and overcome biases in survival estimates, we adopted multievent models (Pradel 2005;Crespin et al. 2008; Pradel and Sanz-Aguilar 2012). We considered two live states, one being individuals alive and at high probability of detection (highly detectable, HD) and the other, individuals alive at low probability of detection (low detectable, LD). A third state is reserved for dead individuals that cannot be observed. Upon initial capture, each individual has a probability (π ) to belong to a given detectability class. This probability relates to the overall encounter history of the individuals, depending on whether individuals have been frequently or sporadically resighted across the study period. Each spring migration period, individuals can be either observed or not (i.e., the events) with a detection probability ( p) conditional on the underlying states (alive and either of high or of low detectability). Thus, birds assigned to the high detectability class with a probability π HD , shared a high probability of detection p HD and, inversely, birds assigned to the low detectability class with a probability (1-π HD ), shared a low probability of detection pLD. Individuals were assumed to survive between years t and t+1 with a survival probability (φ) which could be state-specific (survival probabilities in state HD or LD). Transitions between states were not considered (i.e., no transitions between detectability classes).
Differences in detectability may be linked to unrecognized behavioural features of the birds, features which themselves might be related to their survivorship or fidelity to Friesland. For this reason we integrated detectability classes into the modelling of the apparent survival probability. We also defined two age classes: adults (age ≥ 2 years -noted ad) and subadults (first year -noted 1y). Birds marked as subadults may have different propensities to re-use Friesland and/or to survive the year following marking than adults. Subadults may only be "prospectors" on breeding grounds (e.g., Kentie et al. 2014). According to van Rhijn (1991) only a small proportion of males and females will reproduce when one year old. First year male Ruffs refuel less rapidly (Karlionova et al. 2008) and develop less elaborated nuptial plumage than adults, behave mostly as marginal males and are less involved in displays at leks and have lower mating chances (van Rhijn 1991). We thus considered age as an additive effect. In view of the distinctiveness of males and females, we consider an effect of sex in combination with age and detectability classes.
At last, our most complex model assumed a non-null transience probability for adult male τ (ad males ). The probability of membership of a given detectability class (π ) presented the same temporal trends and variation in adult males, females (sex ad ) and subadults π (sex ad + 1y + t). The apparent survival probabilities (φ) varied independently over years between high and low detectable individuals (H × t) while sex and age effects were modeled as additive. Adult males and females in both detectability class were allowed to have different intercepts (H × sex ad ), but no different trends over years (i.e., low number of female resightings).
Subadults were considered independently φ(H × (sex ad + t) +1y). The detection probabilities of high and low detectable individuals varied independently over years (H × t), the sex effect was modeled as additive with different intercepts for males and females in both detectability p(H ×(t + sex)), but no differences in temporal pattern. Our most complex model was then { τ (ad males ), π (sex ad + 1y + t), φ(H ×(sex ad + t) + 1y), p(H ×(sex + t)) } and constituted our starting point.
Keeping the apparent survival probability time dependent (our parameter of interest), we followed a step-down approach (Lebreton et al. 1992) and consecutively modelled the transience parameter, the probability of detection and the probability of membership in a given detectability class. At last we modelled apparent survival probabilities.

Environmental variables
We tested the influence of environmental covariates on temporal variation of apparent survival probabilities. We considered local weather and regional climate indices as covariates of the survival probability. Foraging conditions in Friesland are crucial for staging Ruffs to combine refuelling and the development of nuptial moult (Jukema et al. 2001). To fuel up, Ruffs rely on soil-dwelling invertebrates (van Rhijn 1991; Verkuil and de Goeij 2003; J. Onrust, unpublished data) made less available by the intensive management of grasslands (Vickery et al. 2001). In absence of comprehensive data on the yearly variations in invertebrate prey availability for different land use categories, we assumed that in this intensively managed landscape foraging conditions for Ruffs will mainly be modulated by weather conditions. Probing for and finding soil-dwelling invertebrates will get harder as the vegetation grows and as the soil dries during the spring season. Hence, we assumed that both dry and early springs might negatively affect refuelling rates and thereby influence survival. Growing Degree Day (GDD), the sum of average daily temperature relative to a temperature threshold of 0 °C , has been shown to be a relevant measure of the onset of spring growth of graminoids in the Netherlands (van Wijk et al. 2012). Here we used the cumulative values of GDD from 1 January to 15 March as a measure of seasonal advancement at the moment that Ruffs started to stage in the study area. To characterize wet and dry springs we used the total amount of precipitation during spring stopover period (from 15 March to 15 May). Both precipitation and temperature data were obtained from Stavoren (52° 52ʹ 56.8₺ N; 5°21ʹ 58.0₺ E) situated in the south of our study area (www.knmi.nl).
To indicate the regional weather pattern over most of the East Atlantic migration flyway, we used a principal component (PC) based North Atlantic Oscillation (Hurrell 2003), provided by NCAR's Climate Analysis Section (http://climatedataguide.ucar.edu/guidance/hurrell-north-atlantic-oscillationnao-index-pc-based). Upon arrival on the breeding grounds, snow cover will make food scarce or unavailable increasing risk of starvation (Meltofte et al. 2007;Buehler and Piersma 2008). We used NAO average values from December through March in the winter preceding breeding (pre-breeding NAO or PBNAO) as a measure of spring snow melt conditions (Stenseth et al. 2003;Kausrud et al. 2008). A small number of Ruffs does not reach Africa but remains in Europe in winter, mainly in the Netherlands, Belgium and the United Kingdom (Prater 1973;Castelijns et al. 1994;Girard 2009;L.E. Schmaltz, unpublished data). These birds could represent a non-negligible proportion of the remaining spring passage population of Ruffs in Friesland. Thus, we included winter NAO values from October to January reflecting mild and wet winter (positive NAO) vs. cold and dry winter (negative NAO) over northern Europe (Stenseth et al. 2003).
For Ruffs wintering in the Sahelian floodplains of West-Africa, dry winters will cause small extents of open water areas to shrink more rapidly and hence force birds to concentrate on fewer feeding and roosting areas (Kone et al. 2002). During the pre-migratory period, this lead to higher competition and hunting pressure (Zwarts et al. 2009). As a covariate we considered the maximum flood extent calculated from aerial pictures made every winter from 2004-5 to 2010-11 from the inner Niger Delta floodplain (L. Zwarts, personal communication), the main wintering area for western migrant Ruffs.
We used principal component analysis (PCA) of the 5 environmental variables to combine their variation in uncorrelated, so-called 'synthetic', covariates which allowed us to reduce the number of environmental covariates and prevent collinearity issues (

Results
Model selection on the transience parameter (τ ), the probability of membership of a given detectability class (π ) and the probability of detection (p).
A transience effect did not improve the fit of our general model to the data and thus was not retained (Table S1 in Electronic Supplementary Material (ESM)). The probability of membership in a high detectability class varied between years, in parallel in subadult and adult males and females over the years (Table S1 in ESM, model selection on π , model 1 π (sex ad + 1y + t)). Subadults had the highest probability to be assigned to the high detectability class, then adult males and finally adult females (Fig. 3.2A). In the multievent approach, the probability to belong to the high detectability class provide an estimate of the actual proportion of high detectable individuals within newly marked birds (Pradel 2005). The probability of membership in a high detectability class showed an increase with  Table  3.1). The error bars show 95 % confidence intervals. time (Fig. 3.2A), which would indicate an increase of the fraction of highly detectable birds within newly captured individuals in the course of this study. However, we acknowledge that this interpretation has to be taken with prudence. As our study is relatively short, intervals of confidence become rapidly larger (i.e., the capture histories of newly marked birds getting shorter) and assignment to either detectability class less reliable with the years.
Detection probabilities varied with years in the same way between males and females (Table S1 in ESM, model selection on p, model 3 p(H + sex + t)), so that detection probabilities were much higher for males ( p -HD = 0.68, range 0.50 -0.75; p -LD = 0.21, range 0.18 -0.33; Fig. 3.2B) than for females ( p -HD = 0.47, range 0.27 -0.62; p -LD = 0.08, range 0.04 -0.15; Fig. 3.2B). In both sexes, the difference in detection probabilities between the two detectability classes was relatively important.

Model selection on the apparent survival probability (φ)
The best fitted parameterization on apparent survival probability indicated interannual variation (Table 3.1, model 8 φ(H × sex + t)). Apparent survival estimates decrease strongly at the end of the study period ( Fig. 3.2, C and D). A linear trend in survival probability was not supported (Table 3.1). Individuals at high and low detectability exhibited quantitative differences in survival, specific in opposite way for each sex (i.e., different intercepts -H × sex). According to the best supported model, the apparent survival probability of males at high detectability was higher than for males at low detectability (Fig. 3.2Dφ -HD males = 0.74, range: 0.51 -0.93; φ -LD males = 0.64, range: 0.37 -0.89), but the inverse was found for females ( Fig. 3.2Cφ -HD females = 0.51, range 0.24 -0.81; φ -LD females = 0.73, range 0.48 -0.93). We only found slight support for age differences in survival. The second best model distinguishing subadults from adult males and females fitted equally well (model 8 vs. model 6, ∆QAIC c = 0.88, Table 3.1).
That variation in apparent annual survival between individuals of high and low detectability were independent was not strongly supported (model 8 vs. model 3, ∆QAIC c ≥ 4.37, Table 3.1). Models considering no difference in survival between individuals at high and low detectability poorly represented our data (∆QAIC c ≥ 10.79, Table 3.1). Finally, our best fitting model (model 8 { π (sex ad + 1y + t), φ(H × sex + t), p(H + sex + t) }, Table 3.1), which included heterogeneity in detection probability, received much better support than an equivalent model under a traditional Cormack-Jolly-Seber model accounting for trapdependence on males (∆QAIC c = 86.21, Table 3.1).  Table S1 in ESM), detection probability ( p) and the probability of membership in the high detectability class (π) were π(sex ad + 1y + t) and p(H + sex + t) and do not vary among models. The overdispersion coefficient is: c = 1.263.

Environmental variables
In the analysis of environmental variables, we retained the two first principal components from the PCA (Eigen values > 1, see Fig. S2 in ESM). The first and second principal component (F1, F2) accounted for 50 % and 30 % of the total variation in the 5 covariates, respectively. Pre-breeding NAO index, winter NAO index and Growing Degree Day were positively correlated and contributed the most to F1 (see Table S2 and Fig. S2 in ESM). Positive values of F1 were thus associated to an early onset of spring and early snow melt on the breeding grounds, and wet and mild winter conditions over northwestern Europe. F2 was positively associated to wet spring conditions on the staging site, but lower flood level in Niger during winter (only relatively to the study period - Table S2 and Fig. S2 in ESM). Hereafter, F1 and F2 were noted respectively as GDD spr /NAOs and PP spr /Flood. Incorporation of the synthetic covariates in our best model (model 8) did not improve the fit to our data and actually received less support than a 'constant' model (Table 3.2). A model which combined the additive effects of both F1 and F2 did not provide better fits.  The detection probability ( p) and the probability of membership in the high detectability class (π) were π(sex ad + 1y + t) and p(H + sex + t) and do not vary among models (cf. Table S1 in ESM). The overdispersion coefficient is: c = 1.263.
Abbreviations: individual detectability class (H), year effect (t); GDD spr /NAOs and PP spr /Flood are our first and second PCA axis based on our 5 selected covariates. Number of parameter (K), Akaike Information Criterion corrected for overdispersion and sample size (QAIC c ), difference in QAIC c score to the best supported model (∆QAIC c -best model in bold), QAIC c weights (w i ).

Discussion
This study examined the year to year variation in apparent survival of western migrant Ruffs of different age and sex staging in Friesland, the Netherlands, during northward migration. Using a capture-mark-resight approach, we found evidence for heterogeneity in the detection probability among colour-marked Ruffs. We used multievent models integrating high and low detectability classes which allowed us to overcome bias and to investigate the biological meaning of heterogeneous detection. Since 2004 apparent survival probability of western migrant Ruffs varied between years. Overall, apparent survival differed according to sex and the detectability of individuals. The males in the high detectability class had a higher apparent survival than the males in the low detectability class, whereas in the females this pattern was reversed. Annual variation in apparent survival probability of western Ruffs could not be explained by the local spring weather conditions, by the broad climate index across the flyway, or by the Inner Niger Delta flood level. During our study, the flood levels of the Inner Niger Delta were relatively high compared to levels during the droughts in the 1980s (L. Zwarts, personal communication) which were responsible of mass mortality events (Zwarts et al. 2009). This might explain why we found no relationship between flood extent and apparent survival within our study period. Factors explaining variation in apparent survival thus remain unknown. The general insights on human-induced habitat changes causing the decline in Afro-Palearctic migrants are there (Vickery et al. 2014), but to fully establish these relationships we need more targeted ecological data.
Survival varied between years and suggest a decline over time. We acknowledge that this pattern to some extent was linked to flag loss. However, the loss of flags likely happened gradually, as the plastic wears out with time. A higher mortality and/or that less bird return to Friesland in the end of our study period may be also involved. The latter would be consistent with previous studies suggesting redistribution of western Ruffs towards more eastern migratory routes (Rakhimberdiev et al. 2010;Verkuil et al. 2012).
During the analysis of apparent annual survival it was necessary to allow for detectability differences between groups of animals. This significantly improved the fit of the models and also uncovered biologically interesting patterns. There are several ideas about how behaviour may affect detectability: (i) the ranging pattern of individuals relative to the study area may consistently differ between birds (e.g., Fearnbach et al. 2012). Although our observation effort covers the core staging area of Ruffs in Friesland, individuals with 'home ranges' edging our study site would have smaller chances to be resighted than birds foraging exclusively within its limits. (ii) the overall distribution of staging Ruffs may have con-centrated in certain areas (see Fig. 1a, L.E. Schmaltz, unpublished data). Therefore the Ruffs were easier to find for the observers thereby increasing overall detectability of the birds. Such interactions between ecologically relevant behaviours of the birds and our way to detect them could explain an increase of the fraction of highly detectable birds within newly captured individuals during the study. (iii) intrinsic behavioural features such as migration strategies may differ between staging Ruffs. In long-distant migrants, a staging area typically harbours a mixture of individuals with different wintering area, timing of migration and breeding latitudes (e.g., Duijns et al. 2012;Franks et al. 2012), and/or staging duration (e.g., Dierschke et al. 2004). For instance, in a population of western gray whales (Eschrichtius robustus) monitored from feeding grounds upon migration (i.e., equivalent to a staging area), detectability differences between individuals were related to their residence time in the study area (Bradford et al. 2006). Our three explanations for detection heterogeneity are non-mutually exclusive.
We conclude that detectability difference should overall be related to the degree of site use, either in space (relatively to the study area) and/or time (fidelity, staging duration).
Among the birds the most reliant on Friesland (i.e., the class of highly detectable individuals), for which apparent survival probability approaches most true survival, adult females survive much less than adult males. In accordance, between 2004 and 2011, the numbers of females captured decreased faster (with 53 %, see Fig. 3.1C) than the numbers of males (with 42 %). A lower female survival is unexpected in Ruffs, as in many other polygynous vertebrate species males usually have a higher mortality than females (Trivers 1972). Together with a loss of females, this hints at female-specific vulnerability to conditions encountered somewhere along the western route. Zwarts et al. (2009) warned, for instance, that the tens of thousands of Ruffs hunted in Mali for commercial purpose are mainly females that face higher risk due to their delayed departure compare to males coinciding with lower water levels and therefore higher hunting pressure (Kone et al. 2002). A selective loss of females would have severe repercussion on recruitment, while perturbation of local operational sex ratio could also expose the population of western Ruffs, already much reduced in size, to higher risk of local extinction (Courchamp et al. 2008;Lee et al. 2011).
Interestingly, unlike females, males in the high detectability class, which we believe are most reliant on Friesland, had a higher apparent survival than the low detectable males that would be less reliant on Friesland. These inverse survival differences according to detectability and sexes are intriguing but the lack of information outside our study area limit further interpretation. Nevertheless it sug gests that males and females with diverse migratory patterns and fidelity to Fries land, might contribute differently to the population dynamics of western Ruffs.
Staging areas as Friesland are the crossroads where very different individuals with respect to origin, destination, behaviour and fidelity pass through. To monitor the survival of a migratory species from a single staging area is certainly challenging: emigration and mortality will always be confounded and the staging popu lation might be very heterogeneous. Our results point at sex-specific behavioural variability among Ruffs staging in Friesland associated with different apparent survival probabilities. Future work is needed to better characterize the apparent survival of these staging Ruffs and investigate the influence of their migratory features. The low survival of highly detectable females and their concomitant decrease in numbers in Friesland might greatly affect the population dynamics of western Ruffs and are in need to be investigated further. To unravel population processes at play, efforts should be made to collect data on female selective harvesting in the Sahel, to study reproductive success and widen our knowledge on migratory connectivity, inter-flyway movements and sex-specific dispersal aptitudes. To achieve such goals, the CMR monitoring should be extended to multiple staging sites, combined to a multi-year satellite telemetry study conducted in both males and females as well as the establishment of throughout breeding monitoring. As Ruffs are a very opportunistic species highly dependent on shallow waterbodies and soft substrates during the entire life cycle (van Rhijn 1991), all efforts to preserve these habitat characteristics throughout its migratory flyway will benefit the species. Such large scale and integrative approach should be a priority to protect and manage this widespread migrating species with weak connectivity in times of rapid global change. For clarity, all model parameters are not presented in each step. Overall model selection is performed starting from our initial general model, and successively a parameter at a time, starting next step with best supported constraint (in bold) on previously reduced parameter. Overdispersion coefficient is: c = 1.263.