When carnivores collide: a review of studies exploring the competitive interactions between bobcats Lynx rufus and coyotes Canis latrans

Competitive interactions and coexistence of sympatric flagship carnivores in Asia

Understanding the competition and coexistence of flagship carnivores is key to creating strategies for their conservation in the face of global carnivore declines. Although studies exploring the dynamics and competition between tigers (Panthera tigris) and leopards (P. pardus) span decades, there is a lack of understanding regarding the factors that influence their coexistence mechanisms on a broad scale, as well as the drivers determining their exploitative and interference competition. We gathered a comprehensive list of research papers among which 36 papers explored the interspecific interactions between tigers and leopards and tested the influence of biotic and abiotic factors on the coexistence mechanisms along three dimensions using multiple response variables regression models; we also tested the influence of ecological drivers determining the exploitative or interference competition between tigers and leopards. Elevation and ungulate density were the most important predictors in regulating the coexistence mechanisms. Tigers and leopards exhibited more positive relations/higher overlaps as elevation increased in the spatial niche. In addition, they showed a higher dietary overlap in the prey-rich regions. We determined that interference competition between tigers and leopards was less frequently observed in habitats with dense tree cover and homogeneous vegetation structures. Meanwhile, studies with multiple metrics would promote the detection of interference competition. Our study provides new insight into the competitive interactions and coexistence mechanisms of tigers and leopards on a broad scale. Policy-makers and managers should pay more attention to the factors of elevation, prey abundance, and habitat structures for the conservation of tigers and leopards.

Assessing mammal population densities in response to urbanization using camera trap distance sampling

Environmental filtering is deemed to play a predominant role in regulating the abundance and distribution of animals during the urbanization process. However, the current knowledge about the effects of urbanization on the population densities of terrestrial mammals is limited. In this study, we compared two invasive mammals (dogs Canis lupus familiaris and cats Felis silvestris) and three indigenous mammals (Siberian weasels Mustela sibirica, Amur hedgehogs Erinaceus amurensis, and Tolai hares Lepus tolai) in response to urbanization using camera trap distance sampling (CTDS) in the rural-urban landscape of Tianjin, China. We used generalized additive mixed models (GAMMs) to test the specific responses of their densities to levels of urbanization. Invasive dogs (2.63 individuals/km2, 95% CI: 0.91-7.62) exhibited similar density estimations to cats (2.15 individuals/km2, 95% CI: 1.31-3.50). Amur hedgehogs were the most abundant species (6.73 individuals/km2, 95% CI: 3.15-14.38), followed by Tolai hares (2.22 individuals/km2, 95% CI: 0.87-5.68) and Siberian weasels (2.15 individuals/km2, 95% CI: 1.06-4.36). The densities of cats, Siberian weasels, and Amur hedgehogs increased with the level of urbanization. The population densities of dogs and cats were only influenced by urban-related variables, while the densities of Siberian weasels and Amur hedgehogs were influenced by both urban-related variables and nature-related variables. Our findings highlight that the CTDS is a suitable and promising method for wildlife surveys in rural-urban landscapes, and urban wildlife management needs to consider the integrated repercussions of urban- and nature-related factors, especially the critical impacts of green space habitats at finer scales.

White-tailed deer (Odocoileus virginianus) fawn survival and the influence of landscape characteristics on fawn predation risk in the Southern Appalachian Mountains, USA

In the Southern Appalachian region of the United States, harvest data has indicated the occurrence of low deer densities while exposing a trend of declining white-tailed deer (Odocoileus virginianus) populations over the past several decades in northern Georgia. A triumvirate of increasing fawn predator populations reside in the Southern Appalachian Mountains including coyotes (Canis latrans), black bears (Ursus americanus) and bobcats (Lynx rufus). This region is also characterized by a homogenous landscape composed of mature forests and sparse understory vegetation, likely lacking adequate cover to offer fawns refugia from predators. Our objectives were to estimate survival and cause-specific mortality rates of fawns while assessing a possible link between mortality risk, intrinsic fawn characteristics (i.e., birth mass, Julian birth date, sibling status), and landscape features within fawn usage areas. During 2018-2020, we radio-collared 71 fawns within the Chattahoochee National Forest of northern Georgia, USA and monitored survival to 12 weeks of age. We observed low fawn survival (cumulative = 0.157, 95% CI = 0.091-0.273; vaginal implant transmitter = 0.196, 95% CI = 0.096-0.403) with predation as the leading cause of all known mortalities (45 of 55 mortalities; 82%) due primarily to coyotes (n = 22), black bears (n = 12), and bobcats (n = 7). Relationships between landscape features and fawn predation risk were minimal with only one informative covariate. Increasing amounts of early successional land cover within fawn usage areas decreased fawn mortality risk within the first 20 days of life, but elevated mortality risk thereafter. All fawns with any amount of early successional land cover in their usage areas died of predation (n = 13) at various time intervals, suggesting limited areas of potential fawning cover may be targeted by predators. However, fawn predation risk seemed to be high regardless of landscape covariates due to the limited number of surviving fawns. Coyote-caused mortality occurred over a longer period at a consistently higher magnitude than all other forms of mortality, indicating possible delayed prey-switching behavior and coyote predation as an important factor of fawn survival. The low recruitment of fawns influenced by high predation rates and homogenous habitat conditions is likely the cause of deer population declines in the region.

Multiscale assessment of habitat selection and avoidance of sympatric carnivores by the endangered ocelot

Habitat selection by animals is a complex, dynamic process that can vary across spatial and temporal scales. Understanding habitat selection is a vital component of managing endangered species. Ocelots (Leopardus pardalis), a medium-sized endangered felid, overlap in their northern range with bobcats (Lynx rufus) and coyotes (Canis latrans), with all three species sharing similar space and resource use. As the potential for competition between these three carnivores is high, understanding differences in habitat use and the effect of these potential competitors on habitat selection of ocelots is essential to conservation. Our objective was to compare habitat selection between species and examine if ocelots avoided areas used by competitors at broad and fine scales. We captured and collared 8 ocelots, 13 bobcats, and 5 coyotes on the East Foundation's El Sauz Ranch and the Yturria San Francisco Ranch in South Texas, USA from 2017 to 2021. We compared 2nd (position of home range) and 3rd (use within the home range) order selection across species and examined whether ocelots avoided areas categorized as high probability of use by bobcats and coyotes across both orders of selection. We found a preference for heterogeneous landscapes by bobcats and coyotes while ocelots were strongly tied to woody cover across both orders. At the 2nd order, ocelots selected areas with higher probability of use by bobcats and showed no response to higher probability of use by coyotes, suggesting ocelots did not avoid either species. However, at the 3rd order, ocelots avoided areas used by coyotes. Ocelots selected for areas of use by bobcats at the 2nd order and 3rd order. Results suggest that at the broader scale, placement of the home range is not affected by the presence of sympatric carnivores, however, at a finer scale, ocelots are avoiding coyotes but not bobcats. Our study emphasizes the importance of woody and herbaceous cover at the broad scale and dense vegetation at the finer scale to sustain ocelots. In addition, we show differing patterns of interspecific avoidance by ocelots across species and scales.

Syntopic species interact with large boreal mammals' response to anthropogenic landscape change

Landscape change alters species' distributions, and understanding these changes is a key ecological and conservation goal. Species-habitat relationships are often modelled in the absence of syntopic species, but niche theory and emerging empirical research suggests heterospecifics should entrain (and statistically explain) variability in distribution, perhaps synergistically by interacting with landscape features. We examined the effects of syntopic species in boreal mammals' relationship to landscape change, using three years of camera-trap data in the western Nearctic boreal forest. Using an information-theoretic framework, we weighed evidence for additive and interactive variables measuring heterospecifics' co-occurrence in species distribution models built on natural and anthropogenic landscape features. We competed multiple hypotheses about the roles of natural features, anthropogenic features, predators, competitors, and species-habitat interaction terms in explaining relative abundance of carnivores, herbivores, and omnivores/scavengers. For most species, models including heterospecifics explained occurrence frequency better than landscape features alone. Dominant predator (wolf) occurrence was best explained by prey, while prey species were explained by apparent competitors and subdominant predators. Evidence for interactions between landscape features and heterospecifics was strong for coyotes and wolves but variable for other species. Boreal mammals' spatial distribution is a function of heterospecific co-occurrence as well as landscape features, with synergistic effects observed for most species. Understanding species' responses to anthropogenic landscape change thus requires a multi-taxa approach that incorporates interspecific relationships, enabling better inference into underlying processes from observed patterns.

Dracula’s ménagerie: A multispecies occupancy analysis of lynx, wildcat, and wolf in the Romanian Carpathians

The recovery of terrestrial carnivores in Europe is a conservation success story. Initiatives focused on restoring top predators require information on how resident species may interact with the re‐introduced species as their interactions have the potential to alter food webs, yet such data are scarce for Europe.

In this study, we assessed patterns of occupancy and interactions between three carnivore species in the Romanian Carpathians. Romania houses one of the few intact carnivore guilds in Europe, making it an ideal system to assess intraguild interactions and serve as a guide for reintroductions elsewhere.

We used camera trap data from two seasons in Transylvanian forests to assess occupancy and co‐occurrence of carnivores using multispecies occupancy models.

Mean occupancy in the study area was highest for lynx (Ψ_winter = 0.76 95% CI: 0.42–0.92; Ψ_autumn = 0.71 CI: 0.38–0.84) and wolf (Ψ_winter = 0.60 CI: 0.34–0.78; Ψ_autumn = 0.81 CI: 0.25–0.95) and lowest for wildcat (Ψ_winter = 0.40 CI: 0.19–0.63; Ψ_autumn = 0.52 CI: 0.17–0.78)

We found that marginal occupancy predictors for carnivores varied between seasons. We also found differences in predictors of co‐occurrence between seasons for both lynx‐wolf and wildcat‐wolf co‐occurrence. For both seasons, we found that conditional occupancy probabilities of all three species were higher when another species was present.

Our results indicate that while there are seasonal differences in predictors of occupancy and co‐occurrence of the three species, co‐occurrence in our study area is high.

Terrestrial carnivore recovery efforts are ongoing worldwide. Insights into interspecific relations between carnivore species are critical when considering the depauperate communities they are introduced in. Our work showcases that apex carnivore coexistence is possible, but dependent on protection afforded to forest habitats and their prey base.

Habitat connectivity and resource selection in an expanding bobcat (Lynx rufus) population

Terrestrial carnivores are among the most imperiled species worldwide, yet some species are resilient and are recovering in human-dominated landscapes after decades or centuries of absence. Bobcat (Lynx rufus) populations were extirpated from much of Midwestern US in the mid-1800's, and are currently expanding and recolonizing their former range. In this study, we investigated multi-scale habitat selection for Ohio's expanding bobcat population, and examined habitat connectivity in order to evaluate the conduits for dispersal statewide. We used citizen observations collected between 1978 and 2019 and logistic regression to evaluate population-level habitat selection, and GPS telemetry data for 20 individuals collected between 2012 and 2014 and a distribution-weighted exponential Resource Selection Function to evaluate individual-level habitat selection within home ranges. At the population level, bobcats selected for higher amounts of forest and pasture (at a 50 km² scale) and herbaceous vegetation (at 15-50 50 km² scales), thus overall heterogeneous forested habitat. At individual (home range) level, bobcats selected for forested habitats with low road density and farther away from high traffic roads; they also showed weak selection for open habitat at the home range level. Male home ranges were significantly greater than female home ranges. Lastly, we used the population-level spatial outputs (i.e. habitat suitability map) to parameterize habitat connectivity models using circuit theory in the program Circuitscape. We tested three relationships between habitat suitability and resistance to movement and used a subset of data on potential dispersing individuals to evaluate which relationship performed best. All three relationships performed almost equally well, and we calculated a weighted averaged connectivity map as our final map. Habitat was highly permeable to movements between core areas of two genetically distinct subpopulations located in southeastern Ohio. We also identified potential dispersal corridors from the core areas to other regions of Ohio dominated by agriculture and suburban development via forested riparian corridors. Overall, our analysis offers new information on habitat selection and connectivity in a rebounding felid population and offers important ecological information for wildlife management strategies. We recommend that the suitability and connectivity models should be periodically updated until the population reaches an equilibrium, and be integrated with data from neighboring states for a comprehensive assessment of a conservation success story.