Estimating site occupancy rates for aquatic plants using spatial sub-sampling designs when detection probabilities are less than one

Estimation of site occupancy rates when detection probabilities are <1 is well established in wildlife science. Data from multiple visits to a sample of sites are used to estimate detection probabilities and the proportion of sites occupied by focal species. In this article we describe how site occupancy methods can be applied to estimate occupancy rates of plants and other sessile organisms. We illustrate this approach and the pitfalls of ignoring incomplete detection using spatial data for 2 aquatic vascular plants collected under the Upper Mississippi River's Long Term Resource Monitoring Program (LTRMP). Site occupancy models considered include: a naïve model that ignores incomplete detection, a simple site occupancy model assuming a constant occupancy rate and a constant probability of detection across sites, several models that allow site occupancy rates and probabilities of detection to vary with habitat characteristics, and mixture models that allow for unexplained variation in detection probabilities. We used information theoretic methods to rank competing models and bootstrapping to evaluate the goodness-of-fit of the final models. Results of our analysis confirm that ignoring incomplete detection can result in biased estimates of occupancy rates. Estimates of site occupancy rates for 2 aquatic plant species were 19–36% higher compared to naive estimates that ignored probabilities of detection <1. Simulations indicate that final models have little bias when 50 or more sites are sampled, and little gains in precision could be expected for sample sizes >300. We recommend applying site occupancy methods for monitoring presence of aquatic species.     

Population size influences amphibian detection probability: implications for biodiversity monitoring programs

Monitoring is an integral part of species conservation. Monitoring programs must take imperfect detection of species into account in order to be reliable. Theory suggests that detection probability may be determined by population size but this relationship has not yet been assessed empirically. Population size is particularly important because it may induce heterogeneity in detection probability and thereby cause bias in estimates of biodiversity. We used a site occupancy model to analyse data from a volunteer-based amphibian monitoring program to assess how well different variables explain variation in detection probability. An index to population size best explained detection probabilities for four out of six species (to avoid circular reasoning, we used the count of individuals at a previous site visit as an index to current population size). The relationship between the population index and detection probability was positive. Commonly used weather variables best explained detection probabilities for two out of six species. Estimates of site occupancy probabilities differed depending on whether the population index was or was not used to model detection probability. The relationship between the population index and detectability has implications for the design of monitoring and species conservation. Most importantly, because many small populations are likely to be overlooked, monitoring programs should be designed in such a way that small populations are not overlooked. The results also imply that methods cannot be standardized in such a way that detection probabilities are constant. As we have shown here, one can easily account for variation in population size in the analysis of data from long-term monitoring programs by using counts of individuals from surveys at the same site in previous years. Accounting for variation in population size is important because it can affect the results of long-term monitoring programs and ultimately the conservation of imperiled species.     

Climate and the range dynamics of species with imperfect detection

Reliable predictions for species range changes require a mechanistic understanding of range dynamics in relation to environmental variation. One obstacle is that most current models are static and confound occurrence with the probability of detecting a species if it occurs at a site. Here we draw attention to recently developed occupancy models, which can be used to examine colonization and local extinction or changes in occupancy over time. These models further account for detection probabilities, which are likely to vary spatially and temporally in many datasets. Occupancy models require repeated presence/absence surveys, for example checklists used in bird atlas projects. As an example, we examine the recent range expansion of hadeda ibises (Bostrychia hagedash) in South African protected areas. Colonization exceeded local extinction in most biomes, and the probability of occurrence was related to local climate. Extensions of the basic occupancy models can estimate abundance or species richness. Occupancy models are an appealing additional tool for studying species' responses to global change.     

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.     

Predicting species distributions from checklist data using site‐occupancy models

Aim (1) To increase awareness of the challenges induced by imperfect detection, which is a fundamental issue in species distribution modelling; (2) to emphasize the value of replicate observations for species distribution modelling; and (3) to show how ‘cheap’ checklist data in faunal/floral databases may be used for the rigorous modelling of distributions by site‐occupancy models. Location Switzerland. Methods We used checklist data collected by volunteers during 1999 and 2000 to analyse the distribution of the blue hawker, Aeshna cyanea (Odonata, Aeshnidae), a common dragonfly in Switzerland. We used data from repeated visits to 1‐ha pixels to derive ‘detection histories’ and apply site‐occupancy models to estimate the ‘true’ species distribution, i.e. corrected for imperfect detection. We modelled blue hawker distribution as a function of elevation and year and its detection probability of elevation, year and season. Results The best model contained cubic polynomial elevation effects for distribution and quadratic effects of elevation and season for detectability. We compared the site‐occupancy model with a conventional distribution model based on a generalized linear model, which assumes perfect detectability (p = 1). The conventional distribution map looked very different from the distribution map obtained using site‐occupancy models that accounted for the imperfect detection. The conventional model underestimated the species distribution by 60%, and the slope parameters of the occurrence–elevation relationship were also underestimated when assuming p = 1. Elevation was not only an important predictor of blue hawker occurrence, but also of the detection probability, with a bell‐shaped relationship. Furthermore, detectability increased over the season. The average detection probability was estimated at only 0.19 per survey. Main conclusions Conventional species distribution models do not model species distributions per se but rather the apparent distribution, i.e. an unknown proportion of species distributions. That unknown proportion is equivalent to detectability. Imperfect detection in conventional species distribution models yields underestimates of the extent of distributions and covariate effects that are biased towards zero. In addition, patterns in detectability will erroneously be ascribed to species distributions. In contrast, site‐occupancy models applied to replicated detection/non‐detection data offer a powerful framework for making inferences about species distributions corrected for imperfect detection. The use of ‘cheap’ checklist data greatly enhances the scope of applications of this useful class of models.     

Multi‐season occupancy models identify biotic and abiotic factors influencing a recovering Arctic Peregrine Falcon Falco peregrinus tundrius population

Critical information for evaluating the effectiveness of management strategies for species of concern include distinguishing seldom occupied (or low‐quality) habitat from habitat that is frequently occupied and thus contributes substantially to population trends. Using multi‐season models that account for imperfect detection and a long‐term (1981–2002) dataset on migratory Arctic Peregrine Falcons Falco peregrinus tundrius nesting along the Colville River, Alaska, we quantified the effects of previous year's productivity (i.e. site quality), amount of prey habitat, topography, climate, competition and year on occupancy dynamics across two spatial scales (nest‐sites, cliffs) during recovery of the population. Initial occupancy probability was positively correlated with area of surrounding prey habitat and height of nest‐sites above the Colville River. Colonization probability was positively correlated with nest height and negatively correlated with date of snowmelt. Local extinction probability was negatively correlated with productivity, area of prey habitat and nest height. Colonization and local extinction probabilities were also positively and negatively correlated, respectively, with year. Our results suggest that nest‐sites (or cliffs) along the Colville River do not need equal protection measures. Nest‐sites and cliffs with historically higher productivity were occupied most frequently and had lower probability of local extinction. These sites were on cliffs high above the river drainage, surrounded by adequate prey habitat and with southerly aspects associated with early snowmelt and warmer microclimates in spring. Protecting these sites is likely to encourage continued occupancy by Arctic Peregrine Falcons along the Colville River and other similar areas. Our findings also illustrate the importance of evaluating fitness parameters along with climate and habitat features when analysing occupancy dynamics, particularly with a long‐term dataset spanning a range of annual climate variation.     

Red-Shouldered Hawk Occupancy Surveys in Central Minnesota,USA

ABSTRACT Forest-dwelling raptors are often difficult to detect because many species occur at low density or are secretive. Broadcasting conspecific vocalizations can increase the probability of detecting forest-dwelling raptors and has been shown to be an effective method for locating raptors and assessing their relative abundance. Recent advances in statistical techniques based on presence—absence data use probabilistic arguments to derive probability of detection when it is < 1 and to provide a model and likelihood-based method for estimating proportion of sites occupied. We used these maximum-likelihood models with data from red-shouldered hawk (Buteo lineatus) call-broadcast surveys conducted in central Minnesota, USA, in 1994–1995 and 2004–2005. Our objectives were to obtain estimates of occupancy and detection probability 1) over multiple sampling seasons (yr), 2) incorporating within-season time-specific detection probabilities, 3) with call type and breeding stage included as covariates in models of probability of detection, and 4) with different sampling strategies. We visited individual survey locations 2–9 times per year, and estimates of both probability of detection (range = 0.28-0.54) and site occupancy (range = 0.81-0.97) varied among years. Detection probability was affected by inclusion of a within-season time-specific covariate, call type, and breeding stage. In 2004 and 2005 we used survey results to assess the effect that number of sample locations, double sampling, and discontinued sampling had on parameter estimates. We found that estimates of probability of detection and proportion of sites occupied were similar across different sampling strategies, and we suggest ways to reduce sampling effort in a monitoring program.     

Spatial dynamics of an invasive bird species assessed using robust design occupancy analysis: the case of the Eurasian collared dove (Streptopelia decaocto) in France

Aim The study of the spatial dynamics of invasive species is a key issue in invasion ecology. While mathematical models are useful for predicting the extent of population expansions, they are not suitable for measuring and characterizing spatial patterns of invasion unless the probability of detection is homogeneous across the distribution range. Here, we apply recently developed statistical approaches incorporating detection uncertainty to characterize the spatial dynamics of an invasive bird species, the Eurasian collared dove (Streptopelia decaocto). Location France. Methods Data on presence/absence of doves were recorded from 1996 to 2004 over 1045 grid cells (28 × 20 km) covering the entire country. Each grid cell included five point counts spaced along a route, which was visited twice a year, allowing for an estimation of detection probability. Each route was assigned to one of six geographical regions. We used robust design occupancy analysis to assess spatial and temporal variation in parameters related to the spatial dynamics of the species. These parameters included occupancy rate, colonization and local extinction probabilities. Our inference approach was based on the selection of the most parsimonious model among competitive models parametrized with conditional probabilities. Results The probability of detecting the presence of doves on a given route was high. However, we found evidence to incorporate detection uncertainty in inference processes about spatial dynamics, since detection probability was neither perfect (i.e. it was < 1), nor constant over space and time. Results showed a clear positive trend in occupancy rate over the study period, increasing from 55% in 1996 to 76% in 2004. In addition, occupancy rate differed among regions (range: 37–79%) and further analysis showed that colonization probability by region was positively related to occupancy rate. Finally, local extinction probability was lower than colonization probability and showed a tendency to decrease over the study period. Main conclusions Our results emphasize the importance of estimating detection probabilities in order to draw proper inferences about the spatial and temporal dynamics of the invasion pattern of the collared dove. In contrast to the perceived spatial dynamics from national atlas surveys, we provide evidence that the range of this species is currently increasing in France. Other results, such as regional specificity in colonization probabilities and time variation in local extinction are consistent with expectations from invasion and metapopulation theory.     

Environmentally constrained null models: site suitability as occupancy criterion

Null models have proven to be an important quantitative tool in the search for ecological processes driving local diversity and species distribution. However, there remains an important concern that different processes, such as environmental conditions and biotic interactions may produce similar patterns in species distributions. In this paper we present an analytical protocol for incorporating habitat suitability as an occupancy criterion in null models. Our approach involves modeling species presence or absence as a function of environmental conditions, and using the estimated site-specific probabilities of occurrence as the likelihood of species occupancy of a site during the generation of "null communities". We validated this approach by showing that type I error is not affected by the use of probabilities as a site occupancy criterion and is robust against a variety of predictive performances of the species-environmental models. We describe the expected differences when contrasting classical and the environmentally constrained null models, and illustrate our approach with a data set of Dutch dune hunting spider assemblages. Together, an environmentally constrained approach to null models will provide a more robust evaluation of species associations by facilitating the distinction between mutually exclusive processes that may shape species distributions and community assembly.     

Capture–recapture models with heterogeneity to study survival senescence in the wild

Detecting senescence in wild populations and estimating its strength raise three challenges. First, in the presence of individual heterogeneity in survival probability, the proportion of high‐survival individuals increases with age. This increase can mask a senescence‐related decrease in survival probability when the probability is estimated at the population level. To accommodate individual heterogeneity we use a mixture model structure (discrete classes of individuals). Second, the study individuals can elude the observers in the field, and their detection rate can be heterogeneous. To account for detectability issues we use capture–mark–recapture (CMR) methodology, mixture models and data that provide information on individuals’ detectability. Last, emigration to non‐monitored sites can bias survival estimates, because it can occur at the end of the individuals’ histories and mimic earlier death. To model emigration we use Markovian transitions to and from an unobservable state. These different model structures are merged together using hidden Markov chain CMR models, or multievent models. Simulation studies illustrate that reliable evidence for survival senescence can be obtained using highly heterogeneous data from non site‐faithful individuals. We then design a tailored application for a dataset from a colony of black‐headed gull Chroicocephalus ridibundus. Survival probabilities do not appear individually variable, but evidence for survival senescence becomes significant only when accounting for other sources of heterogeneity. This result suggests that not accounting for heterogeneity leads to flawed inference and/or that emigration heterogeneity mimics survival heterogeneity and biases senescence estimates.     

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.     

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.     

Can Camera Traps Monitor Komodo Dragons a Large Ectothermic Predator?

Camera trapping has greatly enhanced population monitoring of often cryptic and low abundance apex carnivores. Effectiveness of passive infrared camera trapping, and ultimately population monitoring, relies on temperature mediated differences between the animal and its ambient environment to ensure good camera detection. In ectothermic predators such as large varanid lizards, this criterion is presumed less certain. Here we evaluated the effectiveness of camera trapping to potentially monitor the population status of the Komodo dragon (Varanus komodoensis), an apex predator, using site occupancy approaches. We compared site-specific estimates of site occupancy and detection derived using camera traps and cage traps at 181 trapping locations established across six sites on four islands within Komodo National Park, Eastern Indonesia. Detection and site occupancy at each site were estimated using eight competing models that considered site-specific variation in occupancy (ψ)and varied detection probabilities (p) according to detection method, site and survey number using a single season site occupancy modelling approach. The most parsimonious model [ψ (site), p (site*survey); ω = 0.74] suggested that site occupancy estimates differed among sites. Detection probability varied as an interaction between site and survey number. Our results indicate that overall camera traps produced similar estimates of detection and site occupancy to cage traps, irrespective of being paired, or unpaired, with cage traps. Whilst one site showed some evidence detection was affected by trapping method detection was too low to produce an accurate occupancy estimate. Overall, as camera trapping is logistically more feasible it may provide, with further validation, an alternative method for evaluating long-term site occupancy patterns in Komodo dragons, and potentially other large reptiles, aiding conservation of this species.     

The importance of incorporating imperfect detection in biodiversity assessments: a case study of small mammals in an Australian region

Aim Assessments of biodiversity are an essential requirement of conservation management planning. Species distributional modelling is a popular approach to quantifying biodiversity whereby occurrence data are related to environmental covariates. An important confounding factor that is often overlooked in the development of such models is uncertainty due to imperfect detection. Here, I demonstrate how an analytical approach that accounts for the bias due to imperfect detection can be applied retrospectively to an existing biodiversity survey data set to produce more realistic estimates of species distributions and unbiased covariate relationships. Location Pilbara biogeographic region, Australia. Methods As a component of the Pilbara survey, presence/absence (i.e. undetected) data on small ground‐dwelling mammals were collected. I applied a multiseason occupancy modelling approach to six of the most common species encountered during this survey. Detection and occupancy rates, as well as extinction and colonization probabilities, were determined, and the influence of covariates on these parameters was examined using the multi‐model inference approach. Results Detection probabilities for all six species were considerably lower than 1.0 and varied across time and species. Naïve estimates of occupancy underestimated occupancy rates corrected for species detectability by up to 45%. Seasonal variation in occupancy status was attributed to changes in detection for two of the focal species, while reproductive events explained variation in occupancy in three others. Covariates describing the substrate strongly influenced site occupancy for most of the species modelled. A comparison of the occupancy model with a generalized linear model, assuming perfect detection, showed that the effects of the covariates were underestimated in the latter model. Main conclusions The application of the multiseason occupancy modelling approach to the Pilbara mammal data set demonstrated a powerful framework for examining changes in site occupancy, as well as local colonization and extinction rates of species which are not confounded by variable species detection rates.     

Monitoring occurrence,extinction, and colonization probabilities for gibbon populations

All gibbon species (Family: Hylobatidae) are considered threatened with extinction and recognized on the International Union for Conservation of Nature Red List of Threatened Species. Because gibbons are one of the most threatened families of primates, monitoring their status is now critically important. Long-term monitoring programs applying occupancy approaches, in addition to assessing occurrence probability, improves understanding of other population parameters such as site extinction or colonization probabilities, which elucidate temporal and spatial changes and are therefore important for guiding conservation efforts. In this study, we used multiple season occupancy models to monitor occurrence, extinction, and colonization probabilities for northern yellow-cheeked crested gibbon Nomascus annamensis in three adjacent protected areas in the Central Annamites mountain range, Vietnam. We collected data at 30 listening posts in 2012, 2014, and 2016 using the auditory point count method. Occurrence probabilities were highest in 2012 (0.74, confidence interval [CI]: 0.56–0.87) but slightly lower in 2014 (0.66, CI: 0.51–0.79) and 2016 (0.67, CI: 0.49–0.81). Extinction probabilities during the 2012–2014 and 2014–2016 intervals were 0.26 (0.14–0.44) and 0.25 (0.12–0.44), respectively. Colonization probabilities during 2012–2014 were 0.44 (0.19–0.73) and between 2014 and 2016 was 0.51 (0.26–0.75). Although local site extinctions have occurred, high recolonization probability helped to replenish the unoccupied sites and kept the occurrence probability stable. Long-term monitoring programs which use occurrence probability alone might not fully reveal the true dynamics of gibbon populations. We strongly recommend including multiple season occupancy models to monitor occurrence, extinction, and colonization probabilities in long-term gibbon monitoring programs.     

Dealing with detection error in site occupancy surveys: what can we do with a single survey?

Aim Site occupancy probabilities of target species are commonly used in various ecological studies, e.g. to monitor current status and trends in biodiversity. Detection error introduces bias in the estimators of site occupancy. Existing methods for estimating occupancy probability in the presence of detection error use replicate surveys. These methods assume population closure, i.e. the site occupancy status remains constant across surveys, and independence between surveys. We present an approach for estimating site occupancy probability in the presence of detection error that requires only a single survey and does not require assumption of population closure or independence. In place of the closure assumption, this method requires covariates that affect detection and occupancy.Methods Penalized maximum-likelihood method was used to estimate the parameters. Estimability of the parameters was checked using data cloning. Parametric boostrapping method was used for computing confidence intervals.Important findings The single-survey approach facilitates analysis of historical datasets where replicate surveys are unavailable, situations where replicate surveys are expensive to conduct and when the assumptions of closure or independence are not met. This method saves significant amounts of time, energy and money in ecological surveys without sacrificing statistical validity. Further, we show that occupancy and habitat suitability are not synonymous and suggest a method to estimate habitat suitability using single-survey data.     

A Variational Bayes Approach to the Analysis of Occupancy Models

Detection-nondetection data are often used to investigate species range dynamics using Bayesian occupancy models which rely on the use of Markov chain Monte Carlo (MCMC) methods to sample from the posterior distribution of the parameters of the model. In this article we develop two Variational Bayes (VB) approximations to the posterior distribution of the parameters of a single-season site occupancy model which uses logistic link functions to model the probability of species occurrence at sites and of species detection probabilities. This task is accomplished through the development of iterative algorithms that do not use MCMC methods. Simulations and small practical examples demonstrate the effectiveness of the proposed technique. We specifically show that (under certain circumstances) the variational distributions can provide accurate approximations to the true posterior distributions of the parameters of the model when the number of visits per site (K) are as low as three and that the accuracy of the approximations improves as K increases. We also show that the methodology can be used to obtain the posterior distribution of the predictive distribution of the proportion of sites occupied (PAO).     

Spatial patch occupancy patterns of the Lower Keys marsh rabbit

Reliable estimates of presence or absence of a species can provide substantial information on management questions related to distribution and habitat use but should incorporate the probability of detection to reduce bias. We surveyed for the endangered Lower Keys marsh rabbit (Sylvilagus palustris hefneri) in habitat patches on 5 Florida Key islands, USA, to estimate occupancy and detection probabilities. We derived detection probabilities using spatial replication of plots and evaluated hypotheses that patch location (coastal or interior) and patch size influence occupancy and detection. Results demonstrate that detection probability, given rabbits were present, was <0.5 and suggest that naïve estimates (i.e., estimates without consideration of imperfect detection) of patch occupancy are negatively biased. We found that patch size and location influenced probability of occupancy but not detection. Our findings will be used by Refuge managers to evaluate population trends of Lower Keys marsh rabbits from historical data and to guide management decisions for species recovery. The sampling and analytical methods we used may be useful for researchers and managers of other endangered lagomorphs and cryptic or fossorial animals occupying diverse habitats. © 2011 The Wildlife Society.     

Land use drivers of species re‐expansion: inferring colonization dynamics in Eurasian otters

Aim Land use intensity has been recognized as one of the major determinants of native species declines. The re‐expansion of species previously constrained by habitat degradation has been rarely investigated. Here, we use site occupancy models incorporating imperfect detection to identify the land use drivers of the re‐expansion of the Eurasian otter (Lutra lutra). Location Czech Republic. Methods We applied multi‐season occupancy models to otter presence–non‐detection data collected in three national surveys (1992, 2000, 2006) at 552 sites (11.2 × 12 km grid cells). Model parameters included site occupancy, colonization and extinction probabilities, and detection probability at a sub‐site level. We modelled changes in occupancy over time as a function of agricultural, urban and industrial land use and change in the extent of agricultural land use. Results Under the best fitting model, occupancy was estimated to be 34.6% in 1992, 51.3% in 2000 and 83.7% in 2006. Detection probability was neither perfect nor constant. Occupancy probability in 1992 was negatively related to land use gradients. Colonization was more likely to occur where a reduction in agricultural land was larger. Variation in extinction and colonization rates along land use gradients resulted in increased occupancy in industrial and especially urban landscapes. Conversely, occupancy remained almost unchanged along agricultural gradients. Main conclusions Dynamics of otter expansion were strongly associated with the two main patterns of the rapid environmental transition that has taken place in the Czech Republic since the early 1990s. Results show that a reduction in intensive agricultural land use led to an increase in otter distribution, providing evidence of the impact of agricultural land use on stream ecosystems. Moreover, otters recolonized urban and industrial landscapes, probably as a result of extensive reduction in water pollution from point sources. Our results suggest that active conservation of otter populations should focus on restoration of freshwater habitat at large scales, especially in agricultural landscapes.     

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.