Relationship Between the Abundance of the Endangered Volcano Rabbit Romerolagus Diazi and Vegetation Structure in the Sierra Chichinautzin Mountain Range, Mexico
Total Page:16
File Type:pdf, Size:1020Kb
Relationship between the abundance of the Endangered volcano rabbit Romerolagus diazi and vegetation structure in the Sierra Chichinautzin mountain range, Mexico A RELI R IZO-AGUILAR,JOSÉ A NTONIO G UERRERO,MIRCEA G. HIDALGO-MIHART and A LBERTO G ONZÁLEZ-ROMERO Abstract The volcano rabbit Romerolagus diazi is endemic communities (known locally as zacatonales) characterized to the mountains of central Mexico, where its habitat has mainly by Muhlenbergia spp., Festuca spp. and Jarava ichu been gradually destroyed to make way for agriculture, and associated with trees of the genera Pinus, Alnus and ranching and logging, and by forest fires. The volcano rabbit Quercus (Hoth et al., 1987; Velázquez, 1993). The species’ is categorized as Endangered on the IUCN Red List. habitat has been destroyed gradually for agriculture, We evaluated the relationship between the abundance of ranching and logging and by forest fires (Portales et al., the volcano rabbit and vegetation structure at a small scale 1997). Populations of the volcano rabbit are now at risk (0.25 ha). Using a general linear model we generated a set of and this species is categorized as Endangered, both by the 21 predictive models and proposed the best model as a Mexican government (SEMARNAT, 2010) and on the habitat quality index. Our results suggest that greater height IUCN Red List (AMCELA et al., 2008). and cover of bunchgrasses and the presence of a shrub layer The habitat of a species can be examined at various offer the volcano rabbit the best refuge from predators. The scales, ranging from the local microscale (microhabitat) habitat quality index and the limited available habitat to the regional and global macroscales (Partridge, 1978; documented in our survey indicate that the populations of Delcourt & Delcourt, 1988; Dunning et al., 1992; Morrison volcano rabbits in the study area are more threatened than et al., 2006). Heterogeneous landscapes are a common previously thought. As R. diazi is a habitat specialist it does feature in the geographical distribution of most species not have the option of moving to another habitat type. of mammals, even those with a restricted distribution (Feldhammer et al., 2007). This mosaic of environmental Keywords Endemic, habitat quality index, habitat special- conditions or habitat types is necessary to meet all of the ist, lagomorph, Mexico, pellet count, Romerolagus diazi, requirements for species survival (Hansson, 1979). At the volcano rabbit local scale habitat use can be studied by measuring the population density or abundance of focal species, assuming that the abundance of animals is higher in habitats that are Introduction more suitable for sustaining populations (Partridge, 1978; Duncan et al., 1997). This suggests a unimodal distribution he volcano rabbit Romerolagus diazi is endemic to of population density along a gradient of habitat types central Mexico, restricted to the Ajusco–Chichinautzin T (Velázquez, 1993). mountains and the Popocatépetl, Iztaccíhuatl, Pelado Although the current status of wild populations of and Tláloc volcanoes. Known locally as the teporingo or 2 R. diazi is critical, until now no information has been zacatuche, the species is distributed over c. 386 km reported on their size, structure or demography, and (Velázquez et al., 1996) at elevations of 2,800–4,250 m knowledge of the ecological relationships between this (Fa & Bell, 1990). The volcano rabbit is considered a habitat species and its habitat is scant. The available information on specialist and depends on wild subalpine bunchgrass habitat use is limited to the landscape scale (Velázquez, 1993; Velázquez & Heil, 1996). The only study carried out on the local scale was limited to the area of the Pelado volcano ARELI RIZO-AGUILAR (Corresponding author) and ALBERTO GONZÁLEZ-ROMERO (Fa et al., 1992). Therefore little is known about the habitat Instituto de Ecología, A. C., Carretera antigua a Coatepec No 351, El Haya, relationships of this species at local scales or about its C.P. 91070, Xalapa, Veracruz, Mexico. E-mail [email protected] peripheral populations (Velázquez et al., 1996). JOSÉ ANTONIO GUERRERO Facultad de Ciencias Biológicas, Universidad Autónoma del Estado de Morelos, Cuernavaca, Mexico Here we evaluate habitat characteristics and their relationship to the abundance of the volcano rabbit, and MIRCEA G. HIDALGO-MIHART División Académica de Ciencias Biológicas, Universidad Juárez Autónoma de Tabasco, Villahermosa, Mexico generate a habitat quality index based on this relationship. Received 15 June 2012. Revision requested 31 October 2012. We counted rabbit latrines as an indirect method of Accepted 18 June 2013. First published online 12 May 2014. determining population abundance (Palomares, 2001). Oryx, 2015, 49(2), 360–365 © 2014 Fauna & Flora International doi:10.1017/S0030605313000975 Downloaded from https://www.cambridge.org/core. IP address: 170.106.40.219, on 30 Sep 2021 at 15:37:20, subject to the Cambridge Core terms of use, available at https://www.cambridge.org/core/terms. https://doi.org/10.1017/S0030605313000975 The volcano rabbit 361 FIG. 1 Map of the Sierra Chichinautzin mountain range, with the locations of the 115 sampling sites. The rectangle on the inset shows the location of the main map in Mexico. Study area and has been used to estimate the population abundance of the volcano rabbit elsewhere in its range (Velázquez, 1994). We investigated populations of the volcano rabbit in the At each of the 64 sampling sites we surveyed one 50 × 50 m 1 Sierra Chichinautzin mountain range (Fig. ). This is an quadrat (0.25 ha) for fresh latrines. important area because the populations here may be We selected at random 22 of the 64 sampling points to ff e ectively isolated from other populations of the species describe the structure of the volcano rabbit’s habitat. We that inhabit the Pelado, Tláloc, Popocatépetl and measured seven vegetation attributes: tree height, cover, and Iztaccíhuatl volcanoes, which are considered to be its core diameter at breast height; shrub height and cover; and grass 1996 distribution areas (Velázquez et al., ). Part of the study height and cover. We calculated standard deviation for grass area includes a protected area, the Chichinautzin Biological height, shrub height and tree height, and for diameter at Corridor, which is located in the north-west of the state of breast height as a proxy for habitat complexity. We used a Morelos, on the southern edge of the Federal District clinometer to measure tree height and a telescopic rod to (Mexico City). measure shrub height, in metres to the nearest centimetre. To calculate tree and shrub cover we measured the crown Methods diameter of all trees and shrubs in the quadrat, using an electronic distance measuring tool. From the diameter the Using 1:20,000 aerial photographs from 1995 (National circular area of cover was calculated for each tree and shrub Institute of Statistics, Geography and Informatics, INEGI), and these values were used to calculate percentage cover. we selected 115 sampling points at random above 2,800 m We measured the height of 10 individuals of the two most altitude (using contour lines at the 1:50,000 scale). We abundant species of bunchgrass (Muhlenbergia macroura excluded urban areas (Tres Marías, Huitzilac, Coajomulco, and J. ichu) and the percentage cover was estimated visually Fierro del Toro and Parres). We visited all the sampling by the same person for each quadrat, thereby minimizing points to document habitat availability based on the inter-observer variation. presence of the grasses Muhlenbergia spp. and J. ichu. The variables used to describe a habitat are typically Sixty-four points had suitable habitat for R. diazi and interdependent and cannot be analysed individually 51 points did not (i.e. sites with fir, oak and pine forest, (Cooley & Lohnes, 1971). Therefore we conducted a secondary vegetation, and deforested areas used mainly for principal components analysis, which generates new axes growing oats and potatoes). from the linear combination of the original variables. During June–December 2008 we estimated the abun- The resulting principal components have the advantage dance of the volcano rabbit at the sampling sites by counting of being orthogonal and are considered independent the number of latrines. Following the method of Velázquez variables. For this analysis we used a correlation matrix (1994) a group of 30 or more pellets was defined as a latrine. because the habitat description variables differed in unit and This method is considered reliable for estimating the scale (James & McCulloch, 1990). We selected the principal population abundance of lagomorphs (Palomares, 2001) components with the highest eigenvalues, the sum of which Oryx, 2015, 49(2), 360–365 © 2014 Fauna & Flora International doi:10.1017/S0030605313000975 Downloaded from https://www.cambridge.org/core. IP address: 170.106.40.219, on 30 Sep 2021 at 15:37:20, subject to the Cambridge Core terms of use, available at https://www.cambridge.org/core/terms. https://doi.org/10.1017/S0030605313000975 362 A. Rizo-Aguilar et al. accounted for 90% of the variance. We used the principal and sites with 151–300 latrines were found at 3,060–3,150 m components as new descriptive variables of the habitat to (Fig. 2). construct a habitat quality index. We formulated a set of Of the 22 sites selected for habitat characterization, eight 21 predictive models using different combinations of the (36%) had 1–50 latrines, six (27%) had 51–150 latrines, five major components, which we analysed using generalized (23%) had 151–300 latrines and three (14%) had no latrines. linear models under a Poisson distribution. The indepen- We retained five principal components that together dent variables were the principal components and the accounted for 90% of the total variance in the habitat dependent variable was the abundance of the volcano rabbit data (Table 1). The first component (PC1) was negatively (as indicated by the number of latrines recorded at each associated with tree layer (cover, diameter at breast height, site).