Understanding implications of detection heterogeneity in wildlife abundance estimation

Negative bias in mark-recapture abundance estimators due to heterogeneity in detection (capture) probability is a well-known problem, but we believe most biologists do not understand why heterogeneity causes bias and how bias can be reduced. We demonstrate how heterogeneity creates dependence and bias in mark-recapture approaches to abundance estimation. In comparison, heterogeneity, and hence estimator bias, is not as problematic for distance sampling and mark-resight methods because both techniques estimate detection probabilities based on a known quantity. We show how the introduction of a known number of individuals planted into a study population prior to a mark-recapture survey can reduce bias from heterogeneity in detection probability. We provide examples with simulation and an analysis of motion-sensitive camera data from a study population of introduced eastern wild turkeys (Meleagris gallopavo silvestris) of known size with a subset of telemetered birds. In choosing a method for abundance estimation, careful consideration should be given to assumptions and how heterogeneity in detection probability can be accommodated for each application.     

Pooling robustness in distance sampling: Avoiding bias when there is unmodelled heterogeneity

The pooling robustness property of distance sampling results in unbiased abundance estimation even when sources of variation in detection probability are not modeled. However, this property cannot be relied upon to produce unbiased subpopulation abundance estimates when using a single pooled detection function that ignores subpopulations. We investigate by simulation the effect of differences in subpopulation detectability upon bias in subpopulation abundance estimates. We contrast subpopulation abundance estimates using a pooled detection function with estimates derived using a detection function model employing a subpopulation covariate. Using point transect survey data from a multispecies songbird study, species-specific abundance estimates are compared using pooled detection functions with and without a small number of adjustment terms, and a detection function with species as a covariate. With simulation, we demonstrate the bias of subpopulation abundance estimates when a pooled detection function is employed. The magnitude of the bias is positively related to the magnitude of disparity between the subpopulation detection functions. However, the abundance estimate for the entire population remains unbiased except when there is extreme heterogeneity in detection functions. Inclusion of a detection function model with a subpopulation covariate essentially removes the bias of the subpopulation abundance estimates. The analysis of the songbird point count surveys shows some bias in species-specific abundance estimates when a pooled detection function is used. Pooling robustness is a unique property of distance sampling, producing unbiased abundance estimates at the level of the study area even in the presence of large differences in detectability between subpopulations. In situations where subpopulation abundance estimates are required for data-poor subpopulations and where the subpopulations can be identified, we recommend the use of subpopulation as a covariate to reduce bias induced in subpopulation abundance estimates.     

Monitoring of the Flat-Tailed Horned Lizard With Methods Incorporating Detection Probability

Abstract: Difficulty in monitoring the flat-tailed horned lizard (Phrynosoma mcallii) has led to controversy over its conservation status. The difficulty in detecting this species has discouraged large-scale estimates of abundance and led to uncertainty over whether the species exists in population sizes of sufficient size for long-term persistence. We incorporated detection probability into monitoring of this species using closed mark—recapture and distance-sampling methods. Density estimation from mark—recapture abundance estimates was improved using an estimate of the proportion of time lizards were on the plot. We estimated the probability of detection on the line for distance sampling and adjusted density estimates accordingly. We estimated the populations of the Yuha Basin Management Area in 2002 and the East Mesa Management Area, Imperial County, California, USA, in 2003 to be 25,514 (95% CI 14,444-38,970) and 42,619 (95% CI 23,161-67,639), respectively. Two estimates of detection probability on the line in distance sampling by different methods were 0.45 and 0.65. Density estimates derived from distance analyses for 3 East Mesa Management Area plots and the Yuha Basin Management Area were 1.55 per ha (95% CI 0.64-3.76) and 0.41 per ha (95% CI 0.22-0.7), respectively. These are the first large-scale estimates of abundance and density for P. mcallii.     

A Model-Based Approach for Making Ecological Inference from Distance Sampling Data

Summary .  We consider a fully model-based approach for the analysis of distance sampling data. Distance sampling has been widely used to estimate abundance (or density) of animals or plants in a spatially explicit study area. There is, however, no readily available method of making statistical inference on the relationships between abundance and environmental covariates. Spatial Poisson process likelihoods can be used to simultaneously estimate detection and intensity parameters by modeling distance sampling data as a thinned spatial point process. A model-based spatial approach to distance sampling data has three main benefits: it allows complex and opportunistic transect designs to be employed, it allows estimation of abundance in small subregions, and it provides a framework to assess the effects of habitat or experimental manipulation on density. We demonstrate the model-based methodology with a small simulation study and analysis of the Dubbo weed data set. In addition, a simple ad hoc method for handling overdispersion is also proposed. The simulation study showed that the model-based approach compared favorably to conventional distance sampling methods for abundance estimation. In addition, the overdispersion correction performed adequately when the number of transects was high. Analysis of the Dubbo data set indicated a transect effect on abundance via Akaike's information criterion model selection. Further goodness-of-fit analysis, however, indicated some potential confounding of intensity with the detection function.     

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.     

Estimating abundance and population trends when detection is low and highly variable: A comparison of three methods for the Hermann's tortoise

Assessing population trends is a basic prerequisite to carrying out adequate conservation strategies. Selecting an appropriate method to monitor animal populations can be challenging, particularly for low-detection species such as reptiles. This study compares 3 detection-corrected abundance methods (capture–recapture, distance sampling, and N-mixture) used to assess population size of the threatened Hermann's tortoise. We used a single dataset of 432 adult tortoise observations collected at 118 sampling sites in the Plaine des Maures, southeastern France. We also used a dataset of 520 tortoise observations based on radiotelemetry data collected from 10 adult females to estimate and model the availability (g0) needed for distance sampling. We evaluated bias for N-mixture and capture–recapture, by using simulations based on different values of detection probabilities. Finally, we conducted a power analysis to estimate the ability of the 3 methods to detect changes in Hermann's tortoise abundances. The abundance estimations we obtained using distance sampling and N-mixture models were respectively 1.75 and 2.19 times less than those obtained using the capture–recapture method. Our results indicated that g0 was influenced by temperature variations and can differ for the same temperature on different days. Simulations showed that the N-mixture models provide unstable estimations for species with detection probabilities <0.5, whereas capture–recapture estimations were unbiased. Power analysis showed that none of the 3 methods were precise enough to detect slow population changes. We recommend that great care should be taken when implementing monitoring designs for species with large variation in activity rates and low detection probabilities. Although N-mixture models are easy to implement, we would not recommend using them in situations where the detection probability is very low at the risk of providing biased estimates. Among the 3 methods allowing estimation of tortoise abundances, capture–recapture should be preferred to assess population trends. © 2013 The Wildlife Society.     

Estimating abundance of mountain lions from unstructured spatial sampling

Mountain lions (Puma concolor) are often difficult to monitor because of their low capture probabilities, extensive movements, and large territories. Methods for estimating the abundance of this species are needed to assess population status, determine harvest levels, evaluate the impacts of management actions on populations, and derive conservation and management strategies. Traditional mark–recapture methods do not explicitly account for differences in individual capture probabilities due to the spatial distribution of individuals in relation to survey effort (or trap locations). However, recent advances in the analysis of capture–recapture data have produced methods estimating abundance and density of animals from spatially explicit capture–recapture data that account for heterogeneity in capture probabilities due to the spatial organization of individuals and traps. We adapt recently developed spatial capture–recapture models to estimate density and abundance of mountain lions in western Montana. Volunteers and state agency personnel collected mountain lion DNA samples in portions of the Blackfoot drainage (7,908 km²) in west-central Montana using 2 methods: snow back-tracking mountain lion tracks to collect hair samples and biopsy darting treed mountain lions to obtain tissue samples. Overall, we recorded 72 individual capture events, including captures both with and without tissue sample collection and hair samples resulting in the identification of 50 individual mountain lions (30 females, 19 males, and 1 unknown sex individual). We estimated lion densities from 8 models containing effects of distance, sex, and survey effort on detection probability. Our population density estimates ranged from a minimum of 3.7 mountain lions/100 km² (95% CI 2.3–5.7) under the distance only model (including only an effect of distance on detection probability) to 6.7 (95% CI 3.1–11.0) under the full model (including effects of distance, sex, survey effort, and distance × sex on detection probability). These numbers translate to a total estimate of 293 mountain lions (95% CI 182–451) to 529 (95% CI 245–870) within the Blackfoot drainage. Results from the distance model are similar to previous estimates of 3.6 mountain lions/100 km² for the study area; however, results from all other models indicated greater numbers of mountain lions. Our results indicate that unstructured spatial sampling combined with spatial capture–recapture analysis can be an effective method for estimating large carnivore densities. Published 2012. This article is a U.S. Government work and is in the public domain in the USA.     

Estimating abundance using mark-resight when sampling is with replacement or the number of marked individuals is unknown

Summary .  Although mark–resight methods can often be a less expensive and less invasive means for estimating abundance in long-term population monitoring programs, two major limitations of the estimators are that they typically require sampling without replacement and/or the number of marked individuals available for resighting to be known exactly. These requirements can often be difficult to achieve. Here we address these limitations by introducing the Poisson log and zero-truncated Poisson log-normal mixed effects models (PNE and ZPNE, respectively). The generalized framework of the models allow the efficient use of covariates in modeling resighting rate and individual heterogeneity parameters, information-theoretic model selection and multimodel inference, and the incorporation of individually unidentified marks. Both models may be implemented using standard statistical computing software, but they have also been added to the mark–recapture freeware package Program MARK . We demonstrate the use and advantages of (Z)PNE using black-tailed prairie dog data recently collected in Colorado. We also investigate the expected relative performance of the models in simulation experiments. Compared to other available estimators, we generally found (Z)PNE to be more precise with little or no loss in confidence interval coverage. With the recent introduction of the logit-normal mixed effects model and (Z)PNE, a more flexible and efficient framework for mark–resight abundance estimation is now available for the sampling conditions most commonly encountered in these studies.     

Estimating Abundance From One-Dimensional Passive Acoustic Surveys

ABSTRACT Conventional distance sampling, the most-used method of estimating animal density and abundance, requires ranges to detected individuals, which are not easily measured for vocalizations. However, in some circumstances the sequential pattern of detection of vocalizations along a transect line gives information about the range of detection. Thus, from a one-dimensional acoustic point-transect survey (i.e., records of vocalizations detected or not detected at regularly spaced listening stations) it is possible to obtain a useful estimate of density or abundance. I developed equations for estimation of density for one-dimensional surveys. Using simulations I found that for the method to have little bias when both range of detection and rate of vocalization need to be estimated, stations needed to be spaced at 30–80% of the range of detection and the rate of vocalization should be >0.7. If either the range of detection or rate of vocalization is known, conditions are relaxed, and when both parameters are known the method works well almost universally. In favorable conditions for one-dimensional methods, estimated abundances have overall errors not much larger than those from conventional line-transect distance sampling. The methods appeared useful when applied to real acoustic data from whale surveys. The techniques may also be useful in surveys with nonacoustic detection of animals.     

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.     

Survival, movement, and resource use of the butterfly Parnassius clodius

Abstract.  1. A mark–recapture study was conducted on the American Apollo butterfly Parnassius clodius Menetries during three field seasons (1998–2000) to examine its movement patterns over the course of a season within a sagebrush meadow in Grand Teton National Park, Wyoming, U.S.A. The study examined how resources affected butterfly distribution patterns and used mark–recapture data to gain insight into movement differences between sexes and over time.
2. The average straight-line movement of P. clodius was 202 m day⁻¹, adjusted for sampling effort at different distances. Movement estimates in all 3 years were highly correlated with the average distance between plots sampled.
3. Butterfly abundance was correlated positively with per cent cover of its host plant Dicentra uniflora , but this relationship decreased in importance during the peak of the flight period when individuals may be more interested in finding mates. There was a weak, positive correlation between butterfly abundance and the abundance of its primary nectar source, Eriogonum umbellatum in 1999, but no relationship in 2000.
4. Survival, recapture, and transition probabilities were estimated using open population, capture–recapture models. Survival and recapture probability decreased over the course of each season, while the probability of moving between plots increased. Recapture probability was significantly lower for females than for males among all 3 years, but there was no difference between the sexes in survival rate.     

Accounting for imperfect detection of groups and individuals when estimating abundance

If animals are independently detected during surveys, many methods exist for estimating animal abundance despite detection probabilities <1. Common estimators include double‐observer models, distance sampling models and combined double‐observer and distance sampling models (known as mark‐recapture‐distance‐sampling models; MRDS). When animals reside in groups, however, the assumption of independent detection is violated. In this case, the standard approach is to account for imperfect detection of groups, while assuming that individuals within groups are detected perfectly. However, this assumption is often unsupported. We introduce an abundance estimator for grouped animals when detection of groups is imperfect and group size may be under‐counted, but not over‐counted. The estimator combines an MRDS model with an N‐mixture model to account for imperfect detection of individuals. The new MRDS‐Nmix model requires the same data as an MRDS model (independent detection histories, an estimate of distance to transect, and an estimate of group size), plus a second estimate of group size provided by the second observer. We extend the model to situations in which detection of individuals within groups declines with distance. We simulated 12 data sets and used Bayesian methods to compare the performance of the new MRDS‐Nmix model to an MRDS model. Abundance estimates generated by the MRDS‐Nmix model exhibited minimal bias and nominal coverage levels. In contrast, MRDS abundance estimates were biased low and exhibited poor coverage. Many species of conservation interest reside in groups and could benefit from an estimator that better accounts for imperfect detection. Furthermore, the ability to relax the assumption of perfect detection of individuals within detected groups may allow surveyors to re‐allocate resources toward detection of new groups instead of extensive surveys of known groups. We believe the proposed estimator is feasible because the only additional field data required are a second estimate of group size.     

Abundance estimation of long-diving animals using line transect methods

Line transect sampling is one of the most widely used methods for estimating the size of wild animal populations. An assumption in standard line transect sampling is that all the animals on the trackline are detected without fail. This assumption tends to be violated for marine mammals with surfacing/diving behaviors. The detection probability on the trackline is estimated using duplicate sightings from double-platform line transect methods. The double-platform methods, however, are insufficient to estimate the abundance of long-diving animals because these animals can be completely missed while the observers pass. We developed a more flexible hazard probability model that incorporates information on surfacing/diving patterns obtained from telemetry data. The model is based on a stochastic point process and is statistically tractable. A simulation study showed that the new model provides near-unbiased abundance estimates, whereas the traditional hazard rate and hazard probability models produce considerably biased estimates. As an illustration, we applied the model to data on the Baird's beaked whale (Berardius bairdii) in the western North Pacific.     

Improved methods for estimating abundance and related demographic parameters from mark‐resight data

Over the past decade, there has been much methodological development for the estimation of abundance and related demographic parameters using mark‐resight data. Often viewed as a less‐invasive and less‐expensive alternative to conventional mark recapture, mark‐resight methods jointly model marked individual encounters and counts of unmarked individuals, and recent extensions accommodate common challenges associated with imperfect detection. When these challenges include both individual detection heterogeneity and an unknown marked sample size, we demonstrate several deficiencies associated with the most widely used mark‐resight models currently implemented in the popular capture‐recapture freeware Program MARK. We propose a composite likelihood solution based on a zero‐inflated Poisson log‐normal model and find the performance of this new estimator to be superior in terms of bias and confidence interval coverage. Under Pollock's robust design, we also extend the models to accommodate individual‐level random effects across sampling occasions as a potentially more realistic alternative to models that assume independence. As a motivating example, we revisit a previous analysis of mark‐resight data for the New Zealand Robin (Petroica australis) and compare inferences from the proposed estimators. For the all‐too‐common situation where encounter rates are low, individual detection heterogeneity is non‐negligible, and the number of marked individuals is unknown, we recommend practitioners use the zero‐inflated Poisson log‐normal mark‐resight estimator as now implemented in Program MARK.     

Predicting and correcting bias caused by measurement error in line transect sampling using multiplicative error models

Line transect sampling is one of the most widely used methods for animal abundance assessment. Standard estimation methods assume certain detection on the transect, no animal movement, and no measurement errors. Failure of the assumptions can cause substantial bias. In this work, the effect of error measurement on line transect estimators is investigated. Based on considerations of the process generating the errors, a multiplicative error model is presented and a simple way of correcting estimates based on knowledge of the error distribution is proposed. Using beta models for the error distribution, the effect of errors and of the proposed correction is assessed by simulation. Adequate confidence intervals for the corrected estimates are obtained using a bootstrap variance estimate for the correction and the delta method. As noted by Chen (1998, Biometrics 54, 899-908), even unbiased estimators of the distances might lead to biased density estimators, depending on the actual error distribution. In contrast with the findings of Chen, who used an additive model, unbiased estimation of distances, given a multiplicative model, lead to overestimation of density. Some error distributions result in observed distance distributions that make efficient estimation impossible, by removing the shoulder present in the original detection function. This indicates the need to improve field methods to reduce measurement error. An application of the new methods to a real data set is presented.     

哺乳类动物数量调查中的截线抽样法与逆向截线法

文本介绍了可用于哺乳动物数量调查的一种新方法--截线抽样法,并在此方法的基础上,结合我国动物调查实践,提出了逆向截线法。逆向截线法在贺兰山自然保护区马鹿数量调查中,取得了较好效果。     

Incorporating covariates into standard line transect analyses

An implicit assumption of standard line transect methodology is that detection probabilities depend solely on the perpendicular distance of detected objects to the transect line. Heterogeneity in detection probabilities is commonly minimized using stratification, but this may be precluded by small sample sizes. We develop a general methodology which allows the effects of multiple covariates to be directly incorporated into the estimation procedure using a conditional likelihood approach. Small sample size properties of estimators are examined via simulations. As an example the method is applied to eastern tropical Pacific dolphin sightings data.     

The effect of attractant lures in camera trapping: a case study of population estimates for the Iberian lynx (Lynx pardinus)

Capture–recapture analysis of camera trap data is a conventional method to estimate the abundance of free-ranging wild felids. Due to notorious low detection rates of felids, it is important to increase the detection probability during sampling. In this study, we report the effectiveness of attractants as a tool for improving the efficiency of camera trap sampling in abundance estimation of Iberian lynx. We developed a grid system of camera stations in which stations with and without attractant lures were spatially alternated across known Iberian lynx habitat. Of the ten individuals identified, five were detected at stations with no attractant (blind sets), and nine, at the lured stations. Thirty-eight percent of blind set station’s independent captures and 10?% of lured station’s independent captures resulted in photographs unsuitable for correct individual identification. The total capture probability at lured stations was higher than that obtained at blind set stations. The estimates obtained with blind set cameras underestimated the number of lynxes compared to lured cameras. In our study, it appears that the use of lures increased the efficiency of trail camera captures and, therefore, the accuracy of capture–recapture analysis. The observed failure to detect known individuals at blind set camera stations may violate capture–recapture assumptions and bias abundance estimates.     

Adaptive line transect sampling

Adaptive line transect sampling offers the potential of improved population density estimation efficiency over conventional line transect sampling when populations are spatially clustered. In adaptive sampling, survey effort is increased when areas of high animal density are located, thereby increasing the number of observations. Its disadvantage is that the survey effort required is not known in advance. We develop an adaptive line transect methodology that, by varying the degree of adaptation, allows total effort to be fixed at the design stage. Relative to conventional line transect surveys, it also provides better survey coverage in the event of disruption in survey effort, e.g., due to poor weather. In analysis, sightings from the adaptive sections are downweighted in proportion to the increase in effort. We evaluate the methodology by simulation and report on surveys of harbor porpoise in the Gulf of Maine, in which the approach was compared with conventional line transect sampling.