Exploiting uncertain ecological fieldwork data with multi-event capture--recapture modelling: an example with bird sex assignment

Interspecific synchrony on breeding performance and the role of anthropogenic food subsidies

Synchrony can have important consequences for long-term metapopulations persistence, community dynamics and ecosystems functioning. While the causes and consequences of intra-specific synchrony on population size and demographic rates have received considerable attention only a few factors that may affect inter-specific synchrony have been described. We formulate the hypothesis that food subsidies can buffer the influence of environmental stochasticity on community dynamics, disrupting and masking originally synchronized systems. To illustrate this hypothesis, we assessed the consequences of European policies implementation affecting subsidy availability on the temporal synchrony of egg volume as a proxy of breeding investment in two sympatric marine top predators with differential subsidy use. We show how 7-year synchrony appears on egg volume fluctuations after subsidy cessation suggesting that food subsidies could disrupt interspecific synchrony. Moreover, cross correlation increased after subsidy cessation and environmental buffering seems to act during synchronization period. We emphasize that subsidies dynamics and waste management provide novel insights on the emergence of synchrony in natural populations.

An open spatial capture-recapture model for estimating density, movement, and population dynamics from line-transect surveys

The purpose of many wildlife population studies is to estimate density, movement, or demographic parameters. Linking these parameters to covariates, such as habitat features, provides additional ecological insight and can be used to make predictions for management purposes. Line-transect surveys, combined with distance sampling methods, are often used to estimate density at discrete points in time, whereas capture-recapture methods are used to estimate movement and other demographic parameters. Recently, open population spatial capture-recapture models have been developed, which simultaneously estimate density and demographic parameters, but have been made available only for data collected from a fixed array of detectors and have not incorporated the effects of habitat covariates. We developed a spatial capture-recapture model that can be applied to line-transect survey data by modeling detection probability in a manner analogous to distance sampling. We extend this model to a) estimate demographic parameters using an open population framework and b) model variation in density and space use as a function of habitat covariates. The model is illustrated using simulated data and aerial line-transect survey data for North Atlantic right whales in the southeastern United States, which also demonstrates the ability to integrate data from multiple survey platforms and accommodate differences between strata or demographic groups. When individuals detected from line-transect surveys can be uniquely identified, our model can be used to simultaneously make inference on factors that influence spatial and temporal variation in density, movement, and population dynamics.

Differential survival throughout the full annual cycle of a migratory bird presents a life-history trade-off

Long-distance migrations are among the most physically demanding feats animals perform. Understanding the potential costs and benefits of such behaviour is a fundamental question in ecology and evolution. A hypothetical cost of migration should be outweighed by higher productivity and/or higher annual survival, but few studies on migratory species have been able to directly quantify patterns of survival throughout the full annual cycle and across the majority of a species' range. Here, we use telemetry data from 220 migratory Egyptian vultures Neophron percnopterus, tracked for 3,186 bird months and across approximately 70% of the species' global distribution, to test for differences in survival throughout the annual cycle. We estimated monthly survival probability relative to migration and latitude using a multi-event capture-recapture model in a Bayesian framework that accounted for age, origin, subpopulation and the uncertainty of classifying fates from tracking data. We found lower survival during migration compared to stationary periods (β = -0.816; 95% credible interval: -1.290 to -0.318) and higher survival on non-breeding grounds at southern latitudes (<25°N; β = 0.664; 0.076-1.319) compared to on breeding grounds. Survival was also higher for individuals originating from Western Europe (β = 0.664; 0.110-1.330) as compared to further east in Europe and Asia, and improved with age (β = 0.030; 0.020-0.042). Anthropogenic mortalities accounted for half of the mortalities with a known cause and occurred mainly in northern latitudes. Many juveniles drowned in the Mediterranean Sea on their first autumn migration while there were few confirmed mortalities in the Sahara Desert, indicating that migration barriers are likely species-specific. Our study advances the understanding of important fitness trade-offs associated with long-distance migration. We conclude that there is lower survival associated with migration, but that this may be offset by higher non-breeding survival at lower latitudes. We found more human-caused mortality farther north, and suggest that increasing anthropogenic mortality could disrupt the delicate migration trade-off balance. Research to investigate further potential benefits of migration (e.g. differential productivity across latitudes) could clarify how migration evolved and how migrants may persist in a rapidly changing world.

Decrease in social cohesion in a colonial seabird under a perturbation regime

Social interactions, through influence on behavioural processes, can play an important role in populations’ resilience (i.e. ability to cope with perturbations). However little is known about the effects of perturbations on the strength of social cohesion in wild populations. Long-term associations between individuals may reflect the existence of social cohesion for seizing the evolutionary advantages of social living. We explore the existence of social cohesion and its dynamics under perturbations by analysing long-term social associations, in a colonial seabird, the Audouin’s gull Larus audouinii, living in a site experiencing a shift to a perturbed regime. Our goals were namely (1) to uncover the occurrence of long-term social ties (i.e. associations) between individuals and (2) to examine whether the perturbation regime affected this form of social cohesion. We analysed a dataset of more than 3500 individuals from 25 years of monitoring by means of contingency tables and within the Social Network Analysis framework. We showed that associations between individuals are not only due to philopatry or random gregariousness but that there are social ties between individuals over the years. Furthermore, social cohesion decreased under the perturbation regime. We sustain that perturbations may lead not only to changes in individuals’ behaviour and fitness but also to a change in populations’ social cohesion. The consequences of decreasing social cohesion are still not well understood, but they can be critical for the population dynamics of social species.

Transience effect in capture-recapture studies: The importance of its biological meaning

Capture-recapture (CR) models are an essential tool for estimating demographic parameters in most animal and some plant populations. To avoid drawing incorrect conclusions in any statistical inference, a crucial prerequisite is to assess the goodness of fit of a general model to the data. In CR models, a frequent cause of lack of fit, is the so-called transience effect, which is due to a lower expectation of re-observation of individuals marked for the first time as compared to other individuals. The transience effect may result either from different biological causes or from the sampling procedure. A transience effect is usually treated by distinguishing at least two age-classes in the survival probability, but other approaches may be more suitable. Here we develop a conceptual and analytical framework for including a transience effect in capture-recapture models. We show the implementation of two additional parametrizations that incorporate a transience effect. With these parametrizations, we can directly estimate the "transience probability" defined as the probability that a newly caught individual disappear from the population beyond what is expected based on the behavior of the previously caught individuals in the same sample. Additionally, these parametrizations allow testing biological hypotheses concerning drivers affecting this probability. Results from our case study show differences between parametrizations, with the parametrization most currently used giving different estimates, especially when including covariates. We advocate for a unifying framework for treating a transience effect, that helps clarifying the ideas and terminology, and where the biological reasons should be the rule for choosing the appropriate analytical procedure. This framework will also open new and powerful ways for the detection and exploration of ecological processes such as the costs of the first reproduction or the deleterious effects of some types of marking.

Estimating age-dependent survival when juveniles resemble females: Invasive ring-necked parakeets as an example

Many species only show sexual dimorphism at the age of maturity, such that juveniles typically resemble females. Under these circumstances, estimating accurate age-specific demographic parameters is challenging. Here, we propose a multievent model parameterization able to estimate age-dependent survival using capture-recapture data with uncertainty in age and sex assignment of individuals. We illustrate this modeling approach with capture-recapture data from the ring-necked parakeet Psittacula krameri. We analyzed capture, recapture, and resighting data (439 recaptures/resightings) of 156 ring-necked parakeets tagged with neck collars in Barcelona city from 2003 to 2016 to estimate the juvenile and adult survival rate. Our models successfully estimated the survival probabilities of the different age classes considered. Survival probability was similar between adults (0.83, 95% CI = 0.77-0.87) and juveniles during their second (0.79, 95% CI = 0.58-0.87) and third winter (0.83, 95% CI = 0.65-0.88). The youngest juveniles (1st winter) showed a slightly lower survival (0.57, 95% CI = 0.37-0.79). Among adults, females showed a slightly higher survival than males (0.87, 95% CI = 0.78-0.93; and 0.80, 95% CI = 0.73-0.86, respectively). These high survival figures predict high population persistence in this species and urge management policies. The analysis also stresses the usefulness of multievent models to estimate juvenile survival when age cannot be fully ascertained.

Sex-biased survival contributes to population decline in a long-lived seabird, the Magellanic Penguin

We developed a Hidden Markov mark-recapture model (R package marked) to examine sex-specific demography in Magellanic Penguins (Spheniscus magellanicus). Our model was based on 33 yr of resightings at Punta Tombo, Argentina, where we banded ~44,000 chicks from 1983 to 2010. Because we sexed only 57% of individuals over their lifetime, we treated sex as an uncertain state in our model. Our goals were to provide insight into the population dynamics of this declining colony, to inform conservation of this species, and to highlight the importance of considering sex-specific vital rates in demographic seabird studies. Like many other seabirds, Magellanic Penguins are long-lived, serially monogamous, and exhibit obligate biparental care. We found that the non-breeding-season survival of females was lower than that of males and that the magnitude of this bias was highest for juveniles. Biases in survival accumulated as cohorts aged, leading to increasingly skewed sex ratios. The survival bias was greatest in years when overall survival was low, that is, females fared disproportionality worse when conditions were unfavorable. Our model-estimated survival patterns are consistent with independent data on carcasses from the species' non-breeding grounds, showing that mortality is higher for juveniles than for adults and higher for females than for males. Juveniles may be less efficient foragers than adults are and, because of their smaller size, females may show less resilience to food scarcity than males. We used perturbation analysis of a population matrix model to determine the impact of sex-biased survival on adult sex ratio and population growth rate at Punta Tombo. We found that adult sex ratio and population growth rate have the greatest proportional response, that is, elasticity, to female pre-breeder and adult survival. Sex bias in juvenile survival (i.e., lower survival of females) made the greatest contribution to population declines from 1990 to 2009. Because starvation is a leading cause of morality in juveniles and adults, precautionary fisheries and spatial management in the region could help to slow population decline. Our data add to growing evidence that knowledge of sex-specific demography and sex ratios are necessary for accurate assessment of seabird population trends.

Use of hidden Markov capture-recapture models to estimate abundance in the presence of uncertainty: Application to the estimation of prevalence of hybrids in animal populations

Estimating the relative abundance (prevalence) of different population segments is a key step in addressing fundamental research questions in ecology, evolution, and conservation. The raw percentage of individuals in the sample (naive prevalence) is generally used for this purpose, but it is likely to be subject to two main sources of bias. First, the detectability of individuals is ignored; second, classification errors may occur due to some inherent limits of the diagnostic methods. We developed a hidden Markov (also known as multievent) capture-recapture model to estimate prevalence in free-ranging populations accounting for imperfect detectability and uncertainty in individual's classification. We carried out a simulation study to compare naive and model-based estimates of prevalence and assess the performance of our model under different sampling scenarios. We then illustrate our method with a real-world case study of estimating the prevalence of wolf (Canis lupus) and dog (Canis lupus familiaris) hybrids in a wolf population in northern Italy. We showed that the prevalence of hybrids could be estimated while accounting for both detectability and classification uncertainty. Model-based prevalence consistently had better performance than naive prevalence in the presence of differential detectability and assignment probability and was unbiased for sampling scenarios with high detectability. We also showed that ignoring detectability and uncertainty in the wolf case study would lead to underestimating the prevalence of hybrids. Our results underline the importance of a model-based approach to obtain unbiased estimates of prevalence of different population segments. Our model can be adapted to any taxa, and it can be used to estimate absolute abundance and prevalence in a variety of cases involving imperfect detection and uncertainty in classification of individuals (e.g., sex ratio, proportion of breeders, and prevalence of infected individuals).

Immature survival, fertility, and density dependence drive global population dynamics in a long-lived species

Disentangling the influence of demographic parameters and the role of density dependence on species' population dynamics is a challenge, especially when fractions of the population are unobservable. Additionally, due to the difficulty of gathering data at large spatial scales, most studies ignore the global dynamic of a species, which would integrate local heterogeneity dynamics and remove the noise of dispersal. We developed an integrated population model (IPM) at a global scale to disentangle the main demographic drivers of population dynamics in a long-lived species. We used 28 yr of Audouin's Gull demographic data encompassing 69 local patches (comprising 90% of the world population). Importantly, we took into account the unobservable fraction of non-breeders and also assessed the strength of density dependence for this fraction of the population. As predicted by life histories of long-lived organisms, temporal random variation in survival was highest for immature individuals (1.326, 95% credible interval [CRI] 1.290-1.940) and lowest for adults (0.499, 95% CRI 0.487-0.720). Large temporal fluctuations in the probability of taking a reproductive sabbatical would partly explain the consistency in adult survival, with individuals most likely refraining from breeding when environmental conditions were harsh. Immature survival and fertility were the main drivers of population dynamics during the study period (r²  = 0.83, 0.77-0.87 and 0.73, 0.63-0.79, respectively). We found strong evidence of density dependence, not only due to the number of breeders (r²  = -0.34, -0.43 to -0.24) but also due to individuals on sabbatical (r²  = -0.18, -0.33 to -0.01). From a conservation point of view, the species shows a 5% annual global decrease during the last 10 years, and we propose an update of its conservation status. Even though population dynamics of long-lived organisms are very sensitive to changes in adult survival, we show here that, in the absence of strong environmental perturbations affecting this vital rate, fluctuations in population density are mainly driven by variations in survival of immature individuals and fertility. Integrated models based on long-term monitoring at a global scale may enhance our ecological and evolutionary understanding of how demographic drivers influence population dynamics.

Varying demographic impacts of different fisheries on three Mediterranean seabird species

Fisheries have an enormous economic importance, but reconciling their socio-economic features with the conservation and sustainability of marine ecosystems presents major challenges. Bycatch mortality from fisheries is clearly among the most serious global threats for marine ecosystems, affecting a wide range of top predators. Recent estimates report ca. 200,000 seabirds killed annually by bycatch in European waters. However, there is an urgent need to rigorously estimate actual mortality rates and quantify effects of bycatch on populations. The Mediterranean Sea is one of the most impacted regions. Here, we estimate for the first time both bycatch mortality rates and their population-level effects on three endemic and vulnerable Mediterranean taxa: Scopoli's shearwater, Mediterranean shag, and Audouin's gull, that die in different types of fishing gears: longlines, gillnets and sport trolling, respectively. We use multi-event capture-recapture modelling to estimate crucial demographic parameters, including the probabilities of dying in different fishing gears. We then build stochastic demography models to forecast the viability of the populations under different management scenarios. Longline bycatch was particularly severe for adults of Scopoli's shearwaters and Audouin's gulls (ca. 28% and 23% of total mortality, respectively) and also for immature gulls (ca. 90% of mortality). Gillnets had a lower impact, but were still responsible for ca. 9% of juvenile mortality on shags, whereas sport trolling only slightly influenced total mortality in gulls. Bycatch mortality has high population-level impacts in all three species, with shearwaters having the highest extinction risk under current mortality rates. Different life-history traits and compensatory demographic mechanisms between the three species are probably influencing the different bycatch impact: for shearwaters, urgent conservation actions are required to ensure the viability of their populations. Results will be very useful for guiding future seabird conservation policies and moving towards an ecosystem-based approach to sustainable fisheries management.

Consecutive cohort effects driven by density-dependence and climate influence early-life survival in a long-lived bird

Conditions during early life, including maternal cohort effects, can influence the future fitness of individuals. This may be particularly true for long-distance migrating birds, because, apart from conditions experienced by cohorts during rearing, conditions during early life in regions far from breeding grounds may also influence their population dynamics. Very little is known about the fitness consequences of those conditions experienced by juveniles after independence, especially in wild populations and for long-lived birds. We used multi-event capture-recapture-recovery models and a unique 26-year dataset for the Audouin's gull (Larus audouinii) to assess for the first time whether survival was influenced by early conditions, both during the rearing period (i.e. a maternal cohort effect potentially affected by density dependence) and the first winter (i.e. a cohort effect driven by climate when birds disperse to wintering grounds). Our results show that juvenile survival was highly sensitive to early-life conditions and that survival decreased with stronger density dependence and harsh climate. The two consecutive cohort effects were of similar magnitude and they may represent a selection filter. Thus, early-life conditions had a strong impact on survival, and neglecting this complexity may hinder our understanding on how populations of long-lived animals fluctuate and respond to perturbations.

Multi-event capture-recapture modeling of host-pathogen dynamics among European rabbit populations exposed to myxoma and Rabbit Hemorrhagic Disease Viruses: common and heterogeneous patterns

Host-pathogen epidemiological processes are often unclear due both to their complexity and over-simplistic approaches used to quantify them. We applied a multi-event capture-recapture procedure on two years of data from three rabbit populations to test hypotheses about the effects on survival of, and the dynamics of host immunity to, both myxoma virus and Rabbit Hemorrhagic Disease Virus (MV and RHDV). Although the populations shared the same climatic and management conditions, MV and RHDV dynamics varied greatly among them; MV and RHDV seroprevalences were positively related to density in one population, but RHDV seroprevalence was negatively related to density in another. In addition, (i) juvenile survival was most often negatively related to seropositivity, (ii) RHDV seropositives never had considerably higher survival, and (iii) seroconversion to seropositivity was more likely than the reverse. We suggest seropositivity affects survival depending on trade-offs among antibody protection, immunosuppression and virus lethality. Negative effects of seropositivity might be greater on juveniles due to their immature immune system. Also, while RHDV directly affects survival through the hemorrhagic syndrome, MV lack of direct lethal effects means that interactions influencing survival are likely to be more complex. Multi-event modeling allowed us to quantify patterns of host-pathogen dynamics otherwise difficult to discern. Such an approach offers a promising tool to shed light on causative mechanisms.

Man-induced activities modify demographic parameters in a long-lived species: effects of poisoning and health policies

Recent changes in sanitary policies within the European Union (EU) concerning disposal of carcasses of domestic animals and the increase of non-natural mortality factors, such as illegal poisoning, are threatening European vultures. However, the effects of anthropogenic activities on demographic parameters are poorly studied. Using a long-term study (1994-2011) of the threatened Pyrenean Bearded Vulture Gypaetus barbatus population, we assess the variation in the proportion of breeding pairs, egg-laying dates, clutch size, breeding success, and survival following a sharp reduction in food availability in 2005 due to the application of restrictive sanitary policies decreasing livestock carcass availability. We found a delay in laying dates and a regressive trend in clutch size, breeding success, and survival following policy change. The maintenance of specific supplementary feeding stations for Bearded Vultures probably reduced the negative effects of illegal poisoning and food shortages, which mainly affected subadult survival. A drop in food availability may have produced changes in demographic parameters and an increase in mortality due to an increased exposure to contaminated food. As a result, supplementary feeding as a precautionary measure can be a useful tool to reduce illegal poisoning and declines in demographic parameters until previous food availability scenarios are achieved. This study shows how anthropogenic activities through human health regulations that affect habitat quality can suddenly modify demographic parameters in long-lived species, including those, such as survival, with high sensitivity to population growth rate.

Environmental instability as a motor for dispersal: a case study from a growing population of glossy ibis

Dispersal is a life-history trait directly affecting population dynamics and species range shifts and thus playing a prominent role in the response to climate change. Nonetheless, the relationship between extreme climatic events and dispersal has received little attention in birds. Here we focused on climatic, demographic and individual factors affecting the dispersal propensity of a major glossy ibis population. We performed a capture-resighting analysis on individuals born and observed at Doñana (South-West Spain) over fourteen years. We applied a multiple analytical approach to show that single-site capture-resighting estimates were a reliable index of dispersal propensity from the area. We focused on the emigration of Doñana-born individuals sporadically (transients) and regularly (residents) frequenting their natal area. Droughts during two out of 14 study years caused higher apparent dispersal rates, explaining most of the annual variation in these rates. The age structure of Doñana-born individuals resighted simultaneously locally and in Morocco in one week over the 2010 autumn confirmed that the 2005 drought boosted permanent emigration. As numbers increased steadily during non-drought years since the formation of the colony in 1996 to several thousand pairs, philopatry increased gradually, while transients probability appeared to be related to average breeding success. Age, sex, density, quality of foraging habitat and breeding success in the previous season were not found to directly affect apparent dispersal. Nonetheless, autumn sex ratio gradually switched from male (≈0.68) to female-skewed (≈0.44) by the end of the study period, suggesting that males and females respond differently to high densities reached in recent years. This study demonstrates the importance of extreme climatic events as a powerful motor for spread of species in expansion. Also, it suggests different factors drive emigration of individuals according to their amount of experience in the area (e.g. transients vs residents).

Known unknowns in an imperfect world: incorporating uncertainty in recruitment estimates using multi-event capture-recapture models

Studying the demography of wild animals remains challenging as several of the critical parts of their life history may be difficult to observe in the field. In particular, determining with certainty when an individual breeds for the first time is not always obvious. This can be problematic because uncertainty about the transition from a prebreeder to a breeder state - recruitment - leads to uncertainty in vital rate estimates and in turn in population projection models. To avoid this issue, the common practice is to discard imperfect data from the analyses. However, this practice can generate a bias in vital rate estimates if uncertainty is related to a specific component of the population and reduces the sample size of the dataset and consequently the statistical power to detect effects of biological interest. Here, we compared the demographic parameters assessed from a standard multistate capture-recapture approach to the estimates obtained from the newly developed multi-event framework that specifically accounts for uncertainty in state assessment. Using a comprehensive longitudinal dataset on southern elephant seals, we demonstrated that the multi-event model enabled us to use all the data collected (6639 capture-recapture histories vs. 4179 with the multistate model) by accounting for uncertainty in breeding states, thereby increasing the precision and accuracy of the demographic parameter estimates. The multi-event model allowed us to incorporate imperfect data into demographic analyses. The gain in precision obtained has important implications in the conservation and management of species because limiting uncertainty around vital rates will permit predicting population viability with greater accuracy.

Revisiting the effect of capture heterogeneity on survival estimates in capture-mark-recapture studies: does it matter?

Recently developed capture-mark-recapture methods allow us to account for capture heterogeneity among individuals in the form of discrete mixtures and continuous individual random effects. In this article, we used simulations and two case studies to evaluate the effectiveness of continuously distributed individual random effects at removing potential bias due to capture heterogeneity, and to evaluate in what situation the added complexity of these models is justified. Simulations and case studies showed that ignoring individual capture heterogeneity generally led to a small negative bias in survival estimates and that individual random effects effectively removed this bias. As expected, accounting for capture heterogeneity also led to slightly less precise survival estimates. Our case studies also showed that accounting for capture heterogeneity increased in importance towards the end of study. Though ignoring capture heterogeneity led to a small bias in survival estimates, such bias may greatly impact management decisions. We advocate reducing potential heterogeneity at the sampling design stage. Where this is insufficient, we recommend modelling individual capture heterogeneity in situations such as when a large proportion of the individuals has a low detection probability (e.g. in the presence of floaters) and situations where the most recent survival estimates are of great interest (e.g. in applied conservation).