Metapopulation dynamics in the butterfly Hipparchia semele changed decades before occupancy declined in The Netherlands

The survival of many species in human-dominated, fragmented landscapes depends on metapopulation dynamics, i.e., on a dynamic equilibrium of extinctions and colonizations in patches of suitable habitat. To understand and predict distributional changes, knowledge of these dynamics can be essential, and for this, metapopulation studies are preferably based on long-time-series data from many sites. Alas, such data are very scarce. An alternative is to use opportunistic data (i.e., collected without applying standardized field methods), but these data suffer from large variations in field methods and search intensity between sites and years. Dynamic site-occupancy models offer a general approach to adjust for variable survey effort. These models extend classical metapopulation models to account for imperfect detection of species and yield estimates of the probabilities of occupancy, colonization, and survival of species at sites. By accounting for detection, they fully correct for among-year variability in search effort. As an illustration, we fitted a dynamic site-occupancy model to 60 years of presence-absence data (more precisely, detection-nondetection) of the heathland butterfly Hipparchia semele in The Netherlands. Detection records were obtained from a database containing volunteer-based data from 1950-2009, and nondetection records were deduced from database records of other butterfly species. Our model revealed that metapopulation dynamics of Hipparchia had changed decades before the species' distribution began to contract. Colonization probability had already started to decline from 1950 onward, but this was counterbalanced by an increase in the survival of existing populations, the result of which was a stable distribution. Only from 1990 onward was survival not sufficient to compensate for the further decrease of colonization, and occupancy started to decline. Hence, it appears that factors acting many decades ago triggered a change in the metapopulation dynamics of this species, which ultimately led to a severe decline in occupancy that only became apparent much later. Our study emphasizes the importance of knowledge of changes in survival and colonization of species in modern landscapes over a very long time scale. It also demonstrates the power of site-occupancy modeling to obtain important population dynamics information from databases containing opportunistic sighting records.     

Sampling design trade-offs in occupancy studies with imperfect detection: examples and software.

Researchers have used occupancy, or probability of occupancy, as a response or state variable in a variety of studies (e.g., habitat modeling), and occupancy is increasingly favored by numerous state, federal, and international agencies engaged in monitoring programs. Recent advances in estimation methods have emphasized that reliable inferences can be made from these types of studies if detection and occupancy probabilities are simultaneously estimated. The need for temporal replication at sampled sites to estimate detection probability creates a trade-off between spatial replication (number of sample sites distributed within the area of interest/inference) and temporal replication (number of repeated surveys at each site). Here, we discuss a suite of questions commonly encountered during the design phase of occupancy studies, and we describe software (program GENPRES) developed to allow investigators to easily explore design trade-offs focused on particularities of their study system and sampling limitations. We illustrate the utility of program GENPRES using an amphibian example from Greater Yellowstone National Park, U.S.A.     

Neighborhood and habitat effects on vital rates: expansion of the Barred Owl in the Oregon coast ranges

In this paper, we modify dynamic occupancy models developed for detection-nondetection data to allow for the dependence of local vital rates on neighborhood occupancy, where neighborhood is defined very flexibly. Such dependence of occupancy dynamics on the status of a relevant neighborhood is pervasive, yet frequently ignored. Our framework permits joint inference about the importance of neighborhood effects and habitat covariates in determining colonization and extinction rates. Our specific motivation is the recent expansion of the Barred Owl (Strix varia) in western Oregon, USA, over the period 1990-2010. Because the focal period was one of dramatic range expansion and local population increase, the use of models that incorporate regional occupancy (sources of colonists) as determinants of dynamic rate parameters is especially appropriate. We began our analysis of 21 years of Barred Owl presence/nondetection data in the Tyee Density Study Area (TDSA) by testing a suite of six models that varied only in the covariates included in the modeling of detection probability. We then tested whether models that used regional occupancy as a covariate for colonization and extinction outperformed models with constant or year-specific colonization or extinction rates. Finally we tested whether habitat covariates improved the AIC of our models, focusing on which habitat covariates performed best, and whether the signs of habitat effects are consistent with a priori hypotheses. We conclude that all covariates used to model detection probability lead to improved AIC, that regional occupancy influences colonization and extinction rates, and that habitat plays an important role in determining extinction and colonization rates. As occupancy increases from low levels toward equilibrium, colonization increases and extinction decreases, presumably because there are more and more dispersing juveniles. While both rates are affected, colonization increases more than extinction decreases. Colonization is higher and extinction is lower in survey polygons with more riparian forest. The effects of riparian forest on extinction rates are greater than on colonization rates. Model results have implications for management of the invading Barred Owl, both through habitat alteration and removal.     

Use of Multispecies Occupancy Models to Evaluate the Response of Bird Communities to Forest Degradation Associated with Logging

Forest degradation is arguably the greatest threat to biodiversity, ecosystem services, and rural livelihoods. Therefore, increasing understanding of how organisms respond to degradation is essential for management and conservation planning. We were motivated by the need for rapid and practical analytical tools to assess the influence of management and degradation on biodiversity and system state in areas subject to rapid environmental change. We compared bird community composition and size in managed (ejido, i.e., communally owned lands) and unmanaged (national park) forests in the Sierra Tarahumara region, Mexico, using multispecies occupancy models and data from a 2‐year breeding bird survey. Unmanaged sites had on average higher species occupancy and richness than managed sites. Most species were present in low numbers as indicated by lower values of detection and occupancy associated with logging‐induced degradation. Less than 10% of species had occupancy probabilities >0.5, and degradation had no positive effects on occupancy. The estimated metacommunity size of 125 exceeded previous estimates for the region, and sites with mature trees and uneven‐aged forest stand characteristics contained the highest species richness. Higher estimation uncertainty and decreases in richness and occupancy for all species, including habitat generalists, were associated with degraded young, even‐aged stands. Our findings show that multispecies occupancy methods provide tractable measures of biodiversity and system state and valuable decision support for landholders and managers. These techniques can be used to rapidly address gaps in biodiversity information, threats to biodiversity, and vulnerabilities of species of interest on a landscape level, even in degraded or fast‐changing environments. Moreover, such tools may be particularly relevant in the assessment of species richness and distribution in a wide array of habitats. Uso de Modelos de Ocupación para Múltiples Especies para Evaluar la Respuesta de las Comunidades de Aves a la Degradación de Bosques Asociada con la Tala     

Fragmentation effects on an endangered species across a gradient from the interior to edge of its range

Understanding how habitat fragmentation affects individual species is complicated by challenges associated with quantifying species-specific habitat and spatial variability in fragmentation effects within a species’ range. We aggregated a 29-year breeding survey data set for the endangered marbled murrelet (Brachyramphus marmoratus) from >42,000 forest sites throughout the Pacific Northwest (Oregon, Washington, and northern California) of the United States. We built a species distribution model (SDM) in which occupied sites were linked with Landsat imagery to quantify murrelet-specific habitat and then used occupancy models to test the hypotheses that fragmentation negatively affects murrelet breeding distribution and that these effects are amplified with distance from the marine foraging habitat toward the edge of the species’ nesting range. Murrelet habitat declined in the Pacific Northwest by 20% since 1988, whereas the proportion of habitat comprising edges increased by 17%, indicating increased fragmentation. Furthermore, fragmentation of murrelet habitat at landscape scales (within 2 km of survey stations) negatively affected occupancy of potential breeding sites, and these effects were amplified near the range edge. On the coast, the odds of occupancy decreased by 37% (95% confidence interval [CI] –54 to 12) for each 10% increase in edge habitat (i.e., fragmentation), but at the range edge (88 km inland) these odds decreased by 99% (95% CI 98 to 99). Conversely, odds of murrelet occupancy increased by 31% (95% CI 14 to 52) for each 10% increase in local edge habitat (within 100 m of survey stations). Avoidance of fragmentation at broad scales but use of locally fragmented habitat with reduced quality may help explain the lack of murrelet population recovery. Further, our results emphasize that fragmentation effects can be nuanced, scale dependent, and geographically variable. Awareness of these nuances is critical for developing landscape-level conservation strategies for species experiencing broad-scale habitat loss and fragmentation.     

Transient dynamics of invasive competition: barred owls, spotted owls, habitat, and the demons of competition present

The recent range expansion of Barred Owls (Strix varia) into the Pacific Northwest, where the species now co-occurs with the endemic Northern Spotted Owl (Strix occidentalis caurina), resulted in a unique opportunity to investigate potential competition between two congeneric, previously allopatric species. The primary criticism of early competition research was the use of current species' distribution patterns to infer past processes; however, the recent expansion of the Barred Owl and the ability to model the processes that result in site occupancy (i.e., colonization and extinction) allowed us to address the competitive process directly rather than inferring past processes through current patterns. The purpose of our study was to determine whether Barred Owls had any negative effects on occupancy dynamics of nesting territories by Northern Spotted Owls and how these effects were influenced by habitat characteristics of Spotted Owl territories. We used single-species, multi-season occupancy models and covariates quantifying Barred Owl detections and habitat characteristics to model extinction and colonization rates of Spotted Owl pairs in southern Oregon, USA. We observed a strong, negative association between Barred Owl detections and colonization rates and a strong positive effect of Barred Owl detections on extinction rates of Spotted Owls. We observed increased extinction rates in response to decreased amounts of old forest at the territory core and higher colonization rates when old-forest habitat was less fragmented. Annual site occupancy for pairs reflected the strong effects of Barred Owls on occupancy dynamics with much lower occupancy rates predicted for territories where Barred Owls were detected. The strong Barred Owl and habitat effects on occupancy dynamics of Spotted Owls provided evidence of interference competition between the species. These effects increase the importance of conserving large amounts of contiguous, old-forest habitat to maintain Northern Spotted Owls in the landscape.     

Spatially Explicit Power Analyses for Occupancy‐Based Monitoring of Wolverine in the U.S. Rocky Mountains

Conservation scientists and resource managers often have to design monitoring programs for species that are rare or patchily distributed across large landscapes. Such programs are frequently expensive and seldom can be conducted by one entity. It is essential that a prospective power analysis be undertaken to ensure stated monitoring goals are feasible. We developed a spatially based simulation program that accounts for natural history, habitat use, and sampling scheme to investigate the power of monitoring protocols to detect trends in population abundance over time with occupancy‐based methods. We analyzed monitoring schemes with different sampling efforts for wolverine (Gulo gulo) populations in 2 areas of the U.S. Rocky Mountains. The relation between occupancy and abundance was nonlinear and depended on landscape, population size, and movement parameters. With current estimates for population size and detection probability in the northern U.S. Rockies, most sampling schemes were only able to detect large declines in abundance in the simulations (i.e., 50% decline over 10 years). For small populations reestablishing in the Southern Rockies, occupancy‐based methods had enough power to detect population trends only when populations were increasing dramatically (e.g., doubling or tripling in 10 years), regardless of sampling effort. In general, increasing the number of cells sampled or the per‐visit detection probability had a much greater effect on power than the number of visits conducted during a survey. Although our results are specific to wolverines, this approach could easily be adapted to other territorial species. Poder de Análisis Espacialmente Explícito para el Monitoreo Basado en Ocupación del Glotón (Gulo gulo) en las Montañas Rocallosas de Estados Unidos     

Estimating population impacts via dynamic occupancy analysis of Before-After Control-Impact studies

Estimating environmental impacts on populations is one of the main goals of wildlife monitoring programs, which are often conducted in conjunction with management actions or following natural disturbances. In this study we investigate the statistical power of dynamic occupancy models to detect changes in local survival and colonization from detection-nondetection data, while accounting for imperfect detection probability, in a Before-After Control-Impact (BACI) framework. We simulated impacts on local survival and/or detection probabilities, and asked questions related to: (1) costs and benefits of different analysis models, (2) confounding changes in detection with changes in local survival, (3) sampling design trade-offs, and (4) species with low vs. high rates of turnover. Estimating seasonal effects on local survival and colonization, as opposed to estimating Before-After effects, had little effect on the power to detect changes in local survival. Estimating a parameter that accounted for pretreatment differences in local survival between Control and Impact sites decreased power by 50%, but it was critical to include when such differences existed. When the experimental treatment had a negative impact on species detectability but analysis assumed constant detection, the Type I error rates were dramatically inflated (0.20 0.33). In general, there was low power (< 0.5) to detect a 50% decrease in local survival for all combinations of sites (N = 50 vs. 100), seasons sampled (8 vs. 12), and visits per site per season (4 vs. 6). Unbalanced designs performed worse than balanced designs, with the exception of the case of treatments being implemented in different seasons at different sites. Adding more control sites improved the ability to detect changes in local survival. Surveying more seasons after impact resulted in modest power gains, but at least three seasons before impact were required to successfully implement BACI occupancy studies. Turnover rates had a low impact on power. Occupancy studies conducted in a BACI design offer the opportunity to detect environmental impacts on wildlife populations without the costs of intensive studies. However, given the low power to detect small changes (20%) in local survival, these studies should be used when researchers are confident that major treatment impacts will occur or very large sample sizes are obtainable.     

Using experimental reintroductions to resolve the roles of habitat quality and metapopulation dynamics on patch occupancy in fragmented landscapes

Declines of species in fragmented landscapes can potentially be reversed either by restoring connectivity or restoring local habitat quality. Models fitted to snapshot occupancy data can be used to predict the effectiveness of these actions. However, such inferences can be misleading if the reliability of the habitat and landscape metrics used is unknown. The only way to unambiguously resolve the roles of habitat quality and metapopulation dynamics is to conduct experimental reintroductions to unoccupied patches so that habitat quality can be measured directly from data on vital rates. We, therefore, conducted a 15-year study that involved reintroducing a threatened New Zealand bird to unoccupied forest fragments to obtain reliable data on their habitat quality and reassess initial inferences made by modeling occupancy against habitat and landscape metrics. Although reproductive rates were similar among fragments, subtle differences in adult survival rates resulted in λ (finite rate of increase) estimations of <0.9 for 9 of the 12 fragments that were previously unoccupied. This was the case for only 1 of 14 naturally occupied fragments. This variation in λ largely explained the original occupancy pattern, reversing our original conclusion from occupancy modeling that this occupancy pattern was isolation driven and suggesting that it would be detrimental to increase connectivity without improving local habitat quality. These results illustrate that inferences from snapshot occupancy should be treated with caution and subjected to testing through experimental reintroductions in selected model systems.     

Landscape Characteristics of Fragmented Shrubsteppe Habitats and Breeding Passerine Birds

We examined the influence of local and landscape-level attributes of fragmented habitats in shrubsteppe habitats on the breeding distributions of Sage ( Amphispiza belli ) and Brewer's ( Spizella breweri ) Sparrows, Sage Thrashers ( Oreoscoptes montanus ) Horned Larks ( Eremophila alpestris ), and Western Meadowlarks ( Sturnella neglecta ) in the Snake River Plains of southwestern Idaho. We developed habitat (resource) selection models for each species by combining bird counts conducted from 1991 through 1933 with local vegetation characteristics and landscape attributes derived from satellite imagery. Site selection by shrubsteppe species (Sage and Brewer's Sparrows, and Sage Thrashers) depended on local vegetation cover and landscape features, such as the patch size of shrub habitats or the spatial similarity of sites. Marginal sites for these species (with species present in one of three years) were intermediate between unoccupied (never present) and occupied sites along environmental gradients characterized by increasing size of shrub habitat patches and total shrub cover and by decreasing disturbance. Horned Larks and Western Meadowlarks, typical grassland species, were not sensitive to landscape features, and their occupancy depended on the amount of grassland or shrub cover. In contrast to shrubsteppe species, sites that varied by occupancy rates of Western Meadowlarks did not significantly differ in vegetation or landscape components. Our results demonstrate that fragmentation of shrubsteppe significantly influenced the presence of shrub-obligate species. Because of restoration difficulties, the disturbance of semiarid shrubsteppe may cause irreversible loss of habitat and significant long-term consequences for the conservation of shrub-obligate birds.     

Developing a Monitoring Program for Invertebrates: Guidelines and a Case Study

Abstract:  Invertebrates provide the majority of ecosystem services; thus, it is important that they be inventoried, monitored, and protected. Nevertheless, inventories, monitoring, and management generally focus on vertebrates and flowering plants. Consequently, there are few guidelines or case studies for invertebrates. We present a procedure for developing a monitoring program for species-rich invertebrates that entails (1) characterizing the community; (2) identifying surrogates for biodiversity; and (3) establishing efficient methods to monitor surrogates and any ecologically important or sensitive taxa. We used these procedures, biodiversity-based statistical advances, and a survey of arthropods to develop a monitoring plan for the forests of Shenandoah National Park, Virginia (U.S.A.). Our case study revealed that mixed hardwood and hemlock forests had significantly different compositions of arthropods in their soil and understory strata. Of the 10 orders tested Coleoptera and Hymenoptera were the only two to pass most of the five surrogate tests, and their combination improved predictions of overall arthropod diversity. Because arthropods represent the majority of macroscopic species in most ecosystems, the ability of this assemblage to predict overall arthropod diversity makes it a powerful surrogate. Of the 11 collecting methods used, the beat-sheet method was the most efficient for monitoring this surrogate assemblage. To complement this coarse-filter approach to monitoring at-risk, invasive, or other important taxa (fine filter), we used ordination analyses to match 66 taxa with the methods that most effectively sampled them. Our methods serve as a model for developing an invertebrate monitoring plan and should facilitate linking such monitoring with ecosystem functions and management.     

A Multispecies Framework for Landscape Conservation Planning

Abstract: Rapidly changing landscapes have spurred the need for quantitative methods for conservation assessment and planning that encompass large spatial extents. We devised and tested a multispecies framework for conservation planning to complement single‐species assessments and ecosystem‐level approaches. Our framework consisted of 4 elements: sampling to effectively estimate population parameters, measuring how human activity affects landscapes at multiple scales, analyzing the relation between landscape characteristics and individual species occurrences, and evaluating and comparing the responses of multiple species to landscape modification. We applied the approach to a community of terrestrial birds across 25,000 km² with a range of intensities of human development. Human modification of land cover, road density, and other elements of the landscape, measured at multiple spatial extents, had large effects on occupancy of the 67 species studied. Forest composition within 1 km of points had a strong effect on occupancy of many species and a range of negative, intermediate, and positive associations. Road density within 1 km of points, percent evergreen forest within 300 m, and distance from patch edge were also strongly associated with occupancy for many species. We used the occupancy results to group species into 11 guilds that shared patterns of association with landscape characteristics. Our multispecies approach to conservation planning allowed us to quantify the trade‐offs of different scenarios of land‐cover change in terms of species occupancy.     

Changes in seed dispersal processes and the potential for between-patch connectivity for an arid land daisy

Dispersal is a major and critical process in population biology that has been particularly challenging to study. Animals can have major roles in seed dispersal even in species that do not appear specifically adapted to animal-aided dispersal. This can occur by two processes: direct movement of diaspores by animals and modification of landscape characteristics by animals in ways that greatly influence dispersal. We exploited the production of large, persistent dispersal structures (seed heads, henceforth) by Erodiophyllum elderi (Asteraceae), a daisy from arid Australia, to further understand secondary dispersal. Seed head dispersal on and off animal tracks in eight E. elderi patches was monitored for 9.5 months by periodically recording the location of marked seed heads. Sites were located inside a reserve that excludes sheep but not kangaroos, and in a nearby area with both kangaroos and sheep. The distance moved and likelihood of seed head movement was higher in areas with sheep, and especially along animal tracks. There was clear evidence that seed heads were channeled down animal tracks during large rainfall events. Seed head dispersal away from patches occurred to a limited extent via their physical contact with sheep and potentially via wind dispersal. Thus, the advantages of this study system allowed us to demonstrate the two postulated effects of herbivores on dispersal via direct movement of seed heads, and two distinct indirect effects through landscape modification by herbivores from the creation of animal tracks and the denudation of vegetation.     

Effectiveness of Vegetation Surrogates for Parasitoid Wasps in Reserve Selection

Abstract:  Selecting suitable nature reserves is a continuing challenge in conservation, particularly for target groups that are time-consuming to survey, species rich, and extinction prone. One such group is the parasitoid Hymenoptera, which have been excluded from conservation planning. If basic characteristics of habitats or vegetation could be used as reliable surrogates of specific target taxa, this would greatly facilitate appropriate reserve selection. We identified a range of potential habitat indicators of the species richness of pimpline parasitoid communities (Hymenoptera: Ichneumonidae: Pimplinae, Diacritinae, Poemeniinae) and tested their efficiency at capturing the observed diversity in a group of small woodlands in the agricultural landscape of the Vale of York (United Kingdom). Eight of the 18 vegetation-based reserve-selection strategies were significantly better at parasitoid species inclusion than random selection of areas. The best strategy maximized richness of tree species over the entire reserve network through complementarity. This strategy omitted only 2–3 species more (out of 38 captured in the landscape as a whole) than selections derived from the parasitoid survey data. In general, strategies worked equally well at capturing species richness and rarity. Our results suggest that vegetation data as a surrogate for species richness could prove an informative tool in parasitoid conservation, but further work is needed to test how broadly applicable these indicators may be.     

Community occupancy responses of small mammals to restoration treatments in ponderosa pine forests, northern Arizona, USA

In western North American conifer forests, wildfires are increasing in frequency and severity due to heavy fuel loads that have accumulated after a century of fire suppression. Forest restoration treatments (e.g., thinning and/or burning) are being designed and implemented at large spatial and temporal scales in an effort to reduce fire risk and restore forest structure and function. In ponderosa pine (Pinus ponderosa) forests, predominantly open forest structure and a frequent, low-severity fire regime constituted the evolutionary environment for wildlife that persisted for thousands of years. Small mammals are important in forest ecosystems as prey and in affecting primary production and decomposition. During 2006-2009, we trapped eight species of small mammals at 294 sites in northern Arizona and used occupancy modeling to determine community responses to thinning and habitat features. The most important covariates in predicting small mammal occupancy were understory vegetation cover, large snags, and treatment. Our analysis identified two generalist species found at relatively high occupancy rates across all sites, four open-forest species that responded positively to treatment, and two dense-forest species that responded negatively to treatment unless specific habitat features were retained. Our results indicate that all eight small mammal species can benefit from restoration treatments, particularly if aspects of their evolutionary environment (e.g., large trees, snags, woody debris) are restored. The occupancy modeling approach we used resulted in precise species-level estimates of occupancy in response to habitat attributes for a greater number of small mammal species than in other comparable studies. We recommend our approach for other studies faced with high variability and broad spatial and temporal scales in assessing impacts of treatments or habitat alteration on wildlife species. Moreover, since forest planning efforts are increasingly focusing on progressively larger treatment implementation, better and more efficiently obtained ecological information is needed to inform these efforts.     

Use of ambiguous detections to improve estimates from species distribution models

As large carnivores recover throughout Europe, their distribution needs to be studied to determine their conservation status and assess the potential for human-carnivore conflicts. However, efficient monitoring of many large carnivore species is challenging due to their rarity, elusive behavior, and large home ranges. Their monitoring can include opportunistic sightings from citizens in addition to designed surveys. Two types of detection errors may occur in such monitoring schemes: false negatives and false positives. False-negative detections can be accounted for in species distribution models (SDMs) that deal with imperfect detection. False-positive detections, due to species misidentification, have rarely been accounted for in SDMs. Generally, researchers use ad hoc data-filtering methods to discard ambiguous observations prior to analysis. These practices may discard valuable ecological information on the distribution of a species. We investigated the costs and benefits of including data types that may include false positives rather than discarding them for SDMs of large carnivores. We used a dynamic occupancy model that simultaneously accounts for false negatives and positives to jointly analyze data that included both unambiguous detections and ambiguous detections. We used simulations to compare the performances of our model with a model fitted on unambiguous data only. We tested the 2 models in 4 scenarios in which parameters that control false-positive detections and true detections varied. We applied our model to data from the monitoring of the Eurasian lynx (Lynx lynx) in the European Alps. The addition of ambiguous detections increased the precision of parameter estimates. For the Eurasian lynx, incorporating ambiguous detections produced more precise estimates of the ecological parameters and revealed additional occupied sites in areas where the species is likely expanding. Overall, we found that ambiguous data should be considered when studying the distribution of large carnivores through the use of dynamic occupancy models that account for misidentification.     

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.     

Using interview surveys and multispecies occupancy models to inform vertebrate conservation

Species distribution data are an essential biodiversity variable requiring robust monitoring to inform wildlife conservation. Yet, such data remain inherently sparse because of the logistical challenges of monitoring biodiversity across broad geographic extents. Surveys of people knowledgeable about the occurrence of wildlife provide an opportunity to evaluate species distributions and the ecology of wildlife communities across large spatial scales. We analyzed detection histories of 30 vertebrate species across the Western Ghats biodiversity hotspot in India, obtained from a large-scale interview survey of 2318 people who live and work in the forests of this region. We developed a multispecies occupancy model that simultaneously corrected for false-negative (non-detection) and false-positive (misidentification) errors that interview surveys can be prone to. Using this model, we integrated data across species in composite analyses of the responses of functional species groups (based on disturbance tolerance, diet, and body mass traits) to spatial variation in environmental variables, protection, and anthropogenic pressures. We observed a positive association between forest cover and the occurrence of species with low tolerance of human disturbance. Protected areas were associated with higher occurrence for species across different functional groups compared with unprotected lands. We also observed the occurrence of species with low disturbance tolerance, herbivores, and large-bodied species was negatively associated with developmental pressures, such as human settlements, energy production and mining, and demographic pressures, such as biological resource extraction. For the conservation of threatened vertebrates, our work underscores the importance of maintaining forest cover and reducing deforestation within and outside protected areas, respectively. In addition, mitigating a suite of pervasive human pressures is also crucial for wildlife conservation in one of the world's most densely populated biodiversity hotspots.     

A Framework for Conceptualizing Human Effects on Landscapes and Its Relevance to Management and Research Models

Abstract: The concept of habitat fragmentation is limited in its ability to describe the range of possible landscape configurations created by a variety of disturbances. This limitation is especially problematic in landscapes where human use of the habitat matrix occurs at multiple levels and where habitat modification may be a more important consideration than a simple binary classification of habitat versus nonhabitat. We propose a synthesizing scheme that places intact, variegated, fragmented, and relictual landscape states on a continuum, depending on the degree of habitat destruction. At a second level, the scheme considers the patterns of habitat modification that are imposed on remaining habitats. Management for conservation involves halting and sometimes reversing the trends of habitat destruction and modification. Conservation strategies will differ according to the state of alteration of the landscape, but all strategies include some consideration of the degree of modification of the matrix in determining habitat viability. It is convenient for biologists to assess landscape alteration state in terms of the persistence of large structural elements such as trees. Because animal species use habitats differently, however, they also experience the landscape differently. A landscape considered structurally fragmented by humans may be functionally variegated to other species. Therefore, it is necessary to consider the extent to which the entire landscape, including the matrix, is accessible and utilized by organisms with different spatial scales of resource use.     

Animal dispersal modelling: Handling landscape features and related animal choices

Animal dispersal in a fragmented landscape depends on the complex interaction between landscape structure and animal behavior. To better understand how individuals disperse, it is important to explicitly represent the properties of organisms and the landscape in which they move. A common approach to modelling dispersal includes representing the landscape as a grid of equal sized cells and then simulating individual movement as a correlated random walk. This approach uses a priori scale of resolution, which limits the representation of all landscape features and how different dispersal abilities are modelled.We develop a vector-based landscape model coupled with an object-oriented model for animal dispersal. In this spatially explicit dispersal model, landscape features are defined based on their geographic and thematic properties and dispersal is modelled through consideration of an organism's behavior, movement rules and searching strategies (such as visual cues). We present the model's underlying concepts, its ability to adequately represent landscape features and provide simulation of dispersal according to different dispersal abilities. We demonstrate the potential of the model by simulating two virtual species in a real Swiss landscape. This illustrates the model's ability to simulate complex dispersal processes and provides information about dispersal such as colonization probability and spatial distribution of the organism's path.