Two-species occupancy models: a new parameterization applied to co-occurrence of secretive rails

Impacts of coyote colonization on coastal mammalian predators

Extreme ecosystem modification by humans has caused drastic reductions in populations and ranges of top mammalian predators, while simultaneously allowing synanthropic mesopredator species to expand. These conditions often result in inflated local densities of highly adaptable mesopredators that disrupt trophic dynamics and place unsustainable predation pressure on native prey populations. Colonization of a dominant predator may lead to top-down control of mesopredators and restore trophic balance. Coyotes are a novel colonizer of some coastal barrier islands of eastern North America, offering an opportunity to test how the addition of an apex predator impacts an established guild of mesopredators. To assess their trophic impact, we conducted 75,576 camera trapping hours over an 18-month study period, capturing > 1.5 million images across 108 coastal camera sites. Using two-species occupancy and habitat use models, we found sizeable effects of coyote habitat use on that of red foxes and free-ranging domestic cats, suggesting that coyotes function as apex predators in barrier island ecosystems. In fact, the only factor that determined the spatial pattern of highly ubiquitous red foxes was the sympatric habitat use of the largest carnivore in the food web-coyotes. That 'novel' apex predators can become established in coastal food webs illustrates the highly dynamic nature of conservation challenges for habitats and species at the edge of the sea.

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.

Interplay of competition and facilitation in grazing succession by migrant Serengeti herbivores

Competition, facilitation, and predation offer alternative explanations for successional patterns of migratory herbivores. However, these interactions are difficult to measure, leaving uncertainty about the mechanisms underlying body-size-dependent grazing-and even whether succession occurs at all. We used data from an 8-year camera-trap survey, GPS-collared herbivores, and fecal DNA metabarcoding to analyze the timing, arrival order, and interactions among migratory grazers in Serengeti National Park. Temporal grazing succession is characterized by a "push-pull" dynamic: Competitive grazing nudges zebra ahead of co-migrating wildebeest, whereas grass consumption by these large-bodied migrants attracts trailing, small-bodied gazelle that benefit from facilitation. "Natural experiments" involving intense wildfires and rainfall respectively disrupted and strengthened these effects. Our results highlight a balance between facilitative and competitive forces in co-regulating large-scale ungulate migrations.

American marten occupancy and activity patterns at the southern extent of their range in the eastern United States

The relatively recent rediscovery of an American marten (Martes americana) population that was reintroduced over 30 years ago in southern Vermont provides an opportunity to investigate the relative importance of other mesocarnivores, and forest stand (e.g., DBH, downed logs, vertical structure) and habitat variables to their presence on the Green Mountain National Forest. Marten are state-listed as an endangered species in Vermont and occur there at the southern extent of their range in the eastern United States. We collected detection data from camera surveys in 5 km2 units between 2019 and 2021 (December-April; n = 40 units, 238 cameras). We examined activity patterns and applied an occupancy modeling framework to the detection data to assess the relative importance of covariates at unit and camera levels and assess interactions of marten with other mesocarnivores. We did not find any unit-level occupancy models with significant covariates that were better supported than the base model in the single-season unit-level analysis. Distance to the nearest release site was the covariate most supported for detectability at both spatial scales, and marten occupancy at the camera level was positively influenced by the amount of canopy cover. Two species interaction models did not indicate any positive or negative association beyond random with other mesocarnivores and activity patterns among mesocarnivores had substantial overlap. Marten recovery since the time of the reintroduction appears slow, and even 30 years later, the marten distribution is limited and suggests that dispersal is restricted at some level. We recommend a further investigation of the possible impact of other mesocarnivores to juvenile survival or other vital demographic rate (e.g., recruitment) in marten that were not explicitly measured in this study.

Limited spatiotemporal niche partitioning among mesocarnivores in Gorongosa National Park, Mozambique

Competition drives community composition and structure in many ecosystems. Spatial and temporal niche partitioning, in which competing species divide the environment in space or time, are mechanisms that may allow for coexistence among ecologically similar species. Such division of resources may be especially important for carnivores in African savannas, which support diverse carnivore assemblages. We used camera traps to explore patterns of spatial and temporal niche partitioning among four mesocarnivore species in Mozambique's Gorongosa National Park: large-spotted genet (Genetta maculata), African civet (Civettictis civetta), honey badger (Mellivora capensis) and marsh mongoose (Atilax paludinosus). We applied a multispecies occupancy model to evaluate spatial partitioning among mesocarnivores and to quantify the environmental factors that affect species-specific habitat use, including relative lion (Panthera leo) activity. We also analyzed the temporal activity overlap of the four focal species. We identified species-specific habitat covariates that influenced detection probabilities but found no evidence of spatial or temporal partitioning among mesocarnivores in the study system. Indeed, we found some evidence for spatial co-occurrence between two of our focal species: African civet and marsh mongoose. There may be limited competition among mesocarnivores in this system, perhaps due to niche and diet differentiation among these species and an abundance of resources. While we found limited evidence that lion activity impacts mesocarnivores, ongoing monitoring of intraguild interactions is vital as apex predator populations recover in the system. This study adds to a growing understanding of African mesocarnivore ecology and highlights the importance of understanding these dynamics for effective multispecies conservation and restoration.

Habitat partitioning among sympatric tinamous in semiarid woodlands of central Argentina

Sympatric, phylogenetically related and morphologically similar species that overlap in their distributions at a regional scale display different patterns of co-occurrence at local assemblages. Occurrence of each species at local scales might be the result of interspecific competition for limiting resources. However, these patterns could also arise from species-specific habitat preferences along the abiotic or land use gradients. To assess the role of these mechanisms, we investigated niche partitioning among sympatric tinamou species occurring in semiarid woodlands of central Argentina. We used occupancy models incorporating habitat characteristics and interspecific interactions, while accounting for detectability, to examine the spatial overlap among elegant crested tinamous (Eudromia elegans), brushland tinamous (Nothoprocta cinerascens), and nothura tinamous (Nothura darwinii and N. maculosa) across a wide regional scale. In addition, we investigated time partitioning among these species by estimating the degree of overlap in their daily activity patterns. The regional distribution of the three species was influenced by the gradient in plant productivity and vegetation structure, in agreement with their broad habitat requirements. We also found that the occurrence of each species was presumably affected by the presence of one or two predator species. Models including interactions among tinamou species found weak negative and positive interactions among species pairs, suggesting that co-occurrence patterns were mainly driven by species-specific habitat use rather than interspecific competition. The three species were diurnal, showing two main peaks of activity, and overlapped widely in their overall diel activity, although subordinate species tended to shift their activity patterns to reduce encounters with the dominant tinamou species, suggesting some segregation in this niche dimension. Projected changes in rainfall seasonality and warmer conditions in this region could benefit elegant crested tinamous over the other two species, although climate and land use changes will likely have a negative impact on all tinamou species.

The Minimal Impact of Anthropogenic Disturbances on the Spatial Activities of Leopard Cats in Xinlong, China

The habitat plays a crucial role in ensuring the survival of wildlife. However, the increasing disturbances caused by human activities present a substantial threat to habitats, especially for species such as the leopard cat (Prionailurus bengalensis), which is a significant small predator. Currently, research on leopard cats predominantly focuses on low-altitude regions within its distribution range, leaving plateau areas understudied. To enhance our understanding of the impact of human disturbances on leopard cat habitats, we undertook a study employing infrared camera trappings to monitor leopard cats' activity in Xinlong of southwestern China between 2015 and 2023. We analyzed the spatial distribution and habitat suitability of the leopard cats by utilizing ensemble species distribution models (ESDMs). Moreover, we employed two-species occupancy models to investigate the spatial interaction between leopard cats and human disturbances. The results indicated that (1) the potential suitable habitat area for leopard cats encompassed approximately 1324.93 km2 (14.3%), primarily located along the banks of Yalong river. (2) The distribution of suitable habitat was predominantly influenced by competitors, specifically the yellow-throated marten (YTM), accounting for 52.4% of the influence, as well as environmental factors such as distance to water (DTW) at 12.0% and terrain roughness index (TRI) at 10.0%. Human interference, including cattle presence (4.6%), distance to road (DTD, 4.9%), and distance to settlement (DTS, 3.5%), had a limited impact on the habitat distribution. (3) Within a 5 km radius, habitat suitability increased with proximity to human settlements. (4) Leopard cats exhibited spatial independence from humans and domestic cattle (species interaction factor (SIF) = 1.00) while avoiding domestic horses (SIF = 0.76 ± 0.03). The relatively minor impact of human disturbances in Xinlong could be attributed to traditional cultural practices safeguarding wildlife and the leopard cat's environmental adaptability. We recommend establishing a novel conservation paradigm based on the living dynamics of wildlife communities in Xinlong, thereby offering a more targeted approach to biodiversity preservation in the future.

The Effects of Anthropogenic Disturbances on the Spatiotemporal Patterns of Medium-Large Mammals in Tropical Volcanic Landscapes

A comprehensive understanding of the consequences of human interactions with mammals is a critical factor in supporting and conserving species in landscapes dominated by humans, which are increasingly threatened. This study aimed to identify the spatial and temporal interactions between humans and mammals. A non-parametric statistical approach with kernel density was used to detect human-mammal temporal interactions. The species interaction factor (SIF) was applied to calculate the spatial overlap based on the two-species occupancy detection model. The activity patterns of medium mammals were nocturnal, diurnal, and cathemeral. The human-medium mammal pairs with SIF values that were <1 and statistically significant included the human-long-tailed macaque (Macaca fascicularis) pair, the human-leopard cat (Prionailurus bengalensis) pair, and the human-barking deer (Muntiacus muntjac) pair. Based on their SIF values and the high overlap in their activity times, the human-macaque pairings had a high risk of conflict. Barking deer and leopard cats displayed a coexistence with humans via time-sharing activities. Due to temporal niche variations with human activities, the existence of nocturnal mammals was relatively uninterrupted. This study showed that most mammals are able to adapt spatially and temporally to various human activities. Nonetheless, efforts to mitigate human-wildlife conflict must be maintained, particularly in the case of severely endangered species, such as the Sunda pangolin.

Empirical evidence for effects of invasive American Bullfrogs on occurrence of native amphibians and emerging pathogens

Invasive species and emerging infectious diseases are two of the greatest threats to biodiversity. American Bullfrogs (Rana [Lithobates] catesbeiana), which have been introduced to many parts of the world, are often linked with declines in native amphibians via predation and the spread of emerging pathogens such as amphibian chytrid fungus (Batrachochytrium dendrobatidis [Bd]) and ranaviruses. Although many studies have investigated the potential role of bullfrogs in the decline of native amphibians, analyses that account for shared habitat affinities and imperfect detection have found limited support for clear effects. Similarly, the role of bullfrogs in shaping the patch-level distribution of pathogens is unclear. We used eDNA methods to sample 233 sites in the southwestern USA and Sonora, Mexico (2016-2018) to estimate how the presence of bullfrogs affects the occurrence of four native amphibians, Bd, and ranaviruses. Based on two-species, dominant-subordinate occupancy models fitted in a Bayesian context, federally threatened Chiricahua Leopard Frogs (Rana chiricahuensis) and Western Tiger Salamanders (Ambystoma mavortium) were eight times (32% vs. 4%) and two times (36% vs. 18%), respectively, less likely to occur at sites where bullfrogs occurred. Evidence for the negative effects of bullfrogs on Lowland Leopard Frogs (Rana yavapaiensis) and Northern Leopard Frogs (Rana pipiens) was less clear, possibly because of smaller numbers of sites where these native species still occurred and because bullfrogs often occur at lower densities in streams, the primary habitat for Lowland Leopard Frogs. At the community level, Bd was most likely to occur where bullfrogs co-occurred with native amphibians, which could increase the risk to native species. Ranaviruses were estimated to occur at 33% of bullfrog-only sites, 10% of sites where bullfrogs and native amphibians co-occurred, and only 3% of sites where only native amphibians occurred. Of the 85 sites where we did not detect any of the five target amphibian species, we also did not detect Bd or ranaviruses; this suggests other hosts do not drive the distribution of these pathogens in our study area. Our results provide landscape-scale evidence that bullfrogs reduce the occurrence of native amphibians and increase the occurrence of pathogens, information that can clarify risks and aid the prioritization of conservation actions.

Giant pangolin and white-bellied pangolin observations from a World Heritage site

Pangolins are one of the most threatened mammal groups, as a result of habitat loss and exploitation for their meat, scales, and other body parts. However, there is a lack of quantitative data on pangolin populations; their behaviour and ecology make them challenging to survey. We undertook systematic camera-trap surveys of the 5260 km2 World Heritage Dja Faunal Reserve, Cameroon, sampling 305 sites in eight grids over 28,277 camera-trap days. We recorded 768 images of giant pangolin in 99 independent detections at 57 sites (RAI = 0.35), and 2282 images in 355 detections (RAI = 1.26) of white-bellied pangolin at 137 sites. Ground-dwelling giant pangolins were largely confined to the core of the Reserve. Semi-arboreal white-bellied pangolins were predominantly distributed in the northeast, east and south of the Reserve. Lower occupancy in the west and northwest could partly be due to pressures from human settlements around the Hydromekin Dam and Sud-Cameroun Hévéa rubber plantation. Our study suggests that at the ground-level the two species do not spatially segregate, and both were active throughout the night. We found high diel activity overlap, although there was a significant difference in activity peak times. There was also evidence of white-bellied pangolin possibly exhibiting fine-scale behavioural avoidance of giant pangolin.

Spatial co-occurrence patterns of sympatric large carnivores in a multi-use African system

Interspecific interactions can be a key driver of habitat use, and must be accounted for in conservation planning. However, spatial partitioning between African carnivores, and how this varies with scale, remains poorly understood. Furthermore, most studies have taken place within small or highly protected areas, rather than in the heterogeneous, mixed-use landscapes characteristic of much of modern Africa. Here, we provide one of the first empirical investigations into population-level competitive interactions among an African large carnivore guild. We collected detection/non-detection data for an eastern African large carnivore guild in Tanzania's Ruaha-Rungwa conservation landscape, over an area of ~45,000 km2. We then applied conditional co-occupancy models to investigate co-occurrence between lion, leopard, and African wild dog, at two biologically meaningful scales. Co-occurrence patterns of cheetah and spotted hyaena could not be modelled. After accounting for habitat and detection effects, we found some evidence of wild dog avoidance of lion at the home range scale, and strong evidence of fine-scale avoidance. We found no evidence of interspecific exclusion of leopard by lion; rather, positive associations were observed at both scales, suggesting shared habitat preferences. We found little evidence of leopard habitat use being affected by wild dog. Our findings also reveal some interspecific effects on species detection, at both scales. In most cases, habitat use was driven more strongly by other habitat effects, such as biotic resources or anthropogenic pressures, than by interspecific pressures, even where evidence of the latter was present. Overall, our results help shed light on interspecific effects within an assemblage that has rarely been examined at this scale. We also demonstrate the effectiveness of sign-based co-occurrence modelling to describe interspecific spatial patterns of sympatric large carnivores across large scales. We conclude by discussing the implications of our findings for large carnivore conservation in modern African systems.

A comprehensive multi-state conditional occupancy model for evaluating interactions of non-native and native species

A major challenge in ecology is disentangling interactions of non-native, potentially invasive species on native species. Conditional two-species occupancy models examine the effects of dominant species (e.g., non-native) on subordinate species (e.g., native) while considering the possibility that occupancy of one species may affect occupancy and/ or detection of the other. Although conditional two-species models are useful for evaluating the influence of one species on presence of another, it is possible that species interactions are density dependent. Therefore, we developed a novel two-species occupancy model that incorporates multiple abundance states (i.e., absent, present, abundant) of the native species. We showcase the utility of this model with a case study that incorporates random effects and covariates on both occupancy and detection to help disentangle species interactions given varying occupancy and detection in different abundance states. We use snorkel survey data from the Umpqua basin, Oregon, where it is hypothesized that smallmouth bass Micropterus dolomieu, a non-native piscivore, exclude Umpqua chub Oregonichthys kalawatseti, a small endemic minnow. From our two-species multi-state (2SMS) model, we concluded that average occupancy was low for both fishes, and that when non-native bass were present, overall native chub occupancy in the present (0.18 ± 0.05 SD) and abundant (0.19 ± 0.03) states was higher than when non-natives were absent (0.14 ± 0.02/ 0.08 ± 0.02), indicating the non-native was not excluding the native species. By incorporating a species interaction factor, we found a positive association (6.75 ± 5.54 SD) between native chub and non-native bass. The covariates strongly related to occupancy were elevation, algae, and land cover type (urban and shrub). Detection probability for both species (0.21–0.82) was most strongly related to the covariates day of year, water temperature, gravel substrate, and stream order/ magnitude. Incorporation of detection probability and covariates enabled interpretation of interactions between the two species that may have been missed without their inclusion in the modeling process. Our new 2SMS occupancy model can be used by scientists and managers with a broad range of survey and covariate data to disentangle species interactions problems to help them inform management decisions.

Site occupancy by American martens and fishers in temperate deciduous forests of Québec

Interspecific interactions can mediate site occupancy of sympatric species and can be a key factor in habitat use patterns. American martens (Martes americana) and Fishers (Pekania pennanti) are two sympatric mesocarnivores in eastern North American forests. Due to their larger size, fishers have a competitive advantage over martens. We investigated site occupancy of martens and fishers in temperate deciduous forests of Québec, an environment modified by forest management and climate change. We formulated hypotheses on the spatial distribution of the studied species based on the knowledge of local trappers and on the scientific literature regarding forest cover composition, habitat fragmentation, and competitive relationships. We used a network of 49 camera traps monitored over two fall seasons to document site occupancy by both species. We used two-species site occupancy models to assess habitat use and the influence of fishers on martens at spatial grains of different sizes. None of the habitat variables that we considered explained site occupancy by fishers. Availability of dense old coniferous stands explained the spatial distribution of martens both at the home range grain size and at the landscape grain size. We identified the characteristics of habitat hotspots based on the knowledge of trappers, which highlighted the importance of stand composition, height, age, and canopy closure. The characteristics of habitat hotspots for martens in temperate deciduous forests refine the habitat suitability model for American martens that was originally developed for boreal forests of Québec.

Human density modulates spatial associations among tropical forest terrestrial mammal species

The spatial aggregation of species pairs often increases with the ecological similarity of the species involved. However, the way in which environmental conditions and anthropogenic activity affect the relationship between spatial aggregation and ecological similarity remains unknown despite the potential for spatial associations to affect species interactions, ecosystem function, and extinction risk. Given that human disturbance has been shown to both increase and decrease spatial associations among species pairs, ecological similarity may have a role in mediating these patterns. Here, we test the influences of habitat diversity, primary productivity, human population density, and species' ecological similarity based on functional traits (i.e., functional trait similarity) on spatial associations among tropical forest mammals. Large mammals are highly sensitive to anthropogenic change and therefore susceptible to changes in interspecific spatial associations. Using two-species occupancy models and camera trap data, we quantified the spatial overlap of 1216 species pairs from 13 tropical forest protected areas around the world. We found that the association between ecological similarity and interspecific species associations depended on surrounding human density. Specifically, aggregation of ecologically similar species was more than an order of magnitude stronger in landscapes with the highest human density compared to those with the lowest human density, even though all populations occurred within protected areas. Human-induced changes in interspecific spatial associations have been shown to alter top-down control by predators, increase disease transmission and increase local extinction rates. Our results indicate that anthropogenic effects on the distribution of wildlife within protected areas are already occurring and that impacts on species interactions, ecosystem functions, and extinction risk warrant further investigation.

Stability of patch-turnover relationships under equilibrium and nonequilibrium metapopulation dynamics driven by biogeography

Two controversial tenets of metapopulation biology are whether patch quality and the surrounding matrix are more important to turnover (colonisation and extinction) than biogeography (patch area and isolation) and whether factors governing turnover during equilibrium also dominate nonequilibrium dynamics. We tested both tenets using 18 years of surveys for two secretive wetland birds, black and Virginia rails, during (1) a period of equilibrium with stable occupancy and (2) after drought and arrival of West Nile Virus (WNV), which resulted in WNV infections in rails, increased extinction and decreased colonisation probabilities modified by WNV, nonequilibrium dynamics for both species and occupancy decline for black rails. Area (primarily) and isolation (secondarily) drove turnover during both stable and unstable metapopulation dynamics, greatly exceeding the effects of patch quality and matrix conditions. Moreover, slopes between turnover and patch characteristics changed little between equilibrium and nonequilibrium, confirming the overriding influences of biogeographic factors on turnover.

Effects of free-ranging livestock on occurrence and interspecific interactions of a mammalian community

Mammalian communities inhabiting temperate grasslands are of conservation concern globally, especially in Central Asia, where livestock numbers have dramatically increased in recent decades, leading to overgrazing and land-use change. Yet, how this pervasive presence of livestock herds affects the community of wild mammals remains largely unstudied. We used systematic camera trapping at 216 sites across remote, mountainous areas of the Mongolian Altai Mountains to assess the spatial and temporal patterns of occurrence and the interspecific relationships within a mammalian community that includes different categories of livestock. By adopting a recently proposed multispecies occupancy model that incorporates interspecific correlation in occupancy, we found several statistically strong correlations in occupancy among species pairs, with the majority involving livestock. The sign of such associations was markedly species-dependent, with larger wild species of conservation concern, namely, snow leopard and Siberian ibex, avoiding livestock presence. As predicted, we found evidence of a positive correlation in occupancy between predators and their respective main prey. Contrary to our expectations, a number of intraguild species pairs also showed positive co-occurrence, with no evidence of spatiotemporal niche partitioning. Overall, our study suggests that livestock encroaching into protected areas influences the whole local community of wild mammals. Though pastoralism has coexisted with wildlife for millennia in central Asian grasslands, our findings suggest that policies and practices to decrease the pressure of livestock husbandry on wildlife are needed, with special attention on large species, such as the snow leopard and its wild prey, which seem to be particularly sensitive to this pervasive livestock presence.

Landscape use and co-occurrence pattern of snow leopard (Panthera uncia) and its prey species in the fragile ecosystem of Spiti Valley, Himachal Pradesh

The snow leopard (Panthera uncia) plays a vital role in maintaining the integrity of the high mountain ecosystem by regulating prey populations and maintaining plant community structure. Therefore, it is necessary to understand the role of the snow leopard and its interaction with prey species. Further, elucidating landscape use and co-occurrence of snow leopard and its prey species can be used to assess the differential use of habitat, allowing them to coexist. We used camera trapping and sign survey to study the interactions of snow leopard and its prey species (Siberian Ibex- Capra sibrica and Blue sheep-Pseudois nayaur) in the Spiti valley Himachal Pradesh. Using the occupancy modelling, we examined whether these prey and predator species occur together more or less frequently than would be expected by chance. To understand this, we have used ten covariates considering the ecology of the studied species. Our results suggest habitat covariates, such as LULC16 (barren area), LULC10 (grassland), ASP (aspect), SLP (slope) and DW (distance to water), are important drivers of habitat use for the snow leopard as well as its prey species. Furthermore, we found that the snow leopard detection probability was high if the site was used by its prey species, i.e., ibex and blue sheep. Whereas, in the case of the prey species, the probability of detection was low when the predator (snow leopard) was present and detected. Besides this, our results suggested that both species were less likely to detect together than expected if they were independent (Snow leopard—Ibex, Delta = 0.29, and snow leopard—blue sheep, Delta = 0.28, both the values are <1, i.e., avoidance). Moreover, despite the predation pressure, the differential anti-predation habitat selection and restriction of temporal activities by the prey species when snow leopard is present allows them to co-exist. Therefore, considering the strong link between the habitat use by the snow leopard and its prey species, it is imperative to generate quantitative long-term data on predator-prey densities and the population dynamics of its prey species in the landscape.

Asymmetrical intraguild interactions with coyotes, red foxes, and domestic dogs may contribute to competitive exclusion of declining gray foxes

Species coexistence is governed by availability of resources and intraguild interactions including strategies to reduce ecological overlap. Gray foxes are dietary generalist mesopredators expected to benefit from anthropogenic disturbance, but populations have declined across the midwestern USA, including severe local extirpation rates coinciding with high coyote and domestic dog occurrence and low red fox occurrence. We used data from a large‐scale camera trap survey in southern Illinois, USA to quantify intraguild spatial and temporal interactions among the canid guild including domestic dogs. We used a two‐species co‐occurrence model to make pairwise assessments of conditional occupancy and detection rates. We also estimated temporal activity overlap among species and fit a fixed‐effects hierarchical community occupancy model with the four canid species. We partitioned the posterior distributions to compare gray fox occupancy probabilities conditional on estimated state of combinations of other species to assess support for hypothesized interactions. We found no evidence of broadscale avoidance among native canids and conclude that spatial and temporal segregation were limited by ubiquitous human disturbance. Mean guild richness was two canid species at a site and gray fox occupancy was greater when any combination of sympatric canids was also present, setting the stage for competitive exclusion over time. Domestic dogs may amplify competitive interactions by increasing canid guild size to the detriment of gray foxes. Our results suggest that while human activities can benefit some mesopredators, other species such as gray foxes may serve as bellwethers for habitat degradation with trophic downgrading and continued anthropogenic homogenization.

Co-occurrence models fail to infer underlying patterns of avoidance and aggregation when closure is violated

Advances in multi‐species monitoring have prompted an increase in the use of multi‐species occupancy analyses to assess patterns of co‐occurrence among species, even when data were collected at scales likely violating the assumption that sites were closed to changes in the occupancy state for the target species. Violating the closure assumption may lead to erroneous conclusions related to patterns of co‐occurrence among species. Occurrence for two hypothetical species was simulated under patterns of avoidance, aggregation, or independence, when the closure assumption was either met or not. Simulated populations were sampled at two levels (N = 250 or 100 sites) and two scales of temporal resolution for surveys. Sample data were analyzed with conditional two‐species occupancy models, and performance was assessed based on the proportion of simulations recovering the true pattern of co‐occurrence. Estimates of occupancy were unbiased when closure was met, but biased when closure violations occurred; bias increased when sample size was small and encounter histories were collapsed to a large‐scale temporal resolution. When closure was met and patterns of avoidance and aggregation were simulated, conditional two‐species models tended to correctly find support for non‐independence, and estimated species interaction factors (SIF) aligned with predicted values. By contrast, when closure was violated, models tended to incorrectly infer a pattern of independence and power to detect simulated patterns of avoidance or aggregation that decreased with smaller sample size. Results suggest that when the closure assumption is violated, co‐occurrence models often fail to detect underlying patterns of avoidance or aggregation, and incorrectly identify a pattern of independence among species, which could have negative consequences for our understanding of species interactions and conservation efforts. Thus, when closure is violated, inferred patterns of independence from multi‐species occupancy should be interpreted cautiously, and evidence of avoidance or aggregation is likely a conservative estimate of true pattern or interaction.

Humans Are More Influential Than Coyotes on Mammalian Mesopredator Spatiotemporal Activity Across an Urban Gradient

For mammalian mesopredators, human-dominated landscapes offer a mosaic of risk and reward. While the reward of anthropogenic food resources may attract mesopredators to human-dominated areas, increased mesopredator activity and abundance in these areas may cause interspecific conflict. For smaller-bodied mesopredators, the perceived risk of intraguild predation by larger mesopredators may be enough to drive spatial and temporal avoidance strategies to reduce risk while still benefiting from anthropogenic resources. We evaluated how the spatiotemporal activity of four non-canid mammalian mesopredators – raccoon (Procyon lotor), striped skunk (Mephitis mephitis), opossum (Didelphis virginiana), and domestic cat (Felis catus) – might change in the presence of an intraguild predator, the coyote (Canis latrans), and a “super predator,” humans. We quantified mesopredator activity by deploying camera traps at 110 sites across an urban gradient in the Sacramento Metropolitan Area in central California, USA. We hypothesized that mesopredators would likely change their spatiotemporal activity in response to urban intensity (H1), coyotes (H2), human presence (H3), and if urban intensity might mediate the response to humans and coyotes (H4). We used single-species occupancy models to test how mesopredators responded to different spatial scales of urbanization, as well as the temporal presence of coyotes and humans. Top single-species models then informed two-species conditional occupancy models to evaluate how mesopredators responded to “dominant” coyotes. Finally, we used temporal overlap analyzes to evaluate whether activity patterns of mesopredators changed in response to humans and coyotes across three levels of urban intensity. Mesopredators did not change their spatial or temporal activity across the urban gradient when coyotes were present. Additionally, coyotes did not influence mesopredator temporal activity at 1-3-day scales; however, raccoons and cats may avoid coyotes at finer scales. Humans influenced mesopredator detectability, as cats and coyotes avoided humans after 1 and 2 days respectively, and opossums avoided urban areas when humans were present within 3 days. Coyotes may play a limited role in altering smaller-bodied mesopredator activity, especially when humans are present. While the impacts of human presence and urban features are often linked, mesopredators may perceive the risk each poses differently, and adjust their activity accordingly.

Spatio-Temporal Patterns and Source-Dispersion Modeling Towards Sloth Bear–Human Conflict Management in Central India

The impact of humans on biodiversity, in the form of the spatially extensive occurrence of humans and subsequent habitat degradation, leads to negative interactions between humans and native wildlife. However, knowledge of the spatial and temporal interface between humans and wildlife is necessary to understand the root cause of such negative interactions, yet considerably understudied in the context of human-dominated landscapes in south and south-eastern Asia. We took this opportunity, gaining insights on seasonal spatial interaction and spatio-temporal overlap between sloth bears (Melursus ursinus) and humans, and subsequently predicted the conflict source sites and dispersion (i.e., hotspots) based on the robust geographic profiling (GP) method in the Sanjay Tiger Reserve (STR), a human-dominated landscape of central India. Detection data of sloth bear and human were obtained from camera trap survey conducted for two years (2017–2018) and records of conflict incidents (2009–2019) were collected from forest department. We found that sloth bears can co-occur with humans independently of seasons, based on occupancy models. However, during summer, higher temporal overlap (Δ4 = 0.46) and lower spatial overlap (0.31) were observed between sloth bears and humans. Contrastingly, lower temporal overlap (Δ4 = 0.29) and higher spatial overlap (0.44) were observed between the same two during winter. The activity patterns of sloth bears and humans differed significantly across seasons and within the same species in different seasons. Our findings indicated that significant changes in human activity, especially during summer, increased the likelihood of sloth bear-human interaction and subsequent conflict incidents. The mapping of conflict source and dispersion (with high accuracy) also predicted a greater probability of conflict during summer, compared to winter, and thus showed the successful application of GP models in this field. Also, camera trap data alone were able to predict the occurrence of hotspots, demonstrating the use of camera trap records in the successful prediction of source-dispersion of conflict. This study would be useful for decision-makers to alleviate sloth bear–human conflict based on insights on seasonal variation of spatio-temporal overlap between the two and direct conservation efforts accordingly.

Niche partitioning in a guild of invasive mammalian predators

Predators compete aggressively for resources, establishing trophic hierarchies that influence ecosystem structure. Competitive interactions are particularly important in invaded ecosystems where introduced predators can suppress native prey species. We investigated whether niche partitioning exists within a guild of invasive mammalian predators and determined the consequences for native species. Over 4405 camera-trap days, we assessed interactions among three invasive predators: two apex predators (feral cats Felis catus and ferrets Mustela furo) and a mesopredator (stoats Mustela erminea), in relation to their primary prey (lagomorphs, rodents and birds) and habitat use. Further, we tested for mesopredator release by selectively removing cats and ferrets in a pulse perturbation experiment. We found compelling evidence of niche partitioning; spatiotemporal activity of apex predators maximized access to abundant invasive prey, with ferrets targeting lagomorphs and cats targeting rodents. Mesopredators adjusted their behavior to reduce the risk of interference competition, thereby restricting access to abundant prey but increasing predation pressure on diurnal native birds. Stoats were only recorded at the treatment site after both larger predators were removed, becoming the most frequently detected predator at 6 months post-perturbation. We suggest there is spatial and resource partitioning within the invasive predator guild, but that this is incomplete, and avoidance is achieved by temporal partitioning within overlapping areas. Niche partitioning among invasive predators facilitates coexistence, but simultaneously intensifies predation pressure on vulnerable native species.

Development and validation of an eDNA protocol for monitoring endemic Asian spiny frogs in the Himalayan region of Pakistan

Wildlife monitoring programs are instrumental for the assessment of species, habitat status, and for the management of factors affecting them. This is particularly important for species found in freshwater ecosystems, such as amphibians, as they have higher estimated extinction rates than terrestrial species. We developed and validated two species-specific environmental DNA (eDNA) protocols and applied them in the field to detect the Hazara Torrent Frog (Allopaa hazarensis) and Murree Hills Frog (Nanorana vicina). Additionally, we compared eDNA surveys with visual encounter surveys and estimated site occupancy. eDNA surveys resulted in higher occurrence probabilities for both A. hazarensis and N. vicina than for visual encounter surveys. Detection probability using eDNA was greater for both species, particularly for A. hazarensis. The top-ranked detection model for visual encounter surveys included effects of both year and temperature on both species, and the top-ranked occupancy model included effects of elevation and year. The top-ranked detection model for eDNA data was the null model, and the top-ranked occupancy model included effects of elevation, year, and wetland type. To our knowledge, this is the first time an eDNA survey has been used to monitor amphibian species in the Himalayan region.

Carnivores and their prey in Sumatra: Occupancy and activity in human-dominated forests

Understanding the effect of anthropogenic disturbance, and its interaction with carnivores and their prey, is crucial to support the conservation of threatened carnivores, particularly in rapidly changing landscapes. Based on systematic camera-trap sampling of four protected areas in Riau Province of central Sumatra, we assessed the habitat occupancy and spatiotemporal overlap between people, potential carnivore prey, and four threatened species of medium-sized or large carnivores: Sumatran tigers (Panthera tigris sumatrae), Malayan sun bears (Helarctos malayanus), dholes (Cuon alpinus), and Sunda clouded leopards (Neofelis diardi). To assess spatial overlap of target species, we used single-species occupancy models and applied a Species Interaction Factor (SIF) to conditional two-species occupancy models. We also used kernel density estimation (KDE) to assess temporal overlap among these species. Our habitat use models showed that altitude (elevation) strongly influenced the occupancy of all large carnivores and potential prey species. Except for Sunda clouded leopards, the occurrence of large carnivore species was positively related to the spatial co-occurrence of humans (SIF > 1). In addition, we found that sun bears and dholes both exhibited high spatial overlap with tigers, and that sun bears alone exhibited high temporal overlap with people. Our findings contribute to an improved understanding of the contemporary ecology of carnivores and their prey in rapidly changing, southeast Asian landscapes. Such knowledge is important to the conservation and recovery of large carnivores in conservation hotspots that are increasingly dominated by humans across Sumatra, as well as globally.

Assessing carnivore spatial co-occurrence and temporal overlap in the face of human interference in a semiarid forest

Apex predators drive top-down effects in ecosystems and the loss of such species can trigger mesopredator release. This ecological process has been well documented in human-modified small areas, but for management and conservation of ecological communities, it is important to know which human factors affect apex predator occurrence and which mediate mesopredators release at large scales. We hypothesized that mesopredators would avoid spatial and temporal overlap with the apex predator, the puma; but that human perturbations (i.e., cattle raising and trophy hunting) would dampen top-down effects and mediate habitat use. We installed 16 camera traps in each of 45, 10 × 10 km grid cells in the Caldén forest region of central Argentina resulting in 706 total stations covering 61,611 km² . We used single-season occupancy and two-species co-occurrence models and calculated the species interaction factor (SIF) to explore the contributions of habitat, biotic, and anthropic variables in explaining co-occurrence between carnivore pairs. We also used kernel density estimation techniques to analyze temporal overlap in activity patterns of the carnivore guild. We found that puma habitat use increased with abundance of large prey and with proximity to protected areas. Geoffroy's cats and skunks spatially avoided pumas and this effect was strong and mediated by distance to protected areas and game reserves, but pumas did not influence pampas fox and pampas cat space use. At medium and low levels of puma occupancy, we found evidence of spatial avoidance between three pairs of mesocarnivores. All predators were mostly nocturnal and crepuscular across seasons and mesopredators showed little consistent evidence of changing activity patterns with varying levels of puma occupancy or human interference. We found potential for mesopredator release at large scale, especially on the spatial niche axis. Our results suggest that a combination of interacting factors, in conjunction with habitat features and intervening human activities, may make mesopredator release unlikely or difficult to discern at broad scales. Overall, we believe that promoting the creation of new protected areas linked by small forest patches would likely lead to increased predator and prey abundances, as well as the interactions among carnivores inside and outside of protected areas.

Camera trap reveals the co-occurrence patterns of two sympatric muntjac species in southern Anhui Province, China: No spatial segregation

The competitive relationship and coexistence pattern among close related species have long been one of the hot issues in ecological research. Interspecies interactions can exert important influences on the local distribution of rare species. Black muntjac Muntiacus crinifrons is an endemic species to eastern China, currently restricted to limited regions. In contrast, Chinese muntjac Muntiacus reevesi is the most common and widespread deer in southern China. Both species co-occur in southern Anhui and western Zhejiang Province. Little is known about the interaction of these two sympatric-related species. In this study, to investigate the site use determinants and co-occurrence pattern of the two sympatric muntjac species, we conducted a camera trap survey across about 250 km2 in mountainous area of southern Anhui Province, China. We adopted a multistep approach to incorporate habitat preferences while modeling occupancy and detection. We found that the two species did not separate along elevation gradient (range from 400 m to 1,400 m) as described in previous studies. Results of single-species occupancy models indicated that elevation had positive effects on the site use of both species, while slope had an opposite influence on their site use. Positive effects of elevation on the site use implied that both species try to avoid human interference at low elevations. Significant negative effect of slope on the site use of black muntjac suggested that the species prefer habitat with gentle slope and avoided steep. Co-occurrence models and species interaction factors provided evidence that the two muntjac species had an independent occupancy (ψ BM CM = ψ BM cm, SIF = 1) and exhibited a positive species interaction in detection probability (p BM < r BM CM). Combined with the results of previous studies, we suggested that it was fine differentiation in microhabitats and food resources utilization rather spatial or temporal segregation that allowed the two species co-occurrence. The site use determinants revealed in our study would be useful for the habitat conservation and restoration for the rare black muntjac, and the co-occurrence pattern of the two sympatric muntjac species could provide useful information for deep understanding of the coexistence mechanism among forest-dwelling ungulates.

Perceived risk structures the space use of competing carnivores

Competition structures ecological communities. In carnivorans, competitive interactions are disproportionately costly to subordinate carnivores who must account for the risk of interspecific killing when foraging. Accordingly, missed opportunity costs for meso-carnivores imposed by risk can benefit the smallest-bodied competitors. However, the extent to which the risk perpetuates into spatial partitioning in hierarchically structured communities remains unknown. To determine how risk-avoidance behaviors shape the space-use of carnivore communities, we studied a simple community of carnivores in northern Patagonia, Argentina: pumas (Puma concolor; an apex carnivore), culpeo foxes (Lycalopex culpaeus; a meso-carnivore), and chilla foxes (Lycalopex griseus; a small carnivore). We used multi-species occupancy models to quantify the space use within the carnivore community and giving-up densities to understand the behaviors that structure space use. Notably, we applied an analytical framework that tests whether the actual or perceived risk of predation most strongly influences the space use of subordinate carnivores although accounting for their foraging and vigilance behaviors. We found that there was a dominance hierarchy from the apex carnivore through the meso-carnivore to the subordinate small carnivore, which was reflected in space. Although both meso- and small carnivores exhibited similar predator avoidance behavioral responses to apex carnivores, the habitat associations of apex carnivores only altered meso-carnivore space use. The biases in risk management we observed for meso-carnivores likely translates into stable co-existence of this community of competing carnivores. We believe our analytical framework can be extended to other communities to quantify the spatial-behavioral tradeoffs of risk.

Detection criteria and post-field sample processing influence results and cost efficiency of occupancy-based monitoring

Optimization of occupancy-based monitoring has focused on balancing the number of sites and surveys to minimize field efforts and costs. When survey techniques require post-field processing of samples to confirm species detections, there may be opportunities to further improve efficiency. We used scat-based noninvasive genetic sampling for kit foxes (Vulpes macrotis) in Utah, USA, as a model system to assess post-field data processing strategies, evaluate the impacts of these strategies on estimates of occupancy and associations between parameters and predictors, and identify the most cost-effective approach. We identified scats with three criteria that varied in costs and reliability: (1) field-based identification (expert opinion), (2) statistical-based morphological identification, and (3) genetic-based identification (mitochondrial DNA). We also considered four novel post-field sample processing strategies that integrated statistical and genetic identifications to reduce costly genetic procedures, including (4) a combined statistical-genetic identification, (5) a genetic removal design, (6) a within-survey conditional-replicate design, and (7) a single-genetic-replicate with false-positive modeling design. We considered results based on genetic identification as the best approximation of truth and used this to evaluate the performance of alternatives. Field-based and statistical-based criteria prone to misidentification produced estimates of occupancy that were biased high (˜1.8 and 2.1 times higher than estimates without misidentifications, respectively). These criteria failed to recover associations between parameters and predictors consistent with genetic identification. The genetic removal design performed poorly, with limited detections leading to estimates that were biased high with poor precision and patterns inconsistent with genetic identification. Both statistical-genetic identification and the conditional-replicate design produced occupancy estimates comparable to genetic identification, while recovering the same model structure and associations at cost reductions of 67% and 74%, respectively. The false-positive design had the lowest cost (88% reduction) and recovered patterns consistent with genetic identification but had occupancy estimates that were ˜32% lower than estimated occupancy based on genetic identification. Our results demonstrate that careful consideration of detection criteria and post-field data processing can reduce costs without significantly altering resulting inferences. Combined with earlier guidance on sampling designs for occupancy modeling, these findings can aid managers in optimizing occupancy-based monitoring.

Habitat-dependent effects of predatory spiders on prey frogs in a Neotropical wet forest

In seasonal wet Neotropical forests, many studies have suggested that species-rich terrestrial frog assemblages are regulated bottom-up by the abundance of leaf litter. However, terrestrial frogs are prey to a diverse community of predators, and no studies have tested for top-down effects of predators on this or other anuran assemblages. Here, we used an extensive field dataset to model the relative contribution of food resources, microhabitat resources and predators towards the occupancy and detection of two frog species (Craugastor bransfordii and Oophaga pumilio) at La Selva, Costa Rica. Frog occupancy was most strongly influenced by predatory spiders and secondarily influenced by the abundance of leaf litter. Predators exerted stronger effects on frogs than food resources, and frogs avoided predators more as leaf litter decreased. Detection probability was elevated when predators were present. We found support for bottom-up effects of leaf litter on the terrestrial frog assemblage, but top-down effects by predators exerted stronger effects on frog occupancy and detection. Because predator avoidance varied along a resource gradient, predator and resource effects appear to be dependent, supporting interactions between top-down and bottom-up mechanisms. Climate-driven decreases in leaf litter may drive decreased availability of frog refugia and increased interactions between frogs and predators.

A note on investigating co-occurrence patterns and dynamics for many species, with imperfect detection and a log-linear modeling parameterization

Patterns in, and the underlying dynamics of, species co-occurrence is of interest in many ecological applications. Unaccounted for, imperfect detection of the species can lead to misleading inferences about the nature and magnitude of any interaction. A range of different parameterizations have been published that could be used with the same fundamental modeling framework that accounts for imperfect detection, although each parameterization has different advantages and disadvantages.We propose a parameterization based on log-linear modeling that does not require a species hierarchy to be defined (in terms of dominance) and enables a numerically robust approach for estimating covariate effects.Conceptually, the parameterization is equivalent to using the presence of species in the current, or a previous, time period as predictor variables for the current occurrence of other species. This leads to natural, "symmetric," interpretations of parameter estimates.The parameterization can be applied to many species, in either a maximum likelihood or Bayesian estimation framework. We illustrate the method using camera-trapping data collected on three mesocarnivore species in South Texas.