Small Felids Coexist in Mixed-Use Landscape in the Bolivian Amazon

In the face of global species loss, it is paramount to understand the effects of human activity on vulnerable species, particularly in highly diverse, complex systems. The Greater Madidi Landscape in the Bolivian Amazon includes several biodiverse protected areas that were created with the goal of sustaining healthy and diverse ecosystems while not impeding the livelihoods of local indigenous peoples. In this study, we sought to use camera trap data and single-species occupancy analysis to assess the impacts of different forms of human activity on four species of small felids: ocelots (Leopardus pardalis), margays (Leopardus wiedii), jaguarundi (Herpailurus yagouaroundi), and oncilla (Leopardus tigrinus). We modeled both human variables (proximity to indigenous communities, roads, and tourist camps) and non-human variables (terrain ruggedness, proximity to rivers, canopy height, prey availability, and large cat abundance). Margay occupancy was unaffected by any of these human variables and ocelots showed only weak evidence of being affected by tourism. Ocelots were particularly pervasive throughout the study area and were consistently estimated to have high occupancy probability. We did not obtain sufficient data on jaguarundi or oncilla to reliably model these effects. Our results indicate that small cats successfully coexist both with each other and with the surrounding human activity in this unique landscape, which serves as a model for global protected area management.

First Capture of a Jaguar Using a Minimally Invasive Capture System for GPS Tracking in an Isolated Patch of Atlantic Forest in Southern Brazil

This study presents the first successful capture using GPS tagging of a jaguar (Panthera onca) using a minimally invasive capture system (MICS). We used snare-foot traps and a MICS during two capture campaigns in a fragment of Atlantic Forest in southeastern Brazil. The specimen disarmed snares on different occasions, and capture was only possible with the MICS. The captured jaguar, an estimated 16-year-old adult male, was monitored using a GPS Vertex Plus Iridium collar with an optimal performance of 86% in expected locations. The jaguar's home range (659 km2 by MPC and 174 km2 by 95%K) was within the observed range for the species and the animal was primarily maintained in protected areas. The habitat types most frequently used were native grassland (27.2% of 4798 fixes), marsh (24.8%), and dense lowland forest (24.7%). The use of a MICS for trapping jaguars is a promising technique that shows advantages in terms of efficiency, selectivity, portability, reduced potential risk of injury to animals or trappers, and animal stress compared to other capture methods used for the species.

A review of environmental and anthropogenic variables used to model jaguar occurrence

Jaguars (Panthera onca) are a landscape species of conservation importance and our understanding of environmental and anthropogenic drivers of jaguar occurrence is necessary to improve conservation strategies. We reviewed available literature to simply describe environmental and anthropogenic variables used and found to be significant in occurrence modeling. We reviewed 95 documents published from 1980 to 2021 that focused on jaguar occurrence and that used 39 variable types (21 anthropogenic, 18 environmental) among different techniques, scales, and approaches. In general, these variables included both anthropogenic (roads, land use, human activities, and population) and environmental (climate, vegetation, ecological interactions, topographic, water, and others) factors. Twelve variables were identified as affecting jaguar occurrence overall, eleven at local scale and seven at broad scales (regional and continental). Focusing more specifically on the variables that correlate with occurrence should help researchers to make better predictions in areas without quantitative jaguar data.

A multispecies corridor in a fragmented landscape: Evaluating effectiveness and identifying high-priority target areas

While Misiones, Argentina contains one of the largest remnants of Upper Paraná Atlantic Forest ecoregion, one of the world's biodiversity hotspots, only ~50% of this native forest is protected. Each protected area is at risk of becoming an island of native forest surrounded by a matrix of altered habitats due to ongoing land conversion. In an effort to maximize long-term connectivity between existing protected areas, DeMatteo [1] used a multifaceted cost analysis to determine the optimal location for the region's first multispecies corridor using noninvasive data on jaguars (Panthera onca), pumas (Puma concolor), ocelots (Leopardus pardalis), southern tiger cats (Leopardus guttulus), and bush dogs (Speothos venaticus). This work builds on this framework by integrating new field data that broadens the scope of species-specific data across the region's heterogeneous landscape, which varies in vegetation, disturbance, human proximity, and protective status. In addition, two different land use layers are compared across the distributions of the five carnivores, the overlap in their independent distributions, and their relationship to the multispecies corridor. Interpretation of these land use data to species-specific habitat suitability goes beyond DeMatteo [1], with a subdivision of suitability into marginal and optimal areas. This refined scale allows a reanalysis of key areas in the multispecies corridor, where connectivity was previously defined as at highly-at-risk, allowing for a more directed development of management strategies. These analyses and their interpretation extend beyond northern-central Misiones, as the threats are not unique to this region. The need to develop management strategies that balance human-wildlife needs will continue to grow as humans expand their footprint. The techniques applied in this analysis provide a way to identify key areas that require specific management strategies, either through restoration, protection, or a combination of both.

Ocelot, Leopardus pardalis (Mammalia, Carnivora, Felidae), home range in the Lowland Atlantic Forest of Southeastern Brazil

The ocelot is an important Neotropical mesopredator and information on its spatial ecology remains scarce. Here we estimated the ocelot home range in a remnant of Lowland Atlantic Forest in Southeastern Brazil. The data were collected by camera traps installed at eight known ocelot latrines. We estimated the home range both based on the Minimum Convex Polygon (MCP) and the 95% adaptive Kernel density estimator (95%K) to compare with other published studies. We identified 22 ocelots (adult males = 8; adult females = 12; cubs = 2). Six males were recorded at more than one latrine, while all females were recorded at only one sampling point. In addition to male ocelots being recorded at a large number of points, they showed greater intrasexual spatial overlap as they used the same latrines, suggesting larger home ranges than females. The mean home range size for males was 12.1 ± SE 4.4 km2 (range = 6.2 to 20.8 km2) using MCP, and 19.9 ± SE 9.5 km2 (range = 10.1 to 38.9 km2) applying 95%K. Despite our estimates representing an approximation of the total area used by males, both values are consistent with those reported from other locations. Our data complemented the gradient of vegetation type sampled for ocelots in Atlantic Forest and support the suggestion that this environmental variable and, consequently, its effect on prey availability, influence the home range size of ocelot. Information on population ecology and other spatial ecology data are also presented.

Jaguar (Panthera onca) spatial ecology outside protected areas in the Selva Maya

Jaguars (Panthera onca) are endangered in several countries and a priority species for conservation action. Despite extensive research efforts in Mexico most studies have been associated with natural protected areas far from human habitation. Because protected areas are too few to conserve the jaguar population over the long–term, a landscape approach that includes both protected and unprotected lands is needed. This is the case in Quintana Roo State where an ecological corridor linking two protected areas (Yum Balam and Sian Ka'an) is at risk of disappearing due to tourism–driven activities. Between 2013 and 2015, four male jaguars were captured and monitored using satellite telemetry inside the corridor. The mean home range size (± SD) was 101.5 km2 (± 75.9 km2) for the dry season and 172 km2 (± 107.29 km2) for the rainy season. The mean core area size (± SD) was 17.54 km2 (± 16.21 km2) for the dry season and 29.07 km2 (± 16.19 km2) for the rainy season. No significant seasonal differences were found for home ranges or for core areas. As expected, we observed that jaguars preferred forest or young secondary growth over profusely disturbed areas, using whatever vegetation was available in their home ranges. Although it is not protected, a biological corridor linking Yum Balam and Sian Ka'an still holds its own jaguar population, a population that has learned to coexist with human presence. Conservation actions are recommended at landscape level to maintain what remains of tropical mature forest and to promote the development of long–term secondary growth into close tree canopy

Connectivity conservation at the crossroads: protected areas versus payments for ecosystem services in conserving connectivity for Colombian carnivores

Protected areas (PAs) constitute one of the main tools for global landscape conservation. Recently, payments for environmental services (PES) have attracted interest from national and regional governments and are becoming one of the leading conservation policy instruments in tropical countries. However, the degree to which areas designated for PES overlap with areas that are critical for maintaining species' landscape connectivity is rarely evaluated. We estimated habitat distributions and connectivity for 16 of the 22 mammalian carnivores occurring in the Caribbean region of Colombia, and identified the overlap between existing PAs and areas identified as being important for connectivity for these species. We also evaluated the potential impact of creation of new PAs versus new PES areas on conserving connectivity for carnivores. Our results show that PAs cover only a minor percentage of the total area that is important for maintaining connectivity ( x = 26.8 % ± 20.2   s . d . ). On the other hand, PES, if implemented extensively, could contribute substantially to mammalian carnivores' connectivity ( x = 45.4 % ± 12.8   s . d . ). However, in a more realistic scenario with limited conservation investment in which fewer areas are set aside, a strategy based on implementing new PAs seems superior to PES. We argue that prioritizing designation of new PAs will be the most efficient means through which to maintain connectivity.

Multi‐population puma connectivity could restore genomic diversity to at‐risk coastal populations in California

Urbanization is decreasing wildlife habitat and connectivity worldwide, including for apex predators, such as the puma (Puma concolor). Puma populations along California's central and southern coastal habitats have experienced rapid fragmentation from development, leading to calls for demographic and genetic management. To address urgent conservation genomic concerns, we used double‐digest restriction‐site associated DNA (ddRAD) sequencing to analyze 16,285 genome‐wide single‐nucleotide polymorphisms (SNPs) from 401 pumas sampled broadly across the state. Our analyses indicated support for 4–10 geographically nested, broad‐ to fine‐scale genetic clusters. At the broadest scale, the four genetic clusters had high genetic diversity and exhibited low linkage disequilibrium, indicating that pumas have retained genomic diversity statewide. However, multiple lines of evidence indicated substructure, including 10 finer‐scale genetic clusters, some of which exhibited fixed alleles and linkage disequilibrium. Fragmented populations along the Southern Coast and Central Coast had particularly low genetic diversity and strong linkage disequilibrium, indicating genetic drift and close inbreeding. Our results demonstrate that genetically at risk populations are typically nested within a broader‐scale group of interconnected populations that collectively retain high genetic diversity and heterogenous fixations. Thus, extant variation at the broader scale has potential to restore diversity to local populations if management actions can enhance vital gene flow and recombine locally sequestered genetic diversity. These state‐ and genome‐wide results are critically important for science‐based conservation and management practices. Our nested population genomic analysis highlights the information that can be gained from population genomic studies aiming to provide guidance for the conservation of fragmented populations.

Genetic diversity and population structure for ocelots (Leopardus pardalis) in Costa Rica

Ocelots (Leopardus pardalis) are a wide-ranging felid species, occurring from southern United States to northern Argentina. They occupy various habitats and are usually considered the most abundant wild cat species in the Neotropics. However, genetic studies that include free-ranging Mesoamerican ocelots are rare and generally based on small sample sizes. This is the first conservation genetics study on ocelots in Costa Rica and the second one in Mesoamerica that has conducted a genetic assessment of the species at a countrywide scale. We evaluated genetic diversity and population structure of ocelots using 15 microsatellite loci in 28 successfully genotyped individuals from throughout the country. We also compared genetic diversity of Costa Rican ocelots with that of jaguars (Panthera onca) and pumas (Puma concolor) in the country, and with ocelots in Belize. Genetic diversity of ocelots in Costa Rica was relatively high as measured by rarified allelic richness (AR = 5.50 ± 1.36) and expected heterozygosities (HE = 0.79 ± 0.08). We did not detect patterns of genetic substructure, suggesting high levels of gene flow throughout the country and no strong barriers to movement. As expected, genetic diversity of Costa Rican ocelots was higher than co-occurring jaguars and pumas. Additionally, levels of genetic diversity were slightly higher in Costa Rican ocelots when compared with their counterparts in Belize, confirming the south to north decrease in genetic diversity reported in other studies. Our study provides critical baseline information to understand the status of wild ocelot populations in Costa Rica. Future studies on ocelots and other threatened or keystone species should also integrate genetic monitoring and conservation genetics analysis to properly inform management decisions, guarantee their long-term survival, and improve the resilience of ecosystems.

The differential genetic signatures related to climatic landscapes for jaguars and pumas on a continental scale

Modern and paleoclimate changes may have altered species dynamics by shifting species' niche suitability over space and time. We analyze whether the current genetic structure and isolation of the two large American felids, jaguar (Panthera onca) and puma (Puma concolor), are mediated by changes in climatic suitability and connection routes over modern and paleoclimatic landscapes. We estimate species distribution under 5 climatic landscapes (modern, Holocene, last maximum glaciations [LMG], average suitability, and climatic instability) and correlate them with individuals' genetic isolation through causal modeling on a resemblance matrix. Both species exhibit genetic isolation patterns correlated with LMG climatic suitability, suggesting that these areas may have worked as "allele refuges." However, the jaguar showed higher vulnerability to climate changes, responding to modern climatic suitability and connection routes, whereas the puma showed a continuous and gradual transition of genetic variation. Despite differential responsiveness to climate change, both species are subjected to the climatic effects on genetic configuration, which may make them susceptible to future climatic changes, since these are progressing faster and with higher intensity than changes in the paleoclimate. Thus, the effects of climatic changes should be considered in the design of conservation strategies to ensure evolutionary and demographic processes mediated by gene flow for both species.

Connectivity of Protected Areas: Effect of Human Pressure and Subnational Contributions in the Ecoregions of Tropical Andean Countries

Conservationists recognize the value of protected area (PA) systems, with adequate coverage, ecological representation, connection, and management to deliver conservation benefits. Yet, governments primarily focus on coverage, disregarding quantification of the other criteria. While recent studies have assessed global representation and connectivity, they present limitations due to: (1) limited accuracy of the World Database of Protected Areas used, as governments may report areas that do not meet the IUCN or CBD PA definitions or omit subnational PAs, and (2) failure to include human impacts on the landscape in connectivity assessments. We constructed a validated PA database for Tropical Andean Countries (TAC; Bolivia, Colombia, Ecuador, Perú, and Venezuela) and used the existing Protected-Connected-Land (ProtConn) indicator—incorporating the Global Human Footprint as a spatial proxy for human pressure—to evaluate TAC ecoregions’ representation and connectivity. We found that just 27% of ecoregions in the TAC are both protected and connected on more than 17% of their lands. As we included human pressure, we conclude that previous global ProtConn studies overestimate PA connectivity. Subnational PAs are promising for strengthening the representation of PA systems. If nations seek to meet Aichi target 11, or an upcoming post-2020 30% target, further efforts are needed to implement and report subnational conservation areas and appropriately evaluate PA systems.

The boon and bane of boldness: movement syndrome as saviour and sink for population genetic diversity

BACKGROUND

Many felid species are of high conservation concern, and with increasing human disturbance the situation is worsening. Small isolated populations are at risk of genetic impoverishment decreasing within-species biodiversity. Movement is known to be a key behavioural trait that shapes both demographic and genetic dynamics and affects population survival. However, we have limited knowledge on how different manifestations of movement behaviour translate to population processes. In this study, we aimed to 1) understand the potential effects of movement behaviour on the genetic diversity of small felid populations in heterogeneous landscapes, while 2) presenting a simulation tool that can help inform conservation practitioners following, or considering, population management actions targeting the risk of genetic impoverishment.

METHODS

We developed a spatially explicit individual-based population model including neutral genetic markers for felids and applied this to the example of Eurasian lynx. Using a neutral landscape approach, we simulated reintroductions into a three-patch system, comprising two breeding patches separated by a larger patch of differing landscape heterogeneity, and tested for the effects of various behavioural movement syndromes and founder population sizes. We explored a range of movement syndromes by simulating populations with various movement model parametrisations that range from 'shy' to 'bold' movement behaviour.

RESULTS

We find that movement syndromes can lead to a higher loss of genetic diversity and an increase in between population genetic structure for both "bold" and "shy" movement behaviours, depending on landscape conditions, with larger decreases in genetic diversity and larger increases in genetic differentiation associated with bold movement syndromes, where the first colonisers quickly reproduce and subsequently dominate the gene pool. In addition, we underline the fact that a larger founder population can offset the genetic losses associated with subpopulation isolation and gene pool dominance.

CONCLUSIONS

We identified a movement syndrome trade-off for population genetic variation, whereby bold-explorers could be saviours - by connecting populations and promoting panmixia, or sinks - by increasing genetic losses via a 'founder takes all' effect, whereas shy-stayers maintain a more gradual genetic drift due to their more cautious behaviour. Simulations should incorporate movement behaviour to provide better projections of long-term population viability and within-species biodiversity, which includes genetic diversity. Simulations incorporating demographics and genetics have great potential for informing conservation management actions, such as population reintroductions or reinforcements. Here, we present such a simulation tool for solitary felids.

Habitat use of the ocelot (Leopardus pardalis) in Brazilian Amazon

Amazonia forest plays a major role in providing ecosystem services for human and sanctuaries for wildlife. However, ongoing deforestation and habitat fragmentation in the Brazilian Amazon has threatened both. The ocelot is an ecologically important mesopredator and a potential conservation ambassador species, yet there are no previous studies on its habitat preference and spatial patterns in this biome. From 2010 to 2017, twelve sites were surveyed, totaling 899 camera trap stations, the largest known dataset for this species. Using occupancy modeling incorporating spatial autocorrelation, we assessed habitat use for ocelot populations across the Brazilian Amazon. Our results revealed a positive sigmoidal correlation between remote-sensing derived metrics of forest cover, disjunct core area density, elevation, distance to roads, distance to settlements and habitat use, and that habitat use by ocelots was negatively associated with slope and distance to river/lake. These findings shed light on the regional scale habitat use of ocelots and indicate important species-habitat relationships, thus providing valuable information for conservation management and land-use planning.

Improving estimation of puma (Puma concolor) population density: clustered camera-trapping, telemetry data, and generalized spatial mark-resight models

Obtaining reliable population density estimates for pumas (Puma concolor) and other cryptic, wide-ranging large carnivores is challenging. Recent advancements in spatially explicit capture-recapture models have facilitated development of novel survey approaches, such as clustered sampling designs, which can provide reliable density estimation for expansive areas with reduced effort. We applied clustered sampling to camera-traps to detect marked (collared) and unmarked pumas, and used generalized spatial mark-resight (SMR) models to estimate puma population density across 15,314 km² in the southwestern USA. Generalized SMR models outperformed conventional SMR models. Integrating telemetry data from collars on marked pumas with detection data from camera-traps substantially improved density estimates by informing cryptic activity (home range) center transiency and improving estimation of the SMR home range parameter. Modeling sex of unmarked pumas as a partially identifying categorical covariate further improved estimates. Our density estimates (0.84–1.65 puma/100 km²) were generally more precise (CV = 0.24–0.31) than spatially explicit estimates produced from other puma sampling methods, including biopsy darting, scat detection dogs, and regular camera-trapping. This study provides an illustrative example of the effectiveness and flexibility of our combined sampling and analytical approach for reliably estimating density of pumas and other wildlife across geographically expansive areas.

Scraping marking behaviour of the largest Neotropical felids

BACKGROUND

Details of how, why and in what conditions large felids make scrapes is unknown. Here, we examined the general hypothesis about the use of scrapes for marking proposals, as well as to communicate with other individuals to signalize particular points or areas of interest, by studying scrape-marking behaviour of jaguars and pumas.

METHODS

We surveyed by scrapes between five days and two months mainly during dry season in five study areas from Mexico (El Edén and San Ignacio), Belize (Cockscomb) and Brazil (Angatuba and Serra das Almas), which differed in presence and/or abundance of jaguars and pumas. Paths were slowly walked while searching for scrapes by teams normally composed of two people and tracks were stored in GPS, distinguishing the type of path surveyed (unpaved track roads, trails and cross-country).

RESULTS

We found a total of 269 felid scrapes along 467 km of paths surveyed, obtaining a finding rate of 0.576 scrapes per km. Most scrapes were found in car tracks (0.629 scrapes per km), followed by trails (0.581 scrapes per km), and rarely did we find scrapes in cross country (0.094 scrapes per km). In trails, scrapes were found in a similar frequency in the centre and edge, whereas in car tracks they were mainly found in the edge. There were also clear differences in the position of the scrapes between study areas that differed in presence and/or abundance of pumas and jaguars, with scrapes located mainly in the centre in areas only with pumas, in the centre and in the edge in areas with a similar number of jaguars and pumas, and in the edge in area mainly dominated by jaguars. The remarking rate tended to be higher in one of the areas with only pumas where natural vegetation was scarcer. Felids chose sites mainly covered by leaves and located in paths less wide, clean and rarely used.

DISCUSSION

Scraping was a frequent behaviour in the largest felids of America, although in some areas, scraping behaviour was rare. Scrapes seem to be signalizing some specific areas within territories and data suggest that they are made with the proposal of communication between individuals. It seems that a high scraping behaviour in pumas is not related to the presence of jaguars.

Multi-level, multi-scale resource selection functions and resistance surfaces for conservation planning: Pumas as a case study

The importance of examining multiple hierarchical levels when modeling resource use for wildlife has been acknowledged for decades. Multi-level resource selection functions have recently been promoted as a method to synthesize resource use across nested organizational levels into a single predictive surface. Analyzing multiple scales of selection within each hierarchical level further strengthens multi-level resource selection functions. We extend this multi-level, multi-scale framework to modeling resistance for wildlife by combining multi-scale resistance surfaces from two data types, genetic and movement. Resistance estimation has typically been conducted with one of these data types, or compared between the two. However, we contend it is not an either/or issue and that resistance may be better-modeled using a combination of resistance surfaces that represent processes at different hierarchical levels. Resistance surfaces estimated from genetic data characterize temporally broad-scale dispersal and successful breeding over generations, whereas resistance surfaces estimated from movement data represent fine-scale travel and contextualized movement decisions. We used telemetry and genetic data from a long-term study on pumas (Puma concolor) in a highly developed landscape in southern California to develop a multi-level, multi-scale resource selection function and a multi-level, multi-scale resistance surface. We used these multi-level, multi-scale surfaces to identify resource use patches and resistant kernel corridors. Across levels, we found puma avoided urban, agricultural areas, and roads and preferred riparian areas and more rugged terrain. For other landscape features, selection differed among levels, as did the scales of selection for each feature. With these results, we developed a conservation plan for one of the most isolated puma populations in the U.S. Our approach captured a wide spectrum of ecological relationships for a population, resulted in effective conservation planning, and can be readily applied to other wildlife species.