Acacia woodlots, cattle and dung beetles (Coleoptera: Scarabaeinae) in a Mexican silvopastoral landscape

Matorrales de Acacia, vacas y escarabajos del estiércol (Coleoptera: Scarabaeinae) en un paisaje silvopastoril mexicano

Lucrecia Arellano¹*, Jorge L. León-Cortés², Gonzalo Halffter¹ and Jorge Montero³

¹Instituto de Ecología, A. C., Red de Ecoetología, Carretera Antigua a Coatepec Núm. 351, El Haya, Xalapa, 91070 Veracruz, México. *lucrecia.arellano@inecol.edu.mx

²El Colegio de la Frontera Sur, Departamento de Ecología y Sistemática Terrestre, Carretera Panamericana y Av. Periférico Sur s/n, San Cristóbal de las Casas, 29290 Chiapas, México.

Recibido: 28 agosto 2012
Aceptado: 14 diciembre 2012

Abstract

In landscapes where tropical dry forest was once the dominant vegetation type, traditional silvopastoral systems generate a range of natural and semi-natural habitats, namely fragments of secondary forest and fallow land in various stages of succession; essentially Acacia woodlots. This level of heterogeneity seems to favor the arrival and persistence of a large number of Scarabaeinae species. Changes in dung and carrion beetle assemblages were assessed in a silvopastoral landscape in Chiapas, Mexico by intensive sampling using baited pitfall traps. Total species richness decreased from late successional habitats to early successional habitats and average abundance per site was higher in the intermediate successional stages. The silvopastoral system studied is very dynamic; its components may remain constant or change over time, depending on management. These changing conditions lead to a variable local species richness and allow the permeability of certain species within the landscape. There was high species turnover among successional habitats, generating a gamma diversity similar to that of tropical dry forest landscapes and 2-fold higher than the highest recorded alpha diversity value. Silvopastoral systems can buffer the adverse effects of rapid expansion of open areas and the consequent reduction of tropical dry forest area generated by technified conventional systems.

Key words: fallow land, tropical dry forest, species diversity, scarab, livestock, successional habitats, pasture, Acacia hábitats.

Resumen

Palabras clave: acahual, bosque tropical caducifolio, diversidad de especies, escarabajos, ganado, hábitats sucesionales, potrero, hábitats de Acacia.

Introduction

To meet their needs, humans have used and transformed tropical landscapes through different forms of management (agricultural, livestock and forestry). To clearly understand the current status of conservation in anthropic landscapes, and to predict future trends, a deeper understanding of the biodiversity in actively managed rural landscapes is required (Ballam-Ballote and León-Cortés, 2010). In tropical dry landscapes, the conservation of species diversity depends not only on protecting forested areas, but also on managing the matrix in which they are embedded. Balanced management of the agrosilvopastoral matrix should influence the effects of productive activities on biodiversity and the ecological functions in the landscape (Harvey et al., 2006) and should sustain productivity while maintaining as much original ecosystem integrity as possible (Harvey et al., 2006; Harvey et al., 2008).

Tropical cattle production areas, where farmers use different proportions of fodder trees, shrubs and grasses, produce heterogeneous landscapes that may contribute to the conservation and maintenance of species diversity and ecological processes (Harvey et al., 2006; Giraldo et al., 2011; Rös et al., 2012). In contrast, conventional technified systems are those that feature a monoculture of grass with low or minimal tropical plant diversity: they are treeless and include irrigation systems to maintain the grazing. This type of management has resulted in a rapid expansion of open areas and the consequent reduction of original forest areas. Simplification in the structure and composition of vegetation, along with deforestation and species loss, causes changes in the dynamics of several ecological processes that are crucial to ecosystem functioning (Giraldo et al., 2011). In the anthropic landscapes of southern Mexico, where tropical dry forest was once the dominant vegetation, almost all cattle production activities involve extensive grazing (Gómez et al., 2002) conducted under a traditional agrosilvopastoral management scheme (Jiménez-Ferrer et al., 2007; Ramírez-Marcial et al., 2012). As a result, farms in tropical dry landscapes may include different proportions of grasses and trees: pastures with living fences and/or isolated trees or even forest fallows, where the land is put to a different use during the annual cycle (Nahed et al., 2009). However, some farms remain treeless as a result of management decisions taken by the landowners. The tree cover that is typically retained in silvopastoral landscapes may provide resources and habitats for a variety of organisms (Neita and Escobar, 2012) and may also increase the heterogeneity and functional connectivity of the landscape for certain animal species (Arellano et al. 2008a; Pineda-Diez de Bonilla et al., 2012). Vegetation dynamics and land use modifications can generate both predictable and unpredictable changes in the availability of certain habitats in a successional gradient (pastures, pastures with trees, living fences, secondary forest, managed Acacia habitats, tree and shrubby fallows). The species associated with these habitats could potentially undergo high rates of local extinction as a result of successional processes or disturbances, or may even disappear locally over the course of a few generations (Ballam-Ballote and León-Cortés, 2010); however, these changes could also enable the creation of new successional habitats for other species that depend on livestock activities (León Cortés et al., 2004), such as dung beetles (Scarabaeinae).

Dung and carrion beetles (Coleoptera, Scarabaeinae) are useful indicators of biodiversity in the tropics because their community structure responds rapidly to environmental changes, their biology is relatively well known and they are easy to sample (Favila and Halffter, 1997; McGeoch et al., 2002; Nichols et al., 2007). Our aim was to assess changes in richness, abundance, composition and community structure of dung beetles across a representative range of successional habitats within a silvopastoral landscape in southern Mexico. Our working hypothesis was that species richness and abundance in this landscape would be higher in sites of intermediate successional stage, because these types of heterogeneous habitats would offer conditions suitable for a variety of species with different habitat preferences. We would also expect that this successional stage would support higher proportions of habitat generalists than would be case in early and late successional stages.

Materials and methods

Study area. The landscape is located in the municipality of San Fernando, in Chiapas, Mexico (16°53'01"-16°47'57" N, 93°09'23"93°13'58" W) (see description in Arellano et al., 2008b). The transformation and fragmentation of this tropical dry forest landscape occurred at the turn of the nineteenth century, so the area has been fragmented for at least 150 years.

Site selection. The landscape elements were plotted on a 1:75 000 spatial orthophoto (Inegi, 2001), and digitized into a Geographical Information System (Arc View version 3.2 ©). To characterize the number of trees and shrubs, we randomly established 3 plots of 10 × 10 m at each site. Shrubby tree cover or high (> 2 m) was measured with a spherical densitometer. At each site, a fixed transect of 120 m was established in which we measured tree cover in cardinal point directions (i.e. N, S, E, W). We then selected 9 sites, an average distance of 4 km apart (minimum distance 610 m), within a silvopastoral landscape. The sites were selected to represent different successional stages and management and were classified to 3 successional stages with a cluster analysis, using average linkage as the joining method and Euclidian distance as the distance measure (Fig. 1): i) early (S1), which includes a pasture with no trees or shrubs. ii) Intermediate (S2), including 2 areas of young managed Acacia shrubs (2 and 4 years old, respectively) and a pasture with Acacia pennatula, and iii) late (S3), with a tropical subdeciduous forest, a tropical deciduous forest, a mature tropical deciduous forest fallow, a ravine with tropical dry forest and a mature Acacia forest (Table 1).

Dung beetle sampling. During the period July-August of 2003 and 2004, beetles were sampled along 18 transects (150 m in length), with 6 baited pitfall traps/site (9 transects/year/site). Bait type (150 g of human excrement, cow dung and squid) was alternated so the distance between traps with the same type of bait in each transect was 75 m. Traps remained functional for 48 h and total sampling effort was 2 268 hours.

Data analysis. The final slope of the species accumulation curves was used to determine survey completeness (Hortal and Lobo, 2005). This choice had the advantage of identifying a cut-off value to determine the rate of finding new species in additional samples, which suggests a comparative measure of the survey completeness. As the inventory is being completed, new records of species become rare, so that the slope of the curve decreases. In addition to a decrease in this slope, the sampling effort required to add a significant number of species to the inventory increases, and therefore the balance between extra effort and number of new species becomes less favorable. Therefore, species that are missing from the inventory may be locally rare species that can nonetheless still be found (see Moreno and Halffter, 2000).

The Michaelis-Menten (MM) equation was used as a reliable estimator for dung beetles, since MM commonly used in diversity inventory assessment and is an adequate function for a relatively small-intermediate number of samples (Colwell and Coddington, 1994; Colwell, 2009). For this equation, and using number of individuals as a unit of sampling effort, the survey can be considered sufficiently reliable despite being incomplete, even when the slope becomes approximately < 0.1. Using individuals, one could establish the slope cut-off at 0.05 or 0.01, which would imply a new species every 20 or 100 individuals, respectively. However, if the final slope is much higher than the 0.1 cut-off point, it may be that we are using such a low number of (i) sampling places, (ii) sampling methodologies, or (iii) sampling dates, that we are unable to capture reliable inventories.

Kruskal Wallis tests were used to examine changes in the median and range of richness and abundance per site, while Mann Whitney tests were used to examine changes in the median and range of richness and abundance between early-intermediate successional stage and late successional stage. Early and intermediate stages were grouped together because the early stage included only 1 site. The influence of vegetation variables on species richness and abundance in every site (alpha diversity) was explored using Spearman correlations.

Rank-abundance curves, based on the number of individuals per species, were used to compare the changes in community structure for each habitat on the successional gradient. The equitability values (J') per site were obtained. Ancova analysis was used to compare abundances among sites from rank abundance curves.

To identify indicator species for every successional stage in the silvopastoral landscape, we used INDVAL (Dufrene and Legendre, 1997). Changes in the composition of the dung and carrion beetle communities in different habitats were analyzed according to: a) food relocation method: the proportion of tunnelers, endocoprids and rollers, b) food preference: the proportion of generalists (species for which less than 80% of total individuals were caught in copro- or necro-baited traps) to specialists (species for which more than 80% of the individuals were collected in either copro- or necro-baited traps), c) daily activity: the proportion of nocturnal to diurnal species, d) beetle size: the proportions of large (18-28 mm), medium-sized (9-17.99 mm) and small (< 9 mm) species (Halffter and Arellano, 2002), and e) habitat preference: the proportion of habitat generalists (species for which the individuals were collected in more than one specific habitat) to habitat specialists (species for which more than 80% of total individuals were caught in a specific habitat). The R × C test of independence, using the G-test, was performed to compare the proportion of individuals in each functional trait per habitat (Sokal and Rohlf, 1995).

A Venn diagram was used to show the total number of species, shared species and exclusive species per successional stage. The spatial arrangement of sites based on species composition was analyzed with a non-metric multidimensional scaling ordination procedure (NMDS). Bray-Curtis dissimilarity indices were used for NMDS ordination, since they robustly identify relationships between site dissimilarity and ecological distance among species (Faith et al., 1987). P values were based on 1 000 permutations.

To analyze differences in species composition (using abundances) for successional stages, we applied for Adonis function using R, a multivariate ANOVA based on dissimilarities, and the betadisper function for the homogeneity of groups and beta diversity (Oksanen, 2012) based on a dissimilarity matrix (Bray Curtis index), for successional stages (grouping early and intermediate stages). A test was used to measure differences in the multivariate homoscedasticity among the sites of every successional stage; this test is analogous to a Levene test, for homoscedasticity and is used to measure the dispersal in species composition at every successional stage.

Results

We recorded 1 230 individuals, representing 33 species across the entire landscape (Table 2). In general, the total number of species decreased from the late successional to early successional habitats: there were 28 species in the late successional stages, 18 in the intermediate successional stages and 8 in the early successional stages. Similarly, the average number of species/transects was higher in S3 (with the exception of the ravine) than in the other successional stages (Fig. 2a). The highest abundances were obtained in TDF and in the pasture with Acacia (Fig. 2b). The average number of individuals per site was higher in S2 than in the other successional stages.

The final slope of the species accumulation curve was 0.0058 for the entire landscape; i.e. a new species was added for every 171 individuals collected (Fig. 3). For early (S1), intermediate (S2) and late stage (S3) successional habitats, a new species was collected every 182, 85 and 88 individuals, respectively. According to the Michaelis-Menten estimator (MM), at least 94% of the potential species richness for the entire landscape was collected, representing 88% of species in the early stage, 71% in the intermediate stage, and 82% in the late stage.

No significant differences were detected in local richness (Mann Whitney ( z= -1.49, p= 0.14) or in abundance (Mann Whitney ( z= 0.245, p= 0.80) between early-intermediate successional stages and S3. No significant differences were detected in local abundance (Kruskal-Wallis, H (8, n= 9) = 8, p= 0.43) and in local richness (Kruskal-Wallis, H (8, n= 9) = 7.73, p= 0.46), among sites. The minimum and maximum number of species per site was 8 and 17 species, respectively. We did not find significant correlations for the values of vegetation features among successional stages (Table 3).

The equitability value was 0.67 for S1 while this value in S2 increased with the increasing structural complexity (Fig. 4). In S3, the increase of equitability was apparently related to the increased relative humidity in the wooded sites. These trends might reflect dominant-diversity curves that show the highest slopes in sites with low structural complexity and in wooded sites with low relative humidity (Fig. 5). Significant differences in abundances were detected among the species rank of successional stages (ANCOVA, F= 35.05, p< 0.0001). Canthon (Glaphyrocanthon) delgadoi was distributed in almost all sites and was highly abundant in all successional stages. In early successional stage and in intermediate successional stage, C. indigaceus chiapas was the dominant species, followed by C. cyanellus cyanellus and C. (G.) delgadoi. In late stages, the most abundant species in the forest fragments were Uroxys deavilai and C. (G.) delgadoi, while in the fallows Onthophagus igualensis and Phanaeus endymion were the most abundant (Fig. 5).

According to the INDVAL analysis, the species with the highest indicator values for the early and intermediate successional stages (pasture and pasture with Acacia) was Canthon indigaceus chiapas (IV= 99.4, p= 0.036), a diurnal, coprophagous, median roller and Copris lugubris (IV=88.2, p= 0.05) a nocturnal, coprophagous, median tunneler for the intermediate successional stage; while for the late successional stage Deltochilum scabriusculum (IV= 100, p= 0.006) a nocturnal, coprophagous, large roller and Phanaeus endymion (IV=85.3, p= 0.007) a diurnal, generalist, median tunneler had the highest values. There were no significant values for the other species and their indicator values (IV) were lower than 75.

The proportions of the species categories for every functional trait were independent of the tree cover at each site, with the exception of habitat preferences that were associated with successional stage (x²= 24.11, N= 57, d.f.= 6, p= 0.05). In pastures, there were no forest species or mature Acacia forest species. In tropical deciduous forest, pasture species represented 7% of all the species captured, while young managed Acacia bush (2-year-old) contained 9%, and mature Acacia forests had 12%.

A large number of species were lost between successional habitats: 22 species were lost from S3 to S1, while 12 species were lost from S3 to S2 and 12 species were lost from S2 to S1. Figure 6 shows the total number of species and shared and exclusive species per successional stage. In S1, Euoniticellus intermedius and Scatimus ovatus were exclusive species, in S2 Eurysternus mexicanus and Agamopus lampros were exclusive and in S3, there were 12 exclusive species; e.g. Canthon (Glaphyrocanthon) lucreciae, Copris laeviceps, Phanaeus wagneri, Canthon femoralis, Pedaridium maya and Onthophagus violetae.

The Multivariate ANOVA based on dissimilarities detected significant differences between early-intermediate successional stages and late succesional stages (Pseudo-F= 2.29, p (perm)= 0.038). In Axis 1 of the NMDS, sites included in S3 are separated from the S2 and S1 sites (Fig. 7). This distribution of species is explained by arboreal coverage (DC). Stress type 1 (weak ties) was 0.06.

Discussion

The studied landscape has a variegated structure (see Rös et al., 2012) and dynamic structure, the latter caused by the action of cattle and frequent fires (Table 4). It is formed by remnants of tropical dry forest, Acacia woodlots and grasslands. At landscape level, the diversity of species of Scarabaeinae (33 species) is similar to that of other landscapes where the original vegetation was tropical dry forest (Arellano and Halffter, 2003; Hernández et al., 2003; Harvey et al., 2006). The heterogeneous nature of the silvopastoral landscape and the dynamism of many of its elements allow the entry of foreign species that had abundant populations in other landscapes; e.g. Agamopus lampros, Canthon femoralis, Euoniticellus intermedius, Onthophagus violetae, Pedaridium maya, Phanaeus wagneri, P. demon, etc.

In the pasture land with Acacia shrubs, species abundance was higher than in other sites, and was similar to that recorded in TDF. In contrast, other studies report higher abundance in early successional stages than in late successional stages (Manrique, 2009; Ballam-Ballote and León-Cortés, 2010). In Tizimín, Yucatán, Mexico, Basto-Estrella et al. (2012) found the highest abundance on livestock farms (> 90 000 individuals), perhaps because heliophiles and open site species (Digitonthophagus gazella, Canthon indigaceus chevrolati, Onthophagus landolti) were favored by the continuity and large expanses of open livestock areas in that region.

Contrary to our hypothesis, and the results of other studies with insects (Manrique, 2009; Ballam-Ballote and León-Cortés, 2010), total species richness was not higher in the intermediate successional stages (the ravine and young managed Acacia bush), but instead decreased from late to early successional habitats. For instance, 22 species were lost from S3 to S1, while only 2 species were gained. Structural complexity seems to be an important factor for beetle diversity. Ravines, as in other tropical Mexican landscapes, are sites that allow species with different environmental requirements and biogeographic affinities to colonize silvopastoral landscapes.

Within the landscape, alpha diversity values vary between sites (especially among communities with Acacia) due to the presence of species represented by few individuals. Gamma diversity of the landscape increases in relation to the dry tropical forests, due to heterogeneity and its influence on beta diversity. Intriguingly, the presence of "tourist" species, their number and frequency is highly variable between different successional stages. Therefore, the high number of tourist species is associated with livestock management in the landscape.

When food availability is high, a habitat becomes suitable for larger dung beetle species (Hanski and Cambefort, 1991). In the landscape we studied, small species dominated due in all likelihood to the low availability of dung and the conditions in certain habitats such as the managed young Acacia bush (high densities of trees with many thorns, dry soil, etc.). This pattern is consistent with previous assessments that suggest that open-habitat communities are characterized by highly abundant small species (e.g. Nichols et al., 2007; Caballero et al., 2009; Caballero and León-Cortés, 2012). Throughout this landscape, species that are characteristic of wooded habitats dominate, while those typical of pastures account for a low percentage of the beetle assemblages. This silvopastoral landscape has a high proportion of generalist (53%) and forest species (22%) that have apparently been favored by the landscape structure and for which the functional connectivity is relatively high.

We conclude that the traditional management of tropical dry landscapes has created habitat conditions suitable for the persistence of a variety of species. However, we are concerned that landowners may decide to switch from traditional practices to conventional cattle production techniques with low plant diversity and high dependence on chemical fertilizers and herbicides. Such a transformation would be likely to cause a dramatic change in the abundance of species within the landscape, the local extinction of those species favored by landscape heterogeneity and, in the long term, landscape-scale extinctions.

This silvopastoral landscape allows for a high diversity of beetles while maintaining a very profitable economic activity (livestock production). Acacia pennatula is a native species that increases the profitability of livestock due to its multiple uses, and might -in the long term- allow a potential reestablishment of tropical deciduous forest habitats. A well managed population of Acacia makes it possible to intensify land use while reducing pressure on the forest remnants. Specifically, the use of Acacia farnesiana has been recommended to control erosion and improve soil quality. These trees can reduce peak daytime temperatures near the soil surface, and modify soil physical and chemical conditions, increasing fertility and moisture retention (Cházaro, 1977).

Finally, the persistence of insect species with restricted distributions such as Baronia brevicornis (Lepidoptera: Papilionidae), which is associated with silvopastoral landscapes such as the one we studied (León-Cortés et al., 2004), underlines the importance of maintaining traditionally managed systems in tropical Mesoamerican landscapes. There have been perhaps 10 000 or more insect generations that have been able to persist and interact in these man-made landscape elements, and this alone provides good reason to reconcile landscape management with traditional practices.

Acknowledgments

Manuel Girón, Arcángel Molina, José G. Alonso and José Luis Díaz-Valdivieso helped in the field. Leonardo Delgado, David Edmonds, Bert Kohlmann, Fernando Vaz de Mello and Mario Zunino confirmed our species identifications. Special thanks to Vicente Jiménez and Reinalda Hernández for allowing us to conduct this study on their land. The first author is grateful to SEPCONACYT for a doctoral grant (175987) and El Colegio de la Frontera Sur in San Cristóbal de las Casas for funding and the use of their instalations. We thank the Instituto de Ecología, A. C., for laboratory facilities. Bruce Ferguson, Jorge Macías and José Nahed made helpful comments and suggestions on the manuscript. Keith MacMillan corrected the English.

Literature cited

Arellano, L. and G. Halffter. 2003. Gamma diversity: derived from and a determinant of alpha diversity and beta diversity. An analysis of three tropical landscapes. Acta Zoológica Mexicana (n.s) 90:27-76.         [ Links ]

Arellano, L., J. León-Cortés and O. Ovaskainen. 2008a. Patterns of abundance and movement in relation to landscape structure - a study of a common scarab (Canthon cyanellus cyanellus) in Southern Mexico. Landscape Ecology 23:69-78.         [ Links ]

Arellano, L., J. León-Cortés and G. Halffter. 2008b. Response of dung beetle assemblages and their conservation in remnant natural and modified habitats in southern Mexico. Insect Conservation and Diversity 1:253-262.         [ Links ]

Balam-Ballote, Y. R. and J. L. León-Cortés. 2010. Forest management and biodiversity: a study of an indicator insect group in Southern Mexico. Interciencia 35:526-533.         [ Links ]

Basto-Estrella, G., R. Rodríguez-Vivas, H. Delfin-González and E. Reyes-Novelo.2012. Escarabajos estercoleros (Coleoptera: Scarabaeidae: Scarabaeinae) de ranchos ganaderos de Yucatán, México. Revista Mexicana de Biodiversidad 83:380-386.         [ Links ]

Caballero, U., J. L. León-Cortés and A. Morón Ríos. 2009. Response of rove beetles (Staphylinidae) to various habitat types and change in Southern Mexico. Journal of Insect Conservation 13:67-75.         [ Links ]

Caballero, U. and J. L. León-Cortés. 2012. High diversity beetle assemblages attracted to carrion and dung in threatened tropical oak forests in Southern Mexico. Journal of Insect Conservation 16:537-547.         [ Links ]

Cházaro, B. M. 1977. El Huizache, Acacia pennatula (Sclecht. and Cham.) Benth. Una invasora del centro de Veracruz. Biotica 2:1-18.         [ Links ]

Colwell, R. K. and J. A. Coddington. 1994. Estimating terrestrial biodiversity through extrapolation. In Biodiversity: measurement and estimation, D. L. Hawksworth (ed.). Chapman and Hall. London. p. 75-79.         [ Links ]

Colwell, R. K. 2009. EstimateS: Statistical estimation of species richness and shared species from samples. Version 8.2. User's Guide and application. http://purl.oclc.org/estimates; last access: 14.X.2009.         [ Links ]

Dufrene, M. and P. Legendre. 1997. Species assemblages and indicator species: the need for a flexible asymmetrical approach. Ecological Monographs 67:345-366.         [ Links ]

Faith, D., P. Minchin and L. Belbin.1987. Compositional dissimilarity as a robust measure of ecological distance. Vegetatio 69:57-68.         [ Links ]

Favila, M. E. and G. Halffter. 1997. The use of indicator groups for measuring biodiversity as related to community structure and function. Acta Zoológica Mexicana (n.s.) 72:1-25.         [ Links ]

Ferguson, B. G., T. Alemán-Santillán, G. Jiménez-Ferrer, J. Nahed-Toral, I. Carmona- Muñoz and H. Gómez-Castro. 2007. Desarrollo regional. In Ganadería, desarrollo y ambiente: una visión para Chiapas, T. Alemán-Santillán, B. G. Ferguson and F. J. Medina-Jonapá (eds.). El Colegio de La Frontera Sur, Fundación Produce Chiapas. Tapachula, Chiapas, México. 86-95 p.         [ Links ]

Giraldo, C., F. Escobar, J. N. D. Chara and Z. Calle. 2011. The adoption of silvopastoral systems promotes the recovery of ecological processes regulated by dung beetles in the Colombian Andes. Insect Conservation and Diversity 4:115-122.         [ Links ]

Gómez, C. H., M. A. Tewolde and J. Nahed. 2002. Análisis de los sistemas ganaderos de doble propósito en el centro de Chiapas, México. Latin-American Archives of Animal Production 10:175-183.         [ Links ]

Greenberg, R., P. Bichier and J. Sterling. 1997. Acacia, cattle and migratory birds in southeastern Mexico. Biological Conservation 80:235-247.         [ Links ]

Gutiérrez, J. R. and J. Armesto. 1981. El rol del ganado en la dispersión de las semillas de Acacia caven (Leguminosae). Serie: Ciencia e Investigación Agraria 8:3-8.         [ Links ]

Halffter, G. and L. Arellano. 2002. Response of dung beetle diversity to human-induced changes in a tropical landscape. Biotropica 34:144-154.         [ Links ]

Hanski, I. and I. Cambefort. 1991. Dung Beetles Ecology. Princeton University Press, Princeton, N. Y. 481 p.         [ Links ]

Harvey, C. A., A. Medina, D. M. Sánchez and S. Vílchez. 2006. Patterns of animal diversity in different forms of tree cover in agricultural landscapes. Ecological Applications 16:986-999.         [ Links ]

Harvey, C. A., K. Oliver, R. Chazdon, B. G. Ferguson, B. Finegan, D. M. Griffith, M. Martínez-Ramos, H. Morales, R. Nigh, L. Soto-Pinto, M. van Breugel and M. Wishnie. 2008. Integrating Agricultural Landscapes with Biodiversity Conservation in the Mesoamerican Hotspot. Conservation Biology 22:8-15.         [ Links ]

Hernández, B., J. M. Maes, C. A. Harvey, S. Vílchez, A. Medina and D. Sánchez. 2003. Abundancia y diversidad de escarabajos coprófagos y mariposas diurnas en un paisaje ganadero en el departamento de Rivas, Nicaragua. Agroforestería en las Américas 10:93-102.         [ Links ]

Hortal, J. and J. M. Lobo. 2005. An ED-based protocol for the optimal sampling of biodiversity. Biodiversity and Conservation, 14:2913-2947.         [ Links ]

Inegi (Instituto nacional de estadística, geografía e informática). 2001. Ortofotos E15C59D, E15C59E, escala 1: 75,000. Sistema Nacional Estadístico y de Información Geográfica. México.         [ Links ]

Jiménez-Ferrer, G., J. Nahed and L. Soto (eds.). 2007. Agroforestería pecuaria en Chiapas, México. Ganadería y Ambiente. El Colegio de la Frontera Sur. San Cristóbal de las Casas, Chiapas, México. 76 p.         [ Links ]

León-Cortés, J. L., F. Pérez-Espinosa, L. Marín and A. Molina-Martínez. 2004. Complex habitat requirements and conservation need of the only extant Baroniinae swallowtail butterfly. Animal Conservation 7:1-10.         [ Links ]

Manrique, C. 2009. Composición y abundancia de la comunidad de coleópteros asociada a la vegetación en un gradiente sucesional del Bosque Tropical Caducifolio, en la región de Chamela-Cuixmala, municipio de la Huerta, Jalisco. Tesis, Facultad de Biología, Universidad Michoacana de San Nicolás de Hidalgo. Morelia, Michoacán, México. 64 p.         [ Links ]

McGeoch, M., B. J. van Rensburg and A. Botes. 2002. The verification and application of bioindicators: a case study of dung beetles in a savanna ecosystem. Journal of Applied Ecology 39:661-672.         [ Links ]

Moreno, C. E. and G. Halffter. 2000. Assessing the completeness of bat biodiversity inventories using species accumulation curves. Journal of Applied Ecology 37:149-158.         [ Links ]

Nahed, T. J., P. J. Calderón, J. R. Aguilar, B. Sánchez-Muñoz, J. L. Ruiz-Rojas, Y. Mena, J. M. Castel, F. A. Ruiz, G. Jiménez, J. López-Méndez, G. Sánchez-Moreno and B. Salvatierra. 2009. Aproximación de los sistemas agrosilvopastoriles de tres microrregiones de Chiapas, México, al modelo de producción orgánica. Avances en Investigación Agropecuaria 13:45-58.         [ Links ]

Neita, J. C. and F. Escobar. 2012. The potential value of agroforestry to dung beetle diversity in the wet tropical forests of the Pacific lowlands of Colombia. Agroforestry Systems 85:121-131.         [ Links ]

Nichols, E., T. Larsen, S. Spector, A. L.V. Davis, F. Escobar, M. E. Favila and K. Vulinec. 2007. Global dung beetle response to tropical forest modification and fragmentation: a quantitative literature review and meta-analysis. Biological Conservation 137:1-19.         [ Links ]

Oksanen, J. 2012. Multivariate Analysis of Ecological Communities in R: vegan tutorial. http://CRAN.R-project.org/package=vegan.         [ Links ]

Pineda-Diez E., J. L. León-Cortés and J. L. Rangel-Salazar. 2012. Diversity of bird feeding guilds in relation to habitat heterogeneity and land-use cover in a human-modified landscape in southern Mexico. Journal of Tropical Ecology 28:369-376.         [ Links ]

Ramírez-Marcial, N., M. L. Rueda-Pérez, B. Ferguson and G. Jiménez-Ferrer. 2012. Caracterización del sistema agrosilvopastoril en la depresión central de Chiapas. Avances en Investigación Agropecuaria 16:7-22.         [ Links ]

Rös, M., F. Escobar and G. Halffter. 2012. How dung beetles respond to a human-modified variegated landscape in Mexican cloud forest: a study of biodiversity integrating ecological and biogeographical perspectives. Diversity and Distributions 18:377-384.         [ Links ]

Sokal, R. R. and F. J. Rohlf. 1995. Biometry. 3rd ed. W. H. Freeman and Company, New York. 887 p.         [ Links ]