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.

Modelling nonlinear responses of a desert rodent species to environmental change with hierarchical dynamic generalized additive models

Modeling Community Dynamics Through Environmental Effects, Species Interactions and Movement

Improved state-space models for inference about spatial and temporal variation in abundance from count data

Additive Models with Spatio-Temporal Data

A sequential Monte Carlo algorithm for data assimilation problems in ecology

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

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

Detecting and modelling delayed density-dependence in abundance time series of a small mammal (Didelphis aurita)

Ecological non-linear state space model selection via adaptive particle Markov chain Monte Carlo (AdPMCMC)

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

Hierarchical generalized additive models in ecology: an introduction with mgcv

Integrating counts from rigorous surveys and participatory science to better understand spatiotemporal variation in population processes

A hierarchical model for jointly assessing ecological and anthropogenic impacts on animal demography

Bayesian hierarchical space-time models to improve multispecies assessment by combining observations from disparate fish surveys

Predicting fine-scale distributions and emergent spatiotemporal patterns from temporally dynamic step selection simulations

Time‐varying habitat selection analysis: A model and applications for studying diel, seasonal, and post‐release changes

Addressing temporal trends in survivorship from cross-sectional sampling designs: A modelling framework with applications for megafauna conservation

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

Dynamic structural equation models synthesize ecosystem dynamics constrained by ecological mechanisms

Modelling group dynamic animal movement

Balancing structural complexity with ecological insight in Spatio‐temporal species distribution models