Habitat use of jaguar (Panthera onca) in a tropical forest in northern Quintana Roo, Mexico
Dulce María Ávila-Nájera a, b, Marco Antonio Lazcano-Barrero c, Cuauhtémoc Chávez d, Sergio Pérez-Elizalde e, Barbara Tigar f, Germán David Mendoza a, *
a Departamento de Producción Agrícola y Animal, Universidad Autónoma Metropolitana, Unidad Xochimilco, Calzada del Hueso 1100, Col. Villa Quietud, 04960 Ciudad de México, Mexico
b Universidad Politécnica de Huatusco, Unidad Académica de Biotecnología y Agroindustrial, Av. 1 No. 728, Col. Centro, 94100 Huatusco, Veracruz, Mexico
c Reserva Ecológica El Edén A.C., Teocaltiche 207, SM 45, M4, L3, Fracc. Villas del Sol, 77506 Cancún, Quintana Roo, Mexico
d Departamento de Ciencias Ambientales, CBS Universidad Autónoma Metropolitana, Unidad Lerma, Hidalgo Pte. 46, Col. La Estación Lerma, 52006 Lerma, Estado de México, Mexico
e Colegio de Postgraduados, Campus Montecillo, Carretera México-Texcoco Km. 36.5, 56230 Montecillo, Texcoco, Mexico
f School of Forensic and Applied Sciences, University of Central Lancashire, Preston PR1 2HE, United Kingdom
*Corresponding author: gmendoza@correo.xoc.uam.mx (G.D. Mendoza)
Abstract
The study was carried out in a median stature tropical forest in the north of Quintana Roo during 2008, 2010-2012, using camera traps to record activity and habitat use by jaguar and other vertebrates. Spatial and temporal variations in habitat use, activity patterns, and species associations were analyzed using a chi-squared test (x2). Generalized linear models (GLM) were used to determine which variables influenced the presence of jaguars in the study area. Jaguars were most frequently recorded in acahual (80%) and on roads (66.9%). There were significant differences between years, with selection for particular habitat and trail types in 2008 and 2011 (all χ2 p <0.001). Jaguars were mainly crepuscular-nocturnal (68.5% records). They showed spatial associations with Meleagris ocellata, Dasypus novemncinctus, Mazama temama, Nasua narica, Odocoileus virginianus, Leopardus pardalis, Didelphis spp., and Pecari tajacu. Generalized linear models based on prey, co-predators, and human activity gave the best fit with jaguar occurrence. The felids responded to environmental changes (hurricanes and fires) by modifying their spatial and temporal use of resources.
Keywords:
Acahuales; Associated species; Natural disturbances; Activity patterns; Environmental variables
© 2019 Universidad Nacional Autónoma de México, Instituto de Biología. Este es un artículo Open Access bajo la licencia CC BY-NC-ND
Uso de habitat del jaguar (Panthera onca) en un bosque tropical del norte de Quintana Roo, México
Resumen
El estudio se realizó en una selva tropical mediana del norte de Quintana Roo, México, durante 2008, 2010- 2012. Se empleó la técnica de trampeo fotográfico. Las variables ambientales, los patrones de actividad y las asociaciones de especies se analizaron espacial y temporalmente mediante una prueba de ji cuadrada (x2). Los modelos lineares generalizados (MLG) se usaron para determinar qué variables influyen en la presencia del jaguar en el área de estudio. El jaguar se registró con mayor frecuencia en acahuales (80%) y en caminos anchos (66.9%). Hubo diferencias significativas entre los años, con la selección de hábitat y tipos de senderos en 2008 y 2011 (x2 p <0.001). El jaguar fue principalmente crepuscular-nocturno (68.5%). Se muestran asociaciones espaciales con Dasypus novemncinctus, Didelphis spp., Leopardus pardalis, Mazama temama, Meleagris ocellata, Nasua narica, Odocoileus virginianus y Pecari tajacu. El MLG que incluye presas potenciales, codepredadores y humanos, es el que mejor explica la presencia del jaguar en el área de estudio. Se observó una respuesta a los cambios ambientales (huracanes e incendios), estos modifican el uso espacial y temporal de los recursos.
© 2019 Universidad Nacional Autónoma de México, Instituto de Biología. This is an open access article under the CC BY-NC-ND license
Palabras clave:
Acahuales; Especies asociadas; Perturbaciones naturales; Patrones de actividad; Variables ambientales
Introduction
A habitat is the physical space where a species naturally occurs, and includes all the physical and biological characteristics (environmental factors) necessary to cover its basic needs. However, spatial and temporal variation in habitat conditions generate a selective pressure resulting in the use of different habitat components (Cody, 1985; Krauman, 1999) and this will vary according to the availability, richness, or utility of any resources that offer an advantage towards satisfying the requirements of a species (Hutton, 1985; Johnson, 1980).
Quantifying the abiotic and biotic resources that determine the presence of species in a habitat is complex, because the selection of resources by one species also depends upon the resource-use by the predators, co-predators, and prey species that occur there (McLoughlin et al., 2010; Slagsvold, 1980). In addition, the impact and changes that result from human activity on an ecosystem play an important role in habitat selection and resource-use by many animals (Balestrieri et al., 2009; McLoughlin et al., 2010).
The study and management of large carnivores has proven to be one of the best strategies for conserving their wider ecosystem because of their large home-ranges and extensive areas of activity. Therefore, in order to protect these predators, expansive areas of habitat are protected along with all those species that inhabit the same area. One such carnivore is the jaguar (Panthera onca), which is considered both a keystone and umbrella species, as well as an indicator of the proper functioning of its environment (Manterola et al., 2011; Miller & Rabinowitz, 2002; Núñez et al., 2002). Therefore, understanding the jaguar’s use and selection of resources is of high interest in terms of conserving this species and the ecosystems where it occurs. Despite the importance of the jaguar, its populations have been declining rapidly due to habitat loss in parts of Mexico and Central America (Caso et al., 2008; De la Torre et al., 2017) and it is considered to be a species in danger of extinction (Armella & Yáñez, 2011).
Current knowledge of resource selection by the jaguar suggests that it prefers forested areas near water and avoids disturbed or modified areas (Chávez, 2010; Conde et al., 2010; Cullen, 2006). The habitat components that play an important role in jaguar distribution include vegetation type and habitat continuity (De la Torre et al., 2017; Hidalgo-Mihart et al., 2017), land cover, water availability, prey (Cullen et al., 2013; Scognamillo et al., 2003), and anthropogenic activities (Núñez et al., 2002; Romeu, 1996; Sanderson et al., 2002; Zarza et al., 2007).
Current information about jaguar resource-use only covers a limited part of its geographical range. Studies from Mexico include one from the southern Yucatán Peninsula (Chávez, 2010), and another from the northern Yucatán Peninsula on their population density in the Ría Lagartos Reserve and surrounding areas including the Zapotal, Yum Balam, and Dzilam Reserves (Faller et al., 2007). The Yucatán Peninsula is considered a priority study area for jaguar (Ceballos et al., 2005; Sanderson et al., 2002) because it contains original native vegetation types that give an extensive continuity of habitat and maintain communication with the Maya Forest (De la Torre et al., 2017; Hidalgo-Mihart et al., 2017). Therefore, the aim of this study was to analyze resource-use by jaguars, as well as activity patterns, and information about associated species and environmental variables that influence their occurrence in the tropical forest in the north-eastern Yucatán Peninsula.
Materials and methods
The study was conducted in El Edén Ecological Reserve (EER) and surrounding areas in the municipality of Lázaro Cárdenas, Quintana Roo, in the north-east Yucatán Peninsula, Mexico (21°36’00”-20°34’00” N, 87°06’00”-87°45’00” W; Fig. 1a-e). The EER covers an area of 3,077 ha in the Yalahau region in the northernmost tropical forests of North America, at an elevation of 5-10 m asl (Gómez-Pompa et al., 2003, 2010; Lazcano-Barrero et al., 1992). This encompasses the Yum Balam protected area, which has extensive areas of medium-stature deciduous forest (MSF) and secondary forest or acahual (Navarro et al., 2007). MSF is characterized by trees reaching heights of 10-15 m. The most common species are Manilkara zapota, Lysiloma latisiliqua, Metopium brownei, and Thrinax radiata palm. Acahual is derived from MSF and characterized by trees of 7-10 m in height, typically M. brownei, Bursera simaruba, and Swartzia cubensis. There is also some flooded forest or tintal, dominated by Haematoxylon campechianum, Erythroxylon confusum, and Byrsonima bucidaefolia, plus small fragments of Annonaceae forest and savanna (Schultz, 2003) (Fig. 1).
The study was conducted during July-September 2008, October-December 2010, May-July 2011, and August-December 2012 using camera traps (Cuddeback Expert, Capture, Capture IR, Moultrie, and Wildview). Camera placement followed the CENJAGUAR design (Chávez et al., 2007) with 2 or 3 sampling stations in 9 study plots each measuring 9 km2. The location (longitude and latitude) of each camera station was taken with a Global Positioning System (GPS-Garmin eTrex) and entered into ArcView 8 Geographical Information System (GIS).
The cameras were sited on forest paths or trails, firebreaks, and roads (about 6 m wide) at a distance of 1-3 km from one another. In 2008, there were 9 camera stations in MSF and 18 in acahual. In 2010, there were 10 in MSF and 14 in acahual, in 2011, 10 in MSF and 12 in acahual, and in 2012, 12 in MSF and 22 in acahual (Fig. 1a-d). Despite the variation in the number of cameras over the years, they occupied the same area. The cameras operated 24 h d and were checked and images downloaded once every 15 d.
Photographs were grouped by trap site to allow comparison of trends in the spatial distribution of individual felids, which were identified by their unique coat patterns. Photographic records were considered to be independent events when 1) the same individual was recorded at intervals of more than 30 min apart, 2) it was possible to distinguish between different individuals of the same species in consecutive photos, 3) a photograph showed several different individuals that could be clearly identified, and 4) where members of the same species could not be identified as being different individuals in consecutive photos, a new event was only recorded once in every 3 h (Ávila-Nájera et al., 2016).
Terrestrial or semi-terrestrial mammals and some birds of terrestrial habits that appeared in the photographs were considered to be associated species, and their activity was recorded. For some analyses they were grouped by taxon and weight into 5 classes: large mammals (8-26 kg), medium-sized mammals (3-8 kg), small mammals (1-3 kg), large birds (2-6 kg), and small birds (< 2 kg) (Davis et al., 2010). In addition, photographic records of human activity were included.
Photographs were grouped according to the time of the record: nocturnal (20:00-06:00), crepuscular (06:00-08:00, 18:00-20:00), and diurnal 08:00-18:00 (Monroy-Vilchis et al., 2011), and results were pooled into periods of 2 h (starting at 00:00) in order to detect peak times of activity. All species recorded by the camera traps were classified according to their activity pattern using the following rules: diurnal when ˂ 15% of observations were at night, mainly diurnal when 15-35% of observations at night, nocturnal when > 85% of observations were at night, and mainly nocturnal when 65-85% of observations were at night.
Traps were classified according to vegetation (MSF and acahual) and trail type (whether the camera was by a forest trail, firebreak, or dirt road), and the frequency of capture in the photos was used to determine jaguar activity patterns in each habitat.
Pearson’s chi-squared test (Chernoff & Lehmann, 1954; Siegel & Castellan, 1988) was used to test for differences in jaguar use of vegetation, trail type, and activity pattern between the study years. The test compared the observed activity (from photographs) and expected activity (percentage of cameras sited by each trail type), and the activity time of jaguars (from photographs) and an expected activity pattern of 10/24 h for diurnal/nocturnal and 4/24 h for crepuscular activity. A Monte Carlo simulation was used to estimate the significance levels (p-value; Hope, 1968).
Overall differences in jaguars’ use of resources and activity patterns were tested for evidence of selection of type of vegetation (MSF or acahual), trail type (forest trail, firebreak or dirt road), and time of activity using the equation: H0: pi=p0, where the real proportion is pi for the i-th category, p0 is the expected proportion according to the null hypothesis, at the observed proportion for the i-th category ̂p i. H0 was rejected when the bootstrap interval for pi did not contain p0.
We used 2 × 2 contingency tables to test the null hypothesis that different species’ records in the traps were independent of each other (at α = 0.05; Ludwig & Reynolds, 1988) and using Yates’ correction factor:
Xc2 = (N[|(ad-(bc))|-(N/2)])2/mnrs.
where a = number of sampling units where both species occur, b = number of sampling units where a occurs, but not b, c = number of sampling units where b occurs, but not a, d = number of sampling units where none of the 2 species were found, N= total number of sampling units (N = a + b + c + d); m = a + b; n= c + d; r = a + c; s = b + d. Therefore, if Xc2 > Xc2, where Xt2 is 1 – α for a Χ2 distribution with 2 degrees of freedom, the species are not independent, and the test is rejected.
To test for the overall effect of environmental and resources on jaguar occurrence we used Generalized Linear Models (GLM) (Sileshi, 2008), calculated using the glm.nb function in R. We assumed a negative binomial distribution and considered camera stations 1.5 km or farther apart to be independent sampling units. We also assumed that all vegetation and path types were equally accessible to jaguar and their prey and co-predators recorded in the traps (Table 1).
The following variables were used in the GLMs: vegetation type, horizontal and vertical plant cover, path type, potential prey species, co-predators (puma and ocelot), human activity, activity patterns of potential prey and co-predators, distance to the nearest water body and town, latitude, and longitude. Horizontal plant cover is the percentage of vegetation capable of hiding a subject 15 m away, and vertical plant cover is the visibility through the tree canopy (Griffith & Youtie, 1988). Arc Map 10.3 (GIS) was used to determine the distance to the nearest water body and town. Significance was determined using an asymptotic t-test at α = 0.05 level, with Akaike’s Information Criterion (AIC) to select the best model (Burnham & Anderson, 2002), since this favors simpler models and acts to reduce the number of parameters.
Table 1
List of variables selected according to the Akaike information criterion (AIC) (Burnham & Anderson, 2002) and incorporated into 18 models for each year and 18 models using all the data from the 4 years (90 models in total).
Model name |
Variables selected using AIC |
m0 |
Null |
m2 |
Co-predator activity pattern |
m4 |
Associated species records, adjusted for the presence of humans |
m6 |
Associated species, co-predator records and humans |
m8 |
Associated species records with activity pattern, adjusted for the presence of humans |
m10 |
Activity pattern of associated species and co-predator records, adjusted for the presence of humans |
m12 |
Environmental factors (distance to body water, vegetation types, trail type, longitude, latitude) |
m14 |
Environmental factors and associated species records, adjusted for the presence of humans |
m16 |
Associated species records by taxa and weight group |
m18 |
Associated species records by taxa and weight group and environmental factors |
Results
A total of 181 independent photographic records of jaguar were captured, and a total 12 individual jaguars were identified during the study. There were 96 photographs taken in 2008, 15 in 2010, 12 in 2011, and 58 in 2012, with 7, 5, 5, and 12 individual jaguars identified in each of the study years, respectively.
In some years, there were differences in the observed and predicted frequencies of jaguars in each vegetation and trail type, and by activity pattern (Table 2). Evidence of selective resource-use include significant results for vegetation and trail type in 2008 and 2011 (all at p < 0.001) and activity pattern in 2011 (p < 0.001) and 2012 (p = 0.01) (Table 2).
Table 2
Results of x2 Pearson tests with Monte Carlo estimate for selection vegetation type, trail type, and activity pattern (nocturnal, diurnal, crepuscular) by year on the observed frequency of jaguar records from camera traps in El Edén Ecological Reserve in Quintana Roo, Mexico.
Year of study |
||||||||
Variable |
2008 |
2010 |
2011 |
2012 |
||||
x2 |
p |
x2 |
p |
x2 |
p |
x2 |
p |
|
Vegetation type |
21.9 |
0.00 |
0.31 |
0.86 |
0.159 |
0.00 |
1.76 |
0.41 |
Trail type |
36.9 |
0.00 |
0.61 |
0.89 |
229 |
0.00 |
1.78 |
0.62 |
Activity pattern |
1.40 |
0.71 |
5.96 |
0.11 |
177 |
0.00 |
11.3 |
0.01 |
Most jaguar records were from acahual (80%) with only 20% from MSF. However, the frequency of jaguar records by vegetation type varied depending upon the year of study. There were more records from acahual in 2008, 2010, and 2012 (95%, 62% and 75%, total jaguar records by year of study, respectively) and more from MSF in 2011 (71%). In addition, the observed frequency of records in acahual was higher than expected in 2008 and 2011, but lower than expected in MSF (Table 3).
Jaguars were active both day and night within the study area, although there were differences in activity patterns between the years, with more diurnal records and fewer crepuscular records in 2011. Overall, almost half (43%) of the records were nocturnal (20:00-06:00) with the main peaks of activity occurring between 21:00-22:00 and 02:00-03:00. About a quarter (25.5%) of the records were crepuscular, with most activity occurring between 19:00-20:00 and almost a third (31.3%) of the records occurring during the day (08:00-18:00) with fewer observations between 10:00-11:00 and 13:00-14:00. In 2011, the observed frequently of diurnal records was higher than expected, whilst for crepuscular records it was lower than expected (Table 3). There was no difference in the expected and observed frequency of nocturnal records over the 4 years of study (Table 3).
Table 3
Test of whether jaguar resource-use varied in El Edén Ecological Reserve, Quintana Roo, Mexico according to availability testing for the effect of vegetation type, trail type, and activity pattern with expected values estimated from 10,000 bootstrap simulations (where the test outcome shows Po < Pe, Po > Pe or no difference between Po and Pe (=), with corresponding Bonferroni confidence intervals).
Resource |
Classification type |
Year |
Observed proportion (Po) |
Expected proportion (Pe) |
Confidence interval range |
Test outcome (Po < > or = Pe) |
Vegetation |
Medium-stature forest |
2008 |
0.04 |
0.33 |
0.00-0.12 |
< |
2010 |
0.38 |
0.5 |
0.07-0.69 |
= |
||
2011 |
0.21 |
0.45 |
0.18-0.25 |
< |
||
2012 |
0.25 |
0.34 |
0.11-0.40 |
= |
||
Acahual |
2008 |
0.95 |
0.66 |
0.83-1 |
> |
|
2010 |
0.61 |
0.5 |
0.30-0.92 |
= |
||
2011 |
0.78 |
0.55 |
0.74-0.81 |
> |
||
2012 |
0.75 |
0.65 |
0.59-0.88 |
= |
||
Camera site |
Forest trail |
2008 |
0.04 |
0.37 |
0.00-0.12 |
< |
2010 |
0.38 |
0.5 |
0.07-0.76 |
= |
||
2011 |
0.22 |
0.5 |
0.18-0.26 |
< |
||
2012 |
0.25 |
0.34 |
0.11-0.40 |
= |
||
Firebreak |
2008 |
0.06 |
0.12 |
0.00-0.15 |
= |
|
2010 |
0.15 |
0.18 |
0.00-0.46 |
= |
||
2011 |
0.13 |
0.1 |
0.00-0.16 |
= |
||
2012 |
0.22 |
0.20 |
0.08-0.37 |
= |
||
Road |
2008 |
0.88 |
0.5 |
76.19-0.98 |
> |
|
2010 |
0.46 |
0.31 |
0.15-0.84 |
= |
||
2011 |
0.64 |
0.4 |
0.59-0.69 |
> |
||
2012 |
0.52 |
0.16 |
0.09-0.36 |
> |
||
Activity pattern |
Crepuscular |
2008 |
0.22 |
0.44 |
0.35-0.69 |
= |
2010 |
0.23 |
0.41 |
0.00-0.53 |
= |
||
2011 |
0.18 |
0.41 |
0.15-0.22 |
< |
||
2012 |
0.30 |
0.16 |
0.15-0.47 |
= |
||
Diurnal |
2008 |
0.42 |
0.41 |
0.26-0.58 |
= |
|
2010 |
0.46 |
0.16 |
0.07-0.84 |
= |
||
2011 |
0.63 |
0.41 |
0.58-0.67 |
> |
||
2012 |
0.22 |
0.41 |
0.08-0.37 |
< |
||
Nocturnal |
2008 |
0.34 |
0.41 |
0.19-0.50 |
= |
|
2010 |
0.30 |
0.41 |
0.00-0.69 |
= |
||
2011 |
0.18 |
0.16 |
0.14-0.22 |
= |
||
2012 |
0.47 |
0.41 |
0.30-0.64 |
= |
Figure 2 (a-d) shows the annual relative abundance of other species recorded in the camera traps, grouped according to their taxa and size class, and whether they are co-predators (Fig. 2a) or potential prey (Fig. 2b-d).
There was a significant association (x2 = 3.84, p < 0.05) between jaguar and some other species recorded in the camera traps including potential prey, Aramydes cajanea, Didelphis sp., Meleagris ocellata, Mazama temama, Nasua narica, Odocoileus virginianus, and Pecari tajacu, and the co-predator Leopardus pardalis. When camera trap results were grouped by vegetation type, there were some interesting associations between other species and jaguar. These varied depending upon the year of study. There was an association between jaguar records and traps sited in acahual in 2008 (x2 = 3.84, p < 0.05), as well as for jaguar records and Dasypus novemcinctus and M. temama in 2010, L. pardalis in 2011 and for M. temama, O. virginianus, and Ortalis vetula in 2012. In MSF in 2008, there was a positive association between records of jaguar and Didelphis sp., and for N. narica in 2010, as well as for Leopardus wiedii, N. narica, P. tajacu, and Sciurus sp. in 2012, but there were no species associated with jaguar records from MSF in 2011.
Ninety GLMs were developed using the variables in Table 1 (18 per year and 18 for all 4 years combined). The effect of year was significant in all models, and GLMs based on environment (acahual, roads, and firebreaks) were the best at predicting jaguar presence. However, model 10 (292 AIC) (Table 4) had the lowest AIC value, and used the activity patterns of co-predators, associated species, and humans to predict the presence of jaguar.
Table 4. Model for prey species and co-predators by activity pattern and adjusted for human presence, after variable selection by AIC.
Model 10 |
Estimate |
SE |
z value |
Pr (> |z|) |
(Intercept) |
564 |
173 |
3.26 |
0.00 |
Variable: year |
-0.28 |
0.09 |
-3.29 |
0.00 |
Puma concolor diurnal activity |
0.42 |
0.12 |
3.60 |
0.00 |
Leopardus pardalis crepuscular activity |
0.95 |
0.30 |
3.20 |
0.00 |
Mazama temama crepuscular activity |
-1.77 |
0.88 |
-2.02 |
0.04 |
Odocoileus virginianus diurnal activity |
0.98 |
0.25 |
3.99 |
0.00 |
Urocyon cinereoargenteus diurnal activity |
0.32 |
0.12 |
2.60 |
0.01 |
Aramides cajanea nocturnal activity |
0.37 |
0.17 |
2.10 |
0.04 |
Crax rubra diurnal activity |
-0.79 |
0.28 |
-2.82 |
0.00 |
Dasypus novemncinctus nocturnal activity |
-3.21 |
0.97 |
-3.31 |
0.00 |
Discussion
Although felids have been considered generalist species by many authors, some studies suggest they have specific habitat requirements and need certain components to be present in order to survive in an area (Chávez et al., 2010). In this study, there was evidence that jaguars selected and used resources according to their availability. However, they made higher than expected use of acahual and roads, although this varied by year, and their general presence within EER was defined by the abundance and activity patterns of associated species and co-predators. Jaguars changed their use of vegetation type depending upon the year of study, and this was particularly influenced by other factors such as the presence of humans or prey species. Secondary forest was a preferred habitat, and we believe that this is because it provides both horizontal and vertical vegetation cover, which are essential for jaguars’ permanence in a habitat (Conde et al., 2010).
Existing jaguar distribution data suggest that they inhabit many habitats (Zarza et al., 2007), favoring dense vegetation and areas that provide sufficient water and a diverse and abundant variety of prey (Manterola et al., 2011; Sanderson et al., 2002). They occur primarily in tropical rain forest, followed by tropical dry forest, and are at their lowest densities in arid environments and grasslands (Sanderson et al., 2002). Previous studies from the Yucatán Peninsula show that jaguars have a marked preference for areas with dense vegetation cover (Ceballos et al., 2005; Chávez et al., 2007). In EER, jaguars were recorded on roads, firebreaks, and forest trails, as reported by other authors (Maffei et al., 2004). Most records were from roads, and many authors suggest that jaguars use them to move along or for performing other activities such as scent marking. However, Conde et al. (2010) report a gender bias, with more records of male jaguars from roads, and this may in turn influence the gender of other road users (Maffei et al., 2011). In our study area there was very little traffic, making roads a relatively safe and undisturbed environment, often surrounded by dense vegetation cover, so that their open terrain would facilitate animal movement (Harmsen et al., 2010).
Most jaguar records in the EER were crepuscular or nocturnal, with 2 peaks in nocturnal activity, and crepuscular activity generally occurring around dusk. Similarly, Chávez et al. (2007) found that jaguars in the south of the Yucatán Peninsula (Calakmul) were mainly crepuscular or nocturnal, although peak activity occurred in the early morning. In other parts of their range, jaguars can be active 24 h d, although this varies with locality, and depends upon the availability and activity patterns of their prey (Carrillo, 2000; Scognamillo et al., 2003), and in some case they are mainly diurnal (Álvarez-Castañeda & Patton, 2000; Foster et al., 2013; Harmsen et al., 2009; Maffei et al., 2004, 2011; Rabinowitz & Nottingham, 1986). In the Chaco region of Bolivia, jaguar had 2 peak periods of activity, 03:00-06:00 and 17:00-22:00, but were less active around midnight (Maffei et al., 2004) reflecting the activity of prey such as Mazama sp. and Tayassu sp. (Barrientos & Maffei, 2000). The overall resource-use and nocturnal/crepuscular activity pattern seen in EER generally agrees with previous studies, although the diurnal activity was higher than expected.
The spatial association between jaguar and prey species in EER can be explained by the relative abundance of prey, for example D. novemcinctus, Didelphis sp., M. temama, and N. narica in 2008, and for O. vetula and Sciurus sp. in 2012. However, in 2010 there was a positive association between jaguar and E. barbara, L. wiedii, M. temama, N. narica, and T. tajacu recorded in the camera traps, despite these species not being at high relative abundance. The contingency tables and GLMs suggested that A. cajanea, D. novemcinctus, Didelphis sp., L. pardalis, Meleagris ocellata, M. temama, N. narica, O. virginianus, and P. tajacu were associated with jaguar records in the traps, and these are reported as important prey by Ceballos et al. (2005), Chávez (2010), Chávez et al. (2007), and Davis et al. (2010).
The environmental variables that influence the presence of the jaguar vary along its range alongside its resource-use (Davis et al., 2010). In EER, the year of study had a significant effect on jaguar occurrence in models for all 4 years and the combined models. The photographic catch rates of all species decreased significantly in 2011, following a fire around the reserve, probably inhibiting the activity of all animals nearby. There was also a flood in 2010, which is likely to have had considerable impact on animal mobility, and therefore to have reduced their relative abundance as estimated by the camera traps. The main environmental factor in the models that determined jaguar occurrence was a positive association with acahual, probably because it encompasses many roads, which allow jaguars to move easily in search of prey, as seen in a study from Belize (Davis et al., 2010). Water availability had no impact on jaguar occurrence because EER has large and permanent water sources throughout the year, so it is never a limiting factor.
There was a positive effect of small and large birds and small mammals on jaguar occurrence, and although large birds are a commonly recorded prey item in felid diets, other small vertebrates are not. Models based on the presence of co-predators found that the crepuscular ocelot and nocturnal puma activity had a significant effect on jaguar activity. This agrees with previous studies from this area, and is thought to represent a temporary lack of segregation of predators (Chávez, 2010), also seen in Belize (Davis et al., 2010). However, according to the Akaike value, the principal variables that determined the presence of the jaguar in the study were the activity patterns of their main prey species, such as Odocoileus virginianus (white-tailed deer), Mazama temama (brocket deer), Dasypus novemncinctus (armadillo), and Crax rubra (curassow).
Many researchers have noted the effect of humans on the presence of jaguar (Ávila-Nájera et al., 2011; Ceballos et al., 2005; Zarza et al., 2007), particularly in relation to the decline of prey species, with fewer jaguar occurring in areas of intensive hunting (Chávez, 2010; Montiel et al., 1999). In EER, humans were the most frequently recorded species in the camera traps, and were one of the most important factors that influenced jaguar occurrence in the GLM models, agreeing with similar findings by Davis et al. (2010). Although human activity is an important factor that negatively affects the occurrence and activity of jaguars, it did not seem to influence their presence in this area. This may be because most photographic records were of workers, who pose a minimum threat to jaguar and their prey.
Even though EER is not a large area, it has suitable environmental components including adequate vegetation cover, water availability, and potential prey that allow for the continued presence of the jaguar. There is sufficient primary and secondary MSF within the EER and the surrounding area, as well as sufficient key prey, including ungulates (brocket deer and peccary), opossums, raccoons, and oscillated turkeys. However, these need to be maintained at adequate populations to support sustainable populations of jaguar. Jaguars frequently used roads within the reserve, which received low-level traffic, and were not affected by the presence of humans within the area. However, the long-term survival of jaguars in small areas such as El Edén relies upon the connectivity with other suitable protected areas through the establishment and maintenance of biological corridors. These results provide further evidence to support the planning and management of jaguars in the northern Yucatán Peninsula and adjacent areas.
Acknowledgements
To El Edén Ecological Reserve, the Center of Tropical Research of Veracruz University, and Kathy Cabrero for their support. The study was made possible by funding from Conacyt grant number 211454 and Promep 54310009. To C. Chávez (UAM-PTC-333); Erik J. Torres, Juan Castillo, Alejandro Pacheco, Brady Hollinsgworth, EER staff, and volunteers at Global Vision International for their support in logistics and fieldwork.
References
Álvarez-Castañeda, S. T., & Patton, J. L. (2000). Mamíferos del noroeste de México II. La Paz, Baja California Sur: Centro de Investigaciones Biológicas del Noroeste, S.C.
Armella, A., & Yáñez, M. L. (2011). Mamíferos mexicanos en peligro de extinción. Revista Digital Universitaria, 12, 1−10.
Ávila-Nájera D. M, Chávez, C., Lazcano-Barreto, M. A., Mendoza, G. D., & Pérez-Elizalde, S. (2016). Traslape en patrones de actividad y traslape entre grandes felinos y sus principales presas en el norte de Quintana Roo, México. Therya, 7, 439−448.
Ávila-Nájera, D. M., Rosas-Rosas, O., Tarango-Arámbula, L. A., Martínez-Montoya, J. F., & Santoyo-Brito, E. (2011). Conocimiento, uso y valor cultural de seis presas del jaguar (Panthera onca) y su relación con éste, en San Nicolás de los Montes, San Luis Potosí, México. Revista Mexicana de Biodiversidad, 82, 1020−1028.
Balestrieri, A., Remonti, L., & Prigioni, C. (2009). Habitat selection in a low-density badger Meles meles population: a comparison of radio-tracking and latrine surveys. Wildlife Biology, 15, 442−448.
Barrientos, S., & Maffei, L. (2000). Radio-telemetría de la hurina Mazama gouazoubira en el campamento Cerro Cortado, Izozog, Santa Cruz, Bolivia. In E. Cabrera, C. Mercolli, & R. Resquin (Eds.), Manejo de fauna silvestre en Amazonia y Latinoamérica (pp. 369−372). Florida, Asunción: CITES Paraguay/ Fundación Moises Bertoni/ University of Florida.
Burnham, K. P., & Anderson, D. R. (2002). Model selection and multimodel inference. New York: Springer.
Caso, A., López-González, C., Payan, E., Eizirik, E., Oliveira, T., Leite-Pitman, R. et al. (2008). Panthera onca. The IUCN Red List of Threatened Species 2008: e.T15953A5327466. https://www.iucnredlist.org/
Carrillo, E. (2000). Ecology and conservation of white-lipped peccaries and jaguars in Corcovado National Park, Costa Rica (Ph.D. Thesis). Massachusetts University, Massachusetts.
Ceballos, G., Chávez, C., Zarza, H., & Manterola, C. (2005). Ecología y conservación del jaguar en la región de Calakmul. Conabio. Biodiversitas, 62, 1−7.
Chávez, C. (2010). Ecología y conservación del jaguar (Panthera onca) y puma (Puma concolor) en la región de Calakmul y sus implicaciones para la conservación de la península de Yucatán (Ph.D. Thesis). University of Granada, España.
Chávez, C., Ceballos, G., & Amín, M. (2007). Ecología poblacional del jaguar y sus implicaciones para su conservación en la península de Yucatán. In G. Ceballos, C. Chávez, R. List, & H. Zarza (Eds.), Conservación y manejo del jaguar en México, estudios de caso y perspectivas (pp. 91−100). México D.F.: Conabio/ UNAM/ Alianza WWF-Telcel.
Chernoff, H., & Lehmann, E. L. (1954). The use of maximum likelihood estimates x2 tests for goodness-of-fit. The Annals of Mathematical Statistics, 25, 579−586.
Cody, M. L. (1985). Habitat selection in birds. New York: Academic Press.
Conde, D. A., Colchero, F., Zarza, H., Christensen, J. R., Sexton, J. O., Manterola, C. et al. (2010). Sex matters: modeling male and female habitat differences for jaguar conservation. Biological Conservation, 143, 1980−1988.
Cullen, L. (2006). Jaguars as landscape detectives for the conservation of Atlantic Forest in Brazil (Ph.D. Thesis). Kent University, Canterbury.
Cullen, L., Sana, D. A., Lima, F., Abreu, K., & Uezu, A. (2013). Selection of habitat by the jaguar, Panthera onca (Carnivora: Felidae), in the upper Paraná River, Brazil. Zoologia, 30, 379−387.
Davis, L., Kelly, M., & Stauffer, D. F. (2010). Carnivore co-existence and habitat use in the Mountain Pine Ridge Forest Reserve, Belize. Animal Conservation, 4, 56−65.
De la Torre, A., Núñez, J. M., & Medellín, R. (2017). Spatial requirements of jaguars and pumas in Southern Mexico. Mammal Biology, 82, 52−60.
Faller, J. C., Chávez, C., Johnson, S., & Ceballos, G. (2007). Densidad de la población de jaguares en el noreste de la península de Yucatán. In G. Ceballos, C. Chávez, R. List, & H. Zarza (Eds.), Conservación y manejo del jaguar en México estudios de caso y perspectivas (pp. 111−121). México D.F.: Alianza WWF-Ttelcel/ Conabio/ Conanp/ EcoCiencias, S.C.
Foster, V., Sarmento, P., Sollman, R., Torres, N., Jácomo, A., Negroes, N. et al. (2013). Jaguar and puma activity patterns and predator-prey interactions in four Brazilian biomes. Biotropica, 45, 373−379.
Gómez-Pompa, A., Allen, M., Fedick, S., & Jiménez-Osorio, J. J. (2003). The lowland Maya area. Three millennia at the human-wildland interface. Binghamton, New York: The Haworth Press.
Gómez-Pompa, A., Lazcano-Barrero, M. A., Gómez-Barrero, A., & Macswiney, C. (2010). La Reserva Ecológica El Edén: proyecto de conservación privada. In J. Carabias, J. Sarukhán, J. De la Maza, & C. Galindo (Eds.), Patrimonio natural de México. Cien casos de éxito (pp. 92−93). México D.F.: Comisión Nacional para el Conocimiento y Uso de la Biodiversidad.
Griffith, B., & Youtie, B. A. (1988). Two devices for estimating foliage density and deer hiding cover. Wildlife Society Bulletin, 16, 206−210.
Harmsen, B., Foster, R., Silver, S., Ostro, L., & Doncaster, P. (2009). Spatial and temporal interactions of sympatric jaguars (Panthera onca) and pumas (Puma concolor) in a neotropical forest. Journal of Mammalogy, 90, 612−620.
Harmsen, B., Foster, R., Silver, S., Ostro, L., & Doncaster, P. (2010). Differential use of trails by forest mammals and the implications for camera-trap studies: a case study from Belize. Biotropica, 42, 126−133.
Hidalgo-Mihart, M., Contreras-Moreno, F., De la Cruz, A. J., & Juárez-López, R. (2017). Validation of Calakmul-Laguna de Términos corridor for jaguars Panthera onca in south-eastern Mexico. Oryx. Flora and Fauna International. doi:10.1017/S0030605316001083.
Hope, A. C. (1968). A simplified Monte Carlo significance test procedure. Journal of the Royal Statistical Society, 30, 582−598.
Hutton, L. R. (1985). Habitat selection by nonbreeding migratory land birds. In M. L. Cody (Ed.), Habitat selection in birds (pp. 455−476). Orlando: Academic Press.
Johnson, H. D. (1980). The comparison of usage and availability measurements for evaluating resource preference. Ecology, 61, 65−71.
Krauman, R. P. (1999). Some basic principles of habitat use, grazing behavior of livestock and wildlife. Wildlife and Range Experiment Station Bulletin, 70, 85−90.
Lazcano-Barrero, M. A., March, J., Núñez, H., Ruelas, E., Muñoz-Alonso, A., & Martínez, R. (1992). Inventario faunístico de la Reserva El Edén, Quintana Roo: una prospección. San Cristóbal de Las Casas, Chiapas: ECOSFERA, A. C.
López-González, C., & Miller, B. (2012). Do jaguars (Panthera onca) depend on large prey? Western North American Naturalist, 62, 218−222.
Ludwing, J., & Reynolds, J. (1988). Interspecific association. In J. Ludwing, & J. Reynolds (Eds.), Statistical ecology. A primer on methods and computing (pp. 107−144). New York: Wiley Interscience.
Maffei, L., Cuellar, E., & Noss, A. (2004). One thousand jaguars (Panthera onca) in Bolivia´s chaco? Camera trapping in the Kaa-Iya National Park. Journal of Zoology, London, 262, 295−304.
Maffei, L., Noss, A., Silver, S., & Kelly, J. (2011). Abundance density case study: jaguars in the Americas. In A. O´Connell, J. Nichols, & U. Karanth (Eds.), Camera traps in animal ecology. Methods and analyses (pp. 119−144). Tokio: Springer.
Mandel, M., & Betensky, R. (2008). Simultaneous confidence intervals based on the percentile bootstrap approach. Computational Statistics Data Analysis, 52, 2158−2165.
Manterola, C., Conde, D. A., Colchero, F., Rivera, A., Huerta, E., & Pallares, E. (2011). El jaguar como elemento estratégico para la conservación. México D.F.: Conabio. Colección Corredor Biológico Mesoamericano México. Serie Acciones/ Número 8.
McCullagh, P., & Nelder, J. (1989). Generalized linear models. Florida: Chapman and Hall/ CRC.
McLoughlin, D. P., Morris, D. W., Fortin, D., Wal, E. V., & Contasti, A. L. (2010). Considering ecological dynamics in resource selection functions. Journal of Animal Ecology, 79, 4−12.
Miller, B., & Rabinowitz, A. (2002). ¿Por qué conservar el jaguar? In R. A. Medellín, C. Equihua, C. L. B. Chetkiewizc, P. G. Craeshaw, A. Rabinowitz, K. H. Redford et al. (Eds.), El jaguar en el nuevo milenio (pp. 303−316). México D.F.: Universidad Nacional Autónoma de México/ Wildlife Conservation Society.
Montiel-Ortega, S., Arias-Reyes, L., & Dickinson, F. (1999). La cacería tradicional en el norte de Yucatán: una práctica comunitaria. Revista de Geografía Agrícola, 29, 43−52.
Monroy-Vilchis, O., Zarco-González, M., Rodríguez-Soto, C., Soria-Díaz, L., & Urios, V. (2011). Fototrampeo de mamíferos en la Sierra Nanchititla, México: abundancia relativa y patrón de actividad. Revista Biologia Tropical, 59, 373−383.
Navarro, C. J., Ramolina, J. F., & Pérez, J. J. (2007). El jaguar en Yum Balam y el norte de Quintana Roo. In G. Ceballos, C. Chávez, R. List, & H. Zarza (Eds.), Conservación y manejo del jaguar en México: estudio de caso y perspectivas (pp. 123−132). México D.F.: Conabio/ Alianza WWF- Telcel/ Universidad Nacional Autónoma de México.
Núñez, R., Miller, B., & Lindzey, F. (2002). Ecología del jaguar en la reserva de la biosfera Chamela- Cuixmala, Jalisco, México. In R. A. Medellín, C. Equihua, C. L. B. Chetkiewizc, P. G. Craeshaw, A. Rabinowitz, K. H. Redford et al. (Eds.), El jaguar en el nuevo milenio (pp. 107−126). México D.F.: Universidad Nacional Autónoma de México/ Wildlife Conservation Society.
Rabinowitz, A. R., & Nottingham, B. G. (1986). Ecology and behavior of the jaguar (Panthera onca) in Belize, Central America. Journal of Zoology, 210,149−159.
Ridout, M. S., & Linkie, M. (2009). Estimating overlap of daily activity patterns from camera trap data. Journal of Agricultural, Biological, and Environmental Statistics, 14, 322−337.
Romero-Muñoz, A., Maffei, L., Cuéllar, E., & Noss, A. (2010). Temporal separation between jaguar and puma in the dry forest of southern Bolivia. Journal of Tropical Ecology, 26, 303−311.
Romeu, E. (1996). El Jaguar. Biodiversitas, 2, 1−5.
Sanderson, E. W., Chetkiewicz, C. L. B., Medellín, R. A., Rabinowitz, A., Redford, K., Robinson, J. G. et al.(2002). Un análisis del estado de conservación y distribución de los jaguares a través de sus áreas de distribución. In R. A. Medellín, C. Equihua, C. L. B. Chetkiewizc, P. G. Craeshaw, A. Rabinowitz, K. H. Redford et al. (Eds.), El jaguar en el nuevo milenio (pp. 551−600). México D.F.: Universidad Nacional Autónoma de México/ Wildlife Conservation Society.
Scognamillo, D., Maxit, I. E., Sunquist, M., & Polisar, J. (2003). Coexistence of jaguar (Panthera onca) and puma (Puma concolor) in a mosaic landscape in the Venezuela llanos. Journal of Zoology, 259, 269−279.
Schultz, G. (2003). Structure and diversity of the forests at the El Edén Ecological Reserve. In A. Gómez-Pompa, M. F. Allen, S. Fedick, & J. J. Jiménez-Osorio (Eds.), The lowland Maya area. Three millennia at the human-wildland interface (pp. 91−114). Binghamton, New York: The Haworth Press.
Siegel, S., & Castellan, N. J. (1988). Nonparametric statistics for the behavioral sciences. New York: McGraw-Hill.
Sileshi, G. (2008). The excess-zero problem in soil animal count data choice of appropriate models for statistical inference. Pedobiología, 52, 1−17.
Slagsvold, T. (1980). Habitat selection in birds: on the presence of other bird species with special regard to Turdus pilaris. Journal of Animal Ecology, 49, 523−536.
Zarza, H., Chávez, C., & Ceballos, G. (2007). Uso de hábitat del jaguar a escala regional en un paisaje dominado por actividades humanas en el sur de la península de Yucatán. In G. Ceballos, C. Chávez, R. List, & H. Zarza (Eds.), Conservación y manejo del jaguar en México: estudios de caso y perspectivas (pp. 101−110). México D.F.: Conabio/ Alianza WWF-Telcel/ UNAM.