Integrating multiple genetic detection methods to estimate population density of social and territorial carnivores

Methodological approaches for estimating populations of the endangered dhole Cuon alpinus

Large carnivores are important for maintaining ecosystem integrity and attract much research and conservation interest. For most carnivore species, estimating population density or abundance is challenging because they do not have unique markings for individual identification. This hinders status assessments for many threatened species, and calls for testing new methodological approaches. We examined past efforts to assess the population status of the endangered dhole (Cuon alpinus), and explored the application of a suite of recently developed models for estimating their populations using camera-trap data from India's Western Ghats. We compared the performance of Site-Based Abundance (SBA), Space-to-Event (STE), and Time-to-Event (TTE) models against current knowledge of their population size in the area. We also applied two of these models (TTE and STE) to the co-occurring leopard (Panthera pardus), for which density estimates were available from Spatially Explicit Capture-Recapture (SECR) models, so as to simultaneously validate the accuracy of estimates for one marked and one unmarked species. Our review of literature (n = 38) showed that most assessments of dhole populations involved crude indices (relative abundance index; RAI) or estimates of occupancy and area of suitable habitat; very few studies attempted to estimate populations. Based on empirical data from our field surveys, the TTE and SBA models overestimated dhole population size beyond ecologically plausible limits, but the STE model produced reliable estimates for both the species. Our findings suggest that it is difficult to estimate population sizes of unmarked species when model assumptions are not fully met and data are sparse, which are commonplace for most ecological surveys in the tropics. Based on our assessment, we propose that practitioners who have access to photo-encounter data on dholes across Asia test old and new analytical approaches to increase the overall knowledge-base on the species, and contribute towards conservation monitoring of this endangered carnivore.

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.

Evaluating and integrating spatial capture-recapture models with data of variable individual identifiability

Spatial capture-recapture (SCR) models have become the preferred tool for estimating densities of carnivores. Within this family of models are variants requiring identification of all individuals in each encounter (SCR), a subset of individuals only (generalized spatial mark-resight, gSMR), or no individual identification (spatial count or spatial presence-absence). Although each technique has been shown through simulation to yield unbiased results, the consistency and relative precision of estimates across methods in real-world settings are seldom considered. We tested a suite of models ranging from those only requiring detections of unmarked individuals to others that integrate remote camera, physical capture, genetic, and global positioning system (GPS) data into a hybrid model, to estimate population densities of black bears, bobcats, cougars, and coyotes. For each species, we genotyped fecal DNA collected with detection dogs during a 20-d period. A subset of individuals from each species was affixed with GPS collars bearing unique markings and resighted by remote cameras over 140 d contemporaneous with scat collection. Camera-based gSMR models produced density estimates that differed by <10% from genetic SCR for bears, cougars, and coyotes once important sources of variation (sex or behavioral status) were controlled for. For bobcats, SCR estimates were 33% higher than gSMR. The cause of the discrepancies in estimates was likely attributable to challenges designing a study compatible for species with disparate home range sizes and the difficulty of collecting sufficient data in a timeframe in which demographic closure could be assumed. Unmarked models estimated densities that varied greatly from SCR, but estimates became more consistent in models wherein more individuals were identifiable. Hybrid models containing all data sources exhibited the most precise estimates for all species. For studies in which only sparse data can be obtained and the strictest model assumptions are unlikely to be met, we suggest researchers use caution making inference from models lacking individual identity. For best results, we further recommend the use of methods requiring at least a subset of the population is marked and that multiple data sets are incorporated when possible.

Augustine B, Linden DW, B. Chandler R, Royle JA. 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.

Multiple observation processes in spatial capture-recapture models: How much do we gain?

Population monitoring data may originate from multiple methods and are often sparse and fraught with incomplete information due to practical and economic constraints. Models that can integrate multiple survey methods and are able to cope with incomplete data may help investigators exploit available information more thoroughly. Here, we developed an integrated spatial capture-recapture (SCR) model to incorporate multiple data sources with imperfect individual identification. We contrast inferences drawn from this model with alternate models incorporating only subsets of the data available. Using extensive simulations and an empirical example of multi-method brown bear (Ursus arctos) monitoring data from northern Pakistan, we quantified the benefits of including multiple sources of information in SCR models in terms of parameter precision and bias. Our multiple observation processes SCR model (MOP) yielded a more complete picture of the underlying processes, reduced bias, and led to more precise parameter estimates. Our results suggest that the greatest gains from integrated SCR models can be expected in situations where detection probability is low, a large proportion of detections is not attributable to individuals, and the degree of overlap between individual home ranges is low.