Spatial capture-recapture with random thinning for unidentified encounters

Low signs of territorial behavior in the Eurasian otter during low-water conditions in a Mediterranean river

The Eurasian otter Lutra lutra is a territorial semi-aquatic carnivore usually found at low densities in rivers, coastal areas, and wetlands. Its diet is based on prey associated with aquatic environments. Mediterranean rivers are highly seasonal, and suffer reduced flow during the summer, resulting in isolated river sections (pools) that sometimes can be left with a minimal amount of water, leading to concentrations of food for otters. To our knowledge, this process, which was known to field naturalists, has not been accurately described, nor have otter densities been estimated under these conditions. In this study, we describe the population size and movements of an aggregation of otters in an isolated pool in the Guadiana River in the Tablas de Daimiel National Park (central Spain), which progressively dried out during the spring-summer of 2022, in a context of low connectivity due to the absence of circulating water in the Guadiana and Gigüela rivers. Using non-invasive genetic sampling of 120 spraints collected along 79.4 km of sampling transects and spatial capture-recapture methods, we estimated the otter density at 1.71 individuals/km of river channel length (4.21 individuals/km2) in a progressively drying river pool, up to five times higher than previously described in the Iberian Peninsula. The movement patterns obtained with the spatial capture-recapture model are not quite different from those described in low density, which seems to indicate a wide home range overlap, with low signs of territoriality.

A flexible framework for spatial capture-recapture with unknown identities

Camera traps or acoustic recorders are often used to sample wildlife populations. When animals can be individually identified, these data can be used with spatial capture-recapture (SCR) methods to assess populations. However, obtaining animal identities is often labor-intensive and not always possible for all detected animals. To address this problem, we formulate SCR, including acoustic SCR, as a marked Poisson process, comprising a single counting process for the detections of all animals and a mark distribution for what is observed (eg, animal identity, detector location). The counting process applies equally when it is animals appearing in front of camera traps and when vocalizations are captured by microphones, although the definition of a mark changes. When animals cannot be uniquely identified, the observed marks arise from a mixture of mark distributions defined by the animal activity centers and additional characteristics. Our method generalizes existing latent identity SCR models and provides an integrated framework that includes acoustic SCR. We apply our method to estimate density from a camera trap study of fisher (Pekania pennanti) and an acoustic survey of Cape Peninsula moss frog (Arthroleptella lightfooti). We also test it through simulation. We find latent identity SCR with additional marks such as sex or time of arrival to be a reliable method for estimating animal density.

Wildlife Population Assessment: Changing Priorities Driven by Technological Advances

Advances in technology are having a large effect on the priorities for innovation in statistical ecology. Collaborations between statisticians and ecologists have always been important in driving methodological development, but increasingly, expertise from computer scientists and engineers is also needed. We discuss changes that are occurring and that may occur in the future in surveys for estimating animal abundance. As technology advances, we expect classical distance sampling and capture-recapture to decrease in importance, as camera (still and video) survey, acoustic survey, spatial capture-recapture and genetic methods continue to develop and find new applications. We explore how these changes are impacting the work of the statistical ecologist.

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.

Applications of non-intrusive methods to study the sand cat: a field study in the Sahara Desert

Surveys based on indirect signs and camera trapping are two non-invasive methods extensively used for monitoring elusive mammals. Both approaches can be useful to obtain key information on wildlife in remote areas, since they may allow for the logistically viable design of optimal field frameworks. The sand cat (Felis margarita) is a feline that inhabits the Sahara Desert, the Arabian Peninsula, and western Asian deserts. Its basic ecology is poorly known and the status and impacts of threats are difficult to assess. Some local population declines have been detected, and more research is needed. Based on field surveys carried out in the Atlantic Sahara, we have evaluated the applications of both methods to study this species. Our results show that (a) camera trapping provided reliable data on several key aspects of its ecology, (b) walking surveys to collect feces for molecular data failed completely, and (c) for footprints, identification problems and the marked effects of the absence of optimal substrates and the prevalence of wind are relevant handicaps. Beyond this evaluation, we provide for the first time some key aspects of the ecology of sand cats in the Sahara Desert, including habitat selection, density, diel activity, and predator–prey relationships.

Occupancy data improves parameter precision in spatial capture-recapture models

Population size is one of the basic demographic parameters for species management and conservation. Among different estimation methods, spatially explicit capture-recapture (SCR) models allow the estimation of population density in a framework that has been greatly developed in recent years. The use of automated detection devices, such as camera traps, has impressively extended SCR studies for individually identifiable species. However, its application to unmarked/partially marked species remains challenging, and no specific method has been widely used. We fitted an SCR-integrated model (SCR-IM) to stone marten Martes foina data, a species for which only some individuals are individually recognizable by natural marks, and estimate population size based on integration of three submodels: (1) individual capture histories from live capture and transponder tagging; (2) detection/nondetection or "occupancy" data using camera traps in a bigger area to extend the geographic scope of capture-recapture data; and (3) telemetry data from a set of tagged individuals. We estimated a stone marten density of 0.352 (SD: 0.081) individuals/km2. We simulated four dilution scenarios of occupancy data to study the variation in the coefficient of variation in population size estimates. We also used simulations with similar characteristics as the stone marten case study, comparing the accuracy and precision obtained from SCR-IM and SCR, to understand how submodels' integration affects the posterior distributions of estimated parameters. Based on our simulations, we found that population size estimates using SCR-IM are more accurate and precise. In our stone marten case study, the SCR-IM density estimation increased the precision by 37% when compared to the standard SCR model as regards to the coefficient of variation. This model has high potential to be used for species in which individual recognition by natural markings is not possible, therefore limiting the need to rely on invasive sampling procedures.

Estimating red fox density using non-invasive genetic sampling and spatial capture-recapture modelling

Spatial capture-recapture modelling (SCR) is a powerful tool for estimating density, population size, and space use of elusive animals. Here, we applied SCR modelling to non-invasive genetic sampling (NGS) data to estimate red fox (Vulpes vulpes) densities in two areas of boreal forest in central (2016-2018) and southern Norway (2017-2018). Estimated densities were overall lower in the central study area (mean = 0.04 foxes per km2 in 2016, 0.10 in 2017, and 0.06 in 2018) compared to the southern study area (0.16 in 2017 and 0.09 in 2018). We found a positive effect of forest cover on density in the central, but not the southern study area. The absence of an effect in the southern area may reflect a paucity of evidence caused by low variation in forest cover. Estimated mean home-range size in the central study area was 45 km2 [95%CI 34-60] for females and 88 km2 [69-113] for males. Mean home-range sizes were smaller in the southern study area (26 km2 [16-42] for females and 56 km2 [35-91] for males). In both study areas, detection probability was session-dependent and affected by sampling effort. This study highlights how SCR modelling in combination with NGS can be used to efficiently monitor red fox populations, and simultaneously incorporate ecological factors and estimate their effects on population density and space use.

Evaluating otter reintroduction outcomes using genetic spatial capture-recapture modified for dendritic networks

Monitoring the demographics and genetics of reintroduced populations is critical to evaluating reintroduction success, but species ecology and the landscapes that they inhabit often present challenges for accurate assessments. If suitable habitats are restricted to hierarchical dendritic networks, such as river systems, animal movements are typically constrained and may violate assumptions of methods commonly used to estimate demographic parameters. Using genetic detection data collected via fecal sampling at latrines, we demonstrate applicability of the spatial capture-recapture (SCR) network distance function for estimating the size and density of a recently reintroduced North American river otter (Lontra canadensis) population in the Upper Rio Grande River dendritic network in the southwestern United States, and we also evaluated the genetic outcomes of using a small founder group (n = 33 otters) for reintroduction. Estimated population density was 0.23-0.28 otter/km, or 1 otter/3.57-4.35 km, with weak evidence of density increasing with northerly latitude (β = 0.33). Estimated population size was 83-104 total otters in 359 km of riverine dendritic network, which corresponded to average annual exponential population growth of 1.12-1.15/year since reintroduction. Growth was ≥40% lower than most reintroduced river otter populations and strong evidence of a founder effect existed 8-10 years post-reintroduction, including 13-21% genetic diversity loss, 84%-87% genetic effective population size decline, and rapid divergence from the source population (F ST accumulation = 0.06/generation). Consequently, genetic restoration via translocation of additional otters from other populations may be necessary to mitigate deleterious genetic effects in this small, isolated population. Combined with non-invasive genetic sampling, the SCR network distance approach is likely widely applicable to demogenetic assessments of both reintroduced and established populations of multiple mustelid species that inhabit aquatic dendritic networks, many of which are regionally or globally imperiled and may warrant reintroduction or augmentation efforts.