Effect of detection heterogeneity in occupancy‐detection models: an experimental test of time‐to‐first‐detection methods

Imperfect detection can bias estimates of site occupancy in ecological surveys but can be corrected by estimating detection probability. Time‐to‐first‐detection (TTD) occupancy models have been proposed as a cost–effective survey method that allows detection probability to be estimated from single site visits. Nevertheless, few studies have validated the performance of occupancy‐detection models by creating a situation where occupancy is known, and model outputs can be compared with the truth. We tested the performance of TTD occupancy models in the face of detection heterogeneity using an experiment based on standard survey methods to monitor koala Phascolarctos cinereus populations in Australia. Known numbers of koala faecal pellets were placed under trees, and observers, uninformed as to which trees had pellets under them, carried out a TTD survey. We fitted five TTD occupancy models to the survey data, each making different assumptions about detectability, to evaluate how well each estimated the true occupancy status. Relative to the truth, all five models produced strongly biased estimates, overestimating detection probability and underestimating the number of occupied trees. Despite this, goodness‐of‐fit tests indicated that some models fitted the data well, with no evidence of model misfit. Hence, TTD occupancy models that appear to perform well with respect to the available data may be performing poorly. The reason for poor model performance was unaccounted for heterogeneity in detection probability, which is known to bias occupancy‐detection models. This poses a problem because unaccounted for heterogeneity could not be detected using goodness‐of‐fit tests and was only revealed because we knew the experimentally determined outcome. A challenge for occupancy‐detection models is to find ways to identify and mitigate the impacts of unobserved heterogeneity, which could unknowingly bias many models.     

Conditionally autoregressive models improve occupancy analyses of autocorrelated data: An example with environmental DNA

Site occupancy‐detection models (SODMs) are statistical models widely used for biodiversity surveys where imperfect detection of species occurs. For instance, SODMs are increasingly used to analyse environmental DNA (eDNA) data, taking into account the occurrence of both false‐positive and false‐negative errors. However, species occurrence data are often characterized by spatial and temporal autocorrelation, which might challenge the use of standard SODMs. Here we reviewed the literature of eDNA biodiversity surveys and found that most of studies do not take into account spatial or temporal autocorrelation. We then demonstrated how the analysis of data with spatial or temporal autocorrelation can be improved by using a conditionally autoregressive SODM, and show its application to environmental DNA data. We tested the autoregressive model on both simulated and real data sets, including chronosequences with different degrees of autocorrelation, and a spatial data set on a virtual landscape. Analyses of simulated data showed that autoregressive SODMs perform better than traditional SODMs in the estimation of key parameters such as true‐/false‐positive rates and show a better discrimination capacity (e.g., higher true skill statistics). The usefulness of autoregressive SODMs was particularly high in data sets with strong autocorrelation. When applied to real eDNA data sets (eDNA from lake sediment cores and freshwater), autoregressive SODM provided more precise estimation of true‐/false‐positive rates, resulting in more reasonable inference of occupancy states. Our results suggest that analyses of occurrence data, such as many applications of eDNA, can be largely improved by applying conditionally autoregressive specifications to SODMs.     

Large-scale grid-based detection in occupancy surveys of a threatened small mammal: A comparison of two non-invasive methods

Monitoring the status and trends of wildlife is key to understand how species respond to natural and human-derived threats, and to evaluate and improve conservation planning and management. Large-scale, grid-based assessment of species distribution, abundance, and population trends over time is an important component of biodiversity monitoring. However, such assessments still present important challenges related, for instance, to how the choice of the sampling method may affect species detectability and thus, overall data accuracy. Here, we address this issue, focusing on the Cabrera vole (Microtus cabrerae), a threatened small mammal, listed in the Habitats Directive (Annexes II and IV), hence requiring regular evaluation of its population status and trends. We used occupancy modelling to estimate method-specific detection probability of the species over large-scale, grid-based (10 × 10 km²) surveys relying on two non-invasive sampling techniques: sign surveys and owl pellet analysis. Results provided evidence for a greater cost-effectiveness of sign surveys compared to owl pellet analysis for detecting the species in occupancy surveys, suggesting that large-scale population monitoring of Cabrera voles (or other species also producing easily identifiable signs of their presence) may fairly rely on sign-surveys. Overall, our study supported the view that while owl pellet analysis provides a valuable option when the aim is to assess small mammal assemblages (i.e. multiple species) in a region, other complementary methods may be required to increase the detection probability of certain species that because of their secretive behaviour or rarity remain less predated by owls. We thus argue that the choice of the sampling method should be context-dependent and evaluated based on the study aims, the surveyed area (i.e. local factors), the target species (i.e. life history traits) and the available resources. Based on our results we recommend that researchers and managers explore survey-design trade-offs to ensure the proposed designs have sufficient power to detect real population trends.     

Statistical approaches to account for false‐positive errors in environmental DNA samples

Environmental DNA (eDNA) sampling is prone to both false‐positive and false‐negative errors. We review statistical methods to account for such errors in the analysis of eDNA data and use simulations to compare the performance of different modelling approaches. Our simulations illustrate that even low false‐positive rates can produce biased estimates of occupancy and detectability. We further show that removing or classifying single PCR detections in an ad hoc manner under the suspicion that such records represent false positives, as sometimes advocated in the eDNA literature, also results in biased estimation of occupancy, detectability and false‐positive rates. We advocate alternative approaches to account for false‐positive errors that rely on prior information, or the collection of ancillary detection data at a subset of sites using a sampling method that is not prone to false‐positive errors. We illustrate the advantages of these approaches over ad hoc classifications of detections and provide practical advice and code for fitting these models in maximum likelihood and Bayesian frameworks. Given the severe bias induced by false‐negative and false‐positive errors, the methods presented here should be more routinely adopted in eDNA studies.     

Estimating dung decay rates of roe deer (<Emphasis Type="Italic">Capreolus capreolus</Emphasis>) in different habitat types of a Mediterranean ecosystem: an information theory approach

For elusive species living in concealing habitats (e.g. deer in a forest habitat), indirect methods such as faecal pellet counts are considered more practical means of estimating population density and abundance. Accurate estimation of deer density using the faecal standing crop (FSC) method necessitates the reliable estimation of the mean time to decay of pellet groups present during the survey. Mean time to decay is generally habitat specific, and separate estimations should be made for each habitat type in the study area. In a confined mountainous area of Greece, the habitat-specific mean time to decay of roe deer pellet groups was estimated by locating and marking fresh pellet groups on several dates in the lead up to an FSC survey and returning to the marked signs at the time of the survey to record whether or not each pellet group had survived. Several logistic models were fitted to the data, and estimations were based on a multi-model inference (MMI) approach according to information theory. The highest mean time to decay was estimated in coniferous forests, while mid-ranged values were found in maquis shrubs, and the lowest mean time to decay was observed in open areas. MMI by model averaging, based on Akaike weights, is recommended for making robust parameter estimations and for dealing with uncertainty in model selection.     

Evaluating the interaction of faecal pellet deposition rates and DNA degradation rates to optimize sampling design for DNA‐based mark–recapture analysis of Sonoran pronghorn

Knowledge of population demographics is important for species management but can be challenging in low‐density, wide‐ranging species. Population monitoring of the endangered Sonoran pronghorn (Antilocapra americana sonoriensis) is critical for assessing the success of recovery efforts, and noninvasive DNA sampling (NDS) could be more cost‐effective and less intrusive than traditional methods. We evaluated faecal pellet deposition rates and faecal DNA degradation rates to maximize sampling efficiency for DNA‐based mark–recapture analyses. Deposition data were collected at five watering holes using sampling intervals of 1–7 days and averaged one pellet pile per pronghorn per day. To evaluate nuclear DNA (nDNA) degradation, 20 faecal samples were exposed to local environmental conditions and sampled at eight time points from one to 124 days. Average amplification success rates for six nDNA microsatellite loci were 81% for samples on day one, 63% by day seven, 2% by day 14 and 0% by day 60. We evaluated the efficiency of different sampling intervals (1–10 days) by estimating the number of successful samples, success rate of individual identification and laboratory costs per successful sample. Cost per successful sample increased and success and efficiency declined as the sampling interval increased. Results indicate NDS of faecal pellets is a feasible method for individual identification, population estimation and demographic monitoring of Sonoran pronghorn. We recommend collecting samples >7 days old and estimate that a sampling interval of 4–7 days in summer conditions (i.e. extreme heat and exposure to UV light) will achieve desired sample sizes for mark–recapture analysis while also maximizing efficiency.     

Development and field validation of a regional,management‐scale habitat model: A koala Phascolarctos cinereus case study

Species distribution models have great potential to efficiently guide management for threatened species, especially for those that are rare or cryptic. We used MaxEnt to develop a regional‐scale model for the koala Phascolarctos cinereus at a resolution (250 m) that could be used to guide management. To ensure the model was fit for purpose, we placed emphasis on validating the model using independently‐collected field data. We reduced substantial spatial clustering of records in coastal urban areas using a 2‐km spatial filter and by modeling separately two subregions separated by the 500‐m elevational contour. A bias file was prepared that accounted for variable survey effort. Frequency of wildfire, soil type, floristics and elevation had the highest relative contribution to the model, while a number of other variables made minor contributions. The model was effective in discriminating different habitat suitability classes when compared with koala records not used in modeling. We validated the MaxEnt model at 65 ground‐truth sites using independent data on koala occupancy (acoustic sampling) and habitat quality (browse tree availability). Koala bellows (n = 276) were analyzed in an occupancy modeling framework, while site habitat quality was indexed based on browse trees. Field validation demonstrated a linear increase in koala occupancy with higher modeled habitat suitability at ground‐truth sites. Similarly, a site habitat quality index at ground‐truth sites was correlated positively with modeled habitat suitability. The MaxEnt model provided a better fit to estimated koala occupancy than the site‐based habitat quality index, probably because many variables were considered simultaneously by the model rather than just browse species. The positive relationship of the model with both site occupancy and habitat quality indicates that the model is fit for application at relevant management scales. Field‐validated models of similar resolution would assist in guiding management of conservation‐dependent species.     

Determining Occurrence Dynamics when False Positives Occur: Estimating the Range Dynamics of Wolves from Public Survey Data

Large-scale presence-absence monitoring programs have great promise for many conservation applications. Their value can be limited by potential incorrect inferences owing to observational errors, especially when data are collected by the public. To combat this, previous analytical methods have focused on addressing non-detection from public survey data. Misclassification errors have received less attention but are also likely to be a common component of public surveys, as well as many other data types. We derive estimators for dynamic occupancy parameters (extinction and colonization), focusing on the case where certainty can be assumed for a subset of detections. We demonstrate how to simultaneously account for non-detection (false negatives) and misclassification (false positives) when estimating occurrence parameters for gray wolves in northern Montana from 2007–2010. Our primary data source for the analysis was observations by deer and elk hunters, reported as part of the state’s annual hunter survey. This data was supplemented with data from known locations of radio-collared wolves. We found that occupancy was relatively stable during the years of the study and wolves were largely restricted to the highest quality habitats in the study area. Transitions in the occupancy status of sites were rare, as occupied sites almost always remained occupied and unoccupied sites remained unoccupied. Failing to account for false positives led to over estimation of both the area inhabited by wolves and the frequency of turnover. The ability to properly account for both false negatives and false positives is an important step to improve inferences for conservation from large-scale public surveys. The approach we propose will improve our understanding of the status of wolf populations and is relevant to many other data types where false positives are a component of observations.     

A model‐based solution for observational errors in laboratory studies

Molecular techniques for detecting microorganisms, macroorganisms and infectious agents are susceptible to false‐negative and false‐positive errors. If left unaddressed, these observational errors may yield misleading inference concerning occurrence, prevalence, sensitivity, specificity and covariate relationships. Occupancy models are widely used to account for false‐negative errors and more recently have even been used to address false‐positive errors, too. Current modelling options assume false‐positive errors only occur in truly negative samples, an assumption that yields biased inference concerning detection because a positive sample could be classified as such not because the target agent was successfully detected, but rather due to a false‐positive test result. We present an extension to the occupancy modelling framework that allows false‐positive errors in both negative and positive samples, thereby providing unbiased inference concerning occurrence and detection, as well as reliable conclusions about the efficacy of sampling designs, handling protocols and diagnostic tests. We apply the model to simulated data, showing that it recovers known parameters and outperforms other approaches that are commonly used when confronted with observation errors. We then apply the model to an experimental data set on Batrachochytrium dendrobatidis, a pathogenic fungus that is implicated in the global decline or extinction of hundreds of amphibian species. The model‐based approach we present is not only useful for obtaining reliable inference when data are contaminated with observational errors, but also eliminates the need for establishing arbitrary thresholds or decision rules that have hidden and unintended consequences.     

Improving geographically extensive acoustic survey designs for modeling species occurrence with imperfect detection and misidentification

Acoustic recording units (ARUs) enable geographically extensive surveys of sensitive and elusive species. However, a hidden cost of using ARU data for modeling species occupancy is that prohibitive amounts of human verification may be required to correct species identifications made from automated software. Bat acoustic studies exemplify this challenge because large volumes of echolocation calls could be recorded and automatically classified to species. The standard occupancy model requires aggregating verified recordings to construct confirmed detection/non‐detection datasets. The multistep data processing workflow is not necessarily transparent nor consistent among studies. We share a workflow diagramming strategy that could provide coherency among practitioners. A false‐positive occupancy model is explored that accounts for misclassification errors and enables potential reduction in the number of confirmed detections. Simulations informed by real data were used to evaluate how much confirmation effort could be reduced without sacrificing site occupancy and detection error estimator bias and precision. We found even under a 50% reduction in total confirmation effort, estimator properties were reasonable for our assumed survey design, species‐specific parameter values, and desired precision. For transferability, a fully documented r package, OCacoustic, for implementing a false‐positive occupancy model is provided. Practitioners can apply OCacoustic to optimize their own study design (required sample sizes, number of visits, and confirmation scenarios) for properly implementing a false‐positive occupancy model with bat or other wildlife acoustic data. Additionally, our work highlights the importance of clearly defining research objectives and data processing strategies at the outset to align the study design with desired statistical inferences.     

Dynamic occupancy models for analyzing species' range dynamics across large geographic scales

Large‐scale biodiversity data are needed to predict species' responses to global change and to address basic questions in macroecology. While such data are increasingly becoming available, their analysis is challenging because of the typically large heterogeneity in spatial sampling intensity and the need to account for observation processes. Two further challenges are accounting for spatial effects that are not explained by covariates, and drawing inference on dynamics at these large spatial scales. We developed dynamic occupancy models to analyze large‐scale atlas data. In addition to occupancy, these models estimate local colonization and persistence probabilities. We accounted for spatial autocorrelation using conditional autoregressive models and autologistic models. We fitted the models to detection/nondetection data collected on a quarter‐degree grid across southern Africa during two atlas projects, using the hadeda ibis (Bostrychia hagedash) as an example. The model accurately reproduced the range expansion between the first (SABAP1: 1987–1992) and second (SABAP2: 2007–2012) Southern African Bird Atlas Project into the drier parts of interior South Africa. Grid cells occupied during SABAP1 generally remained occupied, but colonization of unoccupied grid cells was strongly dependent on the number of occupied grid cells in the neighborhood. The detection probability strongly varied across space due to variation in effort, observer identity, seasonality, and unexplained spatial effects. We present a flexible hierarchical approach for analyzing grid‐based atlas data using dynamical occupancy models. Our model is similar to a species' distribution model obtained using generalized additive models but has a number of advantages. Our model accounts for the heterogeneous sampling process, spatial correlation, and perhaps most importantly, allows us to examine dynamic aspects of species ranges.     

Experimental Investigation of Observation Error in Anuran Call Surveys

ABSTRACT Occupancy models that account for imperfect detection are often used to monitor anuran and songbird species occurrence. However, presence—absence data arising from auditory detections may be more prone to observation error (e.g., false-positive detections) than are sampling approaches utilizing physical captures or sightings of individuals. We conducted realistic, replicated field experiments using a remote broadcasting system to simulate simple anuran call surveys and to investigate potential factors affecting observation error in these studies. Distance, time, ambient noise, and observer abilities were the most important factors explaining false-negative detections. Distance and observer ability were the best overall predictors of false-positive errors, but ambient noise and competing species also affected error rates for some species. False-positive errors made up 5% of all positive detections, with individual observers exhibiting false-positive rates between 0.5% and 14%. Previous research suggests false-positive errors of these magnitudes would induce substantial positive biases in standard estimators of species occurrence, and we recommend practices to mitigate for false positives when developing occupancy monitoring protocols that rely on auditory detections. These recommendations include additional observer training, limiting the number of target species, and establishing distance and ambient noise thresholds during surveys.     

Using bioacoustics to enhance the efficiency of spotted owl surveys and facilitate forest restoration

The California spotted owl (Strix occidentalis occidentalis) is an older-forest associated species that resides at the center of forest management planning in the Sierra Nevada and Southern California, USA, which are experiencing increasingly large and severe wildfires and drought-related tree mortality. We leveraged advances in passive acoustic survey technologies to develop an acoustically assisted survey design that could increase the efficiency and effectiveness of project-level surveys for spotted owls, allowing surveys to be completed in a single year instead of in multiple years. We deployed an array of autonomous recording units (ARUs) across a landscape and identified spotted owl vocalizations in the resulting audio using BirdNET. We then evaluated spatio-temporal patterns in spotted owl vocalizations near occupied territories and the ability of a crew naïve to the location of occupied territories to locate spotted owls based on patterns of acoustic detections. After only 3 weeks of acoustic surveys, ≥1 ARU within 750 m of all 17 occupied territories obtained spotted owl detections across ≥2 nights. When active surveys using broadcast calling were conducted near ARUs with spotted owl detections by surveyors naïve to territory occupancy status and locations, surveyors located owls in 93% to 100% of occupied territories with ≤3 surveys. To further improve the efficiency of spotted owl surveys, we developed a statistical model to identify and prioritize areas across the Sierra Nevada for different survey methods (active only, acoustically assisted, no surveys) based on the expected probability of occupancy predicted from remotely sensed measurements of tree height and historical occupancy. Depending on managers' tolerance for false negatives, this model could help identify large areas that might not benefit from surveys based on low expected occupancy probabilities and areas where acoustically assisted surveys might enhance survey effectiveness and efficiency. Collectively, these findings can help managers streamline the survey process and thus increase the pace of forest restoration while minimizing potential near-term adverse effects on California spotted owls.     

Patch occupancy of North American mammals: is patchiness in the eye of the beholder?

Aim Intraspecific variation in patch occupancy often is related to physical features of a landscape, such as the amount and distribution of habitat. However, communities occupying patchy environments typically exhibit non‐random distributions in which local assemblages of species‐poor patches are nested subsets of assemblages occupying more species‐rich patches. Nestedness of local communities implies interspecific differences in sensitivity to patchiness. Several hypotheses have been proposed to explain interspecific variation in responses to patchiness within a community, including differences in (1) colonization ability, (2) extinction proneness, (3) tolerance to disturbance, (4) sociality and (5) level of adaptation to prevailing environmental conditions. We used data on North American mammals to compare the performance of these ‘ecological’ hypotheses and the ‘physical landscape’ hypothesis. We then compared the best of these models against models that scaled landscape structure to ecologically relevant attributes of individual species. Location North America. Methods We analysed data on prevalence (i.e. proportion of patches occupied in a network of patches) and occupancy for 137 species of non‐volant mammals and twenty networks consisting of four to seventy‐five patches. Insular and terrestrial networks exhibited significantly different mean levels of prevalence and occupancy and thus were analysed separately. Indicator variables at ordinal and family levels were included in models to correct for effects caused by phylogeny. Akaike's information criterion was used in conjunction with ordinary least squares and logistic regression to compare hypotheses. Results A patch network's physical structure, indexed using patch area and isolation, received the greatest support among models predicting the prevalence of species on insular networks. Niche breadth (diet and habitat) received the greatest support for predicting prevalence of species occupying terrestrial networks. For both insular and terrestrial systems, physical features (patch area and isolation) received greater support than any of the ecological hypotheses for predicting species occupancy of individual patches. For terrestrial systems, scaling patch area by its suitability to a focal species and by individual area requirements of the species, and scaling patch isolation by species‐specific dispersal ability and niche breadth, resulted in models of patch occupancy that were superior to models relying solely on physical landscape features. For all selected models, unexplained levels of variation were high. Main conclusions Stochasticity dominated the systems we studied, indicating that random events are probably quite important in shaping local communities. With respect to deterministic factors, our results suggest that forces affecting species prevalence and occupancy may differ between insular and terrestrial systems. Physical features of insular systems appeared to swamp ecological differences among species in determining prevalence and occupancy, whereas species with broad niches were disproportionately represented in terrestrial networks. We hypothesize that differential extinction over long time periods in highly variable networks has driven nestedness of mammalian communities on islands, whereas differential colonization over shorter time‐scales in more homogeneous networks probably governed the local structure of terrestrial communities. Our results also demonstrate that integration of a species' ecological traits with physical features of a patch network is superior to reliance on either factor separately when attempting to predict the species' probability of patch occupancy in terrestrial systems.     

Modeling stream fish distributions using interval‐censored detection times

Controlling for imperfect detection is important for developing species distribution models (SDMs). Occupancy‐detection models based on the time needed to detect a species can be used to address this problem, but this is hindered when times to detection are not known precisely. Here, we extend the time‐to‐detection model to deal with detections recorded in time intervals and illustrate the method using a case study on stream fish distribution modeling. We collected electrofishing samples of six fish species across a Mediterranean watershed in Northeast Portugal. Based on a Bayesian hierarchical framework, we modeled the probability of water presence in stream channels, and the probability of species occupancy conditional on water presence, in relation to environmental and spatial variables. We also modeled time‐to‐first detection conditional on occupancy in relation to local factors, using modified interval‐censored exponential survival models. Posterior distributions of occupancy probabilities derived from the models were used to produce species distribution maps. Simulations indicated that the modified time‐to‐detection model provided unbiased parameter estimates despite interval‐censoring. There was a tendency for spatial variation in detection rates to be primarily influenced by depth and, to a lesser extent, stream width. Species occupancies were consistently affected by stream order, elevation, and annual precipitation. Bayesian P‐values and AUCs indicated that all models had adequate fit and high discrimination ability, respectively. Mapping of predicted occupancy probabilities showed widespread distribution by most species, but uncertainty was generally higher in tributaries and upper reaches. The interval‐censored time‐to‐detection model provides a practical solution to model occupancy‐detection when detections are recorded in time intervals. This modeling framework is useful for developing SDMs while controlling for variation in detection rates, as it uses simple data that can be readily collected by field ecologists.     

Spatially explicit estimation of occupancy,detection probability and survey effort needed to inform conservation planning

Aim It is increasingly recognized the importance of accounting for imperfect detection in species distribution modelling and conservation planning. However, the integration of detectability into a spatially explicit frame has received little attention. We aim (1) to show how to develop distribution maps of both detection probability and survey effort required to reliably determine a species presence/absence and (2) to increase awareness of the spatial variation of detection error inherent in studies of species occurrence. Location North‐western Spain. Methods  We registered the presence/absence of the endangered Egyptian vulture (Neophron percnopterus) in 213 surveys performed in 40 of 104 territories once known to be occupied. We model simultaneously both detection probability and occurrence, using site occupancy modelling. With the resulting regression equations, we developed distribution maps of both detection probability and required sampling effort throughout the area. Results Of the studied territories, 72.5% were detected as occupied, but after accounting for imperfect detection, the proportion of sites truly occupied was 79%. Detectability decreased in territories with higher topographical irregularity and increased with both the time of day of the survey and the progress of the season. Spatial distribution of detectability showed a mainly north–south gradient following the distribution of slope in the area. The likelihood of occupancy increased with rockier, less forested surface and less topographical irregularity within the territory. A minimum of five surveys, on average, are needed to assess, with 95% probability, the occupancy status of a site, ranging from ≤ 3 to > 24 visits/territory depending on survey‐ and site‐specific features. Main conclusions Accounting for detectability and its sources of variation allows us to elaborate distribution maps of detectability‐based survey effort. These maps are useful tools to reliably assess (e.g. with 95% probability) occupancy status throughout a landscape and provide guidance for species conservation planning.     

Estimating Detection Probabilities from Sign of Collared Peccary

ABSTRACT Determining presence or absence of collared peccaries (Pecari tajacu) from surveys of sign (tracks and feces) requires information on whether animals in sample units are detected. We estimated detection probabilities of collared peccary from sign surveys using occupancy models. Because it was unlikely that residence status of collared peccary in sampling units remained the same over a survey season, which is a primary assumption of occupancy models, we first determined the time interval for which to pool data. We then examined the influence of rainfall and peccary abundance on detection probabilities. We placed 90 sign stations (25-m-diam circular plots) throughout Chaparral Wildlife Management Area, south Texas, USA. We surveyed plots weekly for the presence or non-presence of collared peccary during 2 11-week sampling seasons in spring and fall 2003. We examined sign data weekly and we pooled the data in intervals from 2 weeks to 5 weeks. Estimates of detection probabilities increased from 1 week to 3 weeks of pooled data and leveled off thereafter. We needed a 3-week time interval to meet the assumption of unchanging residence status. Using sign data pooled in 3-week increments, detection probabilities were influenced by areas that differed in peccary abundance, but they were not influenced by rainfall. Estimates of detection probabilities ranged from 0.43 to 0.77 for 3-week time intervals. Sign surveys and occupancy modeling of data can be used to measure spatial patterns of collared peccary in south Texas as long as multiple 3-week time intervals are sampled.     

Factors affecting detectability of river otters during sign surveys

Sign surveys are commonly used to study and monitor wildlife species but may be flawed when surveys are conducted only once and cover short distances, which can lead to a lack of accountability for false absences. Multiple observers surveyed for river otter (Lontra canadensis) scat and tracks along stream and reservoir shorelines at 110 randomly selected sites in eastern Kansas from January to April 2008 and 2009 to determine if detection probability differed among substrates, sign types, observers, survey lengths, and near access points. We estimated detection probabilities (p) of river otters using occupancy models in Program PRESENCE. Mean detection probability for a 400-m survey was highest in mud substrates (p = 0.60) and lowest in snow (p = 0.18) and leaf litter substrates (p = 0.27). Scat had a higher detection probability (p = 0.53) than tracks (p = 0.18), and experienced observers had higher detection probabilities (p > 0.71) than novice observers (p < 0.55). Detection probabilities increased almost 3-fold as survey length increased from 200 m to 1,000 m, and otter sign was not concentrated near access points. After accounting for imperfect detection, our estimates of otter site occupancy based on a 400-m survey increased >3-fold, providing further evidence of the potential negative bias that can occur in estimates from sign surveys when imperfect detection is not addressed. Our study identifies areas for improvement in sign survey methodologies and results are applicable for sign surveys commonly used for many species across a range of habitats. © 2010 The Wildlife Society     

Designing occupancy surveys and interpreting non‐detection when observations are imperfect

Aim Conservation practitioners use biological surveys to ascertain whether or not a site is occupied by a particular species. Widely used statistical methods estimate the probability that a species will be detected in a survey of an occupied site. However, these estimates of detection probability are alone not sufficient to calculate the probability that a species is present given that it was not detected. The aim of this paper is to demonstrate methods for correctly calculating (1) the probability a species occupies a site given one or more non‐detections, and (2) the number of sequential non‐detections necessary to assert, with a pre‐specified confidence, that a species is absent from a site. Location Occupancy data for a tree frog in eastern Australia serve to illustrate methods that may be applied anywhere species’ occupancy data are used and detection probabilities are < 1. Methods Building on Bayesian expressions for the probability that a site is occupied by a species when it is not detected, and the number of non‐detections necessary to assert absence with a pre‐specified confidence, we estimate occupancy probabilities across tree frog survey locations, drawing on information about where and when the species was detected during surveys. Results We show that the number of sequential non‐detections necessary to assert that a species is absent increases nonlinearly with the prior probability of occupancy, the probability of detection if present, and the desired level of confidence about absence. Main conclusions If used more widely, the Bayesian analytical approaches illustrated here would improve collection and interpretation of biological survey data, providing a coherent way to incorporate detection probability estimates in the design of minimum survey requirements for monitoring, impact assessment and distribution modelling.     

Twenty‐five years of change in southern African passerine diversity: nonclimatic factors of change

We analysed more than 25 years of change in passerine bird distribution in South Africa, Swaziland and Lesotho, to show that species distributions can be influenced by processes that are at least in part independent of the local strength and direction of climate change: land use and ecological succession. We used occupancy models that separate species' detection from species' occupancy probability, fitted to citizen science data from both phases of the Southern African Bird Atlas Project (1987–1996 and 2007–2013). Temporal trends in species' occupancy probability were interpreted in terms of local extinction/colonization, and temporal trends in detection probability were interpreted in terms of change in abundance. We found for the first time at this scale that, as predicted in the context of bush encroachment, closed‐savannah specialists increased where open‐savannah specialists decreased. In addition, the trend in the abundance of species a priori thought to be favoured by agricultural conversion was negatively correlated with human population density, which is in line with hypotheses explaining the decline in farmland birds in the Northern Hemisphere. In addition to climate, vegetation cover and the intensity and time since agricultural conversion constitute important predictors of biodiversity changes in the region. Their inclusion will improve the reliability of predictive models of species distribution.