Abstract:Our integrated modelling method allows for managers to use citizen sightings combined with detection/non‐detection surveys to estimate robust indices of abundance for both high‐ and low‐density, and wide‐spread versus patchily distributed species. Through comparison with previous studies, we highlight how broad‐scale programmes, such as the statewide efforts to estimate species distributions undertaken here, can benefit substantively from integrated models that leverage additional data (here, incidental records) from a larger region of space, and thus capture more landscape heterogeneity than is plausible within formalized surveys alone. Estimating species distribution and abundance is foundational to effective management and conservation. Using an integrated species distribution model that combines presence‐only data from various sources with detection/non‐detection data from structured surveys, we estimated the distribution and expected abundance of three difficult‐to‐monitor mammals of management concern across New York State, namely, coyotes (Canis latrans), bobcats (Lynx rufus) and black bears (Ursus americanus). Three distinct landscape‐scale camera trap surveys provided detection/non‐detection data over 9 years between 2013 and 2021, and we augmented those data with incidental records of our focal species from public repositories. We used an inhomogeneous Poisson point process to construct an integrated model that fit both data types simultaneously. We demonstrate a simple application of spatial point density of all species records in the accessed public databases to inform the thinning process to account for unknown spatial sampling in the presence‐only data, often referred to as the 'magic covariate'. Using this approach, we examine habitat associations and provide spatially explicit estimates in expected abundance across the entirety of New York State for all three focal species. As expected, coyotes were the most widely distributed and abundant species, with a strong positive association with agricultural land uses. Bobcats exhibited low expected abundance throughout the state and showed positive associations with deciduous forest and forest edge, and a negative association with road density. Finally, we observed considerable spatial variation in abundance of black bears with expected abundance increasing in association with various forest cover and composition covariates and decreasing with crop cover. We present insights into habitat associations and spatial variation in abundance, and provide management implications for each of the species of interest. Synthesis and applications. Our integrated modelling method allows for managers to use citizen sightings combined with detection/non‐detection surveys to estimate robust indices of abundance for both high‐ and low‐density, and wide‐spread versus patchily distributed species. Through comparison with previous studies, we highlight how broad‐scale programmes, such as the statewide efforts to estimate species distributions undertaken here, can benefit substantively from integrated models that leverage additional data (here, incidental records) from a larger region of space, and thus capture more landscape heterogeneity than is plausible within formalized surveys alone.

Biogeographic multi-species occupancy models for large-scale survey data

Modeling Complex Species-Environment Relationships Through Spatially-Varying Coefficient Occupancy Models

Scaling multi-species occupancy models to large citizen science datasets

Joint species distribution models with imperfect detection for high-dimensional spatial data

A dynamic multi‐scale occupancy model to estimate temporal dynamics and hierarchical habitat use for nomadic species

spOccupancy: An R package for single‐species, multi‐species, and integrated spatial occupancy models

Predictive multi‐scale occupancy models at range‐wide extents: Effects of habitat and human disturbance on distributions of wetland birds

How do ecologists estimate occupancy in practice?

Spatiotemporal modelling of abundance from multiple data sources in an integrated spatial distribution model

A multispecies hierarchical model to integrate count and distance‐sampling data

A comprehensive multi-state conditional occupancy model for evaluating interactions of non-native and native species

Fast Bayesian inference for large occupancy data sets, using the Polya-Gamma scheme

Should we exploit opportunistic databases with joint species distribution models? Artificial and real data suggest it depends on the sampling completeness

Mapping multiscale breeding bird species distributions across the United States and evaluating their conservation applications

Spatial Clustering of Citizen Science Data Improves Downstream Species Distribution Models

Spatial occupancy models for data collected on stream networks

With super SDMs (machine learning, open access big data, and the cloud) towards more holistic global squirrel hotspots and coldspots

Semiparametric Bivariate Zero-Inflated Poisson Models with Application to Studies of Abundance for Multiple Species

A super SDM (species distribution model) 'in the cloud' for better habitat-association inference with a 'big data' application of the Great Gray Owl for Alaska

Combining Observational Data and Language for Species Range Estimation

Integrating presence‐only and detection/non‐detection data to estimate distributions and expected abundance of difficult‐to‐monitor species on a landscape‐scale