The occupancy-abundance relationship and sampling designs using occupancy to monitor populations of Asian bears

Abstract

Designing a population monitoring program for Asian bears presents challenges associated with their low densities and detectability, generally large home ranges, and logistical or resource constraints. The use of an occupancy-based method to monitor bear populations can be appropriate under certain conditions given the mechanistic relationship between occupancy and abundance. The form of the occupancy–abundance relationship is dependent on species-specific characteristics such as home range size and population density, as well as study area size. To assess the statistical power of tests to detect population change of Asian bears, we conducted a study using a range of scenarios by simulating spatially explicit individual-based capture-recapture data from a demographically open model. Simulations assessed the power to detect changes in population density via changes in site-level occupancy or abundance through time, estimated using a standard occupancy model or a Royle-Nichols model, both with point detectors (representing camera traps). We used IUCN Red List criteria as a guide in selection of two population decline scenarios (20% and 50%), but we chose a shorter time horizon (10 years = 1 bear generation), meaning that declines were steeper than used for IUCN criteria (3 generations). Our simulations detected population declines of 50% with high power (>0.80) and low false positive rates (FPR: incorrectly detecting a decline) (<0.10) when detectors were spaced at > 0.67 times the home range diameter (home-range spacing ratio: HRSR, a measure of spatial correlation), such that bears would tend to overlap no more than two detectors. There was high (0.85) correlation between realized occupancy and N in these scenarios. The FPR increased as the HRSR decreased because of spatial correlation in the occupancy process induced when individual home ranges overlap multiple detectors. The mean statistical power to detect more gradual population declines (20% in 10 years) with HRSR > 0.67 was low for occupancy models 0.22 (maximum power 0.67) and Royle-Nichols models (0.24; maximum power 0.67), suggesting that declines of this magnitude may not be described reliably with 10 years of monitoring. Our results demonstrated that under many realistic scenarios that we explored, false positive rates were unacceptably high. We highlight that when designing occupancy studies, the spacing between point detectors be at least 0.67 times the diameter of the home range size of the larger sex (e.g., males) when the assumptions of the spatial capture-recapture model used for simulation are met.