Inferences about population dynamics from count data using multistate models: a comparison to capture–recapture approaches

Wildlife populations consist of individuals that contribute disproportionately to growth and viability. Understanding a population's spatial and temporal dynamics requires estimates of abundance and demographic rates that account for this heterogeneity. Estimating these quantities can be difficult, requiring years of intensive data collection. Often, this is accomplished through the capture and recapture of individual animals, which is generally only feasible at a limited number of locations. In contrast, N‐mixture models allow for the estimation of abundance, and spatial variation in abundance, from count data alone. We extend recently developed multistate, open population N‐mixture models, which can additionally estimate demographic rates based on an organism's life history characteristics. In our extension, we develop an approach to account for the case where not all individuals can be assigned to a state during sampling. Using only state‐specific count data, we show how our model can be used to estimate local population abundance, as well as density‐dependent recruitment rates and state‐specific survival. We apply our model to a population of black‐throated blue warblers (Setophaga caerulescens) that have been surveyed for 25 years on their breeding grounds at the Hubbard Brook Experimental Forest in New Hampshire, USA. The intensive data collection efforts allow us to compare our estimates to estimates derived from capture–recapture data. Our model performed well in estimating population abundance and density‐dependent rates of annual recruitment/immigration. Estimates of local carrying capacity and per capita recruitment of yearlings were consistent with those published in other studies. However, our model moderately underestimated annual survival probability of yearling and adult females and severely underestimates survival probabilities for both of these male stages. The most accurate and precise estimates will necessarily require some amount of intensive data collection efforts (such as capture–recapture). Integrated population models that combine data from both intensive and extensive sources are likely to be the most efficient approach for estimating demographic rates at large spatial and temporal scales.     

Non‐additive effects of density‐dependent and density‐independent factors on brown trout vital rates

Interactions between density‐dependent and density‐independent processes can lead to variation in both growth and survival rates. Detecting such effects, however, will often require sampling on an individual level and at the appropriate spatial and temporal scale. This study documents substantial variation in survival and growth of stream‐dwelling brown trout Salmo trutta from a small Norwegian stream. The data is based on seasonal capture–recaptures of individually marked trout on fixed stations during eight years. The fish were small‐sized, rarely reaching sizes larger than 20 cm and ages older than seven years. Density varied between 0.2–0.8 fish m^?2. Variation in survival and recapture probabilities was analysed using program MARK. Apparent survival (the probability of being alive and present within the study area) generally decreased with increasing trout density and increasing drought level (measured as lowest observed water flow) during both winter and summer. Further, there was a significant interaction effect between density and water flow, indicating that density‐dependent effects on survival predominated when environmental conditions were benign (no drought), while density‐independent processes were most important under harsh environmental conditions (drought). Observed length‐at‐age during autumn indicated a more or less linear growth trajectory throughout life, and no effect of density, water flow or temperature was found. However, using the individual‐based capture–recapture data to estimated specific growth rate, significant positive effects of water flow and temperature and a negative effect of density were identified. Thus, the capture–recapture data suggest a strong potential for population regulation at the rather low densities found in this stream, and regulation may occur both through effects on survival and growth.     

Estimability of migration survival rates from integrated breeding and winter capture–recapture data

Long‐distance migration is a common phenomenon across the animal kingdom but the scale of annual migratory movements has made it difficult for researchers to estimate survival rates during these periods of the annual cycle. Estimating migration survival is particularly challenging for small‐bodied species that cannot carry satellite tags, a group that includes the vast majority of migratory species. When capture–recapture data are available for linked breeding and non‐breeding populations, estimation of overall migration survival is possible but current methods do not allow separate estimation of spring and autumn survival rates. Recent development of a Bayesian integrated survival model has provided a method to separately estimate the latent spring and autumn survival rates using capture–recapture data, though the accuracy and precision of these estimates has not been formally tested. Here, I used simulated data to explore the estimability of migration survival rates using this model. Under a variety of biologically realistic scenarios, I demonstrate that spring and autumn migration survival can be estimated from the integrated survival model, though estimates are biased toward the overall migration survival probability. The direction and magnitude of this bias are influenced by the relative difference in spring and autumn survival rates as well as the degree of annual variation in these rates. The inclusion of covariates can improve the model's performance, especially when annual variation in migration survival rates is low. Migration survival rates can be estimated from relatively short time series (4–5 years), but bias and precision of estimates are improved when longer time series (10–12 years) are available. The ability to estimate seasonal survival rates of small, migratory organisms opens the door to advancing our understanding of the ecology and conservation of these species. Application of this method will enable researchers to better understand when mortality occurs across the annual cycle and how the migratory periods contribute to population dynamics. Integrating summer and winter capture data requires knowledge of the migratory connectivity of sampled populations and therefore efforts to simultaneously collect both survival and tracking data should be a high priority, especially for species of conservation concern.     

Demographic consequences of changing body size in a terrestrial salamander

Changes in climate can alter individual body size, and the resulting shifts in reproduction and survival are expected to impact population dynamics and viability. However, appropriate methods to account for size‐dependent demographic changes are needed, especially in understudied yet threatened groups such as amphibians. We investigated individual‐ and population‐level demographic effects of changes in body size for a terrestrial salamander using capture–mark–recapture data. For our analysis, we implemented an integral projection model parameterized with capture–recapture likelihood estimates from a Bayesian framework. Our study combines survival and growth data from a single dataset to quantify the influence of size on survival while including different sources of uncertainty around these parameters, demonstrating how selective forces can be studied in populations with limited data and incomplete recaptures. We found a strong dependency of the population growth rate on changes in individual size, mediated by potential changes in selection on mean body size and on maximum body size. Our approach of simultaneous parameter estimation can be extended across taxa to identify eco‐evolutionary mechanisms acting on size‐specific vital rates, and thus shaping population dynamics and viability.     

Open Capture–Recapture Models with Heterogeneity: II. Jolly–Seber Model

Summary Estimation of abundance is important in both open and closed population capture–recapture analysis, but unmodeled heterogeneity of capture probability leads to negative bias in abundance estimates. This article defines and develops a suite of open population capture–recapture models using finite mixtures to model heterogeneity of capture and survival probabilities. Model comparisons and parameter estimation use likelihood‐based methods. A real example is analyzed, and simulations are used to check the main features of the heterogeneous models, especially the quality of estimation of abundance, survival, recruitment, and turnover. The two major advances in this article are the provision of realistic abundance estimates that take account of heterogenetiy of capture, and an appraisal of the amount of overestimation of survival arising from conditioning on the first capture when heterogeneity of survival is present.     

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.     

Movement patterns and study area boundaries: influences on survival estimation in capture–mark–recapture studies

The inability to account for the availability of individuals in the study area during capture–mark–recapture (CMR) studies and the resultant confounding of parameter estimates can make correct interpretation of CMR model parameter estimates difficult. Although important advances based on the Cormack–Jolly–Seber (CJS) model have resulted in estimators of true survival that work by unconfounding either death or recapture probability from availability for capture in the study area, these methods rely on the researcher's ability to select a method that is correctly matched to emigration patterns in the population. If incorrect assumptions regarding site fidelity (non‐movement) are made, it may be difficult or impossible as well as costly to change the study design once the incorrect assumption is discovered. Subtleties in characteristics of movement (e.g. life history‐dependent emigration, nomads vs territory holders) can lead to mixtures in the probability of being available for capture among members of the same population. The result of these mixtures may be only a partial unconfounding of emigration from other CMR model parameters. Biologically‐based differences in individual movement can combine with constraints on study design to further complicate the problem. Because of the intricacies of movement and its interaction with other parameters in CMR models, quantification of and solutions to these problems are needed. Based on our work with stream‐dwelling populations of Atlantic salmon Salmo salar, we used a simulation approach to evaluate existing CMR models under various mixtures of movement probabilities. The Barker joint data model provided unbiased estimates of true survival under all conditions tested. The CJS and robust design models provided similarly unbiased estimates of true survival but only when emigration information could be incorporated directly into individual encounter histories. For the robust design model, Markovian emigration (future availability for capture depends on an individual's current location) was a difficult emigration pattern to detect unless survival and especially recapture probability were high. Additionally, when local movement was high relative to study area boundaries and movement became more diffuse (e.g. a random walk), local movement and permanent emigration were difficult to distinguish and had consequences for correctly interpreting the survival parameter being estimated (apparent survival vs true survival).     

Comparison of photo‐matching algorithms commonly used for photographic capture–recapture studies

Photographic capture–recapture is a valuable tool for obtaining demographic information on wildlife populations due to its noninvasive nature and cost‐effectiveness. Recently, several computer‐aided photo‐matching algorithms have been developed to more efficiently match images of unique individuals in databases with thousands of images. However, the identification accuracy of these algorithms can severely bias estimates of vital rates and population size. Therefore, it is important to understand the performance and limitations of state‐of‐the‐art photo‐matching algorithms prior to implementation in capture–recapture studies involving possibly thousands of images. Here, we compared the performance of four photo‐matching algorithms; Wild‐ID, I3S Pattern+, APHIS, and AmphIdent using multiple amphibian databases of varying image quality. We measured the performance of each algorithm and evaluated the performance in relation to database size and the number of matching images in the database. We found that algorithm performance differed greatly by algorithm and image database, with recognition rates ranging from 100% to 22.6% when limiting the review to the 10 highest ranking images. We found that recognition rate degraded marginally with increased database size and could be improved considerably with a higher number of matching images in the database. In our study, the pixel‐based algorithm of AmphIdent exhibited superior recognition rates compared to the other approaches. We recommend carefully evaluating algorithm performance prior to using it to match a complete database. By choosing a suitable matching algorithm, databases of sizes that are unfeasible to match “by eye” can be easily translated to accurate individual capture histories necessary for robust demographic estimates.     

Cormack–Jolly–Seber model with environmental covariates: A P‐spline approach

In capture–recapture models, survival and capture probabilities can be modelled as functions of time‐varying covariates, such as temperature or rainfall. The Cormack–Jolly–Seber (CJS) model allows for flexible modelling of these covariates; however, the functional relationship may not be linear. We extend the CJS model by semi‐parametrically modelling capture and survival probabilities using a frequentist approach via P‐splines techniques. We investigate the performance of the estimators by conducting simulation studies. We also apply and compare these models with known semi‐parametric Bayesian approaches on simulated and real data sets.     

Parameter Redundancy in Multistate Capture‐Recapture Models

Multistate capture‐recapture models are a powerful tool to address a variety of biological questions concerning dispersal and/or individual variability in wild animal populations. However, biologically meaningful models are often over‐parameterized and consequently some parameters cannot be estimated separately. Identifying which quantities are separately estimable is crucial for proper model selection based upon likelihood tests or information criteria and for the interpretation of the estimates obtained. We show how to investigate parameter redundancy in multistate capture‐recapture models, based on formal methods initially proposed by Catchpole and his associates for exponential family distributions (Catchpole, Freeman and Morgan, 1996. Journal of the Royal Statistical Society Series B 58, 763–774). We apply their approach to three models of increasing complexity.     

Survival costs of within‐ and between‐season mate change in the European blackbird Turdus merula

Many studies of socially monogamous birds discuss the adaptive role of between‐season partner change, but only a handful of them refer to the benefits of pair fidelity in terms of increased survival. Moreover, there are no studies describing the benefits of within‐season mate retention. Our data relating to an urban population of European blackbirds Turdus merula enabled us to test the dependence of survival on pair faithfulness. Because blackbirds divorce within and between seasons, we were able to test the influence of pair faithfulness on their within‐ and between‐season survival and mate fidelity. For this purpose, we used a multievent capture–mark–recapture (MECMR) statistical model, which is based on recapture rates and different pair states (faithful to mate, paired with new partner, or dead). Our study indicated that between‐ and within‐season survival depends on pair states: pair‐bond duration increases survival to the next capture occasion in both sexes. We found that the pair‐bond duration to the current partner increased the chances of being with the same partner during the next breeding occasion, although we failed to find any within‐season pair‐bond influence for females. Our results showed sex differences in mating at the end of the season: females had a much smaller chance of breeding with the current new partner in the next year. This study has demonstrated that within‐ and between‐season survival is dependent on mate retention, and we discuss this in the context of how searching for a new partner could affect the birds’ survival.     

Demographic estimation methods for plants with unobservable life-states

Demographic estimation of vital parameters in plants with an unobservable dormant state is complicated, because time of death is not known. Conventional methods assume that death occurs at a particular time after a plant has last been seen aboveground but the consequences of assuming a particular duration of dormancy have never been tested. Capture–recapture methods do not make assumptions about time of death; however, problems with parameter estimability have not yet been resolved. To date, a critical comparative assessment of these methods is lacking. We analysed data from a 10 year study of Cleistes bifaria, a terrestrial orchid with frequent dormancy, and compared demographic estimates obtained by five varieties of the conventional methods, and two capture–recapture methods. All conventional methods produced spurious unity survival estimates for some years or for some states, and estimates of demographic rates sensitive to the time of death assumption. In contrast, capture–recapture methods are more parsimonious in terms of assumptions, are based on well founded theory and did not produce spurious estimates. In Cleistes, dormant episodes lasted for 1–4 years (mean 1.4, SD 0.74). The capture–recapture models estimated ramet survival rate at 0.86 (SE～0.01), ranging from 0.77–0.94 (SEs≤0.1) in any one year. The average fraction dormant was estimated at 30% (SE 1.5), ranging 16–47% (SEs≤5.1) in any one year. Multistate capture–recapture models showed that survival rates were positively related to precipitation in the current year, but transition rates were more strongly related to precipitation in the previous than in the current year, with more ramets going dormant following dry years. Not all capture–recapture models of interest have estimable parameters; for instance, without excavating plants in years when they do not appear aboveground, it is not possible to obtain independent time‐specific survival estimates for dormant plants. We introduce rigorous computer algebra methods to identify the parameters that are estimable in principle. As life‐states are a prominent feature in plant life cycles, multistate capture–recapture models are a natural framework for analysing population dynamics of plants with dormancy.     

Complete tag loss in capture–recapture studies affects abundance estimates: An elephant seal case study

1. In capture–recapture studies, recycled individuals occur when individuals lose all of their tags and are recaptured as though they were new individuals. Typically, the effect of these recycled individuals is assumed negligible.
2. Through a simulation‐based study of double‐tagging experiments, we examined the effect of recycled individuals on parameter estimates in the Jolly–Seber model with tag loss (Cowen & Schwarz, 2006). We validated the simulation framework using long‐term census data of elephant seals.
3. Including recycled individuals did not affect estimates of capture, survival, and tag‐retention probabilities. However, with low tag‐retention rates, high capture rates, and high survival rates, recycled individuals produced overestimates of population size. For the elephant seal case study, we found population size estimates to be between 8% and 53% larger when recycled individuals were ignored.
4. Ignoring the effects of recycled individuals can cause large biases in population size estimates. These results are particularly noticeable in longer studies.

     

Age‐dependent breeding dispersal and adult survival within a metapopulation of Common Terns Sterna hirundo

Dispersal is increasingly recognized as a process of fundamental importance in population dynamics and other aspects of biology. Concurrently, interest in age‐dependent effects on survival, including actuarial senescence, has increased, especially in studies of long‐lived seabirds. Nevertheless, datasets necessary for studying dispersal and age‐dependent effects are few, as these require simultaneous data collection at two or more sites over many years. We conducted a 22‐year capture‐mark‐recapture study of Common Terns Sterna hirundo at three breeding colonies 10–26 km apart in Buzzards Bay, Massachusetts, USA. All birds in the study were of known age (range 2–28 years, median 7 years, n = 3290) and 77% were of known sex. Estimates of adult recapture, survival and breeding dispersal rates were obtained for all age‐classes from 2 to 20 years. The model that acquired 100% of the QAIC_c (Akaike's Information Criterion adjusted for small sample size and overdispersion) weight in our analysis included age‐specificity in all parameters but no relationship with sex. Our study may be the first to demonstrate age‐specificity in recapture, survival and breeding dispersal rates simultaneously, using a single model. Annual rates of breeding dispersal ranged from <0.01 to 0.27, with a population‐weighted mean of 0.065; they decreased with increasing distance between colony sites and, unexpectedly, increased with age. Breeding dispersal did not increase consistently after years with predation on adults or after an attempt to displace birds from an oiled site. Survival rates did not vary among sites or years. Annual adult survival increased from 0.80 in 2‐year‐old birds to a maximum of approximately 0.88 around age 8 years and then declined to 0.76 at age 20 years, yielding strong evidence for actuarial senescence. The peak annual survival rate of 0.88 is at the low end of other estimates for Common Tern and in the lower part of the range recorded for other terns, but total numbers in the three colonies increased seven‐fold during the study. This was part of a slower increase in the regional population, with net immigration into the study colonies. Our results demonstrate the biological significance of breeding dispersal in local population dynamics and age‐related effects on survival and dispersal from a metapopulation of a long‐lived seabird.     

Updated assessment of hatchery‐reared pallid sturgeon (Forbes & Richardson, 1905) survival in the lower Missouri River

As pallid sturgeon, Scaphirhynchus albus (Forbes & Richardson, 1905), natural reproduction and recruitment remains very minimal in the lower Missouri River from Gavins Point Dam (river kilometer [rkm] 1305.2) to the confluence with the Mississippi River (rkm 0.0), hatchery supplementation and river‐wide monitoring efforts have continued. Annual survival estimates of hatchery‐reared pallid sturgeon stocked in the lower Missouri River were previously estimated during 1994–2008. Low recapture rates prior to 2006 limited the data available to estimate survival, which resulted in considerable uncertainty for the estimate of annual survival of age‐1 fish. Therefore, the objective was to provide more precise estimates of annual survival of pallid sturgeon using five additional years of stocking and sampling. The Cormack‐Jolly‐Seber model structure provided in program MARK was used to estimate the age‐specific survival estimates. Over 135 000 hatchery‐reared pallid sturgeon were released during 1994–2011 and recaptured at a rate of 1.9%, whereby estimates for the annual survival of age‐0 (Ø = 0.048) and >age‐1 (Ø = 0.931) were similar to those previously reported, but the age‐1 (Ø = 0.403) survival estimate was 52% lower. Post hoc analysis using time‐specific survival estimates indicated lower survival for age‐1 fish post‐2003 year classes, relative to the pre‐2002 year classes. An analysis confirms that hatchery‐reared pallid sturgeon continue to survive in the wild. However, low survival during the first 2 years of life is a management concern as efforts are aimed at maximizing genetic diversity and population growth. A follow‐up analysis also demonstrated the variability of capture rates and survival over time, which reinforces the need to continue to monitor and evaluate mark‐recapture data. The mark‐recapture efforts have provided demographic parameter estimates that remain a critical component for species recovery as these data are incorporated into population models.     

Who are we sampling? Apparent survival differs between methods in a secretive species

Survival is a fundamental parameter in population dynamics with increasing importance in the management and conservation strategies of wildlife populations. Survival probability in vertebrates is usually estimated by live‐encounter data obtained by means of physical mark–capture–recapture protocols. Non‐invasive acoustic marking relying on individual‐specific features of signals has been alternatively applied as a marking technique, especially in secretive species. Nevertheless, to date no research has compared survival rate estimates obtained by acoustic and physical marking. We estimated half‐yearly and annual survival and recapture rates of a secretive and threatened passerine, the Dupont's lark Chersophilus duponti, using two separate live‐encounter data sets of males collected simultaneously by physical and acoustic marking in the same study area. The separate analysis of both methods led to different model structures, since transient individuals had to be accounted for in the acoustic marking but not in the physical marking data set. Furthermore, while reencounter probabilities did not differ between methods, survival estimates employing physical marking were lower than those obtained acoustically, especially between the postbreeding and the breeding period when the apparent survival of colour‐banded birds was twice as low as for acoustic marking. The combination of marking methods suggested the existence of different subsets of individuals differentially sampled within the population: whereas colour‐banded males seemed to represent the territorial fraction of the population, both resident and floater individuals were probably detected by acoustic marking. Using traditional mark–recapture methods exclusively could have misled our estimates of survival rates, potentially affecting prospective predictions of population dynamics. Acoustic marking has been poorly applied in mark–recapture studies, but might be a powerful complement to obtain accurate estimates of fundamental demographic parameters such as survival and dispersal.     

An integrated population model for a long‐lived ungulate: more efficient data use with Bayesian methods

We develop an integrated population model for Svalbard reindeer Rangifer tarandus platyrhynchus, and demonstrate how this type of model can be used to extract more information from the data and separate different sources of variability in population estimates. Our model combines individual mark–recapture data with population counts and harvesting data within a Bayesian model framework, and accounts for observation error, environmental and demographic stochasticity, and age structure. From this model we obtain annual estimates of age‐specific population size, survival and fecundity. The model provides estimates of age structure at a finer scale than that found in the census data, and enables us to estimate survival for the period before calves are first caught and marked, i.e. before they enter the individual mark–recapture data. The modeling framework provides an improved approach to studying age‐structured populations that are imperfectly censused and where the demography of only a sample of individuals is known. We use data from independent censuses of the same population to evaluate population estimates obtained from the model, and show that it is successful at correcting for different types of observation error. Based on our model results, we suggest that allocating resources to the collection of supplementary mark–recapture data could improve the reliability of population projections more than making regular population censuses as exhaustive as possible. Our work demonstrates how integrated Bayesian population modeling can be used to increase the amount of information extracted from collections of data, identifying and disentangling sources of variation in individual performance and population size. This represents an important step towards increasing the predictive ability of population growth models for long‐lived species experiencing changes in environmental conditions and harvesting regimes.     

Interactive effects of senescence and natural disturbance on the annual survival probabilities of snail kites

Individuals in wild populations face risks associated with both intrinsic (i.e. aging) and external (i.e. environmental) sources of mortality. Condition‐dependent mortality occurs when there is an interaction between such factors; however, few studies have clearly demonstrated condition‐dependent mortality and some have even argued that condition‐dependent mortality does not occur in wild avian populations. Using large sample sizes (2084 individuals, 3746 re‐sights) of individual‐based longitudinal data collected over a 33 year period (1976–2008) on multiple cohorts, we used a capture–mark–recapture framework to model age‐dependent survival in the snail kite Rostrhamus sociabilis plumbeus population in Florida. Adding to the growing amount of evidence for actuarial senescence in wild populations, we found evidence of senescent declines in survival probabilities in adult kites. We also tested the hypothesis that older kites experienced condition‐dependent mortality during a range‐wide drought event (2000–2002). The results provide convincing evidence that the annual survival probability of senescent kites was disproportionately affected by the drought relative to the survival probability of prime‐aged adults. To our knowledge, this is the first evidence of condition‐dependent mortality to be demonstrated in a wild avian population, a finding which challenges recent conclusions drawn in the literature. Our study suggests that senescence and condition‐dependent mortality can affect the demography of wild avian populations. Accounting for these sources of variation may be particularly important to appropriately compute estimates of population growth rate, and probabilities of quasi‐extinctions.     

Survival estimates of bird species across altered habitats in the tropical Andes

The probability of long‐term persistence of a population is strongly determined by adult survival rates, but estimates of survival are currently lacking for most species of birds in the tropical Andes, a global biodiversity hotspot. We calculated apparent survival rates of birds in the Ecuadorian tropical Andes using a moderately long‐term (11 yr) capture–recapture dataset from three habitats that varied in how much they had been modified by human activities (native forest, introduced forest, and shrubs). We fit mark–recapture models for 28 species with habitat as a covariable. For all species, recapture rates between sampling sessions were low and varied from 0.04 for Rainbow Starfrontlets (Coeligena iris) to 0.41 for Stripe‐headed Brushfinches (Arremon assimilis) when averaged across all occupied habitats. Annual survival rates varied from 0.07 for Black‐crested Warblers (Margarornis squamiger) to 0.75 for Violet‐throated Metaltails (Metallura baroni). We found no significant differences in survival rates either among habitats or species grouped by habitat specialization. Because we found similar survival rates in native forest and human‐modified habitats, our results support those of recent studies concerning the potential value of secondary habitats for the conservation of some species of birds in the tropics. However, our conclusions are tempered by the uncertainty around the estimates of survival rates. Despite the relatively long‐term nature of our study, obtaining survival estimates for bird species in this region was challenging, and either more years of study or modification of field protocols may be needed to obtain more precise survival estimates.     

Using dynamic N‐mixture models to test cavity limitation on northern flying squirrel demographic parameters using experimental nest box supplementation

Dynamic N‐mixture models have been recently developed to estimate demographic parameters of unmarked individuals while accounting for imperfect detection. We propose an application of the Dail and Madsen ( 2011 : Biometrics, 67 , 577–587) dynamic N‐mixture model in a manipulative experiment using a before‐after control‐impact design (BACI). Specifically, we tested the hypothesis of cavity limitation of a cavity specialist species, the northern flying squirrel, using nest box supplementation on half of 56 trapping sites. Our main purpose was to evaluate the impact of an increase in cavity availability on flying squirrel population dynamics in deciduous stands in northwestern Québec with the dynamic N‐mixture model. We compared abundance estimates from this recent approach with those from classic capture–mark–recapture models and generalized linear models. We compared apparent survival estimates with those from Cormack–Jolly–Seber (CJS) models. Average recruitment rate was 6 individuals per site after 4 years. Nevertheless, we found no effect of cavity supplementation on apparent survival and recruitment rates of flying squirrels. Contrary to our expectations, initial abundance was not affected by conifer basal area (food availability) and was negatively affected by snag basal area (cavity availability). Northern flying squirrel population dynamics are not influenced by cavity availability at our deciduous sites. Consequently, we suggest that this species should not be considered an indicator of old forest attributes in our study area, especially in view of apparent wide population fluctuations across years. Abundance estimates from N‐mixture models were similar to those from capture–mark–recapture models, although the latter had greater precision. Generalized linear mixed models produced lower abundance estimates, but revealed the same relationship between abundance and snag basal area. Apparent survival estimates from N‐mixture models were higher and less precise than those from CJS models. However, N‐mixture models can be particularly useful to evaluate management effects on animal populations, especially for species that are difficult to detect in situations where individuals cannot be uniquely identified. They also allow investigating the effects of covariates at the site level, when low recapture rates would require restricting classic CMR analyses to a subset of sites with the most captures.