Return rates in studies of life history evolution: Are biases large?

Studies of life history evolution in passerine birds often depend on examination of annual survival probability of adult birds. Most studies rely on return rates (proportion of marked individuals released in one year that are recaptured in the next year) to estimate annual survival probability. Yet, return rate includes both the probability of survival and the probability of recapturing or resighting the bird in the next time interval. We use numerical estimation to illustrate the increasing bias in return rate as an estimator of annual survival probability as recapture/resighting probability decreases. Recapture/resighting probability is normally assumed to be high and relatively invariant for recapture/resighting studies of color-banded territorial birds. We tested this assumption through examination of 11 color-banding studies of passerines. These studies showed that recapture/resighting probabilities vary strongly and cannot be generalized as high. In short, return rates generally are poor estimators of annual survival probabilities and use of return rates may strongly bias relationships explored in comparative studies or bias results of experiments to test survival costs of reproduction. Recapture/resighting probabilities should be estimated in all studies that attempt to estimate annual survival probabilities.     

Estimating age-specific survival rates of tawny owls--recaptures versus recoveries

We compared estimates of annual survival rates of tawny owls ( Strix aluco ) ringed in southern Finland from several different sampling methods: recoveries of birds ringed as young; recaptures of birds ringed as young; recoveries of birds ringed as adults as well as young; combined recoveries and recaptures of birds ringed as young, and combined recoveries and recaptures of birds ringed as adults and young. From 1979 to 1998, 18 040 young owls were ringed, of which 983 were recaptured as breeders in subsequent years during this period, and 1764 were recovered dead at various locations. In addition, 1751 owls were ringed as adults, of which 612 were later recaptured and 199 were recovered dead. First-year survival rates estimated using only recoveries of birds ringed as young averaged 48%, while apparent survival rates estimated using only recaptures from birds ringed as young averaged 10-13%. Use of combined recapture-recovery models, or supplementary information from recoveries of birds ringed as adults, produced survival estimates of 30-37%. Survival estimates from young-recoveries-only models were biased high, because of violation of the assumption of constant recovery rates with age: birds dying in their first-year were one-third less likely to be found and reported than older birds. In contrast, recaptures-only models confounded emigration with mortality. Despite these differences in mean values, annual fluctuations in estimated first-year survival rates were similar with all models. Estimates of adult survival rates were similar with all models, while those for second-year birds were similar for all models except recaptures-only. These results highlight the potential biases associated with analysing either recaptures or recoveries alone of birds ringed as young, and the benefits of using combined data.     

Transient animals in a resident population of snow geese: Local emigration or heterogeneity?

The estimation of survival rates from analysis of recapture of individually marked animals assumes that all individuals are equally likely to be re-encountered. This assumption is frequently violated in natural populations due to movements to and from the sampling area. We evaluated potential sources of heterogeneity using data from recaptures of 36000 individually marked female lesser snow geese, Anser c. caerulescens , from an expanding population in northern Manitoba, Canada. By stratifying individuals according to marking age and origin (hatched at the colony or not), we assessed the degree to which variation in apparent survival reflected permanent or temporary differences in emigration and effects of handling. In general, for birds ringed as adults, estimated apparent survival rates were significantly lower during the first year after ringing than in subsequent years. By comparing birds ringed as adults (classified by origin) with those ringed as goslings, we were able to demonstrate that these differences are not due to permanent emigration from the colony by transient individuals or heterogeneity of individual capture probability, but more likely reflect differences among individuals in their response to initial marking. Approximately 25% of birds permanently emigrate from the sampling area following marking.     

Transient animals in a resident population of snow geese: local emigration or heterogeneity?

The estimation of survival rates from analysis of recapture of individually marked animals assumes that all individuals are equally likely to be re-encountered. This assumption is frequently violated in natural populations due to movements to and from the sampling area. We evaluated potential sources of heterogeneity using data from recaptures of 36000 individually marked female lesser snow geese, Anser c. caerulescens , from an expanding population in northern Manitoba, Canada. By stratifying individuals according to marking age and origin (hatched at the colony or not), we assessed the degree to which variation in apparent survival reflected permanent or temporary differences in emigration and effects of handling. In general, for birds ringed as adults, estimated apparent survival rates were significantly lower during the first year after ringing than in subsequent years. By comparing birds ringed as adults (classified by origin) with those ringed as goslings, we were able to demonstrate that these differences are not due to permanent emigration from the colony by transient individuals or heterogeneity of individual capture probability, but more likely reflect differences among individuals in their response to initial marking. Approximately 25% of birds permanently emigrate from the sampling area following marking.     

Survival estimation and the effects of dependency among animals

Survival models assume that fates of individuals are independent, yet the robustness of this assumption has been poorly quantified. We examine how empirically derived estimates of the variance of survival rates are affected by dependency in survival probability among individuals. We used Monte Carlo simulations to generate known amounts of dependency among pairs of individuals and analyzed these data with Kaplan-Meier and Cormack-Jolly-Seber models. Dependency significantly increased these empirical variances as compared to theoretically derived estimates of variance from the same populations. Using resighting data from 168 pairs of black brant ( Branta bernicla nigricans ), we used a resampling procedure and program RELEASE to estimate empirical and mean theoretical variances. We estimated that the relationship between paired individuals caused the empirical variance of the survival rate to be 155% larger than the empirical variance for unpaired individuals. Monte Carlo simulations and use of this resampling strategy can provide investigators with information on how robust their data are to this common assumption of independent survival probabilities.     

Survival estimation and the effects of dependency among animals

Survival models assume that fates of individuals are independent, yet the robustness of this assumption has been poorly quantified. We examine how empirically derived estimates of the variance of survival rates are affected by dependency in survival probability among individuals. We used Monte Carlo simulations to generate known amounts of dependency among pairs of individuals and analyzed these data with Kaplan-Meier and Cormack-Jolly-Seber models. Dependency significantly increased these empirical variances as compared to theoretically derived estimates of variance from the same populations. Using resighting data from 168 pairs of black brant ( Branta bernicla nigricans ), we used a resampling procedure and program RELEASE to estimate empirical and mean theoretical variances. We estimated that the relationship between paired individuals caused the empirical variance of the survival rate to be 155% larger than the empirical variance for unpaired individuals. Monte Carlo simulations and use of this resampling strategy can provide investigators with information on how robust their data are to this common assumption of independent survival probabilities.     

Costs of reproduction: Assessing responses to brood size manipulation on life-history and behavioural traits using multi-state capture-recapture models

Costs of reproduction are fundamental trade-offs shaping the evolution of life histories. There has been much interest, discussion and controversy about the nature and type of reproductive costs. The manipulation of reproductive effort (e.g. brood size manipulation) may alter not only life-history traits such as future adult survival rate and future reproductive effort, but also behavioural decisions affecting recapture/resighting and dispersal probabilities. We argue that many previous studies of the costs of reproduction may have erroneously concluded the existence or non-existence of such costs because of their use of local return rates to assess survival. In this paper, we take advantage of the modern multistate capture-recapture methods to highlight how the accurate assessment of the costs of reproduction requires incorporating not only recapture probability, but also behavioural 'state' variables, for example dispersal status and current reproductive investment. The inclusion of state-dependent decisions can radically alter the conclusions drawn regarding the costs of reproduction on future survival or reproductive investment. We illustrate this point by re-analysing data collected to address the question of the costs of reproduction in the collared flycatcher and the great tit. We discuss in some detail the methodological issues and implications of the analytical techniques.     

On the use of secondary capture-recapture samples to estimate temporary emigration and breeding proportions

The use of the Cormack-Jolly-Seber model under a standard sampling scheme of one sample per time period, when the Jolly-Seber assumption that all emigration is permanent does not hold, leads to the confounding of temporary emigration probabilities with capture probabilities. This biases the estimates of capture probability when temporary emigration is a completely random process, and both capture and survival probabilities when there is a temporary trap response in temporary emigration, or it is Markovian. The use of secondary capture samples over a shorter interval within each period, during which the population is assumed to be closed (Pollock's robust design), provides a second source of information on capture probabilities. This solves the confounding problem, and thus temporary emigration probabilities can be estimated. This process can be accomplished in an ad hoc fashion for completely random temporary emigration and to some extent in the temporary trap response case, but modelling the complete sampling process provides more flexibility and permits direct estimation of variances. For the case of Markovian temporary emigration, a full likelihood is required.     

On the use of secondary capture-recapture samples to estimate temporary emigration and breeding proportions

The use of the Cormack-Jolly-Seber model under a standard sampling scheme of one sample per time period, when the Jolly-Seber assumption that all emigration is permanent does not hold, leads to the confounding of temporary emigration probabilities with capture probabilities. This biases the estimates of capture probability when temporary emigration is a completely random process, and both capture and survival probabilities when there is a temporary trap response in temporary emigration, or it is Markovian. The use of secondary capture samples over a shorter interval within each period, during which the population is assumed to be closed (Pollock's robust design), provides a second source of information on capture probabilities. This solves the confounding problem, and thus temporary emigration probabilities can be estimated. This process can be accomplished in an ad hoc fashion for completely random temporary emigration and to some extent in the temporary trap response case, but modelling the complete sampling process provides more flexibility and permits direct estimation of variances. For the case of Markovian temporary emigration, a full likelihood is required.     

Suggestions for future directions for studies of marked migratory landbirds from the perspective of a practitioner in population management and conservation

Significant population declines in landbird species have been documented recently from many areas of the earth, including Europe and North America. Identification of the major causes of these declines and effective management actions to reverse them is difficult, especially for populations of long-distance migrants that winter in tropical areas. Key-factor and sensitivity analyses of critical population parameters in the context of integrated population models provide one promising approach to solving these problems. Key population factors may include breeding productivity, first-year survival, recruitment of young, adult survival and permanent emigration of adults; each of these can be indexed or estimated using data from cooperative ringing programmes, but the usefulness of the indices or estimates is limited by deficiencies in the available data and limitations of the available models. Future methodological directions for ringing studies should include efforts to: (1) develop and implement techniques to distinguish young from adult birds through the first breeding season of the young birds; (2) implement radio-tracking to determine characteristics of dispersal of young birds and transient adults; and (3) implement increased ringing, DNA fingerprinting and stable-isotope analysis to determine correspondence of breeding and winter ranges. Future programme-related directions should include efforts to: (1) integrate multiple methods at individual sites to compare and validate the indices and estimates produced by the different methods; (2) develop cooperative programmes of winter-season mist-netting to generate mark-recapture data to estimate the seasonal components of survival; and (3) develop mutually compatible banding programmes in tropical countries. Future theoretical and analytical directions should include efforts to continue to develop, refine and utilize: (1) key-factor and sensitivity analyses to determine the major causes of population changes; (2) models for dispersal of young birds and transient adults to improve the usefulness of indices of the number of hatch-year and second-year birds; (3) models to determine the proportions of transients in Cormack-Jolly-Seber (CJS) mark-recapture analyses and to eliminate their effects on estimates of survival rate, population size and recruitment of residents; (4) integrated models of population processes that utilize data from multiple methods to provide estimates of first-year survival, recruitment rate of young and permanent emigration rate of adults, parameters that are difficult to obtain from a single method; (5) models to estimate seasonal components of survival to provide insights into the timing and causes of mortality; (6) models incorporating environmental variables and species-specific characteristics as covariates in CJS mark-recapture and key-factor analyses; (7) models for pooling and weighting data obtained from multiple sites in cooperative ringing projects; (8) models for identifying long-term trends in demographic parameters; and (9) techniques for selection of appropriate models. Finally, assumptions implicit in the use of indices of various demographic parameters need to be tested and field techniques need to be improved to increase the numbers of individuals marked and recaptured in order to allow more precise parameter estimation; this will increase the ability to test competing hypotheses of population dynamics from data gathered in ringing programmes.     

Use of the Barker model in an experiment examining covariate effects on first-year survival in Ross's Geese ( Chen rossii ): A case study

The Barker model provides researchers with an opportunity to use three types of data for mark-recapture analyses - recaptures, recoveries, and resightings. This model structure maximizes use of encounter data and increases the precision of parameter estimates, provided the researcher has large amounts of resighting data. However, to our knowledge, this model has not been used for any published ringing studies. Our objective here is to report our use of the Barker model in covariate-dependent analyses that we conducted in Program MARK. In particular, we wanted to describe our experimental study design and discuss our analytical approach plus some logistical constraints we encountered while conducting a study of the effects of growth and parasites on survival of juvenile Ross's Geese. Birds were marked just before fledging, alternately injected with antiparasite drugs or a control, and then were re-encountered during migration and breeding in following years. Although the Barker model estimates seven parameters, our objectives focused on annual survival only, thus we considered all other parameters as nuisance terms. Therefore, we simplified our model structures by maintaining biological complexity on survival, while retaining a very basic structure on nuisance parameters. These analyses were conducted in a two-step approach where we used the most parsimonious model from nuisance parameter analyses as our starting model for analyses of covariate effects. This analytical approach also allowed us to minimize the long CPU times associated with the use of covariates in earlier versions of Program MARK. Resightings made up about 80% of our encounter history data, and simulations demonstrated that precision and bias of parameter estimates were minimally affected by this distribution. Overall, the main source of bias was that smaller goslings were too small to retain neckbands, yet were the birds that we predicted would have the lowest survival probability and highest probability for parasite effects. Consequently, we considered our results conservative. The largest constraint of our study design was the inability to partition survival into biologically meaningful periods to provide insight into the timing and mechanisms of mortality.     

Nonresponse Assessment of a Consumer Price Index

A household budget survey often suffers from a high nonresponse rate and a selective response. The bias that may be introduced in the estimation of budget shares because of this nonresponse can affect the estimate of a consumer price index, which is a weighted sum of partial price index numbers (weighted with the estimated budget shares). The bias is especially important when related to the standard error of the estimate. Because of the impossibility of subsampling nonrespondents to the budget survey, no exact information on the bias can be obtained. To evaluate the nonresponse bias, bounds for this bias are calculated using linear programming methods for several assumptions. The impact on a price index of a high nonresponse rate among people with a high income can also be assessed by using the elasticity with respect to total expenditure. Attention is also given to the possible nonresponse bias in a time series of price index numbers. The possible nonresponse bias is much larger than the standard error of the estimate.     

Mark-resighting analysis of a California gull population

California gulls ( Larus californicus ) of known age and sex were censused on their breeding colony in 1979, 1980 and 1984 through 1993. Ages of 235 males and 196 females ranged from 4 to 27 years. Age classes used in the analysis were limited to 17, 4 through 19, and 20 or more as a final age category because data on gulls over 20 were sparse. Survival declined with age in a way that was parsimoniously modelled with a quadratic function. Other factors, sex and time, did not explain any variation in survival. Resighting depended on age, sex and time. Younger adults skipped breeding more frequently than did older adults, and females skipped breeding more frequently than did males. There was also good evidence for time dependence in resighting probability, but its inclusion in the model occurred at the expense of interpretability and precision. In a data set such as this, resighting probability may assume more importance than a mere 'nuisance parameter'. In this study, resighting history measured attendance at the breeding ground. In turn, attendance rates may be a manifestation of reproductive strategy, which can also have consequences for survival. In this situation, there may be heterogeneity in both survival and resighting probability that is unexplained by the model. While such complexity may well be a nuisance to deal with, it can also point to important biological questions.     

Comparison of survival estimates obtained from three different methods of recapture in the same population of the great tit

Recaptures of marked birds can be used to estimate their survival. We suspect that the various methods of capture might not be equally suited for this purpose. In the long-term study of the great tit at our institute, recaptures were routinely made in three different ways: capture of parents feeding their young at the nest; capture of birds roosting in nest-boxes; mist-netting. We analyzed 20 years of captures from the study site at the Hoge Veluwe to obtain estimates of adult survival and capture rates from captures obtained by each method. The three sets of estimates differ and our analyses suggest that the best method was the capture of breeding birds, while mist-netting was the least suitable.     

Mark-resighting analysis of a California gull population

California gulls ( Larus californicus ) of known age and sex were censused on their breeding colony in 1979, 1980 and 1984 through 1993. Ages of 235 males and 196 females ranged from 4 to 27 years. Age classes used in the analysis were limited to 17, 4 through 19, and 20 or more as a final age category because data on gulls over 20 were sparse. Survival declined with age in a way that was parsimoniously modelled with a quadratic function. Other factors, sex and time, did not explain any variation in survival. Resighting depended on age, sex and time. Younger adults skipped breeding more frequently than did older adults, and females skipped breeding more frequently than did males. There was also good evidence for time dependence in resighting probability, but its inclusion in the model occurred at the expense of interpretability and precision. In a data set such as this, resighting probability may assume more importance than a mere 'nuisance parameter'. In this study, resighting history measured attendance at the breeding ground. In turn, attendance rates may be a manifestation of reproductive strategy, which can also have consequences for survival. In this situation, there may be heterogeneity in both survival and resighting probability that is unexplained by the model. While such complexity may well be a nuisance to deal with, it can also point to important biological questions.     

Test for age-specificity in survival of the common tern

Much effort in life-history theory has been addressed to the dependence of life-history traits on age, especially the phenomenon of senescence and its evolution. Although senescent declines in survival are well documented in humans and in domestic and laboratory animals, evidence for their occurrence and importance in wild animal species remains limited and equivocal. Several recent papers have suggested that methodological issues may contribute to this problem, and have encouraged investigators to improve sampling designs and to analyse their data using recently developed approaches to modelling of capture-mark-recapture data. Here we report on a three-year, two-site, mark-recapture study of known-aged common terns (Sterna hirundo) in the north-eastern USA. The study was nested within a long-term ecological study in which large numbers of chicks had been banded in each year for > 25 years. We used a range of models to test the hypothesis of an influence of age on survival probability. We also tested for a possible influence of sex on survival. The cross-sectional design of the study (one year's parameter estimates) avoided the possible confounding of effects of age and time. The study was conducted at a time when one of the study sites was being colonized and numbers were increasing rapidly. We detected two-way movements between the sites and estimated movement probabilities in the year for which they could be modelled. We also obtained limited data on emigration from our study area to more distant sites. We found no evidence that survival depended on either sex or age, except that survival was lower among the youngest birds (ages 2-3 years). Despite the large number of birds included in the study (1599 known-aged birds, 2367 total), confidence limits on estimates of survival probability were wide, especially for the oldest age-classes, so that a slight decline in survival late in life could not have been detected. In addition, the cross-sectional design of this study meant that a decline in survival probability within individuals (actuarial senescence) could have been masked by heterogeneity in survival probability among individuals (mortality selection). This emphasizes the need for the development of modelling tools permitting separation of these two phenomena, valid under field conditions in which the recapture probabilities are less than one.     

Estimating avian dispersal distances from data on ringed birds

Data from birds ringed as chicks and recaptured during subsequent breeding seasons provide information on avian natal dispersal distances. However, national patterns of ring reports are influenced by recapture rates as well as by dispersal rates. While an extensive methodology has been developed to study survival rates using models that correct for recapture rates, the same is not true for dispersal. Here, we present such a method, showing how corrections for spatial heterogeneity in recapture rate can be built into estimates of dispersal rates if detailed atlas data and ringing totals can be combined with extensive data on birds ringed as chicks and recaptured as breeding adults. We show how the method can be implemented in the software package SURVIV (White, 1992).     

Investigating impacts of complex sampling on latent growth curve modelling

We investigate the impacts of complex sampling on point and standard error estimates in latent growth curve modelling of survey data. Methodological issues are illustrated with empirical evidence from the analysis of longitudinal data on life satisfaction trajectories using data from the British Household Panel Survey, a national representative survey in Great Britain. A multi-process second-order latent growth curve model with conditional linear growth is used to study variation in the two perceived life satisfaction latent factors considered. The benefits of accounting for the complex survey design are considered, including obtaining unbiased both point and standard error estimates, and therefore correctly specified confidence intervals and statistical tests. We conclude that, even for the rather elaborated longitudinal data models that were considered, estimation procedures are affected by variance-inflating impacts of complex sampling.     

Productivity indices and survival rate estimates from MAPS,a continent-wide programme of constant-effort mist-netting in North America

The Monitoring Avian Productivity and Survivorship (MAPS) programme is a cooperative effort to provide annual regional indices of adult population size and post-fledging productivity and estimates of adult survival rates from data pooled from a network of constant-effort mist-netting stations across North America. This paper provides an overview of the field and analytical methods currently employed by MAPS, a discussion of the assumptions underlying the use of these techniques, and a discussion of the validity of some of these assumptions based on data gathered during the first 5 years (1989-1993) of the programme, during which time it grew from 17 to 227 stations. Ageand species-specific differences in dispersal characteristics are important factors affecting the usefulness of the indices of adult population size and productivity derived from MAPS data. The presence of transients, heterogeneous capture probabilities among stations, and the large sample sizes required by models to deal effectively with these two considerations are important factors affecting the accuracy and precision of survival rate estimates derived from MAPS data. Important results from the first 5 years of MAPS are: (1) indices of adult population size derived from MAPS mist-netting data correlated well with analogous indices derived from point-count data collected at MAPS stations; (2) annual changes in productivity indices generated by MAPS were similar to analogous changes documented by direct nest monitoring and were generally as expected when compared to annual changes in weather during the breeding season; and (3) a model using between-year recaptures in Cormack-Jolly-Seber (CJS) mark-recapture analyses to estimate the proportion of residents among unmarked birds was found, for most tropical-wintering species sampled, to provide a better fit with the available data and more realistic and precise estimates of annual survival rates of resident birds than did standard CJS mark-recapture analyses. A detailed review of the statistical characteristics of MAPS data and a thorough evaluation of the field and analytical methods used in the MAPS programme are currently under way.     

Re-analysis of a banding study to test the effects of an experimental increase in bag limits of mourning doves

In 1966-1971, eastern US states with hunting seasons on mourning doves ( Zenaida macroura ) participated in a study designed to estimate the effects of bag limit increases on population survival rates. More than 400 000 adult and juvenile birds were banded and released during this period, and subsequent harvest and return of bands, together with total harvest estimates from mail and telephone surveys of hunters, provided the database for analysis. The original analysis used an ANOVA framework, and resulted in inferences of no effect of bag limit increase on population parameters (Hayne 1975). We used a logistic regression analysis to infer that the bag limit increase did not cause a biologically significant increase in harvest rate and thus the experiment could not provide any insight into the relationship between harvest and annual survival rates. Harvest rate estimates of breeding populations from geographical subregions were used as covariates in a Program MARK analysis and revealed an association between annual survival and harvest rates, although this relationship is potentially confounded by a latitudinal gradient in survival rates of dove populations. We discuss methodological problems encountered in the analysis of these data, and provide recommendations for future studies of the relationship between harvest and annual survival rates of mourning dove populations.