Research Article |
Corresponding author: Miguel Briones-Salas ( miguelbrionessalas@hotmail.com ) Academic editor: Cássio Cardoso Pereira
© 2023 J. Roberto Sosa-López, Nydia Nicté Díaz Bernal, Eugenio Padilla, Miguel Briones-Salas.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Sosa-López JR, Díaz Bernal NN, Padilla E, Briones-Salas M (2023) Analysis of the effects of habitat characteristics, human disturbance and prey on felids presence using long-term community monitoring information. Nature Conservation 53: 279-295. https://doi.org/10.3897/natureconservation.53.104135
|
Predator species are essential for ecosystems as they maintain the ecological integrity of the habitat. Particularly, felids populations have declined globally due to their sensitivity to habitat disturbances. Nevertheless, in Mexico, there are areas protected by indigenous communities to preserve a portion of their territory, benefiting multiple species, including felids. Although the National Commission of Natural Protected Areas of Mexico sponsors a long-term national-wide communal monitoring programme using camera traps, there is not a systematic analysis of the information generated by the programme. We assessed the occurrence of three felids species known to occur in a Zapotec indigenous community conservation area in Oaxaca, Mexico. Specifically, we evaluated how habitat characteristics, human disturbance and prey influence felids’ occurrence across the protected area. None of the variables explained better than the null model the proportion of sites used by Pumas (Puma concolor). Bobcats and Margays favour areas with medium-sized prey. Our study shows the importance of community-based monitoring and information systems (CBMIS) for identifying communal reserve characteristics that contribute to the occupation of carnivores. Further, our results also suggest that management should consider the habitat requirements of felids´ prey. By understanding wildlife habitat use, communal authorities could improve sustainable forest management within the reserves.
biodiversity, CBMIS, conservation, co-occurrence, GLMM, Oaxaca
The combination of natural resources and environmental conditions where organisms can reproduce and survive defines a species´ suitable habitat (
In Mexico, there are six species of felids (Puma, Puma concolor; Jaguarundi, Herpailurus yagouaroundi; Bobcat, Lynx rufus; Ocelot, Leopardus pardalis; Jaguar, Panthera onca; and Margays, Leopardus wiedii). Particularly, Pumas, Bobcats and Margays have been little studied. Studies of habitat use and abundance in pumas stand out (
Since 2008, the Mexican government officially recognised some of these areas through a protection scheme called Voluntary Destined Areas for Conservation (ADVC by its Spanish acronym). The ADVCs must be designated voluntarily by the communities or private owners to get institutional recognition (
Within the ADVCs, one essential activity is wildlife monitoring, which is carried out by locals called “monitores comunitarios” (community monitors). Some authors have questioned this activity for not following a systematised monitoring scheme or lacking analysis of the information generated (e.g.
The National Commission of Natural Protected Areas of Mexico (CONANP by its Spanish acronym) sponsors a national-wide community monitoring programme to increase the knowledge of biodiversity in various country regions by fomenting the use of camera traps to detect wildlife. In Oaxaca, southern Mexico, this programme has been used successfully to assess the status of Jaguar populations (
Here, we analysed the information generated by a long-term wildlife monitoring programme at an ADVC protected by a Zapotec indigenous community to understand habitat use in three felids species (Pumas, Bobcats and Margays). The objective of this study was to evaluate the effect of habitat characteristics, human disturbance and the presence of prey which influence the occurrence of felines in the protected area. We used GLMM models to assess the association between landscape variables related to habitat characteristics, human disturbances and the presence of prey on felids occurrence within the protected area. We hypothesised that felids’ occurrence would be explained by habitat characteristics and prey presence (e.g.
We carried out the study at the Voluntary Conservation Area La Cruz-Corral de Piedra in San Pablo Etla, Oaxaca, Mexico, in a Zapotec ethnic group (17°07' and 17°12'N, 96°39' and 96°48'W; Fig.
The geographical location of the Voluntary Conservation Areas “La Cruz-Corral de Piedra” in the Northern Sierra of Oaxaca, Mexico (red line). Circles and triangles show camera traps installed by community monitors and during our study, respectively. Dark green indicates vegetation dominated by pine, light green by pine-oak and light brown by oak.
The sampling was carried out from January 2013 to June 2017 using 23 camera traps (12 Bushnell Trophy Cam and X8; nine Cuddeback Expert and Capture; two Simmons 119234C). We generated a grid spanning all the ADVC polygons using the ArcGIS programme (v.10.3). We followed previous studies with medium-sized felines to establish the sizes of the cells and generated 26 cells of 1.5 x 1.5 km (
The monitoring programme established by the community allows us to have monthly information regarding the presence/absence of multiple species at each sampling site over time. However, the monitoring programme was not regular for all sampling sites, with maximum sampling periods that ranged from two to 54 months (Suppl. material
We gathered 11 habitat variables associated with each sampling site, following
We assessed shrub layer density, defined as the number of plants per area with > 40 cm and < 2 m in height. Density was estimated using the correction factor proposed by
We also calculated the importance value (IV) of the genus Quercus, the genus Pinus and other tree species (we pooled together plants of the genus: Litesea, Arbutus, Abies, Alnus and Buddleja), as these species were found mostly concentrated in sites, such as glens and creeks, following
Importance Value (IV) = relative density + relative dominance + relative frequency
We assessed seven variables associated with terrain characteristics and human presence: five continuous variables: altitude, slope (assessed with a Clinometer SUUNTO PM-5/360PC); closest distances between the camera trap to any water body (e.g. rivers, dams or waterfalls); closest distances between the camera-trap to any main road within the Reserve (i.e. roads with a width of 7 to 10 m, where cars can circulate); and closest distances to human settlements (i.e. one or more people living in at least one building) and two binary variables: the presence of crags, defined as steep or rugged cliffs; the presence of trails, defined as narrow paths (≤ 2 m wide), formed by the passage of animals or people and located within the vegetation sampling quadrant. For the variables, closest distance to bodies of water, to the main road of the Reserve and human settlements, high-definition topographic maps of San Pablo Etla were reviewed and analysed in the ArcGis programme (version 10.3) generating the Euclidean distance between the sampling stations and the characteristics of the area evaluated.
We identified 16 species as possible prey for the three felids. We classified all potential prey species into three categories, based on size and generated a detection history for each category. In the small-size prey category, we included: Bassariscus astutus, Sylvilagus floridanus, Sylvilagus cunicularis, Conepatus leuconotus, Mephitis macroura, Sciurus aureogaster, Pipilo ocai, Cyanocitta stelleri and Dendrortyx macroura. In the medium-size prey category, we included Eira barbara, Nasua narica, Procyon lotor, Didelphis virginiana and Cuniculus paca. In the large-size prey category, we included Odocoileus virginianus and Dicotyles angulatus.
We used generalised linear mixed models (GLMM) to assess the association between habitat and prey variables with the presence/absence for Pumas, Bobcats and Margays. We defined sampling units as station-month combinations, thus ‘stacking’ the month detection histories, achieving a larger effective sample size. The total sample size for this dataset is n = 531 sites/month combinations for each of the three felid species.
We generated GLMMs with Binomial error and complementary log-log function, as the probability of an event (presence) in our database was small and tended to fit better the data than logistic and probit. We included the presence/absence of Pumas, Bobcats or Margays as a dependent factor, while habitat variables and prey presence/absence were included as fixed factors. We included CameraID and Years as random factors. Including year as a random factor allowed us to meet the assumption of close populations at survey locations since presence/absence in one year is independent of occupancy in the other year. Models were constructed following a forward stepwise procedure, by adding one predictor variable at a time, starting with an empty model (null model) and selecting the variable that provides the best fit to the data according to the Akaike Information Criterion (AIC). We screened sets of predictor variables included in multivariable models to avoid having correlated (r > 0.6) variables together within models. We chose the most parsimonious models using the Akaike Information Criterion (ΔAIC < 2;
We recorded the three focal felids species P. concolor, L. rufus and L. wiedii (Fig.
Images of Puma (Puma concolor: A), Bobcat (Lynx rufus: B) and Margay (Leopardus wiedii: C), obtained by camera trapping in the Voluntary Conservation Area La Cruz-Corral de Piedra in San Pablo Etla, Oaxaca, Mexico.
Images of feline prey: white-tailed deer (Odocoileus virginianus: A), squirrel (Sciurus aureogaster: B), collared peccary (Dicotyles angulatus: C) and white-nosed coati (Nasua narica), obtained by camera trapping in the Voluntary Conservation Area La Cruz -Corral de Piedra in San Pablo Etla, Oaxaca, Mexico.
For the Puma, the presence was best explained by the null model and the model that included medium-size prey (F = 0.69, P = 0.40, Fig.
Best generalised linear mixed models (GLMM) describing the proportion of sampled sites used by the three species felid species. The best five models are presented for each species.
Model | -2log likelihood | AIC | ΔAIC | AICw |
---|---|---|---|---|
Puma | ||||
Null | 1907.21 | 1911.23 | 0.00 | 0.51 |
Medium-size prey | 1907.33 | 1911.35 | 0.12 | 0.48 |
Medium-size prey * Presence of trails | 1921.16 | 1925.18 | 13.95 | 0.00 |
Presence of trails | 1923.71 | 1927.73 | 16.50 | 0.00 |
Medium-size prey + Presence of trails | 1923.89 | 1927.91 | 16.68 | 0.00 |
Bobcat | ||||
Medium-size prey | 2021.75 | 2025.77 | 0.00 | 0.70 |
Medium-size prey + Season | 2024.69 | 2028.72 | 2.95 | 0.16 |
Null | 2025.80 | 2029.82 | 4.05 | 0.09 |
Season | 2027.92 | 2031.94 | 6.17 | 0.03 |
Medium-size prey + Presence of trails | 2028.93 | 2032.95 | 7.18 | 0.02 |
Margay | ||||
Medium-size prey | 2024.41 | 2028.41 | 0.00 | 0.69 |
Medium-size prey + Season | 2026.69 | 2030.71 | 2.30 | 0.22 |
Null | 2094.48 | 2033.48 | 5.07 | 0.05 |
Medium-size prey + Presence of trails | 2031.55 | 2035.57 | 7.16 | 0.02 |
Season | 2031.79 | 2035.79 | 7.38 | 0.02 |
Averaged estimates of the function slopes of variables present in the most parsimonious GLMMs. Estimates of radon factors are shown. Standard errors (SE) and 95% confidence limits (CL) are shown. Tests of significance of variables (F and P) are also given. The asterisk (*) indicates the reference variable.
Variables | Estimate | SE | Lower 95% CL | Upper 95% CL | F (df1, df2) | P |
---|---|---|---|---|---|---|
Puma | ||||||
Null | 0.24 | 0.19 | 0.82 | 1.66 | < 0.001 | |
Medium-size prey | 0.69 (1,529) | 0.40 | ||||
Absence | 0.21 | 0.255 | -0.28 | 0.74 | ||
Presence | 0* | |||||
CameraID | 0.37 | 0.16 | 0.15 | 0.88 | 0.02 | |
Year | 0.05 | 0.05 | 0 | 0.42 | 0.34 | |
Bobcat | 0.43 (1,529) | 0.51 | ||||
Medium-size prey | ||||||
Absence | 0.19 | 0.29 | -0.39 | 0.78 | ||
Presence | 0* | |||||
CameraID | 0.25 | 0.11 | 0.10 | 0.62 | 0.02 | |
Year | 0 | 0.07 | 0 | 0.49 | 0.46 | |
Margay | ||||||
Medium-size prey | 0.20 (1,529) | 0.65 | ||||
Absence | 0.13 | 0.3 | -0.45 | 0.73 | ||
Presence | 0* | |||||
CameraID | 0.27 | 0.12 | 0.11 | 0.66 | 0.02 | |
Year | 0 | 0.11 | 0.21 | 0.60 | 0.69 |
Graphs depicting the relationship between the presence of medium-size prey and Puma (Puma concolor: A), Bobcat (Lynx rufus: B) and Margay (Leopardus wiedii: C). The y axis indicates counts. Model with predictive accuracy by adding one predictor variable at a time, starting with an empty model (null model) and selecting the variable that provides the best fit to the data according to the Akaike Information Criterion (AIC). P-values evaluate the null hypothesis, not model performance. Error bars show standard error. Blue bars refer to absence and red bars refer to presence.
In the study, it was observed that the occurrence of the three felid species was high in areas where their prey species were present, indicating a potential influence of prey on the occupation patterns of these felines. None of the other variables explained the presence of felines; however, there are studies suggesting that human settlements or urbanisation have a negative effect on Pumas, Bobcats and Margays distribution (e.g.
Although the best model suggests that Pumas’ presence was not associated with any of the variables assessed in this study, the second-best model in our study suggests that the occurrence of the Puma was related to the occurrence of medium-size prey, as has been seen in other investigations. (e.g.
Regarding Margay, it is one of the least-studied cats (
Bobcat diet is mainly constituted by lagomorphs and rodents and, to a lesser extent, by opossums, coatis and birds (e.g.
We acknowledge that the study area is very small, considering all three felid species have large home ranges – Margay: 10–21 km2 (
This study highlights that community monitoring (in this case, indigenous monitoring) contributes to scientific knowledge. In this study, information was obtained on the local-scale habitat use of three felids; for one of them, the region represents its southernmost distribution area (L. rufus). Bobcats and Margays favour areas with medium-sized prey. Pumas’ presence did not correlate with the assessed variables. Additionally, the high occurrence of all three felid species in areas with their prey suggests prey influence on their habitat selection. By understanding the relationship between these carnivores, habitat characteristics, human disturbances and the presence of prey, community authorities could improve sustainable forest management. Considering that many indigenous communities around the world protect their natural resources – even without official recognition – (
We thank the communal authorities and residents of San Pablo Etla for access to their lands. To CONANP for providing information on various species of wildlife obtained in San Pablo Etla through the programme “Biological Monitoring of Biodiversity in Protected Natural Areas, Sierra Juárez Mixteca Directorate”. R.E. Galindo contributed with revisions and contributions to the document. MBS thanks the Operation and Promotion of Academic Activities Committee (COFAA) at the IPN. MBS and JRSL thanks the Research Performance Stimulus Program (EDI) at the IPN, as well as the National System of Researchers (SNI) for its recognition and support.
The authors have declared that no competing interests exist.
No ethical statement was reported.
This work was supported by the National Council of Science and Technology (CONACYT) through a scholarship for Master’s studies to NNDB at the Instituto Politécnico Nacional (717791/591006) and the Basic Scientific found a programme (250910 granted to JRSL); and IPN Research and Postgraduate Secretariat (SIP: 20195164 and SIP: 20200713 granted to MBS, and SIP: 20220208 and 20231274 granted to JRSL).
All authors have contributed equally.
J. Roberto Sosa-López https://orcid.org/0000-0002-0120-0704
Eugenio Padilla https://orcid.org/0000-0001-7883-3166
Miguel Briones-Salas https://orcid.org/0000-0003-1413-9946
NNDB, MBS and JRSL conceived and designed the study. NNDB and EP compiled and identified the species for generating the database. NNDB and JRSL performed the statistical analyses. All authors contributed to the interpretation of results and the writing of the manuscript.
Sampling sites and dates on which the camera-traps were installed and Generalized linear mixed models (GLMM) for Puma, Bobcat and Margay
Data type: docx
Explanation note: table S1. Sampling sites and dates on which the camera-traps were installed. table S2. Generalized linear mixed models (GLMM) for Puma. table S3. Generalized linear mixed models (GLMM) for Bobcat. table S4. Generalized linear mixed models (GLMM) for Margay.