Counting the Capital's cats: Estimating drivers of abundance of free-roaming cats with a novel hierarchical model

Free-roaming cats are a conservation concern in many areas but identifying their impacts and developing mitigation strategies requires a robust understanding of their distribution and density patterns. Urban and residential areas may be especially relevant in this process because free-roaming cats are abundant in these anthropogenic landscapes. Here, we estimate the occupancy and density of free-roaming cats in Washington D.C. and relate these metrics to known landscape and social factors. We conducted an extended camera trap survey of public and private spaces across D.C. and analyzed data collected from 1483 camera deployments from 2018 to 2020. We estimated citywide cat distribution by fitting hierarchical occupancy models and further estimated cat abundance using a novel random thinning spatial capture-recapture model that allows for the use of photos that can and cannot be identified to individual. Within this model, we utilized individual covariates that provided identity exclusions between photos of unidentifiable cats with inconsistent coat patterns, thus increasing the precision of abundance estimates. This combined model also allowed for unbiased estimation of density when animals cannot be identified to individual at the same rate as for free-roaming cats whose identifiability depended on their coat characteristics. Cat occupancy and abundance declined with increasing distance from residential areas, an effect that was more pronounced in wealthier neighborhoods. There was noteworthy absence of cats detected in larger public spaces and forests. Realized densities ranged from 0.02 to 1.75 cats/ha in sampled areas, resulting in a district-wide estimate of ~7296 free-roaming cats. Ninety percent of cat detections lacked collars and nearly 35% of known individuals were ear-tipped, indicative of district Trap-Neuter-Return (TNR) programs. These results suggest that we mainly sampled and estimated the unowned cat subpopulation, such that indoor/outdoor housecats were not well represented. The precise estimation of cat population densities is difficult due to the varied behavior of subpopulations within free-roaming cat populations (housecats, stray and feral cats), but our methods provide a first step in establishing citywide baselines to inform data-driven management plans for free-roaming cats in urban environments.

A spatial capture-recapture model for group-living species

Group living in species can have complex consequences for individuals, populations, and ecosystems. Therefore, estimating group density and size is often essential for understanding population dynamics, interspecific interactions, and conservation needs of group-living species. Spatial capture-recapture (SCR) has been used to model both individual and group density in group-living species, but modeling either individual-level or group-level detection results in different biases due to common characteristics of group-living species, such as highly cohesive movement or variation in group size. Furthermore, no SCR method currently estimates group density, individual density, and group size jointly. Using clustered point processes, we developed a cluster SCR model to estimate group density, individual density, and group size. We compared the model to standard SCR models using both a simulation study and a data set of detections of African wild dogs (Lycaon pictus), a group-living carnivore, on camera traps in northern Botswana. We then tested the model's performance under various scenarios of group movement in a separate simulation study. We found that the cluster SCR model outperformed a standard group-level SCR model when fitted to data generated with varying group sizes, and mostly recovered previous estimates of wild dog group density, individual density, and group size. We also found that the cluster SCR model performs better as individuals' movements become more correlated with their groups' movements. The cluster SCR model offers opportunities to investigate ecological hypotheses relating group size to population dynamics while accounting for cohesive movement behaviors in group-living species.

Spatial capture-recapture with random thinning for unidentified encounters

Spatial capture-recapture (SCR) models have increasingly been used as a basis for combining capture-recapture data types with variable levels of individual identity information to estimate population density and other demographic parameters. Recent examples are the unmarked SCR (or spatial count model), where no individual identities are available and spatial mark-resight (SMR) where individual identities are available for only a marked subset of the population. Currently lacking, though, is a model that allows unidentified samples to be combined with identified samples when there are no separate classes of "marked" and "unmarked" individuals and when the two sample types cannot be considered as arising from two independent observation models. This is a common scenario when using noninvasive sampling methods, for example, when analyzing data on identified and unidentified photographs or scats from the same sites.Here we describe a "random thinning" SCR model that utilizes encounters of both known and unknown identity samples using a natural mechanistic dependence between samples arising from a single observation model. Our model was fitted in a Bayesian framework using NIMBLE.We investigate the improvement in parameter estimates by including the unknown identity samples, which was notable (up to 79% more precise) in low-density populations with a low rate of identified encounters. We then applied the random thinning SCR model to a noninvasive genetic sampling study of brown bear (Ursus arctos) density in Oriental Cantabrian Mountains (North Spain).Our model can improve density estimation for noninvasive sampling studies for low-density populations with low rates of individual identification, by making use of available data that might otherwise be discarded.

Improving estimation of puma (Puma concolor) population density: clustered camera-trapping, telemetry data, and generalized spatial mark-resight models

Obtaining reliable population density estimates for pumas (Puma concolor) and other cryptic, wide-ranging large carnivores is challenging. Recent advancements in spatially explicit capture-recapture models have facilitated development of novel survey approaches, such as clustered sampling designs, which can provide reliable density estimation for expansive areas with reduced effort. We applied clustered sampling to camera-traps to detect marked (collared) and unmarked pumas, and used generalized spatial mark-resight (SMR) models to estimate puma population density across 15,314 km² in the southwestern USA. Generalized SMR models outperformed conventional SMR models. Integrating telemetry data from collars on marked pumas with detection data from camera-traps substantially improved density estimates by informing cryptic activity (home range) center transiency and improving estimation of the SMR home range parameter. Modeling sex of unmarked pumas as a partially identifying categorical covariate further improved estimates. Our density estimates (0.84–1.65 puma/100 km²) were generally more precise (CV = 0.24–0.31) than spatially explicit estimates produced from other puma sampling methods, including biopsy darting, scat detection dogs, and regular camera-trapping. This study provides an illustrative example of the effectiveness and flexibility of our combined sampling and analytical approach for reliably estimating density of pumas and other wildlife across geographically expansive areas.

Consequences of severe habitat fragmentation on density, genetics, and spatial capture-recapture analysis of a small bear population

Loss and fragmentation of natural habitats caused by human land uses have subdivided several formerly contiguous large carnivore populations into multiple small and often isolated subpopulations, which can reduce genetic variation and lead to precipitous population declines. Substantial habitat loss and fragmentation from urban development and agriculture expansion relegated the Highlands-Glades subpopulation (HGS) of Florida, USA, black bears (Ursus americanus floridanus) to prolonged isolation; increasing human land development is projected to cause ≥ 50% loss of remaining natural habitats occupied by the HGS in coming decades. We conducted a noninvasive genetic spatial capture-recapture study to quantitatively describe the degree of contemporary habitat fragmentation and investigate the consequences of habitat fragmentation on population density and genetics of the HGS. Remaining natural habitats sustaining the HGS were significantly more fragmented and patchier than those supporting Florida’s largest black bear subpopulation. Genetic diversity was low (A_R = 3.57; H_E = 0.49) and effective population size was small (N_E = 25 bears), both of which remained unchanged over a period spanning one bear generation despite evidence of some immigration. Subpopulation density (0.054 bear/km²) was among the lowest reported for black bears, was significantly female-biased, and corresponded to a subpopulation size of 98 bears in available habitat. Conserving remaining natural habitats in the area occupied by the small, genetically depauperate HGS, possibly through conservation easements and government land acquisition, is likely the most important immediate step to ensuring continued persistence of bears in this area. Our study also provides evidence that preferentially placing detectors (e.g., hair traps or cameras) primarily in quality habitat across fragmented landscapes poses a challenge to estimating density-habitat covariate relationships using spatial capture-recapture models. Because habitat fragmentation and loss are likely to increase in severity globally, further investigation of the influence of habitat fragmentation and detector placement on estimation of this relationship is warranted.

Hidden semi-Markov models reveal multiphasic movement of the endangered Florida panther

Animals must move to find food and mates, and to avoid predators; movement thus influences survival and reproduction, and ultimately determines fitness. Precise description of movement and understanding of spatial and temporal patterns as well as relationships with intrinsic and extrinsic factors is important both for theoretical and applied reasons. We applied hidden semi-Markov models (HSMM) to hourly geographic positioning system (GPS) location data to understand movement patterns of the endangered Florida panther (Puma concolor coryi) and to discern factors influencing these patterns. Three distinct movement modes were identified: (1) Resting mode, characterized by short step lengths and turning angles around 180(o); (2) Moderately active (or intermediate) mode characterized by intermediate step lengths and variable turning angles, and (3) Traveling mode, characterized by long step lengths and turning angles around 0(o). Males and females, and females with and without kittens, exhibited distinctly different movement patterns. Using the Viterbi algorithm, we show that differences in movement patterns of male and female Florida panthers were a consequence of sex-specific differences in diurnal patterns of state occupancy and sex-specific differences in state-specific movement parameters, whereas the differences between females with and without dependent kittens were caused solely by variation in state occupancy. Our study demonstrates the use of HSMM methodology to precisely describe movement and to dissect differences in movement patterns according to sex, and reproductive status.