Using soundscape recordings to estimate bird species abundance, richness, and composition

ABSTRACT Point counts are the most frequently used technique for sampling bird populations and communities, but have well‐known limitations such as inter‐ and intraobserver errors and limited availability of expert field observers. The use of acoustic recordings to survey birds offers solutions to these limitations. We designed a Soundscape Recording System (SRS) that combines a four‐channel, discrete microphone system with a quadraphonic playback system for surveying bird communities. We compared the effectiveness of SRS and point counts for estimating species abundance, richness, and composition of riparian breeding birds in California by comparing data collected simultaneously using both methods. We used the temporal‐removal method to estimate individual bird detection probabilities and species abundances using the program MARK. Akaike's Information Criterion provided strong evidence that detection probabilities differed between the two survey methods and among the 10 most common species. The probability of detecting birds was higher when listening to SRS recordings in the laboratory than during the field survey. Additionally, SRS data demonstrated a better fit to the temporal‐removal model assumptions and yielded more reliable estimates of detection probability and abundance than point‐count data. Our results demonstrate how the perceptual constraints of observers can affect temporal detection patterns during point counts and thus influence abundance estimates derived from time‐of‐detection approaches. We used a closed‐population capture–recapture approach to calculate jackknife estimates of species richness and average species detection probabilities for SRS and point counts using the program CAPTURE. SRS and point counts had similar species richness and detection probabilities. However, the methods differed in the composition of species detected based on Jaccard's similarity index. Most individuals (83%) detected during point counts vocalized at least once during the survey period and were available for detection using a purely acoustic technique, such as SRS. SRS provides an effective method for surveying bird communities, particularly when most species are detected by sound. SRS can eliminate or minimize observer biases, produce permanent records of surveys, and resolve problems associated with the limited availability of expert field observers.     

Comparison of methods for estimating density of forest songbirds from point counts

New analytical methods have been promoted for estimating the probability of detection and density of birds from count data but few studies have compared these methods using real data. We compared estimates of detection probability and density from distance and time-removal models and survey protocols based on 5- or 10-min counts and outer radii of 50 or 100 m. We surveyed singing male Acadian flycatchers (Empidonax virescens), cerulean warblers (Dendroica cerulea), Kentucky warblers (Oporornis formosus), Louisiana waterthrushes (Parkesia motacilla), wood thrushes (Hylocichla mustelina), and worm-eating warblers (Helmitheros vermivorum) in bottomland and upland forest across 5 states in the Central Hardwoods Bird Conservation Region during the breeding season in 2007 and 2008. Detection probabilities differed between distance and time-removal models and species detectabilities were affected differently by year, forest type, and state. Density estimates from distance models were generally higher than from time-removal models, resulting from lower detection probabilities estimated by distance models. We found support for individual heterogeneity (modeled as a finite mixture model) in the time-removal models and that 50-m radius counts generated density estimates approximately twice as high as 100-m radius counts. Users should be aware that in addition to estimating different components of detectability, density estimates derived from distance and time-removal models can be affected by survey protocol because some count durations and plot radii may better meet model assumptions than others. The choice of a method may not affect the use of estimates for relative comparisons (e.g., when comparing habitats) but could affect conclusions when used to estimate population size. We recommend careful consideration of assumptions when deciding on point-count protocol and selection of analysis methods. © 2011 The Wildlife Society.     

Evaluation of three methods to estimate density and detectability from roadside point counts

Roadside point counts are often used to estimate trends of bird populations. The use of aural counts of birds without adjustment for detection probability, however, can lead to incorrect population trend estimates. We compared precision of estimates of density and detectability of whistling northern bobwhites (Colinus virginianus) using distance sampling, independent double-observer, and removal methods from roadside surveys. Two observers independently recorded each whistling bird heard, distance from the observer, and time of first detection at 362 call-count stops in Ohio. We examined models that included covariates for year and observer effects for each method and distance from observer effects for the double-observer and removal methods using Akaike's Information Criterion (AIC). The best model of detectability from distance sampling included observer and year effects. The best models from the removal and double-observer techniques included observer and distance effects. All 3 methods provided precise estimates of detection probability (CV = 2.4–4.4%) with a range of detectability of 0.44–0.95 for a 6-min survey. Density estimates from double-observer surveys had the lowest coefficient of variation (2005 = 3.2%, 2006 = 1.7%), but the removal method also provided precise estimates of density (2005 CV = 3.4%, 2006 CV = 4.8%), and density estimates from distance sampling were less precise (2005 CV = 9.6%, 2006 CV = 7.9%). Assumptions of distance sampling were violated in our study because probability of detecting bobwhites near the observer was <1 or the roadside survey points were not randomly distributed with respect to the birds. Distances also were not consistently recorded by individual members of observer pairs. Although double-observer surveys provided more precise estimates, we recommend using the removal method to estimate detectability and abundance of bobwhites. The removal method provided precise estimates of density and detection probability and requires half the personnel time as double-observer surveys. Furthermore, the likelihood of meeting model assumptions is higher for the removal survey than with independent double-observers. © 2011 The Wildlife Society.     

Modelling occurrence and abundance of species when detection is imperfect

Relationships between species abundance and occupancy are of considerable interest in metapopulation biology and in macroecology. Such relationships may be described concisely using probability models that characterize variation in abundance of a species. However, estimation of the parameters of these models in most ecological problems is impaired by imperfect detection. When organisms are detected imperfectly, observed counts are biased estimates of true abundance, and this induces bias in stated occupancy or occurrence probability. In this paper we consider a class of models that enable estimation of abundance/occupancy relationships from counts of organisms that result from surveys in which detection is imperfect. Under such models, parameter estimation and inference are based on conventional likelihood methods. We provide an application of these models to geographically extensive breeding bird survey data in which alternative models of abundance are considered that include factors that influence variation in abundance and detectability. Using these models, we produce estimates of abundance and occupancy maps that honor important sources of spatial variation in avian abundance and provide clearly interpretable characterizations of abundance and occupancy adjusted for imperfect detection.     

Optimizing surveys for marsh songbirds: does timing matter?

Analysis of data from point counts, a common method for monitoring bird population trends, has evolved to produce estimates of various population parameters (e.g., density, abundance, and occupancy) while simultaneously estimating detection probability. An important consideration when designing studies using point counts is to maximize detection probability while minimizing variation in detection probability both within and between counts. Our objectives were to estimate detection probabilities for three marsh songbirds, including Marsh Wrens (Cistothorus palustris), Swamp Sparrows (Melospiza georgiana), and Yellow‐headed Blackbirds (Xanthocephalus xanthocephalus), as a function of weather covariates and to evaluate temporal variability in detection probability of these three species. We conducted paired, unlimited radius, 10‐min point counts during consecutive morning and evening survey periods for our three focal species at 56 wetlands in Iowa from 20 April to 10 July 2010. Mean detection probabilities ranged from 0.272 (SE = 0.042) for Marsh Wrens to 0.365 (SE = 0.052) for Swamp Sparrows. Time of season was positively correlated with detection probability for Swamp Sparrows, but was negatively correlated with detection probability for Yellow‐headed Blackbirds, suggesting that detection probability increased during the breeding season for Swamp Sparrows and was highest early in the breeding season for Yellow‐headed Blackbirds. Understanding how detection probabilities of marsh songbirds vary throughout the breeding season allows targeted survey efforts that maximize detection probabilities for these species. Furthermore, consistent detection probabilities of marsh songbirds during morning and evening survey periods mean that investigators have more time to conduct surveys for these birds, allowing greater flexibility to increase spatial and temporal replication of surveys that could provide more precise estimates of desired population parameters.     

Evaluation of two methods to estimate and monitor bird populations

Background

Effective management depends upon accurately estimating trends in abundance of bird populations over time, and in some cases estimating abundance. Two population estimation methods, double observer (DO) and double sampling (DS), have been advocated for avian population studies and the relative merits and short-comings of these methods remain an area of debate.

Methodology/Principal Findings

We used simulations to evaluate the performances of these two population estimation methods under a range of realistic scenarios. For three hypothetical populations with different levels of clustering, we generated DO and DS population size estimates for a range of detection probabilities and survey proportions. Population estimates for both methods were centered on the true population size for all levels of population clustering and survey proportions when detection probabilities were greater than 20%. The DO method underestimated the population at detection probabilities less than 30% whereas the DS method remained essentially unbiased. The coverage probability of 95% confidence intervals for population estimates was slightly less than the nominal level for the DS method but was substantially below the nominal level for the DO method at high detection probabilities. Differences in observer detection probabilities did not affect the accuracy and precision of population estimates of the DO method. Population estimates for the DS method remained unbiased as the proportion of units intensively surveyed changed, but the variance of the estimates decreased with increasing proportion intensively surveyed.

Conclusions/Significance

The DO and DS methods can be applied in many different settings and our evaluations provide important information on the performance of these two methods that can assist researchers in selecting the method most appropriate for their particular needs.     

A Comparison of Methods for Estimating Northern Bobwhite Covey Detection Probabilities

Abstract: We compared the time-of-detection and logistic regression methods of estimating probability of detection for northern bobwhite (Colinus virginianus) coveys. Both methods are unusual in that they allow estimation of the total probability of detection (i.e., the product of the probability that a covey is available for detection [i.e., that a covey vocalizes] and detection given availability). The logistic regression method produced an average detection probability of 0.596 (SE = 0.020) and the time-of-detection method produced a detection probability estimate of 0.540 (SE = 0.086), and the 2 estimates were not significantly different. This is the first evaluation of the time-of-detection method with empirical field data. Although the time-of-detection and logistic regression method each have advantages, both can be used under appropriate conditions to improve estimates of bobwhite abundance by allowing for the estimation of detection probabilities. Improved estimates of bobwhite abundance will allow land managers to make more informed management decisions.     

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.     

Use of classical bird census transects as spatial replicates for hierarchical modeling of an avian community

New monitoring programs are often designed with some form of temporal replication to deal with imperfect detection by means of occupancy models. However, classical bird census data from earlier times often lack temporal replication, precluding detection‐corrected inferences about occupancy. Historical data have a key role in many ecological studies intended to document range shifts, and so need to be made comparable with present‐day data by accounting for detection probability. We analyze a classical bird census conducted in the region of Murcia (SE Spain) in 1991 and 1992 and propose a solution to estimating detection probability for such historical data when used in a community occupancy model: the spatial replication of subplots nested within larger plots allows estimation of detection probability. In our study, the basic sample units were 1‐km transects, which were considered spatial replicates in two aggregation schemes. We fit two Bayesian multispecies occupancy models, one for each aggregation scheme, and evaluated the linear and quadratic effect of forest cover and temperature, and a linear effect of precipitation on species occupancy probabilities. Using spatial rather than temporal replicates allowed us to obtain individual species occupancy probabilities and species richness accounting for imperfect detection. Species‐specific occupancy and community size decreased with increasing annual mean temperature. Both aggregation schemes yielded estimates of occupancy and detectability that were highly correlated for each species, so in the design of future surveys ecological reasons and cost‐effective sampling designs should be considered to select the most suitable aggregation scheme. In conclusion, the use of spatial replication may often allow historical survey data to be applied formally hierarchical occupancy models and be compared with modern‐day data of the species community to analyze global change process.     

New technologies in the mix: Assessing N‐mixture models for abundance estimation using automated detection data from drone surveys

1. Reliable estimates of abundance are critical in effectively managing threatened species, but the feasibility of integrating data from wildlife surveys completed using advanced technologies such as remotely piloted aircraft systems (RPAS) and machine learning into abundance estimation methods such as N‐mixture modeling is largely unknown due to the unique sources of detection errors associated with these technologies.
2. We evaluated two modeling approaches for estimating the abundance of koalas detected automatically in RPAS imagery: (a) a generalized N‐mixture model and (b) a modified Horvitz–Thompson (H‐T) estimator method combining generalized linear models and generalized additive models for overall probability of detection, false detection, and duplicate detection. The final estimates from each model were compared to the true number of koalas present as determined by telemetry‐assisted ground surveys.
3. The modified H‐T estimator approach performed best, with the true count of koalas captured within the 95% confidence intervals around the abundance estimates in all 4 surveys in the testing dataset (n = 138 detected objects), a particularly strong result given the difficulty in attaining accuracy found with previous methods.
4. The results suggested that N‐mixture models in their current form may not be the most appropriate approach to estimating the abundance of wildlife detected in RPAS surveys with automated detection, and accurate estimates could be made with approaches that account for spurious detections.

     

Separating Components of the Detection Process With Combined Methods: An Example With Northern Bobwhite

Abstract: There are various methods of estimating detection probabilities for avian point counts. Distance and multiple-observer methods require the sometimes unlikely assumption that all birds in the population are available (i.e., sing or are visible) during a count, but the time-of-detection method allows for the possibility that some birds are unavailable during the count. We combined the dependent double-observer method with the time-of-detection method and obtained field-based estimates of the components of detection probability for northern bobwhite (Colinus virginianus). Our approach was a special case of Pollock's robust capture-recapture design where the probability that a bird does not sing is analogous to the probability that an animal is a temporary emigrant. Top models indicated that observers' detection probabilities were similar (0.78–0.84) if bobwhite were available, but bobwhite only had an approximately 0.61 probability of being available during a 2.5-minute sampling interval. Additionally, observers' detection probabilities increased substantially after the initial encounter with an individual bobwhite (analogous to a trap-happy response on the part of the observer). A simulated data set revealed that the combined method was precise when availability and detection given availability were substantially lower. Combined methods approaches can provide critical information for researchers and land managers to make decisions regarding survey length and personnel requirements for point-count-based surveys.     

Using substitution matrices to estimate probability distributions for biological sequences.

Accurately estimating probabilities from observations is important for probabilistic-based approaches to problems in computational biology. In this paper we present a biologically-motivated method for estimating probability distributions over discrete alphabets from observations using a mixture model of common ancestors. The method is an extension of substitution matrix-based probability estimation methods. In contrast to previous such methods, our method has a simple Bayesian interpretation and has the advantage over Dirichlet mixtures that it is both effective and simple to compute for large alphabets. The method is applied to estimate amino acid probabilities based on observed counts in an alignment and is shown to perform comparably to previous methods. The method is also applied to estimate probability distributions over protein families and improves protein classification accuracy.     

Using mixtures of common ancestors for estimating the probabilities of discrete events in biological sequences

Accurately estimating probabilities from observations is important for probabilistic-based approaches to problems in computational biology. In this paper we present a biologically-motivated method for estimating probability distributions over discrete alphabets from observations using a mixture model of common ancestors. The method is an extension of substitution matrix-based probability estimation methods. In contrast to previous such methods, our method has a simple Bayesian interpretation and has the advantage over Dirichlet mixtures that it is both effective and simple to compute for large alphabets. The method is applied to estimate amino acid probabilities based on observed counts in an alignment and is shown to perform comparably to previous methods. The method is also applied to estimate probability distributions over protein families and improves protein classification accuracy.     

Double-Observer Line Transect Methods: Levels of Independence

Summary .  Double-observer line transect methods are becoming increasingly widespread, especially for the estimation of marine mammal abundance from aerial and shipboard surveys when detection of animals on the line is uncertain. The resulting data supplement conventional distance sampling data with two-sample mark–recapture data. Like conventional mark–recapture data, these have inherent problems for estimating abundance in the presence of heterogeneity. Unlike conventional mark–recapture methods, line transect methods use knowledge of the distribution of a covariate, which affects detection probability (namely, distance from the transect line) in inference. This knowledge can be used to diagnose unmodeled heterogeneity in the mark–recapture component of the data. By modeling the covariance in detection probabilities with distance, we show how the estimation problem can be formulated in terms of different levels of independence. At one extreme, full independence is assumed, as in the Petersen estimator (which does not use distance data); at the other extreme, independence only occurs in the limit as detection probability tends to one. Between the two extremes, there is a range of models, including those currently in common use, which have intermediate levels of independence. We show how this framework can be used to provide more reliable analysis of double-observer line transect data. We test the methods by simulation, and by analysis of a dataset for which true abundance is known. We illustrate the approach through analysis of minke whale sightings data from the North Sea and adjacent waters.     

Comparison of point counts and territory mapping for detecting effects of forest management on songbirds

Point counts are commonly used to assess changes in bird abundance, including analytical approaches such as distance sampling that estimate density. Point‐count methods have come under increasing scrutiny because effects of detection probability and field error are difficult to quantify. For seven forest songbirds, we compared fixed‐radii counts (50 m and 100 m) and density estimates obtained from distance sampling to known numbers of birds determined by territory mapping. We applied point‐count analytic approaches to a typical forest management question and compared results to those obtained by territory mapping. We used a before–after control impact (BACI) analysis with a data set collected across seven study areas in the central Appalachians from 2006 to 2010. Using a 50‐m fixed radius, variance in error was at least 1.5 times that of the other methods, whereas a 100‐m fixed radius underestimated actual density by >3 territories per 10 ha for the most abundant species. Distance sampling improved accuracy and precision compared to fixed‐radius counts, although estimates were affected by birds counted outside 10‐ha units. In the BACI analysis, territory mapping detected an overall treatment effect for five of the seven species, and effects were generally consistent each year. In contrast, all point‐count methods failed to detect two treatment effects due to variance and error in annual estimates. Overall, our results highlight the need for adequate sample sizes to reduce variance, and skilled observers to reduce the level of error in point‐count data. Ultimately, the advantages and disadvantages of different survey methods should be considered in the context of overall study design and objectives, allowing for trade‐offs among effort, accuracy, and power to detect treatment effects.     

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.     

Incorporating capture heterogeneity in the estimation of autoregressive coefficients of animal population dynamics using capture–recapture data

Population dynamic models combine density dependence and environmental effects. Ignoring sampling uncertainty might lead to biased estimation of the strength of density dependence. This is typically addressed using state‐space model approaches, which integrate sampling error and population process estimates. Such models seldom include an explicit link between the sampling procedures and the true abundance, which is common in capture–recapture settings. However, many of the models proposed to estimate abundance in the presence of capture heterogeneity lead to incomplete likelihood functions and cannot be straightforwardly included in state‐space models. We assessed the importance of estimating sampling error explicitly by taking an intermediate approach between ignoring uncertainty in abundance estimates and fully specified state‐space models for density‐dependence estimation based on autoregressive processes. First, we estimated individual capture probabilities based on a heterogeneity model for a closed population, using a conditional multinomial likelihood, followed by a Horvitz–Thompson estimate for abundance. Second, we estimated coefficients of autoregressive models for the log abundance. Inference was performed using the methodology of integrated nested Laplace approximation (INLA). We performed an extensive simulation study to compare our approach with estimates disregarding capture history information, and using R‐package VGAM, for different parameter specifications. The methods were then applied to a real data set of gray‐sided voles Myodes rufocanus from Northern Norway. We found that density‐dependence estimation was improved when explicitly modeling sampling error in scenarios with low process variances, in which differences in coverage reached up to 8% in estimating the coefficients of the autoregressive processes. In this case, the bias also increased assuming a Poisson distribution in the observational model. For high process variances, the differences between methods were small and it appeared less important to model heterogeneity.     

On a log-linear approach to detecting ecological interactions in monitored populations

We discuss a log-linear model for series of regular bird counts taken at a number of survey sites. The model is parameterized in terms of annual growth rates rather than actual indices of abundance, as is more frequently done. This not only permits easy estimation of and inference about these rates, but also allows us to model the effects upon population growth of covariates, such as the local presence of a competitor or predator, which may themselves vary in space and over time. A recursive relationship permits the expected count at a site to be functionally dependent upon the expected count at the previous visit. We discuss the advantages of using this relationship, rather than replacing the latter with their observed counterparts, as has been used previously.     

Mixture models for estimating the size of a closed population when capture rates vary among individuals

We develop a parameterization of the beta-binomial mixture that provides sensible inferences about the size of a closed population when probabilities of capture or detection vary among individuals. Three classes of mixture models (beta-binomial, logistic-normal, and latent-class) are fitted to recaptures of snowshoe hares for estimating abundance and to counts of bird species for estimating species richness. In both sets of data, rates of detection appear to vary more among individuals (animals or species) than among sampling occasions or locations. The estimates of population size and species richness are sensitive to model-specific assumptions about the latent distribution of individual rates of detection. We demonstrate using simulation experiments that conventional diagnostics for assessing model adequacy, such as deviance, cannot be relied on for selecting classes of mixture models that produce valid inferences about population size. Prior knowledge about sources of individual heterogeneity in detection rates, if available, should be used to help select among classes of mixture models that are to be used for inference.     

Estimating red deer abundance in a wide range of management situations in Mediterranean habitats

When human interventions interfere with the natural regulation of wildlife populations by favouring some species, overabundance can emerge. We evaluated different methods of estimating red deer abundance in a wide range of population densities from southern Spain. Distance sampling estimates were used as the reference method across 22 localities and were compared with two kilometric abundance indices (KAIs), four indices based on pellet group counts and two browsing indices (BWIs). The average red deer density estimated by distance sampling was 19.51±3.19 deer per 100 ha, showing a wide range across the study area (0.04–66.77). Distance sampling estimates correlated with the KAIs, pellet group-based index and the BWI. The agreement with distance sampling improved when groups were used instead of individuals in the KAIs, when the minimum pellet group size was fixed at 20 pellets in the dropping counts, and when only highly palatable species were used in the BWI. Thus, several direct and indirect methods can estimate red deer abundance in Mediterranean habitats from Southern Spain with appropriate modifications.