Abstract
Carnivores play key roles in maintaining ecosystem structure and function as well as ecological processes. Understanding how sympatric species coexist in natural ecosystems is a central research topic in community ecology and biodiversity conservation. In this study, we explored intra- and interspecific niche partitioning along spatial, temporal, and dietary niche partitioning between apex carnivores (wolf Canis lupus, snow leopard Panthera uncia, Eurasian lynx Lynx lynx) and mesocarnivores (Pallas’s cat Otocolobus manul, red fox Vulpes vulpes, Tibetan fox Vulpes ferrilata) in Qilian Mountain National Park, China using camera trapping data and DNA metabarcoding sequencing data. Our study showed that apex carnivore species had more overlap temporally (coefficients of interspecific overlap ranging from 0.661 to 0.900) or trophically (Pianka’s index ranging from 0.458 to 0.892), mesocarnivore species had high dietary overlap with each other (Pianka’s index ranging from 0.945 to 0.997), and apex carnivore and mesocarnivore species had high temporal overlap (coefficients of interspecific overlap ranging from 0.497 to 0.855). Large dietary overlap was observed between wolf and snow leopard (Pianka’s index=0.892) and Pallas’s cat and Tibetan fox (Pianka’s index=0.997), suggesting the potential for increased resource competition for these species pairs. We concluded that spatial niche partitioning is likely to key driver in facilitating the coexistence of apex carnivore species, while spatial and temporal niche partitioning likely facilitate the coexistence of mesocarnivore species, and spatial and dietary niche partitioning facilitate the coexistence between apex and mesocarnivore species. Our findings consider partitioning across temporal, spatial and dietary dimensions while examining diverse coexistence patterns of carnivore species in Qilian Mountain National Park, China. These findings will contribute substantially to current understanding of carnivore guilds and effective conservation management in fragile alpine ecosystems.
eLife assessment
This paper provides an important assessment of competition dynamics allowing coexistence of the carnivore guild within a large national park in China. Multiple surveying techniques (camera traps and DNA metabarcoding) provide convincing evidence that spatial segregation represents the main strategy of coexistence, while species have a certain degree of temporal and dietary overlap. Altogether, the manuscript provides information critical to the conservation and management agenda of the park.
Introduction
Carnivores play key roles in maintaining ecosystem structure and function, as well as ecological processes (Ripple et al., 2014). Apex carnivores are often classically specialized hunters occupying top trophic positions that work to suppress the number of herbivores and mesocarnivores through predation, competition and trophic cascades (Ripple et al., 2014; Newsome et al., 2017; Ritchie and Johnson, 2009). Understanding how sympatric carnivore species coexist in natural ecosystems is a central focus in community ecology and biodiversity conservation (Chu et al., 2017). Human activities and climate change are causing large-scale habitat loss and fragmentation, posing significant challenges to carnivore interspecies relationships (Li et al., 2021; Manlick and Pauli, 2020; Ripple et al., 2014). Such relationships are becoming increasingly intricate, leading to greater research efforts aimed at determining the strategies that facilitate coexistence despite intense competition (Kuijper et al., 2016; Smith et al., 2017).
The competitive exclusion principle dictates that species with similar ecological requirements are unable to successfully coexist (Hardin, 1960; Gause, 1934). Thus, carnivores within a guild occupy different ecological niches based on a combination of three niche dimensions, i.e. spatial, temporal, and trophic (Schoener, 1974). Spatially, carnivore species within the same geographic area exhibit distinct distributions that minimize overlap in resource use and competition. For example, carnivores can partition habitats based on habitat feature preferences and availability of prey (De Satgé et al., 2017; Garrote and Pérez De Ayala, 2019; Gołdyn et al., 2003; Strampelli et al., 2023). Temporally, differences in seasonal or daily activity patterns among sympatric carnivores can reduce competitive interactions and facilitate coexistence. For example, carnivores can exhibit temporal segregation in their foraging behaviors, such as diurnal versus nocturnal activity, to avoid direct competition (Finnegan et al., 2021; Nasanbat et al., 2021; Searle et al., 2021). Trophically, carnivore species can diversify their diets to exploit different prey species or sizes, thereby reducing competition for food resources. For example, carnivores can exhibit dietary specialization to optimize their foraging efficiency and minimize competitive pressures (Steinmetz et al., 2021). Currently, research on the niche partitioning of sympatric carnivores primarily focuses on one or two dimensions, with little attention paid to all three (Alexander et al., 2016; Karanth et al., 2017; Li et al., 2022a; Santos et al., 2019; Shao et al., 2021; Strampelli et al., 2023; Tsunoda et al., 2020; Vilella et al., 2020). This limitation restricts the multi-scale understanding of coexistence mechanisms among carnivores.
In recent years, camera trapping and DNA metabarcoding technology has been widely used in wildlife monitoring and research. Camera trapping enables the monitoring of elusive species in remote areas (Boitani, 2016; Palencia et al., 2022), and has become a helpful tool for accumulating large amounts of time-recorded data that can provide detailed information surrounding animal behavior and activity patterns, which is widely used to investigate species interspecies relationships (Frey et al., 2017). For example, Li et al. (2019) used camera trap data to conclude that temporal segregation is a key mechanism for promoting the coexistence of tigers (Panthera tigris) and leopards (P. pardus). DNA metabarcoding provides a noninvasive molecular tool that is more accurate than traditional dietary analysis methods, such as microhistology (Deiner et al., 2017; Newmaster et al., 2013). DNA metabarcording technology has been applied to many studies and has offered profound insight into the ecology, conservation and biological monitoring of rare and endangered species (Deagle et al., 2019; Kartzinel et al., 2015). For example, Shao et al. (2021) found that dietary niche partitioning promoted the coexistence of species in the mountains of southwestern China based on DNA metabarcoding. The combination of camera trapping and DNA metabarcoding can result in robust data for exploring the mechanisms of coexistence surrounding carnivore guilds.
The Qilian Mountains constitute a biodiversity hotspot with one of the richest carnivore assemblages in the world (Di Minin et al., 2016). To understand the competition and coexistence among different carnivore species in this area, we explored the habitat use, activity patterns, and prey item composition of sympatric carnivore species comprised of apex carnivores and mesocarnivores across Qilian Mountain National Park using camera trap data and DNA metabarcoding data. Based on theories surrounding resource partitioning and niche differentiation, as well as studies on interactions and coexistence among carnivorous species (Haswell et al., 2018; Linnell and Strand, 2000), we hypothesized that differentiation along one or more niche axes is beneficial for the coexistence of the carnivorous guild in the Qilian Mountains. We expected that spatial niche differentiation promotes the coexistence of large carnivores in the Qilian Mountain region, as they are more likely than small carnivores to spatially avoid interspecific competition (Davis et al., 2018). Mesocarnivores may coexist either spatially or temporally due to increased interspecific competition for similar prey (Di Bitetti et al., 2010; Donadio and Buskirk, 2006). Dietary niche differentiation may be a significant factor for promoting coexistence between large and mesocarnivore species due to differences in body size (Gómez-Ortiz et al., 2015; Lanszki et al., 2019).
Material and methods
Study sites
The Qilian Mountains laterally span Gansu and Qinghai Provinces in China, located on the northeastern edge of the Qinghai-Tibetan Plateau (Figure 1). Qilian Mountain National Park covers an area of approximately 52,000 km2, with an average elevation of over 3000 m. The area is an alpine ecosystem with a typical plateau continental climate. The average annual temperature is below -4℃ and the average annual rainfall is about 400 mm, with habitats mainly consisting of deserts, grassland, meadows, and wetland (Zheng, 2011). Wildlife present include the wolf (Canis lupus), snow leopard (P. uncia), Eurasian lynx (Lynx lynx), red fox (Vulpes vulpes), Tibetan fox (V. ferrilata), Tibetan brown bear (Ursus arctos), Chinese mountain cat (Felis bieti), wild yak (Bos mutus), blue sheep (Pseudois nayaur), alpine musk-deer (Moschus chrysogaster), Tibetan antelope (Pantholops hodgsonii), Himalayan marmot (Marmota himalayana), woolly hare (Lepus oiostolus), and plateau pika (Ochotona curzoniae), among others (Ma et al., 2021; Xue et al., 2019).
![](https://prod--epp.elifesciences.org/iiif/2/90559%2Fv2%2Fcontent%2F553226v2_fig1.tif/full/max/0/default.jpg)
Locations of camera trap stations and scat collection sites in this study.
Camera-trap monitoring and non-invasive sampling
The study area was subdivided into sample squares of 25 km2 (5×5 km) as a geographical reference for placing camera survey sites and collecting scat samples (Xue et al., 2019). Species occurrence was recorded using camera-trap monitoring (Ltl-6210; Shenzhen Ltl Acorn Electronics Co. Ltd). Two camera traps were placed in each square with a distance of at least 1 km between them. However, due to limitations of terrain, landform, road accessibility and other factors, the number of camera trap in some squares was adjusted in accordance to field conditions. Camera traps were set in areas believed to be important to and heavily used by wildlife, such as the bottoms of cliffs, sides of boulders, valleys and ridges along movement corridors. Taking into account the fact that mammalian communities are sensitive to seasonality, we used camera traps to monitor animals with an extensive survey effort from December 2016 to February 2022, covering the activity of animal species in different seasons, which can reflect the overall distribution of carnivores. We placed a total of 280 infrared cameras at the study site, set them to be active for 4 to 6 months, and considered possible relocation to another position based on animal detection in an effort to improve estimates of the occupancy and detection rates for both common and rare species (Figure 1) (Kays et al., 2020). The camera trap was set to record the time and date on a 24 hr clock when triggered, and to record a 15s video and 1 photo with an interval of 2 minutes between any two consecutive triggers. The sum of camera trap effective days was defined by the total amount of trapping effort during the sampling period, which was calculated from the time the camera was placed in operation to the time the last video or photograph was taken. We visited each camera 2 to 3 times a year to download photos and check batteries.
Non-invasive sampling of scats believed to be of carnivore origin were collected within camera trapping areas. A small portion of scat (approximately 1/3) was broken off and stored in a 15 ml centrifuge tubes with silica desiccant covered by clean filter paper to separate the desiccant from the scat (Janecka et al., 2008). Gloves were replaced between sampling to avoid cross-contamination. Sampling place, date, and sample number were labeled on the tube; GPS coordinates, elevation, and nearby landscape features were recorded on the sample collection sheet (Hacker et al., 2021). A total 480 scat samples were collected from April 2019 to June 2021 (Figure 1).
Data Analysis
Spatial analysis
To investigate the spatial distribution of carnivores, as well as the influence of environmental factors on the site occupancy of species in the study area, we performed single-season, single-species occupancy models to estimate carnivores’ occupancy (ψ) and detection (Pr) probability (Li et al., 2022b; MacKenzie, 2018; Moreno-Sosa et al., 2022). To ensure capture independence, only photo or video records at intervals of 30 min were was included in the data analysis (Li et al., 2020). We created a matrix recording whether each carnivore species was detected (1) or not (0) across several 30-day intervals (that is 0-30, 31-60, 61-90, 91-120, 121-150, >150 days) for each camera location. Based on the previous studies of habitat use of carnivores (Greenspan and Giordano, 2021; Alexander et al., 2016; Gorczynski et al., 2022), we selected terrain, vegetation, biological factors and disturbance to construct the model. Terrain is a fundamental element of wildlife habitat and closely linked to other environmental factors (Chen et al., 2024). Terrain variables include elevation (ele) and roughness index (rix). Vegetation variables include normalized difference vegetation index (ndvi), and provide information on the level of habitat concealment. Biological variables include prey abundance (the number of independent photos of their preferred prey based on dietary analysis in this study, wolf and snow leopard: artiodactyla including livestock; Eurasian lynx and Pallas’s cat: lagomorpha; red fox and Tibetan fox: lagomorpha and rodentia) and reflect habitat preference and distribution patterns of carnivores. Disturbance variables include distance to roads (disrd) and human disturbances (hdis, the number of independent photos of herdsman and livestock) and can provide insight into the habitat selection and behavior patterns of carnivores. In addition, we used elevation, human disturbance and prey as covariates that affect detection rate. Road data were obtained from Open Street Map (OSM, https://www.openstreetmap.org). Others environmental data were obtained from the Resource and Environment Science and Data Center (https://www.resdc.cn). We fitted all possible combinations of covariates and used Akaike’s information criterion (AIC) to rank candidate models, and selected ΔAIC≤2 model as the optimal model. If more than one optimal model resulted, then covariate estimates were obtained by using the equal-weight average.
Carnivore co-occurrence was evaluated using the Sørensen similarity index (value = 0, indicating maximum segregation and value = 1, indicating maximum co-occurrence) based on binary presence absence data within the 5 km × 5 km camera trap grid (Torretta et al., 2021; Sorensen, 1948). Spatial analyses were performed using ArcGIS 10.8 (ESRI Inc.), the “vegan” packages (Oksanen et al., 2019) and “unmarked” package (Fiske and Chandler, 2011) for R studio (version 1.1.463).
Temporal analysis
Estimates of the coefficient of overlap (Δ) for activity patterns were estimated using the non-parametric kernel density method and applying time data obtained by the camera traps. Because the smallest sample had more than 50 records, we used the Δ4 estimator for pairwise comparisons between carnivore activity patterns, and used a smooth bootstrap scheme to generate 1000 resamples with 95% confidence intervals to test the reliability of the overlap value (Ridout and Linkie, 2009). Activity pattern analyses were performed using the “overlap” R package. Values of the Δ4 estimates were calculated relative to 1000 simulated null models of randomized overlap data using the ‘compareCkern’ function in the “activity” R package to test for differences in daily activity patterns (Ridout and Linkie, 2009; Rowcliffe et al., 2014).
Species identification and dietary analysis
Host species and diet were identified using metabarcoding. DNA was extracted using the QIAamp Fast DNA Stool Mini Kit (QIAGEN, Germany) following standard protocols and the MT-RNR1 (12S) and COX1 (cytochrome c oxidase subunit I) gene segments amplified using 12SV5-F/R primer and COX1 primers, respectively (Hacker et al., 2021; Riaz et al., 2011). PCR reaction conditions followed the methods described in Hacker et al. (Hacker et al., 2021). The resulting library was sequenced on an Illumina NovaSeq platform and 250 bp paired-end reads were generated (Guangdong Magigene Biotechnology Co., Ltd. Guangzhou, China).
We used CLC Genomics Workbench v12.0 (QIAGEN, Denmark) to determine the host species as well as the prey consumed by each carnivore by mapping sequence reads to reference sequences of possible prey downloaded from GenBank and BOLD (Barcode of Life Data Systems) with representative haplotypes compiled into one .fasta file. Raw reads were required to have at least 98% similarity across at least 90% of the reference sequence for mapping (Hacker et al., 2021). Species and prey identification were made based on the reference taxa with the highest number of reads mapped and the fewest mismatches. Samples in which species could not be identified were analyzed to ensure the reference file was not incomplete by using the de nova assembly tool in CLC, then blasting the resulting contig sequence with the nucleotide databases in NCBI (https://blast.ncbi.nlm.nih.gov/Blast.cgi). As an additional precaution, the geographic range of the determined host and prey species was researched using the IUCN Red List (https://www.iucnredlist.org/) to ensure that it overlapped with the study site. For complete methods on data parameters and methods used, see Hacker et al. (Hacker et al., 2021).
Dietary data were summarized by the frequency of occurrence of prey species in scats observed. The “bipartite” R package was used to construct food web networks (Dormann, 2011). Dietary diversity for each carnivore host species was assessed by calculating richness and Shannon’s Index (Shannon and Weaver, 1949). Interspecific dietary niche overlap was evaluated using Pianka’s Index (Ojk) (value = 0, no dietary overlap and value = 1, complete dietary overlap) and 95% confidence intervals were obtained by bootstrapping with 1000 resamples via the “spaa” R package (Zhang, 2016). Dietary similarity between any two given carnivore species was assessed by calculating the inversed value of Jaccard’s Index based on binary presence absence data of prey.
Results
Sympatric carnivore identification
Of the 480 scat samples sequenced, those which had no sequencing data, were inconclusive, consisted of non-target species, or host species with low sample sizes (1 Asian badger (Meles leucurus) and 4 upland buzzard (Buteo hemilasius)) were removed (Figure S1). The remaining 404 scat samples were composed of three apex carnivores (49 wolf, 147 snow leopard, 19 Eurasian lynx) and three mesocarnivores (63 Pallas’s cat, 87 red fox, 39 Tibetan fox).
Spatial distribution difference and overlap
A total of 322 camera trap sites were surveyed after relocating infrared cameras that did not capture any target carnivore species. A total of 3 cameras were considered to have failed due to loss. We analyzed data from 319 camera sites and obtained 14,316 independent detections during a total effort of 37,192 effective camera trap days. We recorded wolf in 26 sites, snow leopard in 109 sites, Eurasian lynx in 36 sites, red fox in 92 sites, and Tibetan fox in 34 sites. However, the camera detection rates of Pallas’s cat were too low to analyze species occupancy and daily activity patterns.
According to the findings derived from single-season, single-species occupancy models, the snow leopard demonstrated a notably higher probability of occupancy overall compared to other carnivore species, estimated at 0.437 (Table 1). Conversely, the Eurasian lynx exhibited a lower occupancy probability, estimated at 0.161. Further analysis revealed that the occupancy probabilities of the wolf and Eurasian lynx declined with increasing Normalized Difference Vegetation Index (NDVI) (Table 2, Figure 2). Additionally, wolf occupancy probability displayed a negative relationship with roughness index and a positive relationship with prey availability. Snow leopard occupancy probabilities exhibited a negative relationship with distance to roads and NDVI. In contrast, both red fox and Tibetan fox demonstrated a positive relationship with distance to roads. Moreover, red fox occupancy probability increased with higher human disturbance and greater prey availability.
![](https://prod--epp.elifesciences.org/iiif/2/90559%2Fv2%2Fcontent%2F553226v2_fig2.tif/full/max/0/default.jpg)
The relationship between carnivore species occupancy probability and covariates by the optimal models (ΔAIC≤2). The solid line represents the fitted polynomial regression and the gray area represents 95% confidence intervals.
![](https://prod--epp.elifesciences.org/iiif/2/90559%2Fv2%2Fcontent%2F553226v2_tbl1.tif/full/max/0/default.jpg)
![](https://prod--epp.elifesciences.org/iiif/2/90559%2Fv2%2Fcontent%2F553226v2_tbl1a.tif/full/max/0/default.jpg)
Summary of occupancy rate and detection probability of different species for the optimal models (ΔAIC≤2).
![](https://prod--epp.elifesciences.org/iiif/2/90559%2Fv2%2Fcontent%2F553226v2_tbl2.tif/full/max/0/default.jpg)
Covariates influencing carnivore occupancy rate and detection probability based on the optimal models (ΔAIC≤2).
The detection probabilities of wolf, snow leopard, red fox, and Tibetan fox exhibited an increase with elevation (Table 2). Moreover, there was a positive relationship between the detection probability of Tibetan fox and prey availability. The detection probabilities of snow leopard and Eurasian lynx declined as human disturbance increased.
The Sørensen similarity index (S) ranged from 0.1 to 0.5 (Table 3). Compared with other combinations of apex-mesocarnivore pairs, snow leopard and red fox (S=0.477) had relatively high spatial overlap, while Eurasian lynx and Tibetan fox (S=0.198) had the lowest spatial overlap. Moreover, spatial overlap of apex versus apex carnivores and mesocarnivores versus mesocarnivores was relatively low.
![](https://prod--epp.elifesciences.org/iiif/2/90559%2Fv2%2Fcontent%2F553226v2_tbl3.tif/full/max/0/default.jpg)
Spatial overlap (Sørensen’s index), diel activity overlap (Δ) and dietary overlap (Pianka’s index), as well as confidence intervals for carnivore species.
Daily activity patterns and differences
A total of 1444 independent records were obtained for five carnivore species, consisting of 79 records of wolf, 458 records of snow leopard, 126 records of Eurasian lynx, 421 records of red fox, and 359 records of Tibetan fox. Among apex carnivores (Table 3, Figure 3), the daily activity was similar between snow leopard and Eurasian lynx and their diel activity overlap was close to 1 (Δ4=0.900, P=0.285), their daily activity peak was at 21:00 hr and dawn. However, the wolf had a significantly different daily activity pattern with snow leopard (Δ4=0.676, P<0.001) and Eurasian lynx (Δ4=0.661, P<0.001), and its daily activity peaked happened around 9:00 and 18:00 hr. Tibetan fox and red fox had different activity patterns and peaks (P<0.001). The activity peak for red fox activity peaked at 3:00 and 21:00 hr, while the Tibetan fox had a prolonged active bout between noon and dusk. Temporal activity patterns between apex carnivores and mesocarnivores were significantly different, except for wolf and Tibetan fox (Δ4=855, P=0.118).
![](https://prod--epp.elifesciences.org/iiif/2/90559%2Fv2%2Fcontent%2F553226v2_fig3.tif/full/max/0/default.jpg)
Daily activity patterns of carnivore species. Shades of gray indicate temporal areas of overlap.
Dietary composition, diversity and similarity
A total of 26 unique prey species were identified from 9 taxonomic orders (Figure 4, Table S1). Artiodactyla and lagomorpha were the most frequently detected in the diets of apex carnivores and mesocarnivores, accounting for 32.81% and 70.18% of prey counts, respectively. Blue sheep made up 26.50% of prey counts in apex carnivore diet, while plateau pika made up 67.11% of prey counts in mesocarnivore diet. Livestock were present in 17.98% of apex carnivore diet counts and were present in 4.82% mesocarnivore diet counts.
![](https://prod--epp.elifesciences.org/iiif/2/90559%2Fv2%2Fcontent%2F553226v2_fig4.tif/full/max/0/default.jpg)
The food web of carnivore species (SL-snow leopard, EL-Eurasian lynx, PC-Pallas’s cat, RF-red Fox, TF-Tibetan Fox). The heights of the left bars represent the frequency of occurrence of the taxonomic order of prey species, the middle bars represent the number of samples for each carnivore, and the heights of the right bars represent the frequency of occurrence of prey species in scats. The colors of prey match the taxonomic orders. The connecting line widths represent the prey frequency of occurrence in the diet of each carnivore species.
The dietary niche overlap among all carnivore species can be found in Table 3. Wolf and snow leopard had the highest dietary niche overlap value among apex carnivores (Ojk=0.892). The value of Pianka’s index was generally low between apex carnivores and mesocarnivores, except wolf and red fox (Ojk=0.811). In contrast, observed dietary overlap was greatest among the mesocarnivores, especially Pallas’s cat and Tibetan fox, with a value of 0.997.
Red fox had the greatest richness of prey with a value of 16, while Pallas’s cat and Tibetan fox had the lowest diversity of prey with richness value of 6 (Table S2). Dietary similarity was assessed using inversed Jaccard’s Distances, Table S3). Diets were most similar between wolf and snow leopard with a value of 0.588 and least similar between Pallas’s cat and Tibetan fox with a value of 0.200. All other pairs fell between the values of 0.2 and 0.5.
Discussion
Our study addresses, for the first time, the coexistence patterns of carnivore species present in Qilian Mountain National Park across multiple dimensions of niche partitioning. This work substantially contributes to current understanding of carnivore guilds and offers helpful information for biodiversity conservation at a regional scale. Moreover, our study provides important insights into the potentional mechanisms of niche partitioning among sympatric carnivores, particularly intra- and interspecific relations between apex carnivores and mesocarnivores. Specifically, we found that the overall trend of spatial overlap across carnivores is relatively low, that apex carnivores overlap in time and diet, that mesocarnivores showed a high degree of dietary overlap, and that there was substantial similarity in diel activity patterns between apex carnivores and mesocarnivores. These results indicate that carnivores with similar ecological traits foster co-occurrence by adjusting their daily activity patterns and using differing food resources to minimize competitive interactions.
We found dietary and temporal overlap among apex carnivores, showing that spatial partitioning is responsible for their successful coexistence in this area, consistent with our hypothesis. Wolf and snow leopard had the highest dietary overlap and prey similarity between apex carnivore pairs in this study, showing that the avoidance of space and time plays an important role in their coexistence. Recent evidence proves that habitat preference facilitates the coexistence of wolf and snow leopard (Shrotriya et al., 2022). Their hunting strategies may be impacted by their habitat selection. Solitary snow leopards are more suitable for hiding in habitats with features that favor ambush predators, while wolves hunt in packs (Shrotriya et al., 2022). It is clear that wild ungulates (e.g., blue sheep) constituted the primary proportion of wolf and snow leopard diet, followed by small mammals such as plateau pika, Himalayan marmot (Marmota himalayana), and woolly hare (Lepus oiostolus). In addition, livestock consumption also contributed to the high degree of overlap in their diets (Wang et al., 2014). This supports the optimal foraging theory, in which large predators preferentially select food resources that provide maximum benefit (Brown et al., 1999), but also showed that greater competition for resources is likely to occur between wolf and snow leopard due to their use of the same prey species in this area. This may be especially true in times of habitat stress when resources are poor.
Snow leopard and Eurasian lynx had the highest temporal overlap between apex carnivore pairs in our case, showing that spatial and dietary partitioning facilitate their coexistence. Eurasian lynx is considered an opportunistic predator, and its prey varies among different regions with its primary dietary resource being ungulates and small mammals. For example, Eurasian lynx showed a strong preference for brown hare L. europaeus in Turkey, edible dormice Glis glis in Slovenia and Croatia, and chamois Rupicapra rupicapra or roe deer Capreolus capreolus in Switzerland (Mengüllüoğlu et al., 2018; Krofel et al., 2011; Molinari-Jobin et al., 2007). Varied prey selection may be related to sex, age, population density and season (Mengüllüoğlu et al., 2018; Odden et al., 2006). Our results show that woolly hare make up the majority of the Eurasian lynx diet, followed by blue sheep. Woolly hare is mainly distributed in shrubland, meadow, desert and wetland, while blue sheep tend to choose highly sheltered areas, close to bare rocks and cliffs as habitat. Therefore, prey preferences among snow leopard and Eurasian lynx also contributes to spatial avoidance.
Mesocarnivores had substantial overlap in diet, underscoring food resources as a primary competitive factor that necessitates spatial and temporal partitioning for successful coexistence. This finding aligns with previous research indicating that mesocarnivores use temporal and spatial segregation to reduce competition and the probability of antagonistic interspecific encounters (Ferreiro-Arias et al., 2021; Li et al., 2022a). Differences in habitat preference may lead to spatial niche partitioning among mesocarnivores (Wang et al., 2022). In addition, species can adjust temporal periods of behavior to respond to environmental change, competition, and predation risk (Gallo et al., 2022; Finnegan et al., 2021b; van der Vinne et al., 2019). Pallas’s cat hunts during crepuscular and diurnal periods and inhabit meadow habitat with greater prey abundance (Anile et al., 2021; Greco et al., 2022; Ross et al., 2019). In contrast, red fox is primarily nocturnal and occupies diverse habitats depending on prey abundance (Gołdyn et al., 2003; Pandolfi et al., 1997; Reshamwala et al., 2022). Tibetan fox is a diurnal hunter of the Tibetan plateau, preferring shrub meadow, meadow steppe and alpine meadow steppe (Gong and Hu, 2003). It is worth noting the substantial overlap in diet between Pallas’s cat and Tibetan fox. The dietary overlap between the two was extremely high, with a Pianka’s value close to 1. Dietary analyses showed that pika contributed to more than 85% of their collective diets, with 90% of Pallas’s cat diet being pika. Pika may be an optimal prey item in the area because of size and year-round activity (Ross, 2009).
Previous studies have shown that the Pallas’s cat and Tibetan fox are specialist predators of pikas (Harris et al., 2014; Ross, 2009). However, specialization on pika is facultative in that Pallas’s cat and Tibetan fox can select other prey items when pika availability is low (Harris et al., 2014; Ross, 2009). This was observed in our study, even though dietary diversity was low.
Apex carnivores and mesocarnivores exhibited considerable overlap in time overall, showing that spatial and dietary partitioning may play a large role in facilitating their coexistence. As confirmed by previous research, kit foxes (V. macrotis) successfully coexisted with dominant carnivores by a combination of spatial avoidance and use of alternative resources (Lonsinger et al., 2017). Differences in body mass may play a crucial role in minimizing dietary overlap, effectively reducing interspecific competition between apex and mesocarnivores. Of exception in our study, however, was wolf and red fox, who exhibited more dietary overlap, showing that temporal and spatial avoidance may promote their coexistence. As canid generalist-opportunist species, the wolf and red fox consumed similar prey, albeit the red fox may have obtained livestock and ungulate species via scavenging or by preying on very young individuals (Hacker et al., 2022). In addition, the wolf and red fox had different peak activity times, suggesting temporal segregation as a potentially strong driver of coexistence. Our occupancy model also revealed that mesocarnivores prefer habitats that coincide with human disturbance, possible contributing to spatial niche differentiation with other apex carnivore species like the wolf. Our findings echo other recent research which supports the idea that red foxes can coexist with wolves by exploiting a broader ecological niche (Shrotriya et al., 2022).
Our study highlights the effectiveness of combining camera trapping with DNA metabarcoding for detecting and identifying both cryptic and rare species within a sympatric carnivore guild. This integrated approach allowed us to capture a more comprehensive view of species presence and interactions compared to traditional visual surveys. whereas, it is important to acknowledge the challenges associated with this technique, including the high costs of equipment and the need for specialized training and computational resources to manage and analyze the large volumes of sequence data. Despite these challenges, the benefits of this combined method in improving biodiversity assessments and understanding species coexistence outweigh the drawbacks. However, several restrictions remain for this research. The first limitation involves differences in samples sizes. Although the scat samples of Tibetan fox were relatively low, the accuracy of DNA metabarcoding in informing species presence in diet ensures that data are informative and thus still important for species conservation management decisions (Hacker et al., 2022). Second, the methodology of foraging (e.g., predation or scavenging) and the condition of the prey item (e.g., age or size) cannot be identified in dietary studies (Hacker et al., 2022). Pika is a prime component of diet among mesocarnivores, especially in the diet of the Pallas’s cat and the Tibetan fox. We surmise that the simultaneous dependence on pika led to partial overlap in spatial and activity patterns, resulting in increased potential competitive interactions. Due to the lack of spatial and temporal analysis of Pallas’s cat in our study, further monitoring is needed to develop a comprehensive conservation plan. Despite these limitations, our study provides a foundation from which future studies interested in niche partitioning among carnivores along spatial, temporal, and dietary dimensions, can be modeled.
In summary, our study has shown that the coexistence of carnivore species in the landscapes of Qilian Mountain National Park can be facilitated along three niche axes, with spatial segregation being the most pronounced of the three. Apex carnivore species tended to overlap temporally or trophically, mesocarnivore species had high dietary overlap with each other, and apex carnivore and mesocarnivore species displayed similarity in temporal use. Pika, blue sheep and livestock were found to make up a large proportion of carnivore diet. Resource competition between wolf and snow leopard and the interspecific competition between Pallas’s cat and Tibetan fox were strong in this area. Based on the findings presented above, we recommend targeted efforts for enhanced protection and management in several key areas. First, greater efforts are needed to protect habitat, including the establishment of habitat corridors and optimizing grassland fence layouts to safeguard migration routes. Second, resource competition should be carefully monitored between snow leopards and wolves, as well as between Pallas’s cat and Tibetan foxes. More attention is needed for pika at our study site considering the role they play in the conservation of Pallas’s cat and Tibetan fox populations. Previous poisoning campaigns targeting pikas were initiated out of concern for grassland degradation (Smith and Foggin, 1999). Recent research underscores the importance of pika population health and habitat for the distribution of Pallas’s cats (Greenspan and Giordano, 2021). It’s crucial to note that poisoning campaigns targeting small mammals may be incredibly dangerous for mesocarnivores feeding on them due to secondary poisoning, necessitating caution in such conservation strategies. Further, vigilance regarding the preponderance of smaller prey in predator diets is vital, as this may indicate severe loss of larger prey, which will increase the risk of interference competition (Steinmetz et al., 2021). Third, efforts should focus on restoring vulnerable wild prey populations, strengthening grazing area management, and supporting livelihoods of herders to mitigate livestock predation and consequent human-wildlife conflicts. Regular patrols of protected areas and core habitats are essential, along with educational outreach to conservation staff and herders. Our study corroborates and complements the findings of prior studies on these species and their coexistence mechanisms, offering insights crucial for wildlife conservation in the region.
Acknowledgements
We would like to thank Mr. Jiong He, Yayue Gao, Duifang Ma, Liji Wu, Dazhi Hu and other colleagues of Qilianshan National Park for their generous assistance in the field surveys. We thank Dr. Charlotte Hacker for editing the English text of this manuscript. This work was supported by the Welfare Project of the National Scientific Research Institution (No. CAFYBB2019ZE003 and CAFYBB2018ZD001) and the National Natural Science Foundation of China (No. 32201430 and No. 32101409).
Declaration of interests
The authors declare no competing interests.
References
- 1.On the high trail: examining determinants of site use by the Endangered snow leopard Panthera uncia in Qilianshan, ChinaOryx 50:231–238https://doi.org/10.1017/S0030605315001027
- 2.Density and activity patterns of Pallas’s cats, Otocolobus manul, in central MongoliaWildl Res 48https://doi.org/10.1071/WR20061
- 3.Camera trapping for wildlife research
- 4.The Ecology of Fear: Optimal Foraging, Game Theory, and Trophic InteractionsJournal of Mammalogy 80:385–399https://doi.org/10.2307/1383287
- 5.Coexistence mechanisms of small carnivores in a near-pristine area within the mountains of Southwest ChinaGlobal Ecology and Conservation 49https://doi.org/10.1016/j.gecco.2023.e02777
- 6.Advances in species coexistence theoryBiodiversity Science 25:345–354https://doi.org/10.17520/biods.2017034
- 7.Ecological correlates of the spatial co-occurrence of sympatric mammalian carnivores worldwideEcol Lett 21:1401–1412https://doi.org/10.1111/ele.13124
- 8.Competition and coexistence in a small carnivore guildOecologia 184:873–884https://doi.org/10.1007/s00442-017-3916-2
- 9.Counting with DNA in metabarcoding studies: How should we convert sequence reads to dietary data?Mol Ecol 28:391–406https://doi.org/10.1111/mec.14734
- 10.Environmental DNA metabarcoding: Transforming how we survey animal and plant communitiesMol Ecol 26:5872–5895https://doi.org/10.1111/mec.14350
- 11.Niche partitioning and species coexistence in a Neotropical felid assemblageActa Oecologica 36:403–412https://doi.org/10.1016/j.actao.2010.04.001
- 12.Global priorities for national carnivore conservation under land use changeSci Rep 6https://doi.org/10.1038/srep23814
- 13.Diet, Morphology, and Interspecific Killing in CarnivoraThe American Naturalist 167:524–536https://doi.org/10.1086/501033
- 14.How to be a specialist? Quantifying specialisation in pollination networksNetwork Biology 1:1–20
- 15.Fine-scale coexistence between Mediterranean mesocarnivores is mediated by spatial, temporal, and trophic resource partitioningEcol Evol 11:15520–15533https://doi.org/10.1002/ece3.8077
- 16.“When the felid’s away, the mesocarnivores play”: seasonal temporal segregation in a neotropical carnivore guildMamm Biol 101:631–638https://doi.org/10.1007/s42991-021-00110-9
- 17.“When the felid’s away, the mesocarnivores play”: seasonal temporal segregation in a neotropical carnivore guildMamm Biol 101:631–638https://doi.org/10.1007/s42991-021-00110-9
- 18.. unmarked: An R Package for Fitting Hierarchical Models of Wildlife Occurrence and AbundanceJournal of Statistical Software 43:1–23
- 19.Investigating animal activity patterns and temporal niche partitioning using camera trap data: Challenges and OpportunitiesRemote Sensing in Ecology and Conservation 3https://doi.org/10.1002/rse2.60
- 20.Mammals adjust diel activity across gradients of urbanizationeLife 11https://doi.org/10.7554/eLife.74756
- 21.Spatial segregation between Iberian lynx and other carnivoresAnim Biodiv Conserv :347–354https://doi.org/10.32800/abc.2019.42.0347
- 22.The Struggle for Existence
- 23.Habitat use and diet of the red fox Vulpes vulpes in an agricultural landscape in PolandZeitschrift für Jagdwissenschafte 49:191–200https://doi.org/10.1007/BF02189737
- 24.Feeding interactions in an assemblage of terrestrial carnivores in central MexicoZoological Studies 54https://doi.org/10.1186/s40555-014-0102-7
- 25.The Summer Microhabitat Selection of Tibetan Fox in the Northwest Plateau of SichuanActa Theriologica Sinica 23:266–269
- 26.Human density modulates spatial associations among tropical forest terrestrial mammal speciesGlobal Change Biology 28:7205–7216https://doi.org/10.1111/gcb.16434
- 27.SpatioLJtemporal occurrence and sensitivity to livestock husbandry of Pallas’s cat in the Mongolian AltaiJ Wildl Manag 86https://doi.org/10.1002/jwmg.22150
- 28.A rangewide distribution model for the Pallas’s cat (Otocolobus manul): identifying potential new survey regions for an understudied small catMammalia 85:574–587https://doi.org/10.1515/mammalia-2020-0094
- 29.Dietary diversity and niche-partitioning of carnivores across the Qinghai-Tibetan Plateau of China using DNA metabarcodingJournal of Mammalogy 103:1005–1018https://doi.org/10.1093/jmammal/gyac044
- 30.Regional comparison of snow leopard (panthera uncia) diet using dna metabarcodingBiodiversity and Conservation 30:797–817https://doi.org/10.1007/s10531-021-02118-6
- 31.The Competitive Exclusion PrincipleScience 131:1292–1297
- 32.Evidence that the Tibetan fox is an obligate predator of the plateau pika: Conservation implicationsJournal of Mammalogy 95:0–0https://doi.org/10.1644/14-MAMM-A-021
- 33.Fear, foraging and olfaction: how mesopredators avoid costly interactions with apex predatorsOecologia 187:573–583https://doi.org/10.1007/s00442-018-4133-3
- 34.Population monitoring of snow leopards using noninvasive collection of scat samples: A pilot studyAnimal Conservation 11:401–411https://doi.org/10.1111/j.1469-1795.2008.00195.x
- 35.Spatio-temporal interactions facilitate large carnivore sympatry across a resource gradientProc R Soc B 284https://doi.org/10.1098/rspb.2016.1860
- 36.DNA metabarcoding illuminates dietary niche partitioning by African large herbivoresProc Natl Acad Sci U S A 112:8019–8024https://doi.org/10.1073/pnas.1503283112
- 37.An empirical evaluation of camera trap study design: How many, how long and when?Methods Ecol Evol 11:700–713https://doi.org/10.1111/2041-210X.13370
- 38.Diet of Eurasian lynx Lynx lynx in the northern Dinaric Mountains (Slovenia and Croatia)Acta Theriologica 56:315–322https://doi.org/10.1007/s13364-011-0032-2
- 39.Paws without claws? Ecological effects of large carnivores in anthropogenic landscapesProc R Soc B 283https://doi.org/10.1098/rspb.2016.1625
- 40.Non-linear relationship between body size of terrestrial carnivores and their trophic niche breadth and overlapBasic and Applied Ecology 38:36–46https://doi.org/10.1016/j.baae.2019.06.004
- 41.Spatial co-occurrence and temporal activity patterns of sympatric mesocarnivores guild in Qinling MountainsGlobal Ecology and Conservation 36https://doi.org/10.1016/j.gecco.2022.e02129
- 42.Projected impacts of climate change on snow leopard habitat in Qinghai Province, ChinaEcol Evol 11:17202–17218https://doi.org/10.1002/ece3.8358
- 43.Temporal and Spatial Activity Patterns of Sympatric Wild Ungulates in Qinling Mountains, ChinaAnimals 12https://doi.org/10.3390/ani12131666
- 44.Coexistence of two sympatric flagship carnivores in the human-dominated forest landscapes of Northeast AsiaLandscape Ecology 34:291–305https://doi.org/10.1007/s10980-018-0759-0
- 45.Interference interactions, coLJexistence and conservation of mammalian carnivoresDiversity and Distributions 6:169–176https://doi.org/10.1046/j.1472-4642.2000.00069.x
- 46.The roles of habitat and intraguild predation by coyotes on the spatial dynamics of kit foxesEcosphere 8https://doi.org/10.1002/ecs2.1749
- 47.Camera-trapping survey of the mammal diversity in the Qilian Mountains National Nature Reserve, Gansu ProvinceACTA Theriologica Sinica 41:90–98https://doi.org/10.16829/j.slxb.150418
- 48.Occupancy estimation and modeling: inferring patterns and dynamics of species occurrence, Second edition. ed
- 49.Human disturbance increases trophic niche overlap in terrestrial carnivore communitiesProc Natl Acad Sci U S A 117:26842–26848https://doi.org/10.1073/pnas.2012774117
- 50.Foraging ecology of Eurasian lynx populations in southwest Asia: Conservation implications for a diet specialistEcology and Evolution 8:9451–9463https://doi.org/10.1002/ece3.4439
- 51.Variation in diet, prey selectivity and home-range size of Eurasian lynx Lynx lynx in SwitzerlandWildlife Biology 13:393–405
- 52.Can prey occupancy act as a surrogate for mesopredator occupancy? A case study of ocelot (Leopardus pardalis)Mamm Biol 102:163–175https://doi.org/10.1007/s42991-022-00232-8
- 53.A small neighborhood well-organized: seasonal and daily activity patterns of the community of large and mid-sized mammals around waterholes in the Gobi Desert, MongoliaFront Zool 18https://doi.org/10.1186/s12983-021-00412-1
- 54.Examination of two new technologies to assess the diet of woodland caribou: video recorders attached to collars and DNA barcodingCan J For Res 43:897–900https://doi.org/10.1139/cjfr-2013-0108
- 55.Top predators constrain mesopredator distributionsNat Commun 8https://doi.org/10.1038/ncomms15469
- 56.Diet of Eurasian lynx, Lynx lynx, in the boreal forest of southeastern Norway: The relative importance of livestock and hares at low roe deer densityEuropean Journal of Wildlife Research 52:237–244https://doi.org/10.1007/s10344-006-0052-4
- 57.vegan: Community Ecology Package
- 58.Towards a bestLJpractices guide for camera trapping: assessing differences among camera trap models and settings under field conditionsJournal of Zoology 316:197–208https://doi.org/10.1111/jzo.12945
- 59.Spatial behaviour of the red fox (Vulpes vulpes) in a rural area of central ItalyItalian Journal of Zoology 64:351–358https://doi.org/10.1080/11250009709356222
- 60.On the move: spatial ecology and habitat use of red fox in the Trans-Himalayan cold desertPeerJ 10https://doi.org/10.7717/peerj.13967
- 61.EcoPrimers: Inference of new DNA barcode markers from whole genome sequence analysisNucleic acids research 39https://doi.org/10.1093/nar/gkr732
- 62.Estimating overlap of daily activity patterns from camera trap dataJournal of Agricultural, Biological, and Environmental Statistics 14:322–337https://doi.org/10.1198/jabes.2009.08038
- 63.Status and ecological effects of the world’s largest carnivoresScience 343https://doi.org/10.1126/science.1241484
- 64.Predator interactions, mesopredator release and biodiversity conservationEcol Lett 12:982–998https://doi.org/10.1111/j.1461-0248.2009.01347.x
- 65.Providing an ecological basis for the conservation of the Pallas’s cat (Otocolobus manul). PhD thesis, University of Bristol, UKhttps://doi.org/10.13140/RG.2.1.4338.1521
- 66.The behaviour and ecology of the manul
- 67.Quantifying levels of animal activity using camera-trap dataMethods in Ecology and Evolution 5https://doi.org/10.1111/2041-210X.12278
- 68.Prey availability and temporal partitioning modulate felid coexistence in Neotropical forestsPLoS ONE 14https://doi.org/10.1371/journal.pone.0213671
- 69.Resource Partitioning in Ecological CommunitiesScience
- 70.Temporal partitioning and spatiotemporal avoidance among large carnivores in a human-impacted African landscapePLoS ONE 16https://doi.org/10.1371/journal.pone.0256876
- 71.The Mathematical Theory of CommunicationPhilosophical Review 60
- 72.Prey partitioning and livestock consumption in the world’s richest large carnivore assemblageCurr Biol 31:4887–4897https://doi.org/10.1016/j.cub.2021.08.067
- 73.Feeding Patterns of Three Widespread Carnivores—The Wolf, Snow Leopard, and Red Fox—in the Trans-Himalayan Landscape of IndiaFrontiers in Ecology and Evolution 10https://doi.org/10.3389/fevo.2022.815996
- 74.The plateau pika (Ochotona curzoniae) is a keystone species for biodiversity on the Tibetan plateauAnimal Conservation 2:235–240https://doi.org/10.1017/S1367943099000566
- 75.Fear of the human ‘super predator’ reduces feeding time in large carnivoresProc R Soc B 284https://doi.org/10.1098/rspb.2017.0433
- 76.A method of establishing groups of equal amplitude in plant sociology based on similarity of species content and its application to analyses of the vegetation on Danish commons. Kobenhavn,: Ejnar Minksgaard
- 77.The effects of prey depletion on dietary niches of sympatric apex predators in Southeast AsiaIntegr Zool 16:19–32https://doi.org/10.1111/1749-4877.12461
- 78.Spatial co-occurrence patterns of sympatric large carnivores in a multi-use African systemPLoS ONE 18https://doi.org/10.1371/journal.pone.0280420
- 79.Niche partitioning between sympatric wild canids: the case of the golden jackal (Canis aureus) and the red fox (Vulpes vulpes) in north-eastern ItalyBMC Ecology and Evolution 21https://doi.org/10.1186/s12862-021-01860-3
- 80.Spatio-temporal partitioning facilitates mesocarnivore sympatry in the Stara Planina Mountains, BulgariaZoology (Jena) 141https://doi.org/10.1016/j.zool.2020.125801
- 81.Maximising survival by shifting the daily timing of activityEcol Lett 22:2097–2102https://doi.org/10.1111/ele.13404
- 82.Coexistence of predators in time: Effects of season and prey availability on species activity within a Mediterranean carnivore guildEcol Evol 10:11408–11422https://doi.org/10.1002/ece3.6778
- 83.Spatiotemporal overlap among sympatric Pallas’s cat (Otocolobus manul), Tibetan fox (Vulpes ferrilata) and red fox (V. vulpes) in the source region of the Yangtze RiverBiodiv Sci 30https://doi.org/10.17520/biods.2021365
- 84.Dietary overlap of snow leopard and other carnivores in the Pamirs of Northwestern ChinaChin Sci Bull 59:3162–3168https://doi.org/10.1007/s11434-014-0370-y
- 85.Camera-trapping survey of the mammals and birds in the Qilian Mountain National Park (Qinghai area), ChinaACTA Theriologica Sinica 39:466–475https://doi.org/10.16829/j.slxb.150307
- 86.spaa: SPecies Association Analysis
- 87.The study of Qinghai nature reserve
Article and author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
Copyright
© 2024, Cong et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 171
- downloads
- 15
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.