Habitat Suitability Models and the Shortfall in Conservation Planning for African Vertebrates

Abstract:  Ongoing loss of biodiversity requires identifying large-scale conservation priorities, but the detailed information on the distribution of species required for this purpose is often missing. We present a systematic reserve selection for 1223 African mammals and amphibians in which habitat suitability models are used as estimates of the area occupied by species. In the framework of the World Conservation Union (IUCN) Global Amphibian Assessment and IUCN Global Mammal Assessment, we collected the geographic range (extent of occurrence) and habitat preferences for each species. We used the latter to build species-specific habitat suitability models inside geographic ranges, and for 181 species we verified the models by comparing suitability levels to presence-absence data collected in the field. We then used the suitable areas as estimators of the area of occupancy and compared the results of systematic reserve selection based on geographic ranges to those based on estimated areas of occupancy. Our results showed that the reserve system would need a 30-100% expansion to achieve minimal conservation targets, concentrated in the tropics, where species richness reaches a maximum. Comparative analyses revealed that using geographic ranges, which overestimate the area occupied by species, underestimates the total amount of area that needs to be conserved. The area selected for conservation doubled when we used the estimated area of occupancy in place of the geographic ranges. This happened because the suitable areas potentially occupied by each species overlapped less than their geographic ranges. As a result, any given protected area contained fewer species than predicted by the analysis of ranges. Because species are more specialized than our estimates of distribution based on extent of occurrence suggest, we propose that this is a general effect in systematic conservation planning.     

Integrating habitat-masked range maps with quantifications of prevalence to estimate area of occupancy in IUCN assessments

Estimates of species geographic ranges constitute critical input for biodiversity assessments, including those for the International Union for the Conservation of Nature (IUCN) Red List of Threatened Species. Area of occupancy (AOO) is one metric that IUCN uses to quantify a species’ range, but data limitations typically lead to either under- or overestimates (and unnecessarily wide bounds of uncertainty). Fortunately, existing methods in which range maps and land-cover data are used to estimate the area currently holding habitat for a species can be extended to yield an unbiased range of plausible estimates for AOO. Doing so requires estimating the proportion of sites (currently containing habitat) that a species occupies within its range (i.e., prevalence). Multiplying a quantification of habitat area by prevalence yields an estimate of what the species inhabits (i.e., AOO). For species with intense sampling at many sites, presence–absence data sets or occupancy modeling allow calculation of prevalence. For other species, primary biodiversity data (records of a species’ presence at a point in space and time) from citizen-science initiatives and research collections of natural history museums and herbaria could be used. In such cases, estimates of sample prevalence should be corrected by dividing by the species’ detectability. To estimate detectability from these data sources, extensions of inventory-completeness analyses merit development. With investments to increase the quality and availability of online biodiversity data, consideration of prevalence should lead to tighter and more realistic bounds of AOO for many taxonomic groups and geographic regions. By leading to more realistic and representative characterizations of biodiversity, integrating maps of current habitat with estimates of prevalence should empower conservation practitioners and decision makers and thus guide actions and policy worldwide.     

Predictors of occupancy trend across spatial scale

Many explorations of extinction probability have had a global focus, yet it is unclear whether variables that explain the probability of extinction at large spatial extents are the same as those at small spatial extents. Thus, we used nearly annual presence-absence records for the most recent 40 years of a 110-year data set from Palenque, Mexico, an area with ongoing deforestation, to explore which of >200 species of birds have probabilities of extirpation that are likely to increase. We assessed associations between long-term trends in species presence (i.e., detection in a given year) and body size, geographic range size, diet, dependence on forest cover, taxonomy, and ecological specialization. Our response variable was the estimated slope of a weighted logistic regression for each species. We assessed the relative strength of each predictor by means of a model ranking scheme. Several variables associated with high extinction probability at global extents, such as large body size or small geographic range size, were not associated with occurrence of birds over time at our site. Body size was associated with species loss at Palenque, but occurrence trends of both very large and very small species, particularly the latter, have declined, or the species have been extirpated. We found no association between declining occurrence trend and geographic range size, yet decline correlated with whether a species depends on forest (mean occupancy trend =-0.0380, 0.0263, and 0.0186 for, respectively, species with high, intermediate, or low dependence on forest) and with complex combinations of diet and foraging strata (e.g., occurrence of canopy insectivores and terrestrial omnivores has increased, whereas occurrence of mid-level frugivores and terrestrial granivores has decreased). Our findings emphasize that analyses of local areas are necessary to explicate extirpation risk at various spatial extents.     

The importance of incorporating functional habitats into conservation planning for highly mobile species in dynamic systems

The distribution of mobile species in dynamic systems can vary greatly over time and space. Estimating their population size and geographic range can be problematic and affect the accuracy of conservation assessments. Scarce data on mobile species and the resources they need can also limit the type of analytical approaches available to derive such estimates. We quantified change in availability and use of key ecological resources required for breeding for a critically endangered nomadic habitat specialist, the Swift Parrot (Lathamus discolor). We compared estimates of occupied habitat derived from dynamic presence‐background (i.e., presence‐only data) climatic models with estimates derived from dynamic occupancy models that included a direct measure of food availability. We then compared estimates that incorporate fine‐resolution spatial data on the availability of key ecological resources (i.e., functional habitats) with more common approaches that focus on broader climatic suitability or vegetation cover (due to the absence of fine‐resolution data). The occupancy models produced significantly (P < 0.001) smaller (up to an order of magnitude) and more spatially discrete estimates of the total occupied area than climate‐based models. The spatial location and extent of the total area occupied with the occupancy models was highly variable between years (131 and 1498 km²). Estimates accounting for the area of functional habitats were significantly smaller (2–58% [SD 16]) than estimates based only on the total area occupied. An increase or decrease in the area of one functional habitat (foraging or nesting) did not necessarily correspond to an increase or decrease in the other. Thus, an increase in the extent of occupied area may not equate to improved habitat quality or function. We argue these patterns are typical for mobile resource specialists but often go unnoticed because of limited data over relevant spatial and temporal scales and lack of spatial data on the availability of key resources. Understanding changes in the relative availability of functional habitats is crucial to informing conservation planning and accurately assessing extinction risk for mobile resource specialists.     

Geographic range size and extinction risk assessment in nomadic species

Geographic range size is often conceptualized as a fixed attribute of a species and treated as such for the purposes of quantification of extinction risk; species occupying smaller geographic ranges are assumed to have a higher risk of extinction, all else being equal. However many species are mobile, and their movements range from relatively predictable to‐and‐fro migrations to complex irregular movements shown by nomadic species. These movements can lead to substantial temporary expansion and contraction of geographic ranges, potentially to levels which may pose an extinction risk. By linking occurrence data with environmental conditions at the time of observations of nomadic species, we modeled the dynamic distributions of 43 arid‐zone nomadic bird species across the Australian continent for each month over 11 years and calculated minimum range size and extent of fluctuation in geographic range size from these models. There was enormous variability in predicted spatial distribution over time; 10 species varied in estimated geographic range size by more than an order of magnitude, and 2 species varied by >2 orders of magnitude. During times of poor environmental conditions, several species not currently classified as globally threatened contracted their ranges to very small areas, despite their normally large geographic range size. This finding raises questions about the adequacy of conventional assessments of extinction risk based on static geographic range size (e.g., IUCN Red Listing). Climate change is predicted to affect the pattern of resource fluctuations across much of the southern hemisphere, where nomadism is the dominant form of animal movement, so it is critical we begin to understand the consequences of this for accurate threat assessment of nomadic species. Our approach provides a tool for discovering spatial dynamics in highly mobile species and can be used to unlock valuable information for improved extinction risk assessment and conservation planning.     

Quantifying population declines based on presence‐only records for red‐list assessments

Accurate trend estimates are necessary for understanding which species are declining and which are most in need of conservation action. Imperfect species detection may result in unreliable trend estimates because this may lead to the overestimation of declines. Because many management decisions are based on population trend estimates, such biases could have severe consequences for conservation policy. We used an occupancy‐modeling framework to estimate detectability and calculate nationwide population trends for 14 Swiss amphibian species both accounting for and ignoring imperfect detection. Through the application of International Union for Conservation of Nature Red List criteria to the different trend estimates, we assessed whether ignoring imperfect detection could affect conservation policy. Imperfect detection occurred for all species and detection varied substantially among species, which led to the overestimation of population declines when detectability was ignored. Consequently, accounting for imperfect detection lowered the red‐list risk category for 5 of the 14 species assessed. We demonstrate that failing to consider species detectability can have serious consequences for species management and that occupancy modeling provides a flexible framework to account for observation bias and improve assessments of conservation status. A problem inherent to most historical records is that they contain presence‐only data from which only relative declines can be estimated. A move toward the routine recording of nonobservation and absence data is essential if conservation practitioners are to move beyond this toward accurate population trend estimation.     

Selecting pseudo-absence data for presence-only distribution modeling: How far should you stray from what you know?

An important decision in presence-only species distribution modeling is how to select background (or pseudo-absence) localities for model parameterization. The selection of such localities may influence model parameterization and thus, can influence the appropriateness and accuracy of the model prediction when extrapolating the species distribution across time and space. We used 12 species from the Australian Wet Tropics (AWT) to evaluate the relationship between the geographic extent from which pseudo-absences are taken and model performance, and shape and importance of predictor variables using the MAXENT modeling method. Model performance is lower when pseudo-absence points are taken from either a restricted or broad region with respect to species occurrence data than from an intermediate region. Furthermore, variable importance (i.e., contribution to the model) changed such that, models became increasingly simplified, dominated by just two variables, as the area from which pseudo-absence points were drawn increased. Our results suggest that it is important to consider the spatial extent from which pseudo-absence data are taken. We suggest species distribution modeling exercises should begin with exploratory analyses evaluating what extent might provide both the most accurate results and biologically meaningful fit between species occurrence and predictor variables. This is especially important when modeling across space or time—a growing application for species distributional modeling.     

Examining the relationship between local extinction risk and position in range

Over half of globally threatened animal species have experienced rapid geographic range loss. Identifying the parts of species’ distributions most vulnerable to local extinction would benefit conservation planning. However, previous studies give little consensus on whether ranges decline to the core or edge. We built on previous work by using empirical data to examine the position of recent local extinctions within species’ geographic ranges, address range position as a continuum, and explore the influence of environmental factors. We aggregated point‐locality data for 125 Galliform species from across the Palearctic and Indo‐Malaya into equal‐area half‐degree grid cells and used a multispecies dynamic Bayesian occupancy model to estimate rates of local extinctions. Our model provides a novel approach to identify loss of populations from within species ranges. We investigated the relationship between extinction rates and distance from range edge by examining whether patterns were consistent across biogeographic realm and different categories of land use. In the Palearctic, local extinctions occurred closer to the range edge than range core in both unconverted and human‐dominated landscapes. In Indo‐Malaya, no pattern was found for unconverted landscapes, but in human‐dominated landscapes extinctions tended to occur closer to the core than the edge. Our results suggest that local and regional factors override general spatial patterns of recent local extinction within species’ ranges and highlight the difficulty of predicting the parts of a species’ distribution most vulnerable to threat.     

Testing the abundant center model using range-wide demographic surveys of two coastal dune plants

It is widely accepted that species are most abundant at the center of their geographic ranges and become progressively rarer toward range limits. Although the abundant center model (ACM) has rarely been tested with range-wide surveys, it influences much thinking about the ecology and evolution of species' distributions. We tested ACM predictions using two unrelated but ecologically similar plants, Camissonia cheiranthifolia and Abronia umbellata. We intensively sampled both throughout their one-dimensional distributions within the Pacific coastal dunes of North America, from northern Baja California, Mexico, to southern Oregon, USA. Data from > 1100 herbarium specimens indicated that these limits have been stable for at least the last 100 years. Range-wide field surveys detected C. cheiranthifolia at 87% of 124 sites and A. umbellata at 54% of 113 sites, but site occupancy did not decline significantly toward range limits for either species. Permutation analysis did not detect a significant fit of geographical variation in local density to the ACM. Mean density did not correlate negatively with mean individual performance (plant size or number of seeds/plant), probably because both species occur at low densities. Although size and seeds per plant varied widely, central populations tended to have the highest values for size only. For C. cheiranthifolia, we observed asymmetry in the pattern of variation between the northern and southern halves of the range consistent with the long-standing prediction that range limits are imposed by different ecological factors in different parts of the geographical distribution. However, these asymmetries were difficult to interpret and likely reflect evolutionary differentiation as well as plastic responses to ecological variation. Both density and seeds per plant contributed to variation in seed production per unit area. In C. cheiranthifolia only, sites with highest seed production tended to occur at the range center, as predicted by the ACM and assumed by theory proposing that range limits evolve via antagonism between natural selection and gene flow.     

Assumption‐versus data‐based approaches to summarizing species’ ranges

For conservation decision making, species’ geographic distributions are mapped using various approaches. Some such efforts have downscaled versions of coarse‐resolution extent‐of‐occurrence maps to fine resolutions for conservation planning. We examined the quality of the extent‐of‐occurrence maps as range summaries and the utility of refining those maps into fine‐resolution distributional hypotheses. Extent‐of‐occurrence maps tend to be overly simple, omit many known and well‐documented populations, and likely frequently include many areas not holding populations. Refinement steps involve typological assumptions about habitat preferences and elevational ranges of species, which can introduce substantial error in estimates of species’ true areas of distribution. However, no model‐evaluation steps are taken to assess the predictive ability of these models, so model inaccuracies are not noticed. Whereas range summaries derived by these methods may be useful in coarse‐grained, global‐extent studies, their continued use in on‐the‐ground conservation applications at fine spatial resolutions is not advisable in light of reliance on assumptions, lack of real spatial resolution, and lack of testing. In contrast, data‐driven techniques that integrate primary data on biodiversity occurrence with remotely sensed data that summarize environmental dimensions (i.e., ecological niche modeling or species distribution modeling) offer data‐driven solutions based on a minimum of assumptions that can be evaluated and validated quantitatively to offer a well‐founded, widely accepted method for summarizing species’ distributional patterns for conservation applications.     

Habitat Loss and Changes in the Species-Area Relationship

Abstract: The species-area relationship (SAR) has been used successfully to predict extinction from extent of habitat reduction. These extinction estimates assume that species have uniformly distributed range requirements and a minimum abundance level required for persistence; how many species are lost depends solely on how much habitat is removed, not on where it is removed. We consider another limiting case in which range requirements, rather than abundances, determine extinctions. We used a new method for constructing SARs based on assumptions about geographic ranges of species. Our results show that habitat destruction can change the SAR and consequently the number of species predicted to be lost due to habitat destruction. Our method generates SARs that vary in shape according to the specific distributions of geographic range and occupancy but that have the common feature of being described by a power law with an exponent of <1. When the geographic range of species was included in the SAR, the way habitat was lost became important. Although the SAR before habitat destruction is often used to predict species loss after habitat destruction, assumptions must be clearly stated. To predict the damage caused by habitat loss with our model, it is necessary to know the fraction of aggregated species, the distribution of geographic ranges, the form of habitat destruction, and the sampling protocol. The remaining theoretical challenge is to develop a full theory that links abundance and range.     

Minimum area requirements for an at‐risk butterfly based on movement and demography

Determining the minimum area required to sustain populations has a long history in theoretical and conservation biology. Correlative approaches are often used to estimate minimum area requirements (MARs) based on relationships between area and the population size required for persistence or between species’ traits and distribution patterns across landscapes. Mechanistic approaches to estimating MAR facilitate prediction across space and time but are few. We used a mechanistic MAR model to determine the critical minimum patch size (CMP) for the Baltimore checkerspot butterfly (Euphydryas phaeton), a locally abundant species in decline along its southern range, and sister to several federally listed species. Our CMP is based on principles of diffusion, where individuals in smaller patches encounter edges and leave with higher probability than those in larger patches, potentially before reproducing. We estimated a CMP for the Baltimore checkerspot of 0.7–1.5 ha, in accordance with trait‐based MAR estimates. The diffusion rate on which we based this CMP was broadly similar when estimated at the landscape scale (comparing flight path vs. capture‐mark‐recapture data), and the estimated population growth rate was consistent with observed site trends. Our mechanistic approach to estimating MAR is appropriate for species whose movement follows a correlated random walk and may be useful where landscape‐scale distributions are difficult to assess, but demographic and movement data are obtainable from a single site or the literature. Just as simple estimates of lambda are often used to assess population viability, the principles of diffusion and CMP could provide a starting place for estimating MAR for conservation.     

Extinction Rate Estimates for Plant Populations in Revisitation Studies: Importance of Detectability

Abstract:  Many researchers have obtained extinction-rate estimates for plant populations by comparing historical and current records of occurrence. A population that is no longer found is assumed to have gone extinct. Extinction can then be related to characteristics of these populations, such as habitat type, size, or species, to test ideas about what factors may affect extinction. Such studies neglect the fact that a population may be overlooked, however, which may bias estimates of extinction rates upward. In addition, if populations are unequally detectable across groups to be compared, such as habitat type or population size, comparisons become distorted to an unknown degree. To illustrate the problem, I simulated two data sets, assuming a constant extinction rate, in which populations occurred in different habitats or habitats of different size and these factors affected their detectability. The conventional analysis implicitly assumed that detectability equalled 1 and used logistic regression to estimate extinction rates. It wrongly identified habitat and population size as factors affecting extinction risk. In contrast, with capture-recapture methods, unbiased estimates of extinction rates were recovered. I argue that capture-recapture methods should be considered more often in estimations of demographic parameters in plant populations and communities.     

Responses of pond-breeding amphibians to wildfire: short-term patterns in occupancy and colonization.

Wildland fires are expected to become more frequent and severe in many ecosystems, potentially posing a threat to many sensitive species. We evaluated the effects of a large, stand-replacement wildfire on three species of pond-breeding amphibians by estimating changes in occupancy of breeding sites during the three years before and after the fire burned 42 of 83 previously surveyed wetlands. Annual occupancy and colonization for each species was estimated using recently developed models that incorporate detection probabilities to provide unbiased parameter estimates. We did not find negative effects of the fire on the occupancy or colonization rates of the long-toed salamander (Ambystoma macrodactylum). Instead, its occupancy was higher across the study area after the fire, possibly in response to a large snowpack that may have facilitated colonization of unoccupied wetlands. Naive data (uncorrected for detection probability) for the Columbia spotted frog (Rana luteiventris) initially led to the conclusion of increased occupancy and colonization in wetlands that burned. After accounting for temporal and spatial variation in detection probabilities, however, it was evident that these parameters were relatively stable in both areas before and after the fire. We found a similar discrepancy between naive and estimated occupancy of A. macrodactylum that resulted from different detection probabilities in burned and control wetlands. The boreal toad (Bufo boreas) was not found breeding in the area prior to the fire but colonized several wetlands the year after they burned. Occupancy by B. boreas then declined during years 2 and 3 following the fire. Our study suggests that the amphibian populations we studied are resistant to wildfire and that B. boreas may experience short-term benefits from wildfire. Our data also illustrate how naive presence-non-detection data can provide misleading results.     

Effective Population Size in Red-Cockaded Woodpeckers: Population and Model Differences

Loss of genetic variability in isolated populations is an important issue for conservation biology. Most studies involve only a single population of a given species and a single method of estimating rate of loss. Here we present analyses for three different Red-cockaded Woodpecker ( Picoides borealis ) populations from different geographic regions. We compare two different models for estimating the expected rate of loss of genetic variability, and test their sensitivity to model parameters. We found that the simpler model (Reed et al. 1988) consistently estimated a greater rate of loss of genetic variability from a population than did the Emigh and Pollak (1979) model. The ratio of effective population size (which describes the expected rate of loss of genetic variability) to breeder population size varied widely among Red-cockaded Woodpecker populations due to geographic variation in demography. For this species, estimates of effective size were extremely sensitive to survival parameters, but not to the probability of breeding or reproductive success. Sensitivity was sufficient that error in estimating survival rates in the field could easily mask true population differences in effective size. Our results indicate that accurate and precise demographic data are prerequisites to determining effective population size for this species using genetic models, and that a single estimate of rate of loss of genetic variability is not valid across populations.     

The Quantitative Incidence Function Model and Persistence of an Endangered Butterfly Metapopulation

The incidence function model is derived from a linear first-order Markov chain of the presence or absence of a species in a habitat patch. The model can be parameterized with "snapshot" presence/absence data from a patch network. Using the estimated parameter values the Markov chain can be iterated in the same or in some other patch network to generate quantitative predictions about transient metapopulation dynamics and the stochastic steady state. We tested the ability of the incidence function model to predict patch occupancy using extensive data on an endangered butterfly, the Glanville fritillary ( Melitaea cinxia ) Parameter values were estimated with data collected from a 50-patch network in 1991. In 1993 we surveyed the entire geographic range of the species in Finland, within an area of 50 × 70 km², with 1502 habitat patches (dry meadows) of which 536 were occupied. Model predictions were generated for the 1502 patches and were compared with the observed pattern of occupancy in 1993. The model predicted patch occupancy well in more than half of the study area, but prediction was poor for one quarter of the area, probably because of regional variation in habitat quality and because metapopulations may have been perturbed away from the steady state. The incidence function model provides a practical tool for making quantitative predictions about metapopulation dynamics of species living in fragmented landscapes.     

Scaling range sizes to threats for robust predictions of risks to biodiversity

Assessments of risk to biodiversity often rely on spatial distributions of species and ecosystems. Range‐size metrics used extensively in these assessments, such as area of occupancy (AOO), are sensitive to measurement scale, prompting proposals to measure them at finer scales or at different scales based on the shape of the distribution or ecological characteristics of the biota. Despite its dominant role in red‐list assessments for decades, appropriate spatial scales of AOO for predicting risks of species’ extinction or ecosystem collapse remain untested and contentious. There are no quantitative evaluations of the scale‐sensitivity of AOO as a predictor of risks, the relationship between optimal AOO scale and threat scale, or the effect of grid uncertainty. We used stochastic simulation models to explore risks to ecosystems and species with clustered, dispersed, and linear distribution patterns subject to regimes of threat events with different frequency and spatial extent. Area of occupancy was an accurate predictor of risk (0.81<|r|<0.98) and performed optimally when measured with grid cells 0.1–1.0 times the largest plausible area threatened by an event. Contrary to previous assertions, estimates of AOO at these relatively coarse scales were better predictors of risk than finer‐scale estimates of AOO (e.g., when measurement cells are <1% of the area of the largest threat). The optimal scale depended on the spatial scales of threats more than the shape or size of biotic distributions. Although we found appreciable potential for grid‐measurement errors, current IUCN guidelines for estimating AOO neutralize geometric uncertainty and incorporate effective scaling procedures for assessing risks posed by landscape‐scale threats to species and ecosystems.     

Site‐Occupancy Distribution Modeling to Correct Population‐Trend Estimates Derived from Opportunistic Observations

Abstract: Species’ assessments must frequently be derived from opportunistic observations made by volunteers (i.e., citizen scientists). Interpretation of the resulting data to estimate population trends is plagued with problems, including teasing apart genuine population trends from variations in observation effort. We devised a way to correct for annual variation in effort when estimating trends in occupancy (species distribution) from faunal or floral databases of opportunistic observations. First, for all surveyed sites, detection histories (i.e., strings of detection–nondetection records) are generated. Within‐season replicate surveys provide information on the detectability of an occupied site. Detectability directly represents observation effort; hence, estimating detectablity means correcting for observation effort. Second, site‐occupancy models are applied directly to the detection‐history data set (i.e., without aggregation by site and year) to estimate detectability and species distribution (occupancy, i.e., the true proportion of sites where a species occurs). Site‐occupancy models also provide unbiased estimators of components of distributional change (i.e., colonization and extinction rates). We illustrate our method with data from a large citizen‐science project in Switzerland in which field ornithologists record opportunistic observations. We analyzed data collected on four species: the widespread Kingfisher (Alcedo atthis) and Sparrowhawk (Accipiter nisus) and the scarce Rock Thrush (Monticola saxatilis) and Wallcreeper (Tichodroma muraria). Our method requires that all observed species are recorded. Detectability was <1 and varied over the years. Simulations suggested some robustness, but we advocate recording complete species lists (checklists), rather than recording individual records of single species. The representation of observation effort with its effect on detectability provides a solution to the problem of differences in effort encountered when extracting trend information from haphazard observations. We expect our method is widely applicable for global biodiversity monitoring and modeling of species distributions.     

An evaluation of three statistical methods used to model resource selection

The performance of statistical methods for modeling resource selection by animals is difficult to evaluate with field data because true selection patterns are unknown. Simulated data based on a known probability distribution, though, can be used to evaluate statistical methods. Models should estimate true selection patterns if they are to be useful in analyzing and interpreting field data. We used simulation techniques to evaluate the effectiveness of three statistical methods used in modeling resource selection. We generated 25 use locations per animal and included 10, 20, 40, or 80 animals in samples of use locations. To simulate species of different mobility, we generated use locations at four levels according to a known probability distribution across DeSoto National Wildlife Refuge (DNWR) in eastern Nebraska and western Iowa, USA. We either generated 5 random locations per use location or 10,000 random locations (total) within 4 predetermined areas around use locations to determine how the definition of availability and the number of random locations affected results. We analyzed simulated data using discrete choice, logistic-regression, and a maximum entropy method (Maxent). We used a simple linear regression of estimated and known probability distributions and area under receiver operating characteristic curves (AUC) to evaluate the performance of each method. Each statistical method was affected differently by number of animals and random locations used in analyses, level at which selection of resources occurred, and area considered available. Discrete-choice modeling resulted in precise and accurate estimates of the true probability distribution when the area in which use locations were generated was ≥ the area defined to be available. Logistic-regression models were unbiased and precise when the area in which use locations were generated and the area defined to be available were the same size; the fit of these models improved with increased numbers of random locations. Maxent resulted in unbiased and precise estimates of the known probability distribution when the area in which use locations were generated was small (home-range level) and the area defined to be available was large (study area). Based on AUC analyses, all models estimated the selection distribution better than random chance. Results from AUC analyses, however, often contradicted results of the linear regression method used to evaluate model performance. Discrete-choice modeling was best able to estimate the known selection distribution in our study area regardless of sample size or number of random locations used in the analyses, but we recommend further studies using simulated data over different landscapes and different resource metrics to confirm our results. Our study offers an approach and guidance for others interested in assessing the utility of techniques for modeling resource selection in their study area.     

A Bayesian state-space formulation of dynamic occupancy models

Species occurrence and its dynamic components, extinction and colonization probabilities, are focal quantities in biogeography and metapopulation biology, and for species conservation assessments. It has been increasingly appreciated that these parameters must be estimated separately from detection probability to avoid the biases induced by non-detection error. Hence, there is now considerable theoretical and practical interest in dynamic occupancy models that contain explicit representations of metapopulation dynamics such as extinction, colonization, and turnover as well as growth rates. We describe a hierarchical parameterization of these models that is analogous to the state-space formulation of models in time series, where the model is represented by two components, one for the partially observable occupancy process and another for the observations conditional on that process. This parameterization naturally allows estimation of all parameters of the conventional approach to occupancy models, but in addition, yields great flexibility and extensibility, e.g., to modeling heterogeneity or latent structure in model parameters. We also highlight the important distinction between population and finite sample inference; the latter yields much more precise estimates for the particular sample at hand. Finite sample estimates can easily be obtained using the state-space representation of the model but are difficult to obtain under the conventional approach of likelihood-based estimation. We use R and WinBUGS to apply the model to two examples. In a standard analysis for the European Crossbill in a large Swiss monitoring program, we fit a model with year-specific parameters. Estimates of the dynamic parameters varied greatly among years, highlighting the irruptive population dynamics of that species. In the second example, we analyze route occupancy of Cerulean Warblers in the North American Breeding Bird Survey (BBS) using a model allowing for site-specific heterogeneity in model parameters. The results indicate relatively low turnover and a stable distribution of Cerulean Warblers which is in contrast to analyses of counts of individuals from the same survey that indicate important declines. This discrepancy illustrates the inertia in occupancy relative to actual abundance. Furthermore, the model reveals a declining patch survival probability, and increasing turnover, toward the edge of the range of the species, which is consistent with metapopulation perspectives on the genesis of range edges. Given detection/non-detection data, dynamic occupancy models as described here have considerable potential for the study of distributions and range dynamics.