Estimating spatially variable and density-dependent survival using open-population spatial capture-recapture models

Abstract

Open-population spatial capture-recapture (OPSCR) models use the spatial information contained in individual detections collected over multiple consecutive occasions to estimate occasion-specific density, but also demographic parameters. OPSCR models can also estimate spatial variation in vital rates, but such models are neither widely used nor thoroughly tested. We developed a Bayesian OSPCR model that not only accounts for spatial variation in survival using spatial covariates, but also estimates local density-dependent effects on survival within a unified framework. Using simulations, we show that OPSCR models provide sound inferences on the effect of spatial covariates on survival, including multiple competing sources of mortality, each with potentially different spatial determinants. Estimation of local density-dependent survival was possible but required more data due to the greater complexity of the model. Not accounting for spatial heterogeneity in survival led to positive bias in abundance estimates (up to 10% relative bias). We provide a set of features in R package nimbleSCR that allow computationally efficient fitting of Bayesian OPSCR models with spatially varying survival. The ability to make population-level inferences of spatial variation in survival is an essential step towards a fully spatially-explicit OPSCR model that can disentangle the role of multiple spatial drivers on population dynamics. Open Research statementcode to reproduce the analysis is available on github; https://github.com/Cyril-Milleret/Public/tree/master/SpatialSurvivalOPSCR