Bobcat Populations Research Articles

Dispersal of juvenile bobcats in a recolonizing population

ABSTRACTNatural recolonization of bobcat (Lynx rufus) populations in the midwestern United States presents challenges for managers with limited knowledge of the species’ population dynamics in a highly fragmented agricultural landscape. Dispersal is a component of population dynamics of a recolonizing population, which is likely influenced by landscape features. To better understand population expansion, we examined dispersal characteristics of juvenile (<2 yr) bobcats in recently recolonized south‐central Iowa, USA, from 2003–2009. We radio‐collared and tracked dispersal fates of 61 individuals (34 males, 27 females), with 29 (22 males, 7 females) dispersing by approximately 2 years of age and the rest remaining philopatric. Most (65%) juvenile males dispersed, whereas only 26% of females dispersed. Initiation of dispersal varied, but none occurred in July–August. Dispersal duration ranged from 4–240 days. Average age at dispersal was 16.9 ± 1.1 (SE) months. Dispersal was most prevalent in an east‐west direction and straight‐line distances ranged from 6.6–203.2 km ( = 57.9 km). On average, males dispersed 44 km farther than females. Land cover composition was similar in natal and post‐dispersal core ranges and consisted predominantly (70%) of forest and grassland. Lower abundance of forest and grassland in areas north of the study area may ultimately limit the ability of juvenile bobcats to successfully disperse and expand throughout the Corn Belt of Iowa. © 2019 The Wildlife Society.

Estimating density and detection of bobcats in fragmented midwestern landscapes using spatial capture–recapture data from camera traps

ABSTRACTCamera‐trapping data analyzed with spatially explicit capture–recapture (SCR) models can provide a rigorous method for estimating density of small populations of elusive carnivore species. We sought to develop and evaluate the efficacy of SCR models for estimating density of a presumed low‐density bobcat (Lynx rufus) population in fragmented landscapes of west‐central Illinois, USA. We analyzed camera‐trapping data from 49 camera stations in a 1,458‐km2 area deployed over a 77‐day period from 1 February to 18 April 2017. Mean operational time of cameras was 52 days (range = 32–67 days). We captured 23 uniquely identifiable bobcats 113 times and recaptured these same individuals 90 times; 15 of 23 (65.2%) individuals were recaptured at ≥2 camera traps. Total number of bobcat capture events was 139, of which 26 (18.7%) were discarded from analyses because of poor image quality or capture of only a part of an animal in photographs. Of 113 capture events used in analyses, 106 (93.8%) and 7 (6.2%) were classified as positive and tentative identifications, respectively; agreement on tentative identifications of bobcats was high (71.4%) among 3 observers. We photographed bobcats at 36 of 49 (73.5%) camera stations, of which 34 stations were used in analyses. We estimated bobcat density at 1.40 individuals (range = 1.00–2.02)/100 km 2. Our modeled bobcat density estimates are considerably below previously reported densities (30.5 individuals/100 km 2) within the state, and among the lowest yet recorded for the species. Nevertheless, use of remote cameras and SCR models was a viable technique for reliably estimating bobcat density across west‐central Illinois. Our research establishes ecological benchmarks for understanding potential effects of colonization, habitat fragmentation, and exploitation on future assessments of bobcat density using standardized methodologies that can be compared directly over time. Further application of SCR models that quantify specific costs of animal movements (i.e., least‐cost path models) while accounting for landscape connectivity has great utility and relevance for conservation and management of bobcat populations across fragmented Midwestern landscapes. © 2019 The Wildlife Society.

History matters: contemporary versus historic population structure of bobcats in the New England region, USA

Habitat fragmentation and genetic bottlenecks can have substantial impacts on the health and management of wildlife species by lowering diversity and subdividing populations. Population genetic comparisons across time periods can help elucidate temporal changes in populations and the processes responsible for the changes. Bobcats (Lynx rufus) are wide-ranging carnivores and are currently increasing in abundance across an expanding range. Bobcat populations in New England have fluctuated in the past century in response to changes in their prey base, harvest pressure, and landscape development. We genotyped contemporary (2010–2017) and historic (1952–1964) bobcats from New England and Quebec, Canada at a suite of microsatellite loci and tested for differences in diversity, effective population size, and gene flow. Over 20 generations separated the sampling periods, and the intervening years were marked by drastic changes in land use and species management regimes. We found a general decrease in genetic diversity and differing population genetic structure through time. Effective population size decreased at the end of the historic period, coincident with a spike in harvest, but rebounded to greater numbers in the contemporary period. Our results suggest that bobcat populations in the region are robust, but development and range dynamics may play a significant role in population structure. Our study also highlights the benefits of a historical perspective in interpreting contemporary population genetic data.

High prevalence of Lynx rufus gammaherpesvirus 1 in wild Vermont bobcats.

Gammaherpesviruses (GHVs) are host specific DNA viruses that infect a large range of mammalian species. These viruses preferentially target host lymphocyte cell populations and infection may lead to morbidity or mortality in immunocompromised, co-infected, or non-adapted hosts. In this study, we tested for the presence of Lynx rufus gammaherpesvirus 1 (LruGHV1) in a northeastern United States population of wild bobcats (L. rufus). We estimated prevalence of infection and viral load in infected individuals using quantitative real-time PCR analysis of spleen DNA from 64 Vermont bobcats. We observed an overall prevalence of 64% using this methodology. Bobcat age was significantly positively associated with GHV infection status, and we noted a trend for higher viral loads in young animals, but prevalence and viral load were similar in male and female bobcats. A single LruGHV1 variant was identified from the sequencing of the viral glycoprotein B gene of Vermont bobcats. This gene sequence was 100% similar to that reported in Florida bobcats and slightly variant from other isolates identified in the Western USA. Our work suggests broad geographic distribution and high prevalence of LruGHV1 in bobcat populations across the United States with infection attributes that suggest horizontal transmission of the agent. Geographic differences in viral genotype may reflect historical migration and expansion events among bobcat populations.

Response of two sympatric carnivores to human disturbances of their habitat: the bobcat and coyote

Information about human land uses, producing the least impact on wildlife and their habitats, is crucial to develop management strategies allowing coexistence between human beings and wildlife. We investigated the responses of the bobcat (Lynx rufus) and the coyote (Canis latrans) to different types of human activity by assessing the abundance, habitat use, and activity patterns of these species in response to extraction of firewood, grazing, and crop farming within their habitat. Photographic records of both carnivores were obtained from 18 camera-trap stations, distributed evenly within these three types of habitat disturbances. The highest relative abundance index for both carnivores was obtained at firewood extraction sites, followed by grazing sites, and croplands. The probability of carnivore sightings varied between species by type of disturbance and by season. The probability of bobcat sightings was highest within firewood extraction sites, with no effect of season, whereas that of the coyote was highest during winter, with no effect of type of human disturbance. An analysis of habitat use revealed that both carnivores avoided cropland sites. Activity patterns were similar, which suggests an absence of temporal segregation between these species; moreover, the activity patterns of both carnivores seems to be affected by intensity of human activities. Preference for woodland habitats by both carnivores denotes that firewood extraction is the disturbance that least affects bobcat and coyote populations; therefore, this kind of human activity could be the least impacting type of disturbance for the coexistence between humans and these carnivore species.

Relative Importance of Habitat Characteristics and Interspecific Relations in Determining Terrestrial Carnivore Occurrence

Carnivores act as top-down regulators in terrestrial ecosystems, and their occurrence and relative abundance is a result of complex interactions between food and habitat availability, human pressure (e.g., trapping, hunting, roadkill) and intraguild interactions (competition, predation). Eastern United States has a long history of human impact, which resulted in an altered carnivore community. Specifically, Ohio presents an interesting case for evaluating the relative roles of interspecific relations and habitat characteristics for shaping the carnivore community, as its carnivore community has a unique dynamics and composition: invasive coyote and red fox (Vulpes vulpes), and native bobcat (Lynx rufus), currently recovering and expanding its range, gray fox (Urocyon cinereoargenteus) declining at a fast pace, and the generalist raccoon (Procyon lotor) and Virginian opossum (Didelphis virginiana). We used 50 camera traps to collect presence/absence data in southeastern Ohio, USA. We hypothesized potential interactions between the six carnivores, and used land cover variables, as well as occupancy probabilities of interacting species, to parameterize to single-species occupancy models. We found that landscape composition at three different scales (500m and 1000 m buffer around camera locations, and 3x3 km grid cell) had little effect on species occurrence. We identified strong negative interspecific relations between carnivores, with bobcat occurrence being influenced by presence of coyotes, red fox occurrence by gray foxes, and raccoon occurrence by Virginia possums. While these findings cannot discriminate between habitat partitioning (spatial or temporal) and competition (direct or interference), they lend support to complex dynamics between invasive coyotes and red foxes and recovering (bobcat) and declining (gray fox) native carnivore species. In particular, the negative relation between the apex predator in our system, C. latrans, and L. rufus, raise further questions on whether direct competition from coyotes has the potential to slow bobcat population recovery. In the context of regulated trapping (ongoing for gray fox and proposed season for bobcat), a better understanding of the carnivore intraguild relations can inform management and conservation actions targeted at minimizing the impact of competition on at-risk native species from non-native species.

Genome-wide expression reveals multiple systemic effects associated with detection of anticoagulant poisons in bobcats (Lynx rufus).

Anticoagulant rodenticides (ARs) are indiscriminate toxicants that threaten nontarget predatory and scavenger species through secondary poisoning. Accumulating evidence suggests that AR exposure may have disruptive sublethal consequences on individuals that can affect fitness. We evaluated AR-related effects on genome-wide expression patterns in a population of bobcats in southern California. We identify differential expression of genes involved in xenobiotic metabolism, endoplasmic reticulum stress response, epithelial integrity and both adaptive and innate immune function. Further, we find that differential expression of immune-related genes may be attributable to AR-related effects on leucocyte differentiation. Collectively, our results provide an unprecedented understanding of the sublethal effects of AR exposure on a wild carnivore. These findings highlight potential detrimental effects of ARs on a wide variety of species worldwide that may consume poisoned rodents and indicate the need to investigate gene expression effects of other toxicants added to natural environments by humans.

Urbanization and anticoagulant poisons promote immune dysfunction in bobcats.

Understanding how human activities influence immune response to environmental stressors can support biodiversity conservation across increasingly urbanizing landscapes. We studied a bobcat (Lynx rufus) population in urban southern California that experienced a rapid population decline from 2002-2005 due to notoedric mange. Because anticoagulant rodenticide (AR) exposure was an underlying complication in mange deaths, we aimed to understand sublethal contributions of urbanization and ARs on 65 biochemical markers of immune and organ function. Variance in immunological variables was primarily associated with AR exposure and secondarily with urbanization. Use of urban habitat and AR exposure has pervasive, complex and predictable effects on biochemical markers of immune and organ function in free-ranging bobcats that include impacts on neutrophil, lymphocyte and cytokine populations, total bilirubin and phosphorus. We find evidence of both inflammatory response and immune suppression associated with urban land use and rat poison exposure that could influence susceptibility to opportunistic infections. Consequently, AR exposure may influence mortality and has population-level effects, as previous work in the focal population has revealed substantial mortality caused by mange infection. The secondary effects of anticoagulant exposure may be a worldwide, largely unrecognized problem affecting a variety of vertebrate species in human-dominated environments.

Widespread Anticoagulant Poison Exposure is Linked with Immune Dysregulation and Severe Notoedric Mange in Urban Bobcats

Human activities threaten wildlife with a variety of novel stressors such as exposure to toxicants. Anticoagulant rodenticides (ARs) are toxicants applied worldwide and through bioaccumulation, threaten species that prey on poisoned rodents or their predators. We studied a population of urban bobcats in southern California that declined rapidly from 2002-2005 due to notoedric mange. We first assessed prevalence of AR exposure using blood and liver samples across the population and found widespread exposure (>90%). Death associated with mange was strongly correlated with cumulative first- and second-generation AR exposure. These findings suggested that exposure to both first- and second-generation ARs were an underlying cause of the disease. We next aimed to understand the sublethal immunological and physiological effects of AR exposure in this natural population. We used two approaches: 1) we used a comprehensive suite of health assays (complete blood counts, blood chemistry assessment, and immunological profiling), and 2) we quantified AR-induced differential gene expression in blood for a subset of individuals. We found that sublethal AR exposure, primarily measured as exposure to diphacinone, is associated with hallmark indicators of generalized systemic inflammation that in persistence could promote immune dysfunction. Further, differential gene expression findings supported the results of immunological profiling. Further, a decrease in the expression of genes associated with epithelial maintenance simultaneous to a decrease in gene expression linked with ectoparasitic immune response may explain the link between AR exposure and mange vulnerability. Such indirect effects of sublethal exposure exemplify the challenge of protecting wild populations from common toxicants in human-dominated environments.

The effects of demographic, social, and environmental characteristics on pathogen prevalence in wild felids across a gradient of urbanization.

Transmission of pathogens among animals is influenced by demographic, social, and environmental factors. Anthropogenic alteration of landscapes can impact patterns of disease dynamics in wildlife populations, increasing the potential for spillover and spread of emerging infectious diseases in wildlife, human, and domestic animal populations. We evaluated the effects of multiple ecological mechanisms on patterns of pathogen exposure in animal populations. Specifically, we evaluated how ecological factors affected the prevalence of Toxoplasma gondii (Toxoplasma), Bartonella spp. (Bartonella), feline immunodeficiency virus (FIV), and feline calicivirus (FCV) in bobcat and puma populations across wildland-urban interface (WUI), low-density exurban development, and wildland habitat on the Western Slope (WS) and Front Range (FR) of Colorado during 2009–2011. Samples were collected from 37 bobcats and 29 pumas on the WS and FR. As predicted, age appeared to be positively related to the exposure to pathogens that are both environmentally transmitted (Toxoplasma) and directly transmitted between animals (FIV). In addition, WS bobcats appeared more likely to be exposed to Toxoplasma with increasing intraspecific space-use overlap. However, counter to our predictions, exposure to directly-transmitted pathogens (FCV and FIV) was more likely with decreasing space-use overlap (FCV: WS bobcats) and potential intraspecific contacts (FIV: FR pumas). Environmental factors, including urbanization and landscape covariates, were generally unsupported in our models. This study is an approximation of how pathogens can be evaluated in relation to demographic, social, and environmental factors to understand pathogen exposure in wild animal populations.

Modeling landscape connectivity for bobcats using expert‐opinion and empirically derived models: how well do they work?

AbstractEfforts to retain ecological connectivity have become a conservation priority to permit animal movements within home ranges, allow dispersal between populations and provide opportunities for animals to respond to climate change. We used expert‐opinion and empirically derived models to investigate landscape connectivity at two spatial scales among bobcats Lynx rufus in New Hampshire, USA. Paths of marked bobcats were compared to random movements in the context of program CircuitScape. At the local scale (within home ranges), the empirical model (based on observations and telemetry locations) performed better than the expert‐opinion model. At the regional scale (state of New Hampshire), both models identified urban development as a potential barrier; however, the models differed in predicting how specific natural features (e.g. mountains and large water bodies) and some roads affected bobcat movements. When compared with bobcat population structure based on genetic information, the expert‐opinion model overestimated the influence of roads. Alternatively, the empirical model overestimated the influence of snow. Our findings indicate that the empirically based resistance model was better at describing landscape‐scale effects, whereas the expert‐opinion model provided a good understanding of gene flow at a regional scale. As such, both models may be considered complementary. Bobcats were sensitive to disruptions imposed by habitat fragmentation and thus may be a suitable focal species for evaluating the consequences of land‐use changes on the regional suite of mesocarnivores.

Origin and genetic structure of a recovering bobcat (Lynxrufus) population

Genetic analyses can provide important insights into the demographic processes that underlie recovering populations of mammals of conservation concern such as felid species. To better understand the recent and rapid recovery of bobcats (Lynx rufus (Schreber, 1777)) in Ohio, we analyzed samples from four states in the lower Great Lakes Region using 12 microsatellite DNA loci and a portion of the mtDNA control region. Our results showed that a newly established population of bobcats in the eastern part of Ohio was genetically distinct from a multistate population distributed across Kentucky, southern Ohio, West Virginia, and western Pennsylvania. There was no direct genetic evidence of a bottleneck or inbreeding in this population. A lack of private alleles and only slightly lower levels of allelic richness and heterozygosity compared with its neighbors suggest that the eastern Ohio population likely originated from the migration of relatively large numbers of individuals from a source population rather than re-emerging from an undetected residual population. We recommend that a management plan should define the areas occupied by the two populations in Ohio as separate management units at least for the near future.

Population and genetic outcomes 20years after reintroducing bobcats (Lynx rufus) to Cumberland Island, Georgia USA.

In 1988-1989, 32 bobcats Lynx rufus were reintroduced to Cumberland Island (CUIS), Georgia, USA, from which they had previously been extirpated. They were monitored intensively for 3years immediately post-reintroduction, but no estimation of the size or genetic diversity of the population had been conducted in over 20years since reintroduction. We returned to CUIS in 2012 to estimate abundance and effective population size of the present-day population, as well as to quantify genetic diversity and inbreeding. We amplified 12 nuclear microsatellite loci from DNA isolated from scats to establish genetic profiles to identify individuals. We used spatially explicit capture-recapture population estimation to estimate abundance. From nine unique genetic profiles, we estimate a population size of 14.4 (SE=3.052) bobcats, with an effective population size (N e) of 5-8 breeding individuals. This is consistent with predictions of a population viability analysis conducted at the time of reintroduction, which estimated the population would average 12-13 bobcats after 10years. We identified several pairs of related bobcats (parent-offspring and full siblings), but ~75% of the pairwise comparisons were typical of unrelated individuals, and only one individual appeared inbred. Despite the small population size and other indications that it has likely experienced a genetic bottleneck, levels of genetic diversity in the CUIS bobcat population remain high compared to other mammalian carnivores. The reintroduction of bobcats to CUIS provides an opportunity to study changes in genetic diversity in an insular population without risk to this common species. Opportunities for natural immigration to the island are limited; therefore, continued monitoring and supplemental bobcat reintroductions could be used to evaluate the effect of different management strategies to maintain genetic diversity and population viability. The successful reintroduction and maintenance of a bobcat population on CUIS illustrates the suitability of translocation as a management tool for re-establishing felid populations.

The effects of urbanization on population density, occupancy, and detection probability of wild felids.

Urbanization is a primary driver of landscape conversion, with far-reaching effects on landscape pattern and process, particularly related to the population characteristics of animals. Urbanization can alter animal movement and habitat quality, both of which can influence population abundance and persistence. We evaluated three important population characteristics (population density, site occupancy, and species detection probability) of a medium-sized and a large carnivore across varying levels of urbanization. Specifically, we studied bobcat and puma populations across wildland, exurban development, and wildland-urban interface (WUI) sampling grids to test hypotheses evaluating how urbanization affects wild felid populations and their prey. Exurban development appeared to have a greater impact on felid populations than did habitat adjacent to a major urban area (i.e., WUI); estimates of population density for both bobcats and pumas were lower in areas of exurban development compared to wildland areas, whereas population density was similar between WUI and wildland habitat. Bobcats and pumas were less likely to be detected in habitat as the amount of human disturbance associated with residential development increased at a site, which was potentially related to reduced habitat quality resulting from urbanization. However, occupancy of both felids was similar between grids in both study areas, indicating that this population metric was less sensitive than density. At the scale of the sampling grid, detection probability for bobcats in urbanized habitat was greater than in wildland areas, potentially due to restrictive movement corridors and funneling of animal movements in landscapes influenced by urbanization. Occupancy of important felid prey (cottontail rabbits and mule deer) was similar across levels of urbanization, although elk occupancy was lower in urbanized areas. Our study indicates that the conservation of medium- and large-sized felids associated with urbanization likely will be most successful if large areas of wildland habitat are maintained, even in close proximity to urban areas, and wildland habitat is not converted to low-density residential development.

Abundancia del gato montés (Lynx rufus escuinapae J. A. Allen, 1903) en el Parque Nacional el Cimatario, Querétaro, México

Bobcat conservation status in central Mexico remains unclear. In many localities it is possible that bobcat populations are declining, so it is very important to estimate the densities of this species. In this note we report an abundance of 3 ± 0.78 bobcats and a minimum density of 0.17 ind/km2 at the Parque Nacional El Cimatario, Querétaro. Information on the species is limited and this data are needed to improve our knowledge in fragmented landscapes.

Bobcats (Lynx rufus) as a Model Organism to Investigate the Effects of Roads on Wide-Ranging Carnivores.

We are using bobcats (Lynx rufus) as a model organism to examine how roads affect the abundance, distribution, and genetic structure of a wide-ranging carnivore. First, we compared the distribution of bobcat-vehicle collisions to road density and then estimated collision probabilities for specific landscapes using a moving window with road-specific traffic volume. Next, we obtained incidental observations of bobcats from the public, camera-trap detections, and locations of bobcats equipped with GPS collars to examine habitat selection. These data were used to generate a cost-surface map to investigate potential barrier effects of roads. Finally, we have begun an examination of genetic structure of bobcat populations in relation to major road networks. Distribution of vehicle-killed bobcats was correlated with road density, especially state and interstate highways. Collision models suggested that some regions may function as demographic sinks. Simulated movements in the context of the cost-surface map indicated that some major roads may be barriers. These patterns were supported by the genetic structure of bobcats. The sharpest divisions among genetically distinct demes occurred along natural barriers (mountains and large lakes) and in road-dense regions. In conclusion, our study has demonstrated the utility of using bobcats as a model organism to understand the variety of threats that roads pose to a wide-ranging species. Bobcats may also be useful as one of a group of focal species while developing approaches to maintain existing connectivity or mitigate the negative effects of roads.

Mark-recapture and mark-resight methods for estimating abundance with remote cameras: a carnivore case study.

Abundance estimation of carnivore populations is difficult and has prompted the use of non-invasive detection methods, such as remotely-triggered cameras, to collect data. To analyze photo data, studies focusing on carnivores with unique pelage patterns have utilized a mark-recapture framework and studies of carnivores without unique pelage patterns have used a mark-resight framework. We compared mark-resight and mark-recapture estimation methods to estimate bobcat (Lynx rufus) population sizes, which motivated the development of a new "hybrid" mark-resight model as an alternative to traditional methods. We deployed a sampling grid of 30 cameras throughout the urban southern California study area. Additionally, we physically captured and marked a subset of the bobcat population with GPS telemetry collars. Since we could identify individual bobcats with photos of unique pelage patterns and a subset of the population was physically marked, we were able to use traditional mark-recapture and mark-resight methods, as well as the new “hybrid” mark-resight model we developed to estimate bobcat abundance. We recorded 109 bobcat photos during 4,669 camera nights and physically marked 27 bobcats with GPS telemetry collars. Abundance estimates produced by the traditional mark-recapture, traditional mark-resight, and “hybrid” mark-resight methods were similar, however precision differed depending on the models used. Traditional mark-recapture and mark-resight estimates were relatively imprecise with percent confidence interval lengths exceeding 100% of point estimates. Hybrid mark-resight models produced better precision with percent confidence intervals not exceeding 57%. The increased precision of the hybrid mark-resight method stems from utilizing the complete encounter histories of physically marked individuals (including those never detected by a camera trap) and the encounter histories of naturally marked individuals detected at camera traps. This new estimator may be particularly useful for estimating abundance of uniquely identifiable species that are difficult to sample using camera traps alone.

Anticoagulant rodenticides in urban bobcats: exposure, risk factors and potential effects based on a 16-year study

Anticoagulant rodenticides (ARs) are increasingly recognized as a threat to nontarget wildlife. High exposure to ARs has been documented globally in nontarget predatory species and linked to the high prevalence of an ectoparasitic disease, notoedric mange. In southern California, mange associated with AR exposure has been the proximate cause of a bobcat (Lynx rufus) population decline. We measured AR exposure in bobcats from two areas in southern California, examining seasonal, demographic and spatial risk factors across landscapes including natural and urbanized areas. The long-term study included bobcats sampled over a 16-year period (1997-2012) and a wide geographic area. We sampled blood (N=206) and liver (N=172) to examine exposure ante- and post-mortem. We detected high exposure prevalence (89%, liver; 39%, blood) and for individuals with paired liver and blood data (N=64), 92% were exposed. Moreover, the animals with the most complete sampling were exposed most frequently to three or more compounds. Toxicant exposure was associated with commercial, residential, and agricultural development. Bobcats of both sexes and age classes were found to be at high risk of exposure, and we documented fetal transfer of multiple ARs. We found a strong association between certain levels of exposure (ppm), and between multiple AR exposure events, and notoedric mange. AR exposure was prevalent throughout both regions sampled and throughout the 16-year time period in the long-term study. ARs pose a substantial threat to bobcats, and likely other mammalian and avian predators, living at the urban-wildland interface.

Predicting bobcat abundance at a landscape scale and evaluating occupancy as a density index in central Wisconsin

The abundance of low-density species like carnivores is logistically difficult to directly estimate at a meaningful scale. Predictive distribution models are often used as a surrogate for density estimation. But because density can continue to increase as occupancy asymptotes at 1, occupancy may have little value as an index, and home range expansion in marginal habitat may further confound the association. We sought to estimate bobcat population size at a landscape scale (14,286 km2) in central Wisconsin, which provided an opportunity to relate predicted occurrence to individual space use and population density. We sampled bobcats using motion-sensitive trail cameras at 9 arrays across central Wisconsin. We estimated bobcat site-specific occupancy, and regressed these estimates as linear or asymptotic functions of site-specific density to determine the strength and shape of their association. We subsequently modeled both parameters relative to habitat covariates and repeated the regression process. A linear functional relationship between density and occupancy was most supported when detection parameters were held constant (wi= 0.97, R2 = 0.72) and when detection, occurrence, and density were modeled as a function of habitat covariates (wi= 0.99, R2 = 0.95). This suggests that repeated presence-absence data alone may be an efficient and reliable method for inferring spatial patterns in bobcat density or identifying habitat types with greater density potential in the northern parts of its range. Bobcat occupancy and density were both positively associated with surrounding woody cover and wetland edge density. Our most supported spatially explicit capture-recapture model estimated bobcat abundance as 362 adult individuals (95% CI 272–490) across the study area. © 2015 The Wildlife Society.

Spatially explicit capture–recapture analysis of bobcat (Lynx rufus) density: implications for mesocarnivore monitoring

Context Accurate density estimation is crucial for conservation and management of elusive species. Camera-trapping may provide an efficient method for density estimation, particularly when analysed with recently developed spatially explicit capture–recapture (SECR) models. Although camera-traps are employed extensively to estimate large carnivore density, their use for smaller carnivores has been limited. Moreover, while camera-trapping studies are typically conducted at local scales, the utility of analysing larger-scale patterns by combining multiple camera studies remains poorly known. Aims The goal of the present study was to develop a better understanding of the utility of SECR models and camera-trapping for the estimation of density of small carnivores at local and regional scales. Methods Based on data collected from camera-traps, we used SECR to examine density of bobcats (Lynx rufus) at four study sites in north-central Texas. We then combined our density estimates with previous estimates (from multiple methodologies) across the bobcat’s geographic range, and used linear regression to examine drivers of range-wide density patterns. Key results Bobcat densities averaged 13.2 per 100 km2 across all four study sites, and were lowest at the site in the most heavily modified landscape. Bobcat capture probability was positively related to forest cover around camera-trap sites. At the range-wide scale, 53% of the variation in density was explained by just two factors: temperature and longitude. Conclusions Our results demonstrate the utility of camera-traps, combined with SECR, to generate precise density estimates for mesocarnivores, and reveal the negative effects of landscape disturbance on bobcat populations. The associations revealed in our range-wide analysis, despite variability in techniques used to estimate density, demonstrate how a combination of multiple density estimates for a species can be used for large-scale inference. However, improvement in our understanding of biogeographic density patterns for mesocarnivores could be obtained from a greater number of camera-based density estimates across the range of a species, combined with meta-analytic techniques. Implications Camera-trapping and SECR should be more widely applied to generate local density estimates for many small and medium-sized carnivores, where at least a portion of the individuals are identifiable. If such estimates are more widely obtained, meta-analytic techniques could be used to test biogeographic predictions or for large-scale monitoring efforts.