An evaluation of multistate occupancy models for estimating relative abundance and population trends

Detecting spatiotemporal changes in the abundances of organisms is key to effectively conserving species. While indices of abundance have long been used, there has been a shift toward model-based estimators that account for the detection process. Popular approaches including traditional occupancy models and N-mixture models entail tradeoffs. The traditional occupancy approach requires the researcher coarsen the characterization of abundance to the probability that a site is occupied or unoccupied. Conversely, N-mixture models make use of variation in counts, but perform poorly when individuals have low detectability or move into or out of sites between visits. Multistate occupancy models that differentiate relatively abundant from non-abundant states have the potential to fill this gap but have been underexplored. We conducted a simulation study to test whether multistate occupancy models could capture spatial abundance patterns and detect population declines in the face of low individual detection probability (p ≤ 0.3) and unmodeled heterogeneity (e.g., that arising from individual movement). We considered 10,773 scenarios to examine the effects of differing amounts of heterogeneity as well as alternative study designs, population parameters, and modeling choices. We tracked bias in the proportion of sites estimated to be in the abundant state for single-season models, and power to detect a declining trend across multiple years. We also evaluated data diagnostic metrics to provide guidance to users. Multistate occupancy models were able to differentiate sites with higher abundances from sites with lower abundances when there were at least medium levels of spatial heterogeneity in true abundances. If different sites were randomly selected each year, power to detect even large population declines (65%) was poor (power < 0.8). However, if the same sites were surveyed each year, and a dynamic multistate occupancy was used, multistate occupancy models could detect (power ≥ 0.8) relatively small declines (5-40%) in 20% of scenarios, and frequently detect large declines of 45-60% (mean power = 0.92). Conservation decisions rely on detecting change reliably, rarely needing absolute abundance information. Multistate occupancy models can improve our ability to detect changing abundance while accommodating low individual detection probability and heterogeneity in count monitoring data.